Is Rich and Rare the Common Share? Describing Biodiversity Patterns to Inform Conservation Practices for South American Anurans

Species richness and range size are key features of biogeographic and macroecological analyses, which can yield a first assessment tool to define conservation priorities. Here we combined both features in a simultaneous analysis, based on range-diversity plots, to identify sets of rich-rare (high species richness with restricted ranges) and poor-rare cells (low species richness with restricted ranges). We applied this analysis to the anurans of South America and evaluated the representation of those sets of cells within the protected area system. South American anurans showed high species richness in the Brazilian Atlantic Forest and East Tropical Andes, while regions harboring most of the rare species were concentrated in the Andes and Atlantic Coast from North-Eastern Brazil to River Plate. Based on such patterns, we identified as rich-rare cells the Brazilian Atlantic Forest and Tropical Andes and as poor-rare cells the southern part of Andes and Uruguay. A low fraction of both sets of cells was represented within the protected area system. We show that a simultaneous consideration of species richness and rarity provides a rapid assessment of large-scale biodiversity patterns and may contribute to the definition of conservation priorities.


Introduction
The ongoing biodiversity crisis claims for conservation actions that help mitigate it [1]. Limited availability of resources for conservation demands the definition of priorities [2]. Among these priorities, identifying particular regions of conservation value helps guiding resource allocation and avoiding competition with other potential land-uses [3]. The definition of spatial priorities must be based on an assessment of the geographic distribution of biodiversity features of conservation interest (e.g., numbers of species, phylogenetic and functional distinctiveness or genetic diversity) [4][5][6][7]. This geographic assessment of biodiversity features represents the first step towards informed conservation actions [4,8].
Biogeography and macroecology aim to describe and explain spatial patterns of different biodiversity features. Results from biogeographical and macroecological approaches can provide the primary information for conservation assessments and planning [9]. Nonetheless, additional information (e.g., species' abundances, functional traits, phylogenetic relationships, ecosystem services, and region-specific sociopolitical issues) is still required for conducting detailed conservation actions [6]. Therefore, broadscale conservation assessments constitute only an initial stage of planning from which more detailed assessments and prioritizations (e.g., at smaller spatial scales) can then be conducted [4].
Based on different biodiversity features, various procedures have been used to define spatial conservation priorities [10]. Commonly used criteria include the total number of species, of rare or narrowranged species, and the number of threatened species at a site or region [11,12]. Species richness is a straightforward conservation target owing to its intrinsic significance for biodiversity definition and relative ease of monitoring [13]. Species' rarity, defined either by restricted geographic distribution or low population numbers, is of primary conservation concern owing to its relationship with species' threat status and extinction risk [14], thus the number of rare species is also used to establish spatial priorities.
Despite being relevant to virtually all conservation assessments and plans, species richness and rarity are seldom investigated simultaneously. Instead, they are used individually or comparatively to define different sets of species under study (e.g., all species vs. rare species richness) [12,15]. Moreover, prioritizations based on either richness or rarity have been criticized owing to their emphasis on only one or two aspects of biodiversity (e.g., richness and threat) [16]. At broad spatial scales, species richness and rarity combined can provide a rapid assessment of biodiversity patterns for spatial priority setting. Recently, a macroecological framework that simultaneously considers richness and rarity, measured by range size, was introduced to describe geographic patterns of biodiversity with applications to conservation priority setting [17][18][19][20]. This framework uses primary biological information (i.e., species' presence-absence data) within range-diversity plots that describe the relationship between species richness and range size, providing a straightforward way to identify if species-rich and species-poor regions are composed mainly by rare or common (i.e., geographically restricted vs. widely distributed) species [20].
Here we applied such macroecological framework to describe geographic patterns of biodiversity and identified regions of potential conservation value in South America. Also, we used these regions to evaluate the performance of the protected area system. In doing so, we aimed to inform conservation practices by providing a broad-scale conservation assessment that could be integrated into more detailed prioritizations. For this, we used amphibians as our case study. Amphibians are a particularly threatened vertebrate group, with ca. 32% of species undergoing a combination of hazards including habitat loss, climate change and emergent diseases [21][22][23][24][25]. Their vulnerability is boosted by their overall smaller range sizes, relative to other terrestrial vertebrates [15], which has been shown to correlate with amphibian threat status and extinction risk [25]. Most threatened amphibians are anurans (frogs and toads), which are the most species-rich and geographically widespread amphibian order [26]. Anuran biodiversity, including three of the main conservation aspects (species richness, rarity, and threat), is concentrated over tropical regions, especially in South America [25,26], making this a critical region for amphibian conservation.

Materials and Methods
We mapped the geographic distribution of South American anurans (i.e., extent of occurrence) drawn from the Global Amphibian Assessment [22] onto a 1u61u resolution grid. From this grid, we built a presence-absence matrix (PAM) of 1520 cells by 2437 species. The PAM is a binary matrix representing the basic biogeographical information on species richness and occurrence over a particular region [27]. We followed Arita et al.'s [17] Qr-mode to gather information on the number and mean range size of species present at each cell from the PAM. This procedure extracts data from rows (i.e., grid cells: species richness) by simultaneously considering the information within columns (i.e., species: range size). This information can be summarized and depicted using range-diversity (RD) plots [17,18].
RD plots convey joint information on species richness and range size enabling the exploration of biogeographic patterns of similarity among grid cells, in terms of shared species, and cooccurrence among species [18]. We used the ''by sites'' version of RD plots [17], in which axes represent the mean proportional percell range size and the proportional species richness of cells, in the xand y-axes, respectively. Central tendency of points within the RD plot by sites is determined by the average proportional range size of the whole system depicted by a vertical, dashed line; whereas detailed point dispersion depends on the overall covariance among cells resulting from the number of cells with which each individual cell shares its species [17]. In general, points located to the left of the plot's vertical line depict cells with negative average covariances (i.e., on average, sharing none or few species with most cells) whereas points to the right indicate cells with positive average covariances (i.e., on average, sharing many species with most cells). RD plots and associated parameters were obtained using the script available from [18] for the R statistical language [28].
We applied a quantile approach to define subsets of cells with different combinations of richness and range size values within the RD plot [20]. We defined ''rich-rare'' cells as those lying within the fourth quartile of species richness (i.e., cells with higher proportional richness) and first quartile of range size values (i.e., cells with lower mean range size). Additionally, we defined ''poorrare'' cells as those falling within the first quartiles of species richness (i.e., cells with lower proportional richness) and range size values. This quantile approach is a pragmatic criterion based solely on the frequency distribution of species richness and range size values of the study system, which has been commonly used to define species richness hotspots (e.g., [11,12]) or species' rarity (e.g., [15,29]). We mapped the geographic pattern of rich-rare and poor-rare cells, along with the overall pattern of anuran species richness and per-cell range size within South America.
We used the World Database on Protected Areas [30] to define protected grid-cells throughout South America. Protected cells were defined as those grid-cells with more than 10% of its area covered by protected areas (categories I-IV of IUCN). These protected cells were used to evaluate the performance of the South American protected area system in representing rich-rare and poor-rare cells. Significant spatial congruencies between our richrare and poor-rare cells and protected cells were tested using randomization tests. We quantified the observed spatial overlaps by counting the number of protected cells falling within our defined sets (i.e., observed value). Then we estimated the significance of these spatial overlaps by randomizing the position of protected cells (999 times) and recalculating the spatial overlaps between these and our rich-rare and poor-rare cells (i.e., random values). The proportion of randomized values equal or superior to the observed value was considered as the P-value of the randomization test. The null hypothesis that the observed overlap is lower than expected by chance was tested at the 5% significance level. Statistical analyses were done in R [28].

Results
Anuran species richness varied widely, from one to 161 species among cells (53.83634.20, mean 6 SD). Regions with the highest richness were the eastern tropical Andes, Amazonian basin, and the Brazilian Atlantic Forest (Fig. 1A). Rarity regions (i.e., low mean per-cell range size) concentrated in the tropical Andes and the Atlantic coast from North-Eastern Brazil to River Plate, but also in the southern Andes along the Austral temperate forests. These regions included both species-rich and species-poor cells (Fig. 1B).
Species richness frequency distribution (i.e., the distribution of richness values over all cells) showed a bimodal pattern, peaking at low and intermediate richness values with few cells having high species richness (Fig. 2A). The vast majority of anurans (97.7% of species) occupied less than one quarter of South America and only a few species (0.6% of species) occurred in more than 50% of the continent (mean = 33.6 grid-cells, SD = 110.2). In fact, 35.7% of species occurred in only one grid-cell. This pattern yielded a highly right-skewed range size frequency distribution; with far more species with small range sizes than species with large ranges (Fig. 2B).
The RD plot lacked a clear internal structure, with points arranged towards the lower left corner of the plot but with considerable dispersion along the x-axis (Fig. 3B). Most cells showed positive covariance, indicating that they generally share species (Fig. 3B). Nevertheless, there was a high turnover of species as shown by the Whittaker's beta (b w = 45.27), whose reciprocal value (1/b w = 0.022) equals the average proportional species richness and the average proportional range size in the system [17]. This means that cells contained, on average, 2.2% of species; and that the average anuran species occurred in 2.2% of cells.
The rich-rare set comprised 71 cells harboring more than 79 species and averaging less than 254 cells in per-cell range size. Rich-rare cells were distributed almost exclusively in mountainous areas over the tropical Andes and the Brazilian Atlantic Forest. Conversely, the poor-rare set comprised 255 cells harboring less than 31 species and averaging less than 254 cells in per-cell range size. Poor-rare cells were located mainly over the western-central and southern regions of South America along the southern Andes but also in Uruguay (Fig. 4).
Rich-rare and poor-rare cells were poorly represented by the protected area system within South America (Fig. 4). Less than 25% of cells within both rich-rare and poor-rare sets (21.13% and 16.47% of cells, respectively) contained at least one protected cell within them. This pattern was not significantly different from random placements of the observed number of protected areas within South America (P = 0.463 and P = 0.944 for the rich-rare and poor-rare sets, respectively).

Discussion
We presented the first comprehensive analysis of anuran biodiversity patterns throughout South America simultaneously considering two of the most important biodiversity features: species richness and range size. This biogeographic description allowed us to identify regions within South America that may be relevant for conservation of anuran biodiversity. For instance,  these identified regions can be used as geographical templates where further attention can be directed.
Species richness can be positively related with species rarity when cells occupied by rare species are also species-rich and represent a subset of cells occupied by widespread species [18,20]. Our results showed that rare anuran species occurred at both species-rich and species-poor regions within South America. Regions with the highest species richness, namely the equatorial Andes and the Brazilian Atlantic Forest, tended to have species averaging small range sizes, whereas regions with intermediate species richness like eastern Amazonia and the Cerrado biome (i.e., Brazilian Savanna) harbored species averaging large ranges. Interestingly, regions with low species richness also presented small mean range sizes, especially those at the southern tip and westerncentral South America. These latter findings are in contrast with the common assumption that restricted or endemic species are located mainly in species-rich regions [31,32].
The partial, positive relationship between richness and rarity found for South American anurans could be highly informative to prioritize regions for conservation of this vertebrate group. For instance, geographic congruence between high richness and high rarity regions could guarantee the protection of rare species if only the richest cells were selected for conservation. Indeed, previous studies have supported this richness-rarity relationship. For example, global analyses of terrestrial vertebrates showed that regions with high richness of narrow-ranging species were congruent with species-rich regions, suggesting that rarity could be used as a surrogate for the conservation of overall species richness [15,31]. Moreover, results from those analyses were similar for different taxonomic groups, supporting the role of indicator groups in which richness of rare species from one taxonomic group correlates positively with overall richness of other taxonomic groups (e.g., [31,32]). Accordingly, anurans seem to be a good indicator of biodiversity patterns of other terrestrial vertebrates (e.g., [32]). Therefore, the observed richness-rarity patterns we found could facilitate spatial prioritization for South American vertebrates in general.
Conversely, the observed relationship between species-poor regions and rarity for South American anurans challenges conservation actions. Many cells held low numbers of anuran species with small range sizes. Such cells represent unique assemblages usually not incorporated into common conservation actions focused solely on species richness. Consequently, these species-poor cells or ''coldspots'' [16] may be left out of conservation actions, though they are certainly relevant. Indeed, more comprehensive approaches do consider such species-poor cells [6,33]. For instance, under a complementarity criterion, if a cell holds a low number of species that do not occur elsewhere, this cell will be important in determining the network of cells needed to protect all species [6,7,34]. Yet, our simple approach considering both species richness and range size at large spatial scales would help identifying such coldspots in a straightforward manner.
Identification of poor-rare cells may also point towards particular ecosystems or habitats usually downplayed in most conservation strategies owing to their limited species diversity [16]. Poor-rare cells for anurans in South America concentrated at the southern tip and western-central regions, which are covered by temperate forests and steppes. These vegetation types are indeed poorer in species richness when compared to tropical forests. Nevertheless, the presence of species restricted to those temperate forests and steppes habitats highlight their relevance for biodiversity conservation. In addition, following the reasoning of indicator groups, it is possible that coldspots for one taxonomic group be so for other taxonomic groups too. This possibility could be evaluated in the same manner as the identification of species-rich regions for different taxa [15,29,31].
Recent studies have highlighted the importance of Neotropical ecoregions for the conservation of terrestrial vertebrates, including anurans [35][36][37]. For example, Loyola et al. [35] identified a set of 49 ecoregions that if sufficiently covered by reserves would retain at least 80% of total, endemic and threatened vertebrate species altogether. Some of their priority ecoregions are spatially congruent with our rich-rare and poor-rare cells. For instance, our rich-rare cells were located along the Serra do Mar Coastal forests and Chocó-Darien moist forests, whereas our poor-rare cells fell within the Patagonian steppe, Valdivian temperate forests and Uruguayan savanna. The congruence between Loyola et al. [35] results and ours is not surprising because theirs relied on a complementarity criterion that is bound to select regions with rare species even if they are species poor. Instead, such spatial congruence highlights the advantage of our simple approach for conducting initial assessments of priority regions.
The current protected area system in South America fails to encompass important regions for the conservation of anurans. For instance, large regions with few rare anuran species along the southern Andes in western-central South America and central Uruguay are virtually lacking protected areas. These coldspots would be overlooked by protected area systems if conservation assessments fail to apply more integrative approaches. Urbina-Cardona and Loyola [38] obtained similar findings for the representation of endangered hylid frogs (tree-frogs) within the Neotropical protected areas, showing that most restricted-range tree-frogs have only the periphery of their distributions under protection [38]. Our findings and those of other authors (e.g., [38]) on the scarce protection of South American anuran biodiversity are not surprising. Instead, these results highlight the inefficiency of the ad hoc manner of establishing protected areas to preserve biodiversity of specific taxa [39,40].
Aside from providing objective planning frameworks [9], the application of biogeographic methods and theory allow pattern interpretation and inference of potential processes responsible for such patterns [41,42]. For instance, amphibian biodiversity patterns have been causally related to current environmental conditions (e.g., energy: water and temperature; [43,44]) but also to historical processes (e.g., historical climate oscillations, time for speciation, niche conservatism; [45,46]). Likewise, rich-rare and poor-rare regions identified by our approach can be interpreted as a result of historical processes. For example, processes like lineage origin and speciation followed by extinction within a region or colonization of few lineages followed by low speciation could have given rise to poor-rare assemblages [47]. Conversely, rich-rare regions could be the result of high speciation coupled with low extinction rates [47]. Evaluating the phylogenetic structure of assemblages [48] within rich-rare and poor-rare regions could provide tests for such hypotheses. For example, if the aforementioned processes hold, we could expect phylogenetic clustering (i.e., presence of closely related species) in poor-rare regions and phylogenetic overdispersion (i.e., presence of distantly related species) in rich-rare regions.

Conclusions
Our approach described complex biogeographic patterns of South American anurans to inform conservation prioritization. Other biological data are undoubtedly necessary (e.g., species' traits, phylogenetic relationships) for a proper establishment of conservation priorities, but these are difficult to obtain and not always available. Thus, we believe that our spatial assessment of broad-scale biodiversity patterns of South American anurans will prove useful to inform future conservation planning and practices based on current knowledge for this threatened vertebrate group.