Modelling daisy quorum drive: A short-term bridge across engineered fitness valleys.

Engineered gene-drive techniques for population modification and/or suppression have the potential for tackling complex challenges, including reducing the spread of diseases and invasive species. Gene-drive systems with low threshold frequencies for invasion, such as homing-based gene drive, require initially few transgenic individuals to spread and are therefore easy to introduce. The self-propelled behavior of such drives presents a double-edged sword, however, as the low threshold can allow transgenic elements to expand beyond a target population. By contrast, systems where a high threshold frequency must be reached before alleles can spread-above a fitness valley-are less susceptible to spillover but require introduction at a high frequency. We model a proposed drive system, called "daisy quorum drive," that transitions over time from a low-threshold daisy-chain system (involving homing-based gene drive such as CRISPR-Cas9) to a high-threshold fitness-valley system (requiring a high frequency-a "quorum"-to spread). The daisy-chain construct temporarily lowers the high thresholds required for spread of the fitness-valley construct, facilitating use in a wide variety of species that are challenging to breed and release in large numbers. Because elements in the daisy chain only drive subsequent elements in the chain and not themselves and also carry deleterious alleles ("drive load"), the daisy chain is expected to exhaust itself, removing all CRISPR elements and leaving only the high-threshold fitness-valley construct, whose spread is more spatially restricted. Developing and analyzing both discrete patch and continuous space models, we explore how various attributes of daisy quorum drive affect the chance of modifying local population characteristics and the risk that transgenic elements expand beyond a target area. We also briefly explore daisy quorum drive when population suppression is the goal. We find that daisy quorum drive can provide a promising bridge between gene-drive and fitness-valley constructs, allowing spread from a low frequency in the short term and better containment in the long term, without requiring repeated introductions or persistence of CRISPR elements.

Using early detection data to estimate the date of emergence of an epidemic outbreak.

While the first infection of an emerging disease is often unknown, information on early cases can be used to date it. In the context of the COVID-19 pandemic, previous studies have estimated dates of emergence (e.g., first human SARS-CoV-2 infection, emergence of the Alpha SARS-CoV-2 variant) using mainly genomic data. Another dating attempt used a stochastic population dynamics approach and the date of the first reported case. Here, we extend this approach to use a larger set of early reported cases to estimate the delay from first infection to the Nth case. We first validate our framework by running our model on simulated data. We then apply our model using data on Alpha variant infections in the UK, dating the first Alpha infection at (median) August 21, 2020 (95% interpercentile range across retained simulations (IPR): July 23-September 5, 2020). Next, we apply our model to data on COVID-19 cases with symptom onset before mid-January 2020. We date the first SARS-CoV-2 infection in Wuhan at (median) November 28, 2019 (95% IPR: November 2-December 9, 2019). Our results fall within ranges previously estimated by studies relying on genomic data. Our population dynamics-based modelling framework is generic and flexible, and thus can be applied to estimate the starting time of outbreaks in contexts other than COVID-19.

A stochastic analysis of the interplay between antibiotic dose, mode of action, and bacterial competition in the evolution of antibiotic resistance.

The use of an antibiotic may lead to the emergence and spread of bacterial strains resistant to this antibiotic. Experimental and theoretical studies have investigated the drug dose that minimizes the risk of resistance evolution over the course of treatment of an individual, showing that the optimal dose will either be the highest or the lowest drug concentration possible to administer; however, no analytical results exist that help decide between these two extremes. To address this gap, we develop a stochastic mathematical model of bacterial dynamics under antibiotic treatment. We explore various scenarios of density regulation (bacterial density affects cell birth or death rates), and antibiotic modes of action (biostatic or biocidal). We derive analytical results for the survival probability of the resistant subpopulation until the end of treatment, the size of the resistant subpopulation at the end of treatment, the carriage time of the resistant subpopulation until it is replaced by a sensitive one after treatment, and we verify these results with stochastic simulations. We find that the scenario of density regulation and the drug mode of action are important determinants of the survival of a resistant subpopulation. Resistant cells survive best when bacterial competition reduces cell birth and under biocidal antibiotics. Compared to an analogous deterministic model, the population size reached by the resistant type is larger and carriage time is slightly reduced by stochastic loss of resistant cells. Moreover, we obtain an analytical prediction of the antibiotic concentration that maximizes the survival of resistant cells, which may help to decide which drug dosage (not) to administer. Our results are amenable to experimental tests and help link the within and between host scales in epidemiological models.

Pulled, pushed or failed: the demographic impact of a gene drive can change the nature of its spatial spread.

Understanding the temporal spread of gene drive alleles-alleles that bias their own transmission-through modeling is essential before any field experiments. In this paper, we present a deterministic reaction-diffusion model describing the interplay between demographic and allelic dynamics, in a one-dimensional spatial context. We focused on the traveling wave solutions, and more specifically, on the speed of gene drive invasion (if successful). We considered various timings of gene conversion (in the zygote or in the germline) and different probabilities of gene conversion (instead of assuming 100[Formula: see text] conversion as done in a previous work). We compared the types of propagation when the intrinsic growth rate of the population takes extreme values, either very large or very low. When it is infinitely large, the wave can be either successful or not, and, if successful, it can be either pulled or pushed, in agreement with previous studies (extended here to the case of partial conversion). In contrast, it cannot be pushed when the intrinsic growth rate is vanishing. In this case, analytical results are obtained through an insightful connection with an epidemiological SI model. We conducted extensive numerical simulations to bridge the gap between the two regimes of large and low growth rate. We conjecture that, if it is pulled in the two extreme regimes, then the wave is always pulled, and the wave speed is independent of the growth rate. This occurs for instance when the fitness cost is small enough, or when there is stable coexistence of the drive and the wild-type in the population after successful drive invasion. Our model helps delineate the conditions under which demographic dynamics can affect the spread of a gene drive.

Success of prophylactic antiviral therapy for SARS-CoV-2: Predicted critical efficacies and impact of different drug-specific mechanisms of action.

Repurposed drugs that are safe and immediately available constitute a first line of defense against new viral infections. Despite limited antiviral activity against SARS-CoV-2, several drugs are being tested as medication or as prophylaxis to prevent infection. Using a stochastic model of early phase infection, we evaluate the success of prophylactic treatment with different drug types to prevent viral infection. We find that there exists a critical efficacy that a treatment must reach in order to block viral establishment. Treatment by a combination of drugs reduces the critical efficacy, most effectively by the combination of a drug blocking viral entry into cells and a drug increasing viral clearance. Below the critical efficacy, the risk of infection can nonetheless be reduced. Drugs blocking viral entry into cells or enhancing viral clearance reduce the risk of infection more than drugs that reduce viral production in infected cells. The larger the initial inoculum of infectious virus, the less likely is prevention of an infection. In our model, we find that as long as the viral inoculum is smaller than 10 infectious virus particles, viral infection can be prevented almost certainly with drugs of 90% efficacy (or more). Even when a viral infection cannot be prevented, antivirals delay the time to detectable viral loads. The largest delay of viral infection is achieved by drugs reducing viral production in infected cells. A delay of virus infection flattens the within-host viral dynamic curve, possibly reducing transmission and symptom severity. Thus, antiviral prophylaxis, even with reduced efficacy, could be efficiently used to prevent or alleviate infection in people at high risk.

The French Covid-19 vaccination policy did not solve vaccination inequities: a nationwide study on 64.5 million people.

BACKGROUND: To encourage Covid-19 vaccination, France introduced during the Summer 2021 a 'Sanitary Pass', which morphed into a 'Vaccine Pass' in early 2022. While the sanitary pass led to an increase in Covid-19 vaccination rates, spatial heterogeneities in vaccination rates remained. To identify potential determinants of these heterogeneities and evaluate the French sanitary and vaccine passes' efficacies in reducing them, we used a data-driven approach on exhaustive nationwide data, gathering 141 socio-economic, political and geographic indicators. METHODS: We considered the association between vaccination rates and each indicator at different time points: before the sanitary pass announcement (week 2021-W27), before the sanitary pass came into force (week 2021-W31) and 1 month after (week 2021-W35) and the equivalent dates for the vaccine pass (weeks 2021-W49, 2022-W03 and 2022-W07). RESULTS: The indicators most associated with vaccination rates were the share of local income coming from unemployment benefits, overcrowded households rate, immigrants rate and vote for an 'anti-establishment' candidate at the 2017 Presidential election. These associations increase over time. Consequently, living in a district below the median of such indicator decreases the probability to be vaccinated by about 30% at the end of the studied period, and this probability gradually decreases by deciles of these indicators. CONCLUSIONS: Our analysis reveals that factors related to poverty, immigration and trust in the government are strong determinants of vaccination rate, and that vaccination inequities tended to increase after the introduction of the French sanitary and vaccination passes.

The Effect of Habitat Choice on Evolutionary Rescue in Subdivided Populations.

AbstractEvolutionary rescue is the process by which a population, in response to an environmental change, successfully avoids extinction through adaptation. In spatially structured environments, dispersal can affect the probability of rescue. Here, we model an environment consisting of patches that degrade one after another, and we investigate the probability of rescue by a mutant adapted to the degraded habitat. We focus on the effects of dispersal and of immigration biases. We identify up to three regions delimiting the effect of dispersal on the probability of evolutionary rescue: (i) starting from low dispersal rates, the probability of rescue increases with dispersal; (ii) at intermediate dispersal rates, it decreases; and (iii) at large dispersal rates, it increases again with dispersal, except if mutants are too counterselected in not-yet-degraded patches. The probability of rescue is generally highest when mutant and wild-type individuals preferentially immigrate into patches that have already undergone environmental change. Additionally, we find that mutants that will eventually rescue the population most likely first appear in nondegraded patches. Overall, our results show that habitat choice, compared with the often-studied unbiased immigration scheme, can substantially alter the dynamics of population survival and adaptation to new environments.

Demographic feedbacks can hamper the spatial spread of a gene drive.

This paper is concerned with a reaction-diffusion system modeling the fixation and the invasion in a population of a gene drive (an allele biasing inheritance, increasing its own transmission to offspring). In our model, the gene drive has a negative effect on the fitness of individuals carrying it, and is therefore susceptible of decreasing the total carrying capacity of the population locally in space. This tends to generate an opposing demographic advection that the gene drive has to overcome in order to invade. While previous reaction-diffusion models neglected this aspect, here we focus on it and try to predict the sign of the traveling wave speed. It turns out to be an analytical challenge, only partial results being within reach, and we complete our theoretical analysis by numerical simulations. Our results indicate that taking into account the interplay between population dynamics and population genetics might actually be crucial, as it can effectively reverse the direction of the invasion and lead to failure. Our findings can be extended to other bistable systems, such as the spread of cytoplasmic incompatibilities caused by Wolbachia.

Characterisation of vaccine breakthrough infections of SARS-CoV-2 Delta and Alpha variants and within-host viral load dynamics in the community, France, June to July 2021.

We compared PCR results from SARS-CoV-2-positive patients tested in the community in France from 14 June to 30 July 2021. In asymptomatic individuals, Cq values were significantly higher in fully vaccinated than non-fully vaccinated individuals (effect size: 1.7; 95% CI: 1-2.3; p < 10-6). In symptomatic individuals and controlling for time since symptoms, the difference vanished (p = 0.26). Infections with the Delta variant had lower Cq values at symptom onset than with Alpha (effect size: -3.32; 95% CI: -4.38 to -2.25; p < 10-6).

Detecting the macroevolutionary signal of species interactions.

Species interactions lie at the heart of many theories of macroevolution, from adaptive radiation to the Red Queen. Although some theories describe the imprint that interactions will have over long timescales, we are still missing a comprehensive understanding of the effects of interactions on macroevolution. Current research shows strong evidence for the impact of interactions on macroevolutionary patterns of trait evolution and diversification, yet many macroevolutionary studies have only a tenuous relationship to ecological studies of interactions over shorter timescales. We review current research in this area, highlighting approaches that explicitly model species interactions and connect them to broad-scale macroevolutionary patterns. We also suggest that progress has been made by taking an integrative interdisciplinary look at individual clades. We focus on African cichlids as a case study of how this approach can be fruitful. Overall, although the evidence for species interactions shaping macroevolution is strong, further work using integrative and model-based approaches is needed to spur progress towards understanding the complex dynamics that structure communities over time and space.

Catch Me If You Can: A Spatial Model for a Brake-Driven Gene Drive Reversal.

Population management using artificial gene drives (alleles biasing inheritance, increasing their own transmission to offspring) is becoming a realistic possibility with the development of CRISPR-Cas genetic engineering. A gene drive may, however, have to be stopped. "Antidotes" (brakes) have been suggested, but have been so far only studied in well-mixed populations. Here, we consider a reaction-diffusion system modeling the release of a gene drive (of fitness [Formula: see text]) and a brake (fitness [Formula: see text], [Formula: see text]) in a wild-type population (fitness 1). We prove that whenever the drive fitness is at most 1/2 while the brake fitness is close to 1, coextinction of the brake and the drive occurs in the long run. On the contrary, if the drive fitness is greater than 1/2, then coextinction is impossible: the drive and the brake keep spreading spatially, leaving in the invasion wake a complicated spatiotemporally heterogeneous genetic pattern. Based on numerical experiments, we argue in favor of a global coextinction conjecture provided the drive fitness is at most 1/2, irrespective of the brake fitness. The proof relies upon the study of a related predator-prey system with strong Allee effect on the prey. Our results indicate that some drives may be unstoppable and that if gene drives are ever deployed in nature, threshold drives, that only spread if introduced in high enough frequencies, should be preferred.

Evolutionary dynamics of a quantitative trait in a finite asexual population.

In finite populations, mutation limitation and genetic drift can hinder evolutionary diversification. We consider the evolution of a quantitative trait in an asexual population whose size can vary and depends explicitly on the trait. Previous work showed that evolutionary branching is certain ("deterministic branching") above a threshold population size, but uncertain ("stochastic branching") below it. Using the stationary distribution of the population's trait variance, we identify three qualitatively different sub-domains of "stochastic branching" and illustrate our results using a model of social evolution. We find that in very small populations, branching will almost never be observed; in intermediate populations, branching is intermittent, arising and disappearing over time; in larger populations, finally, branching is expected to occur and persist for substantial periods of time. Our study provides a clearer picture of the ecological conditions that facilitate the appearance and persistence of novel evolutionary lineages in the face of genetic drift.

Evolution of Quantitative Traits under a Migration-Selection Balance: When Does Skew Matter?

Quantitative-genetic models of differentiation under migration-selection balance often rely on the assumption of normally distributed genotypic and phenotypic values. When a population is subdivided into demes with selection toward different local optima, migration between demes may result in asymmetric, or skewed, local distributions. Using a simplified two-habitat model, we derive formulas without a priori assuming a Gaussian distribution of genotypic values, and we find expressions that naturally incorporate higher moments, such as skew. These formulas yield predictions of the expected divergence under migration-selection balance that are more accurate than models assuming Gaussian distributions, which illustrates the importance of incorporating these higher moments to assess the response to selection in heterogeneous environments. We further show with simulations that traits with loci of large effect display the largest skew in their distribution at migration-selection balance.

Host resistance, population structure and the long-term persistence of bubonic plague: contributions of a modelling approach in the Malagasy focus.

Although bubonic plague is an endemic zoonosis in many countries around the world, the factors responsible for the persistence of this highly virulent disease remain poorly known. Classically, the endemic persistence of plague is suspected to be due to the coexistence of plague resistant and plague susceptible rodents in natural foci, and/or to a metapopulation structure of reservoirs. Here, we test separately the effect of each of these factors on the long-term persistence of plague. We analyse the dynamics and equilibria of a model of plague propagation, consistent with plague ecology in Madagascar, a major focus where this disease is endemic since the 1920s in central highlands. By combining deterministic and stochastic analyses of this model, and including sensitivity analyses, we show that (i) endemicity is favoured by intermediate host population sizes, (ii) in large host populations, the presence of resistant rats is sufficient to explain long-term persistence of plague, and (iii) the metapopulation structure of susceptible host populations alone can also account for plague endemicity, thanks to both subdivision and the subsequent reduction in the size of subpopulations, and extinction-recolonization dynamics of the disease. In the light of these results, we suggest scenarios to explain the localized presence of plague in Madagascar.

What we can and cannot learn from SARS-CoV-2 and animals in metagenomic samples from the Huanan market.

While the exact context of the emergence of SARS-CoV-2 remains uncertain, data accumulated since 2020 have provided an increasingly more precise picture of Wuhan's Huanan Seafood Wholesale Market, to which the earliest clusters of human cases of Covid-19 were linked. After the market closed on January 1st 2020, teams from the Chinese Center for Disease Control and Prevention collected environmental samples, and sequenced them. Metagenomic sequencing data from these samples were shared in early 2023. These data confirmed that non-human animals susceptible to SARS-CoV-2 were present in the market before it closed, but also that these animals were located in the side of the market with most human cases, and in a corner with comparatively more SARS-CoV-2-positive environmental samples. The environmental samples were however collected after abundant human-to-human transmission had taken place in the market, precluding any identification of a non-human animal host. Jesse Bloom recently investigated associations between SARS-CoV-2 and non-human animals, concluding that the data failed to indicate whether non-human animals were infected by SARS-CoV-2, despite this being an already acknowledged limitation of the data. Here I explain why a correlation analysis could not confidently conclude which hosts(s) may have shed SARS-CoV-2 in the market, and I rebut the suggestion that such analyses had been encouraged. I show that Bloom's investigation ignores the temporal and spatial structure of the data, which led to incorrect interpretations. Finally, I show that criteria put forward by Bloom to identify the host(s) that shed environmental SARS-CoV-2 would also exclude humans. Progress on the topic of SARS-CoV-2's origin requires a clear distinction between scientific studies and news articles (mis)interpreting them.

Temporal, age, and geographical variation in vaccine efficacy against infection by the Delta and Omicron variants in the community in France, December 2021 to March 2022.

OBJECTIVES: We aimed to quantify how the vaccine efficacy of BNT162b2, messenger RNA-1273, AD26.COV2-S, and ChAdOx1 nCoV-19 against detected infection by the SARS-CoV-2 Delta and Omicron variants varied by time since the last dose, vaccine scheme, age, and geographic areas. METHODS: We analyzed 3,261,749 community polymerase chain reaction tests conducted by private laboratories in France from December 2021 to March 2022 with a test-negative design comparing vaccinated to unvaccinated individuals. RESULTS: Efficacy against detected infection by Delta was 89% (95% confidence interval, 86-91%) at 2 weeks, down to 59% (56-61%) at 26 weeks and more after the second dose. Efficacy against Omicron was 48% (45-51%) at 2 weeks, down to 4% (2-5%) at 16 weeks after the second dose. A third dose temporarily restored efficacy. Efficacy against Omicron was lower in children and the elderly. Geographical variability in efficacy may reflect variability in the ratio of the number of contacts of vaccinated vs unvaccinated individuals. This ratio ranged from 0 to +50% across departments and correlated with the number of restaurants and bars per inhabitant (beta = 15.0 [0.75-29], P-value = 0.04), places that only vaccinated individuals could access in the study period. CONCLUSION: SARS-CoV-2 vaccines conferred low and transient protection against Omicron infection.

Genetic tracing of market wildlife and viruses at the epicenter of the COVID-19 pandemic.

Zoonotic spillovers of viruses have occurred through the animal trade worldwide. The start of the COVID-19 pandemic was traced epidemiologically to the Huanan Wholesale Seafood Market, the site with the most reported wildlife vendors in the city of Wuhan, China. Here, we analyze publicly available qPCR and sequencing data from environmental samples collected in the Huanan market in early 2020. We demonstrate that the SARS-CoV-2 genetic diversity linked to this market is consistent with market emergence, and find increased SARS-CoV-2 positivity near and within a particular wildlife stall. We identify wildlife DNA in all SARS-CoV-2 positive samples from this stall. This includes species such as civets, bamboo rats, porcupines, hedgehogs, and one species, raccoon dogs, known to be capable of SARS-CoV-2 transmission. We also detect other animal viruses that infect raccoon dogs, civets, and bamboo rats. Combining metagenomic and phylogenetic approaches, we recover genotypes of market animals and compare them to those from other markets. This analysis provides the genetic basis for a short list of potential intermediate hosts of SARS-CoV-2 to prioritize for retrospective serological testing and viral sampling.