Skip to main content


TBCRC 002: a phase II, randomized, open-label trial of preoperative letrozole with or without bevacizumab in postmenopausal women with newly diagnosed stage 2/3 hormone receptor-positive and HER2-negative breast cancer



In preclinical studies, the expression of vascular endothelial growth factor (VEGF) in hormone receptor-positive breast cancer is associated with estrogen-independent tumor growth and resistance to endocrine therapies. This study investigated whether the addition of bevacizumab, a monoclonal antibody against VEGF, to letrozole enhanced the antitumor activity of the letrozole in the preoperative setting.


Postmenopausal women with newly diagnosed stage 2 or 3 estrogen and/or progesterone receptor-positive, HER2-negative breast cancer were randomly assigned (2:1) between letrozole 2.5 mg PO daily plus bevacizumab 15 mg/kg IV every 3 weeks (Let/Bev) and letrozole 2.5 mg PO daily (Let) for 24 weeks prior to definitive surgery. Primary objective was within-arm pathologic complete remission (pCR) rate. Secondary objectives were safety, objective response, and downstaging rate.


Seventy-five patients were randomized (Let/Bev n = 50, Let n = 25). Of the 45 patients evaluable for pathological response in the Let/Bev arm, 5 (11%; 95% CI, 3.7–24.1%) achieved pCR and 4 (9%; 95% CI, 2.5–21.2%) had microscopic residual disease; no pCRs or microscopic residual disease was seen in the Let arm (0%; 95% CI, 0–14.2%). The rates of downstaging were 44.4% (95% CI, 29.6–60.0%) and 37.5% (95% CI, 18.8–59.4%) in the Let/Bev and Let arms, respectively. Adverse events typically associated with letrozole (hot flashes, arthralgias, fatigue, myalgias) occurred in similar frequencies in the two arms. Hypertension, headache, and proteinuria were seen exclusively in the Let/Bev arm. The rates of grade 3 and 4 adverse events and discontinuation due to adverse events were 18% vs 8% and 16% vs none in the Let/Bev and Let arms, respectively. A small RNA-based classifier predictive of response to preoperative Let/Bev was developed and confirmed on an independent cohort.


In the preoperative setting, the addition of bevacizumab to letrozole was associated with a pCR rate of 11%; no pCR was seen with letrozole alone. There was additive toxicity with the incorporation of bevacizumab. Responses to Let/Bev can be predicted from the levels of 5 small RNAs in a pretreatment biopsy.

Trial registration

This trial is registered with (Identifier: NCT00161291), first posted on September 12, 2005, and is completed.


Angiogenesis is an integral step in tumor progression [1]. Specifically in hormone receptor-positive breast cancer, enhanced angiogenesis accelerated tumor growth [2], and overexpression of the vascular endothelial growth factor (VEGF) either by the tumor [3, 4] or stromal cells [5] was associated with hormone-independent growth and resistance to endocrine therapies. An interplay between estrogens and angiogenic activity exists [6, 7], and VEGF is one of the genes whose expression is modulated by estrogens. Indeed, the presence of functional estrogen response elements in the VEGF gene [8,9,10] underpins the transcriptional activation of VEGF by estrogens [11]. Tumor-released VEGF recruits stromal cells and promotes a desmoplastic microenvironment; stromal cells, in turn, provide mitogenic and angiogenic growth factors stimulating both tumor and stromal cell growth [4]. VEGF secreted by stromal cells and acting cooperatively with other factors can substitute for estrogens and foster hormone-independent growth of luminal tumors [5]. At the clinical level, in hormone receptor-positive breast cancer, elevated intratumoral levels of VEGF have been associated with suboptimal responses to hormonal therapies and poorer clinical outcomes [12,13,14] lending support to the hypothesis that VEGF and angiogenesis may contribute to resistance to endocrine therapies.

These preclinical studies set the stage for a pilot, single-institution, single-arm study of preoperative letrozole in combination with bevacizumab in postmenopausal women with hormone receptor-positive breast cancer [15]. In that study (n = 25), preoperative letrozole and bevacizumab for 24 weeks resulted in a pCR rate of 12% (n = 3), while the overall objective response rate was 68% (n = 17). The rate of downstaging to stage 0 or stage 1 was 32% [15]. The treatment was overall well tolerated with an 8% discontinuation rate (n = 2) due to adverse events.

The purpose of the current study was to evaluate the benefit of neoadjuvant/preoperative bevacizumab in combination with letrozole in postmenopausal women with hormone receptor-positive breast cancer. The principal hypothesis of the study was that concurrent inhibition of estrogen receptor (ER)-mediated signaling and VEGF-mediated stromal activation and angiogenesis would be more effective than ER blockade alone. The primary objective of the study was to determine the pCR rate with neoadjuvant letrozole for 24 weeks with or without bevacizumab; secondary objectives included assessments of toxicity, overall objective response, and downstaging rates in the two arms. The utility of circulating tumor and endothelial cells as biomarkers of response and pretreatment tumor genomic profiles associated with response or resistance to protocol therapy were explored.


Study design

This was a randomized selection (“pick-the-winner”), open-label, two-treatment arm, phase II study conducted through the Translational Breast Cancer Research Consortium across nine cancer centers in the USA. Bevacizumab was provided by Genentech. This study is registered on the clinical trial website of the US National Cancer Institute (

Patient population

Eligible patients were postmenopausal women, 50 years or older, with newly diagnosed and pathologically confirmed invasive stage 2 and 3 breast cancer (T2, T3, T4a-c, N0-2, and M0), ER- and/or PR-positive, Her2/neu-negative (0 or 1 by immunohistochemistry or non-amplified by in situ hybridization; a 2+ score by immunohistochemistry had to be confirmed as non-amplified by in situ hybridization) who were candidates for mastectomy or breast conserving surgery. Eligible patients had to have measurable disease by mammogram and/or ultrasound; in special cases, dedicated breast magnetic resonance imaging (MRI) could be used. Patients were required to demonstrate controlled blood pressure (< 150/90 mmHg) and a normal ejection fraction with a multigated acquisition (MUGA) scan or echocardiogram at baseline.

Patients with inflammatory (T4d) or metastatic breast cancer were excluded. Other exclusion criteria consisted of uncontrolled endocrine or cardiac disease, bilateral breast cancer, major surgical procedure within 28 days of initiation of therapy, history of abdominal fistula-gastrointestinal perforation-intra-abdominal abscess, or other malignant diseases. Co-administration of other cancer treatments was not allowed, and hormonal replacement therapy had to be discontinued at least 2 weeks before the initiation of study therapy.

Patients with positive axillary lymph nodes, either by clinical examination or by fine needle aspiration, were treated as having at least N1 disease and were not required to undergo sentinel lymph node biopsy prior to initiation of neoadjuvant therapy. Axillary lymph node dissection was performed in those patients at the time of definitive surgery. Sentinel lymph node biopsy before initiation of neoadjuvant therapy for non-palpable nodes was performed at the discretion of the treating physician. All patients gave informed consent. The protocol was reviewed and approved by each participating institution, and the study followed the Declaration of Helsinki and good clinical practice guidelines.

Treatment and dose modification

Patients who met the eligibility criteria and signed informed consent were randomly assigned, 2:1, to daily treatment with letrozole 2.5 mg PO and bevacizumab 15 mg/kg IV every 3 weeks (investigational arm, Let/Bev) or daily treatment with letrozole 2.5 mg PO alone (control arm, Let). The duration of each cycle was 3 weeks.

In the absence of treatment delays due to adverse events, treatment continued until one of the following criteria was met: 24 weeks of active therapy or 9 administrations of bevacizumab, disease progression, intercurrent illness that prevented further administration of treatment, any toxicity that resulted in a treatment delay of > 3 weeks, unacceptable grade 3 or 4 non-hematologic adverse event(s), consent withdrawal, and general or specific changes in the patient’s condition that rendered unacceptable for further continuation of treatment in the judgment of the investigator. No dose modifications were allowed for both bevacizumab and letrozole. The maximum allowable length of treatment interruption for bevacizumab was 6 weeks due to the long elimination half-life of the agent. An interval of at least 4 weeks between the last administration of bevacizumab and definitive surgery was required while letrozole was continued until the day before surgery.

Assessments for safety and efficacy

Patients were monitored for safety and tolerability throughout and after completion of protocol treatment. Toxicity was assessed using the National Cancer Institute Common Terminology Criteria for Adverse Events version 3.0. Physicians’ visits and physical examinations were performed every 6 weeks; laboratory evaluations were performed every 3 weeks in the investigational arm (including a urine protein/creatinine ratio or urine dipstick) and every 6 weeks in the control arm. Study procedures included evaluation of cardiac function at screening and week 24 with echocardiogram or MUGA scan.

At study enrollment, the surgeons were required to state the planned definitive surgery that the patient would have required without preoperative treatment. The type of definitive surgery that the patient underwent after preoperative protocol therapy was recorded.

At 6, 12, ad 18 weeks, patients underwent evaluation of response by physical examination and breast ultrasound. Patients with objective response or stable disease (SD) according to the Response Evaluation Criteria in Solid Tumors (RECIST [16]) continued the same regimen for a maximum of 24 weeks or until progression was observed. Patients with progressive disease (PD) at any time were removed from the protocol and were treated at the investigator’s discretion. At week 24, the evaluation before surgery included a physical evaluation, a breast ultrasound, and a mammogram. Post-surgical adjuvant therapy was administered at the investigator’s discretion.

pCR was defined as the absence of any residual invasive cancer in the resected breast specimen and all sampled ipsilateral lymph nodes (i.e., ypT0/is ypN0) adopting the definition proposed by the FDA [17]; microscopic residual disease was defined as < 5 mm residual invasive cancer in the breast (ypT1mi or ypT1a) and isolated or no tumor cell clusters in the ipsilateral axillary lymph nodes (ypN1mi or ypN0).

Imaging responses were assessed according to RECIST [16] by breast ultrasound or breast MRI in selected cases (whichever modality was consistently obtained throughout the treatment period and most accurately assessed the status of the tumor). Comparisons were made with the same imaging modality in every patient. Pathologic responses were assessed by comparing the maximum cumulative diameter of the target lesion(s) at the time of diagnosis as assessed by imaging studies with the size of the tumor in the final surgical pathology.

Circulating tumor cells

Levels of circulating tumor cells (CTCs) and circulating endothelial cells (CECs) in the peripheral blood prior to treatment (baseline) and during treatment at weeks 6, 18, and 24 were assessed. The number of CTCs in 7.5 mL of blood was measured using the CellSearch system (Veridex LLC, Raritan, NJ, USA) [18]. CTCs were defined as EPCAM-positive nucleated cells expressing cytokeratins but not the leukocyte common antigen CD45 [19]. In a few cases, e.g., when the available blood volume was below 7.5 mL, we performed CTC enumeration using an EPCAM-based immumomagnetic enrichment and flow cytometry procedure [20]. In this assay, CTCs were defined as nucleated EPCAM-positive and CD45-negative. In order to combine the data from the two assays, we reported the CTC values as CTC per mL of blood.

The number of CECs in 50 μL of blood was measured using a four-color flow cytometric assay performed on FACS Calibur (BD Biosciences, San Jose, CA, USA) [21]. CECs were defined as nucleated, CD34-positive, CD31-positive, and CD45-negative. Other CEC-related populations (CD34-positive and CD45-negative) were also evaluated including (1) nucleated and CD146-positive CECs, (2) activated CECs (CD105-positive and CD31-positive), (3) progenitor CECs (CD133-positive and CD31-positive), and (4) CD146-positive and CD31-positive CECs. CEC values were reported as CEC per microliter of blood.

Biopsies of primary tumors

Core biopsies (research biopsies) of the primary breast tumor were obtained at baseline and on week 6 after initiation of therapy. Biopsy samples were obtained using a 14–18 gauge core needle. At least three core biopsies were obtained, two snap frozen individually and a third one was taken for the preparation of paraffin-embedded blocks. Frozen tissues were used for high-throughput genomic analyses after macrodissection.

De-identified fresh frozen tumor tissue biopsy specimens were obtained from the University of Alabama at Birmingham’s Comprehensive Cancer Center Tissue Procurement Shared Facility. The specimens were macrodissected by a board-certified pathologist at the Tissue Procurement Shared Facility to enrich for tumor cell content and remove adjacent uninvolved tissue. The dissected specimens were weighed and transferred to a 15-mL conical tube containing ceramic beads, and RLT Buffer (Qiagen) plus 1% BME was added so that the tube contained 35 μL of buffer for each milligram of tissue. The conical tubes containing tissue, ceramic beads, and buffer were agitated in a MP Biomedicals FastPrep machine at 6.5 m/s for 90 s to homogenize the tissue. The homogenized tissue was stored at − 80 °C. Total RNA was extracted from 350 μL of tissue homogenate (equivalent to 10 mg of tissue) using the Norgen Animal Tissue RNA Purification Kit (Norgen Biotek Corporation). Cell lysate was treated with Proteinase K before it was applied to the column and on-column DNAse treatment was performed according to the manufacturer’s instructions. Total RNA was eluted from the columns and quantified using the Qubit RNA Assay Kit and the Qubit 2.0 fluorometer (Invitrogen). We used 70 ng of the purified total RNA for most samples. Small RNA libraries were created as described previously [22] except that no blocking of any small RNA species was performed. In short, small RNA samples had adaptor oligonucleotide ligated to both ends. cDNA was created by reverse transcription and amplified by 15 cycles of PCR. PCR products were purified and run under extreme denaturing conditions on acrylamide gels. Bands corresponding to libraries with insert sizes of 15–30 base pairs were excised and purified. The libraries were sequenced on a HiSeq 2500 (Illumina).

For qPCR assays, 4 μl of purified RNA for each sample was subjected to reverse transcription using the Universal cDNA Synthesis Kit II (Exiqon). Resulting cDNA was diluted 1:20 and combined with primers and 2X Power SYBR mix (ThermoFisher). The primers were obtained from Exiqon (Supplementary Appendix). The qPCR reactions were run on a QuantStudio 6 Flex under default conditions. CT values were called using the QuantStudio software with default settings.

Statistical considerations

The trial was sized to estimate the within-arm pCR rates to a certain precision; it was not powered to compare arms. Based on a prior pCR rate of 12% in the investigational arm of our pilot trial [15], with n = 50, this trial had a standard error of ± 0.052 and, using the Blyth-Still-Casella method, provided a pCR rate and two-sided confidence interval of 5.4–23.3%. In the control arm, with n = 25 and projected pCR rate of < 1% [23,24,25], no pCR events were anticipated and the probability of observing one pCR event was 22%. The purpose of the control arm was to characterize the pCR response as similar to the prior (historical control) trials with single-agent neoadjuvant letrozole.

Descriptive statistics were used to summarize baseline characteristics of all patients by randomization arm in order to assess relative comparability of the two arms of the study. Mean, standard deviation, and range were summarized for continuous variables, and frequencies (percentages) were calculated for categorical variables. The differences between the pCR and response rates were determined along with 95% exact confidence intervals (CIs). All toxicities were recorded and summarized by calculating frequencies and proportions of toxicity grades for each arm.

A Student’s t test was used to examine the correlation between response and CTC and CEC values at baseline and at each time point for all patients and patients within each arm. Response was categorized into binary variables, 0 for stable and progressive disease and 1 for partial and complete response while CTC and CEC values were tested as continuous variables. Correlation between response and the changes in CTC and CEC numbers between baseline and other time points and the changes between time points were also examined. Correlations with a p value < 0.05 were considered significant.

Genomic data analyses

Raw sequencing reads were analyzed as described previously [26]. Briefly, reads were aligned to the human whole genome (hg19) requiring perfect matches. Features were created by merging overlapping alignments and total read counts reported for each.

To generate a small RNA-based classifier of treatment response, patients were categorized as responders if they had achieved a pathologic treatment response > 30% and non-responders if they had stable or progressive disease. Pathologic response was assessed by comparing the maximum cumulative diameter of the target lesion(s) at the time of diagnosis as assessed by imaging studies with the size of the tumor in the final surgical pathology. Due to the small number of patients who achieved pCR or microscopic residual disease, a genomic classifier on the basis of achievement of pCR or microscopic residual disease could not be generated. Differential expression of feature counts was assessed using DESeq2 [27]. From previously generated full-length RNA-seq data on these samples (data not presented), we had quality control (QC) metrics (fraction of reads mapping to mRNA and cDNA concentration). Significance was assessed using a likelihood ratio test between the full (response variable + QC metrics) and null (QC metrics) models as implemented in DESeq2.

We predicted the agreement between qPCR and sequencing data using previously described methods [26]. We evaluated the relative proportions of 3′ ends of small RNA features, predicting that those with many, equal proportioned 3′ ends would not yield concordant data between qPCR and sequencing measurements. We then selected small RNA features that had low p values from the likelihood ratio test and were predicted to yield concordant qPCR and sequencing measurements. We employed LASSO regression on this subset as implemented in the R package glmnet [28], in order to find small RNA features with optimal ability to classify responder vs. non-responder status. A binomial regression model to predict responder status was created from these optimal classifiers, operating on z-score values from both sequencing and qPCR measurements. Raw CTs from qPCR measurements were normalized using hsa-miR-22-3p. We chose this miRNA based on its robust expression and very low variance across samples in the small RNA sequencing data. Accuracies of model predictions were calculated by summing the true positives and true negatives divided by the total predictions. ROC curve areas and p values were calculated using the R package “verification.”


Patients and dispositions

Patient demographics and tumor characteristics are shown in Table 1. In the study, 75 patients were randomly assigned, 50 in the Let/Bev arm and 25 in the Let arm (2:1 ratio). All patients received at least 1 cycle of therapy. The median age for the patients enrolled in the Let/Bev arm was 61.4 years (range, 50.4 to 81.9) and 65 for the Let arm (range, 50.4 to 86.3). All patients had an ECOG performance status of 0. The protocol arms were well balanced for race, stage, nodal status, and tumor type. The proportion of patients with grade 2 tumors was higher in the Let arm (76% vs. 58% respectively), while no patient with grade 3 tumors was randomized to the Let arm (16% vs. 0% respectively). Randomization was not stratified for any demographic or disease parameter.

Table 1 Patient characteristics

Figure 1 provides the CONSORT diagram for the study. On the Let arm, 92% of the patients (n = 24) completed therapy, whereas 76% (n = 38) completed therapy on the Let/Bev arm. On the Let arm, discontinuations were due to the absence of response, whereas adverse events were the main reasons for discontinuation in the Let/Bev arm. As specified in the protocol, all patients underwent surgery 4 weeks after the last dose of Bev but continued taking Let until the day before surgery.

Fig. 1

CONSORT diagram


Pathologic response

Of the 50 patients enrolled in the Let/Bev arm, 45 patients had available surgical pathology for evaluation of pathologic response (consent withdrawal, n = 2; transition to chemotherapy at the discretion of the treating physician, n = 2; discontinuation due to toxicity, n = 1, Fig. 1). Of the 25 patients in the Let arm, 24 patients had available surgical pathology (1 patient refused surgery upon completion of protocol therapy, Fig. 1).

As seen in Table 2 and Fig. 2a, the rates of pCR, microscopic residual disease, PR, and objective response were numerically higher in the Let/Bev arm vs. Let arm (11.1% vs. none, 8.9%. vs none, 44.4% vs. 37.5%, and 64.4% vs. 37.5% respectively; no formal statistical testing, per protocol). One of the 5 patients who achieved pCR in the Let/Bev arm had residual ductal carcinoma in situ. These data are in agreement with the responses reported previously in a single-arm study of letrozole and bevacizumab in the same patient population [15]. The rates of downstaging (from stage 3/2 to stage 1/0) and lymph node reversion (from positive upon diagnosis to negative at the time of surgery) were also numerically higher in the Let/Bev arm (44.4% vs. 37.5% and 40.0% vs. 26.6% respectively). The rates of stable and progressive disease (disease same as or greater than its original radiological assessment at the time of diagnosis) were numerically lower in the Let/Bev arm vs. Let arm (28.8% vs 50% and 6.6% vs 12.5%, respectively). Of the 6 patients with progressive disease in the trial (3 in each arm), the surgical pathology after protocol therapy in 3 patients (Let arm, n = 2; Bev/Let arm, n = 1) was consistent with invasive lobular carcinoma, grade 1 or 2. Given the diffuse growth pattern of cancer cells and lack of desmoplastic reaction in invasive lobular carcinoma that collectively render its detection by physical examination, imaging, and even gross pathologic evaluation challenging [29], this finding may reflect clinical and radiographic underestimation of the original extent of the disease rather than actual tumor growth..

Table 2 Pathologic response per protocol arm
Fig. 2

Waterfall plot of response by treatment arm. a Pathologic responses were assessed by comparing the maximum cumulative diameter of the target lesion(s) at the time of diagnosis as assessed by imaging studies with the size of the tumor in the final surgical pathology. b Radiological responses were assessed according to RECIST by breast ultrasound or breast MRI in selected cases (whichever modality was consistently obtained throughout the protocol therapy and most accurately assessed the status of the tumor)

pCR and microscopic residual disease were seen only in patients with grade 1 or 2 disease at diagnosis. Additionally, all patients who achieved pCR and microscopic residual disease had stage 2 disease at diagnosis.

Of the 45 patients assigned to the Let/Bev arm with available surgical pathology, 37 were deemed at diagnosis to require mastectomy. Of these patients, 15 (40.5%; 95% CI, 24.8–57.9%) underwent lumpectomy following preoperative letrozole and bevacizumab. Of the 23 patients assigned to the Let arm with available surgical pathology (declined surgery n = 1, definitive surgery without preoperative medical therapy not assessed n = 1), 17 were deemed at diagnosis to require mastectomy. Of these patients, 9 (53%; 95% CI, 27.8–77%) underwent lumpectomy following preoperative letrozole alone.

Radiological response

Forty-nine patients in the Let/Bev arm (one patient withdrew consent before the first assessment of radiological response due to non-protocol related issues) and 25 patients in the Let arm were evaluable for radiological response according to RECIST [16]; results are summarized in Table 3 and Fig. 2b. Complete, partial, and objective radiographic response rates were numerically higher in the Let/Bev arm (10% vs. 4%, 69% vs. 64%, and 79% vs. 68% respectively). Accordingly, the rates of stable and progressive disease were numerically lower in the Let/Bev arm vs. Let arm (21% vs. 28% and none vs. 4%, respectively). The absolute mean percent reduction of the maximum tumor diameter by imaging was − 51.08% (95% CI, − 43.5 to − 58.6) in the Let/Bev arm vs. − 39.2% in the Let arm (95% CI, − 26.4 to − 52.2). Two patients in the Let arm and 3 in the Let/Bev arm discontinued treatment due to inadequate clinical response by physical evaluation and imaging; both patients in the Let arm went for surgery without additional neoadjuvant therapy while in the Let/Bev arm, 2 patients proceeded to chemotherapy after completion of cycles 3 and 9 and 1 patient withdrew consent.

Table 3 Best radiological response per protocol arm


Twenty-five patients in the Let arm and 50 patients in the Let/Bev arm received at least one cycle of therapy and were eligible for toxicity evaluation. The majority of adverse events were grades 1 and 2.

As summarized in Table 4, the most common adverse events observed in more than 5% of the patients in the Let arm deemed by the investigators to be possibly related to the aromatase inhibitor were hot flashes (28%), arthralgias/joint stiffness (24%), fatigue (16%), myalgias/cramps (12%), nausea/vomiting (8%), and night sweats (8%). No grade 4 toxicity was seen in the Let arm and only 2 patients (8%) had grade 3 arthralgias/joint stiffness. No patient discontinued the medication in the letrozole alone arm due to toxicity.

Table 4 Adverse events of any grade seen in > 5% of patients in either protocol arm

Adverse events known to be associated with letrozole were observed in similar frequencies in both arms. As expected, there was additive toxicity in the Let/Bev arm as adverse events known to be associated with bevacizumab were observed. The most common adverse events observed in more than 5% of the patients in the Let/Bev arm deemed by the investigators to be possibly related to the protocol agents were arthralgias/joint stiffness (26%), hot flashes (18%), fatigue (14%), myalgias/cramps (6%) possibly attributable to letrozole; hypertension (32%), hemorrhagic and thrombotic events (16%, includes epistaxis (grade 1, n = 1; grade 2, n = 1), bleeding gingivae (grade 1, n = 1), hematoma (grade 1, n = 1), rectal bleeding (grade 1, n = 1), cerebrovascular ischemia (grade 3, n = 1), hematemesis (grade 3, n = 1), and syncope (grade 3, n = 1)), headache (14%), proteinuria (10%), dyspnea (10%), and rash (8%), possibly attributable to bevacizumab. The most frequent grade 3 adverse events seen in the Let/Bev arm were hypertension (8%), arthralgias/joint stiffness (4%), and proteinuria (2%). Only 1 patient had grade 4 toxicity in the combination arm which was proteinuria (2%). Seven patients discontinued therapy in the Let/Bev arm due to adverse events, but pathology was available after the surgical procedure: 1 patient developed necrosis followed by infection at the site where the dye for sentinel lymph node mapping was injected (protocol therapy discontinued on cycle 5); 3 patients had proteinuria (grades 2, 3, and 4, n = 1 each) and received 6, 6, and 8 cycles, respectively; 1 patient had allergic reaction to the infusion and received 5 cycles; 1 patient had grade 3 fatigue and received 8 cycles; and 1 patient had activation of a Sjögren syndrome and received 4 cycles. In addition, one patient discontinued therapy in the Let/Bev arm due to grade 3 cerebrovascular ischemia but surgical pathology was not available. No healing problems at the time of surgery were seen.

Three serious adverse events were reported in the Let/Bev arm which were classified by the investigators as possibly related to the protocol therapy: 1 patient was admitted with confusion to the hospital due to a possible TIA (patient continued therapy), 1 patient had a neck abscess which required admission to the hospital for IV antibiotics, and 1 patient had grade 3 cerebrovascular ischemia (patient was taken off study). This last patient had presented with expressive aphasia. Thorough workup did not reveal acute infarct, mass, or intracranial hemorrhage. Echocardiogram was consistent with moderate to severe aortic stenosis. Her hospital stay was extended to provide adequate blood pressure control (blood pressure was 220/90 mmHg at presentation). On follow-up, there were no residual neurologic symptoms. No deaths associated with the treatment agents were seen in the trial.

Circulating tumor cells and circulating endothelial cells

CTCs were detected only in 14% of the patients. The mean number of circulating tumor cells was 0.05/mL (SD 0.24, range 0–1.7). Due to the limited number of patients with CTCs, the numbers at baseline and other time points as well as the change in CTCs between time points did not show a statistically significant correlation with response (data not shown).

Samples for CEC analysis at baseline were available in 77% of the patients (Supplementary Table 1). At baseline, the mean number of CD31+ CECs was 20.6 cells/μL (standard deviation [SD] 63.0, range 2.7–489.3); CD146+ CECs, 2.9 cells/μL (SD 3.2, range 0–13.9); activated CECs (CD105+/CD31+cells), 12.9 cells/μL (SD 44.4, range 1.2–296.2); progenitor CECs (CD133+/CD31+), 0.8 cell/μL (SD 1.6, range 0–9.7); and CD146+/CD31+ CECs, 0.7 cell/μL (SD 2.5, range 0–18.5). Absolute numbers and changes in CEC levels did not show significant correlation with response when all patients, regardless of the treatment arm in which they were enrolled, were considered (Supplementary figure 1). However, in the Let/Bev arm, we observed a statistically significant correlation between response and the levels of activated CECs (CD105+/CD31+) at week 18 (p = 0.0432). In addition, changes in the levels of activated CECs (CD105+/CD31+) (p = 0.0038), progenitor CECs (CD133+/CD31+) (p = 0.0120), and CD146+/CD31+ CECs (p = 0.0087) between weeks 6 and 18 significantly correlated to response.

Generation of a small RNA classifier that predicts response to letrozole plus bevacizumab

Macrodissected tumor tissue from research biopsies obtained at the time of diagnosis was successfully sequenced in 26 and 15 patients assigned to Let/Bev and Let arms respectively. Differential gene expression analyses yielded a large number of small RNAs to be significantly associated with response to letrozole/bevacizumab. hsa-mIR-187-3p was the highest ranking differentially expressed small RNA between responders and non-responders to Let/Bev in terms of significance (Supplementary Appendix). Increased hsa-mIR-187-3p expression in breast cancer has been associated with an aggressive, invasive phenotype and has been shown to be an independent predictor of outcome [30]. To overcome the biases associated with small RNA library preparations and the complexities of genomic data analyses, we filtered for small RNAs likely to reproduce the sequencing data when measured by qPCR (see the “Methods” section). We identified a set of five small RNAs as an optimal letrozole/bevacizumab classifier by LASSO regression (Fig. 3a). These include two miRNAs (hsa-miR-141-5p, hsa-miR-449), two small nucleolar RNAs (SNORD51, SNORA21), and fragments of 7SK RNA. hsa-mIR-187-3p was not predicted to generate qPCR data concordant with sequencing measurements.

Fig. 3

a Quantile-quantile plot of small RNA p values for association with letrozole/bevacizumab response. The observed p values from the association of small RNA expression values with letrozole/bevacizumab response (see the “Methods” section) are plotted against the null (uniform) distribution. Small RNAs selected by LASSO regression for an optimal classifier model (see the “Methods” section) are marked by black dots and labeled. b Classifier model performance on three cohorts. Classifier values calculated from qPCR data for each patient tumor sample (x-axis) across three cohorts are plotted on the y-axis (see the “Methods” section). Bar color indicates responder status. The dotted line indicates the common threshold used to calculate model accuracy. The “discovery cohort” is the cohort of patients who received letrozole/bevacizumab in the present study. The “validation cohort” is the cohort of patients who received the same treatment in a separate independent clinical trial [15]. The “letrozole cohort” consists of patients who received letrozole in the present study. The model was trained on sequencing data from the discovery cohort

To confirm the sequencing data for these five candidate classifier small RNAs, we measured them in the same set of samples by qPCR. We found a significant agreement (p < 0.01) between sequencing and qPCR data for all but one (SNORA21, Supplementary Figure 2). To evaluate their performance as classifiers of response to letrozole/bevacizumab in both types of measurement (sequencing and qPCR), we z-scored the sequencing data and created a binomial regression model comprising the five LASSO-selected small RNAs. When we applied the classifier to the sequencing data that it was trained on, the classifier had perfect performance (ROC AUC = 1, p = 7.1 × 10−6). We then applied the classifier to z-scored qPCR data from the same training samples and observed that the classifier still performed well (ROC AUC = 0.934, p = 6.2 10−5, accuracy = 0.885; Fig. 3b). We chose a threshold for calling responders in this data set that maximized model accuracy (Fig. 3b).

As a robust verification of the classifier, we applied it to z-scored qPCR data from breast tumor samples collected from patients that received the same treatment (letrozole/bevacizumab) in a separate independent clinical trial [15]. We successfully sequenced 16 macrodissected tumor samples from the validation cohort. We observed good performance of the classifier in these samples as well (Fig. 3b), with only moderate decrease in performance metrics (ROC AUC = 0.836, p = 9.0 × 10−3, accuracy = 0.875). Lastly, we applied the model to patients from the main trial that received letrozole alone to evaluate the specificity of the classifier for the letrozole/bevacizumab combination, rather than for letrozole alone. We observed non-significant performance of the classifier on patients who received letrozole alone (ROC AUC = 0.53, p = 0.44, accuracy = 0.467; Fig. 3b).


The addition of bevacizumab to letrozole capitalizes on multiple lines of evidence that implicate VEGF and angiogenesis as major mediators of resistance to endocrine therapies in luminal breast cancer [2,3,4,5,6, 13,14,15]. This study met its primary endpoint by determining the pCR rate of neoadjuvant/preoperative letrozole in combination with bevacizumab for 24 weeks to be 11.1% (CI, 3.7–24.1%) (20% including pCR rates and rates of microscopic residual disease). This pCR rate is highly similar to the 12% observed in our pilot study conducted in the same patient population [15] and at least comparable to the pCR rate achieved with neoadjuvant anthracycline-taxane-based chemotherapy [31]. To our knowledge, this pCR rate is the highest reported in endocrine-based neoadjuvant clinical trials, exceeding the < 10% pCR rate reported in other preoperative studies combining endocrine therapy with growth factor pathway inhibitors [32]. In neoadjuvant trials that compared letrozole with or without taselisib and with or without everolimus, the pCR rate in the combination arm was 1.8% [25] and 1.4% [23], respectively. Similarly, the pCR rates of preoperative combinations of endocrine therapy with CDK4/6 inhibitors have been low (0% [33] and 3.8% [34]) as well. No pCR or microscopic residual disease was seen in the control arm, and this result is consistent with multiple prior trials with preoperative aromatase inhibitors. It should be noted that our study was not powered for direct comparisons between arms and a larger study will be required to further confirm these results.

Interestingly, pCRs were confined to patients with stage 2 and well to moderately differentiated tumors. The presence of established angiogenesis in larger tumors promoted and sustained by multiple growth factors may underpin this observation. Indeed, in preclinical studies, anti-VEGF therapies have been shown to preferentially target smaller peripheral blood vessels while centrally located large vessels were spared [35].

Although the toxicities in the investigational arm were not synergistic, the addition of bevacizumab did result in a set of adverse events, known to be associated with bevacizumab, including hypertension, proteinuria, and hemorrhagic and thrombotic events. The nature and frequencies of the bevacizumab-associated adverse events in this study were very similar to the ones reported in other studies of bevacizumab combined with endocrine therapy [36,37,38]. We should acknowledge the uniformly very good performance status at baseline, generally favorable prognosis, and absence of residual toxicities from prior therapies in this population as well as the favorable toxicity profile of letrozole.

The achievement of pCR with preoperative therapy has been associated with favorable long-term outcomes, and the addition of bevacizumab to preoperative chemotherapy has led to higher pCR rates [31, 39,40,41]. Notably though, this benefit has not consistently translated into a long-term disease-free or definitive overall survival advantage [42,43,44,45]. A limitation of our study is the lack of long-term follow-up for recurrence, survival, and toxicities associated with bevacizumab. However, the choice of postoperative therapy was left to the physician’s discretion and many patients with large residual disease or limited treatment effect eventually received adjuvant cytotoxic chemotherapy. Accordingly, whether the higher pCR rate and better overall responses achieved with the addition of bevacizumab translate into better long-term outcomes remains an open question. The controversies regarding antiangiogenic therapy underscore the need to identify a biomarker for optimal patient selection; such a biomarker has remained largely elusive [46,47,48]. To address this gap in the knowledge, in the present study, we explored the utility of circulating tumor and endothelial cells and of a small RNA-based genomic classifier.

The level of CTCs in these patients with early-stage hormone receptor-positive breast cancer was low, in agreement with the literature [49]. The number of CTCs before and during treatment and the changes in levels during therapy were not associated with response. However, the low number of patients with detectable CTCs limits meaningful conclusions. In contrast, CECs were detected in most patients. Interestingly, we observed a correlation between response to letrozole and bevacizumab and the levels of activated CECs at week 18. Changes in the levels of certain CEC subpopulations (activated CECs, progenitor CECs, and CD146+/CD31+ CECs) between weeks 6 and 18 were also significantly associated with response to letrozole and bevacizumab. The significance of these associations remains unclear considering the lack of any correlations at other time points and the limited sample size.

By contrast to the CTCs and CECs, our small RNA-based classifier may constitute a potential biomarker predictive of response to bevacizumab. All small RNA components of the classifier have been implicated in the pathobiology of breast cancer. The hsa-miR-141-5p is the less abundant product from the same pre-miRNA with hsa-miR-141-3p and may in fact represent a surrogate measurement of hsa-miR-141-3p. Overexpression of hsa-miR-141-3p and its family member hsa-miR-200c-3p in MDA-MB-231 breast cancer cells promotes significant migration and invasion and enhances VEGF-A secretion [50]. VEGF-A neutralizing antibodies abrogate these phenotypic consequences conferred by hsa-miR-141-3p and hsa-miR-200c-3p overexpression [50]. hsa-miR-449a targets cysteine-rich intestinal protein 2 (CRIP2) mRNA, a transcription factor and a tumor suppressor [51]. CRIP2 interacts with the NF-κB/p65 to inhibit its DNA-binding ability to the promoter regions of the major pro-angiogenesis cytokines critical for tumor progression, including VEGF [51]. Indeed, increased expression of CRIP2 is associated with impaired tumor angiogenesis [52]. C/D box snoRNAs (like small nucleolar RNA, C/D box 51 (SNORD51)) and H/ACA box snoRNAs (like small nucleolar RNA, H/ACA box 21(SNORA21)) define the target sites for 2′-O-ribose methylation and pseudouridylation on ribosomal RNA (rRNA), respectively [53]. To our knowledge, a direct association between small nucleolar RNAs and VEGFA expression or function has not been reported. The possibility exists though that altered rRNA modifications introduced by the small nucleolar RNAs may modulate the RNA affinities and translational capabilities of the ribosomes and, consequently, prioritize the translation of VEGF mRNA as shown with the rRNA methyl-transferase fibrillarin [54, 55]. Lastly, 7SK snRNA and the La-related protein LARP7 are required for the integrity of the 7SK snRNP complex [56]. This complex sequesters the positive transcription elongation factor b (P-TEFb). Decreased levels of LARP7 and 7SK snRNA redistribute P-TEFb to the transcriptionally active super elongation complex, resulting in accelerated transcription of transcription factors that promote breast cancer epithelial-mesenchymal transition, invasion, and metastasis [56]. Although a direct association between 7SK snRNA or 7SK snRNP and VEGFA expression or function has not been reported, the possibility exists that low levels of 7SK snRNA may allow for upregulated VEGFA transcription. Although our small RNA-based classifier was verified in an independent cohort, incorporation of the classifier as a correlative study in clinical trials with bevacizumab will allow for further refinements and validation. We should note that we also performed RNA-seq on available tumor samples but RNA-seq data analyses did not yield significant classifiers.


In stage 2 and 3 hormone receptor-positive and Her2-negative breast cancer, the achievement of pCR with endocrine therapy combinations, including combinations with PI3K inhibitors and CDK4/6 inhibitors has been an elusive target. In this patient population, the addition of bevacizumab to preoperative letrozole resulted in a pCR rate of 11.1% (95% CI, 3.7–24.1%) at the cost of additive toxicities. We have developed a small RNA-based classifier, which, on the basis of the levels of five small RNAs in a pretreatment biopsy, can select patients more likely to respond to bevacizumab.

Availability of data and materials

All raw and processed sequencing data generated in this study have been submitted to the NCBI Gene Expression Omnibus (GEO; under accession number GSE142308. Analyzed data and the R code used to generate the main and supplementary figures are included as supplementary information files.



Circulating endothelial cells


Confidence interval


Complete response


Circulating tumor cells


Estrogen receptor


Magnetic resonance imaging


Multigated acquisition


Pathologic complete remission


Progressive disease


Partial response


Response Evaluation Criteria in Solid Tumors


Stable disease


Vascular endothelial growth factor


  1. 1.

    Hanahan D, Folkman J. Patterns and emerging mechanisms of the angiogenic switch during tumorigenesis. Cell. 1996;86(3):353–64.

  2. 2.

    Zhang HT, Craft P, Scott PA, Ziche M, Weich HA, Harris AL, et al. Enhancement of tumor growth and vascular density by transfection of vascular endothelial cell growth factor into MCF-7 human breast carcinoma cells. J Natl Cancer Inst. 1995;87(3):213–9.

  3. 3.

    Guo P, Fang Q, Tao HQ, Schafer CA, Fenton BM, Ding I, et al. Overexpression of vascular endothelial growth factor by MCF-7 breast cancer cells promotes estrogen-independent tumor growth in vivo. Cancer Res. 2003;63(15):4684–91.

  4. 4.

    Qu Z, Van Ginkel S, Roy AM, Westbrook L, Nasrin M, Maxuitenko Y, et al. Vascular endothelial growth factor reduces tamoxifen efficacy and promotes metastatic colonization and desmoplasia in breast tumors. Cancer Res. 2008;68(15):6232–40.

  5. 5.

    Pinto MP, Badtke MM, Dudevoir ML, Harrell JC, Jacobsen BM, Horwitz KB. Vascular endothelial growth factor secreted by activated stroma enhances angiogenesis and hormone-independent growth of estrogen receptor-positive breast cancer. Cancer Res. 2010;70(7):2655–64.

  6. 6.

    Losordo DW, Isner JM. Estrogen and angiogenesis: a review. Arterioscler Thromb Vasc Biol. 2001;21(1):6–12.

  7. 7.

    Hyder SM, Stancel GM. Regulation of angiogenic growth factors in the female reproductive tract by estrogens and progestins. Mol Endocrinol (Baltimore). 1999;13(6):806–11.

  8. 8.

    Hyder SM, Nawaz Z, Chiappetta C, Stancel GM. Identification of functional estrogen response elements in the gene coding for the potent angiogenic factor vascular endothelial growth factor. Cancer Res. 2000;60(12):3183–90.

  9. 9.

    Buteau-Lozano H, Ancelin M, Lardeux B, Milanini J, Perrot-Applanat M. Transcriptional regulation of vascular endothelial growth factor by estradiol and tamoxifen in breast cancer cells: a complex interplay between estrogen receptors alpha and beta. Cancer Res. 2002;62(17):4977–84.

  10. 10.

    Mueller MD, Vigne JL, Minchenko A, Lebovic DI, Leitman DC, Taylor RN. Regulation of vascular endothelial growth factor (VEGF) gene transcription by estrogen receptors alpha and beta. Proc Natl Acad Sci U S A. 2000;97(20):10972–7.

  11. 11.

    Applanat MP, Buteau-Lozano H, Herve MA, Corpet A. Vascular endothelial growth factor is a target gene for estrogen receptor and contributes to breast cancer progression. Adv Exp Med Biol. 2008;617:437–44.

  12. 12.

    Foekens JA, Peters HA, Grebenchtchikov N, Look MP, Meijer-van Gelder ME, Geurts-Moespot A, et al. High tumor levels of vascular endothelial growth factor predict poor response to systemic therapy in advanced breast cancer. Cancer Res. 2001;61(14):5407–14.

  13. 13.

    Linderholm B, Grankvist K, Wilking N, Johansson M, Tavelin B, Henriksson R. Correlation of vascular endothelial growth factor content with recurrences, survival, and first relapse site in primary node-positive breast carcinoma after adjuvant treatment. J Clin Oncol. 2000;18(7):1423–31.

  14. 14.

    Linderholm B, Tavelin B, Grankvist K, Henriksson R. Vascular endothelial growth factor is of high prognostic value in node-negative breast carcinoma. J Clin Oncol. 1998;16(9):3121–8.

  15. 15.

    Forero-Torres A, Saleh MN, Galleshaw JA, Jones CF, Shah JJ, Percent IJ, et al. Pilot trial of preoperative (neoadjuvant) letrozole in combination with bevacizumab in postmenopausal women with newly diagnosed estrogen receptor- or progesterone receptor-positive breast cancer. Clin Breast Cancer. 2010;10(4):275–80.

  16. 16.

    Therasse P, Arbuck SG, Eisenhauer EA, Wanders J, Kaplan RS, Rubinstein L, et al. New guidelines to evaluate the response to treatment in solid tumors. European Organization for Research and Treatment of Cancer, National Cancer Institute of the United States, National Cancer Institute of Canada. J Natl Cancer Inst. 2000;92(3):205–16.

  17. 17.

    FDA. Draft guidance for industry: pathologic complete response in neoadjuvant treatment of high-risk early-stage breast cancer: use as an endpoint to support accelerated approval 2012. Available from: Accessed 23 July 2012.

  18. 18.

    Cristofanilli M, Budd GT, Ellis MJ, Stopeck A, Matera J, Miller MC, et al. Circulating tumor cells, disease progression, and survival in metastatic breast cancer. N Engl J Med. 2004;351(8):781–91.

  19. 19.

    Allard WJ, Matera J, Miller MC, Repollet M, Connelly MC, Rao C, et al. Tumor cells circulate in the peripheral blood of all major carcinomas but not in healthy subjects or patients with nonmalignant diseases. Clin Cancer Res. 2004;10(20):6897–904.

  20. 20.

    Garcia JA, Rosenberg JE, Weinberg V, Scott J, Frohlich M, Park JW, et al. Evaluation and significance of circulating epithelial cells in patients with hormone-refractory prostate cancer. BJU Int. 2007;99(3):519–24.

  21. 21.

    Rugo HS, Chien AJ, Franco SX, Stopeck AT, Glencer A, Lahiri S, et al. A phase II study of lapatinib and bevacizumab as treatment for HER2-overexpressing metastatic breast cancer. Breast Cancer Res Treat. 2011;134(1):13–20.

  22. 22.

    Roberts BS, Hardigan AA, Kirby MK, Fitz-Gerald MB, Wilcox CM, Kimberly RP, et al. Blocking of targeted microRNAs from next-generation sequencing libraries. Nucleic Acids Res. 2015;43(21):e145.

  23. 23.

    Baselga J, Semiglazov V, van Dam P, Manikhas A, Bellet M, Mayordomo J, et al. Phase II randomized study of neoadjuvant everolimus plus letrozole compared with placebo plus letrozole in patients with estrogen receptor-positive breast cancer. J Clin Oncol. 2009;27(16):2630–7.

  24. 24.

    Ellis MJ, Tao Y, Luo J, A'Hern R, Evans DB, Bhatnagar AS, et al. Outcome prediction for estrogen receptor-positive breast cancer based on postneoadjuvant endocrine therapy tumor characteristics. J Natl Cancer Inst. 2008;100(19):1380–8.

  25. 25.

    Saura C, Hlauschek D, Oliveira M, Zardavas D, Jallitsch-Halper A, de la Pena L, et al. Neoadjuvant letrozole plus taselisib versus letrozole plus placebo in postmenopausal women with oestrogen receptor-positive, HER2-negative, early-stage breast cancer (LORELEI): a multicentre, randomised, double-blind, placebo-controlled, phase 2 trial. Lancet Oncol. 2019.

  26. 26.

    Roberts BS, Hardigan AA, Moore DE, Ramaker RC, Jones AL, Fitz-Gerald MB, et al. Discovery and validation of circulating biomarkers of colorectal adenoma by high-depth small RNA sequencing. Clin Cancer Res. 2018;24(9):2092–9.

  27. 27.

    Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.

  28. 28.

    Friedman J, Hastie T, Tibshirani R. Regularization paths for generalized linear models via coordinate descent. J Stat Softw. 2010;33(1):1–22.

  29. 29.

    Johnson K, Sarma D, Hwang ES. Lobular breast cancer series: imaging. Breast Cancer Res. 2015;17:94.

  30. 30.

    Mulrane L, Madden SF, Brennan DJ, Gremel G, McGee SF, McNally S, et al. miR-187 is an independent prognostic factor in breast cancer and confers increased invasive potential in vitro. Clin Cancer Res. 2012;18(24):6702–13.

  31. 31.

    von Minckwitz G, Eidtmann H, Rezai M, Fasching PA, Tesch H, Eggemann H, et al. Neoadjuvant chemotherapy and bevacizumab for HER2-negative breast cancer. N Engl J Med. 2012;366(4):299–309.

  32. 32.

    Spring LM, Gupta A, Reynolds KL, Gadd MA, Ellisen LW, Isakoff SJ, et al. Neoadjuvant endocrine therapy for estrogen receptor-positive breast cancer: a systematic review and meta-analysis. JAMA Oncol. 2016;2(11):1477–86.

  33. 33.

    Ma CX, Gao F, Luo J, Northfelt DW, Goetz M, Forero A, et al. NeoPalAna: neoadjuvant palbociclib, a cyclin-dependent kinase 4/6 inhibitor, and anastrozole for clinical stage 2 or 3 estrogen receptor-positive breast cancer. Clin Cancer Res. 2017;23(15):4055–65.

  34. 34.

    Cottu P, D'Hondt V, Dureau S, Lerebours F, Desmoulins I, Heudel PE, et al. Letrozole and palbociclib versus chemotherapy as neoadjuvant therapy of high-risk luminal breast cancer. Ann Oncol. 2018;29(12):2334–40.

  35. 35.

    Davis DW, Inoue K, Dinney CP, Hicklin DJ, Abbruzzese JL, McConkey DJ. Regional effects of an antivascular endothelial growth factor receptor monoclonal antibody on receptor phosphorylation and apoptosis in human 253J B-V bladder cancer xenografts. Cancer Res. 2004;64(13):4601–10.

  36. 36.

    Dickler MN, Barry WT, Cirrincione CT, Ellis MJ, Moynahan ME, Innocenti F, et al. Phase III trial evaluating letrozole as first-line endocrine therapy with or without bevacizumab for the treatment of postmenopausal women with hormone receptor-positive advanced-stage breast cancer: CALGB 40503 (Alliance). J Clin Oncol. 2016;34(22):2602–9.

  37. 37.

    Martin M, Loibl S, Hyslop T, De la Haba-Rodriguez J, Aktas B, Cirrincione CT, et al. Evaluating the addition of bevacizumab to endocrine therapy as first-line treatment for hormone receptor-positive metastatic breast cancer: a pooled analysis from the LEA (GEICAM/2006-11_GBG51) and CALGB 40503 (Alliance) trials. Eur J Cancer. 2019;117:91–8.

  38. 38.

    Martin M, Loibl S, von Minckwitz G, Morales S, Martinez N, Guerrero A, et al. Phase III trial evaluating the addition of bevacizumab to endocrine therapy as first-line treatment for advanced breast cancer: the letrozole/fulvestrant and avastin (LEA) study. J Clin Oncol. 2015;33(9):1045–52.

  39. 39.

    Earl HM, Hiller L, Dunn JA, Blenkinsop C, Grybowicz L, Vallier AL, et al. Efficacy of neoadjuvant bevacizumab added to docetaxel followed by fluorouracil, epirubicin, and cyclophosphamide, for women with HER2-negative early breast cancer (ARTemis): an open-label, randomised, phase 3 trial. Lancet Oncol. 2015;16(6):656–66.

  40. 40.

    Bear HD, Tang G, Rastogi P, Geyer CE Jr, Robidoux A, Atkins JN, et al. Bevacizumab added to neoadjuvant chemotherapy for breast cancer. N Engl J Med. 2012;366(4):310–20.

  41. 41.

    Nahleh ZA, Barlow WE, Hayes DF, Schott AF, Gralow JR, Sikov WM, et al. SWOG S0800 (NCI CDR0000636131): addition of bevacizumab to neoadjuvant nab-paclitaxel with dose-dense doxorubicin and cyclophosphamide improves pathologic complete response (pCR) rates in inflammatory or locally advanced breast cancer. Breast Cancer Res Treat. 2016;158(3):485–95.

  42. 42.

    Nahleh Z, Botrus G, Dwivedi A, Jennings M, Nagy S, Tfayli A. Bevacizumab in the neoadjuvant treatment of human epidermal growth factor receptor 2-negative breast cancer: a meta-analysis of randomized controlled trials. Mol Clin Oncol. 2019;10(3):357–65.

  43. 43.

    Earl HM, Hiller L, Dunn JA, Blenkinsop C, Grybowicz L, Vallier AL, et al. Disease-free and overall survival at 3.5 years for neoadjuvant bevacizumab added to docetaxel followed by fluorouracil, epirubicin and cyclophosphamide, for women with HER2 negative early breast cancer: ARTemis trial. Ann Oncol. 2017;28(8):1817–24.

  44. 44.

    Bear HD, Tang G, Rastogi P, Geyer CE Jr, Liu Q, Robidoux A, et al. Neoadjuvant plus adjuvant bevacizumab in early breast cancer (NSABP B-40 [NRG oncology]): secondary outcomes of a phase 3, randomised controlled trial. Lancet Oncol. 2015;16(9):1037–48.

  45. 45.

    von Minckwitz G, Loibl S, Untch M, Eidtmann H, Rezai M, Fasching PA, et al. Survival after neoadjuvant chemotherapy with or without bevacizumab or everolimus for HER2-negative primary breast cancer (GBG 44-GeparQuinto)dagger. Ann Oncol. 2014;25(12):2363–72.

  46. 46.

    Miles D, Cameron D, Bondarenko I, Manzyuk L, Alcedo JC, Lopez RI, et al. Bevacizumab plus paclitaxel versus placebo plus paclitaxel as first-line therapy for HER2-negative metastatic breast cancer (MERiDiAN): a double-blind placebo-controlled randomised phase III trial with prospective biomarker evaluation. Eur J Cancer. 2017;70:146–55.

  47. 47.

    Schneider BP, Gray RJ, Radovich M, Shen F, Vance G, Li L, et al. Prognostic and predictive value of tumor vascular endothelial growth factor gene amplification in metastatic breast cancer treated with paclitaxel with and without bevacizumab; results from ECOG 2100 trial. Clin Cancer Res. 2013;19(5):1281–9.

  48. 48.

    Jain RK, Duda DG, Willett CG, Sahani DV, Zhu AX, Loeffler JS, et al. Biomarkers of response and resistance to antiangiogenic therapy. Nat Rev Clin Oncol. 2009;6(6):327–38.

  49. 49.

    Lucci A, Hall CS, Lodhi AK, Bhattacharyya A, Anderson AE, Xiao L, et al. Circulating tumour cells in non-metastatic breast cancer: a prospective study. Lancet Oncol. 2012;13(7):688–95.

  50. 50.

    Choi SK, Kim HS, Jin T, Hwang EH, Jung M, Moon WK. Overexpression of the miR-141/200c cluster promotes the migratory and invasive ability of triple-negative breast cancer cells through the activation of the FAK and PI3K/AKT signaling pathways by secreting VEGF-A. BMC Cancer. 2016;16:570.

  51. 51.

    Cheung AK, Ko JM, Lung HL, Chan KW, Stanbridge EJ, Zabarovsky E, et al. Cysteine-rich intestinal protein 2 (CRIP2) acts as a repressor of NF-kappaB-mediated proangiogenic cytokine transcription to suppress tumorigenesis and angiogenesis. Proc Natl Acad Sci U S A. 2011;108(20):8390–5.

  52. 52.

    Shi W, Bruce J, Lee M, Yue S, Rowe M, Pintilie M, et al. MiR-449a promotes breast cancer progression by targeting CRIP2. Oncotarget. 2016;7(14):18906–18.

  53. 53.

    Williams GT, Farzaneh F. Are snoRNAs and snoRNA host genes new players in cancer? Nat Rev Cancer. 2012;12(2):84–8.

  54. 54.

    Bastide A, Karaa Z, Bornes S, Hieblot C, Lacazette E, Prats H, et al. An upstream open reading frame within an IRES controls expression of a specific VEGF-A isoform. Nucleic Acids Res. 2008;36(7):2434–45.

  55. 55.

    Marcel V, Ghayad SE, Belin S, Therizols G, Morel AP, Solano-Gonzalez E, et al. p53 acts as a safeguard of translational control by regulating fibrillarin and rRNA methylation in cancer. Cancer Cell. 2013;24(3):318–30.

  56. 56.

    Ji X, Lu H, Zhou Q, Luo K. LARP7 suppresses P-TEFb activity to inhibit breast cancer progression and metastasis. Elife. 2014;3:e02907.

Download references


The authors wish to acknowledge the Translational Breast Cancer Research Consortium and its supporters, Avon Foundation for Women, the Breast Cancer Research Foundation (BCRF), and Susan G. Komen, for supporting this study. The authors thank the Breast Cancer Research Foundation of Alabama for covering the publication costs.


This study was supported in part by funding from the National Cancer Institute Specialized Program of Research Excellence (SPORE) in Breast Cancer (P50CA089019). Correlative studies were in part supported by Susan G. Komen (IIR13265422, PI: Richard M Myers). The funding agencies were not involved in the design of the study and collection and analysis or in the interpretation of data and writing of the manuscript.

Author information

Protocol writing was done by AFT and AFL. Conception and design were performed by AFT, AFL, RMM, and YL. Provision of study materials or patients was done by NUL, MCL, HSR, SP, RN, AMS, LAC, MNS, and JFD. Collection and assembly of data were performed by all authors. Data analysis and interpretation were performed by CV, BSR, KEV, HSR, YL, WEG, RMM, and AFT. Manuscript writing was done by CV and AFT. Final approval of manuscript was performed by all authors.

Correspondence to Andres Forero-Torres.

Ethics declarations

Ethics approval and consent to participate

The protocol was reviewed and approved by each participating institution, and the study followed the Declaration of Helsinki and good clinical practice guidelines. All patients gave informed consent.

Consent for publication

Not applicable.

Competing interests

C.V. consulting or advisory board: Daiichi-Sankyo, Genentech (unpaid); research funding (paid to the University of Alabama at Birmingham): Roche, Pfizer, Incyte, Novartis, Pharmacyclics, Tracon, Innocrin, H3 Biomedicine, Zymeworks. A.F.T received support from Genentech to conduct clinical and laboratory research (paid to the University of Alabama at Birmingham). A.F.T. is currently an employee of Seattle Genetics. N.U.L.—consulting or advisory board: Puma, Daiichi-Sankyo, Seattle Genetics, Genentech; research funding (paid to institution): Genentech, Pfizer, Novartis, Seattle Genetics, Array Biopharma. The remaining authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Presentations/publications: Data from this work have been partially reported at the ASCO Annual Meeting in 2013 (Abstract # 527).

Supplementary information

Additional file 1. Supplementary Figure 1 Changes in circulating endothelial cells (CECs) during protocol therapy. Enumeration at baseline and at weeks 6, 18, and 24 of circulating endothelial cells (CECs) and related populations in blood samples collected from all patients assigned in both protocol arms. Absolute numbers and changes in CEC levels did not significantly correlate with response.

Additional file 2 Supplementary Figure 2 Comparison of small RNA sequencing data to qPCR. Upper panels. For the 5 small RNAs selected as the optimal classifier by LASSO, we plotted sequencing values (“vsd” for variance stabilized data by DESeq2, x-axis) versus qPCR data (-dCt, y-axis). The R-squared and p-value of the goodness-of-fit by linear regression are provided for each small RNA. Lower panels. Boxplots for each measurement type (qPCR: “-dCt”; sequencing: “vsd”) for each small RNA selected as the optimal classifier by LASSO between non-responders and responders to letrozole/bevacizumab. The p-values are from linear regression for qPCR data (“-dCt” y-axis) and from DESeq2 (“vsd” y-axis) for sequencing data.

Additional file 3. Supplementary Table 1 CEC enumeration in all trial participants.

Additional file 4. Appendix 1. LASSO selected classifiers and qPCR assay information.

Additional file 5. Appendix 2. Summary statistics for all tested small RNA features in sequencing data (sorted by letozole/bevacizumab responder P).

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Vaklavas, C., Roberts, B.S., Varley, K.E. et al. TBCRC 002: a phase II, randomized, open-label trial of preoperative letrozole with or without bevacizumab in postmenopausal women with newly diagnosed stage 2/3 hormone receptor-positive and HER2-negative breast cancer. Breast Cancer Res 22, 22 (2020).

Download citation


  • Preoperative therapy
  • Bevacizumab
  • Letrozole
  • Breast cancer
  • Hormone receptor-positive breast cancer
  • Luminal breast cancer
  • Postmenopausal