鼠疫耶尔森菌低毒种群基因组多样性和选择压力分析
收稿日期: 2025-12-03
修回日期: 2026-02-06
网络出版日期: 2026-02-27
基金资助
国家重点研发计划(2024YFC2310100);国家自然科学基金项目(32500006)
Genomic diversity and selection pressure analyses of low-virulence Yersinia pestis lineages
Received date: 2025-12-03
Revised date: 2026-02-06
Online published: 2026-02-27
Supported by
National Key Research and Development Program of China(2024YFC2310100);National Natural Science Foundation of China(32500006)
鼠疫耶尔森菌(Yersinia pestis)是引发鼠疫的烈性病原体，曾造成三次历史大流行，目前分布于全球多个活跃的自然疫源地，对人类健康与公共安全构成严重威胁。不同型别的鼠疫菌在毒力上存在显著差异，其中 0.PE4 谱系的田鼠型菌株对大型哺乳动物几乎无毒。比较低毒与高毒种群的基因组差异有助于阐明毒力的演化基础，但目前针对低毒种群仍缺乏涵盖多类型变异的系统性群体水平基因组学研究。本研究整合我国历年监测数据及公共数据库中其他国家和地区的菌株基因组，共纳入169株低毒菌株和215株高毒主要谱系代表菌株序列。系统发育与时空分布分析显示，0.PE4低毒种群呈现明显的地理聚集性，其中分布于中国及蒙古国的0.PE4.3亚群进一步分化为3个具有地域特征的三级谱系，并在内蒙古地区形成2个更细分的四级谱系。与高毒主要谱系的比较基因组分析鉴定出低毒种群最近共同祖先及其各亚群分支上固定的变异。其中，81个SNP、19 个Indel和5个大片段缺失为所有低毒菌株所特有，代表了低毒种群最近共同祖先的基因组特征。结合选择压力信号分析，共鉴定出5个受到强选择作用的基因(ail、rovA、tssH、cdiA、alr)，分别与毒力、代谢和适应性相关。本研究重建了全球范围内鼠疫菌低毒种群的精细系统发育拓扑结构，揭示了其进化过程中的关键基因组变异，为鼠疫监测、溯源与毒力演化机制研究提供了重要参考。
张佳怡, 贾艺华, 穆凯, 崔梦楠, 潘海峰, 崔玉军, 武雅蓉 . 鼠疫耶尔森菌低毒种群基因组多样性和选择压力分析[J]. 遗传, 2026 , 48(8) : 828 -841 . DOI: 10.16288/j.yczz.25-318
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.
| [1] | Achtman M, Morelli G, Zhu PX, Wirth T, Diehl I, Kusecek B, Vogler AJ, Wagner DM, Allender CJ, Easterday WR, Chenal-Francisque V, Worsham P, Thomson NR, Parkhill J, Lindler LE, Carniel E, Keim P. Microevolution and history of the plague bacillus, Yersinia pestis. Proc Natl Acad Sci USA, 2004, 101(51): 17837-17842. |
| [2] | Perry RD, Fetherston JD. Yersinia pestis--etiologic agent of plague. Clin Microbiol Rev, 1997, 10(1): 35-66. |
| [3] | Stenseth NC, Atshabar BB, Begon M, Belmain SR, Bertherat E, Carniel E, Gage KL, Leirs H, Rahalison L. Plague: past, present, and future. PLoS Med, 2008, 5(1): e3. |
| [4] | World Health Organization. Plague fact sheet. 2022. https://www.who.int/news-room/fact-sheets/detail/plague. |
| [5] | World Health Organization. Pathogens prioritization: a scientific framework for epidemic and pandemic research preparedness. 2024, 38. |
| [6] | Inglesby TV, Dennis DT, Henderson DA, Bartlett JG, Ascher MS, Eitzen E, Fine AD, Friedlander AM, Hauer J, Koerner JF, Layton M, McDade J, Osterholm MT, O'Toole T, Parker G, Perl TM, Russell PK, Schoch-Spana M, Tonat K. Plague as a biological weapon: medical and public health management. Working group on civilian biodefense. JAMA, 2000, 283(17): 2281-2290. |
| [7] | Yang R, Atkinson S, Chen Z, Cui Y, Du Z, Han Y, Sebbane F, Slavin P, Song Y, Yan Y, Wu Y, Xu L, Zhang C, Zhang Y, Hinnebusch BJ, Stenseth NC, Motin VL. Yersinia pestis and plague: some knowns and unknowns. Zoonoses. 2023, 3(1): 5. |
| [8] | Li YJ, Cui YJ, Hauck Y, Platonov ME, Dai E, Song YJ, Guo ZB, Pourcel C, Dentovskaya SV, Anisimov AP, Yang RF, Vergnaud G. Genotyping and phylogenetic analysis of Yersinia pestis by MLVA: insights into the worldwide expansion of Central Asia plague foci. PLoS One, 2009, 4(6): e6000. |
| [9] | Zhou DS, Tong ZZ, Song YJ, Zhang L, Pei DC, Pang X, Han YP, Li M, Cui BZ, Wang J, Guo ZB, Qi ZZ, Jin LX, Zhai JH, Du ZM, Wang XY, Wang J, Huang PT, Yang RF. Molecular mechanism of biovar conversion of Yersinia pestis and the proposition of a new biovar: microtus. PLA Medical Journal, 2004, 29(3): 211-215. |
| 周冬生, 童宗中, 宋亚军, 张玲, 裴德翠, 庞昕, 韩延平, 李敏, 崔百忠, 王津, 郭兆彪, 祁芝珍, 金丽霞, 翟俊辉, 杜宗敏, 王效义, 汪建, 黄培堂, 杨瑞馥. 鼠疫耶尔森菌生物型变异遗传基础的研究和新生物型--田鼠型的提出. 解放军医学杂志, 2004, 29(3): 211-215. | |
| [10] | Zhang XA, Hai R, Wei JC, Cui ZG, Zhang EM, Song ZZ, Yu DZ. MLVA distribution characteristics of Yersinia pestis in China and the correlation analysis. BMC Microbiol, 2009, 9: 205. |
| [11] | Zhou DS, Tong ZZ, Song YJ, Han YP, Pei DC, Pang X, Zhai JH, Li M, Cui BZ, Qi ZZ, Jin LX, Dai RX, Du ZM, Wang J, Guo ZB, Wang J, Huang PT, Yang RF. Genetics of metabolic variations between Yersinia pestis biovars and the proposal of a new biovar, microtus. J Bacteriol, 2004, 186(15): 5147-5152. |
| [12] | Cui YJ, Yu C, Yan YF, Li DF, Li YJ, Jombart T, Weinert LA, Wang ZY, Guo ZB, Xu LZ, Zhang YJ, Zheng HC, Qin N, Xiao X, Wu MS, Wang XY, Zhou DS, Qi ZZ, Du ZM, Wu HL, Yang XW, Cao HZ, Wang H, Wang J, Yao SS, Rakin A, Li YR, Falush D, Balloux F, Achtman M, Song YJ, Wang J, Yang RF. Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis. Proc Natl Acad Sci USA, 2013, 110(2): 577-582. |
| [13] | Wu YR, Xin YQ, Yang XY, Song K, Zhang QW, Zhao HH, Li CX, Jin Y, Guo Y, Tan YF, Song YJ, Tian HY, Qi ZZ, Yang RF, Cui YJ. Hotspots of genetic change in Yersinia pestis. Nat Commun, 2025, 16(1): 388. |
| [14] | Zheng HY, Yan L, Yang C, Wu YR, Qin JL, Hao TY, Yang DJ, Guo YC, Pei XY, Zhao TY, Cui YJ. Population genomics study of Vibrio alginolyticus. Hereditas(Beijing), 2021, 43(4): 350-361. |
| 郑宏源, 闫琳, 杨超, 武雅蓉, 秦婧靓, 郝彤宇, 杨大进, 郭云昌, 裴晓燕, 赵彤言, 崔玉军. 溶藻弧菌群体基因组学研究. 遗传, 2021, 43(4): 350-361. | |
| [15] | 樊振亚, 罗运珩, 王身荣, 金玲玲, 周祥, 刘俊, 张耀星, 李凤. 布氏田鼠鼠疫不危害人类. 中国地方病防治杂志, 1995, 10(1): 56-57. |
| [16] | Zhang QW, Wang Q, Tian G, Qi ZZ, Zhang XC, Wu XH, Qiu YF, Bi YJ, Yang XY, Xin YQ, He J, Zhou JY, Zeng L, Yang RF, Wang XY. Yersinia pestis biovar Microtus strain 201, an avirulent strain to humans, provides protection against bubonic plague in rhesus macaques. Hum Vaccin Immunother, 2014, 10(2): 368-377. |
| [17] | Chen H, Yang W, Li FZ, Jiang HZ, Xu GR, Li GQ, Yan DL, Luo ZBD, Duan YJ, Deng JY, Zhu XP, Wang LM. The study of harm to Yersinia pestis of Microtus fuscus in Qinghai-Tibet Plateau. Chin J Ctrl Endem Dis, 2006, 21(6): 324-326. |
| 陈虹, 杨文, 李富忠, 蒋和柱, 许光荣, 李光清, 严冬丽, 罗志丹巴, 段勇军, 邓佳云, 祝小平, 汪立茂. 青藏高原青海田鼠型鼠疫菌危害性研究. 中国地方病防治杂志, 2006, 21(6): 324-326. | |
| [18] | 樊振亚, 罗运珩, 李凤, 张春华, 苏晓仙, 王身荣. 鼠疫耶尔森菌锡林郭勒高原型对人致病性试验. 中国地方病防治杂志, 1994, 9(6): 340-342. |
| [19] | Song YJ, Tong ZZ, Wang J, Wang L, Guo ZB, Han YP, Zhang JG, Pei DC, Zhou DS, Qin HO, Pang X, Han YJ, Zhai JH, Li M, Cui BZ, Qi ZZ, Jin LX, Dai RX, Chen F, Li ST, Ye C, Du ZM, Lin W, Wang J, Yu J, Yang HM, Wang J, Huang PT, Yang RF. Complete genome sequence of Yersinia pestis strain 91001, an isolate avirulent to humans. DNA Res, 2004, 11(3): 179-197. |
| [20] | Li YJ, Dai EH, Cui YJ, Li M, Zhang YJ, Wu MS, Zhou DS, Guo ZB, Dai X, Cui BZ, Qi ZZ, Wang ZY, Wang H, Dong XQ, Song ZZ, Zhai JH, Song YJ, Yang RF. Different region analysis for genotyping Yersinia pestis isolates from China. PLoS One, 2008, 3(5): e2166. |
| [21] | Wu YR, Yang C, Mu K, Guo Y, Song YJ, Yang RF, Cui YJ. Insights into Yersinia pestis evolution through rearrangement analysis of 242 complete genomes. Nat Genet, 2025, 57(8): 1994-2003. |
| [22] | Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30(15): 2114-2120. |
| [23] | Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, Lesin VM, Nikolenko SI, Pham S, Prjibelski AD, Pyshkin AV, Sirotkin AV, Vyahhi N, Tesler G, Alekseyev MA, Pevzner PA. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol, 2012, 19(5): 455-477. |
| [24] | Seemann T. Shovill: faster SPAdes assembly of Illumina reads. 2017. https://github.com/tseemann/shovill. |
| [25] | Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics, 2014, 30(14): 2068-2069. |
| [26] | Delcher AL, Salzberg SL, Phillippy AM. Using MUMmer to identify similar regions in large sequence sets. Curr Protoc Bioinformatics, 2003, Chapter 10: Unit 10.3. |
| [27] | Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv: 1303.3997, 2013, 1-3. |
| [28] | Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, Jordan T, Shakir K, Roazen D, Thibault J, Banks E, Garimella KV, Altshuler D, Gabriel S, DePristo MA. From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics, 2013, 43(1110): 11.10.1-11.10.33. |
| [29] | Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res, 1999, 27(2): 573-580. |
| [30] | Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol, 1990, 215(3): 403-410. |
| [31] | Nguyen LT, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol, 2015, 32(1): 268-274. |
| [32] | Rambaut A. FigTree - v1.4.4. 2016. http://tree.bio.ed.ac.uk/software/figtree/. |
| [33] | Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5:an online tool for phylogenetic tree display and annotation. Nucleic Acids Res, 2021, 49(W1): W293-W296. |
| [34] | Tonkin-Hill G, MacAlasdair N, Ruis C, Weimann A, Horesh G, Lees JA, Gladstone RA, Lo S, Beaudoin C, Floto RA, Frost SDW, Corander J, Bentley SD, Parkhill J. Producing polished prokaryotic pangenomes with the Panaroo pipeline. Genome Biol, 2020, 21(1): 180. |
| [35] | Li RQ, Li YR, Kristiansen K, Wang J. SOAP: short oligonucleotide alignment program. Bioinformatics, 2008, 24(5): 713-714. |
| [36] | O'Leary NA, Wright MW, Brister JR, Ciufo S, Haddad D, McVeigh R, Rajput B, Robbertse B, Smith-White B, Ako-Adjei D, Astashyn A, Badretdin A, Bao YM, Blinkova O, Brover V, Chetvernin V, Choi J, Cox E, Ermolaeva O, Farrell CM, Goldfarb T, Gupta T, Haft D, Hatcher E, Hlavina W, Joardar VS, Kodali VK, Li WJ, Maglott D, Masterson P, McGarvey KM, Murphy MR, O'Neill K, Pujar S, Rangwala SH, Rausch D, Riddick LD, Schoch C, Shkeda A, Storz SS, Sun HZ, Thibaud-Nissen F, Tolstoy I, Tully RE, Vatsan AR, Wallin C, Webb D, Wu W, Landrum MJ, Kimchi A, Tatusova T, DiCuccio M, Kitts P, Murphy TD, Pruitt KD. Reference sequence (RefSeq) database at NCBI: current status, taxonomic expansion, and functional annotation. Nucleic Acids Res, 2016, 44(D1): D733-745. |
| [37] | Cantalapiedra CP, Hernández-Plaza A, Letunic I, Bork P, Huerta-Cepas J. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic Scale. Mol Biol Evol, 2021, 38(12): 5825-5829. |
| [38] | Crispell J, Balaz D, Gordon SV. HomoplasyFinder: a simple tool to identify homoplasies on a phylogeny. Microb Genom, 2019, 5(1): e000245. |
| [39] | Wang DP, Zhang YB, Zhang Z, Zhu J, Yu J. KaKs_ Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics, 2010, 8(1): 77-80. |
| [40] | Zhang ZM, Zhou H, Huang XZ, Zhang DY, Zhang JY, Lin Y, Fang LW, Zhang XC, Cui YJ, Wu YR, Li YJ. Study on population genomics of Bacillus anthracis based on multiple types of genetic variations. Hereditas(Beijing), 2025, 47(6): 681-693. |
| 张祖铭, 周豪, 黄学治, 张多悦, 张佳怡, 林昱, 方立崴, 张秀昌, 崔玉军, 武雅蓉, 李艳君. 基于多类型变异的炭疽芽孢杆菌群体基因组学研究. 遗传, 2025, 47(6): 681-693. | |
| [41] | Wu YR, Hao TY, Qian XW, Zhang XL, Song YJ, Yang RF, Cui YJ. Small insertions and deletions drive genomic plasticity during adaptive evolution of Yersinia pestis. Microbiol Spectr, 2022, 10(3): e0224221. |
| [42] | Craig NL. Transposon Tn7. Curr Top Microbiol Immunol, 1996, 204: 27-48. |
| [43] | Young R. Phage lysis: three steps, three choices, one outcome. J Microbiol, 2014, 52(3): 243-258. |
| [44] | Walsh CT. Enzymes in the D-alanine branch of bacterial cell wall peptidoglycan assembly. J Biol Chem, 1989, 264(5): 2393-2396. |
| [45] | Gallegos MT, Schleif R, Bairoch A, Hofmann K, Ramos JL. Arac/XylS family of transcriptional regulators. Microbiol Mol Biol Rev, 1997, 61(4): 393-410. |
| [46] | Felek S, Krukonis ES. The Yersinia pestis Ail protein mediates binding and Yop delivery to host cells required for plague virulence. Infect Immun, 2009, 77(2): 825-836. |
| [47] | Tsang TM, Wiese JS, Felek S, Kronshage M, Krukonis ES. Ail proteins of Yersinia pestis and Y. pseudotuberculosis have different cell binding and invasion activities. PLoS One, 2013, 8(12): e83621. |
| [48] | Aoki SK, Pamma R, Hernday AD, Bickham JE, Braaten BA, Low DA. Contact-dependent inhibition of growth in Escherichia coli. Science, 2005, 309(5738): 1245-1248. |
| [49] | Cathelyn JS, Crosby SD, Lathem WW, Goldman WE, Miller VL. RovA, a global regulator of Yersinia pestis, specifically required for bubonic plague. Proc Natl Acad Sci USA, 2006, 103(36): 13514-13519. |
| [50] | Basler M, Pilhofer M, Henderson GP, Jensen GJ, Mekalanos JJ. Type VI secretion requires a dynamic contractile phage tail-like structure. Nature, 2012, 483(7388): 182-186. |
| [51] | Zhao XN, Cui YJ, Yan YF, Du ZM, Tan YF, Yang HY, Bi YJ, Zhang PP, Zhou L, Zhou DS, Han YP, Song YJ, Wang XY, Yang RF. Outer membrane proteins ail and OmpF of Yersinia pestis are involved in the adsorption of T7-related bacteriophage Yep-phi. J Virol, 2013, 87(22): 12260-12269. |
/
| 〈 |
|
〉 |
www.chinagene.cn
备案号：京ICP备09063187号-4
总访问:,今日访问:,当前在线:
