Coffee is one of the most commonly consumed beverages worldwide. It contains over 1,000 chemical compounds, many of which, including caffeine, diterpenes, and phenolic acids, are known to be bioactive.1 Strong inverse associations of coffee drinking with liver cancer2 and endometrial cancer3,4 have been consistently observed. However, evidence of an association with other cancers, including breast,5 prostate,6 and lung7 is less consistent. In 1991, the International Agency for Research on Cancer (IARC) Monographs Programme found “limited evidence in humans that coffee drinking is carcinogenic in the urinary bladder” and subsequently classified coffee as “possibly carcinogenic” to this site (group 2B; Monographs Volume 51).8 According to IARC, limited evidence means that confounding, bias, and chance cannot be ruled out as potential explanations for the observed positive association. Subsequent studies have weakened rather than strengthened the evidence for a causal association. In 2016, the IARC Monographs Programme reviewed the accumulated evidence and determined that coffee drinking was “not classifiable as to its carcinogenicity to humans” (group 3; Monographs Volume 116).1 The 2016 Monograph also highlighted a need for large prospective studies with sufficient case numbers and detailed information on potential confounding factors, such as tobacco smoking, to better understand the association between coffee drinking and bladder cancer, and to help rule out chance, bias, and confounding as possible explanations for previous positive findings.
In the current study, we analyzed data from the NIH-AARP Study, a large prospective cohort study with 469,047 participants and 6,012 cases of bladder cancer, nearly eight times as many cases as the combined total of the most recent meta-analysis, which included 753 cases and 236,343 participants.9 Our large sample size allowed for stratification by important known risk factors including sex and tobacco smoking. Additionally, more than 15 years of follow-up and detailed information on cigarette smoking use, intensity, and time since quitting permitted an assessment of the impact of reverse causality and residual confounding by smoking, respectively.
METHODS
Study Population and Design
The NIH-AARP Diet and Health Study, which has been described previously,10 launched in 1995–1996 with the mailing of a self-administered questionnaire to 3.5 million AARP members ages 50 to 71 years who resided in one of six US states (California, Florida, Louisiana, New Jersey, North Carolina, and Pennsylvania) or two US metropolitan areas (Atlanta, Georgia and Detroit, Michigan). The baseline questionnaire queried about demographics, health-related behaviors, and dietary intake. The NIH-AARP Diet and Health Study was reviewed and approved by the Special Studies Institutional Review Board of the US National Cancer Institute.
Of the 566,398 participants who satisfactorily completed the questionnaire and provided informed consent, we excluded, in order, proxy-responders (n = 15,760); those with self-reported cancer (except nonmelanoma skin cancer) before baseline (n = 49,318); those with self-reported end-stage renal disease before baseline (n = 997); those with a registry cancer diagnosis before baseline (n = 2,016); those with only a death record for cancer (n = 4,253); those with missing information on coffee intake (n = 2,672); those with extremely low or high caloric intake (n = 4,167), defined as more than two interquartile ranges above the sex-specific 75th percentile or below the 25th percentile of intake; those who died or were diagnosed with bladder cancer on or before the date their questionnaire was received (n = 32); and those with missing information on cigarette smoking (n = 18,136). The resulting analytic cohort consisted of 279,290 men and 189,757 women (N = 469,047).
Cohort Follow-up
Follow-up time accumulated from the date the baseline questionnaire was returned (beginning 25 October 1995) until the date of first urinary system cancer, the date of death, the date the participant moved out of the catchment area, or the end of study follow-up (31 December 2011), whichever came first. Vital status was ascertained by linkage to the Social Security Administration Death Master File and response to mailings. Participant addresses were updated annually in response to change of address requests and by matching cohort participants to the US Post Office National Change of Address database.
Case Ascertainment
We identified incident bladder cancer cases by probabilistic linkage of the NIH-AARP Diet and Health Study cohort to cancer registries of the eight baseline states and three additional states (Arizona, Texas, and Nevada) to which participants were most likely to move during follow-up. A validation study of the eight baseline cancer registries estimated that this approach ascertained about 90% of all cancer cases.11 Bladder cancer cases were defined according to the International Classification of Diseases for Oncology (ICD-O, 3rd edition)12 by anatomic site C67.0–C67.9 and transitional cell (urothelial) morphology (ICD codes 8120, 8122, 8123, or 8130).
Exposure Assessment
The self-administered baseline questionnaire included a 124-item food frequency questionnaire (FFQ) that queried about usual coffee intake in the previous 12 months using 10 predefined frequency categories, ranging from none to ≥6 cups/d. Of those who reported drinking coffee, 96% (n = 452,386) provided information on whether they drank caffeinated or decaffeinated coffee more than half the time. We used responses to these two FFQ items to categorize coffee drinkers into prespecified categories ranging from none to ≥4 cups/d for total coffee and for either caffeinated or decaffeinated coffee more than half the time. For participants with missing information on coffee type, we created a missing category.
The baseline questionnaire also provided information on potential confounding factors. Smoking history was defined as never (<100 cigarettes during lifetime), former, or current smoker. Current and former smokers further reported smoking intensity using six categories of cigarettes/d (1–10, 11–20, 21–30, 31–40, 41–60, and ≥61), and former smokers additionally reported time since quitting using four categories (quit <1 year ago, quit 1–4 years ago, quit 5–9 years ago, and quit ≥10 years ago). We defined former smokers as those who had quit ≥1 year before baseline. Prior methodological studies have shown that the assessment of tobacco smoking via questionnaire is reproducible and valid.13,14 A separate question assessed regular (≥1 year) use of pipes or cigars. Average daily alcohol intake over the past 12 months was calculated from drinks of alcohol from beer, wine, and liquor. Nutrient and food group intakes were estimated using FFQ responses in conjunction with a nutrient and food group database developed using national dietary intake data from the 1994–1996 US Department of Agriculture’s Continuing Survey of Food Intake by Individuals.15 Physical activity during the past 12 months was defined as frequency of activity lasting ≥20 minutes that caused increases in breathing or heart rate or sweating. Body mass index (BMI) was calculated from self-reported height and weight at baseline.
Statistical Methods
We tabulated demographic and lifestyle factors, potentially associated with bladder cancer, by coffee intake separately among men and women because sex is an important risk factor for bladder cancer. We used Cox proportional hazards regression models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for categories of coffee intake and bladder cancer using nondrinkers of coffee as the referent group. We conducted tests for linear trend across categories of coffee intake by assigning participants the midpoint of their coffee intake category and entering this single variable into a separate model. Participants in the ≥4 cups/d category were assigned a value of 5 cups/d. The risk estimate for the continuous value (1 cup/d increase) was based on the original 10-category variable; each category was assigned its midpoint value, and participants in the highest category of intake (i.e., ≥6 cups/d) were assigned a value of 6 cups/d. In our primary analysis, we compared overall models, adjusted for age and sex, and sex-stratified models adjusted for age with those additionally adjusted for tobacco smoking and other potential confounding factors. Only age and cigarette smoking altered risk estimates by more than 5%. For the less than 5% of the cohort that was missing data for a particular covariate, an indicator variable for missing was included in the models. To test for sex heterogeneity and smoking heterogeneity in the association between coffee intake and bladder cancer risk, we compared multivariable models with and without the cross-product terms for either sex or each level of smoking (never, current, or former) and coffee intake (continuous) using the likelihood ratio test. Person-time was used as the underlying time metric; results calculated with age as the underlying time metric were similar. Consistent with the assumption of proportional hazards, inclusion of the time dependent interaction did not improve the model fit (P value for interaction = 0.98). We used SAS software, version 9.3 (SAS Institute, Cary, NC) to conduct our analyses.
To explore the impact of cigarette smoking on risk estimates, we stratified by detailed smoking categories, including categories of smoking dose among current smokers, and categories of smoking dose and time since quitting among former smokers. Next, we further explored the impact of residual confounding by cigarette smoking in a subset of participants (n = 250,716) who provided information on cigarette smoking use and intensity over the lifetime in the 2004–2005 follow-up questionnaire. We compared the HRs from the multivariable-adjusted baseline model to those additionally adjusted for cigarette smoking status and, if applicable, the number of cigarettes per day, on average, during the following age periods: <15, 15–19, 20–24, 25–29, 30–39, 40–49, 50–59, 60–60, and ≥70 years old.
In additional secondary analyses, we estimated the HRs for those who predominantly drank caffeinated coffee or decaffeinated coffee. Caffeinated and decaffeinated coffee variables, as well as an indicator variable for missing caffeine type were included in a single model, and nondrinkers of either type of coffee served as the reference group. Finally, we explored the potential for reverse causation by evaluating associations of coffee drinking with bladder cancer cases that occurred less than 5 years, 5 to less than 10 years, or 10 or more years after baseline.
RESULTS
During more than 6.3 million person-years of follow-up, 6,012 incident cases of bladder cancer (5,088 men and 924 women) were identified. The median age at baseline was 62.7 years and 62.3 years for men and women, respectively. Median follow-up time was 15.5 years. The cohort was predominantly non-Hispanic White (92%) and was well educated (67% of men and 55% of women had at least some college education). A majority of men (70%) and women (54%) had a history of smoking cigarettes; however, most, 57% of men and 38% of women, were former smokers. Coffee intake was positively correlated with smoking intensity among current smokers (Spearman correlation coefficient = 0.18; P < 0.0001). Approximately 91% of men and 88% of women reported drinking coffee and 60% of men and 52% of women reported drinking at least 2 cups/d. At baseline, individuals who reported higher coffee intakes were more likely to be men than women, non-Hispanic White than other races, and ever than never smokers (Table 1), each of which is a known risk factor for bladder cancer.
TABLE 1: Baseline Characteristics of the NIH-AARP Diet and Health Study by Sex and Coffee Intake (N = 469,047)
First, we tested the hypothesis that higher coffee intake is associated with increased risk of bladder cancer (Table 2). Overall, higher coffee intake was associated with greater risk of bladder cancer in age- and sex-adjusted models; the HR for the highest category of coffee intake (≥4 cups/d) relative to who did not drink coffee was 1.91 (95% CI = 1.70, 2.14; P trend < 0.0001). Following adjustment for smoking, a positive, albeit attenuated, association remained between coffee intake and bladder cancer (HR for ≥4 cups/d = 1.21, 95% CI = 1.07, 1.36; P trend = 0.0001). In sex-stratified analyses, adjustment for smoking similarly attenuated the positive association for ≥4 cups/d from a relative risk of 2.04 (95% CI = 1.79, 2.32) to 1.29 (95% CI = 1.13, 1.47) among men and from a relative risk of 1.57 (95% CI = 1.22, 2.02) to 0.98 (95% CI = 0.76, 1.27) among women. Despite the observed differences in risk estimates for men and women, we did not find evidence of effect modification by sex (P for sex heterogeneity = 0.92). Further adjustment for other potential confounders did not meaningfully alter risk estimates (Table 2).
TABLE 2: Hazard Ratios (95% Confidence Intervals) for Bladder Cancer According to Coffee Intake in the NIH-AARP Diet and Health Study by Sex
In analyses stratified by tobacco smoking status (Table 3), we found evidence that risk estimates varied by smoking status (P for smoking heterogeneity = 0.03). We observed evidence of positive associations among both lower intensity (i.e., ≤20 cigarettes/d) current smokers (P trend = 0.01) and higher intensity (i.e., >20 cigarettes/d) current smokers (P trend = 0.02), but we found no evidence of an association between coffee drinking and bladder cancer risk among never smokers (HR for ≥4 cups/d = 0.87, 95% CI = 0.65, 1.17; P trend = 0.84). No clear pattern was observed among former smokers, although most risk estimates were moderately above one. Additionally, the exclusion of participants who ever reported regularly smoking pipes or cigars did not meaningfully alter risk estimates (eTable 1; https://links.lww.com/EDE/B200). Next, we evaluated the association of coffee drinking and bladder cancer in the majority of participants with additional information on lifetime cigarette use. In Table 4, we show that adjustment for both baseline tobacco smoking variables (i.e., cigarette smoking status, intensity, time since quitting among former smokers, and cigar/pipe smoking status) and cigarette smoking use and intensity in each age group, including younger ages, further attenuated positive risk estimates. In the final multivariable-adjusted model, the relative risk for ≥4 cups/d was 1.09 (95% CI = 0.93, 1.29; P trend = 0.16).
TABLE 3: Hazard Ratios (95% Confidence Intervals) for Bladder Cancer According to Coffee Intake in the NIH-AARP Diet and Health Study by Smoking (N = 453,935)a
TABLE 4: Hazard Ratios (95% Confidence Intervals) for Bladder Cancer According to Coffee Intake in the NIH-AARP Diet and Health Study Among Those with Information on Lifetime Cigarette Smoking Intensity (N = 250,716)
Associations were similar for both caffeinated and decaffeinated coffee and were also considerably attenuated after adjustment for tobacco smoking (eTable 2; https://links.lww.com/EDE/B200). Nevertheless, a positive association between the highest level of coffee intake (≥4 cups/d), either caffeinated or decaffeinated, and bladder cancer risk (HR = 1.20, 95% CI = 1.06, 1.35 and HR = 1.18, 95% CI = 1.00, 1.38, respectively) persisted. Finally, our examination of risk estimates over follow-up time suggested slightly stronger associations for cases that were diagnosed within 5 years of baseline, although we found no evidence of a violation of the proportional hazard assumption (P for interaction = 0.98; eTable 3; https://links.lww.com/EDE/B200) indicating that risk estimates did not vary meaningfully by follow-up time.
DISCUSSION
This prospective analysis included, to our knowledge, the largest number of bladder cancer cases to date and examined the association of coffee drinking and bladder cancer risk stratified by cigarette smoking, an important risk factor for bladder cancer.16 Overall, we observed a modest increased risk of bladder cancer among coffee drinkers as compared with those who did not drink coffee. However, this positive association was restricted to ever smokers, particularly current smokers, and we observed no evidence of an association among never smokers. Previous prospective cohort studies, conducted in the United States, Europe, and Asia, with smoking-adjusted analyses have provided little evidence for a positive, independent, dose–response relationship between coffee drinking and risk of bladder cancer,17–24 and two studies observed inverse associations among adults in Japan25 and women in the Netherlands.26 Nevertheless, positive associations in case–control studies,27 which are difficult to interpret due to potential recall and selection bias, as well as suggestive, but imprecise positive trends in some cohort studies20,22,24,26 have sustained speculation that coffee drinking causes bladder cancer. Importantly, the few cohort studies that have conducted sensitivity analyses among never smokers have had very limited case numbers in this subgroup, and risk estimates have been imprecise and inconsistent.20,25,26
Smoking is a strong risk factor for bladder cancer in both men and women in the NIH-AARP cohort with a relative risk estimate of 2.22 (95% CI = 2.03, 2.44) and 4.06 (95% CI = 3.66, 4.50) for former and current smokers as compared with never smokers, respectively.16 In the NIH-AARP cohort and in the overall US population,28 smoking is positively correlated with coffee drinking. One possible explanation is that the cytochrome P450 1A2 (CYP1A2) metabolic pathway is upregulated by both caffeine and compounds in tobacco smoke, including nicotine and polycyclic aromatic hydrocarbons29–31 such that the effect of caffeine is potentially weaker among heavier smokers leading them to drink more coffee than nonsmokers.
Because of this strong positive correlation between coffee drinking and tobacco smoking, associations between coffee drinking and bladder cancer are potentially confounded by smoking even if it is adjusted for in multivariable models. It is, therefore, important to evaluate whether the modest association we observed for coffee drinking and bladder cancer may be due to residual confounding by smoking. Thus, we investigated potential residual confounding using a multifaceted approach.32 First, we found that age- and sex-adjusted relative risks were substantially attenuated following detailed adjustment for baseline tobacco smoking. Our adjustment for smoking was probably imperfect owing to misclassification, which occurs because true tobacco smoke exposure is not accurately captured via self-report for reasons including but not limited to poor recall, differences in smoking practices (e.g., depth of inhalation or number of inhalations), changes during follow-up, or differences in second-hand smoke exposure. To quantify the impact of residual confounding resulting from confounder misclassification, we applied a method developed by Savitz and Barón.33 Assuming a true (i.e., the completely adjusted) relative risk of 1.00, we estimated that 69% of confounding was removed by adjustment for baseline tobacco smoking in our study. We treated the age- and sex-adjusted HR for a 1 cup/d increase in coffee intake (HR = 1.13) as the crude risk estimate and the HR with additional adjustment for smoking (HR = 1.04) as the partially adjusted risk estimate in our percent adjustment calculation ([(1.13 − 1.04)/(1.13 − 1.00) × 100] = 69%). A relatively small degree of tobacco smoking misclassification (i.e., approximately 91% sensitivity and 91% specificity), would result in the modest associations of coffee drinking with bladder cancer risk that we observed in NIH-AARP cohort.33 Strengthening our hypothesis for residual confounding, associations for coffee drinking and bladder cancer were larger in men than in women, and men were nearly twice as likely as women to be ever smokers. Moreover, among ever smokers, men were more than twice as likely as women to have smoked >40 cigarettes/d. Second, in smoking-stratified analyses, we showed that risk estimates for never smokers were close to or even below one, and that risk estimates for former smokers were weaker than those for current smokers, which is consistent with residual confounding by smoking. Finally, we found that additional adjustment for smoking history, among a subgroup of participants with detailed information on lifetime cigarette smoking use and intensity, resulted in further attenuation of risk estimates.
Limitations of this study include its observational design, the use of self-reported data on coffee intake, which is prone to measurement error potentially attenuating risk estimates, and a lack of data on some suspected bladder cancer risk factors. More specifically, the observed positive associations could reflect confounding by factors that were not measured in our cohort including aspects of smoking, like depth of inhalation and second-hand smoke exposure, or occupations in which workers were exposed to bladder cancer carcinogens such as aromatic amines (e.g., rubber workers and painters) or polycyclic aromatic hydrocarbons (e.g., metal, machine, and automobile workers).34 In fact, the relative prevalence of an unmeasured confounder would only have to be slightly higher among coffee drinkers than coffee nondrinkers to explain the observed association.35,36 Other possible explanations include chance or a real positive association.
In conclusion, our study, with more than 6,000 cases, is the largest prospective study of coffee and bladder cancer to date, and we observed little evidence for a strong association. We found a modest positive association between coffee drinking and bladder cancer; however, adjustment for cigarette smoking substantially attenuated risk estimates, and there was no evidence of an association among never smokers suggesting that residual confounding from imperfect measurement of smoking may explain our positive findings. Future pooled or meta-analyses should further explore the association of coffee drinking with bladder cancer among never smokers and confounding by occupational exposures, particularly among men. Coffee is widely consumed throughout the world; thus, it is vital that assessments of its carcinogenicity in humans are based on well-conducted, prospective cohort studies with sufficient case numbers and adequate control for important confounding factors like tobacco smoking.1
ACKNOWLEDGMENTS
We are indebted to the participants in the NIH-AARP Diet and Health Study for their outstanding cooperation. We also thank Sigurd Hermansen and Kerry Grace Morrissey from Westat for study outcomes ascertainment and management and Leslie Carroll at Information Management Services for data support and analysis.
REFERENCES
1. Loomis D, Guyton KZ, Grosse Y, et alInternational Agency for Research on Cancer Monograph Working Group. Carcinogenicity of drinking coffee, mate, and very hot beverages. Lancet Oncol. 2016;17:877–878.
2. Sang LX, Chang B, Li XH, Jiang MConsumption of coffee associated with reduced risk of liver cancer: a meta-analysis. BMC Gastroenterol. 2013;13:34.
3. Je Y, Giovannucci ECoffee consumption and risk of endometrial cancer: findings from a large up-to-date meta-analysis. Int J Cancer. 2012;131:1700–1710.
4. Je Y, Hankinson SE, Tworoger SS, De Vivo I, Giovannucci EA prospective cohort study of coffee consumption and risk of endometrial cancer over a 26-year follow-up. Cancer Epidemiol Biomarkers Prev. 2011;20:2487–2495.
5. Jiang W, Wu Y, Jiang XCoffee and caffeine intake and breast cancer risk: an updated dose-response meta-analysis of 37 published studies. Gynecol Oncol. 2013;129:620–629.
6. Lu Y, Zhai L, Zeng J, et alCoffee consumption and prostate cancer risk: an updated meta-analysis. Cancer Causes Control. 2014;25:591–604.
7. Tang N, Wu Y, Ma J, Wang B, Yu RCoffee consumption and risk of lung cancer: a meta-analysis. Lung Cancer. 2010;67:17–22.
8. Cancer IAfRo. Coffee, Tea, Mate, Methylxanthines and Methylglyoxal. 1991.Lyon;
9. Wu W, Tong Y, Zhao Q, Yu G, Wei X, Lu QCoffee consumption and bladder cancer: a meta-analysis of observational studies. Sci Rep. 2015;5:9051.
10. Schatzkin A, Subar AF, Thompson FE, et alDesign and serendipity in establishing a large cohort with wide dietary intake distributions: the National Institutes of Health-American Association of Retired Persons Diet and Health Study. Am J Epidemiol. 2001;154:1119–1125.
11. Michaud DS, Midthune D, Hermansen S, et alComparison of cancer registry case ascertainment with SEER estimates and self-reporting in a subset of the NIH-AARP Diet and Health Study. J Registry Manag. 2005; 32: 70–75.
12. Fritz AGInternational Classification of Diseases for Oncology, ICD-O. 2000.Geneva, Switzerland;
13. Petitti DB, Friedman GD, Kahn WAccuracy of information on smoking habits provided on self-administered research questionnaires. Am J Public Health. 1981;71:308–311.
14. Assaf AR, Parker D, Lapane KL, McKenney JL, Carleton RAAre there gender differences in self-reported smoking practices? Correlation with thiocyanate and cotinine levels in smokers and nonsmokers from the Pawtucket Heart Health Program. J Womens Health (Larchmt). 2002;11:899–906.
15. Subar AF, Midthune D, Kulldorff M, et alEvaluation of alternative approaches to assign nutrient values to food groups in food frequency questionnaires. Am J Epidemiol. 2000;152:279–286.
16. Freedman ND, Silverman DT, Hollenbeck AR, Schatzkin A, Abnet CCAssociation between smoking and risk of bladder cancer among men and women. JAMA. 2011;306:737–745.
17. Tripathi A, Folsom AR, Anderson KEIowa Women’s Health Study. Risk factors for urinary bladder carcinoma in postmenopausal women. The Iowa Women’s Health Study. Cancer. 2002;95:2316–2323.
18. Stensvold I, Jacobsen BKCoffee and cancer: a prospective study of 43,000 Norwegian men and women. Cancer Causes Control. 1994;5:401–408.
19. Nagano J, Kono S, Preston DL, et alBladder-cancer incidence in relation to vegetable and fruit consumption: a prospective study of atomic-bomb survivors. Int J Cancer. 2000;86:132–138.
20. Mills PK, Beeson WL, Phillips RL, Fraser GEBladder cancer in a low risk population: results from the Adventist Health Study. Am J Epidemiol. 1991;133:230–239.
21. Michaud DS, Spiegelman D, Clinton SK, et alFluid intake and the risk of bladder cancer in men. N Engl J Med. 1999;340:1390–1397.
22. Chyou PH, Nomura AM, Stemmermann GNA prospective study of diet, smoking, and lower urinary tract cancer. Ann Epidemiol. 1993;3:211–216.
23. Ros MM, Bas Bueno-de-Mesquita HB, Büchner FL, et alFluid intake and the risk of urothelial cell carcinomas in the European Prospective Investigation into Cancer and Nutrition (EPIC). Int J Cancer. 2011;128:2695–2708.
24. Kurahashi N, Inoue M, Iwasaki M, Sasazuki S, Tsugane SJapan Public Health Center (JPHC) Study Group. Coffee, green tea, and caffeine consumption and subsequent risk of bladder cancer in relation to smoking status: a prospective study in Japan. Cancer Sci. 2009;100:294–291.
25. Sugiyama K, Sugawara Y, Tomata Y, Nishino Y, Fukao A, Tsuji IThe association between coffee consumption and bladder cancer incidence in a pooled analysis of the Miyagi Cohort Study and Ohsaki Cohort Study. Eur J Cancer Prev. 2017;26:125–130.
26. Zeegers MP, Dorant E, Goldbohm RA, van den Brandt PAAre coffee, tea, and total fluid consumption associated with bladder cancer risk? Results from the Netherlands Cohort Study. Cancer Causes Control. 2001;12:231–238.
27. Zhou Y, Tian C, Jia CA dose-response meta-analysis of coffee consumption and bladder cancer. Prev Med. 2012;55:14–22.
28. Loftfield E, Freedman ND, Dodd KW, et alCoffee drinking is widespread in the united states, but usual intake varies by key demographic and lifestyle factors. J Nutr. 2016;146:1762–1768.
29. Landi MT, Sinha R, Lang NP, Kadlubar FFHuman cytochrome P4501A2. IARC Sci Publ. 1999: 173–95.
30. Kalow W, Tang BKUse of caffeine metabolite ratios to explore CYP1A2 and xanthine oxidase activities. Clin Pharmacol Ther. 1991;50(5 Pt 1):508–519.
31. Gunes A, Dahl MLVariation in CYP1A2 activity and its clinical implications: influence of environmental factors and genetic polymorphisms. Pharmacogenomics. 2008;9:625–637.
32. Guertin KA, Freedman ND, Loftfield E, Graubard BI, Caporaso NE, Sinha RCoffee consumption and incidence of lung cancer in the NIH-AARP Diet and Health Study. Int J Epidemiol. 2016;45:929–939.
33. Savitz DA, Barón AEEstimating and correcting for confounder misclassification. Am J Epidemiol. 1989;129:1062–1071.
34. Cumberbatch MG, Cox A, Teare D, Catto JWContemporary occupational carcinogen exposure and bladder cancer: a systematic review and meta-analysis. JAMA Oncol. 2015;1:1282–1290.
35. Cornfield J, Haenszel W, Hammond EC, Lilienfeld AM, Shimkin MB, Wynder ELSmoking and lung cancer: recent evidence and a discussion of some questions. 1959. Int J Epidemiol. 2009;38:1175–1191.
36. Greenhouse JBCommentary: cornfield, epidemiology and causality. Int J Epidemiol. 2009;38:1199–1201.