- Research article
- Open Access
Protein, iron, and meat consumption and risk for rheumatoid arthritis: a prospective cohort study
Arthritis Research & Therapyvolume 9, Article number: R16 (2007)
A recent prospective study showed that higher consumption of red meat and total protein was associated with increased risk for inflammatory polyarthritis. We therefore prospectively examined the relationship between diet (in particular, protein, iron, and corresponding food sources) and incident rheumatoid arthritis (RA) among 82,063 women in the Nurses' Health Study. From 1980 to 2002, 546 incident cases of RA were confirmed by a connective tissue disease screening questionnaire and medical record review for American College of Rheumatology criteria for RA. Diet was assessed at baseline in 1980 and five additional times during follow up. We conducted Cox proportional hazards analyses to calculate the rate ratio of RA associated with intakes of protein (total, animal, and vegetable) and iron (total, dietary, from supplements, and heme iron) and their primary food sources, adjusting for age, smoking, body mass index, and reproductive factors. The multivariate models revealed no association between RA and any measure of protein or iron intake. In comparisons of highest with lowest quintiles of intake, the rate ratio for total protein was 1.17 (95% confidence interval 0.89–1.54; P for trend = 0.11) and for total iron it was 1.04 (95% confidence interval 0.77–1.41; P for trend = 0.82). Red meat, poultry, and fish were also not associated with RA risk. We were unable to confirm that there is an association between protein or meat and risk for RA in this large female cohort. Iron was also not associated with RA in this cohort.
Rheumatoid arthritis (RA) is associated with both genetic and environmental factors [1–7], but studies of dietary risk factors have been inconclusive . Studies of diet and risk for RA offer the potential to identify modifiable factors and so prevent RA in high-risk patients; they may also provide insights into disease pathogenesis.
Buchanan and Laurent  implicated diets high in protein in the etiology of RA. Furthermore, low-protein diets may improve RA symptoms [10–13]. In ecologic studies, the prevalence of RA is higher in countries with greater consumption of red meat . More recently, Pattison and colleagues  reported the first prospective investigation of red meat and risk for inflammatory polyarthritis (IP) and concluded that higher intakes of both red meat and protein increased the risk for IP, whereas iron – another nutrient component of meat – exhibited no association. The authors acknowledged that it remained unclear whether the observed associations were causative or whether meat consumption was a marker for other lifestyle factors.
To examine this issue further, we prospectively assessed risk for RA in relation to intakes of protein, iron, and meat in women in the Nurses' Health Study (NHS). We examined these intakes with further classifications into animal and vegetable protein; dietary, supplemental, and heme iron; and red meat, poultry, and fish.
Materials and methods
The NHS was established in 1976 when 121,700 female registered nurses (98% white), aged 30–55 years and residing in one of 11 US states, completed and returned the initial NHS mailed questionnaire on their medical history and lifestyle. Every 2 years, follow-up questionnaires have been sent to obtain up-to-date information on risk factors and to identify newly diagnosed diseases. Deaths are reported by family members or by the postal service in response to the follow-up questionnaires. In addition, we use the National Death Index to search for nonrespondents who might have died in the preceding interval. By comparing deaths ascertained from independent sources, we estimate that we have identified at least 98.2% of deaths occurring in the cohort .
The Partners HealthCare Institutional Review Board approved all aspects of this study, and all participants gave informed consent before they were entered into the study.
Ascertainment of rheumatoid arthritis cases
As previously described , self-reports of RA were confirmed using the Connective Tissue Disease Screening Questionnaire  and by medical record review for American College of Rheumatology (ACR) criteria for RA , conducted by two rheumatologists. We confirmed 807 cases of incident RA from 1976 to 2002.
For all analyses, we excluded the following: prevalent RA cases that were diagnosed before June 1980; RA cases with missing date of diagnosis; women who reported RA or connective tissue disease but in whom the diagnosis of RA was not confirmed by medical record review; nonresponders to the semiquantitative Food Frequency Questionnaire (FFQ) in 1980 (the baseline for this analysis); and participants with an unacceptable FFQ (<500 kcal/day or >3,500 kcal/day, accounting for approximately 4% of returned dietary questionnaires). Women were also censored during follow up when they failed to respond to any subsequent biennial questionnaire, because incident RA could not be identified in these cases. Thus, the final group studied included 82,063 women who were followed from 1980 until 2002 and 546 cases of incident RA who met the inclusion criteria, with a total of 1,668,894 person-years of follow up.
Assessment of dietary intake
Dietary intake was assessed in 1980, 1984, 1986, 1990, 1994, and 1998 using a semi-quantitative FFQ. In 1980, a total of 98,462 (81%) of the participants completed the FFQ and the completion rate has remained at about 80% during follow up. The initial FFQ contained 61 food items, but it has been expanded over the years such that 147 foods appeared on the 1998 questionnaire, including nine items for red meat (beef, pork, and lamb), four items for poultry (chicken and turkey), and four items for fish. For each food, participants reported their frequency of consumption of a specified serving size using nine frequency categories, ranging from never to six or more per day.
The validity and reproducibility of the FFQ for nutrients  and foods  have been documented elsewhere. Intakes calculated from the 1980 FFQ were found to be reasonably correlated with those from four 1-week diet records collected over 1 year among 173 NHS participants [20, 22]. The Pearson coefficients were 0.47 for total protein, 0.55 for total iron , and 0.38 for meat .
In this analysis, we examined associations between risk for RA and intakes of the following individual nutrients and components: total protein, animal protein, vegetable protein, total iron, dietary iron (from food sources), supplemental iron (from multivitamins and supplements), and heme iron (the iron with the highest bioavailability). We also examined meat, poultry, and fish (the primary food sources of protein and iron). At the 1998 dietary assessment in this cohort, 19% of protein came from red meat, 14% came from poultry, and 7% from fish. Heme iron also came primarily from the consumption of red meat (28%), poultry (24%), and fish (15%). Supplements contributed 25% of the total iron intake in this cohort.
Assessment of nondietary factors
Age, body mass index (weight [in kilograms] divided by height [in meters]2), and smoking status were updated every 2 years with information from the biennial questionnaires. Other factors were reported once: age at menarche in 1976, total months of breastfeeding for all children in 1986, and regularity of menses from age 20 to 35 years (very regular, usually regular, usually irregular, and very irregular) in 1982.
The number of person-years of follow up was ascertained based on the interval between the date of return of the 1980 questionnaire and the date of diagnosis of RA (as defined in the medical record), death, the end of the study period (1 June 2002), or loss to follow up (defined as no further return of questionnaires) for each participant.
Nutrient and food intakes were categorized into quintiles, and incidence rates for RA were calculated by dividing the number of incident cases by the number of person-years in each quintile of dietary exposure. Rate ratio (RRs) were calculated by dividing the incidence rates in the higher quintiles by the corresponding rate in the reference (lowest) quintile. Age-adjusted and multivariate RRs were estimated using Cox proportional hazards models adjusting for age (continuous variable) and other potential counfounders. We controlled for the following variables because they have either been shown to be associated with RA or were found in this study to be potential confounders: body mass index (categorized as <22, 22 to 24.9, 25 to 29.9, 30 to 34.9, and ≥35 kg/m2), smoking status (never, past, or present), and total lifetime breastfeeding history (nulliparous, parous, and breastfeeding for 0, 1 to 11, ≥12 total months). In addition, we controlled for total energy to reduce measurement error due to general over-reporting or under-reporting of food items . Age at menarche and regularity of menses were not retained as covariates. For all RRs, we calculated the 95% confidence interval (CI). All P values were two-tailed, and P < 0.05 was considered to be statistically significant. Tests for trend were conducted by assigning the median value for each quintile of nutrient and food intake, modeling this variable as a continuous variable.
Nutrient intakes were energy-adjusted using the multivariate residual method . In order to represent the long-term dietary patterns of individual women, our primary analysis used cumulative average food and nutrient intakes from all available dietary questionnaires up to the start of each 2-year interval . For example, the 1980 diet was related to RA incidence during the period from 1980 to 1984; the average of the 1980 and 1984 diets was related to RA incidence during the period from 1984 to 1986; the average of the 1980, 1984, and 1986 diets was related to the RA incidence during the period between 1986 and 1990, and so on, through to 2002.
Age standardized characteristics of the study population in 1990 according to intakes of total protein and heme iron are shown in Table 1. The 1990 time point was chosen because it represents the approximate mid-point of follow up. Body mass index was higher among women in the highest consumption categories of total protein and heme iron. Women with the lowest protein and highest heme iron consumptions were more likely to smoke and, if parous, they were less likely to have breastfed for a total of 12 months or more. Higher total protein intakes were associated with higher heme iron intakes.
In the age-adjusted model, higher total protein intake was associated with greater risk for RA (quintile 5 [89.0 g/day] versus quintile 1 [60.8 g/day]: RR 1.23, 95% CI 0.94–1.61; P for trend = 0.04), but this association was attenuated and the test for trend was no longer significant in the multivariate model (RR 1.17, 95% CI 0.89–1.55; P for trend = 0.12; Table 2). Neither the animal nor vegetable component of protein exhibited any relation to risk for RA. We also did not observe any association with total iron intake (RR 1.00, 95% CI 0.74–1.36 for the highest versus lowest quintile) or with its components of dietary iron, supplemental iron, and heme iron.
No significant associations were observed between the incidence of RA and consumption of total meat, red meat, poultry, or fish (Table 3). For total meat, which included red meat and poultry, the multivariate RR was 0.91 (95% CI 0.67–1.23) in the highest (2.54 servings/day) versus lowest (0.82 servings/day) quintiles of intake. More detailed analyses of individual foods that contribute to each of these major food groups also exhibited no association with RA.
To avoid confounding by indication (for example, dietary changes occurring after RA symptom onset), we also performed analyses in which dietary variables were updated only until the date of first symptom of RA, rather than until the date of RA diagnosis. We also performed lagged analyses such that the dietary intakes associated with RA cases were assessed at least 4 years before the date of diagnosis. In order to account for possible influence of recent dietary intake, we also examined our exposures based on the most recent dietary measures, rather than using long-term average intakes. The results revealed no associations with the nutrient or food exposures.
In this large prospective cohort study involving women, we observed no significant association between protein or iron intakes and risk for RA, including specific analyses of animal and vegetable protein, heme iron, and iron from foods and from supplements. Furthermore, no associations were observed between the primary food sources of these nutrients, namely red meat, poultry, and fish.
Our results differ from those of a nested case-control study  that reported increased risk of IP with greater consumption of protein and red meat. Pattison and coworkers  studied dietary intake and risk for IP between 1993 and 2002, within a prospective population-based study of cancer incidence in Norfolk, England (European Prospective Investigation of Cancer Incidence [EPIC]). In their study they compared 88 patients with IP, identified by linkage with the Norfolk Arthritis Register (a primary care-based inception study of IP), with 167 age-matched and sex-matched control individuals from EPIC who had remained free from IP during the follow-up period. Although the study did not analyze subtypes of protein, animal and vegetable protein, it did analyze the food sources that contribute to each of these categories. The investigators reported an increased risk for IP with greater protein comsumption (>75.3 g/day versus <62.4 g/day: adjusted odds ratio [OR] 2.9, 95% CI 1.1–7.5) and no association with iron. In contrast to our findings, the study by Pattison and coworkers indicated that individuals with the highest level of consumption of red meat (>58.0 g/day versus <25.5 g/day: adjusted OR 1.9, 95% CI 0.9–4.0) and red meat combined with meat products (for instance, sausage and ham; >87.8 g/day versus <49.0 g/day: adjusted OR 2.3, 95% CI 1.1–4.9) were at increased risk for IP.
The discrepancy between the findings of that study and ours could be attributed to methodologic differences. First, the EPIC study assessed dietary intake once, using a 7-day food diary, whereas we used semiquantitative FFQ assessed repeatedly. The FFQ consists of two components : a food list and a frequency response section for individuals to report how often each food was eaten over the previous year. The 7-day food diary consists of a detailed listing of all foods consumed by an individual on 1 day or more . Food intake is recorded by the individual at the time when the foods are eaten, which has the advantages of relying less on memory and permitting direct assessment of portion sizes. In comparison, the FFQ suffers the disadvantages of restrictions imposed by a fixed list of foods, memory, perception of portion sizes, and interpretation of questions. Dietary records provide more precise quantification of foods consumed, but they only reflect short-term diet, because only a limited number of days of diet records are used. Results of validation studies demonstrate greater correlation of blood levels of certain nutrients with 7-day diet diaries than with FFQ findings .
However, our objective was to assess long-term dietary exposures. Therefore, we cumulatively averaged and updated dietary intake assessed six different times over the 22-year period of follow up, which is known to reduce random error in long-term dietary measurement, rather than relying upon one assessment at baseline. Furthermore, results of analyses of more recent diet were consistent with analyses of cumulative diet. Even if absolute measures are not precise, the FFQ is able to rank respondents into higher and lower categories of intake. We energy-adjusted nutrient intakes in order to account for differences due to under-reporting or over-reporting on the FFQ.
Bingham and coworkers  demonstrated a strong association between diet and cancer using 7-day diaries but a modest relationship when the FFQ was used, and they suggested that this pattern might also be seen in other studies analyzing the association of diet and chronic diseases. However, previous studies undertaken in the Nurses' Health Study cohort and others that used the FFQ demonstrated associations between meat and protein and breast cancer, colorectal cancer, lymphoma, coronary heart disease, diabetes, and gout [29–35].
Finally, it is possible that dietary protein intake differs between the USA and the UK. However, comparisons of the median intake of total protein and total iron in the quintiles used in the present study (Table 2) with the tertiles of intake in the EPIC study  demonstrate that the range and categories of intake in the two studies were similar.
A second difference between our study and the EPIC study was that we identified individuals with RA rigorously using the ACR criteria, in which at least four out of seven criteria had to be satisfied in order for a participant to be considered a case. In contrast, the outcome considered by Pattison and colleagues  was the presence of IP, which is defined as inflammation affecting two or more peripheral joints and persisting for 4 weeks or longer. Within 5 years, 60% of IP patients satisfy ACR criteria for RA .
Third, discrepancies between our study and the EPIC study might be related to differences in sex, because our study included women only whereas the EPIC study  included men and women. It is also possible that the discrepant findings resulted from socioeconomic differences; well educated nurses were enrolled in our study, whereas the EPIC cohort included diverse population-based cases and controls.
Strengths of our study include the large number of incident cases of RA, the repeated prospective assessment of exposures, and the lengthy follow-up period. The validation of self-reported RA through medical record review rather than by physical examination is a potential weakness of the study. However, 82% of the RA cases were diagnosed by ACR members, which adds support to the validity of the diagnoses. There is potential for misclassification of RA cases as noncases when diagnosis relies solely upon medical record documentation. Therefore, those women who self-reported RA or other connective tissue diseases in whom the diagnosis of RA was not confirmed by medical record review were excluded from the analyses. It is possible that the null results from this study are due to unmeasured confounding (for example, socioeconomic status), although there are no strong risk factors for RA that could account for attenuation of a true association. Finally, although the participants in the present do not represent a random sample of women living in the USA, it is unlikely that the biologic relationships among these women will differ from those among women in general.
No clear associations were observed between dietary protein, iron, or meat, including red meat, and risk for RA in this large prospective cohort of women.
American College of Rheumatology
Food Frequency Questionnaire
Nepom GT, Byers P, Seyfried C, Healey LA, Wilske KR, Stage D, Nepom BS: HLA genes associated with rheumatoid arthritis. Identification of susceptibility alleles using specific oligonucleotide probes. Arthritis Rheum. 1989, 32: 15-21. 10.1002/anr.1780320104.
Wordsworth BP, Lanchbury JS, Sakkas LI, Welsh KI, Panayi GS, Bell JI: HLA-DR4 subtype frequencies in rheumatoid arthritis indicate that DRB1 is the major susceptibility locus within the HLA class II region. Proc Natl Acad Sci USA. 1989, 86: 10049-10053. 10.1073/pnas.86.24.10049.
Ronningen KS, Spurkland A, Egeland T, Iwe T, Munthe E, Vartdal F, Thorsby E: Rheumatoid arthritis may be primarily associated with HLA-DR4 molecules sharing a particular sequence at residues 67–74. Tissue Antigens. 1990, 36: 235-240.
Symmons DP, Bankhead CR, Harrison BJ, Brennan P, Barrett EM, Scott DG, Silman AJ: Blood transfusion, smoking, and obesity as risk factors for the development of rheumatoid arthritis: results from a primary care-based incident case-control study in Norfolk, England. Arthritis Rheum. 1997, 40: 1955-1961. 10.1002/art.1780401106.
Stolt P, Bengtsson C, Nordmark B, Lindblad S, Lundberg I, Klareskog L, Alfredsson L: Quantification of the influence of cigarette smoking on rheumatoid arthritis: results from a population based case-control study, using incident cases. Ann Rheum Dis. 2003, 62: 835-841. 10.1136/ard.62.9.835.
Karlson EW, Lee IM, Cook NR, Manson JE, Buring JE, Hennekens CH: A retrospective cohort study of cigarette smoking and risk of rheumatoid arthritis in female health professionals. Arthritis Rheum. 1999, 42: 910-917. 10.1002/1529-0131(199905)42:5<910::AID-ANR9>3.0.CO;2-D.
Criswell LA, Merlino LA, Cerhan JR, Mikuls TR, Mudano AS, Burma M, Folsom AR, Saag KG: Cigarette smoking and the risk of rheumatoid arthritis among postmenopausal women: results from the Iowa Women's Health Study. Am J Med. 2002, 112: 465-471. 10.1016/S0002-9343(02)01051-3.
Choi HK: Dietary risk factors for rheumatic diseases. Curr Opin Rheumatol. 2005, 17: 141-146. 10.1097/01.cco.0000152664.87204.3c.
Buchanan WW, Laurent RM: Rheumatoid arthritis: an example of ecological succession?. Can Bull Med Hist. 1990, 7: 77-91.
Skoldstam L: Fasting and vegan diet in rheumatoid arthritis. Scand J Rheumatol. 1986, 15: 219-221.
Kjeldsen-Kragh J, Haugen M, Borchgrevink CF, Laerum E, Eek M, Mowinkel P, Hovi K, Forre O: Controlled trial of fasting and one-year vegetarian diet in rheumatoid arthritis. Lancet. 1991, 338: 899-902. 10.1016/0140-6736(91)91770-U.
Kjeldsen-Kragh J, Haugen M, Forre O, Laache H, Malt UF: Vegetarian diet for patients with rheumatoid arthritis: can the clinical effects be explained by the psychological characteristics of the patients?. Br J Rheumatol. 1994, 33: 569-575. 10.1093/rheumatology/33.6.569.
Nenonen MT, Helve TA, Rauma AL, Hanninen OO: Uncooked, lactobacilli-rich, vegan food and rheumatoid arthritis. Br J Rheumatol. 1998, 37: 274-281. 10.1093/rheumatology/37.3.274.
Grant WB: The role of meat in the expression of rheumatoid arthritis. Br J Nutr. 2000, 84: 589-595.
Pattison DJ, Symmons DP, Lunt M, Welch A, Luben R, Bingham SA, Khaw KT, Day NE, Silman AJ: Dietary risk factors for the development of inflammatory polyarthritis: evidence for a role of high level of red meat consumption. Arthritis Rheum. 2004, 50: 3804-3812. 10.1002/art.20731.
Rich-Edwards JW, Corsano KA, Stampfer MJ: Test of the National Death Index and Equifax Nationwide Death Search. Am J Epidemiol. 1994, 140: 1016-1019.
Costenbader KH, Feskanich D, Mandl LA, Karlson EW: Smoking intensity, duration, and cessation, and the risk of rheumatoid arthritis in women. Am J Med. 2006, 119: e1-e9. 10.1016/j.amjmed.2005.09.053.
Karlson EW, Sanchez-Guerrero J, Wright EA, Lew RA, Daltroy LH, Katz JN, Liang MH: A connective tissue disease screening questionnaire for population studies. Ann Epidemiol. 1995, 5: 297-302. 10.1016/1047-2797(94)00096-C.
Arnett FC, Edworthy SM, Bloch DA, McShane DJ, Fries JF, Cooper NS, Healey LA, Kaplan SR, Liang MH, Luthra HS, et al: The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988, 31: 315-324. 10.1002/art.1780310302.
Willett WC, Sampson L, Stampfer MJ, Rosner B, Bain C, Witschi J, Hennekens CH, Speizer FE: Reproducibility and validity of a semiquantitative food frequency questionnaire. Am J Epidemiol. 1985, 122: 51-65.
Salvini S, Hunter DJ, Sampson L, Stampfer MJ, Colditz GA, Rosner B, Willett WC: Food-based validation of a dietary questionnaire: the effects of week-to-week variation in food consumption. Int J Epidemiol. 1989, 18: 858-867. 10.1093/ije/18.4.858.
Willett WC, Sampson L, Browne ML, Stampfer MJ, Rosner B, Hennekens CH, Speizer FE: The use of a self-administered questionnaire to assess diet four years in the past. Am J Epidemiol. 1988, 127: 188-199.
Willett W, Stampfer M: Implications of total energy intake for epidemiologic analyses. Nutritional epidemiology. Edited by: Willett W. 1998, Oxford: Oxford University Press, 273-301. 2
Hu FB, Stampfer MJ, Manson JE, Rimm E, Colditz GA, Rosner BA, Hennekens CH, Willett WC: Dietary fat intake and the risk of coronary heart disease in women. N Engl J Med. 1997, 337: 1491-1499. 10.1056/NEJM199711203372102.
Willett W: Food Frequency Methods. Nutritional epidemiology. 2
Bingham SA, Gill C, Welch A, Cassidy A, Runswick SA, Oakes S, Lubin R, Thurnham DI, Key TJ, Roe L, et al: Validation of dietary assessment methods in the UK arm of EPIC using weighed records, and 24-hour urinary nitrogen and potassium and serum vitamin C and carotenoids as biomarkers. Int J Epidemiol. 1997, 26 Suppl 1 (): S137-S151. 10.1093/ije/26.suppl_1.S137.
Day N, McKeown N, Wong M, Welch A, Bingham S: Epidemiological assessment of diet: a comparison of a 7-day diary with a food frequency questionnaire using urinary markers of nitrogen, potassium and sodium. Int J Epidemiol. 2001, 30: 309-317. 10.1093/ije/30.2.309.
Bingham SA, Luben R, Welch A, Wareham N, Khaw KT, Day N: Are imprecise methods obscuring a relation between fat and breast cancer?. Lancet. 2003, 362: 212-214. 10.1016/S0140-6736(03)13913-X.
Cho E, Chen WY, Hunter DJ, Stampfer MJ, Colditz GA, Hankinson SE, Willett WC: Red meat intake and risk of breast cancer among premenopausal women. Arch Intern Med. 2006, 166: 2253-2259. 10.1001/archinte.166.20.2253.
Zhang S, Hunter DJ, Rosner BA, Colditz GA, Fuchs CS, Speizer FE, Willett WC: Dietary fat and protein in relation to risk of non-Hodgkin's lymphoma among women. J Natl Cancer Inst. 1999, 91: 1751-1758. 10.1093/jnci/91.20.1751.
Giovannucci E, Rimm EB, Stampfer MJ, Colditz GA, Ascherio A, Willett WC: Intake of fat, meat, and fiber in relation to risk of colon cancer in men. Cancer Res. 1994, 54: 2390-2397.
Choi HK, Atkinson K, Karlson EW, Willett W, Curhan G: Purine-rich foods, dairy and protein intake, and the risk of gout in men. N Engl J Med. 2004, 350: 1093-1103. 10.1056/NEJMoa035700.
van Dam RM, Willett WC, Rimm EB, Stampfer MJ, Hu FB: Dietary fat and meat intake in relation to risk of type 2 diabetes in men. Diabetes Care. 2002, 25: 417-424. 10.2337/diacare.25.3.417.
Zhang C, Schulze MB, Solomon CG, Hu FB: A prospective study of dietary patterns, meat intake and the risk of gestational diabetes mellitus. Diabetologia. 2006, 49: 2604-2613. 10.1007/s00125-006-0422-1.
Hu FB, Bronner L, Willett WC, Stampfer MJ, Rexrode KM, Albert CM, Hunter D, Manson JE: Fish and omega-3 fatty acid intake and risk of coronary heart disease in women. JAMA. 2002, 287: 1815-1821. 10.1001/jama.287.14.1815.
Wiles N, Symmons DP, Harrison B, Barrett E, Barrett JH, Scott DG, Silman AJ: Estimating the incidence of rheumatoid arthritis: trying to hit a moving target?. Arthritis Rheum. 1999, 42: 1339-1346. 10.1002/1529-0131(199907)42:7<1339::AID-ANR6>3.0.CO;2-Y.
We would like to acknowledge all of the nurses who participate in this study and also Gideon Aweh, for programming assistance. We would also like to thank the support for this research by grants CA87979, R01 AR42630, P60 AR47782 and R0149880 from the National Institutes of Health.
The authors declare that they have no competing interests.
EB-G contributed to the study concept and design, data collection and analyses, and manuscript writing and editing. DF contributed to the data analyses and statistical support, as well as to the manuscript writing and editing. FB contributed to the concept and design, and to the manuscript writing and editing. LAM contributed to data collection and manuscript writing and editing. EWK obtained the funding and contributed to the concept and design, data collection and analyses, and manuscript writing and editing.