CN102495214B - 鞘脂激活蛋白原作为效应标志物在检测多氯联苯雄性生殖毒性的应用 - Google Patents

Abstract

鞘脂激活蛋白原作为效应标志物在检测多氯联苯雄性生殖毒性的应用，涉及鞘脂激活蛋白原。鞘脂激活蛋白原是多功能的糖蛋白，分子量为70kD，存在于一些生物体液中。鞘脂激活蛋白原的筛选方法：用不同剂量的多氯联苯对小鼠进行暴露后收集样本；样本差异表达基因的筛选；对筛选的差异表达基因进行功能聚类，获得与生殖功能相关的差异表达基因；利用实时荧光定量PCR技术和免疫印迹法在mRNA转录水平和蛋白表达水平上对获得的差异表达基因随着PCB暴露剂量以及暴露时间的增加而变化进行检测分析；进行精子参数的检测分析；对效应标志物、体内多氯联苯水平以及精液指标之间的相关系数通过斯皮尔曼相关分析进行统计。

Description

鞘脂激活蛋白原作为效应标志物在检测多氯联苯雄性生殖毒性的应用

技术领域

本发明涉及鞘脂激活蛋白原，尤其是涉及鞘脂激活蛋白原作为效应标志物在检测多氯联苯雄性生殖毒性的应用。

背景技术

多氯联苯(Polychlorinated Biphenyls，PCBs)是一类人工合成的有机物。根据氯原子取代数和取代位置的不同共有209种同族异构体。由于PCBs具有良好的化学惰性、抗热性、不可燃性和高介电常数等优点，早期曾用于热介质、特殊润滑油、可塑剂、涂料、防尘剂、油墨添加剂、杀虫剂及复写纸等的制造和在电容器、变压器等电力设备中作为绝缘油。随着PCBs的不断应用，在环境中的积累不断增加。而且PCBs一旦进入环境，极难降解，在整个生态系统中循环传输，从而造成“全球性的环境污染”。由于PCBs脂溶性较高，易于在脂肪组织中蓄积，沿食物链逐级富集和放大，在含脂量的脏器中含量明显要高，在生物体内的富集系数高达10⁵。人类暴露于PCBs主要是通过摄食含有PCBs的鱼类肉类和乳制品等，PCBs进入体内会通过血液转运到肝脏、肌肉和脂肪，并在那逐渐累积，产生积累性中毒，严重影响人类健康。

现有的研究表明，雄性生殖系统是PCBs作用的重要靶器官之一，PCBs对生精细胞、支持细胞和间质细胞均可产生明显的毒性效应。但现有的结果因为：1)动物实验染毒剂量远远高于生物体和自然环境中的暴露剂量，无法估计真实的人群效应；2)PCBs各种睾丸毒性效应产生的分子机制模糊不清，大量人群流行病学调查得到的提示信息没有后续实验的支持；3)细胞、动物和人群3个层次的研究结果契合度不高。由于上述原因，目前尚未找到环境水平PCBs暴露对雄性动物生殖系统是否产生毒性作用的检测方法。

生物标志物(biomarker)是用于监测和评价能够导致生物有机体的生物化学和生理学改变的化学污染物。生物标志物在亚个体和个体水平上既可以测定污染物暴露水平，也可以测定污染物效应的生理和生化指标。生物标志物是生物体受到严重损害之前，在不同生物学水平(分子、细胞、个体等)上因受环境污染物影响而异常化的信号指标，它可以对严重毒性伤害提供早期警报。生物标志物从功能上一般分为：接触(暴露)生物标志物；效应生物标志物；敏感性生物标志物([8]Biological markers in environmental health research.Committee onBiological Markers of the National Research Council[J].Environ Health Perspect，1987，74：3-9)。作为生物标志物一般必须具有以下一些特性：1)具有一定的敏感性，敏感性应高于一般生物检测指标，低剂量下就可测出，可微量操作；2)具有反应的时间效应。反应要有一定的稳定时间，同时要快速；3)效应标志物在分子和生化水平上的效应要与高级生物学水平上的效应(如生长、繁殖)紧密相联，各级水平上的效应要有因果关系。

化学物对生态系统的毒性效应是经过基因→分子→细胞→组织→器官→个体→群体→生态系统这个途径产生影响的。其中基因和生物大分子是化学物最早产生作用的靶标。污染物的直接化学检测虽然数值准确、针对性强，但工作量大，而且无法对污染物的生态毒理学效应做出预测。而污染物造成生物死亡、生长受阻或繁殖受影响、最终导致生态系统结构破坏，这些已是污染物造成的晚期影响。因此，迫切需要发展快速有效的生态毒理学研究手段和技术对污染物的影响做出更为准确的生态毒理学预测和早期报警。相对于传统的研究方法和技术寻找污染物的作用靶标，基因芯片技术具有诸多的优势：1)基因是化学物产生毒性效应的第一级靶分子，较为敏感，也比较不易受到其他因素的影响，是用于检测痕量污染的首选指标([9]Afshari，C.A.，E.F.Nuwaysir，and J.C.Barrett，Application of complementary DNAmicroarray technology to carcinogen identification，toxicology，and drug safety evaluation[J].Cancer Res，1999，59(19)：4759-4760)；2)基因芯片提供的信息多，符合综合分析和评价的要求([10]Bartosiewicz，M.，et al.，Development of a toxicological gene array and quantitativeassessment ofthis technology[J].Arch Biochem Biophys，2000，376(1)：66-73)。但直接利用基因芯片评价化学物的毒性，预测有毒化学物对生物体的毒性效应，因操作程序复杂，所需费用高。因此，利用基因芯片技术对污染物作用靶标进行筛查，结合其它生物技术进行验证，获得化学污染物的效应生物标志物成为污染物毒性效应检测的首选方法。

鞘脂激活蛋白原(Prosaposin，psap)作为一种多功能的糖蛋白，存在于一些生物体液中，例如精液、乳汁、脑脊液([11]Kishimoto，Y.，M.Hiraiwa，and J.S.O′Brien，Saposins：structure，function，distribution，and molecular genetics[J].J Lipid Res，1992，33(9)：1255-1267)。对小鼠的Psap基因进行靶向突变能够引起雄性生殖系统的一系列异常，包括睾丸重量下降、精子数量减少以及前列腺、精囊腺和附睾的形态学变化([12]Morales，C.R.，et al.，Targeted disruption ofthe mouse prosaposin gene affects the development of the prostate gland and other malereproductive organs[J].JAndrol，2000，21(6)：765-775)。基于Psap与生殖终点的生物学相关性，Psap有可能作为检测PCBs雄性生殖毒性的效应标志物。

发明内容

本发明的目的是提供鞘脂激活蛋白原作为效应标志物在检测多氯联苯雄性生殖毒性的应用。

所述鞘脂激活蛋白原(Prosaposin，简称Psap)是一种多功能的糖蛋白，分子量为70kD，存在于一些生物体液中。

所述鞘脂激活蛋白原可采用以下方法筛选：

1)用不同剂量的多氯联苯对小鼠进行暴露后收集样本；所述不同剂量的多氯联苯可采用以下剂量：0.5μg/kg、5μg/kg、50μg/kg、500μg/kg；所述收集样本，可在暴露50天后收集样本；所述样本为附睾和肝脏样本。

2)采用基因芯片技术进行样本差异表达基因的筛选，具体方法如下：

对步骤1)所收集的正常对照样本不同剂量的多氯联苯处理的附睾样本提取的总RNA，以RNA为模板合成cDNA；利用获得的cDNA合成cRNA，片段化荧光标记的cRNA；标记的cRNA经片段化处理后用于基因芯片杂交，基因芯片进行杂交、洗脱、染色及检测；对获得的检测数据进行归一化处理，使不同的样本数据具有可比性；

3)通过生物信息学统计和分析软件，实现对筛选的差异表达基因进行功能聚类，获得与生殖功能相关的差异表达基因，具体方法如下：

对步骤2)获得的数据进行统计学分析，利用微阵列显著性分析算法(Significance Analysisof Microarray，SAM)分析差异表达基因的显著性，利用分子注释系统(Molecular AnnotationSystem，MAS)数据库，对差异表达基因以进行检索，并建立PCB暴露于附睾功能变化间的相互关系；

4)利用实时荧光定量PCR技术和免疫印迹法(Western Blot)在mRNA转录水平和蛋白表达水平上对获得的差异表达基因随着PCB暴露剂量以及暴露时间的增加而变化进行检测分析，具体方法如下：

将步骤2)提取的总RNA利用特异性引物进行实时荧光PCR，定量Psap的mRNA转录水平，所述特异性引物具有如下序列：

引物1：5’-GCCAGAGGGCAGGAGCATT-3’

引物2：5’-CTGACCCAGGGACAGCAACA-3’

样品在缓冲液中匀浆后，依次进行聚丙烯酰胺凝胶电泳和电转移，利用特异性抗体检测Psap的蛋白丰度。

所述特异性抗体为鞘脂激活蛋白原(羊抗兔)多克隆抗体。

5)进行精子参数的检测分析，具体方法如下：

将步骤1)取得的附睾尾部剪开，使用Chemicon公司HTF培养基(MR-070-D)37℃孵育30min，让精子游出，吸取10μl新鲜精子悬液于载玻片上，40倍物镜下观察精子活力，统计活动精子的总百分率，另取10μl精子悬液制备精子涂片，固定后1％伊红染色，统计每200个精子中畸形精子的百分比。

6)对效应标志物、体内多氯联苯水平以及精液指标之间的相关系数通过斯皮尔曼相关分析进行统计，具体方法如下：通过气象色谱-电子捕获检测器检测步骤1)取得肝脏中PCB的富集，将PCB的富集与步骤4)的Psap蛋白丰度，以及步骤5)的精子参数进行斯皮尔曼相关性分析。

本发明的原理是利用基因芯片检测PCB对附睾表达谱的影响，结合生物信息学工具寻找附睾中对PCB暴露敏感的基因。通过实验研究将上述敏感基因的表达与PCB在生物体内的富集以及生殖终点建立相关性，以筛选和验证PCB暴露的特异性标志物。

本发明提供了一种筛选PCB暴露标志物以及检测和评价PCB暴露的生殖毒性的方法。本发明涉及的生物学样品为非侵入性样品，对检测对象损伤小。与已有的PCB标志物相比，本发明所述的生物标志物不仅能够表征检测对象的PCB暴露，同是能够反映检测对象的生殖功能变化。

附图说明

图1为Real-TimePCR验证Psap在附睾中的表达。在图1中，横坐标为多氯联苯(μg/kg)，纵坐标为相对表达量；a为cDNA微阵列，b为实时定量PCR；“*”表示与对照组有显著差异。

图2为Western blot检测在不同暴露浓度下Psap在附睾中的表达。

图3为Psap在附睾中的相对表达值。在图3中，横坐标为多氯联苯(μg/kg)，纵坐标为相对表达量；不同上标a，b，c表示相互间有显著差异。

图4为Western blot检测在不同暴露浓度下Psap在精子中的表达。

图5为Psap在精子中的相对表达值。在图5中，横坐标为多氯联苯(μg/kg)，纵坐标为相对表达量；不同上标a，ab，c表示相互间有显著差异。

图6为Western blot检测500μg/KgAroclor 1254暴露不同时间后Psap在精子中的表达。

图7为Psap在对照组和暴露组精子中相对表达的比值。在图7中，横坐标为暴露天数，纵坐标为暴露组/对照组；不同上标a，a，ab，c表示相互间有显著差异。

具体实施方式

以下实施例将结合附图对本发明作进一步的说明。

1.实验动物的暴露

以C57BL/6品系健康青春期雄性小鼠为受试动物，平均体重15g，由厦门大学实验动物中心提供(SPF级)。实验温度控制在24±1℃，昼夜周期为12:12h(早上08:00开灯)，实验期间小鼠自由取食饮水。将雄小鼠随机分组，每组12只。根据环境中PCBs的浓度设计染毒浓度：0.5、5、50和500μg/kg。每3天对小鼠灌胃染毒一次，灌胃剂量为每克体重5μL，对照组给予每克体重5μL生理盐水，暴露持续50天。取小鼠附睾用于步骤7或迅速放入液氮中保存。

用于观察时间效应的雄性C57BL/6品系小鼠随机分成2组，分别暴露于对照溶剂和500μg/kg Aroclor 1254，饲养条件和操作方法如前所述。分别在暴露第7、14、28、50天取样。

选用3个浓度梯度即0，5和500μg/kg暴露组的小鼠样品进行基因表达谱微阵列分析。

2.总RNA提取

2.1提取总RNA

用Invitrogen公司的TRIzol提取附睾总RNA。

2.2总RNA质量检测

1)紫外检测：

检测260nm和OD280nm处的吸光度，OD260/OD280比值为1.8～2.0，260nm处一个OD＝40μg/mL RNA。

2)电泳检测：

1.5％琼脂糖凝胶电泳后，紫外灯下可见清晰明亮的28s、18s条带，28S rRNA条带的亮度应约为18S rRNA的1.5～2.0倍。

3.基因表达谱微阵列样品的制备

根据晶芯

cRNA扩增标记试剂盒的要求完成基因表达谱微阵列样品的制备。

3.1合成1st-strand cDNA

合成1st-strand cDNA，使用晶芯

cRNA扩增标记试剂盒。

1)取总RNA0.1～5μg调整体积到10μl。

2)加入晶芯

基因表达谱外参(Cat.No.360030)1μl。

3)加入T7-Oligo(dT)Primer 1μl，调整体积到12μl。

4)轻轻混匀溶液，70℃反应10min，反应结束后立即放入冰浴骤冷5min。

5)配制反转录Master Mix，每个20μl反应体系包含如下成分：10×First-strand Buffer 2μl，dNTP Mix 4μl，RNase Inhibitor 1μl，CbcScript 1μl，轻轻混匀，短暂离心将溶液收集于管底。上一步骤中反应完的每个样品管中加入上述Master Mix 8μl。

6)轻轻混匀溶液，42℃反应2h。PCR仪的热盖温度42℃，反应完后冰浴5min。

3.2合成2nd-strand cDNA

合成2nd-strand cDNA，使用晶芯cRNA扩增标记试剂盒。反应系统见表1，16℃反应2h。产物置于冰上钝化。

表1合成2nd-strand cDNA的反应体系

3.3纯化cDNA

使用

Extract II(MN公司，Cat.No.740609.250)试剂盒纯化cDNA

3.4cRNA的合成

将纯化后的cDNA真空浓缩至16μl，加入到已准备好的0.2ml离心管中；加入24μl体外转录MasterMix，轻轻混匀，37℃反应4～14h。

体外转录Master Mix包含以下成分：NTP Mix 16μl，T710×Reaction Buffer 4μl，T7 EnzymeMix 4μl。

3.5纯化cRNA

使用

RNA Clean-up(Cat.No.MN-740948.250)试剂盒纯化cRNA

3.6cRNA定量

使用紫外分光光度计进行定量cRNA。1OD260＝40ng/μl。琼脂糖变性胶电泳检测得到cRNA smear片段范围在250～5000nt，平均大小接近1150nt。

3.7cRNA反转录

使用晶芯

cRNA扩增标记试剂盒反转录cRNA。

取cRNA纯化产物2μg，调整体积到7.5μl，加入到0.2ml无核酸酶离心管中。加入4μl Random Primer，混匀，65℃反应10min，冰浴5min。反应完的每个样品管中加入8.5μl二次反转录Master Mix，轻轻混匀，25℃反应10min，37℃反应1.5h。反应结束后加入Terminate Solution 5μl，混匀。65℃反应10min，冰浴5min后加入Neutralize Solution 1μl，混匀。

二次反转录Master Mix包含以下成分：4×CbcScript II Buffer 5μl，0.1M DTT 2μl，CbcScriptII 1.5μl。

3.8cRNA反转录产物纯化和定量

使用

Extract II(MN公司，Cat.No.740609.250)试剂盒纯化反转录产物。使用紫外分光光度计定量纯化后的DNA溶液。1OD260＝40ng/μl。

3.9样品的标记

将纯化后反转录得到的DNA真空浓缩到14μl。加入4μl Random Primer，混匀。95℃反应3min，冰浴5min。依次加入下列试剂：5×Klenow Buffer 5μl，Cy5-dCTP(或Cy3-dCTP)1μl  (GE Healthcare)，Klenow Fragment 1.2μl。混匀，37℃反应1.5h，70℃反应5min，冰浴5min。使用

Extract II(MN公司，Cat.No.740609.250)试剂盒纯化产物。使用紫外分光光度计定量纯化后的DNA溶液。1OD260＝40ng/μl。

4.基因表达谱微阵列样品的杂交、洗脱、检测

标记后的样品用于基因表达谱微阵列样品的杂交。基因表达谱微阵列样品的杂交、洗脱、检测的一切操作均按博奥公司推荐的条件进行。

5.基因表达谱微阵列结果的分析

用微阵列显著性分析算法(Significance Analysis of Microarray，SAM)分析差异表达基因的显著性，仅有通过SAM分析表现出显著性差异(p＜0.05)且变化倍数不小于2倍的数据被用于下游的生物信息学分析。

为了进一步了解Arolcor 1254暴露对小鼠附睾功能的影响，利用分子注释系统(Molecularannotation system，MAS)对表现出显著性差异的基因进行分子注释，检索相关基因的功能，以及基因分类学(Gene Ontology，GO)注释。

以变化倍数≥3为阈值，同时包括GO注释：cell component为extracelluar space或extracelluar region的差异表达基因为候选生物标志物。其中上调的基因9个，下调的基因16个。表2给出部分潜在标志物基因列表，从表2可以看出，这些蛋白包括分泌性的蛋白、多肽、受体以及酶类。

表2部分潜在标志物基因列表

6.实时荧光定量PCR(Real-time PCR)

6.1cDNA模板的合成

总RNA的提取和检测参照步骤3。将提取的总RNA用TaKaRa 1^st srand合成试剂盒合成cDNA模板。

6.2标准品的制备

以上述cDNA模板为模板加入引物对，通过聚合酶链式反应扩增，其中使用的引物对是：

引物1：5’-GCCAGAGGGCAGGAGCATT-3’

引物2：5’-CTGACCCAGGGACAGCAACA-3’

将扩增完毕的PCR产物进行1.5％的琼脂糖凝胶电泳，将PCR产物条带切下，用PCR参悟纯化柱提取纯化PCR产物，用紫外分光光度法，计测260nm出OD值。

6.3实时荧光定量聚合酶链式反应Real-time PCR检测

反应总体积为25μl，反应体系含10倍系列稀释的标准品或同法制备的样品2μl，浓度为10pmol/μl的引物个1μl，Sybr Green realtime PCR master mix(Toyobo)12.5μl，用灭菌双蒸水补足25μl。循环条件为：95℃预变性2min；95℃，15s变性，60℃，30s退火，72℃，30s延伸，扩增40个循环。在每个循环的延伸步骤进行荧光检测，分别得到标准品和样品的扩增曲线和Ct值。

6.4结果的计算与分析

利用样品扩增后的Ct值在标准曲线上对应起始拷贝数，数值用REST软件得出。结果表明Psap的mRNA水平在PCB暴露后表现出显著的下降，与微阵列的结果表现出一致性(图1)。

7.精子参数检测

新鲜取得的附睾尾部剪开，使用Chemicon公司HTF培养基(MR-070-D)37℃孵育30min，让精子充分游出。迅速吸取10μl新鲜精子悬液于载玻片上，40倍物镜下观察精子活力，统计活动精子的总百分率。另取10μl精子悬液制备精子涂片，固定后1％伊红染色，统计每200个精子中畸形精子的百分比。

8.Psap蛋白表达的检测

使用RIPA提取样品组织的总蛋白，使用Thermo Scientific BCA蛋白定量分析试剂盒对总蛋白进行定量。用Western blot法检测Psap蛋白在附睾中的表达。附睾蛋白提取物μg，用10％SDS PAGE电泳分离后电转移到PVDF膜上，5％脱脂牛奶封闭后，加入兔抗鼠Psap多克隆抗体(工作浓度1∶2000)，4℃过夜，加入HRP标记的羊抗兔抗体(工作浓度1∶5000)与37℃反应一h，ECL显色后观察结果。

图2、图3显示PCB暴露后小鼠附睾Psap蛋白的表达随着PCB暴露剂量变化的趋势。结果显示附睾中Psap蛋白的表达表现出剂量依赖性的下调。进一步考察Psap在精子中的表达，发现Psap在精子中的表达也呈现出剂量依赖的下降(图4、图5)。

图6、图7显示PCB暴露对Psap在精子中表达的时间效应。随着暴露时间的延长，精子中的Psap表达表现出递减的趋势，在28天以后表现出显著的下降。

9.PCB残留的检测

9.1微波萃取

称取1.0g肝脏样品，加入PCB 103为内标，微波辅助萃取。在微波功率300W，350psi及100℃的条件下，丙酮/正己烷萃取15min。

9.2净化与浓缩

萃取液经无水硫酸钠后，转移至吹脱管中进行浓缩。氮吹浓缩。浓缩液依次通过中性氧化铝、活化硅胶、无水硫酸钠，丙酮：正己烷淋洗，氮吹近干后，用丙酮：正己烷定容至1mL。

9.3气相色谱分析检测

取上述定容样品1μ进样，按照下列色谱条件进行检测。

载气：高纯N2(＞99.999％)；进样口温度：280。℃；检测器(ECD)温度：300。℃；柱压力：68、95kPa；尾吹60mL/min。升温程序：初始温度85℃，保持2m in，再以30℃/min升温至180℃，保持2min；然后以20℃/min升温至200℃，保持3min。然后以3℃/min升温至230℃，保持10min；然后以5℃/min升温至250℃，保持8min。无分流进样。获得PCB的气相色谱图，对小峰和干扰峰实行积分抑制，最终选取11个与基底不重合且峰型较明显的特征峰来表征整体的浓度

10.Psap表达与多氯联苯暴露对雄性哺乳动物产生生殖毒性的相关性分析

效应标志物Psap、体内PCB水平以及精液指标之间的相关性采用Spearman相关分析进行统计，结果见表3。结果显示精子中Psap表达量与体内PCB水平呈显著的负相关，Psap表达量与精子活力和精子畸形率均有显著的相关性。

表3Psap表达与肝脏PCB浓度、精子活力和畸形率的Spearman相关系数。

“*”表示相关性显著(p＜0.05)；“**”表示相关性极显著(p＜0.01)

Claims (1)

1.鞘脂激活蛋白原在制备多氯联苯雄性生殖毒性检测效应标志物中的应用，所述鞘脂激活蛋白原是一种多功能的糖蛋白，分子量为70kD，存在于一些生物体液中；

所述鞘脂激活蛋白原采用以下方法筛选：

1）用不同剂量的多氯联苯对小鼠进行暴露后收集样本；所述不同剂量的多氯联苯采用以下剂量：0.5μg/kg、5μg/kg、50μg/kg、500μg/kg；所述收集样本，在暴露50天后收集样本；所述样本为附睾和肝脏样本；

2）采用基因芯片技术进行样本差异表达基因的筛选，具体方法如下：

对步骤1）所收集的正常对照样本不同剂量的多氯联苯处理的附睾样本提取的总RNA，以RNA为模板合成cDNA；利用获得的cDNA合成cRNA，片段化荧光标记的cRNA；标记的cRNA经片段化处理后用于基因芯片杂交，基因芯片进行杂交、洗脱、染色及检测；对获得的检测数据进行归一化处理，使不同的样本数据具有可比性；

3）通过生物信息学统计和分析软件，实现对筛选的差异表达基因进行功能聚类，获得与生殖功能相关的差异表达基因，具体方法如下：

对步骤2）获得的数据进行统计学分析，利用微阵列显著性分析算法分析差异表达基因的显著性，利用分子注释系统数据库，对差异表达基因以进行检索，并建立PCB暴露于附睾功能变化间的相互关系；

4）利用实时荧光定量PCR技术和免疫印迹法在mRNA转录水平和蛋白表达水平上对获得的差异表达基因随着PCB暴露剂量以及暴露时间的增加而变化进行检测分析，具体方法如下：

将步骤2）提取的总RNA利用特异性引物进行实时荧光PCR，定量Psap的mRNA转录水平，所述特异性引物具有如下序列：

引物1：5’-GCCAGAGGGCAGGAGCATT-3’

引物2：5’-CTGACCCAGGGACAGCAACA-3’

样品在缓冲液中匀浆后，依次进行聚丙烯酰胺凝胶电泳和电转移，利用特异性抗体检测Psap的蛋白丰度；

所述特异性抗体为鞘脂激活蛋白原多克隆抗体；

5）进行精子参数的检测分析，具体方法如下：

将步骤1）取得的附睾尾部剪开，使用Chemicon公司HTF培养基37℃孵育30min，让精子游出，吸取10μl新鲜精子悬液于载玻片上，40倍物镜下观察精子活力，统计活动精子的总百分率，另取10μl精子悬液制备精子涂片，固定后1%伊红染色，统计每200个精子中畸形精子的百分比；

6）对效应标志物、体内多氯联苯水平以及精液指标之间的相关系数通过斯皮尔曼相关分析进行统计，具体方法如下：通过气象色谱-电子捕获检测器检测步骤1）取得肝脏中PCB的富集，将PCB的富集与步骤4）的Psap蛋白丰度，以及步骤5）的精子参数进行斯皮尔曼相关性分析。

Images