- Systematic Review Protocol
- Open Access
The effects of exotic seaweeds on native benthic assemblages: variability between trophic levels and influence of background environmental and biological conditions
© Bulleri et al.; licensee BioMed Central Ltd. 2012
Received: 9 February 2012
Accepted: 13 June 2012
Published: 23 July 2012
Biological invasions are among the most severe threats to marine biodiversity. The impacts of introduced seaweeds on native macroalgal assemblages have been thoroughly reviewed. In contrast, no attempt has been made to synthesize the available information on the effects of exotic seaweeds on other trophic levels. In addition, it has not been clarified whether the effects of introduced seaweeds on native assemblages vary according to background physical and biological conditions.
This protocol provides details of our proposed method to carry out a systematic review aiming to identify and synthesize existing knowledge to answer the following primary questions: a) how does the impact of the presence of exotic seaweeds on native primary consumers (across trophic levels) compare in magnitude and extent to that observed on native primary producers (same trophic level)?; b) does the intensity of the effects of the presence of exotic seaweeds on native benthic ecosystems vary along a gradient of human disturbance (i.e. from urban/industrial areas to extra-urban areas to pristine areas)?
Biological invasions are an important component of global change, posing major threats to marine biodiversity . The introduction of non-indigenous species can alter the structure and functioning of natural ecosystems, with potential repercussions for their ability to provide goods and services to humans [2, 3]. Thus, assessing the impacts of introduced species on native assemblages is among the top priorities for ecologists.
In the marine environment, the spread of introduced benthic macroalgae can lead to their complete domination of substrata, strongly affecting native assemblages and biodiversity. The impacts of introduced seaweeds on native macroalgal assemblages have been thoroughly reviewed [4, 5]. In contrast, no attempt has been made to synthesize the available information on the effects of exotic seaweeds on other trophic levels. This is at odds with mounting evidence indicating that exotic seaweeds can be consumed by guilds of native herbivores, which, in some cases, receive physiological damage through the ingestion of these novel sources of food [6–8]. For instance, the red pigment caulerpin, the most abundant secondary metabolite of the green exotic macroalga, Caulerpa racemosa, enters food chains and accumulates in the fish tissues . Significant correlations among caulerpin tissue load, fish condition factor and hepatosomatic index suggest a possible detrimental effect of the dietary exposure to C. racemosa on Diplodus sargus. Under such circumstances, evaluating the magnitude of the impacts of introduced macroalgae at upper trophic levels (i.e. primary consumers) and how they relate to impacts at the same trophic level (i.e. on native macroalgal assemblages) is of paramount importance for the management of invasions. Addressing bottom-up effects associated with invasions of primary producers will also help to refine current invasion theories that envision the success of invasive species in the introduced range as the result of the loss of consumers or parasites (e.g. enemy-release hypothesis ).
In addition, we can predict that the effects of introduced seaweeds on native assemblages will vary according to background physical and biological conditions. Biological (e.g. richness, relative abundance and evenness of species) and physical attributes (e.g. availability of resources; disturbance regimes) of recipient ecosystems are known to play a key role in determining the susceptibility to invasion of natural communities [10–13]. These same attributes may also determine the severity of the impacts of invaders that become established. In stressful environments, species are often limited by unfavorable environmental conditions, rather than competition . Thus, we can predict that species in heavily degraded environments will not be limited by competition and are likely to be weakly affected by the presence of novel species. In well preserved, species-rich environments, competition is expected to be stronger, and resource uptake by exotic seaweeds could generate large negative effects on resident species.
Here, we propose to fill these two major gaps in our understanding of the impacts of biological invasions, taking advantage of the availability of a large body of empirical studies that have been carried out in marine environments.
Objective of the review
How does the impact of exotic seaweeds on native primary consumers (across trophic levels) compare in magnitude and extent to that observed on native primary producers (same trophic level)?
Do the effects exotic seaweeds on native benthic ecosystems vary in intensity along a gradient of human disturbance (i.e. from urban/industrial areas to extra-urban areas to pristine areas)?
Definitions of the components of primary question a
Outcomes (response variables)
Native seaweed and herbivore assemblages
Introduction, establishment and spread of exotic seaweeds
Species/taxa diversity (e.g richness, Simpson/Shannon/Berger-Parker diversity indices), species/taxa evenness (e.g. Pielou evenness index), total community cover, total community biomass, total biomass, mortality, density, single species/taxa cover, single species/taxa biomass, individual growth, individual size, community productivity, per capita reproductive output
Invaded vs non-invaded units Unit: organisms/individual/plots/ treatments/areas/sites/locations/regions Invaded unit vs Invader removal unit Unit: organisms/individual/plots/ treatments/areas/sites/locations/regions Control unit (no Invader) vs Invader transplanted unit Unit: organisms/individual/plots/ treatments/areas/sites/locations/regions
Definitions of the components of primary question b
Outcomes (response variables)
Native benthic ecosystems
Introduction, establishment and spread of exotic seaweeds
Species/taxa diversity (e.g richness, Simpson/Shannon/Berger-Parker diversity indices), species/taxa evenness (e.g. Pielou evenness index), total community cover, total community biomass, single species/taxa cover, single species/taxa biomass, mortality, density, individual size, individual growth, per capita reproductive output, community productivity, nutrient cycling (e.g. nutrient uptake rate [N, P], nitrification rates), metabolism (e.g. individual/community metabolic rate “B”;, oxygen consumption, production of primary or secondary metabolites), decomposition (e.g. abundance of particulate organic carbon [POC] or nitrogen [PON]), carbon flux, respiration, sediment stabilisation (e.g. eroding pulse pressure, sediment mass erosion), sediment mixing (e.g. sediment mixing depth, bioturbation potential “BP”; ), resilience (transient differences in spatial/temporal patterns of abundance of single or total species/taxa or composition of assemblages in presence vs absence of the invasive species), temporal (in)stability (e.g. temporal variance in abundance of single species/taxa or total community; temporal (dis)similarity indices), resistance (lack of differences in spatial/temporal patterns of abundance of single or total species/taxa or composition of assemblages in presence vs absence of the invasive species), abundance of epiphytes, richness (i.e. number) of epiphytes, invasibility (presence or abundance of other invasive species)
Invaded unit vs non-invaded unit Unit: organisms/individual/plots/ treatments/areas/sites/locations/regions Invaded unit vs Invader removal unit Unit: organisms/individual/plots/ treatments/areas/sites/locations/regions Control unit (no Invader) vs Invader transplanted unit Unit: organisms/individual/plots/ treatments/areas/sites/locations/regions
The strength of consumer pressure and competitive interactions are known to vary with latitude (i.e. decrease towards higher latitudes; ). In addition, the number of established exotic species increases with latitude . Thus, we will ask whether variation in the intensity of the effects of exotic seaweeds on primary producers and primary consumers depends upon latitude. In addition, we ask if latitude influences the comparison in magnitude and extent of the effects of exotic seaweeds between primary producers (same trophic level) and primary consumers (across trophic level).
Environmental conditions vary with latitude (i.e. solar radiation, temperature, productivity). Thus, we will ask whether variation in the intensity of the effects of exotic seaweeds along disturbance gradients depends upon latitude.
Exotic seaweeds can be present at different levels of abundance. Thus, we will ask whether the intensity of the effects of exotic seaweeds varies according to their abundance.
Scope of search
We will search the following databases for relevant literature and data:
ISI Web of Science
Bibliographies of all reviews identified as relevant during assessment of their full text will be searched for further material. In the case of papers reporting incomplete information, we will contact the authors to ask whether they might help in the retrieval of missing information. In order to maximize efficient use of limited resources, grey literature (reports, theses and non-published research works) will not be included. However, a test for publication bias and the potential impact of excluding grey literature will be carried out. References to the grey literature reported by peer-reviewed studies will be examined and included when matching the standards required .
Strings of search terms
Alga*, macroalga* Seaweed*
Alien*, invasi*, introduced, allochthonous, nonindigenous, non-indigenous, “non native*”, non-native*, exotic*
marine, brackish, estuar*, coastal,shallow, sea*, aquatic, maritime, lagoon*, pelagic, benth*, demersal, shore*, intertidal, subtidal, ocean*, bay, cove,
Impact*, effect*, influence, consequence*, food-web*, “food web*”, ecosystem*, biomass, biodiversity,“biological diversity”, communit*, richness, diversity, abundance, evenness, cover, density, “reproductive capacity”, mortality, growth, assemblage*, producti*, decomposition, “nutrient cycl*”, oxygen, carbon, flux, respiration, “ecosystem matabolism”, “sediment stabilisation”, epiphyte*, “sediment mixing”, resilience, stability, resistance, invasibility,
An asterisk (*) indicates a ‘wildcard’, which allows the database or search engine to look for multiple words that have different endings, e.g. estuar* captures [estuary OR estuaries OR estuarine]. Quotation marks (“”) around two words restricts the search to where that phrase occurs.
For these searches, the individual key words will be combined using ‘Boolean’ operators (restricting or expanding the searches using ‘AND’ and ‘OR’ to maximise the efficiency of searching). The individual terms or phrases listed within each of the columns (Exposure1, Exposure 2, Population, Outcome) will be combined by ‘OR’ within each category. These sets (one from each category) will be combined by ‘AND’ to generate the final set of results. Records from each database will be then integrated to yield the final list of papers. At the beginning, we will include all of the outcomes reported in Table 3, but the outcome list may be revised after having checked what is available in the literature.
Records of the search strategy used will be maintained to ensure repeatability and transparency and modifications will be made to the review protocol as necessary.
Relevant studies will be selected through a 3-step process: 1) scanning article titles, 2) reading of the abstract and 3) reading of full text. This hierarchical process will speed up the retrieval of relevant studies whilst minimizing the probability of missing relevant information.
Studies will be deemed as relevant when including the subjects specified as Population, Exposure, Comparator and Outcome, as defined in Tables 1 and 2. We will begin by searching for studies examining a broad range of outcomes. Due to limited resources, we may later decide to exclude studies examining some of these outcomes. This decision will be based on the number of studies found for different outcomes, not on the apparent magnitude or direction of effects.
Only studies written in English will be considered, but there will be no restriction in regard to the date of publication or location. The decision not to consider studies written in languages other than English is justified by an efficient use of available resources. Information on the latitude, distance from major urban or industrial centers (either raw data or synthetic indexes of human pressures), and abundance of the exotic seaweed will be considered for the secondary questions. We will provide details of the reason why some studies have been rejected.
Both field and lab studies will be initially retained, as the “study setting” will be later considered as a possible effect-modifier (see Section Potential effect modifiers and reasons for heterogeneity).
Articles from meeting abstracts or book chapters will be assessed for relevance following the same criteria stated above. Original data or relevant information thereof will be sought directly from the authors.
The first two steps will be carried out by two investigators working together. Repeatability of the article selection process will be determined through the assessment of a subset of the literature database. Both investigators will work independently at the third step of the process (reading of full-text), and repeatability will be assessed via kappa analysis. If there are significant discrepancies in the assessments of the investigators, these will be discussed and the inclusion criteria will be amended for clarity, if necessary.
Potential effect modifiers and reasons for heterogeneity
It is likely that we will find additional sources of heterogeneity among studies. Outcomes may also depend on the:
– study type (lab, field, observational, experimental)
– scale of the study (spatial, temporal)
– design (replication, controls, confounding variables, BACI)
– execution (independence of samples, independence of treatments, randomization, sampling techniques and protocols)
– habitat type
– removal versus addition of the invader
– geographical region (Europe vs non Europe)
– morphological groups of invaders
Study quality assessment
Study type: Comparison between invaded and non-invaded plots / Comparison between invaded and invader removal plots/ Comparison between non-invaded and invader transplanted plots
Distance from major urban or industrial centers
Study type: lab*/ field (*studies done in the lab will be excluded for the primary question b)
Study approach: observational /experimental
Spatial scale - extent: <1 km2 / 1–100 km2 / spanning area >100 km2
Spatial scale – grain (size of plots): plots <1 m / plots >1 m / individual /sub-individual (i.e. portions of an individual)
Description of the hierarchical design – space
# of spatial scales included: 1 / 2 / >2
Temporal scale – extent
Temporal scale – grain
Description of the hierarchical design – time
Number of temporal scales included: 1 / 2 / >2
Design - replication: Unreplicated / Replicated
Design - controls/reference sites: None / Procedural controls / Unmanipulated controls / Both
Design – CI
how many Control sites?
Design - BA
how many times Before?
Design - BACI
if Beyond BACI: how many Control sites?
if Beyond BACI: how many Impacted sites?
if Beyond BACI: how many 'Before' times?
Design – confounding: Confounded/Not confounded
Execution - specify the sampling method
Execution - sample independence: Not independent/Probably independent/ Independent
Execution - randomisation (allocation of sampling units): Neither/Haphazard/ Random
Confounding variables: Certainly present/Likely present/Not present
Studies that are considered of insufficient quality will be excluded from subsequent analysis (see section on confounding variables). However, the results of the quality assessment will be made available for each study (including studies excluded due to poor quality) in an appendix to the final review. The impact of excluding and including studies from the analysis that are categorized as different quality will be assessed in the review.
Determining whether replication has been carried out appropriately will be somewhat subjective, but the following criteria will be considered:
– do the controls appear to be spatially/temporally independent of the affected areas?
– do the replicates appear to be independent of one another in space and time (are they interspersed geographically, spread sufficiently in time, etc.)
– are controls and affected areas sufficiently similar (consider, for example habitat and substrate type, degree of exposure, salinity, proximity to human activities).
If clear evidence is available in the paper that all of these criteria are met, then the controls should be considered valid, if there is evidence of failure to meet any one of these criteria, the controls should be considered invalid and if it is not possible to make a proper assessment based on the information provided, then the study should be classified as ‘unclear’ in this regard.
If the authors refer to the replicates being assigned randomly, or make reference to the use of a random number table, they will be classified as being randomly assigned. If the authors refer to the allocation as haphazard, or make reference to a procedure such as throwing a quadrat over their shoulder, allocation will be classified as 'haphazard'. All others methods will be classified as 'other' and described so that their susceptibility to bias can be assessed.
A descriptive synthesis of the number and type of studies retained will be made available.
We will record the presence of factors that may have caused the observed changes other than the presence of exotic seaweeds. These might include pulse or chronic sources of disturbance. For example, the occurrence of accidents due to human activities (e.g. oil spill, nuclear waste) or natural extreme events (e.g. storms; run-off due to heavy rains) during the execution of the study could mask the effects of exotic seaweeds. Details of these events will be taken into account to evaluate whether a study will have to be retained or rejected within the analysis. This, although clearly reported within our review, will be somewhat subjective.
For each paper, details of the confounding variables will be given when present. Confounding variables will be formally accounted for if possible; otherwise, studies will have to be rejected.
In the case of observational (i.e. not experimental) studies, we will carry out sensitivity tests according to the following categories used to score studies:
Control - 1 Impacted
≥ 2 Control - ≥ 1 Impacted
Before - 1 After
≥ 2 Before - ≥ 2 Impacted
Control - 1 Impacted : 1 Before – 1After
≥ 2 Control - ≥ 1 Impacted: 1 Before – 1After
Control - 1 Impacted : ≥ 2 Before - ≥ 2 Impacted
The meta-analysis will be run including all the studies and repeated separately for each category of study to check how the category of the study can influence the results of the analysis. Of course, these analyses will be carried out only for those categories for which a sufficient number of studies is available. Categories will be reported in the summary Excel spreadsheet.
Papers suffering from one or more of the following major flaws could be rejected:
– Design – confounding: confounded
– Execution – sample independence: not independent
– Execution – treatment independence: not independent
– Randomization: no randomization in the allocation of experimental units to different treatments or, in the case of observational studies, of control sites.
– Specific details of the reasons why a study has been deemed as suffering from one or more of the flaws mentioned above will be provided.
Data extraction strategy
We will extract data on the response variables listed within the Outcomes column in Table 1 and 2. These will include response variables describing the response at the level of native assemblages, populations and individuals. Specifically, we will extract means, measures of variability (standard errors, standard deviations, confidence limits) and sample sizes for plots or sites where the invasive species was present or absent. In addition we will extract information on invader abundance, and other potential effect-modifiers (see Section Potential effect modifiers and reasons for heterogeneity) and confounding variables (see Section Confounding variables).
We will also extract geographical coordinates. This will enable us to calculate the distance of each site from major urban or industrial centers, and use it as a proxy for human disturbance levels (relevant for primary question b). This will also allow us to assess the level of human pressure at each study site, by referring to the map of human impacts developed by Halpern et al. . In addition, data on latitudes will be relevant to secondary questions a) and b). Authors will be contacted and asked to provide geographical coordinates, when these cannot be easily extracted from the papers.
Data will be extracted from tables or graphs presented in the selected papers. Extraction from graphs will be carried out by means of image analysis software, such as IMAGEJ or DATATHIEF. If data cannot be directly extracted from papers, the authors of the original study will be approached by e-mail, and asked to provide either raw data or relevant information (e.g. means, standard deviation/variance, sample size). Papers will be excluded when authors do not provide necessary information. In such cases, and in cases where only partial data is available, detailed information will be provided.
Data synthesis and presentation
For each study examined, a synthetic Excel spreadsheet will be generated in order to promote availability of relevant information. Data will be analyzed by means of meta-analysis techniques. It is likely that we will use random- or mixed-effect models. However, technical aspects, such as effect size measure (e.g. response ratio, Hedges’ d), nature of the factor involved (fixed or random) will be dependent upon the type and quality of data. For instance, in the case in which we will use mixed models, the error structure of effect sizes might be modeled by means of fixed-effect model. Likewise, analyses might be run on effect sizes calculated for each of the relevant response variables reported by each experiment or on average effect sizes calculated using equal weight for all the individual effect size that are reported per experiment. These technical aspects have been widely described in the literature [23–26].
Results of the meta-analysis will be reported graphically using standard approaches. Forest plots will be used to represent effect sizes, and funnel plots and normal q-q plots will be used to assess whether there is any publication bias.
Sub-group analyses will be performed to assess the effects of potential effect-modifiers.
All authors participated in the design of the study. FB drafted the manuscript. All authors read and approved the final manuscript.
This research has received funding from the European Community’s Seventh Framework Programme (FP7/2007-2013) under Grant Agreement No. 266445 for the project Vectors of Change in Oceans and Seas Marine Life, Impact on Economic Sectors (VECTORS).
- Gray JS: Marine biodivesity: patterns, threats and conservation needs. Biodiv Conserv 1997, 6: 153–157. 10.1023/A:1018335901847View ArticleGoogle Scholar
- Williamson MH: Biological invasions. Chapman and Hall, London; 1996.Google Scholar
- Pimentel D, Lach L, Zuniga R, Morrison D: Environmental and economic costs of nonindigenous species in the United States. Bioscience 2000, 50: 53–65. 10.1641/0006-3568(2000)050[0053:EAECON]2.3.CO;2View ArticleGoogle Scholar
- Smith JE, Williams SL: A global review of the distribution, taxonomy and impacts of introduced seaweeds. Annu Rev Ecol Syst 2007, 38: 327–359. 10.1146/annurev.ecolsys.38.091206.095543View ArticleGoogle Scholar
- Thomsen MS, Wenberg T, Tuya F, Silliman BR: Evidence for impacts of nonindigenous macroalgae: a meta-analysis of experimental field studies. J Phycol 2009, 45: 812–819. 10.1111/j.1529-8817.2009.00709.xView ArticleGoogle Scholar
- Trowbridge CD, Todd CD: Host-plant change in marine specialist herbivores: ascoglossan sea slugs on introduced macroalgae. Ecol Monogr 2001, 71: 219–243. 10.1890/0012-9615(2001)071[0219:HPCIMS]2.0.CO;2View ArticleGoogle Scholar
- Box A, Deudero S, Sureda A, Blanco A, Alòs J, Terrados J, Grau AM, Riera F: Diet and physiological responses of Spondyliosoma cantharus (Linnaeus, 1758) to the Caulerpa racemosa var. cylindracea invasion. J Exp Mar Biol Ecol 2009, 380: 11–19. 10.1016/j.jembe.2009.08.010View ArticleGoogle Scholar
- Terlizzi A, Felline S, Lionetto MG, Caricato R, Perfetti V, Cutignano A, Mollo E: Detrimental physiological effects of the invasive alga Caulerpa racemosa on the Mediterranean white seabream Diplodus sargus. Aquat Biol 2011, 12: 109–117. 10.3354/ab00330View ArticleGoogle Scholar
- Kean RM, Crawley MJ: Exotic plant invasions and the enemy release hypothesis. Trends Ecol Evol 2002, 17: 164–170. 10.1016/S0169-5347(02)02499-0View ArticleGoogle Scholar
- Hobbs RJ, Huenneke LF: Disturbance, diversity, and invasion: implications for conservation. Cons Biol 1992, 6: 324–337. 10.1046/j.1523-1739.1992.06030324.xView ArticleGoogle Scholar
- Davis MA, Grime JP, Thompson K: Fluctuating resources in plant communities: a general theory of invisibility. J Ecol 2000, 88: 528–534. 10.1046/j.1365-2745.2000.00473.xView ArticleGoogle Scholar
- Stachowicz JJ, Fried H, Osman RW, Whitlatch RB: Biodiversity, invasion resistance, and marine ecosystem function: reconciling pattern and process. Ecology 2002, 83: 2575–2590. 10.1890/0012-9658(2002)083[2575:BIRAME]2.0.CO;2View ArticleGoogle Scholar
- Bulleri F, Balata D, Bertocci I, Tamburello L, Benedetti-Cecchi L: The seaweed Caulerpa racemosa on Mediterranean rocky reefs: from passenger to driver of ecological change. Ecology 2010, 91: 2205–2212. 10.1890/09-1857.1View ArticleGoogle Scholar
- Bertness M, Callaway RM: Positive interactions in communities. Trends Ecol Evol 1994, 9: 191–193. 10.1016/0169-5347(94)90088-4View ArticleGoogle Scholar
- Allen PA, Brown JH, Gillooly GF: Global biodiversity, biochemical kinetics and the energetic-equivalence rule. Science 2002, 297: 1545–1548. 10.1126/science.1072380View ArticleGoogle Scholar
- Solan M, Cardinale BJ, Downing AL, Engelhardt KAM, Ruesink JL, Srivastava DS: Extinction and ecosystem function in the marine benthos. Science 2004, 306: 1177–1180. 10.1126/science.1103960View ArticleGoogle Scholar
- Pennings SC, Silliman BR: Linking biogeography and community ecology: Latitudinal variation in plant-herbivore interaction strength. Ecology 2005, 86: 2310–2319. 10.1890/04-1022View ArticleGoogle Scholar
- Sax DF: Latitudinal gradients and geographic ranges of exotic species: implications for biogeography. J Biogeogr 2001, 28: 139–150.View ArticleGoogle Scholar
- Vilà M, Espinar JL, Hejda M, Hulme PE, Jarosik V, Maron JL, Pergl J, Schaffner U, Sun Y, Pysek P: Ecological impacts of alien plants: a meta-analysis of their effects on species, communities and ecosystems. Ecol Lett 2011, 14: 702–708. 10.1111/j.1461-0248.2011.01628.xView ArticleGoogle Scholar
- Underwood AJ: Experiments in ecology. Cambridge University Press, Cambridge; 1997.Google Scholar
- Winer BJ, Brown DR, Michels KM: Statistical principles in experimental 410 design. McGraw-Hill, New York; 1991.Google Scholar
- Halpern BS, Walbridge S, Selkoe KA, Kappel CV, Micheli F, D’Agrosa C, Bruno JF, Casey KS, Ebert C, Fox HE, Fujita R, Heinemann D, Lenihan HS, Madin EMP, Perry MT, Selig ER, Spalding M, Steneck R, Watson R: A global map of human impact on marine ecosystems. Science 2008, 319: 948–952. 10.1126/science.1149345View ArticleGoogle Scholar
- Gurevitch J, Hedges LV: Statistical issues in ecological meta-analysis. Ecology 1999, 80: 1142–1149. 10.1890/0012-9658(1999)080[1142:SIIEMA]2.0.CO;2View ArticleGoogle Scholar
- Hedges LV, Gurevitch J, Curtis P: The meta-analysis using response ratios in experimental ecology. Ecology 1999, 80: 1150–1156. 10.1890/0012-9658(1999)080[1150:TMAORR]2.0.CO;2View ArticleGoogle Scholar
- Osenberg CW, Sarnelle O, Cooper SD, Holt RD: Resolving ecological questions through meta-analysis: goals, metrics, and models. Ecology 1999, 80: 1105–1117. 10.1890/0012-9658(1999)080[1105:REQTMA]2.0.CO;2View ArticleGoogle Scholar
- Rosenberg MS, Adams DC, Gurevitch J: Metawin: Statistical software for meta-analysis. Sinauer Associates, Massachusetts; 2000.Google 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.