Phylodynamic theory of persistence, extinction and speciation of rapidly adapting pathogens
Abstract
Rapidly evolving pathogens like influenza viruses can persist by changing their antigenic properties fast enough to evade the adaptive immunity, yet they rarely split into diverging lineages. By mapping the multi-strain Susceptible-Infected-Recovered model onto the traveling wave model of adapting populations, we demonstrate that persistence of a rapidly evolving, Red-Queen-like state of the pathogen population requires long-ranged cross-immunity and sufficiently large population sizes. This state is unstable and the population goes extinct or 'speciates' into two pathogen strains with antigenic divergence beyond the range of cross-inhibition. However, in a certain range of evolutionary parameters, a single cross-inhibiting population can exist for times long compared to the time to the most recent common ancestor (TMRCA) and gives rise to phylogenetic patterns typical of influenza virus. We demonstrate that the rate of speciation is related to fluctuations of TMRCA and construct a 'phase diagram' identifying different phylodynamic regimes as a function of evolutionary parameters.
Data availability
Computer programs used for numerical simulations and analysis have been made publicly available athttps://github.com/neherlab/FluSpeciation
Article and author information
Author details
Funding
Simons Foundation (326844)
- Boris I Shraiman
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Katia Koelle, Emory University, United States
Version history
- Received: December 7, 2018
- Accepted: September 14, 2019
- Accepted Manuscript published: September 18, 2019 (version 1)
- Version of Record published: October 23, 2019 (version 2)
Copyright
© 2019, Yan et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,512
- views
-
- 388
- downloads
-
- 30
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Computational and Systems Biology
- Physics of Living Systems
Locomotion is a fundamental behavior of Caenorhabditis elegans (C. elegans). Previous works on kinetic simulations of animals helped researchers understand the physical mechanisms of locomotion and the muscle-controlling principles of neuronal circuits as an actuator part. It has yet to be understood how C. elegans utilizes the frictional forces caused by the tension of its muscles to perform sequenced locomotive behaviors. Here, we present a two-dimensional rigid body chain model for the locomotion of C. elegans by developing Newtonian equations of motion for each body segment of C. elegans. Having accounted for friction-coefficients of the surrounding environment, elastic constants of C. elegans, and its kymogram from experiments, our kinetic model (ElegansBot) reproduced various locomotion of C. elegans such as, but not limited to, forward-backward-(omega turn)-forward locomotion constituting escaping behavior and delta-turn navigation. Additionally, ElegansBot precisely quantified the forces acting on each body segment of C. elegans to allow investigation of the force distribution. This model will facilitate our understanding of the detailed mechanism of various locomotive behaviors at any given friction-coefficients of the surrounding environment. Furthermore, as the model ensures the performance of realistic behavior, it can be used to research actuator-controller interaction between muscles and neuronal circuits.
-
- Physics of Living Systems
Termites build complex nests which are an impressive example of self-organization. We know that the coordinated actions involved in the construction of these nests by multiple individuals are primarily mediated by signals and cues embedded in the structure of the nest itself. However, to date there is still no scientific consensus about the nature of the stimuli that guide termite construction, and how they are sensed by termites. In order to address these questions, we studied the early building behavior of Coptotermes gestroi termites in artificial arenas, decorated with topographic cues to stimulate construction. Pellet collections were evenly distributed across the experimental setup, compatible with a collection mechanism that is not affected by local topography, but only by the distribution of termite occupancy (termites pick pellets at the positions where they are). Conversely, pellet depositions were concentrated at locations of high surface curvature and at the boundaries between different types of substrate. The single feature shared by all pellet deposition regions was that they correspond to local maxima in the evaporation flux. We can show analytically and we confirm experimentally that evaporation flux is directly proportional to the local curvature of nest surfaces. Taken together, our results indicate that surface curvature is sufficient to organize termite building activity and that termites likely sense curvature indirectly through substrate evaporation. Our findings reconcile the apparently discordant results of previous studies.