Alpine areas of the tropical Andes constitute the largest of all tropical alpine regions worldwide. They experience a particularly harsh climate, and they are fragmented into tropical alpine islands at various spatial scales. These factors generate unique patterns of continental insularity, whose impacts on biodiversity remain to be examined precisely. By reviewing existing literature and by presenting unpublished data on beta-diversity and endemism for a wide array of taxonomic groups, we aimed at providing a clear, overall picture of the isolation-biodiversity relationship in the tropical alpine environments of the Andes. Our analyses showed that (1) taxa with better dispersal capacities and wider distributions (e.g., grasses and birds) were less restricted to alpine areas at local scale; (2) similarity among communities decreased with spatial distance between isolated alpine areas; and (3) endemism reached a peak in small alpine areas strongly isolated from main alpine islands. These results pinpoint continental insularity as a powerful driver of biodiversity in the tropical High Andes. A combination of human activities and warming is expected to increase the effects of continental insularity in the next decades, especially by amplifying the resistance of the lowland matrix that surrounds tropical alpine islands.
Tropical alpine environments (TAEs) represent as much as 10% of alpine systems worldwide and display unique environmental characteristics among alpine regions (Rundel et al., 1994; Körner, 2003). Among them, low seasonality, inverted rainfall gradients above 3000–3500 m, and the high frequency of freeze-thaw cycles (see Anthelme and Dangles, 2012) set up a peculiar pool of environmental conditions that have shaped the evolution of species assemblages substantially different from those that live in extratropical alpine systems (Hedberg and Hedberg, 1979; Ramsay and Oxley, 1997; Navas, 2003; Jacobsen et al., 2010). For instance, plant growth forms like giant rosettes, giant cushions (Fig. 1, parts A and C), and tussock grasses are more abundant under the tropics than anywhere else in alpine regions (Hedberg and Hedberg, 1979; Körner, 2003). Another specificity of TAE is that they develop at higher elevation than any other alpine ecosystem (Körner, 2003), maybe with the exception of dry subtropical regions (Halloy, 1991; Macek et al., 2012). This argument has been commonly cited to generate a higher isolation of alpine systems under the tropics, with fragmented TAE constituting an archipelago of continental islands (Smith and Young, 1987; Luteyn, 1999; Hughes and Eastwood, 2006; Sklenář et al., 2014).
Yet, the effects of such pronounced continental insularity on biodiversity patterns and dynamics are not straightforward. Some authors argue that continental insularity is one of the main drivers of the high proportion of endemism observed in TAE (Luteyn, 1999), especially because topographical barriers may be more effective than in extratropical regions (Janzen, 1967; Richter et al., 2009). Other authors tend to reduce its influence because the lowlands located between TAE are not only a matrix resistant to the dispersion of alpine taxa, but they may act instead as a source of colonization during climatic fluctuations (e.g., Smith and Young, 1987). The extent to which continental insularity explains patterns of diversity in tropical alpine systems in comparison with oceanic insularity is a topical question that has been poorly addressed so far and certainly deserves more attention given the high levels of biodiversity found in TAE and their high sensitivity to current global changes (Luteyn, 1999; Bradley et al., 2006).
The tropical High Andes represent an ideal study region to examine the effects of continental insularity on biodiversity patterns for several reasons. First, they are particularly representative of TAEs as they gather more than 90% of their area worldwide (Jacobsen, 2008a). Second, they shelter a wide latitudinal range of environmental conditions from the humid páramos in the north to the dry puna in the south. Third, they display both highly fragmented alpine areas, as in Colombia and Venezuela (Vuilleumier, 1970), and continuous alpine areas in the central Andes, as found in the Peruvian and Bolivian altiplano (Fig. 2). Also, by encompassing Colombia, Venezuela, Ecuador, Peru, Bolivia, North Chile, and North Argentina, tropical Andes are the largest tropical mountain range worldwide, providing opportunities for tropical alpine biodiversity to diverge not only because of habitat fragmentation, but also because of the effects of geographical distance along a continuous topographical unit (Hughes and Eastwood, 2006; Sklenář et al., 2014).
Continental insularity in the tropical Andes not only is a matter of single isolated summits but also is effective at a coarser, regional scale, with entire mountain ranges being separated from the main cordillera (e.g., Cordillera del Condor in Ecuador, Cordillera Apolobamba in Bolivia, and Sierra Nevada de Santa Marta in Colombia). The most important interruption of Andean tropical alpine regions is the Huancabamba deflection in northern Peru (Fig. 2; Lauer, 1968; Simpson, 1975). There, both the Eastern and the Central Cordillera are entirely interrupted and the only mountain chain providing continuous elevation above 2150 m a.s.l. is the very dry mountain range of the Western Cordillera, which may have contributed to increasing contemporary isolation between the southern and the northern tropical Andes, thus being likely a driver of dissimilarities in floristic elements between these two regions (Smith and Cleef, 1988; Sklenář et al., 2011). However, current patterns of community dissimilarity between Central and Northern Andes are also indebted to more ancient biogeographical processes, including the effects of a former lowland corridor invaded by wetlands, termed the “Western Andean Portal” at this location (Antonelli and Sanmartín, 2011).
The tropical High Andes display three main features—harsh environmental conditions, large spatial extent, and strong habitat fragmentation—that are expected to greatly influence the organization of biodiversity at both local and regional scales. By analyzing variations in beta-diversity of different plant and animal groups at various spatial scales, the objective of our study is to describe general patterns in the organization and endemism of natural communities in the alpine areas of the tropical High Andes. Beta-diversity among sampling sites is indeed an important criterion for obtaining adequate representation of regional biodiversity (Kattan et al., 2006), an issue that has been poorly considered at the scale of the tropical Andes (Kessler et al., 2011). We specifically intended to answer the following questions: (1) What proportions of the taxa that live in tropical alpine highlands are strictly restricted to this environment? (2) How does spatial distance between pairs of highlands affect community similarity? (3) How does the spatial configuration of tropical highlands influence the number of endemic species? It is not the objective of this paper to discuss the biogeographical/historical reasons for the present patterns of biodiversity found in tropical alpine islands of the Andes. Rather, our “ecological” approach, by providing an original insight from a different angle, may be used by biogeographers to provide a synthetic overview on these patterns.
The study area was limited to the tropical Andes sensu stricto, from Colombia (Sierra Nevada de Santa Marta; 11°N) to the Southern tropics (23°S; Fig. 1), thus excluding the tropical alpine regions of Central America. Most of the literature considers that the lower limit of TAE in the Andes occur between 3000 m and 3500 m a.s.l. (Vuilleumier, 1970; Smith and Cleef, 1988; Rundel et al., 1994; Luteyn, 1999; Körner, 2003). We arbitrary used the elevation 3200 m as the lower limit of TAE, and focused on open environments in our review. However, given that Andean forest in some places develop at higher elevation, we eliminated data collected in mountain forests above 3200 m. As an exception, distribution data for terrestrial arthropods were only available from 3400 m upward (Moret, 2005, unpublished data), which sets the lower limit of alpine areas for this group.
Our database was extracted from the literature and unpublished data, and each study contributed one or several analyses (Table 1). Data in the literature had to fulfil two criteria: (1) to be located in tropical alpine regions of the Andes, and (2) to provide clear indices of beta-diversity, or at least raw data, to be considered in our meta-analyses. Among plants, information was provided for spermatophytes (wetland, terrestrial, and the specific group of grasses), pteridophytes, and bryophytes. Among animals, information was provided for aquatic (trichopteran) and terrestrial arthropods (carabids), frogs, reptiles, birds, and mammals. Most data were provided at species level, but some relevant data at genera level were also included in the analyses (Table 1).
We then performed three types of analyses. First, at local/country scale we calculated the proportion of taxa restricted to (or precinctive of) alpine areas as the total number of taxa found only above the lower elevational limit of the alpine area divided by the total number of taxa found in the alpine area. Depending on available data, this calculation was applied along specific elevation transects or along elevation transects at country scale (see Table 1).
Second, at regional scale we plotted the community similarity vs. geographical distance relationship among pairs of alpine areas separated by a non-alpine lowland matrix for six taxa in Bolivia (3) and in Ecuador (3; see number of sites and matrices in Appendix Table A1). Community similarity was quantified using the Sørensen's index of similarity S between two alpine communities as follows:http://rsb.info.nih.gov/ij). When the alpine area was clearly dominated by a single summit, we chose this summit as the starting/ending point for the spatial distance. However, when no single summit was easily identified (e.g., Nudo de Cajas, Ecuador) or when several summits of similar importance coexisted in the same area (e.g., Antisana and Cotopaxi, Ecuador), then we selected the center or the equidistance between the summits of the alpine area to perform our measurement.
Third, we compiled data on “strict” endemic species of several alpine areas for three studies that provided particularly clean data (plants and birds in Colombia; Arthropods in Ecuador; 14, 15, and 15 sites, respectively, see Appendix Table A2). Strict endemism occurred as soon as one species was restricted to only one alpine area in the whole study area. Strict endemism was put in relation with the spatial extent of each alpine area (ranging from 0.73 km2 to 3487 km2; see Appendix Table A2) and with the distance to nearest large alpine area (minimum 200 km2 for plants and birds; 100 km2 for carabid beetles), two of the best predictors of strict endemism as described in the literature (Vuilleumier, 1970; Simpson, 1975). Note that for beetles (Moret, 2005, unpublished data) the data set comes from Ecuador where alpine areas are much less fragmented than in Colombia (data on birds and plants; Vuilleumier, 1970; Simpson, 1975). Therefore, to make taxonomic groups comparable, we lifted the lower limit of alpine areas of arthropods to 4000 m, which resulted in a similar pattern of habitat fragmentation.
The different nomenclatures used in each of the five countries sampled may generate lists of species composition that are not fully comparable (e.g., the genus Deyeuxia [Poaceae] in Bolivia is described as Calamagrostis in Peru). However, since our analyses did not compare species composition between countries, but rather within each country or each alpine site, this divergence did not affect our data.
The similarity—geographic distance relationship was fitted with linear regressions and mean values of similarity were compared across two-sample T-tests. The proportion of animal vs. plant precinctive taxa was also inferred across two-sample T-tests. The relation between strict endemism and the two variables “area” and “distance to nearest large alpine area” was tested for each taxonomic group with stepwise regression (alpha to enter = 0.15; alpha to exit = 0.15; R2 adjusted when two explaining variables). Results are presented in the form of mean ± standard error. Extrapolation provided by three-dimensional surface plot for strict endemism used the distance method (distance power: 2). All analyses were performed using MINITAB 15.
BIODIVERSITY SPECIFIC TO ALPINE AREAS (FIG. 3)
The overall mean proportion of taxa precinctive to alpine areas was 0.35 ± 0.03, ranging from 0.09 (mammals and pteridophytes in Bolivia) to 0.71 (frogs in Bolivia). No significant differences were observed between animal and plant taxa (two-samples T-tests: p > 0.05). Among animals, the most precinctive taxa were amphibians/reptiles (0.71 in Bolivia and 0.67 in Peru) and terrestrial arthropods (0.62 in Ecuador). In comparison, mammals and birds displayed lower precinctive scores, from 0.09 to 0.33. Among plants, the six most precinctive scores concerned spermatophytes, with proportions ranging from 0.50 to 0.56. The lowest value of precinctive taxa among spermatophytes was obtained for Poaceae (0.26). Pteridophytes were less restricted to alpine areas, but their proportion of precinctive taxa were highly variable (average value: 0.26 ± 0.05, ranging from 0.09 to 0.45). Bryophytes were overall the less precinctive taxon and had an average value of 0.22 ± 0.01 (three studies).
Local studies provided a larger proportion of precinctive taxa than studies at the scale of a whole country. This was observed (1) for amphibians where the three most precinctive studies were supported by local data, (2) for birds (local study: 0.33, study at country level: 0.12), and (3) for pteridophytes, with the only study at local scale providing higher precinctive proportion than any other study at country scale (0.45 versus an average value of 0.22). We observed no clear effect of latitude or country on precinctive scores (one-way ANOVAs on both variables: p > 0.05).
COMMUNITY SIMILARITY BETWEEN ALPINE AREAS: EFFECTS OF DISTANCE (FIG. 4)
The overall similarity between communities from disjointed alpine areas was highest for amphibians, terrestrial plants, and reptiles (mean value of Sørensen's indices: 0.41 ± 0.07, 0.35 ± 0.03, and 0.34 ± 0.06, respectively); intermediate in wetland plants (0.22 ± 0.03); and lowest in mammals and arthropods (0.14 ± 0.04 and 0.16 ± 0.02, respectively). Similarity between communities reduced significantly with spatial distance between alpine areas for each taxonomic group (R2 ranging from 0.18 to 0.53; p ≤ 0.05), except for amphibians (R2 = 0.05; p > 0.05). Among plants, average similarity between terrestrial plants exceeded similarity between wetland plants (0.35 ± 0.03 vs. 0.22 ± 0.03), but pairwise distance correlated better with similarity for terrestrial plants than for wetland plants (R2 = 0.53 and R2 = 0.28, respectively). Among animals, the best relationship between distance and community similarity was obtained with reptiles and carabids (R2 = 0.47 and R2 = 0.38, respectively; p ≤ 0.001), whereas distance explained 18% of variation in the similarity observed in mammal communities (p ≤ 0.05).
DRIVERS OF STRICT ENDEMISM IN ISOLATED ALPINE AREAS (FIG. 5)
The average proportion of strict endemics in selected alpine areas was 0.32 ± 0.09 for carabid beetles (terrestrial arthropods) in Ecuador, 0.21 ± 0.05 for plants in Colombia, and 0.13 ± 0.05 for birds in Colombia (see Appendix Table A2). In each group, the maximum proportion of strict endemism was observed (1) in the smallest alpine areas, and (2) at the most remote sites from nearest large alpine areas. Both variables explained 33% of the strict endemism in plants (stepwise regression; R2 adjusted = 0.33, p ≤ 0.05). However, only the “distance to nearest large alpine area” term, and not “area,” had a significant effect on the observed number of strict endemics in carabids and birds (stepwise regressions; arthropods: R2 = 0.61, p ≤⃒ 0.001; birds: R2 = 0.47, p ≤ 0.01). In each group, a secondary peak of strict endemism was observed in small alpine areas with short distance to the nearest large alpine area. In plants and birds, this peak referred to data from the Sierra Nevada de Santa Marta (Colombia), while in carabids it was sustained by data from the páramo ‘Llanganates’; (Ecuador).
CONTINENTAL INSULARITY AT VARIOUS SCALES EXPLAINS BIODIVERSITY PATTERNS
In contrast with patterns of plant species richness, which are poorly related to current degree of isolation and extent of tropical alpine areas of the Andes (see review in Smith and Young, 1987), our data demonstrate that beta-diversity and endemism in the tropical alpine islands of the Andes are certainly dependent on the current degree of continental insularity, at various spatial scales.
At local/country scale, elevation shifts from the Andean forests to alpine areas (3000–3500 m) resulted in the existence of alpine communities with more than one-third of precinctive taxa (Fig. 3), thus revealing a strong environmental barrier between these two types of ecosystems. Contrastingly to other alpine environments, TAEs experience unusual abiotic constraints such as high frequency of freeze-thaw cycles, edaphic and atmospheric aridity overall stronger than in other alpine regions due to the absence of protective snow cover, and/or the inversion of precipitation gradients (see review in Anthelme and Dangles, 2012). In some taxonomic groups (e.g., spermatophytes), these constraints are easily noticeable in the development of peculiar growth forms, such as giant rosettes of giant cushion-like species (Fig. 1, parts A and C), which are well adapted to these specific constraints (Monasterio and Sarmiento, 1991; Kleier and Rundel, 2009; Anthelme et al., 2012), and by physiological traits particularly resistant to low temperature and water stress (e.g., Rundel et al., 2003; Macek et al., 2009; Sklenář et al., 2010; Almeida et al., 2013). Lower precinctive scores in pteridophytes and bryophytes may be indebted to the higher dispersal capacities of these taxa combined with an overall high tolerance and/or resistance to desiccation (Sakai and Larcher, 1987; Muñoz et al., 2004; Anthelme et al., 2011), which may permit species to occupy wider vertical distribution ranges. However, higher endemism in bryophytes observed in the northern Andes than in the central Andes (Churchill, 2009) suggests that continental insularity also affects the diversity of this group. Concerning animals, while taxon richness generally decreases with elevation (Jacobsen, 2008b; Herzog et al., 2011, and reference therein), the highland fauna, however, is not simply an attenuated version of the adjacent lowland fauna and many taxa are specific to high-elevation sites (e.g., the frogs of the genera Atelopus [Navas, 2003], the rodent Thomasomys vulcani [Tirira, 2007], the bird Theristicus melanopis [Ridgely and Greenfield, 2006]; Fig. 1, part B). Moreover, the lower proportion of precinctive taxa in birds and mammals in comparison with amphibians, reptiles, and arthropods is likely due to group-specific differences in (1) mobility (Brown and Lomolino, 1998) and (2) a tolerance to greater climatic constraints at higher elevation, both in favor of endotherm mammals and amphibians (Laurance et al., 2011). From this viewpoint, carabids considered in our study are not fully representative of terrestrial arthropods, which may be expected to display overall lower proportions of precinctive taxa (Gonzalez and Engel, 2004). However, morphological adaptations such as wing atrophy or smaller body size are known from a number of tropical alpine insects, increasing their possibility to adapt to local conditions by finding sheltered microhabitats, but limiting their dispersion to adjacent lowlands (Sømme, 1989). Specialization to high elevation habitats can also occur in highly mobile organisms such as birds, for which ca. 15–20 species are limited to Polylepis woodlands above 3500 m a.s.l. in the Peruvian and Bolivian puna (e.g., Oreomanes fraseri, which feed exclusively on branches and trunks of Polylepis; Fjeldsa, 1993). However, from a methodological viewpoint, we must keep in mind that variation in the area sampled has a substantial influence on the proportion of precinctive taxa. The reduction of this proportion at country scale might be the result of a “compensation effect” with taxa precinctive on a single elevation gradient becoming ubiquitous at country scale.
At regional scale, isolation between pairs of islands has frequently been cited to increase beta-diversity in the tropical High Andes (Smith and Young, 1987; Luteyn, 1999; Hughes and Eastwood, 2006). Our data corroborate this opinion for a large pool of taxonomic groups, with distance among tropical alpine islands being negatively correlated with the similarity between alpine communities (except for amphibians). In plants, terrestrial communities are without surprise more similar than communities restricted to humid ecosystems (wetland plants), which are much highly scattered in the High Andes. But the reduced influence of distance on the community similarity of the high Andean wetlands is an original result, which may deserve further investigation (existence of a relatively homogeneous wetland plant community at the scale of the tropical Andes?). Interestingly, when considering only the Poaceae family among terrestrial plants (which is one of the families with the most efficient propagule dispersion capacities across wind dispersal), similarity is much higher than in other studies on spermatophytes. Similarly, wind pollination capacities have been shown to reduce the effects of isolation by distance in the High Andes (Schmidt-Lebuhn et al., 2007). Accordingly, seed dispersal and pollination capacity (by wind) are certainly major drivers of plant beta-diversity in isolated alpine islands of the Andes, as initially suggested by Hedberg (1969). Among others, it explains why one of the most isolated alpine areas in the tropical Andes, the Sumaco Volcano in Ecuadorian Amazonia, is largely dominated by Poaceae and more generally speaking by plants with high seed dispersal capacities (Løjtnant and Molau, 1983).
Among animals, the overall pattern of similarity between pairs of communities in the Bolivian puna may appear counterintuitive as (non-volant) mammals display a stronger similarity among alpine islands than amphibians and reptiles, which generally have lower dispersal capacities. As proposed by Tarifa et al. (2007), such pattern can be explained by (1) a reduced species pool of amphibians and reptiles in the dry puna, (2) the occurrence of a few cosmopolitan species, and (3) a better adaptation of mammals to specific niches in alpine areas. However, distance-similarity relationships may differ in more humid, high-elevation locations of the tropical Andes, such as in the Colombian and Ecuadorian páramos, which are known to shelter site-specific amphibian communities (Navas, 2003; Ron et al., 2012). Also, species of the frog genera Telmatobius (see Fig. 1, part D) are known to be highly endemic from specific locations, with species turnover taking place within a few tens or, at most, hundreds of kilometers (De La Riva et al., 2005). Overall, our knowledge of the factors influencing community similarity among alpine areas is likely to be influenced by the size of the species pool considered and also strongly limited by the patchiness of available information in the tropical Andes.
In addition to variation in elevation at local scale and distance between pairs of alpine islands, two other variables—“distance to nearest large alpine area” and to a lesser extent “alpine area”—were positively correlated with strict endemism for the three studied groups (arthropods, plants, birds). This finding supports our argumentation that the current degree of isolation of continental insularity in the tropical Andes is a strong driver of the organization of beta-diversity (Vuilleumier, 1970; Simpson, 1975).
CONTINENTAL INSULARITY DEPENDS ON THE RESISTANCE OF THE MATRIX
The existence of secondary peaks of strict endemism in relatively small areas such as the Sierra Nevada de Santa Marta, which is separated from the main Cordillera by a very low lowland matrix (Fig. 5, parts A and B), suggests that lowland matrices likely have a key influence on beta-diversity patterns. This is supported by various observations on TAE plant communities in East Africa (Hedberg, 1969), Ecuador (Løjtnant and Molau, 1983), and Colombia (Sklenář and Balslev, 2005), and on arthropod communities in Ecuador (Moret, 2005). In contrast to oceanic islands where the quality of the matrix is homogenous, continental insularity implies different degrees of matrix resistances among isolated alpine areas. As shown above, the resistance may primarily be dependent on elevation, with lowlands being barriers, but also opportunities for colonization, especially under the effects of large-scale climatic fluctuations (Smith and Young, 1987; Sklenář et al., 2011). Nevertheless, other factors may influence the resistance of the continental matrix, among which are human activities such as the rapidly advancing agricultural frontier at higher elevations (Hedberg, 1969; Sklenář and Ramsay, 2001; Morales and Sarmiento, 2002). Furthermore, this resistance differs among groups, as for plants it is usually lethal landing outside of TAE during dispersal event, while for animals it does not have to be.
MULTIPLE DRIVERS OF CONTINENTAL INSULARITY
The purpose of this study was to examine the effects of contemporary isolation by distance, area, and variation in elevation on the organization of biodiversity in the tropical Andes. However, apart from biogeographical and historical drivers, whose role on biodiversity organization in the High Andes has been amply documented (e.g., Fjeldsa, 1993; Richter et al., 2009; Sklenář et al., 2011), other factors would certainly deserve more attention when trying to provide a synthetic explanation of observed diversity patterns in the tropical Andes. For example, local climatic variations, especially those observed between western and eastern sides of large Andean summits, are potential drivers of isolation (e.g., see dry and isolated “Arenal Grande” on the western side of Chimborazo, Ecuador; Moret, 2009; Sklenář and Ramsay, 2001; Kessler et al., 2011). Also, the presence of large lakes (Lake Titicaca; ancient lake Minchin) was cited as a potential factor of isolation in the altiplano (Simpson, 1975). As said before, human activities are expected to modify the degree of isolation of alpine areas by making the surrounding matrix more resistant to migration of natural species, but also by promoting the transport and introduction of species adapted to such disturbances (e.g., invasive species; Pauchard et al., 2009).
From this viewpoint, the secondary peaks of strict endemism that have been observed in several small alpine areas, despite their short distance to the nearest large alpine area (Fig. 5), deserve a special attention. For example, the Llanganates area in central Ecuador displays an outstanding 50% of strict endemic Carabid beetles above 4000 m (Moret, unpublished data), which cannot be explained by the existence of a topographical barrier. In this case, future studies should consider a combination of other factors, such as wet microclimate, absence of anthropogenic perturbation, specific soil properties, and most of all absence of Pleistocene and Holocene volcanic activity.
Finally, because of their high elevation, tropical alpine regions of the Andes are probably one of the terrestrial ecosystems that will face the highest warming up to 2100 (Bradley et al., 2006). In this context, species that occur in alpine areas may be especially prone to extinction under the effects of climate change because (1) they inhabit environments with a relatively low total colonizable area and have nowhere to migrate upward (Gosling et al., 2009), and (2) the rapid shift of the upper limits of alpine areas to higher elevation increases the isolation of contiguous viable populations (Jørgensen et al., 2011; Larsen et al., 2011; Laurance et al., 2011; Velásquez-Tibatá et al., 2013). It must be taken into account, however, that species at higher elevations may be able to reduce their extinction risk by having possibly larger vertical distribution ranges (e.g., Herzog et al., 2013). Further research in the line of that presented in this paper is urgently needed to better predict the response of the biodiversity of high mountains to rapid anthropogenic changes in the tropical Andes.
We warmly thank Carlos Boada, Michael Kessler, Santiago Ron, the Herbario Nacional de Bolivia (LPB), and the Missouri Botanical Garden for providing us with unpublished information that consistently improved our database. We are also grateful to James Juvik for inviting Fabien Anthelme to the international conference “Vulnerable Islands in the Sky: Science and Management of Tropical Island Alpine & Sub-alpine Ecosystems” held in Hawai‘i in August 2012. Macek was supported by MSMT LM2010009 CzechPolar.