Latitudinal patterns in stabilizing density dependence of forest communities.

Numerous studies have shown reduced performance in plants that are surrounded by neighbours of the same species1,2, a phenomenon known as conspecific negative density dependence (CNDD)3. A long-held ecological hypothesis posits that CNDD is more pronounced in tropical than in temperate forests4,5, which increases community stabilization, species coexistence and the diversity of local tree species6,7. Previous analyses supporting such a latitudinal gradient in CNDD8,9 have suffered from methodological limitations related to the use of static data10-12. Here we present a comprehensive assessment of latitudinal CNDD patterns using dynamic mortality data to estimate species-site-specific CNDD across 23 sites. Averaged across species, we found that stabilizing CNDD was present at all except one site, but that average stabilizing CNDD was not stronger toward the tropics. However, in tropical tree communities, rare and intermediate abundant species experienced stronger stabilizing CNDD than did common species. This pattern was absent in temperate forests, which suggests that CNDD influences species abundances more strongly in tropical forests than it does in temperate ones13. We also found that interspecific variation in CNDD, which might attenuate its stabilizing effect on species diversity14,15, was high but not significantly different across latitudes. Although the consequences of these patterns for latitudinal diversity gradients are difficult to evaluate, we speculate that a more effective regulation of population abundances could translate into greater stabilization of tropical tree communities and thus contribute to the high local diversity of tropical forests.

Tropical tree ectomycorrhiza are distributed independently of soil nutrients.

Mycorrhizae, a form of plant-fungal symbioses, mediate vegetation impacts on ecosystem functioning. Climatic effects on decomposition and soil quality are suggested to drive mycorrhizal distributions, with arbuscular mycorrhizal plants prevailing in low-latitude/high-soil-quality areas and ectomycorrhizal (EcM) plants in high-latitude/low-soil-quality areas. However, these generalizations, based on coarse-resolution data, obscure finer-scale variations and result in high uncertainties in the predicted distributions of mycorrhizal types and their drivers. Using data from 31 lowland tropical forests, both at a coarse scale (mean-plot-level data) and fine scale (20 × 20 metres from a subset of 16 sites), we demonstrate that the distribution and abundance of EcM-associated trees are independent of soil quality. Resource exchange differences among mycorrhizal partners, stemming from diverse evolutionary origins of mycorrhizal fungi, may decouple soil fertility from the advantage provided by mycorrhizal associations. Additionally, distinct historical biogeographies and diversification patterns have led to differences in forest composition and nutrient-acquisition strategies across three major tropical regions. Notably, Africa and Asia's lowland tropical forests have abundant EcM trees, whereas they are relatively scarce in lowland neotropical forests. A greater understanding of the functional biology of mycorrhizal symbiosis is required, especially in the lowland tropics, to overcome biases from assuming similarity to temperate and boreal regions.

Mycorrhizal feedbacks influence global forest structure and diversity.

One mechanism proposed to explain high species diversity in tropical systems is strong negative conspecific density dependence (CDD), which reduces recruitment of juveniles in proximity to conspecific adult plants. Although evidence shows that plant-specific soil pathogens can drive negative CDD, trees also form key mutualisms with mycorrhizal fungi, which may counteract these effects. Across 43 large-scale forest plots worldwide, we tested whether ectomycorrhizal tree species exhibit weaker negative CDD than arbuscular mycorrhizal tree species. We further tested for conmycorrhizal density dependence (CMDD) to test for benefit from shared mutualists. We found that the strength of CDD varies systematically with mycorrhizal type, with ectomycorrhizal tree species exhibiting higher sapling densities with increasing adult densities than arbuscular mycorrhizal tree species. Moreover, we found evidence of positive CMDD for tree species of both mycorrhizal types. Collectively, these findings indicate that mycorrhizal interactions likely play a foundational role in global forest diversity patterns and structure.

Long-term seedling and small sapling census data from the Barro Colorado Island 50 ha Forest Dynamics Plot, Panama.

Tropical forests are well known for their high woody plant diversity. Processes occurring at early life stages are thought to play a critical role in maintaining this high diversity and shaping the composition of tropical tree communities. To evaluate hypothesized mechanisms promoting tropical tree species coexistence and influencing composition, we initiated a census of woody seedlings and small saplings in the permanent 50 ha Forest Dynamics Plot (FDP) on Barro Colorado Island (BCI), Panama. Situated in old-growth, lowland tropical moist forest, the BCI FDP was originally established in 1980 to monitor trees and shrubs ≥1 cm diameter at 1.3 m above ground (dbh) at ca. 5-year intervals. However, critical data on the dynamics occurring at earlier life stages were initially lacking. Therefore, in 2001 we established a 1-m2 seedling plot in the center of every 5 × 5 m section of the BCI FDP. All freestanding woody individuals ≥20 cm tall and <1 cm dbh (hereafter referred to as seedlings) were tagged, mapped, measured, and identified to species in 19,313 1-m2 seedling plots. Because seedling dynamics are rapid, we censused these seedling plots every 1-2 years. Here, we present data from the 14 censuses of these seedling plots conducted between the initial census in 2001 to the most recent census, in 2018. This data set includes nearly 1 M observations of ~185,000 individuals of >400 tree, shrub, and liana species. These data will permit spatially-explicit analyses of seedling distributions, recruitment, growth, and survival for hundreds of woody plant species. In addition, the data presented here can be linked to openly-available, long-term data on the dynamics of trees and shrubs ≥1 cm dbh in the BCI FDP, as well as existing data sets from the site on climate, canopy structure, phylogenetic relatedness, functional traits, soil nutrients, and topography. This data set can be freely used for non-commercial purposes; we request that users of these data cite this data paper in all publications resulting from the use of this data set.

RecruitNet: A global database of plant recruitment networks.

Plant recruitment interactions (i.e., what recruits under what) shape the composition, diversity, and structure of plant communities. Despite the huge body of knowledge on the mechanisms underlying recruitment interactions among species, we still know little about the structure of the recruitment networks emerging in ecological communities. Modeling and analyzing the community-level structure of plant recruitment interactions as a complex network can provide relevant information on ecological and evolutionary processes acting both at the species and ecosystem levels. We report a data set containing 143 plant recruitment networks in 23 countries across five continents, including temperate and tropical ecosystems. Each network identifies the species under which another species recruits. All networks report the number of recruits (i.e., individuals) per species. The data set includes >850,000 recruiting individuals involved in 118,411 paired interactions among 3318 vascular plant species across the globe. The cover of canopy species and open ground is also provided. Three sampling protocols were used: (1) The Recruitment Network (RN) protocol (106 networks) focuses on interactions among established plants ("canopy species") and plants in their early stages of recruitment ("recruit species"). A series of plots was delimited within a locality, and all the individuals recruiting and their canopy species were identified; (2) The paired Canopy-Open (pCO) protocol (26 networks) consists in locating a potential canopy plant and identifying recruiting individuals under the canopy and in a nearby open space of the same area; (3) The Georeferenced plot (GP) protocol (11 networks) consists in using information from georeferenced individual plants in large plots to infer canopy-recruit interactions. Some networks incorporate data for both herbs and woody species, whereas others focus exclusively on woody species. The location of each study site, geographical coordinates, country, locality, responsible author, sampling dates, sampling method, and life habits of both canopy and recruit species are provided. This database will allow researchers to test ecological, biogeographical, and evolutionary hypotheses related to plant recruitment interactions. There are no copyright restrictions on the data set; please cite this data paper when using these data in publications.

Distribution of biomass dynamics in relation to tree size in forests across the world.

Tree size shapes forest carbon dynamics and determines how trees interact with their environment, including a changing climate. Here, we conduct the first global analysis of among-site differences in how aboveground biomass stocks and fluxes are distributed with tree size. We analyzed repeat tree censuses from 25 large-scale (4-52 ha) forest plots spanning a broad climatic range over five continents to characterize how aboveground biomass, woody productivity, and woody mortality vary with tree diameter. We examined how the median, dispersion, and skewness of these size-related distributions vary with mean annual temperature and precipitation. In warmer forests, aboveground biomass, woody productivity, and woody mortality were more broadly distributed with respect to tree size. In warmer and wetter forests, aboveground biomass and woody productivity were more right skewed, with a long tail towards large trees. Small trees (1-10 cm diameter) contributed more to productivity and mortality than to biomass, highlighting the importance of including these trees in analyses of forest dynamics. Our findings provide an improved characterization of climate-driven forest differences in the size structure of aboveground biomass and dynamics of that biomass, as well as refined benchmarks for capturing climate influences in vegetation demographic models.

Demographic composition, not demographic diversity, predicts biomass and turnover across temperate and tropical forests.

The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.

The relative abundance of languages: Neutral and non-neutral dynamics.

Credible estimates suggest that a large number of the nearly 7000 languages in the world could go extinct this century, a prospect with profound cultural, socioeconomic, and political ramifications. Despite its importance, we still have little predictive theory for language dynamics and richness. Critical to the language extinction problem, however, is to understand the dynamics of the number of speakers of languages, the dynamics of language abundance distributions (LADs). Many regional LADs are very similar to the bell-shaped distributions of relative species abundance predicted by neutral theory in ecology. Using the tenets of neutral theory, here we show that LADs can be understood as an equilibrium or disequilibrium between stochastic rates of origination and extinction of languages. However, neutral theory does not fit some regional LADs, which can be explained if the number of speakers has grown systematically faster in some languages than others, due to cultural factors and other non-neutral processes. Only the LADs of Australia and the United States, deviate from a bell-shaped pattern. These deviations are due to the documented higher, non-equilibrium extinction rates of low-abundance languages in these countries.

Local canopy disturbance as an explanation for long-term increases in liana abundance.

Canopy disturbance explains liana abundance and distribution within tropical forests and thus may also explain the widespread pattern of increasing liana abundance; however, this hypothesis remains untested. We used a 10-year study (2007-2017) of 117,100 rooted lianas in an old-growth Panamanian forest to test whether local canopy disturbance explains increasing liana abundance. We found that liana density increased 29.2% and basal area 12.5%. The vast majority of these increases were associated with clonal stem proliferation following canopy disturbance, particularly in liana-dense, low-canopy gaps, which had far greater liana increases than did undisturbed forest. Lianas may be ecological niche constructors, arresting tree regeneration in gaps and thus creating a high-light environment that favours sustained liana proliferation. Our findings demonstrate that liana abundance is increasing rapidly and their ability to proliferate via copious clonal stem production in canopy gaps explains much of their increase in this and possibly other tropical forests.

Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide.

Arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) associations are critical for host-tree performance. However, how mycorrhizal associations correlate with the latitudinal tree beta-diversity remains untested. Using a global dataset of 45 forest plots representing 2,804,270 trees across 3840 species, we test how AM and EcM trees contribute to total beta-diversity and its components (turnover and nestedness) of all trees. We find AM rather than EcM trees predominantly contribute to decreasing total beta-diversity and turnover and increasing nestedness with increasing latitude, probably because wide distributions of EcM trees do not generate strong compositional differences among localities. Environmental variables, especially temperature and precipitation, are strongly correlated with beta-diversity patterns for both AM trees and all trees rather than EcM trees. Results support our hypotheses that latitudinal beta-diversity patterns and environmental effects on these patterns are highly dependent on mycorrhizal types. Our findings highlight the importance of AM-dominated forests for conserving global forest biodiversity.

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.

Genome assemblies for two Neotropical trees: Jacaranda copaia and Handroanthus guayacan.

The lack of genomic resources for tropical canopy trees is impeding several research avenues in tropical forest biology. We present genome assemblies for two Neotropical hardwood species, Jacaranda copaia and Handroanthus (formerly Tabebuia) guayacan, that are model systems for research on tropical tree demography and flowering phenology. For each species, we combined Illumina short-read data with in vitro proximity-ligation (Chicago) libraries to generate an assembly. For Jacaranda copaia, we obtained 104X physical coverage and produced an assembly with N50/N90 scaffold lengths of 1.020/0.277 Mbp. For H. guayacan, we obtained 129X coverage and produced an assembly with N50/N90 scaffold lengths of 0.795/0.165 Mbp. J. copaia and H. guayacan assemblies contained 95.8% and 87.9% of benchmarking orthologs, although they constituted only 77.1% and 66.7% of the estimated genome sizes of 799 and 512 Mbp, respectively. These differences were potentially due to high repetitive sequence content (>59.31% and 45.59%) and high heterozygosity (0.5% and 0.8%) in each species. Finally, we compared each new assembly to a previously sequenced genome for Handroanthus impetiginosus using whole-genome alignment. This analysis indicated extensive gene duplication in H. impetiginosus since its divergence from H. guayacan.

The interspecific growth-mortality trade-off is not a general framework for tropical forest community structure.

Resource allocation within trees is a zero-sum game. Unavoidable trade-offs dictate that allocation to growth-promoting functions curtails other functions, generating a gradient of investment in growth versus survival along which tree species align, known as the interspecific growth-mortality trade-off. This paradigm is widely accepted but not well established. Using demographic data for 1,111 tree species across ten tropical forests, we tested the generality of the growth-mortality trade-off and evaluated its underlying drivers using two species-specific parameters describing resource allocation strategies: tolerance of resource limitation and responsiveness of allocation to resource access. Globally, a canonical growth-mortality trade-off emerged, but the trade-off was strongly observed only in less disturbance-prone forests, which contained diverse resource allocation strategies. Only half of disturbance-prone forests, which lacked tolerant species, exhibited the trade-off. Supported by a theoretical model, our findings raise questions about whether the growth-mortality trade-off is a universally applicable organizing framework for understanding tropical forest community structure.

Demographic trade-offs predict tropical forest dynamics.

Understanding tropical forest dynamics and planning for their sustainable management require efficient, yet accurate, predictions of the joint dynamics of hundreds of tree species. With increasing information on tropical tree life histories, our predictive understanding is no longer limited by species data but by the ability of existing models to make use of it. Using a demographic forest model, we show that the basal area and compositional changes during forest succession in a neotropical forest can be accurately predicted by representing tropical tree diversity (hundreds of species) with only five functional groups spanning two essential trade-offs-the growth-survival and stature-recruitment trade-offs. This data-driven modeling framework substantially improves our ability to predict consequences of anthropogenic impacts on tropical forests.

Lightning is a major cause of large tree mortality in a lowland neotropical forest.

The mortality rates of large trees are critical to determining carbon stocks in tropical forests, but the mechanisms of tropical tree mortality remain poorly understood. Lightning strikes thousands of tropical trees every day, but is commonly assumed to be a minor agent of tree mortality in most tropical forests. We use the first systematic quantification of lightning-caused mortality to show that lightning is a major cause of death for the largest trees in an old-growth lowland forest in Panama. A novel lightning strike location system together with field surveys of strike sites revealed that, on average, each strike directly kills 3.5 trees (> 10 cm diameter) and damages 11.4 more. Given lightning frequency data from the Earth Networks Total Lightning Network and historical total tree mortality rates for this site, we conclude that lightning accounts for 40.5% of the mortality of large trees (> 60 cm diameter) in the short term and probably contributes to an additional 9.0% of large tree deaths over the long term. Any changes in cloud-to-ground lightning frequency due to climatic change will alter tree mortality rates; projected 25-50% increases in lightning frequency would increase large tree mortality rates in this forest by 9-18%. The results of this study indicate that lightning plays a critical and previously underestimated role in tropical forest dynamics and carbon cycling.

Testing for changes in biomass dynamics in large-scale forest datasets.

Tropical forest responses to climate and atmospheric change are critical to the future of the global carbon budget. Recent studies have reported increases in estimated above-ground biomass (EAGB) stocks, productivity, and mortality in old-growth tropical forests. These increases could reflect a shift in forest functioning due to global change and/or long-lasting recovery from past disturbance. We introduce a novel approach to disentangle the relative contributions of these mechanisms by decomposing changes in whole-plot biomass fluxes into contributions from changes in the distribution of gap-successional stages and changes in fluxes for a given stage. Using 30 years of forest dynamic data at Barro Colorado Island, Panama, we investigated temporal variation in EAGB fluxes as a function of initial EAGB (EAGBi ) in 10 × 10 m quadrats. Productivity and mortality fluxes both increased strongly with initial quadrat EAGB. The distribution of EAGB (and thus EAGBi ) across quadrats hardly varied over 30 years (and seven censuses). EAGB fluxes as a function of EAGBi varied largely and significantly among census intervals, with notably higher productivity in 1985-1990 associated with recovery from the 1982-1983 El Niño event. Variation in whole-plot fluxes among census intervals was explained overwhelmingly by variation in fluxes as a function of EAGBi , with essentially no contribution from changes in EAGBi distributions. The high observed temporal variation in productivity and mortality suggests that this forest is very sensitive to climate variability. There was no consistent long-term trend in productivity, mortality, or biomass in this forest over 30 years, although the temporal variability in productivity and mortality was so strong that it could well mask a substantial trend. Accurate prediction of future tropical forest carbon budgets will require accounting for disturbance-recovery dynamics and understanding temporal variability in productivity and mortality.

Seed-to-seedling transitions exhibit distance-dependent mortality but no strong spacing effects in a Neotropical forest.

Patterns of seed dispersal and seed mortality influence the spatial structure of plant communities and the local coexistence of competing species. Most seeds are dispersed in proximity to the parent tree, where mortality is also expected to be the highest, because of competition with siblings or the attraction of natural enemies. Whereas distance-dependent mortality in the seed-to-seedling transition was often observed in tropical forests, few studies have attempted to estimate the shape of the survival-distance curves, which determines whether the peak of seedling establishment occurs away from the parent tree (Janzen-Connell pattern) or if the peak attenuates but remains at the parent location (Hubbell pattern). In this study, we inferred the probability density of seed dispersal and two stages of seedling establishment (new recruits, and seedlings 20 cm or taller) with distance for 24 tree species present in the 50-ha Forest Dynamics Plot of Barro Colorado Island, Panama. Using data from seed traps, seedling survey quadrats, and tree-census records spanning the 1988-2014 period, we fit hierarchical Bayesian models including parameters for tree fecundity, the shape of the dispersal kernel, and overdispersion of seed or seedling counts. We combined predictions from multiple dispersal kernels to obtain more robust inferences. We find that Hubbell patterns are the most common and Janzen-Connell patterns are very rare among those species; that distance-dependent mortality may be stronger in the seed stage, in the early recruit stage, or comparable in both; and that species with larger seeds experience less overall mortality and less distance-dependent mortality. Finally, we describe how this modeling approach could be extended at a community scale to include less abundant species.

Pervasive decreases in living vegetation carbon turnover time across forest climate zones.

Forests play a major role in the global carbon cycle. Previous studies on the capacity of forests to sequester atmospheric CO2 have mostly focused on carbon uptake, but the roles of carbon turnover time and its spatiotemporal changes remain poorly understood. Here, we used long-term inventory data (1955 to 2018) from 695 mature forest plots to quantify temporal trends in living vegetation carbon turnover time across tropical, temperate, and cold climate zones, and compared plot data to 8 Earth system models (ESMs). Long-term plots consistently showed decreases in living vegetation carbon turnover time, likely driven by increased tree mortality across all major climate zones. Changes in living vegetation carbon turnover time were negatively correlated with CO2 enrichment in both forest plot data and ESM simulations. However, plot-based correlations between living vegetation carbon turnover time and climate drivers such as precipitation and temperature diverged from those of ESM simulations. Our analyses suggest that forest carbon sinks are likely to be constrained by a decrease in living vegetation carbon turnover time, and accurate projections of forest carbon sink dynamics will require an improved representation of tree mortality processes and their sensitivity to climate in ESMs.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix).

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

Biodiversity recovery of Neotropical secondary forests.

Old-growth tropical forests harbor an immense diversity of tree species but are rapidly being cleared, while secondary forests that regrow on abandoned agricultural lands increase in extent. We assess how tree species richness and composition recover during secondary succession across gradients in environmental conditions and anthropogenic disturbance in an unprecedented multisite analysis for the Neotropics. Secondary forests recover remarkably fast in species richness but slowly in species composition. Secondary forests take a median time of five decades to recover the species richness of old-growth forest (80% recovery after 20 years) based on rarefaction analysis. Full recovery of species composition takes centuries (only 34% recovery after 20 years). A dual strategy that maintains both old-growth forests and species-rich secondary forests is therefore crucial for biodiversity conservation in human-modified tropical landscapes.