刘晓美 范亚光 姜勇 向剑 王继先 孙志娟 任冠华 姚树祥 常润生 赵永成 乔友林 周清华
肺癌是世界范围内首要癌症死亡原因[1]。依据第3次全国恶性肿瘤回顾性死因抽样调查结果,我国肺癌死亡率居恶性肿瘤死亡率首位(男性和女性年龄标化死亡率分别为39.1/105、16.7/105)[2]。自20世纪70年代发现我国云南锡矿矿工(下文简称云锡矿)肺癌发病率、死亡率高且已被确认为职业癌。姚树祥等[3]研究发现,在1954年-2002年间,云锡矿男性矿工肺癌年平均发病率为187.7/105,而有30年以上井下工作史人群的发病率则高达1,694.9/105,年平均死亡率为161.0/105,远高于当地普通男性居民肺癌标化死亡率(26.6/105)。
吸烟被公认为肺癌首要危险因素,此外环境烟草烟雾、氡、职业暴露等其它因素也增加肺癌发病危险[4-6]。研究[7]已证明云锡矿工肺癌死亡率高与生产环境中的氡、砷等职业暴露因素有关。本研究拟利用云锡矿工职业高危人群前瞻性队列研究资料,探讨云锡矿工肺癌危险因素及交互作用。
1.1 研究人群 研究对象来自于1992年-1999年参加“云锡矿工肺癌早期标志物”队列研究的成员,其纳入标准为:①年龄为40岁及以上;②至少有10年井下工作史或冶炼史或两者相加满10年;③无恶性肿瘤史;④知情同意。本研究中,研究对象至少参加过1次1992年-1999年年度性肺癌筛查。
1.2 肺癌病例的确认和死亡 1992年-1999年,每年对队列成员进行一次痰细胞学和X线胸片筛查。肺癌病例确认途径主要为:①肺癌筛查阳性者住院进一步检查后确诊;②因症状于云锡总医院就医而确诊;③根据云锡肿瘤登记系统提供的资料确诊;④个别病例根据随访结果死因推断。
1.3 暴露情况
1.3.1 累积氡暴露量资料 根据云锡矿井作业环境中氡监测浓度,按照国际通用计算方法计算。累积氡暴露量的常用单位为工作水平月(working level month, WLM),1 WLM相当于在1 WL下暴露1个工作月(每年工作2,000 h/12 m≈170 h)[8]。
累积氡暴露量计算:按照氡暴露健康效应委员会(Committee on Health Effects of Exposure to Radon, BEIR VI)给出的暴露量-年龄-暴露持续时间(exposure-ageduration, EAD)模型计算[9](暴露截止时间为1995年12月31日)。模型表达式如下所示:
W5-14、W15-24、W25+分别表示发病前5年-14年、15年-24年及≥25年间的累积氡暴露量,θ15-24、θ25+表示发病前15年-24年及≥25年间氡暴露危险的相对贡献。其具体参数为:θ5-14=1.00,θ15-24=0.72,θ25+=0.44。累积氡暴露量W(WLM)=W5-14+θ15-24W15-24+θ25+W25+。
1.3.2 氡暴露率 氡暴露率计算方法有多种,本文的计算方法为:肺癌潜伏期外的累积氡暴露量/潜伏期外的暴露时间(以年为单位)。
1.3.3 砷暴露量 以砷暴露指数(index of arsenic exposure months, IAEM)估算矿工砷暴露量,计算方法:按照不同年代作业环境中砷浓度(mg/m3)与累积月之乘积做累加。
1.3.4 吸烟量 用吸烟指数(包/日*年,pack-years)表示,为每天吸纸烟包数乘以吸烟年数。对于旱烟或水烟,按每1.25 g旱烟或水烟为一支纸烟计算[10]。
1.3.5 年龄、性别、教育水平、既往肺部疾病史 上述信息均在研究对象进入队列时填写基线调查表获得。既往肺部疾病由医生询问获得,此外,矽肺必须有云锡劳动保护研究所开具尘肺诊断证明书才能确认。具体分组情况见表1。
1.4 统计学分析 计算不同人群特征和职业暴露肺癌死亡率,利用Cox回归模型分析云锡矿工肺癌危险因素,分析氡暴露人群中初始氡暴露年龄和暴露率对肺癌危险的影响,分析不同吸烟、氡暴露水平下肺癌死亡的危险。基于相加模型和相乘模型分析吸烟和累积氡暴露量的交互作用:①在Cox回归模型中纳入乘积项分析吸烟和累积氡暴露量之间的相乘交互作用;②利用Rothman提出的针对Logistic或Cox回归模型的三个评价相加交互作用的指标评价吸烟与累积氡暴露量的相加交互作用[11],分别为①相对超危险度比(relative excess risk due to interaction,RERI);②归因比(attributable proportion due to interaction, AP);③交互作用指数(synergy index, S)。三个指标及其95%CI的计算根据Li的方法[12]。如果两因素无相加交互作用,RERI和AP的可信区间应包含0,S的可信区间应包含1,以双侧P<0.05为差异有统计学意义。
表 1 云锡矿工队列人群特征及肺癌死亡分布Tab 1 Characteristics and lung cancer mortality of Yunnan tin miner cohort
2.1 队列特征 1992年-1999年,进入队列并参加过肺癌年度性筛查的云锡矿工共9,295人,截止到2001年12月31日,确诊肺癌死亡病例443人。表1显示不同人群特征和职业暴露肺癌死亡率,从肺癌死亡率比可以看出,年龄、性别、教育水平、吸烟、职业性氡砷暴露、既往肺疾病史均为肺癌的可能危险因素。表2为云锡矿工肺癌危险因素的多因素Cox回归分析结果,表明年龄、吸烟量、职业性氡砷暴露量、慢性支气管炎史是此队列人群的独立危险因素,而教育是保护性因素。
2.2 初始氡暴露年龄、氡暴露率与肺癌的关系 前瞻性队列人群中共9,295人,其中有氡暴露者7,406人,表3显示了7,406名氡暴露矿工肺癌危险与初始氡暴露年龄和氡暴露率的Cox多因素回归分析结果:在分析中,当初始氡暴露年龄和氡暴露率为分类变量时,初始氡暴露年龄、氡暴露率与肺癌无显著的剂量效应关系;而当初始暴露年龄和氡暴露率为连续变量时,氡暴露率是肺癌死亡危险的保护性因素,调整肺癌死亡危险比(hazard ratio, HR)=0.98(0.96-0.99)。
2.3 不同吸烟、氡暴露水平下肺癌死亡的危险 由表4可见,相同的吸烟水平下,肺癌的危险随累积氡暴露量的增大而增大,在不吸烟者中,低水平氡暴露、高水平氡暴露的HR分别为5.75(95%CI: 1.26-26.12)和11.66(95%CI: 3.24-42.05)。在相同的累积氡暴露水平内,在高剂量组中肺癌的危险随吸烟量的增大而增大。与无吸烟和氡暴露组相比,吸烟量>43.7 pys且氡暴露>118.44 WLM的HR高达17.15(95%CI: 5.33-55.16)。
2.4 累积氡暴露量与吸烟对肺癌危险的交互作用分析 基于相加模型和相乘模型分析吸烟和累积氡暴露量的交互作用,结果如表5所示。在Cox回归模型中纳入乘积项分析吸烟和累积氡暴露量之间的相乘交互作用,乘积项系数β为0.09(HR=1.10, 95%CI: 0.57-2.12),无统计学意义,尚不能认为二者对肺癌危险存在单纯的相乘交互作用。
利用Rothman提出的针对Cox回归模型的三个评价相加交互作用的指标评价吸烟与累积氡暴露量的相加交互作用,结果如表5所示,RERI、AP、S分别为0.92(0.11-1.73)、0.25(-0.05-0.55)和1.52(0.85-2.70),其95%区间分别大于0、包括0、包括1,说明累积氡暴露量与吸烟对云锡矿工肺癌的死亡有相加交互作用。图1为累积氡暴露量与吸烟相加交互作用示意图。
表 2 云锡矿工肺癌危险因素的多因素分析Tab 2 Multiple analysis for lung cancer risk among Yunnan tin miners
表 3 氡暴露云锡矿工肺癌危险与初始氡暴露年龄和氡暴露率的相关分析Tab 3 Relationship between lung cancer risk and age at first radon exposure, rate of radon exposure among radon-exposed miners
表 4 按吸烟量分层后肺癌调整危险比与累积氡暴露量的相关分析Tab 4 Relation of cumulative radon exposure to lung cancer risk by smoking
表 5 累积氡暴露量和吸烟对肺癌危险的交互作用分析Tab 5 Interaction of cumulative radon and smoking on lung cancer risk
图 1 累积氡暴露量和吸烟交互作用示意图Fig 1 Additive interaction of cumulative radon and smoking on lung cancer risk
国内外研究结果表明,矿工肺癌发病与多种暴露因素有关。姚树祥等[13]曾报道云锡矿工的主要危险因素为职业性因素,如氡、砷、粉尘等;同时还与一些非职业性危险因素有关,如吸烟、慢性支气管炎史、文化水平等。我们的研究除分析各危险因素对云锡矿工肺癌死亡危险的影响外,还对初始氡暴露年龄、氡暴露率及吸烟和氡暴露交互作用对肺癌死亡危险的影响进行了分析。不同人群特征及职业暴露水平下肺癌死亡率的差别显示年龄、性别、文化水平、吸烟、职业性氡砷暴露、肺疾病既往史等因素与肺癌死亡有关,继而多因素Cox回归分析结果显示,年龄、慢性支气管炎史、吸烟、职业性氡砷暴露是肺癌的独立危险因素,教育水平为保护性因素。
Xuan等[14]、Lubin等[15]研究发现肺癌危险与初始氡暴露年龄有关。本研究利用Cox回归模型分析了初始氡暴露年龄与肺癌危险的关系,分别把初始氡暴露年龄按分层变量和连续性定量变量进行分析,肺癌死亡危险与初始氡暴露年龄间均无统计学意义。其他研究[16]也有类似报道,初始氡暴露年龄与肺癌危险的关系有待进一步深入研究。用同样的方法分析肺癌危险与氡暴露率的关系,当暴露率为连续性变量时,肺癌危险与氡暴露率之间存在逆剂量率效应(HR=0.98, 95%CI: 0.96-0.99),这与之前的报道结果一致[17,18]。
本研究对吸烟和累积氡暴露量对肺癌危险的作用作了重点分析。吸烟被公认为是引起肺癌最重要的危险因素,为分析累积氡暴露量对肺癌的影响,同时排除吸烟的影响,本研究在校正混杂因素后采用分层分析,结果显示随累积氡暴露量的增加肺癌危险递增;无氡暴露吸烟者肺癌危险要远小于有氡暴露不吸烟者。Taylor等[19]曾在研究中指出,云锡矿工肺癌病因中,砷暴露<吸烟<氡暴露,本文结果与上述结果一致。为进一步研究吸烟和累积氡暴露对肺癌危险的交互作用,本文又对二者之间的相加交互作用和相乘交互作用进行分析,结果显示二者之间存在相加交互作用,相乘交互作用无统计学意义。Xuan等[14]就吸烟者和非吸烟者的研究发现:吸烟和氡暴露存在介于相加和相乘作用之间的亚相乘作用,Lubin等[20]的研究也支持二者的关系介于相加和相乘模型之间。Breslow等[21]提出利用广义相对危险度模型分析因素间的包括相加、相乘的一系列交互作用形式,本研究中仅对相加和相乘交互作用进行了分析,我们将继续对吸烟和氡暴露的交互作用进行深入分析。
本研究的优点在于:累积氡暴露量的计算较以往的计算方法更科学;研究为对云锡矿工队列的随访,较以往云锡肺癌危险研究随访时间更长,随访的肺癌病例更多,此外职业性氡砷暴露、吸烟、肺疾病等信息是在基线调查中获得,均在肺癌发病、死亡之前。缺点在于:云锡矿工肺癌危险因素较多,在分析过程中尽管采用多因素分析方法,但在分析过程中不能保证控制得当。
总之,在研究职业性肺癌的危险因素时要注意危险因素种类的复杂性。通过本研究结果发现,控制矿工肺癌的关键是控制肺癌危险因素,通过通风除尘、湿式作业等途径降低作业环境中氡、砷以及粉尘的浓度;控制吸烟等不良嗜好;尽量减少职业性暴露机会;严格控制矿工的纳入标准等措施的实施,从而以最大的可能控制和降低职业性肺癌的发病率和死亡率。