skip to main content


Title: Interactions between all pairs of neighboring trees in 16 forests worldwide reveal details of unique ecological processes in each forest, and provide windows into their evolutionary histories
When Darwin visited the Galapagos archipelago, he observed that, in spite of the islands’ physical similarity, members of species that had dispersed to them recently were beginning to diverge from each other. He postulated that these divergences must have resulted primarily from interactions with sets of other species that had also diverged across these otherwise similar islands. By extrapolation, if Darwin is correct, such complex interactions must be driving species divergences across all ecosystems. However, many current general ecological theories that predict observed distributions of species in ecosystems do not take the details of between-species interactions into account. Here we quantify, in sixteen forest diversity plots (FDPs) worldwide, highly significant negative density-dependent (NDD) components of both conspecific and heterospecific between-tree interactions that affect the trees’ distributions, growth, recruitment, and mortality. These interactions decline smoothly in significance with increasing physical distance between trees. They also tend to decline in significance with increasing phylogenetic distance between the trees, but each FDP exhibits its own unique pattern of exceptions to this overall decline. Unique patterns of between-species interactions in ecosystems, of the general type that Darwin postulated, are likely to have contributed to the exceptions. We test the power of our null-model method by using a deliberately modified data set, and show that the method easily identifies the modifications. We examine how some of the exceptions, at the Wind River (USA) FDP, reveal new details of a known allelopathic effect of one of the Wind River gymnosperm species. Finally, we explore how similar analyses can be used to investigate details of many types of interactions in these complex ecosystems, and can provide clues to the evolution of these interactions.  more » « less
Award ID(s):
1831952
NSF-PAR ID:
10304280
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; more » ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; « less
Editor(s):
Pascual, Mercedes
Date Published:
Journal Name:
PLOS Computational Biology
Volume:
17
Issue:
4
ISSN:
1553-7358
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Aim

    Spatially explicit protections of coastal habitats determined on the current distribution of species and ecosystems risk becoming obsolete in 100 years if the movement of species ranges outpaces management action. Hence, a critical step of conservation is predicting the efficacy of management actions in future. We aimed to determine how foundational, habitat‐building species will respond to climate change in Fiji.

    Location

    The Republic of Fiji.

    Methods

    We develop species distribution models (SDMs) using MaxEnt, General Additive Models and Boosted Regression Trees and publicly available data from the Global Biodiversity Information Facility to predict changes in distribution of suitable habitat for mangrove forests, coral habitat, seagrass meadows and critical fisheries invertebrates under several IPCC climate change scenarios in 2070 or 2100. We then overlay predicted distribution models onto existing Fijian protected area network to assess whether today's conservation measures will afford protection to tomorrow's distributions.

    Results

    We found that mangrove suitability is projected to decrease along the Coral Coast and increase northward towards the Yasawa Islands due to precipitation changes. The response of seagrass meadows was predicted to be inconsistent and dependent on the climate scenario. Meanwhile, suitability for coral reefs was not predicted to decline significantly overall. The mangrove crabScylla serrata, an important resource for fisherwomen in Fiji, is projected to increase in habitat suitability while economically important sea cucumber species will have highly variable responses to climate change.

    Main conclusions

    Species distribution models are a critical tool for conservation managers, as linking spatial distribution data with future climate change scenarios can aid in the creation and resiliency of protected area programmes. New protected area designations should consider the future distribution of species to maximize benefits to those taxa.

     
    more » « less
  2. All life on earth is linked by a shared evolutionary history. Even before Darwin developed the theory of evolution, Linnaeus categorized types of organisms based on their shared traits. We now know these traits derived from these species’ shared ancestry. This evolutionary history provides a natural framework to harness the enormous quantities of biological data being generated today. The Open Tree of Life project is a collaboration developing tools to curate and share evolutionary estimates (phylogenies) covering the entire tree of life (Hinchliff et al. 2015, McTavish et al. 2017). The tree is viewable at https://tree.opentreeoflife.org, and the data is all freely available online. The taxon identifiers used in the Open Tree unified taxonomy (Rees and Cranston 2017) are mapped to identifiers across biological informatics databases, including the Global Biodiversity Information Facility (GBIF), NCBI, and others. Linking these identifiers allows researchers to easily unify data from across these different resources (Fig. 1). Leveraging a unified evolutionary framework across the diversity of life provides new avenues for integrative wide scale research. Downstream tools, such as R packages developed by the R OpenSci foundation (rotl, rgbif) (Michonneau et al. 2016, Chamberlain 2017) and others tools (Revell 2012), make accessing and combining this information straightforward for students as well as researchers (e.g. https://mctavishlab.github.io/BIO144/labs/rotl-rgbif.html). Figure 1. Example linking phylogenetic relationships accessed from the Open Tree of Life with specimen location data from Global Biodiversity Information Facility. For example, a recent publication by Santorelli et al. 2018 linked evolutionary information from Open Tree with species locality data gathered from a local field study as well as GBIF species location records to test a river-barrier hypothesis in the Amazon. By combining these data, the authors were able test a widely held biogeographic hypothesis across 1952 species in 14 taxonomic groups, and found that a river that had been postulated to drive endemism, was in fact not a barrier to gene flow. However, data provenance and taxonomic name reconciliation remain key hurdles to applying data from these large digital biodiversity and evolution community resources to answering biological questions. In the Amazonian river analysis, while they leveraged use of GBIF records as a secondary check on their species records, they relied on their an intensive local field study for their major conclusions, and preferred taxon specific phylogenetic resources over Open Tree where they were available (Santorelli et al. 2018). When Li et al. 2018 assessed large scale phylogenetic approaches, including Open Tree, for measuring community diversity, they found that synthesis phylogenies were less resolved than purpose-built phylogenies, but also found that these synthetic phylogenies were sufficient for community level phylogenetic diversity analyses. Nonetheless, data quality concerns have limited adoption of analyses data from centralized resources (McTavish et al. 2017). Taxonomic name recognition and reconciliation across databases also remains a hurdle for large scale analyses, despite several ongoing efforts to improve taxonomic interoperability and unify taxonomies, such at Catalogue of Life + (Bánki et al. 2018). In order to support innovative science, large scale digital data resources need to facilitate data linkage between resources, and address researchers' data quality and provenance concerns. I will present the model that the Open Tree of Life is using to provide evolutionary data at the scale of the entire tree of life, while maintaining traceable provenance to the publications and taxonomies these evolutionary relationships are inferred from. I will discuss the hurdles to adoption of these large scale resources by researchers, as well as the opportunities for new research avenues provided by the connections between evolutionary inferences and biodiversity digital databases. 
    more » « less
  3. The phenology of critical biological events in aquatic ecosystems are rapidly shifting due to climate change. Growing variability in phenological cues can increase the likelihood of trophic mismatches, causing recruitment failures in commercially, culturally, and recreationally important fisheries. We tested for changes in spawning phenology of regionally important walleye (Sander vitreus) populations in 194 Midwest US lakes in Minnesota, Michigan, and Wisconsin spanning 1939-2019 to investigate factors influencing walleye phenological responses to climate change and associated climate variability, including ice-off timing, lake physical characteristics, and population stocking history. Data from Wisconsin and Michigan lakes (185 and 5 out of 194 total lakes, respectively) were collected by the Wisconsin Department of Natural Resources (WDNR) and the Great Lakes Indian Fish and Wildlife Commission (GLIFWC) through standardized spring walleye mark-recapture surveys and spring tribal harvest season records. Standardized spring mark-recapture population estimates are performed shortly after ice-off, where following a marking event, a subsequent recapture sampling event is conducted using nighttime electrofishing (typically AC – WDNR, pulsed-DC – GLIFWC) of the entire shoreline including islands for small lakes and index stations for large lakes (Hansen et al. 2015) that is timed to coincide with peak walleye spawning activity (G. Hatzenbeler, WDNR, personal communication; M. Luehring, GLIFWC, personal communication; Beard et al. 1997). Data for four additional Minnesota lakes were collected by the Minnesota Department of Natural Resources (MNDNR) beginning in 1939 during annual collections of walleye eggs and broodstock (Schneider et al. 2010), where date of peak egg take was used to index peak spawning activity. For lakes where spawning location did not match the lake for which the ice-off data was collected, the spawning location either flowed into (Pike River) or was within 50 km of a lake where ice-off data were available (Pine River) and these ice-off data were used. Following the affirmation of off-reservation Ojibwe tribal fishing rights in the Ceded Territories of Wisconsin and the Upper Peninsula of Michigan in 1987, tribal spearfishers have targeted walleye during spring spawning (Mrnak et al. 2018). Nightly harvests are recorded as part of a compulsory creel survey (US Department of the Interior 1991). Using these records, we calculated the date of peak spawning activity in a given lake-year as the day of maximum tribal harvest. Although we were unable to account for varying effort in these data, a preliminary analysis comparing spawning dates estimated using tribal harvest to those determined from standardized agency surveys in the same lake and year showed that they were highly correlated (Pearson’s correlation: r = 0.91, P < 0.001). For lakes that had walleye spawning data from both agency surveys and tribal harvest, we used the data source with the greatest number of observation years. Ice-off phenology data was collected from two sources – either observed from the Global Lake and River Ice Phenology database (Benson et al. 2000)t, or modeled from a USGS region-wide machine-learning model which used North American Land Data Assimilation System (NLDAS) meteorological inputs combined with lake characteristics (lake position, clarity, size, depth, hypsography, etc.) to predict daily water column temperatures from 1979 - 2022, from which ice-off dates could be derived (https://www.sciencebase.gov/catalog/item/6206d3c2d34ec05caca53071; see Corson-Dosch et al. 2023 for details). Modeled data for our study lakes (see (Read et al. 2021) for modeling details), which performed well in reflecting ice phenology when compared to observed data (i.e., highly significant correlation between observed and modeled ice-off dates when both were available; r = 0.71, p < 0.001). Lake surface area (ha), latitude, and maximum depth (m) were acquired from agency databases and lake reports. Lake class was based on a WDNR lakes classification system (Rypel et al. 2019) that categorized lakes based on temperature, water clarity, depth, and fish community. Walleye stocking history was defined using the walleye stocking classification system developed by the Wisconsin Technical Working Group (see also Sass et al. 2021), which categorized lakes based on relative contributions of naturally-produced and stocked fish to adult recruitment by relying heavily on historic records of age-0 and age-1 catch rates and stocking histories. Wisconsin lakes were divided into three groups: natural recruitment (NR), a combination of stocking and natural recruitment (C-ST), and stocked only (ST). Walleye natural recruitment was indexed as age-0 walleye CPE (number of age-0 walleye captured per km of shoreline electrofished) from WDNR and GLIFWC fall electrofishing surveys (see Hansen et al. 2015 for details). We excluded lake-years where stocking of age-0 fish occurred before age-0 surveys to only include measurements of naturally-reproduced fish. 
    more » « less
  4. Abstract Island biogeography has classically focused on abiotic drivers of species distributions. However, recent work has highlighted the importance of mutualistic biotic interactions in structuring island floras. The limited occurrence of specialist pollinators and mycorrhizal fungi have been found to restrict plant colonization on oceanic islands. Another important mutualistic association occurs between nearly 15,000 plant species and nitrogen-fixing (N-fixing) bacteria. Here, we look for evidence that N-fixing bacteria limit establishment of plants that associate with them. Globally, we find that plants associating with N-fixing bacteria are disproportionately underrepresented on islands, with a 22% decline. Further, the probability of N-fixing plants occurring on islands decreases with island isolation and, where present, the proportion of N-fixing plant species decreases with distance for large, but not small islands. These findings suggest that N-fixing bacteria serve as a filter to plant establishment on islands, altering global plant biogeography, with implications for ecosystem development and introduction risks. 
    more » « less
  5. null (Ed.)
    ABSTRACT Miniature insects must overcome significant viscous resistance in order to fly. They typically possess wings with long bristles on the fringes and use a clap-and-fling mechanism to augment lift. These unique solutions to the extreme conditions of flight at tiny sizes (<2 mm body length) suggest that natural selection has optimized wing design for better aerodynamic performance. However, species vary in wingspan, number of bristles (n) and bristle gap (G) to diameter (D) ratio (G/D). How this variation relates to body length (BL) and its effects on aerodynamics remain unknown. We measured forewing images of 38 species of thrips and 21 species of fairyflies. Our phylogenetic comparative analyses showed that n and wingspan scaled positively and similarly with BL across both groups, whereas G/D decreased with BL, with a sharper decline in thrips. We next measured aerodynamic forces and visualized flow on physical models of bristled wings performing clap-and-fling kinematics at a chord-based Reynolds number of 10 using a dynamically scaled robotic platform. We examined the effects of dimensional (G, D, wingspan) and non-dimensional (n, G/D) geometric variables on dimensionless lift and drag. We found that: (1) increasing G reduced drag more than decreasing D; (2) changing n had minimal impact on lift generation; and (3) varying G/D minimally affected aerodynamic forces. These aerodynamic results suggest little pressure to functionally optimize n and G/D. Combined with the scaling relationships between wing variables and BL, much wing variation in tiny flying insects might be best explained by underlying shared growth factors. 
    more » « less