点击显示 收起
1 From the Division of Epidemiology and Community Health, School of Public Health, University of Minnesota, Minneapolis, MN (JAN, LMS, and DRJ); the Department of Nutrition, University of Oslo, Oslo, Norway (DRJ); the Department of Epidemiology, German Institute of Human Nutrition Potsdam-Rehbruecke, Nuthetal, Germany (MBS); the Department of Pathology, University of Vermont, Burlington, VT (NSJ); the Department of Medicine, College of Physicians and Surgeons, Columbia University, New York, NY (RGB); and the Department of Public Health Sciences, Wake Forest University, Winston-Salem, NC (AGB)
2 Supported by training grant T32 HL07779 (to JAN) and contracts N01-HC-95159 through N01-HC-95166 from the National Heart, Lung, and Blood Institute, and General Clinical Research Center grant M01-RR00645 from the National Center for Research Resources. 3 Reprints not available. Address correspondence to JA Nettleton, University of Minnesota, Division of Epidemiology and Community Health, School of Public Health, 1300 South 2nd Street, Suite 300, Minneapolis, MN 55454. E-mail: nettleton{at}epi.umn.edu.
ABSTRACT
Background: The association between diet and cardiovascular disease (CVD) may be mediated partly through inflammatory processes and reflected by markers of subclinical atherosclerosis.
Objective: We investigated whether empirically derived dietary patterns are associated with coronary artery calcium (CAC) and common and internal carotid artery intima media thickness (IMT) and whether prior information about inflammatory processes would increase the strength of the associations.
Design: At baseline, dietary patterns were derived with the use of a food-frequency questionnaire, and inflammatory biomarkers, CAC, and IMT were measured in 5089 participants aged 45–84 y, who had no clinical CVD or diabetes, in the Multi-Ethnic Study of Atherosclerosis. Dietary patterns based on variations in C-reactive protein, interleukin-6, homocysteine, and fibrinogen concentrations were created with reduced rank regression (RRR). Dietary patterns based on variations in food group intake were created with principal components analysis (PCA).
Results: The primary RRR (RRR 1) and PCA (PCA factor 1) dietary patterns were high in total and saturated fat and low in fiber and micronutrients. However, the food sources of these nutrients differed between the dietary patterns. RRR 1 was positively associated with CAC [Agatston score >0: OR (95% CI) for quartile 5 compared with quartile 1 = 1.34 (1.05, 1.71); ln(Agatston score = 1): P for trend = 0.023] and with common carotid IMT [1.0 mm: OR (95% CI) for quartile 5 compared with quartile 1 = 1.33 (0.99, 1.79); ln(common carotid IMT): P for trend = 0.006]. PCA 1 was not associated with CAC or IMT.
Conclusion: The results suggest that subtle differences in dietary pattern composition, realized by incorporating measures of inflammatory processes, affect associations with markers of subclinical atherosclerosis.
Key Words: Dietary patterns • principal components analysis • reduced rank regression • carotid artery intima media thickness • coronary artery calcium
INTRODUCTION
Coronary artery calcium (CAC) and carotid artery intima media thickness (IMT) are considered early indicators of atherosclerotic disease (1-3) and have been associated with cardiovascular disease (CVD) events in large prospective cohort studies (4-7). Inflammatory biomarkers, such as C-reactive protein (CRP), interleukin-6 (IL-6), homocysteine, and fibrinogen, are associated with the development of overt CVD (8-13). Dietary intake also plays a role in the development of CVD, and studies have reported significant relations between dietary factors and subclinical atherosclerosis and inflammation. Plasma vitamin C (14) and dietary linolenic acid (15) have been shown to be inversely associated, whereas dietary saturated fat (16) has been shown to be positively associated, with CAC. Fiber and whole-grain intake have been shown to be associated with reduced progression of carotid artery IMT (17) or coronary artery stenosis (18). Diets rich in combinations of these individual components have also shown inverse associations with subclinical atherosclerosis markers (19-21). Randomized, multifactorial lifestyle interventions also support a beneficial role for dietary patterns in relation to the progression of atherosclerosis (22, 23). Similar associations between inflammatory biomarkers and nutrients, foods, and dietary patterns also exist (24-29), which suggests a potential intermediate pathway between diet and markers of atherosclerosis and, thus, CVD.
Dietary pattern analysis facilitates the examination of food and nutrient synergies in relation to disease, whereas traditional reductive approaches may not capture such synergies and associations (30, 31). Recently, Hoffmann et al (32) proposed the application of reduced rank regression (RRR) in dietary pattern analysis. In contrast with the commonly used principal components analysis (PCA), which produces a linear combination of food groups (a factor) that maximally explains variation in food group intake, RRR produces a linear combination of food groups that maximally explains variation in investigator-specified intermediate response variables. Thus, the technique allows hypotheses regarding pathways (represented by the response variables) between diet and disease to be evaluated (33). Although PCA has been successful, it may not achieve the full potential of dietary pattern analysis because it relies solely on intercorrelations among dietary variables, which may not optimally represent the diet qualities most relevant to specific disease etiology (32).
In the present study we used RRR to derive dietary patterns on the basis of their ability to maximally explain variation in biochemical markers related to inflammation (ie, CRP, IL-6, fibrinogen, and homocysteine). A comparison of the food group composition of behavior-based PCA patterns with the food group composition of inflammation-based RRR patterns will provide an interesting perspective on the value of the 2 methods, ie how much information is lost if patterns are based on behavior alone and do not take into account inflammatory processes? On the basis of previous successes of PCA in this (26) and other (34) cohorts, we hypothesized that both PCA and RRR techniques would produce dietary patterns associated with CAC and carotid artery IMT. However, by incorporating inflammatory intermediates in the RRR analysis, we hypothesized that the RRR method would produce dietary patterns more strongly associated with CAC and IMT than dietary patterns derived by PCA.
SUBJECTS AND METHODS
Subjects
The primary objectives of the Multi-Ethnic Study of Atherosclerosis (MESA) are to study the characteristics of subclinical CVD and the risk factors that predict progression to clinically overt CVD or progression of subclinical disease (35). The MESA population consists of an ethnically diverse group of 6814 white, black, Hispanic, and Asian men and women aged 45–84 y, who were recruited from 6 field centers across the United States (Baltimore City and County, MD; Chicago, IL; Forsyth County, NC; New York, NY; Los Angeles County, CA; and St Paul, MN). All participants were free of clinical CVD at baseline and gave informed consent. The MESA protocols were approved by each field center's institutional review board.
The current cross-sectional investigation includes data from 5089 participants enrolled in MESA, including 2407 men and 2682 women aged 45–84 y. This sample was chosen after the exclusion of individuals who were currently taking oral steroid or antiinflammatory asthma medications (n = 134), who had diabetes mellitus defined according to the American Diabetes Association 2003 criteria (n = 919) (36), and who provided insufficient or implausible dietary information (n = 630: 167 because of extreme energy intake and 463 because of implausible responses or improperly completed questionnaires; 26).
Carotid intima media thickness and coronary calcium
Thickness of the intima media layers of the common and internal carotid arteries was determined from images obtained by high-resolution B-mode ultrasonography (Logiq 700 ultrasound machine; GE Medical Systems, Waukesha, WI). Calculation of IMT was performed at the central MESA ultrasound reading center (Tufts–New England Medical Center, Boston, MA). In the present study, high IMT was defined as measures above the 75th percentile of the sample distribution, separately for each internal carotid IMT (1.3 mm) and common carotid IMT (1.0 mm).
Computed tomography (CT) of the coronary arteries was performed with cardiac-gated (at 80% of the R-R interval) electron beam scanners at 3 centers (Imatron C-150; Imatron Inc, San Francisco, CA) (37) and with a prospective electrocardiogram-triggered scan acquisition at 50% of the R-R interval with multidetector scanners (38) at the remaining 3 centers. These scanners are comparable in their ability to measure calcium (39). All scans were read at the central MESA CT reading center (Harbor-University of California Medical Center; Los Angeles, CA). CAC scores (Agatston scores) were determined by blinded CT analysts. Presence of CAC was defined as an Agatston score >0.
CRP, fibrinogen, interleukin-6, and homocysteine
CRP, fibrinogen, IL-6, and homocysteine concentrations were measured in blood samples collected at baseline, processed (40), and stored at –80 °C until analyzed. CRP and fibrinogen antigen were measured in plasma with a BNII nephelometer (N High Sensitivity CRP and N Antiserum to Human Fibrinogen; Dade Behring Inc, Deerfield, IL). IL-6 concentrations were measured by ultrasensitive enzyme-linked immunosorbent assay (Qantikine HS Human IL-6 Immunoassay; R&D Systems, Minneapolis, MN). Plasma homocysteine was measured by polarization immunoassay with an IMx Analyzer (IMx Homocysteine Assay; Axis Biochemicals ASA; Oslo, Norway). Analytic CVs were 3.6%, 2.6%, 6.3%, and 4.5% for CRP, fibrinogen antigen, IL-6, and homocysteine, respectively.
Dietary assessment
At the baseline examination, each participant completed a self-administered, modified-Block style food-frequency questionnaire (FFQ) (26, 41, 42). The FFQ was provided in English, Spanish, or Chinese with staff assistance available on request. Responses were converted to servings/d and categorized into 47 food groups.
Dietary pattern analysis
Dietary patterns were derived from 47 food groups with the use of 2 different empirical methods: PCA (26) and RRR (32) (SAS PROC FACTOR with varimax rotation and SAS PROC PLS with the RRR method option, respectively). In the RRR analysis, the specified response variables were CRP, IL-6, fibrinogen, and homocysteine (natural log values due to non-normal distribution of original values). For the PCA analysis, a 4-pattern solution was chosen after the evaluation of a scree plot of the eigenvalues and after the interpretability of the final solution was considered. For the RRR analysis, a 4-pattern solution was used in accordance with the number of response variables specified (32). The 4-factor PCA solution explained 26.2% of the variation in food group intake and 1.5% of the variation in response variables. The 4-factor RRR solution explained 12.6% of the variation in food group intake and 5.0% of the variation in response variables. To present data concisely and because other dietary patterns were not significantly associated with the outcomes, the results are presented for only the primary factor (dietary pattern) derived by the PCA and the primary factor derived by RRR.
A score was calculated for each participant for each PCA and RRR dietary pattern as a sum of the 47 food groups, each weighted according to the factor loadings. Participant scores were then categorized into quintiles separately for each dietary pattern. Thus, for each dietary pattern, quintile 5 was composed of persons whose diets conformed most closely to that particular pattern.
Assessment of other relevant variables
Information on demographics, education, medication use, smoking history, and physical activity were collected at the baseline examination with a combination of self-administered and interviewer-administered questionnaires. Waist circumference (cm) was measured at the umbilicus, and body mass index (BMI; in kg/m2) was calculated from measured height and weight.
Statistical analysis
All analyses were conducted with SAS version 9.1 (SAS Institute Inc, Cary, NC). The percentage variance in total dietary pattern score explained by each food group was calculated as follows (32):
RESULTS
Dietary patterns, foods, and nutrients
Pearson correlation coefficients between the dietary patterns and food groups that explained most of the variance in dietary pattern scores are presented in Table 1. As reported previously (26), the first PCA factor, "fats and processed meats" (PCA 1), was high in fats and oils, processed meats, fried potatoes, and desserts. RRR 1 was somewhat similar to PCA 1 and also high in processed meats and fats and oils, but other unique contributors to this pattern were nondiet sodas (high consumption) and soy foods and dark-yellow and cruciferous vegetables (low consumption). Percentage variance in scores explained by all food groups for each diet pattern are given in Appendixes A and B . PCA 1 and RRR 1 scores were significantly correlated (0.41; P < 0.001); however, the percentage agreement of quintile classification between PCA 1 and RRR 1 was low (32.3%).
View this table:
TABLE 1. Correlations between empirically derived dietary patterns and food groups, nutrients, and biomarkers of inflammation in 5089 men and women from the Multi-Ethnic Study of Atherosclerosis1
View this table:
Appendix A. Food group statistics for primary dietary patterns empirically derived by principal components analysis (PCA) in 5089 men and women from the Multi-Ethnic Study of Atherosclerosis
View this table:
Appendix B. Food group statistics for primary dietary patterns empirically derived by reduced rank regression (RRR) in 5089 men and women from the Multi-Ethnic Study of Atherosclerosis
Partial correlations between dietary patterns and nutrients were consistent with the food groups largely composing the patterns (Table 1). PCA 1 was positively correlated with total, saturated, monounsaturated, and trans fats but negatively correlated with fiber, carbohydrate, and several vitamins. RRR 1 was also positively correlated with saturated and trans fats, but correlations with total fat were not as strong as for PCA 1, and RRR 1 positively correlated with carbohydrate intake. However, carbohydrate contributors were of low quality, as evidenced by the negative correlations between RRR 1 and fiber, vitamin C, folate, and ß-carotene.
Dietary patterns and CRP, IL-6, fibrinogen, and homocysteine concentrations
Both PCA 1 and RRR 1 were significantly and positively correlated with CRP, IL-6, fibrinogen, and homocysteine (Table 1), although the magnitudes of the correlations with RRR 1 were approximately twice as large as those with PCA 1. The percentages of total biomarker variance (RRR response variables) explained by each PCA 1 (0.23%) and RRR 1 (3.0%) also differed to a similar degree. PCA 1 explained 0.55%, 0.36%, 0.002%, and 0.01% of the variation in CRP, IL-6, fibrinogen, and homocysteine, respectively, whereas for RRR 1, these values were considerably higher (5.3%, 4.0%, 2.4%, and 0.34%, respectively).
Dietary patterns and participant characteristics
The distribution of participant characteristics across quintiles of PCA 1 was different from that of RRR 1 (Table 2). For example, educational attainment was positively associated with PCA 1 scores, whereas educational attainment was negatively associated with RRR 1 scores. The directions of the associations for age, sex, and percentage Hispanic also differed between the 2 patterns (P for trend < 0.001 for all). However, both patterns were associated positively with inactive leisure, waist circumference, BMI, and smoking prevalence and negatively with supplement use (P for trend < 0.001 for all).
View this table:
TABLE 2. Descriptive characteristics of quintiles (Q) 1, 3, and 5 for primary dietary patterns derived by principal components analysis (PCA) and reduced rank regression (RRR) in 5089 men and women from the Multi-Ethnic Study of Atherosclerosis1
Dietary patterns and markers of subclinical atherosclerosis
The PCA-derived dietary pattern was not associated significantly with internal or common IMT or the presence of CAC after adjustment for demographic and lifestyle confounders (Table 3). However, RRR 1 was associated positively with the presence of CAC (Agatston score >0), independent of demographic and lifestyle confounders (P for trend = 0.011; Table 3). For persons in the fifth quintile, the odds of having CAC were 34% higher than for persons in the first quintile (OR: 1.34; 95% CI:.05, 1.71). Formal tests of interaction between dietary patterns and race, sex, current smoking status, waist circumference, or BMI were not statistically significant. The estimate for risk of high common carotid IMT based on RRR 1 scores in all participants was similar in magnitude to that of CAC > 0 (OR: 1.33), but the 95% CIs were wide (0.99, 1.79). If potential mechanistic variables were added to the main multivariable model, risk estimates for CAC > 0 were attenuated—modestly after adjustment for waist circumference [quintile 5 versus quintile 1 (OR: 1.28; 95% CI: 1.00, 1.64)] and more notably after adjustment for systolic blood pressure, HDL cholesterol, LDL cholesterol, triacylglycerols, insulin, and glucose [quintile 5 versus quintile 1 (OR: 1.19; 95% CI: 0.93, 1.53)], which suggests that any relations between this particular dietary pattern and CAC are mediated through the effects of diet on traditional CVD risk factors other than inflammatory processes. Results of regression analysis of the continuous variables [mean internal IMT (mm), common IMT (mm), and CAC (Agatston score)] across dietary pattern quintiles were similar to those for the dichotomized analyses. Mean CAC (Agatston score) and mean common carotid IMT increased across quintiles of RRR 1 (P for trend = 0.023 and 0.006 for CAC and common carotid IMT, respectively), although the range of values across quintiles was rather narrow (2% increase in common carotid IMT and 29% increase in CAC; Table 4). Again, PCA 1 was not associated with mean CAC or IMT (data not shown). As a supplementary sensitivity analysis, we also evaluated multivariable-adjusted associations between continuously modeled exposures (dietary patterns) and outcomes (log mean common and internal carotid IMT and Agatston score). Again, results showed a significant positive association between RRR 1 and CAC (ß ± SE: 0.064 ± 0.031; P = 0.039) and a nearly significant positive association between RRR 1 and common carotid IMT (0.002 ± 0.001; P = 0.054), but not internal carotid IMT (0.0007 ± 0.003; P = 0.82). PCA 1 showed no significant associations with outcomes (P = 0.43–0.98).
View this table:
TABLE 3. Odds ratios (and 95% CIs) for risk of subclinical atherosclerosis according to quintiles (Q) of dietary pattern score for primary dietary patterns empirically derived with principal components analysis (PCA) and reduced rank regression (RRR) in 5089 white, black, Hispanic, and Chinese men and women from the Multi-Ethnic Study of Athersclerosis1
View this table:
TABLE 4. Intima media thickness (IMT) of the internal carotid artery and common carotid artery (CAC) and coronary artery calcium (CAC) scores according to quintiles (Q) of dietary pattern score for primary dietary patterns empirically derived by principal components analysis (PCA) and reduced rank regression (RRR) in 5089 white, black, Hispanic, and Chinese men and women from the Multi-Ethnic Study of Atherosclerosis1
Food groups and markers of subclinical atherosclerosis
No single food group was able to explain the associations ascribed to the RRR dietary pattern. Although age- and energy-adjusted models showed significant associations between food groups and markers of subclinical atherosclerosis in the expected directions [quintile 5 versus quintile 1: processed meat consumption (OR: 1.38; 95% CI: 1.13, 1.68); cruciferous vegetable consumption (OR: 0.82; 95% CI: 0.67, 0.99)], values were attenuated and 95% CIs included the null after full adjustment for demographic and lifestyle characteristics [processed meat consumption: (OR: 1.11; 95% CI: 0.89, 1.39); cruciferous vegetable consumption (OR: 1.10; 95% CI: 0.87, 1.39)].
DISCUSSION
In this cross-sectional analysis, a dietary pattern based on variation in CRP, IL-6, fibrinogen, and homocysteine (RRR 1) was significantly, although modestly, associated with IMT and CAC. However, a dietary pattern based on variation in food group intake (PCA 1) was not significantly associated with these subclinical disease markers. The RRR dietary pattern was high in fats and oils, processed meats, and soda and low in soy and vegetable food groups. Consistent with the hypothesized ill effects of diets rich in saturated and trans fats and low in fiber and micronutrients, RRR 1 scores were positively associated with the presence of CAC, marginally associated with high common carotid IMT, and associated with continuous measures of common carotid IMT and CAC. The results were independent of demographic and lifestyle confounders, but not of classic CVD risk factors, nor of CRP, IL-6, fibrinogen, or homocysteine (by design).
The associations noted in relation to our RRR dietary pattern are similar to those reported in the one other study that investigated the relation between empirically derived dietary patterns and IMT (20). In a cluster analysis in a cohort of white women, Millen et al (20) reported that women falling into dietary pattern clusters that were high in total and saturated fats and low in micronutrients and fiber had an increased risk of carotid artery atherosclerosis relative to women in a low-fat, high-fiber, and micronutrient-rich "heart healthy" dietary pattern cluster. Because our RRR dietary pattern was low in fiber and antioxidant vitamins but high in saturated fat, our findings are also consistent with the studies reporting the beneficial (14, 17, 18) and deleterious (saturated fat) (16) effects of these nutrients in relation to IMT, CAC, or both. Although our PCA dietary pattern also shared similar characteristics, it was not as discriminating as was the RRR dietary pattern in terms of inflammatory biomarker concentrations, which possibly explained the absence of associations between PCA 1 and IMT or CAC.
RRR applied to the nutritional discipline is currently fairly novel. One of our goals was to add to this smaller body of research by comparing the associations between IMT or CAC and an RRR dietary pattern with the associations between IMT or CAC and a dietary pattern derived by a more common method, PCA. We hoped that such a comparison would increase knowledge of the RRR method and also further our understanding of diet-disease relations. We found the RRR dietary pattern to be associated with CAC and common carotid IMT, whereas the PCA dietary pattern was not significantly associated with these markers. A comparison of these 2 dietary patterns may aid in identifying specific food groups with biological effects relevant to inflammation and the development of atherosclerosis. Several nutrients showed similar correlations with RRR 1 and PCA 1 (saturated fat, trans fat, fiber, and antioxidant vitamins), and intakes of the food groups fats and oils and processed meats were important determinants of both RRR 1 and PCA 1 dietary pattern scores. However, the consumption of cruciferous and dark-yellow vegetables, soy foods, and soda notably differed between RRR 1 and PCA 1 (as well as other food groups not among the highest score contributors; Appendixes A and B). It is possible that such subtle differences in dietary pattern composition are biologically relevant to the atherosclerotic process. For example, Adams et al (43) recently showed that diets high in green and yellow vegetables (freeze-dried peas, green beans, broccoli, corn, and carrots) inhibited aortic atherosclerosis in a mouse model of atherosclerosis. Soy protein-feeding studies in monkeys fed otherwise atherogenic diets showed reduced atherosclerosis compared with milk proteins (44, 45), which may have been due partly to modified monocyte-endothelial cell interactions in the early stages of atherosclerosis (46). Although soda consumption has not been evaluated previously in relation to atherosclerosis, its association with obesity and type 2 diabetes is well known (47, 48); thus, a negative effect on atherosclerosis would not be unexpected. However, in the present study, individual food groups were not associated with IMT or CAC, which suggests that all facets of the dietary pattern act synergistically. One important role of dietary pattern analysis is to guide further investigations at the food and nutrient level (30). Although it is likely that the effects of single constituents are partly dependent on the overall dietary context, randomized clinical studies are needed to disentangle the effects of single components from those of the larger dietary pattern.
The limitations of this study deserve discussion. First, this cross-sectional analysis may be confounded by recent dietary changes that are not reflective of lifetime dietary practices, which are influential in the development of atherosclerosis. Similarly, these cross-sectional data do not allow for mechanistic hypotheses to be tested, ie, "does diet impact the atherosclerotic process by influencing inflammatory processes? " We can only conjecture that all factors (diet, inflammation, and atherosclerotic disease) may be interrelated. Second, as a result of the number of statistical tests conducted, noted associations may have been significant by chance alone. Indeed, the ORs were small but were robust to adjustment for multiple confounders and were consistent across strata of several characteristics. Third, we caution that ORs overestimate the relative prevalence when outcomes are common, as was the case for CAC presence (48% of cohort). Fourth, it is possible that the collection of biomarkers used in our RRR procedure did not adequately represent inflammatory processes relevant in the development of atherosclerosis. It is possible that the use of different biomarkers may have produced larger risk estimates, such as those reported in relation to coronary artery disease when more traditional CVD risk factors were used (33). However, we chose these biomarkers because of general interest in novel disease pathways. Last, the relation between diet and atherosclerosis may be mediated through pathways other than inflammation, and likewise, the relation between diet and CVD events may be mediated through pathways other than those reflected by carotid artery IMT and CAC (49).
Because RRR patterns are not based on dietary behavior, the applicative value of the resulting patterns in real-world settings remains largely unknown. In contrast, PCA dietary patterns reflect real-world food usage patterns but do not take into account biological measures of health (eg, inflammatory markers). Nevertheless, PCA-derived dietary patterns have consistently shown relevant associations with many CVD-related outcomes (34), including subclinical biomarkers whose concentrations are typically unknown by participants (25, 26). The unresolved question then becomes whether it is more important to identify a dietary pattern that reflects how people really eat and how that pattern is associated with disease outcomes, or alternatively, to identify an artificial dietary pattern (collection of foods) that is associated with disease outcomes. Likely, a valuable role for each method exists, depending on the question at hand—eg, PCA for planning behavioral interventions and RRR for etiologic investigations.
Dietary pattern research continues to be an area of interest because of the potential for food and nutrient synergies (30), utility of dietary patterns for making public health recommendations (50), and values of dietary pattern change (51-54). Our study is unique in that we were able to study a large multiethnic cohort, with a well-characterized measure of overall dietary intake and validated measures of inflammation and subclinical atherosclerosis. We found that a dietary pattern derived by RRR and based on variations in CRP, IL-6, fibrinogen, and homocysteine was positively associated with CAC and common carotid IMT, but a dietary pattern derived by PCA was not significantly associated with CAC or carotid artery IMT. Food group differences between these patterns did not solely explain the difference in relations to outcomes, which suggests the importance of the overall dietary pattern and method of dietary pattern derivation.
ACKNOWLEDGMENTS
We thank the other investigators, the staff, and the participants of MESA for their valuable contributions. A full list of participating MESA investigators and institutions can be found at http://www.mesa-nhlbi.org.
The authors' responsibilities were as follows—JAN: analytic design, data analysis, and manuscript preparation; LMS: critical review of the manuscript; MBS: data analysis and critical review of the manuscript; NSJ, RGB, and AGB: data acquisition and manuscript review; DRJ: data acquisition, critical review of the manuscript, and data analysis. None of the authors had a conflict of interest to report.
REFERENCES