Demographic determinants of acute gastrointestinal illness in Canada: a population study
© Majowicz et al; licensee BioMed Central Ltd. 2007
Received: 26 January 2007
Accepted: 18 July 2007
Published: 18 July 2007
Gastrointestinal illness is an important global public health issue, even in developed countries, where the morbidity and economic impact are significant. Our objective was to evaluate the demographic determinants of acute gastrointestinal illness in Canadians.
We used data from two population-based studies conducted in select communities between 2001 and 2003. Together, the studies comprised 8,108 randomly selected respondents; proxies were used for all respondents under 12 years and for respondents under 19 years at the discretion of the parent or guardian. Using univariate and multivariate logistic regression, we evaluated the following demographic determinants: age, gender, cultural group, and urban/rural status of the respondent, highest education level of the respondent or proxy, number of people in the household, and total annual household income. Two-way interaction terms were included in the multivariate analyses. The final multivariate model included income, age, gender, and the interaction between income and gender.
After adjusting for income, gender, and their interaction, children under 10 years had the highest risk of acute gastrointestinal illness, followed by young adults aged 20 to 24 years. For males, the risk of acute gastrointestinal illness was similar across all income levels, but for females the risk was much higher in the lowest income category. Specifically, in those with total annual household incomes of less than $20,000, the odds of acute gastrointestinal illness were 2.46 times higher in females than in males.
Understanding the demographic determinants of acute gastrointestinal illness is essential in order to identify vulnerable groups to which intervention and prevention efforts can be targeted.
Gastrointestinal illness (GI) remains an important global public health issue [1, 2]. In developed countries, although GI tends to be self-limiting and mild, the associated morbidity and economic impact are significant [3–5]. To address this, numerous countries have estimated the incidence and burden of GI in the community via population-based studies [6–12], including two studies conducted recently in Canadian communities [13, 14]. These population-based studies collect information on gastrointestinal symptoms experienced by a random sample of study area residents, the severity of those symptoms, and demographic and other information. Since the main purpose of these studies is to estimate disease incidence and burden, demographic determinants of illness are often investigated as a secondary objective only, and are not necessarily thoroughly explored.
Understanding the relationships between GI and determinants of health in the general population, however, is essential to identify vulnerable groups to which intervention and prevention efforts can be targeted. Therefore, the objective of this study was to investigate the demographic determinants of acute gastrointestinal illness (AGI) in Canadians using available data from population-based studies.
At the time of this analysis, the Public Health Agency of Canada (formerly Health Canada) had conducted two studies designed to ascertain the burden and distribution of self-reported AGI in defined Canadian populations. One study was conducted in Hamilton, Ontario, Canada from February 2001 to February 2002, and one was conducted in three communities in the province of British Columbia from June 2002 to June 2003; these studies used the same methodology and core survey tool, and have been described in detail elsewhere [13, 14]. Briefly, both were retrospective, cross-sectional telephone surveys administered to randomly selected residents of the study area. A random sample of households was chosen using a commercial database of residential telephone numbers. One individual from each household was then randomly selected to participate by identifying the individual in the household with the next birthday. Proxy respondents were used for all individuals under 12 years of age, and for individuals 12 to 18 years of age at the discretion of the parent or guardian.
Respondents were asked whether they had experienced any vomiting or diarrhea in the 28 days prior to the interview. Cases were those respondents who reported vomiting or diarrhea in the four weeks prior to the interview, excluding those who reported that their vomiting or diarrhea was due to a chronic condition including pregnancy, medication use, colitis, diverticulitis, Crohn's disease, irritable bowel syndrome, or other chronic condition. Respondents who did not report vomiting or diarrhea, as well as those whose symptoms were due to chronic conditions, were included in the non-case group. A broad case definition for AGI was deliberately chosen to ensure high sensitivity and case capture.
Ethical approval for these studies was obtained from one or more of the following boards: the Research Ethics Board of St. Joseph's Hospital (Hamilton, Ontario, Canada), McMaster University (Hamilton, Ontario, Canada), the Human Subjects Committee of the University of Guelph (Guelph, Ontario, Canada), and the University of British Columbia Behavioural Research Ethics Board (Vancouver, British Columbia, Canada). The response rates for the surveys were 36.6%  and 44.3% , which are in the range of response rates for similar studies [7, 9, 15, 16]. Together, the studies comprised 8,108 respondents (Hamilton, 3,496; British Columbia, 4,612).
Frequency (n), percent (%), odds ratios (OR), 95% confidence intervals (CI) for odds ratios and P-values from the univariate analysis of the relationship between demographic determinants and acute gastrointestinal illness in randomly selected residents of two Canadian study areas, 2001 to 2003 (N = 8,108)
CI for OR
Number of People in Household
Three or more b
No high school diploma b
High school diploma
College/trade school diploma
University, graduate, or professional diploma
North American b
Native North American/Aboriginal
>$20 000 to <$40 000
>$40 000 to <$60 000
>$60 000 to <$80 000
Logistic regression was used to determine how the risk of AGI related to demographic variables. To examine whether the relationship between each demographic factor and the risk of AGI varied between the two study areas, we fit (separately for each demographic factor) multivariate models with the demographic factor, study area, and their interaction as independent variables (results not shown). However, no significant interactions were found and data from the two studies were combined for all further analyses. Univariate models were then fit for each demographic variable (Table 1). To explore various multivariate models and ultimately choose a final model, all demographic variables and their two-way interactions were used in stepwise selection, with the criteria for entry and exit into the model being that the P-value of the score chi-square test be less than 0.05. Individuals with missing data for a given variable were excluded from any models in which that variable was present. To assess whether any variables in the final model were subject to confounding by any variables that had been omitted from the final model, each omitted variable was re-introduced individually (results not shown). The impact on the sign, magnitude, and significance of each of the original coefficients was examined; a change from significant to non-significant (or vice versa) at P = 0.05, or a change in the resulting odds ratio of ± 0.5, was considered biologically significant enough to retain the variable in the final model as a confounder.
All statistical analyses were carried out in SAS version 9.1 (SAS Institute Inc., Cary, NC, USA, 2002–2003). Likelihood ratios were used to compare models and Wald tests were used for tests of global hypotheses and tests involving individual parameters .
Results of the univariate analysis are shown (Table 1). The risk of AGI was significantly associated with the number of people in the household (P = 0.005), culture (P =< 0.001), income (P = 0.033), age (P < 0.001), and gender (P = 0.001), but not with urban/rural status (P = 0.080) or education (P = 0.150). Specifically, the risk of AGI in the past four weeks increased significantly as the number of people in the household increased.
The odds of AGI for respondents who identified themselves as Asian were only 0.37 times higher (i.e. 2.69 times lower) than respondents who identified themselves as North American (P < 0.001). Respondents with household incomes between $40,000 and $60,000 had odds of AGI that were 0.76 times higher (i.e. 1.32 times lower) than respondents with household incomes less than $20,000 (P = 0.047). A significantly higher risk of AGI was observed in children less than 10 years (P < 0.001), and young adults 20 to 24 years (P = 0.003), compared to those 25 to 64 years. The odds of AGI in females were 1.28 times higher than males (P = 0.001).
Adjusted odds ratios (OR), 95% confidence interval (CI) for odds ratios and P-values from the final multivariate model of the relationship between demographic determinants and acute gastrointestinal illness in randomly selected residents of two Canadian study areas, 2001 to 2003 (N = 5,732)
CI for OR
>$20 000 to <$40 000
>$40 000 to <$60 000
>$60 000 to <$80 000
>$20 000 to <$40 000 (Female)
>$40 000 to <$60 000 (Female)
>$60 000 to <$80 000 (Female)
>$80 000 (Female)
It is interesting to compare the odds of AGI in males and females with the same income. For instance, the odds ratio for females earning $60,000–$80,000 versus males with the same income was 1.38 (1.74/1.26). This difference is not significant (P = 0.107). The difference between the odds of AGI in males and females was statistically significant only at the lowest income level (P < 0.001). In those with total annual household incomes of less than $20,000, the odds of AGI were 2.46 times higher in females than in males.
In Figure 2, the differential effect of income on the risk of AGI for males and females is clearly visible. For males, the risk of AGI was approximately constant across all income levels. For females, however, the risk was much higher in those with total household incomes less than $20,000.
This study investigated the demographic determinants of AGI in select Canadian-based populations. Adjusting for gender and income, children under 10 years had the highest risk of AGI, followed by those aged 20 to 24 years. Adjusting for age, in households with annual incomes under $20,000, females were significantly more likely to have AGI than males. These results highlight that, although AGI is typically self-limiting and is rarely a cause of mortality in developed countries, it still represents a particular health risk in specific sub-populations.
Interestingly, there was no interaction between the place in which the study occurred and each of the demographic variables assessed. Thus, it is possible that the associations observed here are consistent across two distinct geographic areas of Canada (southern Ontario, British Columbia), such that the increased risk in low income females and children may be due to underlying factors unrelated to geography. Additionally, income, age, and gender were significantly associated with the risk of AGI whether or not the following variables were controlled for: total number of people in the household, the urban/rural status of the respondent, education, and cultural group, suggesting that the observed associations are not confounded by these variables.
In the univariate analysis, we observed a significantly higher risk of AGI in children under 10 years, and young adults between 20 and 24 years, as compared to adults aged 25 to 64 years. As noted above, even when income, gender, and the income-gender interaction were accounted for, children under 10 years and young adults between 20 and 24 years remained at a higher risk of AGI. In children, this increased risk likely reflects an increased susceptibility due to immune status. In young adults, this increased risk may reflect behavioural factors. Our findings are somewhat consistent with one study from the Netherlands, which reported highest incidences in children and the elderly , and with other international studies which report higher rates in children [8, 9, 18, 19]. Interestingly, we did not observe an increased risk of AGI in the elderly, as has been observed elsewhere  and which is biologically plausible due to decreasing immune function. Although individuals over the age of 65 years were well represented in our data (n = 1,346), it is likely that the sampling method inadvertently selected for healthier individuals in this age category, since residents of institutions (including nursing homes) were not included in the sampling frame. Thus, despite the lack of observed risk in the elderly in this study, it is possible that elderly Canadians may be at increased risk for AGI.
Past studies report higher rates of AGI in females than males [6, 13, 14, 20, 21], which may be due to an increased exposure to infectious causes of AGI via foodborne or person-to-person transmission. Given that the kitchen serves as a reservoir for many foodborne causes of AGI , the higher risk observed in females may reflect the fact that, in Canada, females generally do more food preparation than males and thereby incur greater exposure to pathogens. This is supported by a study from England and Wales which found female gender a significant risk for infection with Campylobacter jejuni infection , which is primarily a foodborne pathogen . Additionally, although the relationship between gender and the risk of AGI observed here was unconfounded by the number of people in the household, we were unable to account for the ages of any household residents other than the respondent. In a recent Australian study , women between the ages of 25 and 64 years both with and without AGI were compared; the study found that 18% of women with AGI had at least one child under five years of age in their household, compared to 5% of women without AGI. Thus, it is possible that females are at an increased risk due to the presence of and interaction with young children within the household.
It is possible that the higher rate of AGI observed in females may be due to recall bias (where females may be more likely to take note of and thus report AGI than males), or to reporting bias (where females may be more comfortable discussing and thus more likely to report AGI than males). However, such explanations are unlikely given the plausibility for greater exposure to infectious causes of AGI in females. In addition, if reporting bias was the reason for the higher rate of AGI observed in females, we would expect higher rates in females consistently across age groups and income levels, which was not observed here.
There is a growing body of literature linking health to income or income inequality, with lower household income associated with an increased risk of morbidity and mortality [25–28]. In Canada, income and good health are positively associated [26, 29], and low socio-economic status is associated with poor health , despite the existence of universally insured health services. However, the specific association between income and the risk of AGI has not been previously explored.
Here, we found that total annual household income was associated with the risk of AGI in females, but not males. In males, the risk of AGI was consistent regardless of income. In females, a higher risk occurred in those in the lowest income category (total annual household income of less than $20,000). In this category, the odds of AGI for females were 2.5 times higher than the odds for males, regardless of age. Unfortunately, no literature exists which provides explicit reasons for this observation. Specific hypotheses may include different occupational risk settings in low income males versus females, or an increased susceptibility in females (due perhaps to increased foodborne exposure or increased exposure to infected children) that is exacerbated by low income living conditions. Further research evaluating reasons for this apparent increased risk in low-income females is warranted. Regardless of the cause, this finding calls for targeting information and interventions to this segment of the population, potentially via local public health outreach programs.
The results of the multivariate analysis can be interpreted to yield stratum-specific odds ratios; for example, the odds of AGI in a female child under 10 years of age who resides in a low income household (i.e. total annual household income of less than $20,000) is 0.554, and the odds of AGI in a male child in the same age and income category is 0.225, yielding an odds ratio of 2.46. Examining the risk of AGI in females across income categories showed that those in households with total annual incomes of $20,000 to $40,000 (OR = 0.566), $40,000 to $60,000 (OR = 0.530), $60,000 to $80,000 (OR = 0.706), and over $80,000 (OR = 0.424) all had a lower risk of AGI than females in low income households.
Low response rate in the original studies (36.6% in Hamilton and 44.3% in B.C) was the main limitation of this study, and is a limitation typical of such telephone surveys. Other similar studies report response rates ranging from 27% to 71% [6, 7, 9, 15, 16] Non-response in our data likely related to the subject-matter of the questionnaire. If the nature of the relationship between AGI and demographic determinants in respondents is different than in the non-respondents, then the results presented here will be biased. Future studies should attempt to minimize non-response to mitigate this potential source of bias.
Item non-response is also a concern in surveys. Here, urban/rural status was missing for 9% of respondents, and total annual household income was missing for 25% of respondents. Thus our final multivariate model, which included income, used data from 5,732 respondents. This may be of concern if the likelihood of responding to specific questions differed across income levels. For example, if non-response to the income question is greater in those with lower incomes, and low income is a risk factor for AGI, the odds ratio for low versus high income will be a biased underestimate of the true risk of illness. To address this, we also examined models without income (results not shown). Since such analyses yielded consistent conclusions with respect to the other demographic variables, the impact of income non-response on the results presented here is likely minor.
In our study, the education variable measured the highest level of education attained by either the respondent (for those over 18 years) or their proxy (for those 18 years and younger), in an attempt to capture the education of the person guiding the behaviour of the respondent. Thus, for those 18 years and under, we assumed that the proxy (parent or legal guardian) was more instrumental in guiding the behaviour of the respondent than the respondent themselves. Although this is likely true for very young respondents, it may be less true as the age of the respondent approaches 18 years. Thus for teenage respondents, we may have over-estimated the education level of the person guiding the behaviour of the respondent. However, in our analysis, we found no interaction between education and age, suggesting that the lack of association between education and the risk of AGI observed here is the same regardless of the age of the respondent. In any case, our findings do not negate the need for future in-depth analyses of specific components of education (e.g. food handling or hygiene training) that may decrease the risk of AGI.
Cases of AGI in this analysis were those who reported vomiting or diarrhea in the past four weeks, excluding those whose symptoms were due to a chronic condition. In these data, no attempt was made to differentiate infectious AGI from other causes such as food allergies or intolerances, over-indulgence of drugs or alcohol, or other causes of AGI. Although respondents reported what they believed to be the cause of their illness, this information was not used to exclude cases since the validity of these self-diagnosed causes was highly variable. Thus our case definition of AGI, although highly sensitive for infectious GI (i.e. includes most true cases of infectious GI), should not be considered specific for infectious GI (i.e. includes some non-infectious cases of AGI).
Another possible limitation of this study is that the data were collected via telephone interview. Thus, the results presented here may not be applicable to those without telephones, such as the homeless, those in institutions, or who are incarcerated. Lastly, it is debated whether analyses such as these should be weighted by the number of persons in the household . We repeated our analyses (not shown) with weighting and found no difference in the substantive conclusions.
In Canada, it appears that children less than 10 years, young adults 20 to 24 years, and females in households with annual incomes under $20,000 are at an increased risk for AGI. In children, this increased risk may reflect an increased susceptibility to gastrointestinal infections due to immune status, and in young adults, this increased risk may be due to behavioural factors. In low income females, however, the specific reasons for this increased risk are unclear, and further research is needed. Understanding these relationships between AGI and determinants of health in the population is necessary to guide intervention and prevention efforts. These results suggest that children, young adults, and low income females should be targeted by public health programs aimed at decreasing the incidence of AGI in Canada.
The authors thank the members of Hamilton Population Study research team and the British Columbia Population Study research team for their earlier work on the studies from which the data for this analysis were obtained. Funding for this research was provided by the Public Health Agency of Canada and the Natural Sciences and Engineering Research Council of Canada (NSERC).
- Bern C, Martines J, de Zoysa I, Glass RI: The magnitude of the global problem of diarrhoeal disease: a ten-year update. Bulletin of the World Health Organization. 1992, 70 (6): 705-714.PubMed CentralPubMed
- Kosek M, Bern C, Guerrant RL: The global burden of diarrhoeal disease, as estimated from studies published between 1992 and 2000. Bulletin of the World Health Organization. 2003, 81 (3): 197-204.PubMed CentralPubMed
- Mead PS, Slutsker L, Dietz V, McCaig LF, Bresee JS, Shapiro C, Griffin PM, Tauxe RV: Food-related illness and death in the United States. Emerging Infectious Diseases. 1999, 5 (5): 607-625.PubMed CentralView ArticlePubMed
- Helms M, Vastrup P, Gerner-Smidt P, Molbak K: Short and long term mortality associated with foodborne bacterial gastrointestinal infections: registry based study. British Medical Journal. 2003, 326 (7385): 357-10.1136/bmj.326.7385.357.PubMed CentralView ArticlePubMed
- Majowicz SE, McNab WB, Sockett PN, Henson S, Doré K, Edge VL, Buffett MC, Fazil A, Read S, McEwen SA, Stacey D, Wilson JB: Burden and cost of gastroenteritis in a Canadian community. Journal of Food Protection. 2006, 69 (3): 651-659.PubMed
- Hoogenboom-Verdegaal A, de Jong JC, During M, Hoogenveen R, Hoekstra JA: Community-based study of the incidence of gastrointestinal diseases in The Netherlands. Epidemiology and Infection. 1994, 112 (3): 481-487.PubMed CentralView ArticlePubMed
- Wheeler JG, Sethi D, Cowden JM, Wall PG, Rodrigues LC, Tompkins DS, Hudson MJ, Roderick PJ: Study of infectious intestinal disease in England: rates in the community, presenting to general practice, and reported to national surveillance. British Medical Journal. 1999, 318 (7190): 1046-1050.PubMed CentralView ArticlePubMed
- de Wit MA, Koopmans MP, Kortbeek LM, Wannet WJ, Vinje J, van Leusden F, Bartelds AI, van Duynhoven YT: Sensor, a population-based cohort study on gastroenteritis in the Netherlands, incidence and etiology. American Journal of Epidemiology. 2001, 154 (7): 666-674. 10.1093/aje/154.7.666.View ArticlePubMed
- Herikstad H, Yang S, van Gilder TJ, Vugia D, Hadler J, Blake P, Deneen V, Shiferaw B, Angulo F: A population-based estimate of the burden of diarrhoeal illness in the United States: FoodNet, 1996–7. Epidemiology and Infection. 2002, 129 (1): 9-17. 10.1017/S0950268801006628.PubMed CentralPubMed
- Kuusi M, Aavitsland P, Gondrosen B, Kapperud G: Incidence of gastroenteritis in Norway – a population-based survey. Epidemiology and Infection. 2003, 131 (1): 591-597. 10.1017/S0950268803008744.PubMed CentralView ArticlePubMed
- Scallan E, Fitzgerald M, Collins C, Crowley D, Daly L, Devine M, Igoe D, Quigley T, Robinson T, Smyth B: Acute gastroenteritis in northern Ireland and the Republic of Ireland: a telephone survey. Communicable Disease and Public Health. 2004, 7 (1): 61-67.PubMed
- Hall G, Kirk MD, Becker N, Gregory JE, Unicomb L, Millard G, Stafford R, Lalor K, the OzFoodNet Working Group: Estimating foodborne gastroenteritis, Australia. Emerging Infectious Diseases. 2005, 11 (8): 1257-1264.PubMed CentralView ArticlePubMed
- Majowicz SE, Doré K, Flint J, Edge VL, Read S, Buffet MC, McEwen SA, McNab WB, Stacey D, Sockett PN, Wilson JB: Magnitude and distribution of acute, self-reported gastrointestinal illness in a Canadian community. Epidemiology and Infection. 2004, 132 (4): 607-617. 10.1017/S0950268804002353.PubMed CentralView ArticlePubMed
- Thomas MK, Majowicz SE, MacDougall L, Sockett PN, Kovacs SJ, Fyfe M, Edge VL, Doré K, Flint JA, Henson S, Jones AQ: Population distribution and burden of acute gastrointestinal illness in British Columbia, Canada. BMC Public Health. 2006
- Tijssen JGP: Representativeness and response rates from the Domestic/International Gastroenterology Surveillance Study (DIGEST). Scandinavian Journal of Gastroenterology. 1999, 231: 15-19. 10.1080/003655299750025228.View Article
- de Wit MA, Koopmans MP, Kortbeek LM, van Leeuwen WJ, Bartelds AIM, van Duynhoven YT: Gastroenteritis in sentinel general practices, The Netherlands. Emerging Infectious Diseases. 2001, 7 (1): 82-91.PubMed CentralView ArticlePubMed
- Hosmer DW, Lemeshow S: Applied Logistic Regression. 2000, New York, John Wiley & Sons, Inc, 2View Article
- Monto AS, Koopman JS: The Tecumseh study. American Journal of Epidemiology. 1980, 112 (3): 323-333.PubMed
- Palmer S, Houston H, Lervy B, Ribeiro D, Thomas P: Problems in the diagnosis of foodborne infection in general practice. Epidemiology and Infection. 1996, 117 (3): 479-484.PubMed CentralView ArticlePubMed
- Stanghellini V: Three-month prevalence rates of gastrointestinal symptoms and the influence of demographic factors: results from the Domestic/International Gastroentorology Surveillance Study (DIGEST). Scandinavian Journal of Gastroenterology. 1999, 231 (Suppl): 20-28. 10.1080/003655299750025237.View Article
- Hall GV, Kirk MD, Ashbolt R, Stafford R, Lalor K: Frequency of infectious intestinal illness in Australia, 2002: regional, seasonal and demographic variation. Epidemiology and Infection. 2006, 134 (1): 111-118. 10.1017/S0950268805004656.PubMed CentralView ArticlePubMed
- Kagan LJ, Aiello AE, Larson E: The role of the home environment in the transmission of infectious diseases. Journal of Community Health. 2002, 27 (4): 247-267. 10.1023/A:1016378226861.View ArticlePubMed
- Gillespie IA, O'Brien SJ, Frost JA, Tam C, Tompkins D, Neal KR, Syed Q, Farthing MJG, The Campylobacter Sentinel Surveillance Scheme Collaborators: Investigating vomiting and/or bloody diarrhoea in Campylobacter jejuni infection. Journal of Medical Microbiology. 2006, 55: 741-746. 10.1099/jmm.0.46422-0.View ArticlePubMed
- Chin J: Control of communicable diseases manual. 2000, Washington, D.C, American Public Health Association, 17
- Mustard CA, Derksen S, Berthelot JM, Wolfson M, Roos LL: Age-specific education and income gradients in morbidity and mortality in a Canadian province. Social Science & Medicine. 1997, 45: 383-397. 10.1016/S0277-9536(96)00354-1.View Article
- Humphries KH, van Doorslaer E: Income-related health inequality in Canada. Social Science & Medicine. 2000, 50 (5): 663-671. 10.1016/S0277-9536(99)00319-6.View Article
- Berthelot J, Wilkins R, Ng E: Trends in mortality by neighbourhood income in urban Canada from 1971 to 1996. Statistics Canada Health Reports. 2002, 13: 45-
- Holtgrave DR, Crosby RA: Social capital, poverty, and income inequality as predictors of gonorrhoea, syphilis, chlamydia and AIDS case rates in the United States. Sexually Transmitted Infections. 2003, 79 (1): 62-64. 10.1136/sti.79.1.62.PubMed CentralView ArticlePubMed
- Jerrett M, Eyles J, Cole D: Socioeconomic and environmental covariates of premature mortality in Ontario. Social Science & Medicine. 1998, 47 (1): 33-49. 10.1016/S0277-9536(98)00008-2.View Article
- Mustard CA, Frohlich N: Socioeconomic status and the health of the population. Med Care. 1995, 33 (12 Suppl): DS43-54. 10.1097/00005650-199512001-00007.PubMed
- Czaja R, Blair J: Designing surveys: a guide to decisions and procedures. 2005, Thousand Oaks, California, Pine Forge Press, 2
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/7/162/prepub
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.