
==== Front
Environ Evid
Environ Evid
Environmental Evidence
2047-2382
BioMed Central London

257
10.1186/s13750-022-00257-z
Systematic Review Protocol
What are the effects of herbivore diversity on tundra ecosystems? A systematic review protocol
http://orcid.org/0000-0002-8120-5248
Barrio Isabel C. isabel@lbhi.is

1
http://orcid.org/0000-0002-3495-5839
Barbero-Palacios Laura 1
http://orcid.org/0000-0003-0014-0073
Kaarlejärvi Elina 2
http://orcid.org/0000-0002-0633-5595
Speed James D. M. 3
http://orcid.org/0000-0003-0329-6882
Heiðmarsson Starri 4
http://orcid.org/0000-0002-8994-9305
Hik David S. 5
http://orcid.org/0000-0003-4280-8350
Soininen Eeva M. 6
1 grid.432856.e 0000 0001 1014 8912 Faculty of Environmental and Forest Sciences, Agricultural University of Iceland, Árleyni 22, Keldnaholt, 112 Reykjavík, Iceland
2 grid.7737.4 0000 0004 0410 2071 Organismal and Evolutionary Research Programme, University of Helsinki, FI-00014 Helsinki, Finland
3 grid.5947.f 0000 0001 1516 2393 Department of Natural History, NTNU University Museum, Norwegian University of Science and Technology, 7491 Trondheim, Norway
4 grid.435368.f 0000 0001 0660 3759 Icelandic Institute of Natural History, Borgum við Norðurslóð, 600 Akureyri, Iceland
5 grid.61971.38 0000 0004 1936 7494 Department of Biological Sciences, Simon Fraser University, Burnaby, BC V5A 1S6 Canada
6 grid.10919.30 0000000122595234 Department of Arctic and Marine Biology, UiT – the Arctic University of Norway, 9037 Tromsö, Norway
31 1 2022
31 1 2022
2022
11 117 9 2021
12 1 2022
© The Author(s) 2022
2022
https://creativecommons.org/licenses/by/4.0/ Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
Background

Changes in the diversity of herbivore communities can strongly influence the functioning of northern ecosystems. Different herbivores have different impacts on ecosystems because of differences in their diets, behaviour and energy requirements. The combined effects of different herbivores can in some cases compensate each other but lead to stronger directional changes elsewhere. However, the diversity of herbivore assemblages has until recently been a largely overlooked dimension of plant–herbivore interactions. Given the ongoing environmental changes in tundra ecosystems, with increased influx of boreal species and changes in the distribution and abundance of arctic herbivores, a better understanding of the consequences of changes in the diversity of herbivore assemblages is needed. This protocol presents the methodology that will be used in a systematic review on the effects of herbivore diversity on different processes, functions and properties of tundra ecosystems.

Methods

This systematic review builds on an earlier systematic map on herbivory studies in the Arctic that identified a relatively large number of studies assessing the effects of multiple herbivores. The systematic review will include primary field studies retrieved from databases, search engines and specialist websites, that compare responses of tundra ecosystems to different levels of herbivore diversity, including both vertebrate and invertebrate herbivores. We will use species richness of herbivores or the richness of functional groups of herbivores as a measure of the diversity of the herbivore assemblages. Studies will be screened in three stages: title, abstract and full text, and inclusion will follow clearly identified eligibility criteria, based on their target population, exposure, comparator and study design. The review will cover terrestrial Arctic ecosystems including the forest-tundra ecotone. Potential outcomes will include multiple processes, functions and properties of tundra ecosystems related to primary productivity, nutrient cycling, accumulation and dynamics of nutrient pools, as well as the impacts of herbivores on other organisms. Studies will be critically appraised for validity, and where studies report similar outcomes, meta-analysis will be performed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13750-022-00257-z.

Keywords

Herbivore assemblage
Browsing
Grazing
Defoliation
Ecosystem function
Plant–herbivore interaction
Species richness
Icelandic Research Fund217754 Barrio Isabel C. http://dx.doi.org/10.13039/100010661 Horizon 2020 Framework Programme 869471 issue-copyright-statement© The Author(s) 2022
==== Body
pmcBackground

Herbivores in tundra ecosystems include organisms varying in size, from large ungulates, like muskoxen or reindeer, to small rodents, birds and invertebrates [1]. Arctic herbivores also vary in other traits that influence how they interact with plants and their environment, such as migratory behaviour, home range size, or their digestive physiology [2]. Given these differences, the diet choices and energy requirements of herbivores differ strongly [3, 4] and so do their impacts on ecosystems. Further, the effects of herbivores on processes, functions and properties of tundra ecosystems can be direct or indirect. For example, trampling by large herbivores can influence soil structure and soil biota directly, but also indirectly through changes in plant abundance and community composition [5]. To add to this complexity, the combined effects of different herbivore assemblages can lead to opposite outcomes (Fig. 1a). For example, in some cases the combined effects of large and small mammalian herbivores on vegetation are stronger than would be predicted for each group of herbivores alone [6–8]. In other cases, the effects of different herbivores can compensate each other if herbivores consume different plant species, leading to little to no changes in plant community composition [3, 9]. However, herbivore diversity has until recently been a largely overlooked dimension of plant–herbivore interactions in tundra ecosystems [10, 11], and we know little about how functionally different herbivore assemblages will affect these systems.Fig. 1 The diversity of the herbivore community (here expressed as richness of species or groups of herbivores) can have different impacts on processes, functions or properties of tundra ecosystems. a Taking as an example the effects of herbivores on vegetation, if the herbivores present in the herbivore community have similar diet we can expect additive effects on vegetation, where the combined effects of herbivores promote directional changes in plant community composition. On the other hand, if herbivores have different diets, we can expect compensatory effects, where the effects of different herbivores on plant species may balance each other by reducing competition among plant species and lead to stabilizing effects on plant community composition. b We hypothesize that herbivore diversity will influence ecosystem processes, functions or properties directly (e.g. consuming biomass, trampling) and indirectly (via effects on plant community composition as in panel a). Through both direct and indirect effects, we expect the combined effects of multiple herbivores on processes, functions and properties to be greater compared to the effects of single herbivores, although the slope of this relationship can change depending on the ecosystem process, function or property and the herbivore assemblage being considered. At very high herbivore diversity this relationship may saturate (dashed line) but we expect that in tundra ecosystems (shaded area), with relatively low numbers of species of herbivores, the relationship will not saturate

Tundra ecosystems are highly sensitive to environmental changes [12], and trophic dynamics are already responding to climatic changes [13, 14]. It has been suggested that the extent of ecosystem change in higher latitude systems at the end of the Pleistocene was determined by the loss of key species of herbivores [15], shifting in many cases vegetation from a grass-dominated steppe to a moss-dominated tundra [16]. Contemporary herbivore communities in tundra ecosystems are species-poor, compared to those present in the Pleistocene [17], but they can still drive ecosystem shifts [18, 19]. It has been hypothesized that a more diverse herbivore assemblage would have higher potential to influence ecosystem dynamics and drive transitions between vegetation states [17], so we are crucially in need of a better understanding of the status of herbivore assemblages and their role in tundra ecosystems. Knowledge on how the effects of different herbivores combine is thus particularly relevant to guide conservation and sustainable grazing management in tundra ecosystems that are grazed by wild and domestic herbivores [20].

We will use a systematic review to assess the effects of herbivore diversity on the functioning of tundra ecosystems. We will focus on multiple ecosystem processes, functions and properties related to primary productivity, nutrient cycling and nutrient pools, as well as the impacts of herbivores on other organisms, to comprehensively assess the role of herbivore diversity in the functioning of tundra ecosystems [21, 22]. By synthesizing the results of studies that evaluate the effects of different herbivore assemblages, we will be able to assess the influence of herbivore diversity on ecosystem functioning. As a measure of diversity, we will use the richness of species or groups of herbivores. We hypothesize that changes in the processes, functions or properties of tundra ecosystems will be affected by herbivore diversity (Fig. 1b). Specifically, we expect that increased diversity of herbivores will enhance ecosystem functioning through the complementarity of different herbivores, particularly in tundra ecosystems because they host relatively low herbivore diversity [2].

The topic for this systematic review was identified following the work on a systematic map [23] authored by a large number of scientists working on Arctic herbivory. These scientists represent the main stakeholder group for the topic of the present systematic review. Many of these scientists also work for conservation and environmental management agencies, thus broadening the scope to include other relevant stakeholders for the topic of the systematic review. Similar to the process followed by Soininen et al. [23], we will issue an open call for collaboration through the Herbivory Network (https://herbivory.lbhi.is/), to ensure openness and engagement of key participants in the systematic review.

Objective of the review

The main objective of the systematic review is to synthesize available evidence on the effects of herbivore diversity on tundra ecosystems. A recent systematic map on the effects of herbivory on Arctic vegetation [23] identified a substantial number of studies (98) that investigated the effect of several herbivores on vegetation responses, indicating that a systematic review on this topic is possible. In the present review, we will extend the effects of herbivores to include non-vegetation functions, processes and properties in tundra ecosystems. Thus, we expect the number of studies to be larger than those identified by the systematic map, although this number will in turn be limited by the studies from which we can extract data and by our definition of the comparator (see components of the primary question and eligibility criteria below).

Primary question. What are the effects of herbivore diversity on processes, functions and properties of tundra ecosystems?

Components of the primary questionPopulation terrestrial Arctic ecosystems

Exposure herbivory (including disturbance and fertilization effects of herbivores)

Comparator contrasting levels of herbivore diversity (species richness or richness of functional groups of herbivores)

Outcome measured ecological processes, functions and properties in response to herbivory

Methods

Searching for articles

The search string for this systematic review will be the same as in a previously published systematic map [23, 24]. That search string was broad enough to include all potentially relevant studies for the present systematic review because it focused on the region/system (terrestrial Arctic ecosystems) and exposure (herbivory), which are common to the systematic map and the present systematic review. Yet, the search string did not pose any restriction on the outcome or comparator. In contrast to the systematic map, where the outcomes included the responses of plants to herbivory and the comparator was no herbivory or alternative levels of herbivory [24], the present systematic review focuses on the responses of ecological processes, functions and properties to herbivory, where the comparator is different levels of herbivore diversity. The search for the systematic map was conducted in February 2019 and retrieved 3200 records, which were filtered based on the eligibility criteria defined for the systematic map. Since the outcomes of interest of this systematic review are broader and include also other ecosystem components and a different comparator, we defined new eligibility criteria for the systematic review (see eligibility criteria below).

The search string comprises two substrings, one targeted at delimiting the study region and the other targeted at the exposure element (see details in [22]). Our full search string (formatted for Web of Science) is:

(arctic OR subarctic OR tundra) AND (herbivor* OR graz* OR browser OR browsing OR grubb* OR trampl* OR defolia* OR ((invertebrate OR insect) AND (gall* OR mining OR miner)))

The search will be conducted in English in global search sources, and English together with relevant local languages (Russian, French, Finnish, Swedish, Norwegian, Icelandic, and Danish) in searches from local/regional sources (Additional file 1).

Publications will be searched in the following global search sources:Scopus (article title, abstract and keyword search) with no further limitations applied.

Web of Science (topic search), including all databases: Web of Science Core Collection, KCI-Korean Journal database, MEDLINE, Russian Science Citation Index and SciELO Citation Index.

Google Scholar: title search standardized so that search history is not taken into account. We will only include the first 300 search results as recommended by Haddaway et al. [25].

Details for institutional subscriptions for the final searches will be reported. In addition, we will search in local and specialist databases for grey literature using the list of sources from [23] (Additional file 1). We will complement the bibliographic database searches by checking the reference lists in relevant articles to identify other potentially relevant articles (i.e. “snowballing” [26]). We will apply no time or document type restrictions for the search, but we will exclude publications for which we cannot access the full text, either in electronic or paper form.

Assessing the specificity and sensitivity of the search

The specificity (minimizing the proportion of irrelevant studies returned by the search) and sensitivity (finding all relevant studies) of the search string for the systematic review was assessed in scoping exercises during protocol development [24]. Further, the comprehensiveness of the search was tested against an independent test-list of articles that the protocol development team identified as relevant to answer the question of the systematic review (Additional file 2). In May 2021, searches using the search string returned 2135 hits in the Scopus database, 3037 within Web of Science and ‘about 4650’ hits on Google Scholar. All of the articles included in the test-list (20 articles) were included in the results from this search. We screened the titles of the records retrieved from Scopus to check the specificity of the search string. Following the same procedure as in Soininen et al. [24] we subsampled the first 500 records sorted alphabetically by first author name. Based on their titles, we excluded 52% of the records, mainly because they were not focusing on terrestrial ecosystems (71% of excluded documents). This result is similar to the systematic map (46% of documents were excluded; [24]), where the specificity of the search string was deemed adequate, as no additional search terms would screen away studies conducted in non-terrestrial ecosystems.

We assessed the sensitivity of the search string by comparing the records retrieved in Scopus by the full search string (combining the two substrings: region- and exposure-specific) to the substrings separately. A total of 115,964 records were identified by the region-specific search string that were not included in the full search string, while 226,264 records identified by the exposure-specific search string were not retrieved by the full search string. From each of these, the first 1000 records ordered alphabetically by first author name were screened for relevance, based on title, abstract and full text, to identify potentially relevant records that could have been missed by the combined search string. None of these records were deemed as potentially relevant.

Article screening and study eligibility criteria

Screening process

Studies will be evaluated for inclusion in three stages: title, abstract and full text screening, based on the eligibility criteria presented below. When uncertain about whether a study should pass to the next stage of screening, reviewers will be inclusive. Reviewers will not take part in the critical appraisal of studies they have authored or co-authored. A list of studies excluded at full text stage with reasons for exclusion will be provided as an additional file to the systematic review.

During the development of the systematic review protocol, we tested the repeatability of the screening process at the abstract stage. Four co-authors (EMS, JDMS, LBP, ICB) screened the abstracts of 100 publications (from a list of alphabetically ordered records previously assessed as relevant by one of the authors based on their title). The observers unequivocally agreed on the classification of 58 studies, either for inclusion (15) or exclusion (43). From the remaining 42 studies, 32 corresponded to cases where three observers agreed and one disagreed, while 10 cases corresponded to cases where two observers suggested inclusion and two suggested exclusion. Of those 42 studies, 13 focused on a single herbivore species and it was not clear whether areas without herbivores were included in the study (i.e. no eligible comparator); 7 studies referred to simulated herbivory, 7 were paleo-ecological studies where the data on herbivores was unclear. In 6 cases observers disagreed on whether the studies assessed effects of herbivores, 2 studies focused on effects on the herbivore itself, three potentially referred to non-Arctic systems (e.g. using locality names that were identified by some reviewers as non-Arctic but not by others) and 4 had study-specific issues. These disagreements were discussed, and the eligibility criteria were refined. Based on these discussions and following the recommendations of Foo et al. [27] we built a decision tree to guide the screening process at the title and abstract stage and at the full text stage (Fig. 2). Two observers (ICB, LBP) screened 100 additional abstracts using this decision tree and their agreement in including or excluding documents increased from 82 to 93%.Fig. 2 Decision tree for the title, abstract and full text screening stages

In the systematic review, we will include an assessment of the repeatability of the screening process, where 10% of the articles at each stage (title, abstract and full text) will be assessed by two reviewers. Cohen’s Kappa (κ) statistic will be calculated to assess the consistency among reviewers [28]. After double screening 10% of the articles, reviewers will discuss and reconcile the disagreements and adjust the inclusion criteria accordingly. If deemed necessary, double screening of another subset of articles will be performed.

Eligibility criteria

Studies will be screened for eligibility based on their target population, exposure, comparator and outcome, as well as for the type of study design.

Eligible population (terrestrial Arctic ecosystems)

To be included, studies must focus on Arctic terrestrial ecosystems. As in Soininen et al. [23] we define the Arctic using the southern limit of the subarctic as defined by the Conservation of Arctic Flora and Fauna [29] to delimit the study region. To further limit the scope to tundra and ecotone forests, we will exclude studies in areas South of the subarctic zone, based on the definition used in the Arctic Biodiversity Assessment [30]. Studies conducted clearly in boreal forests or other non-Arctic terrestrial habitats will be excluded. We will extract geographical coordinates provided in the text, maps or place names, and we will exclude studies where it is not possible to extract this information. Articles reporting studies conducted at several locations including both Arctic and non-Arctic sites will be included, but only for the part conducted within the Arctic provided that it is possible to separate data from Arctic and non-Arctic sites; otherwise, the article will be excluded.

Eligible exposure (herbivory)

We will include studies that assess effects of terrestrial herbivores on tundra ecosystems (i.e., grazing, browsing, trampling and other types of disturbance, such as fertilizing, digging or grubbing). Studies have to report the identity of the herbivores, irrespective of taxonomic resolution; for example, studies may report the species of herbivore (e.g. Rangifer tarandus) or a broader group of herbivores (e.g. small mammalian herbivores). We will limit the scope to include only multicellular herbivores; we will thus exclude grazing by unicellular organisms. We will include studies simulating herbivory, as long as they clearly connect their experimental manipulation to the activities of an identifiable herbivore or group of herbivores; for example, we will include studies adding or removing reindeer faeces, studies clipping plant biomass simulating biomass consumption by lemmings and voles, or those mimicking grubbing disturbances by geese.

Eligible comparator (contrasting herbivore diversity)

As a measure of herbivore diversity, we will use the richness of herbivores in the herbivore assemblage, either as taxonomic units (species or subspecies) or broader groups of species (e.g. small mammals, geese). We acknowledge that richness is a coarse measure of diversity that only takes into account the presence or absence of species (or groups) and not their relative abundances. The inclusion of abundance data for herbivore populations would provide a more nuanced understanding of the role of herbivores in the functioning of ecosystems [31], but reliable information is limited to a few species and only for some regions [30]. To be included, studies must assess the effect of different herbivore assemblages by comparing areas exposed to different numbers of herbivore species or groups of species (or no herbivore species). We will thus include studies that exclude all herbivores, provided that there is a contrast to areas where some herbivore species are present. However, we will exclude studies where the change in the herbivore assemblage is due to changes in the relative abundance of some herbivore species or groups but not in its diversity; for example, we will exclude studies that assess the impacts of changes in abundance of lemmings and voles, unless there is a change in the number of species in the herbivore assemblage (e.g. populations of one species go locally extinct). We will place no restrictions on the type of comparison (e.g. experimental treatments, changes across spatial or temporal gradients).

Eligible outcome (changes in ecological processes, functions and properties of terrestrial Arctic ecosystems)

We will include studies that assess the effects of herbivory on processes, functions and properties of Arctic terrestrial ecosystems, including their effects on other organisms. We will place no restriction on the ecosystem processes, functions and properties reported, but we will exclude studies that only report aspects of the herbivore itself and not its impacts on other ecosystem components (e.g. population dynamics of the focal herbivore, habitat selection studies). We will also exclude studies that do not report an outcome and those that present no analysable primary data for the outcome. We will therefore exclude reviews and book chapters, unless they contain original, primary data. However, we will include, at title and abstract screening, any corrections and erratum to published papers, as these can include relevant information for our review results; such documents will be identified as partly redundant.

Eligible type of study design

We will include primary field studies (observational or manipulative) comparing ecosystem processes, functions and properties in areas and/or time periods with different levels of herbivore diversity. Remote sensing studies will be included but not modelling studies as these do not represent direct effects of herbivores. We will not include greenhouse experiments because they restrict access of natural herbivores, but we will include common garden experiments where herbivores have free access to the experimental areas.

Redundancy

We will exclude studies that report data that are reported in another study. In such cases, we will include the study presenting the longest time series or the greatest number of replicates. For example, if a study presents a spatial or temporal subset of a larger dataset presented in another study, when possible, we will include the larger study. We will assess this first by checking the references cited in the methods section for potential overlap, and by checking studies that were conducted at the same location following Soininen et al. [23].

Study validity assessment

We will assess the validity of studies fulfilling all the eligibility criteria described above. Studies will be categorised as having high, medium or low susceptibility to bias based on seven criteria proposed in the new prototype CEE critical appraisal tool [32]. These criteria include an assessment of the risk of bias due to the presence of confounding variables, post-intervention or exposure biases, misclassified comparison variables (in the case of observational studies) or performance biases (in the case of experimental studies), detection biases, and risks of outcome reporting or assessment biases. Analyses in the systematic review will be conducted with and without studies with high susceptibility of bias to assess their influence on synthesis results. Study validity will be critically appraised by the reviewers conducting the data coding process based on the seven criteria (Additional files 3 and 4). During protocol development we assessed the repeatability of the critical appraisal of study validity as part of the assessment of the data coding and extraction processes (see data coding strategy below) and refined the criteria for study validity assessment as needed. In the systematic review, a random subset of at least 10% of studies will be appraised by a second reviewer (same as for the data extraction process). Information about study validity assessment will be provided for all included studies as an additional file to the systematic review.

Data coding and extraction strategy

For each eligible article, we will extract information on the variables described in Table 1. One article can contain several studies, when separate parts of the article differ in terms of outcome, environmental context or methodological approach. Information will be recorded separately for each study. If an article reports an outcome separately for different habitats, sites, etc. we will consider them as separate studies. If studies report repeated annual measurements, we will extract data from the last measurement. In case of studies reporting repeated measurements over the growing season, we will extract data from the peak of the growing season (i.e. late July-early August in Arctic terrestrial ecosystems). When a study reports several comparisons of herbivore diversity levels, for example in studies using size-selective exclosures, pairwise comparisons will be extracted as separate studies. We will conduct our analyses using species richness of herbivores when studies report sufficient taxonomic detail, but also grouping species into relevant functional groups of herbivores [2].Table 1 List of variables to code from the studies, including key sources of heterogeneity known to influence the effects of herbivores on tundra ecosystems

Topic	Coding variable	Variable description	Source	
Article ID	author_list	List of authors	P	
	title	Title of the publication	P	
	year	Year of publication	P	
	journal	Journal or publishing house	P	
	language	Language of the publication	P	
	study_ID	Unique ID number for study	C	
	searchable_pdf	Whether the document allows for automated text searches	P	
Study location	country	Country or region name	P	
	locality	Site name describing the locality, as specified in the study	P	
	coordinate_source	Whether coordinates are provided in the text, based on maps, figures or place names	P	
	coordinates_N, coordinates_E	Geographic coordinates in decimal degrees	P	
	elevation	Elevation (m a.s.l.)	P	
	elevation_DEM	Elevation (m a.s.l.) extracted from digital elevation model	D	
Study type	study_design	Whether the study involves experimental manipulations, observational approaches or whether the study design is unclear	C	
	diversity_contrast	Approach used to create or assess the difference in herbivore diversity, such as exclosures or spatial contrasts	C	
	size_selective_exclosures	Does the study use different exclosures that prevent access of different sized herbivores?	C	
	spatial_extent	Size of the study area	C	
	spatial_resolution	Spatial scale at which the outcome is reported for each study	C	
	temporal_resolution	Interval between measurements of the outcome	C	
	sampling_frequency	Frequency of measurements	C	
	year_start	Year when the study started	P	
	year_end	Year when the study finished	P	
Exposure and comparator	herbivore_data	How is the composition of the herbivore assemblages assessed?	C	
	herbivore_ID_higher_diversity	Species (or group) list in the areas with higher herbivore diversity	P	
	herbivore_ID_lower_diversity	Species (or group) list in the areas with lower herbivore diversity	P	
Outcome	measured_response_variable	What was measured in areas exposed to different herbivore diversity?	P/C	
	reported_units	Units in which the response variable is reported; "unitless" if there is no unit	P	
	value_higher_diversity	Value of the outcome in the higher diversity area	P	
	value_lower_diversity	Value of the outcome in the lower diversity area	P	
	value_type	Does the value refer to the mean, median, etc	P	
	variability_higher_diversity	Variability of the outcome in the higher diversity area	P	
	variability_lower_diversity	Variability of the outcome in the lower diversity area	P	
	variability_type	Type of variability measurement (SE, confidence intervals (CI), etc.)	P	
	sample_size_higher_diversity	Number of observations in the higher diversity area	P	
	sample_size_lower diversity	Number of observations in the lower diversity area	P	
	effect_size	Value of the comparison between areas with different diversity	P	
	effect_size_variability	Variability of the comparison	P	
	effect_size_variability_type	Type of variability measurement (SE, confidence intervals (CI), etc.)	P	
	effect_size_direction	Direction of the comparison, or unclear	C	
	effect_size_sample_size	Number of observations for the comparison if reported	P	
	effect_size_type	Type of comparison (e.g. model estimate)	C	
	effect_size_comments	Any additional comments	C	
	statistical_test	Type of statistical test reported (e.g. paired t-test)	P	
	statistical_test_value	Value of the test statistic	P	
	statistical_test_df	Degrees of freedom of the statistic	P	
	statistical_test_p	p-value of the statistical test	P	
	outcome_comments	Any additional comments	C	
	more_data	Is there data available in the paper that is not extractable in its current form?	C	
Study validity assessment	bias_risk_criterion1	Biases due to uncontrolled confounding variables that influence both the herbivore diversity levels and the response variable	C	
	bias_risk_criterion2	Biases arising from systematic differences in the selection of areas into the study or analyses after treatment	C	
	bias_risk_criterion3or4	Biases arising from misclassification or measurement of contrasts of herbivore diversity (observational studies) or treatments (experimental studies)	C	
	bias_risk_criterion5	Biases arising from systematic differences in measurements of outcomes	C	
	bias_risk_criterion6	Biases in reporting of study findings	C	
	bias_risk_criterion7	Biases due to error in applied statistical methods	C	
	bias_risk_comments	Please indicate any relevant comments; if there were confounding variables please report them here	C	
Context	distance_to_treeline	Distance (km) to southern border of arctic subzone E	D	
	distance_from_coast	Distance (km) to the coast	D	
	bioclimatic_zone	Bioclimatic zone A to E, or "other"	D	
	temperature	Temperature related variable(s) extracted from WorldClim	D	
	precipitation	Precipitation related variable(s) extracted from Worldclim	D	
	growing_season	Duration of growing season (days)	D	
	productivity	Value of NDVI (vegetation greenness)	D	
	recent_warming	Extent of recent warming	D	
	recent_greening	Extent of recent greening	D	
	extent_of_recent__change	Extent of recent change in growing season length	D	
	soil_chemistry	Soil chemistry	C	
	soil_texture	Broad categories of soil texture	C	
	soil_moisture	Soil moisture as described in the study	C	
	soil_type	Soil type as reported by the authors	P	
	soil_type_D	Soil type	D	
	permafrost	Presence of permafrost	C	
	permafrost_D	Presence of permafrost	D	
	habitat_type	Habitat types using the broad categories defined in CAVM	C	
	habitat_type_D	Habitat types	D	
	management_focus	Whether the study is framed within a management context	C	
	conservation_focus	Whether the study is framed within a conservation context	C	
The coding variables are grouped into broader topics. Source indicates where the data is extracted from: P for publication, D for digital spatial data layers and C classified by the reviewers based on information available in the publication. More detailed descriptions, examples of the variables and references for the digital spatial data (D) are provided in Additional file 3

To be included, studies will have to provide quantitative data on the outcome variables, although this information may take different forms, including estimates of means, measures of variation (standard deviation, standard error, confidence intervals) and sample size for the different levels of herbivore diversity, values for the comparison (effect size) or results of a statistical test for the comparison between different levels of herbivore diversity providing enough information to enable calculation of effect sizes. For example, measures of variability should be reported by the studies or be calculable from the information provided in the study [33]. When possible, we will convert these comparisons to raw values, so that we can compare absolute effects rather than relative changes. If the study provides more than one type of data (for example the study reports both the statistical test for the comparison and the mean and variance for the levels of herbivore diversity being compared), all data will be extracted. For the outcome variables we will extract data from the text, tables and graphs, using image analysis software (ImageJ; [34]) when needed. When raw data are provided, we will calculate summary statistics. Data on potential effect modifiers will also be extracted (see next section). In some cases, it may be useful to ask authors of relevant studies for access to unpublished primary data or to ask for confirmation of missing or unclear information. During the protocol development, we created an Excel sheet with drop-down menus and open fields to ensure consistency during the data coding and data extraction process. This coding template includes additional information per variable (definition, potential values) and examples that help clarify the decision-making process (Additional file 3).

During protocol development we tested our data coding and extraction strategy to evaluate whether it was possible to extract the proposed variables and whether we had excluded potentially relevant variables or categories presented in the studies. Three co-authors (ICB, LBP, EK) coded four studies each from among the studies that had passed the title, abstract and full text screening stages, representing different challenges in the data extraction process. Three of these studies were common to all observers to test for consistency in coding, and one was different to identify additional issues where the data coding strategy had to be refined. While most of the variables coded by all co-authors were consistent, we identified several minor issues and reviewed the coding template accordingly. During the systematic review data from 10% publications will be recorded by two reviewers to assess the repeatability of the data extraction process. The extracted data records will be made available as additional files.

Potential effect modifiers and reasons for heterogeneity

Extraction of meta-data from the studies will include data regarding key sources of heterogeneity (Table 1). The list of potential variables was identified based on expert knowledge and discussions with relevant stakeholders in the protocol development team, and included variables known to influence the effects of herbivores on tundra ecosystems, such as the geographical location of the study (latitude and longitude, distance to treeline), climate, soil type, vegetation or habitat type, proximity to human activities or human management of the herbivores (e.g. hunting, herding) or the ecosystem (e.g. protected areas), the presence and identity of predators and the identity of the herbivores in the herbivore assemblage. These potential modifiers will be used in meta-analyses to account for differences between studies. In most cases, these variables will be extracted from existing data layers or coded in the data extraction process (Additional file 3) to ensure that the information on these variables is extracted in a consistent way across studies (see e.g. [23]).

Data synthesis and presentation

In the systematic review, we will describe the review process and the evidence base, focusing on the different ecosystem processes, functions and properties of terrestrial ecosystems in the Arctic that have been measured in the literature. The outcomes of this systematic review are purposely open-ended and will provide a sense of the scope and volume of the research that has been conducted on herbivore diversity and tundra ecosystems. This will in itself be a valuable output of the systematic review and the overview of outcomes measured by herbivore diversity studies will be summarized using descriptive statistics, tables and figures.

A narrative synthesis of the data will describe the quality of the results and the findings of studies taking into account the study validity assessment and will provide an overview of the evidence. Where enough studies (at least 5) report similar outcomes, meta-analysis will be undertaken [33], including in the models the effect of potential modifiers when possible. We will include studies that report outcomes for contrasting levels of herbivore diversity. Data will thus be often presented as average values and measures of variance for two groups separately. In such cases, we will calculate as an effect size the standardized mean difference (e.g. Hedges’ g), as the difference between group means standardized by the pooled standard deviation. Alternatively, when the outcomes refer to rates of change between the different groups, we may use the response ratio as a measure of effect size [33]. Ultimately, the type of effect size calculated will depend on how the outcomes are reported in these studies. In all cases, effect sizes will be standardized and weighted based on inverse variance. We will take into account the potential for non-independence of effect sizes arising from articles reporting several studies by including article ID as a random factor in our meta-analyses. These mixed-effects models allow assessing whether the heterogeneity in effect sizes arises from the effects of moderator variables or random sources by using inverse variance weights adjusted for the estimated random effects component [33]. Where possible, publication bias and sensitivity analysis will also be conducted using diagnostic plots, like funnel plots of effect size vs study size and the trim-and-fill method [35].

We will present the results using the following illustrations:A flow diagram illustrating the inclusion/exclusion process including the number of papers retained at each stage of the process (ROSES diagram; [36]).

A geographic map showing the locations of the studies included in the systematic review

Scatterplots, plots of means or barplots showing the relationship between herbivore assemblages and different ecosystem processes, functions and properties (mean effect sizes and variance)

Forest plots will summarize meta-analysis outputs, separating by groups where at least five studies are included in each group.

For effective communication of results we will publish the final dataset in full on an interactive map server, where readers can explore and filter the results and visualize results (see for example [23]). Finally, we will include as an additional file to the systematic review a check-list of reporting standards for systematic reviews (i.e. ROSES form), as we do for this systematic review protocol (Additional file 5).

Supplementary Information

Additional file 1.List of local and specialist databases for grey literature searches.

Additional file 2. Test-list of articles used to assess the comprehensiveness of the search.

Additional file 3. Data coding template.

Additional file 4. Criteria for study validity assessment.

Additional file 5. ROSES form for systematic review protocol.

Acknowledgements

This study is a contribution of the Herbivory Network (http://herbivory.lbhi.is), a UArctic Thematic Network.

Authors' contributions

Conceptualization and funding acquisition (ICB, EMS, JDMS, EK, DSH); methodology (ICB, LBP, EMS, JDMS, EK); all authors contributed to writing and revisions of the manuscript (ICB, EMS, JDMS, EK, DSH, LBP, SH). All authors read and approved the final manuscript.

Funding

The project was funded by the Icelandic Research Fund (Grant Nr. 217754) and the European Union’s Horizon 2020 programme (CHARTER project, Grant Agreement Nr. 869471). Funding for open access publication was provided by the Agricultural University of Iceland. The funding bodies had no influence in the design of the study and collection, analysis and interpretation of data.

Availability of data and materials

All data generated or analysed during this study will be included in the published article and its supplementary information files.

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Barrio IC Bueno CG Gartzia M Soininen EM Christie KS Speed JDM Biotic interactions mediate patterns of herbivore diversity in the Arctic Glob Ecol Biogeogr 2016 25 9 1108 1118 10.1111/geb.12470
Barrio IC, Bueno CG, Gartzia M, Soininen EM, Christie KS, Speed JDM, et al. Biotic interactions mediate patterns of herbivore diversity in the Arctic. Glob Ecol Biogeogr. 2016;25(9):1108–18.10.1111/geb.12470
2. Speed JDM Skjelbred IÅ Barrio IC Martin MD Berteaux D Bueno CG Trophic interactions and abiotic factors drive functional and phylogenetic structure of vertebrate herbivore communities across the Arctic tundra biome Ecography 2019 42 1 12 10.1111/ecog.04347
Speed JDM, Skjelbred IÅ, Barrio IC, Martin MD, Berteaux D, Bueno CG, et al. Trophic interactions and abiotic factors drive functional and phylogenetic structure of vertebrate herbivore communities across the Arctic tundra biome. Ecography. 2019;42:1–12. 10.1111/ecog.04347.10.1111/ecog.04347
3. Ritchie M Olff H Olff H Brown V Drent R Herbivore diversity and plant dynamics: compensatory and additive effects Herbivores: between plants and predators 1999 Oxford Blackwell Science 175 204
Ritchie M, Olff H. Herbivore diversity and plant dynamics: compensatory and additive effects. In: Olff H, Brown V, Drent R, editors. Herbivores: between plants and predators. Oxford: Blackwell Science; 1999. p. 175–204.
4. Klein DR Bay C Resource partitioning by mammalian herbivores in the high Arctic Oecologia 1994 97 439 450 10.1007/BF00325880 28313731
Klein DR, Bay C. Resource partitioning by mammalian herbivores in the high Arctic. Oecologia. 1994;97:439–50.28313731 10.1007/BF00325880
5. Tuomi M Väisänen M Ylänne H Brearley FQ Barrio IC Bråthen KA Stomping in silence: conceptualizing trampling effects on soils in polar tundra Funct Ecol 2021 35 306 317 10.1111/1365-2435.13719
Tuomi M, Väisänen M, Ylänne H, Brearley FQ, Barrio IC, Bråthen KA, et al. Stomping in silence: conceptualizing trampling effects on soils in polar tundra. Funct Ecol. 2021;35:306–17.10.1111/1365-2435.13719
6. Olofsson J Hulme PE Oksanen L Suominen O Importance of large and small mammalian herbivores for the plant community structure in the forest tundra ecotone Oikos 2004 106 2 324 334 10.1111/j.0030-1299.2004.13224.x
Olofsson J, Hulme PE, Oksanen L, Suominen O. Importance of large and small mammalian herbivores for the plant community structure in the forest tundra ecotone. Oikos. 2004;106(2):324–34.10.1111/j.0030-1299.2004.13224.x
7. Ravolainen VT Bråthen KA Yoccoz NG Nguyen JK Ims RA Complementary impacts of small rodents and semi-domesticated ungulates limit tall shrub expansion in the tundra J Appl Ecol 2014 51 1 234 241 10.1111/1365-2664.12180
Ravolainen VT, Bråthen KA, Yoccoz NG, Nguyen JK, Ims RA. Complementary impacts of small rodents and semi-domesticated ungulates limit tall shrub expansion in the tundra. J Appl Ecol. 2014;51(1):234–41. 10.1111/1365-2664.12180.10.1111/1365-2664.12180
8. Petit Bon M Gunnarsdotter Inga K Jónsdóttir IS Utsi TA Soininen EM Bråthen KA Interactions between winter and summer herbivory affect spatial and temporal plant nutrient dynamics in tundra grassland communities Oikos 2020 129 8 1229 1242 10.1111/oik.07074
Petit Bon M, Gunnarsdotter Inga K, Jónsdóttir IS, Utsi TA, Soininen EM, Bråthen KA. Interactions between winter and summer herbivory affect spatial and temporal plant nutrient dynamics in tundra grassland communities. Oikos. 2020;129(8):1229–42. 10.1111/oik.07074.10.1111/oik.07074
9. Van Ruijven J De Deyn GB Raaijmakers CE Berendse F Van der Putten WH Interactions between spatially separated herbivores indirectly alter plant diversity Ecol Lett 2004 8 1 30 37 10.1111/j.1461-0248.2004.00688.x
Van Ruijven J, De Deyn GB, Raaijmakers CE, Berendse F, Van der Putten WH. Interactions between spatially separated herbivores indirectly alter plant diversity. Ecol Lett. 2004;8(1):30–7. 10.1111/j.1461-0248.2004.00688.x.10.1111/j.1461-0248.2004.00688.x
10. Burkepile DE Parker JD Recent advances in plant-herbivore interactions F1000Research. 2017 6 119 10.12688/f1000research.10313.1 28232868
Burkepile DE, Parker JD. Recent advances in plant-herbivore interactions. F1000Research. 2017;6:119.28232868 10.12688/f1000research.10313.1
11. Burkepile DE Hay ME Herbivore species richness and feeding complementarity affect community structure and function on a coral reef Proc Natl Acad Sci 2008 105 42 16201 16206 10.1073/pnas.0801946105 18845686
Burkepile DE, Hay ME. Herbivore species richness and feeding complementarity affect community structure and function on a coral reef. Proc Natl Acad Sci. 2008;105(42):16201–6.18845686 10.1073/pnas.0801946105
12. Post E Forchhammer MC Bret-Harte MS Callaghan TV Christensen TR Elberling B Ecological dynamics across the Arctic associated with recent climate change Science 2009 325 5946 1355 1358 10.1126/science.1173113 19745143
Post E, Forchhammer MC, Bret-Harte MS, Callaghan TV, Christensen TR, Elberling B, et al. Ecological dynamics across the Arctic associated with recent climate change. Science. 2009;325(5946):1355–8.19745143 10.1126/science.1173113
13. Juhasz CC Shipley B Gauthier G Berteaux D Lecomte N Direct and indirect effects of regional and local climatic factors on trophic interactions in the Arctic tundra J Anim Ecol 2020 89 3 704 715 10.1111/1365-2656.13104 31538330
Juhasz CC, Shipley B, Gauthier G, Berteaux D, Lecomte N. Direct and indirect effects of regional and local climatic factors on trophic interactions in the Arctic tundra. J Anim Ecol. 2020;89(3):704–15.31538330 10.1111/1365-2656.13104
14. Speed JDM Chimal JAB Martin MD Barrio IC Vuorinen KEM Soininen EM Will borealization of Arctic tundra herbivore communities be driven by climate warming or vegetation change? Glob Chang Biol 2021 27 24 6568 6577 10.1111/gcb.15910 34592044
Speed JDM, Chimal JAB, Martin MD, Barrio IC, Vuorinen KEM, Soininen EM. Will borealization of Arctic tundra herbivore communities be driven by climate warming or vegetation change? Glob Chang Biol. 2021;27(24):6568–77.34592044 10.1111/gcb.15910
15. Barnosky AD Lindsey EL Villavicencio NA Bostelmann E Hadly EA Variable impact of late-Quaternary megafaunal extinction in causing ecological state shifts in North and South America Proc Natl Acad Sci 2016 113 4 856 861 10.1073/pnas.1505295112 26504219
Barnosky AD, Lindsey EL, Villavicencio NA, Bostelmann E, Hadly EA. Variable impact of late-Quaternary megafaunal extinction in causing ecological state shifts in North and South America. Proc Natl Acad Sci. 2016;113(4):856–61.26504219 10.1073/pnas.1505295112
16. Zimov SA Chuprynin VI Oreshko AP Chapin FS III Reynolds JF Chapin M Steppe-tundra transition: a herbivore-driven biome shift at the end of the Pleistocene Am Nat 1995 146 5 765 794 10.1086/285824
Zimov SA, Chuprynin VI, Oreshko AP, Chapin FS III, Reynolds JF, Chapin M. Steppe-tundra transition: a herbivore-driven biome shift at the end of the Pleistocene. Am Nat. 1995;146(5):765–94.10.1086/285824
17. Olofsson J Post E Effects of large herbivores on tundra vegetation in a changing climate, and implications for rewilding Philos Trans R Soc London 2018 373 20170437 10.1098/rstb.2017.0437 30348880
Olofsson J, Post E. Effects of large herbivores on tundra vegetation in a changing climate, and implications for rewilding. Philos Trans R Soc London. 2018;373:20170437.30348880 10.1098/rstb.2017.0437
18. Egelkraut D Kardol P De Long JR Olofsson J The role of plant–soil feedbacks in stabilizing a reindeer-induced vegetation shift in subarctic tundra Funct Ecol 2018 32 8 1959 1971 10.1111/1365-2435.13113
Egelkraut D, Kardol P, De Long JR, Olofsson J. The role of plant–soil feedbacks in stabilizing a reindeer-induced vegetation shift in subarctic tundra. Funct Ecol. 2018;32(8):1959–71.10.1111/1365-2435.13113
19. Vindstad OPL Jepsen JU Ek M Pepi A Ims RA Can novel pest outbreaks drive ecosystem transitions in northern-boreal birch forest? J Ecol 2019 107 3 1141 1153 10.1111/1365-2745.13093
Vindstad OPL, Jepsen JU, Ek M, Pepi A, Ims RA. Can novel pest outbreaks drive ecosystem transitions in northern-boreal birch forest? J Ecol. 2019;107(3):1141–53.10.1111/1365-2745.13093
20. Speed JDM Austrheim G Kolstad AL Solberg EJ Long-term changes in northern large-herbivore communities reveal differential rewilding rates in space and time PLoS ONE 2019 14 5 e0217166 10.1371/journal.pone.0217166 31112579
Speed JDM, Austrheim G, Kolstad AL, Solberg EJ. Long-term changes in northern large-herbivore communities reveal differential rewilding rates in space and time. PLoS ONE. 2019;14(5):e0217166.31112579 10.1371/journal.pone.0217166
21. Duffy JE Cardinale BJ France KE McIntyre PB Thébault E Loreau M The functional role of biodiversity in ecosystems: incorporating trophic complexity Ecol Lett 2007 10 6 522 538 10.1111/j.1461-0248.2007.01037.x 17498151
Duffy JE, Cardinale BJ, France KE, McIntyre PB, Thébault E, Loreau M. The functional role of biodiversity in ecosystems: incorporating trophic complexity. Ecol Lett. 2007;10(6):522–38.17498151 10.1111/j.1461-0248.2007.01037.x
22. Wang L Delgado-Baquerizo M Wang D Isbell F Liu J Feng C Diversifying livestock promotes multidiversity and multifunctionality in managed grasslands Proc Natl Acad Sci 2019 116 13 201807354 10.1073/pnas.1807354116
Wang L, Delgado-Baquerizo M, Wang D, Isbell F, Liu J, Feng C, et al. Diversifying livestock promotes multidiversity and multifunctionality in managed grasslands. Proc Natl Acad Sci. 2019;116(13):201807354.10.1073/pnas.1807354116
23. Soininen EM Barrio IC Bjørkås R Björnsdóttir K Ehrich D Hopping K Location of studies and evidence of effects of herbivory on Arctic vegetation: a systematic map Environ Evid 2021 10 25 10.1186/s13750-021-00240-0
Soininen EM, Barrio IC, Bjørkås R, Björnsdóttir K, Ehrich D, Hopping K, et al. Location of studies and evidence of effects of herbivory on Arctic vegetation: a systematic map. Environ Evid. 2021;10:25.10.1186/s13750-021-00240-0
24. Soininen EM Barrio I Jepsen JU Ehrich D Ravolainen VT Speed JDM Evidence of effects of herbivory on Arctic vegetation: a systematic map protocol Environ Evid 2018 7 1 23 10.1186/s13750-018-0135-1
Soininen EM, Barrio I, Jepsen JU, Ehrich D, Ravolainen VT, Speed JDM. Evidence of effects of herbivory on Arctic vegetation: a systematic map protocol. Environ Evid. 2018;7(1):23. 10.1186/s13750-018-0135-1.10.1186/s13750-018-0135-1
25. Haddaway NR Collins AM Coughlin D Kirk S The role of Google Scholar in evidence reviews and its applicability to grey literature searching PLoS ONE 2015 10 9 e0138237 10.1371/journal.pone.0138237 26379270
Haddaway NR, Collins AM, Coughlin D, Kirk S. The role of Google Scholar in evidence reviews and its applicability to grey literature searching. PLoS ONE. 2015;10(9):e0138237.26379270 10.1371/journal.pone.0138237
26. Livoreil B Glanville J Haddaway NR Bayliss H Bethel A De Lachapelle FF Systematic searching for environmental evidence using multiple tools and sources Environ Evid 2017 6 1 23 10.1186/s13750-017-0099-6
Livoreil B, Glanville J, Haddaway NR, Bayliss H, Bethel A, De Lachapelle FF, et al. Systematic searching for environmental evidence using multiple tools and sources. Environ Evid. 2017;6(1):23.10.1186/s13750-017-0099-6
27. Foo YZ O’Dea RE Koricheva J Nakagawa S Lagisz M A practical guide to question formation, systematic searching and study screening for literature reviews in ecology and evolution Methods Ecol Evol. 2021 12 9 1705 1720 10.1111/2041-210X.13654
Foo YZ, O’Dea RE, Koricheva J, Nakagawa S, Lagisz M. A practical guide to question formation, systematic searching and study screening for literature reviews in ecology and evolution. Methods Ecol Evol. 2021;12(9):1705–20.10.1111/2041-210X.13654
28. Cohen J A coefficient of agreement for nominal scales Educ Psychol Meas 1960 20 1 37 46 10.1177/001316446002000104
Cohen J. A coefficient of agreement for nominal scales. Educ Psychol Meas. 1960;20(1):37–46.10.1177/001316446002000104
29. CAFF [Conservation of Arctic Flora and Fauna]. Boundary for Conservation of Arctic Flora and Fauna (CAFF) working group of the Arctic Council. Arctic Council; 2017.
30. CAFF [Conservation of Arctic Flora and Fauna]. Arctic Biodiversity Assessment. Status and trends in Arctic biodiversity. Akureyri; 2013. http://www.arcticbiodiversity.is/
31. Jia S Wang X Yuan Z Lin F Ye J Hao Z Global signal of top-down control of terrestrial plant communities by herbivores Proc Natl Acad Sci U S A 2018 115 24 6237 6242 10.1073/pnas.1707984115 29848630
Jia S, Wang X, Yuan Z, Lin F, Ye J, Hao Z, et al. Global signal of top-down control of terrestrial plant communities by herbivores. Proc Natl Acad Sci U S A. 2018;115(24):6237–42.29848630 10.1073/pnas.1707984115
32. Konno K, Livoreil B, Pullin A. Collaboration for environmental evidence critical appraisal tool version 0.3 (prototype). 2021. https://environmentalevidence.org/cee-critical-appraisal-tool/
33. Harrison F Getting started with meta-analysis Methods Ecol Evol 2011 2 1 1 10 10.1111/j.2041-210X.2010.00056.x
Harrison F. Getting started with meta-analysis. Methods Ecol Evol. 2011;2(1):1–10. 10.1111/j.2041-210X.2010.00056.x.10.1111/j.2041-210X.2010.00056.x
34. Schneider CA Rasband WS Eliceiri KW NIH Image to ImageJ: 25 years of image analysis Nat Methods 2012 9 7 671 675 10.1038/nmeth.2089 22930834
Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9(7):671–5.22930834 10.1038/nmeth.2089
35. Duval S Tweedie R Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis Biometrics 2000 56 455 463 10.1111/j.0006-341X.2000.00455.x 10877304
Duval S, Tweedie R. Trim and fill: a simple funnel-plot-based method of testing and adjusting for publication bias in meta-analysis. Biometrics. 2000;56:455–63.10877304 10.1111/j.0006-341X.2000.00455.x
36. Haddaway NR Macura B Whaley P Pullin AS ROSES Reporting standards for Systematic Evidence Syntheses: pro forma, flow-diagram and descriptive summary of the plan and conduct of environmental systematic reviews and systematic maps Environ Evid 2018 7 1 4 11 10.1186/s13750-018-0121-7
Haddaway NR, Macura B, Whaley P, Pullin AS. ROSES Reporting standards for Systematic Evidence Syntheses: pro forma, flow-diagram and descriptive summary of the plan and conduct of environmental systematic reviews and systematic maps. Environ Evid. 2018;7(1):4–11.10.1186/s13750-018-0121-7
