Community Phylogenetics: Assessing Tree Reconstruction Methods and the Utility of DNA Barcodes

Origins of the central Macaronesian psyllid lineages (Hemiptera; Psylloidea) with characterization of a new island radiation on endemic Convolvulus floridus (Convolvulaceae) in the Canary Islands

A molecular survey of native and adventive psyllids in the central Macaronesian islands provides the first comprehensive phylogenetic assessment of the origins of the psyllid fauna of the Canary and Madeira archipelagos. We employ a maximum likelihood backbone constraint analysis to place the central Macaronesian taxa within the Psylloidea mitogenome phylogeny. The native psyllid fauna in these central Macaronesian islands results from an estimated 26 independent colonization events. Island host plants are predicted by host plants of continental relatives in nearly all cases and six plant genera have been colonized multiple times (Chamaecytisus, Convolvulus, Olea, Pistacia, Rhamnus, and Spartocytisus) from the continent. Post-colonization diversification varies from no further cladogenesis (18 events, represented by a single native taxon) to modest in situ diversification resulting in two to four native taxa and, surprisingly, given the diverse range of islands and habitats, only one substantial species radiation with more than four native species. Specificity to ancestral host plant genera or family is typically maintained during in situ diversification both within and among islands. Characterization of a recently discovered island radiation consisting of four species on Convolvulus floridus in the Canary Islands shows patterns and rates of diversification that reflect island topographic complexity and geological dynamism. Although modest in species diversity, this radiation is atypical in diversification on a single host plant species, but typical in the primary role of allopatry in the diversification process.

Launching insectphylo.org; a new hub facilitating construction and use of synthesis molecular phylogenies of insects

The Holy Grail of an Insect Tree of Life can only be 'discovered' through extensive collaboration among taxon specialists, phylogeneticists and centralized frameworks such as Open Tree of Life, but insufficient effort from stakeholders has so far hampered this promising approach. The resultant unavailability of synthesis phylogenies is an unfortunate situation given the numerous practical usages of phylogenies in the near term and against the backdrop of the ongoing biodiversity crisis. To resolve this issue, we establish a new online hub that centralizes the collation of relevant phylogenetic data and provides the resultant synthesis molecular phylogenies. This is achieved through key developments in a proposed pipeline for the construction of a species-level insect phylogeny. The functionality of the framework is demonstrated through the construction of a highly supported, species-comprehensive phylogeny of Diptera, built from integrated omics data, COI DNA barcodes, and a compiled database of over 100 standardized, published Diptera phylogenies. Machine-readable forms of the phylogeny (and subsets thereof) are publicly available at insectphylo.org, a new public repository for species-comprehensive phylogenies for biological research.

Detection of Microsporidia in Pollinator Communities of a Mediterranean Biodiversity Hotspot for Wild Bees

Insect pollination is crucial for the maintenance of natural and managed ecosystems but the functioning of this ecosystem service is threatened by a worldwide decline of pollinators. Key factors in this situation include the spread and interspecific transmission of pathogens worldwide through the movement of managed pollinators. Research on this field has been mainly conducted in some particular species, while studies assessing the interspecific transmission of pathogens at a community level are scarce. However, this information is pivotal to design strategies to protect pollinators. Herein, we analysed the prevalence of two common microsporidia pathogens of managed honey bees (Nosema ceranae and N. apis) in bee communities of semiarid Mediterranean areas from the Southeast of the Iberian Peninsula. Our results confirm the ability of N. ceranae to disperse across wild bee communities in semiarid Mediterranean ecosystems since it was detected in 36 Apoidea species (39% of the sampling; for the first time in nine genera). The prevalence of the pathogen did not show any phylogenetic signal which suggests a superfamily host range of the pathogen or that wild bees may be acting only as vectors of N. ceranae. In addition, N. apis was detected in an Eucera species, which is the second time it has been detected by molecular techniques in a host other than the honey bee. Our study represents the primary assessment of the prevalence of microsporidia at community level in Mediterranean areas and provides outstanding results on the ability of Nosema pathogens to spread across the landscape.

Phylogenetic signal of sub-arctic beetle communities

Postglacial dispersal and colonization processes have shaped community patterns in sub-Arctic regions such as Churchill, Manitoba, and Canada. This study investigates evolutionary community structure within the beetle (Coleoptera) families of Churchill and tests whether biological traits have played a role in governing colonization patterns from refugial and southerly geographic regions. This study quantifies sub-Arctic beetle phylogenetic community structure for each family using the net relatedness index (NRI) and nearest taxon index (NTI), calculated using publicly available data from the Barcode of Life Data Systems (BOLD); compares patterns across families with different traits (habitat, diet) using standard statistical analysis (ANOVA) as well as phylogenetic generalized least squares (PGLS) using a family-level beetle phylogeny obtained from the literature; and compares community structure in Churchill with a region in southern Canada (Guelph, Ontario). These analyses were also repeated at a genus level. The dominant pattern detected in our study was that aquatic families were much better represented in Churchill compared to terrestrial families, when compared against richness sampled from across Canada and Alaska. Individually, most families showed significant phylogenetic clustering in Churchill, likely due to the strong environmental filtering present in Arctic environments. There was no significant difference in phylogenetic structure between Churchill and Guelph but with a trend toward stronger clustering in the North. Fungivores were significantly more overdispersed than other feeding modes, predators were significantly more clustered, and aquatic families showed significantly stronger clustering compared to terrestrial. This study contributes to our understanding of the traits and processes structuring insect biodiversity and macroecological trends in the sub-Arctic.

Plastome-based phylogeny improves community phylogenetics of subtropical forests in China

Phylogenetic trees have been extensively used in community ecology. However, how the phylogeny construction affects ecological inferences is poorly understood. In this study, we constructed three different types of phylogenetic trees (a synthetic-tree generated using V.PhyloMaker, a barcode-tree generated using rbcL+matK+trnH-psbA, and a plastome-tree generated from plastid genomes) that represented an increasing level of phylogenetic resolution among 580 woody plant species from six forest dynamic plots in subtropical evergreen broadleaved forests of China. We then evaluated the performance of each phylogeny in estimations of community phylogenetic structure, turnover and phylogenetic signal in functional traits. As expected, the plastome-tree was most resolved and most supported for relationships among species. For local phylogenetic structure, the three trees showed consistent results with Faith's PD and MPD; however, only the synthetic-tree produced significant clustering patterns using MNTD for some plots. For phylogenetic turnover, contrasting results between the molecular trees and the synthetic-tree occurred only with nearest neighbor distance. The barcode-tree agreed more with the plastome-tree than the synthetic-tree for both phylogenetic structure and turnover. For functional traits, both the barcode-tree and plastome-tree detected phylogenetic signal in maximum height, but only the plastome-tree detected signal in leaf width. This is the first study that uses plastid genomes in large-scale community phylogenetics. Our results highlight the improvement of plastome-trees over barcode-trees and synthetic-trees for the analyses studied here. Our results also point to the possibility of type I and II errors in estimation of phylogenetic structure and turnover and detection of phylogenetic signal when using synthetic-trees.

Phylogenetic uncertainty and the inference of patterns in community ecology and comparative studies

Progress in phylogenetic community ecology is often limited by the availability of phylogenetic information and the lack of appropriate methods and solutions to deal with this problem. We estimate the effect of the lack of phylogenetic information on the relations among taxa measured by commonly used phylogenetic metrics in comparative studies and community ecology, namely: Blomberg's K phylogenetic signal, Faith's Phylogenetic Diversity (PD), Mean Phylogenetic Distance (MPD) and Mean Nearest Taxon Distance (MNTD). To overcome this problem, we tested two possible solutions: Polytomic trees and Operational trees. Our results show that the effects on K values strongly depended on the level of phylogenetic signal. In the case of the community metrics, the effects were insensitive to the patterns of species distribution in the communities. Community metrics tended to be overestimated with both Polytomic and Operational trees, but the overestimation was higher with Polytomic trees. PD and MPD metrics were less biased than MNTD metric. We show that the lack of phylogenetic resolution is not necessarily problematic for all analyses and that its effect will depend on the chosen metric and on the solutions used to deal with the problem. Based on our results, we suggest that ecologists should prefer the Operational tree solution to remove polytomies in the phylogenetic tree and take careful consideration while designing experiments, and analyzing and interpreting the results of phylogenetic metrics.

Building a Robust, Densely-Sampled Spider Tree of Life for Ecosystem Research

Phylogenetic relatedness is a key diversity measure for the analysis and understanding of how species and communities evolve across time and space. Understanding the nonrandom loss of species with respect to phylogeny is also essential for better-informed conservation decisions. However, several factors are known to influence phylogenetic reconstruction and, ultimately, phylogenetic diversity metrics. In this study, we empirically tested how some of these factors (topological constraint, taxon sampling, genetic markers and calibration) affect phylogenetic resolution and uncertainty. We built a densely sampled, species-level phylogenetic tree for spiders, combining Sanger sequencing of species from local communities of two biogeographical regions (Iberian Peninsula and Macaronesia) with a taxon-rich backbone matrix of Genbank sequences and a topological constraint derived from recent phylogenomic studies. The resulting tree constitutes the most complete spider phylogeny to date, both in terms of terminals and background information, and may serve as a standard reference for the analysis of phylogenetic diversity patterns at the community level. We then used this tree to investigate how partial data affect phylogenetic reconstruction, phylogenetic diversity estimates and their rankings, and, ultimately, the ecological processes inferred for each community. We found that the incorporation of a single slowly evolving marker (28S) to the DNA barcode sequences from local communities, had the highest impact on tree topology, closely followed by the use of a backbone matrix. The increase in missing data resulting from combining partial sequences from local communities only had a moderate impact on the resulting trees, similar to the difference observed when using topological constraints. Our study further revealed substantial differences in both the phylogenetic structure and diversity rankings of the analyzed communities estimated from the different phylogenetic treatments, especially when using non-ultrametric trees (phylograms) instead of time-stamped trees (chronograms). Finally, we provide some recommendations on reconstructing phylogenetic trees to infer phylogenetic diversity within ecological studies.

Distribution and Community Assembly of Trees Along an Andean Elevational Gradient

Highlighting patterns of distribution and assembly of plants involves the use of community phylogenetic analyses and complementary traditional taxonomic metrics. However, these patterns are often unknown or in dispute, particularly along elevational gradients, with studies finding different patterns based on elevation. We investigated how patterns of tree diversity and structure change along an elevation gradient using taxonomic and phylogenetic diversity metrics. We sampled 595 individuals (36 families; 53 genera; 88 species) across 15 plots along an elevational gradient (2440-3330 m) in Ecuador. Seventy species were sequenced for the rbcL and matK gene regions to generate a phylogeny. Species richness, Shannon-Weaver diversity, Simpson's Dominance, Simpson's Evenness, phylogenetic diversity (PD), mean pairwise distance (MPD), and mean nearest taxon distance (MNTD) were evaluated for each plot. Values were correlated with elevation and standardized effect sizes (SES) of MPD and MNTD were generated, including and excluding tree fern species, for comparisons across elevation. Taxonomic and phylogenetic metrics found that species diversity decreases with elevation. We also found that overall the community has a non-random phylogenetic structure, dependent on the presence of tree ferns, with stronger phylogenetic clustering at high elevations. Combined, this evidence supports the ideas that tree ferns have converged with angiosperms to occupy the same habitat and that an increased filtering of clades has led to more closely related angiosperm species at higher elevations.

For common community phylogenetic analyses, go ahead and use synthesis phylogenies

Should we build our own phylogenetic trees based on gene sequence data, or can we simply use available synthesis phylogenies? This is a fundamental question that any study involving a phylogenetic framework must face at the beginning of the project. Building a phylogeny from gene sequence data (purpose-built phylogeny) requires more effort, expertise, and cost than subsetting an already available phylogeny (synthesis-based phylogeny). However, we still lack a comparison of how these two approaches to building phylogenetic trees influence common community phylogenetic analyses such as comparing community phylogenetic diversity and estimating trait phylogenetic signal. Here, we generated three purpose-built phylogenies and their corresponding synthesis-based trees (two from Phylomatic and one from the Open Tree of Life, OTL). We simulated 1,000 communities and 12,000 continuous traits along each purpose-built phylogeny. We then compared the effects of different trees on estimates of phylogenetic diversity (alpha and beta) and phylogenetic signal (Pagel's λ and Blomberg's K). Synthesis-based phylogenies generally yielded higher estimates of phylogenetic diversity when compared to purpose-built phylogenies. However, resulting measures of phylogenetic diversity from both types of phylogenies were highly correlated (Spearman's ρ  > 0.8 in most cases). Mean pairwise distance (both alpha and beta) is the index that is most robust to the differences in tree construction that we tested. Measures of phylogenetic diversity based on the OTL showed the highest correlation with measures based on the purpose-built phylogenies. Trait phylogenetic signal estimated with synthesis-based phylogenies, especially from the OTL, was also highly correlated with estimates of Blomberg's K or close to Pagel's λ from purpose-built phylogenies when traits were simulated under Brownian motion. For commonly employed community phylogenetic analyses, our results justify taking advantage of recently developed and continuously improving synthesis trees, especially the Open Tree of Life.

The use of DNA barcodes to estimate phylogenetic diversity in forest communities of southern China

To elucidate potential ecological and evolutionary processes associated with the assembly of plant communities, there is now widespread use of estimates of phylogenetic diversity that are based on a variety of DNA barcode regions and phylogenetic construction methods. However, relatively few studies consider how estimates of phylogenetic diversity may be influenced by single DNA barcodes incorporated into a sequence matrix (conservative regions vs. hypervariable regions) and the use of a backbone family-level phylogeny. Here, we use general linear mixed-effects models to examine the influence of different combinations of core DNA barcodes (rbcL, matK, ITS, and ITS2) and phylogeny construction methods on a series of estimates of community phylogenetic diversity for two subtropical forest plots in Guangdong, southern China. We ask: (a) What are the relative influences of single DNA barcodes on estimates phylogenetic diversity metrics? and (b) What is the effect of using a backbone family-level phylogeny to estimate topology-based phylogenetic diversity metrics? The combination of more than one barcode (i.e., rbcL + matK + ITS) and the use of a backbone family-level phylogeny provided the most parsimonious explanation of variation in estimates of phylogenetic diversity. The use of a backbone family-level phylogeny showed a stronger effect on phylogenetic diversity metrics that are based on tree topology compared to those that are based on branch lengths. In addition, the variation in the estimates of phylogenetic diversity that was explained by the top-rank models ranged from 0.1% to 31% and was dependent on the type of phylogenetic community structure metric. Our study underscores the importance of incorporating a multilocus DNA barcode and the use of a backbone family-level phylogeny to infer phylogenetic diversity, where the type of DNA barcode employed and the phylogenetic construction method used can serve as a significant source of variation in estimates of phylogenetic community structure.

Preserving the tree of life of the fish family Cyprinidae in Africa in the face of the ongoing extinction crisis

Our understanding of how the phylogenetic tree of fishes might be affected by the ongoing extinction risk is poor. This is due to the unavailability of comprehensive DNA data, especially for many African lineages. In addition, the ongoing taxonomic confusion within some lineages, e.g., Cyprinidae, makes it difficult to contribute to the debate on how the fish tree of life might be shaped by extinction. Here, we combine COI sequences and taxonomic information to assemble a fully sampled phylogeny of the African Cyprinidae and investigate whether we might lose more phylogenetic diversity (PD) than expected if currently threatened species go extinct. We found evidence for phylogenetic signal in extinction risk, suggesting that some lineages might be at higher risk than others. Based on simulated extinctions, we found that the loss of all threatened species, which approximates 37% of total PD, would lead to a greater loss of PD than expected, although highly evolutionarily distinct species are not particularly at risk. Pending the reconstruction of an improved multi-gene phylogeny, our results suggest that prioritizing high-EDGE species (evolutionary distinct and globally endangered species) in conservation programmes, particularly in some geographic regions, would contribute significantly to safeguarding the tree of life of the African Cyprinidae.

Flower-visitor communities of an arcto-alpine plant-Global patterns in species richness, phylogenetic diversity and ecological functioning

Pollination is an ecosystem function of global importance. Yet, who visits the flower of specific plants, how the composition of these visitors varies in space and time and how such variation translates into pollination services are hard to establish. The use of DNA barcodes allows us to address ecological patterns involving thousands of taxa that are difficult to identify. To clarify the regional variation in the visitor community of a widespread flower resource, we compared the composition of the arthropod community visiting species in the genus Dryas (mountain avens, family Rosaceae), throughout Arctic and high-alpine areas. At each of 15 sites, we sampled Dryas visitors with 100 sticky flower mimics and identified specimens to Barcode Index Numbers (BINs) using a partial sequence of the mitochondrial COI gene. As a measure of ecosystem functioning, we quantified variation in the seed set of Dryas. To test for an association between phylogenetic and functional diversity, we characterized the structure of local visitor communities with both taxonomic and phylogenetic descriptors. In total, we detected 1,360 different BINs, dominated by Diptera and Hymenoptera. The richness of visitors at each site appeared to be driven by local temperature and precipitation. Phylogeographic structure seemed reflective of geological history and mirrored trans-Arctic patterns detected in plants. Seed set success varied widely among sites, with little variation attributable to pollinator species richness. This pattern suggests idiosyncratic associations, with function dominated by few and potentially different taxa at each site. Taken together, our findings illustrate the role of post-glacial history in the assembly of flower-visitor communities in the Arctic and offer insights for understanding how diversity translates into ecosystem functioning.

Expedited assessment of terrestrial arthropod diversity by coupling Malaise traps with DNA barcoding 1

Monitoring changes in terrestrial arthropod communities over space and time requires a dramatic increase in the speed and accuracy of processing samples that cannot be achieved with morphological approaches. The combination of DNA barcoding and Malaise traps allows expedited, comprehensive inventories of species abundance whose cost will rapidly decline as high-throughput sequencing technologies advance. Aside from detailing protocols from specimen sorting to data release, this paper describes their use in a survey of arthropod diversity in a national park that examined 21 194 specimens representing 2255 species. These protocols can support arthropod monitoring programs at regional, national, and continental scales.

Trait and phylogenetic diversity provide insights into community assembly of reef-associated shrimps (Palaemonidae) at different spatial scales across the Chagos Archipelago

Coral reefs are the most biodiverse marine ecosystem and one of the most threatened by global climate change impacts. The vast majority of diversity on reefs is comprised of small invertebrates that live within the reef structure, termed the cryptofauna. This component of biodiversity is hugely understudied, and many species remain undescribed. This study represents a rare analysis of assembly processes structuring a distinct group of cryptofauna, the Palaemonidae, in the Chagos Archipelago, a reef ecosystem under minimal direct human impacts in the central Indian Ocean. The Palaemonidae are a diverse group of Caridae (infraorder of shrimps) that inhabit many different niches on coral reefs and are of particular interest because of their varied habitat associations. Phylogenetic and trait diversity and phylogenetic signal were used to infer likely drivers of community structure. The mechanisms driving palaemonid community assembly and maintenance in the Chagos Archipelago showed distinct spatial patterns. At local scales, among coral colonies and among reefs fringing individual atolls, significant trait, and phylogenetic clustering patterns suggest environmental filtering may be a dominant ecological process driving Palaemonidae community structure, although local competition through equalizing mechanisms may also play a role in shaping the local community structure. Importantly, we also tested the robustness of phylogenetic diversity to changes in evolutionary information as multi‐gene phylogenies are resource intensive and for large families, such as the Palaemonidae, are often incomplete. These tests demonstrated a very modest impact on phylogenetic community structure, with only one of the four genes (PEPCK gene) in the phylogeny affecting phylogenetic diversity patterns, which provides useful information for future studies on large families with incomplete phylogenies. These findings contribute to our limited knowledge of this component of biodiversity in a marine locality as close to undisturbed by humans as can be found. It also provides a rare evaluation of phylogenetic diversity methods.

Taxon sampling effects on the quantification and comparison of community phylogenetic diversity

Ecologists are increasingly making use of molecular phylogenies, especially in the fields of community ecology and conservation. However, these phylogenies are often used without full appreciation of their underlying assumptions and uncertainties. A frequent practice in ecological studies is inferring a phylogeny with molecular data from taxa only within the community of interest. These "inferred community phylogenies" are inherently biased in their taxon sampling. Despite the importance of comprehensive sampling in constructing phylogenies, the implications of using inferred community phylogenies in ecological studies have not been examined. Here, we evaluate how taxon sampling affects the quantification and comparison of community phylogenetic diversity using both simulated and empirical data sets. We demonstrate that inferred community trees greatly underestimate phylogenetic diversity and that the probability of incorrectly ranking community diversity can reach up to 25%, depending on the dating methods employed. We argue that to reach reliable conclusions, ecological studies must improve their taxon sampling and generate the best phylogeny possible.

Host Specificity in Subarctic Aphids

Plants and herbivorous (or parasitic) insects form the majority of macroscopic life. The specificity of interaction between host plant and parasitic insect depends on the adaptations of both the host and the parasite. Over time, these interactions evolve and change as a result of an 'arms race' between host and parasite, and the resulting species-specific adaptations may be maintained, perpetuating these interactions across speciation events. This can lead to specialisation between species or clades. With speciation and species sorting over time, complex interactions evolve. Here, we elucidate a three-tier method to test these interactions using the aphids (Hemiptera: Aphididae) and plants of Churchill (Manitoba, Canada) as a model system. We analyzed these interactions by testing for three patterns in host specificity: monophagy, phylogenetic clustering, and cophylogeny. We defined monophagy strictly as one species feeding exclusively upon a single host plant species (an association likely driven by arms races in morphology, chemical resistance/tolerance, and visual appearance) and observed this in 7 of 22 aphid species. In all the remaining 'polyphagous' cases, there was a strong trend toward monophagy (80% of individuals were found on a single host plant species). Second, we observed two separate examples of phylogenetic clustering where groups of closely related aphid species fed upon individual plant species. Finally, we found no support for cophylogenetic relationships where both aphids and plants cospeciate to form congruent phylogenetic trees (evidence of coadaptation through an ongoing arms race). One explanation for uncovering species-specific interactions in a recently deglaciated, subarctic locality is that the species involved in the associations moved north together. Testing different levels of specificity in the most predominant species-species interactions on the planet will allow us to elucidate these patterns accurately and gives us insight into where to direct future research.

DNA Barcoding of an Assembly of Montane Andean Butterflies (Satyrinae): Geographical Scale and Identification Performance

DNA barcoding is a technique used primarily for the documentation and identification of biological diversity based on mitochondrial DNA sequences. Butterflies have received particular attention in DNA barcoding studies, although varied performance may be obtained due to different scales of geographic sampling and speciation processes in various groups. The montane Andean Satyrinae constitutes a challenging study group for taxonomy. The group displays high richness, with more of 550 species, and remarkable morphological similarity among taxa, which renders their identification difficult. In the present study, we evaluated the effectiveness of DNA barcodes in the identification of montane Andean satyrines and the effect of increased geographical scale of sampling on identification performance. Mitochondrial sequences were obtained from 104 specimens of 39 species and 16 genera, collected in a forest remnant in the northwest Andes. DNA barcoding has proved to be a useful tool for the identification of the specimens, with a well-defined gap and producing clusters with unambiguous identifications for all the morphospecies in the study area. The expansion of the geographical scale with published data increased genetic distances within species and reduced those among species, but did not generally reduce the success of specimen identification. Only in Forsterinaria rustica (Butler, 1868), a taxon with high intraspecific variation, the barcode gap was lost and low support for monophyly was obtained. Likewise, expanded sampling resulted in a substantial increase in the intraspecific distance in Morpho sulkowskyi (Kollar, 1850); Panyapedaliodes drymaea (Hewitson, 1858); Lymanopoda obsoleta (Westwood, 1851); and Lymanopoda labda Hewitson, 1861; but for these species, the barcode gap was maintained. These divergent lineages are nonetheless worth a detailed study of external and genitalic morphology variation, as well as ecological features, in order to determine the potential existence of cryptic species. Even including these cases, DNA barcoding performance in specimen identification was 100% successful based on monophyly, an unexpected result in such a taxonomically complicated group.

Forecasting pollination declines through DNA barcoding: the potential contributions of macroecological and macroevolutionary scales of inquiry

While pollinators are widely acknowledged as important contributors to seed production in plant communities, we do not yet have a good understanding of the importance of pollinator specialists for this ecosystem service. Determination of the prevalence of pollinator specialists is often hindered by the occurrence of cryptic species and the limitations of observational data on pollinator visitation rates, two areas where DNA barcoding of pollinators and pollen can be useful. Further, the demonstrated adequacy of pollen DNA barcoding from historical records offers opportunities to observe the effects of pollinator loss over longer timescales, and phylogenetic approaches can elucidate the historical rates of extinction of specialist lineages. In this Viewpoint article, we review how advances in DNA barcoding and metabarcoding of plants and pollinators have brought important developments to our understanding of specialization in plant-pollinator interactions. We then put forth several lines of inquiry that we feel are especially promising for providing insight on changes in plant-pollinator interactions over space and time. Obtaining estimates of the effects of reductions in specialists will contribute to forecasting the loss of ecosystem services that will accompany the erosion of plant and pollinator diversity.