Cluster of liver cancer and immigration: A geographic analysis of incidence data for Ontario 1998–2002
© Chen et al; licensee BioMed Central Ltd. 2008
Received: 26 March 2008
Accepted: 02 June 2008
Published: 02 June 2008
Liver cancer is not common in Canada in general; however, clustering of the disease causes a concern. We conducted a spatial analysis to determine the geographic variation of liver cancer and its association with the proportion of immigration in Ontario. Liver cancer incidence data between 1998 and 2002 were obtained from the Ontario Cancer Registry. The Canadian Community Health Survey (CCHS) in 2001 provided information on potential risk factors.
Age standardized incidence ratios (SIR) for liver cancer and prevalence of potential risk factors were calculated for each of 35 health regions. The SIRs for liver cancer varied across the 35 health regions (p < 0.01). Toronto and York health regions had a significantly higher SIR than other regions, indicated by the Scan method (p < 0.001). Poisson models with and without random effects were fitted to determine independent ecological contributors. After adjustment for sex, age and spatial location, the proportion of immigrants remained a significant determinant. Smoking, alcohol drinking, physical activity, education, income, obesity and diabetes did not substantially explain the geographic variation of liver cancer in Ontario.
Immigration is an important reason for the clustering of liver cancer in Ontario. More attention should be paid to areas with a high proportion of immigrants.
Liver cancer is the sixth most common cancer worldwide and the third most common cause of death from cancer. Liver cancer incidence has been increasing in developed countries [2–6] including Canada, but 82% of cases (and deaths) are in developing countries. Hepatitis B and C viruses are major risk factors for liver cancer, and more than 75% of cases worldwide, and 85% of cases in developing countries, are caused by these two viruses[1, 8, 9]. Liver cancer affects Asian/Pacific Islander Americans disproportionately in the U.S.[10, 11]. In Canada, the age-standardized mortality rate for liver cancer is most significantly elevated for birthplace of Asia and to a lesser degree for Europe. However, a study from Houston, Texas, suggests that the increase in liver cancer risk is not from immigration but represents a true rise among the native born population.
Ontario, located in the centre part of Canada, has a population of 12 million or 38% of the total Canadian population. The sex and age standardized incidence rate (/100,000) of liver cancer was 3.18 in 2004 compared to 2.18 in. In this analysis, we used the data collected by the Ontario Cancer Registry in Ontario for the period from 1998 to 2002 and the 2001 Canadian Community Health Survey (CCHS). We examined the geographic distribution of liver cancer incidence and assessed the impacts of area level measures of immigrant, smoking, alcohol use, obesity, and diabetes in Ontario, Canada.
Average and range of liver cancer incidence and population characteristics of 35 health regions in Ontario, Canada, 2000–2001
Minimum – maximum
Coefficient of variance
Incidence of liver cancer (/100,000)
Current smoking (%)
Alcohol drinking (%)
Physically active (%)
Reduced fruit and vegetable intake (%)
High education (%)
Low income (%)
Overall spatial correlation for liver cancer across 35 health regions in Ontario – Hierarchical Poisson models
Variance (SE) for RREa
Moran's Ib (p value)
Geary's Cb (p value)
+ selected factorsc
Regression coefficients (β) and standard errors (SEβ) for liver cancer in Ontario – Poisson regression analysis
Sex × age interaction
Proportion of immigrants
Our analysis demonstrated that the incidence of liver cancer in Ontario was not randomly distributed geographically. There was a significant variation among health regions and one high risk cluster included Toronto, York and Peel health regions. This geographic variation in liver cancer incidence did not result from differences in age and sex. Although elderly people and men had increased risk of liver cancer, the observed geographic variation remained when results were adjusted for age and sex.
We found that the proportion of immigrants was a main determinant for the geographic variation of liver cancer in Ontario. After the proportion of immigrants was taken into consideration, both Moran's I and Geary's C were no longer significant. Rosenblatt et al. compared the incidence of primary liver cancer in Chinese, Japanese, and Filipino migrants to the U.S. with that of their descendants and U.S. born Whites. The incidence was higher for men born in Asia than Asian men born in the U.S., and U.S. Whites had the lowest rate. Luo et al. examined the incidence of cancer among Chinese immigrants in Alberta, and found that the overall cancer incidence was lower among immigrants, but the incidence of liver cancer was much higher compared to Canadian-born Alberta residents. Globally, sub-Saharan Africa, eastern and southeastern Asia, and Melanesia have high rates of liver cancer.
Although it is known that hepatitis B and C viruses are major risk factors for liver cancer[1, 8, 9], the cause of rising incidence of liver cancer in both the U.S. and Canada remains unclear. Hepatitis C virus (HCV) infection can explain a substantial proportion of liver cancer incidence in the U.S. in recent years[16, 17] and might also explain the differences between Danish and U.S. trends in liver cancer incidence. In addition to drug use and blood transfusions, immigration can also be an major source of HCV. Recent immigrants from hepatitis B virus (HBV) endemic areas and their descendants are at high risk of chronic HBV infection, and HBV-related liver cancer, and it is expected to be similar for HCV.
Canadian census data show an in increase in immigration from Africa and Asia in recent decades. Immigrants from Africa and Asia made up about 17% of the foreign-born population in 1981, increasing to 28% in 1996 and 42% in 2001. Concurrently, immigrants from Europe made up a decreasing proportion of the foreign-born population, beginning at 67% in 1981 and dropping to 55% in 1996 and then 42% in 2001. There is a possibility that the increased incidence rate of liver cancer is related to the increased proportion of immigrants from Asia and Africa during the past several decades.
It has been documented previously that increased incidence of liver cancer is associated with smoking[20, 21], heavy alcohol drinking[22, 23] and obesity In this analysis, we did not find that smoking, alcohol drinking, obesity, diabetes, fruit and vegetable consumption, physical activity, income and education at the health region level explained a significant part of the geographic variation of liver cancer incidence in Ontario.
The study has potential limitations. Within health regions, there are clusters for people with different ethnic backgrounds. In the city of Toronto for example, there is an increased proportion of immigrants from Asia who live in the Scarborough area. Smaller area units may result in finer results. This is related to the scale effects of the modifiable areal unit problem[25, 26]. However, in this study there were no smaller regional data available. Although the association between immigration and liver cancer is important, there is no conclusion of causation from this study. We did not have information on important risk factors such as HBV or HCV infection. Other potential risk factors data were collected in 2001 while the cancer incidence was for the period from 1998 to 2002. We can only assume that the current exposure correlates the previous exposure, which may or may not be true. Although immigration status did not change over time, the proportion of immigrants in health regions might change. In addition, information on race and ethnicity is not available in the Canadian Cancer Registry data, and therefore, we do not know if these immigrants are more likely to come from areas where HBV and/or HCV are endemic. In this analysis, immigrant proportion explained most of the variation of region random effects, which decreased from 0.080 to 0.049. However, the residual variation remained significant even after other covariates were included in the model, indicating that there are some unidentified important factors for the regional difference. Heavy alcohol drinking is an important risk factor for liver cancer; however, the quantity and frequency of alcohol consumption could not be easily adjusted for at the ecological level.
The significant geographic variations in the incidence of liver cancer in Ontario could largely be explained by differences in the proportion of immigrants. The prevalence of smoking, alcohol use, obesity and diabetes as well as physical activity, fruit and vegetable consumption, education and income were not significant contributors to the geographic variation. Further studies are needed to examine the reasons for the association between immigration and liver cancer, such as the prevalence of HBV or HCV infections among immigrants and non-immigrants. The study provides useful information that may help planning and distribution of resources. The results may also underscore the importance of efforts to promote screening and treatment programs for liver cancer in the regions with a high risk.
Cancer incidence data were obtained from the Ontario Cancer Registry (OCR). The OCR has accumulated data for liver cancer since the 1960s. In this analysis we used data collected during the period from 1998 to 2002. Cancers were coded according to the International Classification of Diseases Second Edition (ICDO-2). Eligible cases for inclusion in this analysis were all individuals from 35 health regions in the province with primary liver cancer (ICDO-3 C22.0 excluding 9590:9989). We obtained number of cases and population for each region by gender and age groups from the OCR. We excluded cases under age of 15 years. Hepatocellular carcinoma is the most common type of primary liver cancer in Canada, accounting for about 85% of all cases. Prevalence of potential cancer risk factors were calculated based on data from the CCHS conducted from September 2000 to November 2001 (cycle 1.1). The target population of the survey was household residents aged 12 years or older in all the ten provinces and three territories in Canada, excluding individuals living on Indian Reserves or Crown lands, clientele of institutions, full-time members of the Canadian Armed Forces and residents of certain remote regions. The survey targeted approximately 98% of Canadians aged 12 years or more and included more than 130,000 subjects. The questionnaire was administered using computer-assisted interviewing. Sampling units selected from a telephone list frame were interviewed from centralized centres, and those selected from an area frame were mainly interviewed by decentralized field interviewers. In all selected dwellings, a knowledgeable household member was asked to supply basic demographic information on all residents of the dwelling. One member of the household was then selected for a more in-depth interview. The survey included questions related to health status, health care utilization and health determinants. To match with the cancer data, only data from people aged 15+ years living in Ontario were used in the analysis. In the public use micro data file, some regions were combined with a total of 35 health regions in Ontario. Based on the health region map, an adjacency matrix was created manually for 35 collapsed regions to reflect the neighborhood structure and used in the spatial analysis. Geographic centroids for each health region were calculated based on the map polygon.
The CCHS covered information on potential cancer risk factors. For our study, we obtained the proportion of individuals in each health region in Ontario with daily smoking, regular alcohol use, leisure time physical activity, obesity and diabetes. Current smokers were respondents who had smoked at least 100 cigarettes during their lifetime and reported smoking cigarettes every day or almost every day at the time of the survey. A regular alcohol drinker was defined as a person who drank alcoholic beverages at least once a month during the past twelve month. Active subjects were those that had total daily energy expenditure values of 3+ kcal/kg/day for leisure time physical activity. Obesity was defined as body mass index equal to or larger than 30 kg/m2. People were considered to have diabetes if they answered the following question with an affirmative response: "Do you have diabetes diagnosed by a health professional?" Other variables included in the analysis were income, education and stress.
The standardized incidence ratio (SIR) for live cancer was calculated for each of the 35 health regions and mapped to visually show the geographic distribution. To detect the specific locations of either high rate or low rate clusters with a minimum of assumptions about cluster size, and to evaluate their statistical significance, we employed the spatial scan statistic proposed by Kulldorff[27, 28]. This method imposes an infinite number of circles on the map, with circle centroids at any of the region coordinates and with different radius. We tested the difference in cancer incidence between inside and outside the circle. Each circle is a potential cluster that is least likely to have occurred by chance. Moran's I and Geary's C were calculated for the distribution of incidence ratio of liver cancer to assess the overall spatial correlation between neighbouring regions[29, 30] Moran's I ranges from -1 to 1, with a positive/negative signage representing positive/negative spatial autocorrelation and zero representing no spatial autocorrelation. Geary's C ranges from 0 to +2, with zero being a strong positive spatial autocorrelation, through to 2 which represents a strong negative spatial autocorrelation.
For both Moran's I (or Geary's C) and the scan method describe above, p-values were obtained by using Monte Carlo hypothesis testing with an adjustment for multiple testing. Hierarchical Poisson models were fitted to explore the association between aggregated factors and liver cancer incidence with health region treated as a clustering factor. The aggregated factors included alcohol use, smoking, fruits and vegetable consumption, obesity, immigrant status, income and education as well as geographic location of the health region (longitude/latitude was used as a proxy for the location of health regions). Both random effects (health regions) and fixed effects (alcohol use, smoking, exercise, diabetes, obesity, fruit and vegetable consumption, education and income) were analyzed simultaneously. A conditional autocorrelation model (CAR) was constructed to examine the impacts of those factors on liver cancer with both spatial correlation and regional random effects being considered. All the analyses were conducted by using SAS (SAS Institute Inc., Cary, NC), SaTScan package http://www.satscan.org/ and WINBUGS software http://www.mrc-bsu.cam.ac.uk/bugs. Maps were created using ArcGIS package http://www.esri.com/software/arcgis/.
- CAR model:
conditional autocorrelation model
Canadian Community Health Survey
Hepatitis B virus
hepatitis C virus
International Classification of Diseases Second Edition
leisure time physical activity
Ontario Cancer Registry
standardized incidence ratio
The work is partially supported by fund from the Public Health Agency of Canada. The cancer data were provided by the Ontario Cancer Registry. The authors thank Robert Semenciw for his assistance in conducting the project and his helpful comments on the manuscript.
- Parkin DM, Bray F, Ferlay J, Pisani P: Global cancer statistics, 2002. CA Cancer J Clin. 2005, 55 (2): 74-108.PubMedView ArticleGoogle Scholar
- El-Serag HB, Davila JA, Petersen NJ, McGlynn KA: The continuing increase in the incidence of hepatocellular carcinoma in the United States: an update. Ann Intern Med. 2003, 139 (10): 817-823.PubMedView ArticleGoogle Scholar
- El-Serag HB, Mason AC: Rising incidence of hepatocellular carcinoma in the United States. N Engl J Med. 1999, 340 (10): 745-750. 10.1056/NEJM199903113401001.PubMedView ArticleGoogle Scholar
- McGlynn KA, Tsao L, Hsing AW, Devesa SS, Fraumeni JF: International trends and patterns of primary liver cancer. Int J Cancer. 2001, 94 (2): 290-296. 10.1002/ijc.1456.PubMedView ArticleGoogle Scholar
- Capocaccia R, Sant M, Berrino F, Simonetti A, Santi V, Trevisani F: Hepatocellular carcinoma: trends of incidence and survival in Europe and the United States at the end of the 20th century. Am J Gastroenterol. 2007, 102 (8): 1661-70; quiz 1660, 1671. 10.1111/j.1572-0241.2007.01337.x.PubMedView ArticleGoogle Scholar
- McGlynn KA, Tarone RE, El-Serag HB: A comparison of trends in the incidence of hepatocellular carcinoma and intrahepatic cholangiocarcinoma in the United States. Cancer Epidemiol Biomarkers Prev. 2006, 15 (6): 1198-1203. 10.1158/1055-9965.EPI-05-0811.PubMedView ArticleGoogle Scholar
- Dyer Z, Peltekian K, van Zanten SV: Review article: the changing epidemiology of hepatocellular carcinoma in Canada. Aliment Pharmacol Ther. 2005, 22 (1): 17-22. 10.1111/j.1365-2036.2005.02504.x.PubMedView ArticleGoogle Scholar
- Donato F, Boffetta P, Puoti M: A meta-analysis of epidemiological studies on the combined effect of hepatitis B and C virus infections in causing hepatocellular carcinoma. Int J Cancer. 1998, 75 (3): 347-354. 10.1002/(SICI)1097-0215(19980130)75:3<347::AID-IJC4>3.0.CO;2-2.PubMedView ArticleGoogle Scholar
- Tanaka Y, Kurbanov F, Mano S, Orito E, Vargas V, Esteban JI, Yuen MF, Lai CL, Kramvis A, Kew MC, Smuts HE, Netesov SV, Alter HJ, Mizokami M: Molecular tracing of the global hepatitis C virus epidemic predicts regional patterns of hepatocellular carcinoma mortality. Gastroenterology. 2006, 130 (3): 703-714. 10.1053/j.gastro.2006.01.032.PubMedView ArticleGoogle Scholar
- Chang ET, Keegan TH, Gomez SL, Le GM, Clarke CA, So SK, Glaser SL: The burden of liver cancer in Asians and Pacific Islanders in the Greater San Francisco Bay Area, 1990 through 2004. Cancer. 2007, 109 (10): 2100-2108. 10.1002/cncr.22642.PubMedPubMed CentralView ArticleGoogle Scholar
- Choe JH, Chan N, Do HH, Woodall E, Lim E, Taylor VM: Hepatitis B and liver cancer beliefs among Korean immigrants in Western Washington. Cancer. 2005, 104 (12 Suppl): 2955-2958. 10.1002/cncr.21518.PubMedPubMed CentralView ArticleGoogle Scholar
- ElSaadany S, Giulivi A: Epidemiology of hepatocellular carcinoma in Canada, 1995: analysis of death certificates. Chronic Dis Can. 2006, 27 (3): 125-129.PubMedGoogle Scholar
- Kulkarni K, Barcak E, El-Serag H, Goodgame R: The impact of immigration on the increasing incidence of hepatocellular carcinoma in the United States. Aliment Pharmacol Ther. 2004, 20 (4): 445-450. 10.1111/j.1365-2036.2004.02090.x.PubMedView ArticleGoogle Scholar
- Rosenblatt KA, Weiss NS, Schwartz SM: Liver cancer in Asian migrants to the United States and their descendants. Cancer Causes Control. 1996, 7 (3): 345-350. 10.1007/BF00052940.PubMedView ArticleGoogle Scholar
- Luo W, Birkett NJ, Ugnat AM, Mao Y: Cancer Incidence Patterns Among Chinese Immigrant Populations in Alberta. Journal of Immigrant Health. 2004, 6 (1): 41-48. 10.1023/B:JOIH.0000014641.68476.2d.PubMedView ArticleGoogle Scholar
- El-Serag HB, Mason AC: Risk factors for the rising rates of primary liver cancer in the United States. Arch Intern Med. 2000, 160 (21): 3227-3230. 10.1001/archinte.160.21.3227.PubMedView ArticleGoogle Scholar
- Davila JA, Morgan RO, Shaib Y, McGlynn KA, El-Serag HB: Hepatitis C infection and the increasing incidence of hepatocellular carcinoma: a population-based study. Gastroenterology. 2004, 127 (5): 1372-1380. 10.1053/j.gastro.2004.07.020.PubMedView ArticleGoogle Scholar
- Jepsen P, Vilstrup H, Tarone RE, Friis S, Sorensen HT: Incidence rates of hepatocellular carcinoma in the U.S. and Denmark: recent trends. Int J Cancer. 2007, 121 (7): 1624-1626. 10.1002/ijc.22860.PubMedView ArticleGoogle Scholar
- McQuillan GM, Townsend TR, Fields HA, Carroll M, Leahy M, Polk BF: Seroepidemiology of hepatitis B virus infection in the United States. 1976 to 1980. Am J Med. 1989, 87 (3A): 5S-10S. 10.1016/0002-9343(89)90523-8.PubMedView ArticleGoogle Scholar
- Zhu K, Moriarty C, Caplan LS, Levine RS: Cigarette smoking and primary liver cancer: a population-based case-control study in US men. Cancer Causes Control. 2007, 18 (3): 315-321. 10.1007/s10552-006-0105-8.PubMedView ArticleGoogle Scholar
- Mizoue T, Tokui N, Nishisaka K, Nishisaka S, Ogimoto I, Ikeda M, Yoshimura T: Prospective study on the relation of cigarette smoking with cancer of the liver and stomach in an endemic region. Int J Epidemiol. 2000, 29 (2): 232-237. 10.1093/ije/29.2.232.PubMedView ArticleGoogle Scholar
- Makimoto K, Higuchi S: Alcohol consumption as a major risk factor for the rise in liver cancer mortality rates in Japanese men. Int J Epidemiol. 1999, 28 (1): 30-34. 10.1093/ije/28.1.30.PubMedView ArticleGoogle Scholar
- Ioannou GN, Weiss NS, Kowdley KV: Relationship between transferrin-iron saturation, alcohol consumption, and the incidence of cirrhosis and liver cancer. Clin Gastroenterol Hepatol. 2007, 5 (5): 624-629. 10.1016/j.cgh.2007.01.008.PubMedView ArticleGoogle Scholar
- Qian Y, Fan JG: Obesity, fatty liver and liver cancer. Hepatobiliary Pancreat Dis Int. 2005, 4 (2): 173-177.PubMedGoogle Scholar
- Willis A, Krewski D, Jerrett M, Goldberg MS, Burnett RT: Selection of ecologic covariates in the American Cancer Society study. J Toxicol Environ Health A. 2003, 66 (16-19): 1563-1589. 10.1080/15287390306425.PubMedView ArticleGoogle Scholar
- Klliott P: Spatial Epidemiology. Geographical analysis of communicable disease data. Edited by: Atkinson PJMA. 2000, Oxford , Oxford Unviersity Press, 253-266.Google Scholar
- Kulldorff M: A spatial scan statistics. Communications in Statistics: Theory and Methods. 1997, 26: 1487-1496.Google Scholar
- Kulldorff M: Spatial scan statistics: models, calculations, and applications. Scan statistics and applications. Edited by: Glaz J, Balakrishnan N. 1999, Boston , Birkhèauser, 303-319.View ArticleGoogle Scholar
- Moran PA: Notes on continuous stochastic phenomena. Biometrika. 1950, 37 (1-2): 17-23. 10.1093/biomet/37.1-2.17.PubMedView ArticleGoogle Scholar
- Geary RC: The contiguity ratio and statistical mapping. The Incorporated Statistician. 1954, 5: 115-145. 10.2307/2986645.View ArticleGoogle Scholar
- Banerjee S, Carlin BP, Gelfand AE: Hierarchical modeling and analysis for spatial data. Monographs on statistics and applied probability ; 101. 2004, Boca Raton, Fla. , Chapman & Hall/CRC, xvii, 452 p.-Google Scholar
- Lawson A, Browne WJ, Vidal Rodeiro CL: Disease mapping with WinBUGS and MLwiN. 2003, Hoboken, NJ , J. Wiley, xiii, 277 p.-View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.