The evolution of antibiotic susceptibility and resistance during the formation of Escherichia colibiofilms in the absence of antibiotics
- Jabus G Tyerman†1, 2, 6,
- José M Ponciano†3,
- Paul Joyce2, 4, 5,
- Larry J Forney1, 2 and
- Luke J Harmon1, 2Email author
© Tyerman et al.; licensee BioMed Central Ltd. 2013
Received: 18 June 2012
Accepted: 11 January 2013
Published: 28 January 2013
Explanations for bacterial biofilm persistence during antibiotic treatment typically depend on non-genetic mechanisms, and rarely consider the contribution of evolutionary processes.
Using Escherichia coli biofilms, we demonstrate that heritable variation for broad-spectrum antibiotic resistance can arise and accumulate rapidly during biofilm development, even in the absence of antibiotic selection.
Our results demonstrate the rapid de novo evolution of heritable variation in antibiotic sensitivity and resistance during E. coli biofilm development. We suggest that evolutionary processes, whether genetic drift or natural selection, should be considered as a factor to explain the elevated tolerance to antibiotics typically observed in bacterial biofilms. This could be an under-appreciated mechanism that accounts why biofilm populations are, in general, highly resistant to antibiotic treatment.
Bacteria that form biofilms have been shown to be highly resistant to antimicrobial therapy [1–3] and contribute to the chronic nature of many bacterial infections in humans because cells in biofilms are highly resistant to antibiotic treatment (e.g. [4, 5]). Developing effective treatments for biofilm-related infections requires an understanding of the processes that lead biofilms to persist in the face of antimicrobial treatment .
There are multiple hypotheses to explain biofilm persistence during antibiotic treatment [7–10]. Physical factors, like diffusion limitation, may prevent antibiotic concentrations from reaching inhibitory or lethal levels within biofilms [11, 12]. However, several studies report biofilm persistence despite substantial antibiotic diffusion (e.g. [13, 14], reviewed in [1, 15]). Another hypothesis posits that non-genetic phenotypic heterogeneity, including the plastic expression of phenotypes that are resistant to antibiotics, may occur in biofilms. For example, resource gradients may lead cells to experience different microenvironments and thus express distinct phenotypes in different parts of a biofilm [1, 9]. In particular, a “persister” phenotype has been hypothesized to arise in biofilms and provide immunity against antibiotics . These persister cells are postulated to have the ability to reconstitute biofilms upon release from antibiotic threat [16–19].
Alternatively, bacteria cultured as biofilms may evolve heritable variation for resistance to antibiotics de novo. We know that evolutionary change can occur rapidly within tens to hundreds of generations , which are time scales relevant to medical treatment of infectious diseases . Furthermore, bacteria in biofilms have huge population sizes so that many new mutations will arise over relatively short time scales. It is possible that antibiotic resistance might arise in bacterial biofilms through straightforward population genetic processes. We suggest that, given enough time, variation in antibiotic resistance may arise in biofilms even in the absence of antibiotic selection. This could happen due to the accumulation of neutral variation or as a result of selection for phenotypes that by happenstance are correlated with antibiotic resistance (e.g. ).
Evolutionary change in biofilms is plausible since recent studies have shown tremendous phenotypic variation among cells isolated from biofilms growing in patients [23–25] and demonstrated heritable variation in traits within experimentally cultured biofilm populations [26–32]. However, few hypotheses about antibiotic resistance in biofilms invoke evolutionary change as an alternative explanation for biofilm persistence during antibiotic treatment (but see [7, 8, 27]). Evolution could involve mutations that convey resistance to single antibiotics (specialized resistance) or to whole suites of antibiotics (broad-spectrum resistance ). The occurrence of variants resistant to antibiotics may provide an “insurance effect” [30, 32, 33] by creating subpopulations of cells that can survive or even proliferate should the biofilm come under antibiotic assault. Resistant variants could thus facilitate reconstitution of bacterial biofilm populations following cessation of antibiotic treatment.
Documenting whether genetic variation for antibiotic resistance arises during the course of biofilm development is an important first step to exploring biofilm persistence in the face of antibiotic treatment. As we suggest, evolutionary explanations predict that genetic variation in the susceptibility to antibiotics will arise in biofilms, and that the frequency of antibiotic resistant cells will increase through time. Here we demonstrate the rapid evolution of heritable variation for broad-spectrum antibiotic resistance during the course of biofilm development by E. coli. We tested three hypotheses in this work: (1) heritable variation for antibiotic resistance evolves during biofilm development; (2) this variation includes both resistant and susceptible mutants to a range of antibiotics; and (3) phenotypic variation in biofilms increases through time. We found evidence for both (1) and (2) but not (3). Variation in antibiotic resistance emerged within 15 days of biofilm growth, a time frame that is consistent with many common bacterial infections. Thus, these findings have important implications for the development of treatments for bacteria that form biofilms during infection.
Results and discussion
Summary of the antibiotics used in this study
Ampicillin with Sulbactam
Semi-synthetic (rifamycin group)
Cationic basic protein
Results of planned contrasts following one-way MANOVA on mean ZOI across biofilm replicates and time
Antibiotics that differ significantly
Ancestor vs. 15 days
Ancestor vs. 30 days
CFP75↓, GM10↓, PB300↓, SAM20↓, TE30↓
Ancestor vs. 60 days
15 days vs. 30 days
15 days vs. 60 days
CIP5↓, E15↓, GM10↓, K30↓, NA30↓, PB300↓, S10↓, TE30↓
30 days vs. 60 days
CIP5↓, GM10↓, K30↓, PB300↓, S10↓, SAM20↓, TE30↓
Biofilm 1 vs. biofilm 2
Biofilm 1 vs. biofilm 3
Biofilm 2 vs. biofilm 3
Our results suggest that antibiotic resistance and susceptibility can rapidly evolve in biofilms over relatively short time scales (<15 days), which begs the question of how these rates of mutation accumulation compare to those in well-mixed liquid cultures that support exponential growth. Such a comparison is difficult because estimating the “mutation rate” in the spatially structured bacterial cells of biofilms is problematic. Mutation rates are almost always calculated and compared on a “per-generation” basis (e.g. [35, 36]), but rates of bacterial cell division in biofilms vary widely depending on location within the biofilm matrix. This variation in cellular growth rates is a consequence of nutrient depletion and the creation of strong gradients of substrates, electron acceptors and other resources within the spatially structured environment of biofilms [37, 38]. These gradients cause growth rates to vary tremendously within biofilms, such that cells deep within the biofilm matrix may not divide at all . Because of this, mutation frequency cannot be expressed in the same terms, i.e., per generation, as in well-mixed liquid cultures; nor can one calculate a meaningful population-wide average growth rate for cells in biofilms. One can imagine applying models that account for differential growth in biofilms (e.g., ), and then using current data to calculate mutation rates that can be compared to rates in well-mixed cultures. However, such calculations require data about mutation rates in in non-growing bacterial cells that is largely lacking, so direct and simple comparisons between biofilms and well-mixed cultures are not possible at this time.
The evolution of antibiotic resistance and susceptibility in bacterial biofilms involves the interaction between mutation, selection, genetic drift, and spatial structure [26, 40]. The data presented here cannot determine the importance of these multiple explanatory factors. It seems likely that evolution in biofilms typically occurs under conditions contrary to what is typically assumed in standard population genetics theory (e.g. strong selection and weak mutation) and rather involves strong mutational mechanisms typical in bacteria under stress  coupled with weak selection (see also ). Future work combining spatially explicit models for biofilm growth (e.g. [40, 42]) with model-based estimates of mutation rates and effect sizes for bacteria (e.g. ) would provide more insight into the details of evolution in biofilms.
These data show the rapid de novo evolution of heritable variation in antibiotic sensitivity and resistance during E. coli biofilm development. We suggest that evolutionary processes, whether genetic drift or natural selection, should be considered as a factor to explain the elevated tolerance to antibiotics typically observed in bacterial biofilms. We do not yet know whether evolution of antibiotic resistance requires high rates of mutation as can arise in biofilms (e.g. ) or can be explained by normal mutation rates in bacteria. In either case, biofilms quickly evolve high levels of variation in antibiotic resistance. We hypothesize that rare, highly resistant variants may allow biofilms to regrow following antibiotic treatment. This mechanism is an important potential explanation for why biofilm populations are, in general, highly resistant to antibiotic treatment.
Strain, media and growth conditions
Bacterial biofilms were grown as described by Ponciano et al.. Escherichia coli K12 MG1655 was grown in minimal salts media (M9) augmented with vitamins and trace elements with 0.05% glucose as the carbon source. The inoculum for flow cells was prepared by inoculating 10 ml of minimal medium with a scraping from a -80°C freezer stock and incubating the culture for 24 h at 37°C. Biofilms were cultured in flow-cells that had been sterilized by flowing 5% bleach for >24 h, followed by rinsing with minimal medium for 24 h. A 100 ul inoculum was introduced to each flow cell using a sterile syringe and needle. Bacterial cells were allowed to settle for 6 h before the flow was restarted (with a mean hydraulic retention time of 2.5 h). The biofilms were cultured for 15, 30 or 60 days prior to sampling them through a port on the upper surface of the flow-cell using a syringe and needle. Each sample was vortexed for 1 minute, then serially diluted in saline and plated on minimal medium solidified with agar. Ten randomly chosen clones were obtained from each of three replicate biofilms sampled at four times: 0, 15, 30 or 60 days. The 0d samples are referred to as “ancestors”. All clones were grown overnight in minimal medium and archived as glycerol stock cultures at -80°C.
We determined the sensitivity of ancestral and biofilm clones (15, 30 and 60 days old) to 12 antibiotics using the Kirby-Bauer disk method. The antibiotics (Table 1) were selected to target a range of cellular processes. Individual clones were grown in minimal medium for 24 h (final optical density at 600 nm = 0.15-0.2) and spread on Mueller-Hinton agar using sterile cotton swabs to form a lawn. After allowing the plate to dry for about 10 minutes, antibiotic-infused disks (Sensi-Disks, BD, New Jersey) were placed on the plates, which were then incubated for 18 h at 37°C. We photographed plates from a standard distance and measured the zones of (growth) inhibition (ZOI) for each antibiotic disk using ImageJ software available for download from NIH (http://rsbweb.nih.gov/ij/).
Antibiotic resistance was quantified as the diameter (in pixels) of the ZOI around the antibiotic-infused disks. Susceptible clones had relatively larger ZOI, while resistant clones had relatively smaller ZOI. For each bacterial clone, we replicated the resistance score for each antibiotic three times by growing three independent cultures from the frozen stock of that clone. For each independent replicate, we used the mean ZOI from two antibiotic disks. For each antibiotic disk, we scored the ZOI as the mean of three arbitrarily drawn diameters across the ZOI. Thus, each resistance score represents a mean of 3 × 2 × 3 = 18 individual measurements. Finally, at each time point we measured two replicates of the ancestor as a control. At each sampling time we standardized scores by dividing each by the mean score of the control to reduce variation introduced by day-to-day fluctuations in media (i.e., agar thickness, dryness, concentration, etc.).
To test the hypothesis that heritable variation for antibiotic resistance arose during biofilm development, we carried out a one-way MANOVA across all biofilms and time points simultaneously (Table 1). This MANOVA used ZOI diameter across all antibiotics as a response variable, and a concatenated variable of biofilm identity and time as predictor variable (the 10 levels of the predictor variable were then: ancestral line at 0 days, Biofilm replicate 1 at 15, 30, and 60 days, Biofilm replicate 2 at 15, 30, and 60 days, and biofilm replicate 3 at 15, 30, and 60 days). We also subjected these data to a set of nine planned contrasts (all paired comparisons of the ancestral population and all biofilms from 15, 30, and 60 days, as well as all paired comparisons between the three biofilm replicates; see Table 2). To account for an inflated Type I error associated with multiple comparisons, we computed the conservative simultaneous confidence intervals for each contrast .
We identified sensitive and resistant forms, respectively, as clones whose mean ZOI was more than two standard deviations above or below the mean of the ancestor. To test for increasing variation through time, we used linear regression to compare both total multivariate phenotypic variation [disparity] among clones and the number of resistant or sensitive clones to biofilm age.
All analyses were conducted using R (version 2.12.2 ).
We thank A. Spangler and M. Yarger for technical assistance, and H.J. La for advice on cultivation of biofilms. Research reported in this publication was supported by an Institutional Development Award (IDeA) to LJF from the National Institute of General Medical Sciences of the National Institutes of Health under grant number P20RR16448, a grand from the National Science Foundation (NSF DEB-0919499), and a grant from Proctor & Gamble Co
- Costerton JW, Stewart PS, Greenberg EP: Bacterial biofilms: a common cause of persistent infections. Science. 1999, 284: 1318-1322. 10.1126/science.284.5418.1318.PubMedView ArticleGoogle Scholar
- Costerton JW: Microbial ecology comes of age and joins the general ecology community. Proc Natl Acad Sci U S A. 2004, 49: 16983-16984.View ArticleGoogle Scholar
- Fux CA, Costerton JW, Stewart PS, Stoodley P: Survival strategies of infectious biofilms. Trends Microbiol. 2005, 13: 34-40. 10.1016/j.tim.2004.11.010.PubMedView ArticleGoogle Scholar
- Thomassen MJ, Boxerbaum B, Demko CA, Kuchenbrod PJ, Dearborn DG, Wood RE: Inhibitory effect of cystic fibrosis serum on pseudomonas phagocytosis by rabbit and human alveolar macrophages. Pediatr Res. 1979, 13: 1085-1088. 10.1203/00006450-197909000-00030.PubMedView ArticleGoogle Scholar
- Lang BJ, Arron SD, Ferris SW, Hebert PC, MacDonald NE: Multiple combination bactericidal antibiotic testing for patients with cystic fibrosis infected with multiresistant strains of Pseudomonas aeruginosa. Am J Respir Crit Care Med. 2000, 162: 2241-2245.PubMedView ArticleGoogle Scholar
- Smith EE, Buckley DG, Wu Z, Saenphimmachak C, Hoffman LR, D'Argenio DA, Miller SI, Ramsey BW, Speert DP, Moskowitz SM, Burns JL, Kaul R, Olson MV: Genetic adaptation by Pseudomonas aeruginosa to the airways of cystic fibrosis patients. Proc Natl Acad Sci U S A. 2006, 103: 8487-8492. 10.1073/pnas.0602138103.PubMed CentralPubMedView ArticleGoogle Scholar
- Conibear TC, Collins SL, Webb JS: Role of mutation in Pseudomonas aeruginosa biofilm development. PLoS One. 2009, 14: e6289-View ArticleGoogle Scholar
- Cloete TE: Resistance mechanisms of bacteria to antimicrobial compounds. Int Biodeter & Biodegrad. 2003, 51: 277-282. 10.1016/S0964-8305(03)00042-8.View ArticleGoogle Scholar
- Mah TC, O’Toole GA: Mechanisms of biofilm resistance to antimicrobial agents. Trends Microbiol. 2001, 9: 34-39. 10.1016/S0966-842X(00)01913-2.PubMedView ArticleGoogle Scholar
- Stewart PS, Costerton JW: Antibiotic resistance of bacteria in biofilms. Lancet. 2001, 358: 135-138. 10.1016/S0140-6736(01)05321-1.PubMedView ArticleGoogle Scholar
- Cheema MS, Rassing JE, Marriot C: The diffusion characteristics of antibiotics in mucus glycoprotein gels. J Pharm Pharmacol. 1986, 38: 53-10.1111/j.2042-7158.1986.tb14282.x.View ArticleGoogle Scholar
- Gordon CA, Hodges NA, Marriott C: Antibiotic interaction and diffusion through alginate and exopolysaccharide of cystic fibrosis-derived Pseudomonas aeruginosa. J Antimicrob Chemother. 1988, 22: 667-674. 10.1093/jac/22.5.667.PubMedView ArticleGoogle Scholar
- Anderl JN, Franklin MJ, Stewart PS: Role of antibiotic penetration limitation in Klebsiella pneumoniae biofilm resistance to ampicillin and ciprofloxacin. Antimicrob Agents Chemother. 2000, 44: 1818-1824. 10.1128/AAC.44.7.1818-1824.2000.PubMed CentralPubMedView ArticleGoogle Scholar
- Stewart PS, Rayner J, Roe F, Rees WM: Biofilm penetration and disinfection efficacy of alkaline hypochlorite and chlorosulfamates. J App Microbiol. 2001, 91: 525-532. 10.1046/j.1365-2672.2001.01413.x.View ArticleGoogle Scholar
- Stewart PS: Mechanisms of antibiotic resistance in bacterial biofilms. Int J Med Microbiol. 2002, 292: 107-113. 10.1078/1438-4221-00196.PubMedView ArticleGoogle Scholar
- Lewis K: Persister cells, dormancy and infectious disease. Nature Rev Microbiol. 2007, 5: 48-56. 10.1038/nrmicro1557.View ArticleGoogle Scholar
- Dhar N, McKinney D: Microbial phenotypic heterogeneity and antibiotic tolerance. Curr Opinion in Microbiol. 2007, 10: 30-38. 10.1016/j.mib.2006.12.007.View ArticleGoogle Scholar
- Lewis K: Riddle of biofilm resistance. Antimicrob Agents Chemother. 2001, 45: 999-1007. 10.1128/AAC.45.4.999-1007.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Roberts ME, Stewart PS: Modeling protection from antimicrobial agents in biofilms through the formation of persister cells. Microbiol. 2005, 151: 75-80. 10.1099/mic.0.27385-0.View ArticleGoogle Scholar
- Hairston NG, Ellner SP, Geber MA, Yoshida T, Fox JA: Rapid evolution and the convergence of ecological and evolutionary time. Ecol Lett. 2005, 8: 1114-1127. 10.1111/j.1461-0248.2005.00812.x.View ArticleGoogle Scholar
- Nesse RM, Stearns SC: The great opportunity: evolutionary applications to medicine and public health. Evolutionary Applications. 2008, 1: 28-48. 10.1111/j.1752-4571.2007.00006.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Wrande M, Roth JR, Hughes D: Accumulation of mutants in “ageing” bacterial colonies is due to growth under selection, not stress-induced mutagenesis. Proc Natl Acad Sci U S A. 2008, 105: 11863-11868. 10.1073/pnas.0804739105.PubMed CentralPubMedView ArticleGoogle Scholar
- Foweraker JE, Laughton CR, Brown DFJ, Bilton D: Phenotypic variability of Pseudomonas aeruginosa in sputa from patients with acute infective exacerbation of cystic fibrosis and its impact on the validity of antimicrobial susceptibility testing. J Antimicrobl Chemotherapy. 2005, 55: 921-927. 10.1093/jac/dki146.View ArticleGoogle Scholar
- Martínez-Solano L, Macia MD, Fajardo A, Oliver A, Martinez JL: Chronic Pseudomonas aeruginosa infection in chronic obstructive pulmonary disease. Clin Infect Dis. 2008, 47: 1526-1533. 10.1086/593186.PubMedView ArticleGoogle Scholar
- Mena A, Smith EE, Burns JL, Speert DP, Moskowitz SM, Perez JL, Oliver A: Genetic adaptation of Pseudomonas aeruginosa to the airways of cystic fibrosis patients is catalyzed by hypermutation. J Bacteriol. 2008, 190: 7910-7917. 10.1128/JB.01147-08.PubMed CentralPubMedView ArticleGoogle Scholar
- Ponciano JM, La H-J, Joyce P, Forney LJ: Evolution of Diversity in Spatially Structured E. coli populations. Appl Env Microbiol. 2009, 75: 6047-6054. 10.1128/AEM.00063-09.View ArticleGoogle Scholar
- Allegrucci M, Sauer K: Characterization of colony morphology variants isolated from Streptococcus pneumoniae biofilms. J Bac. 2007, 189: 2030-2038. 10.1128/JB.01369-06.View ArticleGoogle Scholar
- Bantinaki E, Kassen R, Knight CG, Robinson Z, Spiers AJ, Rainey PB: Adaptive divergence in experimental populations of Pseudomonas fluorescens. III. Mutational origins of wrinkly spreader diversity. Genetics. 2007, 176: 441-453. 10.1534/genetics.106.069906.PubMed CentralPubMedView ArticleGoogle Scholar
- Yarwood JM, Paquette KM, Tikh IB, Volper EM, Greenberg EP: Generation of virulence factor variants in Staphylococcus aureus biofilms. J Bac. 2007, 189: 7961-7967. 10.1128/JB.00789-07.View ArticleGoogle Scholar
- Boles BR, Singh PK: Endogenous oxidative stress produces diversity and adaptability in biofilm communities. Proc Natl Acad Sci U S A. 2008, 105: 12503-12508. 10.1073/pnas.0801499105.PubMed CentralPubMedView ArticleGoogle Scholar
- Boles BR, Thoendel M, Singh PK: Self-generated diversity produces “insurance effects” in biofilm communities. Proc Natl Acad Sci U S A. 2004, 101: 16630-16635. 10.1073/pnas.0407460101.PubMed CentralPubMedView ArticleGoogle Scholar
- Waite RD, Struthers JK, Dowson CG: Spontaneous sequence duplication within an open reading frame of the pneumococcal type 3 capsule locus causes high-frequency phase variation. Mol Microbiol. 2001, 42: 1223-1232.PubMedView ArticleGoogle Scholar
- Yachi S, Loreau M: Biodiversity and ecosystem productivity in a fluctuating environment: the insurance hypothesis. Proc Natl Acad Sci U S A. 1999, 96: 1463-1468. 10.1073/pnas.96.4.1463.PubMed CentralPubMedView ArticleGoogle Scholar
- Ortez JH: Disk diffusion testing. Manual of Antimicrobial Susceptibility Testing. Edited by: Coyle MB. 2005, Washington, DC: American Society for MicrobiologyGoogle Scholar
- Rosche WA, Foster PL: Determining Mutation Rates in Bacterial Populations. Methods. 2000, 20: 4-17. 10.1006/meth.1999.0901.PubMed CentralPubMedView ArticleGoogle Scholar
- Martinez JL, Baquero F: Mutation Frequencies and Antibiotic Resistance. Antimicrob Agents Chemother. 2000, 44: 1771-1777. 10.1128/AAC.44.7.1771-1777.2000.PubMed CentralPubMedView ArticleGoogle Scholar
- Stewart PS, Franklin MJ: Physiological heterogeneity in biofilms. Nat Rev Microbiol. 2008, 6: 199-210. 10.1038/nrmicro1838.PubMedView ArticleGoogle Scholar
- Stewart PW: Diffusion in biofilms. J Bacteriol. 2003, 185: 1485-1491. 10.1128/JB.185.5.1485-1491.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Perez-Osorio AC, Williamson KS, Franklin MJ: Heterogeneous rpoS and rhlR mRNA Levels and 16S rRNA/rDNA (rRNA Gene) Ratios within Pseudomonas aeruginosa Biofilms Sampled by Laser Capture Microdissection. J Bacteriol. 2010, 192: 2991-3000. 10.1128/JB.01598-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Eastman JM, Harmon LJ, La H-J, Joyce P, Forney LJ: The onion model, a simple neutral model for the evolution of diversity in bacterial biofilms. J Evol Biol. 2011, 11: 2496-2504.View ArticleGoogle Scholar
- Tenaillon O, Toupance B, Le Nagard H, Taddei F, Godelle B: Mutators, population size, adaptive landscape and the adaptation of asexual populations of bacteria. Genetics. 1999, 152: 485-493.PubMed CentralPubMedGoogle Scholar
- Chang I, Gilbert ES, Eliashberg E, Keasling JD: A three-dimensional, stochastic simulation of biofilm growth and transport-related factors that affect structure. Microbiology. 2003, 149: 2859-2871. 10.1099/mic.0.26211-0.PubMedView ArticleGoogle Scholar
- Perfeito L, Fernandes L, Mota C, Gordo I: Adaptive Mutations in Bacteria: High Rate and Small Effects. Science. 2007, 317: 813-815. 10.1126/science.1142284.PubMedView ArticleGoogle Scholar
- Barrick JE, Yu DS, Yoon SH, Jeong H, Oh TK, Schneider D, Lenski RE, Kim FJ: Genome evolution and adaptation in a long-term experiment with Escherichia coli. Nature. 2009, 461: 1243-1247. 10.1038/nature08480.PubMedView ArticleGoogle Scholar
- Johnson RA, Wichern DW: Applied Multivariate Analysis. 2002, New Jersey: Prentice HallGoogle Scholar
- Development Core Team R: R: A language and environment for statistical computing. 2012, Vienna, Austria: R Foundation for Statistical Computing, URL http://www.R-project.orgGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.