While Botrytis cinerea causes gray mold on many plants, its close relative, Botrytis fabae, is host-specifically infecting predominantly faba bean plants. To explore the basis for its narrow host range, a gapless genome sequence of B. fabae strain G12 (BfabG12) was generated. The BfabG12 genome encompasses 45.0 Mb, with 16 chromosomal telomere-to-telomere contigs that show high synteny and sequence similarity to the corresponding B. cinerea B05.10 (BcB0510) chromosomes. Compared to BcB0510, it is 6% larger, due to many AT-rich regions containing remnants of transposable elements, but encodes fewer genes (11,420 vs. 11,707), due to losses of chromosomal segments with up to 20 genes. The coding capacity of BfabG12 is further reduced by nearly 400 genes that had been inactivated by mutations leading to truncations compared to their BcB0510 orthologues. Several species-specific gene clusters for secondary metabolite biosynthesis with stage-specific expression were identified. Comparison of the proteins secreted during infection revealed high similarities, including 17 phytotoxic proteins that were detected in both species. Our data indicate that evolution of the host-specific B. fabae occurred from an ancestral pathogen with wide host range similar to B. cinerea and was accompanied by losses and degeneration of genes, thereby reducing its pathogenic flexibility.

AB5-type toxins are a diverse family of protein toxins composed of an enzymatic active (A) subunit and a pentameric delivery (B) subunit. Salmonella enterica serovar Typhi\'s typhoid toxin features two A subunits, CdtB and PltA, in complex with the B subunit PltB. Recently, it was shown that S. Typhi encodes a horizontally acquired B subunit, PltC, that also assembles with PltA/CdtB to produce a second form of typhoid toxin. S. Typhi therefore produces two AB5 toxins with the same A subunits but distinct B subunits, an evolutionary twist that is unique to typhoid toxin. Here, we show that, remarkably, the Salmonella bongori species independently evolved an analogous capacity to produce two typhoid toxins with distinct B subunits. S. bongori\'s alternate B subunit, PltD, is evolutionarily distant from both PltB and PltC and outcompetes PltB to form the predominant toxin. We show that, surprisingly, S. bongori elicits similar levels of CdtB-mediated intoxication as S. Typhi during infection of cultured human epithelial cells. This toxicity is exclusively due to the PltB toxin, and strains lacking pltD produce increased amounts of PltB toxin and exhibit increased toxicity compared to the wild type, suggesting that the acquisition of the PltD subunit potentially made S. bongori less virulent toward humans. Collectively, this study unveils a striking example of convergent evolution that highlights the importance of the poorly understood \"two-toxin\" paradigm for typhoid toxin biology and, more broadly, illustrates how the flexibility of A-B interactions has fueled the evolutionary diversification and expansion of AB5-type toxins.
OBJECTIVE: Typhoid toxin is an important Salmonella Typhi virulence factor and an attractive target for therapeutic interventions to combat typhoid fever. The recent discovery of a second version of this toxin has substantial implications for understanding S. Typhi pathogenesis and combating typhoid fever. In this study, we discover that a remarkably similar two-toxin paradigm evolved independently in Salmonella bongori, which strongly suggests that this is a critical aspect of typhoid toxin biology. We observe significant parallels between how the two toxins assemble and their capacity to intoxicate host cells during infection in S. Typhi and S. bongori, which provides clues to the biological significance of this unusual toxin arrangement. More broadly, AB5 toxins with diverse activities and mechanisms are essential virulence factors for numerous important bacterial pathogens. This study illustrates the capacity for novel A-B interactions to evolve and thus provides insight into how such a diverse arsenal of toxins might have emerged.

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.

Tellurite resistance gene clusters have been identified in numerous pathogenic bacteria, including clinical isolates of Escherichia coli. The rareness of tellurium in host organisms and the noncontaminated environment raises a question about the true functionality of tellurite resistance gene clusters in pathogenesis and their possible contribution to bacterial fitness. The study aims to point out the beneficial effects of the tellurite resistance gene cluster of pathogenic bacteria to survive in ROS-rich environments. Here, we analysed the bacterial response to oxidative stress conditions with and without tellurite resistance gene clusters, which are composed of terWY1XY2Y3 and terZABCDEF genes. By measuring the levels of protein carbonylation, lipid peroxidation, and expression changes of oxidative stress genes upon oxidative stress, we propose a tellurite resistance gene cluster contribution to the elimination of oxidative damage, potentially increasing fitness and resistance to reactive oxygen species during macrophage attack. We have shown a different beneficial effect of various truncated versions of the tellurite resistance gene cluster on cell survival. The terBCDEF genes increased the survival of E. coli strain MC4100 by 13.21%, terW and terZABCDEF by 10.09%, and terWY1XY2Y3 and terZABCDEF by 25.57%, respectively. The ability to survive tellurite treatment is the most significant at 44.8% in wild clinical strain KL53 compared to laboratory strain E. coli MC4100 due to a complete wild-type plasmid presence.

Pathogen evolution of drug resistance often occurs in a stepwise manner via the accumulation of multiple mutations that in combination have a non-additive impact on fitness, a phenomenon known as epistasis. The evolution of resistance via the accumulation of point mutations in the DHFR genes of Plasmodium falciparum (Pf) and Plasmodium vivax (Pv) has been studied extensively and multiple studies have shown epistatic interactions between these mutations determine the accessible evolutionary trajectories to highly resistant multiple mutations. Here, we simulated these evolutionary trajectories using a model of molecular evolution, parameterised using Rosetta Flex ddG predictions, where selection acts to reduce the target-drug binding affinity. We observe strong agreement with pathways determined using experimentally measured IC50 values of pyrimethamine binding, which suggests binding affinity is strongly predictive of resistance and epistasis in binding affinity strongly influences the order of fixation of resistance mutations. We also infer pathways directly from the frequency of mutations found in isolate data, and observe remarkable agreement with the most likely pathways predicted by our mechanistic model, as well as those determined experimentally. This suggests mutation frequency data can be used to intuitively infer evolutionary pathways, provided sufficient sampling of the population.

The major resistance gene BvCR4 recently bred into sugar beet hybrids provides a high level of resistance to Cercospora leaf spot caused by the fungal pathogen Cercospora beticola. The occurrence of pathogen strains that overcome BvCR4 was studied using field trials in Switzerland conducted under natural disease pressure. Virulence of a subset of these strains was evaluated in a field trial conducted under elevated artificial disease pressure. We created a new C. beticola reference genome and mapped whole genome sequences of 256 isolates collected in Switzerland and Germany. These were combined with virulence phenotypes to conduct three separate genome-wide association studies (GWAS) to identify candidate avirulence genes. We identified a locus associated with avirulence containing a putative avirulence effector gene named AvrCR4. All virulent isolates either lacked AvrCR4 or had nonsynonymous mutations within the gene. AvrCR4 was present in all 74 isolates from non-BvCR4 hybrids, whereas 33 of 89 isolates from BvCR4 hybrids carried a deletion. We also mapped genomic data from 190 publicly available US isolates to our new reference genome. The AvrCR4 deletion was found in only one of 95 unique isolates from non-BvCR4 hybrids in the United States. AvrCR4 presents a unique example of an avirulence effector in which virulent alleles have only recently emerged. Most likely these were selected out of standing genetic variation after deployment of BvCR4. Identification of AvrCR4 will enable real-time screening of C. beticola populations for the emergence and spread of virulent isolates.

We have extended previously published sets of simple sequence repeat markers for Synchytrium endobioticum, selected to be polymorphic for the German-standard isolates of pathotypes P1, P2, P6, P8, and P18. These markers also complement the extensive published information on DNA polymorphisms for the mitogenomes of Synchytrium endobioticum. This extended set of 35 markers representing 73 alleles differentiated 51 isolates from Europe and North America into three large, well-separated clusters and subclusters using dendrogram analysis, principal coordinates analysis (PCoA), and population substructure analysis using STRUCTURE 2.3.4 software. This suggests a limited number of introgressions of the wart disease pathogen into current potato growing areas, followed by recombination and admixture of populations through human activities. The new markers extend the published marker sets and are useful tools for future analyses of population structure and dynamics in Synchytrium endobioticum, which are necessary to understand the biology of the interaction between the pathogen and its potato host and to develop future control strategies.

Major histocompatibility complex (MHC) genes (referred to as human leukocyte antigen or HLA in humans) are a key component of vertebrate immune systems, coding for proteins which present antigens to T-cells. These genes are outstanding in their degree of polymorphism, with important consequences for human and animal health. The polymorphism is thought to arise from selection pressures imposed by pathogens on MHC allomorphs, which differ in their antigen-binding capacity. However, the existing theory has not considered MHC selection in relation to the formation of immune memory. In this paper, we argue that this omission limits our understanding of the evolution of MHC polymorphism and its role in disease. We review recent evidence that has emerged from the massive research effort related to the SARS-CoV-2 pandemics, and which provides new evidence for the role of MHC in shaping immune memory. We then discuss why the inclusion of immune memory within the existing theory may have non-trivial consequence for our understanding of the evolution of MHC polymorphism. Finally, we will argue that neglecting immune memory hinders our interpretation of empirical findings, and postulate that future studies focusing on pathogen-driven MHC selection would benefit from stratifying the available data according to the history of infection (and vaccination, if relevant).

OBJECTIVE: Polymicrobial infections are common. In chronic infections, the different pathogens may repeatedly interact, which could spur evolutionary dynamics with pathogens adapting to one another. Here, we explore the potential of Staphylococcus aureus to adapt to its competitor Pseudomonas aeruginosa. These two pathogens frequently co-occur, and P. aeruginosa is seen as the dominant species being able to displace S. aureus. We studied three different S. aureus strains and found that all became quickly resistant to inhibitory compounds secreted by P. aeruginosa. Our experimental evolution revealed strains-specific adaptations with three main factors contributing to resistance evolution: (i) overproduction of staphyloxanthin, a molecule protecting from oxidative stress; (ii) the formation of small colony variants also protecting from oxidative stress; and (iii) alterations of membrane transporters possibly reducing toxin uptake. Our results show that species interactions can change over time potentially favoring species co-existence, which in turn could affect disease progression and treatment options.

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.