High exposure of global tree diversity to human pressure.

Safeguarding Earth's tree diversity is a conservation priority due to the importance of trees for biodiversity and ecosystem functions and services such as carbon sequestration. Here, we improve the foundation for effective conservation of global tree diversity by analyzing a recently developed database of tree species covering 46,752 species. We quantify range protection and anthropogenic pressures for each species and develop conservation priorities across taxonomic, phylogenetic, and functional diversity dimensions. We also assess the effectiveness of several influential proposed conservation prioritization frameworks to protect the top 17% and top 50% of tree priority areas. We find that an average of 50.2% of a tree species' range occurs in 110-km grid cells without any protected areas (PAs), with 6,377 small-range tree species fully unprotected, and that 83% of tree species experience nonnegligible human pressure across their range on average. Protecting high-priority areas for the top 17% and 50% priority thresholds would increase the average protected proportion of each tree species' range to 65.5% and 82.6%, respectively, leaving many fewer species (2,151 and 2,010) completely unprotected. The priority areas identified for trees match well to the Global 200 Ecoregions framework, revealing that priority areas for trees would in large part also optimize protection for terrestrial biodiversity overall. Based on range estimates for >46,000 tree species, our findings show that a large proportion of tree species receive limited protection by current PAs and are under substantial human pressure. Improved protection of biodiversity overall would also strongly benefit global tree diversity.

Hill-Chao numbers allow decomposing gamma multifunctionality into alpha and beta components.

Biodiversity-ecosystem functioning (BEF) research has provided strong evidence and mechanistic underpinnings to support positive effects of biodiversity on ecosystem functioning, from single to multiple functions. This research has provided knowledge gained mainly at the local alpha scale (i.e. within ecosystems), but the increasing homogenization of landscapes in the Anthropocene has raised the potential that declining biodiversity at the beta (across ecosystems) and gamma scales is likely to also impact ecosystem functioning. Drawing on biodiversity theory, we propose a new statistical framework based on Hill-Chao numbers. The framework allows decomposition of multifunctionality at gamma scales into alpha and beta components, a critical but hitherto missing tool in BEF research; it also allows weighting of individual ecosystem functions. Through the proposed decomposition, new BEF results for beta and gamma scales are discovered. Our novel approach is applicable across ecosystems and connects local- and landscape-scale BEF assessments from experiments to natural settings.

Dimensions of invasiveness: Links between local abundance, geographic range size, and habitat breadth in Europe's alien and native floras.

Understanding drivers of success for alien species can inform on potential future invasions. Recent conceptual advances highlight that species may achieve invasiveness via performance along at least three distinct dimensions: 1) local abundance, 2) geographic range size, and 3) habitat breadth in naturalized distributions. Associations among these dimensions and the factors that determine success in each have yet to be assessed at large geographic scales. Here, we combine data from over one million vegetation plots covering the extent of Europe and its habitat diversity with databases on species' distributions, traits, and historical origins to provide a comprehensive assessment of invasiveness dimensions for the European alien seed plant flora. Invasiveness dimensions are linked in alien distributions, leading to a continuum from overall poor invaders to super invaders-abundant, widespread aliens that invade diverse habitats. This pattern echoes relationships among analogous dimensions measured for native European species. Success along invasiveness dimensions was associated with details of alien species' introduction histories: earlier introduction dates were positively associated with all three dimensions, and consistent with theory-based expectations, species originating from other continents, particularly acquisitive growth strategists, were among the most successful invaders in Europe. Despite general correlations among invasiveness dimensions, we identified habitats and traits associated with atypical patterns of success in only one or two dimensions-for example, the role of disturbed habitats in facilitating widespread specialists. We conclude that considering invasiveness within a multidimensional framework can provide insights into invasion processes while also informing general understanding of the dynamics of species distributions.

Long-term Consequences on Soil Fungal Community Structure: Monoculture Planting and Natural Regeneration.

Understanding the regeneration and succession of belowground communities, particularly in forests, is vital for maintaining ecosystem health. Despite its importance, there is limited knowledge regarding how fungal communities change over time during ecosystem development, especially under different forest restoration strategies. In this study, we focused on two restoration methods used in northern Japan: monoculture planting and natural regeneration. We examined the responses of the fungal community to monoculture plantations (active tree planting) and naturally regenerated (passive regeneration) forests over a 50-year chronosequence, using natural forests as a reference. Based on DNA metabarcoding, we assessed the richness of fungal Operational Taxonomic Units (OTUs) and their dissimilarity. Our findings revealed that soil fungal richness remained stable after natural regeneration but declined in monoculture plantations, from 354 to 247 OTUs. While the compositional dissimilarity of fungal assemblages between monoculture plantations and natural forests remained consistent regardless of the time since tree planting, it significantly decreased after natural regeneration, suggesting recovery to a state close to the reference level. Notably, the composition of key functional fungal groups-saprotrophic and ectomycorrhizal- has increasingly mirrored that of natural forests over time following passive natural regeneration. In summary, our study suggests that monoculture plantations may not be effective for long-term ecosystem function and service recovery because of their limited support for soil fungal diversity. These results underscore the importance of natural regeneration in forest restoration and management strategies.

Geography-Dependent Horizontal Gene Transfer from Vertebrate Predators to Their Prey.

Horizontal transfer (HT) of genes between multicellular animals, once thought to be extremely rare, is being more commonly detected, but its global geographic trend and transfer mechanism have not been investigated. We discovered a unique HT pattern of Bovine-B (BovB) LINE retrotransposons in vertebrates, with a bizarre transfer direction from predators (snakes) to their prey (frogs). At least 54 instances of BovB HT were detected, which we estimate to have occurred across time between 85 and 1.3 Ma. Using comprehensive transcontinental sampling, our study demonstrates that BovB HT is highly prevalent in one geographical region, Madagascar, suggesting important regional differences in the occurrence of HTs. We discovered parasite vectors that may plausibly transmit BovB and found that the proportion of BovB-positive parasites is also high in Madagascar where BovB thus might be physically transported by parasites to diverse vertebrates, potentially including humans. Remarkably, in two frog lineages, BovB HT occurred after migration from a non-HT area (Africa) to the HT hotspot (Madagascar). These results provide a novel perspective on how the prevalence of parasites influences the occurrence of HT in a region, similar to pathogens and their vectors in some endemic diseases.

Liver Resection for Hepatocellular Carcinoma With Tumor Thrombus in the Inferior Vena Cava or Right Atrium: A Large-scale Multicenter Survey Conducted in Japan.

OBJECTIVE: To clarify the short and long-term postoperative outcomes and surgical indications for patients accompanied by hepatocellular carcinoma with tumor thrombus (TT) in the inferior vena cava (IVC) or right atrium (RA). BACKGROUND: These patients are known to have an extremely poor prognosis; however, the postoperative outcomes have not been fully verified because of the rarity of this disease. METHODS: We contacted 211 specialized centers in Japan and collected data on liver resection for hepatocellular carcinoma with TT in the IVC or RA from centers with experience performing surgery for such patients. The patient characteristics, operative procedures, and surgical outcomes were then analyzed. RESULTS: A total of 119 patients from 23 institutions were enrolled; 49 patients had TT in the IVC below the diaphragm (type I), 42 had TT in the IVC above the diaphragm (type II), and 28 had TT entering the RA (type III). The severity and frequency of postoperative complications did not differ among the 3 groups. There was one surgery-related death in the type III group. The median survival times were 2.47 years in the type I group, 1.77 years in the type II group, and 1.02 years in the type III group. Multivariate analysis identified an indocyanine green retention rate at 15 minutes >15% and ≥3 tumors as prognostic factors affecting survival, whereas the use of cardiopulmonary bypass and ≥3 tumors were risk factors for recurrence. CONCLUSIONS: As the postoperative prognosis of patients with type I or type II disease and of patients with no risk factors is relatively good, surgery should be considered for these patient populations.

Propensity Score Matching Analysis of the Safety of Completion Total Pancreatectomy for Remnant Pancreatic Tumors Versus that of Initial Total Pancreatectomy for Primary Pancreatic Tumors.

BACKGROUND: The safety and feasibility of completion total pancreatectomy (TP) for remnant pancreatic neoplasms remain controversial and are rarely compared with that of initial TP. Thus, we aimed to compare the safety of these two procedures inducing a pancreatic state. METHODS: Patients who underwent TP for pancreatic neoplasms between 2006 and 2018 at our institution were included in this study. Tumor pathologies were classified into three subgroups according to survival curves. We used 1:1 propensity score matching (PSM) to analyze age, sex, Charlson Comorbidity Index, and tumor stage. Finally, we analyzed the primary outcome Clavien-Dindo classification (CDC) grade, risks of other safety-related outcomes, and the survival rate of patients with invasive cancer. RESULTS: Of 54 patients, 16 underwent completion TP (29.6%) and 38 (70.4%) underwent initial TP. Before PSM analysis, age and Charlson Comorbidity Index were significantly higher, and T category and stage were significantly lower for the completion TP group. Upon PSM analysis, these two groups were equivalent in CDC grade [initial TP vs. completion TP: 71.4% (10/14) vs. 78.6% (11/14); p = 0.678] and other safety-related outcomes. Additionally, while the overall survival and recurrence-free survival of patients with invasive cancer were not significantly different between these two groups, the T category and stage tended to be remarkably severe in the initial TP group. CONCLUSIONS: PSM analysis for prognostic factors showed that completion TP and initial TP have similar safety-related outcomes that can be used as a decision-making reference in the surgery of pancreatic tumors.

Enhancing the structural diversity between forest patches-A concept and real-world experiment to study biodiversity, multifunctionality and forest resilience across spatial scales.

Intensification of land use by humans has led to a homogenization of landscapes and decreasing resilience of ecosystems globally due to a loss of biodiversity, including the majority of forests. Biodiversity-ecosystem functioning (BEF) research has provided compelling evidence for a positive effect of biodiversity on ecosystem functions and services at the local (α-diversity) scale, but we largely lack empirical evidence on how the loss of between-patch ß-diversity affects biodiversity and multifunctionality at the landscape scale (γ-diversity). Here, we present a novel concept and experimental framework for elucidating BEF patterns at α-, ß-, and γ-scales in real landscapes at a forest management-relevant scale. We examine this framework using 22 temperate broadleaf production forests, dominated by Fagus sylvatica. In 11 of these forests, we manipulated the structure between forest patches by increasing variation in canopy cover and deadwood. We hypothesized that an increase in landscape heterogeneity would enhance the ß-diversity of different trophic levels, as well as the ß-functionality of various ecosystem functions. We will develop a new statistical framework for BEF studies extending across scales and incorporating biodiversity measures from taxonomic to functional to phylogenetic diversity using Hill numbers. We will further expand the Hill number concept to multifunctionality allowing the decomposition of γ-multifunctionality into α- and ß-components. Combining this analytic framework with our experimental data will allow us to test how an increase in between patch heterogeneity affects biodiversity and multifunctionality across spatial scales and trophic levels to help inform and improve forest resilience under climate change. Such an integrative concept for biodiversity and functionality, including spatial scales and multiple aspects of diversity and multifunctionality as well as physical and environmental structure in forests, will go far beyond the current widely applied approach in forestry to increase resilience of future forests through the manipulation of tree species composition.

Drivers of community assembly change during succession in wood-decomposing beetle communities.

The patterns of successional change of decomposer communities is unique in that resource availability predictably decreases as decomposition proceeds. Saproxylic (i.e. deadwood-dependent) beetles are a highly diverse and functionally important decomposer group, and their community composition is affected by both deadwood characteristics and other environmental factors. Understanding how communities change with faunal succession through the decomposition process is important as this process influences terrestrial carbon dynamics. Here, we evaluate how beta-diversity of saproxylic beetle communities change with succession, as well as the effects of different major drivers of beta-diversity, such as deadwood tree species, spatial distance between locations, climate and forest structure. We studied spatial beta-diversity (i.e. dissimilarity of species composition between deadwood logs in the same year) of saproxylic beetle communities over 8 years of wood decomposition. Our study included 379 experimental deadwood logs comprising 13 different tree species in 30 forest stands in Germany. We hypothesized that the effects of tree species dissimilarity, measured by phylogenetic distance, and climate on beta-diversity decrease over time, while the effects of spatial distance between logs and forest structure increase. Observed beta-diversity of saproxylic beetle communities increased over time, whereas standardized effects sizes (SES; based on null models) of beta-diversity decreased indicating higher beta-diversity than expected during early years. Beta-diversity increased with increasing phylogenetic distance between tree species and spatial distance among regions, and to a lesser extent with spatial distance within regions and differences in climate and forest structure. Whereas effects of space, climate and forest structure were constant over time, the effect of phylogenetic distance decreased. Our results show that the strength of the different drivers of saproxylic beetle community beta-diversity changes along deadwood succession. Beta-diversity of early decay communities was strongly associated with differences among tree species. Although this effect decreased over time, beta-diversity remained high throughout succession. Possible explanations for this pattern include differences in decomposition rates and fungal communities between logs or the priority effect of early successional communities. Our results suggest that saproxylic beetle diversity can be enhanced by promoting forests with diverse tree communities and structures.

Intrinsic Factors Associated with Dietary Toxin Quantity and Concentration in the Nuchal Glands of a Natricine Snake Rhabdophis Tigrinus.

The snake Rhabdophis tigrinus sequesters cardiotonic steroids, bufadienolides (BDs), from ingested toads and stores them in the nuchal glands as defensive toxins. It has previously been shown that there are individual differences in the total quantity of BDs stored in the nuchal glands of adult R. tigrinus and that BD quantities and profiles of R. tigrinus exhibit geographic variation. However, no previous study has examined the total quantity of BDs as a percentage of body mass (relative BD quantity) and the concentration of BDs in the nuchal gland fluid (BD gland concentration). In addition, intrinsic factors that are associated with relative BD quantity and BD concentration have not been examined within a single population. We collected 158 adult snakes from an area of central Japan from May to October and analyzed their BD quantities by UV analysis. We assessed individual differences in BD quantity, relative BD quantity and BD gland concentration. We found that 1) in approximately 60% of the 158 individuals, the BD gland concentration was greater than 50%; 2) body length and body condition are positively correlated with relative BD quantity and BD gland concentration; 3) even in a single population, individual differences of BD quantity are large, and are greater in females than in males; and 4) relative BD quantity and BD gland concentration of females during the gestation season are lower than those during the non-gestation season.

DO Toxic Invasive Prey Become a Toxin Source for Native Consumers?

Toxic organisms can become food that potentially harms consumers. When these organisms become invasive species, the harm often turns to a serious threat that disrupts native ecosystems. On the other hand, there are consumers that can exploit toxic organisms for food and sequester intact toxins from them for the consumers' own chemical defense. Therefore, it can be expected that toxic invasive prey can become a toxin source for native consumers. Here, we focused on the relationship between toads, which are one of the major toxic invasive organisms and possess bufadienolides (BDs), and Rhabdophis snakes, which sequester BDs from toads. On Sado Island, Japan, R. tigrinus is native, but no toads had inhabited this island until Bufo japonicus formosus was introduced as a domestic invasive species in 1963 and 1964. At present, invasive toads are distributed only in the southwestern part of the island. We collected a total of 25 and 24 R. tigrinus from areas allopatric and sympatric with toads, respectively. Then, we investigated the possession of BDs and the BD profile of these snakes. We found that only R. tigrinus sympatric with toads possessed BDs, whereas all snakes allopatric with toads lacked BDs. Based on the characteristics of the BD profile, the toxin source was identified as B. j. formosus. Our findings show that a new case of impact caused by toxic invasive species, i.e., "toxin supply to native consumers from invasive prey", could occur.

Dramatic dietary shift maintains sequestered toxins in chemically defended snakes.

Unlike other snakes, most species of Rhabdophis possess glands in their dorsal skin, sometimes limited to the neck, known as nucho-dorsal and nuchal glands, respectively. Those glands contain powerful cardiotonic steroids known as bufadienolides, which can be deployed as a defense against predators. Bufadienolides otherwise occur only in toads (Bufonidae) and some fireflies (Lampyrinae), which are known or believed to synthesize the toxins. The ancestral diet of Rhabdophis consists of anuran amphibians, and we have shown previously that the bufadienolide toxins of frog-eating species are sequestered from toads consumed as prey. However, one derived clade, the Rhabdophis nuchalis Group, has shifted its primary diet from frogs to earthworms. Here we confirm that the worm-eating snakes possess bufadienolides in their nucho-dorsal glands, although the worms themselves lack such toxins. In addition, we show that the bufadienolides of R. nuchalis Group species are obtained primarily from fireflies. Although few snakes feed on insects, we document through feeding experiments, chemosensory preference tests, and gut contents that lampyrine firefly larvae are regularly consumed by these snakes. Furthermore, members of the R. nuchalis Group contain compounds that resemble the distinctive bufadienolides of fireflies, but not those of toads, in stereochemistry, glycosylation, acetylation, and molecular weight. Thus, the evolutionary shift in primary prey among members of the R. nuchalis Group has been accompanied by a dramatic shift in the source of the species' sequestered defensive toxins.

The latitudinal gradient in plant community assembly processes: A meta-analysis.

Beta(ß)-diversity, or site-to-site variation in species composition, generally decreases with increasing latitude, and the underlying processes driving this pattern have been challenging to elucidate because the signals of community assembly processes are scale-dependent. In this meta-analysis, by synthesising the results of 103 studies that were distributed globally and conducted at various spatial scales, we revealed a latitudinal gradient in the detectable assembly processes of vascular plant communities. Variations in plant community composition at low and high latitudes were mainly explained by geographic variables, suggesting that distance decay and dispersal limitations causing spatial aggregation are influential in these regions. In contrast, variation in species composition correlated most strongly with environmental variables at mid-latitudes (20-30°), reflecting the importance of environmental filtering, although this unimodal pattern was not statistically significant. Importantly, our analysis revealed the effects of different spatial scales, such that the correlation with spatial variables was stronger at smaller sampling extents, and environmental variables were more influential at larger sampling extents. We concluded that plant communities are driven by different community assembly processes in distinct biogeographical regions, suggesting that the latitudinal gradient of biodiversity is created by a combination of multiple processes that vary with environmental and species size differences.

Cryptic diversity and phylogeography of the Rhabdophis nuchalis group (Squamata: Colubridae).

Previous studies, have found that the rapid uplift of the Tibetan plateau accelerated the diversification of species. However, there are few relevant biogeographic data for the Colubridae in this region. We conducted a comprehensive study of the Rhabdophis nuchalis Group, which presently contains four nominal species, R. nuchalis, R. pentasupralabialis, R. leonardi, and R. chiwen. Building upon previous studies with specimens we have recently examined, greater interspecific and intraspecific diversity has been revealed. Here we address three questions: (1) Do the intraspecific differences represent only geographic variation within lineages, or are there cryptic species? (2) What are the interspecific relationships among members of the R. nuchalis Group? (3) What has been the biogeographic history of this species group? To resolve these questions we used four mitochondrial gene sequences and one nuclear sequence to investigate the molecular phylogenetic and geographic relationships among populations. Our molecular analysis reveals cryptic species diversity within the R. nuchalis Group, and seven clades were identified in the analysis. Ancestral area estimation suggests that the R. nuchalis Group originated in the Hengduan Mountains approximately 6.24 Mya and expanded its range northward to the Qinling-Daba Mountains. The Sichuan Basin appears to have been a barrier to migration. Species divergence seems to have been related to the rapid uplift of the Qinghai-Tibet Plateau.

A Conceptual Framework to Integrate Biodiversity, Ecosystem Function, and Ecosystem Service Models.

Global biodiversity and ecosystem service models typically operate independently. Ecosystem service projections may therefore be overly optimistic because they do not always account for the role of biodiversity in maintaining ecological functions. We review models used in recent global model intercomparison projects and develop a novel model integration framework to more fully account for the role of biodiversity in ecosystem function, a key gap for linking biodiversity changes to ecosystem services. We propose two integration pathways. The first uses empirical data on biodiversity-ecosystem function relationships to bridge biodiversity and ecosystem function models and could currently be implemented globally for systems and taxa with sufficient data. We also propose a trait-based approach involving greater incorporation of biodiversity into ecosystem function models. Pursuing both approaches will provide greater insight into biodiversity and ecosystem services projections. Integrating biodiversity, ecosystem function, and ecosystem service modeling will enhance policy development to meet global sustainability goals.

Implications of scale dependence for cross-study syntheses of biodiversity differences.

Biodiversity studies are sensitive to well-recognised temporal and spatial scale dependencies. Cross-study syntheses may inflate these influences by collating studies that vary widely in the numbers and sizes of sampling plots. Here we evaluate sources of inaccuracy and imprecision in study-level and cross-study estimates of biodiversity differences, caused by within-study grain and sample sizes, biodiversity measure, and choice of effect-size metric. Samples from simulated communities of old-growth and secondary forests demonstrated influences of all these parameters on the accuracy and precision of cross-study effect sizes. In cross-study synthesis by formal meta-analysis, the metric of log response ratio applied to measures of species richness yielded better accuracy than the commonly used Hedges' g metric on species density, which dangerously combined higher precision with persistent bias. Full-data analyses of the raw plot-scale data using multilevel models were also susceptible to scale-dependent bias. We demonstrate the challenge of detecting scale dependence in cross-study synthesis, due to ubiquitous covariation between replication, variance and plot size. We propose solutions for diagnosing and minimising bias. We urge that empirical studies publish raw data to allow evaluation of covariation in cross-study syntheses, and we recommend against using Hedges' g in biodiversity meta-analyses.

The function-dominance correlation drives the direction and strength of biodiversity-ecosystem functioning relationships.

Community composition is a primary determinant of how biodiversity change influences ecosystem functioning and, therefore, the relationship between biodiversity and ecosystem functioning (BEF). We examine the consequences of community composition across six structurally realistic plant community models. We find that a positive correlation between species' functioning in monoculture versus their dominance in mixture with regard to a specific function (the "function-dominance correlation") generates a positive relationship between realised diversity and ecosystem functioning across species richness treatments. However, because realised diversity declines when few species dominate, a positive function-dominance correlation generates a negative relationship between realised diversity and ecosystem functioning within species richness treatments. Removing seed inflow strengthens the link between the function-dominance correlation and BEF relationships across species richness treatments but weakens it within them. These results suggest that changes in species' identities in a local species pool may more strongly affect ecosystem functioning than changes in species richness.

Grand challenges in biodiversity-ecosystem functioning research in the era of science-policy platforms require explicit consideration of feedbacks.

Feedbacks are an essential feature of resilient socio-economic systems, yet the feedbacks between biodiversity, ecosystem services and human wellbeing are not fully accounted for in global policy efforts that consider future scenarios for human activities and their consequences for nature. Failure to integrate feedbacks in our knowledge frameworks exacerbates uncertainty in future projections and potentially prevents us from realizing the full benefits of actions we can take to enhance sustainability. We identify six scientific research challenges that, if addressed, could allow future policy, conservation and monitoring efforts to quantitatively account for ecosystem and societal consequences of biodiversity change. Placing feedbacks prominently in our frameworks would lead to (i) coordinated observation of biodiversity change, ecosystem functions and human actions, (ii) joint experiment and observation programmes, (iii) more effective use of emerging technologies in biodiversity science and policy, and (iv) a more inclusive and integrated global community of biodiversity observers. To meet these challenges, we outline a five-point action plan for collaboration and connection among scientists and policymakers that emphasizes diversity, inclusion and open access. Efforts to protect biodiversity require the best possible scientific understanding of human activities, biodiversity trends, ecosystem functions and-critically-the feedbacks among them.

New Insights Into Dietary Toxin Metabolism: Diversity in the Ability of the Natricine Snake Rhabdophis tigrinus to Convert Toad-Derived Bufadienolides.

The Japanese natricine snake Rhabdophis tigrinus sequesters cardiotonic steroids, bufadienolides (BDs), from ingested toads in the nuchal glands as defensive toxins. A previous study showed that R. tigrinus in captivity converts dietary BDs when it sequesters them. However, it is unknown whether the dietary BDs are actually converted and the modified products accumulated under natural conditions. It is also unknown to what extent the BD profile of ingested toads is reflected in that of the snake. We collected 123 snakes from throughout Japan, analyzed their BD profiles by liquid chromatography/mass spectrometry, and identified 15 BDs from R. tigrinus by nuclear magnetic resonance analyses. We also compared their BD profiles using hierarchical cluster analysis (HCA). HCA exhibited two main clusters associated with their collection locations: eastern and western regions of the Japanese main islands. These results, coupled with previous findings on the BDs of Japanese toads, suggest that 1) R. tigrinus converts toad-derived BDs into other compounds under natural conditions; 2) there are both universal and regionally-specific conversions of dietary BDs by R. tigrinus; and 3) geographic variation in toad BD profiles is partially reflected in the variation of snake BD profiles.

Biodiversity increases the resistance of ecosystem productivity to climate extremes.

It remains unclear whether biodiversity buffers ecosystems against climate extremes, which are becoming increasingly frequent worldwide. Early results suggested that the ecosystem productivity of diverse grassland plant communities was more resistant, changing less during drought, and more resilient, recovering more quickly after drought, than that of depauperate communities. However, subsequent experimental tests produced mixed results. Here we use data from 46 experiments that manipulated grassland plant diversity to test whether biodiversity provides resistance during and resilience after climate events. We show that biodiversity increased ecosystem resistance for a broad range of climate events, including wet or dry, moderate or extreme, and brief or prolonged events. Across all studies and climate events, the productivity of low-diversity communities with one or two species changed by approximately 50% during climate events, whereas that of high-diversity communities with 16-32 species was more resistant, changing by only approximately 25%. By a year after each climate event, ecosystem productivity had often fully recovered, or overshot, normal levels of productivity in both high- and low-diversity communities, leading to no detectable dependence of ecosystem resilience on biodiversity. Our results suggest that biodiversity mainly stabilizes ecosystem productivity, and productivity-dependent ecosystem services, by increasing resistance to climate events. Anthropogenic environmental changes that drive biodiversity loss thus seem likely to decrease ecosystem stability, and restoration of biodiversity to increase it, mainly by changing the resistance of ecosystem productivity to climate events.