候选基因多态性与慢性苯中毒的关联分析研究苯主要通过其代谢产物引致接触者的血液毒性和遗传毒性。慢性苯中毒的发生存在着遗传易感性。代谢酶基因、DNA损伤修复基因、细胞周期调控基因的多态性与慢性苯中毒遗传易感性的关系,本课题组已进行了多次单独的研究探讨。我们在研究中注意到以下几个问题：①以往的研究结果是建立在对苯的代谢、DNA损伤修复、细胞周期调控分别进行分析的基础上得到的。由于苯中毒的发生与以上3个方面均可能有关,如果仅分析某一机制的作用,则可能使研究结果产生偏倚。②单个SNP的作用可能是非常微弱的,研究者在进行统计分析时,需要考虑SNP的交互作用。但是,当样本量固定,而研究的SNP个数逐渐增多,在考虑交互作用时,则估计的参数呈指数形式增加,进行Logistic回归会导致非常大的估计偏差。③在所研究的候选基因多态中,哪些对苯中毒发生的影响比较大,哪些有潜在的影响,目前还不明确。应该选取重要的变量来预测苯中毒发生的概率,从而提高预测的准确性。针对上述问题,本次研究在职业流行病学调查和实验室研究的基础上,综合分析了候选基因多态与慢性苯中毒的关联,包括基因的主效应、基因与基因的交互作用、基因与环境因素的交互作用。分析策略主要为：①在单个位点分析时控制多重检验的阳性错误发生率(FDR),并筛选SNP用于Logistic回归模型,以分析主效应和2阶交互效应；②控制吸烟等混杂因素的影响,分析单体型与慢性苯中毒之间的关联；③采用数据挖掘中的多因子降维法(MDR)进一步分析高阶交互作用；④采用随机森林(Random Forest)分析各个SNP影响苯中毒发生的重要性,并根据重要性评分,筛选相对重要的SNP以利今后进一步研究。为探讨候选基因多态性与慢性苯中毒的关联,我们开展了病例-对照研究,病例组为152名苯中毒工人,对照组为152名调查时正在从事接苯作业而没有慢性苯中毒表现的工人。采用Dosmeci等的接触评估方法,对研究对象的累积接触水平进行了评估,并调查了他们吸烟、饮酒等情况。职业流行病学调查结果表明,病例组和对照组的一般特征如性别、民族、年龄、工龄和累积接触评分在两组人群中的分布无统计学差异,表明病例组和对照组均衡可比。根据苯中毒的发生机制,检测了3类基因的SNP。代谢酶基因：CYP2E1、MPO, NQO1、GSTT1、GSTM1、GSTP1、EPHX1、EPHX2、UGT1A6、UGT1A7、SULT1A1、CYP1A1、CYP2D6,共33个SNP。DNA损伤修复基因：hMTH1、hOGG1、hMYH、XPD、APE1、XRCC1、ADPRT、XRCC2、XRCC3、共14个SNP。细胞周期调控基因：p53, p21、mdm2、gadd45α、p14ARF,共13个SNP。通过排除未检出变异基因型的以及严重偏离HWE的SNP,最后共有51个SNP用于关联研究(代谢酶基因、DNA损伤修复酶基因、细胞周期调控基因分别有33个、10个、8个SNP)由于SNP个数较多,在进行Logistic回归时考虑交互作用的策略如下：首先筛选出单变量检验P值小于0.05的SNP作为自变量,环境因素为协变量；再用Forward法建立Logistic回归模型分析各SNP的主效应、2阶交互作用；最后,将以往研究得出的无主效应但有交互作用的SNP及其交互项也纳入到Logistic回归模型(Backward法)中,分析其在新的模型中是否还存在交互作用。Logistic回归分析结果如下。①未检测到主效应的有以下40个SNP。代谢酶基因：CYP2E1所有3个SNP; GSTT1缺失多态；GSTM1缺失多态；MPO的所有4个SNP; NQO1 rs113134、rs1800566; EPHX1 rs2234922、rs1051741、rs2854451、rs3738047; EPHX2 rs751141; UGT1A6的所有7个SNP; UGT1A7 rs11692021; SULT1A1 rs9282861; CYP1A1 rs4646421、rs4646422、rs1048943。DNA损伤修复基因：XPD rs13181; APE1 rs1130409; XRCC1 rs25487; ADPRT的rs1136410; XRCC3的rs861539。细胞周期调控基因：p53的所有3个SNP；p14ARF rs3731217、rs3088440; gadd45αrs581000、rs532446;②检测出有主效应的有如下11个SNP。代谢酶基因：GSTP1 rs947894、EPHX1 rs1051740; CYP1A1rs4646903：CYP2D6*10 rsl065852和rsl 135840。DNA损伤修复基因：hMTH1 rs4866、hOGG1 rs1052133、hMYH rs3219489、XPD rs1799793、XRCC1 rs1799782。细胞周期调控基因：p21 rs1059234。③检出的交互作用有如下10种：CYP2E1rs3813867与EPHX1 rs3738047; EPHX1 rs3738047与饮酒；GSTP1 rs947894与饮酒;CYPlAl 4646903与CYP2D6 rs1135840; hMTH1 rs4866与XRCC1 rs1799782; hOGG1 rs1052133与hMYH rs3219489; hMYH rs3219489与吸烟；XRCC1 rs1799782与APE1 rs1130409; APE1 rs1130409与饮酒;hOGG1rs1052133与XPD rs1799793.与以往研究结果不同的是：①本次研究新检测出了有主效应的2个SNP：GSTP1 rs947894和hMYH rs3219489。②本次研究新检测出的交互作用有4种：CYP2E1 rs3813867与EPHX1 rs3738047; XRCC1 rs1799782与APE1 rs1130409;CYPlAl rs4646903与CYP2D6 rs1135840; hMTH1 rs4866与XRCC1 rs1799782。③以往研究中有交互作用,在本次研究未得到验证的有：NQOlrs1800566与吸烟或饮酒的交互作用。本次研究的单体型分析考虑了吸烟、饮酒、苯接触强度等协变量的影响。与以往研究相同的是,均检出CYP2D6*10的单体型与慢性苯中毒存在关联,携带CC单体型的个体与携带TC单体型相比,易感性增加。而以往研究的单体型与苯中毒存在关联,本次研究没有得到验证的有：EPHX、UGT1A6、CYP1A1、XRCC1等单体型。对代谢酶基因多态性研究的MDR模型中,我们发现了1个3阶交互作用,即CYP1Al rs4646903、CYP2D6 rsl065852、CYP2D6 rsll35840的3因子组合。有3种组合被归入高危组：①CYP2D6 rs1135840的CC基因型、CYPlAl rs4646903的TC+CC基因型、CYP2D6 rsl065852的CT+CC基因型的个体；②CYP2D6rs1135840的CC基因型、CYPlAl rs4646903的TT基因型、CYP2D6 rsl065852的CT+CC基因型的个体；③CYP2D6 rs1135840的CG+GG基因型、CYPlAlrs4646903的TT基因型、CYP2D6 rsl065852的CT+CC基因型的个体。其它5种组合属于低危组。对DNA损伤修复、细胞周期调控基因多态性研究的MDR模型中,未发现3阶以上的高阶交互作用。hOGG1 rs1052133和XPD rs17997933因子组合,是预测效果最佳的2阶交互作用,hOGG1 rs1052133的GG基因型和XPD rs1799793的GG基因型的个体属于高危组,其它3种基因型组合属于低危组。对全部3类基因多态性研究的MDR模型中,最佳的因子模型与对代谢酶基因多态性研究的MDR模型一致。随机森林法对代谢酶、DNA损伤修复酶和细胞周期调控基因多态性影响慢性苯中毒的重要性排序结果显示,重要性从大到小排列的前25个SNP或环境因素依次为：EPHX1 rs1051740、hOGG1 rs1052133、CYP2D6 rs1065852、CYP1A1 rs4646903、CYP2D6 rs1135840、p21rs1059234、XPD rs1799793、p53 rs17878362、hMYH rs3219489、hMTH1 rs4866、GSTP1 rs947894、EPHX1 rs1051741、CYP2E196-bp插入、MPO rs7208693、EPHX1 rs3738047、XRCC1 rs25487、SULT1A1rs9282861、吸烟、UGT1A6 rs6759892、UGT1A7 rs11692021、XRCC1 rs1779782. NQO1 rs1800566、XPD rs13181、UGT1A6 rs2070959、GSTT1。重要性评分最高的EPHX1 rs1051740是如何通过交互作用对慢性苯中毒产生影响的,还需要进一步研究。除了p53 rs17878362外,排在第2至第11位的9个SNP中,均在Logistic回归分析中检测到,且大部分无主效应,但存在交互作用。对于Logistic回归分析中由于样本量的限制而未能检测到效应的SNP,随机森林的结果能进一步检测它们的重要性,从而提示哪些SNP最有可能与慢性苯中毒有关联,以筛选出这些相对重要的SNP,再加大样本量进行研究根据本次研究结果,发现对代谢酶基因、DNA损伤修复基因、细胞周期调控基因的多态性与慢性苯中毒关联的综合研究结果与以往的研究结果基本相同,这证明我们已有的研究结论是比较可靠的。本次研究得到的单体型关联分析的结果、高阶交互作用分析的结果、以及各候选基因多态性相对重要性的比较,是对已有研究结果的进一步补充,有利于进一步阐明候选基因多态性与慢性苯中毒的关系。综合本次研究结果,发现与慢性苯中毒发生关系密切的代谢酶基因、DNA损伤修复基因、细胞周期调控基因的多态性与个体慢性苯中毒存在关联。基因-基因、基因-环境交互作用是影响个体慢性苯中毒遗传易感性的重要方式。这些研究结果可为更好地做好接苯工人的健康监护工作和筛选有效的易感性生物标志提供理论依据。本次研究结果中,也有很多SNP未发现有主效应或交互作用,一方面可能这些SNP的作用确实比较微弱或者缺乏,另一方是可能是样本量较小的原因。因此本课题组的下一步工作将是进行SNP的筛选(例如,在随机森林方法对各SNP重要性排序的提示下,依靠生物学意义和统计学变量筛选方法开展筛选工作),以及扩大样本量进一步研究。