Five challenges in evolution and infectious diseases

Mutation Effects on Structure and Dynamics: Adaptive Evolution of the SARS-CoV-2 Main Protease

The main protease of SARS-CoV-2 (Mpro) plays a critical role in viral replication; although it is relatively conserved, Mpro has nevertheless evolved over the course of the COVID-19 pandemic. Here, we examine phenotypic changes in clinically observed variants of Mpro, relative to the originally reported wild-type enzyme. Using atomistic molecular dynamics simulations, we examine effects of mutation on protein structure and dynamics. In addition to basic structural properties such as variation in surface area and torsion angles, we use protein structure networks and active site networks to evaluate functionally relevant characters related to global cohesion and active site constraint. Substitution analysis shows a continuing trend toward more hydrophobic residues that are dependent on the location of the residue in primary, secondary, tertiary, and quaternary structures. Phylogenetic analysis provides additional evidence for the impact of selective pressure on mutation of Mpro. Overall, these analyses suggest evolutionary adaptation of Mpro toward more hydrophobicity and a less-constrained active site in response to the selective pressures of a novel host environment.

The mechanism shaping the logistic growth of mutation proportion in epidemics at population scale

Virus evolution is a common process of pathogen adaption to host population and environment. Frequently, a small but important fraction of virus mutations are reported to contribute to higher risks of host infection, which is one of the major determinants of infectious diseases outbreaks at population scale. The key mutations contributing to transmission advantage of a genetic variant often grow and reach fixation rapidly. Based on classic epidemiology theories of disease transmission, we proposed a mechanistic explanation of the process that between-host transmission advantage may shape the observed logistic curve of the mutation proportion in population. The logistic growth of mutation is further generalized by incorporating time-varying selective pressure to account for impacts of external factors on pathogen adaptiveness. The proposed model is implemented in real-world data of COVID-19 to capture the emerging trends and changing dynamics of the B.1.1.7 strains of SARS-CoV-2 in England. The model characterizes and establishes the underlying theoretical mechanism that shapes the logistic growth of mutation in population.

Modelling: Understanding pandemics and how to control them

New disease challenges, societal demands and better or novel types of data, drive innovations in the structure, formulation and analysis of epidemic models. Innovations in modelling can lead to new insights into epidemic processes and better use of available data, yielding improved disease control and stimulating collection of better data and new data types. Here we identify key challenges for the structure, formulation, analysis and use of mathematical models of pathogen transmission relevant to current and future pandemics.

Phylogenetic and Ancestral Sequence Reconstruction of SARS-CoV-2 Reveals Latent Capacity to Bind Human ACE2 Receptor

SARS-CoV-2 is a unique event, having emerged suddenly as a highly infectious viral pathogen for human populations. Previous phylogenetic analyses show its closest known evolutionary relative to be a virus detected in bats (RaTG13), with a common assumption that SARS-CoV-2 evolved from a zoonotic ancestor via recent genetic changes (likely in the Spike protein receptor-binding domain or RBD) that enabled it to infect humans. We used detailed phylogenetic analysis, ancestral sequence reconstruction, and in situ molecular dynamics simulations to examine the Spike-RBD's functional evolution, finding that the common ancestral virus with RaTG13, dating to no later than 2013, possessed high binding affinity to the human ACE2 receptor. This suggests that SARS-CoV-2 likely possessed a latent capacity to bind to human cellular targets (though this may not have been sufficient for successful infection) and emphasizes the importance of expanding efforts to catalog and monitor viruses circulating in both human and non-human populations.

Challenges in modeling the emergence of novel pathogens

The emergence of infectious agents with pandemic potential present scientific challenges from detection to data interpretation to understanding determinants of risk and forecasts. Mathematical models could play an essential role in how we prepare for future emergent pathogens. Here, we describe core directions for expansion of the existing tools and knowledge base, including: using mathematical models to identify critical directions and paths for strengthening data collection to detect and respond to outbreaks of novel pathogens; expanding basic theory to identify infectious agents and contexts that present the greatest risks, over both the short and longer term; by strengthening estimation tools that make the most use of the likely range and uncertainties in existing data; and by ensuring modelling applications are carefully communicated and developed within diverse and equitable collaborations for increased public health benefit.

Genotypic variation in an ecologically important parasite is associated with host species, lake and spore size

Genetic variation in parasites has important consequences for host–parasite interactions. Prior studies of the ecologically important parasite Metschnikowia bicuspidata have suggested low genetic variation in the species. Here, we collected M. bicuspidata from two host species (Daphnia dentifera and Ceriodaphnia dubia) and two regions (Michigan and Indiana, USA). Within a lake, outbreaks tended to occur in one host species but not the other. Using microsatellite markers, we identified six parasite genotypes grouped within three distinct clades, one of which was rare. Of the two main clades, one was generally associated with D. dentifera, with lakes in both regions containing a single genotype. The other M. bicuspidata clade was mainly associated with C. dubia, with a different genotype dominating in each region. Despite these associations, both D. dentifera- and C. dubia-associated genotypes were found infecting both hosts in lakes. However, in lab experiments, the D. dentifera-associated genotype infected both D. dentifera and C. dubia, but the C. dubia-associated genotype, which had spores that were approximately 30% smaller, did not infect D. dentifera. We hypothesize that variation in spore size might help explain patterns of cross-species transmission. Future studies exploring the causes and consequences of variation in spore size may help explain patterns of infection and the maintenance of genotypic diversity in this ecologically important system.

Quantifying the Evolutionary Constraints and Potential of Hepatitis C Virus NS5A Protein

RNA viruses, such as hepatitis C virus (HCV), influenza virus, and SARS-CoV-2, are notorious for their ability to evolve rapidly under selection in novel environments. It is known that the high mutation rate of RNA viruses can generate huge genetic diversity to facilitate viral adaptation. However, less attention has been paid to the underlying fitness landscape that represents the selection forces on viral genomes, especially under different selection conditions. Here, we systematically quantified the distribution of fitness effects of about 1,600 single amino acid substitutions in the drug-targeted region of NS5A protein of HCV. We found that the majority of nonsynonymous substitutions incur large fitness costs, suggesting that NS5A protein is highly optimized. The replication fitness of viruses is correlated with the pattern of sequence conservation in nature, and viral evolution is constrained by the need to maintain protein stability. We characterized the adaptive potential of HCV by subjecting the mutant viruses to selection by the antiviral drug daclatasvir at multiple concentrations. Both the relative fitness values and the number of beneficial mutations were found to increase with the increasing concentrations of daclatasvir. The changes in the spectrum of beneficial mutations in NS5A protein can be explained by a pharmacodynamics model describing viral fitness as a function of drug concentration. Overall, our results show that the distribution of fitness effects of mutations is modulated by both the constraints on the biophysical properties of proteins (i.e., selection pressure for protein stability) and the level of environmental stress (i.e., selection pressure for drug resistance).

IMPORTANCE Many viruses adapt rapidly to novel selection pressures, such as antiviral drugs. Understanding how pathogens evolve under drug selection is critical for the success of antiviral therapy against human pathogens. By combining deep sequencing with selection experiments in cell culture, we have quantified the distribution of fitness effects of mutations in hepatitis C virus (HCV) NS5A protein. Our results indicate that the majority of single amino acid substitutions in NS5A protein incur large fitness costs. Simulation of protein stability suggests viral evolution is constrained by the need to maintain protein stability. By subjecting the mutant viruses to selection under an antiviral drug, we find that the adaptive potential of viral proteins in a novel environment is modulated by the level of environmental stress, which can be explained by a pharmacodynamics model. Our comprehensive characterization of the fitness landscapes of NS5A can potentially guide the design of effective strategies to limit viral evolution.

The Expectations and Challenges of Wildlife Disease Research in the Era of Genomics: Forecasting with a Horizon Scan-like Exercise

The outbreak and transmission of disease-causing pathogens are contributing to the unprecedented rate of biodiversity decline. Recent advances in genomics have coalesced into powerful tools to monitor, detect, and reconstruct the role of pathogens impacting wildlife populations. Wildlife researchers are thus uniquely positioned to merge ecological and evolutionary studies with genomic technologies to exploit unprecedented "Big Data" tools in disease research; however, many researchers lack the training and expertise required to use these computationally intensive methodologies. To address this disparity, the inaugural "Genomics of Disease in Wildlife" workshop assembled early to mid-career professionals with expertise across scientific disciplines (e.g., genomics, wildlife biology, veterinary sciences, and conservation management) for training in the application of genomic tools to wildlife disease research. A horizon scanning-like exercise, an activity to identify forthcoming trends and challenges, performed by the workshop participants identified and discussed 5 themes considered to be the most pressing to the application of genomics in wildlife disease research: 1) "Improving communication," 2) "Methodological and analytical advancements," 3) "Translation into practice," 4) "Integrating landscape ecology and genomics," and 5) "Emerging new questions." Wide-ranging solutions from the horizon scan were international in scope, itemized both deficiencies and strengths in wildlife genomic initiatives, promoted the use of genomic technologies to unite wildlife and human disease research, and advocated best practices for optimal use of genomic tools in wildlife disease projects. The results offer a glimpse of the potential revolution in human and wildlife disease research possible through multi-disciplinary collaborations at local, regional, and global scales.

pMHC Structural Comparisons as a Pivotal Element to Detect and Validate T-Cell Targets for Vaccine Development and Immunotherapy-A New Methodological Proposal

The search for epitopes that will effectively trigger an immune response remains the "El Dorado" for immunologists. The development of promising immunotherapeutic approaches requires the appropriate targets to elicit a proper immune response. Considering the high degree of HLA/TCR diversity, as well as the heterogeneity of viral and tumor proteins, this number will invariably be higher than ideal to test. It is known that the recognition of a peptide-MHC (pMHC) by the T-cell receptor is performed entirely in a structural fashion, where the atomic interactions of both structures, pMHC and TCR, dictate the fate of the process. However, epitopes with a similar composition of amino acids can produce dissimilar surfaces. Conversely, sequences with no conspicuous similarities can exhibit similar TCR interaction surfaces. In the last decade, our group developed a database and in silico structural methods to extract molecular fingerprints that trigger T-cell immune responses, mainly referring to physicochemical similarities, which could explain the immunogenic differences presented by different pMHC-I complexes. Here, we propose an immunoinformatic approach that considers a structural level of information, combined with an experimental technology that simulates the presentation of epitopes for a T cell, to improve vaccine production and immunotherapy efficacy.

A general framework for modelling the impact of co-infections on pathogen evolution

Theoretical models suggest that mixed-strain infections, or co-infections, are an important driver of pathogen evolution. However, the within-host dynamics of co-infections vary enormously, which complicates efforts to develop a general understanding of how co-infections affect evolution. Here, we develop a general framework which condenses the within-host dynamics of co-infections into a few key outcomes, the most important of which is the overall R0 of the co-infection. Similar to how fitness is determined by two different alleles in a heterozygote, the R0 of a co-infection is a product of the R0 values of the co-infecting strains, shaped by the interaction of those strains at the within-host level. Extending the analogy, we propose that the overall R0 reflects the dominance of the co-infecting strains, and that the ability of a mutant strain to invade a population is a function of its dominance in co-infections. To illustrate the utility of these concepts, we use a within-host model to show how dominance arises from the within-host dynamics of a co-infection, and then use an epidemiological model to demonstrate that dominance is a robust predictor of the ability of a mutant strain to save a maladapted wild-type strain from extinction (evolutionary emergence).

Bridging the gap between HIV epidemiology and antiretroviral resistance evolution: Modelling the spread of resistance in South Africa

The scale-up of antiretroviral therapy (ART) in South Africa substantially reduced AIDS-related deaths and new HIV infections. However, its success is threatened by the emergence of resistance to non-nucleoside reverse-transcriptase inhibitors (NNRTI). The MARISA (Modelling Antiretroviral drug Resistance In South Africa) model presented here aims at investigating the time trends and factors driving NNRTI resistance in South Africa. MARISA is a compartmental model that includes the key aspects of the local HIV epidemic: continuum of care, disease progression, and gender. The dynamics of NNRTI resistance emergence and transmission are then added to this framework. Model parameters are informed using data from HIV cohorts participating in the International epidemiology Databases to Evaluate AIDS (IeDEA) and literature estimates, or fitted to UNAIDS estimates. Using this novel approach of triangulating clinical and resistance data from various sources, MARISA reproduces the time trends of HIV in South Africa in 2005–2016, with a decrease in new infections, undiagnosed individuals, and AIDS-related deaths. MARISA captures the dynamics of the spread of NNRTI resistance: high levels of acquired drug resistance (ADR, in 83% of first-line treatment failures in 2016), and increasing transmitted drug resistance (TDR, in 8.1% of ART initiators in 2016). Simulation of counter-factual scenarios reflecting alternative public health policies shows that increasing treatment coverage would have resulted in fewer new infections and deaths, at the cost of higher TDR (11.6% in 2016 for doubling the treatment rate). Conversely, improving switching to second-line treatment would have led to lower TDR (6.5% in 2016 for doubling the switching rate) and fewer new infections and deaths. Implementing drug resistance testing would have had little impact. The rapid ART scale-up and inadequate switching to second-line treatment were the key drivers of the spread of NNRTI resistance in South Africa. However, even though some interventions could have substantially reduced the level of NNRTI resistance, no policy including NNRTI-based first line regimens could have prevented this spread. Thus, by combining epidemiological data on HIV in South Africa with biological data on resistance evolution, our modelling approach identified key factors driving NNRTI resistance, highlighting the need of alternative first-line regimens.

Resistance to non-nucleoside reverse transcriptase inhibitors (NNRTI) threatens the long-term success of antiretroviral therapy (ART) roll-out in South Africa. We developed a compartmental model integrating the local HIV epidemiology with biological mechanisms of drug resistance. A first dimension of the model accounts for the continuum of care: infection, diagnosis, first-line treatment with suppression or failure, and second-line treatment. Other dimensions include: disease progression (CD4 counts), gender, and acquisition and transmission of NNRTI resistance. Whenever possible, we informed the parameters using the data available from local cohorts. Other parameters were informed using literature or UNAIDS estimates. The model captured the rise of NNRTI resistance during the period. We assessed the impact of counter-factual scenarios reflecting alternative countrywide policies during the period 2005 to 2016, considering either increasing ART coverage, improving management of treatment failure, broadening ART eligibility, or implementing drug resistance testing before ART initiation. We identified key drivers of the NNRTI resistance epidemic: large-scale ART roll-out and insufficient monitoring of first-line treatment failure. The model also suggested that no policy including NNRTI-based first line regimens could have prevented the spread of NNRTI resistance.

Digest: Little evidence exists for a virulence-transmission trade-off

Research predicting the impact and spread of infectious disease has been heavily influenced by the idea of an evolutionary trade-off between a pathogen's virulence and its transmission rate. In a meta-analysis of the key underlying relationships, Acevedo et al. (2019) highlight the surprising lack of empirical evidence for this influential hypothesis.

Fogarty International Center collaborative networks in infectious disease modeling: Lessons learnt in research and capacity building

Due to a combination of ecological, political, and demographic factors, the emergence of novel pathogens has been increasingly observed in animals and humans in recent decades. Enhancing global capacity to study and interpret infectious disease surveillance data, and to develop data-driven computational models to guide policy, represents one of the most cost-effective, and yet overlooked, ways to prepare for the next pandemic. Epidemiological and behavioral data from recent pandemics and historic scourges have provided rich opportunities for validation of computational models, while new sequencing technologies and the 'big data' revolution present new tools for studying the epidemiology of outbreaks in real time. For the past two decades, the Division of International Epidemiology and Population Studies (DIEPS) of the NIH Fogarty International Center has spearheaded two synergistic programs to better understand and devise control strategies for global infectious disease threats. The Multinational Influenza Seasonal Mortality Study (MISMS) has strengthened global capacity to study the epidemiology and evolutionary dynamics of influenza viruses in 80 countries by organizing international research activities and training workshops. The Research and Policy in Infectious Disease Dynamics (RAPIDD) program and its precursor activities has established a network of global experts in infectious disease modeling operating at the research-policy interface, with collaborators in 78 countries. These activities have provided evidence-based recommendations for disease control, including during large-scale outbreaks of pandemic influenza, Ebola and Zika virus. Together, these programs have coordinated international collaborative networks to advance the study of emerging disease threats and the field of computational epidemic modeling. A global community of researchers and policy-makers have used the tools and trainings developed by these programs to interpret infectious disease patterns in their countries, understand modeling concepts, and inform control policies. Here we reflect on the scientific achievements and lessons learnt from these programs (h-index = 106 for RAPIDD and 79 for MISMS), including the identification of outstanding researchers and fellows; funding flexibility for timely research workshops and working groups (particularly relative to more traditional investigator-based grant programs); emphasis on group activities such as large-scale modeling reviews, model comparisons, forecasting challenges and special journal issues; strong quality control with a light touch on outputs; and prominence of training, data-sharing, and joint publications.

Coxiella burnetii in wild mammals: A systematic review

Coxiella burnetii is a multi-host bacterium that causes Q fever in humans, a zoonosis that is emerging worldwide. The ecology of C. burnetii in wildlife is still poorly understood and the influence of host, environmental and pathogen factors is almost unknown. This study gathers current published information on different aspects of C. burnetii infection in wildlife, even in species with high reservoir potential and a high rate of interaction with livestock and humans, in order to partially fill the existing gap and highlight future needs. Exposure and/or infection by C. burnetii has, to date, been reported in 109 wild mammal species. The limited sample size of most of the existing studies could suggest an undervalued prevalence of C. burnetii infection. Knowledge on the clinical outcome of C. burnetii infection in wildlife is also very limited, but currently includes reproductive failure in waterbuck (Kobus ellipsiprymnus), roan antelope (Hippotragus niger), dama gazelle (Nanger dama) and water buffalo (Bubalus bubalis) and placentitis in the Pacific harbor seal (Phoca vitulina richardsi), Steller sea lion (Eumetopias jubatus) and red deer (Cervus elaphus). The currently available serological tests need to be optimised and validated for each wildlife species. Finally, there is a huge gap in the research on C. burnetii control in wildlife, despite of the increasing evidence that wildlife is a source of C. burnetii for both livestock and humans.

Opportunities and challenges for modelling epidemiological and evolutionary dynamics in a multihost, multiparasite system: Zoonotic hybrid schistosomiasis in West Africa

Multihost multiparasite systems are evolutionarily and ecologically dynamic, which presents substantial trans-disciplinary challenges for elucidating their epidemiology and designing appropriate control. Evidence for hybridizations and introgressions between parasite species is gathering, in part in line with improvements in molecular diagnostics and genome sequencing. One major system where this is becoming apparent is within the Genus Schistosoma, where schistosomiasis represents a disease of considerable medical and veterinary importance, the greatest burden of which occurs in sub-Saharan Africa. Interspecific hybridizations and introgressions bring an increased level of complexity over and above that already inherent within multihost, multiparasite systems, also representing an additional source of genetic variation that can drive evolution. This has the potential for profound implications for the control of parasitic diseases, including, but not exclusive to, widening host range, increased transmission potential and altered responses to drug therapy. Here, we present the challenging case example of haematobium group Schistosoma spp. hybrids in West Africa, a system involving multiple interacting parasites and multiple definitive hosts, in a region where zoonotic reservoirs of schistosomiasis were not previously considered to be of importance. We consider how existing mathematical model frameworks for schistosome transmission could be expanded and adapted to zoonotic hybrid systems, exploring how such model frameworks can utilize molecular and epidemiological data, as well as the complexities and challenges this presents. We also highlight the opportunities and value such mathematical models could bring to this and a range of similar multihost, multi and cross-hybridizing parasites systems in our changing world.

A complete categorization of multiscale models of infectious disease systems

Modelling of infectious disease systems has entered a new era in which disease modellers are increasingly turning to multiscale modelling to extend traditional modelling frameworks into new application areas and to achieve higher levels of detail and accuracy in characterizing infectious disease systems. In this paper we present a categorization framework for categorizing multiscale models of infectious disease systems. The categorization framework consists of five integration frameworks and five criteria. We use the categorization framework to give a complete categorization of host-level immuno-epidemiological models (HL-IEMs). This categorization framework is also shown to be applicable in categorizing other types of multiscale models of infectious diseases beyond HL-IEMs through modifying the initial categorization framework presented in this study. Categorization of multiscale models of infectious disease systems in this way is useful in bringing some order to the discussion on the structure of these multiscale models.

Mutation and recombination in pathogen evolution: Relevance, methods and controversies

Mutation and recombination drive the evolution of most pathogens by generating the genetic variants upon which selection operates. Those variants can, for example, confer resistance to host immune systems and drug therapies or lead to epidemic outbreaks. Given their importance, diverse evolutionary studies have investigated the abundance and consequences of mutation and recombination in pathogen populations. However, some controversies persist regarding the contribution of each evolutionary force to the development of particular phenotypic observations (e.g., drug resistance). In this study, we revise the importance of mutation and recombination in the evolution of pathogens at both intra-host and inter-host levels. We also describe state-of-the-art analytical methodologies to detect and quantify these two evolutionary forces, including biases that are often ignored in evolutionary studies. Finally, we present some of our former studies involving pathogenic taxa where mutation and recombination played crucial roles in the recovery of pathogenic fitness, the generation of interspecific genetic diversity, or the design of centralized vaccines. This review also illustrates several common controversies and pitfalls in the analysis and in the evaluation and interpretation of mutation and recombination outcomes.

The within-host population dynamics of Mycobacterium tuberculosis vary with treatment efficacy

Background

Combination therapy is one of the most effective tools for limiting the emergence of drug resistance in pathogens. Despite the widespread adoption of combination therapy across diseases, drug resistance rates continue to rise, leading to failing treatment regimens. The mechanisms underlying treatment failure are well studied, but the processes governing successful combination therapy are poorly understood. We address this question by studying the population dynamics of Mycobacterium tuberculosis within tuberculosis patients undergoing treatment with different combinations of antibiotics.

Results

By combining very deep whole genome sequencing (~1000-fold genome-wide coverage) with sequential sputum sampling, we were able to detect transient genetic diversity driven by the apparently continuous turnover of minor alleles, which could serve as the source of drug-resistant bacteria. However, we report that treatment efficacy has a clear impact on the population dynamics: sufficient drug pressure bears a clear signature of purifying selection leading to apparent genetic stability. In contrast, M. tuberculosis populations subject to less drug pressure show markedly different dynamics, including cases of acquisition of additional drug resistance.

Conclusions

Our findings show that for a pathogen like M. tuberculosis, which is well adapted to the human host, purifying selection constrains the evolutionary trajectory to resistance in effectively treated individuals. Nonetheless, we also report a continuous turnover of minor variants, which could give rise to the emergence of drug resistance in cases of drug pressure weakening. Monitoring bacterial population dynamics could therefore provide an informative metric for assessing the efficacy of novel drug combinations.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1196-0) contains supplementary material, which is available to authorized users.

Genomic Analysis of Viral Outbreaks

Genomic analysis is a powerful tool for understanding viral disease outbreaks. Sequencing of viral samples is now easier and cheaper than ever before and can supplement epidemiological methods by providing nucleotide-level resolution of outbreak-causing pathogens. In this review, we describe methods used to answer crucial questions about outbreaks, such as how they began and how a disease is transmitted. More specifically, we explain current techniques for viral sequencing, phylogenetic analysis, transmission reconstruction, and evolutionary investigation of viral pathogens. By detailing the ways in which genomic data can help us understand viral disease outbreaks, we aim to provide a resource that will facilitate the response to future outbreaks.

Antagonistic within-host interactions between plant viruses: molecular basis and impact on viral and host fitness

Double infections of related or unrelated viruses frequently occur in single plants, the viral agents being inoculated into the host plant simultaneously (co-infection) or sequentially (super-infection). Plants attacked by viruses activate sophisticated defence pathways which operate at different levels, often at significant fitness costs, resulting in yield reduction in crop plants. The occurrence and severity of the negative effects depend on the type of within-host interaction between the infecting viruses. Unrelated viruses generally interact with each other in a synergistic manner, whereas interactions between related viruses are mostly antagonistic. These can incur substantial fitness costs to one or both of the competitors. A relatively well-known antagonistic interaction is cross-protection, also referred to as super-infection exclusion. This type of interaction occurs when a previous infection with one virus prevents or interferes with subsequent infection by a homologous second virus. The current knowledge on why and how one virus variant excludes or restricts another is scant. Super-infection exclusion between viruses has predominantly been attributed to the induction of RNA silencing, which is a major antiviral defence mechanism in plants. There are, however, presumptions that various mechanisms are involved in this phenomenon. This review outlines the current state of knowledge concerning the molecular mechanisms behind antagonistic interactions between plant viruses. Harmful or beneficial effects of these interactions on viral and host plant fitness are also characterized. Moreover, the review briefly outlines the past and present attempts to utilize antagonistic interactions among viruses to protect crop plants against destructive diseases.

Opportunities and challenges of Integral Projection Models for modelling host-parasite dynamics

Epidemiological dynamics are shaped by and may in turn shape host demography. These feedbacks can result in hard to predict patterns of disease incidence. Mathematical models that integrate infection and demography are consequently a key tool for informing expectations for disease burden and identifying effective measures for control.

A major challenge is capturing the details of infection within individuals and quantifying their downstream impacts to understand population‐scale outcomes. For example, parasite loads and antibody titres may vary over the course of an infection and contribute to differences in transmission at the scale of the population. To date, to capture these subtleties, models have mostly relied on complex mechanistic frameworks, discrete categorization and/or agent‐based approaches.

Integral Projection Models (IPMs) allow variance in individual trajectories of quantitative traits and their population‐level outcomes to be captured in ways that directly reflect statistical models of trait–fate relationships. Given increasing data availability, and advances in modelling, there is considerable potential for extending this framework to traits of relevance for infectious disease dynamics.

Here, we provide an overview of host and parasite natural history contexts where IPMs could strengthen inference of population dynamics, with examples of host species ranging from mice to sheep to humans, and parasites ranging from viruses to worms. We discuss models of both parasite and host traits, provide two case studies and conclude by reviewing potential for both ecological and evolutionary research.

Siderophore production and the evolution of investment in a public good: An adaptive dynamics approach to kin selection

Like many other bacteria, Pseudomonas aeruginosa sequesters iron from the environment through the secretion, and subsequent uptake, of iron-binding molecules. As these molecules can be taken up by other bacteria in the population than those who secreted them, this is a form of cooperation through a public good. Traditionally, this problem has been studied by comparing the relative fitnesses of siderophore-producing and non-producing strains, but this gives no information about the fate of strains that do produce intermediate amounts of siderophores. Here, we investigate theoretically how the amount invested in this form of cooperation evolves. We use a mechanistic description of the laboratory protocols used in experimental evolution studies to describe the competition and cooperation of the bacteria. From this dynamical model we derive the fitness following the adaptive dynamics method. The results show how selection is driven by local siderophore production and local competition. Because siderophore production reduces the growth rate, local competition decreases with the degree of relatedness (which is a dynamical variable in our model). Our model is not restricted to the analysis of small phenotypic differences and allows for theoretical exploration of the effects of large phenotypic differences between cooperators and cheats. We predict that an intermediate ESS level of cooperation (molecule production) should exist. The adaptive dynamics approach allows us to assess evolutionary stability, which is often not possible in other kin-selection models. We found that selection can lead to an intermediate strategy which in our model is always evolutionarily stable, yet can allow invasion of strategies that are much more cooperative. Our model describes the evolution of a public good in the context of the ecology of the microorganism, which allows us to relate the extent of production of the public good to the details of the interactions.

The adaptive evolution of virulence: a review of theoretical predictions and empirical tests

Why is it that some parasites cause high levels of host damage (i.e. virulence) whereas others are relatively benign? There are now numerous reviews of virulence evolution in the literature but it is nevertheless still difficult to find a comprehensive treatment of the theory and data on the subject that is easily accessible to non-specialists. Here we attempt to do so by distilling the vast theoretical literature on the topic into a set of relatively few robust predictions. We then provide a comprehensive assessment of the available empirical literature that tests these predictions. Our results show that there have been some notable successes in integrating theory and data but also that theory and empiricism in this field do not ‘speak’ to each other very well. We offer a few suggestions for how the connection between the two might be improved.

Modelling challenges in context: lessons from malaria, HIV, and tuberculosis

Malaria, HIV, and tuberculosis (TB) collectively account for several million deaths each year, with all three ranking among the top ten killers in low-income countries. Despite being caused by very different organisms, malaria, HIV, and TB present a suite of challenges for mathematical modellers that are particularly pronounced in these infections, but represent general problems in infectious disease modelling, and highlight many of the challenges described throughout this issue. Here, we describe some of the unifying challenges that arise in modelling malaria, HIV, and TB, including variation in dynamics within the host, diversity in the pathogen, and heterogeneity in human contact networks and behaviour. Through the lens of these three pathogens, we provide specific examples of the other challenges in this issue and discuss their implications for informing public health efforts.

Seven challenges in modeling vaccine preventable diseases

Vaccination has been one of the most successful public health measures since the introduction of basic sanitation. Substantial mortality and morbidity reductions have been achieved via vaccination against many infections, and the list of diseases that are potentially controllable by vaccines is growing steadily. We introduce key challenges for modeling in shaping our understanding and guiding policy decisions related to vaccine preventable diseases.

Eight challenges in phylodynamic inference

The field of phylodynamics, which attempts to enhance our understanding of infectious disease dynamics using pathogen phylogenies, has made great strides in the past decade. Basic epidemiological and evolutionary models are now well characterized with inferential frameworks in place. However, significant challenges remain in extending phylodynamic inference to more complex systems. These challenges include accounting for evolutionary complexities such as changing mutation rates, selection, reassortment, and recombination, as well as epidemiological complexities such as stochastic population dynamics, host population structure, and different patterns at the within-host and between-host scales. An additional challenge exists in making efficient inferences from an ever increasing corpus of sequence data.

Seven challenges for modelling indirect transmission: vector-borne diseases, macroparasites and neglected tropical diseases

Many of the challenges which face modellers of directly transmitted pathogens also arise when modelling the epidemiology of pathogens with indirect transmission – whether through environmental stages, vectors, intermediate hosts or multiple hosts. In particular, understanding the roles of different hosts, how to measure contact and infection patterns, heterogeneities in contact rates, and the dynamics close to elimination are all relevant challenges, regardless of the mode of transmission. However, there remain a number of challenges that are specific and unique to modelling vector-borne diseases and macroparasites. Moreover, many of the neglected tropical diseases which are currently targeted for control and elimination are vector-borne, macroparasitic, or both, and so this article includes challenges which will assist in accelerating the control of these high-burden diseases. Here, we discuss the challenges of indirect measures of infection in humans, whether through vectors or transmission life stages and in estimating the contribution of different host groups to transmission. We also discuss the issues of “evolution-proof” interventions against vector-borne disease.