The Role of Socio-Demographic Factors in the Coverage of Breast Cancer Screening: Insights From a Quantile Regression Analysis

Background: In Flanders, breast cancer (BC) screening is performed in a population-based breast cancer screening program (BCSP), as well as in an opportunistic setting. Women with different socio-demographic characteristics are not equally covered by BC screening. Objective: To evaluate the role of socio-demographic characteristics on the lowest 10th and highest 90th quantile levels of BC screening coverage. Methods: The 2017 neighborhood-level coverage rates of 8,690 neighborhoods with women aged 50–69 and eligible for BCSP and opportunistic screening were linked to socio-demographic data. The association between socio-demographic characteristics and the coverage rates of BCSP and opportunistic screening was evaluated per quantile of coverage using multivariable quantile regression models, with specific attention to the lowest 10th and highest 90th quantiles. Results: The median coverage in the BCSP was 50%, 33.5% in the 10th quantile, and 64.5% in the 90th quantile. The median coverage of the opportunistic screening was 12, 4.2, and 24.8% in the 10th and 90th quantile, respectively. A lower coverage of BCSP was found in neighborhoods with more foreign residents and larger average household size, which were considered indicators for a lower socioeconomic status (SES). However, a higher average personal annual income, which was considered an indicator for a higher SES, was also found in neighborhoods with lower coverage of BCSP. For these neighborhoods, that have a relatively low and high SES, the negative association between the percentage of foreign residents, average household size, and average personal annual income and the coverage in the BCSP had the smallest regression coefficient and 95% confidence interval (CI) values were −0.75 (95% CI: −0.85, −0.65), −13.59 (95% CI: −15.81, −11.37), and −1.05 (95% CI: −1.18, −0.92), respectively, for the 10th quantile. The neighborhoods with higher coverage of opportunistic screening had a relatively higher average personal annual income, with the largest regression coefficient of 1.72 (95% CI: 1.59, 1.85) for the 90th quantile. Conclusions: Women from relatively low and high SES neighborhoods tend to participate less in the BCSP, whereas women with a relatively high SES tend to participate more in opportunistic screening. For women from low SES neighborhoods, tailored interventions are needed to improve the coverage of BCSP.


INTRODUCTION
Worldwide, breast cancer is the most common cause of cancer death in women (1). In 2018 the global age-standardized incidence and mortality rates of breast cancer were 54.4 and 11.6 per 100,000 women, respectively (1). Randomized controlled trials have confirmed that mammography screening can reduce the risk of breast cancer mortality by 20% for women aged 50-70 who attend this screening (2). The European guideline for quality assurance in breast cancer screening and diagnosis suggests to strive for 70% screening coverage in order to have a significant effect on breast cancer burden in the population. However, this percentage is not obtained in many countries where a breast cancer screening program (BCSP) has been established. In Europe, the mean screening coverage is about 50% (range 28-92%), meaning that a large proportion of women who are eligible for screening are not covered by the population breast cancer screening programs. This does not mean that these women do not receive any form of BC screening, since opportunistic BC screening exists in many countries (3).
To understand the reasons for this low coverage, several studies assessed the determinants of screening coverage. The main reported determinants of low screening coverage were related to low socioeconomic status (SES): lower income level, having an immigrant background, and being a single parent (4)(5)(6)(7). The conclusions of these studies are mainly derived from linear regression modeling, in which the change in the mean coverage is estimated as a function of the explanatory variables. However, the determinants of the lowest and highest screening coverage are more of interest to policymakers than the determinants of the average coverage. In order to improve the screening coverage of BCSP, it is important to know the characteristics of the women with the lowest coverage in the BCSP, as well as those of the women with the highest coverage in the opportunistic screening. These insights can assist policy makers to implement specific actions to increase the screening coverage of BCSP, especially for those groups that are not covered by opportunistic screening either.
The aim of this study was therefore to evaluate the effect of socio-demographic factors and indicators for SES at the lower end of screening coverage in the BCSP and the higher end of opportunistic screening coverage. Data of screening coverage in Flanders, Belgium were analyzed with a quantile regression model (8).

Screening in Flanders
Flanders, the largest and most populated region in Belgium (around 6 million inhabitants), is among the regions with the highest BC incidence in Europe, despite the early implementation of a BCSP in 2001, which offers biennial screening to women between 50 and 69 years (1,9). Women in Flanders can also choose to be screened for BC in opportunistic screening outside the BCSP. Where the BCSP is fully reimbursed by health insurance, opportunistic screening does not have a systematic quality control such as daily quality checks of mammography equipment, double reading and case-based feedback to readers (10). Although these two screening strategies coexist in Flanders, the combined coverage rate is only 64.1%, with a coverage rate of 50.0% for the BCSP and 14.1% for the opportunistic screening (11).

Data Source
For this analysis, publicly available data of 2017 from the Center for Cancer Detection and regional authorities in Flanders were linked at a statistical sector level. The statistical sector level is the smallest level at which administrative information is systematically collected in Flanders and is comparable to the neighborhood in literature (12). We will therefore use the term neighborhood hereafter. From the Center for Cancer Detection data were received of the screening invitations by the BCSP and attendances at the BCSP. The attendances to the opportunistic screening were identified from reimbursement data of the InterMutualist Agency (IMA). The data of opportunistic screening was sent to the Center for Cancer Detection regularly. These data were collected in a collaboration between the Belgian Cancer Registry and the Center for Cancer Detection, as described by their respective statutes (13,14). From the regional authorities data on socio-demographic characteristics were retrieved (15). Linkage was performed in a protected environment and provided at an aggregated neighborhood level to mask all the individual level information. To prevent identification of any individual, neighborhoods were excluded with less than five women screened. Due to that, 8,690 of the 9,490 neighborhoods in Flanders that provided the data of screening coverage rate in the BCSP and the opportunistic screening were included in this study. For these reasons, ethical approval was not needed for this study because the above-noted measures were taken for privacy protection.

Study Variables
The primary outcome in this analysis was the coverage rate of screening in the BCSP, which was defined as the percentage of women aged 50-69 screened in the BCSP. The secondary outcome in this analysis was the coverage rate of opportunistic screening which was defined as the percentage of women aged 50-69 screened in the opportunistic screening. The socio-demographic variables per neighborhood which can characterize the study population from different angles were used as covariates are listed in Table 1 and were defined as follows: (1) Population density: the number of residents per km²; (2) Same address as last year: the percentage of residents with the same address as last year; (3) Single parents: the percentage of unmarried residents that are single and live together with at least one child; (4) Married resident with child(ren) living at home: the percentage of married residents that live together with at least one child; (5) Unmarried cohabiting resident with child(ren) living at home: the percentage of unmarried residents that live together with at least one child; (6) The percentage of foreign residents: the percentage of residents without Belgian nationality; (7) Average personal annual income: the quotient of the total net taxable income and number of residents on January 1 of the tax year; (8) Average household size: the average number of persons in a household.

Statistical Analysis
Women were categorized into two groups, the group that was screened in the BCSP and the group that was screened in the opportunistic screening. In the first step of the analysis, the included neighborhoods in the analysis were ranked based on the coverage rates in the BCSP and in the opportunistic screening, respectively, and categorized in 9 quantiles (Q10-Q90). Then the socio-demographic variables were described per quantile. As we have aggregated data, these descriptions reflect the median estimates per quantile.
For each socio-demographic variable, to evaluate the deviations per quantile from the median regression, which is an indicator for the central tendency of the data, multivariable quantile regression models were performed. The quantile regression does not make any assumptions on the distribution of the observations neither on the residuals of the regression model (Supplementary Material Figure 1). The significance of the quantile regression coefficients was tested first to determine if the coefficients were significantly different from 0, then they were tested to determine if the coefficients were significantly different from the median regression coefficient. The quantile regression was considered to provide more information if both tests were significant (8,16). All statistical tests were two-sided and considered statistically significant at 0.05. All analyses were performed in R 4.0.2.

Description of the Coverage and Neighborhoods
Of the 8,690 included neighborhoods in the analysis, the median coverage in the BCSP was 50%. The coverage in the 10th quantile (Q10) was 33.5% and in the 90th quantile (Q90), it was 64.5%. The median percentage of single-parents, the median percentage of foreign residents, and the median average personal annual income decreased with increasing quantiles of BCSP coverage: for Q10-Q90 of the BCSP coverage, the median per quantile decreased from 3.8 to 3.0%, 11.1 to 3.9%, and e21,100 to 19,500, respectively ( Table 1). The median coverage of the opportunistic screening was 12%, and for the Q10-Q90 this increased from 4.2 to 24.8%. The median average personal annual income increased with increasing quantiles of the coverage of opportunistic screening from e18,700 to 21,700 for the Q10 to Q90 ( Table 2).

Determinants of the Coverage in the BCSP
A significant difference between the coefficients of the quantile regression and the median regression was observed for the population density, the percentage of foreign residents, the average household size, and the average personal annual income (Figure 1) (Supplementary Material Figure 2). These four determinants were all negatively associated with the coverage of BCSP with median regression coefficients of −0.97 (95%   Table 3). At Q10, the statistically significant association of the quantile regression of the percentage of foreign residents, average household size, and average personal annual income was stronger than the median regression, and the difference was statistically significant (

Determinants of Coverage of the Opportunistic Screening
The coefficients of the quantile regression for most of the determinants were significantly different from the median regression (Figure 2) (Supplementary Material Figure 3).    Table 4).

DISCUSSION
We found that a larger percentage of foreign residents and a larger average household size as indicators for lower SES were found in neighborhoods with a lower coverage in the BCSP. However, a higher average personal annual income as an indicator for a higher SES was also found in neighborhoods with a lower coverage in the BCSP. For these neighborhoods, the negative association between SES and coverage in the BCSP was stronger than in the neighborhoods with a relatively middle SES level. The neighborhoods with a higher coverage of opportunistic screening had a relatively higher average personal annual income and the positive association between SES and coverage of the opportunistic screening was stronger than in neighborhoods with a relatively lower SES level. The median coverage in the BCSP for all neighborhoods was 50%. The median coverage for neighborhoods in the lowest and the highest quantile was 33.5 and 64.5%, respectively. The coverage in the BCSP in Flanders was lower than in countries like the Netherlands and the United Kingdom where BC screening is mainly performed in the BCSP, but it was close to countries like France and Germany where both BCSP and opportunistic screening are performed for BC screening (17). As the median coverage of the opportunistic screening was 12%, with 4.2% coverage in the neighborhoods in the lowest quantile and 24.8% coverage in the neighborhoods in the highest quantile, the combined coverage in the BCSP and the opportunistic screening was relatively close to the BC screening coverage in other Western European countries (3,17).
We observed that women in neighborhoods with a lower SES, as indicated by a crowded housing condition (18,19) or being an immigrant with a foreign nationality (20), tended to participate less in the BCSP than women in neighborhoods with middle level SES. The negative association between these variables and the coverage in the BCSP was previously also found in other European countries (21,22). An explanation for this phenomenon is that women with a relatively low SES have a possible lack of health literacy, leaving them less informed on the benefit of BCSP (23). Other reasons may include language barriers for women with an immigration background (24), and the fact that the hardship of life for women with relatively low SES may require more time to work than high SES women and reduce the attention to health care (25).
Interestingly, we also found that women in neighborhoods with a higher income as indicated by a higher average personal annual income also tend to participate less in the BCSP than the women in neighborhoods with middle level income. However, these women tend to participate more in the opportunistic screening than the women in neighborhoods with middle level income. Studies in France, where BCSP and opportunistic screening coexist as in Flanders, also found that women in the most affluent group tend to participate more in opportunistic Numbers in bold: The 10th (Q10) to the 90th (Q90) quantile regression coefficient that was significantly different from zero at the 5% significance level. The 95%CI in parentheses. F statistics below the parentheses. *P < 0.05, **P < 0.01, ***P < 0.001.
Frontiers in Public Health | www.frontiersin.org screening (26). A possible explanation for this phenomenon is that women who belong to the neighborhoods with a relatively high SES may prefer the opportunistic screening which has a more flexible time schedule and more personalized service than the BCSP (2, 27).

STRENGTH AND LIMITATION
The strength of this study is that the validity of the screening coverage data was warranted by the administrative database of the screening program in Flanders. The recall bias in selfreported surveys was avoided in our data (28). Moreover, the determinants of screening coverage were evaluated for the full distribution of the coverage in the BCSP and the opportunistic screening. The heterogeneous effect of the determinants of screening coverage in different quantiles can provide a more specific target for the potential interventions to improve the screening coverage rate. A limitation of this study was that all data were aggregated at a neighborhood level. It was therefore not possible to explore the variation within the neighborhoods due to the aggregated nature of the data. Moreover, the aggregated data in the current data do not contain factors that have a potential impact on the implementation of the screening, such as women's contract information with general practitioners which might play an important role in the implementation of opportunistic screening. However, studies have shown that residents who live in the same neighborhood share similar SES and health behavior (29,30) and the association between the determinants in our study and the screening coverage was consistent with the literature. Therefore, the findings in this study can demonstrate the disproportional impact of SES on women who have different breast cancer screening rate and can inform the resource allocation in policymaking. Another potential limitation of this study was the exclusion of neighborhoods that have less than five screened women, which was mainly related to the small size of a neighborhood. This measure was taken by the breast cancer screening program administrators with the aim to protect privacy and is mandatory by the Flemish government. However, we think that the effect of this exclusion was limited due to the fact that we included 92% of all neighborhoods (31).

CONCLUSION
Women in neighborhoods that have a relatively low SES, that are characterized by being an immigrant with a foreign nationality, and having a large average household size, as well as women in neighborhoods that have a relatively high SES, that are characterized by a high average personal annual income, participate less frequently in the BCSP. On the other hand, women who belong to neighborhoods with a relatively high SES tend to participate more in the opportunistic screening. Tailored intervention that aims to increase the coverage of BCSP should pay more attention to women in these neighborhoods.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
Ethical review and approval was not required for the study on human participants in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

AUTHOR CONTRIBUTIONS
LD, MJG, GB, and GV conceived the study, performed the data analysis, and wrote the manuscript. SJ was involved in the data curation and review of the manuscript. KV, MG, HD, and PM wrote the manuscript and verified the analysis. All authors discussed and interpreted the results and contributed to the final manuscript.