Fixation times of de novo and standing beneficial variants in subdivided populations.

The rate at which beneficial alleles fix in a population depends on the probability of and time to fixation of such alleles. Both of these quantities can be significantly impacted by population subdivision and limited gene flow. Here, we investigate how limited dispersal influences the rate of fixation of beneficial de novo mutations, as well as fixation time from standing genetic variation. We investigate this for a population structured according to the island model of dispersal allowing us to use the diffusion approximation, which we complement with simulations. We find that fixation may take on average fewer generations under limited dispersal than under panmixia when selection is moderate. This is especially the case if adaptation occurs from de novo recessive mutations, and dispersal is not too limited (such that approximately FST<0.2). The reason is that mildly limited dispersal leads to only a moderate increase in effective population size (which slows down fixation), but is sufficient to cause a relative excess of homozygosity due to inbreeding, thereby exposing rare recessive alleles to selection (which accelerates fixation). We also explore the effect of metapopulation dynamics through local extinction followed by recolonization, finding that such dynamics always accelerate fixation from standing genetic variation, while de novo mutations show faster fixation interspersed with longer waiting times. Finally, we discuss the implications of our results for the detection of sweeps, suggesting that limited dispersal mitigates the expected differences between the genetic signatures of sweeps involving recessive and dominant alleles.

Assessing the best time interval between doses in a two-dose vaccination regimen to reduce the number of deaths in an ongoing epidemic of SARS-CoV-2.

The SARS-CoV-2 pandemic is a major concern all over the world and, as vaccines became available at the end of 2020, optimal vaccination strategies were subjected to intense investigation. Considering their critical role in reducing disease burden, the increasing demand outpacing production, and that most currently approved vaccines follow a two-dose regimen, the cost-effectiveness of delaying the second dose to increment the coverage of the population receiving the first dose is often debated. Finding the best solution is complex due to the trade-off between vaccinating more people with lower level of protection and guaranteeing higher protection to a fewer number of individuals. Here we present a novel extended age-structured SEIR mathematical model that includes a two-dose vaccination schedule with a between-doses delay modelled through delay differential equations and linear optimization of vaccination rates. By maintaining the minimum stock of vaccines under a given production rate, we evaluate the dose interval that minimizes the number of deaths. We found that the best strategy depends on an interplay between the vaccine production rate and the relative efficacy of the first dose. In the scenario of low first-dose efficacy, it is always better to vaccinate the second dose as soon as possible, while for high first-dose efficacy, the best strategy of time window depends on the production rate and also on second-dose efficacy provided by each type of vaccine. We also found that the rate of spread of the infection does not affect significantly the thresholds of the best window, but is an important factor in the absolute number of total deaths. These conclusions point to the need to carefully take into account both vaccine characteristics and roll-out speed to optimize the outcome of vaccination strategies.

Percolation across households in mechanistic models of non-pharmaceutical interventions in SARS-CoV-2 disease dynamics.

Since the emergence of the novel coronavirus disease 2019 (COVID-19), mathematical modelling has become an important tool for planning strategies to combat the pandemic by supporting decision-making and public policies, as well as allowing an assessment of the effect of different intervention scenarios. A proliferation of compartmental models were developed by the mathematical modelling community in order to understand and make predictions about the spread of COVID-19. While compartmental models are suitable for simulating large populations, the underlying assumption of a well-mixed population might be problematic when considering non-pharmaceutical interventions (NPIs) which have a major impact on the connectivity between individuals in a population. Here we propose a modification to an extended age-structured SEIR (susceptible-exposed-infected-recovered) framework, with dynamic transmission modelled using contact matrices for various settings in Brazil. By assuming that the mitigation strategies for COVID-19 affect the connections among different households, network percolation theory predicts that the connectivity among all households decreases drastically above a certain threshold of removed connections. We incorporated this emergent effect at population level by modulating home contact matrices through a percolation correction function, with the few additional parameters fitted to hospitalisation and mortality data from the city of São Paulo. Our model with percolation effects was better supported by the data than the same model without such effects. By allowing a more reliable assessment of the impact of NPIs, our improved model provides a better description of the epidemiological dynamics and, consequently, better policy recommendations.

Magnetoreception in multicellular magnetotactic prokaryotes: a new analysis of escape motility trajectories in different magnetic fields.

Magnetotactic microorganisms can be found as unicellular prokaryotes, as cocci, vibrions, spirilla and rods, and as multicellular organisms. Multicellular magnetotactic prokaryotes are magnetotactic microorganisms composed by several magnetotactic bacteria organized almost in a spherical helix, and one of the most studied is Candidatus Magnetoglobus multicellularis. Several studies have shown that Ca. M. multicellularis displays forms of behavior not well explained by magnetotaxis. One of these is escape motility, also known as "ping-pong" motion. Studies done in the past associated the "ping-pong" motion to some magnetoreceptive behavior, but those studies were never replicated. In the present manuscript a characterization of escape motility trajectories of Ca. M. multicellularis was done for several magnetic fields, considering that this microorganism swims in cylindrical helical trajectories. It was observed that the escape motility can be separated into three phases: (I) when the microorganism jumps from the drop border, (II) where the microorganism moves almost perpendicular to the magnetic field and (III) when the microorganism returns to the drop border. The total time of the whole escape motility, the time spent in phase II and the displacement distance in phase I decreases when the magnetic field increases. Our results show that the escape motility has several characteristics that depend on the magnetic field and cannot be understood by magnetotaxis, with a magnetoreceptive mechanism being the best explanation.