Host genotype and genetic diversity shape the evolution of a novel bacterial infection

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.

Bdelloid rotifers deploy horizontally acquired biosynthetic genes against a fungal pathogen

Coevolutionary antagonism generates relentless selection that can favour genetic exchange, including transfer of antibiotic synthesis and resistance genes among bacteria, and sexual recombination of disease resistance alleles in eukaryotes. We report an unusual link between biological conflict and DNA transfer in bdelloid rotifers, microscopic animals whose genomes show elevated levels of horizontal gene transfer from non-metazoan taxa. When rotifers were challenged with a fungal pathogen, horizontally acquired genes were over twice as likely to be upregulated as other genes - a stronger enrichment than observed for abiotic stressors. Among hundreds of upregulated genes, the most markedly overrepresented were clusters resembling bacterial polyketide and nonribosomal peptide synthetases that produce antibiotics. Upregulation of these clusters in a pathogen-resistant rotifer species was nearly ten times stronger than in a susceptible species. By acquiring, domesticating, and expressing non-metazoan biosynthetic pathways, bdelloids may have evolved to resist natural enemies using antimicrobial mechanisms absent from other animals.

Emergence failure of early epidemics: A mathematical modeling approach

Epidemic or pathogen emergence is the phenomenon by which a poorly transmissible pathogen finds its evolutionary pathway to become a mutant that can cause an epidemic. Many mathematical models of pathogen emergence rely on branching processes. Here, we discuss pathogen emergence using Markov chains, for a more tractable analysis, generalizing previous work by Kendall and Bartlett about disease invasion. We discuss the probability of emergence failure for early epidemics, when the number of infected individuals is small and the number of the susceptible individuals is virtually unlimited. Our formalism addresses both directly transmitted and vector-borne diseases, in the cases where the original pathogen is 1) one step-mutation away from the epidemic strain, and 2) undergoing a long chain of neutral mutations that do not change the epidemiology. We obtain analytic results for the probabilities of emergence failure and two features transcending the transmission mechanism. First, the reproduction number of the original pathogen is determinant for the probability of pathogen emergence, more important than the mutation rate or the transmissibility of the emerged pathogen. Second, the probability of mutation within infected individuals must be sufficiently high for the pathogen undergoing neutral mutations to start an epidemic, the mutation threshold depending again on the basic reproduction number of the original pathogen. Finally, we discuss the parameterization of models of pathogen emergence, using SARS-CoV1 as an example of zoonotic emergence and HIV as an example for the emergence of drug resistance. We also discuss assumptions of our models and implications for epidemiology.

Incomplete immunity in a natural animal-microbiota interaction selects for higher pathogen virulence

Incomplete immunity in recovered hosts is predicted to favor more virulent pathogens upon re-infection in the population.1 The microbiota colonizing animals can generate a similarly long-lasting, partial immune response, allowing for infection but dampened disease severity.2 We tracked the evolutionary trajectories of a widespread pathogen (Pseudomonas aeruginosa), experimentally passaged through populations of nematodes immune-primed by a natural microbiota member (P. berkeleyensis). This bacterium can induce genes regulated by a mitogen-activated protein kinase (MAPK) signaling pathway effective at conferring protection against pathogen-induced death despite infection.3 Across host populations, this incomplete immunity selected for pathogens more than twice as likely to kill as those evolved in non-primed (i.e., naive) or immune-compromised (mutants with a knockout of the MAPK ortholog) control populations. Despite the higher virulence, pathogen molecular evolution in immune-primed hosts was slow and constrained. In comparison, evolving pathogens in immune-compromised hosts were characterized by substantial genomic differentiation and attenuated virulence. These findings directly attribute the incomplete host immunity induced from microbiota as a significant force shaping the virulence and evolutionary dynamics of novel infectious diseases.

Pathogen evolution: Protective microbes act as a double-edged sword

Vaccines and infection can sometimes cause incomplete immunity, which allows for pathogen re-infection with decreased disease severity but also contributes to the evolution of pathogen virulence. A new study demonstrates that incomplete immunity from resident protective microbes results in similar evolutionary trajectories.

Genetic variants associated with hantavirus infection in a reservoir host are related to regulation of inflammation and immune surveillance

The immunogenetics of wildlife populations influence the epidemiology and evolutionary dynamic of the host-pathogen system. Profiling immune gene diversity present in wildlife may be especially important for those species that, while not at risk of disease or extinction themselves, are host to diseases that are a threat to humans, other wildlife, or livestock. Hantaviruses (genus: Orthohantavirus) are globally distributed zoonotic RNA viruses with pathogenic strains carried by a diverse group of rodent hosts. The marsh rice rat (Oryzomys palustris) is the reservoir host of Orthohantavirus bayoui, a hantavirus that causes fatal cases of hantavirus cardiopulmonary syndrome in humans. We performed a genome wide association study (GWAS) using the rice rat "immunome" (i.e., all exons related to the immune response) to identify genetic variants associated with infection status in wild-caught rice rats naturally infected with their endemic strain of hantavirus. First, we created an annotated reference genome using 10× Chromium Linked Reads sequencing technology. This reference genome was used to create custom baits which were then used to target enrich prepared rice rat libraries (n = 128) and isolate their immunomes prior to sequencing. Top SNPs in the association test were present in four genes (Socs5, Eprs, Mrc1, and Il1f8) which have not been previously implicated in hantavirus infections. However, these genes correspond with other loci or pathways with established importance in hantavirus susceptibility or infection tolerance in reservoir hosts: the JAK/STAT, MHC, and NFκB. These results serve as informative markers for future exploration and highlight the importance of immune pathways that repeatedly emerge across hantavirus systems. Our work aids in creating cross-species comparisons for better understanding mechanisms of genetic susceptibility and host-pathogen coevolution in hantavirus systems.

Intrahost evolution of the gut microbiota

A massive number of microorganisms, belonging to different species, continuously divide inside the guts of animals and humans. The large size of these communities and their rapid division times imply that we should be able to watch microbial evolution in the gut in real time, in a similar manner to what has been done in vitro. Here, we review recent findings on how natural selection shapes intrahost evolution (also known as within-host evolution), with a focus on the intestines of mice and humans. The microbiota of a healthy host is not as static as initially thought from the information measured at only one genomic marker. Rather, the genomes of each gut-colonizing species can be highly dynamic, and such dynamism seems to be related to the microbiota species diversity. Genetic and bioinformatic tools, and analysis of time series data, allow quantification of the selection strength on emerging mutations and horizontal transfer events in gut ecosystems. The drivers and functional consequences of gut evolution can now begin to be grasped. The rules of this intrahost microbiota evolution, and how they depend on the biology of each species, need to be understood for more effective development of microbiota therapies to help maintain or restore host health.

Microbes are potential key players in the evolution of life histories and aging in Caenorhabditis elegans

Microbes can have profound effects on host fitness and health and the appearance of late-onset diseases. Host-microbe interactions thus represent a major environmental context for healthy aging of the host and might also mediate trade-offs between life-history traits in the evolution of host senescence. Here, we have used the nematode Caenorhabditis elegans to study how host-microbe interactions may modulate the evolution of life histories and aging. We first characterized the effects of two non-pathogenic and one pathogenic Escherichia coli strains, together with the pathogenic Serratia marcescens DB11 strain, on population growth rates and survival of C. elegans from five different genetic backgrounds. We then focused on an outbred C. elegans population, to understand if microbe-specific effects on the reproductive schedule and in traits such as developmental rate and survival were also expressed in the presence of males and standing genetic variation, which could be relevant for the evolution of C. elegans and other nematode species in nature. Our results show that host-microbe interactions have a substantial host-genotype-dependent impact on the reproductive aging and survival of the nematode host. Although both pathogenic bacteria reduced host survival in comparison with benign strains, they differed in how they affected other host traits. Host fertility and population growth rate were affected by S. marcescens DB11 only during early adulthood, whereas this occurred at later ages with the pathogenic E. coli IAI1. In both cases, these effects were largely dependent on the host genotypes. Given such microbe-specific genotypic differences in host life history, we predict that the evolution of reproductive schedules and senescence might be critically contingent on host-microbe interactions in nature.

Multi-Host Pathogen Staphylococcus aureus-Epidemiology, Drug Resistance and Occurrence in Humans and Animals in Poland

Staphylococcus aureus is a drug resistant pathogen with zoonotic potential commonly isolated from humans and animals. The aim of this study was to compare the occurrence of drug resistance, resistance genes, sequence types (STs), and genotypes of S. aureus isolated from humans, livestock, and wildlife in eastern Poland. A high percentage of isolates resistant to penicillin (63%), erythromycin (39%), clindamycin (37%), tetracycline (31%), and methicillin (MRSA-19%), an intermediate resistant to vancomycin (VISA-13%), and a multidrug resistant (MDR-39%) was obtained. Multilocus sequence typing analysis showed the presence of 35 different STs (with dominance ST 15, ST 45, ST 7, and ST 582 in human, and ST 398 and ST 8139 in porcine and cattle isolates, respectively), including 9 new ones that had never been reported before (ST 8133-8141). Identical genotypic patterns were detected among porcine and cattle isolates as well as from humans and cattle. A high percentage of MDR, MRSA, and VISA in humans and livestock combined with the presence of the same genotypes among S. aureus isolated from human and cattle indicates the circulation of strains common in the region and their zoonotic potential. There is a need to develop new strategies to counteract this phenomenon according to the One Health policy.

Interspecific host competition and parasite virulence evolution

Virulence, the harm to hosts caused by parasite infection, can be selected for by several ecological factors acting synergistically or antagonistically. Here, we focus on the potential for interspecific host competition to shape virulence through such a network of effects. We first summarize how host natural mortality, body mass changes, population density and community diversity affect virulence evolution. We then introduce an initial conceptual framework highlighting how these host factors, which change during host competition, may drive virulence evolution via impacts on life-history trade-offs. We argue that the multi-faceted nature of both interspecific host competition and virulence evolution still requires consideration and experimentation to disentangle contrasting mechanisms. It also necessitates a differential treatment for parasites with various transmission strategies. However, such a comprehensive approach focusing on the role of interspecific host competition is essential to understand the processes driving the evolution of virulence in a tangled bank.

Virulence evolution during a naturally occurring parasite outbreak

Virulence, the degree to which a pathogen harms its host, is an important but poorly understood aspect of host-pathogen interactions. Virulence is not static, instead depending on ecological context and potentially evolving rapidly. For instance, at the start of an epidemic, when susceptible hosts are plentiful, pathogens may evolve increased virulence if this maximizes their intrinsic growth rate. However, if host density declines during an epidemic, theory predicts evolution of reduced virulence. Although well-studied theoretically, there is still little empirical evidence for virulence evolution in epidemics, especially in natural settings with native host and pathogen species. Here, we used a combination of field observations and lab assays in the Daphnia-Pasteuria model system to look for evidence of virulence evolution in nature. We monitored a large, naturally occurring outbreak of Pasteuria ramosa in Daphnia dentifera, where infection prevalence peaked at ~ 40% of the population infected and host density declined precipitously during the outbreak. In controlled infections in the lab, lifespan and reproduction of infected hosts was lower than that of unexposed control hosts and of hosts that were exposed but not infected. We did not detect any significant changes in host resistance or parasite infectivity, nor did we find evidence for shifts in parasite virulence (quantified by host lifespan and number of clutches produced by hosts). However, over the epidemic, the parasite evolved to produce significantly fewer spores in infected hosts. While this finding was unexpected, it might reflect previously quantified tradeoffs: parasites in high mortality (e.g., high predation) environments shift from vegetative growth to spore production sooner in infections, reducing spore yield. Future studies that track evolution of parasite spore yield in more populations, and that link those changes with genetic changes and with predation rates, will yield better insight into the drivers of parasite evolution in the wild.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10682-022-10169-6.

Time Course RNA-seq Reveals Soybean Responses against Root-Lesion Nematode and Resistance Players

Pratylenchus brachyurus causes serious damage to soybean production and other crops worldwide. Plant molecular responses to RLN infection remain largely unknown and no resistance genes have been identified in soybean. In this study, we analyzed molecular responses to RLN infection in moderately resistant BRSGO (Chapadões-BRS) and susceptible TMG115 RR (TMG) Glycine max genotypes. Differential expression analysis revealed two stages of response to RLN infection and a set of differentially expressed genes (DEGs) in the first stage suggested a pattern-triggered immunity (PTI) in both genotypes. The divergent time-point of DEGs between genotypes was observed four days post-infection, which included the activation of mitogen-activated protein kinase (MAPK) and plant-pathogen interaction genes in the BRS, suggesting the occurrence of an effector-triggered immunity response (ETI) in BRS. The co-expression analyses combined with single nucleotide polymorphism (SNP) uncovered a key element, a transcription factor phytochrome-interacting factor (PIF7) that is a potential regulator of moderate resistance to RLN infection. Two genes for resistance-related leucine-rich repeat (LRR) proteins were found as BRS-specific expressed genes. In addition, alternative splicing analysis revealed an intron retention in a myo-inositol oxygenase (MIOX) transcript, a gene related to susceptibility, may cause a loss of function in BRS.

Transgenerational genomic analyses reveal allelic oscillation and purifying selection in a gut parasite Nosema ceranae

Standing genetic variation is the predominant source acted on by selection. Organisms with high genetic diversity generally show faster responses toward environmental change. Nosema ceranae is a microsporidian parasite of honey bees, infecting midgut epithelial cells. High genetic diversity has been found in this parasite, but the mechanism for the parasite to maintain this diversity remains unclear. This study involved continuous inoculation of N. ceranae to honey bees. We found that the parasites slowly increased genetic diversity over three continuous inoculations. The number of lost single nucleotide variants (SNVs) was balanced with novel SNVs, which were mainly embedded in coding regions. Classic allele frequency oscillation was found at the regional level along the genome, and the associated genes were enriched in apoptosis regulation and ATP binding. The ratio of synonymous and non-synonymous substitution suggests a purifying selection, and our results provide novel insights into the evolutionary dynamics in microsporidian parasites.

Trade-offs in defence to pathogen species revealed in expanding nematode populations

Many host organisms live in polymicrobial environments and must respond to a diversity of pathogens. The degree to which host defences towards one pathogen species affect susceptibility to others is unclear. We used a panel of Caenorhabditis elegans nematode isolates to test for natural genetic variation in fitness costs of immune upregulation and pathogen damage, as well as for trade-offs in defence against two pathogen species, Staphylococcus aureus and Pseudomonas aeruginosa. We examined the fitness impacts of transient pathogen exposure (pathogen damage and immune upregulation) or exposure to heat-killed culture (immune upregulation only) by measuring host population sizes, which allowed us to simultaneously capture changes in reproductive output, developmental time and survival. We found significant decreases in population sizes for hosts exposed to live versus heat-killed S. aureus and found increased reproductive output after live P. aeruginosa exposure, compared with the corresponding heat-killed challenge. Nematode isolates with relatively higher population sizes after live P. aeruginosa infection produced fewer offspring after live S. aureus challenge. These findings reveal that wild C. elegans genotypes display a trade-off in defences against two distinct pathogen species that are evident in subsequent generations.

Genetic diversity and disease: The past, present, and future of an old idea

Why do infectious diseases erupt in some host populations and not others? This question has spawned independent fields of research in evolution, ecology, public health, agriculture, and conservation. In the search for environmental and genetic factors that predict variation in parasitism, one hypothesis stands out for its generality and longevity: genetically homogeneous host populations are more likely to experience severe parasitism than genetically diverse populations. In this perspective piece, I draw on overlapping ideas from evolutionary biology, agriculture, and conservation to capture the far-reaching implications of the link between genetic diversity and disease. I first summarize the development of this hypothesis and the results of experimental tests. Given the convincing support for the protective effect of genetic diversity, I then address the following questions: (1) Where has this idea been put to use, in a basic and applied sense, and how can we better use genetic diversity to limit disease spread? (2) What new hypotheses does the established disease-diversity relationship compel us to test? I conclude that monitoring, preserving, and augmenting genetic diversity is one of our most promising evolutionarily informed strategies for buffering wild, domesticated, and human populations against future outbreaks.

Population size impacts host-pathogen coevolution

Ongoing host-pathogen interactions are characterized by rapid coevolutionary changes forcing species to continuously adapt to each other. The interacting species are often defined by finite population sizes. In theory, finite population size limits genetic diversity and compromises the efficiency of selection owing to genetic drift, in turn constraining any rapid coevolutionary responses. To date, however, experimental evidence for such constraints is scarce. The aim of our study was to assess to what extent population size influences the dynamics of host-pathogen coevolution. We used Caenorhabditus elegans and its pathogen Bacillus thuringiensis as a model for experimental coevolution in small and large host populations, as well as in host populations which were periodically forced through a bottleneck. By carefully controlling host population size for 23 host generations, we found that host adaptation was constrained in small populations and to a lesser extent in the bottlenecked populations. As a result, coevolution in large and small populations gave rise to different selection dynamics and produced different patterns of host-pathogen genotype-by-genotype interactions. Our results demonstrate a major influence of host population size on the ability of the antagonists to co-adapt to each other, thereby shaping the dynamics of antagonistic coevolution.

Increased Pathogenicity of the Nematophagous Fungus Drechmeria coniospora Following Long-Term Laboratory Culture

Domestication provides a window into adaptive change. Over the course of 2 decades of laboratory culture, a strain of the nematode-specific fungus Drechmeria coniospora became more virulent during its infection of Caenorhabditis elegans. Through a close comparative examination of the genome sequences of the original strain and its more pathogenic derivative, we identified a small number of non-synonymous mutations in protein-coding genes. In one case, the mutation was predicted to affect a gene involved in hypoxia resistance and we provide direct corroborative evidence for such an effect. The mutated genes with functional annotation were all predicted to impact the general physiology of the fungus and this was reflected in an increased in vitro growth, even in the absence of C. elegans. While most cases involved single nucleotide substitutions predicted to lead to a loss of function, we also observed a predicted restoration of gene function through deletion of an extraneous tandem repeat. This latter change affected the regulatory subunit of a cAMP-dependent protein kinase. Remarkably, we also found a mutation in a gene for a second protein of the same, protein kinase A, pathway. Together, we predict that they result in a stronger repression of the pathway for given levels of ATP and adenylate cyclase activity. Finally, we also identified mutations in a few lineage-specific genes of unknown function that are candidates for factors that influence virulence in a more direct manner.

Trends in the molecular epidemiology and population genetics of emerging Sporothrix species

Sporothrix (Ophiostomatales) comprises species that are pathogenic to humans and other mammals as well as environmental fungi. Developments in molecular phylogeny have changed our perceptions about the epidemiology, host-association, and virulence of Sporothrix. The classical agent of sporotrichosis, Sporothrix schenckii, now comprises several species nested in a clinical clade with S. brasiliensis, S. globosa, and S. luriei. To gain a more precise view of outbreaks dynamics, structure, and origin of genetic variation within and among populations of Sporothrix, we applied three sets of discriminatory AFLP markers (#3 EcoRI-GA/MseI-TT, #5 EcoRI-GA/MseI-AG, and #6 EcoRI-TA/MseI-AA) and mating-type analysis to a large collection of human, animal and environmental isolates spanning the major endemic areas. A total of 451 polymorphic loci were amplified in vitro from 188 samples, and revealed high polymorphism information content (PIC = 0.1765-0.2253), marker index (MI = 0.0001-0.0002), effective multiplex ratio (E = 15.1720-23.5591), resolving power (Rp = 26.1075-40.2795), discriminating power (D = 0.9766-0.9879), expected heterozygosity (H = 0.1957-0.2588), and mean heterozygosity (Havp  = 0.000007-0.000009), demonstrating the effectiveness of AFLP markers to speciate Sporothrix. Analysis using the program structure indicated three genetic clusters matching S. brasiliensis (population 1), S. schenckii (population 2), and S. globosa (population 3), with the presence of patterns of admixture amongst all populations. AMOVA revealed highly structured clusters (PhiPT = 0.458-0.484, P < 0.0001), with roughly equivalent genetic variability within (46-48 %) and between (52-54 %) populations. Heterothallism was the exclusive mating strategy, and the distributions of MAT1-1 or MAT1-2 idiomorphs were not significantly skewed (1:1 ratio) for S. schenckii (χ2 = 2.522; P = 0.1122), supporting random mating. In contrast, skewed distributions were found for S. globosa (χ2 = 9.529; P = 0.0020) with a predominance of MAT1-1 isolates, and regional differences were highlighted for S. brasiliensis with the overwhelming occurrence of MAT1-2 in Rio de Janeiro (χ2 = 14.222; P = 0.0002) and Pernambuco (χ2 = 7.364; P = 0.0067), in comparison to a higher prevalence of MAT1-1 in the Rio Grande do Sul (χ2 = 7.364; P = 0.0067). Epidemiological trends reveal the geographic expansion of cat-transmitted sporotrichosis due to S. brasiliensis via founder effect. These data support Rio de Janeiro as the centre of origin that has led to the spread of this disease to other regions in Brazil. Our ability to reconstruct the source, spread, and evolution of the ongoing outbreaks from molecular data provides high-quality information for decision-making aimed at mitigating the progression of the disease. Other uses include surveillance, rapid diagnosis, case connectivity, and guiding access to appropriate antifungal treatment.