Pneumococcal Colonization and Virulence Factors Identified Via Experimental Evolution in Infection Models

Spontaneous mutations and mutational responses to penicillin treatment in the bacterial pathogen Streptococcus pneumoniae D39

Bacteria with functional DNA repair systems are expected to have low mutation rates due to strong natural selection for genomic stability. However, our study of the wild-type Streptococcus pneumoniae D39, a pathogen responsible for many common diseases, revealed a high spontaneous mutation rate of 0.02 per genome per cell division in mutation-accumulation (MA) lines. This rate is orders of magnitude higher than that of other non-mutator bacteria and is characterized by a high mutation bias in the A/T direction. The high mutation rate may have resulted from a reduction in the overall efficiency of selection, conferred by the tiny effective population size in nature. In line with this, S. pneumoniae D39 also exhibited the lowest DNA mismatch-repair (MMR) efficiency among bacteria. Treatment with the antibiotic penicillin did not elevate the mutation rate, as penicillin did not induce DNA damage and S. pneumoniae lacks a stress response pathway. Our findings suggested that the MA results are applicable to within-host scenarios and provide insights into pathogen evolution.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00220-6.

Identification of two distinct phylogenomic lineages and model strains for the understudied cystic fibrosis lung pathogen Burkholderia multivorans

Burkholderia multivorans is the dominant Burkholderia pathogen recovered from lung infection in people with cystic fibrosis. However, as an understudied pathogen there are knowledge gaps in relation to its population biology, phenotypic traits and useful model strains. A phylogenomic study of B. multivorans was undertaken using a total of 283 genomes, of which 73 were sequenced and 49 phenotypically characterized as part of this study. Average nucleotide identity analysis (ANI) and phylogenetic alignment of core genes demonstrated that the B. multivorans population separated into two distinct evolutionary clades, defined as lineage 1 (n=58 genomes) and lineage 2 (n=221 genomes). To examine the population biology of B. multivorans, a representative subgroup of 77 B. multivorans genomes (28 from the reference databases and the 49 novel short-read genome sequences) were selected based on multilocus sequence typing (MLST), isolation source and phylogenetic placement criteria. Comparative genomics was used to identify B. multivorans lineage-specific genes - ghrB_1 in lineage 1 and glnM_2 in lineage 2 - and diagnostic PCRs targeting them were successfully developed. Phenotypic analysis of 49 representative B. multivorans strains showed considerable inter-strain variance, but the majority of the isolates tested were motile and capable of biofilm formation. A striking absence of B. multivorans protease activity in vitro was observed, but no lineage-specific phenotypic differences were demonstrated. Using phylogenomic and phenotypic criteria, three model B. multivorans CF strains were identified, BCC0084 (lineage 1), BCC1272 (lineage 2a) and BCC0033 lineage 2b, and their complete genome sequences determined. B. multivorans CF strains BCC0033 and BCC0084, and the environmental reference strain, ATCC 17616, were all capable of short-term survival within a murine lung infection model. By mapping the population biology, identifying lineage-specific PCRs and model strains, we provide much needed baseline resources for future studies of B. multivorans.

Metabolic Network Models of the Gardnerella Pangenome Identify Key Interactions with the Vaginal Environment

Gardnerella is the primary pathogenic bacterial genus present in the polymicrobial condition known as bacterial vaginosis (BV). Despite BV's high prevalence and associated chronic and acute women's health impacts, the Gardnerella pangenome is largely uncharacterized at both the genetic and functional metabolic levels. Here, we used genome-scale metabolic models to characterize in silico the Gardnerella pangenome metabolic content. We also assessed the metabolic functional capacity in a BV-positive cervicovaginal fluid context. The metabolic capacity varied widely across the pangenome, with 38.15% of all reactions being core to the genus, compared to 49.60% of reactions identified as being unique to a smaller subset of species. We identified 57 essential genes across the pangenome via in silico gene essentiality screens within two simulated vaginal metabolic environments. Four genes, gpsA, fas, suhB, and psd, were identified as core essential genes critical for the metabolic function of all analyzed bacterial species of the Gardnerella genus. Further understanding these core essential metabolic functions could inform novel therapeutic strategies to treat BV. Machine learning applied to simulated metabolic network flux distributions showed limited clustering based on the sample isolation source, which further supports the presence of extensive core metabolic functionality across this genus. These data represent the first metabolic modeling of the Gardnerella pangenome and illustrate strain-specific interactions with the vaginal metabolic environment across the pangenome. IMPORTANCE Bacterial vaginosis (BV) is the most common vaginal infection among reproductive-age women. Despite its prevalence and associated chronic and acute women's health impacts, the diverse bacteria involved in BV infection remain poorly characterized. Gardnerella is the genus of bacteria most commonly and most abundantly represented during BV. In this paper, we use metabolic models, which are a computational representation of the possible functional metabolism of an organism, to investigate metabolic conservation, gene essentiality, and pathway utilization across 110 Gardnerella strains. These models allow us to investigate in silico how strains may differ with respect to their metabolic interactions with the vaginal-host environment.

Influence of Streptococcus pneumoniae Within-Strain Population Diversity on Virulence and Pathogenesis

The short generation time of many bacterial pathogens allows the accumulation of de novo mutations during routine culture procedures used for the preparation and propagation of bacterial stocks. Taking the major human pathogen Streptococcus pneumoniae as an example, we sought to determine the influence of standard laboratory handling of microbes on within-strain genetic diversity and explore how these changes influence virulence characteristics and experimental outcomes. A single culture of S. pneumoniae D39 grown overnight resulted in the enrichment of previously rare genotypes present in bacterial freezer stocks and the introduction of new variation to the bacterial population through the acquisition of mutations. A comparison of D39 stocks from different laboratories demonstrated how changes in bacterial population structure taking place during individual culture events can cumulatively lead to fixed, divergent change that profoundly alters virulence characteristics. The passage of D39 through mouse models of infection, a process used to standardize virulence, resulted in the enrichment of high-fitness genotypes that were originally rare (<2% frequency) in D39 culture collection stocks and the loss of previously dominant genotypes. In the most striking example, the selection of a <2%-frequency genotype carrying a mutation in sdhB, a gene thought to be essential for the establishment of lung infection, was associated with enhanced systemic virulence. Three separately passaged D39 cultures originating from the same frozen stocks showed considerable genetic divergence despite comparable virulence. IMPORTANCE Laboratory bacteriology involves the use of high-density cultures that we often assume to be clonal but that in reality are populations consisting of multiple genotypes at various abundances. We have demonstrated that the genetic structure of a single population of a widely used Streptococcus pneumoniae strain can be substantially altered by even short-term laboratory handling and culture and that, over time, this can lead to changes in virulence characteristics. Our findings suggest that caution should be applied when comparing data generated in different laboratories using the same strain but also when comparing data within laboratories over time. Given the dramatic reductions in the cost of next-generation sequencing technology in recent years, we advocate for the frequent sampling and sequencing of bacterial isolate collections.

The Yin and Yang of Pneumolysin During Pneumococcal Infection

Pneumolysin (PLY) is a pore-forming toxin produced by the human pathobiont Streptococcus pneumoniae, the major cause of pneumonia worldwide. PLY, a key pneumococcal virulence factor, can form transmembrane pores in host cells, disrupting plasma membrane integrity and deregulating cellular homeostasis. At lytic concentrations, PLY causes cell death. At sub-lytic concentrations, PLY triggers host cell survival pathways that cooperate to reseal the damaged plasma membrane and restore cell homeostasis. While PLY is generally considered a pivotal factor promoting S. pneumoniae colonization and survival, it is also a powerful trigger of the innate and adaptive host immune response against bacterial infection. The dichotomy of PLY as both a key bacterial virulence factor and a trigger for host immune modulation allows the toxin to display both "Yin" and "Yang" properties during infection, promoting disease by membrane perforation and activating inflammatory pathways, while also mitigating damage by triggering host cell repair and initiating anti-inflammatory responses. Due to its cytolytic activity and diverse immunomodulatory properties, PLY is integral to every stage of S. pneumoniae pathogenesis and may tip the balance towards either the pathogen or the host depending on the context of infection.

IL-6 Prevents Lung Macrophage Death and Lung Inflammation Injury by Inhibiting GSDME- and GSDMD-Mediated Pyroptosis during Pneumococcal Pneumosepsis

Streptococcus pneumoniae is a leading bacterial cause of a wide range of infections, and pneumococcal pneumosepsis causes high mortality in hosts infected with antibiotic-resistant strains and those who cannot resolve ongoing inflammation. The factors which influence the development and outcome of pneumosepsis are currently unclear. IL-6 is critical for maintaining immune homeostasis, and we determined that this cytokine is also essential for resisting pneumosepsis, as it inhibits macrophage pyroptosis and pyroptosis-related inflammation injury in the lung. IL-6 affected infection outcomes in mice and exerted a protective role, primarily via macrophages. We further found that IL-6 deficiency led to increased lung macrophage death and aggravated lung inflammation, and that exogenous administration of IL-6 protein could decrease macrophage death and alleviate lung tissue inflammation. IL-6 also protected Streptococcus pneumoniae-induced lung macrophage death and lung inflammation injury by inhibiting gasdermin E (GSDME)- and gasdermin D (GSDMD)-mediated pyroptosis. Together, these data reveal a novel mechanism for the development of pneumosepsis and the critical protective role of IL-6. These findings may assist in the early identification and treatment of pneumococcal pneumosepsis.

IMPORTANCE Pneumococcal pneumonia has been a significant cause of morbidity and mortality throughout human history. Failing to control pneumococcal pneumonia and resolve ongoing inflammation in a host can cause sepsis, namely pneumococcal pneumosepsis, and death ensues. Few theories have suggested an optimally therapeutic option for this infectious disease. The interleukin-6 (IL-6, a cytokine featuring pleiotropic activity) theory, proposed here, implies that IL-6 acts as a protector against pneumococcal pneumosepsis. IL-6 prevents lung macrophage death and lung inflammation injury by inhibiting a caspase-3-GSDME-mediated switch from apoptosis to pyroptosis and inhibiting caspase-1-GSDMD-mediated classic pyroptosis during pneumococcal pneumosepsis. Thus, IL-6 is an important determinant for controlling bacterial invasion and a homeostatic coordinator of pneumococcal pneumosepsis. This study clarifies a novel mechanism of occurrence and development of pneumonia and secondary sepsis following a Streptococcus pneumoniae infection. It is important for the early identification and treatment of pneumococcal pneumosepsis.

The glycerol-3-phosphate dehydrogenases GpsA and GlpD constitute the oxidoreductive metabolic linchpin for Lyme disease spirochete host infectivity and persistence in the tick

We have identified GpsA, a predicted glycerol-3-phosphate dehydrogenase, as a virulence factor in the Lyme disease spirochete Borrelia (Borreliella) burgdorferi: GpsA is essential for murine infection and crucial for persistence of the spirochete in the tick. B. burgdorferi has a limited biosynthetic and metabolic capacity; the linchpin connecting central carbohydrate and lipid metabolism is at the interconversion of glycerol-3-phosphate and dihydroxyacetone phosphate, catalyzed by GpsA and another glycerol-3-phosphate dehydrogenase, GlpD. Using a broad metabolomics approach, we found that GpsA serves as a dominant regulator of NADH and glycerol-3-phosphate levels in vitro, metabolic intermediates that reflect the cellular redox potential and serve as a precursor for lipid and lipoprotein biosynthesis, respectively. Additionally, GpsA was required for survival under nutrient stress, regulated overall reductase activity and controlled B. burgdorferi morphology in vitro. Furthermore, during in vitro nutrient stress, both glycerol and N-acetylglucosamine were bactericidal to B. burgdorferi in a GlpD-dependent manner. This study is also the first to identify a suppressor mutation in B. burgdorferi: a glpD deletion restored the wild-type phenotype to the pleiotropic gpsA mutant, including murine infectivity by needle inoculation at high doses, survival under nutrient stress, morphological changes and the metabolic imbalance of NADH and glycerol-3-phosphate. These results illustrate how basic metabolic functions that are dispensable for in vitro growth can be essential for in vivo infectivity of B. burgdorferi and may serve as attractive therapeutic targets.

Lyme disease (borreliosis) is the most common tick-borne disease in the Northern hemisphere and its prevalence is increasing. Borrelia burgdorferi, the etiological agent of Lyme disease, is an enzootic pathogen that alternates between a tick vector and vertebrate host. Humans are considered an incidental host after transmission of B. burgdorferi following the bite of an infected tick. The mechanisms by which B. burgdorferi persists in the Ixodid tick, transmits to a vertebrate host and establishes infection are not well understood. Therefore, identifying virulence factors and uncovering the pathogenic strategies in the spirochete remain important to address the public health concerns of Lyme disease. In this study, we identify an enzyme involved in three-carbon metabolism, GpsA, as a new virulence factor with an effect on persistence in ticks. GpsA and GlpD, another enzyme, constitute a bidirectional metabolic node connecting lipid biosynthesis and glycolysis, which serves as the linchpin for regulating carbon utilization for B. burgdorferi throughout its enzootic cycle. Disruption of this node causes a lethal metabolic imbalance revealing a potential therapeutic target for the treatment of Lyme disease.

Lysine245 plays a crucial role in stability and function of glycerol 3-phosphate dehydrogenase (Gpd1) in Saccharomyces cerevisiae

Glycerol 3-phosphate dehydrogenase (Gpd1 isoform) catalyzes the rate limiting step of glycerol synthesis and is a critical component of the osmo-responsive machinery in yeast. The three-dimensional structure of the enzyme is similar to the enzyme from many other organisms, including humans. A recent study with the human enzyme has proposed K120 (K152 in yeast) to be in the correct orientation for catalysis; K204 (K245 in yeast) is out of plane and is not a participant in the catalytic cycle. The current work was carried out to establish the role of K245 in the catalytic cycle of yeast Gpd1. K245A mutant was found to possess lower catalytic activity. Osmotically stressed cells expressing Gpd1 (K245A) showed no change in intracellular glycerol as compared with wild-type cells which showed ~60% increase. Fluorescence microscopy, native polyacrylamide gel electrophoresis (PAGE) analysis, fluorescence spectroscopy, and Thioflavin T spectrofluorimetry showed a relatively unstable, aggregation- and degradation-prone conformation for the mutant. In silico studies showed an aggregation "hotspot" around K245. This study establishes the requirement of K245 for conformational stability and functional adaptation of Gpd1 in Saccharomyces cerevisiae.