Association of diabetes and cancer mortality in American Indians: the Strong Heart Study

Purpose The metabolic abnormalities that accompany diabetes mellitus are associated with an increased risk of many cancers. These associations, however, have not been well studied in American Indian populations, which experience a high prevalence of diabetes. The Strong Heart Study is a population-based, prospective cohort study with extensive characterization of diabetes status. Methods Among a total cohort of 4,419 participants who were followed for up to 20 years, 430 cancer deaths were identified. Results After adjusting for sex, age, education, smoking status, drinking status, and body mass index, participants with diabetes at baseline showed an increased risk of gastric (HR 4.09; 95 % CI 1.42–11.79), hepatocellular (HR 2.94; 95 % CI 1.17–7.40), and prostate cancer mortality (HR 3.10; 95 % CI 1.22–7.94). Further adjustment for arsenic exposure showed a significantly increased risk of all-cause cancer mortality with diabetes (HR 1.27; 95 % CI 1.03–1.58). Insulin resistance among participants without diabetes at baseline was associated with hepatocellular cancer mortality (HR 4.70; 95 % CI 1.55–14.26). Conclusions Diabetes mellitus, and/or insulin resistance among those without diabetes, is a risk factor for gastric, hepatocellular, and prostate cancer in these American Indian communities, although relatively small sample size suggests cautious interpretation. Additional research is needed to evaluate the role of diabetes and obesity on cancer incidence in American Indian communities as well as the importance of diabetes prevention and control in reducing the burden of cancer incidence and mortality in the study population.


Introduction
A growing body of evidence links diabetes and chronic hyperglycemia with the development of several common cancers, including pancreas, liver, colon, and breast & Lyle G. Best lbest@restel.com cancers [1][2][3][4][5]. In a meta-analysis of 12 studies from Europe, Asia, and the USA, diabetes was also associated with increased all-cause cancer mortality (pooled relative risk 1.16, 95 % CI 1.03-1.30) [6]. Biologic plausibility for a direct causal association between cancer and hyperglycemia was first suggested in the early 1900s by Warburg et al. based on the increased reliance of cancer cells on glucose and aerobic glycolysis (the ''Warburg Effect'') [7]. This effect has recently been further corroborated by the detection of increased glucose metabolism in malignant tissue [8] and the demonstration of these metabolic changes in individual cancer cells [9]. The increased availability of glucose may enhance the survival of malignant cells in their early stages or provide a growth advantage during later stages. Insulin, and insulin-like growth factors (IGFs) which are increased during much of the course of type II diabetes [10,11], could also increase cancer risk.
Cancer has become a leading cause of death among American Indians [12,13], who have the highest rates of type 2 diabetes in the USA [14][15][16]. Little is known, however, about the associations between hyperglycemic conditions and cancer among American Indians. In 2003, the Strong Heart Study, a population-based prospective cohort study in American Indian communities of Arizona, Oklahoma, and North and South Dakota, failed to identify an association between diabetes and all-cause cancer mortality over 9 years of follow-up [17]. However, that initial study was limited by a small number of cases of cancer mortality as well as a relatively short period of follow-up. In addition, specific cancers were not evaluated in that report [17].
The Strong Heart Study has continued to follow older American Indians since 1989, contributing important information on the relationship of diabetes and hyperglycemic conditions to complex diseases [18][19][20][21]. Although the study was primarily funded to investigate cardiovascular disease and its risk factors, cancer deaths have also been identified as part of the mortality surveillance system. The objective of the current study is to determine the association of diabetes, hyperglycemia, and insulin resistance with cancer mortality among American Indians participating in the Strong Heart Study. Based on the evidence from other populations, we hypothesized that diabetes would be associated with an increase in total cancer mortality, as well as with the following cancer-specific mortalities: gastric, colorectal, breast, renal cell, liver, and pancreas cancer.

Study population
From 1989 to 1991, men and women 45-75 years of age from 13 American Indian communities were invited to participate in the Strong Heart Study. In Arizona and Oklahoma, every eligible person was invited, whereas in North/South Dakota, a cluster sampling technique was used. The baseline participation rate was 62 %, with a final sample of 4,549 participants. All analysis was based on baseline data collected in the first examination of this cohort in 1989-1991. For this study, we excluded participants missing information on diabetes status (n = 90), body mass index (BMI) (n = 20), smoking status (n = 8), alcohol consumption (n = 8), and education level (n = 4), leaving 4,419 participants for analysis. The SHS protocol was approved by multiple institutional review boards and by the participating communities. All participants provided oral and written informed consent.

Data collection
Study visits were performed by trained and certified examiners following a standard protocol and included a questionnaire (sociodemographic factors, environmental exposures, and medical history), a physical examination (height, weight, and blood pressure), and bio-specimen collection (blood, urine). Participants were asked to fast for 12 h before blood collection. Plasma glucose was measured by a hexokinase method at MedStar Health Research Institute (MHRI), Hyattsville, MD. A 2-h, 75-g oral glucose tolerance test was performed in participants without diabetes medication. The oral glucose tolerance test was not performed if the participant received renal dialysis (or a renal transplant) or if the fasting glucose level was greater than 225 mg/dL as determined by Acucek II (Baxter Healthcare, Grand Prairie, Texas). HbA1c was measured by a high-performance liquid chromatography method at the laboratory of the National Institute of Diabetes and Digestive and Kidney Diseases Epidemiology and Clinical Research Branch, Phoenix, Arizona. Plasma insulin was measured at MHRI by radioimmunoassay (Linco, St. Louis, Missouri).

Diabetes definition
Diabetes was defined as a fasting glucose level C126 mg/ dL or the use of insulin or oral hypoglycemic medication. Those on renal dialysis or with a kidney transplant and that responded positively to the question ''Has a medical person ever told you that you had diabetes?'' were also classified as having diabetes. The date of diagnosis was either the date of the baseline examination resulting in diagnostic criteria or the self-reported date of diagnosis from the baseline questionnaire. Because type 1 diabetes is rare, we assumed essentially all cases in our study have type 2 diabetes. Impaired fasting glucose (IFG) was defined as fasting glucose C100 and \126 mg/dL in the absence of diabetes as defined above. Insulin resistance was estimated in non-diabetic participants using the homeostasis model assessment to quantify insulin resistance (HOMA-IR) = [fasting plasma insulin (mU/L) 9 fasting plasma glucose (mmol/L)]/22.5. In addition to HOMA-IR, we also evaluated plasma insulin.

Cancer mortality follow-up
Follow-up for mortality was complete for 99.8 % of the study population. Death certificates were obtained from the State Departments of Health and the underlying cause of death used in analysis. If the death certificate indicated that an autopsy had been performed, the medical examiner's report was obtained. Death certificate codes were recorded according to the International Classification of Diseases, 9th Revision (ICD-9). In addition to overall cancer, we evaluated the following specific cancers: esophagus and stomach (ICD-9 150-151), colon and rectum (153-154), liver and intrahepatic bile ducts (157) (referred to from now on as liver cancer), gallbladder and extrahepatic bile ducts (156), bronchus and lung (162.2-162.9), breast (174), prostate (185), kidney (189.0), and lymphatic and hematopoietic tissue (200-208). The follow-up interval was calculated from the date of baseline examination to the date of death or 31 December 2008, whichever occurred first. The median follow-up time among participants who did not develop cancer was 17.2 years.
We conducted and have presented the analyses for all cancers, except for those which had too few deaths to properly analyze (brain, skin, and reproductive organs). Analyses for HbA1c and HOMA-IR are only shown for cancers that we found a statistically significant association with diabetes or a suggestive trend for a cancer previously associated with diabetes such as colorectal cancer. For the other cancers, we confirmed that none of them was associated with HbA1c or HOMA-IR (data not shown).

Statistical methods
Statistical analyses were conducted with Stata version 12.1 (StataCorp LP, College Station, Texas). The figures were conducted with R version 3.0.1 (R project for statistical computing). The prospective association between baseline diabetes status (as a dichotomous trait, as well as a categorical trait (i.e., no IFG, IFG, or diabetes)) and cancer mortality (overall and site specific) was assessed using Cox proportional hazards models with age as timescale and individual starting follow-up times treated as staggered entries. The nonparametric underlying baseline hazards were allowed to differ by study region using the strata option. By adding study region as strata, the model adjusts for region, but no association is estimated for region and region does not need to meet the assumption of proportionality of hazards. Initial Cox proportional hazards models accounted for age, sex, and region (model 1). Model 2 also adjusted for BMI, and in the case of breast cancer for baseline menopausal status (pre/post), hormone replacement therapy (current/past/never use), and parity (0/ 1-2/3-4/C5). Model 3 further adjusted for education, smoking (never, former, current), and drinking status (never drinker, former drinker, light current drinker (\4 drinks/week), moderate current drinkers (4-12 drinks/ week), heavy drinkers ([12 drinks/week)) [22]. The assumption of hazards proportionality was evaluated visually based on the smoothed association between age and scaled Schoenfeld residuals, with no major departures from proportionality.
To evaluate the relationship between HbA1c and cancer mortality, we modeled HbA1c as a continuous variable and derived adjusted hazard ratios of cancer mortality comparing the 80th versus the 20th percentiles of HbA1c. Additionally, HbA1c levels were separately modeled using restricted cubic splines with knots at the 10th (4.6 %), 50th (5.6 %), and 90th (10.7 %) percentile of its distribution.
HOMA-IR values were log-transformed, and their association with cancer mortality was evaluated using Cox proportional hazards models, using the same levels of adjustment as in models 1-3 but restricted to participants without diabetes. Plasma insulin was also log-transformed and analyzed similarly in a separate model. Subgroup analyses (e.g., by sex, age groups, or smoking status) were not performed due to limited sample size, especially for specific cancer types. We performed several sensitivity analyses. First, we further adjusted for cigarette pack-years as well as for arsenic in the subset of participants with information on pack-years and urine arsenic concentrations available (n = 3,737). Arsenic is an established carcinogen that has been related to cancer mortality [23] and diabetes [24] in the Strong Heart Study. Second, we further adjusted for time since the diagnosis of diabetes. Third, to account for competing risks by causes of death other than cancer, we estimated proportional hazards regression models according to the method of Fine and Gray [25]. Findings from all sensitivity analyses were consistent with those reported (data not shown).

Results
Tables 1 and 2 show participant characteristics stratified by cancer mortality and diabetes status, respectively. In this cohort of American Indians, there were 13.7 % participants with IFG and 45.7 % with diabetes at baseline. The prevalence of diabetes in this cohort was previously reported as 43 % for men and 52 % for women; a random sample (n = 311) of non-participants from this population showed a similar proportion of participants and non-participants with diabetes (40 and 38 %, respectively) [14]. During follow-up, 187 men and 243 women died from cancer, mainly from lung and prostate cancer in men and lung and breast cancer in women.
While confidence intervals are wide, due to small sample size, after multivariate adjustment for age, sex, BMI, education, drinking status, and smoking status (Table 3), the hazard ratio (95 % CI) for overall cancer mortality was 1.19 (0.97-1.45) for those with diabetes compared to those without. The corresponding hazard ratios (95 % CI) for gastric, liver, and prostate cancers were 4.09 (1.42-11.79), 2.94 (1.17-7.40), and 3.11 (1.22-7.94), respectively. Similar results were observed in models that further adjusted for log-transformed urine arsenic and for number of packyears (n = 3,737), although the effect for overall cancer became significant in that analysis (HR 1.27; 95 % CI 1.03-1.58), and a suggestive, but not statistically significant, association was observed for colorectal cancer (HR 1.92; 95 % CI 0.92-4.00) (data not shown in tables). In the analysis of breast cancer, additional adjustment for parity and menopausal status did not affect our results. Similarly, adjustment for physical activity did not reveal significant changes in the risk of specific cancers, compared with model 3.
When modeling the dose-response relationship between HbA1C and cancer mortality, none of the associations were statistically significant (Fig. 1). For gastric cancer, the nonsignificant association was positive and fairly linear. For prostate cancer, the nonsignificant association was nonlinear. Total cancer was not associated with HbA1C (HR per IQR in HbA1c levels: 0.99; 95 % CI 0.85, 1.15).
Results for HOMA-IR and cancer (Table 4) showed only a positive association between HOMA-IR values and liver cancer mortality (HR 4.70 per log unit; 95 % CI 1. 55-14.26). Unexpectedly, we also observed a negative association between HOMA-IR values and prostate cancer mortality. Similar findings were observed with plasma insulin (Table 4). HOMA-IR was not associated with total cancer mortality (HR per log unit of HOMA-IR was 0.92; 95 % CI 0.73, 1.18).

Discussion
In this population-based cohort study with almost 20 years of follow-up, we found that baseline type 2 diabetes was associated with a borderline significant increase risk of overall cancer mortality among American Indians from Arizona, Oklahoma, and North/South Dakota. The association became significant in a subset of the population with information available on arsenic exposure, after adjustment for this variable. Hyperglycemia, as measured by HbA1c, was also nonsignificantly associated with increased cancer risk. We found no clear association between baseline insulin resistance and overall cancer mortality, although HOMA-IR levels were positively associated with liver cancer and inversely associated with prostate cancer mortality.
The present study provides further evidence for increased risk of specific cancer mortality, including gastric and liver cancer, with baseline diabetes, although cautious interpretation is warranted, given the relatively small sample size. These findings are consistent with our a priori hypothesis based on reports conducted in Australia, East  Asia, North America, and Europe [26][27][28][29][30]. Also consistent with our a priori hypothesis, we found a nonsignificant, but increased mortality risk estimate for colorectal cancer. Hyperglycemia may contribute to the initial mutagenic events initiating neoplasia, via an increase in reactive oxygen species [31]. Additionally, two subsequent independent mechanisms enhancing carcinogenesis could promote cancer development: (1) the increased availability of glucose coupled with the Warburg effect and (2) the availability of endogenous insulin and insulin-like growth factors (IGFs), which are increased before and during the early course of type II diabetes and have a synergistic effect on tumor growth. Interestingly, Laron dwarfism, a syndrome characterized by a lack of functional IGF-1 receptors, confers a marked reduction in cancer risk [32]. In epidemiologic studies, treatment with insulin [33,34] or insulin secretagogues [35] has been associated with an increased risk of some cancers, including breast, colon, pancreas, and prostate. Also, agents that decrease insulin resistance and thus decrease insulin levels (metformin) have been associated with a decreased risk of stomach, liver, pancreatic, and colorectal cancers [36], although there is recent conflicting evidence [37]. It is unclear whether hyperglycemia could show effects for certain, specific cancers and not for others.
In our study, HbA1c, HOMA-IR, and plasma insulin per se were associated with an increased risk of liver cancer mortality among those without diabetes, a finding that is consistent with previous epidemiologic studies [38][39][40].
An extensive review of meta-analyses points out various weaknesses in previous studies associating specific cancers with diabetes, but finds the most robust evidence for breast, intrahepatic cholangiocarcinoma, colorectal, and endometrial cancers [41]. This is consistent with our findings, while allowing for potential overlap in the classification of hepatocellular versus cholangiocarcinoma. Interestingly, the Beckwith-Wiedemann syndrome, with known elevations of insulin and IGFs during the first years of life, is associated with very high incidence of cancer in childhood, especially hepatoblastoma and Wilm's tumor [42]. Another Fig. 1 Distribution of glycated hemoglobin (% HbA1c) and hazard ratios for selected cancer mortality by % HbA1c mechanism that may increase liver cancer risk of those with diabetes is the impaired immune function that potentially interferes with neoplastic surveillance and encourages chronic viral infection, such as hepatitis B and hepatitis C.
The association between diabetes and gastric cancer mortality in our study is consistent with findings from several meta-analyses published between 2011 and 2013 showing an increased risk of gastric cancer in individuals with diabetes. A recent review has considered the multiple interacting factors that may play a role in this association [43]. These include possibly higher infection rates or prevalence of Helicobacter pylori infection, potentially increased salt intake, or increased detection rates among those with diabetes requiring routine medical care for a chronic condition.
There appear to be differential effects of elevated levels of glucose and insulin during the prodromal/early stages of diabetes on the development of prostate cancer, compared with the influence of fully developed diabetes on the risk of prostate cancer mortality [44]. Evidence from a metaanalysis [45] and two large cohort studies [46,47] showed lower prostate cancer incidence in individuals with longstanding diabetes compared to those without diabetes. The present study found that increased levels of insulin and insulin resistance (typical of subclinical diabetes) among those without diabetes were associated with lower risk of prostate cancer death. After the development of prostate cancer, however, the presence of diabetes seems to confer an increased risk of prostate cancer death. A recently published meta-analysis found an increased risk of prostate cancer mortality among those with diabetes [48]; however, a large collaborative study including a total of 2,217 prostate cancer deaths [49] found a nonsignificant decreased risk of mortality in individuals with diabetes. The present study provides further evidence for diabetes as a risk factor for prostate cancer death. Overall, the relationship of diabetes and its physiologic changes with prostate cancer incidence and mortality is complex [50,51]. Future research is needed to more accurately consider the time course of these physiologic changes, as well as the social and medical influences on the relationship between prostate cancer development and associated mortality.
Interpretation of this epidemiologic evidence should be undertaken with caution considering the multiple and possibly countervailing influences of factors difficult to track over the relevant exposure period. These could include the direct effects of hyperglycemia as well as secondary physiologic responses to hyperglycemia (e.g., insulin and the effects of the insulin-like growth factor system) and treatment of diabetes. Although we adjusted for obesity, there could also be effects of obesity per se. The evaluation of the role of obesity on cancer is very complex, especially in a study that evaluates cancer mortality instead of cancer incidence. Indeed, at baseline, participants who later died of cancer were less likely to be obese compared to those without a subsequent cancer death. It is well known that the obesity of those with diabetes tends to ameliorate with advancing age; insulin levels decline with time as beta-cell function fails [52,53]. There is also evidence that the influence of the IGF system declines with age in a parallel fashion to that of insulin  [54]. It is challenging to confidently model or adjust for these multiple factors and their possible interactions. There are a number of effects, by which diabetes could increase the likelihood of death among those who have already developed cancer, such as increased cardiovascular complications and increased susceptibility to infections. The Strong Heart Study was initially designed to study the prevalence of known or suspected cardiovascular disease risk factors in American Indians and to assess their influence on cardiovascular morbidity and mortality (http://stron gheart.ouhsc.edu/). However, the long follow-up of Strong Heart Study participants has allowed us to secondarily evaluate the influence of some of these risk factors on cancer mortality. For this reason, we do not have information on certain confounders that are of obvious interest (e.g., family history of the studied cancers, stage at diagnosis, treatment, and some environmental exposures), making residual confounding likely. It was not possible to evaluate for potential confounding due to aspirin or non-steroidal anti-inflammatory agents because the only information available was on sporadic use. Since the analysis is based on cancer mortality, conclusions may have also been biased by differences in healthcare utilization, or access to more specialized services by those with diabetes, compared to those without diabetes. Finally, we lack incidence data, and particularly among tumors with high survival (e.g., breast or colorectal), the power to detect significant associations is limited by the small number of deaths from those cancers.
The strengths of this study include the prospective analysis, the availability of relevant covariates (including arsenic, an established carcinogen that was present in drinking water at levels above 10 lg/L, the current US EPA standard) [55] to adjust for confounding, and a population with a unique genetic background.
In a population with a high prevalence of diabetes and obesity and with unique environmental exposures, our data provide additional evidence of heightened risk of cancer mortality associated with diabetes and confirmation of increased risk of mortality due to stomach and liver cancers. Additional research is needed to evaluate the role of diabetes and obesity on cancer incidence in American Indian communities as well as the importance of diabetes prevention and control in reducing the burden of cancer incidence and mortality in the study population.