|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: KDY GHD ZZS ZMS. Performed the experiments: KDY. Analyzed the data: KDY JJL GHD JW. Contributed reagents/materials/analysis tools: JW ZZS ZMS. Wrote the paper: KDY JJL.
To investigate the relationship between age and axillary lymph node (LN) involvement in Chinese breast cancer patients, and to replicate a recently identified piecewise relationship between age and LN involvement.
A dataset, consisting of 3,715 patients (with complete information on study variables) with operable breast cancer consecutively surgically treated between 1996 and 2006, was derived from the database of Shanghai Cancer Hospital. Univariate and multivariate logistic regression were employed to analyze the relationship between age and LN. We subsequently performed a similar analysis on another dataset including 1,832 consecutive patients treated between 2007 and 2008 to replicate our findings in the first dataset.
A U-shaped relationship (previously observed in two European populations) between age and LN status failed to be replicated in our dataset of Chinese patients. Instead, we observed a linear rather than piecewise relationship. After multivariate adjustment, the linear relationship was still present. Moreover, the interaction between age and LN involvement was not modified by tumor size. The odds of LN involvement decreased by 1.5% for each year increase in age (OR 0.985, 95% CI 0.979–0.991, P<0.001). Breast cancer subtypes were also associated with LN status. Proportions of basal-like and ERBB2+ subtypes decreased with increasing age. The observations in the first dataset were successfully replicated in a second independent dataset.
We confirmed a straightforward but not piecewise relationship between age and LN status in Chinese patients. The different pattern between Chinese and European elderly patients should be considered when making clinical decisions.
Despite relatively lower incidence of breast cancer in Asian countries than in the west , China has documented a 20–30% increase in urban areas over the past decade , . The disease status of axillary lymph nodes (LNs) is the most significant prognostic factor for patients with breast cancer  . At present, although the use of sentinel LNs identification and sampling procedure could be reliably performed in selected early stage patients by a trained multidisciplinary team , axillary LN dissection (ALND) is still the mainstay of the surgical management of breast cancer in China .
The determination of factors associated with axillary LN involvement may help us predict or identify a subgroup of patients with a low rate of node involvement. Regarding the relationship between age and LN status, some investigator suggested the older patients were associated with an increased probability of nodes involvement , while others proposed that tumors of elderly patients were biologically more favorable  and had a decrease in LN involvement . In addition, Wildiers et al. recently reported a large sample size study (referred to as the Leuven study) examining the relationship between patient age and risk of axillary LN involvement . Interestingly, they identified a piecewise effect of age on LN positivity. Because of these discordant results, we here decided to perform a large scale retrospective analysis for two aims. One aim was to validate the piecewise U-shaped relationship between age and LN status, which had been replicated in two European patient populations; the other was to explore the relationship between age and LN status in Chinese breast cancer patients, which has been seldom studied.
The study subjects were from the Breast Malignancy Database established by the Department of Surgery, Shanghai Cancer Hospital of Fudan University in Shanghai, China. The database has more than 10,000 case records including detailed clinicopathologic data and follow-up results. The information of this database has been reported elsewhere . All patients gave written informed consent for their information to be stored in the hospital database and used for research, and this study was approved by the Ethical Committee of Shanghai Cancer Hospital of Fudan University. The first dataset we used was derived from the entire database of the period between 1996 and 2006, in which period there were total 6,030 consecutive patients with breast malignancies. Patients selected for the present retrospective study fulfilled the following inclusion criteria: (1) female patients diagnosed with ipsilateral invasive breast cancer. Patients with breast carcinoma in situ (with or without microinvasion) and breast sarcoma were excluded; (2) all these patients received primary surgery in our hospital; (3) the pathological examination of patients' tumor specimens was performed in the Department of Pathology in our hospital; (4) patients were operable without any evidence of metastases at diagnosis. Preoperative evaluation and examination had been described previously ; (5) patients receiving neoadjuvant systemic therapy (chemotherapy and/or hormone therapy) and preoperative irradiation were excluded. A preoperative diagnostic biopsy was allowed. As a result, 4,263 patients met all these criteria (see Figure 1 for patient selection).
In order to validate the results of the Leuven study, we chose the following variables for analysis: age at diagnosis (year), pathological tumor size (the largest diameter), tumor grade, pathology (invasive ductal cancer vs. invasive lobular cancer vs. other invasive cancers), involvement of axillary LN (number of positive nodes), estrogen receptor (ER) and progesterone receptor (PR) status, and HER2/neu status. Age at diagnosis and tumor sizes were recorded as continuous variables, the other variables treated as categorical data. Axillary LN status was analyzed either as a binary variable (positive or negative) or as a ternary variable according to the number of involved LNs (0, 1–3, or ≥4). Determination of tumor grading, ER, PR and HER2 status was performed by pathologists in the Department of Pathology in our hospital. The work was done according to established procedures which had been described elsewhere .
According to molecular profile of breast cancer subtypes , , the immunohistochemical (IHC) surrogates for four major subtypes had been developed: luminal A (ER+ and/or PR+, HER2−), luminal B (ER+ and/or PR+, HER2+), ERBB2+ (ER−, PR−, HER2+), and basal-like (ER−, PR−, HER2-, also known as triple-negative) . In the present study, we analyzed the status of ER, PR and HER2 receptors in a combined manner. The luminal A and luminal B breast cancer were combined as the ‘luminal-like subtype’ (ER+ and/or PR+, any HER2 status).
Of the 4,263 eligible breast cancer patients, some of them had missing values of the studied variables, 4.3% of tumor size (n=182), 12.5% of combined receptors status (n=532), and 1.4% of axillary LN status (n=58). Of note, 35.7% cases (n=1,523) were not evaluated for tumor histological grade (most breast carcinomas of special histologic types such as mucinous, tubular, medullary, and papillary carcinoma were not evaluated). Because of such a high proportion, we did not exclude those cases with missing values of tumor grade. Finally, 3,715 patients (87.1%) with complete information of age, tumor size, pathology, LN status, and ER, PR and HER2 status were selected for analysis.
The second dataset for validation was also derived from the entire database, and the patients corresponded to the breast cancer cases of the period between 2007 and 2008 (Figure 1). During that period, more than 3,000 consecutive breast malignancy patients received primary surgical treatment in our department. A total 3,112 cases were entered into our database. Among these, 2,172 patients fulfilled the inclusion criteria applied for patient selection in the first dataset.
The IHC data (ER/PR/HER2) in the second dataset is currently being checked and replenished, and only half of patients (n=997, 45.9%) in the second dataset are available for ER, PR and HER2. Therefore, we excluded the independent variable “breast cancer subtype” from further analysis. Percentages of other missing data on tumor size, tumor grade and LN involvement were 15.0%, 20.4% and 3.0%, respectively. Finally, we selected 1,832 (84.3%) patients with complete information (except tumor grade and subtype) for analysis. Similar to the situation in the first dataset, we performed multivariate logistic regression to predict LN involvement using three independent variables (age, tumor size and pathology).
We regressed on age using nonparametric logistic regression based on smoothing (LOWESS) method  to explore the relationship pattern in our study subjects. In descriptive statistics, associations between categorical variables were tested using Pearson's χ2 test. The odds ratios (OR) with 95% confidence interval (CI) for relationship between each variable and LN involvement (yes or no) were calculated using logistic regression. We employed multivariate logistic regression to predict LN involvement (method: backward stepwise, likelihood ratio). For the first dataset, the regression model was established based on four basic independent variables (age, tumor size, pathology and subtype). Similarly, for the second dataset, the model was established based on three basic independent variables (age, tumor size and pathology). Breast cancer subtype was chosen for modeling in the first dataset, but not in the second dataset. Hosmer-Lemeshow test was performed to assess goodness-of-fit of model. A P-value less than or equal to 0.05 was considered statistically significant. Statistical analysis was performed using Stata/SE version 10.0 (Stata, College Station, TX) and SPSS Software version 12.0 (SPSS, Chicago, IL, USA).
The basic characteristics of the patients from the first dataset (1996–2006, n=3,715) are shown in Table 1. Since more than half of breast cancer patients in China are premenopausal and the patients older than 80 years are relatively few , we arbitrarily divided the study subjects into five groups according to age: younger than 40 years, 40–49 years, 50–59 years, 60–69 years, and 70 years and older.
As the Leuven study previously demonstrated a U-shaped relationship between age and LN involvement , we first tried to replicate that result in our population. However, in our study, the LOWESS smoothing plot displayed a linear rather than a piecewise effect of age on LN status (Figure 2A). Age was negatively related to LN involvement in the whole group (OR per year 0.985, 95% CI 0.979–0.990 Table 2). Other independent variables including tumor size (OR per centimeter 1.244, 95% CI 1.182–1.310), pathology (OR 0.625, 95% CI 0.551–0.709), ER (OR of positive 0.645, 95% CI 0.566–0.735) and PR status (OR of positive 0.820, 95% CI 0.720–0.933), and HER2 status (OR of positive 1.363, 95% CI 1.167–1.591) had significant relationships with LN involvement. When combining ER, PR and HER2 as variables of breast cancer subtype, ERBB2+ and basal-like subtypes, both contributed to the risk of LN involvement (OR 1.773, 95% CI 1.543–2.037).
In the second dataset, the univariate effect of age on LN involvement using LOWESS smoothing method also displayed a similar linear relationship that was observed in the first study (Figure 2B). Age was negatively related to LN involvement (OR per year 0.986, 95% CI 0.978–0.994 Table 2). Other independent variables associated with LN involvement included tumor size (OR per centimeter 1.244, 95% CI 1.182–1.310) and pathology (OR 0.565, 95% CI 0.449–0.711). When we combined the first dataset and the second dataset together, the LOWESS smoothing plot showed the same trend, i.e., that older women had a lower probability of LN involvement (Figure 2C).
After adjustment of other predictors, the linear relationship between age and LN involvement was still present. As the results from the first dataset showed, the odds of LN involvement decreased by 1.5% for each year increase in age (OR 0.985, 95% CI 0.979–0.991), increased by 25–30% for each centimeter increase in tumor size (OR 1.262, 95% CI 1.197–1.331), and increased in ERBB2+ and basal-like subtypes (OR 1.303, 95% CI 1.203–1.411). Special histologic type (such as mucinous, tubular, medullary, and papillary carcinoma) also displayed a decreased risk of LN involvement compared with invasive ductal cancer. Similarly, the reverse relationship between age and LN involvement was successfully replicated.
Since the increase in age had a negative effect on LN involvement, and no piecewise effect was observed, we therefore did not develop piecewise models for prediction of node involvement. One logistic regression model worked well for patients of any age. We regressed a model using four independent variables (age, size, pathology and subtype) from the first dataset and three independent variables (age, size, and pathology) from the second dataset, respectively. The goodness-of-fit tests suggested that both models were good fits. Hazard ratios of each independent variable are presented in Table 2. The probability of lymph node involvement was estimated as eL/(1+eL), where the value of L was derived by multivariate logistic regression analysis. The dependent variable was “lymph node involvement”, and the independent variables for the first dataset model were age, pathology, size, and subtype; the independent variables for the second dataset model were age, pathology, and size. The values of age and size were continuous. For pathology, invasive ductal=1, invasive lobular=2, and other invasive carcinomas=3; for subtype, luminal-like=1, ERBB+=2, and basal-like=3.
Formula for the first model:
Formula for the second model:
The Leuven study revealed an interaction between age and tumor size on the frequency of LN involvement. It was demonstrated that the piecewise effect for age on the LN involvement was clearer in small tumors but not in tumors larger than 3.5 cm . However, in our study, we observed that the effect of age on LN positivity was not affected by tumor size through subgroup analysis (Table 3). Despite subdivision of the study subjects according to tumor size (≤2.0 cm, 2.1–3.5 cm, ≥3.6 cm), the reverse relationship between age and LN status did change either in the first dataset or in the combined dataset. Furthermore, the predicted probability of LN involvement was negatively associated with age in each subgroup (Figure 3A) in the first dataset, which was replicated in the second dataset (Figure 3B).
A piecewise U-shaped relationship between age and LN involvement has been recently reported for two European breast cancer populations . The present study attempted to validate it in our study subjects of Chinese breast cancer patients. However, our results displayed a linear rather than U-shaped relationship. The linear relationship observed in our first dataset was subsequently replicated in a second independent dataset. Successful validation strengthens the reliability of our results, indicating an increase in age related to a continuous decrease in LN involvement.
There could be a number of reasons responsible for the inconsistence between our study and the Leuven study. First, selection bias may be a causal factor. Results from the National Surgical Adjuvant Breast Project-B04 trial had demonstrated that ALND did not improve survival in patients with clinically negative axillary nodes . Subsequent clinical trials demonstrated that there was no significant survival difference between ALND and no axillary dissection in the selected patients with older age small tumor, and clinically negative axilla LNs –. Therefore, in clinical practice, a certain number of early stage elderly patients received lumpectomy or simple mastectomy only without ALND. These patients were excluded from analysis due to no pathological node status available. As a result, a higher percentage of LN involvement was observed in the older population. This conservative therapeutic strategy seems to be more prevalent in Europe –. In China, however, only a slight proportion of elderly patients were subjected to simple mastectomy or lumpectomy without ALND . The dataset which we used included all the consecutive patients who underwent surgery during the studied period. Every operable elderly patient was included despite her disease stage. In the Leuven study, however, whether all elderly patients during the studied period were recruited is unknown. Second, the epidemiologic difference of breast cancer between the east and the west could influence the outcomes as well. In most European and American countries, the peak of breast cancer incidence was observed in patients >70–75 years old. In contrast, the incidence peak emerged in the middle-age group (45–60 years) in Chinese population . Distinct age distribution of breast cancer patients between the two populations, as well as ethnic heterogeneity between Europeans and Asians, might make the nonrepeatability of the piecewise effect of age on LN involvement.
Our results were in line with the results in other studies. Previous studies suggested that age was negatively associated with node involvement, and tumor of older patients frequently had ER+ phenotype , –. Regarding intrinsic subtypes, we found decreasing proportions of basal-like subtype (an aggressive phenotype) and ERBB2+ subtype (another aggressive phenotype) with increasing age. Basal-like breast cancer is associated with higher grade, younger age, higher probability of node involvement and poor prognosis , , . Luminal-like tumors, however, have a good biologic behavior and are more frequently observed in older patients , , . Additionally, a special type, mucinous breast carcinoma, is more likely to occur in elderly patients. Mucinous carcinomas have substantially less nodal involvement, have higher expression rates of ER and/or PR, and have a lower S-phase fraction, compared with infiltrating ductal carcinoma . The features mentioned above may determine the indolent nature of breast cancers in elderly women, implying that advancing age is associated with more favorable tumor biology and less involved LNs. Nowadays, sentinel lymph node biopsy (SLNB) has become an alternative procedure of ALND. SLNB presents a lower risk of significant operative morbidity but similar benefit, and has emerged as the milestone advance in the surgical management of early-stage primary breast cancer , . In our hospital, the proportion of SLNB in the early breast cancer is relatively low (approximately 10%). Considering the low number of SLNB performed in our hospital and the low false-negative rate (<5%), we did not think SLNB would cause obvious bias of this study.
Some limitations of our study should be acknowledged. First, the first dataset had missing values of tumor grade, and the second dataset included no information on ER/PR/HER2 status. High percentage of missing values could result in increased likelihood of biased results. We predicted the probability of LN involvement using partial independent variables; a stable, robust linear relationship was repeatedly observed in different models. Second, the studied period of the first dataset was fairly long. There should be unavoidable biases in surgical procedure, histopathology evaluation, and ER/PR/HER2 detection. Third, the second dataset for validation, although independent, was also derived from the entire database, which could have increased the probability of successful validation. Fourth, in our Chinese database, the proportion of ≥80 years-old patients was low, and there were only 79 women with age ≥80 years, even after combining these two datasets. The small proportion of elderly and the under-representation of the ≥80 years-old group might explain the lack of reproducibility.
So far, most studies showed a linear rather than piecewise relationship between age and LN involvement. Some observations in the Leuven study are not straightforward to explain. Typically, although the authors offered the possible biologic reasons for their findings, such as the possibility that suppressed cellular immunity in the elderly might offset the favorable tumor biology, this explanation was not consistent with the observation that the increased rate of axillary node involvement in older patients versus younger patients was only limited to patients with smaller tumors . We stress that the novel piecewise relationship should be further verified in other populations.
In conclusion, we confirmed a straightforward but not piecewise relationship between age and LN status in Chinese breast cancer patients. The U-shaped relationship observed in European breast cancer patients does not appear to be applicable to breast cancer patients in China. The ethnic differences in LN involvement in elderly patients should be considered when making clinical decisions, as well as when establishing global clinical practice guidelines.
Competing Interests: The authors have declared that no competing interests exist.
Funding: This research is supported by grants from the National Basic Research Program of China (2006CB910501), 2009 Youth Foundation of Shanghai Public Health Bureau, 2009 Youth Foundation of Shanghai Medical College, and the National Natural Science Foundation of China (30971143, 30972936). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.