aSchool of Geography and the Environment, University of Oxford, Oxford OX1 3QY, United Kingdom;bCentre for Macroecology, Evolution and Climate, Department of Biology, University of Copenhagen, Copenhagen, Denmark;

Significance

Biogeographic theory builds upon a long history of analyzing species-diversity patterns of remote islands, but no previous studies have attempted to investigate corresponding patterns in functional traits on islands. Our analyses of functional diversity (FD) for spiders and beetles in the Azorean archipelago reveal that FD increases with species richness, which, in turn scales with island area regardless of the taxa and distributional group considered (endemics, natives, and exotics). Our results also support the hypothesis that each group contributes to FD in proportion to their species richness and that, being dominant, exotic species have significantly extended the realized trait space of the Azorean islands. Further analyses in other archipelagos are needed to establish whether our findings are representative of oceanic islands.

Abstract

Analyses of species-diversity patterns of remote islands have been crucial to the development of biogeographic theory, yet little is known about corresponding patterns in functional traits on islands and how, for example, they may be affected by the introduction of exotic species. We collated trait data for spiders and beetles and used a functional diversity index (FRic) to test for nonrandomness in the contribution of endemic, other native (also combined as indigenous), and exotic species to functional-trait space across the nine islands of the Azores. In general, for both taxa and for each distributional category, functional diversity increases with species richness, which, in turn scales with island area. Null simulations support the hypothesis that each distributional group contributes to functional diversity in proportion to their species richness. Exotic spiders have added novel trait space to a greater degree than have exotic beetles, likely indicating greater impact of the reduction of immigration filters and/or differential historical losses of indigenous species. Analyses of species occurring in native-forest remnants provide limited indications of the operation of habitat filtering of exotics for three islands, but only for beetles. Although the general linear (not saturating) pattern of trait-space increase with richness of exotics suggests an ongoing process of functional enrichment and accommodation, further work is urgently needed to determine how estimates of extinction debt of indigenous species should be adjusted in the light of these findings.

In comparison with less-isolated settings, the biotas of oceanic island archipelagos are species poor, disharmonic, endemic-rich, and particularly sensitive to disturbance (1). Founded by strongly dispersal-filtered sets of colonists, some of which have since diversified, their precontact biotas represent the dynamic outcome of immigration, speciation, and extinction processes. Following human contact, they have experienced waves of anthropogenic species immigration, ecosystem transformation, indigenous species population reductions, and extinction (1). We may sustain such broad generalizations by reference to many case studies but, in so doing, disregard some intense debates within island biogeography: for instance, whether prehuman contact biotas can be regarded as equilibrial, what rules govern disassembly of precontact assemblages and the assembly of the novel combinations of indigenous and exotic species, how long it will take for the many ecological adjustments to play out before the tally of extinctions of declining indigenous species is clear, or indeed whether the net effect of anthropogenic interference (land-use change, species introductions, etc.) may be higher levels of species diversity (1, 2).

Although traditional island biogeographic analyses sensu MacArthur and Wilson (3) have identified notable differences between the emergent diversity patterns of endemic, other native, and exotic species across islands (4⇓–6), they have so far failed to resolve the questions posed above. This failure may in part be because they mostly rely upon analyses in which species are treated as equivalent (but see, e.g., ref. 7). In contrast, case-study evidence suggests that the introduction of particular species possessing novel trait characteristics (e.g., mammalian predators, large herbivores) can have dramatic and disproportionate impacts on island ecosystems and assemblages (1, 8, 9). Such observations suggest that analyses based on functional traits (“components of an organism’s phenotype that influence ecosystem level processes”; see ref. 10, p. 742) may provide especially valuable insights for island biogeographic theory (and perhaps vice versa).

Herein, we assess the scaling of functional diversity (FD) with island area and species richness (SR), testing for nonrandomness in the contributions of differing distributional groups of species. For exploratory purposes, we used the following distributional categories: endemics, natives (excluding endemics), indigenous (endemics plus natives), exotics, and all species, noting that the overlaps between categories renders certain tests nonindependent. Our study system comprises the spiders and beetles of the Azorean archipelago, a uniquely well-specified oceanic island system. The Azores were originally almost fully forested, but, since human colonization (ca. A.D. 1440), native-forest cover has been reduced to <5% (11). Spiders are top predators among invertebrates and present an array of potentially important traits related to habitat structure and resource use. Beetles are the most diverse insect Azorean order (as they are globally), include a wide range of feeding ecologies and are prevalent in most habitats. Current estimates are that 22 spider species are endemic, 17 native, and 85 exotic, with respective values for beetles being 65, 134, and 325 (12, 13). It is likely that several endemic species are already globally extinct (14) and analyses of the Azorean “extinction debt” based on calibrated species-area models predict further extinctions to follow (2). However, the logic of those extinction-debt analyses is based on reductions in habitat area of native forest without adjustment for the role of exotics.

Several FD indices have been developed for the analysis of interspecific variation in functional traits for whole assemblages (15, 16). However, there have been few applications of these methods to island faunas (but see, e.g., ref. 8). More specifically, we are not aware of any work quantifying at island scale how functional-trait space of invertebrates has been changed, and with what consequences, by anthropogenic alteration of oceanic islands. Our goal is thus to explore the absolute and relative contribution of endemic, native, and exotic invertebrate species to functional-trait space at island and archipelagic scales in a system of nine oceanic islands spanning from 17 to 757 km2 in area. We also test for nonrandomness of FD contributions at finer resolution within the remaining native-forest enclaves.

To estimate FD, we use FRic (a measure of functional richness) (17) calculated from input data based on seven traits for spiders and three traits for beetles. Given the broad differences in ecology (encompassing trophic status, flight ability, etc.) between these groups, the traits selected necessarily differ: the spider traits relate to size, feeding behavior, and habitat use whereas the beetle traits relate to size, feeding guild, and flight capabilities (Materials and Methods). We test for evidence of trait-space saturation with increasing island area and species richness for each distributional group and for all species (Materials and Methods and SI Appendix, Table S1). We also test for the contribution made by, and degree of overlap between, each distributional group by means of null-model simulations. We thus address a question of fundamental island- and conservation-biogeographical interest (e.g., refs. 18 and 19): Are the exotic species occupying and thus competing within essentially the same functional space as indigenous groups, or are they serving to expand the functional biodiversity of the archipelago?

Island species–area relationships (ISARs) typically show a positive trend, with slopes (z values) indicative of the process(es) establishing species richness and composition patterns (1, 5, 20). As we move from speciation-dominated systems (e.g., oceanic islands) to immigration–extinction dynamics (e.g., continental-shelf islands) and then to systems of little dispersal-limitation (e.g., inland islands), we generally observe lower values of z, as is also the case when comparing narrowly endemic species (higher z) to nonendemic natives (lower z) within oceanic archipelagos (5, 20). Recent work on functional traits has shown that FD, when measured as functional richness, is typically strongly related to SR (10, 17). Combining these insights, we next consider the implications for patterns of FD in an oceanic archipelago. As remote oceanic islands have strongly colonization-filtered (“disharmonic”) biotas, we predict that, in the absence of large-scale trophic collapse driven by ecosystem transformers (9), anthropogenic introductions should add FD disproportionately with reference to native species and endemics and should therefore significantly expand the occupied functional-trait “space” of the archipelago. Because of the expected strong scaling between FD and area, we predict that smaller islands, possessing the fewest indigenous species, will see the greatest proportional gains in FD through the introduction of exotic species even if (compared with larger islands) comparatively few exotic species manage to arrive (less economic activity, less transport to the island) and establish (smaller human settlements) on those small islands. Higher proportional gains in smaller islands are expected because, with increasing richness (larger islands, more introductions), the chances of each new arrival adding truly novel trait combinations to the system diminishes. Therefore, the expansion of FD space by the addition of exotic species is predicted to decrease with island size. Based on the foregoing, we would predict that (i) log–log ISARs should increase in slope (z values) from exotics to natives and especially to endemics whereas (ii) the equivalent island functional diversity–area relationships (IFDARs) should increase in slope from exotics to indigenous species, with little difference between natives and endemics, as the latter groups experienced the same colonization filter in reaching the archipelago. By contrast, the rank order of “intercepts” (c values) may not be so readily predicted as a function of distributional groups, as much depends on the age and isolation of the islands and the extent of human transportation (1, 5).

Our analyses provide no support for the saturation of trait space at island level: Exotic species add FD approximately in proportion to increasing SR. This result holds for both spiders and beetles although exotic spiders appear to contribute more functional novelty than exotic beetles. Our predictions for the form of ISARs and IFDARs met with only equivocal support (for spider ISARs and beetle IFDARs). Finer scale analyses, restricted to species occurring in native-forest remnants, provide some intriguing indications of habitat filtering on particular islands for beetles, but not for spiders. Our results indicate that anthropogenic habitats and species introductions are combining to permit the persistence in these highly altered oceanic islands of much higher SR than might be anticipated from traditional applications of island theory while processes connected with habitat loss nonetheless continue to drive the attrition of the indigenous fauna (2, 14, 21).

Results

We first established that the significant log–log ISARs previously established for both spiders and beetles (2) hold, with some variation, for the distributional subsets (Fig. 1 and SI Appendix, Tables S2 and S3). For spiders, endemics display a steep ISAR, and natives a flat (nonsignificant) ISAR, but when the two categories are combined as “indigenous” species, the resulting ISAR is indistinguishable in slope from that of the exotic category whereas exotics have a higher intercept, reflecting their greater SR (Fig. 1 and SI Appendix, Table S3A). For beetles, the distributional subsets have indistinguishable slopes but differing intercepts [analysis of covariance (ANCOVA)] (SI Appendix, Table S3B). Having established that the traits we used are mostly complementary rather than redundant with one another (SI Appendix, Tables S1 and S4), we tested for the equivalent island functional diversity–area relationships (IFDARs). The IFDARs show that FD generally increases with island area, with the exception of the endemics subsets, for which no significant relationship was obtained (Fig. 1 and SI Appendix, Table S3A). For spiders, the slopes of the IFDARs for different subsets are indistinguishable, but the intercepts differ, whereas the between-group differences for beetles involve both intercept and slope (Fig. 1 and SI Appendix, Table S3B). Thus, in general, SR and FD each scale with island area, consistent with classic island theory (1, 3, 20) while statistical support for our specific predictions regarding ISAR and IFDAR form was limited to spider ISARs and beetle IFDARs.

Comparison of linear and polynomial model fits for the FD–SR relationship provides a test for saturation in trait space with increasing SR across the nine islands (eight for endemics) (Fig. 2). Linear models are preferred [based on small sample size-corrected Akaike information criterion (AICc) values], with the following exceptions: for endemic spiders and native beetles, neither model is adequate (F statistic tests, P > 0.05); for endemic beetles, the ΔAICc between the linear and polynomial models is <2, indicating equivalence in fit, and, for native spiders, the fitted polynomial indicates increasing FD with higher SR rather than a tendency toward asymptote. Comparison with a sensitivity analysis indicates the patterns for endemics to be equivocal and unstable (Fig. 2 and SI Appendix, Fig. S1); in addition, a separate sensitivity test established that linear models are preferred for each distributional group of spiders upon reduction of the number of traits from seven (main analyses) down to three (as for beetles) (SI Appendix, SI Materials and Methods, and Table S5).

The relationship between functional richness values (FRic, ref. 17) and species richness (SR) for Azorean spiders (Upper) and beetles (Lower) for each distributional category, where each data point represents an island (n = 9, except for endemics, where n = 8 as Corvo Island was excluded). Linear (full lines) and polynomial (dotted lines) fits and their respective ΔAICc and R2 values are given. Bold lines indicate that the model is significant (F statistic test P < 0.05) and has the lowest ΔAICc value. In all cases, apart from endemic beetles, this is also the single “best” model and the preferred model (based on ΔAICc). Note that the polynomial model for endemic beetles is unrealistic in indicating negative FRic values, meaning that, even in this case, the linear model is preferable on grounds of ecological realism. FRic was computed based on six PCoA axes (81% total inertia) and seven PCoA axes (92%) for spiders and beetles, respectively.

We next used a simulation approach to create null distributions of expected FD per island for each observed SR value. We report two such simulations (SI Appendix, Table S6). First, based on the entire archipelagic species pool, we found no deviation from the null expectation for either taxon for any distributional grouping; i.e., observed FD is simply a function of SR. Our second simulation was restricted to species occurring in the native forest (SONF), to assess the outcome of species assembly processes in the habitats where most remaining threatened endemic species occur. In this case, spider values did not depart from expected, but significantly lower than expected values were observed for six islands for native beetles (Faial, Pico, S. Jorge, S. Maria, S. Miguel, and Terceira) and for three islands for exotic beetles (S. Maria, S. Miguel, and Terceira). Lower values imply that, for the given SR, the FRic value is lower than expected were these species a random draw from the archipelagic SONF pool.

To further interrogate the trait-space occupancy and overlap among endemics, natives, and exotics, we undertook analyses of the whole archipelagic species pool (rather than per island) and estimated FRic as before, but using a slightly lower proportion of the derived trait data (Materials and Methods and SI Appendix, Table S7 and Fig. S2). For spiders, 99.74% of the trait space is occupied by exotic species, 36.24% by natives, and 20.64% by endemics (Table 1). The endemic trait space is nested within that of natives, and natives within the exotics. By contrast, endemic (83.1%) and native (89.8%) beetles occupy far larger proportions relative to exotics (95.8%) although, again, the space occupancy is largely nested. The null simulations indicate that each distributional category contributes to the overall trait space as expected as a simple function of their SR. There are two marginal exceptions: First, endemic spiders occupy less functional space than expected and, second, the overlap between endemic and exotic spiders is less than expected. One corollary is that spider mean body size increases from endemic to native to exotic spiders (although only endemic body size is statistically distinct) whereas, for beetles, the three distributional categories are not distinguishable by body size (SI Appendix, SI Results, Body Length). These results indicate that the emergent patterns in FD can be related back to the underlying functional-trait data (compare SI Appendix, Fig. S2). Moreover, even when analyzing a single, common trait, there are differences between the two taxa in the relative contribution to FD made by endemics, natives, and exotics. The patterns in body size are also broadly as expected as larger bodied spiders are relatively unlikely to reach such remote islands naturally by ballooning, but may readily do so by human transportation (cargo in boats and planes). In contrast, beetles are less likely to demonstrate such a strong body size-related filter due to their active flight capability.

Functional-space occupancy and overlap per distributional group, for all Azorean species (n values in SI Appendix, Table S2) based on FRic estimated from the first four PCoA axes using a Gower distance matrix

As the within-island habitat occupancy of the different groups may be crucial to the ecological interpretation (e.g., ref. 21), we undertook a further randomized simulation analysis of the species occurring in native forest, analyzing the degree of overlap in FD between exotic SONF and indigenous SONF and how this overlap compares with that expected if the exotic SONF were a random draw from the overall set of exotic species. This test showed only one deviation from the null expectation across the nine islands (Table 2). Only for beetles on Terceira Island is the overlap marginally lower than expected by chance, implying a degree of functional divergence between exotic SONF and indigenous SONF (Table 2). Thus, based on the traits assessed here, the assemblage of exotic species that has been observed in native-forest habitat in each island is (with that one exception) essentially a random subset of the exotics pool and not a pool of species filtered with respect to their similarity with indigenous species.

Deviations expressed as standardized effect sizes (SESs) from the null expectation (given observed SR) for the functional overlap between exotic species occurring in native forest (SONF) and indigenous SONF for Azorean spiders and beetles species, based on the PCoA analyses reported in Table 1

Discussion

For most distributional groups, we obtained a strong positive relationship between island area and species richness (SR), with exotic SR exhibiting a similar but elevated ISAR to indigenous species, while FD also increases in proportion to both island area and SR across the nine islands of the archipelago. Thus, the principal effect of increasing SR by the addition of exotics is to generate a corresponding increase in FD, with no evidence of saturation. Whereas the ISARs vary mostly in their intercepts, the form of the IFDARs indicates that exotic spiders contribute an increasingly large proportion of FD with increasing island log-area. The pattern for beetles differs, with indigenous species showing comparable FD to exotics across the larger islands (Fig. 1).

Our analyses of functional space occupancy and overlap show that exotic spiders have contributed much more to the increase in functional trait-space occupancy than have exotic beetles, both in absolute terms and in proportion to increasing SR (Table 1, Fig. 1, and SI Appendix, Fig. S3). Exotic spiders have added new taxonomic diversity to the Azores. For example, 13 of the 26 spider families now found on the archipelago are wholly comprised of exotics (13), including species bringing completely new foraging strategies to the regional functional space. For beetles, a rather lower proportion—18 of 54 families—are represented only by exotics (13), consistent with their making a lesser contribution to functional space across the archipelago. In summary, the most striking emergent pattern is that, for both taxa, FD scales up with the area and number of species per island, indicating that the present-day Azorean fauna appears to be largely unsaturated in terms of functional-trait space and consequently remains vulnerable to further colonization by newly introduced exotics (22).

These observations of expanding “biodiversity” contrast with claims (2, 14) of past and impending extinctions as a result of habitat alteration, yet these scenarios may be easily reconcilable with reference to the massive anthropogenic changes since colonization ca. A.D. 1440. The native habitats, on which many indigenous arthropods are dependent (23), have been reduced across the archipelago to <5% of their original cover, and an array of new agricultural, forestry, urban, and semiurban habitats have been created in their place. Concurrently, large numbers of exotic species have been introduced by anthropogenic means, many of which occur wholly or predominantly in novel, anthropogenic habitats (24). The resulting ecological relationships between indigenous and exotic species must span the full gamut from closely co-occurring and interacting, to nonexistent, depending on the ecology and habitat occupancy of the particular species under consideration.

Previous work has shown evidence of source–sink dynamics operating between native and anthropogenic habitats in these islands. This process operates in both directions, with some endemic spiders and very few endemic beetles (see, e.g., ref. 23) spilling into novel sink habitats surrounding native forest whereas intensive pasture is the main habitat source of exotic species invading native forest (25). The success (i.e., abundance and spread) of such exotic “invasive” species has been found to vary between soil and canopy habitats. In particular, the canopies of Azorean endemic trees seem to be a specialized habitat, and few exotic insects and spiders maintain “source” populations there (12, 25, 26). These findings suggest that there may be assembly rules determined by a combination of functional traits and characteristics of the habitat template at local scales of analysis. In the present coarser-scaled analyses, we found no evidence for saturation at the island scale, and, when analyzing overlap between exotic and indigenous species occurring in native forest (SONF), we again found no deviation from expected, apart from the single case of beetles on Terceira (for which functional divergence was found). We did find slightly more extensive indications of lower FD than expected as a function of SR for native SONF in six islands and exotic SONF for three islands for beetles (Table 2 and SI Appendix, Table S6). These findings are intriguing. However, when endemic and native species are recombined to represent the whole indigenous assemblage of SONF, no departure from random is evident. The three islands showing lower than expected FD for exotic SONF are S. Maria, S. Miguel, and Terceira. These islands are the three richest islands in terms of exotic species and, on the basis of this test, would appear to show some evidence of saturation of native forest for exotics. Alternative explanations are possible; for example, these results might reflect the existence of large numbers of urban beetle species, or of recent arrivals that have yet to reach and enter the remaining native-forest fragments (24). It thus appears that multiscale analyses will be necessary, assessing species abundances and population trajectories, to determine how exotic species are influencing survival prospects of indigenous species in the remaining forest patches.

We recognize that FD values can be sensitive to properties such as the number and range of traits considered and that the outcome of certain FD analyses can also be sensitive to how the species pool is defined (15, 27). In the present analysis, we necessarily began by using the archipelagic species pool, as we lack a more realistic basis for defining a regional species pool. It is also important to note that we cannot be certain of the number, or trait selectivity, of extinctions of indigenous species that may have occurred without record, and which may have impacted on the functional trait-space occupancy and FD of indigenous groups (2, 13, 14) (SI Appendix, SI Materials and Methods). Moreover, a fuller picture of trait occupancy requires the use of additional FD indices that in essence interrogate the occupancy of the convex hull space analyzed herein (e.g., functional evenness or functional divergence) (17), and which make use of species-abundance data: unfortunately such data are lacking at a whole-island scale.

To conclude, it is apparent that exotic species have significantly expanded the realized trait space and FD of the Azorean islands, generating emergent FD patterns strikingly consistent with those of the indigenous fauna (compare ref. 4), while remaining largely dependent on the extensive anthropogenic habitats that have replaced the native forests. These findings suggest that, whereas island area strongly constrains both species richness and FD, the diversity levels observed are not indicative of a fixed area-dependent species carrying capacity. Rather, the capacity attainable has increased as the biological isolation of the islands has been progressively diminished by human agency. Although there is no evidence of saturation of the anthropogenic habitat space, endemic species are threatened in their native-forest enclaves (2), in which a degree of habitat filtering is apparent, at least for beetles. Further analyses of the functional biogeography, for other taxa and archipelagos, at multiple scales of analysis, are needed to establish whether the findings reported herein are representative of oceanic islands (compare ref. 28). In addition, long-term monitoring data are required to determine the extent to which indigenous species of the Azores can persist within the now much larger pool of species present on the islands.

Materials and Methods

In all cases, we undertook and reported separate analyses for spiders and for beetles, conducting statistical analyses implemented within the R programming environment (29) and the software FDiversity (30). Where significance tests are reported, the critical value used was <0.05.

Study Area.

The Azores are volcanic, oceanic islands spanning the mid-Atlantic ridge, roughly 1,600 km from the European mainland. Humid evergreen broadleaf laurel forest (laurisilva) covered most of the land area before human colonization ca. A.D. 1440. Around 95% of these forests have since been destroyed by humans (11). The current major vegetation habitats across the archipelago comprise native forest, exotic forest, seminatural pasture, and intensively managed pasture (2, 21, 23, 25).

Functional Traits and Distributional Data.

Data on the occurrence of species in each island and categorizations as endemic, native, or exotic were drawn from the latest compilations (13). We collated and/or recorded information for each spider species for (i) body length (average between females and males), (ii) web use, (iii) web architecture, (iv) foraging strategy, (v) prey range, (vi) vertical stratification, and (vii) circadian activity; and for each beetle species (i) body length, (ii) feeding guild, and (iii) wing morphology (see SI Appendix, SI Materials and Methods for details). Except for body length, all traits were coded as nominal variables. For analyses focused on species occurring in native forest (hereafter SONF), we used our distributional data (SI Appendix, SI Materials and Methods) to denote exotic SONF and indigenous SONF. As some of the exotic SONF may have been “tourists” or “sink populations,” we avoided terming this group “invasive species.”

Computing FD.

Functional diversity (FD) was computed using the functional richness index FRic (17), which measures the multidimensional volume of trait space defined by using a convex hull volume. FRic does not require abundance data (which are unavailable in this study) and is not sensitive to species splitting; it is also effective for discriminating assembly rules (e.g., ref. 15) (see SI Appendix, SI Materials and Methods for details). First, Gower’s distance was used to calculate multivariate distances between species based on the raw trait data. To avoid negative eigenvalues, the Gower’s distance matrix was subject to primary square-root transformation (31). These distance data were then subject to principal coordinates analysis (PCoA)s and the resulting PCoA axes were used as new composite functional traits. FRic is then based on the minimum convex hull volume that contains all of the species selected for analysis (e.g., all endemics, exotic species on Faial Island, etc.). As the index requires that there are more species than there are traits (17) and in certain categories we have small numbers of species, we used only the first six (spiders, 81% total inertia) and seven (beetles, 92% total inertia) PCoA axes to estimate FRic for the island-by-island analyses. For endemics, as Corvo island has only three endemic spider and four endemic beetle species, we ran these analyses without Corvo (for subsidiary analyses including Corvo, see SI Appendix, SI Results). Due to computational restrictions, for analyses at the archipelago level (where very large numbers of species are involved), we used only the first four PCoA axes. We calculated FD (i.e., FRic) for five species groups: (a) endemics, (b) natives (i.e., excluding endemics), (c) indigenous (i.e., a plus b), (d) exotics, and (e) all species.

Relationships with Island Area.

We estimated the log–log island species–area relationship (ISAR) and the equivalent island functional diversity–area relationship (IFDAR) for each distributional category. The significance of both linear and quadratic models was assessed using the F statistic. Differences between the slopes and intercept between categories were estimated by analysis of covariance (ANCOVA), first considering endemics, native nonendemic, and exotics and, second, indigenous and exotic species.

Testing Redundancy.

Potential presence of saturation/redundancy, for each grouping, was analyzed by means of the interisland SR–FD relationship. We used linear modeling and evaluated the relevance of a quadratic term (polynomial) to account for potential curvilinear relationships (i.e., indications of saturation) (32, 33). The significance of both linear and quadratic models was assessed using the F statistic. The additional information provided by the quadratic model was assessed on the basis of percentage of explained variance (R2) and the small sample size-corrected AICc. The most parsimonious model has the lowest AICc, and thus ΔAICc = 0. Models with ΔAICc between 0 and 2 were considered to have equivalent empirical support, and models with ΔAICc > 2 were considered to have less support (34). Second, to distinguish whether FD was simply a product of SR, we used a simulations approach to create null distributions of FD for each observed SR value. Keeping SR constant for each island and category, we randomly selected species from the global pool of species without distinction between endemic, native, and exotic species. We repeated this procedure 999 times to produce a distribution of null FD values. The final null distribution included the 999 random values plus the observed value. Based on this null distribution, we computed the standardized effect size (SES) for each island to measure deviation of the observed FD from the null distribution. The P value associated with our null hypotheses, i.e., no deviation from a random assortment for species, was estimated per island. To avoid type I errors, we applied a false discovery rate analysis (FDR) (35) to account for multiple test comparisons.

Estimating Functional Space and Overlap.

Overlap between endemic, native, and exotic functional space at the archipelago scale was estimated following ref. 36. The multidimensional volume (i.e., convex hull) for each set of species was estimated, allowing quantification of the intersection between two (or more) species assemblages. The overlap is maximal when the two functional spaces are equally similar or one is nested in the other and minimal (zero) when assemblages do not intersect in their functional spaces. Due to computational restrictions, analyses for the entire archipelagic pool of species were run using only the first four PCoA axes (see Computing FD), summarizing 71% and 77% of the total inertia present in the initial Gower distance matrix for spiders and beetles, respectively, and for each pairwise comparison of endemics, natives, and exotics. To determine whether outcomes are simply a function of SR variation between categories, two different hypotheses were tested: (i) Does the functional space for each category differ significantly from a random expectation for the given SR? (ii) Does functional pairwise overlap differ significantly from a random expectation for a given difference in SR? Expected distributions were obtained by 999 randomizations of the label of the distributional categories for each species (i.e., endemics, natives, and exotics) keeping SR constant. Statistical significance of the observed values was estimated by calculating the P value from the null distributions.

Functional Overlap Between Exotic and Indigenous Species Occurring in Native-Forest Habitat.

To test whether the exotic species co-occurring with indigenous species in native-forest habitat are a random draw from the overall set of exotic species, we undertook a further test focusing only on species occurring in native forest (SONF). Expected distributions of the overlap between exotic SONF and indigenous SONF were obtained by measuring the overlap between indigenous SONF and random selections of all exotic species 999 times, keeping richness of the exotic SONF per island in the native habitat constant. Statistical significance was estimated as above.

(2008) Do non-native species invasions lead to biotic homogenization at small scales? The similarity and functional diversity of habitats compared for alien and native components of Mediterranean floras. Divers Distrib14:774–785.

(2010) in Canopy habitat area effect on the arthropod species densities in the Azores: Pondering the contribution of tourist species and other life histories. Terrestrial Arthropods Of Macaronesia: Biodiversity, Ecology And Evolution, eds SerranoARM, BorgesPAV, BoieiroM, OromíP (Sociedade Portuguesa de Entomologia, Lisbon), pp 89–114.

(2010) What does species richness tell us about functional trait diversity? Predictions and evidence for responses of species and functional trait diversity to land-use change. Glob Ecol Biogeogr19:423–431.

Researchers report biparental inheritance of mitochondrial DNA in 17 members of three unrelated multigeneration families, paving the way for insights into alternative mechanisms for the treatment of inherited mitochondrial diseases.

Researchers report a machine-learning approach to identify land plants at risk of extinction, suggesting that the approach can be used to guide policies aimed at allocating resources for biodiversity conservation.

A study explores how cats groom fur using fine structures called papillae on the surface of the tongue and presents a biologically inspired hairbrush to remove allergens from cat fur and apply medications on cat skin.