Genetic diversity among offspring from archived Salmonella enterica ssp. enterica serovar typhimurium (Demerec Collection): in search of survival strategies

Genomic population structure associated with repeated escape of Salmonella enterica ATCC14028s from the laboratory into nature

Salmonella enterica serovar Typhimurium strain ATCC14028s is commercially available from multiple national type culture collections, and has been widely used since 1960 for quality control of growth media and experiments on fitness (“laboratory evolution”). ATCC14028s has been implicated in multiple cross-contaminations in the laboratory, and has also caused multiple laboratory infections and one known attempt at bioterrorism. According to hierarchical clustering of 3002 core gene sequences, ATCC14028s belongs to HierCC cluster HC20_373 in which most internal branch lengths are only one to three SNPs long. Many natural Typhimurium isolates from humans, domesticated animals and the environment also belong to HC20_373, and their core genomes are almost indistinguishable from those of laboratory strains. These natural isolates have infected humans in Ireland and Taiwan for decades, and are common in the British Isles as well as the Americas. The isolation history of some of the natural isolates confirms the conclusion that they do not represent recent contamination by the laboratory strain, and 10% carry plasmids or bacteriophages which have been acquired in nature by HGT from unrelated bacteria. We propose that ATCC14028s has repeatedly escaped from the laboratory environment into nature via laboratory accidents or infections, but the escaped micro-lineages have only a limited life span. As a result, there is a genetic gap separating HC20_373 from its closest natural relatives due to a divergence between them in the late 19^th century followed by repeated extinction events of escaped HC20_373.

Clades of closely related bacteria exist in nature. Individual isolates from such clades are often distinguishable by genomic sequencing because genomic sequence differences can be acquired over a few years due to neutral drift and natural selection. The evolution of laboratory strains is often largely frozen, physically due to storage conditions and genetically due to long periods of storage. Thus, laboratory strains can normally be readily distinguished from natural isolates because they show much less diversity. However, laboratory strain ATCC14028s shows modest levels of sequence diversity because it has been shipped around the world to multiple laboratories and is routinely used for analyses of laboratory evolution. Closely related natural isolates also exist, but their genetic diversity is not dramatically greater at the core genome level. Indeed, many scientists doubt that such isolates are natural, and interpret them as undetected contamination by the laboratory strain. We present data indicating that ATCC14028s has repeatedly escaped from the laboratory through inadvertent contamination of the environment, infection of technical staff and deliberate bioterrorism. The escapees survive in nature long enough that some acquire mobile genomic elements by horizontal gene transfer, but eventually they go extinct. As a result, even extensive global databases of natural isolates lack closely related isolates whose ancestors diverged from ATCC14028s within the last 100 years.

A Combination of Metagenomic and Cultivation Approaches Reveals Hypermutator Phenotypes within Vibrio cholerae-Infected Patients

Vibrio cholerae can cause a range of symptoms, from severe diarrhea to asymptomatic infection. Previous studies using whole-genome sequencing (WGS) of multiple bacterial isolates per patient showed that V. cholerae can evolve modest genetic diversity during symptomatic infection. To further explore the extent of V. cholerae within-host diversity, we applied culture-based WGS and metagenomics to a cohort of both symptomatic and asymptomatic cholera patients from Bangladesh. While metagenomics allowed us to detect more mutations in symptomatic patients, WGS of cultured isolates was necessary to detect V. cholerae diversity in asymptomatic carriers, likely due to their low V. cholerae load. Using both metagenomics and isolate WGS, we report three lines of evidence that V. cholerae hypermutators evolve within patients. First, we identified nonsynonymous mutations in V. cholerae DNA repair genes in 5 out of 11 patient metagenomes sequenced with sufficient coverage of the V. cholerae genome and in 1 of 3 patients with isolate genomes sequenced. Second, these mutations in DNA repair genes tended to be accompanied by an excess of intrahost single nucleotide variants (iSNVs). Third, these iSNVs were enriched in transversion mutations, a known hallmark of hypermutator phenotypes. While hypermutators appeared to generate mostly selectively neutral mutations, nonmutators showed signs of convergent mutation across multiple patients, suggesting V. cholerae adaptation within hosts. Our results highlight the power and limitations of metagenomics combined with isolate sequencing to characterize within-patient diversity in acute V. cholerae infections, while providing evidence for hypermutator phenotypes within cholera patients.

IMPORTANCE Pathogen evolution within patients can impact phenotypes such as drug resistance and virulence, potentially affecting clinical outcomes. V. cholerae infection can result in life-threatening diarrheal disease or asymptomatic infection. Here, we describe whole-genome sequencing of V. cholerae isolates and culture-free metagenomic sequencing from stool of symptomatic cholera patients and asymptomatic carriers. Despite the typically short duration of cholera, we found evidence for adaptive mutations in the V. cholerae genome that occur independently and repeatedly within multiple symptomatic patients. We also identified V. cholerae hypermutator phenotypes within several patients, which appear to generate mainly neutral or deleterious mutations. Our work sets the stage for future studies of the role of hypermutators and within-patient evolution in explaining the variation from asymptomatic carriage to symptomatic cholera.

Evaluations of CRC2631 toxicity, tumor colonization, and genetic stability in the TRAMP prostate cancer model

Conventional cancer chemotherapies are not fully efficacious and do not target tumors, leading to significant treatment-related morbidities. A number of genetically attenuated cancer-targeting bacteria are being developed to safely target tumors in vivo. Here we report the toxicological, tumor-targeting, and efficacy profiles of Salmonella enterica serovar Typhimurium CRC2631 in a syngeneic and autochthonous TRAMP model of aggressive prostate cancer. CRC2631 preferentially colonize primary and metastatic tumors in the TRAMP animals. In addition, longitudinal whole genome sequencing studies of CRC2631 recovered from prostate tumor tissues demonstrate that CRC2631 is genetically stable. Moreover, tumor-targeted CRC2631 generates an anti-tumor immune response. Combination of CRC2631 with checkpoint blockade reduces metastasis burden. Collectively, these findings demonstrate a potential for CRC2631 in cancer immunotherapy strategies.

Comparative stress response to food preservation conditions of ST19 and ST213 genotypes of Salmonella enterica serotype Typhimurium

The replacement of the most prevalent Salmonella enterica genotypes has been documented worldwide. Here we tested the hypothesis that the current prevalent sequence type ST213 of serotype Typhimurium in Mexico has a higher resistance to stressful food preservation conditions than the displaced sequence ST19. ST19 showed higher cell viability percentages than ST213 in osmotic (685 mM NaCl) and acidic (pH 3.5) stress conditions and in combination with refrigeration (4 °C) and ambient (≈22 °C) temperatures. Both genotypes showed the same poststress recovery growth. ST213 formed biofilm and filamentous cells (FCs) under stress, whereas ST19 did not. ST213 cells also showed higher motility. The capacity of ST213 to form FCs may explain its lower viability percentages when compared with ST19, i.e., ST213 cells divided less under stress conditions, but FCs had the same recovery capacity of ST19 cells. ST213 presented a higher unsaturated/saturated fatty acids ratio (0.5-0.6) than ST19 (0.2-0.5), which indicates higher membrane fluidity. The transcript levels of the rpoS gene were similar between genotypes under the experimental conditions employed. Biofilm formation, the generation of FCs, cell motility and membrane modification seem to make ST213 more resistant than ST19 to food preservation environments.

Salmonella Bacterial Monotherapy Reduces Autochthonous Prostate Tumor Burden in the TRAMP Mouse Model

Attenuated Salmonella typhimurium injected in the circulatory system of mammals selectively targets tumors. Using weekly intraperitoneal injections of attenuated Salmonella strain CRC2631, we tested for regression and/or inhibition of tumor development in the TRAMP prostate tumor mouse model, which utilizes SV40 early region expression for autochthonous formation of prostate tumors that progress into metastatic, poorly differentiated prostatic carcinomas in an immunocompetent murine model. Thirteen weekly intraperitoneal administrations of 10⁵–10⁷ CFU CRC2631 into 10 week old mice were well tolerated by the TRAMP model. Sacrifice and histological analysis of TRAMP prostates at 22 weeks indicated that Salmonella monotherapy at administrated levels decrease visible tumor size (>29%) but did not significantly inhibit previously described SV40 expression-driven TRAMP tumor progression to undifferentiated carcinomas when histologically examined. In conclusion, this work demonstrates baseline results for CRC2631 Salmonella monotherapy using the immunocompetent TRAMP prostate tumor model in preparation for study of combination therapies that resolve autochthonously generated TRAMP prostate tumors, further reduce tumor size, or inhibit prostate tumor progression.

Strains, Mechanism, and Perspective: Salmonella-Based Cancer Therapy

Recently, investigation of bacterial-based tumor therapy has regained focus due to progress in molecular, cellular, and microbial biology. Many bacteria such as Salmonella, Listeria, Escherichia, and Clostridium have proved to have tumor targeting and in some cases even tumor-destroying phenotypes. Furthermore, bacterial clinical treatments for cancer have been improved by combination with other therapeutic methods such as chemotherapeutic drugs and radioactive agents. Synthetic biology techniques have also driven the development of new bacterial-based cancer therapies. However, basic questions about the mechanisms of bacterial-mediated tumor targeting and destruction are still being elucidated. In this review, we focus on three tumor-therapeutic Salmonella models, the most intensively studied bacterial genus in this field. One of these Salmonella models is our Salmonella enterica serovar Typhimurium LT2 derived strain CRC2631, engineered to minimize toxicity but maximize tumor-targeting and destruction effects. The other two are VNP20009 and A1-R. We compare the means by which these therapeutic candidate strain models were selected for study, their tumor targeting and tumor destruction phenotypes in vitro and in vivo, and what is currently known about the mechanisms by which they target and destroy tumors.

Predation and selection for antibiotic resistance in natural environments

Genes encoding resistance to antibiotics appear, like the antibiotics themselves, to be ancient, originating long before the rise of the era of anthropogenic antibiotics. However, detailed understanding of the specific biological advantages of antibiotic resistance in natural environments is still lacking, thus limiting our efforts to prevent environmental influx of resistance genes. Here, we propose that antibiotic-resistant cells not only evade predation from antibiotic producers but also take advantage of nutrients released from cells that are killed by the antibiotic-producing bacteria. Thus, predation is potentially an important mechanism for driving antibiotic resistance during slow or stationary phase of growth when nutrients are deprived. This adds to explain the ancient nature and widespread occurrence of antibiotic resistance in natural environments unaffected by anthropogenic antibiotics. In particular, we suggest that nutrient-poor environments including indoor environments, for example, clean rooms and intensive care units may serve as a reservoir and source for antibiotic-producing as well as antibiotic-resistant bacteria.

The role of China in the global spread of the current cholera pandemic

Epidemics and pandemics of cholera, a severe diarrheal disease, have occurred since the early 19th century and waves of epidemic disease continue today. Cholera epidemics are caused by individual, genetically monomorphic lineages of Vibrio cholerae: the ongoing seventh pandemic, which has spread globally since 1961, is associated with lineage L2 of biotype El Tor. Previous genomic studies of the epidemiology of the seventh pandemic identified three successive sub-lineages within L2, designated waves 1 to 3, which spread globally from the Bay of Bengal on multiple occasions. However, these studies did not include samples from China, which also experienced multiple epidemics of cholera in recent decades. We sequenced the genomes of 71 strains isolated in China between 1961 and 2010, as well as eight from other sources, and compared them with 181 published genomes. The results indicated that outbreaks in China between 1960 and 1990 were associated with wave 1 whereas later outbreaks were associated with wave 2. However, the previously defined waves overlapped temporally, and are an inadequate representation of the shape of the global genealogy. We therefore suggest replacing them by a series of tightly delineated clades. Between 1960 and 1990 multiple such clades were imported into China, underwent further microevolution there and then spread to other countries. China was thus both a sink and source during the pandemic spread of V. cholerae, and needs to be included in reconstructions of the global patterns of spread of cholera.

Cholera is a life-threatening, diarrheal disease caused by the bacterium Vibrio cholerae. After a long interregnum of decades without epidemics, the seventh cholera pandemic spread globally since 1961, causing considerable morbidity and mortality. Our analysis of published and newly sequenced genomes provides details on genetic groupings within V. cholerae, so-called clades, that have developed during the recent pandemic spread of these bacteria, and, in some cases, persisted to modern times. We reconstructed some of the pathways taken by the current pandemic since its origins in Indonesia, and show that both South Asia and East Asia are important pathogenic reservoirs and sources of international transmissions.

Laboratory divergence of Methylobacterium extorquens AM1 through unintended domestication and past selection for antibiotic resistance

Background

A common assumption of microorganisms is that laboratory stocks will remain genetically and phenotypically constant over time, and across laboratories. It is becoming increasingly clear, however, that mutations can ruin strain integrity and drive the divergence or “domestication” of stocks. Since its discovery in 1960, a stock of Methylobacterium extorquens AM1 (“AM1”) has remained in the lab, propagated across numerous growth and storage conditions, researchers, and facilities. To explore the extent to which this lineage has diverged, we compared our own “Modern” stock of AM1 to a sample archived at a culture stock center shortly after the strain’s discovery. Stored as a lyophilized sample, we hypothesized that this Archival strain would better reflect the first-ever isolate of AM1 and reveal ways in which our Modern stock has changed through laboratory domestication or other means.

Results

Using whole-genome re-sequencing, we identified some 29 mutations – including single nucleotide polymorphisms, small indels, the insertion of mobile elements, and the loss of roughly 36 kb of DNA - that arose in the laboratory-maintained Modern lineage. Contrary to our expectations, Modern was both slower and less fit than Archival across a variety of growth substrates, and showed no improvement during long-term growth and storage. Modern did, however, outperform Archival during growth on nutrient broth, and in resistance to rifamycin, which was selected for by researchers in the 1980s. Recapitulating selection for rifamycin resistance in replicate Archival populations showed that mutations to RNA polymerase B (rpoB) substantially decrease growth in the absence of antibiotic, offering an explanation for slower growth in Modern stocks. Given the large number of genomic changes arising from domestication (28), it is somewhat surprising that the single other mutation attributed to purposeful laboratory selection accounts for much of the phenotypic divergence between strains.

Conclusions

These results highlight the surprising degree to which AM1 has diverged through a combination of unintended laboratory domestication and purposeful selection for rifamycin resistance. Instances of strain divergence are important, not only to ensure consistency of experimental results, but also to explore how microbes in the lab diverge from one another and from their wild counterparts.

Genetic characterization of atypical Citrobacter freundii

The ability of a bacterial population to survive in different niches, as well as in stressful and rapidly changing environmental conditions, depends greatly on its genetic content. To survive such fluctuating conditions, bacteria have evolved different mechanisms to modulate phenotypic variations and related strategies to produce high levels of genetic diversity. Laboratories working in microbiological diagnosis have shown that Citrobacter freundii is very versatile in its colony morphology, as well as in its biochemical, antigenic and pathogenic behaviours. This phenotypic versatility has made C. freundii difficult to identify and it is frequently confused with both Salmonella enterica and Escherichia coli. In order to determine the genomic events and to explain the mechanisms involved in this plasticity, six C. freundii isolates were selected from a phenotypic variation study. An I-CeuI genomic cleavage map was created and eight housekeeping genes, including 16S rRNA, were sequenced. In general, the results showed a range of both phenotypes and genotypes among the isolates with some revealing a greater similarity to C. freundii and some to S. enterica, while others were identified as phenotypic and genotypic intermediary states between the two species. The occurrence of these events in natural populations may have important implications for genomic diversification in bacterial evolution, especially when considering bacterial species boundaries. In addition, such events may have a profound impact on medical science in terms of treatment, course and outcomes of infectious diseases, evading the immune response, and understanding host-pathogen interactions.

Insights from genomic comparisons of genetically monomorphic bacterial pathogens

Some of the most deadly bacterial diseases, including leprosy, anthrax and plague, are caused by bacterial lineages with extremely low levels of genetic diversity, the so-called 'genetically monomorphic bacteria'. It has only become possible to analyse the population genetics of such bacteria since the recent advent of high-throughput comparative genomics. The genomes of genetically monomorphic lineages contain very few polymorphic sites, which often reflect unambiguous clonal genealogies. Some genetically monomorphic lineages have evolved in the last decades, e.g. antibiotic-resistant Staphylococcus aureus, whereas others have evolved over several millennia, e.g. the cause of plague, Yersinia pestis. Based on recent results, it is now possible to reconstruct the sources and the history of pandemic waves of plague by a combined analysis of phylogeographic signals in Y. pestis plus polymorphisms found in ancient DNA. Different from historical accounts based exclusively on human disease, Y. pestis evolved in China, or the vicinity, and has spread globally on multiple occasions. These routes of transmission can be reconstructed from the genealogy, most precisely for the most recent pandemic that was spread from Hong Kong in multiple independent waves in 1894.