Hawaiʻi Gets a Hit by New Birds on the Block

Credit: Samuel B Case

Editor's Introduction

Structure, spatial dynamics, and stability of novel seed dispersal mutualistic networks in Hawaiʻi

annotated by

On the Hawaiian island of O’ahu, local plants rely on birds to eat and disperse their seeds. However, native birds have been replaced with new, invasive species. Do these invasive birds still distribute the seeds of native plants, or does their diet favor invasive plant species? Can these non-native neighbors maintain the complex, yet stable relationships once shared between native plants and animals despite not being closely related to the original birds of the island? This paper investigates how well invasive bird species of O’ahu adapt into their new roles as the primary seed dispersers of the island.

Paper Details

Original title
Structure, spatial dynamics, and stability of novel seed dispersal mutualistic networks in Hawaiʻi
Original publication date
Reference
Vol. 364, Issue 6435, pp. 78-82
Issue name
Science
DOI
10.1126/science.aau8751
Topics

Abstract

Increasing rates of human-caused species invasions and extinctions may reshape communities and modify the structure, dynamics, and stability of species interactions. To investigate how such changes affect communities, we performed multiscale analyses of seed dispersal networks on Oʻahu, Hawaiʻi. Networks consisted exclusively of novel interactions, were largely dominated by introduced species, and exhibited specialized and modular structure at local and regional scales, despite high interaction dissimilarity across communities. Furthermore, the structure and stability of the novel networks were similar to native-dominated communities worldwide. Our findings suggest that shared evolutionary history is not a necessary process for the emergence of complex network structure, and interaction patterns may be highly conserved, regardless of species identity and environment. Introduced species can quickly become well integrated into novel networks, making restoration of native ecosystems more challenging than previously thought.

Report

High rates of human-caused species invasions and extinctions are a ubiquitous feature of the Anthropocene (12). As a result, “novel communities” have emerged, characterized by a reshuffling of species, changes in species interactions, and, in some cases, alteration or disruption of ecosystem services maintained by these interactions (34). Mutualistic plant-animal networks are particularly susceptible to species loss (5) and invasions (467), increasing the vulnerability of species and communities to further perturbations (4). Previous studies have focused on native-dominated communities in which few or no invasive species occur and mutualistic partners have interacted for prolonged periods of time, developing complex and often coevolved interactions (89). By contrast, the architecture and stability of novel interaction networks across spatial scales and how they compare to native-dominated communities remain virtually unknown. This knowledge gap hampers our ability to forecast and mitigate the impacts of extinctions and invasions on ecosystem functions.

Here we address these gaps by examining the structure, dynamics, and stability to perturbations of multiple novel communities in the Hawaiian archipelago and compare our results with networks from communities worldwide. Hawaiʻi provides an opportunity to investigate the consequences of an extreme scenario of loss of native species and their replacement by non-native species. Most native Hawaiian forest plants are bird-dispersed, yet no native dispersers remain in most ecosystems (1011). Thus, seed dispersal is almost entirely dependent on a handful of introduced vertebrate dispersers, nearly all of which are birds (1011). Oʻahu, in particular, is among the areas most affected by extinctions and biological invasions in the world (12): All its native frugivores are extinct.

To what extent are introduced species integrated into seed dispersal networks (SDNs), and do introduced dispersers replace extinct native animals? To investigate these questions, we examined interactions based on 3278 fecal samples from 21 bird species [tables S1 to S3 and (13)] collected over 3 years at seven sites encompassing broad environmental variation across Oʻahu (fig. S1 and table S1). We identified 109,424 viable seeds, representing 1792 seed dispersal events (presence of viable seeds in a sample). Oʻahu’s SDN included 15 bird and 44 plant species connected by 112 distinct links (Fig. 1). Most birds (86.7%) and plants (65.9%) are not native to Hawaiʻi; introduced plants accounted for 93.3% of dispersal events, and there was no interaction between a native bird and a native plant. Proportions of introduced species varied from 60.0 to 100% for birds and 50.0 to 100% for plants, two local networks consisted entirely of introduced species, and the number of species and links was highly variable across sites (table S4). We found that 59.0% of fecal samples contained seeds (table S4), but only 0.22% of interactions (n = 4 events) involved native birds (two species not specialized for fruit consumption). Thus, although introduced birds are critical for seed dispersal in the ecosystem, they are primarily dispersing introduced plants (only 6.7% of interactions involved native plants).

 

HawaiiSeeds_1

Fig. 1 Structure of the island-wide seed dispersal network on Oʻahu and illustration of two emblematic interactions. (A and B) The novel network was nested [specialist species interacting with proper subsets of partners of the most generalist species (wNODF) = 48.67; 95% confidence interval (CI) = 34.24 to 46.66] (A) and modular [subsets of species interacting preferentially with each other, forming modules of highly connected species (QW) = 0.24; 95% CI = 0.07 to 0.09] (B). Species and links from distinct modules are depicted by different colors (blue, orange, and green), and gray links are interactions connecting modules. For a list of interacting species, see fig. S2. (C) Japanese white-eye feeding on Pipturus albidus, the most commonly consumed native plant. (D) A red-billed leiothrix feeding on Clidemia hirta, the most widely consumed and widespread introduced plant. [Illustration credit: P. Lorenzo]

Question 

What type of network patterns are formed between introduced bird species and native and introduced plant species?

Methods

The authors identified 44 different plant species and 15 different birds in O'ahu. Figure 1 depicts each bird (left column) and each plant (right column) as rectangles. The lines connecting a bird to a plant represent a seed dispersal event.

Results A

Networks such as seed dispersal networks consist of specialist species (which interacts with only a few, select species) and generalist species (which interacts with a broad range of other species). When specialist species interact with one (or a few) of the same species that the generalists interact with, it generates a nested organization of the interactions. Metaphorically, nestedness can be compared to a Russian Doll, where the diet of specialist species (smaller dolls) fit within the diet of the more generalists species (larger dolls). In Figure 1A, the specialists species are depicted as very thin rectangles to represent the few interactions they have with other species, whereas the generalist species are bigger rectangles to encompass the many interactions they have with other species.

Results B

A module is a group of species that preferentially interacts with each other instead of species outside the module. In Figure 1B, each module is represented by a distinct color.

Connectance measures the proportion of interactions taking place in a network out of the total amount of possible interactions. The low connectance observed means that not all possible interactions actually occur in the ecosystem.

We assessed species interaction patterns via complex network analyses and used four complementary metrics known to vary geographically and reflect community-level responses to major drivers of biodiversity patterns, such as productivity, climatic seasonality, and historical climatic stability [e.g., (1416)]. A network is an interaction matrix for which each row i is a plant species and each column j is a bird species, with intersections aij describing interaction intensity. The statistical significance of the observed topological patterns was assessed by contrasting observed values for each metric with the confidence interval from null models (13). Like other mutualistic networks, SDNs in native-dominated communities typically have consistent structures: (i) low connectance—not all possible interactions are realized; (ii) high specialization—few supergeneralist species exist and most species interact with a few partners in a complementary way; (iii) nested topology—specialist species tend to interact with subsets of partners of the most generalist species; and (iv) modular structure—subsets of species interacting preferentially with each other, forming modules of highly connected species (1720).

Novel insular communities are predicted to have low specialization because of niche broadening (21) and interaction release (22). For example, both fleshy-fruited plants and frugivores on islands tend to have wide niches owing to resource limitation (4). Consequently, high connectance and nonmodular structures are expected, because both are linked to low specialization [e.g., (1423)]. For nestedness, contrasting predictions exist because low specialization can either lead to non-nested topology, owing to random partners associating, or to nested topology, driven by species’ relative abundances, which defines probabilities of species encountering one another (2024). In contrast to theoretical predictions, we found that O‘ahu networks were nonrandom and had highly complex structures at local (site-specific) and regional (island-wide) scales. The regional network had low connectance, moderate specialization, and nested and modular topologies, with three distinct modules (Fig. 1, fig. S2, and table S4). At the local scale, networks had low to intermediate connectance and, unlike the regional network, were not nested. Similar to the regional network, six of seven local networks were specialized and modular, presenting three or four modules (Fig. 2, fig. S3, and table S4). We found that despite all interactions being novel and primarily involving introduced species, networks were structurally complex and notably similar between scales (local versus regional) and across sites. Furthermore, partner sharing (how distinct species share resources) in SDNs on Oʻahu is structured in a complementary way among bird and plant species, giving rise to distinct modules in which certain birds and plants interact preferentially. The emergence of such structures indicates that these novel SDNs largely reproduce the well-known patterns exhibited in mutualistic networks (18) and that SDN structure is highly conserved, regardless of variation in plant and bird communities. Given the low generalization in our novel insular networks, interaction release (22) either is not occurring or may occur in the form of consumption of more food types (e.g., insects, fruits, and nectar), rather than increased diversity within a specific resource type (e.g., greater number of species of fruits).

 

HawaiiSeeds_2

Fig. 2 Local SDNs on Oʻahu. Each of the seven networks includes all birds (left side of network) and plants (right side of network) consumed in a specific site. Blue, orange, green, and yellow depict modules with species interacting more among themselves than with other species, as identified by Beckett’s algorithm (33). All local networks but MTK were modular, presenting three or four modules. Line thickness indicates frequency of interactions. For a list of interacting species, see fig. S3. EKA, ʻĒkahanui; KAH, Kahanahāiki; MOA, Moanalua; MTK, Mount Kaʻala; PAH, Pahole; TAN, Tantalus; WAI, Waimea Valley.

Question 

What types of networks are found at each local site in O’ahu?

Methods

For each location, the birds are listed on the left column and plants are listed on the right column. Similar to Figure 1, a seed dispersal event is represented by a line between a bird and a plant. The thicker the line, the more frequent the connected bird and plant interacted. Each module has a distinct color.

Results A

The modules were determined using an algorithm that can identify groups in a network that interact more with each other than other participants of the network. 

Results B

Six out of seven O’ahu sites had modular networks.

Several studies suggest that the phylogenetic relationships of species contribute to structuring mutualistic networks [e.g., (2526)], which is an expected consequence of coevolved interactions among species interacting for prolonged (evolutionary) periods of time (8). Here we show that the interaction patterns recurrently identified in native-dominated networks also emerge in novel mutualistic networks composed of species with little or no shared evolutionary history. This result indicates that prolonged shared evolutionary history is not necessary for the emergence of complex network structure. We should note, however, that preexisting adaptations of introduced birds for frugivory and fruits for bird dispersal are necessary for their integration into novel networks. Furthermore, the presence of nested structure at regional, but not local, scale indicates the critical importance of spatial scale to understanding network patterns and their underlying processes. The wider variety of partners used at the larger scale (regional network) corresponds to the “fundamental niche,” whereas the subset of partners found at local scales indicates that local populations have much more restricted “realized niches” (2728). Therefore, not all species use available resources in the same way across all sites. By sampling across large spatial scales, researchers may be evaluating species’ fundamental niches and not population-level realized niches. Therefore, processes operating at different spatial scales may be overlooked or confounded (2728).

Most networks have been studied primarily as static entities at single sites, despite the importance of multiscale studies for understanding the processes underlying network structure and for evaluating the generalizability of network patterns (29). To examine interaction dynamics across sites and to test their association with environmental variables, we calculated the dissimilarity (interaction turnover) between pairs of networks, using data limited to species present in the networks. Highest dissimilarity occurs when two sites share no interactions. We decomposed this metric into two components: species turnover (βST—the proportion of interactions that are not shared owing to differences in species composition between two networks) and linkage turnover [βOS, also called rewiring—the proportion of interactions unique to a single network despite the occurrence of both partners in both networks (30)].

We found high interaction dissimilarity among sites owing to both changes in species composition and rewiring. This suggests high flexibility of birds and plants to switch partners, which is a major characteristic of highly successful invasive species (31). Interaction turnover across sites was high [interaction dissimilarity (βWN) = 0.57 ± 0.11, mean ± SE; n = 21 pairwise sites; Fig. 3 and table S5], indicating that, on average, only 43% of interactions were shared between sites despite the most common bird and plant species occurring at all sites (tables S2 and S3). Surprisingly, only 53% of the interaction dissimilarity was due to differences in species composition among sites (βST = 0.30 ± 0.09), whereas 47% was because pairs of species that interacted in one site did not interact in another site where they co-occurred (βOS = 0.27 ± 0.07; fig. S4). This indicates that, in addition to its influence on the structure of mutualistic networks [i.e., nestedness; (32)], partner switching is a major component of the spatial dynamics of novel networks. High interaction dissimilarity has also been reported in specialized, native-dominated pollination networks, even between spatially close networks (33). Thus, plant-animal networks appear to have distinct links (high interaction rewiring) even when the same species are present in both sites, irrespective of whether networks are dominated by native or introduced species.

Figure 3

Fig. 3 Interaction dissimilarity between each pair of sites on Oʻahu. Interaction dissimilarity (βWN) was decomposed into its two components: species turnover (βST) and linkage turnover among species shared by pairwise sites (i.e., rewiring, βOS).

Question 

Are there similar plant-bird interactions occurring between different sites of O’ahu? 

Methods

The overall interaction dissimilarity (βWN), which measures the percentage of interactions two sites do not share in common, was calculated using the Whittaker’s equation. This equation essentially considers the number of interactions shared between two sites (a) and the number of interactions unique to each individual site (b and c). A βWN value can range from zero (meaning that both sites completely share interaction patterns) to 1 (1 meaning the sites do not share any interactions). 

βWN
Results βOS 

The graph depicts the rewiring as a black bar, representing the portion of the interaction dissimilarity which was due to different interactions between the two sites despite both areas having the same species of plants and animals. 

Results βST

Species turnover (shown as the grey bar) represents the portion of the interaction dissimilarity that is caused by two sites having different species of plants and animals, and so as a consequence, different interactions.

Abiotic factors had a greater effect than biotic factors on the overall interaction dissimilarity and the dissimilarity caused by species turnover between sites, whereas interaction rewiring was not influenced by any factor examined (tables S6 to S11). Specifically, interaction dissimilarity and the dissimilarity caused by species turnover were influenced by elevation and rainfall, but not by percent of introduced plant species (tables S6 to S9). This suggests that the environment indirectly influences interactions via effects on species distributions, including the distribution of introduced species. However, the lack of association between rewiring and examined factors indicates that birds and plants in the system are highly flexible and can switch partners, irrespective of abiotic conditions and the identity of species in the community.

Lastly, we compared O‘ahu SDNs with native-dominated SDNs around the world and found that O‘ahu’s novel networks resemble the structure and stability of native-dominated networks. We assembled and analyzed a dataset of 42 avian SDNs encompassing a broad geographical range, with data from islands (n = 17) and continents (n = 25) in tropical (n = 18) and nontropical (n = 24) areas (table S12). Although some of the other SDNs in the analyses included introduced species [e.g., (734)], SDNs on O‘ahu present an extreme case of dominance by introduced species (>50%), coupled with extinction of all native frugivorous birds. For these 42 networks and the seven on O‘ahu, we calculated a set of weighted (for 26 networks where frequency of interaction was reported) and binary (for all 42 networks) descriptors of network structure. We also estimated robustness (stability to species loss) of each network as the rate of secondary extinction expected under the simulated loss of network partners, assuming a species goes extinct when all connected partners are lost (3536). We estimated robustness of animals to the extirpation of plants (assuming bottom-up control) and robustness of plants to the extirpation of animals (top-down control). We simulated two scenarios, one in which order of extirpation was random and another—more extreme—scenario in which order was from the most generalist to the most specialist species. After using a null model correction on each metric to account for variation in sampling intensity and network dimensions across studies (14), we compared the 95% confidence intervals for the O‘ahu networks with the global dataset. We found that specialization, modularity, nestedness, and the simulated robustness in all scenarios to species loss of the O‘ahu networks overlapped with the range of values observed in other networks. These results held true for both weighted and binary data and when O‘ahu’s networks were compared to subsets of networks from tropical and nontropical islands and continents. The only exceptions were that the specialization and weighted modularity observed in O‘ahu networks were lower than those in networks from nontropical continental areas (Fig. 4 and table S13).

 

Figure 4 comparison of structure and stability

Fig. 4 Structure and stability of 42 SDNs from islands and continents in tropical and nontropical communities worldwide in comparison to novel networks on Oʻahu. Significant difference (*) occurs when the 95% confidence interval of a metric for the seven sites in Oʻahu (gray shaded area) does not overlap the intervals for non-Oʻahu networks (colored bars). H2′, complementary specialization; wNODF and NODF, nestedness; Qw and Qb, modularity; ranP, ranA, degP, and degA, network robustness to the sequential extinction of animals (A) and plants (P) by random (ran) or from the most-generalist to the most-specialist species (deg). The latter is calculated only for binary data.

Question 

How does the structure and stability of seed dispersal networks in O’ahu compare to seed dispersal networks in other habitats around the world?

Methods

Data of seed dispersal networks from different habitats were calculated into weighted (top) or binary (bottom) networks. On the y- axis are measurements of different metrics of network interactions (such as nestedness, specialization, and modularity) and the x -axis organizes the different habitats into islands, continents, or both.

Confidence Interval 

Each graph has a gray bar representing the 95% confidence interval for data collected from O’ahu. This bar depicts the range of values estimated to include the true value for the network metric measured in O’ahu. The colored bars represent the range of values measured from the non-O’ahu habitats. When the colored bar doesn’t overlap the gray bar, that means there is a significant difference between the measurements taken from O’ahu and non-O’ahu locations.

Results

The colored bars overlap the gray bars for all of the different network measurements except for modularity (Qw) and nestedness (wNODF) under weighted networks. The seed dispersal networks of O’ahu share the majority of network structures of other habitats around the world.

Most SDNs from communities around the world have been described as specialized, nested, and modular [e.g., (19)], and the variation in such structures reveals the responses of species interactions to biotic and abiotic factors at both small (37) and large scales (141534). Here we show that O‘ahu’s novel networks notably resemble the structure and stability of native-dominated networks elsewhere. This high degree of similarity between novel and native-dominated networks suggests that the processes that structure interactions in such communities are largely independent of species identity and that ecological filtering occurs over relatively short (ecological) time, leading to functionally similar sets of players as compared with systems that have long evolutionary histories. Yet, because filtering depends on the pool of species introduced, novel networks may have an incomplete set of roles fulfilled. For example, in Hawai‘i, large frugivorous birds are absent, resulting in a lack of dispersal of large native fruits (38). Therefore, functional characteristics (e.g., beak, seed, and fruit sizes) and species abundance (39) may be more important in the structure of mutualistic networks than species identity, supporting the role of ecological fitting (40). Thus, further investigation on the influence of functional traits and abundances on novel networks may shed light on the ultimate mechanisms driving network structure and species roles.

By studying novel networks across scales and comparing them with native-dominated networks worldwide, we identify several key considerations. First, sampling across scales is critical for testing generalizability of patterns and identifying the underlying processes (e.g., abiotic or biotic) structuring networks. Thus, explicitly examining multiple spatial scales is an essential next step toward advancing the understanding of processes that define specialization and shape ecological networks (41). We also predict that the patterns described here are more likely to be found in other isolated ecosystems, such as oceanic islands or isolated habitat patches, which are more prone to species invasions than less-isolated ecosystems. Second, our results show that introduced dispersers incompletely fulfill species roles lost by O‘ahu’s extreme scenario of plant and bird loss and introductions. Although these introduced birds on O‘ahu are the only dispersers of native plants, they disperse a much higher proportion of seeds from invasive plants; therefore, their presence is a “double-edged sword” for conservation. The flexibility of birds and plants for partner switching and the fact that novel networks may be highly robust to species removal should be considered in restoration efforts. These efforts would benefit from initiatives that increase use of restoration sites by targeted frugivores and their consumption of native fruits. This would include outplanting commonly consumed native plants (e.g., Pipturus albidus) within plant restoration areas, removing commonly consumed introduced plants in sites with high densities of native fruits, and attracting (e.g., via playback) specific frugivores to restoration sites. The dramatic changes that have occurred in Hawaiian ecosystems provide an opportunity to better understand, anticipate, and mitigate the impacts of widespread and increasing biological invasions and species extinction, while also determining how network complexity develops.

 

Supplementary Materials

www.sciencemag.org/content/364/6435/78/suppl/DC1

Materials and Methods

Figs. S1 to S5

Tables S1 to S14

References (4382)

http://www.sciencemag.org/about/science-licenses-journal-article-reuse

This is an article distributed under the terms of the Science Journals Default License.

 

References and Notes

1. S. L. Lewis, M. A. Maslin, Nature 519, 171–180 (2015).

2. C. Hui, D. M. Richardson, Invasion Dynamics (Oxford Univ. Press, 2017).

3. J. F. Brodie et al., Trends Ecol. Evol. 29, 664–672 (2014).

4. A. Traveset, D. M. Richardson, Annu. Rev. Ecol. Evol. Syst. 45, 89–113 (2014).

5. M. J. O. Pocock, D. M. Evans, J. Memmott, Science 335, 973–977 (2012).

6. C. E. Mitchell et al., Ecol. Lett. 9, 726–740 (2006).

7. R. H. Heleno, J. A. Ramos, J. Memmott, Biol. Invasions 15, 1143–1154 (2013).

8. J. N. Thompson, The Geographic Mosaic of Coevolution (Univ. of Chicago Press, 2005).

9. J. Bascompte, P. Jordano, J. M. Olesen, Science 312, 431–433 (2006).

10. J. T. Foster, S. K. Robinson, Conserv. Biol. 21, 1248–1257 (2007).

11. C. G. Chimera, D. R. Drake, Biotropica 42, 493–502 (2010).

12. P. J. Conry, R. Cannarella, “Hawaii statewide assessment of forest conditions and trends: 2010” (Hawaii Department of Land and Natural Resources, Division of Forestry and Wildlife, 2010).

13. Materials and methods are available as supplementary materials.

14. B. Dalsgaard et al., Ecography 40, 1395–1401 (2017).

15. M. Schleuning et al., Curr. Biol. 22, 1925–1931 (2012).

16. E. Sebastián-González, B. Dalsgaard, B. Sandel, P. R. Guimarães Jr., Glob. Ecol. Biogeogr. 24, 293–303 (2015).

17. J. Bascompte, P. Jordano, Annu. Rev. Ecol. Evol. Syst. 38, 567–593 (2007).

18. D. P. Vázquez, N. Blüthgen, L. Cagnolo, N. P. Chacoff, Ann. Bot. 103, 1445–1457 (2009).

19. A. de Almeida, S. B. Mikich, Oikos 127, 184–197 (2018).

20. J. Vizentin-Bugoni, P. K. Maruyama, C. S. de Souza, F. Ollerton, A. R. Rech, M. Sazima, in Ecological Networks in the Tropics, W. Dáttilo, V. Rico-Gray, Eds. (Springer, 2018), pp. 73–91.

21. R. H. MacArthur, J. M. Diamond, J. R. Karr, Ecology 53, 330–342 (1972).

22. A. Traveset et al., Nat. Commun. 6, 6376 (2015).

23. A. M. Martín González et al., Glob. Ecol. Biogeogr. 24, 1212–1224 (2015).

24. A. Krishna, J. P. R. Guimaraes Jr., P. Jordano, J. Bascompte, Oikos 117, 1609–1618 (2008).

25. E. L. Rezende, J. E. Lavabre, P. R. Guimarães, P. Jordano, J. Bascompte, Nature 448, 925–928 (2007).

26. R. S. Vitória, J. Vizentin-Bugoni, L. D. S. Duarte, Oikos 127, 561–569 (2018).

27. N. Blüthgen, F. Menzel, N. Blüthgen, BMC Ecol. 6, 9 (2006).

28. V. Devictor et al., J. Appl. Ecol. 47, 15–25 (2010).

29. W. Dáttilo, V. Rico-Gray, Eds., Ecological Networks in the Tropics (Springer, 2018).

30. T. Poisot, E. Canard, D. Mouillot, N. Mouquet, D. Gravel, Ecol. Lett. 15, 1353–1361 (2012).

31. H. A. Mooney, E. E. Cleland, Proc. Natl. Acad. Sci. U.S.A. 98, 5446–5451 (2001).

32. F. Zhang, C. Hui, J. S. Terblanche, Ecol. Lett. 14, 797–803 (2011).

33. D. W. Carstensen, M. Sabatino, K. Trøjelsgaard, L. P. C. Morellato, PLOS ONE 9, e112903 (2014).

34. M. Nogales et al., Glob. Ecol. Biogeogr. 25, 912–922 (2016).

35. J. Memmott, N. M. Waser, M. V. Price, Proc. R. Soc. London Ser. B 271, 2605–2611 (2004).

36. E. Burgos et al., J. Theor. Biol. 249, 307–313 (2007).

37. C. N. Kaiser-Bunbury et al., Nature 542, 223–227 (2017).

38. S. Culliney, L. Pejchar, R. Switzer, V. Ruiz-Gutierrez, Ecol. Appl. 22, 1718–1732 (2012).

39. A. González-Castro, S. Yang, M. Nogales, T. A. Carlo, AoB Plants 7, 1–10 (2015).

40. D. H. Janzen, Oikos 45, 308–310 (1985).

41. C. F. Dormann, J. Fründ, H. M. Schaefer, Annu. Rev. Ecol. Evol. Syst. 48, 559–584 (2017).

42. J. Vizentin-Bugoni et al., Dataset from: Structure, spatial dynamics, and stability of novel seed dispersal mutualistic networks in Hawai’i, Version 1, Harvard Dataverse (2019); https://doi.org/10.7910/DVN/CMRVAK.