Competitive assembly of South Pacific invasive ant communities
© Lester et al; licensee BioMed Central Ltd. 2009
Received: 01 October 2008
Accepted: 24 January 2009
Published: 24 January 2009
The relative importance of chance and determinism in structuring ecological communities has been debated for nearly a century. Evidence for determinism or assembly rules is often evaluated with null models that randomize the occurrence of species in particular locales. However, analyses of the presence or absence of species ignores the potential influence of species abundances, which have long been considered of major importance on community structure. Here, we test for community assembly rules in ant communities on small islands of the Tokelau archipelago using both presence-absence and abundance data. We conducted three sets of analyses on two spatial scales using three years of sampling data from 39 plots on 11 islands.
First, traditional null model tests showed support for negative species co-occurrence patterns among plots within islands, but not among islands. A plausible explanation for this result is that analyses at larger spatial scales merge heterogeneous habitats that have considerable effects on species occurrences. Second, analyses of ant abundances showed that samples with high ant abundances had fewer species than expected by chance, both within and among islands. One ant species, the invasive yellow crazy ant Anoplolepis gracilipes, appeared to have a particularly strong effect on community structure correlated with its abundance. Third, abundances of most ant species were inversely correlated with the abundances of all other ants at both spatial scales. This result is consistent with competition theory, which predicts species distributions are affected by diffuse competition with suites of co-occurring species.
Our results support a pluralistic explanation for ant species abundances and assembly. Both stochastic and deterministic processes interact to determine ant community assembly, though abundance patterns clearly drive the deterministic patterns in this community. These deterministic patterns were observed at two spatial scales. Results indicate that abundance-based null models may be more sensitive in detecting non-random patterns in community assembly than species co-occurrences analyses.
Opinions on the processes governing community assembly are polarized [1–3]. One school of thought maintains that ecological communities are assembled deterministically according to 'assembly rules', generated by biotic interactions [4–6]. An opposing view asserts that communities are assembled stochastically by chance dispersal of species with life history characteristics suited to local environmental conditions [7–9]. The relative importance of chance and determinism in structuring ecological communities has been fiercely debated for nearly a century and we are far from a resolution to the debate.
Diamond  hypothesized that interspecific interactions result in community assembly rules. He suggested that some pairs of species never coexist, either by themselves or as part of a larger combination. Evidence for assembly rules are often evaluated with null models that randomize the occurrence of species in particular locales. However, analyses of the presence or absence of species ignore the potential influence of species abundances, which have long been considered of major importance on community structure [10–13]. Abundance effects may be particularly important in determining the distribution of invasive species . Less abundant invasive species likely have weaker ecological effects on community assembly than more abundant species. Nevertheless, the effects of species abundances have largely been ignored when testing for community assembly rules.
Patterns in community structure, and their underlying mechanisms, can differ among spatial scales , as can our understanding of the invasiveness of introduced species. For example, investigations over larger spatial scales often incorporate a variety of habitats, some of which may not be suitable for an invasive species and thus provide refuges for native species. Despite evidence that community assembly rules can be scale-dependent [16, 17], null model tests for community assembly rules are typically assessed at single scales .
Here, we test for community assembly rules in ant communities on islands in the South Pacific, by conducting abundance-based null model analyses at two spatial scales. Ant communities on Islands in the South Pacific are often composed of introduced species, which are an ideal group for the study of processes governing community assembly. Interspecific aggression , chemical warfare , dominance hierarchies and competitive displacement [21–23] have been documented in numerous ant species inhabiting different parts of the globe. In addition, many invasive ants are considered superior competitors . Most islands in the South Pacific are volcanic in origin, so their flora and fauna are derived from over-water dispersal. Unlike most other island colonists, most ant species have colonized islands comparatively recently, with the aide of human travellers, and the process of island colonization and ant community assembly continues . Because South Pacific ant communities are often comprised almost entirely of introduced species, they provide an ideal opportunity to study how ecological communities are assembled. Unfortunately, they might also provide a template for the future of many other ecological communities, as invasive species continue their global spread.
We censused ant communities in 39 plots on 11 islands and used null models to test for community assembly rules on two spatial scales (i.e. within and among islands). First, we randomised the presence of ant species among samples to test for negative co-occurrence patterns. Second, we randomised the presence of individual ants among samples to test whether ant species richness was inversely related to total ant abundances. Third, we randomised the presence of individual ants among samples to test whether the abundances of each ant species was inversely related to the abundance of all other ant species.
Study site, species and field sampling
Ant communities on islands in the Southern Pacific Ocean typically have few native species and are often comprised entirely of species that were brought their by human travellers . Only 28 ant species have been recorded on Tokelau and of those only two are Pacific endemics, with no species endemic to Tokelau .
Censuses were conducted during three successive visits to the archipelago (November 2002, November–December 2004 and June–July 2005). A standard sampling design was used throughout [29–31]. Survey plots consisted of 15 × 15 m quadrats located at least 40 m apart in forested areas. The number of quadrats per island ranged from between two and nine over the period of 2002–2005. Four to five pitfall traps were haphazardly situated throughout the plot and one-third filled with Gault's solution, which is an insect killing agent and preservative . Placement of traps was haphazard (rather than random) so as to avoid traps being placed beside or on individual ant nests, which could provide biased estimates of ant abundances. The pitfall traps were plastic containers 9 cm tall, 7.5 cm diameter at the top, tapered to 5 cm diameter at the bottom and placed flush with the ground surface. Traps were only left out for 24 hrs due to the speed at which they occasionally accumulated abundant ant species, such as Anoplolepis gracilipes. In total we placed a total of 184 pitfall traps, within 39 plots from 11 islands. Analyses at the "plot" scale refer to the pooled number and abundance of each species from all pitfall traps within each quadrat. Analyses at the "island" scale refer to the pooled number and abundance of each species from all quadrats within each island, sampled over the period of 2002–2005.
To test for negative co-occurrence patterns, we used standard analyses in the freely available software package EcoSim . In this analysis, there were 17 rows of species and either 11 or 39 columns for islands or plots, respectively. Negative co-occurrence patterns were quantified with three indices: the C-score index, the checkerboard score and the number of unique species combinations. The C-score is obtained by calculating the number of checkerboard units cu for each species pair. For example, at the island scale of our analyses: cu = (o i - s)(o j - s), where o i is the total number of islands occupied by species i, o j is the total number of islands occupied by species j, and s is the number of islands occupied by both species . A single C-score, which describes community-level species co-occurrences, is obtained by averaging cu values of all species pairs. The checkerboard score is obtained by tallying the number of species pairs that never co-occur with one another. The number of unique species combinations is a count of all species pairs that co-occurred with one another.
To test for non-random co-occurrence patterns, the observed values of each index were compared to values generated in 5000 simulations of the observed species × plot presence-absence matrix using fixed row and column sums. The sequential swap technique was used to generate random permutations of the observed matrix. Miklós & Podani  show that this method biases against support for non-random co-occurrence patterns and provides a more conservative test of the co-occurrence assembly rule. If communities are structured deterministically according to assembly rules, observed communities should have higher C-scores, fewer unique species combinations, and more checkerboard species pairs than expected under the null model. Among-island and within-island analyses were conducted separately.
Competition theory predicts that as the total number of individuals in an area increases, competition for resources leads to the exclusion of competitively subordinate species . We tested this prediction by evaluating whether areas with higher numbers of ants had lower numbers of ant species. To test for an inverse relationship between species richness and overall ant abundances, we began by obtaining expected values of species richness for each sample locale (i.e. plot and island). Expected values of species richness were obtained using rarefaction, which estimates species richness on a per-individual, rather than a per-area, basis [2, 36, 37]. To generate rarefaction curves we used the programme RAREFACT 1.0 . Different numbers of ants were randomly sampled from the total pool encountered during sampling using this computer simulation, which estimates species richness for a given number of individual ants sampled. Separate rarefaction curves were generated for each spatial scale. The total pool of ants sampled across the archipelago was used to calculate expected species richness values for each island and the total pool found on each island was used to calculate expected species richness values for each plot.
This application of rarefaction yields expected values of species richness if individuals were randomly distributed among samples. Differences between observed (O) and expected (E) species richness values (O - E) therefore indicate whether samples sites have higher or lower richness than expected by chance. If ants compete for resources, sites with greater total ant abundances will contain fewer species than expected, resulting in negative O - E values. Conversely, sites with lower overall ant abundances will contain more species than expected, resulting in positive O - E values. However, differences between observed and expected values increase passively with the magnitude of expected values. To remove this confounding effect we divided the difference by expected values [(O-E)E-1], which results in an unbiased estimates of deviations from expected species richness values.
General linear models were used to establish whether differences from expected species richness varied with total ant abundances. Standardised differences from expected species richness values were included in the model as a dependent variable. Total ant abundances were treated as a covariate. Separate analyses were conducted for each scale (i.e. plot or island). In the plot-scale analysis, island was included in the model as a random factor to account for the independence problem generated by multiple plots occurring on the same island. Data were log10 transformed to conform to normality and homoscedasticity assumptions at the island-scale. Data conformed to assumptions without transformation at the plot-scale. Analyses were conducted in SPSS .
A significant effect of expected ant abundances would suggest stochastic processes shape community assembly. Support for community assembly rules would require negative relationships between the observed abundance of each species and the summed abundance of all other ant species. Ant species was included as a random factor in the island-scale analysis to account for the independence problem associated with including more than one value for each species. Both ant species and island were included as random factors in the plot-scale analysis to account for the independence problem associated with including more than one value for each species and each island. Separate analyses were conducted at each scale using the generalized linear model procedure with a Poisson distribution in SPSS .
Ant species with abundances found in plots and islands.
Anoplolepis gracilipes (Smith)
Paratrechina longicornis (Latr.)
Paratrechina vaga (Forel)
Cardiocondyla nuda (Mayr)
Monomorium floricola (Jerdon)
Monomorium liliuokalanii Forel
Pheidole fervens Smith
Pheidole oceanica Mayr
Pheidole sexspinosa Mayr
Pheidole umbonata Mayr
Rogeria stigmatica Emery
Strumigenys sp. 1
Tetramorium bicarinatum (Nyl.)
Tetramorium lanuginosum Mayr
Tetramorium simillimum (Smith)
Tetramorium tonganum Mayr
Anochetus graeffei Mayr
Patterns of species co-occurrences were scale-dependent. No evidence for co-occurrence assembly rules were found at the analysis scale of among islands. The observed C-score (4.368) was not statistically different from null model expectations among islands (4.369, 0.006; mean, σ; P = 0.485). Similar results were obtained for the number of checkerboard species pairs (observed = 13, expected = 15.997, 9.985; P = 0.882) and the number of unique species combinations (observed = 11, expected = 10.999, 0.028; P = 0.998). In contrast, within-island analyses showed evidence for co-occurrence assembly rules. The observed c-score (41.007) was significantly greater than null model expectations (40.214, 0.151; P= 0.023). While the observed number of checkerboard species pairs (32) did not significantly differ from expected numbers (34.099, 11.918; P = 0.760), the observed number of unique species combinations (35) was lower than null model expectations (39.623, 0.376; P < 0.001).
The ant communities in the Tokelau Archipelago appear to be assembled deterministically. However, evidence for community assembly rules differed between abundance-based and community-based analyses. Species co-occurrences were distributed at random among islands, but evidence for negative co-occurrence patterns were observed within islands. The number of species present in samples at both spatial scales differed from randomized expectations; differences between observed and expected species richness values were negatively correlated with total ant abundances, indicating that areas with high ant abundances supported a reduced number of ant species. Population abundances of most ant species were correlated with both randomized expectations and the summed abundance of all other ant species. Therefore, consistent support for community assembly rules were found in abundance based analyses, but scale dependent results were found in analyses of species co-occurrences.
Co-occurrence analyses showed variable support for community assembly rules. Within islands, species tended to co-occur less frequently than expected by chance. Negative co-occurrence patterns are common in community ecology  and are commonly attributed to inter-specific competition. However, our results differ substantially from other analyses involving other invasive ants. Both red imported fire ants (Solenopsis invicta) and Argentine ants (Linepithema humile) have caused the disassembly of ant communities in other locales [17, 41]. Our study involved several invasive species including A. gracilipes, which is thought to have similar effects on ant communities as S. invicta and L. humile . Conversely, between island analyses indicated that species co-occurrences were distributed at random. This result is consistent with several other studies which found that community assembly rules can be scale-dependent [16, 42]. Spatial scale thus appears to have a strong influence on the outcome of co-occurrence patterns. While there may be any number of hypotheses for this result, a plausible explanation is that analyses at larger spatial scales merge heterogeneous habitats that have considerable effects on species occurrences. Different species may have different habitat requirements. Determining the appropriate scale for co-occurrence analysis is probably community and habitat specific, and represents a difficult task.
Species richness analyses showed that the number of species present in samples declined with the total number of ants present. This result is consistent with theoretical predictions concerning how groups of competing species interact. Competition from groups of similar species may increase local extinction rates and decrease local colonization rates [5, 43]. Islands and plots with large numbers of ants housed fewer species than expected by chance, whereas samples with smaller numbers of ants housed more species. Although several previous studies have documented similar patterns in island communities (see Burns ; and references therein) null model analyses that attempt to link community level patterns of abundance to species are rare. Therefore, the generality of the abundance-species richness relationship observed here remains unclear.
Rarefaction can overestimate species richness when spatial distributions are clumped . If species interactions (e.g., competition) cause the spatial segregation of species distributions, the number of species present in samples should be less than predicted by rarefaction. We found that deviations from expected richness were negatively related to total ant abundances, suggesting that densely populated areas are subject to a greater influence of deterministic processes.
Patterns of the abundance of most ant species were negatively associated with the abundances of all other species. This result is consistent with competition theory, which predicts species distributions are affected by diffuse competition with suites of co-occurring species . Patterns of ant abundances were also associated with null model expectations, indicating a stochastic component to ant community assembly. Random dispersal events and/or stochastic population dynamics may therefore help to determine ant abundance patterns, supporting recent work highlighting the effects of 'neutral' processes on spatial patterns of biodiversity . When both relationships are viewed jointly, overall results support a pluralistic explanation for ant population abundances. Under this view, both stochastic and deterministic processes interact to determine ant community assembly. A similar conclusion has been reached by others examining ants on islands (e.g. Cole ), though a novel feature of our results is how abundance clearly drives deterministic patterns in this community.
Why patterns indicative of deterministic processes were consistently observed in abundance-based analyses, yet inconsistently in analyses of species co-occurrences is unclear. Perhaps the qualitative nature of binary data was not sensitive enough to detect evidence for assembly rules at all scales. Suites of analyses across a variety of scales using both quantitative and qualitative data may provide more accurate tests for assembly rules.
One ant species, the yellow crazy ant Anoplolepis gracilipes, appeared to have a particularly strong effect on community structure. When present in low abundances it appeared to have little to no effect on ant communities. However, in high abundance, A. gracilipes was associated with reductions in the number of co-occurring species and their abundance. Similarly, on Christmas Island A. gracilipes has its strongest influence on the communities when in high densities . High abundances may help invasive ants such as A. gracilipes obtain resources through numerical superiority, but are also associated with increased aggression . Our experimental work on these islands has demonstrated that when in high abundance, A. gracilipes often dominates preferred resources and only those species with different foraging modes or food preference can co-occur with them [31, 47]. Spatial variation in the distribution of such dominant species also likely creates room for subordinate species, and the establishment of these subordinate species may have a large stochastic element . Whatever the mechanisms for co-occurrence or exclusion, our results are consistent with competition theory, which predicts that interspecific competition and community structuring is abundance-dependent.
Overall results showed that ant communities in the Tokelau archipelago are assembled deterministically. Although some support for stochastic processes was observed, most of our results are consistent with the hypothesis that competition structures community assembly, both within and among islands. Because South Pacific ant communities are often comprised entirely of introduced species, they may provide a template for the future, given the continuing global spread of invasive species.
We thank John Haywood, Nate Sanders, Dan Simberloff and three anonymous reviewers for commenting on this manuscript. This project was funded by the Royal Society of New Zealand Marsden Fund, The Pacific Conservation and Development Trust, the US National Geographic Society, and Victoria University of Wellington.
- Roughgarden J: The structure and assembly of communities. Perspectives in ecological theory. Edited by: Roughgarden J, May RM, Levin SA. 1989, Cambridge: Cambridge University Press, 203-226.
- Gotelli NJ, Graves GR: Null models in ecology. 1996, Washington DC: Smithsonian Press
- Crawley MJ: The structure of plant communities. Plant Ecology. Edited by: Crawley MJ. 1997, Oxford: Blackwell Scientific Publications, 475-531.
- Clements FE: Plant succession: an analysis of the development of vegetation. Publication of the Carnegie Institute, Washington. 1916, 242: 1-512.
- Diamond JM: Assembly of species communities. Ecology and evolution of communities. Edited by: Cody ML, Diamond JM. 1975, Cambridge: Harvard University Press, 342-444.
- Silvertown J, Poulton P, Johnston E, Edwards G, Heard M, Biss PM: The Park Grass Experiment 1856–2006: Its contribution to ecology. Journal of Ecology. 2006, 94: 801-814. 10.1111/j.1365-2745.2006.01145.x.View Article
- Gleason HA: The structure and development of the plant association. Bulletin of the Torrey Botanical Club. 1917, 44: 463-481. 10.2307/2479596.View Article
- Connor EF, Simberloff D: The assembly of species communities: chance or competition?. Ecology. 1979, 60: 1132-1140. 10.2307/1936961.View Article
- Hubbell SP: The unified neutral theory of biodiversity and biogeography. 2001, Princeton: Princeton University Press
- Lack D: The natural regulation of animal numbers. 1954, London: Oxford University Press
- Sutherland J: From Individual Behaviour to Population Ecology. 1996, Oxford: Oxford University Press
- Harte J, Conlisk E, Ostling A, Green JL, Smith AB: A theory of spatial structure in ecological communities at multiple spatial scales. Ecological Monographs. 2005, 75: 179-197. 10.1890/04-1388.View Article
- Ulrich W, Zalewski M: Abundance and co-occurrence patterns of core and satellite species of ground beetles on small lake islands. Oikos. 2006, 114: 338-348. 10.1111/j.2006.0030-1299.14773.x.View Article
- Colautti RI, MacIsaac HJ: A neutral terminology to define 'invasive' species. Diversity and Distributions. 2004, 10: 135-141. 10.1111/j.1366-9516.2004.00061.x.View Article
- Levin SA: The problem with pattern and scale in ecology. Ecology. 1992, 73: 1943-1967. 10.2307/1941447.View Article
- Gotelli NJ, Ellison AM: Assembly rules for New England ant assemblages. Oikos. 2002, 99: 591-599. 10.1034/j.1600-0706.2002.11734.x.View Article
- Sanders NJ, Gotelli NJ, Heller NE, Gordon DM: Community disassembly by an invasive species. Proceedings of the Natural Academy of Sciences, USA. 2003, 100: 2474-2477. 10.1073/pnas.0437913100.View Article
- Burns KC: Patterns in the assembly of an island plant community. Journal of Biogeography. 2007, 34: 760-768. 10.1111/j.1365-2699.2006.01625.x.View Article
- Cole BJ: Assembly of mangrove ant communities: patterns of geographical distribution. Journal of Animal Ecology. 1983, 52: 339-347. 10.2307/4557.View Article
- Andersen AN, Blum MS, Jones TH: Venom alkaloids in Monomorium "rothsteini" Forel repel other ants: is this the secret to success by Monomorium in Australian ant communities. Oecologia. 1991, 88: 157-160. 10.1007/BF00320805.View Article
- Holway DA: Effect of Argentine ant invasions on ground-dwelling arthropods in northern California riparian woodlands. Oecologia. 1998, 116: 252-258. 10.1007/s004420050586.View Article
- Palmer TM, Stanton ML, Young TP, Goheen JR, Pringle RM, Karban R: Breakdown of an ant-plant mutualism follows the loss of large herbivores from an African Savanna. Science. 2008, 319: 192-195. 10.1126/science.1151579.View ArticlePubMed
- Savolainen P, Vespäläinen K: Niche differentiation of ant species within territories of the wood ant Formica polyctena. Oikos. 1989, 56: 3-16. 10.2307/3566082.View Article
- Holway DA, Lach L, Suarez AV, Tsutsui ND, Case TJ: The causes and consequences of ant invasions. Annual Review of Ecology and Systematics. 2002, 33: 181-233. 10.1146/annurev.ecolsys.33.010802.150444.View Article
- Wilson EO, Taylor RW: The ants of Polynesia (Hymenoptera: Formicidae). Pacific Insects Monograph. 1967, 14: 1-109.
- Mueller-Dombois D, Fosberg FR: Vegetation of the tropical Pacific islands. 1998, New York: Springer-VerlagView Article
- Parham BEV: The vegetation of the Tokelau Islands with special reference to the plants of Nukunonu Atoll. New Zealand Journal of Botany. 1971, 9: 576-609.View Article
- Abbott KL, Sarty M, Lester PJ: The ants of Tokelau. New Zealand Journal of Zoology. 2006, 33: 157-164.View Article
- Lester PJ, Tavite A: Long-legged ants (Anoplolepis gracilipes) have invaded Tokelau, changing the composition and dynamics of ant and invertebrate communities. Pacific Science. 2004, 58: 391-401. 10.1353/psc.2004.0031.View Article
- Abbott KL, Greaves SNJ, Ritchie PA, Lester PJ: Behaviourally and genetically distinct populations of an invasive ant provide insight into invasion history and impacts on a tropical ant community. Biological Invasions. 2007, 9: 453-463. 10.1007/s10530-006-9052-2.View Article
- Sarty M, Abbott KL, Lester PJ: Community level impacts of an ant invader and food mediated coexistence. Insectes Sociaux. 2007, 54: 166-173. 10.1007/s00040-007-0927-8.View Article
- Walker AK, Crosby TK: The preparation and curation of insects. 1988, Auckland: DSIR Information Series 163 Entomology Division
- Gotelli NJ, Entsminger GL: EcoSim: Null models software for ecology. Version 7. 2007, Acquired Intelligence Inc. & Kesey-Bear. Jericho, VT 05465, USA, [http://garyentsminger.com/ecosim/index.htm]
- Stone L, Roberts A: Competitive exclusion, or species aggregation? An aid in deciding. Oecologia. 1992, 91: 419-424. 10.1007/BF00317632.View Article
- Miklós I, Podani J: Randomization of presence-absence matrices: comments and new algorithms. Ecology. 2004, 85: 86-92. 10.1890/03-0101.View Article
- Sanders HL: Marine benthic diversity: a comparative study. American Naturalist. 1968, 102: 243-282. 10.1086/282541.View Article
- Gotelli NJ, Colwell RK: Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecology Letters. 2001, 4: 379-391. 10.1046/j.1461-0248.2001.00230.x.View Article
- RAREFACT 1.0. [http://people.hofstra.edu/j_b_bennington/research/paleoecology/rarefact.html]
- SPSS: SPSS for Windows, version 16.0.1. 2007, Chicago: Lead Technologies;
- Gotelli NJ, McCabe DJ: Species co-occurrence: A meta-analysis of J.M. Diamond's assembly rules model. Ecology. 2002, 83: 2091-2096.View Article
- Gotelli NJ, Arnett AE: Biogeographic effects of red fire ant invasion. Ecology Letters. 2000, 3: 257-261. 10.1046/j.1461-0248.2000.00138.x.View Article
- Sanders NJ, Gotelli NJ, Wittman SE, Ratchford JS, Ellison AM, Jules ES: Assembly rules for ant communities across spatial scales and habitats. Journal of Biogeography. 2007, 34: 1632-1641. 10.1111/j.1365-2699.2007.01714.x.View Article
- Keddy PA: Competition. 1989, Chapman and Hall, London, UKView Article
- Collins MD, Simberloff D: Rarefaction and non-random spatial dispersion patterns. Environmental and Ecological Statistics. 2008
- O'Dowd DJ, Green PT, Lake PS: Invasional 'meltdown' on an oceanic island. Ecology Letters. 2003, 6: 812-817. 10.1046/j.1461-0248.2003.00512.x.View Article
- Sagata K, Lester PJ: Behavioural plasticity associated with propagule size, resources, and the invasion success of the Argentine ant Linepithema humile. Journal of Applied Ecology. 2009, 46: 19-27. 10.1111/j.1365-2664.2008.01523.x.View Article
- Sarty M, Abbott KL, Lester PJ: Habitat complexity facilitates coexistence in a tropical ant community. Oecologia. 2006, 149: 465-473. 10.1007/s00442-006-0453-9.View ArticlePubMed
- Andersen AN: Not enough niches: non-equilibrial processes promoting species coexistence in diverse ant communities. Austral Ecology. 2008, 33: 211-220. 10.1111/j.1442-9993.2007.01810.x.View Article
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.