- Research article
- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
Mortality after infection with methicillin-resistant Staphylococcus aureus(MRSA) diagnosed in the community
© Delaney et al; licensee BioMed Central Ltd. 2008
Received: 15 June 2007
Accepted: 31 January 2008
Published: 31 January 2008
Outbreak reports suggest that community-acquired methicillin-resistant Staphylococcus aureus (MRSA) infections can be life-threatening. We conducted a population based cohort study to assess the magnitude of mortality associated with MRSA infections diagnosed in the community.
We used the United Kingdom's General Practice Research Database (GPRD) to form a cohort of all patients with MRSA diagnosed in the community from 2001 through 2004 and up to ten patients without an MRSA diagnosis. The latter were frequency-matched with the MRSA patients on age, GPRD practice and diagnosis date. All patients were older than 18 years, had no hospitalization in the 2 years prior to cohort entry and medical history information of at least 2 years prior to cohort entry. The cohort was followed up for 1 year and all deaths and hospitalizations were identified. Hazard ratios of all-cause mortality were estimated using the Cox proportional hazards model adjusted for patient characteristics.
The cohort included 1439 patients diagnosed with MRSA and 14,090 patients with no MRSA diagnosis. Mean age at cohort entry was 70 years in both groups, while co-morbid conditions were more prevalent in the patients with MRSA. Within 1 year, 21.8% of MRSA patients died as compared with 5.0% of non-MRSA patients. The risk of death was increased in patients diagnosed with MRSA in the community (adjusted hazard ratio 4.1; 95% confidence interval: 3.5–4.7).
MRSA infections diagnosed in the community are associated with significant mortality in the year after diagnosis.
Methicillin-resistant Staphylococcus aureus (MRSA) is increasingly implicated in potentially lethal hospital-acquired infections [1, 2]. Some of the MRSA strains that are found in patients in the community have been acquired during hospitalization or a visit to a hospital. However, MRSA strains likely originating outside of hospitals and different from those found in hospitals have been identified [3–5]. People who are infected by MRSA in the community have different risk factors than those infected in hospitals [1, 2], particularly as common risk factors in hospitals may be rare or absent in the community.
While hospital-acquired MRSA infections can be fatal [1, 2], cases of severe and life-threatening MRSA infections from the community have also been reported [6–8] and these case reports suggest that the prognosis of community-acquired MRSA infections may be poor . However, hospital-based studies do not capture patients from the community whose infection does not require hospital care, and thus such studies or case reports cannot evaluate the impact of this infection on mortality. To date, no population-based study has evaluated the prognosis and mortality of MRSA patients diagnosed in the community.
We thus conducted a large cohort study to assess mortality associated with MRSA infections diagnosed in the community in comparison with a population-based sample of patients free of this infection.
We used the General Practice Research Database (GPRD) between 1 January 2001 and 31 December 2005 for this study. The GPRD contains the diagnostic, testing and prescribing records of approximately 3.2 million patients from over 400 general practices in the United Kingdom (UK) as recorded by their GPs. The GPRD is extensively used for medical research [10, 11] and has been used previously for prognostic work in cardiovascular disease , as well as infectious disease research [13, 14]. GPRD data is routinely audited for quality  and is considered to be of high validity for medical research.
Cohort definition and follow-up
We conducted a frequency-matched cohort study. The cohort consisted of all patients with MRSA diagnosed between 1 January 2001 and 31 December 2004, and a population-based sample of patients without MRSA. We defined MRSA to be the diagnosis 'Methicillin-resistant Staphylococcus aureus positive' (READ code 4JP..00) entered into the patient's medical record. The date of the MRSA diagnosis was the MRSA patient's cohort entry date. For each patient with MRSA, we randomly selected up to 10 patients free of MRSA that were registered in the GPRD at their corresponding MRSA patient's cohort entry date and that were matching on practice, cohort entry date and age (± 2 years). We assigned to these MRSA-free patients their corresponding MRSA patient's cohort entry date. Patients of the cohort, both with and without MRSA, had to be at least 18 years of age at their cohort entry date, had been registered in the GPRD for at least 2 years prior to their cohort entry date and had no hospitalizations recorded in the GPRD in the 2 years before cohort entry date. Frequency matching on cohort entry date is particularly important to control for calendar time effects, because the rate of MRSA diagnosis in the GPRD was not constant over the time period of our study . This time variation in risk could have been a potential source of bias in our estimates of effect . By frequency matching on practice we also ensured that patients with and without an MRSA diagnosis came from the same geographic region of the UK, thus indirectly controlling for factors such as socio-economic status and geographical variations in disease frequency.
Follow-up was from day 1 to day 365 after cohort entry date. Patients were censored if they transferred out of practice for a reason other than death or were free of the outcome 365 days after cohort entry. Patients who died on day 0 were excluded. It was possible for a patient without MRSA to be diagnosed with MRSA during follow-up.
The outcome, death, was defined as a transfer out of the GPRD with the reason being 'Death' recorded in the patient database of the GPRD. We used the date of the transfer out as the date of death. As the GPRD is based on GP records and deaths are a variable that practices in the GPRD are audited on , death should be a well-recorded variable in this database.
For a secondary outcome, to assess the level of morbidity caused by MRSA infections diagnosed in the community, we considered a composite outcome of either death or hospitalization. For this outcome, the date of outcome was either the date of hospitalization or the date of death; whichever came first. The recording of reason for hospitalization is poor in the GPRD, but the hospital visit itself is well recorded [10, 11]. Therefore, we were not able to separate hospitalization by reason and used 'any cause' hospitalization. The recording of hospitalization information also made it impossible to assess whether or not patients diagnosed with MRSA had serious infectious complications such as bacteraemia.
In addition to age and sex, we defined a series of baseline covariates that could have an impact on the prognosis of individuals diagnosed with MRSA in the community. We included the presence of cancer, heart disease, renal failure, autoimmune diseases and diabetes, all identified using medical records in the GPRD in the 2 years prior to cohort entry date. We also considered use of antibiotics and oral prednisone in the year before cohort entry date.
We used the Kaplan-Meier method to estimate the cumulative mortality in patients with and without MRSA over the one-year follow-up. The Cox proportional hazards model was used to estimate the hazard ratio of death associated with MRSA infection, adjusting for covariates . Proportionality of hazards was assessed using graphical methods.
We conducted several sensitivity analyses. First, to determine the maximum possible effect of censoring on the results, we used both worst-case and best-case imputation, appropriate techniques because of the small number of censored observations . For worst-case imputation, we assumed that all MRSA patients who were censored die immediately and that all patients without MRSA lived until the end of follow-up. For best-case imputation, we assumed that all censored MRSA patients lived until the end of follow-up and all censored MRSA-free patients died immediately. Second, as an additional control for overall health status, we considered past antibiotic use as a marker of susceptibility to infection. Third, we determined the effect of MRSA diagnosis among patients with none of the co-morbid conditions in our study to estimate the effect of MRSA among healthy patients. We used SAS version 9.1.3 in all analyses.
We obtained ethics approval for this study from the Scientific and Ethical Advisory Group of the GPRD and the McGill University Health Center Research Ethics Board.
Baseline characteristics of patients diagnosed with MRSA (1439 subjects; mean ± SD age 70.3 ± 18.3 years) and a matched population sample free of MRSA (14,090 subjects; mean ± SD age 69.9 ± 18.2 years) from GPRD (2001–2004)
Patients with MRSA
Matched MRSA-free patients
Co-morbidity in the 2 years prior to cohort entry
Peripheral vascular disease
Drug exposure in the 1 year prior to cohort entry
Oral prednisone use
All-cause mortality and hospitalization or all-cause mortality within a year in patients diagnosed with MRSA compared with matched patients free of MRSA from the GPRD (2001–2004) The hazard ratio is adjusted for all variables in Table 1. Death or hospitalization patients were followed until first event of either kind and only counted once.
Person-years of follow-up
Rate of event per 100 years
Adjusted HR (95% CI)
Patients with MRSA
Patients without MRSA
Death or hospitalization
Patients with MRSA
Patients without MRSA
The sensitivity analysis for informative censoring produced an adjusted HR of 4.91 (95% CI: 4.30–5.61) for worst-case imputation, and an adjusted HR of 2.69 (95% CI: 2.36–3.07) for best-case imputation. Prior antibiotic use had a small effect on mortality once we adjusted for other factors with an adjusted HR of 1.30 (95% CI: 1.14–1.49). Finally, among the subset of patients with no co-morbidities, the effect of MRSA infection on mortality was comparable to that obtained for all patients with an adjusted HR of 4.86 (95% CI: 4.10–5.77).
In this study, we demonstrated an increase in mortality within a year of a diagnosis with an MRSA infection in the community. This increased mortality cannot be completely explained by underlying co-morbid conditions as we adjusted for the major risk factors of death. The four-fold increase in mortality persisted across several sensitivity analyses. In addition, we found a risk increase of similar magnitude with the combined outcome of death and hospitalization.
Patients with nosocomial MRSA bacteraemia [18, 19] have been reported as having as great as a ten-fold increased risk of death . However, these extremely severe infections are not representative of all nosocomial MRSA infections, and likely not representative of community-acquired MRSA infections, many of which do not require hospitalization . As expected, our observed death rate is much lower than in these studies.
In this study we observe an average of 360 MRSA infections diagnosed in the community per year among the adult GPRD population during the study period. This is about 13 cases per 100,000 adults in the GPRD population, which is similar to the incidence found in other studies such as that by Fridkin et al . This is especially true if we assume that some MRSA infections are diagnosed in hospital emergency rooms and thus not counted in our study. Fridkin et al also found a significant proportion of community-associated MRSA infections occurring in patients above 65 years of age .
In this study we considered a large cohort of patients diagnosed with MRSA infections in the community and a large group of matching disease-free patients all selected from a database that is representative of the UK population . The death rate among the disease-free patients is broadly consistent for this age and sex distribution with that seen in UK vital statistics . Using the GPRD allows us to study patients with all levels of disease severity encountered by GPs, including those not treated in hospitals and thus not part of any hospital-based investigation. The study design (frequency matched cohort study) ensured that the distribution of key covariates (general practice, age and calendar time) was balanced between patients with and without MRSA at baseline. By matching on general practice we indirectly matched on unmeasured factors such as socio-economic status, area of residence and nursing home care, which are usually common to the patient population of a given general practice. These demographic factors have been shown to be important in the epidemiology of MRSA . By matching for practice we control for the potential confounding nature of these variables.
There are nevertheless some limitations to our study. In particular, we could not distinguish site or severity of individual MRSA infections. As we could not examine antibiotics given in hospital, it was not feasible to use post-diagnosis treatment either to separate infection site or infection source as different therapy is recommended for community-acquired MRSA . However, our observed mortality of MRSA infections diagnosed by general practitioners in the community rules out that these are the extremely severe infections seen in hospitals . On the other hand, these results imply that individuals with MRSA diagnoses in the GPRD are not mostly asymptomatic carriers but include clinically significant infections: a conclusion supported by other work in the GPRD on MRSA infections . Moreover, we strongly believe that currently, general practitioners in the UK do not routinely screen and thus diagnose MRSA in asymptomatic individuals in the community.
We could not exclude the possibility that some MRSA infections were acquired in hospitals and not in the community, because we had no information on potential transmission from visits to or work in hospitals. We could not look at the spread of MRSA from close contacts  or clusters in workplaces or schools . Therefore, we describe here the prognosis of MRSA infections detected by general practitioners, but we cannot demonstrate that these infections arose independently in the community. However, the high mortality rate underlines the importance of their control. To contribute to achieving this, we recently proposed the appropriate use of antibiotics .
The increased mortality observed for MRSA infections may be a result of a health state that predisposes to MRSA infections or residual confounding by this health state. Frailty and residing in a nursing home are examples of potential confounders that we were unable to adjust for in this study. If the mortality associated with MRSA is a manifestation of underlying susceptibility to infections, then past antibiotic use could serve as an estimate of the degree of bias in our study. Past antibiotic use should not be associated with distant outcomes but will show patients who have a history of infections. While we found a small increase in mortality predicted by antibiotic prescriptions given in the year prior to cohort entry, it was much too small to be a major confounder of the observed increase in mortality. Therefore, some degree of susceptibility to infection may explain a small part of the increased death rate. However, because of the small size of this effect we can infer that the risk of death in this study is not principally a marker of an underlying susceptibility to infections.
We did not report on the effect of antibiotic therapy on mortality after MRSA diagnosis for three reasons. One, the GPRD does not report drug exposures during hospitalization. This would result in the most severe cases of MRSA infection being misclassified as unexposed. This could lead to a falsely protective effect of antibiotic use if it was only recorded among the healthiest patients. Two, the use of antibiotics would be a dynamic treatment regime and more advanced statistical methods would be required to obtain a valid estimate . Three, the GPRD recording system does not allow us to determine whether an antibiotic prescription is being given for the original infection or for a second infection that developed independently.
While previous work in this database has shown differences between classes of antibiotics and the development of MRSA diagnosed in the community , we did not see any strong link between antibiotic class and prognosis. Therefore, we considered antibiotics only as a single drug class for this analysis as the principal reason for including antibiotics as a covariate was to control for possible differences in susceptibility to infections between patients with and without MRSA.
The broad use of sensitivity analysis [26, 27] is a key strength of our study. The results of these analyses indicate that the increase in mortality after a diagnosis of MRSA is robust to changes in assumptions. The use of the READ code 4JP..00 for MRSA infection has been discussed previously  and patients with this code were more similar in characteristics to post-operative MRSA wound infections than to carriers. The very high mortality rate seen in patients with MRSA in the present study, independent of baseline health status, is also highly suggestive of infections rather than carrier status. However, carrier status may be associated with lower mortality and hospitalization rate in comparison to clinically relevant infections. Thus, including this code may have led to an underestimation of the impact of a diagnosis of MRSA in the community on hospitalization and mortality.
The largest remaining source of possible bias in this study is the possibility of unmeasured confounding or, more likely, residual confounding. While non-parametric approaches could be used to reduce residual confounding in this study, they are unlikely to explain effects of the size that we observe . It is not likely that such large effects could be explained by anything other than many extremely strong confounders  based on estimates of the effect of confounding. We also avoided introducing bias from either the misallocation of person time  or from adjusting for variables potentially in the causal pathway between MRSA infection and mortality  by defining confounders using only the baseline characteristics of the subjects. Adjusting for factors (such as secondary infections) that arise after the diagnosis of MRSA infection could lead to overadjustment  and could bias the estimates towards the null.
In conclusion, our study suggests that MRSA can be a potentially serious infection in the community leading to increased mortality. Further research will allow a better understanding of MRSA diagnosed in the community and the prognosis of patients with this diagnosis. This is essential for effective prevention of MRSA arisen from the community, which will not only reduce hospital admissions with infections, spread into and within hospitals, and transmission in the community, but also likely will prevent deaths in the community.
The authors want to gratefully acknowledge Dr Tim Williams and the GPRD research team who provided feedback on MRSA endpoints in the GPRD. This study was funded by the Canadian Institutes of Health Research (CIHR) and the Canadian Foundation for Innovation. VSL is the recipient of a fellowship from CIHR, PB is supported by a new investigator career award from CIHR and SS is the recipient of a Distinguished Investigator Award from CIHR.
- Grundmann H, Aires-de-Sousa M, Boyce J, Tiemersma E: Emergence and resurgence of methicillin-resistant Staphylococcus aureus as a public-health threat. Lancet. 2006, 368: 874-885. 10.1016/S0140-6736(06)68853-3.View ArticlePubMedGoogle Scholar
- Zetola N, Francis JS, Nuermberger EL, Bishai WR: Community-acquired meticillin-resistant Staphylococcus aureus: an emerging threat. Lancet Infect Dis. 2005, 5: 275-286. 10.1016/S1473-3099(05)70112-2.View ArticlePubMedGoogle Scholar
- Robinson DA, Kearns AM, Holmes A, Morrison D, Grundmann H, Edwards G, O'Brien FG, Tenover FC, McDougal LK, Monk AB, Enright MC: Re-emergence of early pandemic Staphylococcus aureus as a community-acquired meticillin-resistant clone. Lancet. 2005, 365: 1256-1258. 10.1016/S0140-6736(05)74814-5.View ArticlePubMedGoogle Scholar
- Fridkin SK, Hageman JC, Morrison M, Sanza LT, Como-Sabetti K, Jernigan JA, Harriman K, Harrison LH, Lynfield R, Farley MM, Active Bacterial Core Surveillance Program of the Emerging Infections Program Network: Methicillin-resistant Staphylococcus aureus disease in three communities. N Engl J Med. 2005, 352: 1436-1444. 10.1056/NEJMoa043252.View ArticlePubMedGoogle Scholar
- Moran GJ, Amii RN, Abrahamian FM, Talan DA: Methicillin-resistant Staphylococcus aureus in community-acquired skin infections. Emerg Infect Dis. 2005, 11: 928-930.View ArticlePubMedPubMed CentralGoogle Scholar
- Maltezou HC, Giamarellou H: Community-acquired methicillin-resistant Staphylococcus aureus infections. Int J Antimicrob Agents. 2006, 27: 87-96. 10.1016/j.ijantimicag.2005.11.004.View ArticlePubMedGoogle Scholar
- Francis JS, Doherty MC, Lopatin U, Johnston CP, Sinha G, Ross T, Cai M, Hansel NN, Perl T, Ticehurst JR, Carroll K, Thomas DL, Nuermberger E, Bartlett JG: Severe community-onset pneumonia in healthy adults caused by methicillin-resistant Staphylococcus aureus carrying the Panton-Valentine leukocidin genes. Clin Infect Dis. 2005, 40: 100-107. 10.1086/427148.View ArticlePubMedGoogle Scholar
- Gilbert M, MacDonald J, Gregson D, Siushansian J, Zhang K, Elsayed S, Laupland K, Louie T, Hope K, Mulvey M, Gillespie J, Nielsen D, Wheeler V, Louie M, Honish A, Keays G, Conly J: Outbreak in Alberta of community-acquired (USA300) methicillin-resistant Staphylococcus aureus in people with a history of drug use, homelessness or incarceration. CMAJ. 2006, 175: 149-154.View ArticlePubMedPubMed CentralGoogle Scholar
- Miller LG, Quan C, Shay A, Mostafaie K, Bharadwa K, Tan N, Matayoshi K, Cronin J, Tan J, Tagudar G, Bayer AS: A prospective investigation of outcomes after hospital discharge for endemic, community-acquired methicillin-resistant and -susceptible Staphylococcus aureus skin infection. Clin Infect Dis. 2007, 44: 483-492. 10.1086/511041.View ArticlePubMedGoogle Scholar
- Walley T, Mantgani A: The UK General Practice Research Database. Lancet. 1997, 350: 1097-1099. 10.1016/S0140-6736(97)04248-7.View ArticlePubMedGoogle Scholar
- Van Staa TP, Abenhaim L: The quality of information recorded on a UK database of primary care records: a study of hospitalization due to hypoglycemia and other conditions. Pharmacoepidemiol Drug Saf. 1994, 3: 15-21. 10.1002/pds.2630030106.View ArticleGoogle Scholar
- Ruigomez A, Johansson S, Wallander MA, Garcia Rodriguez LA: Predictors and prognosis of paroxysmal atrial fibrillation in general practice in the UK. BMC Cardiovasc Disord. 2005, 5: 20-10.1186/1471-2261-5-20.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneider-Lindner V, Delaney JA, Dial S, Dascal A, Suissa S: Antibacterial drugs and the risk of infection with methicillin resistant Staphylococcus aureus (MRSA) in the community, United Kingdom. Emerg Infect Dis. 2007, 13: 994-1000.View ArticlePubMedPubMed CentralGoogle Scholar
- Dial S, Delaney JA, Schneider V, Suissa S: Proton pump inhibitor use and the risk of treated Clostridium difficile associated disease in the community. CMAJ. 2006, 175: 745-748.View ArticlePubMedPubMed CentralGoogle Scholar
- Card TR, Solaymani-Dodaran M, Hubbard R, Logan RF, West J: Is an internal comparison better than using national data when estimating mortality in longitudinal studies?. J Epidemiol Community Health. 2006, 60: 819-821. 10.1136/jech.2005.041202.View ArticlePubMedPubMed CentralGoogle Scholar
- Cox DR: Regression models and life tables (with discussion). J R Stat Soc B. 1972, 34: 187-220.Google Scholar
- Touloumi G, Pocock SJ, Babiker AG, Darbyshire JH: Impact of missing data due to selective dropouts in cohort studies and clinical trials. Epidemiology. 2002, 13: 347-355. 10.1097/00001648-200205000-00017.View ArticlePubMedGoogle Scholar
- Wyllie DH, Crook DW, Peto TE: Mortality after Staphylococcus aureus bacteraemia in two hospitals in Oxfordshire, 1997–2003: cohort study. BMJ. 2006, 333: 281-10.1136/bmj.38834.421713.2F.View ArticlePubMedPubMed CentralGoogle Scholar
- Guilarde AO, Turchi MD, Martelli CM, Primo MG: Staphylococcus aureus bacteraemia: incidence, risk factors and predictors for death in a Brazilian teaching hospital. J Hosp Infect. 2006, 63: 330-336. 10.1016/j.jhin.2006.02.011.View ArticlePubMedGoogle Scholar
- National Statistics: Deaths: age and sex, numbers and rates, 1976 onwards (England and Wales): Population Trends 130. [http://www.statistics.gov.uk/STATBASE/ssdataset.asp?vlnk=9552]
- Kerttula AM, Lyytikäinen O, Salmenlinna S, Vuopio-Varkila J: Changing epidemiology of methicillin-resistant Staphylococcus aureus in Finland. J Hosp Infect. 2004, 58: 109-114. 10.1016/j.jhin.2004.05.019.View ArticlePubMedGoogle Scholar
- Micek ST, Dunne M, Kollef MH: Pleuropulmonary complications of Panton-Valentine leukocidin-positive community-acquired methicillin-resistant Staphylococcus aureus: importance of treatment with antimicrobials inhibiting exotoxin production. Chest. 2005, 128: 2732-2738. 10.1378/chest.128.4.2732.View ArticlePubMedGoogle Scholar
- Moran GJ, Krishnadasan A, Gorwitz RJ, Fosheim GE, McDougal LK, Carey RB, Talan DA, EMERGEncy ID Net Study Group: Methicillin-resistant S. aureus infections among patients in the emergency department. N Engl J Med. 2006, 355: 666-674. 10.1056/NEJMoa055356.View ArticlePubMedGoogle Scholar
- Lo WT, Lin WJ, Tseng MH, Lu JJ, Chu ML, Wang CC: Nasal carriage of a single clone of community-acquired methicillin-resistant Staphylococcus aureus among kindergarten attendees in northern Taiwan. BMC Infect Dis. 2007, 7: 51-10.1186/1471-2334-7-51.View ArticlePubMedPubMed CentralGoogle Scholar
- Moodie EM, Richardson TS, Stephens DA: Demystifying optimal dynamic treatment regimes. Biometrics. 2007, 63: 447-455. 10.1111/j.1541-0420.2006.00686.x.View ArticlePubMedGoogle Scholar
- Greenland S: Basic methods for sensitivity analysis of biases. Int J Epidemiol. 1996, 25: 1107-1116. 10.1093/ije/25.6.1107-a.View ArticlePubMedGoogle Scholar
- Schneeweiss S: Sensitivity analysis and external adjustment for unmeasured confounders in epidemiologic database studies of therapeutics. Pharmacoepidemiol Drug Saf. 2006, 15: 291-303. 10.1002/pds.1200.View ArticlePubMedGoogle Scholar
- Benedetti A, Abrahamowicz M: Using generalized additive models to reduce residual confounding. Stat Med. 2004, 23: 3781-3801. 10.1002/sim.2073.View ArticlePubMedGoogle Scholar
- McCandless LC, Gustafson P, Levy A: Bayesian sensitivity analysis for unmeasured confounding in observational studies. Stat Med. 2007, 26: 2331-2347. 10.1002/sim.2711.View ArticlePubMedGoogle Scholar
- Suissa S: Immortal time bias in observational studies of drug effects. Pharmacoepidemiol Drug Saf. 2007, 16: 241-249. 10.1002/pds.1357.View ArticlePubMedGoogle Scholar
- Rothman KJ, Greenland S: Modern Epidemiology. 1998, Philadelphia, PA: Lippincott Williams & Wilkins, 125-2Google Scholar
- Suissa S, Garbe E: Primer: administrative health databases in observational studies of drug effects – advantages and disadvantages. Nat Clin Pract Rheumatol. 2007, 3: 725-732. 10.1038/ncprheum0652.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1741-7015/6/2/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.