Exploring geographical differences in the incidence of colorectal cancer in the Norwegian Women and Cancer Study: a population-based prospective study
Received 3 March 2019
Accepted for publication 5 July 2019
Published 8 August 2019 Volume 2019:11 Pages 669—682
Checked for plagiarism Yes
Review by Single-blind
Peer reviewer comments 2
Editor who approved publication: Professor Henrik Toft Sørensen
Sunday Oluwafemi Oyeyemi,1 Tonje Braaten,1 Edoardo Botteri,2,3 Paula Berstad,2 Kristin Benjaminsen Borch1
1Department of Community Medicine, Uit-the Arctic University of Norway, Tromsø, Norway; 2Department of Bowel Cancer Screening, Cancer Registry of Norway, Oslo, Norway; 3Norwegian National Advisory Unit for Women’s Health, Women’s Clinic, Oslo University Hospital, Oslo, Norway
Purpose: Norway has experienced an unexplained, steep increase in colorectal cancer (CRC) incidence in the last half-century, with large differences across its counties. We aimed to determine whether geographical distribution of lifestyle-related CRC risk factors can explain these geographical differences in CRC incidence in Norwegian women.
Methods: We followed a nationally representative cohort of 96,898 women with self-reported information on lifestyle-related CRC risk factors at baseline and at follow-up 6–8 years later in the Norwegian Women and Cancer Study. We categorized Norwegian counties into four county groups according to CRC incidence and used Cox proportional hazard models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for risk factors. We used the Karlson, Holm, and Breen (KHB) method of mediation analysis to investigate the extent to which the risk factors accounted for the observed differences in CRC incidence between counties.
Results: During an average of 15.5 years of follow-up, 1875 CRC cases were diagnosed. Height (HR=1.12; 95% CI 1.08, 1.17 per 5 cm increase); being a former smoker who smoked ≥10 years (HR=1.34; 95% CI 1.15, 1.57); or being a current smoker who has smoked for ≥10 years (HR=1.28; 95% CI 1.12, 1.46) relative to never smokers was associated with increased CRC risk. Duration of education >12 years (HR=0.78; 95% CI 0.69, 0.87) vs ≤12 years, and intake of vegetables and fruits >300 g (HR=0.90; 95% CI 0.80, 0.99) vs ≤300 g per day were associated with reduced CRC risk. However, these risk factors did not account for the differences in CRC risk between geographical areas of low and high CRC incidence. This was further confirmed by the KHB method using baseline and follow-up measurements (b=0.02, 95% CI −0.02, 0.06, p=0.26).
Conclusion: Lifestyle-related CRC risk factors did not explain the geographical variations in CRC incidence among Norwegian women. Possible residual explanations may lie in heritable factors.
Keywords: lifestyle, diet, risk factors, colorectal cancer, women, NOWAC study
Colorectal cancer (CRC) is the second most common malignancy in women globally,1 and the second leading cause of cancer-related death in high-income countries.2 Norway has experienced an unexplained, steep increase in the incidence of CRC in both men and women in the last half-century.3,4 From 1957–61 to 2012–16, incidence rates among Norwegian women increased from 21 to 54 per 100,000 person-years for colon cancer, and from 9 to 20 per 100,000 person-years for rectal cancer.5 The CRC incidence rates among women in Norway are currently among the highest in the world,6 having almost tripled from 1957–61 to 2012–16, and surpassing the rates in other Nordic countries with apparently similar lifestyles. So far, the reasons for this steep increase have been elusive. Moreover, differences in CRC incidence vary over 10-fold across countries,7 which may be ascribed to variations in dietary and environmental exposures, coupled with genetic susceptibility.8 CRC incidence also varies within Norway, with a more than 20 per 100,000 person-years difference between areas of high and low CRC incidence.9,10 The factors responsible for this geographical heterogeneity are yet to be determined, and knowledge of these factors could be useful to guide screening strategies and health policy.
Therefore, this study aimed to determine whether the geographical distribution of lifestyle-related CRC risk factors can explain the geographical differences in CRC incidence, using the Norwegian Women and Cancer (NOWAC) Study.
Materials and methods
The NOWAC Study is a nationwide, representative prospective cohort study which started in 1991.11 The full detail of the cohort profile has been described previously.11,12 Summarily, the study consists of over 172,000 women who were recruited over three different time periods: 1991–92, 1996–97, and 2003–04. Potential participants aged 30–70 years were randomly selected from the Norwegian Central Population Register (Statistics Norway) and received a questionnaire by mail that collected information on their lifestyle and health status at enrollment (baseline questionnaire). Similar follow-up questionnaires were sent to the same women about 6–8 years later. All women who agreed to participate completed and returned the questionnaires with written informed consent. The NOWAC Study was approved by the Regional Committee for Medical Research Ethics and the Norwegian Data Inspectorate.11
NOWAC participants who were enrolled in 1991–92, 1996–97, and 2003–04 and completed a food frequency questionnaire (FFQ) in 1998, 1996–97, and 2003–04, respectively, were eligible for inclusion in the present study. Those who were enrolled in 1991–92 completed an FFQ in 1998 because an FFQ was not included in the 1991–92 questionnaire. Thus, we used the 1998 information as baseline for the participants enrolled in 1991–92. This represented 101,321 participants who completed a baseline questionnaire with dietary information between 1996 and 2004. We subsequently excluded women who died or emigrated (n=14) prior to the start of follow-up, and all cases of prevalent cancer except non-melanoma skin cancer (n=4,414). This resulted in a final study sample of 96,893 women. Follow-up information was available for 68,626 (70.8%) of these women.
Assessment of CRC risk factors
Information on age, physical activity, height, weight, duration of education, alcohol intake, smoking status and intensity (pack-years), annual household income, hormone replacement therapy use, oral contraceptive use, and dietary habits (daily intake of red meat, processed meat, fish, fruits and vegetables, fiber, calcium, vitamin D, and milk) were taken from the NOWAC questionnaire. Physical activity was reported on a validated 10-point scale, on which 1 was “very low” and 10 was “very high”. This is a global (ie, all-inclusive) physical activity score that has been found valid to rank the physical activity of women in the NOWAC Study.13 The validated, self-reported height and weight measurements from the questionnaires were used to compute body mass index (BMI).14 Information on the duration of education and alcohol intake was obtained from the questionnaire, while information on smoking status and smoking intensity (pack-years) were combined into one variable of smoking history. Information on annual household income, hormone replacement therapy use, and oral contraceptive use were also extracted from the NOWAC questionnaire. The FFQ includes foods that are common in Norway and has been validated.15,16
Assessment of county of residence and creation of county groups by CRC incidence
County of residence at baseline was accessed through linkage to the Norwegian Central Population Register (Statistics Norway). There were 19 counties in Norway at the time of data collection (Figure 1). We used percentiles of CRC incidence rate (Table 1) to categorize the counties into four groups. The intent was to compare the lowest 10% to the highest 10% to discern possible differences in lifestyle-related CRC risk factors. However, we raised the limit of the low-incidence counties to the 15th percentile to allow for more cases of CRC in this group. Thus, we grouped counties from 0 to 15th percentile as low-incidence counties (Oppland, Sør Trøndelag, and Telemark); 15–50th as mid-low-incidence counties (Hedmark, Hordaland, Oslo, Møre and Romsdal, Nord-Trøndelag, Vest-Agder, and Buskerud); 50–90th as mid-high-incidence counties (Rogaland, Akershus, Aust-Agder, Vestfold, Østfold, Finnmark, and Troms); and 90–100th as high-incidence counties (Nordland, Sogn and Fjordane).
Figure 1 Map of Norway showing the 19 counties and regions.
Table 1 Basic parameters and endpoints in the 19 counties of Norway in the Norwegian Women and Cancer Study
We also conducted sensitivity analyses in which, we grouped participants by region of residence (Oslo, East, South, West, Middle, and North)18 and by rural/urban area of residence. Urban residence was defined as living in a “dense area” with a maximum distance of 50 m between houses, except for public areas or natural barriers, and inhabited by at least 200 persons.19
CRC incidence, emigration, and death
Participants diagnosed with primary colon or rectal cancer were ascertained through linkage to the Cancer Registry of Norway. We used the International Statistical Classification of Diseases and Related Health Problems, Tenth Edition (ICD-10), which uses code C18 for colon and C19-20 for rectal cancer. The county of residence, date of emigration, and date of death were ascertained via linkage to the Norwegian Central Population Register (Statistics Norway).
We carried out an initial analysis using the baseline data to assess the CRC risk factors for multi-collinearity. This initial analysis included height (continuous, in meters); physical activity (dichotomized into inactive (1–5) and active (6–10)); BMI (<20.0, 20.0–24.9, 25.0–29.9, and ≥30.0 kg/m2); duration of education (≤12 and >12 years); alcohol intake (0, ≤3.0, >3.0–10.0, and >10.0 g/day); smoking history (never, former smoker of <10 years, former smoker of ≥10 years, current smoker of <10 years, current smoker of ≥10 years); annual household income in Norwegian kroner (NOK) (low: <300,000 NOK, medium: 300–600,000 NOK, and high: >600,000 NOK); hormone replacement therapy use (never/ever); and oral contraceptive use (never/ever). All the dietary variables were dichotomized along their median values: red meat intake (0, ≤15, >15 g/day); processed meat intake (0, ≤70, >70 g/day); fish intake (0–90, >90 g/day); fruit and vegetable intake (0–300, >300 g/day); fiber (0–21, >21 g/day); calcium intake from food (0–700, >700 mg/day); vitamin D intake (0–6, >6 µg/day); and milk intake (0, ≤170, >170 g/day). Where possible, we used the median values (50th percentile) to split the variables into categories, as the median values are more robust and undistorted by outliers.20
We present descriptive statistics at baseline as mean values (±standard errors, SEs) or percentages. We used Cox proportional hazard regression models with age as the time scale to estimate the hazard ratios (HRs) and 95% confidence intervals (CIs) for the associations between the county groups (low-, mid-low-, mid-high-, and high-incidence counties), risk factors, and CRC incidence. Follow-up time was defined as the period in years between age at baseline and age at diagnosis of incident cancer, death, emigration, or age at the end of follow-up (31 December 2016), whichever came first.
We assessed predefined possible interaction effects between physical activity versus BMI, smoking history, alcohol intake, and dietary factors, respectively. We also checked for interaction effects between duration of education and BMI, smoking history, alcohol intake, and dietary factors, respectively. We tested for multi-collinearity between calcium versus milk and vitamin D intake, respectively; red meat versus processed meat intake; and fiber versus fruit and vegetable intake. We excluded milk because of high collinearity with calcium and >25% missing values in the variable. We repeated the baseline analyses following exclusion of cancers diagnosed in the first 2 years of follow-up to control for possible reverse causality. Sensitivity analyses were carried out by region of residence, and area of residence (rural/urban).
Mediation analysis using Karlson, Holm, and Breen (KHB) method of decomposition
We used the KHB method of mediation analysis21 to investigate the extent to which the CRC risk factors (mediating variables) account for the observed difference in CRC incidence between individual counties. The KHB method provides decomposition of the total effects of counties on CRC incidence into direct and indirect effects.21 The basic outputs from the KHB method include three models: the reduced model, the full model, and the difference (model). The reduced model describes the estimated effect of the counties with no mediating variables in the model (total effect). The full model describes the estimated effect of counties with all mediating variables in the model (direct effect). The difference between these two models represents the indirect effect. The indirect effect is interpreted as the mediation effect. The KHB method assumes a normal distribution of the indirect effect, and this assumption has been shown to be legitimate in large samples such as the NOWAC Study.22 We fitted the KHB models using the data collected at baseline and then used the multiply imputed data.
Multiple imputation and repeated measurements analyses
Multiple imputation using chained equations was used to handle missing data, under the assumption that this data was missing at random.23 The missing values were replaced by multiply imputed values from 20 duplicate datasets. We created 20 duplicates datasets from the imputation simulation to reduce sampling variability.24 We included all the CRC risk factors used in the analyses and the Nelson–Aalen cumulative hazard estimator as predictors in the imputation model.25,26 We used Rubin’s rules to combine the estimates from the 20 imputed datasets to estimate HRs and corresponding 95% CIs.27 The KHB method also computes the total, direct, and indirect effects for each imputed dataset and combines the estimates using Rubin’s rules.
We used baseline information up to the point when follow-up information was available on physical activity, BMI, alcohol intake, smoking history, hormone replacement therapy use, and all dietary intakes. We then used the follow-up information until death, emigration, or the end of the study, whichever occurred first.
All the analyses and multiple imputations were done in Stata version 15.0 (StataCorp, College Station, TX, USA). Figure 1 is produced using GraphPad Prism 8 (GraphPad Software, San Diego, CA). All statistical analyses were two-sided, and p-values were considered statistically significant at a level of <0.05.
During an average of 15.5 years of follow-up and 1.5 million person-years, 1875 CRC cases (1276 [68%] colon cancers and 599 [32%] rectal cancers) were diagnosed in the study sample. The counties of lowest and highest crude incidence rates were Oppland, and Sogn and Fjordane, respectively (Table 1).
The median age at baseline was 51 years, while the median age at diagnosis of CRC was 66 years (range 43–89). When looking at county groups, low-incidence counties had a higher proportion of physically active women compared to high-incidence counties (46% vs 41%) at baseline. Similarly, the low-incidence counties had a higher proportion of women with a longer duration of education (38% vs 25%), never smokers (38% vs 34%), high annual household income (12% vs 5%), hormone replacement therapy use (34% vs 30%), and oral contraceptive use (53% vs 43%), compared to high-incidence counties. Conversely, high-incidence counties had higher proportion of women with overweight (33% vs 31%), obese (10% vs 9.6%), ever smokers (64% vs 60%), and low annual household income (48% vs 36%), compared to low-incidence counties (Table 2).
Table 2 Selected participant characteristics by county group at study enrollment (baseline) in the Norwegian Women and Cancer Study
The variables with the highest proportion of missing values at baseline were physical activity (9.5%), annual household income (7.3%), and duration of education (5.8%). At follow-up, 38% of the women had missing values on physical activity, and approximately 30% had missing information on BMI, alcohol intake, smoking history, hormone replacement therapy use, and dietary intakes. There was no substantial change in the characteristic features of the study sample between the imputed and the complete-case dataset (Table S1).
The multivariable-adjusted model of repeated measurements showed that the high-incidence county group had an HR of 1.37 (95% CI 1.13–1.66) relative to the low-incidence county group (Figure 2), which was similar to the unadjusted estimate (Table S2). Height (HR=1.12; 95% CI 1.08, 1.17 per 5 cm increase), being a former smoker who smoked ≥10 years (HR=1.34; 95% CI 1.15, 1.57), or a current smoker who had been smoking ≥10 years (HR=1.28; 95% CI 1.12, 1.46), compared to never smokers, were significantly associated with a higher CRC risk. Duration of education >12 years (HR=0.78; 95% CI 0.69, 0.87) compared to ≤12 years, and daily fruit and vegetable intake >300 g (HR=0.90; 95% CI 0.80, 0.99) compared to ≤300 g, were associated with decreased CRC risk (Figure 2).
Figure 2 Multivariable hazard ratios (HRs) and 95% confidence intervals (CIs) of factors associated with colorectal cancer (CRC) incidence at baseline and follow-up with chained multiple imputations, in the Norwegian Women and Cancer study.
No substantial difference was seen after excluding those who were diagnosed with CRC during the first 2 years of follow-up (data not shown). Sensitivity analyses by region showed no differences in the HR estimates for CRC risk factors, nor were any statistically significant differences seen in the HR estimates for the regions before and after multivariable adjustment. This was also the case in sensitivity analyses that used rural/urban area of residence (Table S2).
The KHB analysis showed the extent to which the mediating variables (CRC risk factors) account for the difference in CRC incidence between the low-incidence county group (reference) and that of other county groups. At baseline, the log odds of having CRC in the high-incidence county group were 0.41 higher than those in the low-incidence county group (Table 3). After adjusting for mediating factors, the effect of living in the high-incidence county group reduced to 0.39, leaving an indirect effect of 0.02 (b=0.02; 95% CI −0.02, 0.06, p=0.26). This shows that the differences in CRC incidence between the low- and high-incidence county groups are not significantly mediated by the combined effects of the investigated CRC risk factors (Table 3). The mediation analysis results in the imputed dataset were similar to the baseline results. We conducted a sensitivity analysis using the 19 counties individually (without grouping), which also showed that the combined effects of the risk factors did not significantly mediate the variations in CRC incidence across counties (data not shown).
Table 3 Decomposition of total effects of county groups into direct and indirect effects using the Karlson, Holm, and Breen method at baseline and follow-up in the Norwegian Women and Cancer Study
In this large cohort of Norwegian women, we found that county-level differences in CRC incidence were not explained by differences in lifestyle-related CRC risk factors. This was demonstrated by two different approaches: Cox proportional hazards models and the relatively new KHB method of decomposition.
The lifestyle-related CRC risk factors significantly associated with CRC incidence in our cohort of women included height, smoking history, duration of education, and fruit and vegetable intake. Our results showed that these factors, together with other CRC risk factors, did not significantly explain the differences in the CRC incidence between the counties. CRC risk in county groups remained statistically the same before and after adjusting for risk factors. These results remained consistent when using baseline data, as well as when using repeated measurements with multiple imputation. Our findings suggest that there are other important or unmeasured risk factors that are responsible for the differences in CRC incidence between Norwegian counties.
Previous international studies have rationalized that variations in CRC incidence in different areas of a country are due to different, but overlapping, contributory factors, such as rural–urban disparities, socioeconomic status (SES), ease of access to health care, public health campaigns, unique social and lifestyle risk factors, differences in exposure to risk factors, such as in dietary customs and ethnic variations in food preparation, and different exposures to unknown risk factors.28–31 Some studies have indicated that rural–urban disparities confer an increased risk of CRC in rural areas32,33 and suggested that the relationship may be mediated through screening behavior.32,33 Other studies have reported that the increased risk may simply reflect the socioeconomic differences between rural and urban communities.34 Other studies found a higher risk of CRC in urban areas.34–36 These findings differ by country and time period of assessment, and differences in the definition of rural/urban areas may mask the relationship between this variable and CRC risk.35 There is currently no national CRC screening program in Norway, which could expound on some of the geographical differences in the present population.
Education and household income are often used as proxy indicators of SES. We found a significant inverse association between duration of education and CRC risk, while we found no such association with annual household income. Results of previous similar studies regarding SES have been inconsistent. A recent review showed that, in the United States and Canada, low SES groups have a higher CRC incidence than high SES groups (RR from 1.0 to 1.5), while these findings were mostly reversed (RR from 0.3 to 0.9) in Europe.30 Nonetheless, education, and not necessarily income, may be a better predictor of a healthy lifestyle.37,38
Cigarette smoking has been associated with increased incidence of CRC, and our data further suggest that the risk remains even among former smokers. A meta-analysis of 106 observational studies concluded that smokers have an increased risk of developing CRC compared to never smokers (RR 1.18, 95% CI 1.11–1.25).39 Height was also associated with increased CRC risk in our study sample. This finding is in agreement with two recent systematic reviews of prospective studies, which posited a potential causal association of adult attained height with the risk of CRC.40,41 Our study found a significant inverse association between fruit and vegetable intake and CRC risk, which is in concurrence with the findings in the European Prospective Investigation into Cancer and Nutrition (EPIC) study.42
In our study, participants in the low-incidence county group were more physically active, had a longer duration of education, were more often never smokers, and had a higher fruit and vegetable intake. These are markers of a generally healthy lifestyle, and the reduced CRC risk observed in this county group may be a reflection of this lifestyle. Notwithstanding, these factors failed to account for the risk differences between low- and high-incidence county groups.
Occurrence of exposure to established risk factors for cancer has been reported to vary geographically within some countries. For instance, the prevalence of obesity varies within Finland,43 while the use of hormone replacement therapy is more likely in women living in urban areas of Denmark.44 Therefore, it is plausible that the risk of CRC could vary in different counties or areas due to different prevalences of exposure to established CRC risk factors. However, since these established risk factors did not account for the observed risk differences in CRC between the counties in the present study, considerable uncertainty remains about what is responsible for these differences. This may be a partial reflection of the incomplete understanding of the carcinogenesis of CRC,34 although the unexplained risk differences could also come from unmeasured risk factors. A large Scandinavian study, which combined cohorts of twins from Sweden, Denmark, and Finland, demonstrated that genetically inheritable factors account for 35% of the CRC cases, while non-shared environmental factors account for 60%, and shared environmental factors the remaining 5%.45 Thus, a possible explanation for our observed differences in risk between high- and low-incidence county groups probably lies more in genetically inherited factors. The well-described CRC-related inheritable syndromes (such as hereditary nonpolyposis colon cancer (HNPCC) and familial adenomatous polyposis (FAP)), where inheritance is highly penetrant, only account for about 3–5% of the inherited cases of CRC.46
The main limitations of this study are the unmeasured established CRC risk factors. This includes family history of CRC and its precursors (such as adenomatous polyps), as genetically inherited factors can increase the likelihood of CRC oncogenesis.45,46 Our study lacks information on the use of aspirin and other non-steroidal anti-inflammatory drugs, the regular use of which has been associated with reduced CRC risk.47,48 The lack of information on these factors may have confounded our study. The county of residency used in this study was captured only at baseline; thus, some of the participants could have changed their county of residence in the course of the study. However, most women at the age of our cohort would have settled down at a county on a long-term basis. We lack the power to explore the CRC risk in each county or in each county group separately. Most variables in our study are self-reported and therefore are saddled with the errors inherent with self-reported measurements. However, most of these variables, such as physical activity, duration of education, BMI, alcohol intake, and dietary habits, have been validated with good results.12–16
The strengths of our study include the prospective and population-based design, with a large sample size of participants who were randomly recruited and are representative of Norwegian women between 30 and 70 years at recruitment,12 information on important risk factors, and the high quality of the national cancer registry with almost 100% completeness.49 The NOWAC Study has been shown to have almost the same observed cumulative incidence rates for all cancer sites as that of the national figures.11,12 We used repeated measurements of variables to account for changes in these variables over time in order to lower the risk of measurement error. We used chained multiple imputation to deal with missing data, and thus maximize the number of participants, and by extension, the number of CRC cases included in the analyses.
The lifestyle-related CRC risk factors that we investigated did not account for the risk differences between the areas of low and high incidence of CRC. A possible explanation lies in inheritable factors. Thus, the family history of CRC cases may be especially important in determining the appropriate preventive screening strategy in areas of high incidence.
CRC, colorectal cancer; NOWAC, Norwegian Women and Cancer Study; FFQ, Food frequency questionnaire; BMI, body mass index; NOK, Norwegian kroner; KHB method, Karlson, Holm, and Breen method.
Ethical approval and informed consent
The Norwegian Women and Cancer Study was approved by the Regional Committee for Medical Research Ethics and the Norwegian Data Inspectorate (P REK NORD 141/2008). All participants gave written informed consent.
To access the data supporting the findings presented, kindly contact the person in charge of the NOWAC Study - https://site.uit.no/nowac/contact-information/.
We thank the staff and participants in the Norwegian Women and Cancer Study for the time and effort they put forth to advancing medical science. Many thanks to Rolf Wynn, Dolley Charles, Runa Borgund Barnung, and Jan Håkon Rudolfsen for their constructive critiques during the work on this study. The publication of this manuscript was funded by the publication fund of UiT-The Arctic University of Norway. The funder had no role in the study design, data analysis, preparation of the manuscript, and decision to publish.
SOO and KBB conceived the study idea. All authors contributed to the data analysis. SOO organised the writing and wrote the initial draft. All authors contributed toward drafting and critically revising the article, gave final approval of the version to be published, and agreed to be accountable for all aspects of the work.
The authors report no conflicts of interest in this work.
2. WHO. World Health Organisation 2018. Top 10 causes of deaths in high-income countries in 2016. 2018. Available from: https://www.who.int/en/news-room/fact-sheets/detail/the-top-10-causes-of-death.
3. NORDCAN. Association of the Nordic Cancer Registries 2017. Cancer statistical fact sheet: Colorectal cancer in Norway. 2017. Available from: http://www-dep.iarc.fr/NORDCAN/English/StatsFact.asp?cancer=590&country=578.
4. Larsen IK, Bray F. Trends in colorectal cancer incidence in Norway 1962–2006: an interpretation of the temporal patterns by anatomic subsite. Int J Cancer. 2010;126(3):721–732.
5. Cancer Registry of Norway: facts on colorectal cancer “Tykk- og endetarmskreft”. 2017. Available from: https://www.kreftregisteret.no/Generelt/Om-kreft/Tykk-og-endetarmskreft/.
7. Global Burden of Disease Cancer C, Fitzmaurice C, Allen C, et al. Global, regional, and national cancer incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 cancer groups, 1990 to 2015: a systematic analysis for the Global Burden of Disease Study. JAMA Oncol. 2017;3(4):524–548. doi:10.1001/jamaoncol.2016.5688.
8. Finlay A, Macrae M. Colorectal cancer: epidemiology, risk factors, and protective factors. UpToDate. 2018. Available from: https://www.uptodate.com/contents/colorectal-cancer-epidemiology-risk-factors-and-protective-factors.
9. Cancer Registry of Norway: women in Western Norway are most affected by bowel cancer “Kvinner på Vestlandet rammes oftest av tarmkreft”. 2017. Available from: https://www.kreftregisteret.no/Generelt/Nyheter/kvinner-pa-vestlandet-rammes-oftest-av-tarmkreft/.
10. Larsen I, Møller B, Johannesen TB, et al. Cancer in Norway 2016 - Cancer Incidence, Mortality, Survival and Prevalence in Norway. Oslo: Cancer Registry of Norway; 2017.
11. Lund E, Dumeaux V, Braaten T, et al. Cohort profile: the Norwegian Women and Cancer Study—NOWAC—Kvinner og kreft. Int J Epidemiol. 2008;37(1):36–41. doi:10.1093/ije/dym137
12. Eiliv L, Merethe K, Tonje B, et al. External validity in a population-based national prospective study – the Norwegian Women and Cancer Study (NOWAC). Cancer Causes Control. 2003;14(10):1001–1008.
13. Borch KB, Ekelund U, Brage S, Lund E. Criterion validity of a 10-category scale for ranking physical activity in Norwegian women. (Research) (Report). Int J Behav Nutr Phys Act. 2012;9:2. doi:10.1186/1479-5868-9-2
14. Skeie G, Mode N, Henningsen M, Borch KB. Validity of self-reported body mass index among middle-aged participants in the Norwegian Women and Cancer study. Clin Epidemiol. 2015;7:313–323. doi:10.2147/CLEP.S83839
15. Hjartaker A, Andersen LF, Lund E. Comparison of diet measures from a food-frequency questionnaire with measures from repeated 24 hr dietary recalls. The Norwegian Women and Cancer Study. Public Health Nutr. 2007;10(10):1094–1103. doi:10.1017/S1368980007702872
16. Parr CL, Veierod MB, Laake P, Lund E, Hjartaker A. Test-retest reproducibility of a food frequency questionnaire (FFQ) and estimated effects on disease risk in the Norwegian Women and Cancer Study (NOWAC). Nutr J. 2006;5:4. doi:10.1186/1475-2891-5-4
17. Brenner H, Chen C. The colorectal cancer epidemic: challenges and opportunities for primary, secondary and tertiary prevention. Br J Cancer. 2018;119(7):785–792. doi:10.1038/s41416-018-0264-x
18. Statista-The Statistics Portal. Population of Norway 2017, by region. 2017. Available from: https://www.statista.com/statistics/586908/population-in-norway-by-region/.
19. Statistics Norway. Population and land area in urban settlements. 2017. Available from: https://www.ssb.no/en/befolkning/statistikker/beftett/aar.
20. Manikandan S. Measures of central tendency: median and mode. J Pharmacol Pharmacother. 2011;2(3):214–215. doi:10.4103/0976-500X.83300
21. Kohler U, Karlson K, Holm A. Comparing Coefficients of Nested Nonlinear Probability Models. The Stata Journal 2011;11.
22. MacKinnon DP, Lockwood CM, Hoffman JM, West SG, Sheets V. A comparison of methods to test mediation and other intervening variable effects. Psychol Methods. 2002;7(1):83–104.
23. White IR, Carlin JB. Bias and efficiency of multiple imputation compared with complete-case analysis for missing covariate values. Stat Med. 2010;29(28):2920–2931. doi:10.1002/sim.3944
24. Sterne JAC, White IR, Carlin JB, et al. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393–b2393. doi:10.1136/bmj.b902
25. White IR, Royston P. Imputing missing covariate values for the Cox model. Stat Med. 2009;28(15):1982–1998. doi:10.1002/sim.3618
26. White IR, Royston P, Wood AM. Multiple imputation using chained equations: issues and guidance for practice. Stat Med. 2011;30(4):377–399. doi:10.1002/sim.4067
27. Rubin DB. Multiple imputation after 18+ years. J Am Stat Assoc. 1996;91(434):473–489. doi:10.1080/01621459.1996.10476908
28. Singh H, Nugent Z, Decker K, Demers A, Samadder J, Torabi M. Geographic variation and factors associated with colorectal cancer incidence in Manitoba. 2018. 2018;108(5–6):7.
29. Brooke HL, Talbäck M, Martling A, Feychting M, Ljung R. Socioeconomic position and incidence of colorectal cancer in the Swedish population. Cancer Epidemiol. 2016;40:188–195. doi:10.1016/j.canep.2016.01.004
30. Aarts MJ, Lemmens VEPP, Louwman MWJ, Kunst AE, Coebergh JWW. Socioeconomic status and changing inequalities in colorectal cancer? A review of the associations with risk, treatment and outcome. Eur J Cancer. 2010;46(15):2681–2695. doi:10.1016/j.ejca.2010.04.026
31. Doubeni CA, Laiyemo AO, Major JM, et al. Socioeconomic status and the risk of colorectal cancer: an analysis of more than a half million adults in the National Institutes of Health-AARP Diet and Health Study. Cancer. 2012;118(14):3636–3644. doi:10.1002/cncr.26677
32. Zahnd WE, James AS, Jenkins WD, et al. Rural–urban differences in cancer incidence and trends in the United States. Cancer Epidemiol Biomarkers Prev. 2018;27(11):1265. doi:10.1158/1055-9965.EPI-17-0430
33. Kinney AY, Harrell J, Slattery M, Martin C, Sandler RS. Rural-urban differences in colon cancer risk in blacks and whites: the North Carolina Colon Cancer Study. J Rural Health. 2006;22(2):124–130. doi:10.1111/j.1748-0361.2006.00020.x
34. Sharp L, Donnelly D, Hegarty A, et al. Risk of several cancers is higher in urban areas after adjusting for socioeconomic status. Results from a two-country population-based study of 18 common cancers. J Urban Health. 2014;91(3):510–525. doi:10.1007/s11524-013-9846-3
35. Hawley ST, Chang S, Risser D, Zhang Q. Colorectal cancer incidence and mortality in Texas 1990–1992: a comparison of rural classifications. J Rural Health. 2002;18(4):536–546.
36. Coughlin SS, Richards TB, Thompson T, et al. Rural/nonrural differences in colorectal cancer incidence in the United States, 1998–2001. Cancer. 2006;107(S5):1181–1188. doi:10.1002/cncr.22015
37. Hahn RA, Truman BI. Education improves public health and promotes health equity. Int J Health Serv. 2015;45(4):657–678. doi:10.1177/0020731415585986
38. Lutz W, Kebede E. Education and health: redrawing the preston curve. Popul Dev Rev. 2018;44(2):343–361. doi:10.1111/padr.12141
39. Botteri E, Iodice S, Bagnardi V, Raimondi S, Lowenfels AB, Maisonneuve P. Smoking and colorectal cancer: a meta-analysis. JAMA. 2008;300(23):2765–2778. doi:10.1001/jama.2008.839
40. Khankari NK, Shu X-O, Wen W, et al. Association between Adult Height and Risk of Colorectal, Lung, and Prostate Cancer: results from Meta-analyses of Prospective Studies and Mendelian Randomization Analyses. PLoS Med. 2016;13(9):e1002118–e1002118. doi:10.1371/journal.pmed.1002118
41. Abar L, Vieira AR, Aune D, et al. Height and body fatness and colorectal cancer risk: an update of the WCRF-AICR systematic review of published prospective studies. Eur J Nutr. 2018;57(5):1701–1720. doi:10.1007/s00394-017-1557-1
42. van Duijnhoven FJB, Bueno-De-Mesquita HB, Ferrari P, et al. Fruit, vegetables, and colorectal cancer risk: the European Prospective Investigation into Cancer and Nutrition. Am J Clin Nutr. 2009;89(5):1441–1452. doi:10.3945/ajcn.2008.27120
43. Lahti-Koski M, Taskinen O, Similä M, et al. Mapping geographical variation in obesity in Finland. Eur J Public Health. 2008;18(6):637–643. doi:10.1093/eurpub/ckn089
44. Løkkegaard E, Lidegaard Ø, Nørgaard Møller L, Agger C, Helms Andreasen A, Jørgensen T. Hormone replacement therapy in Denmark, 1995–2004. Acta Obstet Gynecol Scand. 2007;86(11):1342–1351. doi:10.1080/00016340701505523
45. Lichtenstein P, Holm NV, Verkasalo PK, et al. Environmental and heritable factors in the causation of cancer — analyses of cohorts of twins from Sweden, Denmark, and Finland. N Engl J Med. 2000;343(2):78–85. doi:10.1056/NEJM200007133430201
46. Burt R. Inheritance of colorectal cancer. Drug Discov Today. 2007;4(4):293–300. doi:10.1016/j.ddmec.2008.05.004
47. Rothwell PM, Wilson M, Elwin C-E, et al. Long-term effect of aspirin on colorectal cancer incidence and mortality: 20-year follow-up of five randomised trials. Lancet. 2010;376(9754):1741–1750. doi:10.1016/S0140-6736(10)61543-7
48. Giovannucci E, Egan KM, Hunter DJ, et al. Aspirin and the risk of colorectal cancer in women. N Engl J Med. 1995;333(10):609–614. doi:10.1056/NEJM199509073331001
49. Larsen IK, Småstuen M, Johannesen TB, et al. Data quality at the Cancer Registry of Norway: an overview of comparability, completeness, validity and timeliness. Eur J Cancer. 2009;45(7):1218–1231. doi:10.1016/j.ejca.2008.10.037
Table S1 Comparison of the complete-case and imputed dataset, the Norwegian Women and Cancer study
Table S2 Hazard ratios (HRs) and 95% confidence intervals (CIs) before and after multivariable risk adjustment at baseline (complete-case analysis) in the Norwegian Women and Cancer study
This work is published by Dove Medical Press Limited, and licensed under a Creative Commons Attribution License. The full terms of the License are available at http://creativecommons.org/licenses/by/4.0/. The license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.Download Article [PDF]