# Influence of taxonomic resolution on mutualistic network properties

Estelle Renaud, Emmanuelle Baudry, Carmen Bessa-Gomes

► To cite this version:

Estelle Renaud, Emmanuelle Baudry, Carmen Bessa-Gomes. Influence of taxonomic resolution on mutualistic network properties. Ecology and Evolution, 2020, 10 (7), pp.3248-3259. 10.1002/ece3.6060. hal-04241364

HAL Id: hal-04241364

https://agroparistech.hal.science/hal-04241364v1

Submitted on 13 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Influence of taxonomic resolution on mutualistic network properties

Estelle Renaud ID |

Emmanuelle Baudry |

Carmen Bessa-Gomes

Ecologie Systématique Evolution, CNRS,
AgroParisTech, Université Paris-Saclay,
Orsay, France

**Correspondence**

Estelle Renaud, Université Paris-Saclay, 362 rue du Doyen André Guinier, 91400 Orsay, France.

Email: estelle.renaud@ens-lyon.fr

**Funding information**

ENS de Lyon

## Abstract

1. Ecologists are increasingly interested in plant-pollinator networks that synthesize in a single object the species and the interactions linking them within their ecological context. Numerous indices have been developed to describe the structural properties and resilience of these networks, but currently, these indices are calculated for a network resolved to the species level, thus preventing the full exploitation of numerous datasets with a lower taxonomic resolution. Here, we used datasets from the literature to study whether taxonomic resolution has an impact on the properties of plant-pollinator networks.

2. For a set of 41 plant-pollinator networks from the literature, we calculated nine network index values at three different taxonomic resolutions: species, genus, and family. We used nine common indices assessing the structural properties or resilience of networks: nestedness (estimated using the nestedness index based on overlap and decreasing fill [NODF], weighted NODF, discrepancy [BR], and spectral radius [SR]), connectance, modularity, robustness to species loss, motifs frequencies, and normalized degree.

3. We observed that modifying the taxonomic resolution of these networks significantly changes the absolute values of the indices that describe their properties, except for the spectral radius and robustness. After the standardization of indices measuring nestedness with the Z-score, three indices—NODF, BR, and SR for binary matrices—are not significantly different at different taxonomic resolutions. Finally, the relative values of all indices are strongly conserved at different taxonomic resolutions.

4. We conclude that it is possible to meaningfully estimate the properties of plant-pollinator interaction networks with a taxonomic resolution lower than the species level. We would advise using either the SR or robustness on untransformed data, or the NODF, discrepancy, or SR (for weighted networks only) on Z-scores. Additionally, connectance and modularity can be compared between low taxonomic resolution networks using the rank instead of the absolute values.

## KEYWORDS

connectance, modularity, nestedness, networks, plant-pollinator, robustness, taxonomic resolution

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

© 2020 The Authors. *Ecology and Evolution* published by John Wiley & Sons Ltd.
# 1 | INTRODUCTION

The study of species interactions has always been central in ecology. Such interactions have historically been examined by focusing on two interacting species, but in recent years, the marked increase in the amount of biological information available and the development of novel approaches and tools have placed a new focus on the study of interaction networks (Proulx, Promislow, & Phillips, 2005). Ecological networks may provide important insights that cannot be gained when species are studied in isolation. They currently play a central role in key aspects of ecological theory such as the long-standing question of the relationship between complexity and stability in ecosystems (Montoya, Pimm, & Solé, 2006; Thébault & Fontaine, 2010) or the interplay between interspecific competition and ecological niche (Bastolla et al., 2009). Ecological networks are also powerful tools for applied ecology, as they can be used to monitor the impact of biological perturbations on an ecosystem or the efficiency of restoration programs (Kaiser-Bunbury & Blüthgen, 2015; Kaiser-Bunbury et al., 2017; Memmott, 2009).

Most studies on ecological networks have focused on three main categories of networks defined according to the type of species and their interactions: food webs, parasitoid host interaction networks, and more recently, mutualist interaction networks (Ings et al., 2009). In this paper, we concentrate on the case of the mutualistic networks linking plants and pollinators, which have attracted particular attention in recent years. Indeed, pollinators have an essential ecological function, namely the pollination function, which is threatened in many parts of the world by the sharp decline in pollinators on account of the many threats that they face (Goulson, Nicholls, Botías, & Rotheray, 2015). Such a decline in pollinator populations may harm both wild biodiversity and agricultural productivity (Garibaldi et al., 2013).

The use of a network makes it possible to synthesize in a single object the species and interactions linking them and thus constitute the community of species (Delmas et al., 2019). It thus becomes possible to use the many methods developed to study ecological networks to describe their structure and properties using different indices (Lau, Borrett, Baiser, Gotelli, & Ellison, 2017). One structural characteristic that has received particular attention in the study of plant-pollinator networks is nestedness (Bascompte, Jordano, Melián, & Olesen, 2003; Table 1). A nested network is characterized by the extent to which interactions of less-connected species form subsets of the interactions of more-connected species. Other frequently examined structural characteristics of mutualistic networks are connectance, the proportion of realized interactions among all possible ones, and modularity, that is, the extent to which linked interactions between pollinators and plants are organized into delimited modules, as well as motifs, which are subnetworks representing the interactions between a given number of taxa (Milo et al., 2002). These properties have been associated with the ecosystem's resilience to perturbations (Soares, Ferreira, & Lopes, 2017). It has, for example, been shown that high levels of connectivity, modularity, and nestedness promote both the structural and dynamic stability of mutualist interaction networks (Vanbergen, Woodcock, Heard, & Chapman, 2017).

A large number of datasets on plants and their pollinators have been collected to date. However, given the large number of pollinator species potentially present in a community, as well as the relative difficulty in identifying some of these pollinators at the species level, a significant portion of the collected datasets has a taxonomic resolution lower than the species level. For a given research effort, there is therefore a trade-off between the quantity of possible identifications and the taxonomic accuracy of these identifications, which makes it difficult to produce large or numerous sets of data identified down to the species level. An extreme point in this regard is the datasets provided by citizen science programs for pollinators (Toomey & Domroese, 2013) such as the Spipoll program in France (Deguines, Julliard, Flores, & Fontaine, 2012), which generally allow very large datasets to be collected, although their taxonomic accuracy does not generally extend to the species level (Dickinson, Zuckerberg, & Bonter, 2010; Kremen, Ullman, & Thorp, 2011).

Currently, network analyses are performed on networks with varying levels of taxonomic precision, which makes comparisons between studies or even sites of the same studies potentially invalid, because we do not know how taxonomic resolution influences the indices of those networks, nor how they should be interpreted. If possible, it would, however, be interesting to use network analyses on such datasets in order to fully exploit the information contained therein and allow comparisons with other studies. Here, we sought to establish whether taxonomic resolution has an influence on the architecture and properties of a mutualistic network estimated using several indices. We used a set of 41 plant-pollinator networks from the literature and compared their index values at three different taxonomic resolutions: species, genus, and family. We showed that for a given network, changing the taxonomic resolution usually significantly changes the value of most indices. We also show that after the standardization (with the Z-score, using null models) of the indices measuring nestedness, these three indices are no longer differed significantly at different taxonomic resolutions. We also used another normalization measure for one nestedness index (NODF) called NODFc and show that this measure is robust to a lower taxonomic resolution (Song, Rohr, & Saavedra, 2017). Additionally, we showed that among the set of 41 networks, the relative value of a given network for a given index is well conserved across different taxonomic resolutions, particularly between the species and genus levels.

# 2 | MATERIALS AND METHODS

## 2.1 | Overview

We used plant-pollinator networks from the literature (Vázquez, Goldberg, & Naik, 2003) determined to the level of species. For each species-level network, we deduced the equivalent network at the genus and family levels. We then calculated several indices commonly used to estimate mutualistic network properties for each of these networks and then compared their values across taxonomic
resolutions. Data manipulation and analysis were conducted with the R language (R version 3.2.3, 2015-12-10). The script used for those results is accessible here: https://gitlab.com/EstelleRenaud/taxonomic_influence_network_properties

## 2.2 | Network indices

We selected frequently used indices that describe various properties of interaction networks, namely nestedness, connectance, modularity, motifs, and robustness. Given the particular importance

of generalist pollinator species in maintaining plant-pollinator networks (Martín González, Dalsgaard, & Olesen, 2010), we also added one index calculated at the species level, that is, the normalized degree.

The characteristics of these indices are summarized in Table 1.

To estimate nestedness, no unique index has been established to date as consensual, which led us to use four indices: nestedness index based on overlap and decreasing fill (NODF; Almeida-Neto, Guimarães, Guimarães, Loyola, & Ulrich, 2008), spectral radius (SR; Staniczenko, Kopp, & Allesina, 2013), discrepancy (BR; Brualdi & Sanderson, 1999), and NODF for weighted matrices (WNODF;

<table><caption>TABLE 1 Commonly used network indices</caption><thead><tr><th>Index</th><th>Matrix<sup>a</sup></th><th>Nature of the index, per network</th><th>Representation</th></tr></thead><tbody><tr><td>Nestedness</td><td></td><td></td><td></td></tr><tr><td>NODF</td><td>B</td><td>Extent to which interactions of less-connected species form subsets of the interactions of more-connected species</td><td></td></tr><tr><td>BR</td><td>B</td><td></td><td></td></tr><tr><td>SR</td><td>B</td><td></td><td></td></tr><tr><td>WNODF</td><td>W</td><td></td><td></td></tr><tr><td>Connectance</td><td>B</td><td>Proportion of realized interactions among all possible ones</td><td>High Low</td></tr><tr><td>Modularity</td><td>W</td><td>Extent to which interactions between pollinators and plants are organized into delimited modules</td><td></td></tr><tr><td>Robustness</td><td>B</td><td>Speed at which plant taxa disappear as pollinator taxa disappear</td><td></td></tr><tr><td>Normalized degree</td><td>B</td><td>Connectance of each taxa (this is the only index calculated per taxa and not per network)</td><td></td></tr><tr><td>Motif frequency</td><td>B</td><td>Frequency of each of the 17 kinds of motifs that can link up to 5 taxa between them</td><td>1 3 6 15</td></tr></tbody></table>

<figure><img src="image_2.png" /></figure>

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

Abbreviations: BR, discrepancy; NODF, nestedness index based on overlap and decreasing fill; SR_Bin, spectral radius calculated on binary (absence/ presence) matrices; SR_Qua, spectral radius calculated on weighted (abundance) matrices; WNODF, NODF calculated on weighted matrices.

<footnote><sup>a</sup>B and W indicate an index calculated on binary and weighted matrices, respectively.</footnote>
Almeida-Neto & Ulrich, 2011). NODF and BR are indices for binary matrices, while WNODF is adapted to weighted matrices and SR can be used for both. The values for NODF, WNODF, BR, and SR were calculated using the *Falcon* package (Beckett, Boulton, & Williams, 2014) for R. Furthermore, BR was calculated using the method presented in Brualdi and Sanderson (Brualdi & Sanderson, 1999): it is the minimal number of differences with a perfectly nested matrix with the same size, number of links, and column (or row) sums as the real network. The SR of a network is thus the largest of its matrix eigenvalues (Staniczenko et al., 2013).

We used five additional indices. Four of them—connectance, robustness, motifs, and normalized degree—are calculated on the presence/absence matrices, whereas modularity is calculated on frequency matrices. Network connectance was calculated as the sum of links divided by the number of cells in the matrix. Network modularity was measured according to the Beckett algorithm DIRTLPAwb+ (Beckett, 2016), which aims to estimate the modularity of the network using three steps. The first uses label propagation to obtain a locally maximized modularity (bottom-up); the second agglomerates the modules found in the first step if it allows for an increased modularity; the third repeats these steps until modularity can no longer be increased. DIRTLPAwb+ then randomizes the initial labeling of nodules multiple times and returns the result with the greatest modularity score. Modularity itself was then calculated as the modularity M proposed by Newman (Newman, 2006). Following Burgos et al. (2007), network robustness was measured as the area under the attack tolerance curve, defined as the speed at which plant taxa disappear as pollinator taxa disappear. Basically, as pollinator taxa disappear, plants that rely exclusively on them (according to the network) also disappear, thus creating a curve in which the percentage of remaining plants depends on the percentage of remaining pollinators. This method assumes that preferences are static, that is, that plants that rely on one given pollinator taxon will not be able to switch to another pollinator if it disappears. Normalized degree is calculated for each taxon in a network as the sum of the links of that taxon scaled by its number of potential partners. Because the normalized degree index produces one value per taxon, for each matrix, we chose to characterize matrices by their quartile values in order to accommodate differently skewed distributions of the normalized degrees between networks. The R package *bipartite* (Dormann, Gruber, & Fründ, 2008) was used to calculate the normalized degree, robustness, modularity, and connectance of the network. Motifs were compared using the frequency of the 17 different motifs that involve up to five different taxa, which we calculated using the *mcount* function of the *bmotif* package (Simmons et al., 2019). The frequency of each motif was calculated as the number of times a given motif occurs in the network, divided by the number of motifs of the same size (involving the same number of taxa) that occurs in the network.

Additionally, because most nestedness indices are known to be sensitive to the size (number of rows and columns) and fill (number of nonzeros) of the input matrix (Rodríguez-Gironés & Santamaría,

2006)—two properties that are modified when the taxonomic resolution is changed—we performed standardization with Z-scores for nestedness indices (Ulrich & Almeida-Neto, 2012). These are calculated as the difference between the observed index value and the value expected under a null model divided by the standard deviation under this null model. Z-scores were obtained by calculating 500 null models of each matrix and comparing the resulting mean value to that calculated for the matrix. As the null models take into account the size of the matrix, this minimizes the possible effect of size on the index values. We used the *Falcon* package to calculate the Z-scores. For the binary indices (BR, NODF, SR), we followed Bascompte et al. (2003) and used null models obtained from 500 iterations of the DD (degreeprobable-degreeprobable) model. This model is intermediate in terms of constraints on row and column totals (part of the class termed “PP,” proportional to both row totals and column totals; Strona, Ulrich, & Gotelli, 2018), with one extreme being the fixed-fixed model that is susceptible to type II errors (Gotelli, 2000) and the other the equiprobable-equiprobable model that is susceptible to type I error (Wright, Patterson, Mikkelson, Cutler, & Atmar, 1997). This model has statistically determined elements following the degree distribution of the initial matrix as $p_{ij} = 1/2*(d_j/r + k_j/c)$, where $p_{ij}$ is the probability of assigning a 1 to the ith row and jth column, $d_j$ is the column degree of the jth column, $k_j$ is the row degree of the ith row, and r and c are the respective numbers of rows and columns. For the weighted indices (WNODF, SR), we used two kinds of null models, as no null model has been established as more suited to WNODF or SR yet: The first set of null matrices is obtained from 500 iterations of the row and column total average model (introduced in the Falcon software) that averages two matrices: a matrix created conserving the row totals and redistributing a random portion of that total to each element of a given row, and a matrix following the same principle with the column totals. The second kind is Patefield's historical r2dtable model, implemented with the null model function (option “r2dtable”) of the bipartite R package (Dormann et al., 2008). We also generated 500 matrices under that model.

Finally, because Z-scores for NODF have been criticized for their sensitivity to connectance and number of taxa (Song et al., 2017), we used the normalization proposed by Song et al. as $NODF_c = NODF_n/(C*\log(S))$, where C is the connectance, S is the geometric mean of plants and pollinators in the network, and $NODF_n = NODF/\max(NODF)$, where $\max(NODF)$ is the maximal NODF value that could be attained in a network with the same number of rows, columns, and links as the original network. Max(NODF) was calculated using the maxnodf R package (Hoeppke, 2019).

## 2.3 | Pollination networks

We extracted all plant-pollinator interaction networks from the Interaction Web Database (Vázquez et al., 2003). All networks
were issued from previously published data (Table S1). We only kept matrices for which taxa determination was possible using the taxize package; that is, valid taxonomic names resolved at the genus or species level. In some cases, we replaced old taxonomic names by a current valid synonym. We also only kept matrices that dealt with several families, which left us with a dataset of 41 matrices, 10 of which were binary (presence/absence) matrices. The remaining 31 were weighted according to the frequency of the visitation or a proxy of that frequency. The number of taxa in the matrices varied from seven to 135 for plants, and 12 to 144 for pollinators.

We then used the taxize package (Chamberlain & Szöcs, 2013; version 0.9.0) from R to extract from the taxonomic information supplied by the authors the taxonomic affiliation from the superior ranks. Only the identification from the species, genus, family, and order ranks was retained, as these were the ranks most often known for all observations. The database GBIF (GBIF, 2018) was used as a reference.

We transformed each of the 41 previously described matrices into interaction matrices determined at the species level by keeping only the observations (within each network) for which both the plant and pollinator were determined to the species level. From these species-level matrices, we deduced the genus-level and then the family-level matrices.

## 2.4 | Statistical tests

To examine the influence of the taxonomic level on the structure of a given matrix, we compared the values of the indices for Species-level matrices and Genus-level matrices, Genus-level matrices and Family-level matrices, and Species-level matrices and Family-level matrices, using a one-way analysis of variance. Post hoc tests were performed with a Bonferroni correction, using the built-in pairwise.t. test R function, with the "paired" option. We also performed the same analyses after standardizing nestedness values using Z-scores.

To investigate whether an index was useful for comparing different observed matrices, we performed a nonparametric correlation test (cor.test on R) to calculate both the value and significance of Spearman's rho for a given index in Species-level matrices, Genus-level matrices, and Family-level matrices. This allowed us to test whether the relative ranks of this index's values were significantly correlated between one taxonomical level and another.

## 3 | RESULTS

For each of the species-, genus-, and family-determined matrices, we obtained a set of values per index for the species-determined matrices, and another set for genus- and family-determined matrices. We compared these using a one-way analysis of variance. Most indices show a significant effect of taxonomic level on their value (Table 2 and Figure 1; results for normalized degrees and motifs frequencies are presented in the Appendix A). Robustness and SR (both on binary and on weighted matrices) are the only indices without

<table><caption>**TABLE 2** Results of one-way analysis of variance comparing index values at three different taxonomic resolutions</caption><thead><tr><th>Indices</th><th>df</th><th>F</th><th>p-Value</th><th>Post hoc test results</th></tr></thead><tbody><tr><td>Nestedness</td><td></td><td></td><td></td><td></td></tr><tr><td>NODF</td><td>2, 120</td><td>15.03</td><td>1.496e-05</td><td>All levels differ significantly</td></tr><tr><td>BR</td><td>2, 120</td><td>4.2035</td><td>.0172</td><td>All levels differ significantly</td></tr><tr><td>SR_Bin</td><td>2, 120</td><td>1.1573</td><td>.3178</td><td></td></tr><tr><td>SR_Qua</td><td>2, 90</td><td>0.5765</td><td>.5639</td><td></td></tr><tr><td>WNODF</td><td>2, 90</td><td>11.626</td><td>3.229e-05</td><td>All levels differ significantly</td></tr><tr><td>Connectance</td><td>2, 120</td><td>20.237</td><td>2.67e-08</td><td>All levels differ significantly</td></tr><tr><td>Modularity</td><td>2, 90</td><td>5.2213</td><td>.007155</td><td>All levels differ significantly</td></tr><tr><td>Robustness</td><td>2, 120</td><td>0.0831</td><td>.9204</td><td></td></tr></tbody></table>

Note: NODF, BR, SR_Bin, robustness, and connectance are calculated on presence/absence networks. SR_Qua, WNODF, and modularity are calculated on abundance-based networks. The results for the normalized degree are presented in the Appendix A.

Abbreviations: BR, discrepancy; NODF, nestedness index based on overlap and decreasing fill; SR_Bin, spectral radius calculated on binary (absence/presence) matrices; SR_Qua, spectral radius calculated on weighted (abundance) matrices; WNODF, NODF calculated on weighted matrices.

any significant influence of taxonomic level on their values. This indicates that two matrices cannot be directly compared with most indices if they are not at the same level of taxonomic resolution. The same is true with family versus genus comparisons, and family versus species comparisons.

We used the Z-score (with two kinds of null models for the weighted indices) to take into account the difference in the matrix fills and sizes caused by the change in taxonomic resolutions, as well as another normalization by the maximal NODF (noted as NODF<sub>c</sub>). We compared nestedness Z-scores and NODF<sub>c</sub> values from one taxonomic level to another using a one-way analysis of variance test. After this standardization, only the WNODF showed a significant effect of taxonomic level on its Z-score value (Figure 2 and Table 3), with a higher level of significance using the r2dtable null model than the RTCA. NODF<sub>c</sub> showed no influence of the taxonomic level on its values.

While most untransformed indices for species-, genus-, and family-determined matrices are significantly different, their values seem very strongly correlated among them. Using the Spearman's rho calculation, we observed that all tested index ranks showed a high positive correlation (Figure 3; for the results of the normalized degrees and motifs frequencies, see the Appendix A), indicating that the ranks of indices are well conserved at different taxonomic resolutions. All indices show a rank correlation superior to .8 between species and genus ranks, except for one motif frequency and certain
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /></figure>

<figure><img src="image_14.png" /></figure>

**FIGURE 1** Untransformed index value distribution for species (S)-, genus (G)-, and family (F)-determined matrices. From left to right and top to bottom: NODF, BR, SR_Bin, Robustness, Connectance, WNODF, SR_Qua, Modularity. BR, discrepancy; NODF, nestedness index based on overlap and decreasing fill; SR_Bin, spectral radius calculated on binary (absence/presence) matrices; WNODF, NODF calculated on weighted matrices. 0.001 < *** < 0.01 < ** < 0.05; \*Next to the index name reflects the result of the ANOVA test, on the graph itself reflects the results of the post hoc paired t test, with a Bonferroni correction

quartile values of the normalized degree. This strong correlation becomes weaker for most indices as the taxonomic levels increase, particularly for comparisons with family versus species. Among the indices, the strongest options seem to be modularity and SR for nestedness, as they show a strong correlation (superior to .8) between ranks at the species, genus, and family levels for all possible comparisons. On the contrary, motif frequencies seem to be less reliable, as they show lower values for correlation between ranks than other indices, especially at the family versus species level, with most correlations being below .6 and some even below .35 (see the Appendix A for the full values).

4 | DISCUSSION

Using plant-pollinator interaction networks from the literature, we showed that modifying the taxonomic resolution of these networks significantly changes the absolute values of the indices that describe their properties, except for two indices, namely the SR (both for

binary and quantitative matrices) and robustness to species loss. If a standardization of the indices measuring nestedness is performed using the Z-score, then three indices—NODF, BR, and SR for both binary and weighted matrices—are not significantly different at different taxonomic resolutions. Finally, the ranks of all indices are strongly conserved at different taxonomic resolutions, particularly between the species and genus levels.

For nestedness, we observed for both NODF and BR that the absolute values of these indices, but not the associated Z-scores, are strongly modified by the change in taxonomic resolution. This result is in agreement with the work of Almeida-Neto et al. (2008) who showed that NODF and BR are markedly affected by the matrix fill, a parameter that is modified when the taxonomic resolution is changed. However, this effect is no longer present after standardization using the Z-score. Nevertheless, we observed that the absolute values of the SR for both binary and quantitative matrices are not modified by the change in taxonomic resolution, in accordance with the results of Strona and Fattorini (2014), which show an absence of relationship between the SR values and the filling of the matrices.
<figure><img src="image_15.png" /><figcaption>**FIGURE 2** Z-score distribution for nestedness indices for species-, genus-, and family-determined matrices. BR, discrepancy; NODF, nestedness index based on overlap and decreasing fill; NODFc, NODF normalized according to Song et al.'s method (see Materials and Methods); SR_Bin, spectral radius calculated on binary (absence/presence) matrices; SR_Qua, spectral radius calculated on weighted (abundance) matrices; WNODF, NODF calculated on weighted matrices. 0.001 < *** < 0.01 < ** < 0.05; \*Next to the index name reflects the result of the ANOVA test, on the graph itself reflects the results of the post hoc paired t test, with a Bonferroni correction. (RTCA) and (r2dtable) specify the results obtained through two different null models</figcaption></figure>

<table><caption>**TABLE 3** Results of one-way analysis of variance comparing Z-score values (as well as NODF normalized according to Song et al. (2017)) at three different taxonomic resolutions</caption><thead><tr><th>Indices</th><th>df</th><th>F</th><th>p-Value</th><th>Post hoc test results</th></tr></thead><tbody><tr><td>Nestedness</td><td></td><td></td><td></td><td></td></tr><tr><td>NODF</td><td>2, 120</td><td>0.5698</td><td>.5672</td><td></td></tr><tr><td>BR</td><td>2, 120</td><td>0.4745</td><td>.6233</td><td></td></tr><tr><td>SR_Bin</td><td>2, 120</td><td>2.6813</td><td>.07258</td><td></td></tr><tr><td>SR_Qua (RTCA)</td><td>2, 90</td><td>0.0729</td><td>.9297</td><td></td></tr><tr><td>SR_Qua (r2dtable)</td><td>2, 90</td><td>0.3565</td><td>.7002</td><td></td></tr><tr><td>WNODF (RTCA)</td><td>2, 90</td><td>3.0987</td><td>.04995</td><td>Genus and species do not differ significantly, while the two other comparisons are significantly different</td></tr><tr><td>WNODF (r2dtable)</td><td>2, 90</td><td>52.971</td><td>&lt;2.2e-16</td><td>All levels differ significantly</td></tr><tr><td>NODF<sub>c</sub></td><td>2, 104</td><td>0.9219</td><td>.4011</td><td></td></tr></tbody></table>

Note: (RTCA) and (r2dtable) indicate which null model was used.
<figure><img src="image_16.png" /><figcaption>**FIGURE 3** Correlation strength for all taxonomic levels and all indices. All correlations are significant. BR, discrepancy; NODF, nestedness index based on overlap and decreasing fill; SR_Bin, spectral radius calculated on binary (absence/presence) matrices; SR_Qua, spectral radius calculated on weighted (abundance) matrices; VS, versus; WNODF, NODF calculated on weighted matrices</figcaption></figure>

We also showed that the normalization method proposed by Song et al. (2017) offers an NODF index robust to taxonomic resolution, which is in line with their own conclusions that NODF<sub>c</sub> is independent from network number of rows, columns, and number of links, making it remarkably relevant to compare networks across studies or spatial gradients. Our results are also in agreement with an as yet unpublished study by Hemprich-Bennett, Oliveira, Comber, Rossiter, and Clare, (2018) who also found that absolute measures of most of the metrics they tested (which includes NODF, robustness, connectivity, but not SR or BR) vary according to the taxonomical level of the networks (both observed networks and networks deduced from metabarcoding data).

The main objective of our study was to determine whether it is possible to meaningfully estimate indices describing the characteristics of plant-pollinator interaction networks with a taxonomic resolution lower than the species. Our results suggest that it is indeed the case. To estimate nestedness, our results suggest that only the absolute values of SR indices are minimally impacted by changes in taxonomic resolution and should therefore probably be preferred when the objective is to compare nestedness levels for networks with a lower resolution than the species. Alternatively, it is possible to use NODF and BR after standardization using the Z-score. For the other properties of the networks that we examined, namely connectivity, modularity, normalized degree, and robustness to species loss, the absolute values of the indices cannot be directly compared at resolutions lower than the species level, but it is still possible to rank networks according to their values for these indices, because such ranks are well preserved when the level of taxonomic resolution changes. Motif frequencies do not present a unique pattern of sensitivity to taxonomic resolution. Indeed, some motif frequencies are significantly influenced by taxonomic resolution, while others are not. However, they all show a good preservation of the ranks between species and genus networks. We would advise not to use motif frequencies at a family level, though, as the correlation between ranks gets rather low (sometimes as low as 0.2).

Note, however, that whereas the taxonomic resolution lower than the species seems to allow us to characterize the properties of plant-pollinator interaction networks, it may make it more difficult to interpret these properties. One of the main objectives of the measurement of network properties is to make or test inferences about their underlying mechanisms. For example, Junker et al. (2013) showed that sets of plant traits such as phenology, floral reflectance, and morphology can predict plant-pollinator interactions and thus network structure. Similarly, Klumpers, Stang, and Klinkhamer (2019) showed that size matching between the pollinator proboscis length and the nectar tube depth is important in shaping plant-pollinator interactions. Such conclusions would be more difficult to reach when working above the species level. In the future, working on these levels would require careful consideration: Can functional traits be extended to the whole genus in that particular case? If this is not possible, then working on these levels could thus deprive us of a significantly explanatory variable. This means that while genus- and family-level networks are usable and interpretable, they still entail a loss of information for future studies. For this reason, future studies need to consider the gain in network explicitness versus the loss of information before choosing to work at the genus or family level.

Our results support the relevance of citizen science for ecological research. The major strengths of citizen science programs lie in their ability to conduct studies at large geographic scales and on private properties, which are usually impossible to perform with traditional field research (Dickinson et al., 2010), although these are often at the price of a lower taxonomical precision. Here, we showed that datasets with a taxonomic resolution lower than the species level can be used to estimate the properties of networks assembled at the same resolution, even if it is lower than the
species. However, plant-pollinator interaction data produced by citizen science are probably characterized by relatively low sampling completeness, because detecting all the species interactions is extremely labor-intensive (Chacoff et al., 2012), which can have an effect on the estimated properties of the networks. For the indices that we studied, Rivera-Hutinel, Bustamante, Marín, and Medel (2012) showed that nestedness, modularity, and robustness to species loss are little affected by sampling completeness, whereas connectance is very sensitive to low sampling. In conclusion, sets of plant-pollinator networks produced by citizen science, frequently characterized by low taxonomic resolution and low sampling efforts, are probably best analyzed by calculating their nestedness with SR (or NODF and BR after standardization using the Z-score) and their robustness with species loss, and then ranking them according to their modularity.

Our work confirms that we can use protocols with only genus- or family-level data and still use network-level analyses of plant-pollinator interactions. An interesting complement would be to study the same question for other kinds of mutualistic networks such as ant-plant networks or even for other kinds of interaction networks such as food webs.

## ACKNOWLEDGMENTS

We thank the ENS de Lyon, who funded E.R. during the project.

## AUTHORS' CONTRIBUTIONS

CBG and EB conceived the ideas and designed methodology. ER collected the data. ER, CBG, and EB analyzed the data. ER led the writing of the manuscript. All authors contributed critically to the drafts and gave final approval for publication.

## DATA AVAILABILITY STATEMENT

The networks used in this analysis are available at the Interaction Web Database (https://www.nceas.ucsb.edu/interactionweb/resources.html#plant_pollinator), which collects datasets from authors publishing their work. Some of them, however, were not used due to the difficulty of treatment. The list of articles used in our study is available in Supplementary Materials. The complete script used to obtain the result and diverse networks is accessible here: https://gitlab.com/EstelleRenaud/taxonomic_influence_network_properties

## ORCID

Estelle Renaud https://orcid.org/0000-0003-0963-2648

## REFERENCES

* Almeida-Neto, M., Guimarães, P., Guimarães, P. R., Loyola, R. D., & Ulrich, W. (2008). A consistent metric for nestedness analysis in ecological systems: Reconciling concept and measurement. *Oikos*, *117*(8), 1227-1239. https://doi.org/10.1111/j.0030-1299.2008.16644.x

* Almeida-Neto, M., & Ulrich, W. (2011). A straightforward computational approach for measuring nestedness using quantitative matrices. *Environmental Modelling & Software*, *26*(2), 173-178. https://doi.org/10.1016/j.envsoft.2010.08.003

* Bascompte, J., Jordano, P., Melián, C. J., & Olesen, J. M. (2003). The nested assembly of plant-animal mutualistic networks. *Proceedings of the National Academy of Sciences of the United States of America*, *100*(16), 9383-9387. https://doi.org/10.1073/pnas.1633576100

* Bastolla, U., Fortuna, M. A., Pascual-García, A., Ferrera, A., Luque, B., & Bascompte, J. (2009). The architecture of mutualistic networks minimizes competition and increases biodiversity. *Nature*, *458*(7241), 1018-1020. https://doi.org/10.1038/nature07950

* Beckett, S. J. (2016). Improved community detection in weighted bipartite networks. *Royal Society Open Science*, *3*(1), 140536. https://doi.org/10.1098/rsos.140536

* Beckett, S. J., Boulton, C. A., & Williams, H. T. P. (2014). FALCON: A software package for analysis of nestedness in bipartite networks. *F1000Research*, *3*, 185. https://doi.org/10.12688/f1000resea rch.4831.1

* Brualdi, R. A., & Sanderson, J. G. (1999). Nested species subsets, gaps, and discrepancy. *Oecologia*, *119*(2), 256-264. https://doi.org/10.1007/s004420050784

* Burgos, E., Ceva, H., Perazzo, R. P. J., Devoto, M., Medan, D., Zimmermann, M., & Maria Delbue, A. (2007). Why nestedness in mutualistic networks? *Journal of Theoretical Biology*, *249*(2), 307-313. https://doi.org/10.1016/j.jtbi.2007.07.030

* Chacoff, N. P., Vázquez, D. P., Lomáscolo, S. B., Stevani, E. L., Dorado, J., & Padrón, B. (2012). Evaluating sampling completeness in a desert plant-pollinator network. *Journal of Animal Ecology*, *81*(1), 190-200. https://doi.org/10.1111/j.1365-2656.2011.01883.x

* Chamberlain, S. A., & Szöcs, E. (2013). taxize: Taxonomic search and retrieval in R. *F1000Research*, *2*, 191. https://doi.org/10.12688/f1000 research.2-191.v2

* Deguines, N., Julliard, R., de Flores, M., & Fontaine, C. (2012). The whereabouts of flower visitors: Contrasting land-use preferences revealed by a country-wide survey based on citizen science. *PLoS ONE*, *7*(9), e45822. https://doi.org/10.1371/journal.pone.0045822

* Delmas, E., Besson, M., Brice, M.-H., Burkle, L. A., Dalla Riva, G. V., Fortin, M.-J., ... Poisot, T. (2019). Analysing ecological networks of species interactions. *Biological Reviews*, *94*(1), 16-36. https://doi.org/10.1111/brv.12433

* Dickinson, J. L., Zuckerberg, B., & Bonter, D. N. (2010). Citizen science as an ecological research tool: Challenges and benefits. *Annual Review of Ecology, Evolution, and Systematics*, *41*(1), 149-172. https://doi.org/10.1146/annurev-ecolsys-102209-144636

* Dormann, C. F., Gruber, B., & Fründ, J. (2008). Introducing the bipartite package: Analysing ecological networks. *Rnews*, *8*. Retrieved from http://www.uni-goettingen.de/de/document/download/96729 eb9d30a6f2dc4403df15854305c.pdf/Rnews2008,8_8-11_open.pdf

* Garibaldi, L. A., Steffan-Dewenter, I., Winfree, R., Aizen, M. A., Bommarco, R., Cunningham, S. A., ... Klein, A. M. (2013). Wild pollinators enhance fruit set of crops regardless of honey bee abundance. *Science*, *339*(6127), 1608-1611. https://doi.org/10.1126/science.1230200

* GBIF (2018). Retrieved from https://www.gbif.org/

* Gotelli, N. J. (2000). Null model analysis of species co-occurrence patterns. *Ecology*, *81*(9), 2606-2621. https://doi.org/10.1890/0012-9658(2000)081[2606: NMAOSC]2.0.CO;2

* Goulson, D., Nicholls, E., Botías, C., & Rotheray, E. L. (2015). Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. *Science*, *347*(6229), 1255957. https://doi.org/10.1126/science ce.1255957

* Hemprich-Bennett, D. R., Oliveira, H. F. M., Comber, S. C. L., Rossiter, S. J., & Clare, E. L. (2018). Assessing the impact of taxon resolution on network structure, with implication for comparative ecology. *BioRxiv*, 357376. https://doi.org/10.1101/357376

* Hoeppke, C. (2019). CHoeppke/maxnodf. R. Retrieved from https://github.com/CHoeppke/maxnodf (Original work published 8 August 2018)
Ings, T. C., Montoya, J. M., Bascompte, J., Blüthgen, N., Brown, L., Dormann, C. F., ... Woodward, G. (2009). Review: Ecological networks - beyond food webs. *Journal of Animal Ecology*, 78(1), 253-269.
https://doi.org/10.1111/j.1365-2656.2008.01460.x

Junker, R. R., Blüthgen, N., Brehm, T., Binkenstein, J., Paulus, J., Martin Schaefer, H., & Stang, M. (2013). Specialization on traits as basis for the niche-breadth of flower visitors and as structuring mechanism of ecological networks. *Functional Ecology*, 27(2), 329-341. https://doi.org/10.1111/1365-2435.12005

Kaiser-Bunbury, C. N., & Blüthgen, N. (2015). Integrating network ecology with applied conservation: A synthesis and guide to implementation. *AoB Plants*, 7, 1-15. https://doi.org/10.1093/aobpla/plv076

Kaiser-Bunbury, C. N., Mougal, J., Whittington, A. E., Valentin, T., Gabriel, R., Olesen, J. M., & Blüthgen, N. (2017). Ecosystem restoration strengthens pollination network resilience and function. *Nature*, 542(7640), 223-227. https://doi.org/10.1038/nature21071

Klumpers, S. G. T., Stang, M., & Klinkhamer, P. G. L. (2019). Foraging efficiency and size matching in a plant-pollinator community: The importance of sugar content and tongue length. *Ecology Letters*, 22(3), 469-479. https://doi.org/10.1111/ele.13204

Kremen, C., Ullman, K. S., & Thorp, R. W. (2011). Evaluating the quality of citizen-scientist data on pollinator communities. *Conservation Biology*, 25(3), 607-617. https://doi.org/10.1111/j.1523-1739.2011.01657.x

Lau, M. K., Borrett, S. R., Baiser, B., Gotelli, N. J., & Ellison, A. M. (2017). Ecological network metrics: Opportunities for synthesis. *Ecosphere*, 8(8), e01900. https://doi.org/10.1002/ecs2.1900

Martín González, A. M., Dalsgaard, B., & Olesen, J. M. (2010). Centrality measures and the importance of generalist species in pollination networks. *Ecological Complexity*, 7(1), 36-43. https://doi.org/10.1016/j.ecocom.2009.03.008

Memmott, J. (2009). Food webs: A ladder for picking strawberries or a practical tool for practical problems? *Philosophical Transactions of the Royal Society B: Biological Sciences*, 364(1524), 1693-1699. https://doi.org/10.1098/rstb.2008.0255

Milo, R., Shen-Orr, S., Itzkovitz, S., Kashtan, N., Chklovskii, D., & Alon, U. (2002). Network motifs: Simple building blocks of complex networks. *Science*, 298(5594), 824-827. https://doi.org/10.1126/science.298.5594.824

Montoya, J. M., Pimm, S. L., & Solé, R. V. (2006). Ecological networks and their fragility. *Nature*, 442(7100), 259-264. https://doi.org/10.1038/nature04927

Newman, M. E. (2006). Modularity and community structure in networks. *Proceedings of the National Academy of Sciences of the United of America*, 103(23), 8577-8582. https://doi.org/10.1073/pnas.060160210

Proulx, S. R., Promislow, D. E. L., & Phillips, P. C. (2005). Network thinking in ecology and evolution. *Trends in Ecology & Evolution*, 20(6), 345-353. https://doi.org/10.1016/j.tree.2005.04.004

Rivera-Hutinel, A., Bustamante, R. O., Marín, V. H., & Medel, R. (2012). Effects of sampling completeness on the structure of plant-pollinator networks. *Ecology*, 93(7), 1593-1603. https://doi.org/10.1890/11-1803.1

Rodríguez-Gironés, M. A., & Santamaría, L. (2006). A new algorithm to calculate the nestedness temperature of presence-absence matrices. *Journal of Biogeography*, 33(5), 924-935. https://doi.org/10.1111/j.1365-2699.2006.01444.x

Simmons, B. I., Sweering, M. J. M., Schillinger, M., Dicks, L. V., Sutherland, W. J., & Clemente, R. D. (2019). bmotif: A package for motif analyses of bipartite networks. *Methods in Ecology and Evolution*, 10(5), 695-701. https://doi.org/10.1111/2041-210X.13149

Soares, R. G. S., Ferreira, P. A., & Lopes, L. E. (2017). Can plant-pollinator network metrics indicate environmental quality? *Ecological Indicators*, 78, 361-370. https://doi.org/10.1016/j.ecolind.2017.03.037

Song, C., Rohr, R. P., & Saavedra, S. (2017). Why are some plant-pollinator networks more nested than others? *The Journal of Animal Ecology*, 86(6), 1417-1424. https://doi.org/10.1111/1365-2656.12749

Staniczenko, P. P. A., Kopp, J. C., & Allesina, S. (2013). The ghost of nestedness in ecological networks. *Nature Communications*, 4, 1391. https://doi.org/10.1038/ncomms2422

Strona, G., & Fattorini, S. (2014). On the methods to assess significance in nestedness analyses. *Theory in Biosciences*, 133(3), 179-186. https://doi.org/10.1007/s12064-014-0203-1

Strona, G., Ulrich, W., & Gotelli, N. J. (2018). Bi-dimensional null model analysis of presence-absence binary matrices. *Ecology*, 99(1), 103-115. https://doi.org/10.1002/ecy.2043

Thébault, E., & Fontaine, C. (2010). Stability of ecological communities and the architecture of mutualistic and trophic networks. *Science*, 329(5993), 853-856. https://doi.org/10.1126/science.1188321

Toomey, A. H., & Domroese, M. C. (2013). Can citizen science lead to positive conservation attitudes and behaviors? *Human Ecology Review*, 20(1), 50-62.

Ulrich, W., & Almeida-Neto, M. (2012). On the meanings of nestedness: Back to the basics. *Ecography*, 35(10), 865-871. https://doi.org/10.1111/j.1600-0587.2012.07671.x

Vanbergen, A. J., Woodcock, B. A., Heard, M. S., & Chapman, D. S. (2017). Network size, structure and mutualism dependence affect the propensity for plant-pollinator extinction cascades. *Functional Ecology*, 31(6), 1285-1293. https://doi.org/10.1111/1365-2435.12823

Vázquez, D., Goldberg, J., & Naik, R. (2003). *Interaction web database*. Retrieved from https://www.nceas.ucsb.edu/interactionweb/resources.html#plant_pollinator

Wright, D. H., Patterson, B. D., Mikkelson, G. M., Cutler, A., & Atmar, W. (1997). A comparative analysis of nested subset patterns of species composition. *Oecologia*, 113(1), 1-20. https://doi.org/10.1007/s004420050348

## SUPPORTING INFORMATION

Additional supporting information may be found online in the Supporting Information section.

**How to cite this article:** Renaud E, Baudry E, Bessa-Gomes C. Influence of taxonomic resolution on mutualistic network properties. *Ecol Evol.* 2020;00:1-13. https://doi.org/10.1002/ece3.6060
APPENDIX A

<table><caption>TABLE A1 Results of one-way analysis of variance comparing normalized degree (ND) values and motif frequencies at three different taxonomic resolutions</caption><thead><tr><th>Indices</th><th>df</th><th>F</th><th>p-Value</th><th>Post hoc test result</th></tr></thead><tbody><tr><td>Normalized degrees for pollinators—lower limit of 1st quartile</td><td>2, 120</td><td>5.7164</td><td>.004252</td><td>All levels differ significantly</td></tr><tr><td>ND_high_2</td><td>2, 120</td><td>10.92</td><td>4.396e-05</td><td>All levels differ significantly</td></tr><tr><td>ND_high_3</td><td>2, 120</td><td>13.129</td><td>6.977e-06</td><td>All levels differ significantly</td></tr><tr><td>ND_high_4</td><td>2, 120</td><td>18.155</td><td>1.293e-07</td><td>All levels differ significantly</td></tr><tr><td>ND_high_5</td><td>2, 120</td><td>9.2779</td><td>.0001792</td><td>All levels differ significantly</td></tr><tr><td>ND_low_1</td><td>2, 120</td><td>20.594</td><td>2.046e-08</td><td>All levels differ significantly</td></tr><tr><td>ND_low_2</td><td>2, 120</td><td>15.191</td><td>1.316e-06</td><td>All levels differ significantly</td></tr><tr><td>ND_low_3</td><td>2, 120</td><td>11.922</td><td>1.893e-05</td><td>All levels differ significantly</td></tr><tr><td>ND_low_4</td><td>2, 120</td><td>18.947</td><td>7.06e-08</td><td>All levels differ significantly</td></tr><tr><td>ND_low_5</td><td>2, 120</td><td>13.857</td><td>3.85e-06</td><td>All levels differ significantly</td></tr><tr><td>Motif 1</td><td>2, 120</td><td>1</td><td>.3709</td><td></td></tr><tr><td>Motif 2</td><td>2, 120</td><td>1.1355</td><td>.3247</td><td></td></tr><tr><td>Motif 3</td><td>2, 120</td><td>1.1355</td><td>.3247</td><td></td></tr><tr><td>Motif 4</td><td>2, 120</td><td>0.112</td><td>.8942</td><td></td></tr><tr><td>Motif 5</td><td>2, 120</td><td>2.1753</td><td>.118</td><td></td></tr><tr><td>Motif 6</td><td>2, 120</td><td>9.3845</td><td>.0001634</td><td>All levels differ significantly</td></tr><tr><td>Motif 7</td><td>2, 120</td><td>2.5995</td><td>.07849</td><td></td></tr><tr><td>Motif 8</td><td>2, 120</td><td>0.1919</td><td>.8256</td><td></td></tr><tr><td>Motif 9</td><td>2, 120</td><td>0.2043</td><td>.8155</td><td></td></tr><tr><td>Motif 10</td><td>2, 120</td><td>1.7185</td><td>.1837</td><td></td></tr><tr><td>Motif 11</td><td>2, 120</td><td>9.5841</td><td>.0001376</td><td>All levels differ significantly</td></tr><tr><td>Motif 12</td><td>2, 120</td><td>5.4115</td><td>.005621</td><td>All levels differ significantly</td></tr><tr><td>Motif 13</td><td>2, 120</td><td>0.0374</td><td>.9633</td><td></td></tr><tr><td>Motif 14</td><td>2, 120</td><td>0.9314</td><td>.3968</td><td></td></tr><tr><td>Motif 15</td><td>2, 120</td><td>3.4154</td><td>.03609</td><td>Species and genus do not differ significantly, while the two other comparisons are significantly different</td></tr><tr><td>Motif 16</td><td>2, 120</td><td>4.8719</td><td>.00924</td><td>Species and genus do not differ significantly, while the two other comparisons are significantly different</td></tr><tr><td>Motif 17</td><td>2, 120</td><td>3.4329</td><td>.0355</td><td>All levels differ significantly</td></tr></tbody></table>

Note: Results are compiled for pollinator ("high") and plant ("low") taxa, for each bound of the quartile values (1–5).
<table><caption>TABLE A2 Results of Spearman's rho correlation tests comparing the relative values of all networks between two taxonomic levels, for each index</caption><thead><tr><th rowspan="2">Index</th><th colspan="2">Species vs Genus</th><th colspan="2">Genus vs Family</th><th colspan="2">Species vs Family</th></tr><tr><th>Rho estimate</th><th>p-Value</th><th>Rho estimate</th><th>p-Value</th><th>Rho estimate</th><th>p-Value</th></tr></thead><tbody><tr><td>NODF</td><td>.914</td><td>&lt;2.2e-16</td><td>.751</td><td>1.611e-07</td><td>.617</td><td>2.62e-05</td></tr><tr><td>WNODF</td><td>.926</td><td>1.541e-08</td><td>.6375</td><td>.0001612</td><td>.467</td><td>.008709</td></tr><tr><td>BR</td><td>.959</td><td>&lt;2.2e-16</td><td>.797</td><td>4.478e-10</td><td>.678</td><td>1.093e-06</td></tr><tr><td>SR Binary</td><td>.962</td><td>&lt;2.2e-16</td><td>.956</td><td>&lt;2.2e-16</td><td>.903</td><td>&lt;2.2e-16</td></tr><tr><td>SR Weighted</td><td>.988</td><td>&lt;2.2e-16</td><td>.960</td><td>&lt;2.2e-16</td><td>.939</td><td>&lt;2.2e-16</td></tr><tr><td>Connectance</td><td>.968</td><td>&lt;2.2e-16</td><td>.948</td><td>2.705e-12</td><td>.754</td><td>1.224e-08</td></tr><tr><td>Modularity</td><td>.952</td><td>&lt;2.2e-16</td><td>.900</td><td>1.343e-07</td><td>.863</td><td>3.964e-07</td></tr><tr><td>Robustness</td><td>.921</td><td>&lt;2.2e-16</td><td>.793</td><td>6.482e-10</td><td>.650</td><td>4.219e-06</td></tr><tr><td>ND_high_1</td><td>.994</td><td>&lt;2.2e-16</td><td>.828</td><td>2.357e-11</td><td>.823</td><td>3.734e-11</td></tr><tr><td>ND_high_2</td><td>.904</td><td>5.834e-16</td><td>.819</td><td>6.084e-11</td><td>.681</td><td>9.452e-07</td></tr><tr><td>ND_high_3</td><td>.933</td><td>&lt;2.2e-16</td><td>.776</td><td>2.533e-09</td><td>.712</td><td>1.786e-07</td></tr><tr><td>ND_high_4</td><td>.914</td><td>&lt;2.2e-16</td><td>.824</td><td>3.773e-11</td><td>.724</td><td>8.565e-08</td></tr><tr><td>ND_high_5</td><td>.886</td><td>1.436e-14</td><td>.693</td><td>5.144e-07</td><td>.634</td><td>8.632e-06</td></tr><tr><td>ND_low_1</td><td>.923</td><td>&lt;2.2e-06</td><td>.741</td><td>2.984e-08</td><td>.641</td><td>6.3e-06</td></tr><tr><td>ND_low_2</td><td>.890</td><td>6.742e-15</td><td>.762</td><td>7.313e-09</td><td>.726</td><td>7.695e-08</td></tr><tr><td>ND_low_3</td><td>.940</td><td>&lt;2.2e-16</td><td>.779</td><td>1.947e-09</td><td>.666</td><td>2.013e-06</td></tr><tr><td>ND_low_4</td><td>.936</td><td>&lt;2.2e-16</td><td>.853</td><td>1.453e-12</td><td>.721</td><td>1.039e-07</td></tr><tr><td>ND_low_5</td><td>.910</td><td>&lt;2.2e-16</td><td>.734</td><td>4.871e-08</td><td>.605</td><td>2.791e-05</td></tr><tr><td>Motif 1</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr><tr><td>Motif 2</td><td>.935</td><td>&lt;2.2e-16</td><td>.641</td><td>6.319e-06</td><td>.575</td><td>8.341e-05</td></tr><tr><td>Motif 3</td><td>.935</td><td>&lt;2.2e-16</td><td>.641</td><td>6.319e-06</td><td>.575</td><td>8.341e-05</td></tr><tr><td>Motif 4</td><td>.935</td><td>&lt;2.2e-16</td><td>.756</td><td>1.356e-07</td><td>.685</td><td>1.742e-06</td></tr><tr><td>Motif 5</td><td>.859</td><td>&lt;2.2e-16</td><td>.661</td><td>4.746e-06</td><td>.633</td><td>1.404e-05</td></tr><tr><td>Motif 6</td><td>.906</td><td>&lt;2.2e-16</td><td>.702</td><td>3.201e-07</td><td>.558</td><td>.0001514</td></tr><tr><td>Motif 7</td><td>.930</td><td>&lt;2.2e-16</td><td>.614</td><td>1.989e-05</td><td>.578</td><td>7.461e-05</td></tr><tr><td>Motif 8</td><td>.927</td><td>&lt;2.2e-16</td><td>.793</td><td>6.482e-10</td><td>.735</td><td>4.315e-08</td></tr><tr><td>Motif 9</td><td>.937</td><td>&lt;2.2e-16</td><td>.703</td><td>8.524e-07</td><td>.631</td><td>1.542e-05</td></tr><tr><td>Motif 10</td><td>.921</td><td>&lt;2.2e-16</td><td>.402</td><td>.009617</td><td>.399</td><td>0.01027</td></tr><tr><td>Motif 11</td><td>.913</td><td>&lt;2.2e-16</td><td>.418</td><td>.006511</td><td>.285</td><td>0.07137</td></tr><tr><td>Motif 12</td><td>.953</td><td>&lt;2.2e-16</td><td>.708</td><td>2.208e-07</td><td>.608</td><td>2.464e-05</td></tr><tr><td>Motif 13</td><td>.900</td><td>&lt;2.2e-16</td><td>.585</td><td>8.009e-05</td><td>.435</td><td>.004839</td></tr><tr><td>Motif 14</td><td>.788</td><td>3.833e-08</td><td>.617</td><td>1.742e-05</td><td>.519</td><td>.0005086</td></tr><tr><td>Motif 15</td><td>.807</td><td>6.773e-09</td><td>.742</td><td>2.875e-08</td><td>.508</td><td>.0006918</td></tr><tr><td>Motif 16</td><td>.815</td><td>8.911e-11</td><td>.830</td><td>2.023e-11</td><td>.644</td><td>5.625e-06</td></tr><tr><td>Motif 17</td><td>.933</td><td>&lt;2.2e-16</td><td>.606</td><td>2.661e-05</td><td>.560</td><td>.0001424</td></tr></tbody></table>

Abbreviation: BR, discrepancy; ND, normalized degree compiled for pollinator ("high") and plant ("low") taxa, for each bound of the quartile values (1-5); NODF, nestedness index based on overlap and decreasing fill; SR_Bin, spectral radius calculated on binary (absence/presence) matrices; SR_Qua, spectral radius calculated on weighted (abundance) matrices; WNODF, NODF calculated on weighted matrices.