Risk assessment model for invasive breast cancer in Hong Kong women.

No risk assessment tool is available for identifying high risk population of breast cancer (BCa) in Hong Kong. A case-control study including 918 BCa cases and 923 controls was used to develop the risk assessment model among Hong Kong Chinese women.Each participant received an in-depth interview to obtain their lifestyle and environmental risk factors. Least absolute shrinkage and selection operator (LASSO) selection model was used to select the optimal risk factors (LASSO-model). A risk score system was constructed to evaluate the cumulative effects of selected factors. Bootstrap simulation was used to test the internal validation of the model. Model performance was evaluated by receiver-operator characteristic curves and the area under the curve (AUC).Age, number of parity, number of BCa cases in 1st-degree relatives, exposure to light at night, and sleep quality were the common risk factors for all women. Alcohol drinking was included for premenopausal women; body mass index, age at menarche, age at 1st give birth, breast feeding, using of oral contraceptive, hormone replacement treatment, and history of benign breast diseases were included for postmenopausal women. The AUCs were 0.640 (95% CI, 0.598-0.681) and 0.655 (95% CI, 0.621-0.653) for pre- and postmenopausal women, respectively. Further subgroup evaluation revealed that the model performance was better for women aged 50 to 70 years or ER-positive.This BCa risk assessment tool in Hong Kong Chinese women based on LASSO selection is promising, which shows a slightly higher discriminative accuracy than those developed in other populations.

Introduction

Breast cancer (BCa) is the most common cancer among Hong Kong women. About 6% of Hong Kong women (1 in 17) will develop BCa during their lives. Although this lifetime risk is half of that (1 in 8) in the United States, the BCa incidence rate of Hong Kong women has almost doubled during recent 3 decades. It has been reported that the BCa incidence rates of recent Chinese generations in Singapore and Taiwan (with similar economic development as Hong Kong) were even higher than those in the United States. Hence, there is an urgent need for developing efficient prevention strategies against BCa among Chinese women of Hong Kong.

The 1st step for BCa prevention is to accurately identify individuals at an increased BCa risk. Gail model is one of the most popular BCa risk assessment and prediction tools that has been widely used clinically to predict risk for individual women with Caucasian and African ethnic origin,[ ] although it was originally designed to identify women at an increased risk for entry onto chemoprevention trials and has limited discriminatory power (about 60% for the values of receiver-operator characteristic [ROC] curves).[ ] There are several other models developed to predict BCa risk in different populations by using similar indicators, such as the model from International Breast Cancer Intervention Study (IBIS) which focuses on familial BCa and the Mayo BBD-to-BCa model which focuses on women with history of benign breast diseases. Up till now, only 2 studies assessed the risk of BCa (Shanghai and Nanjing) in mainland China[ ] by including a limited number of SNPs and some risk factors, while the discriminatory power was also about 60% in these studies. The risk prediction model among Chinese women in Singapore also shows a similar accuracy. Furthermore, different fertility pattern, lifestyle, and environmental factors between mainland and Hong Kong also makes the direct application of prediction tools based on mainland Chinese to Hong Kong women a challenge.

Some established risk factors are known to vary by menopausal status ; however, menopausal status is absent in most risk assessment models. Further, some potential risk factors identified in recent years (e.g., exposure to light at night [LAN], sleep disturbance, shift work with circadian disruption) [ ] were not addressed in the previous risk assessment models, which may lead to underestimated discriminatory power of these models. This study aims to construct a BCa risk assessment model specific to Hong Kong Chinese women, by including menopausal status and additional environmental risk factors that have never been considered in previous risk models.

Methods

Study subjects and specification of risk factors

The study protocol was approved by both the Clinical Research Ethics Committees of Joint Chinese University of Hong Kong – New Territories East Cluster (Joint CUHK-NTEC CREC) and the Kowloon West Cluster (KWC CREC). It was a hospital-based case–control study and the detailed process of subjects recruitment was described previously. In brief, eligible cases were consecutively identified from the Department of Surgery or Clinical Oncology of 3 hospitals in Hong Kong during the period November 2011 to March 2015. All included cases were Chinese women, aged between 20 and 84 with newly diagnosed primary BCa (ICD-10 code 50), and confirmed by histology. Specifically, we obtained the diagnostic results from either the tissue pathological or biopsy report for the cases receiving breast tissue dissection or biopsy, respectively. The response rate of cases was 90.8%, and main reasons for the nonresponses were no interest or poor medical condition. Each case was frequency matched in 5-year age groups by a control patient selected from the same hospital where the cases came from, with a response rate of 93.2%. Participants were excluded if they had a physician-diagnosed cancer history at any site.

After obtaining the written informed consent, face-to-face interviews were conducted by trained interviewers. A standardized questionnaire was used to obtain the information on demographic data, menstrual and reproductive history (such as age at menarche, menopausal status, parity, breast feeding, age at 1st give birth, and hormone replacement therapy), history of benign breast diseases, family cancer history, occupational history, smoking habit, alcohol drinking, self-report exposure to LAN (1, dark; 2, few bright; 3, little bright; and 4, bright), and self-report sleep quality (1, good; 2, common; 3, poor; and 4, poor with sleep pill). The clinical information (including estrogen receptor [ER]) for the cases and controls was extracted from the hospital medical records.

Statistical analyses

We transformed relevant continuous variables to the categorical scale for the use of clinical practice. In the development of the risk model, we initially performed univariate logistic regression analyses separately by menopausal status to identify potential predictors. To include more variables in the initial selection, predictors with P ≤ 0.10 were considered for further evaluation. Then, LASSO model was performed to select the potential predictors to be further used by shrinking the coefficients toward zero through setting a constraint on the sum of the absolute standardized coefficients. Shrinkage estimates with LASSO method in fitting model provided an important way for adjusting model's overfitting and preventing extreme predictions. The factors selected by LASSO selection were the predictors for the final model (LASSO-model). Third, multivariate logistic regression analyses were performed to generate the coefficients of all predictors from LASSO selection. We then assessed the effect of LASSO-model using a risk score analysis on the basis of a linear combination of the selected predictors weighted by their coefficients, using a formula as following:  , where S is the risk score, x the selected predictors, and n is the number of selected predictors. To facilitate result demonstration, the final risk score was multiplied by 10.

Ten-fold cross-validation was used to evaluate the model's internal validity by partitioning the original sample into a training set to train the model and a test set to evaluate it. The model performance was evaluated by ROC curves and the area under the curve (AUC) was used to classify the BCa cases and controls. Hosmer–Lemeshow fit test was used to assess the agreement between observed and model-predicted proportions of BCa events. The difference of AUCs was tested by a nonparametric approach developed by DeLong et al.

An optional model (OPT-model) which included the established risk factors of age, age at menarche, age at 1st give birth, number of BCa cases in 1st-degree relatives, and history of benign breast diseases was also constructed to compare the discrimination accuracy with LASSO-model. Their coefficients derived from multivariate logistic regression analysis were used as their category weightings. The discriminative power of LASSO-model was also evaluated among women in different age groups and different ER status of BCa. All statistical analyses were 2-sided and performed with Stata software (Version 11.2; StataCorp LP, TX).

Results

A total of 923 BCa cases and 918 age-matched controls were included in the final analysis. The characteristics of these subjects were summarized in Table 1. The mean age was 56.0 ± 11.8 years old in cases and 53.8 ± 11.8 years old in controls. BMI, age at menarche, number of parity, age at 1st give birth, ever breast feeding, hormone replace treatment, BCa history among 1st-degree relatives, and exposure to LAN were significantly different between cases and controls. After stratified by menopausal status, the alcohol drinking rate was significantly higher in premenopausal cases of BCa. The distribution of number of parity, breast feeding, and hormone replace treatment were significantly different between cases and controls among postmenopausal women (Supplemental Table 1). There were 784 BCa patients with data on ER status, and among them, 601 (76.7%) cases were ER positive (H-score > 50, or Allred score ≥3).

Table 1

Selected characteristics between breast cancer cases and controls of Hong Kong women.

The potential predictors selected by univariate logistic regression are presented in Table 2 and supplemental Table 2. These potential predictors were refined by LASSO selection (Table 2), which were generated by multivariate logistic regression model that were comparable to the OPT-model (supplemental Table 2). Eventually, 6 risk factors were included in the final LASSO model for premenopausal women and 12 factors were included for postmenopausal women. Age, number of parity, case number of BCa in first-degree relatives, LAN, and sleep quality were the common predictors for both pre- and postmenopausal women. The coefficient for most variables was similar in both groups of women, except for number of parity (0.137 and −0.184 for pre- and postmenopausal women, respectively). In addition to these factors, alcohol drinking was included for premenopausal women; BMI, age at menarche, age at 1st give birth, ever breast feeding, ever using of oral contraceptive, hormone replacement treatment, and history of benign breast diseases were included for postmenopausal women.

Table 2

Predictors for breast cancer by LASSO selection and their coefficient calculated by multivariate logistic regression model.

Figure 1 demonstrates the results of risk score analysis from the LASSO-model incorporating with a linear combination of the risk predictors weighted by their coefficients. The means of the risk scores was 7.43 ± 4.63 for controls and 9.74 ± 5.05 for BCa cases among premenopausal women, and 10.44 ± 5.32 for controls and 13.60 ± 5.97 for cases among postmenopausal women. The discrimination accuracy of this score system was determined by the ROC curves (Fig. 2), given the AUC of 0.640 (95% CI, 0.598–0.681) for premenopausal women and 0.655 (95% CI, 0.623–0.686) for postmenopausal women, respectively. Hosmer–Lemeshow fit test showed good agreement between observed and model-predicted proportions of BCa in both pre- and postmenopausal women (P = 0.302 and 0.848, respectively).

Figure 1

Distribution of risk scores calculated with a linear combination of the least absolute shrinkage and selection operator (LASSO) selected predictors weighted by multivariate logistic regression coefficient. (A) Premenopausal women; (B) postmenopausal women.

Figure 2

Discrimination performance of LASSO-model for breast cancer and comparison with OPT-model by receiver-operator characteristic analyses among pre- and postmenopausal Hong Kong women (AUC). (A) Premenopausal women; (B) postmenopausal women. AUC = area under the curve, LASSO = least absolute shrinkage and selection operator, OPT-model = optional model, ROC = receiver-operator characteristic.

A 10-fold cross-validation method was applied to split all the samples randomly into 10 partitions, and taking 9-fold as training set, and another 1-fold as the validation set, 1000 times were repeated to evaluate the OPT-model's ability. The risk score calculated by the coefficients of the predictors was evaluated with each random datasets, and the results were similar with the observed data (average AUC: 0.621 and 0.632 for pre- and postmenopausal women, respectively). These simulation results showed a robust internal consistence between the estimated effects of predictors based on the original results and the bootstrapped results.

Compared to the OPT-model, the discrimination accuracy of LASSO-model improved significantly (Fig. 2). The AUCs were increased from 0.586 to 0.640 among premenopausal women (P = 0.011) and from 0.621 to 0.655 among postmenopausal women (P = 0.006), respectively. The comparison of discrimination accuracy among women with different ages is presented in Table 3. This score system showed better accuracy among postmenopausal women aged between 50 and 70 years old, and among premenopausal women aged between 40 and 50 years. For ER positive BCa, the accuracy was higher in postmenopausal women with AUC of 0.663 (95% CI, 0.628–0.696) (Supplemental Table 3).

Table 3

Prediction accuracy of risk scores among different age groups in Hong Kong women.

Discussion

There was no risk assessment tool specifically for identifying Hong Kong women at increased risk of developing BCa prior to this study. We developed a risk assessment model for invasive BCa of Hong Kong women. By adding modifiable risk factors (e.g., exposure to LAN and sleep quality) that had never been addressed in previous risk models, the discrimination accuracy significantly improved from 0.586 to 0.640 among premenopausal women and from 0.621 to 0.655 among postmenopausal women, and a better model was demonstrated for the postmenopausal women aged between 50 and 70 years and for the ER positive BCa. Overall, the performance of this model is slightly superior to other risk models in which the discriminatory power is generally around 60%.

In our study, only 5 of 12 predictors in the risk model of postmenopausal women were enrolled in the risk model of premenopausal women, and the effect of parity was opposite between pre- and postmenopausal women (coefficient were 0.137 and −0.184, respectively). These results suggest that the influence of menopause is rather comprehensive than a single event for the initiation of BCa. Hormones produced by the ovaries play central roles in breast tissue development, maintenance, and tumorigenesis.[ ] However, menopause changes the existing circumstance which might affect the entire risk profile of BCa. For example, obesity is associated with BCa incidence, but it showed different risk patterns according to the women's menopausal status.[ ] The study of KoBRCAT observed the same phenomenon as our study. The major risk factors and their effects in KoBRCAT were different after dividing the participants into 2 subgroups by age of 50 years old. Our study is able to attribute this variety in terms of risk factors to the etiology of BCa incidences before and after menopause, and thus provides scientific basis for decision makers to resettle the proper role of menopause in BCa risk assessment and prediction.

Gail model was used to predict the individual risk for developing BCa in western population. However, although recalibrations have been made, it does not perform well in Asian women. The recalibrated Gail model performed good calibration in the study of Singapore Breast Screening Programme (model name, Gail-SBSP) and Seoul Breast Cancer Study (model name, KoBRCAT), but the discriminative power was limited with the AUC of about 0.60.[ ] After involving modifiable risk factors of BMI, oral contraceptive usage, and exercise, the discriminative power of KoBRCAT improved significantly to 0.63 for women <50 years old and 0.65 for women ≥50 years old, respectively. In the current study, we enrolled more predictors in the risk model including BMI, breast feeding history, LAN, sleep quality and alcohol drinking, and the discriminative accuracy was significantly improved.

It is the 1st time that LAN exposure was enrolled into the BCa risk assessment model. There is a theoretical link between LAN and BCa occurrence. Long-term exposure to LAN can suppress the secretion of hormone melatonin and cause the circadian disruption. Melatonin appears to protect against cancer development, and decreased secretion of melatonin may induce continuous production of estrogen and alter the function of ER. There is no epidemiological report which focused on LAN and BCa risk. However, as an extreme scenario of exposure to LAN, night shift work showed positive association with BCa risk. Except for the working circumstance with bright light, mechanism study showed a significantly suppressed of nocturnal melatonin secretion by very low illuminance of blue-appearing light. In this study, the brightness might be very low in bedroom during the sleep time and we did not assess the light levels objectively. But this result warrantees further studies on LAN exposure and the risk of BCa.

Three studies reported the relationship between sleep quality and risk of BCa, [ ] but no association was found; however, we observed a negative relation between poor sleep quality and BCa risk, and the underlying reason remains unclear. We tried to exclude those with chronic comorbidities (may potentially influence sleep quality) from both cases and controls, and results remain unchanged. Further studies with prospective study design are needed to investigate the role of sleep quality and BCa occurrence.

BMI, which was used as the indicator for obesity, was included into our postmenopausal women model as well as in the KoBRCAT model for women aged older than 50 years. BMI was not selected into the model for premenopausal women in our study or for young women in KoBRCAT, although it was reported to have a protective effect among premenopausal women.[ ] Considering potential differences in the manifestation of weight gain between Asian and western women, waist circumference may be a better indicator that should be addressed in future studies.

Alcohol drinking has long been shown to be associated with an increased risk of BCa. Every 10 g of alcohol consumption per day might increase the risk up to 12%.[ ] In this study, alcohol drinking was involved into the risk assessment model for premenopausal women with high weight in the score system. However, it was excluded from the model for postmenopausal women. The possible reason is that the drinking rate is higher in young women than old women in Hong Kong, according to a recent study of BCa in Hong Kong.

Our BCa risk assessment model showed better performance for ER positive BCa and women aged between 50–70 years. It is easy to understand the better accuracy for ER positive BCa, because most selected factors are directly or indirectly associated with previous estrogen levels. The lower discriminative power for women older than age 70 years might be the results of a weaker association of risk factors with BCa risk among older women and fewer participants.

There are several strengths of this study. First, we constructed the risk assessment model by menopausal status which might more clearly reflect the risk changes before and after stopping menstruation. Second, several modifiable risk factors were involved in the risk model, which might be used for the primary prevention of BCa. Third, our study is the 1st to provide a potential value of using LAN exposure to predict the risk of BCa. Nevertheless, several limitations should be mentioned. First, this study used a case–control study design and recall bias might be a concern. We compared information on a special group of 117 patients (not including in this study) who were handled as BCa initially but finally confirmed to be noncases, and found that the prevalence of various risk factors were slightly lower than the true cases, which suggested the information bias may not be a major methodological issue. And we also conducted test–retest reliability among 25% cases and controls after 6 weeks of the initial interview. The kappa value was 0.62 which indicates a relatively high reliability and low misclassification of the potential variables. Therefore, the issues of residual confounding effect and reliability of the potential risk predictors in our study should be low. Second, the study design did not allow us to capture the changes of some modifiable risk factors such as BMI and sleep quality through menopause. These changes need to be investigated in further study with cohort study design. Third, we recruited controls from hospitals rather than from general population, and some exposures may differ between hospital patients and the general population. To reduce the potential bias caused by using hospital-based controls, we excluded controls with breast-related diseases and recruited controls with a broad spectrum of diagnosis. Fourth, germline mutation of BRCA1/2 status was unavailable to us; however, this should not influence the risk assessment significantly as the prevalence of germline BRCA1/2 mutations is very low (1%) in Hong Kong. Also genetic variants and mammographic density were unavailable for this study. They are promising factors for BCa prediction, although limited contribution was observed in previous studies. Fifth, we did not collect potential risk factors’ changes as the subject progress through menopause, because the potential misclassification from the recall of previous environmental exposures in different time points may be substantial, especially for the nonhabitual exposures. Finally, all BCa patients came from 3 public hospitals that may not represent all cases in the entire population; however, our case samples are comparable to those obtained from the Hong Kong Cancer Registry in which a wide coverage of approximate 95% of the Hong Kong general population is indicated regarding age and histological subtypes.

In conclusion, considering the rapid increase of BCa incidence in recent decades, a better risk assessment model specific to Hong Kong women is urgently needed for assessing and predicting BCa risk. We developed a novel risk assessment model by including menopausal status and emerging environmental risk factors that have never been explored in previous models, and demonstrated better discriminative accuracy than previous risk models in other populations; hence, results from our study have added new scientific evidence to the current literature in terms of risk prediction on BCa. We expect that newly developed model shall be used for the screening of high risk population of BCa and contribute to primary and secondary prevention of BCa in Hong Kong.

Acknowledgements

The authors thank Research Grants Council of Hong Kong (Grant number 474811) for the support. The authors also thank Miss Yin-shan Magdalene Leung, Hung-kuen Ivy Hsu, Kit-ping Apple Kwok, and Angela Chan for their assistance in patients’ recruitment and data collection.