- Research article
- Open Access
- Open Peer Review
A worksite intervention to reduce the cardiovascular risk: proposal of a study design easy to integrate within Italian organization of occupational health surveillance
© Mastrangelo et al.; licensee BioMed Central. 2015
- Received: 29 July 2014
- Accepted: 8 January 2015
- Published: 21 January 2015
Despite the substantial amount of knowledge on effectiveness of worksite health promotion (WHP) in reducing cardiovascular disease (CVD) risk, WHP programs are not systematically applied in Italy. The aim was to design an intervention easy to integrate within the Italian organization of workplace health surveillance.
We used the “pretest-posttest design”. Workers were employed in multiple occupations and resident in Veneto region, Italy. Occupational physicians (OPs) performed all examinations, including laboratory evaluation (capillary blood sampling and measure of glycaemia and cholesterolemia with portable devices), during the normal health surveillance at worksite. CVD risk was computed based on sex, age, smoking habit, diabetes, systolic pressure and cholesterol level. After excluding those with <40 years of age, missing consent, CVD diagnosis or current therapy for CVD, missing information, CVD risk <5%, out of 5,536 workers 451 underwent the intervention and 323 male workers were re-examined at 1 year. CVD risk was the most compelling argument for changing lifestyle. The counseling was based on the individual risk factors. Individuals examined at posttest were a small fraction of the whole (6% = 323/5,536). In these workers we computed the ratio pretest/posttest of proportions (such as percent of subjects with cardiovascular risk >5%) as well as the exact McNemar significance probability or the exact test of table symmetry.
CVD risk decreased by 24% (McNemar p = 0.0000) after the intervention; in a sensitivity analysis assuming that all subjects lost to follow-up kept their pretest cardiovascular risk value, the effect (−18%) was still significant (symmetry p < 0.0000). Each prevented CVD case was expected to cost about 5,700 euro.
The present worksite intervention promoted favorable changes of CVD risk that were reasonably priced and consistent across multiple occupations.
- Worksite health promotion
- Cardiovascular disease
- Cardiovascular risk
- Pretest-posttest design
The worksite has been proposed by the World Health Organization (WHO) as a priority setting for health promotion in the 21st century: The worksite directly influences the physical, mental, economic and social well-being of workers and in turn the health of their families, communities and society. It offers an ideal setting and infrastructure to support the promotion of health of a large audience .
Worksite health promotion programs originated in US from executive fitness programs that were created in the years after World War II. Initiated by business leaders who endorsed the benefits of a healthful lifestyle, the number of in-house corporate programs grew steadily throughout the 1970s. During the next decades, employer benefits began to focus on management of prevalent chronic conditions (obesity, diabetes, heart disease, cancer, and depression) instead of focusing on fitness and were increasingly offered to employees of all job levels . Work health programs (WHPs) carried out in the past decade showed, in particular, promising results in contrasting the modifiable risk factors of cardiovascular diseases (CVDs) defined as: 1) physical inactivity, 2) tobacco use, 3) hypertension, 4) dyslipidemia, 5) poor diet, 6) hyperglycemia, and 7) elevated psychological stress . Several review papers have been recently issued on WHPs and cardiac rehabilitation, authored by leading authors from the United States , Canada , Brazil , Europe , India  and Japan . Although the delivery models, level of development/utilization, and legislative support varied among countries, these reviews clearly indicated that worksite health and wellness are important lifestyle intervention strategies and should be viewed as integral components of global healthcare with respect to combating CVDs .
Little has been done in Italy. Even the specific notion of WHP and lifestyle modification interventions is unknown within the Italian legal system, particularly in the recent set of rules for health and safety in workplaces contained in Legislative Decree 9/4/2008 No. 81 (updated in 2013) . A barrier is probably the fact that in Italy common diseases are entrusted to the general practitioner (public service), while occupational diseases are assigned to the occupational physician (private service). Employers are reluctant to support additional financial costs to improve their employees’ health, because this task falls to public health care organizations. On the other hand, employed individuals are unable to attend primary care services during the working day and may not wish to utilize their “citizen time” (time spent outside work) for this; in addition, males are less likely than their female counterparts to schedule annual health checks, seek medical advice, or attend educational meetings .
Despite the substantial amount of knowledge on effectiveness of WHPs, these interventions are not systematically applied in Italy. We therefore aimed to design an educational intervention that could be easy to integrate within the Italian current organization of occupational health surveillance and reasonably priced.
Quasi-experimental study designs, often described as pre-post intervention studies or before-and-after studies, are common in the medical literature. We used a quasi-experimental design and precisely the “one-group pretest-posttest design” .
The key outcome was reducing the cardiovascular risk over the next 10 years. The latter was computed with an algorithm − proposed by the European Society of Cardiology  and acknowledged by the Istituto Superiore di Sanità (Italian National Institute of Health) , combining information on sex, age, smoking habit, diabetes, blood pressure (mmHg), and blood cholesterol level (mg/dl). The modifiable risk factors (smoking, blood pressure, cholesterol), along with physical inactivity and alcohol intake were targeted by the intervention. There are multiple primary outcomes and this is, therefore, a multi-faceted worksite intervention to promote favorable changes in cardiovascular disease risk factors.
Study size was set at about 5,000 workers based on the available budget rather than power considerations. Workers were employed in a wide range of occupational sectors (private and public businesses, industry and services) and were resident in various provinces (Padova, Verona, Vicenza) of Veneto, Northeastern region of Italy. An intervention across multiple occupational groups from separate geographic communities could increase confidence that the intervention was responsible for a change in outcome. The ultimate rationale of this procedure was, therefore, assessing consistency of results.
Lastly, an unstructured interview (qualitative method) was conducted by a trained occupational physician (OP) among a small group of information-rich workers to answer the questions: «How did the intervention have that effect?»; and «What was the reaction of participants to the intervention?» .
Sampling frame and running the study
A snowball sampling method was used to select these workers. Two authors (LM and GM) chose the occupational physicians based on their scientific interests; and OPs chose the companies where they had the best relationships with both employers and employees. An invitation letter was sent to the management of these companies, explaining the aim and methods of the study; the companies willing to cooperate constituted the final sample of 5,536 workers.
Before investigation, OPs were trained on counseling techniques, mainly focused on diagnosing the worker’s motivational state to change risky behaviors. Each OP was given a fixed incentive (20 euro) for each worker examined.
All investigations were performed by OPs during the normal health surveillance and took place in the worksite. A computer aided interviewing software (Microsoft Access) was set up to store the data. Information was collected on: lifestyle factors (physical activity, cigarette smoking, alcohol consumption); past medical history (particularly, occurrence of CVDs, diabetes and obesity) and whether the subject was under therapy for diabetes, hypertension and hypercholesterolemia. At physical examination, blood pressure was measured twice at 4–5 minutes distance, always right arm and worker in upright position; the lowest value was used in the statistical analysis. The laboratory evaluation was performed in the workplace itself, collecting specimens of capillary blood and measuring glycaemia and cholesterolemia with portable devices. The procedure for collecting capillary blood specimens by fingerstick was that recommended by Centers for Disease Control . The cardiovascular risk over the next 10 years (CVD risk) was computed with the above algorithm and scored in classes (<5%; 5-10%; 10-15%; 15-20%; 20-30%; >30%). After 12 months, workers were re-examined with the same protocol.
The exams began after a letter of information to the supervisory body for workplace safety and health of the relevant Local Health Authority. The project was run from January 2011 up to December 2012. Clearance by Ethics Committee was not necessary because the study was a mandatory activity deliberated by Veneto Region with a formal act (Regional Decree n. 2008 3 Aug 2010). All workers signed an informed consent at enrollment.
The original 5,536 workers were divided in several subsets during the course of the study. We did not consider 2,062 subjects aged less than 40 years, 65 individuals not giving the consent, 36 patients already affected by CVD, and 537 under current therapy for hypertension, hypercholesterolemia and diabetes. 2,836 subjects older than 40 years without past CVD history or current therapy for diabetes, hypertension, and hypercholesterolemia underwent laboratory evaluation. 59 workers with missing information on one or more of the six components (sex, age, smoking, diabetes, blood pressure and cholesterol) used to estimate the cardiovascular risk, and 2,326 workers with a CVD risk below 5% were excluded. The remaining 451 underwent the educational intervention. Out of the latter, 330 workers (323 males and 7 females) with a CVD risk >5% were re-examined at 1 year, while 121 (26.8% = 121/451) were lost to follow-up. All analyses were carried out in the 323 males because seven subjects could not be used to arrive at any conclusions regarding female gender.
There were two aspects: motivation and education. The most compelling argument for changing lifestyle was the estimated risk of CVDs in the next 10 years. Then subjects received an individualized counseling based on the presence of risk factors. Physical activity was generally mistaken with “exercise” (activity that is planned, structured, repetitive, and purposeful). In agreement with World Health Organization , workers were recommended to do at least 150 minutes of moderate-intensity aerobic physical activity throughout the week (example 30 minutes 5 times/week). For diet, recommendations were to limit energy intake from total fats and shift from saturated fats to unsaturated fats, increase consumption of fruits, vegetables, legumes, whole grains and nuts, limit the intake of free sugars and limit salt consumption from all sources . Subjects with hypertension and/or hypercholesterolemia or hyperglycemia were interviewed about their attitude towards lifestyle change; whenever they could not cope to recommendations they were referred to their general practitioners for medical therapy, even when the CVD risk was lower than 20%. Likewise, most smokers with CVD risk >5% were addressed to receive an anti-smoking counseling from counselors with educational competence. In other words, we used an aggressive approach that combined both a primary and secondary prevention.
In order to determine whether the intervention had the intended effect we calculated proportions with the factor (such as, proportions of smokers, before and after), the ratio between proportions (point estimates and confidence intervals) and the exact McNemar significance probability (for 2 × 2 tables). For outcome variables with multiple discrete levels (k × k tables), we performed an exact test of table symmetry.
A sensitivity analysis was carried out, performing an exact test of table symmetry on 451 subjects that included 121 subjects (26.8%) lost to follow-up. The latter contributed to the analysis assuming that their pretest value of cardiovascular risk remained unchanged at 1 year.
We coded a binary variable (delta) that was 1 if pretest CVD risk was higher that posttest CVD risk, and 0 otherwise. A low value of delta seemingly indicates a worst impact of the intervention. Since figures became too sparse in the subset of 330 workers undergoing intervention, occupational categories were merged in four groups: “basic metals” (original category); “other industries” (multiple categories); hospital workers (original category); other service workers (multiple categories). Using delta as outcome we fitted two models of logistic regression where the predictors were age, gender and work sectors (model 1); or age, gender, posttest smoking, posttest blood cholesterol, posttest systolic blood pressure and work sectors (model 2). In all models the work sector with the lowest value of delta was the reference. Odds ratios (OR) with 95% confidence interval (CI) and p-value were calculated with Stata 13 (Stata Corporation, College Station, Texas, USA).
A statistical process  had determined a detectable characteristic (CVD risk) associated with an increased chance of experiencing future unwanted outcomes. By identifying risk factors before the occurrence of the event, we developed targeted interventions to mitigate their impact. In order to obtain the prevented cases of CVD expected by the intervention, we multiplied in each class of risk the median CVD risk by the number of subjects. The sum of the latter values were the cases expected at pretest (A) or posttest (B). The number of CVD cases potentially prevented by the intervention was the difference (A – B). The cost outcome analysis was obtained by dividing overall cost by the number of potentially prevented cases.
Number of subjects, percentage of males, mean and standard deviation (sd) of age by occupational category in the whole study and in subjects enrolled in the pretest-posttest study
Age (years) mean ± sd
Age (years) mean ± sd
41.5 ± 20.5
53.3 ± 12.9
46.3 ± 9.6
41.3 ± 8.8
45.3 ± 14.2
46.0 ± 4.2
38.3 ± 12.5
34.6 ± 8.4
40.6 ± 8.1
52.4 ± 2.5
38.2 ± 9.2
54.3 ± 3.1
36.5 ± 8.9
52.0 ± 7.2
41.0 ± 9.6
51.3 ± 4.6
40.4 ± 8.5
50.1 ± 3.5
40.8 ± 9.6
55.2 ± 3.5
40.8 ± 9.3
51.5 ± 4.8
46.3 ± 8.6
54.0 ± 4.8
40.7 ± 9.4
51.6 ± 1.7
47.4 ± 7.8
54.6 ± 4.4
43.6 ± 8.0
53.4 ± 2.0
43.0 ± 9.9
52.8 ± 4.2
42.0 ± 8.8
52.3 ± 3.5
42.1 ± 9.2
53.0 ± 3.8
42.3 ± 9.2
53.0 ± 4.0
Proportions with the factor (at pretest and posttest), point estimates and confidence intervals (95% CI) for the ratio between proportions and the exact McNemar significance p-value, criterion for positivity in the footnote
Ratio post/pre (95% CI)
McNemar exact p-value
Proportion with factor
Physical activity (PA)
1.46 (1.32 - 1.62)
Cigarette smoking (CS)
0.84 (0.78 - 0.92)
Alcohol drinking (AD)
0.86 (0.79 - 0.94)
Systolic blood pressure (SBP)
0.83 (0.74 - 0.92)
Blood cholesterol (BC)
0.85 (0.78 - 0.93)
Cardiovascular risk (CR)
0.76 (0.71 - 0.80)
In 108 posttest smokers cigarettes consumption decreased after the intervention; for example, heavy smokers (>20 cigarettes/day) were 43 before and 32 after the intervention. These changes were statistically significant (symmetry exact significance probability = 0.0135).
Sensitivity analysis on the risk of cardiovascular disease: pretest and posttest comparison and exact test of table symmetry in 451 subjects, including 121 lost to follow-up that were presumed to keep the pretest values
Work sectors merged in the analysis of 323 subjects undergoing the pretest-posttest study: number of subjects (N) and percent (%) with delta = 1 (pretest cardiovascular risk > posttest cardiovascular risk), Odds ratio (OR) with 95% confidence interval (95% CI) and p-value for work sectors according to three models of logistic regression (see footnote)
N (% with delta = 1)
Other service activities
Expected number of cardiovascular disease cases based on stratification of cardiovascular risk at pretest and posttest
Cardiovascular risk %
No. of subjects
No. of subjects
An experienced occupational physician interviewed a small sample of workers already known as “opinion leaders” in their respective groups; the latter reported that both workers and employers perceived the intervention as useful. When the occupational physician was told to express his personal view, he answered: “we have gained esteem of workers”.
Thus, quantitative statistically significant results and some qualitative evidence, together, suggested that the intervention had been effective.
Governmental agencies and private sector groups are working hard to help employers to improve the health of their employees in an efficient, integrated, and cost-effective way. The objective is clear; it is the “how to” that is difficult .
At present, the accepted gold standard for the evaluation of interventions in health care is the randomized controlled trial (RCT). The medical literature reports several RCTs on workplace health promotion programs. In a recent meta-analysis, a surprising observation is that studies with poor methodological quality reported an average effect size 2.9-fold larger than good-quality studies. Analyses stratifıed by outcome showed the same result for sickness absence, work productivity, and work ability. This might indicate publication bias: poor-quality studies are more frequently published if they show a greater effect . Another reason could be the fact that RCTs conducted in the worksite may be affected by a threat to internal validity that occurs when the intervention delivered to one group “diffuses” to another (contamination threat). This can easily happen when the intervention is educational in nature, since workers naturally share information with one another. A contamination is undesirable for an evaluation because it reduces the differences observed between the intervention and control groups . Therefore, the logistic requirements of RCTs often cause them to be unfeasible, especially for single smaller worksites.
Given these limitations, the Cochrane Effective Practice and Organization of Care Group endorses three alternative methodologies for evaluating population interventions: (1) the non-RCT, (2) the controlled before-and-after study, and (3) the interrupted time series design. In a non-RCT, individuals or groups are allocated to experimental conditions using a nonrandom method. While nonrandom allocation may be more convenient in some circumstances, it increases the probability that unmeasured characteristics that may influence the outcomes introduces a systematic bias that could artificially exaggerate, or reduce, true intervention effects .
It has been recently suggested that the optimal study design for a workplace health promotion program may be a quasi-experimental design in which medical cost data are collected for several years before the program and participants and nonparticipants are matched through propensity scoring .
Another approach, whose advantages and methodological limitations have been recently discussed, is multiple baseline design. It involves conducting multiple time-series in multiple populations, each of which receives the intervention at a different point in time . Like RCTs, the multiple baseline design can demonstrate that a change in behavior has occurred, the change is a result of the intervention, and the change is significant. Especially important practical advantages over the RCT are that, first, this design requires fewer population groups and, second, communities may act as their own controls . As explained in Methods, the present study was conducted in multiple occupational categories even though, because of time constraints, we could not stagger the intervention and all categories were examined concurrently.
Individuals examined at posttest was a small fraction of the whole (6% = 323/5,536). This decreased the cost of prevention (about 5,700 euro for each prevented cases of cardiovascular disease) but involved a before-and-after design of the study. The latter is a non-experimental approach that must be used with caution, because of circumstances that threaten the ability to correctly infer whether the intervention had the desired effect. When the basis for choosing the intervention group is a greater apparent need for the intervention, an alternative explanation for the apparent success of the safety initiative is “regression-to the-mean” . In the present study, the intervention included only subjects with a high risk of cardiovascular disease. Thus, part of any decrease observed may have nothing to do with the intervention itself. Rather, CVD risk could be simply fluctuating closer to the average (year-to-year fluctuations). Strictly speaking, however, this consideration applies when one group is being examined and one outcome is being evaluated. As explained in Methods, the present study is a multiple risk factor intervention conducted in multiple occupational categories. There was a consistent performance of all indicators of cardiovascular risk (Table 2). Despite their original heterogeneity, work sectors were not found to influence the posttest risk of cardiovascular disease after taking into account the changes in modifiable risk factors produced by the intervention (Table 4). On the other hand, the characteristics of the intervention group could be altered when enough people drop out of the study (dropout threat) . In the present study the before measurements were available; even in the extreme assumption that all dropouts kept their initial value of cardiovascular risk, a significant decrease of cardiovascular risk was observed at posttest (Table 3). Overall, these pieces of evidence might increase confidence that the intervention was responsible for the change in the outcome.
Cardiovascular risk can be viewed as a surrogate endpoint to investigate the primary event (cause-specific mortality). Adoption of surrogate criteria must, however, be regarded with some caution because the link between surrogate and primary event is not always linear; furthermore risk factor changes could not be maintained in the long term .
The concept of the health promoting workplace is becoming increasingly relevant as more private and public organizations recognize that future success in a globalizing marketplace can only be achieved with a healthy, qualified and motivated workforce … For nations, the development of HPW will be a pre-requisite for sustainable social and economic development . In this context, health promotion activities fall into the mission of OPs, who should be already trained during the course of education.
We tried to estimate a rough cost of the intervention. Regarding the financial impact of WHP programs, an extensive review of the literature  and the major WHP study on cardiac risk factors  showed a positive return on investment, demonstrating that such programs seem to pay for themselves.
The results of this multi-faceted worksite intervention across multiple occupational groups from several geographic communities consistently converged on the evidence of a decreased risk of cardiovascular disease after an educational intervention. The intervention was reasonably priced and easy to integrate within the current organization of occupational health surveillance in Italy.
Authors thank the Italian Ministry of Health and the Veneto Region (Regional Decree n. 2008 of August 3, 2010) that supported the study. The sponsors of the study had no role in study design, data collection, data analysis, data interpretation, or manuscript preparation. The corresponding author took final responsibility for the decision to submit for publication.
Authors wish also to thank Dr. Paolo Civitelli, Ser.T, ULSS 4, Veneto Region, for holding the course “The motivational counseling for prevention of cardiovascular diseases in workplaces”.
- Occupational health. Workplace health promotion. WHO; 2014. [http://www.who.int/occupational_health/topics/workplace/en/]
- Sparling PB. Worksite health promotion: principles, resources, and challenges. Prev Chronic Dis. 2010;7(1) [http://www.cdc.gov/pcd/issues/2010/jan/09_0048.htm]
- Arena R. Lifestyle modification interventions and cardiovascular health: global perspectives on worksite health and wellness and cardiac rehabilitation. Prog Cardiovasc Dis. 2014;56(5):473–5.View ArticlePubMedGoogle Scholar
- Cahalin LP, Myers J, Kaminsky L, Briggs P, Forman DE, Patel MJ, et al. Current trends in reducing cardiovascular risk factors in the United States: focus on worksite health and wellness. Prog Cardiovasc Dis. 2014;56(5):476–83.View ArticlePubMedGoogle Scholar
- Després JP, Alméras N, Gauvin L. Worksite health and wellness programs: Canadian achievements & prospects. Prog Cardiovasc Dis. 2014;56(5):484–92.View ArticlePubMedGoogle Scholar
- Cipriano Jr G, Neves LM, Cipriano GF, Chiappa GR, Borghi-Silva A. Cardiovascular disease prevention and implications for worksite health promotion programs in Brazil. Prog Cardiovasc Dis. 2014;56(5):493–500.View ArticlePubMedGoogle Scholar
- Guazzi M, Faggiano P, Mureddu GF, Faden G, Niebauer J, Temporelli PL. Worksite health and wellness in the European union. Prog Cardiovasc Dis. 2014;56(5):508–14.View ArticlePubMedGoogle Scholar
- Babu AS, Madan K, Veluswamy SK, Mehra R, Maiya AG. Worksite health and wellness programs in India. Prog Cardiovasc Dis. 2014;56(5):501–7.View ArticlePubMedGoogle Scholar
- Okamura T, Sugiyama D, Tanaka T, Dohi S. Worksite wellness for the primary and secondary prevention of cardiovascular disease in Japan: the current delivery system and future directions. Prog Cardiovasc Dis. 2014;56(5):515–21.View ArticlePubMedGoogle Scholar
- Legislative decree 9/4/2008 No. 81 updated in October 2013. [http://www.lavoro.gov.it/SicurezzaLavoro/MS/Normativa/Documents/TU%2081-08%20-%20Ed.%20Dicembre%202014.pdf]
- Deeks A, Lombard C, Michelmore J, Teede H. The effects of gender and age on health related behaviors. BMC Public Health. 2009;30(9):213.View ArticleGoogle Scholar
- Harris AD, McGregor JC, Perencevich EN, Furuno JP, Zhu J, Peterson DE, et al. The Use and Interpretation of quasi-experimental studies in medical informatics. J Am Med Inform Assoc. 2006;13:16–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Conroy RM, Pyörala K, Fitzgerald AP, Sans S, Menotti A, De Backer G, et al. Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE project. Eur Heart J. 2003;24(11):987–1003.View ArticlePubMedGoogle Scholar
- Istituto Superiore della Sanità. Il Progetto Cuore. Carta del Rischio Cardiovascolare. 2013. [http://www.cuore.iss.it/valutazione/carte.asp]
- Robson LS, Shannon HS, Goldenhar LM, Hale AR. Guide to evaluating the effectiveness of strategies for preventing work injuries: how to show whether a safety intervention really works. Cincinnati OH: U.S. Department of Health and Human Services, Public Health Service, Centers for Disease Control and Prevention, National Institute for Occupational Safety and Health; 2001. p. 119.Google Scholar
- Centers for Disease Control. [http://www.cdc.gov/nceh/lead/publications/1997/pdf/c2.pdf]
- Global strategy on diet, physical activity and health. WHO; 2011. [http://www.who.int/dietphysicalactivity/pa/en/]
- Rongen A, Robroek SJ, van Lenthe FJ, Burdorf A. Workplace health promotion: a meta-analysis of effectiveness. Am J Prev Med. 2013;44(4):406–15.View ArticlePubMedGoogle Scholar
- O’Donnell MP. What is the ROI of workplace health promotion? The answer just got simpler by making the question more complicated. Am J Health Promot. 2014;28(6):IV–V.View ArticleGoogle Scholar
- Hawkins NG, Sanson-Fisher RW, Shakeshaft A, D’Este C, Green LW. The multiple baseline design for evaluating population-based research. Am J Prev Med. 2007;33(2):162–8.View ArticlePubMedGoogle Scholar
- Ebrahim S, Taylor F, Ward K, Beswick A, Burke M, Smith GD. Multiple risk factor interventions for primary prevention of coronary heart disease. The Cochrane Library. 2011; Suppl 1 [http://0-onlinelibrary.wiley.com.brum.beds.ac.uk/doi/10.1002/14651858.CD001561.pub3/pdf]
- Milani RV, Lavie CJ. Impact of worksite wellness intervention on cardiac risk factors and one-year health care costs. Am J Cardiol. 2009;104:1389–92.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.