Emergence and Spread of Basal Lineages of Yersinia pestis during the Neolithic Decline.

Between 5,000 and 6,000 years ago, many Neolithic societies declined throughout western Eurasia due to a combination of factors that are still largely debated. Here, we report the discovery and genome reconstruction of Yersinia pestis, the etiological agent of plague, in Neolithic farmers in Sweden, pre-dating and basal to all modern and ancient known strains of this pathogen. We investigated the history of this strain by combining phylogenetic and molecular clock analyses of the bacterial genome, detailed archaeological information, and genomic analyses from infected individuals and hundreds of ancient human samples across Eurasia. These analyses revealed that multiple and independent lineages of Y. pestis branched and expanded across Eurasia during the Neolithic decline, spreading most likely through early trade networks rather than massive human migrations. Our results are consistent with the existence of a prehistoric plague pandemic that likely contributed to the decay of Neolithic populations in Europe.

Mutation rates and adaptive variation among the clinically dominant clusters of Mycobacterium abscessus.

Mycobacterium abscessus (Mab) is a multidrug-resistant pathogen increasingly responsible for severe pulmonary infections. Analysis of whole-genome sequences (WGS) of Mab demonstrates dense genetic clustering of clinical isolates collected from disparate geographic locations. This has been interpreted as supporting patient-to-patient transmission, but epidemiological studies have contradicted this interpretation. Here, we present evidence for a slowing of the Mab molecular clock rate coincident with the emergence of phylogenetic clusters. We performed phylogenetic inference using publicly available WGS from 483 Mab patient isolates. We implement a subsampling approach in combination with coalescent analysis to estimate the molecular clock rate along the long internal branches of the tree, indicating a faster long-term molecular clock rate compared to branches within phylogenetic clusters. We used ancestry simulation to predict the effects of clock rate variation on phylogenetic clustering and found that the degree of clustering in the observed phylogeny is more easily explained by a clock rate slowdown than by transmission. We also find that phylogenetic clusters are enriched in mutations affecting DNA repair machinery and report that clustered isolates have lower spontaneous mutation rates in vitro. We propose that Mab adaptation to the host environment through variation in DNA repair genes affects the organism's mutation rate and that this manifests as phylogenetic clustering. These results challenge the model that phylogenetic clustering in Mab is explained by person-to-person transmission and inform our understanding of transmission inference in emerging, facultative pathogens.

Gene evolutionary trajectories in Mycobacterium tuberculosis reveal temporal signs of selection.

Genetic differences between different Mycobacterium tuberculosis complex (MTBC) strains determine their ability to transmit within different host populations, their latency times, and their drug resistance profiles. Said differences usually emerge through de novo mutations and are maintained or discarded by the balance of evolutionary forces. Using a dataset of ∼5,000 strains representing global MTBC diversity, we determined the past and present selective forces that have shaped the current variability observed in the pathogen population. We identified regions that have evolved under changing types of selection since the time of the MTBC common ancestor. Our approach highlighted striking differences in the genome regions relevant for host­pathogen interaction and, in particular, suggested an adaptive role for the sensor protein of two-component systems. In addition, we applied our approach to successfully identify potential determinants of resistance to drugs administered as second-line tuberculosis treatments.

Defense Heterogeneity in Host Populations Gives Rise to Pathogen Diversity.

AbstractHost organisms can harbor microbial symbionts that defend them from pathogen infection in addition to the resistance encoded by the host genome. Here, we investigated how variation in defenses, generated from host genetic background and symbiont presence, affects the emergence of pathogen genetic diversity across evolutionary time. We passaged the opportunistic pathogen Pseudomonas aeruginosa through populations of the nematode Caenorhabditis elegans varying in genetic-based defenses and prevalence of a protective symbiont. After 14 passages, we assessed the amount of genetic variation accumulated in evolved pathogen lineages. We found that diversity begets diversity. An overall greater level of pathogen whole-genome and per-gene genetic diversity was measured in pathogens evolved in mixed host populations compared with those evolved in host populations composed of one type of defense. Our findings directly demonstrate that symbiont-generated heterogeneity in host defense can be a significant contributor to pathogen genetic variation.

Sequence capture identifies fastidious chytrid fungi directly from host tissue.

The chytrid fungus Batrachochytrium dendrobatidis (Bd) was discovered in 1998 as the cause of chytridiomycosis, an emerging infectious disease causing mass declines in amphibian populations worldwide. The rapid population declines of the 1970s-1990s were likely caused by the spread of a highly virulent lineage belonging to the Bd-GPL clade that was introduced to naïve susceptible populations. Multiple genetically distinct and regional lineages of Bd have since been isolated and sequenced, greatly expanding the known biological diversity within this fungal pathogen. To date, most Bd research has been restricted to the limited number of samples that could be isolated using culturing techniques, potentially causing a selection bias for strains that can grow on media and missing other unculturable or fastidious strains that are also present on amphibians. We thus attempted to characterize potentially non-culturable genetic lineages of Bd from distinct amphibian taxa using sequence capture technology on DNA extracted from host tissue and swabs. We focused our efforts on host taxa from two different regions that likely harbored distinct Bd clades: (1) wild-caught leopard frogs (Rana) from North America, and (2) a Japanese Giant Salamander (Andrias japonicus) at the Smithsonian Institution's National Zoological Park that exhibited signs of disease and tested positive for Bd using qPCR, but multiple attempts failed to isolate and culture the strain for physiological and genetic characterization. We successfully enriched for and sequenced thousands of fungal genes from both host clades, and Bd load was positively associated with number of recovered Bd sequences. Phylogenetic reconstruction placed all the Rana-derived strains in the Bd-GPL clade. In contrast, the A. japonicus strain fell within the Bd-Asia3 clade, expanding the range of this clade and generating additional genomic data to confirm its placement. The retrieved ITS locus matched public barcoding data from wild A. japonicus and Bd infections found on other amphibians in India and China, suggesting that this uncultured clade is widespread across Asia. Our study underscores the importance of recognizing and characterizing the hidden diversity of fastidious strains in order to reconstruct the spatiotemporal and evolutionary history of Bd. The success of the sequence capture approach highlights the utility of directly sequencing pathogen DNA from host tissue to characterize cryptic diversity that is missed by culture-reliant approaches.

Paleomicrobiology: Diagnosis and Evolution of Ancient Pathogens.

The last century has witnessed progress in the study of ancient infectious disease from purely medical descriptions of past ailments to dynamic interpretations of past population health that draw upon multiple perspectives. The recent adoption of high-throughput DNA sequencing has led to an expanded understanding of pathogen presence, evolution, and ecology across the globe. This genomic revolution has led to the identification of disease-causing microbes in both expected and unexpected contexts, while also providing for the genomic characterization of ancient pathogens previously believed to be unattainable by available methods. In this review we explore the development of DNA-based ancient pathogen research, the specialized methods and tools that have emerged to authenticate and explore infectious disease of the past, and the unique challenges that persist in molecular paleopathology. We offer guidelines to mitigate the impact of these challenges, which will allow for more reliable interpretations of data in this rapidly evolving field of investigation.

Death is overrated: the potential role of detection in driving virulence evolution.

A common assumption in the evolution of virulence theory literature is that pathogens transmit better when they exploit their host more heavily, but by doing so, they impose a greater risk of killing their host, thus truncating infectious periods and reducing their own opportunities for transmission. Here, I derive an equation for the magnitude of this cost in terms of the infection fatality rate, and in doing so, I show that there are many cases where mortality costs are too small to plausibly constrain increases in host exploitation by pathogens. I propose that pathogen evolution may often be constrained by detection costs, whereby hosts alter their behaviour when infection is detectable, and thus reduce pathogen opportunities for onward transmission. I then derive an inequality to illustrate when mortality costs or detection costs impose stronger constraints on pathogen evolution, and I use empirical data from the literature to demonstrate that detection costs are frequently large in both human and animal populations. Finally, I give examples of how evolutionary predictions can change depending on whether costs of host exploitation are borne out through mortality or detection.

Out of Asia? Expansion of Eurasian Lyme borreliosis causing genospecies display unique evolutionary trajectories.

Vector-borne pathogens exist in obligate transmission cycles between vector and reservoir host species. Host and vector shifts can lead to geographic expansion of infectious agents and the emergence of new diseases in susceptible individuals. Three bacterial genospecies (Borrelia afzelii, Borrelia bavariensis, and Borrelia garinii) predominantly utilize two distinct tick species as vectors in Asia (Ixodes persulcatus) and Europe (Ixodes ricinus). Through these vectors, the bacteria can infect various vertebrate groups (e.g., rodents, birds) including humans where they cause Lyme borreliosis, the most common vector-borne disease in the Northern hemisphere. Yet, how and in which order the three Borrelia genospecies colonized each continent remains unclear including the evolutionary consequences of this geographic expansion. Here, by reconstructing the evolutionary history of 142 Eurasian isolates, we found evidence that the ancestors of each of the three genospecies probably have an Asian origin. Even so, each genospecies studied displayed a unique substructuring and evolutionary response to the colonization of Europe. The pattern of allele sharing between continents is consistent with the dispersal rate of the respective vertebrate hosts, supporting the concept that adaptation of Borrelia genospecies to the host is important for pathogen dispersal. Our results highlight that Eurasian Lyme borreliosis agents are all capable of geographic expansion with host association influencing their dispersal; further displaying the importance of host and vector association to the geographic expansion of vector-borne pathogens and potentially conditioning their capacity as emergent pathogens.

Host shift and natural long-distance dispersal to an oceanic island of a host-specific parasite.

Parasite dispersal and host-switching may be better understood by knowing when they occurred. We estimated when the ancestor of a parasite of great reed warblers (Acrocephalus arundinaceus) dispersed to the Seychelles and began infecting the endemic Seychelles warbler (A. sechellensis). We used mitochondrial genomes and published molecular divergence rates to estimate the date of divergence between mitochondrial haplotypes of the parasite Haemoproteus nucleocondensis (lineage GRW01) in the great reed warbler and the Seychelles warbler. We also constructed a time-calibrated phylogeny of the hosts and their relatives to determine when the ancestor of the Seychelles warbler dispersed to the Seychelles. The two GRW01 lineages diverged ca 20-451 kya, long after the ancestor of the Seychelles warbler colonized the Seychelles ca 1.76-4.36 Mya. GRW01 rarely infects other species despite apparent opportunity. Humans were likely not involved in the dispersal of this parasite because humans settled the Seychelles long after the parasite diverged from its mainland relative. Furthermore, introduced birds are unlikely hosts of GRW01. Instead, the ancestor of GRW01 may have dispersed to the Seychelles with an errant migrating great reed warbler. Our results indicate that even specialized parasites can naturally disperse long distances to become emerging infectious diseases.

Cholera Dynamics and the Emergence of Pandemic Vibrio cholerae.

Cholera is a severe diarrheal disease caused by the aquatic bacterium Vibrio cholerae. Interestingly, to date, only one major clade has emerged to cause pandemic disease in humans: the clade that encompasses the strains from the O1 and O139 serogroups. In this chapter, we provide a comprehensive perspective on the virulence factors and mobile genetic elements (MGEs) associated with the emergence of pandemic V. cholerae strains and highlight novel findings such as specific genomic background or interactions between MGEs that explain their confined distribution. Finally, we discuss pandemic cholera dynamics contextualizing them within the evolution of the bacterium.

Honey bee virus causes context-dependent changes in host social behavior.

Anthropogenic changes create evolutionarily novel environments that present opportunities for emerging diseases, potentially changing the balance between host and pathogen. Honey bees provide essential pollination services, but intensification and globalization of honey bee management has coincided with increased pathogen pressure, primarily due to a parasitic mite/virus complex. Here, we investigated how honey bee individual and group phenotypes are altered by a virus of concern, Israeli acute paralysis virus (IAPV). Using automated and manual behavioral monitoring of IAPV-inoculated individuals, we find evidence for pathogen manipulation of worker behavior by IAPV, and reveal that this effect depends on social context; that is, within versus between colony interactions. Experimental inoculation reduced social contacts between honey bee colony members, suggesting an adaptive host social immune response to diminish transmission. Parallel analyses with double-stranded RNA (dsRNA)-immunostimulated bees revealed these behaviors are part of a generalized social immune defensive response. Conversely, inoculated bees presented to groups of bees from other colonies experienced reduced aggression compared with dsRNA-immunostimulated bees, facilitating entry into susceptible colonies. This reduction was associated with a shift in cuticular hydrocarbons, the chemical signatures used by bees to discriminate colony members from intruders. These responses were specific to IAPV infection, suggestive of pathogen manipulation of the host. Emerging bee pathogens may thus shape host phenotypes to increase transmission, a strategy especially well-suited to the unnaturally high colony densities of modern apiculture. These findings demonstrate how anthropogenic changes could affect arms races between human-managed hosts and their pathogens to potentially affect global food security.

Pathogen transmission modes determine contact network structure, altering other pathogen characteristics.

Pathogen traits can vary greatly and heavily impact the ability of a pathogen to persist in a population. Although this variation is fundamental to disease ecology, little is known about the evolutionary pressures that drive these differences, particularly where they interact with host behaviour. We hypothesized that host behaviours relevant to different transmission routes give rise to differences in contact network structure, constraining the space over which pathogen traits can evolve to maximize fitness. Our analysis of 232 contact networks across mammals, birds, reptiles, amphibians, arthropods, fish and molluscs found that contact network topology varies by contact type, most notably in networks that are representative of fluid-exchange transmission. Using infectious disease model simulations, we showed that these differences in network structure suggest pathogens transmitted through fluid-exchange contact types will need traits associated with high transmissibility to successfully proliferate, compared to pathogens that transmit through other types of contact. These findings were supported through a review of known traits of pathogens that transmit in humans. Our work demonstrates that contact network structure may drive the evolution of compensatory pathogen traits according to transmission strategy, providing essential context for understanding pathogen evolution and ecology.

Panzootic chytrid fungus exploits diverse amphibian host environments through plastic infection strategies.

While some pathogens are limited to single species, others can colonize many hosts, likely contributing to the emergence of novel disease outbreaks. Despite this biodiversity threat, traits associated with host niche expansions are not well understood in multihost pathogens. Here, we aimed to uncover functional machinery driving multihost invasion by focusing on Batrachochytrium dendrobatidis (Bd), a pathogen that infects the skin of hundreds of amphibians worldwide. We performed a meta-analysis of Bd gene expression using data from published infection experiments and newly generated profiles. We analysed Bd transcriptomic landscapes across the skin of 14 host species, reconstructed Bd isolates phylogenetic relationships, and inferred the origin and evolutionary history of differentially expressed genes under a phylogenetic framework comprising other 12 zoosporic fungi. Bd displayed plastic infection strategies when challenged by hosts with different disease susceptibility. Our analyses identified sets of differentially expressed genes under host environments with similar infection outcome. We stressed nutritional immunity and gene silencing as important processes required to overcome challenging skin environments in less susceptible hosts. Overall, Bd genes expressed during amphibian skin exploitation have arisen mainly via gene duplications with great family expansions, increasing the gene copy events previously described for this fungal species. Finally, we provide a comprehensive gene data set that can be used to further examine eco-evolutionary hypotheses for this host-pathogen system. Our study supports the idea that host environments exert contrasting selective pressures, such that gene expression plasticity could be one of the evolutionary keys leading to the success of multihost pathogens.

Rapid sequence evolution driven by transposable elements at a virulence locus in a fungal wheat pathogen.

BACKGROUND: Plant pathogens cause substantial crop losses in agriculture production and threaten food security. Plants evolved the ability to recognize virulence factors and pathogens have repeatedly escaped recognition due rapid evolutionary change at pathogen virulence loci (i.e. effector genes). The presence of transposable elements (TEs) in close physical proximity of effector genes can have important consequences for gene regulation and sequence evolution. Species-wide investigations of effector gene loci remain rare hindering our ability to predict pathogen evolvability. RESULTS: Here, we performed genome-wide association studies (GWAS) on a highly polymorphic mapping population of 120 isolates of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe. We identified a major locus underlying significant variation in reproductive success of the pathogen and damage caused on the wheat cultivar Claro. The most strongly associated locus is intergenic and flanked by genes encoding a predicted effector and a serine-type endopeptidase. The center of the locus contained a highly dynamic region consisting of multiple families of TEs. Based on a large global collection of assembled genomes, we show that the virulence locus has undergone substantial recent sequence evolution. Large insertion and deletion events generated length variation between the flanking genes by a factor of seven (5-35 kb). The locus showed also strong signatures of genomic defenses against TEs (i.e. RIP) contributing to the rapid diversification of the locus. CONCLUSIONS: In conjunction, our work highlights the power of combining GWAS and population-scale genome analyses to investigate major effect loci in pathogens.

Comparative Omics Analysis of Historic and Recent Isolates of Bordetella pertussis and Effects of Genome Rearrangements on Evolution.

Despite high vaccination coverage, pertussis is increasing in many industrialized countries, including the Czech Republic. To better understand Bordetella pertussis resurgence, we analyzed historic strains and recent clinical isolates by using a comparative omics approach. Whole-genome sequencing showed that historic and recent isolates of B. pertussis have substantial variation in genome organization and form separate phylogenetic clusters. Subsequent RNA sequence analysis and liquid chromatography with mass tandem spectrometry analyses showed that these variations translated into discretely separated transcriptomic and proteomic profiles. When compared with historic strains, recent isolates showed increased expression of flagellar genes and genes involved in lipopolysaccharide biosynthesis and decreased expression of polysaccharide capsule genes. Compared with reference strain Tohama I, all strains had increased expression and production of the type III secretion system apparatus. We detected the potential link between observed effects and insertion sequence element-induced changes in gene context only for a few genes.

Rapid diversification of Pseudomonas aeruginosa in cystic fibrosis lung-like conditions.

Chronic infection of the cystic fibrosis (CF) airway by the opportunistic pathogen Pseudomonas aeruginosa is the leading cause of morbidity and mortality for adult CF patients. Prolonged infections are accompanied by adaptation of P. aeruginosa to the unique conditions of the CF lung environment, as well as marked diversification of the pathogen into phenotypically and genetically distinct strains that can coexist for years within a patient. Little is known, however, about the causes of this diversification and its impact on patient health. Here, we show experimentally that, consistent with ecological theory of diversification, the nutritional conditions of the CF airway can cause rapid and extensive diversification of P. aeruginosa Mucin, the substance responsible for the increased viscosity associated with the thick mucus layer in the CF airway, had little impact on within-population diversification but did promote divergence among populations. Furthermore, in vitro evolution recapitulated traits thought to be hallmarks of chronic infection, including reduced motility and increased biofilm formation, and the range of phenotypes observed in a collection of clinical isolates. Our results suggest that nutritional complexity and reduced dispersal can drive evolutionary diversification of P. aeruginosa independent of other features of the CF lung such as an active immune system or the presence of competing microbial species. We suggest that diversification, by generating extensive phenotypic and genetic variation on which selection can act, may be a key first step in the development of chronic infections.

Natural selection drives population divergence for local adaptation in a wheat pathogen.

Evolution favors the emergence of locally-adapted optimum phenotypes that are likely to differ across a wide array of environmental conditions. The emergence of favorable adaptive characteristics is accelerated in agricultural pathogens due to the unique properties of agro-ecosystems. We performed a QST - FST comparison using 164 strains of Parastagonospora nodorum sampled from eight global field populations to disentangle the predominant evolutionary forces driving population divergence in a wheat pathogen. We used digital image analysis to obtain quantitative measurements of growth rate and melanization at different temperatures and under different fungicide concentrations in a common garden experiment. FST measures were based on complete genome sequences obtained for all 164 isolates. Our analyses indicated that all measured traits were under selection. Growth rates at 18 °C and 24 °C were under stabilizing selection (QST < FST), while diversifying selection (QST > FST) was the predominant evolutionary force affecting growth under fungicide and high temperature stress. Stabilizing selection (QST < FST) was the predominant force affecting melanization across the different environments. Melanin production increased at 30 °C but was negatively correlated with higher growth rates, consistent with a trade-off under heat stress. Our results demonstrate that global populations of P. nodorum possess significant evolutionary potential to adapt to changing local conditions, including warmer temperatures and applications of fungicides.

On the probability of strain invasion in endemic settings: Accounting for individual heterogeneity and control in multi-strain dynamics.

Pathogen evolution is an imminent threat to global health that has warranted, and duly received, considerable attention within the medical, microbiological and modelling communities. Outbreaks of new pathogens are often ignited by the emergence and transmission of mutant variants descended from wild-type strains circulating in the community. In this work we investigate the stochastic dynamics of the emergence of a novel disease strain, introduced into a population in which it must compete with an existing endemic strain. In analogy with past work on single-strain epidemic outbreaks, we apply a branching process approximation to calculate the probability that the new strain becomes established. As expected, a critical determinant of the survival prospects of any invading strain is the magnitude of its reproduction number relative to that of the background endemic strain. Whilst in most circumstances this ratio must exceed unity in order for invasion to be viable, we show that differential control scenarios can lead to less-fit novel strains invading populations hosting a fitter endemic one. This analysis and the accompanying findings will inform our understanding of the mechanisms that have led to past instances of successful strain invasion, and provide valuable lessons for thwarting future drug-resistant strain incursions.

An evolutionary framework for host shifts - jumping ships for survival.

Host jumping is a process by which pathogens settle in new host groups. It is a cornerstone in the evolution of pathogens, as it leads to pathogen diversification. It is unsurprising that host jumping is observed in facultative pathogens, as they can reproduce even if they kill their hosts. However, host jumps were thought to be rare in obligate biotrophic pathogens, but molecular phylogenetics has revealed that the opposite is true. Here, I review some concepts and recent findings and present several hypotheses on the matter. In short, pathogens evolve and diversify via host jumps, followed by radiation, specialisation and speciation. Host jumps are facilitated by, for example, effector innovations, stress, compatible pathogens and physiological similarities. Host jumping, subsequent establishment, and speciation takes place rapidly - within centuries and millennia rather than over millions of years. If pathogens are unable to evolve into neutral or mutualistic interactions with their hosts, they will eventually be removed from the host population, despite balancing trade-offs. Thus, generally, plant pathogens only survive in the course of evolution if they jump hosts. This is also reflected by the diversity patterns observed in many genera of plant pathogens, where it leads to a mosaic pattern of host groups over time, in which the original host group becomes increasingly obscure.

Controlled fire use in early humans might have triggered the evolutionary emergence of tuberculosis.

Tuberculosis (TB) is caused by the Mycobacterium tuberculosis complex (MTBC), a wildly successful group of organisms and the leading cause of death resulting from a single bacterial pathogen worldwide. It is generally accepted that MTBC established itself in human populations in Africa and that animal-infecting strains diverged from human strains. However, the precise causal factors of TB emergence remain unknown. Here, we propose that the advent of controlled fire use in early humans created the ideal conditions for the emergence of TB as a transmissible disease. This hypothesis is supported by mathematical modeling together with a synthesis of evidence from epidemiology, evolutionary genetics, and paleoanthropology.