Regulatory Functions of PurR in Yersinia pestis: Orchestrating Diverse Biological Activities.

The bacterium Yersinia pestis has developed various strategies to sense and respond to the complex stresses encountered during its transmission and pathogenic processes. PurR is a common transcriptional regulator of purine biosynthesis among microorganisms, and it modulates the transcription level of the pur operon to suppress the production of hypoxanthine nucleotide (IMP). This study aims to understand the functions and regulatory mechanisms of purR in Y. pestis. Firstly, we constructed a purR knockout mutant of Y. pestis strain 201 and compared certain phenotypes of the null mutant (201-ΔpurR) and the wild-type strain (201-WT). The results show that deleting purR has no significant impact on the biofilm formation, growth rate, or viability of Y. pestis under different stress conditions (heat and cold shock, high salinity, and hyperosmotic pressure). Although the cytotoxicity of the purR knockout mutant on HeLa and 293 cells is reduced, the animal-challenging test found no difference of the virulence in mice between 201-ΔpurR and 201-WT. Furthermore, RNA-seq and EMSA analyses demonstrate that PurR binds to the promoter regions of at least 15 genes in Y. pestis strain 201, primarily involved in purine biosynthesis, along with others not previously observed in other bacteria. Additionally, RNA-seq results suggest the presence of 11 potential operons, including a newly identified co-transcriptional T6SS cluster. Thus, aside from its role as a regulator of purine biosynthesis, purR in Y. pestis may have additional regulatory functions.

mStrain: strain-level identification of Yersinia pestis using metagenomic data.

Motivation: High-resolution target pathogen detection using metagenomic sequencing data represents a major challenge due to the low concentration of target pathogens in samples. We introduced mStrain, a novel Yesinia pestis strain/lineage-level identification tool that utilizes metagenomic data. mStrain successfully identified Y. pestis at the strain/lineage level by extracting sufficient information regarding single-nucleotide polymorphisms (SNPs), which can therefore be an effective tool for identification and source tracking of Y. pestis based on metagenomic data during a plague outbreak. Definition: . Strain-level identification: Assigning the reads in the metagenomic sequencing data to an exactly known or most closely representative Y. pestis strain. Lineage-level identification: Assigning the reads in the metagenomic sequencing data to a specific lineage on the phylogenetic tree. canoSNPs: The unique and typical SNPs present in all representative strains. Ancestor/derived state: An SNP is defined as the ancestor state when consistent with the allele of Yersinia pseudotuberculosis strain IP32953; otherwise, the SNP is defined as the derived state. Availability and implementation: The code for running mStrain, the test dataset, and instructions for running the code can be found at the following GitHub repository: https://github.com/xwqian1123/mStrain.

Genomic epidemiological analysis of county-scale Yersinia pestis spread pattern over 50 years in a Southwest Chinese prefecture.

Plague, one of the most devastating infectious diseases in human history, is caused by the bacterium Yersinia pestis. Since the 1950s, the Dehong Dai-Jingpo Autonomous Prefecture (DH) in Yunnan Province, China, has recorded plague outbreaks that have resulted in 1,153 human cases and 379 deaths. The genetic diversity and transmission characteristics of Y. pestis strains in this region remain unknown. Here, we performed high-resolution genomic epidemiological analysis of 175 Y. pestis strains isolated from five counties and 19 towns in DH between 1953 and 2007. Phylogenetic analysis revealed that most DH strains were located in lineage 1.ORI2, which could be further subdivided into seven sub-phylogroups (SPG1-SPG7). The dominant sub-phylogroups of Y. pestis in DH varied during different periods and presented a population shift. Genomic evidence showed that plague might have emerged from the southwest of DH (e.g., Longchuan or Ruili counties) or its bordering countries, and subsequently spread to the northeast in multiple waves between 1982 and 2007. Our study infers a fine-scale phylogeny and spread pattern of the DH Y. pestis population, which extends our knowledge regarding its genetic diversity and provides clues for the future prevention and control of plague in this region.

Genomic diversity of Yersinia pestis from Yunnan Province, China, implies a potential common ancestor as the source of two plague epidemics.

Plague, caused by Yersinia pestis, is a zoonotic disease that can reemerge and cause outbreaks following decades of latency in natural plague foci. However, the genetic diversity and spread pattern of Y. pestis during these epidemic-silent cycles remain unclear. In this study, we analyze 356 Y. pestis genomes isolated between 1952 and 2016 in the Yunnan Rattus tanezumi plague focus, China, covering two epidemic-silent cycles. Through high-resolution genomic epidemiological analysis, we find that 96% of Y. pestis genomes belong to phylogroup 1.ORI2 and are subdivided into two sister clades (Sublineage1 and Sublineage2) characterized by different temporal-spatial distributions and genetic diversity. Most of the Sublineage1 strains are isolated from the first epidemic-silent cycle, while Sublineage2 strains are predominantly from the second cycle and revealing a west to east spread. The two sister clades evolved in parallel from a common ancestor and independently lead to two separate epidemics, confirming that the pathogen responsible for the second epidemic following the silent interval is not a descendant of the causative strain of the first epidemic. Our results provide a mechanism for defining epidemic-silent cycles in natural plague foci, which is valuable in the prevention and control of future plague outbreaks.

Population genomics implies potential public health risk of two non-toxigenic Vibrio cholerae lineages.

Diarrheal cases caused by non-toxigenic Vibrio cholerae have been reported globally. Lineages L3b and L9, characterized as ctxAB-negative and tcpA-positive (CNTP), pose the highest risk and have caused long-term epidemics in different regions worldwide. From 2001 to 2018, two waves (2001-2012 and 2013-2018) of epidemic caused by non-toxigenic V. cholerae occurred in the developed city of Hangzhou, China. In this study, through the integrated analysis of 207 genomes of Hangzhou isolates from these two waves (119 and 88) and 1573 publicly available genomes, we showed that L3b and L9 lineages together caused the second wave as had happened in the first wave, but the dominant lineage shifted from L3b (first wave: 69%) to L9 (second wave: 50%). We further found that the genotype of a key virulence gene, tcpF, in the L9 lineage during the second wave shifted to type I, which may have enhanced bacterial colonization in humans and potentially promoted the pathogenic lineage shift. Moreover, we found that 21% of L3b and L9 isolates had changed to predicted cholera toxin producers, providing evidence that gain of complete CTXφ-carrying ctxAB genes, rather than ctxAB gain in pre-CTXφ-carrying isolates, led to the transition. Taken together, our findings highlight the possible public health risk associated with L3b and L9 lineages due to their potential to cause long-term epidemics and turn into high-virulent cholera toxin producers, which necessitates a more comprehensive and unbiased sampling in further disease control efforts.

Small Insertions and Deletions Drive Genomic Plasticity during Adaptive Evolution of Yersinia pestis.

The life cycle of Yersinia pestis has changed a lot to adapt to flea-borne transmission since it evolved from an enteric pathogen, Yersinia pseudotuberculosis. Small insertions and deletions (indels), especially frameshift mutations, can have major effects on phenotypes and contribute to virulence and host adaptation through gene disruption and inactivation. Here, we analyzed 365 Y. pestis genomes and identified 2,092 genome-wide indels on the core genome. As recently reported in Mycobacterium tuberculosis, we also detected "indel pockets" in Y. pestis, with average complexity scores declining around indel positions, which we speculate might also exist in other prokaryotes. Phylogenic analysis showed that indel-based phylogenic tree could basically reflect the phylogenetic relationships of major phylogroups in Y. pestis, except some inconsistency around the Big Bang polytomy. We observed 83 indels arising in the trunk of the phylogeny, which played a role in accumulation of pseudogenes related to key metabolism and putatively pathogenicity. We also discovered 32 homoplasies at the level of phylogroups and 7 frameshift scars (i.e., disrupted reading frame being rescued by a second frameshift). Additionally, our analysis showed evidence of parallel evolution at the level of genes, with sspA, rpoS, rnd, and YPO0624, having enriched mutations in Brazilian isolates, which might be advantageous for Y. pestis to cope with fluctuating environments. The diversified selection signals observed here demonstrates that indels are important contributors to the adaptive evolution of Y. pestis. Meanwhile, we provide potential targets for further exploration, as some genes/pseudogenes with indels we focus on remain uncharacterized. IMPORTANCE Yersinia pestis, the causative agent of plague, is a highly pathogenic clone of Yersinia pseudotuberculosis. Previous genome-wide SNP analysis provided few adaptive signatures during its evolution. Here by investigating 365 public genomes of Y. pestis, we give a comprehensive overview of general features of genome-wide indels on the core genome and their roles in Y. pestis evolution. Detection of "indel pockets," with average complexity scores declining around indel positions, in both Mycobacterium tuberculosis and Y. pestis, gives us a clue that this phenomenon might appear in other bacterial genomes. Importantly, the identification of four different forms of selection signals in indels would improve our understanding on adaptive evolution of Y. pestis, and provide targets for further physiological mechanism researches of this pathogen. As evolutionary research based on genome-wide indels is still rare in bacteria, our study would be a helpful reference in deciphering the role of indels in other species.