The principal complication of bone loss is fragility fracture. Approximately 44 million US adults aged 50 years and older (80 % women) have low bone mass that puts them at risk for fragility fractures( 1 ). In 2005, more than 2 million new fragility fractures occurred in the USA, with direct health care costs totalling more than $17 billion( Reference Ray, Chan and Thamer 2 ). Roughly 50 % of those with hip fractures will never walk again without assistance and 25 % will require long-term care( Reference Riggs and Melton 3 ).
Protein-rich foods are often associated with higher levels of satiety than carbohydrate-rich foods( Reference Halton and Hu 4 ). Owing to this beneficial effect, diets that include protein-rich foods are often advocated for weight loss and weight maintenance( Reference Halton and Hu 4 ). However, dietary protein has unclear effects on bone health( Reference Ginty 5 ). Short-term intervention studies using purified protein supplements have shown that 1 mg Ca is, on average, lost in the urine for every 1 g increase in protein intake( Reference Kerstetter and Allen 6 ). However, high-protein diets appear to increase intestinal Ca absorption, while low-protein diets decrease it( Reference Hunt, Johnson and Fariba Roughead 7 ).
Observational epidemiological studies are inconsistent on the subject of dietary protein and the risk of fractures. In ecological studies of animal protein intake and rates of hip fracture, the highest rates were observed in Western countries where dietary protein intake, particularly from animal sources, is also highest( Reference Abelow, Holford and Insogna 8 , Reference Frassetto, Todd and Morris 9 ). Consumption of animal proteins, high in acidifying amino acids such as cysteine and methionine, could potentially increase the risk of bone loss( Reference Thorpe, Mojtahedi and Chapman-Novakofski 10 , Reference Sellmeyer, Stone and Sebastian 11 ). A few large prospective cohort studies have also associated higher protein intake (particularly from animal sources) with higher incidence of hip fracture( Reference Feskanich, Willett and Stampfer 12 – Reference Meyer, Pedersen and Loken 14 ), while others have found the opposite association( Reference Munger, Cerhan and Chiu 15 – Reference Promislow, Goodman-Gruen and Slymen 17 ), and one found no association( Reference Key, Appleby and Spencer 18 ). Aside from differences in measurement methods, this inconsistency may be due in part to racial/ethnic differences. Recent studies have suggested that the association between protein intake and bone mineral density (BMD) may differ by race and/or ethnicity( Reference Henderson, Price and Cole 19 – Reference Chiu, Lan and Yang 24 ), which could reflect differences in dietary patterns and sources of protein.
Diagnosis of osteoporosis and clinical guidelines for its management are often based on BMD of the lumbar spine, total hip or femoral neck( Reference Schousboe, Shepherd and Bilezikian 25 , Reference Papaioannou, Morin and Cheung 26 ). Measures of volumetric or three-dimensional BMD can be obtained from computed tomography (CT) scans, while traditional dual-energy X-ray absorptiometry scans provide two-dimensional areal BMD. Unlike dual-energy X-ray absorptiometry scans, quantitative CT scans can also provide separate measures for cortical and trabecular compartments of bone, aspects of geometry and architecture that may be related to fracture risk( Reference Schousboe, Shepherd and Bilezikian 25 , Reference Cauley, Blackwell and Zmuda 27 , Reference Roux, Wegrzyn and Arlot 28 ).
In the present study, we examined the association between dietary protein and volumetric bone mineral density (vBMD) of the lumbar spine obtained from quantitative CT, specifically the average of lumbar vBMD in the second through the fourth lumbar vertebrae (primary outcome) and vBMD Z-score generated by subtracting the mean and dividing the difference by the standard deviation for each observation based on sex- and race/ethnicity-specific distributions (secondary outcome) among participants of the Multi-Ethnic Study of Atherosclerosis (MESA). We also evaluated whether relationships differed by source of protein intake (animal or vegetable), sex and race/ethnicity (non-Hispanic white, Chinese-American, African-American and Hispanic-American).
Methods
Participants
The methods of the MESA study have been described in detail elsewhere( Reference Bild, Bluemke and Burke 29 ). In brief, the MESA cohort was initially examined between July 2000 and August 2002 (baseline). Participants were recruited from six US communities. The study population consisted of 6814 men and women aged 45–84 years who were free of CVD at their baseline examination. Approximately 38 % of the cohort was non-Hispanic white, 11 % Chinese-American, 28 % African-American and 23 % Hispanic-American. We used data from a subsample of MESA participants who participated in the MESA Abdominal Aortic Calcium Study. Between August 2002 and September 2005, participants were recruited to the MESA Abdominal Aortic Calcium Study during their follow-up visits in the following five MESA communities: Chicago, IL; Forsyth County, NC; Los Angeles County, CA; New York, NY; St Paul, MN. In total, 1967 men and women who participated in the MESA Abdominal Aortic Calcium Study completed CT scanning. Of those, data from fourteen participants were excluded because they lacked data on bone density of the second through fourth lumbar vertebrae; 235 participants who reported consuming >25 104 kJ/d or < 2510 kJ/d were excluded; and an additional sixty were excluded because they provided implausible data on physical activity (>24 h/d). After these exclusions, data from 1658 participants were included in the present analysis. The MESA study was conducted according to the guidelines of the Declaration of Helsinki, and all MESA procedures were approved by the institutional review boards of participating institutions. Written informed consent was obtained from all MESA participants during primary data collection. This secondary data analysis was approved by the institutional review board of Tulane University.
Bone density measurement
The measurement of vBMD using CT scans of the abdomen in MESA has been reported in a previous study( Reference Hyder, Allison and Wong 30 ). Scanning was performed using an electron-beam CT scanner (Chicago, IL; New York, NY; and Los Angeles, CA, USA; Imatron C-150; General Electric Medical Systems) or a multi-detector CT system that utilised helical scanning with reconstruction in 5 mm thick cuts and 350 mm field of view (New York, NY; Forsyth County, NC; and St Paul, MN field centres, USA; Siemens, Inc., General Electric Medical system). During the scanning, CT parameters were directly converted to equivalent vBMD in mg/cm3( Reference Cann 31 ). Using the Image Analysis QCT 3D PLUS software program (Image Analysis), vBMD was measured in a virtual 10 mm-thick slice from each vertebral bone (L2–L4) using software-directed, automated placement of the region of interest in the anterior one-half to one-third of the vertebral body, where it encompassed a large area of trabecular or cancellous bone, excluded cortical bone, and excluded the basivertebral plexus. A trained reader adjusted placement of the reading region and excluded vertebrae if fractures, metastatic lesions, osteophytes or benign focal lesions were detected. All scans were brightness adjusted with a standard phantom. The scans were read by a blinded, experienced reader at the MESA CT reading centre. Previous studies have demonstrated the comparability and reproducibility of image data collected using these scanners( Reference Carr, Nelson and Wong 32 , Reference Detrano, Anderson and Nelson 33 ).
Dietary assessment
All participants were asked to complete a 120-item FFQ at baseline, on average 3 years before the CT scan, to assess their usual intake of specific foods and beverages over the past year. The FFQ utilised the Block format and was patterned after the FFQ used in the Insulin Resistance Atherosclerosis Study, which has been validated among non-Hispanic white, African-American and Hispanic-American, and modified to include Chinese foods( Reference Nettleton, Rock and Wang 34 , Reference Block, Woods and Potosky 35 ). For each food item, the consumption frequency (times/d, week or month) and serving size (small, medium or large) were recorded. Frequency options included nine responses ranging from ‘rare or never’ to ‘ ≥ 2 times/d’ for each food item( Reference Block, Woods and Potosky 35 ). Related line items were combined to form forty-seven different food groups( Reference Block, Woods and Potosky 35 ). Daily macronutrient and micronutrient intakes were estimated by multiplying frequency and serving size (age- and sex-specific and portion-size gram weights) for each food/beverage consumed by the nutrient content of that food or beverage (Nutrition Data Systems for Research; University of Minnesota; Minneapolis). Dietary protein intake, including total, animal and vegetable protein, was calculated in g and also as a percentage of total energy intake (% energy). In studies examining the reproducibility of a single FFQ, correlations ranged between 0·5 and 0·7 for nutrient intakes measured at intervals of 1 to 10 years( Reference Willett 36 ).
Additional variables
Baseline data collection also included a medical history, standardised medical examination, laboratory tests and anthropometric measurements. Age, sex, race/ethnicity, age at menopause, educational level, current cigarette smoking and current use of hormone replacement therapy were collected using interviewer and self-administered questionnaires. We calculated BMI at baseline as weight (kg) divided by height (m) squared. Physical activity was assessed as metabolic equivalents (MET), using a semi-quantitative questionnaire adapted from the Cross-Cultural Activity Participation Study( Reference Lutsey, Jacobs and Kori 37 ). Leisure-time physical activity was calculated as the sum of min/week of moderate to vigorous activities (walking, sports, dance and conditioning) multiplied by the activity's individual MET value. Likewise, a sedentariness score was calculated as the sum of MET-min/week of sitting, reclining, knitting, sewing, reading, driving a car and watching television. Estimated glomerular filtration rate was calculated using the CKD-EPI (Chronic Kidney Disease Epidemiology Collaboration) equation( Reference Levey, Stevens and Schmid 38 ).
Statistical analysis
For each baseline characteristic, the mean value or corresponding percent of study participants was calculated by quartile of total protein intake as a percentage of total energy. The statistical significance of differences among quartiles was examined using χ2 tests for categorical variables and generalised linear regression analyses for continuous variables.
There was evidence that the association between dietary protein intake and BMD differed by race/ethnicity and sex (P for interaction = 0·03 for protein (as % energy) and P for interaction = 0·02 for protein (in g)). Thus, we categorised participants into groups based on sex- and race/ethnicity-specific quartiles of dietary protein intake for analyses.
We used substitution modelling to investigate the association of ‘substituting’ an energy equivalent amount of dietary protein (total, animal and/or vegetable protein) for fat with lumbar spine vBMD by including dietary protein and dietary carbohydrate intakes in the same multivariable model with total energy intake. Tests for linear trend were conducted for adjusted means of vBMD across the quartiles of energy intake from dietary protein (for total, animal and vegetable protein). Adjusted mean vBMD was computed by quartile of sex- and race/ethnicity-specific protein intake using linear regression models. In multivariable-adjusted models, a variety of potential confounders including age (years), BMI (kg/m2), physical activity (MET-min/week), sedentariness score (MET-min/week), current smoking (yes/no), educational level ( ≤ high school v. ≥ college), hormone replacement therapy use, age at menopause (years), and intakes of total energy (kJ/d), dietary carbohydrate (% energy), Ca (mg), P (mg), Mg (mg) and alcohol (g/d) were included. We examined the influence of renal function by including estimated glomerular filtration rate, and generated a propensity score to represent Ca, Mg and P intakes as a single covariate. In order to limit the influence of heteroscedasticity or potential outliers in the present analysis, robust regression models were employed to examine the impact of quartile of dietary protein intake (as % energy), one standard deviation changes in protein intake (as % energy), and absolute protein intake (in g) on vBMD Z-score among the different racial/ethnic groups. We tested for interactions between dietary protein intake and hormone replacement therapy as well as Ca intake. All P values were two-sided, and statistical significance was defined as P< 0·05. Multiple comparisons were controlled by using the false discovery rate method( Reference Benjamini and Hochberg 39 ). We used SAS (version 9.2; SAS Institute, Inc.) for all analyses.
Results
The present analysis included 1658 participants, of whom 801 or 48 % were women. Of these women, 317 (39·6 %) were non-Hispanic white, 104 (13·0 %) were Chinese-American, 173 (21·6 %) were African-American and 207 (25·8 %) were Hispanic-American. Of men, 375 (43·8 %), 126 (14·7 %), 133 (15·5 %) and 223 (26·0 %) were non-Hispanic white, Chinese-American, African-American and Hispanic-American, respectively. The mean age among women was 63 years and that among men was 62 years. Baseline characteristics are presented by quartile of dietary intake of total protein in Table 1. Those in the highest quartile of protein intake were slightly younger, more often female and less likely to smoke cigarettes. Non-Hispanic whites and African-Americans tended to consume less protein, while roughly half of Chinese were in the highest quartile of protein intake. The amount of vegetable protein and fat intakes (as % energy) was consistent across the quartiles, but intake of animal protein was higher and carbohydrate was lower among those with higher total protein intake.
MET, metabolic equivalents; vBMD, volumetric bone mineral density; HRT, hormone replacement therapy.
* The ranges for the quartiles 1–4 of dietary total protein intake were 6·1–13·6, 13·7–15·7, 15·8–17·9 and 18·0–33·5, respectively.
† Cold cereal is a part of whole grain.
‡ Trabecular vBMD of Lumbar L2–4 vertebrae.
§ Female only.
We identified a significant interaction by which race/ethnicity modified the effect of dietary protein on vBMD; therefore, race/ethnicity-specific quartiles of protein intake (% energy) were used to assess the association between vBMD and protein intake in men and women separately. These quartiles are presented in Table 2.
Q, quartile; TP, total protein; AP, animal protein; VP, vegetable protein.
No significant relationships were observed among men of any race/ethnicity (Table 3). For women, there was no significant interaction between dietary protein and hormone replacement therapy use or Ca intake, in association with lumbar spine vBMD. Non-Hispanic white women with higher vegetable protein intake had a higher mean vBMD (P for trend = 0·03), after adjustment for age, BMI, physical activity, sedentariness score, current smoking, educational level, hormone therapy use, age at menopause, and intakes of total energy, dietary carbohydrate, animal protein, Ca, P, Mg and alcohol (Table 3). The associations with regard to total and animal protein intakes were not significant among non-Hispanic white women.
Q, quartile; TP, total protein; AP, animal protein; VP, vegetable protein.
* Upper and lower limits of each race/ethnicity- and sex-specific quartile are shown in Table 2.
† From the trend test for the median of bone mineral density in each quartile.
‡ Adjusted for age, BMI, physical activity, sedentariness score, smoking, educational level, hormone therapy use (if applicable), age at menopause (if applicable), and intakes of total energy, dietary carbohydrate as a percentage of energy, Ca, P, Mg and alcohol.
§ TP model additionally adjusted for vegetable protein intake as a percentage of energy.
∥ TP model additionally adjusted for animal protein intake as a percentage of energy.
We further examined whether the associations observed among women could be due to outliers or heterogeneity by employing robust regression methods. These results were consistent with the adjusted robust regression analyses examining the association between protein intake and vBMD Z-score (Table 4). Compared with the lowest quartile of vegetable protein intake, non-Hispanic white women in the highest quartile had significantly greater lumbar spine vBMD Z-scores after adjusting for potential confounders. Moreover, this result was also consistent when continuous dietary vegetable protein intake was used in robust regression models; a greater intake of vegetable protein was significantly associated with a higher vBMD (for each 1 % greater protein intake lumbar spine vBMD Z-score was 3·87 % greater, P= 0·01) among non-Hispanic white women. In contrast, there were no consistently significant associations between dietary protein (total, animal or vegetable protein) intake and vBMD in other sex and racial/ethnic groups. For example, Hispanic-American women with higher animal protein intake had a higher mean vBMD (P for trend = 0·04), and Chinese-American women with higher animal protein intake had a lower mean vBMD (P for trend = 0·02) (Table 3). These associations were no longer significant in robust regression analyses, regardless of whether the categorical (Table 4) or continuous (data not shown) dietary intake of animal protein was used in the models. Additional models without substitution of fat for protein or the exchange of protein for carbohydrate did not alter the results (data not shown). Modelling absolute protein intakes as one standard deviation difference (in g) was not associated with lumbar spine vBMD (see online supplementary Appendix Table S1), and there was no alteration of the results after further adjustment for estimated glomerular filtration rate or propensity score representing Ca, P and Mg (see online supplementary Appendix Tables S2–S8).
* From adjusted robust regression analyses.
† Adjusted for age, BMI, physical activity, sedentariness score, smoking, educational level, hormone therapy use, age at menopause, and intakes of total energy, dietary carbohydrate as a percentage of energy, Ca, P, Mg and alcohol.
‡ Total protein model additionally adjusted for vegetable protein intake as a percentage of energy.
§ Total protein model additionally adjusted for animal protein intake as a percentage of energy.
We also made adjustment for multiple comparisons to limit the potential for false-positive findings. After adjustment for multiple comparisons, the observed associations did not reach conventional levels of significance.
Discussion
In the present multi-ethnic cohort of middle-aged men and women from the MESA cohort, sex and race/ethnicity jointly modified the association of dietary protein with lumbar spine trabecular vBMD. No detrimental effects of high-protein diets on vBMD were observed. In contrast, the most consistent association was that of a higher vegetable protein intake substituting for fat with a greater lumbar spine vBMD among non-Hispanic white women, independent of age, BMI, smoking, physical activity, postmenopausal hormone use, age at menopause, and intakes of total energy, carbohydrate, Ca, P, Mg and alcohol. However, this association was not statistically significant after adjustment for multiple comparisons.
It is still unknown whether the relationship between protein intake and bone health is modified by Ca intake. We did not observe a significant interaction between dietary Ca and protein intake in our sample. In a randomised controlled trial among 342 healthy men and women (mean age 71 years; mean total protein intake 79·1 g/d), Dawson-Hughes & Harris( Reference Dawson-Hughes and Harris 40 ) reported that the association of total protein intake with a 3-year change in total and femoral neck BMD was significantly different between groups with and without Ca/vitamin D supplementation, but not for lumbar spine BMD. A Ca balance study among fifteen postmenopausal, mostly white, women (mean age 61 years) showed that urinary Ca loss or indicators of bone metabolism did not differ between low-protein (mean total protein intake 68 g/d) and high-protein (mean total protein intake 117 g/d) diets( Reference Roughead, Johnson and Lykken 41 ). In contrast, a prospective cohort study among 960 older men and women has shown that increasing protein consumption appeared to be more beneficial to women with lower Ca intakes (mean age 71 years)( Reference Promislow, Goodman-Gruen and Slymen 17 ).
Our findings with regard to total and animal protein intakes are consistent with several studies that examined this association among postmenopausal women( Reference Thorpe, Mojtahedi and Chapman-Novakofski 10 , Reference Wang, Luz Villa and Marcus 22 , Reference Cooper, Atkinson and Hensrud 42 ). In contrast, Teegarden et al. ( Reference Teegarden, Lyle and McCabe 43 ) and Quintas et al. ( Reference Quintas, Ortega and Lopez-Sobaler 44 ) found a positive association of total protein intake with lumbar spine BMD among non-Hispanic white women, but their populations were younger and more physically active. A study conducted by Whiting et al. ( Reference Whiting, Boyle and Thompson 45 ) has also shown a significant correlation between total protein intake and BMD in men but did not control for important potential confounders, such as age and physical activity. In a cross-sectional survey, Rapuri et al. ( Reference Rapuri, Gallagher and Haynatzka 46 ) reported that women in the highest and second highest quartiles of total protein intake had higher lumbar spine BMD, compared with their counterparts in the lowest quartile; however, this association was not significant in their longitudinal analysis. Darling et al. ( Reference Darling, Millward and Torgerson 47 ) conducted a meta-analysis of randomised controlled trials and found that participants on protein supplementation had a slightly higher lumbar spine areal BMD compared with those on placebo supplementation, with a weighted mean difference of 0·02 (95 % CI 0·00, 0·04) g/cm2; however, only three small studies (overall n 110) were included, and sex or racial/ethnic effects were not evaluated. The present results suggest that BMD is not adversely affected by high-protein diets, and build on these prior findings by suggesting that the relationship between dietary protein and bone density may differ by race/ethnicity and protein source, which may partially account for conflicting results.
With regard to the initial findings of vegetable protein intake and bone health in non-Hispanic white women, a cross-sectional study examining the association among 560 white women enrolled in a managed care organisation has shown the same trend: women in a higher tertile of vegetable protein intake had a higher BMD at each skeletal site (P for trend: 0·03 for hip, 0·10 for lumbar spine and 0·04 for whole body) than did women consuming less vegetable protein( Reference Beasley, Ichikawa and Ange 48 ). Among elderly Chinese women, no association between vegetable protein intake and lumbar spine BMD was identified in one study( Reference Lau, Kwok and Woo 49 ), while another( Reference Chiu, Lan and Yang 24 ) found that vegetable protein intake was positively associated with lumbar BMD in postmenopausal Taiwanese vegetarian women, a population with very different dietary characteristics than the present study population. However, none of the aforementioned results were adjusted for the number of comparisons. Although the basis for these ethnic differences is not clear, the findings as a whole support a neutral or beneficial role of dietary protein on body mineral density.
Mechanisms by which protein intake may influence bone health include increasing Ca re-absorption or regulating plasma insulin-like growth factor-1 that increases bone formation( Reference Ginty 5 ). Isoflavones that are usually contained in plant-based protein diets may influence bone health, although their protective effects as dietary supplements have not been consistently shown in data from randomised controlled trials( Reference Ricci, Cipriani and Chiaffarino 50 , Reference Liu, Ho and Su 51 ). In contrast, protein intake releases protons during the oxidation of sulphur-containing amino acids, and thus lowers pH and increases bone dissolution( Reference Thorpe, Mojtahedi and Chapman-Novakofski 10 , Reference Sellmeyer, Stone and Sebastian 11 , Reference Rizzoli and Bonjour 52 ). Because animal protein is rich in acidifying amino acids such as cysteine and methionine, whereas vegetable protein is not, the effect of protein intake on BMD could be different by protein source. However, the differences between racial/ethnic groups could also be due to differences in physiological characteristics and the amount and duration of protein intake.
The results should be interpreted with caution as conclusions from the present study are limited by the following factors. First, this is an observational study so residual confounding cannot be excluded. Second, dietary assessment methods usually include measurement error that might affect both accuracy and precision of study results. Measurement error in the assessment of protein intake is likely to be non-differential because diet was measured before data collection of lumbar spine vBMD. Third, the type 2 error rate may be increased in some sex- and race/ethnicity-specific subgroups due to reduced sample sizes. However, we conducted power analysis and found at least 80 % power to detect an association in each racial/ethnic group with the exception of 76 % power among Chinese-American women. Fourth, we lacked data on Vitamin D supplementation, medicines for osteoporosis, vBMD on other bone sites, and biochemical markers of bone health. Finally, our findings may not be directly applicable to the elderly because elderly adults require more protein consumption than those in late middle age (mean age of the study population 62 years) in order to maintain N equilibrium, heal wounds or pressure ulcers, fight infections, or build new tissue to replace normal losses. Recent evidence suggests that the minimum requirement of dietary protein is 1·0 g/kg body weight among elders, when compared with 0·8 g/kg recommended for middle-aged adults( Reference Chernoff 53 ). The findings may also not apply to athletes due mainly to increased requirements for protein with high levels of physical activity( Reference Phillips and Van Loon 54 ).
The present study also has several strengths that lend confidence to these findings. First, we examined vBMD in a well-characterised sample of multiple ethnicities, which enabled us to investigate and compare the association in participants of diverse racial and ethnic origins. Second, MESA followed a stringent protocol, and collected measures using standardised procedures to increase the measurement validity and reliability. Third, validated questionnaires were used for the assessment of habitual dietary intake and confounding variables. Fourth, we adjusted for multiple comparisons to better control type 1 error rates. In addition, information on a variety of potential confounders, including physical activity, BMI, smoking, alcohol consumption, hormone replacement therapy and other dietary factors, was available for adjustment.
In summary, data from the present large, multi-ethnic, population-based study suggest that a higher level of protein intake, when substituted for fat, is not associated with poor bone health among late middle-aged men and women. On the contrary, a higher intake of vegetable protein may be beneficial for bone health among middle-aged non-Hispanic white women. Differences in the relationship between protein source and race/ethnicity of study populations may in part explain inconsistent findings reported previously. Given the prevalence of low bone mass among older women, these results have important public health implications. Additional studies with more detailed dietary assessments are needed to confirm these findings and to examine bone formation or reabsorption markers that may clarify the potential mechanisms for racial/ethnic differences in bone health.
Supplementary material
To view supplementary material for this article, please visit http://dx.doi.org/10.1017/S0007114514002220
Acknowledgements
The authors thank the other investigators, the staff and the participants of the MESA study for their valuable contribution. A full list of participating MESA investigators and institutions can be found at http://www.mesa-nhlbi.org
The present study was supported by contracts N01-HC-95159 through N01-HC-95166 and N01-HC-95169 from the National Heart, Lung, and Blood Institute (NHLBI) for the collection of data in the MESA study. L. A. B. was supported in part by K08 HL091108 from the NHLBI and in part by P20RR017659 from the National Center for Research Resources (NCRR) through the Tulane Hypertension and Renal Center of Excellence. Both NHLBI and NCRR are components of the National Institutes of Health (NIH). The contents herein are solely the responsibility of the authors and do not necessarily represent the official view of NHLBI, NCRR or NIH.
The authors’ responsibilities are as follows: T. H. and L. A. B. contributed to the study concept and design and analysis and interpretation of the data; D. R. J. and M. H. C. contributed to the acquisition of the data; N. J. R., J. A. N., J. A. H., L. M. S., J. H., D. R. J., M. H. C. and L. A. B. provided critical revision of the manuscript for important intellectual content; T. H. contributed to the drafting of the manuscript; L. A. B. had full access to all of the data in the study, and responsibility for the integrity of the data and the accuracy of the data analysis.
The authors declare that they have no conflicts of interest.