- Research article
- Open Access
- Open Peer Review
Milk consumption in relation to incidence of nasopharyngeal carcinoma in 48 countries/regions
BMC Cancervolume 15, Article number: 994 (2015)
Decreasing trends of nasopharyngeal carcinoma (NPC) incidence have been consistently reported in endemic populations but the etiology of NPC remains unclear. The objective of our study was to assess the international and local (Hong Kong) correlations of milk and dairy products per capita consumption with NPC incidence.
We conducted an ecological study in 48 countries/regions. Age standardized incidence rates of NPC were obtained from the Cancer Incidence in Five Continents. Dairy product consumption and Human Development Index were obtained from the Food and Agriculture Organization of the United Nations and the United Nations Development Programme. Spearman correlation, multivariate analysis and time-lagged analysis were performed.
The negative correlations between milk consumption and decreased age standardized incidence rates of NPC were observed in the 48 countries/regions adjusting for Human Development Index in endemic countries/regions. In Hong Kong, multivariate analysis, after adjusting for other potential confounders, including salted fish, cigarette, vegetable consumption and socioeconomic status, showed consistently negative and significant correlations between milk consumption and NPC incidence (The strongest coefficient (β) was observed at 10-year lag in males [β = −0.439; P < 0.01] and in females [β = −0.258; P < 0.01]).
Our study showed the correlations on milk consumption per capita and against lower risk of NPC in 48 countries/regions and in Hong Kong. These hypothesis-generating results could support further studies on individual exposures and the disease.
Nasopharyngeal carcinoma (NPC) is rare worldwide, but much more common in Southeast Asia . One of the highest age-standardized incidence rates (ASIRs) of NPC was observed in Hong Kong (12.5 per 105 person-years in males, 2012) , and similar ASIRs were observed in nearby cities including Zhongshan and Guangzhou, southern China . By contrast, the ASIRs for most part of the world are less than 0.5 per 105 person-years . In addition to the marked geographic incidence pattern of NPC, decreasing trends of NPC incidence have been consistently reported in endemic populations . In Hong Kong, for example, NPC was the 8th most common cancer in 2012; however, the ASIR in males and females combined decreased by 68 % during 1983 to 2012 . The endemic regions of NPC were marked by rapid economic growths and these growths might be associated with the observed declines , but the underlying causes are unclear.
Economic development may have a link with the changes of lifestyle, including decreased consumption of preserved foods and increased consumption of various food . Cantonese-style salted fish was a staple preserved foods in southern China because they were cheap. It is considered as a major NPC risk factor and has been rated as a Group 1 carcinogen by the International Agency for Research on Cancer . However, our recent ecological study did not find a correlation between the decline in salted fish consumption per capita and the corresponding changes in NPC incidence in 8 regions with varied NPC risk, including Hong Kong which has the highest risk . The decreasing trend of NPC incidence may be due to other risk or protective factors.
The World Cancer Research Fund suggests that approximately 30–40 % of cancer cases are potentially preventable through changes in food consumption patterns . Milk consumption has greatly increased in Asia, though is still far lower than Western countries . Differences in milk consumption and in the incidence of NPC among populations worldwide are distinct. The limited evidence of milk consumption and NPC came from case-control studies, and majority of these in endemic regions were conducted two decades ago before the changes of dietary patterns with rapid economic development. The correlation between milk consumption and NPC remains unclear, and new evidence is needed. Case-control or cohort studies take some years to deliver results and should be conducted to test the hypothesis generated by ecological studies. Ecological study is the most cost-effective approach to assess the population disease impact from changes over several decades in the consumption of food items that are suspected to be harmful or protective, and can provide quick albeit preliminary evidence at low cost. A notable example is that the hypothesis of the correlation between higher milk consumption and increased risk of prostate cancer was generated from ecological studies , and has been later confirmed by further studies (case-control and cohort studies, and meta-analyses) which provide stronger evidence .
We therefore assessed the international and local (Hong Kong) correlations of milk and dairy products per capita consumption with NPC incidence.
A total of 48 countries/regions with comprehensive and reliable cancer and food statistics data were included.
NPC incidence statistics
ASIR for NPC (C11, ICD-10: The International Classification of Diseases 10th Revision) were obtained from Cancer Incidence in Five Continents, published by the International Agency for Research on Cancer, 1998–2002 (Vol. IX) . All incidence data were adjusted to the Segi standard population . The data collection procedures are shown in Appendix 1.
Milk and dairy product consumption per capita in 48 countries/regions
Per capita milk consumption data from 1968 to 2002 were obtained from two international sources. Food supply data, including those for “milk (whole)”, “cheese”, “butter and ghee”, and “fresh cream”, were obtained from the Food Supply database maintained by the Food and Agriculture Organization of the United Nations (FAO) (http://faostat.fao.org/). Milk excluding butter per capita consumption was estimated from import, export, and change in stock volumes provided in the Food Balance Sheet of the FAO according to the formula in the International Farm Comparison Network dairy report 2004 [(import − export − change in stock) ÷ total population of the different countries/regions]. The 1950–2010 population data were obtained from the World Population Prospects: the 2012 Revision (http://esa.un.org/wpp/unpp/panel_population.html) .
Indicators of the social economic status of the 48 countries/regions were obtained from the Human Development Index (HDI) published by the United Nations, which includes national income, education and life expectancy.
Salted fish, tobacco, vegetables consumption per capita, and Gross Domestic Product (GDP) per capita in multivariate analysis
Hong Kong has the highest ASIR among the 48 countries/regions, and has data on other risk factors. The data collection procedures have been described previously . Briefly, per capita consumption of salted fish, tobacco and vegetables was estimated using the Service Centre on Trade Statistics, Census and Statistics Department, Hong Kong Special Administrative Region (HKSAR) Government, the Hong Kong Council on Smoking and Health via the General Household Survey and Thematic Household Survey 2008, and the Department Annual Report of the Agriculture, Fisheries and Conservation Department (procedures included in Appendix 2). The per capita GDP was obtained from the 2009 Gross Domestic Product (Yearly), Census and Statistics Department, HKSAR Government.
Spearman correlation analysis was used to examine the correlation between the ASIR of NPC and milk and dairy product consumption per capita (1998–2002). Six additional sets of correlation were calculated for consumption with the following time lags: 5-year (1993–1997), 10-year (1988–1992), 15-year (1983–1987), 20-year (1978–1982), 25-year (1973–1977) and 30-year (1968–1972).
Regression analysis was used to examine the correlation of milk consumption with the log NPC incidence adjusting for HDI in order to comply with the assumptions underlying outliers. In our analysis,
where β0 is the intercept, β1 and β2 represent parameters in the multiplicative model for each milk and dairy products, and (Milk & dairy product) i denotes each type of milk and dairy product, i from 1 to 5.
Each histopathological diagnosis of NPC has its unique epidemiologic profile: keratinizing squamous carcinoma is the main histological subtype in low-risk areas, while non-keratinizing undifferentiated carcinoma is the predominant subtype in endemic areas . Countries/regions were stratified by low- and high-risk countries/regions of NPC to indirectly examine the correlation between milk consumption and histology-specific NPC incidence according to the GLOBOCAN 2012  which showed the countries/regions with the 20 highest ASIR of NPC in males in the world. Eight countries/regions in our analysis, including Hong Kong (ASIR of NPC: 17.82/105 per person-years), Malaysia (11.91/105), Mainland China (11.65/105), Singapore (10.92/105), Philippines (5.66/105), Algeria (5.08/105), Tunisia (4.28/105) and Thailand (2.94/105) were among the 20 countries/regions with the highest ASIR of NPC (Appendix 3). Our literature review showed that non-keratinizing undifferentiated carcinoma are the predominant subtype in these 8 high-risk countries/regions. The other 40 countries/regions with keratinizing squamous carcinoma as the predominant subtype were classified as low-risk (Appendix 4).
More detailed analysis was conducted for Hong Kong, which has the highest NPC risk and has reliable NPC and risk factors data. Multivariate linear regression analysis was performed between milk excluding butter (which was the only dairy item with available consumption data) consumption per capita and NPC incidence in Hong Kong, with adjustment for other recognized NPC risk factors (salted fish, cigarette smoking, vegetables consumption and socioeconomic status). Repeated analysis using various time lags (range: 0, 5, 10 years) of milk excluding butter consumption per capita prior to the ASIR of NPC was performed to test the robustness of estimates based on different assumptions in the latency period of cancer development subsequent to relevant exposures.
All analyses were performed using SPSS 20.0. p < 0.05 was considered statistically significant.
Ethics approval was not required as data in our study were secondary data from Cancer Incidence in Five Continents series published by the International Agency for Research on Cancer (http://ci5.iarc.fr/) and the Food Supply database maintained by the Food and Agriculture Organization of the United Nations (FAO) (http://faostat.fao.org/). These data were not directly collected from humans or animals.
Correlation between NPC incidence and milk consumption in 48 countries/regions
Spearman rank correlation coefficients (ρ) for the a) cross-sectional and b) time-lagged analysis of the ASIR of NPC by different milk and dairy products are shown in Table 1.
Cross-sectional analysis from 1998 to 2002 (Table 1)
The strongest negative correlation was observed between the ASIR of NPC and fresh cream (CRF) in females (ρ = −0.63; p < 0.001). The second strongest significant, but moderately negative, correlation was observed in milk, excluding butter (MEBC) in males (ρ = −0.55; p < 0.001). In addition, negative and significant correlations were also observed in the categories of milk (whole) (MWF), cheese (CHF), and butter and ghee (BGF), with Spearman coefficients of − 0.35 in males and −0.29 in females,−0.40 in males and −0.42 in females, and −0.36 in males and −0.44 in females, respectively.
Time-lagged analysis in seven lagged models (Table 1)
The strongest negative correlations between CRF and NPC were shown in a 10-year lagged model in males (ρ = −0.71; p < 0.001) and in females (ρ = −0.75; p < 0.001). The significant coefficients increased from the 5-year lag, were maximal at the 10-year lag in males and in females, and declined to the lowest values at the 30-year lag. MEBC, MWF, and BGF had similar patterns for their negative and significant correlations, where the Spearman coefficients increased monotonically. CHF showed a U-shape with significant coefficients in males, which became minimal at the 15-year lag in males and increased and peaked at the longest lag, but the pattern was unclear in females.
Figure 1 shows the negative nonlinear correlation between the ASIR of NPC and per capita consumption of milk excluding butter (MEBC) in 48 countries/regions. The other dairy items showed similar patterns (figures not shown, but are included in Appendix 5). The ASIR of NPC ranged from 0.16/105 person-years (Ecuador) to 17.82/105 (HKSAR) in males (>100-fold difference), and from 0.02/105 (Ecuador) to 6.72/105 (HKSAR) in females (>300-fold difference). Per capita consumption of milk excluding butter varied markedly across the world (>60-fold difference), from about 10 kg/capita/year in China Mainland with lower milk consumption to about 640 kg/capita/year in New Zealand (Appendix 6).
Table 2 shows that the standardized coefficients (sβ) were consistently significant and negative between different milk and dairy product consumption per capita and the ASIR of NPC without adjustment. The strongest negative and significant correlations were observed in MEBC (sβc = −0.499 in males, sβc = −0.500 in females; p < 0.001).
After adjusting for HDI, the negative and significant correlations (adjusted standardized coefficients [asβ]) in MEBC (asβa = −0.552 in males, asβa = −0.534 in females; p < 0.01) and MWF (asβa = −0.344 in males, p < 0.05 and asβa = −0.344 in females; p < 0.05) were found consistently, though the negative correlations between per capita consumption of other dairy products and the ASIR of NPC were attenuated and no longer statistically significant.
Correlation between NPC incidence and milk consumption in high- and low-risk regions
Table 3 shows that the standardized coefficients in high-risk and in low-risk regions with adjustment for HDI. Negative correlations were observed in high-risk regions, but were not statistically significant. In low-risk regions, negative but insignificant coefficients were also observed in different dairy products with adjustment except for CRF with asβl = −0.444 in males; p = 0.04, and −0.375 in females; p = 0.08 which was marginally significant.
Multivariate analysis in Hong Kong
Table 4 shows that negative and significant correlations between per capita consumption of milk excluding butter and ASIR of NPC were observed in Hong Kong at three lagged periods in model 1 (adjusted for per capita consumption of vegetables, salted fish and cigarette smoking) and model 2 (model 1 with additional adjustments, for GDP). The strongest negative correlation in model 1 was observed at 10-year lag in males (the adjusted unstandardized coefficient [auβ] = −0.454; P < 0.05) and in females (auβ = −0.232; P < 0.01), and in model 2 at 10-year lag in males (auβ = −0.439; P < 0.01) and in females (auβ = −0.258; P < 0.01).
The present study showed negative and significant correlations between ASIR of NPC and milk consumption in 48 countries/regions. Countries/regions have different histopathological subtypes of NPC of varied etiology. This may be an important reason underlying the varying patterns of risk factors found across different countries/regions . Our detailed analysis in 8 high-risk countries/regions with predominantly (over 85 %) Type 3 histopathological diagnosis of NPC (non-keratinizing undifferentiated carcinoma) showed negative but non-significant correlations between increased milk consumption and decreased NPC incidence adjusting for HDI. In Hong Kong with the highest NPC risk, we found negative and significant correlations after adjusting for additional potential confounders, including salted fish, cigarette, vegetables consumption and GDP per capita.
Strengths and limitations
Our study had several strengths, including the most comparable and reliable data, a global perspective (48 countries/regions) and a more detailed analysis on Hong Kong data with predominantly (over 95 %) Type 3 of histopathological diagnosis of NPC (non-keratinizing undifferentiated carcinoma) adjusting for potential confounders.
However, it also had several limitations. Although the FAO data are a cost-effective source for ecological comparisons and the best available data for a specific period in most countries/regions, these data did not represent real consumption in individuals (only the milk and dairy products available in the food supply). The major limitation is the ecological fallacy due to unmeasurable individual exposures and confounders. There may be residual confounding. The factors that affect consumption at population level are not the same at individual level. At population level, there are many determinant of various nature and studies based on individual exposure factors are needed. We took account important and known factors by adjusting for the specific time-dependent covariates (socioeconomic status) in the global analysis, with additional adjustment of major confounders, including salted fish, cigarette and vegetables consumption and socioeconomic status for Hong Kong analysis. Race and histology of NPC, are important factors for NPC studies, and should be taken into account in ecological analysis, but race- and histology-specific incidence data are not available in most cancer registries, suggesting that race- and histology-specific incidence data should be collected routinely to enable more in-depth analysis. In addition, temporal sequence was unclear in the cross-sectional analysis to establish causality. Our time-lagged analysis was performed in 5-year to 30-year lagged periods. The results showed that negative and significant coefficients generally increased with longer time lag but with inconsistent patterns for different dairy products, warranting further investigations using longitudinal and time-series analysis for specific dairy products.
We found negative and significant correlations between milk and dairy products per capita consumption and NPC incidence in 48 countries/regions and in Hong Kong. These hypothesis-generating results could support further studies on individual exposures and the disease.
adjusted standardized regression coefficients
adjusted unstandardized regression coefficients
age-standardized incidence rates
butter & ghee
Food and Agriculture Organization of the United Nations
Human Development Index
The International Classification of Diseases 10th Revision
Gross Domestic Product
Hong Kong Special Administrative Region
milk excluding butter
Spearman rank correlation coefficients
standardized regression coefficients
Feng BJ. Descriptive, Environmental and Genetic Epidemiology of Nasopharyngeal Carcinoma. In: Busson P, editor. Nasopharyngeal Carcinoma: Keys for Translational Medicine and Biology. New York: Springer; 2013. p. 23.
Hong Kong Cancer Registry, Hospital Authority. Nasopharyngeal Cancer in 2012. http://www3.ha.org.hk/cancereg/npc_2012.pdf.
Silva IS. Cancer Epidemiology: Principles and Methods. International Agency for Research on Cancer. Lyon: World Health Organization; 1999.
Lee AW, Ng W, Chan Y, Sze H, Chan C, Lam T. The battle against nasopharyngeal cancer. Radiother Oncol. 2012;104(3):272–8.
Jeannel D, Bouvier G, Hubert A. Nasopharyngeal carcinoma: an epidemiological approach to carcinogenesis. Cancer Surv. 1999;33:125–55.
Athyros VG, Tziomalos K, Katsiki N, Gossios TD, Giouleme O, Anagnostis P, et al. The impact of smoking on cardiovascular outcomes and comorbidities in statin-treated patients with coronary artery disease: a post hoc analysis of the GREACE study. Curr Vasc Pharmacol. 2013;11(5):779–84.
Lau HY, Leung CM, Chan YH, Lee AWM, Kwong DLW, Lung ML, et al. Secular trends of salted fish consumption and nasopharyngeal carcinoma: a multi-jurisdiction ecological study in 8 regions from 3 continents. BMC Cancer. 2013;13(1):298.
World Cancer Research Fund / American Institute for Cancer Research. Food, Nutrition, Physical Activity, and the Prevention of Cancer: a Global Perspective. Washington DC: AICR; 2007.
Davoodi H, Esmaeili S, Mortazavian A. Effects of Milk and Milk Products Consumption on Cancer: A Review. Compr Rev Food Sci Food Saf. 2013;12(3):249–64.
Ganmaa D, Li XM, Wang J, Qin LQ, Wang PY, Sato A. Incidence and mortality of testicular and prostatic cancers in relation to world dietary practices. Int J Cancer. 2002;98(2):262–7.
Aune D, Rosenblatt DAN, Chan DS, Vieira AR, Vieira R, Greenwood DC, et al. Dairy products, calcium, and prostate cancer risk: a systematic review and meta-analysis of cohort studies. Am J Clin Nutr. 2015;87–117.
Ahmad OB, Boschi-Pinto C, Lopez AD, Murray CJ, Lozano R, Inoue M. Age standardization of rates: a new WHO standard. Geneva: World Health Organization; 2001.
United Nations, Department of Economic and Social Affairs, Population Division (2013). World Population Prospects: The 2012 Revision. http://esa.un.org/wpp/unpp/panel_population.html. Accessed 17 March 2014.
Nicholls J, Niedobitek G. Histopathological Diagnosis of Nasopharyngeal Carcinoma: Looking beyond the Blue Book. In: Busson P, editor. Nasopharyngeal Carcinoma. New York: Springer; 2013. p. 10–22.
Ferlay J, Soerjomataram I, Ervik M, Dikshit R, Eser S, Mathers C, et al. GLOBOCAN 2012 v1.0, Cancer Incidence and Mortality Worldwide: IARC CancerBase No. 11 [Internet]. Lyon, France: International Agency for Research on Cancer; 2013. Available from: http://globocan.iarc.fr. Accessed 17 March 2014.
Lin JH, Jiang CQ, Ho SY, Zhang WS, Mai ZM, Xu L, et al. Smoking and nasopharyngeal carcinoma mortality: a cohort study of 101,823 adults in Guangzhou, China. BMC Cancer. 2015;15(1):906.
This study supported by Hong Kong RGC Area of Excellence Scheme (Grant Number AoE/M-06/08), and by World Cancer Research Fund UK (WCRF UK) and Wereld Kanker Onderzoek Fonds (WCRF NL), as part of the WCRF International grant programme (grant 2011/460).
The authors have no conflict of interests.
ZMM participated in data collection, conducted data analysis, and drafted the manuscript. CML participated in data collection, and helped to draft the manuscript and revised it critically for important intellectual content. JX, KPC and CMW contributed to the analysis and interpretation of the data, performed the statistical analysis, and revised the manuscript critically for important intellectual content. MLL helped to draft the manuscript and revised it critically for important intellectual content. THL conceived of the study, participated in its design and coordination, and helped to draft the manuscript and revised it critically for important intellectual content. All authors read and approved the final manuscript.
There were four categories of countries/regions according to the different methods of data collection using the CI5, Vol. IX:
For the 16 countries (Australia, Brazil, France, Germany, India, Italy, Japan, China, Malaysia, Poland, Portugal, Spain, Switzerland, Thailand, Turkey, and United Kingdom), each with one or more cancer registries operating, the mean value of the ASIR from different cancer registries was calculated and used as the representative rate for the country. In Italy, for example, 22 ASIRs of NPC in 22 cancer registries were averaged, resulting in one mean value to represent the ASIR of NPC in Italy.
For 12 countries (Algeria, Egypt, Tunisia, Uganda, Zimbabwe, Argentina, Chile, Colombia, Ecuador, Peru, Pakistan, and the Philippines), only one cancer registry was operating in each. The ASIR of the single cancer registry was used to estimate the NPC incidence for the country.
The US has many cancer registries by different states. Among them, the registry for all racial groups in the Surveillance, Epidemiology, and End Results (SEER) Program of the National Cancer Institute (9 registries), which collected data in longer periods, was used as representative for the US.
For the remaining 19 countries/regions (Costa Rica, Canada, Cyprus, Israel, Korea, Kuwait, Austria, Bulgaria, Denmark, Finland, Hong Kong Special Administrative Region (HKSAR), Iceland, Ireland, Malta, The Netherlands, Norway, Singapore, Sweden and New Zealand), the cancer registries collected data for the whole country/region. The incidence rates for these countries/regions were used as provided.
In some countries (Zimbabwe, Israel, and Kuwait), analysis by ethnic groups was supported by some cancer registries. In our analysis, Zimbabwe included the African race as the representative population, Israel included all ethnic groups (both Jews and non-Jews), and Kuwait included Kuwaitis.
In a separate analysis, annual NPC ASIRs in Hong Kong from 1983 to 2009 were obtained from the Hong Kong Cancer Registry with the 1966 world standard population as reference.