Rapid evolution of SARS-CoV-2 challenges human defenses.

The race between pathogens and their hosts is a major evolutionary driver, where both reshuffle their genomes to overcome and reorganize the defenses for infection, respectively. Evolutionary theory helps formulate predictions on the future evolutionary dynamics of SARS-CoV-2, which can be monitored through unprecedented real-time tracking of SARS-CoV-2 population genomics at the global scale. Here we quantify the accelerating evolution of SARS-CoV-2 by tracking the SARS-CoV-2 mutation globally, with a focus on the Receptor Binding Domain (RBD) of the spike protein determining infection success. We estimate that the > 820 million people that had been infected by October 5, 2021, produced up to 1021 copies of the virus, with 12 new effective RBD variants appearing, on average, daily. Doubling of the number of RBD variants every 89 days, followed by selection of the most infective variants challenges our defenses and calls for a shift to anticipatory, rather than reactive tactics involving collaborative global sequencing and vaccination.

Detection of SARS-CoV-2 variants requires urgent global coordination.

OBJECTIVES: To assess the efforts deployed by different nations and territories in sequencing SARS-CoV-2 isolates, thus enabling detection of variants, known and novel, of concern. METHODS: The sources of over one million full genome sequences of SARS-CoV-2 virus available in the COVID-19 virus Mutation Tracker (CovMT) were analyzed to determine the number of variants in the RBD region of the genome determining infectivity detected in the various nations and territories. RESULTS: The number of detected variants increased as the square root of sequencing effort by nations. Eight nations have contributed 79% of all SARS-CoV-2 isolates that have been sequenced, with two-thirds of all unique variants, adding to 1118 RBD variants, reported by five nations. The median number of sequenced isolates required to detect, on average, one novel RBD variant is 24.05, which is a threshold achieved by 70 nations. CONCLUSIONS: Many developing nations have not contributed any sequences due to lack of capacity. This poses a risk of dangerous virus variants in these under-sampled regions spreading globally before being detected. A collaborative program to sequence SARS-CoV-2 isolates, and other pathogens of concern, is needed to monitor, track, and control the pandemic.

No Evidence for Temperature-Dependence of the COVID-19 Epidemic.

The pandemic of the COVID-19 extended from China across the north-temperate zone, and more recently to the tropics and southern hemisphere. The hypothesis that COVID-19 spread is temperature-dependent was tested based on data derived from nations across the world and provinces in China. No evidence of a pattern between spread rates and ambient temperature was found, suggesting that the SARS-CoV-2 is unlikely to behave as a seasonal respiratory virus.

Multiple stressor effects on coral reef ecosystems.

Global climate change has profound implications on species distributions and ecosystem functioning. In the coastal zone, ecological responses may be driven by various biogeochemical and physical environmental factors. Synergistic interactions can occur when the combined effects of stressors exceed their individual effects. The Red Sea, characterized by strong gradients in temperature, salinity, and nutrients along the latitudinal axis provides a unique opportunity to study ecological responses over a range of these environmental variables. Using multiple linear regression models integrating in situ, satellite and oceanographic data, we investigated the response of coral reef taxa to local stressors and recent climate variability. Taxa and functional groups responded to a combination of climate (temperature, salinity, air-sea heat fluxes, irradiance, wind speed), fishing pressure and biogeochemical (chlorophyll a and nutrients - phosphate, nitrate, nitrite) factors. The regression model for each species showed interactive effects of climate, fishing pressure and nutrient variables. The nature of the effects (antagonistic or synergistic) was dependent on the species and stressor pair. Variables consistently associated with the highest number of synergistic interactions included heat flux terms, temperature, and wind speed followed by fishing pressure. Hard corals and coralline algae abundance were sensitive to changing environmental conditions where synergistic interactions decreased their percentage cover. These synergistic interactions suggest that the negative effects of fishing pressure and eutrophication may exacerbate the impact of climate change on corals. A high number of interactions were also recorded for algae, however for this group, synergistic interactions increased algal abundance. This study is unique in applying regression analysis to multiple environmental variables simultaneously to understand stressor interactions in the field. The observed responses have important implications for understanding climate change impacts on marine ecosystems and whether managing local stressors, such as nutrient enrichment and fishing activities, may help mitigate global drivers of change.

Marine heatwaves reveal coral reef zones susceptible to bleaching in the Red Sea.

As the Earth's temperature continues to rise, coral bleaching events become more frequent. Some of the most affected reef ecosystems are located in poorly monitored waters, and thus, the extent of the damage is unknown. We propose the use of marine heatwaves (MHWs) as a new approach for detecting coral reef zones susceptible to bleaching, using the Red Sea as a model system. Red Sea corals are exceptionally heat-resistant, yet bleaching events have increased in frequency. By applying a strict definition of MHWs on >30 year satellite-derived sea surface temperature observations (1985-2015), we provide an atlas of MHW hotspots over the Red Sea coral reef zones, which includes all MHWs that caused major coral bleaching. We found that: (a) if tuned to a specific set of conditions, MHWs identify all areas where coral bleaching has previously been reported; (b) those conditions extended farther and occurred more often than bleaching was reported; and (c) an emergent pattern of extreme warming events is evident in the northern Red Sea (since 1998), a region until now thought to be a thermal refuge for corals. We argue that bleaching in the Red Sea may be vastly underrepresented. Additionally, although northern Red Sea corals exhibit remarkably high thermal resistance, the rapidly rising incidence of MHWs of high intensity indicates this region may not remain a thermal refuge much longer. As our regionally tuned MHW algorithm was capable of isolating all extreme warming events that have led to documented coral bleaching in the Red Sea, we propose that this approach could be used to reveal bleaching-prone regions in other data-limited tropical regions. It may thus prove a highly valuable tool for policymakers to optimize the sustainable management of coastal economic zones.

Bayesian inference for psychology. Part I: Theoretical advantages and practical ramifications.

Bayesian parameter estimation and Bayesian hypothesis testing present attractive alternatives to classical inference using confidence intervals and p values. In part I of this series we outline ten prominent advantages of the Bayesian approach. Many of these advantages translate to concrete opportunities for pragmatic researchers. For instance, Bayesian hypothesis testing allows researchers to quantify evidence and monitor its progression as data come in, without needing to know the intention with which the data were collected. We end by countering several objections to Bayesian hypothesis testing. Part II of this series discusses JASP, a free and open source software program that makes it easy to conduct Bayesian estimation and testing for a range of popular statistical scenarios (Wagenmakers et al. this issue).

Bayesian inference for psychology. Part II: Example applications with JASP.

Bayesian hypothesis testing presents an attractive alternative to p value hypothesis testing. Part I of this series outlined several advantages of Bayesian hypothesis testing, including the ability to quantify evidence and the ability to monitor and update this evidence as data come in, without the need to know the intention with which the data were collected. Despite these and other practical advantages, Bayesian hypothesis tests are still reported relatively rarely. An important impediment to the widespread adoption of Bayesian tests is arguably the lack of user-friendly software for the run-of-the-mill statistical problems that confront psychologists for the analysis of almost every experiment: the t-test, ANOVA, correlation, regression, and contingency tables. In Part II of this series we introduce JASP ( http://www.jasp-stats.org ), an open-source, cross-platform, user-friendly graphical software package that allows users to carry out Bayesian hypothesis tests for standard statistical problems. JASP is based in part on the Bayesian analyses implemented in Morey and Rouder's BayesFactor package for R. Armed with JASP, the practical advantages of Bayesian hypothesis testing are only a mouse click away.

What Are the Odds? Modern Relevance and Bayes Factor Solutions for MacAlister's Problem From the 1881 Educational Times.

In 1881, Donald MacAlister posed a problem in the Educational Times that remains relevant today. The problem centers on the statistical evidence for the effectiveness of a treatment based on a comparison between two proportions. A brief historical sketch is followed by a discussion of two default Bayesian solutions, one based on a one-sided test between independent rates, and one on a one-sided test between dependent rates. We demonstrate the current-day relevance of MacAlister's original question with a modern-day example about the effectiveness of an educational program.

Default "Gunel and Dickey" Bayes factors for contingency tables.

The analysis of R×C contingency tables usually features a test for independence between row and column counts. Throughout the social sciences, the adequacy of the independence hypothesis is generally evaluated by the outcome of a classical p-value null-hypothesis significance test. Unfortunately, however, the classical p-value comes with a number of well-documented drawbacks. Here we outline an alternative, Bayes factor method to quantify the evidence for and against the hypothesis of independence in R×C contingency tables. First we describe different sampling models for contingency tables and provide the corresponding default Bayes factors as originally developed by Gunel and Dickey (Biometrika, 61(3):545-557 (1974)). We then illustrate the properties and advantages of a Bayes factor analysis of contingency tables through simulations and practical examples. Computer code is available online and has been incorporated in the "BayesFactor" R package and the JASP program ( jasp-stats.org ).

A unimodal species response model relating traits to environment with application to phytoplankton communities.

In this paper we attempt to explain observed niche differences among species (i.e. differences in their distribution along environmental gradients) by differences in trait values (e.g. volume) in phytoplankton communities. For this, we propose the trait-modulated Gaussian logistic model in which the niche parameters (optimum, tolerance and maximum) are made linearly dependent on species traits. The model is fitted to data in the Bayesian framework using OpenBUGS (Bayesian inference Using Gibbs Sampling) to identify according to which environmental variables there is niche differentiation among species and traits. We illustrate the method with phytoplankton community data of 203 lakes located within four climate zones and associated measurements on 11 environmental variables and six morphological species traits of 60 species. Temperature and chlorophyll-a (with opposite signs) described well the niche structure of all species. Results showed that about 25% of the variance in the niche centres with respect to chlorophyll-a were accounted for by traits, whereas niche width and maximum could not be predicted by traits. Volume, mucilage, flagella and siliceous exoskeleton are found to be the most important traits to explain the niche centres. Species were clustered in two groups with different niches structures, group 1 high temperature-low chlorophyll-a species and group 2 low temperature-high chlorophyll-a species. Compared to group 2, species in group 1 had larger volume but lower surface area, had more often flagella but neither mucilage nor siliceous exoskeleton. These results might help in understanding the effect of environmental changes on phytoplankton community. The proposed method, therefore, can also apply to other aquatic or terrestrial communities for which individual traits and environmental conditioning factors are available.

Generalized linear mixed models can detect unimodal species-environment relationships.

Niche theory predicts that species occurrence and abundance show non-linear, unimodal relationships with respect to environmental gradients. Unimodal models, such as the Gaussian (logistic) model, are however more difficult to fit to data than linear ones, particularly in a multi-species context in ordination, with trait modulated response and when species phylogeny and species traits must be taken into account. Adding squared terms to a linear model is a possibility but gives uninterpretable parameters. This paper explains why and when generalized linear mixed models, even without squared terms, can effectively analyse unimodal data and also presents a graphical tool and statistical test to test for unimodal response while fitting just the generalized linear mixed model. The R-code for this is supplied in Supplemental Information 1.