Abstract
Background: Although unopposed estrogen exposure is considered the main driver of endometrial carcinogenesis, factors associated with states of insulin resistance and hyperinsulinemia are independently associated with endometrial cancer risk. We used dietary insulin load and insulin index scores to represent the estimated insulin demand of overall diets and assessed their association with endometrial cancer risk in the prospective Nurses' Health Study.
Methods: We estimated incidence rate ratios (RR) and 95% confidence intervals (CI) for risk of invasive endometrial cancer using Cox proportional hazards models. Between the baseline dietary questionnaire (1980) and 2010, we identified a total of 798 incident-invasive epithelial endometrial adenocarcinomas over 1,417,167 person-years of follow-up.
Results: Dietary insulin scores were not associated with overall risk of endometrial cancer. Comparing women in the highest with the lowest quintile, the multivariable-adjusted RRs of endometrial cancer were 1.07 (95% CI, 0.84–1.35) for cumulative average dietary insulin load and 1.03 (95% CI, 0.82–1.31) for cumulative average dietary insulin index. Findings did not vary substantially by alcohol consumption, total dietary fiber intake, or body mass index and/or physical activity (Pheterogeneity ≥ 0.10).
Conclusions: Intake of a diet predicted to stimulate a high postprandial insulin response was not associated with endometrial cancer risk in this large prospective study. Considering the complex interplay of diet, lifestyle, and genetic factors contributing to the hyperinsulinemic state, dietary measures alone may not sufficiently capture absolute long-term insulin exposure.
Impact: This study is the first to investigate dietary insulin scores in relation to endometrial cancer risk. Cancer Epidemiol Biomarkers Prev; 23(8); 1512–20. ©2014 AACR.
This article is featured in Highlights of This Issue, p. 1443
Introduction
Two major risk factors for endometrial cancer include obesity and use of estrogen-only postmenopausal hormone therapy (HT). These established risk factors, in addition to several others [e.g., parity, oral contraceptive (OC) use, smoking], correspond with alterations in sex steroid hormones such that the predominant driver of endometrial carcinogenesis is thought to be a relative increase in exposure to estrogen unopposed by progesterone (1, 2). However, endometrial cancer risk factors such as polycystic ovarian syndrome, type 2 diabetes, and the accumulation of visceral fat also are associated with insulin resistance and hyperinsulinemia (3–8). Up to one third of non-diabetic endometrial cancer patients exhibit insulin resistance based on an assessment of basal insulinemia and glycemia (9). Epidemiologic evidence suggests that individuals with higher levels of insulin or C-peptide, a stable marker of insulin secretion (10, 11), and those with evidence of metabolic syndrome are at increased endometrial cancer risk (12–18). Factor analysis of various prediagnostic plasma hormones, binding proteins, and cytokines identified an “insulin resistance/metabolic syndrome” component associated with postmenopausal endometrial cancer risk (19). Finally, insulin induces proliferation of endometrial cancer cell lines in vitro (20, 21), suggesting that increased insulin signaling may promote endometrial cancer development and/or progression.
Western diets, which consist of high-fat and highly processed carbohydrate-rich foods, are substantially more insulinogenic (i.e., greater insulin secretion per gram of food) compared with traditional diets based on less refined foods (22). Type, amount, and digestibility of dietary carbohydrate intake have direct physiologic effects on circulating insulin levels (22–24), which are highly correlated with postprandial blood glycemia (r = 0.70; P < 0.001; ref. 22). Given the putative link between insulin signaling and endometrial tumor growth, frequent consumption of foods associated with elevated insulin or blood glucose response has been hypothesized to increase endometrial cancer risk.
Epidemiologic studies have investigated dietary carbohydrate quality (glycemic index, GI) and/or a measure of both carbohydrate quality and quantity (glycemic load, GL) as surrogates of insulin and blood glucose levels with respect to endometrial cancer risk. Except for the Prostate, Lung, Colorectal, and Ovarian (PLCO) Cancer Screening Trial (25), prior studies observed nonsignificant elevations in endometrial cancer risk among women in the highest category of GL compared with the lowest. This translated into a modest, but significant approximately 20% elevation in risk associated with a high GL diet in a recent meta-analysis (26). However, although postprandial blood glycemia from carbohydrate consumption is highly correlated with circulating insulin levels, protein and fat can induce insulin secretion without raising blood glucose (22). Thus, quantifying the postprandial insulin response for various food items, including those with low or no carbohydrate content, may address the insulin hypothesis more directly. In this analysis, we used novel dietary insulin index (II) and insulin load (IL) scores developed for the Nurses' Health Study (NHS) cohort to investigate prospectively whether diets high in insulinogenic foods are associated with endometrial cancer risk.
Materials and Methods
Study population
The NHS is an ongoing prospective cohort following 121,700 female registered nurses ages 30 to 55 years from 11 U.S. states at enrollment in 1976. At baseline and biennially thereafter, participants completed self-administered questionnaires providing detailed information on anthropometric, lifestyle, menstrual, and reproductive factors. Participants also report their medical history from which newly diagnosed cancers and other diseases are identified. Follow-up of the cohort is high, with >90% of total possible person-years. Vital status was ascertained through next-of-kin, the U.S. Postal Service, and the National Death Index. These methods have identified an estimated 98% of deaths in the cohort (27). Completion of the self-administered questionnaire was considered to imply informed consent. The NHS protocol was approved by the Human Research Committee of Brigham and Women's Hospital (Boston, MA).
Case ascertainment
On each questionnaire, women reported whether they had been diagnosed with endometrial cancer during the previous 2 years. We then sought permission to obtain the relevant medical records and pathology reports. Study physicians blinded to questionnaire data reviewed the documents to verify and establish the date of diagnosis. As <10% of NHS endometrial cancer cases were diagnosed with serous, clear cell, or other rare histologic types, we restricted our analysis to cases with invasive endometrial adenocarcinoma (stages IA–IV) diagnosed between return of the 1980 dietary questionnaire and June 1, 2010. Over 1,417,167 person-years of follow-up, we identified 798 incident invasive epithelial endometrial adenocarcinomas.
Dietary assessment
Dietary information was first assessed in 1980 when a 61 semiquantitative food frequency questionnaire (FFQ) was sent to participants. Expanded questionnaires inquiring about annual average intake of over 130 individual foods and 22 beverages were sent in 1984, 1986, and every 4 years thereafter. FFQs asked, on average, how often during the previous year participants consumed each food or beverage using standard portion sizes. Nine responses were possible, ranging from never or less than once per month to 6 or more times per day. Nutrient intakes were calculated by multiplying the portion size of a single serving of each food by its reported frequency of intake, then multiplying the amount consumed by the nutrient content of the food, and summing the nutrient contributions of all food items using the U.S. Department of Agriculture food composition data (28). Reproducibility of the FFQ has been reported elsewhere (29–31).
II values for individual food items were obtained from published estimates (31 foods; refs. 22, 32) or provided by Dr. Jennie Brand-Miller of the University of Sydney, Australia (73 foods). U.S. food samples were shipped to the laboratory in Sydney for testing as previously described (22). The food II value was calculated by dividing the area under the insulin response curve for 1,000 kilojoules of the reference food (glucose). The II value for each test food represents the mean response of 11 to 13 subjects. On the basis of these data, we built an II database for a large number of foods queried on our FFQs. II values for individual foods ranged from 2.0 (butter) to 120.0 (jelly beans). For foods not directly tested, II values were either assigned values from similar foods (meats, fish, poultry, fats, and oils) or calculated using an adjustment based on percentage of carbohydrates/1,000 kilojoule (cereal grains, baked goods, sweets and snack foods, beverages, fruits, and dairy foods), an inverse water ratio (nonstarchy vegetables), or inverse dietary fiber content ratio (starchy vegetables).
We calculated the average dietary IL for each participant by multiplying the II value of each food and beverage by its energy content, and summing values for all items reported [∑(food II × kilocalories per serving × servings per day)]. Each unit of dietary IL represents the equivalent insulin response generated by 1 kilocalorie of glucose. The dietary II for the overall diet, which is the weighted mean of II values for each of the component foods, was calculated by dividing dietary IL by total energy intake [IL/∑(kilocalories per serving × servings per day)]. For secondary analyses, we calculated average dietary IL and II scores excluding the estimated insulin response to alcoholic beverages (beer, wine, and liquor).
Bao et al. assessed validity of dietary IL in predicting the actual insulin response (33). Healthy participants (n = 10 or 11 for each meal) consumed 13 different isoenergetic (2,000 kilojoules) mixed meals of varying macronutrient content. Despite limited statistical power, dietary IL strongly correlated with the observed postprandial insulin response (r = 0.78, P = 0.002), suggesting that dietary IL may be a valid measure of actual insulin response to composite meals.
Assessment of non-dietary factors
On the questionnaire administered in 1976, women reported current weight, height, smoking history, age at menarche, OC use, pregnancy history, menopausal status, age at menopause, use of postmenopausal HT, and personal history of other diseases. OC use and parity were asked on each questionnaire until 1984, whereas most other covariate data (except height and age at menarche) were updated on each biennial questionnaire. Type of postmenopausal hormone used (i.e., estrogen alone or estrogen with progesterone) was first asked in 1978 and on all subsequent questionnaires. Weight at age 18 was reported on the 1980 questionnaire. First-degree family history of endometrial cancer was collected in 1996 and 2008. Family history of colorectal cancer was first collected in 1982, updated in 1988 and biennially thereafter. Starting in 1980, participants were asked hours per week of overall moderate and vigorous physical activity. Hours spent on specific leisure-time activities were first reported in 1986 and updated every 2 to 4 years. From the frequency, duration, and intensity of reported activities, we calculated total metabolic equivalent (MET) hours of activity per week (34, 35). Body mass index (BMI) was calculated as weight in kilograms divided by the square of height in meters (kg/m2). Pack-years of smoking were calculated by multiplying the duration and dose of smoking; 1 pack-year is equivalent to having smoked 20 cigarettes per day for 1 year.
Statistical analysis
We used the 1980 dietary questionnaire as baseline. We sequentially excluded women who did not respond to the 1980 FFQ, left an extensive number of items blank (>10) or reported implausible total energy intake (<500 or >3,500 kcal/day; N = 29,232), had been diagnosed with cancer before 1980 (other than nonmelanoma skin cancer; N = 3,678), were born before 1920 or had an unknown date of birth (N = 58), or had undergone hysterectomy with or without oophorectomy or experienced menopause due to pelvic irradiation as of the 1980 questionnaire (N = 20,625). Because obesity is an important risk factor for endometrial cancer, we also excluded women missing BMI (N = 138; these women could reenter the analysis in subsequent cycles once information on BMI was available) for a total analytic population of 67,969 women at baseline. Person-time of follow-up for each eligible participant started with the return date of the baseline questionnaire and ended with date of endometrial cancer diagnosis, diagnosis of another cancer (except nonmelanoma skin cancer), hysterectomy, pelvic irradiation, death from any cause, or end of follow-up (June 1, 2010), whichever occurred first.
Cox proportional hazards regression models estimated incidence rate ratios (RR) and 95% confidence intervals (CI) with age in months and 2-year questionnaire cycle as the time scale. Average dietary II, dietary IL, dietary GI, dietary GL, and total dietary fiber intake were energy-adjusted by the residual method (36). This addresses the question of whether the composition of the diet, independent of total energy intake, is most relevant to endometrial cancer risk. Participants were classified according to dietary insulin score tertiles or quintiles based on the distribution observed within the cohort. Women in the lowest quantile of each score served as the reference group in analyses.
To assess the importance of timing, we analyzed dietary insulin scores in three ways: (i) using 1980 baseline values to reflect past exposure, (ii) simple update, to assess most recent dietary exposure, and (iii) cumulatively updating the average dietary insulin scores with each successive questionnaire data to reflect long-term exposure. To illustrate the cumulative average method, we used the average of dietary IL from 1980 and 1984 to predict endometrial cancer risk from 1984 to 1986. Similarly, we averaged IL scores from 1980, 1984, and 1986 to predict risk from 1986 to 1988. Factors included in the final multivariable analyses are listed in table footnotes. We additionally assessed dietary GI (quintiles), dietary GL (quintiles), total dietary fiber intake (quintiles), BMI at age 18 (continuous), physical activity (<3, 3 to <9, 9 to <18, 18 to <27, 27+ METs/week, unknown), alcohol consumption (none, <5, 5–14.9, 15+ grams/day), coffee consumption (continuous), caffeine intake (quintiles), and personal history of diabetes (no, yes) as potential confounders. Because results remained similar, these factors were not included in final models.
We performed a number of sensitivity analyses. We used the 1984 dietary questionnaire as baseline because the FFQ included more food items than the 1980 FFQ. To minimize the possibility that preclinical disease influenced dietary assessment, we examined associations after incorporating a 4-year time lag between exposure assessment and endometrial cancer diagnosis (e.g., 1980 exposure to predict disease risk from 1984 to 1986). We repeated analyses in which the dietary insulin scores were not energy-adjusted by the residual method and with dietary insulin scores that excluded the alcohol component. We additionally conducted an analysis restricted to women who did not report a personal history of diabetes to reduce dietary changes occurring as a result of a diabetes diagnosis.
We used the Wald test to assess linear trends of dietary insulin score variables assigned the median value for each intake category. Likelihood ratio tests assessed heterogeneity by BMI (<30 vs. ≥30 kg/m2) and total physical activity (high vs. low, using the median as the cutpoint) strata. We additionally assessed heterogeneity by alcohol intake (drinkers vs. nondrinkers) and total dietary fiber intake (above vs. below the median). We tested the proportional hazards assumption by comparing the log likelihood of models with and without interaction terms between age and the dietary insulin score trend variables. The proportional hazards assumption was not violated (P ≥ 0.14). P values were based on two-sided tests and considered statistically significant at P < 0.05. All analyses were conducted using SAS version 9.3 (SAS Institute Inc.).
Results
Characteristics of women according to selected quintiles of dietary IL and II at the midpoint of follow-up (1994) are shown in Table 1. Women in higher categories of the dietary insulin scores were slightly older, had a higher intake of carbohydrates and total dietary fiber, but consumed less protein, total fat, alcohol, and caffeine, and were less likely to have ever smoked or taken OCs. Among women that had used OCs, duration of use was shorter among those in the highest categories of the dietary insulin scores. Parity, menopausal status, BMI, total physical activity, and family history of endometrial or colorectal cancer did not differ by categories of the dietary insulin scores. Cumulative average dietary insulin scores were highly correlated with each other (rs = 0.96) and with cumulative average dietary GL (rs = 0.75 with dietary IL; rs = 0.73 with dietary II), but only modestly correlated with cumulative average GI (rs = 0.37 and 0.39, respectively).
. | Dietary IL . | Dietary II . | ||||
---|---|---|---|---|---|---|
. | Quintile 1 (217–621; n = 9,558) . | Quintile 3 (668–706; n = 9,556) . | Quintile 5 (755–1,279; n = 9,795) . | Quintile 1 (13.9–38.8; n = 9,448) . | Quintile 3 (41.7–44.1; n = 9,568) . | Quintile 5 (47.0–80.9; n = 9,729) . |
Age (y)a | 59.3 (6.9) | 59.5 (7.2) | 60.6 (7.3) | 59.3 (6.9) | 59.4 (7.1) | 60.7 (7.3) |
Age at menarche (y) | 12.6 (1.4) | 12.5 (1.4) | 12.5 (1.3) | 12.6 (1.4) | 12.5 (1.4) | 12.5 (1.3) |
Ever use of OCs, % | 52 | 50 | 47 | 52 | 50 | 47 |
Duration of OC useb (mo) | 53.8 (48.5) | 51.4 (46.9) | 50.0 (47.9) | 53.7 (48.4) | 51.3 (46.5) | 50.0 (47.9) |
Parous, % | 93 | 94 | 94 | 93 | 95 | 94 |
Number of pregnanciesc | 3.1 (1.5) | 3.2 (1.5) | 3.1 (1.4) | 3.1 (1.5) | 3.2 (1.5) | 3.1 (1.4) |
Postmenopausal, % | 80 | 80 | 81 | 80 | 79 | 81 |
Age at menopaused (y) | 49.9 (3.7) | 50.1 (3.5) | 50.1 (3.5) | 49.9 (3.7) | 50.0 (3.6) | 50.1 (3.5) |
Ever use of postmenopausal hormonesd, % | 50 | 50 | 53 | 50 | 50 | 52 |
Never smoker, % | 30 | 45 | 50 | 29 | 45 | 50 |
Past smoker, % | 46 | 41 | 40 | 47 | 41 | 40 |
Current smoker, % | 24 | 14 | 10 | 23 | 14 | 10 |
BMI at age 18 (kg/m2) | 21.1 (3.4) | 21.2 (3.4) | 21.1 (3.4) | 21.2 (3.4) | 21.2 (3.3) | 21.2 (3.4) |
Current BMI (kg/m2) | 26.1 (5.2) | 26.6 (5.3) | 26.0 (5.0) | 26.0 (5.2) | 26.6 (5.3) | 26.1 (5.0) |
Total physical activity (METs/week) | 19.7 (25.7) | 19.3 (25.0) | 20.1 (23.1) | 19.9 (26.0) | 19.2 (24.5) | 19.7 (23.1) |
Dietary II | 35.7 (3.2) | 42.9 (1.1) | 50.1 (3.4) | 35.5 (3.1) | 42.9 (0.7) | 50.3 (3.2) |
Dietary IL | 570 (48) | 687 (11) | 810 (57) | 571 (50) | 687 (18) | 808 (58) |
GL | 85.8 (17.5) | 109 (14) | 131 (18) | 85.7 (17.4) | 109 (14) | 131 (17) |
GI | 50.4 (3.8) | 52.8 (2.9) | 54.7 (3.5) | 50.4 (3.7) | 52.8 (2.9) | 54.7 (3.5) |
Carbohydrates (grams/day) | 175 (30) | 210 (24) | 241 (26) | 175 (30) | 211 (24) | 241 (26) |
Protein (grams/day) | 75.9 (15.4) | 74.9 (11.9) | 70.6 (12.4) | 76.0 (15.4) | 74.9 (12.0) | 70.5 (12.4) |
Total fat (grams/day) | 59.7 (12.4) | 52.3 (9.6) | 42.7 (9.6) | 59.5 (12.5) | 52.3 (9.7) | 42.6 (9.5) |
Total dietary fiber (grams/day) | 17.7 (6.8) | 18.8 (5.5) | 19.5 (5.5) | 17.6 (6.4) | 18.9 (5.6) | 19.5 (5.6) |
Alcohol intake (grams/day) | 13.0 (15.4) | 3.7 (6.0) | 1.9 (4.0) | 13.5 (15.4) | 3.5 (5.6) | 1.8 (3.7) |
Caffeine intake (mg/day) | 309 (231) | 266 (216) | 223 (206) | 303 (229) | 267 (216) | 228 (209) |
Caloric intake (kcal/day) | 1,711 (550) | 1,777 (524) | 1,702 (508) | 1,707 (546) | 1,777 (524) | 1,699 (507) |
Personal history of diabetes, % | 6 | 6 | 4 | 6 | 6 | 4 |
Personal history of hypertension, % | 35 | 33 | 33 | 35 | 33 | 33 |
Family history of endometrial cancer, % | 3 | 3 | 3 | 3 | 3 | 3 |
Family history of colon cancer, % | 13 | 15 | 14 | 13 | 15 | 14 |
. | Dietary IL . | Dietary II . | ||||
---|---|---|---|---|---|---|
. | Quintile 1 (217–621; n = 9,558) . | Quintile 3 (668–706; n = 9,556) . | Quintile 5 (755–1,279; n = 9,795) . | Quintile 1 (13.9–38.8; n = 9,448) . | Quintile 3 (41.7–44.1; n = 9,568) . | Quintile 5 (47.0–80.9; n = 9,729) . |
Age (y)a | 59.3 (6.9) | 59.5 (7.2) | 60.6 (7.3) | 59.3 (6.9) | 59.4 (7.1) | 60.7 (7.3) |
Age at menarche (y) | 12.6 (1.4) | 12.5 (1.4) | 12.5 (1.3) | 12.6 (1.4) | 12.5 (1.4) | 12.5 (1.3) |
Ever use of OCs, % | 52 | 50 | 47 | 52 | 50 | 47 |
Duration of OC useb (mo) | 53.8 (48.5) | 51.4 (46.9) | 50.0 (47.9) | 53.7 (48.4) | 51.3 (46.5) | 50.0 (47.9) |
Parous, % | 93 | 94 | 94 | 93 | 95 | 94 |
Number of pregnanciesc | 3.1 (1.5) | 3.2 (1.5) | 3.1 (1.4) | 3.1 (1.5) | 3.2 (1.5) | 3.1 (1.4) |
Postmenopausal, % | 80 | 80 | 81 | 80 | 79 | 81 |
Age at menopaused (y) | 49.9 (3.7) | 50.1 (3.5) | 50.1 (3.5) | 49.9 (3.7) | 50.0 (3.6) | 50.1 (3.5) |
Ever use of postmenopausal hormonesd, % | 50 | 50 | 53 | 50 | 50 | 52 |
Never smoker, % | 30 | 45 | 50 | 29 | 45 | 50 |
Past smoker, % | 46 | 41 | 40 | 47 | 41 | 40 |
Current smoker, % | 24 | 14 | 10 | 23 | 14 | 10 |
BMI at age 18 (kg/m2) | 21.1 (3.4) | 21.2 (3.4) | 21.1 (3.4) | 21.2 (3.4) | 21.2 (3.3) | 21.2 (3.4) |
Current BMI (kg/m2) | 26.1 (5.2) | 26.6 (5.3) | 26.0 (5.0) | 26.0 (5.2) | 26.6 (5.3) | 26.1 (5.0) |
Total physical activity (METs/week) | 19.7 (25.7) | 19.3 (25.0) | 20.1 (23.1) | 19.9 (26.0) | 19.2 (24.5) | 19.7 (23.1) |
Dietary II | 35.7 (3.2) | 42.9 (1.1) | 50.1 (3.4) | 35.5 (3.1) | 42.9 (0.7) | 50.3 (3.2) |
Dietary IL | 570 (48) | 687 (11) | 810 (57) | 571 (50) | 687 (18) | 808 (58) |
GL | 85.8 (17.5) | 109 (14) | 131 (18) | 85.7 (17.4) | 109 (14) | 131 (17) |
GI | 50.4 (3.8) | 52.8 (2.9) | 54.7 (3.5) | 50.4 (3.7) | 52.8 (2.9) | 54.7 (3.5) |
Carbohydrates (grams/day) | 175 (30) | 210 (24) | 241 (26) | 175 (30) | 211 (24) | 241 (26) |
Protein (grams/day) | 75.9 (15.4) | 74.9 (11.9) | 70.6 (12.4) | 76.0 (15.4) | 74.9 (12.0) | 70.5 (12.4) |
Total fat (grams/day) | 59.7 (12.4) | 52.3 (9.6) | 42.7 (9.6) | 59.5 (12.5) | 52.3 (9.7) | 42.6 (9.5) |
Total dietary fiber (grams/day) | 17.7 (6.8) | 18.8 (5.5) | 19.5 (5.5) | 17.6 (6.4) | 18.9 (5.6) | 19.5 (5.6) |
Alcohol intake (grams/day) | 13.0 (15.4) | 3.7 (6.0) | 1.9 (4.0) | 13.5 (15.4) | 3.5 (5.6) | 1.8 (3.7) |
Caffeine intake (mg/day) | 309 (231) | 266 (216) | 223 (206) | 303 (229) | 267 (216) | 228 (209) |
Caloric intake (kcal/day) | 1,711 (550) | 1,777 (524) | 1,702 (508) | 1,707 (546) | 1,777 (524) | 1,699 (507) |
Personal history of diabetes, % | 6 | 6 | 4 | 6 | 6 | 4 |
Personal history of hypertension, % | 35 | 33 | 33 | 35 | 33 | 33 |
Family history of endometrial cancer, % | 3 | 3 | 3 | 3 | 3 | 3 |
Family history of colon cancer, % | 13 | 15 | 14 | 13 | 15 | 14 |
NOTE: Values are mean (SD) or percentages and are standardized to the age distribution of the study population.
aValue is not age standardized.
bAmong ever OC users.
cAmong parous women.
dAmong postmenopausal women.
Past, recent, and long-term measures of dietary insulin scores were not associated with endometrial cancer risk in age- or multivariable-adjusted analyses. Women in the highest quintile of cumulative average dietary IL had an RR of 1.07 (95% CI, 0.84–1.35; Ptrend = 0.51; Table 2) compared with women in the lowest quintile after adjusting for established and suspected endometrial cancer risk factors. The comparable RR for cumulative average dietary II was 1.03 (95% CI, 0.81–1.31; Ptrend = 0.59). We observed similar associations using the 1984 FFQ as baseline. Results remained essentially unchanged after incorporating a 4-year time lag between exposure assessment and diagnosis. Dietary insulin scores were not associated with endometrial cancer risk after additionally adjusting for dietary GI, dietary GL, total dietary fiber intake, BMI at age 18, total physical activity, personal history of diabetes, hypertension, caffeine intake, alcohol consumption, or coffee consumption. Similar results were obtained when the dietary insulin scores were not energy-adjusted by the residual method, when the alcohol component was removed from the scores, or when analyses were restricted to women without a personal history of diabetes (data not shown). Estimates did not differ by alcohol drinking status (Pheterogeneity ≥ 0.71).
. | Quintile 1 . | Quintile 2 . | Quintile 3 . | Quintile 4 . | Quintile 5 . | Ptrend . |
---|---|---|---|---|---|---|
Dietary IL, RR (95% CI) | ||||||
Baseline (1980) | ||||||
Cases, N | 140 | 173 | 163 | 156 | 166 | |
Person-years | 280,015 | 285,065 | 289,085 | 285,454 | 277,547 | |
Age-adjusted | 1.00 (ref) | 1.25 (1.00–1.57) | 1.15 (0.92–1.45) | 1.14 (0.90–1.43) | 1.22 (0.97–1.53) | 0.21 |
Multivariablea | 1.00 (ref) | 1.17 (0.93–1.47) | 1.03 (0.82–1.30) | 1.00 (0.79–1.26) | 1.03 (0.82–1.30) | 0.78 |
Simple update | ||||||
Cases, N | 139 | 160 | 179 | 165 | 155 | |
Person-years | 282,179 | 283,090 | 285,420 | 284,570 | 281,908 | |
Age-adjusted | 1.00 (ref) | 1.18 (0.94–1.48) | 1.31 (1.05–1.63) | 1.20 (0.95–1.50) | 1.07 (0.85–1.35) | 0.50 |
Multivariablea | 1.00 (ref) | 1.08 (0.86–1.36) | 1.16 (0.93–1.46) | 1.08 (0.86–1.36) | 0.99 (0.79–1.26) | 0.99 |
Cumulative average | ||||||
Cases, N | 131 | 169 | 150 | 186 | 162 | |
Person-years | 275,154 | 285,596 | 289,088 | 288,057 | 279,272 | |
Age-adjusted | 1.00 (ref) | 1.27 (1.01–1.60) | 1.12 (0.88–1.41) | 1.39 (1.11–1.74) | 1.21 (0.96–1.52) | 0.08 |
Multivariablea | 1.00 (ref) | 1.14 (0.90–1.43) | 0.97 (0.76–1.23) | 1.22 (0.97–1.54) | 1.07 (0.84–1.35) | 0.51 |
Dietary II, RR (95% CI) | ||||||
Baseline (1980) | ||||||
Cases, N | 144 | 161 | 171 | 166 | 156 | |
Person-years | 276,530 | 285,550 | 287,460 | 287,715 | 279,912 | |
Age-adjusted | 1.00 (ref) | 1.10 (0.88–1.38) | 1.16 (0.92–1.44) | 1.14 (0.91–1.43) | 1.10 (0.87–1.38) | 0.38 |
Multivariablea | 1.00 (ref) | 1.04 (0.83–1.31) | 1.06 (0.85–1.33) | 1.03 (0.82–1.29) | 0.98 (0.77–1.23) | 0.81 |
Simple update | ||||||
Cases, N | 132 | 167 | 169 | 171 | 159 | |
Person-years | 280,524 | 282,222 | 284,101 | 286,213 | 284,108 | |
Age-adjusted | 1.00 (ref) | 1.28 (1.02–1.61) | 1.31 (1.04–1.65) | 1.31 (1.04–1.64) | 1.15 (0.91–1.45) | 0.23 |
Multivariablea | 1.00 (ref) | 1.17 (0.93–1.47) | 1.18 (0.93–1.49) | 1.19 (0.94–1.50) | 1.07 (0.84–1.35) | 0.56 |
Cumulative average | ||||||
Cases, N | 133 | 168 | 156 | 181 | 160 | |
Person-years | 271,761 | 285,369 | 287,736 | 287,945 | 284,357 | |
Age-adjusted | 1.00 (ref) | 1.24 (0.99–1.56) | 1.13 (0.90–1.43) | 1.36 (1.08–1.70) | 1.14 (0.90–1.44) | 0.18 |
Multivariablea | 1.00 (ref) | 1.09 (0.86–1.37) | 1.00 (0.79–1.26) | 1.19 (0.94–1.50) | 1.03 (0.81–1.31) | 0.59 |
. | Quintile 1 . | Quintile 2 . | Quintile 3 . | Quintile 4 . | Quintile 5 . | Ptrend . |
---|---|---|---|---|---|---|
Dietary IL, RR (95% CI) | ||||||
Baseline (1980) | ||||||
Cases, N | 140 | 173 | 163 | 156 | 166 | |
Person-years | 280,015 | 285,065 | 289,085 | 285,454 | 277,547 | |
Age-adjusted | 1.00 (ref) | 1.25 (1.00–1.57) | 1.15 (0.92–1.45) | 1.14 (0.90–1.43) | 1.22 (0.97–1.53) | 0.21 |
Multivariablea | 1.00 (ref) | 1.17 (0.93–1.47) | 1.03 (0.82–1.30) | 1.00 (0.79–1.26) | 1.03 (0.82–1.30) | 0.78 |
Simple update | ||||||
Cases, N | 139 | 160 | 179 | 165 | 155 | |
Person-years | 282,179 | 283,090 | 285,420 | 284,570 | 281,908 | |
Age-adjusted | 1.00 (ref) | 1.18 (0.94–1.48) | 1.31 (1.05–1.63) | 1.20 (0.95–1.50) | 1.07 (0.85–1.35) | 0.50 |
Multivariablea | 1.00 (ref) | 1.08 (0.86–1.36) | 1.16 (0.93–1.46) | 1.08 (0.86–1.36) | 0.99 (0.79–1.26) | 0.99 |
Cumulative average | ||||||
Cases, N | 131 | 169 | 150 | 186 | 162 | |
Person-years | 275,154 | 285,596 | 289,088 | 288,057 | 279,272 | |
Age-adjusted | 1.00 (ref) | 1.27 (1.01–1.60) | 1.12 (0.88–1.41) | 1.39 (1.11–1.74) | 1.21 (0.96–1.52) | 0.08 |
Multivariablea | 1.00 (ref) | 1.14 (0.90–1.43) | 0.97 (0.76–1.23) | 1.22 (0.97–1.54) | 1.07 (0.84–1.35) | 0.51 |
Dietary II, RR (95% CI) | ||||||
Baseline (1980) | ||||||
Cases, N | 144 | 161 | 171 | 166 | 156 | |
Person-years | 276,530 | 285,550 | 287,460 | 287,715 | 279,912 | |
Age-adjusted | 1.00 (ref) | 1.10 (0.88–1.38) | 1.16 (0.92–1.44) | 1.14 (0.91–1.43) | 1.10 (0.87–1.38) | 0.38 |
Multivariablea | 1.00 (ref) | 1.04 (0.83–1.31) | 1.06 (0.85–1.33) | 1.03 (0.82–1.29) | 0.98 (0.77–1.23) | 0.81 |
Simple update | ||||||
Cases, N | 132 | 167 | 169 | 171 | 159 | |
Person-years | 280,524 | 282,222 | 284,101 | 286,213 | 284,108 | |
Age-adjusted | 1.00 (ref) | 1.28 (1.02–1.61) | 1.31 (1.04–1.65) | 1.31 (1.04–1.64) | 1.15 (0.91–1.45) | 0.23 |
Multivariablea | 1.00 (ref) | 1.17 (0.93–1.47) | 1.18 (0.93–1.49) | 1.19 (0.94–1.50) | 1.07 (0.84–1.35) | 0.56 |
Cumulative average | ||||||
Cases, N | 133 | 168 | 156 | 181 | 160 | |
Person-years | 271,761 | 285,369 | 287,736 | 287,945 | 284,357 | |
Age-adjusted | 1.00 (ref) | 1.24 (0.99–1.56) | 1.13 (0.90–1.43) | 1.36 (1.08–1.70) | 1.14 (0.90–1.44) | 0.18 |
Multivariablea | 1.00 (ref) | 1.09 (0.86–1.37) | 1.00 (0.79–1.26) | 1.19 (0.94–1.50) | 1.03 (0.81–1.31) | 0.59 |
aMultivariable model adjusted for age at menarche (<12, 12, >12 years of age, unknown), age at last birth and number of pregnancies (nulliparous, <30 years and 1–2 children, 30+ years and 1–2 children, <30 years and 3–4 children, 30+ years and 3–4 children, 5 or more children), OC use (never use, <13, 13–36, 37–72, >72 months), age at menopause (<45, 45–46, 47–48, 49–50, 51–52, >52 years), menopause and HT use status (premenopausal, postmenopausal-never HT user, -past HT use, -current E only, -current E + P use, unknown), smoking status and pack-years of smoking (never smoker, past/<20 pack-years, past/21–40 pack-years, past/>40 pack-years, current/<20 pack-years, current/21–40 pack-years, current/>40 pack-years, unknown), family history of colorectal cancer (no, yes), family history of endometrial cancer (no, yes), total caloric intake (continuous), hypertension (no, yes), and BMI (continuous). E refers to estrogen-only hormone therapy use; E + P refers to estrogen plus progesterone hormone therapy use.
Because higher dietary fiber content lowers insulin response to foods (37), we stratified our analysis by total dietary fiber intake (above vs. below median intake). Cumulative average dietary insulin scores were not associated with endometrial cancer risk among women with higher total dietary fiber intake (Ptrend ≥ 0.10). Among women who consumed less total dietary fiber, those in the highest tertiles of cumulative average dietary insulin scores had significant 35% to 40% increased risks of endometrial cancer (Ptrend = 0.03) compared with women in the lowest tertiles. However, heterogeneity by total dietary fiber intake was not significant (Pheterogeneity ≥ 0.77).
Food intake among insulin-resistant individuals stimulates greater insulin production to overcome the decreased responsiveness of target tissues (38). Obese or less active individuals are more likely to exhibit insulin resistance (39) and, therefore, may experience the greatest detriment from a high insulinogenic diet. We examined whether the association with cumulative average dietary insulin scores differed by BMI and total physical activity. We did not observe a positive association between cumulative average dietary insulin scores and endometrial cancer risk among obese women or among women who reported total physical activity less than the median of the study population. Likewise, no association was observed in the group of women expected to have the highest proportion of insulin-resistant individuals (obese and low physical activity; Table 3). We observed a marginal increased risk associated with the highest compared with lowest tertile of cumulative average dietary II among nonobese women with high levels of physical activity (RR, 1.37; 95% CI, 0.99–1.90; Ptrend = 0.06; Pheterogeneity = 0.26). However, the same trend was not observed using the 1984 FFQ as baseline or in the 4-year lagged analysis, suggesting the observation is due to chance. Finally, because postmenopausal HT use may alter insulin responsiveness (12, 40), we conducted an analysis restricted to women who never used postmenopausal hormones. Cumulative average dietary insulin scores were not associated with endometrial cancer risk among these women (data not shown).
. | . | Dietary IL (tertiles) . | . | Dietary II (tertiles) . | . | ||||
---|---|---|---|---|---|---|---|---|---|
. | . | RRa (95% CI) . | . | RRa (95% CI) . | . | ||||
. | Cases, N . | 1 . | 2 . | 3 . | Ptrend . | 1 . | 2 . | 3 . | Ptrend . |
BMI <30 kg/m2 | 500 | 1.00 | 0.94 (0.75–1.17) | 1.04 (0.84–1.30) | 0.82 | 1.00 | 1.06 (0.85–1.33) | 1.19 (0.95–1.49) | 0.15 |
BMI > = 30 kg/m2 | 298 | 1.00 | 1.01 (0.75–1.36) | 1.10 (0.81–1.49) | 0.57 | 1.00 | 1.09 (0.81–1.46) | 1.06 (0.78–1.44) | 0.71 |
High physical activity | 341 | 1.00 | 0.88 (0.67–1.16) | 1.07 (0.82–1.40) | 0.72 | 1.00 | 1.07 (0.81–1.40) | 1.19 (0.91–1.56) | 0.22 |
Low physical activity | 414 | 1.00 | 1.01 (0.79–1.29) | 1.00 (0.77–1.28) | 0.97 | 1.00 | 1.04 (0.81–1.33) | 1.01 (0.79–1.31) | 0.96 |
BMI < 30 and high PA | 244 | 1.00 | 1.01 (0.73–1.41) | 1.15 (0.83–1.59) | 0.45 | 1.00 | 1.14 (0.81–1.59) | 1.37 (0.99–1.90) | 0.06 |
Intermediate group | 329 | 1.00 | 0.81 (0.61–1.06) | 0.83 (0.63–1.09) | 0.17 | 1.00 | 0.97 (0.74–1.28) | 0.89 (0.67–1.18) | 0.38 |
BMI ≥ 30 and low PA | 182 | 1.00 | 1.34 (0.89–2.01) | 1.35 (0.88–2.05) | 0.17 | 1.00 | 1.27 (0.86–1.89) | 1.22 (0.81–1.83) | 0.35 |
. | . | Dietary IL (tertiles) . | . | Dietary II (tertiles) . | . | ||||
---|---|---|---|---|---|---|---|---|---|
. | . | RRa (95% CI) . | . | RRa (95% CI) . | . | ||||
. | Cases, N . | 1 . | 2 . | 3 . | Ptrend . | 1 . | 2 . | 3 . | Ptrend . |
BMI <30 kg/m2 | 500 | 1.00 | 0.94 (0.75–1.17) | 1.04 (0.84–1.30) | 0.82 | 1.00 | 1.06 (0.85–1.33) | 1.19 (0.95–1.49) | 0.15 |
BMI > = 30 kg/m2 | 298 | 1.00 | 1.01 (0.75–1.36) | 1.10 (0.81–1.49) | 0.57 | 1.00 | 1.09 (0.81–1.46) | 1.06 (0.78–1.44) | 0.71 |
High physical activity | 341 | 1.00 | 0.88 (0.67–1.16) | 1.07 (0.82–1.40) | 0.72 | 1.00 | 1.07 (0.81–1.40) | 1.19 (0.91–1.56) | 0.22 |
Low physical activity | 414 | 1.00 | 1.01 (0.79–1.29) | 1.00 (0.77–1.28) | 0.97 | 1.00 | 1.04 (0.81–1.33) | 1.01 (0.79–1.31) | 0.96 |
BMI < 30 and high PA | 244 | 1.00 | 1.01 (0.73–1.41) | 1.15 (0.83–1.59) | 0.45 | 1.00 | 1.14 (0.81–1.59) | 1.37 (0.99–1.90) | 0.06 |
Intermediate group | 329 | 1.00 | 0.81 (0.61–1.06) | 0.83 (0.63–1.09) | 0.17 | 1.00 | 0.97 (0.74–1.28) | 0.89 (0.67–1.18) | 0.38 |
BMI ≥ 30 and low PA | 182 | 1.00 | 1.34 (0.89–2.01) | 1.35 (0.88–2.05) | 0.17 | 1.00 | 1.27 (0.86–1.89) | 1.22 (0.81–1.83) | 0.35 |
Abbreviation: PA, physical activity.
aMultivariable model adjusted for age at menarche (<12, 12, >12 years of age, unknown), age at last birth and number of pregnancies (nulliparous, <30 years and 1–2 children, 30+ years and 1–2 children, <30 years and 3–4 children, 30+ years and 3–4 children, 5 or more children), OC use (never use, <13, 13–36, 37–72, >72 months), age at menopause (<45, 45–46, 47–48, 49–50, 51–52, >52 years), menopause and HT use status (premenopausal, postmenopausal-never HT user, -past HT use, -current E only, -current E + P use, unknown), smoking status and pack-years of smoking (never smoker, past/<20 pack-years, past/21–40 pack-years, past/>40 pack-years, current/<20 pack-years, current/21–40 pack-years, current/>40 pack-years, unknown), family history of colorectal cancer (no, yes), family history of endometrial cancer (no, yes), total caloric intake (continuous), hypertension (no, yes), and BMI (continuous). All P values for interaction ≥0.10. E refers to estrogen-only hormone therapy use; E + P refers to estrogen plus progesterone hormone therapy use.
Discussion
We used novel dietary insulin scores to assess the impact of diets high in insulinogenic foods on endometrial cancer risk. We did not find evidence for increased endometrial cancer risk associated with dietary IL or II in this large cohort of U.S. women. Further, the scores were not associated with elevated risk among heavy and/or less active individuals, subgroups of women who collectively experience greater postprandial insulin response due to a higher proportion of insulin-resistant individuals (39).
Prior studies assessed endometrial cancer risk using dietary GL and GI scores to reflect the postprandial glucose response of carbohydrate-containing foods and to act as surrogates of insulin response to different diets. With the exception of PLCO, which observed significant reduced risk among women who consumed high GL diets (RR, 0.63; 95% CI, 0.46–0.84 for highest vs. lowest quartile; ref. 25), most studies report nonsignificant elevations in endometrial cancer risk with GL or GI (26). Cohort studies observed elevated risks of 15% to 36% with high GL diets (pooled RR, 1.22; 95% CI, 1.09–1.37 for highest vs. lowest category), including the NHS (RR, 1.29; 95% CI, 0.99–1.67; ref. 41), but not with GI (pooled RR, 1.00; 95% CI, 0.87–1.14; ref. 26). In contrast, case–control studies observed significant increased endometrial cancer risk with GI (pooled OR, 1.56; 95% CI, 1.21–2.02), but the relationship with high GL diets was weaker (pooled OR, 1.14; 95% CI, 0.91–1.44) compared with cohort studies (26).
To address the insulin hypothesis more directly, we used novel dietary insulin scores to quantify the postprandial insulin response. The food II, on which the scores were based, was developed under highly standardized conditions (22). Dietary IL accurately predicted insulin response to mixed meals in healthy volunteers, demonstrating validity of the index (33). Moreover, in the NHS and Health Professionals Follow-up Study cohorts, dietary insulin scores positively correlated with plasma triglyceride levels, particularly among obese individuals (42), confirming the ability to predict an expected biologic response associated with insulin resistance. Contrary to our hypothesis, we did not find evidence of increased endometrial cancer risk associated with the dietary insulin scores. Given that imprecise FFQ assessments may have biased results toward the null, we examined long-term exposure using cumulative average dietary scores, which dampen variation due to measurement error. The scores remained unrelated to endometrial cancer risk.
Insulin resistance/metabolic syndrome is an emerging risk factor for endometrial cancer (14–19). Fasting insulin, a marker of insulin resistance (43), was associated with overall endometrial cancer risk in a Canadian population-based case–control study (44), and in subgroups of non-HT users and/or overweight/obese women in a Swedish population-based case–control study (45) and the Women's Health Initiative Observational Study (12). Because excess energy is a major contributor to insulin resistance and sustained elevated insulin levels (39), we hypothesized that insulin-resistant individuals (obese/inactive) who consume an insulinogenic diet would be at highest risk for endometrial cancer. However, we did not observe significant associations between dietary insulin scores and endometrial cancer risk among these individuals (Ptrend ≥ 0.17; Pheterogeneity ≥ 0.19).
Our study benefits from the large, prospective design in which dietary and risk factor information was updated every 2 to 4 years over follow-up. The prospective nature and high follow-up rate minimized recall and selection biases, whereas repeated assessments reduced measurement error. Known and suspected risk factors were controlled to minimize residual confounding, although confounding by unmeasured or unknown exposures cannot be ruled out. Given our extensive data, we conducted a number of sensitivity analyses to rule out chance findings or a potential effect of subclinical disease on dietary intake.
We recognize that our study has several limitations. Although dietary insulin scores were developed to assess total postprandial insulin response to food intake, the scores do not account for overall composition or frequency of food intake, which may affect insulin response. In developing the database, acute postprandial insulin responses were measured in young lean individuals (22), which may not reflect the absolute response in heavier, middle-aged women. However, the method is valid if the relative insulin response to each food is comparable between the two groups. It is important to consider that in addition to dietary factors that influence insulin secretion acutely, elevated insulin levels are determined by a host of factors including anthropometric characteristics throughout the lifecourse, physical activity, dietary factors that influence insulin resistance, and genetic predisposition (39, 46–50). Thus, the acute postprandial insulin response alone may not sufficiently reflect absolute exposure to insulin levels. Finally, our cohort is not representative of the general U.S. population. By design, the target population was chosen to maximize internal validity through enhanced reporting of exposures and health outcomes. As underlying biology is likely very similar across ethnic groups and social class, observations made in our cohort may apply more broadly.
In summary, our results do not support an association between dietary factors specifically related to the postprandial insulin response and endometrial cancer risk. An integrated surrogate measure of hyperinsulinemia that incorporates dietary, lifestyle, and genetic factors may provide a more useful indicator for studies investigating the insulin hypothesis in disease risk.
Disclosure of Potential Conflicts of Interest
No potential conflicts of interest were disclosed.
Disclaimer
The opinions or assertions contained in this work are the private views of the authors and are not considered as official or reflecting the views of the NIH. Certain data used in this publication were obtained from the Connecticut Department of Public Health (DPH). The authors assume full responsibility for analyses and interpretation of these data.
Authors' Contributions
Conception and design: J. Prescott, I. De Vivo
Development of methodology: Y. Bao, E.L. Giovannucci, I. De Vivo
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): S.E. Hankinson, I. De Vivo
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): J. Prescott, Y. Bao, A.N. Viswanathan, E.L. Giovannucci, S.E. Hankinson, I. De Vivo
Writing, review, and/or revision of the manuscript: J. Prescott, Y. Bao, A.N. Viswanathan, E.L. Giovannucci, S.E. Hankinson, I. De Vivo
Study supervision: I. De Vivo
Acknowledgments
The authors thank the NHS participants and cohort staff for their valuable contributions as well as the following state registries: AL, AZ, AR, CA, CO, CT, DE, FL, GA, ID, IL, IN, IA, KY, LA, ME, MD, MA, MI, NE, NH, NJ, NY, NC, ND, OH, OK, OR, PA, RI, SC, TN, TX, VA, WA, WY. This study was approved by the Human Investigations Committee of the DPH.
Grant Support
The NHS project, E.L. Giovannucci, and S.E. Hankinson are supported by the NIH (P01 CA87969). I. De Vivo (R01 CA082838) and J. Prescott (U54 CA155626) were supported by funds from the NIH.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.