Skip to main content

Do body mass index trajectories affect the risk of type 2 diabetes? A case–control study

Abstract

Background

Although obesity is a well-studied risk factor for diabetes, there remains an interest in whether “increasing body mass index (BMI),” “high BMI per se,” or both are the actual risk factors for diabetes. The present study aimed to retrospectively compare BMI trajectories of individuals with and without diabetes in a case–control design and to assess whether increasing BMI alone would be a risk factor.

Methods

Using comprehensive health check-up data measured over ten years, we conducted a case–control study and graphically drew the trajectories of BMIs among diabetic patients and healthy subjects, based on coefficients in fitted linear mixed-effects models. Patient group was matched with healthy control group at the onset of diabetes with an optimal matching method in a 1:10 ratio. Simple fixed-effects models assessed the differences in increasing BMIs over 10 years between patient and control groups.

Results

At the time of matching, the mean ages in male patients and controls were 59.3 years [standard deviation (SD) = 9.2] and 57.7 years (SD = 11.2), whereas the mean BMIs were 25.0 kg/m2 (SD = 3.1) and 25.2 kg/m2 (SD = 2.9), respectively. In female patients and controls, the mean ages were 61.4 years (SD = 7.9) and 60.1 years (SD = 9.6), whereas the mean BMIs were 24.8 kg/m2 (SD = 3.5) and 24.9 kg/m2 (SD = 3.4), respectively. The simple fixed-effects models detected no statistical significance for the differences of increasing BMIs between patient and control groups in males (P = 0.19) and females (P = 0.67). Sudden increases in BMI were observed in both male and female patients when compared with BMIs 1 year prior to diabetes onset.

Conclusions

The present study suggested that the pace of increasing BMIs is similar between Japanese diabetic patients and healthy individuals. The increasing BMI was not detected to independently affect the onset of type 2 diabetes.

Peer Review reports

Background

The worldwide population of individuals with diabetes is more than 382 million people [1], and potentially, 20.5 million of these are in Japan [2]. Because the subsequent ischemic heart disease, cerebral infarction, retinopathy, nephropathy, and neuropathy are life threatening to numerous people, the prevention of diabetes is of great importance. Among many environmental factors associated with diabetes, obesity is a particularly well-known risk factor for type 2 diabetes [35]. Notably, studies indicate that East Asian individuals with body mass indexes (BMIs) within normal range are much more likely to develop diabetes compared with Caucasians and African Americans with those BMIs [6].

Several studies have indicated that individuals with high BMIs are at a high risk of developing diabetes,[713]; however, due to their cross-sectional or 2–3 time-point cohort designs, these reports have not distinguished the risk of “increasing BMI” on diabetes onsets from that of “high BMI per se.” The present study from Japanese large-scale clinical data was driven by the interest in comparing BMI trajectories in diabetic patients and healthy individuals and by the question of whether increasing BMI alone would cause diabetes.

Methods

Ethics statement

This study was approved by the Ethical Review Board of the Faculty of Medicine, University of Yamanashi (approval number: H22 No.620). Approval was based on the ethical guidelines and regulations of the Declaration of Helsinki. The Japanese guidelines permit the use of data from medical check-ups without consent if the data are anonymous. Hence, informed consent was not required in this study because the data were obtained anonymously from observational check-up data. In addition, in Japan, subjects are not required to approve or disapprove the use of anonymous information from medical check-up data for research. All data were analyzed anonymously, and thorough care was taken to prevent assessed individuals from suffering any disadvantages.

Study subjects

Data were from 64,762 individuals who underwent comprehensive health check-ups during the period from April 1, 1999 to March 31, 2009 at the Yamanashi Koseiren Health Care Center. The criteria for inclusion were undergoing at least two health check-ups over the 5-year period from April 1999 to March 2004. During the period, glycated hemoglobin A1c (HbA1c) and fasting plasma glucose (FPG) levels of the study subjects should not meet the Japanese diagnostic criteria for prediabetes of HbA1c ≥ 6.5 % or FPG ≥ 126 mg/dL (7.0 mmol/L) [14]. In addition to the criteria, the subjects were eligible if they underwent at least one check-up during the period from April 2004 to March 2009.

Patient and control definitions

Patients were defined as subjects who were newly diagnosed with diabetes at other hospitals outside or who exhibited HbA1c ≥ 6.5 % or FPG ≥ 126 mg/dL during the period from April 2004 to March 2009. Control healthy subjects, whose HbA1c and FPG levels did not exceed the diagnostic criteria for diabetes within the 10-year survey period, were matched to diabetic patients for BMI and age at a 1:10 ratio at the time of diabetes onset of patients (time point 0). These operations were all performed separately for males and females. The diagnostic criteria were also based on the guidelines for epidemiologic research on diabetes by “Committee Report Regarding Diabetes Classifications and Diagnostic Criteria (internationally standardized version)” [15]. Matching was performed using an optimal matching [16].

Definitions of time points

The year of diabetes onset in the patient group was set as time point 0. Each year was then measured with reference to time point 0. Thus, 1 year before onset was time point −1, and 2 years before onset was time point −2. This was continued to time point −9. Accordingly, time points are expressed in years.

Clinical parameters

Clinical parameters included in our analysis comprised sex, date of birth, height, weight, HbA1c, and FPG levels. Age was calculated from date of birth. BMI was calculated as weight in kilograms divided by square of height in meters.

Statistical analysis

We calculated mean and standard deviation (SD) values for age, weight, height, BMI, HbA1c, and FPG levels for each group at the matched time point. Student’s t test compared baseline characteristics between patient and control groups. As the major focus of this study, BMI trajectory was subjected to regression analysis, using a linear mixed-effects model that included time points, group categories, and interaction terms as explanatory categorical variables. The group categories included diabetic patients and healthy subjects, and the interaction terms were defined between time points and group categories. The past histories of the patient and control groups were reviewed, and mean BMIs were estimated using the linear mixed-effects model, from which a trajectory of each group was drawn, stratified for sex. With the aim of taking individual variability into consideration, a random effect for individuals was set for the intercept of the linear mixed-effects models. To assess whether increasing BMI affects the risk for diabetes, P value for the difference in increasing BMIs between patient and control groups was evaluated in another simple fixed-effects model with interaction terms that categorizes patient or control group and time points 0 or −9. All statistical analyses were performed using SAS version 9.3 (Cary, NC, USA). All reporteP values are two-sided; P values of <0.05 were considered to be statistically significant.

Results

Figure 1 shows how eligible and matched control subjects were selected. Table 1 represents the baseline characteristics in patient and control groups. The mean ages of males at the matching time point 0 were 59.3 (SD = 9.2) years in the patient group and 57.7 (SD = 11.2) years in the control group (P = 0.003), whereas the mean BMIs were25.0 (SD = 3.1) kg/m2 in the patient group and 25.2 (SD = 2.9) kg/m2 in the control group (P = 0.22), respectively. In females, the mean ages were 61.4 (SD = 7.9) years in the patient group and 60.1 (SD = 9.6) years in the control group (P = 0.06), whereas the mean BMIs were 24.8 (SD = 3.5) kg/m2 in the patient group and 24.9 (SD = 3.4) kg/m2 in the control group (P = 0.67), respectively.

Fig. 1
figure 1

Flowchart of study individuals

Table 1 Characteristics of diabetic patients and healthy subjects at the matching time of diabetes onset

Table 2 and Figs. 2 and 3 show the trajectories in patient and control groups. In males, the mean BMI was consistently higher in the control group than in the patient group. In females, BMI also tended to be higher in the control group than in the patient group. In contrast to males, a lesser increase in BMI was observed over the nine years in the female patient group than in the female control group. Sudden increases in BMI were observed in both male and female patients from one year before diabetes onset. Statistically significant differences in BMI as interaction terms were observed between the patient and control groups at −1 (P = 0.01), −2 (P = 0.01), −3 (P = 0.045), and −4 (P = 0.0002) time points among males and at −4 (P = 0.03), and −8 (P = 0.04) time points among females. The simple fixed-effects models, which examined the major question of this study, represented P of 0.19 in males and 0.67 in females between the increasing BMIs in patient and control groups.

Table 2 Estimated mean BMIs before onsets of diabetes based on coefficients in mixed-effects models in patients and matched healthy subjects
Fig. 2
figure 2

Retrospective BMI trajectories of patient and control groups by a linear mixed-effects model: Males

Fig. 3
figure 3

Retrospective BMI trajectories of patient and control groups by a linear mixed-effects model: Females

Discussion

The results of this study suggest that the retrospective trajectories of BMIs in patients with newly-acquired diabetes are similar to those in healthy individuals. To the best of our knowledge, this is the first study comparing BMI trajectories between Japanese individuals with and without diabetes using repeated data.

Vistisen et al. described large changes in BMI in patients with diabetes in the UK without using matching case to control [17], although typical BMI for the onset of type 2 diabetes greatly differed between individuals of East Asian descent and individuals of European descent [18]. Therefore, the present data may be specific to a type of Japanese diabetic patients and hence may result in limited generalizability. No previous reports described sudden increases in BMI prior to diabetes onset. One Japanese study, although performed over two time points, reported that the risk for diabetes onset increased when either BMI or body weight increased [19]. In a previous study, Kodama et al. demonstrated that individuals of East Asian descent with non-overweight or slightly overweight BMIs also have a high risk of developing type 2 diabetes [6]. These studies are consistent with our results presenting increasing but almost normal-ranged BMIs among the diabetic patients.

The most likely mechanism for the diabetic onsets in patient group is insulin resistance and impaired insulin secretion in Japanese individuals with moderately high and increasing BMI. According to a previous systematic review, individuals of East Asian descent, which includes Japanese people, have higher insulin sensitivity than people from Western countries, whereas the insulin secretory ability is low only in East Asian individuals with diabetes. It has been reported that when the body weight of an individual with impaired insulin secretion increases, the risk of developing type 2 diabetes also increases [6]. It appears that this mechanism is involved in diabetes onset in patients in the present study. In another study that identified a polymorphism involved in the onset and progression of type 2 diabetes, the KCNJ15 polymorphism was observed in many non-overweight diabetic patients of Japan, whose ability to secrete insulin tended to decrease over time [20]. Although the patients and healthy subjects in this study had not undergone genetic testing, increased body weight in such individuals with decreased insulin secretory ability may have increased insulin resistance [21, 22] and it may have increased their susceptibility to developing diabetes. Another pathway of diabetic onset may be attributed to the observed fluctuation in BMI prior to the onset of diabetes in this study. Because literature suggests that weight fluctuation increases the risk of type 2 diabetes in American women [23], the observed fluctuation among the patient group may have affected their insulin resistance [24]. Alternatively, the observed fluctuation may reflect the small sample sizes of the male and female patient groups.

Literature is available on the sensitivity and specificity for an HbA1c threshold of 6.5 %. Ito reports a sensitivity of 53.4 % and a specificity of 94.6 % for an HbA1c threshold of 6.5 % in a Japanese population [25]. In a Chinese population, Bao et al. report a sensitivity of 53.7 % and a specificity of 97.4 % for an HbA1c threshold of 6.5 % [26]. Because Choi et al. report a sensitivity of 52 % and a specificity of 97 % for an HbA1c threshold of 6.2 % in a Korean population, the sensitivity for an HbA1c threshold of 6.5 % should be lower than 52 % [27]. Furthermore, the latest data from the American National Health and Nutrition Examination Survey describes a sensitivity of approximately 44 % and a specificity of approximately 98 % at an HbA1c threshold of 6.5 % [28]. On the other hand, data from a Japanese population has reported a sensitivity of 86.5 % and a specificity of 87.3 % at the FPG threshold of 115 mg/dL (6.4 mmol/L) [29]. Thus, the present epidemiologic-research-based diagnosis for diabetes by the combination of HbA1c ≥ 6.5 % or FPG ≥ 126 mg/dL [15] would have achieved an acceptable sensitivity and specificity.

The strengths of this study include the fact that BMI trajectories of both patients and controls were matched for the timing of diabetes onset, the BMI, and the age. Thus, this study is considered to have excluded the effect of having high BMI per se on the risk of developing diabetes. Additionally, we used chronological data in which BMIs were measured at three or more time points. Although subjects were not identical at each time point, it is noteworthy that the data depicted the trajectories over 10 time points.

A limitation of this study is that, in addition to the interviews with the patients regarding their newly-acquired diabetes, diabetes was partly diagnosed based on HbA1c and FPG levels rather than on a confirmed diagnosis from secondary tests or glucose tolerance tests. Although we used criteria equivalent to parts of the diagnostic criteria used by physicians in Japan [14], and guidelines for epidemiologic research on diabetes allow for the use of the criteria applied in the present study [15], it is possible that some subjects were incorrectly classified as patient group subjects. Results in Tables 1 and 2 show older ages and larger weight gains from the previous year in the patient group than the control group. Because people who are either older or are gaining weight are more likely to undergo health check-ups and be diagnosed with diabetes [30], the misclassification may have been differential and hence may have produced an associated bias toward broadening the difference in BMI trajectories between patient and control groups. The second limitation of our study was that all subjects underwent a paid comprehensive health check-up. Therefore, these results may not be directly applicable to the entire Japanese population. Another limitation was the fact that because this study used retrospective data from the time of disease onset (time point 0), the population was not identical at each observational time point. However, this had probably only little effect, as no marked generational changes were seen over the period from 1999 to 2009 in attitudes toward health, diet, and exercise habits of Japanese people. Furthermore, although adjustment by multivariate regression for confounding factors could be employed in this study, we have exploited simple illustrations to compare trajectories, matching patients and healthy subjects with confounding factors.

This study suggests that the body weight of some patients with type 2 diabetes increases immediately before disease onset. Therefore, healthy but pre-diabetic individuals whose body weight has been increasing for the past few years should receive proactive guidance regarding diet and exercise from their physicians, and they should be encouraged to control their body weight. The gradually increasing mean BMI of individuals with and without diabetes in our study reached approximately 25 at the final measurement. This suggests that even if BMI is within normal range (18.5–25), physicians should be recommended to strongly advise East Asian adults not to gain any weight.

Conclusions

The present study from large-scale clinical data of Japanese residents suggests that the pace of increasing BMI in diabetic patients is similar to that in healthy individuals. The results also indicate that the BMIs of the Japanese patients do not markedly increase prior to the onset of type 2 diabetes.

Abbreviations

BMI:

Body mass index

HbA1c:

Glycated hemoglobin A1c

FPG:

Fasting plasma glucose

SD:

Standard deviation

References

  1. International Diabetes Federation. Diabetes Atlas 6th Edition. 2013. Available: http://www.idf.org/diabetesatlas. Accessed 29 June 2015.

  2. Headquarters of Cancer Control, Ministry of Health, Labour and Welfare. Outline of the results of the 2012 National Health and Nutrition Survey. 2013. Available: http://www.mhlw.go.jp/file/04-Houdouhappyou-10904750-Kenkoukyoku-Gantaisakukenkouzoushinka/0000032813.pdf. Japanese. Accessed 29 June 2015.

    Google Scholar 

  3. Mokdad AH, Ford ES, Bowman BA, Dietz WH, Vinicor F, Bales VS, et al. Prevalence of obesity, diabetes, and obesity-related health risk factors, 2001. JAMA. 2003;289:76–9.

    Article  PubMed  Google Scholar 

  4. Helmrich SP, Ragland DR, Leung RW, Paffenbarger Jr RS. Physical activity and reduced occurrence of non-insulin-dependent diabetes mellitus. New Eng J Med. 1991;325:147–52.

    Article  CAS  PubMed  Google Scholar 

  5. Nguyen NT, Nguyen XM, Lane J, Wang P. Relationship between obesity and diabetes in a US adult population: findings from the National Health and Nutrition Examination Survey, 1999–2006. Obes Surg. 2011;21:351–5.

    Article  PubMed  Google Scholar 

  6. Kodama K, Tojjar D, Yamada S, Toda K, Patel CJ, Butte AJ. Ethnic differences in the relationship between insulin sensitivity and insulin response. A systematic review and meta-analysis. Diabetes Care. 2013;36:1789–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Ma RC, Chan JC. Type 2 diabetes in East Asians: similarities and differences with populations in Europe and the United States. Ann NY Acad Sci. 2013;1281:64–91.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Tirosh A, Shai I, Afek A, Dubnov-Raz G, Ayalon N, Gordon B, et al. Adolescent BMI trajectory and risk of diabetes versus coronary disease. New Eng J Med. 2011;364:1315–25.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Hwang LC, Chen CJ, Lin BJ. Obesity and changes in body weight related to 10-year diabetes incidence in women in Taiwan (1993–2003). Asia Pac J Clin Nutr. 2007;16:677–82.

    CAS  PubMed  Google Scholar 

  10. Ishikawa-Takata K, Ohta T, Moritaki K, Gotou T, Inoue S. Obesity, weight change and risks for hypertension, diabetes and hypercholesterolemia in Japanese men. Eur J Clin Nutr. 2002;56:601–7.

    Article  CAS  PubMed  Google Scholar 

  11. Mishra GD, Carrigan G, Brown WJ, Barnett AG, Dobson AJ. Short-term weight change and the incidence of diabetes in midlife - Results from the Australian Longitudinal Study on Women’s Health. Diabetes Care. 2007;30:1418–24.

    Article  PubMed  Google Scholar 

  12. Oguma Y, Sesso HD, Paffenbarger RS, Lee IM. Weight change and risk of developing type 2 diabetes. Obes Res. 2005;13:945–51.

    Article  PubMed  Google Scholar 

  13. Koh-Banerjee P, Wang YF, Hu FB, Spiegelman D, Willett WC, Rimm EB. Changes in body weight and body fat distribution as risk factors for clinical diabetes in US men. Am J Epidemiol. 2004;159:1150–9.

    Article  PubMed  Google Scholar 

  14. Kuzuya T, Nakagawa S, Satoh J, Kanazawa Y, Iwamoto Y, Kobayashi M, et al. Report of the Committee on the classification and diagnostic criteria of diabetes mellitus. Diabet Res Clin Pract. 2002;55:65–85.

    Article  Google Scholar 

  15. The Committee of the Japan Diabetes Society on the diagnostic criteria of diabetes mellitus. Committee Report Regarding Diabetes Classifications and Diagnostic Criteria (internationally standardized version). J Jpn Diabet Soc. 2012;55:485–504. Japanese.

    Google Scholar 

  16. Rosenbaum PR. Optimal matching for observational studies. J Am Stat Assoc. 1989;84:1024–32.

    Article  Google Scholar 

  17. Vistisen D, Witte DR, Tabák AG, Herder C, Brunner EJ, Kivimäki M, et al. Patterns of obesity development before the diagnosis of type 2 diabetes: The Whitehall II cohort study. PLoS Med. 2014;11, e1001602.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Sone H, Yoshimura Y, Ito H, Ohashi Y, Yamada N. Energy intake and obesity in Japanese patients with type 2 diabetes. Lancet. 2004;363:248–9.

    Article  PubMed  Google Scholar 

  19. Nanri A, Mizoue T, Takahashi Y, Matsushita Y, Noda M, Inoue M, et al. Association of weight change in different periods of adulthood with risk of type 2 diabetes in Japanese men and women: the Japan Public Health Center-Based Prospective Study. J Epidemiol Commun H. 2011;65:1104–10.

    Article  Google Scholar 

  20. Okamoto K, Iwasaki N, Nishimura C, Doi K, Noiri E, Nakamura S, et al. Identification of < i > KCNJ15</i > as a Susceptibility Gene in Asian Patients with Type 2 Diabetes Mellitus. Am J Hum Genet. 2010;86:54–64.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Okuno A, Tamemoto H, Tobe K, Ueki K, Mori Y, Iwamoto K, et al. Troglitazone increases the number of small adipocytes without the change of white adipose tissue mass in obese Zucker rats. J Clin Invest. 1998;101:1354.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Yamauchi T, Kamon J, Waki H, Murakami K, Motojima K, Komeda K, et al. The mechanisms by which both heterozygous peroxisome proliferator-activated receptor γ (PPARγ) deficiency and PPARγ agonist improve insulin resistance. J Biol Chem. 2001;276:41245–54.

    Article  CAS  PubMed  Google Scholar 

  23. Field AE, Manson JE, Laird N, Williamson DF, Willett WC, Colditz GA. Weight cycling and the risk of developing type 2 diabetes among adult women in the United States. Obes Res. 2004;12:267–74.

    Article  PubMed  Google Scholar 

  24. Kahn BB, Flier JS. Obesity and insulin resistance. J Clin Invest. 2000;106:473.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Ito C. Evidence for diabetes mellitus criteria in 2010 using HbA1c. Diabetol Int. 2013;4:9–15.

    Article  Google Scholar 

  26. Bao Y, Ma X, Li H, Zhou M, Hu C, Wu H, et al. Glycated haemoglobin A1c for diagnosing diabetes in Chinese population: cross sectional epidemiological survey. BMJ. 2010;340.

  27. Choi SH, Kim TH, Lim S, Park KS, Jang HC, Cho NH. Hemoglobin A1c as a Diagnostic Tool for Diabetes Screening and New-Onset Diabetes Prediction A 6-year community-based prospective study. Diabetes Care. 2011;34:944–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Buell C, Kermah D, Davidson MB. Utility of A1C for diabetes screening in the 1999–2004 NHANES population. Diabetes Care. 2007;30:2233–5.

    Article  CAS  PubMed  Google Scholar 

  29. Miyazaki M, Kubo M, Kiyohara Y, Okubo K, Nakamura H, Fujisawa K, et al. Comparison of diagnostic methods for diabetes mellitus based on prevalence of retinopathy in a Japanese population: the Hisayama Study. Diabetologia. 2004;47:1411–5.

    Article  CAS  PubMed  Google Scholar 

  30. Liang W, Shediac-Rizkallah MC, Celentano DD, Rohde C. A population-based study of age and gender differences in patterns of health-related behaviors. Am J Prev Med. 1999;17:8–17.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This study was supported by a Grant-in-Aid for Young Scientists (B) (20596879) and Research (C) (15 K08730) from the Ministry of Education, Culture, Sports, Science and Technology of Japan to HY. On completion of this study, we would like to extend our deepest gratitude to all the staff members at the Yamanashi Koseiren Health Care Center who readily provided data for our study.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Hiroshi Yokomichi or Zentaro Yamagata.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ZY conceived the study. KS designed the study. YY and MT carried out the laboratory investigations. YM and HY conducted the statistical analyses. AT, MS, and RS provided input on analytical strategy. YM and HY drafted the manuscript with assistance from ZY. MM and SM agreed with results, discussion, and conclusion and contributed to the writing of the manuscript. ZY is guarantor of the manuscript. All authors contributed to drafts of the manuscript and have read and approved the final manuscript.

Authors’ information

Yoshihiko Mano and Hiroshi Yokomichi are the co-first authors.

Yoshihiko Mano and Hiroshi Yokomichi contributed equally to this work.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mano, Y., Yokomichi, H., Suzuki, K. et al. Do body mass index trajectories affect the risk of type 2 diabetes? A case–control study. BMC Public Health 15, 718 (2015). https://doi.org/10.1186/s12889-015-2073-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-015-2073-y

Keywords