[an error occurred while processing this directive]

Hereditas(Beijing) ›› 2026, Vol. 48 ›› Issue (8): 828-841.doi: 10.16288/j.yczz.25-318

• Research Article • Previous Articles     Next Articles

Genomic diversity and selection pressure analyses of low-virulence Yersinia pestis lineages

Jiayi Zhang1,2(), Yihua Jia2,3, Kai Mu2, Mengnan Cui2, Haifeng Pan1, Yujun Cui1,2(), Yarong Wu2()   

  1. 1 School of Public Health, Anhui Medical University, Hefei 230032, China
    2 State Key Laboratory of Pathogen and Biosecurity, Academy of Military Medical Sciences, Beijing 100071, China
    3 School of Pharmacy, Suzhou Medical College, Soochow University, Suzhou 215000, China
  • Received:2025-12-03 Revised:2026-02-06 Online:2026-08-20 Published:2026-02-27
  • Contact: Yujun Cui, Yarong Wu E-mail:2345010442@stu.ahmu.edu.cn;cuiyujun.new@gmail.com;wuyarong525@126.com
  • Supported by:
    National Key Research and Development Program of China(2024YFC2310100);National Natural Science Foundation of China(32500006)

Abstract:

Yersinia pestis, the causative agent of plague, is a highly virulent pathogen that has caused three historical pandemics and currently persists in multiple active natural foci worldwide, posing a serious threat to human health and public safety. Distinct phylogenetic lineages of Y. pestis exhibit differences in virulence, among which the Microtus biovar within the 0.PE4 lineage is nearly avirulent to large mammals. Comparative genomics between low-virulence and highly virulent populations is therefore key to elucidating the evolutionary basis of virulence. However, comprehensive population-genomic analysis of the low-virulence lineage remains limited, particularly those integrating multiple types of genomic variations. In this study, we analyzed 169 low-virulence strains and 215 representative high-virulence strains collected from long-term national surveillance data in China and public databases. Phylogenetic and spatiotemporal analyses revealed geographic clustering within the low-virulence 0.PE4 lineage. The 0.PE4.3 subclade circulating in China and Mongolia was further subdivided into three regionally associated tertiary lineages, with two more finely resolved quaternary lineages identified specifically in Inner Mongolia. Comparative genomic analyses with representative high-virulence genomes identified mutations fixed in the most recent common ancestor (MRCA) of the low-virulence population and its subclades. Among these, 81 single nucleotide polymorphisms (SNPs), 19 insertion and deletions (Indels), and 5 large fragment losses are shared by all low-virulence strains, defining the genomic features of their MRCA. Together with selection pressure analysis, we identified five genes (ail, rovA, tssH, cdiA, and alr) under strong positive selection, which are involved in virulence, metabolism, and adaptation. Collectively, this study reconstructs the phylogenetic topology of the global low-virulence Y. pestis population and identifies key genomic variations that occurred during its evolutionary process, providing valuable insights for fine-scale tracing and identifying potential molecular targets for elucidating virulence mechanisms.

Key words: Yersinia pestis, virulence, low-virulence lineage, phylogenetic analysis, comparative genomics, selection pressure