Reconstruction of Microbial Haplotypes by Integration of Statistical and Physical Linkage in Scaffolding

DNA sequencing technologies provide unprecedented opportunities to analyze within-host evolution of microorganism populations. Often, within-host populations are analyzed via pooled sequencing of the population, which contains multiple individuals or “haplotypes.” However, current next-generation sequencing instruments, in conjunction with single-molecule barcoded linked-reads, cannot distinguish long haplotypes directly. Computational reconstruction of haplotypes from pooled sequencing has been attempted in virology, bacterial genomics, metagenomics, and human genetics, using algorithms based on either cross-host genetic sharing or within-host genomic reads. Here, we describe PoolHapX, a flexible computational approach that integrates information from both genetic sharing and genomic sequencing. We demonstrated that PoolHapX outperforms state-of-the-art tools tailored to specific organismal systems, and is robust to within-host evolution. Importantly, together with barcoded linked-reads, PoolHapX can infer whole-chromosome-scale haplotypes from 50 pools each containing 12 different haplotypes. By analyzing real data, we uncovered dynamic variations in the evolutionary processes of within-patient HIV populations previously unobserved in single position-based analysis.     

From within-host interactions to epidemiological competition: a general model for multiple infections

Many hosts are infected by several parasite genotypes at a time. In these co-infected hosts, parasites can interact in various ways thus creating diverse within-host dynamics, making it difficult to predict the expression and the evolution of virulence. Moreover, multiple infections generate a combinatorial diversity of cotransmission routes at the host population level, which complicates the epidemiology and may lead to non-trivial outcomes. We introduce a new model for multiple infections, which allows any number of parasite genotypes to infect hosts and potentially coexist in the population. In our model, parasites affect one another''s within-host growth through density-dependent interactions and by means of public goods and spite. These within-host interactions determine virulence, recovery and transmission rates, which are then integrated in a transmission network. We use analytical solutions and numerical simulations to investigate epidemiological feedbacks in host populations infected by several parasite genotypes. Finally, we discuss general perspectives on multiple infections.     

Antigenic Diversity,Transmission Mechanisms,and the Evolution of Pathogens

Pathogens have evolved diverse strategies to maximize their transmission fitness. Here we investigate these strategies for directly transmitted pathogens using mathematical models of disease pathogenesis and transmission, modeling fitness as a function of within- and between-host pathogen dynamics. The within-host model includes realistic constraints on pathogen replication via resource depletion and cross-immunity between pathogen strains. We find three distinct types of infection emerge as maxima in the fitness landscape, each characterized by particular within-host dynamics, host population contact network structure, and transmission mode. These three infection types are associated with distinct non-overlapping ranges of levels of antigenic diversity, and well-defined patterns of within-host dynamics and between-host transmissibility. Fitness, quantified by the basic reproduction number, also falls within distinct ranges for each infection type. Every type is optimal for certain contact structures over a range of contact rates. Sexually transmitted infections and childhood diseases are identified as exemplar types for low and high contact rates, respectively. This work generates a plausible mechanistic hypothesis for the observed tradeoff between pathogen transmissibility and antigenic diversity, and shows how different classes of pathogens arise evolutionarily as fitness optima for different contact network structures and host contact rates.     

Microevolutionary analysis of Clostridium difficile genomes to investigate transmission

Background

The control of Clostridium difficile infection is a major international healthcare priority, hindered by a limited understanding of transmission epidemiology for these bacteria. However, transmission studies of bacterial pathogens are rapidly being transformed by the advent of next generation sequencing.

Results

Here we sequence whole C. difficile genomes from 486 cases arising over four years in Oxfordshire. We show that we can estimate the times back to common ancestors of bacterial lineages with sufficient resolution to distinguish whether direct transmission is plausible or not. Time depths were inferred using a within-host evolutionary rate that we estimated at 1.4 mutations per genome per year based on serially isolated genomes. The subset of plausible transmissions was found to be highly associated with pairs of patients sharing time and space in hospital. Conversely, the large majority of pairs of genomes matched by conventional typing and isolated from patients within a month of each other were too distantly related to be direct transmissions.

Conclusions

Our results confirm that nosocomial transmission between symptomatic C. difficile cases contributes far less to current rates of infection than has been widely assumed, which clarifies the importance of future research into other transmission routes, such as from asymptomatic carriers. With the costs of DNA sequencing rapidly falling and its use becoming more and more widespread, genomics will revolutionize our understanding of the transmission of bacterial pathogens.     

An out-of-body experience: the extracellular dimension for the transmission of mutualistic bacteria in insects

Across animals and plants, numerous metabolic and defensive adaptations are a direct consequence of symbiotic associations with beneficial microbes. Explaining how these partnerships are maintained through evolutionary time remains one of the central challenges within the field of symbiosis research. While genome erosion and co-cladogenesis with the host are well-established features of symbionts exhibiting intracellular localization and transmission, the ecological and evolutionary consequences of an extracellular lifestyle have received little attention, despite a demonstrated prevalence and functional importance across many host taxa. Using insect–bacteria symbioses as a model, we highlight the diverse routes of extracellular symbiont transfer. Extracellular transmission routes are unified by the common ability of the bacterial partners to survive outside their hosts, thereby imposing different genomic, metabolic and morphological constraints than would be expected from a strictly intracellular lifestyle. We emphasize that the evolutionary implications of symbiont transmission routes (intracellular versus extracellular) do not necessarily correspond to those of the transmission mode (vertical versus horizontal), a distinction of vital significance when addressing the genomic and physiological consequences for both host and symbiont.     

'SEEDY' (Simulation of Evolutionary and Epidemiological Dynamics): An R Package to Follow Accumulation of Within-Host Mutation in Pathogens

Genome sequencing is an increasingly common component of infectious disease outbreak investigations. However, the relationship between pathogen transmission and observed genetic data is complex, and dependent on several uncertain factors. As such, simulation of pathogen dynamics is an important tool for interpreting observed genomic data in an infectious disease outbreak setting, in order to test hypotheses and to explore the range of outcomes consistent with a given set of parameters. We introduce ‘seedy’, an R package for the simulation of evolutionary and epidemiological dynamics (http://cran.r-project.org/web/packages/seedy/). Our software implements stochastic models for the accumulation of mutations within hosts, as well as individual-level disease transmission. By allowing variables such as the transmission bottleneck size, within-host effective population size and population mixing rates to be specified by the user, our package offers a flexible framework to investigate evolutionary dynamics during disease outbreaks. Furthermore, our software provides theoretical pairwise genetic distance distributions to provide a likelihood of person-to-person transmission based on genomic observations, and using this framework, implements transmission route assessment for genomic data collected during an outbreak. Our open source software provides an accessible platform for users to explore pathogen evolution and outbreak dynamics via simulation, and offers tools to assess observed genomic data in this context.     

Crossing the scale from within-host infection dynamics to between-host transmission fitness: a discussion of current assumptions and knowledge

The progression of an infection within a host determines the ability of a pathogen to transmit to new hosts and to maintain itself in the population. While the general connection between the infection dynamics within a host and the population-level transmission dynamics of pathogens is widely acknowledged, a comprehensive and quantitative understanding that would allow full integration of the two scales is still lacking. Here, we provide a brief discussion of both models and data that have attempted to provide quantitative mappings from within-host infection dynamics to transmission fitness. We present a conceptual framework and provide examples of studies that have taken first steps towards development of a quantitative framework that scales from within-host infections to population-level fitness of different pathogens. We hope to illustrate some general themes, summarize some of the recent advances and—maybe most importantly—discuss gaps in our ability to bridge these scales, and to stimulate future research on this important topic.     

The population genomics of within-host Mycobacterium tuberculosis

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.Subject terms: Genetic variation, Population genetics     

Transmission of Equine Influenza Virus during an Outbreak Is Characterized by Frequent Mixed Infections and Loose Transmission Bottlenecks

The ability of influenza A viruses (IAVs) to cross species barriers and evade host immunity is a major public health concern. Studies on the phylodynamics of IAVs across different scales – from the individual to the population – are essential for devising effective measures to predict, prevent or contain influenza emergence. Understanding how IAVs spread and evolve during outbreaks is critical for the management of epidemics. Reconstructing the transmission network during a single outbreak by sampling viral genetic data in time and space can generate insights about these processes. Here, we obtained intra-host viral sequence data from horses infected with equine influenza virus (EIV) to reconstruct the spread of EIV during a large outbreak. To this end, we analyzed within-host viral populations from sequences covering 90% of the infected yards. By combining gene sequence analyses with epidemiological data, we inferred a plausible transmission network, in turn enabling the comparison of transmission patterns during the course of the outbreak and revealing important epidemiological features that were not apparent using either approach alone. The EIV populations displayed high levels of genetic diversity, and in many cases we observed distinct viral populations containing a dominant variant and a number of related minor variants that were transmitted between infectious horses. In addition, we found evidence of frequent mixed infections and loose transmission bottlenecks in these naturally occurring populations. These frequent mixed infections likely influence the size of epidemics.     

An integrated model of Plasmodium falciparum dynamics

The within-host and between-host dynamics of malaria are linked in myriad ways, but most obviously by gametocytes, the parasite blood forms transmissible from human to mosquito. Gametocyte dynamics depend on those of non-transmissible blood forms, which stimulate immune responses, impeding transmission as well as within-host parasite densities. These dynamics can, in turn, influence antigenic diversity and recombination between genetically distinct parasites. Here, we embed a differential-equation model of parasite-immune system interactions within each of the individual humans represented in a discrete-event model of Plasmodium falciparum transmission, and examine the effects of human population turnover, parasite antigenic diversity, recombination, and gametocyte production on the dynamics of malaria. Our results indicate that the local persistence of P. falciparum increases with turnover in the human population and antigenic diversity in the parasite, particularly in combination, and that antigenic diversity arising from meiotic recombination in the parasite has complex differential effects on the persistence of founder and progeny genotypes. We also find that reductions in the duration of individual human infectivity to mosquitoes, even if universal, produce population-level effects only if near-absolute, and that, in competition, the persistence and prevalence of parasite genotypes with gametocyte production concordant with data exceed those of genotypes with higher gametocyte production. This new, integrated approach provides a framework for investigating relationships between pathogen dynamics within an individual host and pathogen dynamics within interacting host and vector populations.     

Secondary Analysis of the NCI-60 Whole Exome Sequencing Data Indicates Significant Presence of Propionibacterium acnes Genomic Material in Leukemia (RPMI-8226) and Central Nervous System (SF-295, SF-539, and SNB-19) Cell Lines

The NCI-60 human tumor cell line panel has been used in a broad range of cancer research over the last two decades. A landmark 2013 whole exome sequencing study of this panel added an exceptional new resource for cancer biologists. The complementary analysis of the sequencing data produced by this study suggests the presence of Propionibacterium acnes genomic sequences in almost half of the datasets, with the highest abundance in the leukemia (RPMI-8226) and central nervous system (SF-295, SF-539, and SNB-19) cell lines. While the origin of these contaminating bacterial sequences remains to be determined, observed results suggest that computational control for the presence of microbial genomic material is a necessary step in the analysis of the high throughput sequencing (HTS) data.     

Temporal dynamics of SARS-CoV-2 mutation accumulation within and across infected hosts

Analysis of SARS-CoV-2 genetic diversity within infected hosts can provide insight into the generation and spread of new viral variants and may enable high resolution inference of transmission chains. However, little is known about temporal aspects of SARS-CoV-2 intrahost diversity and the extent to which shared diversity reflects convergent evolution as opposed to transmission linkage. Here we use high depth of coverage sequencing to identify within-host genetic variants in 325 specimens from hospitalized COVID-19 patients and infected employees at a single medical center. We validated our variant calling by sequencing defined RNA mixtures and identified viral load as a critical factor in variant identification. By leveraging clinical metadata, we found that intrahost diversity is low and does not vary by time from symptom onset. This suggests that variants will only rarely rise to appreciable frequency prior to transmission. Although there was generally little shared variation across the sequenced cohort, we identified intrahost variants shared across individuals who were unlikely to be related by transmission. These variants did not precede a rise in frequency in global consensus genomes, suggesting that intrahost variants may have limited utility for predicting future lineages. These results provide important context for sequence-based inference in SARS-CoV-2 evolution and epidemiology.     

Targeted sequence capture as a powerful tool for evolutionary analysis1

Next-generation sequencing technologies (NGS) have revolutionized biological research by significantly increasing data generation while simultaneously decreasing the time to data output. For many ecologists and evolutionary biologists, the research opportunities afforded by NGS are substantial; even for taxa lacking genomic resources, large-scale genome-level questions can now be addressed, opening up many new avenues of research. While rapid and massive sequencing afforded by NGS increases the scope and scale of many research objectives, whole genome sequencing is often unwarranted and unnecessarily complex for specific research questions. Recently developed targeted sequence enrichment, coupled with NGS, represents a beneficial strategy for enhancing data generation to answer questions in ecology and evolutionary biology. This marriage of technologies offers researchers a simple method to isolate and analyze a few to hundreds, or even thousands, of genes or genomic regions from few to many samples in a relatively efficient and effective manner. These strategies can be applied to questions at both the infra- and interspecific levels, including those involving parentage, gene flow, divergence, phylogenetics, reticulate evolution, and many more. Here we provide a brief overview of targeted sequence enrichment, and emphasize the power of this technology to increase our ability to address a wide range of questions of interest to ecologists and evolutionary biologists, particularly for those working with taxa for which few genomic resources are available.     

Deciphering conjugative plasmid permissiveness in wastewater microbiomes

Wastewater treatment plants (WWTPs) are designed to robustly treat polluted water. They are characterized by ceaseless flows of organic, chemical and microbial matter, followed by treatment steps before environmental release. WWTPs are hotspots of horizontal gene transfer between bacteria via conjugative plasmids, leading to dissemination of potentially hazardous genetic material such as antimicrobial resistance genes (AMRGs). While current focus is on the threat of AMRGs spreading and their environmental maintenance, conjugative plasmid transfer dynamics within and between bacterial communities still remains largely uncharted. Furthermore, current in vitro methods used to assess conjugation in complex microbiomes do not include in situ behaviours of recipient cells, resulting in partial understanding of transfers. We investigated the in vitro conjugation capacities of WWTP microbiomes from inlet sewage and outlet treated water using the broad‐host range IncP‐1 conjugative plasmid, pKJK5. A thorough molecular approach coupling metagenomes to 16S rRNA DNA/cDNA amplicon sequencing was established to characterize microbiomes using the ecological concept of functional response groups. A broad diversity of recipient bacterial phyla for the plasmid was observed, especially in WWTP outlets. We also identified permissive bacteria potentially able to cross WWTPs and engage in conjugation before and after water treatment. Bacterial activity and lifestyle seem to influence conjugation extent, as treated water copiotrophs were the most represented strategist amongst transconjugants. Correlation analysis highlighted possible plasmid transmission routes into communities between the sewage to the environment, with identification of keystone members (e.g., Arcobacter) potentially involved in cross‐border exchanges between distant Gram‐positive and Gram‐negative phyla.     

The Genealogical Population Dynamics of HIV-1 in a Large Transmission Chain: Bridging within and among Host Evolutionary Rates

Transmission lies at the interface of human immunodeficiency virus type 1 (HIV-1) evolution within and among hosts and separates distinct selective pressures that impose differences in both the mode of diversification and the tempo of evolution. In the absence of comprehensive direct comparative analyses of the evolutionary processes at different biological scales, our understanding of how fast within-host HIV-1 evolutionary rates translate to lower rates at the between host level remains incomplete. Here, we address this by analyzing pol and env data from a large HIV-1 subtype C transmission chain for which both the timing and the direction is known for most transmission events. To this purpose, we develop a new transmission model in a Bayesian genealogical inference framework and demonstrate how to constrain the viral evolutionary history to be compatible with the transmission history while simultaneously inferring the within-host evolutionary and population dynamics. We show that accommodating a transmission bottleneck affords the best fit our data, but the sparse within-host HIV-1 sampling prevents accurate quantification of the concomitant loss in genetic diversity. We draw inference under the transmission model to estimate HIV-1 evolutionary rates among epidemiologically-related patients and demonstrate that they lie in between fast intra-host rates and lower rates among epidemiologically unrelated individuals infected with HIV subtype C. Using a new molecular clock approach, we quantify and find support for a lower evolutionary rate along branches that accommodate a transmission event or branches that represent the entire backbone of transmitted lineages in our transmission history. Finally, we recover the rate differences at the different biological scales for both synonymous and non-synonymous substitution rates, which is only compatible with the ‘store and retrieve’ hypothesis positing that viruses stored early in latently infected cells preferentially transmit or establish new infections upon reactivation.     

Insights from Genomics into Bacterial Pathogen Populations

Bacterial pathogens impose a heavy burden of disease on human populations worldwide. The gravest threats are posed by highly virulent respiratory pathogens, enteric pathogens, and HIV-associated infections. Tuberculosis alone is responsible for the deaths of 1.5 million people annually. Treatment options for bacterial pathogens are being steadily eroded by the evolution and spread of drug resistance. However, population-level whole genome sequencing offers new hope in the fight against pathogenic bacteria. By providing insights into bacterial evolution and disease etiology, these approaches pave the way for novel interventions and therapeutic targets. Sequencing populations of bacteria across the whole genome provides unprecedented resolution to investigate (i) within-host evolution, (ii) transmission history, and (iii) population structure. Moreover, advances in rapid benchtop sequencing herald a new era of real-time genomics in which sequencing and analysis can be deployed within hours in response to rapidly changing public health emergencies. The purpose of this review is to highlight the transformative effect of population genomics on bacteriology, and to consider the prospects for answering abiding questions such as why bacteria cause disease.     

Bridging scales in the evolution of infectious disease life histories: application

Within- and between-host disease processes occur on the same timescales, therefore changes in the within-host dynamics of parasites, resources, and immunity can interact with changes in the epidemiological dynamics to affect evolutionary outcomes. Consequently, studies of the evolution of disease life histories, that is, infection-age-specific patterns of transmission and virulence, have been constrained by the need for a mechanistic understanding of within-host disease dynamics. In a companion paper (Day et al. 2011), we develop a novel approach that quantifies the relevant within-host aspects of disease through genetic covariance functions. Here, we demonstrate how to apply this theory to data. Using two previously published datasets from rodent malaria infections, we show how to translate experimental measures into disease life-history traits, and how to quantify the covariance in these traits. Our results show how patterns of covariance can interact with epidemiological dynamics to affect evolutionary predictions for disease life history. We also find that the selective constraints on disease life-history evolution can vary qualitatively, and that "simple" virulence-transmission trade-offs that are often the subject of experimental investigation can be obscured by trade-offs within one trait alone. Finally, we highlight the type and quality of data required for future applications.     

Backward bifurcations and multiple equilibria in epidemic models with structured immunity

Many disease pathogens stimulate immunity in their hosts, which then wanes over time. To better understand the impact of this immunity on epidemiological dynamics, we propose an epidemic model structured according to immunity level that can be applied in many different settings. Under biologically realistic hypotheses, we find that immunity alone never creates a backward bifurcation of the disease-free steady state. This does not rule out the possibility of multiple stable equilibria, but we provide two sufficient conditions for the uniqueness of the endemic equilibrium, and show that these conditions ensure uniqueness in several common special cases. Our results indicate that the within-host dynamics of immunity can, in principle, have important consequences for population-level dynamics, but also suggest that this would require strong non-monotone effects in the immune response to infection. Neutralizing antibody titer data for measles are used to demonstrate the biological application of our theory.