Observational Dose-Response Meta-Analysis Methods May Bias Risk Estimates at Low Consumption Levels: The Case of Meat and Colorectal Cancer

Abstract

Observational studies of foods and health are susceptible to bias, particularly from confounding between diet and other lifestyle factors. Common methods for deriving dose-response meta-analysis (DRMA) may contribute to biased or overly certain risk estimates. We used DRMA models to evaluate the empirical evidence for colorectal cancer (CRC) association with unprocessed red meat (RM) and processed meats (PM), and the consistency of this association for low and high consumers under different modeling assumptions. Using the Global Burden of Disease project’s systematic reviews as a start, we compiled a data set of studies of PM with 29 cohorts contributing 23,522,676 person-years and of 23 cohorts for RM totaling 17,259,839 person-years. We fitted DRMA models to lower consumers only [consumption < United States median of PM (21 g/d) or RM (56 g/d)] and compared them with DRMA models using all consumers. To investigate impacts of model selection, we compared classical DRMA models against an empirical model for both lower consumers only and for all consumers. Finally, we assessed if the type of reference consumer (nonconsumer or mixed consumer/nonconsumer) influenced a meta-analysis of the lowest consumption arm. We found no significant association with consumption of 50 g/d RM using an empirical fit with lower consumption (relative risk [RR] 0.93 (0.8–1.02) or all consumption levels (1.04 (0.99–1.10)), while classical models showed RRs as high as 1.09 (1.00–1.18) at 50g/day. PM consumption of 20 g/d was not associated with CRC (1.01 (0.87–1.18)) when using lower consumer data, regardless of model choice. Using all consumption data resulted in association with CRC at 20g/day of PM for the empirical models (1.07 (1.02–1.12)) and with as little as 1g/day for classical models. The empirical DRMA showed nonlinear, nonmonotonic relationships for PM and RM. Nonconsumer reference groups did not affect RM (P = 0.056) or PM (P = 0.937) association with CRC in lowest consumption arms. In conclusion, classical DRMA model assumptions and inclusion of higher consumption levels influence the association between CRC and low RM and PM consumption. Furthermore, a no-risk limit of 0 g/d consumption of RM and PM is inconsistent with the evidence.