Skip to main content

Main menu

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
  • Submit
  • About
    • Editorial Board
    • PNAS Staff
    • FAQ
    • Accessibility Statement
    • Rights and Permissions
    • Site Map
  • Contact
  • Journal Club
  • Subscribe
    • Subscription Rates
    • Subscriptions FAQ
    • Open Access
    • Recommend PNAS to Your Librarian

User menu

  • Log in
  • My Cart

Search

  • Advanced search
Home
Home
  • Log in
  • My Cart

Advanced Search

  • Home
  • Articles
    • Current
    • Special Feature Articles - Most Recent
    • Special Features
    • Colloquia
    • Collected Articles
    • PNAS Classics
    • List of Issues
  • Front Matter
    • Front Matter Portal
    • Journal Club
  • News
    • For the Press
    • This Week In PNAS
    • PNAS in the News
  • Podcasts
  • Authors
    • Information for Authors
    • Editorial and Journal Policies
    • Submission Procedures
    • Fees and Licenses
  • Submit
Research Article

Experimental dispersal reveals characteristic scales of biodiversity in a natural landscape

Rachel M. Germain, Sharon Y. Strauss, and Benjamin Gilbert
  1. aDepartment of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada M5S 3G5;
  2. bDepartment of Evolution and Ecology, University of California, Davis, CA 95616

See allHide authors and affiliations

PNAS April 25, 2017 114 (17) 4447-4452; first published April 17, 2017; https://doi.org/10.1073/pnas.1615338114
Rachel M. Germain
aDepartment of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada M5S 3G5;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: rgermain@zoology.ubc.ca
Sharon Y. Strauss
bDepartment of Evolution and Ecology, University of California, Davis, CA 95616
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Benjamin Gilbert
aDepartment of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, Canada M5S 3G5;
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  1. Edited by David Tilman, University of Minnesota, St. Paul, MN, and approved March 28, 2017 (received for review September 14, 2016)

  • Article
  • Figures & SI
  • Info & Metrics
  • PDF
Loading

Significance

Biological communities differ dramatically in numbers and identities of species, a pattern that could be explained by many mechanisms that each vary with spatial scale. Testing how ecological mechanisms transition among scales is key to understanding the maintenance of diversity but is infeasible in most systems. In a natural plant community, we experimentally enhanced seed dispersal over scales ranging from 1 m to 10 km using a seed vacuum. Our results indicate that pervasive dispersal limitation constrains local communities at and above 100 m, causing communities to contain half as many species as the environment is capable of supporting. Our results suggest that many species at risk for extirpation in this global biodiversity hotspot could be restored through managed dispersal.

Abstract

Ecological theory posits that dispersal among habitat patches links local communities and is a key “regional” process that maintains biological diversity. However, manipulations required to experimentally test regional processes are infeasible for most systems, and thus more work is needed to detect the scales at which regional processes manifest and their overall effect on diversity. In a Californian grassland, a hotspot for global biodiversity, we used a seed vacuum to increase dispersal at spatial scales varying from 1 m to 10 km while maintaining a realistic spatial structure of species pools and environmental conditions. We found that dispersal limitation has a profound influence on diversity; species richness increased with the spatial scale of seed mixing, doubling in plots that received seed from large (≥5 km) compared with small (≤5 m) scales. This increase in diversity corresponded to an increase in how well species distributions were explained by environmental conditions, from modest at small scales (R2 = 0.34) to strong at large scales (R2 = 0.52). Responses to the spatial scale of seed mixing were nonlinear, with no differences below 5 m or above 5 km. Nonlinearities were explained by homogeneity of environmental conditions below 5 m and by a lack of additional variation in the species pool above 5 km. Our approach of manipulating natural communities at different spatial scales reveals (i) nonlinear transitions in the importance of environmental sorting and dispersal, and (ii) the negative effects of dispersal limitation on local diversity, consistent with previous research suggesting that large numbers of species are headed toward regional extinction.

  • dispersal limitation
  • McLaughlin Natural Reserve
  • metacommunity
  • seed addition
  • spatial scale

The problem of pattern and scale is the central problem in ecology, unifying population biology and ecosystem science, and marrying basic and applied ecology.

S. A. Levin (1992)

The processes that structure ecological populations, biodiversity, and ecosystem properties transition in importance across spatial scales (1, 2). As the spatial scale of observation increases, the range of environments sampled (1, 3) and the geographic distance separating localities (4) become increasingly important in shaping species distributions. As a consequence, the relative spatial scaling of different ecological processes is thought to underlie some of the most important patterns in ecology, such as species–area (5) and biodiversity–ecosystem function relationships (6). Because identifying the important scales is challenging, ecologists often compare ecological patterns among local and regional scales to simplify theoretical (7⇓–9) and empirical research (10, 11). However, how closely local and regional delineations match up with the actual scaling of ecological processes is rarely known (2, 12). Quantifying the spatial scaling of these processes promises to enrich our understanding of the mechanisms that maintain diversity, yet remains elusive even in biodiversity hotspots that require this information for conservation decisions (13).

A major challenge to testing how ecological processes transition among spatial scales in natural communities is that regional processes, unlike local processes, are not often amenable to experimental manipulation in the field due to the inability to move most communities of organisms. Manipulative tests of local and regional diversity are typically performed in mesocosm experiments on communities constructed using simplified environments, species pools, or dispersal patterns (14). Because it is often unclear how such simplifications affect experimental outcomes (15), mesocosm experiments allow essential tests of the range of potential outcomes under different sets of experimental conditions, but cannot capture the importance of processes that occur in nature. Finding new ways to combine the power of manipulative field experiments with the biological realism of natural landscape structure is necessary to address fundamental questions, such as: how widespread is dispersal limitation in a community and at what spatial scales does it manifest? How strongly do species sort along environmental gradients in the absence of dispersal limitation? And, even more basically, at what spatial scales are local and regional communities most appropriately defined?

The above questions can be tested by experimentally removing dispersal limitation through the homogenization of species pools, that is, by redistributing species equally across habitat patches while maintaining their abundance distributions. Theory makes distinct predictions for how such homogenization would affect species diversity, and whether it would increase or decrease the probability that species are found in environmentally suitable habitat patches (9, 16, 17). This latter concept, referred to here as “species–environment associations,” is commonly measured by the amount of variation in species distributions that is explained by environmental covariates (18). If dispersal chronically limits species movement and the landscape is patchy and heterogeneous, then homogenization of species pools would increase both species richness (19, 20) and the strength of species–environment associations (16) (SI Appendix, Fig. S1A). These changes occur as species access suitable localities that were previously inaccessible (9), causing the average number of species in a site to increase and the variation in species distributions explained by the environment to similarly increase.

Alternatively, if species are not dispersal-limited and already well-matched to environmental conditions (SI Appendix, Fig. S1B), the strength of species–environment associations should decrease even as local richness increases due to mass effects (16, 17), meaning that species are increasingly found in environments to which they are not suited. When species–environment associations are neutral with respect to species identity, species–environment associations should not change whether species have limited dispersal (SI Appendix, Fig. S1C; increased species richness) or are already dispersing evenly among communities (SI Appendix, Fig. S1D; no change in species richness) (8, 16). In the metacommunity ecology literature, these four potential outcomes of increasing dispersal would correspond to species sorting, mass effects, and neutral dynamics (with or without dispersal limitation), respectively (9, 16). The approach of homogenizing species pools and tracking diversity responses to local conditions is similar to that taken by many experiments that use mesocosm (e.g., ref. 11) or field (e.g., refs. 19⇓–21) manipulations with artificial species pools (14, 15, 22), but has yet to be implemented in the field with natural spatial structure in species pools, abundance distributions, and environments.

We tested the spatial scales at which dispersal and environmental heterogeneity impose constraints on species distributions and diversity using a “hay transfer” technique (23) in serpentine annual plant communities in California. This technique involved vacuuming the seed bank and other loose material from field plots, pooling the material among plots, and then redistributing it. Our “species-pooling” treatments captured five spatial scales (∼1 m, 5 m, 100 m, 5 km, and 10 km; SI Appendix, Fig. S2), with the 100-m treatment receiving a mix of seeds from all plots within 100 m, for example. In the annual plant communities studied, pooling the seed bank redistributes all plants during their dormant stage. In the following growing season, we collected data on species occupancy and environmental conditions to answer four questions: (i) How common is dispersal limitation in a metacommunity? (ii) How strongly do species sort along environmental gradients once potential dispersal limitation is experimentally removed? (iii) At what spatial scales do signals of dispersal limitation and environmental sorting manifest? (iv) Do these signals correspond to the natural structure of species pools and environmental conditions? Our surveys were conducted at peak flowering, after the filtering effects of competition and environment, which can be strongest at the seedling stage (24, 25), had time to take effect.

Our serpentine study system is ideal for testing questions of the spatial scaling of diversity for two reasons. First, the high occurrence of annual species and patchy distribution of serpentine habitat among a nonserpentine matrix is well-suited to testing spatial questions (26) with the hay transfer method (23). Second, our study took place in the California Floristic Province, a global biodiversity hotspot that is threatened by the pervasive invasion of European grasses (26, 27). Serpentine soils are hypothesized to act as spatial refugia for native species to escape the direct competitive effects of invasive grasses (28), but less often considered are the indirect negative effects of invaders on diversity through the isolation of habitat patches (13, 29).

Given that invasive grasses and other human impacts [road building, fire suppression (30)] have increased fragmentation of native plant communities in the past 200 y, we hypothesized that diversity is constrained by limited dispersal and that, if this dispersal constraint were removed, species would sort deterministically according to their environmental niche requirements (i.e., increased species–environment associations; SI Appendix, Fig. S1A). Environmental conditions were described by two composite variables—principal components analysis (PCA) axes—that summarized soil chemistry, fertility, moisture, and site topography. The strength of species–environment associations was the pseudo-R2 of logistic regressions that predicted the occurrences of all species with environmental conditions (PCA axes) used as predictors (see Data analysis in Materials and Methods and SI Appendix, Fig. S4 A–C for an example of data that underlie pseudo-R2 results). The degree of dispersal limitation, and the scales at which its effects manifest, are important to understanding biodiversity in this global hotspot and to projecting community stability over longer timescales (13).

Results

Increasing the spatial scale of species pooling caused a net increase of 12 species in plots (F4,112 = 69.4, P < 0.001; Fig. 1A). This increase was sigmoidal rather than linear, with no significant differences in species richness between the two smallest (1 m and 5 m; P = 0.72) or the two largest (5 km and 10 km; P = 0.34) spatial scales, but highly significant (all P < 0.001) differences among small (≤5 m), intermediate (100 m), and large (≥5 km) spatial scales. As a result, species richness in our plots doubled (F4,112 = 69.9, P < 0.001) from 10 species in the 1-m treatment to 21 species in the 10-km treatment (SI Appendix, Fig. S3A; all measures are for 0.75- × 0.75-m plots). Species richness postmanipulation was composed of species initially present at sites (not dispersal limited) plus those gained by species pooling (dispersal limited; SI Appendix, Fig. S3B).

Fig. 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 1.

Effect of spatial scale of species pooling on (A) net change in species richness (mean ± SE), relative to an unmanipulated control plot (C3 in SI Appendix, Fig. S2B) and (B) the strength of species–environment associations across a community of species (pseudo-R2 values); observed patterns correspond to prediction in SI Appendix, Fig. S1A. Points with the same letter were not significantly different in a multiple comparisons test. Net changes in species richness are driven by new species gained, which translates into a doubling in species richness when species are pooled at the 1-m versus at the 10-km scale (SI Appendix, Fig. S3).

To explore whether this increase in species richness corresponded with an increase or decrease in species–environment associations, we tested whether environmental conditions predicted the occurrences of all 73 annual species observed in our study. Environmental conditions were highly correlated and thus were summarized by the first and second axis scores of a principal components analysis (Materials and Methods). Axis 1 primarily summarized soil chemistry (i.e., Ca/Mg, Olsen-P, X-K, NO3-N, and organic matter), whereas axis 2 summarized elevation and soil moisture (SI Appendix, Fig. S5B). We found the proportion of variance explained (pseudo-R2) by our models increased from 0.34 to 0.52 with increasing spatial scale of species pooling (Fig. 1B); the sigmoidal response of species richness to the spatial scale of species pooling was closely mirrored by the sigmoidal increase in variation explained by the environment (Fig. 1). Supplementary analyses of species–environment associations with percentage of cover data, rather than presence/absence data, produced broadly similar increases with increasing spatial scale of species pooling (SI Appendix, Fig. S9 and discussion in SI Appendix, Supplementary Methods, Results, and Discussion).

In addition to plot-level species richness (alpha diversity), we tested the effects of our species-pooling treatments on other components of diversity. As with species richness, the composition of species in plots responded sigmoidally to the spatial scale of species pooling (Fig. 2A), as estimated using the axis 1 and 2 scores of a principal coordinates analysis (PCoA) (Materials and Methods). These axes were significantly associated with 42 plant species (P < 0.05), including many endemic serpentine-associated species (e.g., Clarkia gracilis, Collinsia sparsiflora, and Navarretia jepsonii). Beta diversity (i.e., among-plot turnover in species composition, measured with Jaccard’s dissimilarity index) decreased with the scale of species pooling from 0.81 to 0.56 between the 1-m to 10-km spatial scaling treatments (Fig. 2B). The size of the regional species pool (gamma diversity), estimated as the total number of unique species that emerged across plots, was not affected by the species-pooling treatments, with a mean pool size of 58.0 ± 0.5 SE across treatments (x-intercept; SI Appendix, Fig. S6). At a common gamma diversity, alpha and beta diversity are inversely related, as compositional dissimilarity among sites decreases as a consequence of species being found in more sites (31). Site occupancy (the number of sites occupied by each species) also increased with the spatial scale of species pooling (SI Appendix, Fig. S6), as did community evenness due to the increasing occurrences of regionally rare species (shallower slope in SI Appendix, Fig. S6). Not all 73 species were observed in all treatments; on average, each treatment had 1.2 species (±0.58 SE) that were not observed in any other treatment.

Fig. 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 2.

Effect of spatial scale of species pooling on (A) species composition and (B) the compositional dissimilarity of plots (beta diversity). Species composition shows the mean and SE of the first axis scores from a PCoA with Jaccard’s dissimilarity index (see biplot in SI Appendix, Fig. S5A); results with second-axis scores are qualitatively similar. Points with the same letter were not significantly different in a multiple comparisons test. Beta diversity patterns mirror changes in species richness, as increasing local (alpha) diversity with no change in regional (gamma) diversity reduces among-plot compositional dissimilarity.

Differences in the naturally occurring spatial structure of environmental conditions and species pools provide further evidence of how regional processes structure local patterns of diversity. We first tested for differences in environmental conditions among sites at spatial scales corresponding to our experimental setup (SI Appendix, Fig. S2) to characterize the spatial scaling of the environment. Environmental conditions differed significantly among groups of sites within halves (100-km scale; F = 14.05, P = 0.001), but differences between halves of the reserve (5-km scale; F = 247.68, P = 0.001) were larger, with limited overlap in environmental conditions (Fig. 3A). Differences among reserve halves were driven by axis 2 (% soil moisture, elevation) but not axis 1 (soil chemistry, fertility). For those environmental variables that were quantified with multiple within-site measurements (5-m scale), there were no statistically significant differences among plots that occurred at the same site (all P > 0.3; SI Appendix, Table S2). Despite these average changes in environmental conditions with distance, even distant sites (≥10 km apart) frequently shared similar environmental conditions (Fig. 3A and SI Appendix, Figs. S7 and S8), supporting the finding that dispersal limitation at these scales can have important consequences.

Fig. 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Fig. 3.

Natural differences among halves of the reserve in (A) environmental conditions (from a PCA) and (B) species composition (from a PCoA). The contours outline the bivariate normal distributions of each half of the reserve; dashed contours/upward triangles are sites that occur in the northern half, and solid contours/downward triangles are sites that occur in the southern half. Environmental conditions (P = 0.001) but not species (P = 0.175) significantly differed among halves (Results).

In contrast to environmental conditions, the composition of species significantly differed at the 100-m scale (F = 1.96, P = 0.001) but not between halves of the reserve (5-km scale; F = 1.26, P = 0.168). Despite the lack of a significant difference in species composition, the two halves of the reserve had only 62% (37/60) of species in common, with 20% (12/60) of species unique to the northern half and 18% (11/60) unique to the southern half; these numbers were calculated from the 1-m treatment and sum to less than 73 because not all species were observed in all treatments. The species unique to one-half of the reserve were significantly less associated with the axis 1 and 2 PCoA scores than species that were common to both halves (P = 0.023), which potentially reconciles the high compositional similarity among halves of the reserve despite the presence of unique species.

Discussion

Understanding the scale dependency of ecological processes is central to protecting biodiversity and a core interest of ecology research (1, 2, 32-34), yet experimental demonstrations have remained elusive for natural communities (2). By experimentally manipulating the spatial scale of seed dispersal in a biodiversity hotspot, our study resulted in two key findings: (i) compared with natural patterns, species richness and the strength of environmental matching doubled when seeds were distributed at large spatial scales, and (ii) nonlinear diversity responses to the spatial scale of species pooling suggests that there are distinct spatial transitions in the forces that structure diversity. These findings provide insight into how metacommunities and the distributions of their constituent species are structured, with implications for conservation.

In the absence of species pooling, species in our study showed modest associations with environmental conditions (blue point in Fig. 1B), consistent with the findings of several observational studies (e.g., refs. 34, 35). However, this type of observational evidence alone cannot identify the mechanisms underlying the pervasive lack of environmental associations, such as dispersal limitation, neutrality, and mass effects (SI Appendix, Fig. S1). The experimental enhancement of dispersal via seed pool manipulations allowed us to discriminate among potential mechanisms and supports a model in which dispersal limitation precludes the colonization of otherwise suitable sites (SI Appendix, Fig. S1A). Consistent with these results, a subsequent analysis showed that species with more effective dispersal modes (wind and vertebrate dispersed) tended to occupy more sites in the absence of species pooling (SI Appendix, Fig. S10 and Supplementary Methods, Results, and Discussion).

Similar evidence of the role of dispersal limitation has been obtained for individual species, such as the canary violet (4); we extend this perspective beyond individual species to an entire plant community and show a near doubling in the strength of species–environment associations that explains over half of the variation in how communities of species are distributed (Fig. 1B). An additional feature of our community-level seed additions is that they allowed natural spatial structure in both species composition and relative abundances, rather than fixed amounts of seed added for all sites and all species (19). Doing so most closely reflects natural variation in propagule pressure, an important determinant of dispersal limitation and community structure (20, 36). Although we cannot refute the possibility that the observed responses are transient, evidence from related seed addition research has shown that initial diversity effects are persistent, even in competitive environments (37), and increases in species–environment associations would not be expected in a transient scenario. Additionally, our results were insensitive to the exclusion of species with small populations, which are those most likely to be “sink” populations (SI Appendix, Table S3); we discuss these points at length in SI Appendix, Supplementary Methods, Results, and Discussion.

To what extent differences in species composition among local communities can be explained by distance or environment has been a major area of ecological research (38, 39), stimulated by the development of statistical methods to partition their relative influences from observational data (18, 40). However, variance-partitioning methods can produce biased estimates of spatial and environmental components of variation, as well as inflated model fits (41), and should be interpreted as rough estimates. Moreover, partitioning approaches cannot isolate the interactive effects of dispersal limitation and environmental components (18), yet it is these interactive effects that determine the persistence of species within landscapes (13, 29). In our study, we find that environmental conditions explain 32% of variation in species occupancy patterns in the absence of species pooling and an additional 20% when distance effects are experimentally removed (Fig. 1B). This increase when distance effects are removed quantifies the interactive effect of distance and environmental conditions, or the degree to which dispersal limitation prevents species from accessing suitable habitats.

Nearly all components of diversity that were examined showed consistent, nonlinear responses to the spatial scale of species pooling, including species richness (alpha diversity), species–environment associations, and compositional turnover (beta diversity). At small spatial scales (<5 m), conditions are homogenous, and most species are able to disperse these distances in one to few generations (42); our species-pooling manipulations had negligible effects at this scale (blue and green points, Figs. 1 and 2 and SI Appendix, Fig. S3). At intermediate (100 m) and larger (≥5 km) scales, our manipulations allowed seed to cross the nonserpentine grassland matrix to reach distant serpentine sites, resulting in increased species richness and species–environment associations consistent with dispersal limitation. Although environmental conditions also tended to be more different as distance increased, many sites were similar even when separated by ∼10 km (Fig. 3 and SI Appendix, Figs. S7 and S8), making dispersal across these distances important for accessing distant sites with similar environments.

Together, our evidence reveals three “characteristic scales” [sensu Levin (2)] or “domains of scale” [sensu Wiens (1)] among which strengths of ecological processes appear to differ. When seed pooling occurs at local scales ≤5 m, communities experience relatively homogenous environmental conditions, but the make-up of communities is severely limited by dispersal, such that species–environment associations frequently fail to emerge. As the scale of species pooling increases to 100 m and beyond, the increased environmental difference in new sites appears to be outweighed by the large numbers of similar sites (SI Appendix, Figs. S7 and S8) and the release from dispersal limitation. The net result is that species are better able to exploit the increased range of environmental heterogeneity, increasing species–environment associations. However, as the scale of species pooling increases beyond 5 km, dispersal limitation no longer alters community responses to continually increasing environmental heterogeneity.

The lack of differences among 5-km and 10-km scales was surprising given their differences in spatial extent and environmental conditions, particularly soil moisture and elevation (as summarized by PCA axis 2, Fig. 3A); there are two potential, nonmutually exclusive explanations for this trend. First, the high similarity in the composition of species pool among halves of the reserve (the 5-km scale, Fig. 3B) indicates that mixing seed among halves (10-km scale) would introduce little variation to the species pool, consistent with the lack of responses that we observed. Although we did observe a number of species that were unique to each half, these species were so rare that they did not contribute strongly to compositional variation among sites or the strengthening of species–environment associations in response to seed mixing. Second, it is possible that most species in the reserve are more specialized to soil fertility and biochemistry than to differences in elevation and soil moisture (Fig. 3, axis 1 and 2, respectively). Because the reserve halves did not significantly differ in axis 1 scores, specialization along this axis did not generate among-half differences in species composition.

Our finding that local communities are half as speciose as the environment is capable of supporting sheds light on the past and potential future of this global biodiversity hotspot. Californian landscapes no longer resemble those of 200 y ago; European invaders now form a matrix of unsuitable habitat (28) and have restricted native species to harsh refuge environments, such as serpentine (13, 26). Our findings of low species richness and high regional rarity (SI Appendix, Fig. S6) are consistent with those predicted for areas experiencing regional “extinction debts”—the delayed extinction of species from a region through chronic reductions in colonization rates (13). The prognosis for native diversity under this scenario is bleak because old invaders are adapting to tolerate harsh serpentine conditions (43) and new invaders that thrive in these harsh environments are spreading rapidly [e.g., barbed goatgrass (Aegilops triuncialis)]. The net effect is a shrinking and loss of native plant refuge patches that threaten the regional stability of native plant diversity.

In this scenario and others like it, a more regionally focused approach to conservation is needed to preserve many types of local communities and their constituent species, understanding that local persistence relies on colonization from adjacent patches (44). Management at the wrong scale may miss important covariation in species pools and environments, generating mismatches between species and local conditions (1, 4). In our system, management at a spatial extent of 5 km is likely to be most effective, particularly through assisted dispersal or creation of habitat patches that act as stepping stones between otherwise isolated patches. Although the exact scale of spatial delineations will differ among study systems depending on landscape structure, habitat matrix permeability, and dispersal distances of focal organisms, we have offered a clear example of how these critical scales can be identified. In other systems, approaches such as genetic analyses may be required to better understand characteristic scales of dispersal before assessing scales of environmental turnover.

The forces that dictate how species are distributed across landscapes have long fascinated ecologists. In the serpentine system, we have identified scale-specific processes that structure plant communities and how this spatial scaling is explained by the accumulation of species and environments across space. These results provide insights that are specific to an area that is considered a model system for biodiversity research (45) and more generally provide an experimental test of predictions from a well-developed body of theory (e.g., ref. 9). Broadly applying such approaches can better match conservation actions with ecological processes and promises to advance our understanding of one of the longest-standing challenges in ecology.

Materials and Methods

Study System.

The field experiment was conducted at the 2,800-ha McLaughlin Natural Reserve (nrs.ucdavis.edu/McL/) in northern California (38.8739° N, 122.4317° W). The region has a Mediterranean climate, featuring cool, wet winters (November–March) and hot, dry summers (April–October) with ∼750 mm of annual rainfall. The landscape is largely composed of chaparral, oak woodland, and grassland meadow habitat; common herbivores include mule deer, jack rabbits, and pocket gophers. The reserve lies on the San Andreas Fault and has unique soil chemistry owing to the emergence and erosion of the Earth’s mantle into the serpentine ultramafic soils that characterize the region. Serpentine soils have Ca/Mg ratios <1, as well as low levels of essential nutrients, high heavy metal content, and poor soil moisture retention.

We focused specifically on serpentine meadow habitat and observed 113 species in our plots from a potential pool of the 310 species that occur in all habitat types at the reserve; 73 of the 113 species were annual plants (SI Appendix, Table S1). The four most common species observed in unmanipulated plots were Vulpia microstachys (89% of sites), Hemizonia congesta (79%), Plantago erecta (69%), and Lasthenia californica (52%). Sixteen of the 113 total species could not be identified or classified as having an annual or perennial life history and were thus excluded from all analyses; however, these individuals occurred only at single sites and are thus unlikely to have large effects on diversity patterns.

Experimental Setup.

In May 2013, we surveyed the reserve for 30 serpentine meadow sites. Site locations were chosen such that they could be hierarchically grouped at five spatial scales, which resulted in three groups of five sites each at the northwestern and southeastern ends of the reserve (SI Appendix, Fig. S2). Within each site, we flagged and GPS-located eight 0.75- × 0.75-m plots arranged in a 2 × 4 block of plots, with plots separated by 1 m of bare ground. Each plot in a block was randomly assigned to receive a different experimental treatment (discussed below). In total, there were 240 plots (eight plots × 30 sites) and five spatial scales for comparison [at the level of the plot (1 m), a block of plots at a site (5 m), a group of sites (100 m), a reserve half (5 km), and a whole reserve (10 km)] (SI Appendix, Fig. S2). In late July 2013, after all winter annual species had senesced and the majority of summer annual species had set seed, we harvested all seed and standing vegetation from seven of the eight plots at each site using garden shears and a powerful gas-powered leaf vacuum (Stihl BG86); the eighth plot (C2; SI Appendix, Fig. S1B) was left unaltered to evaluate any unintended effects of the vacuuming procedure on plant diversity. All collected materials were stored outside in paper bags to allow natural heat stratification until they could be processed (<6 wk).

All collected materials were reintroduced to the seven vacuumed plots per site in two control and five species-pooling treatments (SI Appendix, Fig. S2B). The control treatments were as follows: a “vacuum without replacement” treatment (C1) to identify individuals that were left behind following vacuuming, and a “vacuum without movement” treatment (C3), where the collected material was homogenized at the plot level and redistributed back onto the source plot. The five spatial scaling treatments involved the pooling, homogenization, and redistribution of material collected from a single plot (1 m; blue plot in SI Appendix, Fig. S2B), multiple plots at a site (5 m; green plot), 5 sites of a single group (100 m; yellow plot), 15 sites from the same half of the reserve (5 km; orange plot), and all 30 sites across the entire reserve (10 km; red plot). Because treatments were nested within sites, and plots within sites are highly similar in environments and species composition, all treatments received a similar regional pool in terms of species richness and relative abundances. The redistributed material was secured to each recipient plot with twine.

We surveyed the plot-level community structure and corresponding environmental parameters during peak biomass in the following growing season. Plots were surveyed April 20 to May 2 in 2014, using percentage of cover estimates of each species because small-statured annuals can occur at densities up to 5,500 individuals per square meter (46); additional surveys were conducted later in the season to confirm the identities of late-flowering species. We surveyed the innermost 0.5 × 0.5 m2 of each plot to account for any edge effects in our analyses. We measured plot-level percentage of soil moisture content, understory photosynthetically active radiation (PAR) in full sun, and slope inclination, as well as site-level elevation, slope aspect, hillside slope steepness, and soil depth. We also performed site-level soil fertility analyses (NO3-N, Olsen-P, X-K, X-Na, X-Ca, X-Mg, pH, cation exchange capacity, organic matter; University of California at Davis Analytical Lab) on soil samples collected and pooled between four plots per site; X-Ca and X-Mg were converted to a ratio of Ca/Mg. Although site-level estimates of environmental conditions preclude finer-scale estimates, those variables for which we do have plot-level estimates showed no difference among plots that occur at the same site (all P > 0.35; SI Appendix, Table S2). This means that, although site-level estimates preclude estimates of within-site (among-plot) error in most of our environmental variables, all treatments and controls are subject to the same error because each is represented within all sites, and plots within sites do not significantly differ in environmental conditions (or species composition, Fig. 3).

Data Analysis.

A presence/absence matrix was created from percentage of cover estimates of all 73 annual species in 29 of the 30 sites that were sampled; data from one site were lost due to a corrupt data file. We used presence/absence data because it is most appropriate for tracking gains and losses of species in response to manipulations of species pools and is most comparable to other studies (14, 29).

We tested the responses of five components of diversity to our species pool manipulations: net changes in species richness, species richness, regional site occupancy, species composition (PCoA axis 1 and 2 scores), and beta diversity. The first four components were tested using separate linear mixed-effects models in the “lmerTest” R package, each with spatial species-pooling treatment as a fixed factor and plot nested within site as a random factor. The nested random effects account for the nonindependence of our nested experimental design in terms of both error structure and degrees of freedom (47). Because species-pooling treatment emerged as a significant predictor in all analyses, we used post hoc Tukey tests using the “multcomp” package to identify treatment levels that differed significantly from each other. Net changes in species richness were estimated by subtracting the control plot C3 (SI Appendix, Fig. S2B) presence/absence matrix from each species-pooling treatment matrix; doing so allowed us to decompose net changes into gains and losses of individual species (SI Appendix, Fig. S3B). Species composition was the first two axis scores of a PCoA using a Jaccard dissimilarity matrix (R package “vegan”; SI Appendix, Fig. S5A).

The last component, beta diversity, was calculated as a matrix of Jaccard dissimilarity coefficients among all pairwise site combinations and separately for each species-pooling treatment. We then performed a permutational multivariate analyses of variance (“adonis” function in R package vegan) to test differences in Jaccard dissimilarity among treatments, constrained to account for treatments being nested within sites. This analysis is analogous to a univariate analysis of variance, except expanded to handle multivariate response variables. Because the overall model was significant, we ran adonis with a Bonferroni correction on all pairwise treatment combinations to identify treatments that significantly differed.

Before testing the responses of species–environment associations to our species-pooling treatments, we used complementary multivariate methods in the R package vegan to reduce the dimensionality of the environmental variables that were measured. First, we used variance inflation factors (VIFs) to confirm that multicollinearity was low (VIFs < 10) and performed a constrained correspondence analysis (CCA), forward-selecting environmental variables using the R function “ordistep” to identify the subset of environmental variables that significantly influenced species occurrences (all but soil depth). Second, we ran a PCA on the subset of environmental variables identified as meaningful in the CCA, with the axis 1 and 2 scores of the PCA summarizing 51.4% of the variation in the among-site environmental dataset (SI Appendix, Fig. S5B).

We used the first (PCA1) and second (PCA2) PCA axis scores as composite environmental variables to examine how the species-pooling treatments affected the strength of species–environment associations. We performed a generalized linear model, separately for each species-pooling treatment, with species occurrences as the binomially distributed response variable and poly(PCA1,2)*species + poly(PCA2,2)*species as predictor variables. The poly(x,2) function in the R “stats” package calculates orthogonal first- and second-order polynomials of PCA1 and PCA2 to detect linear and quadratic relationships; species present in sites of intermediate environmental conditions but absent at the extremes would be best characterized by quadratic relationships. “Species” was a fixed factor in these analyses to facilitate model convergence and because we were directly interested in the amount of variation explained (pseudo-R2 [1-{residual deviance/null deviance}]) by species-specific differences in responses to environmental conditions (47). We then compared the amount of variation explained by composite environmental conditions among species-pooling treatments (18), with the prediction that pseudo-R2 values should increase with the spatial scale of species pooling. To further examine species-level responses (i.e., how strongly the occurrences of each species were associated with environmental conditions), we performed separate generalized linear models for each species and species-pooling treatment, testing the additive effects of poly(PCA1,2) and poly(PCA2,2). We visualized how the distribution of the pseudo-R2 values of the 73 species shifted across species-pooling treatments. We then visualized the distribution of pseudo-R2 values across species for each of the species-pooling treatments.

To quantify similarity among sites in environments and species that occurred naturally at distinct spatial scales, we performed permutational multivariate analyses of variance (adonis function in R package vegan). The response variables were matrices of Euclidian (for environment) and Jaccard (for species) distances among sites, and each matrix was tested for differences among groups of five sites (100-m scale) nested within each half of the reserve (5-km scale). The species distance matrix was created from plots that received the 1-m treatment only, as this treatment most closely reflected natural unmanipulated species distributions.

The control plots were used to assess the presence of any unintended effects of the hay transfer manipulation that were unrelated to the species-pooling treatment. We found no difference in species richness or composition between the unmanipulated (C2 in SI Appendix, Fig. S2) plots and our 1-m treatment plots (all P > 0.998), indicating that the hay transfer manipulation did not affect local diversity. Similarly, there was no difference among plots that received their own hay back (C3) and the 1-m treatment plots that received hay from a different adjacent plot (all P > 0.807); thus, the removal and replacement of hay did not affect plot diversity. As a result, we report only the results from the 1-m treatment plots. The removal without replacement control plot (C1) had significantly lower species richness than the unmanipulated control plot (C2) and all treatment plots (all P < 0.002), indicating that our seed vacuum was effective.

We incorporated percentage-of-cover data to assess whether our species-pooling manipulations caused small transient “sink” populations to establish, leading to biased diversity estimates. To do this, we revisited the original percentage-of-cover matrix from which the presence/absence matrix was generated and converted any percentage-of-cover entries less than or equal to four cutoff values (0.05, 0.10, 0.25, and 0.5%) to zero. For small-statured annual plants, 0.5% cover roughly translates into 10 individual plants; values greater than 0.5% risks removing small but stable populations. We then converted those four percentage-of-cover matrices to presence/absences and reanalyzed our species richness and species–environment association results (SI Appendix, Table S3).

Acknowledgments

We thank S. Harrison, B. Anacker, C. Kohler, C. Lee, E. Case, S. Copeland, P. Arenas, and L. Sergison for advice and/or assistance in the field; S.Y.S. and R.M.G. laboratory members for project feedback; and the editor, two anonymous reviewers, M. McPeek, and J. Stinchcombe for comments on an earlier version of the manuscript. This work was supported by the National Sciences and Engineering Research Council (NSERC)’s Canada Graduate Scholarship and Michael Smith Foreign Study Supplement, as well as Sigma Xi (R.M.G.), and by the California Agricultural Experiment Station (S.Y.S.).

Footnotes

  • ↵1To whom correspondence should be addressed. Email: rgermain{at}zoology.ubc.ca.
  • Author contributions: R.M.G., S.Y.S., and B.G. designed research; R.M.G. performed research; R.M.G. and B.G. analyzed data; and R.M.G., S.Y.S., and B.G. wrote the paper.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1615338114/-/DCSupplemental.

References

  1. ↵
    1. Wiens JA
    (1989) Spatial scaling in ecology. Funct Ecol 3:385–397.
    .
    OpenUrl
  2. ↵
    1. Levin SA
    (1992) The problem of pattern and scale in ecology: The Robert H. MacArthur Award Lecture. Ecology 73:1943–1967.
    .
    OpenUrlCrossRef
  3. ↵
    1. Bell G,
    2. Lechowicz MJ
    (1991) The ecology and genetics of fitness in forest plants. I. Environmental heterogeneity measured by explant traits. J Ecol 79:663–685.
    .
    OpenUrl
  4. ↵
    1. Pinto SM,
    2. MacDougall AS
    (2010) Dispersal limitation and environmental structure interact to restrict the occupation of optimal habitat. Am Nat 175:675–686.
    .
    OpenUrlCrossRefPubMed
  5. ↵
    1. Rosenzweig ML
    (1995) Species Diversity in Space and Time (Cambridge Univ Press, Cambridge, UK).
    .
  6. ↵
    1. Tonkin JD,
    2. Death RG
    (2013) Scale dependent effects of productivity and disturbance on diversity in streams. Fundam Appl Limnol 182:283–295.
    .
    OpenUrl
  7. ↵
    1. MacArthur RH,
    2. Wilson EO
    (1967) The Theory of Island Biogeography (Princeton Univ Press, Princeton, NJ).
    .
  8. ↵
    1. Hubbell SP
    (2001) The Unified Neutral Theory of Biodiversity and Biogeography (Princeton Univ Press, Princeton, NJ).
    .
  9. ↵
    1. Leibold MA, et al.
    (2004) The metacommunity concept: A framework for multi-scale community ecology. Ecol Lett 7:601–613.
    .
    OpenUrlCrossRef
  10. ↵
    1. Cadotte MW
    (2007) Competition-colonization trade-offs and disturbance effects at multiple scales. Ecology 88:823–829.
    .
    OpenUrlCrossRefPubMed
  11. ↵
    1. Chase JM
    (2010) Stochastic community assembly causes higher biodiversity in more productive environments. Science 328:1388–1391.
    .
    OpenUrlAbstract/FREE Full Text
  12. ↵
    1. Ricklefs RE
    (2008) Disintegration of the ecological community. Am Nat 172:741–750.
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Gilbert B,
    2. Levine JM
    (2013) Plant invasions and extinction debts. Proc Natl Acad Sci USA 110:1744–1749.
    .
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Grainger TN,
    2. Gilbert B
    (2016) Dispersal and diversity in experimental metacommunities: Linking theory and practice. Oikos 125:1213–1223.
    .
    OpenUrl
  15. ↵
    1. Logue JB,
    2. Mouquet N,
    3. Peter H,
    4. Hillebrand H, Metacommunity Working Group
    (2011) Empirical approaches to metacommunities: A review and comparison with theory. Trends Ecol Evol 26:482–491.
    .
    OpenUrlCrossRefPubMed
  16. ↵
    1. Pulliam HR
    (2000) On the relationship between niche and distribution. Ecol Lett 3:349–361.
    .
    OpenUrlCrossRef
  17. ↵
    1. Mouquet N,
    2. Loreau M
    (2003) Community patterns in source-sink metacommunities. Am Nat 162:544–557.
    .
    OpenUrlCrossRefPubMed
  18. ↵
    1. Borcard D,
    2. Legendre P,
    3. Drapeau P
    (1992) Partialling out the spatial component of ecological variation. Ecology 73:1045–1055.
    .
    OpenUrlCrossRef
  19. ↵
    1. Tilman D
    (1997) Community invasibility, recruitment limitation, and grassland biodiversity. Ecology 78:81–92.
    .
    OpenUrlCrossRef
  20. ↵
    1. Myers JA,
    2. Harms KE
    (2009) Seed arrival, ecological filters, and plant species richness: A meta-analysis. Ecol Lett 12:1250–1260.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Germain RM, et al.
    (2013) Spatial variability in plant predation determines the strength of stochastic community assembly. Am Nat 182:169–179.
    .
    OpenUrl
  22. ↵
    1. Cadotte MW
    (2006) Dispersal and species diversity: A meta-analysis. Am Nat 167:913–924.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Coiffait-Gombault C,
    2. Buisson E,
    3. Dutoit T
    (2010) Hay transfer promotes establishment of Mediterranean steppe vegetation on soil disturbed by pipeline construction. Restor Ecol 19:214–222.
    .
    OpenUrl
  24. ↵
    1. Goldberg DE,
    2. Turkington R,
    3. Olsvig-Whittaker L,
    4. Dyer AR
    (2001) Density dependence in an annual plant community: Variation among life history stages. Ecol Monogr 71:423–446.
    .
    OpenUrlCrossRef
  25. ↵
    1. Baldeck CA, et al.
    (2013) Habitat filtering across tree life stages in tropical forest communities. Proc R Soc B Biol Sci R Soc 280:20130548.
    .
    OpenUrl
  26. ↵
    1. Harrison S
    (1999) Local and regional diversity in a patchy landscape: Native, alien, and endemic herbs on serpentine. Ecology 80:70–80.
    .
    OpenUrlCrossRef
  27. ↵
    1. Strauss SY,
    2. Webb CO,
    3. Salamin N
    (2006) Exotic taxa less related to native species are more invasive. Proc Natl Acad Sci USA 103:5841–5845.
    .
    OpenUrlAbstract/FREE Full Text
  28. ↵
    1. Anacker BL
    (2014) The nature of serpentine endemism. Am J Bot 101:219–224.
    .
    OpenUrlAbstract/FREE Full Text
  29. ↵
    1. Hanski I,
    2. Ovaskainen O
    (2002) Extinction debt at extinction threshold. Conserv Biol 16:666–673.
    .
    OpenUrlCrossRef
  30. ↵
    1. Szaro RC,
    2. Johnston DW
    1. Pearson S,
    2. Turner M,
    3. Gardner R,
    4. O’Neill R
    (1996) An organism-based perspective of habitat fragmentation. Biodiversity in Managed Landscapes: Theory and Practice, eds Szaro RC, Johnston DW (Oxford Univ Press, New York), pp 77–95.
    .
  31. ↵
    1. Chase JM,
    2. Kraft NJB,
    3. Smith KG,
    4. Velland M,
    5. Inouye B
    (2011) Using null models to disentangle variation in community dissimilarity from variation in α-diversity. Ecosphere 2:1–11.
    .
    OpenUrl
  32. ↵
    1. Karlson RH,
    2. Cornell HV
    (1998) Scale-dependent variation in local vs. regional effects on coral species richness. Ecol Monogr 68:259–274.
    .
    OpenUrlCrossRef
    1. McGill BJ
    (2010) Ecology. Matters of scale. Science 328:575–576.
    .
    OpenUrlAbstract/FREE Full Text
  33. ↵
    1. Borcard D,
    2. Legendre P,
    3. Avois-Jacquet C,
    4. Tuomisto H
    (2004) Dissecting the spatial structure of ecological data at multiple scales. Ecology 85:1826–1832.
    .
    OpenUrlCrossRef
  34. ↵
    1. Ramette A,
    2. Tiedje JM
    (2007) Multiscale responses of microbial life to spatial distance and environmental heterogeneity in a patchy ecosystem. Proc Natl Acad Sci USA 104:2761–2766.
    .
    OpenUrlAbstract/FREE Full Text
  35. ↵
    1. Levine JM
    (2001) Local interactions, dispersal, and native and exotic plant diversity along a California stream. Oikos 95:387–408.
    .
    OpenUrl
  36. ↵
    1. Foster BL,
    2. Tilman D
    (2003) Seed limitation and the regulation of community structure in oak savanna grassland. J Ecol 91:999–1007.
    .
    OpenUrlCrossRef
  37. ↵
    1. Tuomisto H,
    2. Ruokolainen K,
    3. Yli-Halla M
    (2003) Dispersal, environment, and floristic variation of western Amazonian forests. Science 299:241–244.
    .
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Legendre P, et al.
    (2009) Partitioning beta diversity in a subtropical broad-leaved forest of China. Ecology 90:663–674.
    .
    OpenUrlCrossRefPubMed
  39. ↵
    1. Peres-Neto PR,
    2. Legendre P,
    3. Dray S,
    4. Borcard D
    (2006) Variation partitioning of species data matrices: Estimation and comparison of fractions. Ecology 87:2614–2625.
    .
    OpenUrlCrossRefPubMed
  40. ↵
    1. Bennet JR,
    2. Gilbert B
    (2010) Partitioning variation in ecological communities: Do the numbers add up? J Appl Ecol 47:1071–1082.
    .
    OpenUrl
  41. ↵
    1. Thomson FJ,
    2. Moles AT,
    3. Auld TD,
    4. Kingsford RT
    (2011) Seed dispersal distance is more strongly correlated with plant height than with seed mass. J Ecol 99:1299–1307.
    .
    OpenUrlCrossRef
  42. ↵
    1. Harrison S,
    2. Rice K,
    3. Maron J
    (2001) Habitat patchiness promotes invasion by alien grasses on serpentine soil. Biol Conserv 100:45–53.
    .
    OpenUrlCrossRef
  43. ↵
    1. Hanski I,
    2. Ovaskainen O
    (2000) The metapopulation capacity of a fragmented landscape. Nature 404:755–758.
    .
    OpenUrlCrossRefPubMed
  44. ↵
    1. Harrison S,
    2. Rajakaruna R
    , eds (2011) Serpentine: The Evolution and Ecology of a Model System (University of California Press, Berkeley).
    .
  45. ↵
    1. Bartolome JW
    (1979) Germination and seedling establishment in California annual grassland. J Ecol 67:273–281.
    .
    OpenUrlCrossRef
  46. ↵
    1. Bolker BM, et al.
    (2009) Generalized linear mixed models: A practical guide for ecology and evolution. Trends Ecol Evol 24:127–135.
    .
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Article Alerts
Email Article

Thank you for your interest in spreading the word on PNAS.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Experimental dispersal reveals characteristic scales of biodiversity in a natural landscape
(Your Name) has sent you a message from PNAS
(Your Name) thought you would like to see the PNAS web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Spatial scale of ecological processes
Rachel M. Germain, Sharon Y. Strauss, Benjamin Gilbert
Proceedings of the National Academy of Sciences Apr 2017, 114 (17) 4447-4452; DOI: 10.1073/pnas.1615338114

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Request Permissions
Share
Spatial scale of ecological processes
Rachel M. Germain, Sharon Y. Strauss, Benjamin Gilbert
Proceedings of the National Academy of Sciences Apr 2017, 114 (17) 4447-4452; DOI: 10.1073/pnas.1615338114
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Mendeley logo Mendeley

Article Classifications

  • Biological Sciences
  • Ecology
Proceedings of the National Academy of Sciences: 114 (17)
Table of Contents

Submit

Sign up for Article Alerts

Jump to section

  • Article
    • Abstract
    • Results
    • Discussion
    • Materials and Methods
    • Acknowledgments
    • Footnotes
    • References
  • Figures & SI
  • Info & Metrics
  • PDF

You May Also be Interested in

Smoke emanates from Japan’s Fukushima nuclear power plant a few days after tsunami damage
Core Concept: Muography offers a new way to see inside a multitude of objects
Muons penetrate much further than X-rays, they do essentially zero damage, and they are provided for free by the cosmos.
Image credit: Science Source/Digital Globe.
Water from a faucet fills a glass.
News Feature: How “forever chemicals” might impair the immune system
Researchers are exploring whether these ubiquitous fluorinated molecules might worsen infections or hamper vaccine effectiveness.
Image credit: Shutterstock/Dmitry Naumov.
Venus flytrap captures a fly.
Journal Club: Venus flytrap mechanism could shed light on how plants sense touch
One protein seems to play a key role in touch sensitivity for flytraps and other meat-eating plants.
Image credit: Shutterstock/Kuttelvaserova Stuchelova.
Illustration of groups of people chatting
Exploring the length of human conversations
Adam Mastroianni and Daniel Gilbert explore why conversations almost never end when people want them to.
Listen
Past PodcastsSubscribe
Panda bear hanging in a tree
How horse manure helps giant pandas tolerate cold
A study finds that giant pandas roll in horse manure to increase their cold tolerance.
Image credit: Fuwen Wei.

Similar Articles

Site Logo
Powered by HighWire
  • Submit Manuscript
  • Twitter
  • Facebook
  • RSS Feeds
  • Email Alerts

Articles

  • Current Issue
  • Special Feature Articles – Most Recent
  • List of Issues

PNAS Portals

  • Anthropology
  • Chemistry
  • Classics
  • Front Matter
  • Physics
  • Sustainability Science
  • Teaching Resources

Information

  • Authors
  • Editorial Board
  • Reviewers
  • Subscribers
  • Librarians
  • Press
  • Cozzarelli Prize
  • Site Map
  • PNAS Updates
  • FAQs
  • Accessibility Statement
  • Rights & Permissions
  • About
  • Contact

Feedback    Privacy/Legal

Copyright © 2021 National Academy of Sciences. Online ISSN 1091-6490