This article has Open Peer Review reports available.
Clinical, epidemiological, and spatial characteristics of Vibrio parahaemolyticus diarrhea and cholera in the urban slums of Kolkata, India
© Suman Kanungo et al.; licensee BioMed Central Ltd. 2012
Received: 17 June 2012
Accepted: 26 September 2012
Published: 28 September 2012
There is not much information on the differences in clinical, epidemiological and spatial characteristics of diarrhea due to V. cholerae and V. parahaemolyticus from non-coastal areas. We investigated the differences in clinical, epidemiological and spatial characteristics of the two Vibrio species in the urban slums of Kolkata, India.
The data of a cluster randomized cholera vaccine trial were used. We restricted the analysis to clusters assigned to placebo. Survival analysis of the time to the first episode was used to analyze risk factors for V. parahaemolyticus diarrhea or cholera. A spatial scan test was used to identify high risk areas for cholera and for V. parahaemolyticus diarrhea.
In total, 54,519 people from the placebo clusters were assembled. The incidence of cholera (1.30/1000/year) was significantly higher than that of V. parahaemolyticus diarrhea (0.63/1000/year). Cholera incidence was inversely related to age, whereas the risk of V. parahaemolyticus diarrhea was age-independent. The seasonality of diarrhea due to the two Vibrio species was similar. Cholera was distinguished by a higher frequency of severe dehydration, and V. parahaemolyticus diarrhea was by abdominal pain. Hindus and those who live in household not using boiled or treated water were more likely to have V. parahaemolyticus diarrhea. Young age, low socioeconomic status, and living closer to a project healthcare facility were associated with an increased risk for cholera. The high risk area for cholera differed from the high risk area for V. parahaemolyticus diarrhea.
We report coexistence of the two vibrios in the slums of Kolkata. The two etiologies of diarrhea had a similar seasonality but had distinguishing clinical features. The risk factors and the high risk areas for the two diseases differ from one another suggesting different modes of transmission of these two pathogens.
Vibrio parahaemolyticus diarrhea and cholera diarrhea (due to Vibrio cholerae O1 and less commonly V. cholerae O139) are both public health concerns. V. parahaemolyticus is a halophilic etiologic agent of diarrheal disease, having an ability to produce outbreaks of gastroenteritis . Recently, V. parahaemolyticus of specific serotypes were associated with the outbreaks in several parts of the world with the earliest cases being reported from Kolkata, India in 1996 [2, 3]. The bacterium was first identified as a cause of seafood-borne illness in Japan in 1950, when 272 individuals became ill and 20 died after the consumption of semidried juvenile sardines . It is associated with three major syndromes of clinical illness: gastroenteritis, wound infections, and septicemia. The most common syndrome is gastroenteritis; the symptoms include diarrhea with abdominal cramps, nausea, vomiting, headache, and low-grade fever . Between 1988 and 1997, a review of infections found that 88% of patients with V parahaemolyticus gastroenteritis and 91% of patients with V. parahaemolyticus primary septicemia had a known food history reported eating raw oysters . Thus, consumption of crustacean and molluscan shellfish has commonly been implicated in the transmission of V. parahaemolyticus.
Cholera infection is often mild or without symptoms, assuming 90% of the infections are asymptomatic , but sometimes it can be severe. Approximately one in 20 infected persons has severe disease characterized by profuse watery diarrhea, vomiting, and leg cramps. In these persons, rapid loss of body fluids leads to dehydration, acute renal failure and shock. Without treatment, death can occur within hours. The epidemiological patterns of cholera depend largely on environmental factors including sanitary conditions and social aspects, prior immune status of the population at risk, and the inherent properties of the vibrios themselves. Until the mid-1980s, humans were thought to be the only reservoir of V. cholerae. It is now believed that the organism has a free-living cycle and is a natural resident of aquatic environs . A significant marine reservoir of V. cholerae is plankton, and the bacterium attaches primarily to zooplankton, specifically copepods . It is now believed that the transmission of the vibrios might occur through water without fecal contamination, and the evidences suggest natural reservoir of two vibrios is the aquatic environment.
Although the reservoir of both the organisms (V. parahaemolyticus and V. cholerae) is believed to be marine environment, they have the ability to infect people in areas far from coastal areas. However, there is not much information on the differences in clinical, epidemiological and spatial characteristics of diarrhea due to V. cholerae and V. parahaemolyticus from non-coastal areas. This paper describes these characteristics in the urban slums of, a non-coastal area, Kolkata, India, where a cohort of population was under uniform surveillance for diarrhea.
The study area and data
The study was conducted in urban slum communities in Kolkata, the capital of the state of West Bengal. Kolkata, the third largest city in India, has 14 million inhabitants living within an area of 1,450 km2, making it one of the world’s most densely populated cities. The Kolkata Municipal Corporation consists of 141 civic administrative units called wards, with each ward having an office responsible for public health supervised by a medical officer. The study site comprises three contiguous wards (29, 30, and 33) with about 100,000 residents. These residents live in homes tightly-spaced together along winding sewage-littered pathways, and they rely on shared toilets and drinking-water . The area has a high population density. Through the years, extensive subletting has resulted in overcrowding as more and more people are squeezed into available housing . Sufficient water supply and sanitary facilities are unavailable in the area. Several households share one or two latrines and water taps. Most sewage is collected in open drainage gutters which tend to overflow during the rainy season, flooding adjacent homes. Kolkata has three seasons, the cool dry months from November to February, the hot dry period from March to May, and the monsoon season from June to October. Seventy percent of the people are Hindus, and the rest are predominately Muslims.
This study used the data of a cluster-randomized, double-blind, placebo-controlled trial of a killed oral cholera vaccine . The clusters were dwellings, which were randomly assigned to receive either vaccine or placebo, so that individuals living in the same dwelling (cluster) received the same agent (vaccine or placebo). A dwelling was defined as a hut, a group of huts, or a multistory building with several households using shared water pipes, bathrooms, and latrines as assigned by the Kolkata Municipal Corporation. There were 3933 clusters of which 1966 clusters were assigned to receive vaccine and 1967 clusters were assigned to receive placebo. Residents were eligible to receive a study intervention if they were aged one year or older and were not pregnant.
Nine clinics were established in the community to undertake diarrhea surveillance in addition to the two referral hospitals (Infectious Diseases Hospital and B.C. Roy Children’s Hospital). Private medical practitioners were encouraged to refer patients with diarrhea to the clinics. Patients from the study area were identified by use of household identification cards and a computerized database. Study physicians recorded pertinent clinical details on a structured clinical data form. Rectal swabs were obtained from all patients presenting with history of loose stools and transported in Cary-Blair media to a laboratory at the National Institute of Cholera and Enteric Diseases (NICED) within 8 h of specimen collection. At the laboratory, rectal swabs were examined for V. parahaemolyticus and V. cholerae by use of conventional methods .
The disease surveillance data of four years from January 1, 2007 to December 31, 2010 were used in this analysis. A diarrheal visit was defined as a visit by a patient who had, in the 24 h before presentation, three or more loose or liquid stools; or at least one loose or liquid stool with blood; or, if one to two or an indeterminate number of loose or liquid stools were reported, at least some evidence of dehydration, according to WHO criteria . The onset of a diarrheal visit was the day on which the patient first experienced loose or liquid stools. Diarrheal visits for which the date of onset was less than or equal to 7 days from the date of discharge for the previous visit were grouped into the same diarrheal episode. A cholera episode was one in which V. cholerae O1 or O139 was isolated from a fecal specimen during any component diarrheal visit; an episode of V. parahaemolyticus diarrhea was one in which V. parahaemolyticus was isolated.
To detect potential geographic areas of high risk for the diseases, the spatial scan test has been widely used in recent times [15–18]. The spatial scan test is also suitable for uneven geographic distribution of cases and population density . We used the spatial scan test implemented through SaTScan®  to identify unique non-random spatial clusters that are higher risk for cholera and for V. parahaemolyticus. We assumed that the incidences of two diseases followed a Poisson distribution. Under the null hypothesis, the incidence of disease in a particular location is proportional to the number of residents in that location . Using SaTScan®, we estimated the probability that the frequency of disease at each peak surpasses that expected by chance. We set the space limitations to 50% population at risk, which allowed us to scan for both large and small clusters of disease risk. We took into account the observed number of cases inside and outside the circle when calculating the highest likelihood for each circle. This circle was the most probable cluster and had a rate that was the least likely to happen by chance alone. The statistical significance of possible clusters was calculated using 999 Monte Carlo simulations . Purely spatial analysis was performed using circular windows. The output from SaTScan® was imported into ArcGIS (Version 9.2, California, USA) to map significant (p<0.01) clusters of higher risk.
This is an observational population based study in which the population was classified in clusters of residential dwelling units. The study participants were geographically identified based on their residence, and this spatial component was included in the study. Several demographic and socioeconomic variables that were thought to be independently associated with the risk for the diseases were evaluated in the study. However, this is an exploratory analysis of a secondary data source. To avoid distortions due to the impact of the cholera vaccination on the outcomes, we limited our analysis to residential dwellings assigned to placebo. Comparison of individual characteristics was performed using generalized estimating equations (GEE) with logit link function adjusting for the design effect of the clusters used to allocated subjects to the two agents in the randomized trial . To analyze the risk factors of V. parahaemolyticus diarrhea or cholera, we used survival analyses of the time to the first episode of the disease, censoring the follow-up of individuals who died or migrated out . We fitted unadjusted and covariate adjusted Cox proportional hazard regression models verifying first that the proportionality assumption was satisfied for all independent variables [23–25]. In both covariate-unadjusted and covariate-adjusted analyses of the risk for V. parahaemolyticus diarrhea or cholera, we accounted for the design effect induced by cluster randomization by use of robust sandwich variance estimates . These estimates enabled inferences about risk of the disease at the individual level, adjusting for the design effect. Final adjusted risk estimates were obtained from the model significant at p<0·05 in a forward selection algorithm.
The study received approval from the Health Ministry Screening Committee of the Government of India, Scientific Advisory Committee and Institutional Ethics Committee of NICED, Institutional Review board of International Vaccine Institute and also from the Secretariat Committee for Research Involving Human Subjects, World Health Organization Geneva, Switzerland.
Number of episodes (incidence rate/1000/year) of V. parahaemolyticus diarrhea and cholera in the study area during 2007–2010
Population at Jan 1, 2007
V. parahaemolyticusDiarrhea No. of episodes (incidence rate/1000/year)
Cholera No. of episodes (incidence rate/1000/year)
Clinical symptoms of V. parahaemolyticus diarrhea and cholera in the study area, 2007–2010. No (%)
V. parahaemolyticus diarrhea (n=137)
Blood in stool
Abnormal mental state
Socio demographic characteristics of the diarrhea cases infected by the two vibrio species in the study area, 2007–2010
Mean (SD) age (years)
No. (%) of children aged less than 5 years
No. (%) of males
No. (%) of Hindus
No. (%) of individuals living in a household with literate household head
No. (%) of individuals living in a household with household head having more than 5 years of schooling
No. (%) of individuals living in a household with important economic contributor having stable occupation1)
No. (%) of individuals living in a household using safe toilet2)
No. (%) of individuals living in a household using safe water source (tap at household level)
No. (%) of individuals living in a household using safe water source (tap or hand pump at household level)
No. (%) of individuals living in a household using boiled or filtered water
No. (%) of individuals living in a household always wash hands with soap and water after defecation
No. (%) of individuals living in a household having specific place for waste disposal
No. (%) of individuals living in their own house
No. (%) of individuals living in a household owning refrigerator
No. (%) of individuals living in a household owning motorbike
No. (%) of individuals living in a household owning telephone
No. (%) of individuals living in a household owning television
No. (%) of individuals living in a household owning all luxury items3)
No. (%) of individuals living in a household owning at least one luxury item3)
Mean (SD) monthly household expenditure (Indian Rupee)
No. (%) of individuals with high monthly household expenditure4)
Mean (SD) monthly per-capita expenditure5) of house hold (Indian Rupee)
No. (%) of individuals having high monthly per-capita expenditure4) of household
Mean (SD) Mean (SD) of the household size
Mean (SD) distance (m) from household to the nearest health clinic
No. (%) of individuals living in household longer distance to the nearest health clinic6)
Mean (SD) distance (m) from household to the nearest water body
No. (%) of individuals living in household longer distance to the nearest water body 7)
No. (%) of individuals living in Ward 33
Predictors of the risk of V. parahaemolyticus diarrhea in the study area, 2007–2010
Individuals living in a household using boiled or filtered water
Predictors of the risk of cholera in the study area, 2007–2010
Individuals 5 years and above
Individuals living in a household with important economic contributor having stable occupation1)
Individuals living in a household using safe toilet2)
Individuals living in a household owning at least one luxury item3)
Longer distance from the household to the nearest health clinic6)
This paper reports coexistence of V. parahaemolyticus and V. cholerae in the slums of Kolkata, India. Literature suggests foods frequently incriminated in V. parahaemolyticus infections are raw or inadequately cooked seafood and foods contaminated by seafood materials . However, the transmission and epidemiology of V. parahaemolyticus infections in the study area may be different because seafood is never eaten raw and freshwater fish is preferred over seawater fish by the local people.3 Contamination by seawater fish at the fish market, and secondary contamination of other foods in the kitchen by V. parahaemolyticus-contaminated fish brought from markets as well as cutting and dressing of sea food especially prawn (affecting women more) may be the routes of transmission of the V. parahaemolyticus in the study area [27–29].
The clinical symptoms and seasonality for both the diseases are almost identical, thus it is difficult to diagnose the patients without microbiological test of the patients’ fecal specimens. However, some degree of differentiation of these patients may be made by looking for severity of diarrhoea and stomach pain. The incidence of cholera was higher than that of V. parahaemolyticus diarrhoea in the study area. The incidence rates reported here are certainly underestimates of the true rates, because the cases were detected only through augmented passive surveillance in project health clinics. The higher risk for cholera for people living closer to the project health facilities also suggests that some patients living far from the project health facilities might not have sought their care for diarrhea from the project health clinics. The proportion of patients with diarrhea in the study area seeking treatment in those clinics was not assessed.
It was interesting to note that the Hindus, the dominant religious group, were at greater risk for V. parahaemolyticus diarrhea compared to Muslims. The Hindus are relatively affluent in that society in comparison to Muslims. A study in coastal Vietnam  showed that more affluent members of the community, assessed by their higher professional status, better living conditions, and possession of luxury objects, were more frequently infected by V. parahaemolyticus compared to less affluent individuals. A possible explanation for that finding could be that only the more affluent members of the community can afford to include fish in their diet, which is thought to be the source of infection.
The cluster of significantly higher risk for cholera was observed in a part of Ward 29 where the density of population is very high. A study observed that the density of refuse dumps is an important environmental predictor of cholera . The refuse dump was considered as an index of basic sanitation in that study. It is reasonable to believe that an area with higher density of population may create increased density of refuse dumps, and the resulting breakdown of sanitation can lead to a higher risk for the cholera in the area. Likewise, the cause of clustering of the V. parahaemolyticus diarrhea in ward 30 may be explained by the fact that it has predominantly Hindu population and was thereby more prone to V. parahaemolyticus diarrhea (as mentioned above). However, further in-depth studies are needed to find the cause of spatial risk for the two organisms in the study area.
There is not much information on the differences in clinical, epidemiological and spatial characteristics of diarrhea due to Vibrio Cholerae and V parhaemolyticus from non-coastal areas, and we observed some distinctive risk factors and spatial patterns of risk for diarrhea due to cholera and V. parahemolyticus suggesting different modes of transmission of these two pathogens. This information may be helpful for the health policy makers. However, further research is needed to delineate the modes of transmission of the two vibrios.
This work was supported by the Bill & Melinda Gates Foundation through the Diseases of the Most Impoverished Program and the Cholera Vaccine Initiative. Additional funding was provided by the Swedish International Development Cooperation Agency and the Governments of South Korea and, Sweden. Shantha Biotechnics donated vaccine and placebo for the study. We are grateful to the people of Kolkata where our study is being undertaken and to our field staff who provided valuable support in our study.
- Tantillo GM, Fontanarosa M, Di Pinto A, Musti M: Updated perspectives on emerging vibrios associated with human infections. Lett Appl Microbiol. 2004, 39: 117-126. 10.1111/j.1472-765X.2004.01568.x.View ArticlePubMedGoogle Scholar
- Okuda J, Ishibashi M, Hayakawa E, Nishino T, Takeda Y, Mukhopadhyay AK, Garg S, Bhattacharya SK, Nair GB, Nishibuchi M: Emergence of a unique O3:K6 clone of Vibrio parahaemolyticus in Calcutta, India, and isolation of strains from the same clonal group from Southeast Asian travelers arriving in Japan. J Clin Microbiol. 1997, 35 (12): 3150-3155.PubMedPubMed CentralGoogle Scholar
- Nair GB, Ramamurthy T, Bhattacharya SK, Dutta B, Takeda Y, Sack DA: Global dissemination of Vibrio parahaemolyticus serotype O3:K6 and its serovariants. Clin Microbiol Rev. 2007, 20 (1): 39-48. 10.1128/CMR.00025-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Fujino T, Okuno Y, Nakada D, Aoyama A, Fukai K, Mukai T, Uebo T: On the bacteriological examination of Shirasu food poisoning. Med J Osaka Univ. 1953, 4: 299-304.Google Scholar
- Honda T, Iida T: The pathogenicity of Vibrio parahaemolyticus and the role of the thermostable direct haemolysin and related haemolysins. Rev. Med. Microbiol. 1993, 4: 106-113. 10.1097/00013542-199304000-00006.View ArticleGoogle Scholar
- Daniels NA, MacKinnon L, Bishop R, Altekruse S, Ray B, Hammond RM, Thompson S, Wilson S, Bean NH, Griffin PM, Slutsker L: Vibrio parahaemolyticus infections in the United States, 1973–1998. J Infect Dis. 2000, 181: 1661-1666. 10.1086/315459.View ArticlePubMedGoogle Scholar
- Longini IM, Azhar N, Ali M, Yunus M, Shenvi N, Clemens J: Controlling endemic cholera with oral vaccine. PLoS Med. 2007, 4 (11): e336-10.1371/journal.pmed.0040336.View ArticlePubMedPubMed CentralGoogle Scholar
- Feachem RG: Environmental aspects of cholera epidemiology. II. Occurrences and survival of Vibrio cholerae in the environment. Trop Dis Bull. 1981, 78: 865-880.PubMedGoogle Scholar
- Huq A, Colwell RR: Vibrios in the marine and estuarine environments. J Mar Biotechnol. 1995, 3: 60-63.Google Scholar
- Sur D, Deen JL, Manna B, Niyogi SK, Deb AK, Kanungo S, Sarkar BL, Kim DR, Danovaro-Holliday MC, Holliday K, Gupta VK, Ali M, von Seidlein L, Clemens JD, Bhattacharya SK: The burden of cholera in the slums of Kolkata, India: data from a prospective, community based study. Arch Dis Child. 2005, 90: 1175-1181. 10.1136/adc.2004.071316.View ArticlePubMedPubMed CentralGoogle Scholar
- Thomas FC: Calcutta: The human face of poverty. 1997, London: Penguin Books Ltd, 3-166.Google Scholar
- Sur D, Lopez AL, Kanungo S, Paisley A, Manna B, Ali M, Niyogi SK, Park JK, Sarkar BL, Puri MK, Kim DR, Deen J, Holmgren J, Carbis R, Rao R, Van NT, Donner A, Ganguli NK, Nair GB, Bhattacharya SK, Clemens JD: Efficacy and safety of a modified killed whole-cell oral cholera vaccine in India: an interim analysis of a cluster-randomized, double-blind, placebo-controlled trial. Lancet. 2009, 374: 1694-1702. 10.1016/S0140-6736(09)61297-6.View ArticlePubMedGoogle Scholar
- Bopp C, Ries A, Wells J: Laboratory methods for the diagnosis of epidemic dysentery and cholera. 1999, Atlanta, GA: Centers for Disease Control and PreventionGoogle Scholar
- World Health Organization: The treatment of diarrhea: a manual for physicians and other senior health workers, 4th revision. 2005, Geneva: WHOGoogle Scholar
- Kulldorff M, Feuer EJ, Miller BA, Freedman LS: Breast cancer clusters in the northeast United States: a geographic analysis. Am J Epidemiol. 1997, 146: 161-170. 10.1093/oxfordjournals.aje.a009247.View ArticlePubMedGoogle Scholar
- Kulldorff M, Athas WF, Feurer EJ, Miller BA, Key CR: Evaluating cluster alarms: a space-time scan statistic and brain cancer in Los Alamos, New Mexico. Am J Public Health. 1998, 88: 1377-1380. 10.2105/AJPH.88.9.1377.View ArticlePubMedPubMed CentralGoogle Scholar
- Hjalmars U, Kulldorff M, Gustafsson G, Nagarwalla N: Childhood leukaemia in Sweden: using GIS and a spatial scan statistic for cluster detection. Stat Med. 1996, 15: 707-715. 10.1002/(SICI)1097-0258(19960415)15:7/9<707::AID-SIM242>3.0.CO;2-4.View ArticlePubMedGoogle Scholar
- Gosselin P, Lebel G, Rivest S, Douville-Fradet M: The integrated system for public health monitoring of West Nile virus (ISPHM-WNV): a real-time GIS for surveillance and decision-making. Int J Health Geogr. 2005, 4: 21-10.1186/1476-072X-4-21.View ArticlePubMedPubMed CentralGoogle Scholar
- Sheehan TJ, DeChello LM, Kulldorff M, Gregorio DI, Gershman S, Mroszczyk M: The geographic distribution of breast cancer incidence in Massachusetts 1988 to 1997, adjusted for covariates. Int J Health Geogr. 2004, 3: 17-10.1186/1476-072X-3-17.View ArticleGoogle Scholar
- Mostashari F, Kulldorff M, Hartman JJ, Miller JR, Kulasekera V: Dead bird clusters as an early warning system for West Nile virus activity. Emerg Infect Dis. 2003, 9: 641-646. 10.3201/eid0906.020794.View ArticlePubMedPubMed CentralGoogle Scholar
- Liang K, Zeger S: Longitudinal data analysis using generalized linear models. Biometrika. 1986, 73: 13-22. 10.1093/biomet/73.1.13.View ArticleGoogle Scholar
- Miller R: Survival analysis. 1981, New York: WileyGoogle Scholar
- Reid N, Crepeau H: Influence functions for proportional hazards regression. Biometrika. 1985, 72: 1-9. 10.1093/biomet/72.1.1.View ArticleGoogle Scholar
- Lin D, Wei L, Ying Z: Checking the Cox model with cumulative sums of Martingale-based residuals. Biometrika. 1993, 80: 557-572. 10.1093/biomet/80.3.557.View ArticleGoogle Scholar
- Lin DY, Wei LJ: The robust inference for the proportional hazards model. J Am Stat Assoc. 1989, 84: 1074-1078. 10.1080/01621459.1989.10478874.View ArticleGoogle Scholar
- Tuyet DT, Thiem VD, von Seidlein L, Chowdhury A, Park E, Canh DG, Chien BT, Tung TV, Naficy A, Rao MR, Ali M, Lee H, Sy TH, Nichibuchi M, Clemens J, Trach DD: Clinical, epidemiological, and socioeconomic analysis of an outbreak of Vibrio parahaemolyticus in Khanh Hoa Province, Vietnam. J Infect Dis. 2002, 186: 1615-1620. 10.1086/345731.View ArticlePubMedGoogle Scholar
- Pal SC, Sircar BK, Nair GB, Deb BC: Epidemiology of bacterial diarrhoeal diseases in India with special reference to Vibrio parahaemolyticus infections. Bacterial diarrhoeal disease. Edited by: Takeda Y, Miwatani T. 1984, Tokyo: KTK Scientific Publishers, 65-73.Google Scholar
- Sarkar BL, Nair GB, Sircar BK, Pal SC: Incidence and level of Vibrio parahaemolyticus associated with freshwater plankton. Appl Environ Microbiol. 1983, 46 (1): 288-290.PubMedPubMed CentralGoogle Scholar
- Sarkar BL, Nair GB, Banerjee AK, Pal SC: Seasonal distribution of Vibrio parahaemolyticus in freshwater environs and in association with freshwater fishes in Calcutta. Appl Environ Microbiol. 1985, 49 (1): 132-136.PubMedPubMed CentralGoogle Scholar
- Osei FB, Duker AA: Spatial dependency of V. cholera prevalence on open space refuse dumps in Kumasi, Ghana: a spatial statistical modelling. Int J Health Geogr. 2008, 7: 62-10.1186/1476-072X-7-62.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/12/830/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.