自2013年3月,我国长江三角洲首次发现人感染H7N9禽流感病例,当月福建省也报告了本地首例人感染H7N9禽流感病例,截至目前疫情经历了多个流行阶段,没有发现明显停止迹象[1]。本研究于2014年11月 — 2017年10月在福建省连续收集75个病例,探测福建省人感染H7N9禽流感发病的可能危险因素,以期为制定和调整防控策略提供科学指导。
1 对象与方法 1.1 对象(1)病例组 连续地收集2014年11月 — 2017年10月,现住址为福建省的人感染H7N9禽流感实验室确诊病例,含外地报告本地,不含港澳台,不含本地报告外地病例,病例中明确为活禽从业人员不予纳入病例组。(2)对照组 运用成组设计的病例对照研究,采用随机数字法选取与病例在同一小区或自然村居住时间 ≥ 6个月;年龄差异在 ± 3岁以内的社区非病例为对照。(3)样本量估计 参照外省同类研究[2],对照组暴露率26.7 %,OR值3.4,双侧α = 0.05,β = 0.10,运用成组设计样本量计算公式得两组分别需要60个样本,测算公式:
先征得受访对象知情同意,再采用自制的问卷开展调查,主要收集人口学特征、卫生行为习惯、患慢性病情况和禽类暴露情况等。洗手情况判定:经常洗手(≥ 5次/d)、偶尔洗手(3~4次/d)、很少洗手(1~2次/d)。吸烟情况:吸烟者根据世界卫生组织(World Health Organization,WHO)定义[3]指一生中连续或累积吸烟 ≥ 6个月者;现在吸烟[3]指符合“吸烟者”定义,且在发病前30 d内吸过烟的人;现已戒烟指每天连续吸烟至少6个月,调查时已不再吸烟 > 30 d;从不吸烟是指未达到WHO定义吸烟者以外的人员。
1.3 统计分析采用Epi Data 3.1软件双录入,用SAS 9.2软件开展相应的统计分析及假设检验。将单因素分析筛选出的危险因素进行无序多因素logistic回归分析,拟合回归方程。检验水准取双侧α = 0.05,P < 0.05,回归模型因素筛选α 出 = α入 = 0.05,P < 0.05,拟合优度判断采用一致性指数法(index of concordance,C)。
2 结 果 2.1 基本情况(表1)2014年11月 — 2017年10月,成功收集75个病例,300个对照的流行病学资料,覆盖所有设区市和平潭综合实验区,占同期86例病例的87.21 %。病例与对照组在年龄中位数、年龄组构成、性别构成、城乡构成和体质指数(body mass index,BMI)构成方面均无统计学差异。在是否本人回答问卷方面病例组要低于对照(病例组49.3 %,对照组88.00 %,χ2 = 56.64,P < 0.01),但病例组中非本人回答和本人回答的问卷在患慢性病、市场接触情况和家养禽类等信息无统计学差异( χ2 = 1.60,P = 0.21;χ2 = 1.25,P = 0.26;χ2 = 0.35,P = 0.55)。
2.2 单因素分析 2.2.1 基础疾病情况(表2)BMI异常(BMI ≥ 28或 < 18.5)、患慢性病、患肺部疾病、糖尿病、慢性肝脏疾病、慢性支气管炎病例组高于对照,有统计学意义,可能是发病的危险因素。患高血压病2组没有统计学意义。
2.2.2 市场暴露情况(表3)暴露活禽市场(含停留、未停留和农贸市场活禽摊点)提示可能是危险因素。暴露活禽摊点为开放性、设宰杀点、设置机械脱毛和正在机械脱毛提示可能是危险因素,而半密闭活禽摊点暴露2组无统计学差异(P = 1.00)。在活禽市场摊点直接禽类或摊点物品也可能是危险因素。暴露接触统一屠宰的鲜禽(白条禽)2组之间没有统计学差异,提示可能不是危险因素。
2.2.3 家庭禽类暴露情况(表4)暴露家养禽、暴露家养鸡、鸭提示可能是危险因素。直接接触家养禽和其中的捕捉家禽、给家禽喂药,非直接接触家养禽和其中给家禽投食、触摸家禽笼舍、打扫清理家禽禽舍提示可能是危险因素。而家中出现病死禽、在家宰杀市场买回的活禽、围观宰杀过程、家中购买摊点现杀的活禽和接触接触家中买回的鲜杀活禽、接触家中冰冻禽两组间没有统计学差异。家中购买生鲜禽(白条禽)、接触过购买的生鲜禽(白条禽)和家中购买过冰冻禽病例组比例要明显低于对照组。
2.2.4 其他禽类暴露和卫生与行为习惯(表5)暴露于野鸟、野鸟排泄物、暴露禽类职业人员可能是危险因素。经常洗手、常用肥皂/洗手液洗手可能是疾病的保护因素。吸烟率对照高于病例组(χ2 = 12.27,P < 0.01)。暴露活禽加工餐馆、访问有家养禽家庭和接触过发热病人2组间没有统计学差异。
2.3 多因素逐步logistic回归分析(表6)根据单因素分析结果,以人感染H7N9禽流感发病作为因变量,将单因素分析有统计学意义的变量作为自变量。采用多因素逐步logistic回归分析,模型构建(α出 = α入 = 0.05,P < 0.05)。其中8个因素仍有意义纳入回归模型,分别为:患肺部疾病(X 1)、糖尿病(X2)、暴露活禽市场(X3)、暴露使用机械脱毛摊点(X4)、经常洗手(X5)、触摸家禽笼舍(X6)、暴露野鸟排泄物(X7)和现在吸烟(X8)。通过模型拟合得到回归方程:
本研究以福建省人感染H7N9禽流感确诊病例为病例组,以社区人口为对照研究发病危险因素,具有地域特色,病例和对照组在性别、年龄组构成、城乡构成和BMI构成方面都均衡可比。研究对象涵盖所有设区市和平潭综合实验区,病例组占同期报告确诊病例数的87.21 %,具有一定的代表性。
逐步多因素logistic回归分析得到,暴露活禽市场、暴露使用机械脱毛摊点都是人感染H7N9禽流感发病的独立危险因素,暴露统一屠宰的鲜禽(白条禽)市场暴露没有提示风险,与国内同类研究结果一致[4]。因此活禽市场仍是感染H7N9禽流感重要来源[2, 4 – 7],活禽市场在疾病传播过程中起到了病毒汇集、增殖、储存、杂交、变异和传播放大等作用[8 – 9]。主要依据如下:1是活禽市场禽类来源渠道复杂、多种活禽混合在同一市场销售容易造成交叉污染,活禽市场外环境和禽类标本不仅检出H7N9还存在H5N1、H9N2等多种禽类病毒亚型[10 – 11],病毒间存在杂交和变异可能。2是H7N9病毒已出现珠江三角和长三角病毒2个分支[12],说明病毒不断与当地其他禽类感病毒亚型重配。3是广东地区已在活禽市场或病例中分离到突变的高致病H7N9禽流感病毒[13 – 14],部分病毒耐药位点也出现突变提示对现有的抗病毒治疗药物出现耐药可能[14 – 15]。4是人间H7N9疫情暴发存在接触活禽市场时间依赖关系[16],季节性干预策略可以大大降低疫情暴发的规模。5是大部分病例有活禽市场暴露史[2, 5, 7 – 8, 16],从病例和活禽市场标本分离的病毒存在高度同源[17 – 18]。本研究发现的暴露使用机械脱毛摊点除接触传播外可能与机械脱毛比传统手工脱毛方式更容易产生气溶胶有一定关系,美国疾病预防控制中心(Centers for Disease Control and Prevention,CDC)和埃默里大学模拟体外的自然感染,验证H7N9(LPAI)在人类支气管上皮细胞(Calu-3)中都具有高度传染性[19]。因此,活禽市场仍是人感染H7N9禽流感重要来源,努力减少活禽市场接触对于防止长期流行是至关重要,规范管理和关闭活禽市场是控制疫情的有效措施[16, 20],但考虑到人民群众消费习惯短时间内较难改变,因此加强市场管控是当前重要的替代手段,保持活禽市场清洁,落实零存栏,加强活禽摊点建设做好脱毛机与销售窗口物理隔离,加强禽类摊点从业人员防护措施。并加大宣传教育力度,逐步推进禽类统一宰杀白条上市,减少活禽暴露等对防控人感染H7N9禽流感非常必要。
暴露触摸家禽笼舍是福建省人感染H7N9禽流感发病的独立危险因素。本研究发现福建省家庭禽类养殖达23.67 %,高于国内其他研究水平(14.12 %)[4],触摸家禽笼舍可能与直接到携带病毒的禽类和受病毒污染的环境有关,接触禽笼舍与清理打扫卫生容易形成扬尘和病毒气溶胶等增加感染发病风险。本研究还发现暴露野鸟排泄物是可能感染的危险因素与国内研究不同[4, 21],该因素可能具有福建省地域特点。1是福建省暴露于野鸟等野禽的比例为31.33 % 要高于国内其他省份(14.71 %)[4],可能与福建省山区多森林覆盖率高,气候适宜野鸟生存[22],增加暴露概率。2 是H7N9禽流感病毒最早就起源于野禽再传播给家禽[23]。福建省广大农村地区家禽养殖较高,且存在散养现象,增加野禽和家禽混合概率,可能会继续演绎野禽病毒感染到家禽的可能性,应进一步加强散养禽类管理,加强病原学监测等工作。暴露禽类职业人员也提示是危险因素,国内研究发现禽类从业人员H7N9禽流感抗体滴度要高于普通人群,提示禽类从业人员是高危人群,应采取必要的防护措施[4]。
患肺部疾病、患糖尿病是人感染H7N9禽流感发病的危险因素,与国内研究患慢性病中的肺部疾病、糖尿病和肝脏疾病可能会使机体更容易感染流感发病风险的结果相似[4, 24],因此要做好患高危慢性病等重点人群防控工作。患高血压病并不会增加发病风险,与国内同类研究相同[2, 4 – 5, 21]。经常洗手(OR = 0.07,95 % CI = 0.04~0.16)是人感染H7N9禽流感的保护因素,国内同类研究也得出相似结果,认为肥皂或洗手液能够使H5N1流感病毒失活[4]。但吸烟降发病风险也可能是研究中出现的混杂偏倚,主要分析如下:(1)国内相关研究尚未很好证实吸烟在人感染H7N9禽流感发病中所起的保护作用和保护机制;(2)国内同类研究发现调查过程中病人可能隐瞒某些不良卫生习惯,如吸烟等[4]。因此,正在吸烟是否为人感染H7N9禽流感发病的保护因素,研究不够充分,机制尚不明确,仍需要进一步研究。
本研究具有一定的局限性:(1)病例对照研究中纳入75例病例,其他在本研究开始之前的病例和部分无法配合调查的病例未纳入,可能存在选择偏倚,研究结果外推时受限。(2)本研究病例组与对照组在年龄、性别和现住址等因素上均无统计学差异,而前期研究年龄、性别和农村可能是发病的影响因素,因此在危险因素筛选中可能将上述因素排除在外而影响研究结果。(3)相关信息收集都在发病之后的一段时间进行,可能存在回忆偏倚。(4)是病例组中本人回答率和对照组存在差异。虽然通过分析病例组中非本人回答和本人回答的问卷在患慢性病、市场接触情况和家养禽类等信息没有统计学差异。但非本人回答的信息可能会存在其他细节方面的信息偏倚。(5)病例和对照样本量还是偏小,部分人群中频率低的因素,受到统计检验效能影响可能没能纳入模型,制约了关联性分析。(6)能存在混在偏倚,本研究采用成组设计的病例对照研究,尽管采用多因素分析模型处理,但仍有可能无法完全排除可能的混杂因素,因此在结论外推时仍需谨慎。
[1] | World Health Organization. Influenza at the human-animal interface: summary and assessment as of 30 October to 7 December 2017[EB/OL]. (2017 – 12 – 31)[2017 Dec 31]. http://www.who.int/influenza/human_animal_interface/Influenza_Summary_IRA_HA_interface_12_07_2017.pdf?ua=1. |
[2] | Li J, Chen J, Yang G, et al. Case-control study of risk factors for human infection with avian influenza A (H7N9) virus in Shanghai, China, 2013[J]. Epidemiol Infect, 2015, 143(9): 1826–1832. DOI:10.1017/S0950268814003264 |
[3] | 何权瀛, 高莹慧. 关于吸烟问题若干名词定义[J]. 中华结核和呼吸杂志, 2009, 32(1): 56. DOI:10.3760/cma.j.issn.1001-0939.2009.01.008 |
[4] | 任瑞琦. 人感染H7N9禽流感流行病学特征及发病危险因素研究[D]. 中国疾病预防控制中心, 2015. |
[5] | Ai J, Huang Y, Xu K, et al. Case-control study of risk factors for human infection with influenza A(H7N9) virus in Jiangsu province, China, 2013[J]. Euro Surveill, 2013, 18(26): 20510. DOI:10.2807/1560-7917.ES2013.18.26.20510 |
[6] | 王富良, 孙向珏, 蔡维未, 等. 2013 — 2016年浙江省杭州市萧山区人感染H7N9禽流感病例流行病学特征分析[J]. 中国人兽共患病学报, 2017, 33(3): 208–211. DOI:10.3969/j.issn.1002-2694.2017.03.004 |
[7] | 陈恩富, 柴程良, 孙继民, 等. 浙江省人感染H7N9禽流感流行特征与防控对策[J]. 中国公共卫生, 2013(5): 625–627. |
[8] | 陈金堃, 董恒进, 傅利军, 等. 浙江省绍兴市人感染H7N9禽流感流行病学特征和外环境监测结果[J]. 疾病监测, 2015(2): 122–125. |
[9] | 王素春, 蒋文明, 侯广宇, 等. 全面取代活禽市场的必要性和可行性分析[J]. 中国动物检疫, 2017, 34(6). DOI:10.3969/j.issn.1005-944X.2017.06.011 |
[10] | 邓立权, 范明, 黄飚, 等. 吉林省活禽市场外环境禽流感污染情况调查[J]. 中国公共卫生, 2017(3): 493–495. |
[11] | 邹丽容, 毛小晓, 武婕, 等. 2012 — 2015年广东省活禽市场外环境禽流感病毒污染状况研究[J]. 激光生物学报, 2016(3): 276–282. DOI:10.3969/j.issn.1007-7146.2016.03.015 |
[12] | Kile JC, Ren R, Liu L, et al. Update: increase in human infections with novel Asian lineage avian influenza A(H7N9) viruses during the fifth epidemic – China, October 1, 2016 – August 7, 2017[J]. MMWR Morb Mortal Wkly Rep, 2017, 66(35): 928–932. DOI:10.15585/mmwr.mm6635a2 |
[13] | Su W, Cheng KL, Chu D, et al. Genetic analysis of H7N9 highly pathogenic avian influenza virus in Guangdong, China, 2016 – 2017[J]. J Infect, 2018, 76(1): 93–96. DOI:10.1016/j.jinf.2017.09.001 |
[14] | Qi W, Jia W, Liu D, et al. Emergence and adaptation of a novel highly pathogenic H7N9 influenza virus in birds and humans from a 2013 human-infecting low-pathogenic ancestor[J]. Journal of Virology, 2018, 92(2): e00921–17. |
[15] | Deng Y, Li C, Han J, et al. Phylogenetic and genetic characterization of a 2017 clinical isolate of H7N9 virus in Guangzhou, China during the fifth epidemic wave[J]. Sci China Life Sci, 2017, 60(12): 1331–1339. DOI:10.1007/s11427-017-9152-1 |
[16] | Teng Y, Bi D, Guo X, et al. Contact reductions from live poultry market closures limit the epidemic of human infections with H7N9 influenza[J]. J Infect, 2018, 76(3): 295–304. DOI:10.1016/j.jinf.2017.12.015 |
[17] | Wu J, Zou L, Ni H, et al. Serologic screenings for H7N9 from three sources among high-risk groups in the early stage of H7N9 circula-tion in Guangdong province, China[J]. Virol J, 2014, 11–184. DOI:10.1186/1743-422X-11-184 |
[18] | Lu J, Wu J, Zeng X, et al. Continuing reassortment leads to the genetic diversity of influenza virus H7N9 in Guangdong, China[J]. J Virol, 2014, 88(15): 8297–8306. DOI:10.1128/JVI.00630-14 |
[19] | Creager HM, Zeng H, Pulit-Penaloza JA, et al. In vitro exposure system for study of aerosolized influenza virus[J]. Virology, 2017, 500: 62–70. DOI:10.1016/j.virol.2016.10.007 |
[20] | 阮菁如, 黄佳峰, 钟春燕, 等. 思明区人感染H7N9禽流感流行特征及防控[J]. 海峡预防医学杂志, 2017, 23(1): 25–27. |
[21] | Zhou L, Ren R, Ou J, et al. Risk factors for influenza A (H7N9) disease in China, a matched case control study, October 2014 to April 2015[J]. Open Forum Infect Dis, 2016, 3(3): w182. DOI:10.1093/ofid/ofw182 |
[22] | 关玉贤. 福建省森林覆盖率影响因素分析及对策研究[J]. 林业勘察设计, 2015(1): 10–13. DOI:10.3969/j.issn.1004-2180.2015.01.003 |
[23] | Gao R, Cao B, Hu Y, et al. Human infection with a novel avian-origin influenza A (H7N9) virus[J]. N Engl J Med, 2013, 368(20): 1888–1897. DOI:10.1056/NEJMoa1304459 |
[24] | Zhong Y, Qin Y, Yu H, et al. Avian influenza virus infection risk in humans with chronic diseases[J]. Sci Rep, 2015, 5: 8971. DOI:10.1038/srep08971 |