关键词：h9亚型禽流感 中国大陆地区 进化分析 

摘要：目的从分子水平研究中国大陆地区H9亚型禽流感病毒的起源、进化和地理扩散情况,为禽流感的防控提供科学依据。方法从全球公共流感病毒数据库收集中国大陆地区的H9亚型禽流感病毒HA基因序列,通过对下载的序列进行比对和聚类筛选,得到1997—2018年180株含有准确采样信息的序列,使用BEAST等生物信息学软件分析病毒选择压力、构建贝叶斯法时间进化树以及重建病毒进化时空动态。结果序列的核苷酸差异在0.2%~20.7%之间,氨基酸的差异在0.0%~15.8%之间。编码区的核苷酸替代率分别为:0.540、0.607和1.853。非同义突变率与同义突变率的比值的平均值为0.1789,但是在7个位点上比值大于1。HA基因的平均进化替代率为2.476×10^-3替换/位点/年,最早的公共祖先出现在1947.457年。构建的基因进化树的拓扑结构分为多个分支,根节点最有可能的位置是广东省,后验概率为21.7%。基于目前数据的时空动态分析显示,1959—2018年间病毒从广东传播到安徽,再到广西、江苏,再以江苏为中心向周边多个省份传播,最后在多个省份之间出现传播。其中江苏传出的6条路径具有统计学支持。结论H9亚型禽流感病毒出现的时间可能比国内最早报导的早很多,在中国多个省份之间传播扩散情况复杂,江苏省是传播扩散的一个重要中心。

实用预防医学杂志要求:

{1}非直接引文注释，注释前应加“参见”；非引用原始资料时，应先注明原始作品之相关信息，再以“转引自”为引领词注明转引之文献详细信息。

{2}本刊倡导严谨、求实的学风，欢迎有真知灼见、有理有据、论证严密、合乎语法规范的学术论文。

{3}文中图要有图名，图、表、文字切忌相互重复。文中数据一律采用法定计量单位。

{4}注释：主要用于对文章篇名、作者及文内某一专题内容作必要的解释或说明。序号在注释处右上角用阿拉伯数字①②等表示。

{5}正文中的数据组不能同时采用图、表、文字描述，可选用其一，选择次序为：首先用图、次之用表、最后用文字。

注：因版权方要求，不能公开全文，如需全文，请咨询杂志社