论文部分内容阅读
目的:1.分析孕早期PAHs宫内暴露及代谢酶（CYP450、GST）基因多态性与胎停育的相关性。2.分析PAHs暴露对端粒长度（TL）的影响,并探讨端粒长度与胎停育的关系。3.探讨影响胎停育的环境、遗传因素及交互作用。方法:采用病例对照的流行病学方法。连续收集2019年3月-2019年12月在山西医科大学第一医院产科进行产前检查并确诊的胎停育妇女为病例组。选取同期同医院正常早孕但自愿终止妊娠的妇女为对照组,共纳入114名病例和139名对照。采用面对面问卷调查的方式,收集研究对象一般情况、既往疾病史、不良妊娠史、孕期生活行为方式等。经负压吸引术获得流产的绒毛组织标本,用生理盐水冲洗去除残留的血迹和杂质,低温运送至实验室,-80℃保存。选取一定量的绒毛组织,按比例（1:9）加入生理盐水,使用匀浆器进行组织匀浆,随后离心取上清液。采用双抗夹心酶联免疫吸附试验（ELISA）检测组织上清液PAH-DNA加合物含量。采用酚氯仿法提取组织DNA,通过多重PCR扩增和高通量测序法进行基因分型。采用荧光定量PCR（Q-PCR）法检测端粒相对长度。采用Epidata 3.1软件进行问卷录入。正态分布资料采用均数±标准差（x±s）表示,非正态分布资料采用中位数M(P25～P75)表示。采用SPSS 23.0进行t检验、卡方检验、Logistic回归分析、Mann-Whitney U检验和Kendall’s tau-b相关性分析等。采用SAS 9.2进行Armitage趋势检验。采用R 3.6.3软件进行哈迪-温伯格平衡（Hardy-Weinberg equilibrium,HWE）检验和单体型分析。采用GMDR 0.9软件进行基因-基因、基因-环境交互作用分析。结果:1.PAH-DNA加合物检出情况:病例组、对照组PAH-DNA加合物含量分别为466.26（424.55～507.53）pg/ml,446.51（384.14～506.19）pg/ml。病例组PAH-DNA加合物水平高于对照组（Z=-2.10,P=0.036）。2.PAH-DNA加合物水平与胎停育的相关性:调整协变量后,与PAH-DNA加合物低水平组（Q1）相比,Q2、Q3和Q4组胎停育的发生风险分别为3.79、3.59和3.05倍（调整后OR=3.79,95%CI:1.67-8.58,P=0.001;调整后OR=3.59,95%CI:1.59-8.12,P=0.002;调整后OR=3.05,95%CI:1.33-7.00,P=0.008）。3.哈迪-温伯格平衡检验（Hardy-Weinberg equilibrium test）:纳入分析的11个基因位点均符合Hardy-Weinberg遗传平衡检验（P>0.05）。4.两组各位点基因型、等位基因频率比较:胎停育与GSTM1 rs1065411基因位点的基因型有关联（P=0.005）。病例组与对照组GSTM1 rs1065411位点等位基因频率分布不同（P<0.001）,胎停育组等位基因C所占比率较高。其余位点基因型、等位基因与胎停育无统计学关联（P>0.05）。5.各基因位点与胎停育易感性的Logistic回归分析:调整协变量后,GSTM1rs1065411基因位点与胎停育有关联（P<0.05）。其中GC、CC基因型携带者胎停育发生的风险分别为GG基因型的4.82倍、3.29倍（调整后OR=4.82,95%CI:1.74-13.42,P=0.003;OR=3.29,95%CI:1.55-7.00,P=0.002）。GSTM1 rs449856基因位点TA基因型携带者胎停育的发生风险为TT基因型的0.21倍（调整后OR=0.21,95%CI:0.06-0.71,P=0.044）。6.CYP450、GST基因与胎停育发病的单体型关联分析:GSTM1 rs1065411、rs449856基因位点创建的单体型结果显示（调整协变量）,单体型G-T、G-A胎停育的发病风险分别为C-T的1.50、4.00倍（调整后OR=1.50,95%CI:1.10-2.10,P=0.012;调整后OR=4.00,95%CI:1.40-11.00,P=0.010）。7.基因-基因交互作用分析:调整协变量,所测得CYP450各基因位点间不存在交互作用（P>0.05）。GSTM1 rs1065411、GSTM3 rs7483与GSTP1 rs1695基因位点多态性间存在交互作用,该模型训练集平衡精度为0.6472,测试集平衡精度为0.6159,交叉验证一致性为10/10,P=0.011。所测得CYP450与GST基因位点间不存在交互作用（P>0.05）。8.基因-环境交互作用分析:调整协变量,CYP450基因位点多态性与PAH-DNA加合物之间不存在交互作用（P>0.05）。GSTM1 rs1065411基因多态性与PAH-DNA加合物间存在交互作用,该模型训练集平衡精度为0.6620,测试集平衡精度为0.6432,交叉验证一致性为10/10,P=0.001。9.端粒相对长度检出情况:病例组、对照组端粒相对长度差异具有统计学意义（Z=-4.02,P<0.001）,对照组端粒相对长度大于病例组。10.端粒相对长度与胎停育发病风险的关联性分析:根据端粒相对长度的中位数分为两组。调整协变量后显示,端粒较短组（M2）胎停育的发生风险为端粒较长组（M1）的2.11倍（调整后OR=2.11,95%CI:1.22-3.63,P=0.007）。11.PAHs暴露与端粒长度的相关性分析:Kendall’s tau-b相关性分析显示,PAH-DNA加合物水平与端粒长度间不存在相关性（r=0.08,P=0.170）。12.PAHs暴露对端粒长度影响的Logistic回归分析:调整协变量后,PAH-DNA加合物为Q3水平时造成端粒长度变短的风险为低水平Q1组的2.67倍（OR=2.67,95%CI:1.28-5.58）,差异具有统计学意义（P=0.009）。13.影响胎停育的多因素Logistic回归分析:将第一、二部分分析结果有统计学意义的因素纳入Logistics回归分析,结果显示有不良妊娠史（OR=3.08,95%CI:1.61-5.91）、PAH-DNA加合物Q2、Q3水平（OR=2.93,95%CI:1.32-6.52;OR=2.77,95%CI:1.23-6.25）、GSTM1 rs1065411位点GC、CC基因型（OR=3.37,95%CI:1.46-7.80;OR=2.57,95%CI:1.40-4.69）、端粒相对长度较短（OR=2.05,95%CI:1.19-3.56）均为胎停育的危险因素。结论:1.PAHs高暴露、GSTM1 rs1065411基因位点变异可能是胎停育的危险因素。2.PAHs暴露、GST基因多态性间存在基因-基因、基因-环境交互作用,可共同影响胎停育的发生。3.PAHs高暴露可能会造成端粒长度的磨损,而端粒长度较短与胎停育有关。