Tropical tree mortality has increased with rising atmospheric water stress.

Evidence exists that tree mortality is accelerating in some regions of the tropics1,2, with profound consequences for the future of the tropical carbon sink and the global anthropogenic carbon budget left to limit peak global warming below 2 °C. However, the mechanisms that may be driving such mortality changes and whether particular species are especially vulnerable remain unclear3-8. Here we analyse a 49-year record of tree dynamics from 24 old-growth forest plots encompassing a broad climatic gradient across the Australian moist tropics and find that annual tree mortality risk has, on average, doubled across all plots and species over the last 35 years, indicating a potential halving in life expectancy and carbon residence time. Associated losses in biomass were not offset by gains from growth and recruitment. Plots in less moist local climates presented higher average mortality risk, but local mean climate did not predict the pace of temporal increase in mortality risk. Species varied in the trajectories of their mortality risk, with the highest average risk found nearer to the upper end of the atmospheric vapour pressure deficit niches of species. A long-term increase in vapour pressure deficit was evident across the region, suggesting that thresholds involving atmospheric water stress, driven by global warming, may be a primary cause of increasing tree mortality in moist tropical forests.

Key Substitutions in the Spike Protein of SARS-CoV-2 Variants Can Predict Resistance to Monoclonal Antibodies, but Other Substitutions Can Modify the Effects.

Mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants can compromise the effectiveness of therapeutic antibodies. Most clinical-stage therapeutic antibodies target the spike receptor binding domain (RBD), but variants often have multiple mutations in several spike regions. To help predict antibody potency against emerging variants, we evaluated 25 clinical-stage therapeutic antibodies for neutralization activity against 60 pseudoviruses bearing spikes with single or multiple substitutions in several spike domains, including the full set of substitutions in B.1.1.7 (alpha), B.1.351 (beta), P.1 (gamma), B.1.429 (epsilon), B.1.526 (iota), A.23.1, and R.1 variants. We found that 14 of 15 single antibodies were vulnerable to at least one RBD substitution, but most combination and polyclonal therapeutic antibodies remained potent. Key substitutions in variants with multiple spike substitutions predicted resistance, but the degree of resistance could be modified in unpredictable ways by other spike substitutions that may reside outside the RBD. These findings highlight the importance of assessing antibody potency in the context of all substitutions in a variant and show that epistatic interactions in spike can modify virus susceptibility to therapeutic antibodies. IMPORTANCE Therapeutic antibodies are effective in preventing severe disease from SARS-CoV-2 infection (COVID-19), but their effectiveness may be reduced by virus variants with mutations affecting the spike protein. To help predict resistance to therapeutic antibodies in emerging variants, we profiled resistance patterns of 25 antibody products in late stages of clinical development against a large panel of variants that include single and multiple substitutions found in the spike protein. We found that the presence of a key substitution in variants with multiple spike substitutions can predict resistance against a variant but that other substitutions can affect the degree of resistance in unpredictable ways. These findings highlight complex interactions among substitutions in the spike protein affecting virus neutralization and, potentially, virus entry into cells.

Tropical tree growth sensitivity to climate is driven by species intrinsic growth rate and leaf traits.

A better understanding of how climate affects growth in tree species is essential for improved predictions of forest dynamics under climate change. Long-term climate averages (mean climate) drive spatial variations in species' baseline growth rates, whereas deviations from these averages over time (anomalies) can create growth variation around the local baseline. However, the rarity of long-term tree census data spanning climatic gradients has so far limited our understanding of their respective role, especially in tropical systems. Furthermore, tree growth sensitivity to climate is likely to vary widely among species, and the ecological strategies underlying these differences remain poorly understood. Here, we utilize an exceptional dataset of 49 years of growth data for 509 tree species across 23 tropical rainforest plots along a climatic gradient to examine how multiannual tree growth responds to both climate means and anomalies, and how species' functional traits mediate these growth responses to climate. We show that anomalous increases in atmospheric evaporative demand and solar radiation consistently reduced tree growth. Drier forests and fast-growing species were more sensitive to water stress anomalies. In addition, species traits related to water use and photosynthesis partly explained differences in growth sensitivity to both climate means and anomalies. Our study demonstrates that both climate means and anomalies shape tree growth in tropical forests and that species traits can provide insights into understanding these demographic responses to climate change, offering a promising way forward to forecast tropical forest dynamics under different climate trajectories.

Climate shapes and shifts functional biodiversity in forests worldwide.

Much ecological research aims to explain how climate impacts biodiversity and ecosystem-level processes through functional traits that link environment with individual performance. However, the specific climatic drivers of functional diversity across space and time remain unclear due largely to limitations in the availability of paired trait and climate data. We compile and analyze a global forest dataset using a method based on abundance-weighted trait moments to assess how climate influences the shapes of whole-community trait distributions. Our approach combines abundance-weighted metrics with diverse climate factors to produce a comprehensive catalog of trait-climate relationships that differ dramatically-27% of significant results change in sign and 71% disagree on sign, significance, or both-from traditional species-weighted methods. We find that (i) functional diversity generally declines with increasing latitude and elevation, (ii) temperature variability and vapor pressure are the strongest drivers of geographic shifts in functional composition and ecological strategies, and (iii) functional composition may currently be shifting over time due to rapid climate warming. Our analysis demonstrates that climate strongly governs functional diversity and provides essential information needed to predict how biodiversity and ecosystem function will respond to climate change.

Solar radiation and functional traits explain the decline of forest primary productivity along a tropical elevation gradient.

One of the major challenges in ecology is to understand how ecosystems respond to changes in environmental conditions, and how taxonomic and functional diversity mediate these changes. In this study, we use a trait-spectra and individual-based model, to analyse variation in forest primary productivity along a 3.3 km elevation gradient in the Amazon-Andes. The model accurately predicted the magnitude and trends in forest productivity with elevation, with solar radiation and plant functional traits (leaf dry mass per area, leaf nitrogen and phosphorus concentration, and wood density) collectively accounting for productivity variation. Remarkably, explicit representation of temperature variation with elevation was not required to achieve accurate predictions of forest productivity, as trait variation driven by species turnover appears to capture the effect of temperature. Our semi-mechanistic model suggests that spatial variation in traits can potentially be used to estimate spatial variation in productivity at the landscape scale.

Scale dependence of canopy trait distributions along a tropical forest elevation gradient.

Average responses of forest foliar traits to elevation are well understood, but far less is known about trait distributional responses to elevation at multiple ecological scales. This limits our understanding of the ecological scales at which trait variation occurs in response to environmental drivers and change. We analyzed and compared multiple canopy foliar trait distributions using field sampling and airborne imaging spectroscopy along an Andes-to-Amazon elevation gradient. Field-estimated traits were generated from three community-weighting methods, and remotely sensed estimates of traits were made at three scales defined by sampling grain size and ecological extent. Field and remote sensing approaches revealed increases in average leaf mass per unit area (LMA), water, nonstructural carbohydrates (NSCs) and polyphenols with increasing elevation. Foliar nutrients and photosynthetic pigments displayed little to no elevation trend. Sample weighting approaches had little impact on field-estimated trait responses to elevation. Plot representativeness of trait distributions at landscape scales decreased with increasing elevation. Remote sensing indicated elevation-dependent increases in trait variance and distributional skew. Multiscale invariance of LMA, leaf water and NSC mark these traits as candidates for tracking forest responses to changing climate. Trait-based ecological studies can be greatly enhanced with multiscale studies made possible by imaging spectroscopy.

Variation in leaf wettability traits along a tropical montane elevation gradient.

Leaf wetting is often considered to have negative effects on plant function, such that wet environments may select for leaves with certain leaf surface, morphological, and architectural traits that reduce leaf wettability. However, there is growing recognition that leaf wetting can have positive effects. We measured variation in two traits, leaf drip tips and leaf water repellency, in a series of nine tropical forest communities occurring along a 3300-m elevation gradient in southern Peru. To extend this climatic gradient, we also assembled published leaf water repellency values from 17 additional sites. We then tested hypotheses for how these traits should vary as a function of climate. Contrary to expectations, we found that the proportion of species with drip tips did not increase with increasing precipitation. Instead, drip tips increased with increasing temperature. Moreover, leaf water repellency was very low in our sites and the global analysis indicated high repellency only in sites with low precipitation and temperatures. Our findings suggest that drip tips and repellency may not solely reflect the negative effects of wetting on plant function. Understanding the drivers of leaf wettability traits can provide insight into the effects of leaf wetting on plant, community, and ecosystem function.

Predicting trait-environment relationships for venation networks along an Andes-Amazon elevation gradient.

Understanding functional trait-environment relationships (TERs) may improve predictions of community assembly. However, many empirical TERs have been weak or lacking conceptual foundation. TERs based on leaf venation networks may better link individuals and communities via hydraulic constraints. We report measurements of vein density, vein radius, and leaf thickness for more than 100 dominant species occurring in ten forest communities spanning a 3,300 m Andes-Amazon elevation gradient in Peru. We use these data to measure the strength of TERs at community scale and to determine whether observed TERs are similar to those predicted by physiological theory. We found strong support for TERs between all traits and temperature, as well weaker support for a predicted TER between maximum abundance-weighted leaf transpiration rate and maximum potential evapotranspiration. These results provide one approach for developing a more mechanistic trait-based community assembly theory.

Quantifying ecological memory in plant and ecosystem processes.

The role of time in ecology has a long history of investigation, but ecologists have largely restricted their attention to the influence of concurrent abiotic conditions on rates and magnitudes of important ecological processes. Recently, however, ecologists have improved their understanding of ecological processes by explicitly considering the effects of antecedent conditions. To broadly help in studying the role of time, we evaluate the length, temporal pattern, and strength of memory with respect to the influence of antecedent conditions on current ecological dynamics. We developed the stochastic antecedent modelling (SAM) framework as a flexible analytic approach for evaluating exogenous and endogenous process components of memory in a system of interest. We designed SAM to be useful in revealing novel insights promoting further study, illustrated in four examples with different degrees of complexity and varying time scales: stomatal conductance, soil respiration, ecosystem productivity, and tree growth. Models with antecedent effects explained an additional 18-28% of response variation compared to models without antecedent effects. Moreover, SAM also enabled identification of potential mechanisms that underlie components of memory, thus revealing temporal properties that are not apparent from traditional treatments of ecological time-series data and facilitating new hypothesis generation and additional research.

Deviation from symmetrically self-similar branching in trees predicts altered hydraulics, mechanics, light interception and metabolic scaling.

The West, Brown, Enquist (WBE) model derives symmetrically self-similar branching to predict metabolic scaling from hydraulic conductance, K, (a metabolism proxy) and tree mass (or volume, V). The original prediction was Kα V(0.75). We ask whether trees differ from WBE symmetry and if it matters for plant function and scaling. We measure tree branching and model how architecture influences K, V, mechanical stability, light interception and metabolic scaling. We quantified branching architecture by measuring the path fraction, Pf : mean/maximum trunk-to-twig pathlength. WBE symmetry produces the maximum, Pf = 1.0. We explored tree morphospace using a probability-based numerical model constrained only by biomechanical principles. Real tree Pf ranged from 0.930 (nearly symmetric) to 0.357 (very asymmetric). At each modeled tree size, a reduction in Pf led to: increased K; decreased V; increased mechanical stability; and decreased light absorption. When Pf was ontogenetically constant, strong asymmetry only slightly steepened metabolic scaling. The Pf ontogeny of real trees, however, was 'U' shaped, resulting in size-dependent metabolic scaling that exceeded 0.75 in small trees before falling below 0.65. Architectural diversity appears to matter considerably for whole-tree hydraulics, mechanics, photosynthesis and potentially metabolic scaling. Optimal architectures likely exist that maximize carbon gain per structural investment.

Quantifying the timescales over which exogenous and endogenous conditions affect soil respiration.

Understanding how exogenous and endogenous factors and above-ground-below-ground linkages modulate carbon dynamics is difficult because of the influences of antecedent conditions. For example, there are variable lags between above-ground assimilation and below-ground efflux, and the duration of antecedent periods are often arbitrarily assigned. Nonetheless, developing models linking above- and below-ground processes is crucial for estimating current and future carbon dynamics. We collected data on leaf-level photosynthesis (Asat ) and soil respiration (Rsoil ) in different microhabitats (under shrubs vs under bunchgrasses) in the Sonoran Desert. We evaluated timescales over which endogenous and exogenous factors control Rsoil by analyzing data in the context of a semimechanistic temperature-response model of Rsoil that incorporated effects of antecedent exogenous (soil water) and endogenous (Asat ) conditions. For both microhabitats, antecedent soil water and Asat significantly affected Rsoil , but Rsoil under shrubs was more sensitive to Asat than that under bunchgrasses. Photosynthetic rates 1 and 3 d before the Rsoil measurement were most important in determining current-day Rsoil under bunchgrasses and shrubs, respectively, indicating a significant lag effect. Endogenous and exogenous controls are critical drivers of Rsoil , but the relative importance and the timescale over which each factor affects Rsoil depends on above-ground vegetation and ecosystem structure characteristics.

Inclusion of vein traits improves predictive power for the leaf economic spectrum: a response to Sack et al. (2013).

Our model for the worldwide leaf economics spectrum (LES) based on venation networks (Blonder et al., 2011, 2013) was strongly criticized by Sack et al. (2013) in this journal. Here, we show that the majority of criticisms by Sack et al. are based on mathematical and conceptual misunderstandings. Using empirical data from both our original study as well as others in the literature, we show support for our original hypothesis, that venation networks provide predictive power and conceptual unification for the LES. In an effort to reconcile differing viewpoints related to the role of leaf venation traits for the LES, we highlight several lines of further investigation.

An empirical assessment of tree branching networks and implications for plant allometric scaling models.

Several theories predict whole-tree function on the basis of allometric scaling relationships assumed to emerge from traits of branching networks. To test this key assumption, and more generally, to explore patterns of external architecture within and across trees, we measure branch traits (radii/lengths) and calculate scaling exponents from five functionally divergent species. Consistent with leading theories, including metabolic scaling theory, branching is area preserving and statistically self-similar within trees. However, differences among scaling exponents calculated at node- and whole-tree levels challenge the assumption of an optimised, symmetrically branching tree. Furthermore, scaling exponents estimated for branch length change across branching orders, and exponents for scaling metabolic rate with plant size (or number of terminal tips) significantly differ from theoretical predictions. These findings, along with variability in the scaling of branch radii being less than for branch lengths, suggest extending current scaling theories to include asymmetrical branching and differential selective pressures in plant architectures.

Differential daytime and night-time stomatal behavior in plants from North American deserts.

Night-time stomatal conductance (g(night)) occurs in many ecosystems, but the g(night) response to environmental drivers is relatively unknown, especially in deserts. Here, we conducted a Bayesian analysis of stomatal conductance (g) (N=5013) from 16 species in the Sonoran, Chihuahuan, Mojave and Great Basin Deserts (North America). We partitioned daytime g (g(day)) and g(night) responses by describing g as a mixture of two extreme (dark vs high light) behaviors. Significant g(night) was observed across 15 species, and the g(night) and g(day) behavior differed according to species, functional type and desert. The transition between extreme behaviors was determined by light environment, with the transition behavior differing between functional types and deserts. Sonoran and Chihuahuan C(4) grasses were more sensitive to vapor pressure difference (D) at night and soil water potential (Ψ(soil)) during the day, Great Basin C(3) shrubs were highly sensitive to D and Ψ(soil) during the day, and Mojave C(3) shrubs were equally sensitive to D and Ψ(soil) during the day and night. Species were split between the exhibition of isohydric or anisohydric behavior during the day. Three species switched from anisohydric to isohydric behavior at night. Such behavior, combined with differential D, Ψ(soil) and light responses, suggests that different mechanisms underlie g(day) and g(night) regulation.

SARS-CoV-2 BA.1 variant is neutralized by vaccine booster-elicited serum but evades most convalescent serum and therapeutic antibodies.

The rapid spread of the highly contagious Omicron variant of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) along with its high number of mutations in the spike gene has raised alarms about the effectiveness of current medical countermeasures. To address this concern, we measured the neutralization of the Omicron BA.1 variant pseudovirus by postvaccination serum samples after two and three immunizations with the Pfizer/BioNTech162b2 SARS-CoV-2 mRNA (Pfizer/BNT162b2) vaccine, convalescent serum samples from unvaccinated individuals infected by different variants, and clinical-stage therapeutic antibodies. We found that titers against the Omicron variant were low or undetectable after two immunizations and in many convalescent serum samples, regardless of the infecting variant. A booster vaccination increased titers more than 30-fold against Omicron to values comparable to those seen against the D614G variant after two immunizations. Neither age nor sex was associated with the differences in postvaccination antibody responses. We also evaluated 18 clinical-stage therapeutic antibody products and an antibody mimetic protein product obtained directly from the manufacturers. Five monoclonal antibodies, the antibody mimetic protein, three antibody cocktails, and two polyclonal antibody preparations retained measurable neutralization activity against Omicron with a varying degree of potency. Of these, only three retained potencies comparable to the D614G variant. Two therapeutic antibody cocktails in the tested panel that are authorized for emergency use in the United States did not neutralize Omicron. These findings underscore the potential benefit of mRNA vaccine boosters for protection against Omicron and the need for rapid development of antibody therapeutics that maintain potency against emerging variants.

Functional susceptibility of tropical forests to climate change.

Tropical forests are some of the most biodiverse ecosystems in the world, yet their functioning is threatened by anthropogenic disturbances and climate change. Global actions to conserve tropical forests could be enhanced by having local knowledge on the forests' functional diversity and functional redundancy as proxies for their capacity to respond to global environmental change. Here we create estimates of plant functional diversity and redundancy across the tropics by combining a dataset of 16 morphological, chemical and photosynthetic plant traits sampled from 2,461 individual trees from 74 sites distributed across four continents together with local climate data for the past half century. Our findings suggest a strong link between climate and functional diversity and redundancy with the three trait groups responding similarly across the tropics and climate gradient. We show that drier tropical forests are overall less functionally diverse than wetter forests and that functional redundancy declines with increasing soil water and vapour pressure deficits. Areas with high functional diversity and high functional redundancy tend to better maintain ecosystem functioning, such as aboveground biomass, after extreme weather events. Our predictions suggest that the lower functional diversity and lower functional redundancy of drier tropical forests, in comparison with wetter forests, may leave them more at risk of shifting towards alternative states in face of further declines in water availability across tropical regions.

Venation networks and the origin of the leaf economics spectrum.

The leaf economics spectrum describes biome-invariant scaling functions for leaf functional traits that relate to global primary productivity and nutrient cycling. Here, we develop a comprehensive framework for the origin of this leaf economics spectrum based on venation-mediated economic strategies. We define a standardized set of traits - density, distance and loopiness - that provides a common language for the study of venation. We develop a novel quantitative model that uses these venation traits to model leaf-level physiology, and show that selection to optimize the venation network predicts the mean global trait-trait scaling relationships across 2548 species. Furthermore, using empirical venation data for 25 plant species, we test our model by predicting four key leaf functional traits related to leaf economics: net carbon assimilation rate, life span, leaf mass per area ratio and nitrogen content. Together, these results indicate that selection on venation geometry is a fundamental basis for understanding the diversity of leaf form and function, and the carbon balance of leaves. The model and associated predictions have broad implications for integrating venation network geometry with pattern and process in ecophysiology, ecology and palaeobotany.

Greater seasonal carbon gain across a broad temperature range contributes to the invasive potential of Phalaris arundinacea (Poaceae; reed canary grass) over the native sedge Carex stricta (Cyperaceae).

PREMISE OF THE STUDY: Most invasive plants grow faster and produce more biomass than the species that they displace, but physiological mechanisms leading to invasive success are poorly understood. To foster novel control approaches, our goal was to determine whether the grass Phalaris arundinacea possessed superior physiological strategies that contributed to its success over native sedges. METHODS: Data for spring, summer, and autumn diel gas-exchange, leaf morphology, and nitrogen content for plants of P. arundinacea and Carex stricta in water-saturated, drained, and periodically flooded sites in northern Indiana, USA, were compared with similar data for plants in a greenhouse. KEY RESULTS: Phalaris arundinacea had higher maximum CO(2) assimilation (A) across a broad range of temperatures, greater summer/autumn net carbon gain, higher water use efficiencies, larger leaf areas per shoot, and higher specific leaf areas than did C. stricta. Species differences in gas-exchange data were similar in the greenhouse. However, long-term flooding reduced A for P. arundinacea. Greater declines in leaf A and nitrogen content from July to October compared to P. arundinacea were suggestive of earlier leaf senescence for C. stricta. CONCLUSIONS: We propose that superior daily and seasonal carbon gain, especially during rhizome carbohydrate storage in the summer and autumn, contribute to the success of invasive P. arundinacea over C. stricta. This advantage may be enhanced by frequent summer/autumn heat waves. The poor performance of P. arundinacea during long-term flooding is consistent with C. stricta's dominance in water-saturated soil, implying that water management strategies could be crucial to controlling P. arundinacea.

SARS-CoV-2 Delta Variant Displays Moderate Resistance to Neutralizing Antibodies and Spike Protein Properties of Higher Soluble ACE2 Sensitivity, Enhanced Cleavage and Fusogenic Activity.

The SARS-CoV-2 B.1.617 lineage variants, Kappa (B.1.617.1) and Delta (B.1.617.2, AY) emerged during the second wave of infections in India, but the Delta variants have become dominant worldwide and continue to evolve. Here, we compared B.1.617 variants for neutralization resistance by convalescent sera, mRNA vaccine-elicited sera, and therapeutic neutralizing antibodies using a pseudovirus neutralization assay. B.1.617.1, B.1.617.2, and AY.1 pseudoviruses showed a modest 1.5- to 4.4-fold reduction in neutralization by convalescent sera and vaccine-elicited sera. In comparison, similar modest reductions were also observed for C.37, P.1, R.1, and B.1.526 pseudoviruses, but 7- and 16-fold reductions for vaccine-elicited and convalescent sera, respectively, were seen for B.1.351 pseudoviruses. Among twenty-three therapeutic antibodies tested, four antibodies showed either complete or partial loss of neutralization against B.1.617.2 pseudoviruses and six antibodies showed either complete or partial loss of neutralization against B.1.617.1 and AY.1 pseudoviruses. Our results indicate that the current mRNA-based vaccines will likely remain effective in protecting against B.1.617 variants. Finally, the P681R substitution confers efficient cleavage of B.1.617 variants' spike proteins and the spike of Delta variants exhibited greater sensitivity to soluble ACE2 neutralization, as well as fusogenic activity, which may contribute to enhanced spread of Delta variants.

SARS-COV-2 Delta variant displays moderate resistance to neutralizing antibodies and spike protein properties of higher soluble ACE2 sensitivity, enhanced cleavage and fusogenic activity.

The SARS-CoV-2 B.1.617 lineage variants, Kappa (B.1.617.1) and Delta (B.1.617.2, AY) emerged during the second wave of infections in India, but the Delta variants have become dominant worldwide and continue to evolve. The spike proteins of B.1.617.1, B.1.617.2, and AY.1 variants have several substitutions in the receptor binding domain (RBD), including L452R+E484Q, L452R+T478K, and K417N+L452R+T478K, respectively, that could potentially reduce effectiveness of therapeutic antibodies and current vaccines. Here we compared B.1.617 variants, and their single and double RBD substitutions for resistance to neutralization by convalescent sera, mRNA vaccine-elicited sera, and therapeutic neutralizing antibodies using a pseudovirus neutralization assay. Pseudoviruses with the B.1.617.1, B.1.617.2, and AY.1 spike showed a modest 1.5 to 4.4-fold reduction in neutralization titer by convalescent sera and vaccine-elicited sera. In comparison, similar modest reductions were also observed for pseudoviruses with C.37, P.1, R.1, and B.1.526 spikes, but seven- and sixteen-fold reduction for vaccine-elicited and convalescent sera, respectively, was seen for pseudoviruses with the B.1.351 spike. Four of twenty-three therapeutic neutralizing antibodies showed either complete or partial loss of neutralization against B.1.617.2 pseudoviruses due to the L452R substitution, whereas six of twenty-three therapeutic neutralizing antibodies showed either complete or partial loss of neutralization against B.1.617.1 pseudoviruses due to either the E484Q or L452R substitution. Against AY.1 pseudoviruses, the L452R and K417N substitutions accounted for the loss of neutralization by four antibodies and one antibody, respectively, whereas one antibody lost potency that could not be fully accounted for by a single RBD substitution. The modest resistance of B.1.617 variants to vaccine-elicited sera suggest that current mRNA-based vaccines will likely remain effective in protecting against B.1.617 variants, but the therapeutic antibodies need to be carefully selected based on their resistance profiles. Finally, the spike proteins of B.1.617 variants are more efficiently cleaved due to the P681R substitution, and the spike of Delta variants exhibited greater sensitivity to soluble ACE2 neutralization, as well as fusogenic activity, which may contribute to enhanced spread of Delta variants.