通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法

摘要

本发明公开了一种通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，包括培育SD大鼠并获取睾丸、测定T水平、提取总RNA纯化合成cDNA并进行qPCR反应、获取总蛋白样品进行Western蛋白质条带密度分析、睾丸切片并染色、转录组测序分析、分析差异表达基因(DEGs)功能、计算并比较不同睾丸细胞类型的基因百分比等步骤，通过该方法可以对双酚A对睾丸间质细胞发育影响的机理进行研究判断，同时步骤简单高效。

权利要求

1.一种通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，其特征在于：包括以下步骤，

步骤一，将SD大鼠置于23±2℃的恒温、55±5%的恒湿环境中，进行12小时的光/暗循环，并允许SD大鼠自由获取食物和水，将BPA溶解在乙醇中，然后稀释在生育酚剥脱玉米油中，直至乙醇的最终浓度为0.5%，用于胃内给药；将SD大鼠随机分为三组，分别设置为对照组、低剂量组和高剂量组，每组8只，从出生后第21天开始，每天分别灌胃0、0.1或10 mg/kg BW剂量的BPA，持续14天；在出生第56天处死SD大鼠；收集躯干血，在4℃下以250g的转速离心15min，以收集血清，并在‑80℃下储存；将其中一个睾丸浸入液氮中冷冻，并在‑80℃下保存；将另一个睾丸固定在4%多聚甲醛中；

步骤二，使用IMMULITE®2000试剂盒，通过免疫化学发光法测定血清T水平；采用酶联免疫吸附试验（ELISA）试剂盒测定血清黄体生成素（LH）、卵泡刺激素（FSH）、双氢睾酮（DHT）、生长激素（GH）和胰岛素样生长因子1（IGF‑1）的浓度；T、LH、FSH、DHT、GH和IGF‑1试剂盒的灵敏度分别为0.15 ng/ml、0.2 mIU/ml、0.1 mIU/ml、0.5 pg/ml、0.1 ng/ml和3 ng/ml；

步骤三，通过定量PCR（qPCR）检测睾丸间质细胞类固醇合成基因和支持细胞标记基因的转录水平，睾丸间质细胞类固醇合成基因包括Cyp11a1, Cyp17a1, Star, Hsd3b1, Hsd17b3, Scarb1, Lhcgr, Akr1c14, Srd5a1，支持细胞标记基因为Sox9和 Dhh，使用RNeasy®微型试剂盒从睾丸中提取总RNA，使用NanoDrop ND‑1000测量RNA的数量和纯度，然后通过逆转录系统合成cDNA；在含有7.5ul LightCycler®480 SYBR®Green I Master、

0.75ul反向引物、0.75ul正向引物、10ul总量为1ng的 RNA和4ul去酶水的容器中进行PCR反应；将核糖体蛋白S16（RPS16）作为qPCR结果的标准内参物；

步骤四，用RIPA裂解缓冲液对睾丸组织进行均质和裂解获得总蛋白样品，使用增强型BCA蛋白质分析试剂盒测量蛋白质浓度，取60µg总蛋白样品，在含有10%十二烷基硫酸钠的聚丙烯酰胺凝胶中进行电泳分离，然后转移到聚偏氟乙烯膜上；用含5%BSA的TBST封闭2h后，在4℃下将聚偏氟乙烯膜与一抗一起培育过夜，所述一抗是针对于CYP11A1, CYP17A1, STAR, SCARB1, LHCGR, HSD3B1, HSD17B3, INSL3, SOX9, DHH, SRD5A1, IRS1，p‑IRS1，AKT1，p‑AKT1，CREB，p‑CREB及b‑ACTIN的抗体；然后，将聚偏氟乙烯膜与山羊抗小鼠IgG（H+L）HRP结合物或山羊抗兔IgG（H+L）在室温下培育2小时；通过Image J 软件进行蛋白质条带密度分析，将β‑肌动蛋白作为内部对照值；

步骤五，选择CYP11A1作为识别睾丸间质细胞的生物标记物，先将睾丸固定，在乙醇和二甲苯中梯度脱水，将各SD大鼠睾丸切成4部分后包埋于石蜡块中，然后在每个睾丸中随机取相同大小的部分作为组织阵列包埋在石蜡中，切割获得横截面厚度为5微米的切片，并将切片粘附在载玻片上；完成脱蜡和复水后，用含10%H2O2的甲醇孵育10min以阻断内源性过氧化物酶，用PBS冲洗切片；切片在PH6.0的0.01M柠檬酸钠中加热，并在室温下冷却；用二抗宿主血清进行封闭后，切片在4℃下与CYP11A1兔克隆抗体孵育过夜，用PBS洗涤，然后在室温下与Dylight 488山羊抗兔IgG（H+L）孵育1小时；在用PBS洗涤后，向切片中添加含有DAPI的抗荧光淬灭剂，以显示胞质为绿色的CYP11A1阳性细胞和蓝色细胞核；

步骤六，从步骤三所得部分总RNA，委托LC Bio Technology Co对对照组和低剂量组的睾丸RNA样本进行RNA‑seq测序；通过生物分析仪2100对完整性RIN值>7.0的RNA进行评估，并通过变性琼脂糖凝胶电泳进行确认；将片段化的RNA通过逆转录酶的作用下合成cDNA；序列由Illumina NovaseqTM 6000平台完成匹配，采用R‑package‑edgeR软件包筛选差异倍数>1.5，P值<0.05的差异表达基因（DEGs）；采用基因本体（GO）、京都基因和基因组百科全书（KEGG）分析DEGs功能；

步骤七，将10种不同睾丸细胞类型的100个最特异标记基因与DEGs对照，睾丸细胞类型为精原细胞，精母细胞，精子细胞，支持细胞，管周myoid细胞，Leydig细胞，血管内皮细胞，间充质细胞，巨噬细胞和淋巴细胞；计算并比较10种不同睾丸细胞类型的匹配基因百分比。

2.根据权利要求1所述的通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，其特征在于：所述步骤五中，在用含10%H2O2的甲醇阻断内源性过氧化物酶后，用PBS冲洗切片3次，每次持续5min。

3.根据权利要求1所述的通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，其特征在于：所述步骤五中，切片在4℃下与CYP11A1兔克隆抗体培育过夜后，用PBS洗涤切片3次，每次5分钟。

4.根据权利要求1所述的通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，其特征在于：所述步骤五中，还包括睾丸间质细胞的统计，在20X物镜下对每个区域中的CYP11A1阳性细胞进行计数，在至少四只SD大鼠的睾丸切片上统计多个区域。

说明书

通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测

试方法

技术领域

[0001] 本发明涉及生物医药技术领域。

背景技术

[0002] 1981年，双酚A，即2,2‑双（4‑羟基苯基）丙烷，缩写为BPA，首次作为聚碳酸酯塑料制造单体合成。从那时起，它已被广泛用于各种消费品，包括与食品密切接触的产品，如塑料婴儿奶瓶、食品罐内衬和牙科修复材料。除了来自一般环境，也就是水、空气和土壤的潜在暴露外，BPA也可能直接通过受污染的食品进入人体。调查表明，90%的美国普通人群的尿液样本中可检测到BPA，平均血清浓度为2.5‑10.95 ng/ml。同时最近的一项研究发现，中国3至11岁儿童尿液中的总双酚含量在0.43至31.5 ng/ml之间。

[0003] 目前有研究发现，双酚A被认为是一种内分泌干扰物，其作用类似于雌激素，并具有干扰雄激素合成的作用。双酚A对睾丸间质细胞发育的影响似乎取决于剂量和暴露时间。例如，子宫内暴露于高剂量（40‑400mg/kg BW/day）的双酚A会抑制胎儿睾丸间质细胞和支持细胞的分化，并导致大鼠体内睾酮（T）的产生降低。然而，在新生小鼠第0‑9天使用剂量高达97mg/kg BW的BPA不会影响成年的生殖功能，而在小鼠青春期，即出生后第21‑35天，使用极低剂量（2.4μg/kg BW/day）的BPA暴露会降低35天时的血清LH和T水平。此外，从妊娠第12天到出生后第21天的双酚A暴露，成年后睾丸间质液中的T水平降低，这意味着围产期是BPA暴露的敏感时期，成年大鼠直接暴露于双酚A 两周也会降低睾丸重量和T水平，这意味着双酚A不仅会影响睾丸间质细胞的发育，还会影响其功能。尽管如此，大多数研究使用相对较高剂量的双酚A（1‑100毫克/千克体重/天），但是不清楚双酚A作用的潜在机制，也没有很好的检测方法进行检测，并且睾丸间质细胞对于男性生殖系统的发育和男性生殖功能的维持至关重要。在大鼠中，成年睾丸间质细胞（ALC）由青春期的睾丸间质干细胞（SLC）发育而来，经过4个发育阶段：SLC、PLC（祖细胞LC）、ILC（未成熟LC）和ALC，成年鼠每个睾丸中总共形成了约2500万个ALC。出生后21‑35天，是PLC向ILC发展的时期，也是ALC发育中最关键的阶段之一，也是BPA暴露最敏感的阶段。然而，在这一阶段，BPA暴露的长期影响以及研究机理还不清楚，目前没有很好的判断研究方式。

发明内容

[0004] 针对现有技术的不足，本发明提供了一种通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，通过该方法可以对双酚A对睾丸间质细胞发育影响的机理进行研究判断，同时步骤简单高效。

[0005] 为实现上述目的，本发明提供了一种通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法，包括以下步骤，

[0006] 步骤一，将SD大鼠置于23±2℃的恒温、55±5%的恒湿环境中，进行12小时的光/暗循环，并允许SD大鼠自由获取食物和水，将BPA溶解在乙醇中，然后稀释在生育酚剥脱玉米油中，直至乙醇的最终浓度为0.5%，用于胃内给药；将SD大鼠随机分为三组，分别设置为对照组、低剂量组和高剂量组，每组8只，从出生后第21天开始，每天分别灌胃0、0.1或10 mg/kg BW剂量的BPA，持续14天；在出生第56天处死SD大鼠；收集干血，在4℃下以250g的转速离心15min，以收集血清，并在‑80℃下储存；将其中一个睾丸浸入液氮中冷冻，并在‑80℃下保存；将另一个睾丸固定在4%多聚甲醛中；

[0007] 步骤二，使用IMMULITE®2000试剂盒，通过免疫化学发光法测定血清T水平；采用酶联免疫吸附试验（ELISA）试剂盒测定血清黄体生成素（LH）、卵泡刺激素（FSH）、双氢睾酮（DHT）、生长激素（GH）和胰岛素样生长因子1（IGF‑1）的浓度；T、LH、FSH、DHT、GH和IGF‑1试剂盒的灵敏度分别为0.15 ng/ml、0.2 mIU/ml、0.1 mIU/ml、0.5 pg/ml、0.1 ng/ml和3 ng/ml；

[0008] 步骤三，通过定量PCR（qPCR）检测睾丸间质细胞类固醇合成基因和支持细胞标记基因的转录水平，使用RNeasy®微型试剂盒从睾丸中提取总RNA，使用NanoDrop ND‑1000测量RNA的数量和纯度，然后通过逆转录系统合成cDNA；在含有7.5ul LightCycler®480 SYBR®Green I Master、0.75ul反向引物、0.75ul正向引物、10ul总量为1ng的 RNA和4ul去酶水的容器中进行PCR反应；将核糖体蛋白S16（RPS16）作为qPCR结果的标准内参物；

[0009] 步骤四，用RIPA裂解缓冲液对睾丸组织进行均质和裂解获得总蛋白样品，使用增强型BCA蛋白质分析试剂盒测量蛋白质浓度，取60µg总蛋白样品，在含有10%十二烷基硫酸钠的聚丙烯酰胺凝胶中进行电泳分离，然后转移到聚偏氟乙烯膜上；用含5%BSA的TBST封闭2h后，在4℃下将聚偏氟乙烯膜与一抗一起培育过夜；然后，将聚偏氟乙烯膜与山羊抗小鼠IgG（H+L）HRP结合物或山羊抗兔IgG（H+L）在室温下培育2小时；通过Image J 软件进行蛋白质条带密度分析，将β‑肌动蛋白作为内部对照值；

[0010] 步骤五，选择CYP11A1作为识别睾丸间质细胞的生物标记物，先将睾丸固定，在乙醇和二甲苯中梯度脱水，将各SD大鼠睾丸切成4部分后包埋于石蜡块中，然后在每个睾丸中随机取相同大小的部分作为组织阵列包埋在石蜡中，切割获得横截面厚度为5微米的切片，并将切片粘附在载玻片上；完成脱蜡和复水后，用含10%H2O2的甲醇孵育10min以阻断内源性过氧化物酶，用PBS冲洗切片；切片在PH6.0的0.01M柠檬酸钠中加热，并在室温下冷却；用二抗宿主血清进行封闭后，切片在4℃下与CYP11A1兔克隆抗体孵育过夜，用PBS洗涤，然后在室温下与Dylight 488山羊抗兔IgG（H+L）孵育1小时；在用PBS洗涤后，向切片中添加含有DAPI的抗荧光淬灭剂，以显示胞质为绿色的CYP11A1阳性细胞和蓝色细胞核；

[0011] 步骤六，采用LC Bio Technology Co对对照组和低剂量组的睾丸RNA样本进行RNA‑seq分析；通过生物分析仪2100对完整性RIN值>7.0的RNA进行评估，并通过变性琼脂糖凝胶电泳进行确认；将片段化的RNA通过逆转录酶的作用下合成cDNA；序列由Illumina NovaseqTM 6000平台完成匹配，采用R‑package‑edgeR软件包筛选差异倍数>1.5，P值<0.05的差异表达基因（DEGs）；采用基因本体（GO）、京都基因和基因组百科全书（KEGG）分析DEGs功能；

[0012] 步骤七，将10种不同睾丸细胞类型的100个最特异标记基因与DEGs对照；计算并比较10种不同睾丸细胞类型的匹配基因百分比。

[0013] 作为本发明的进一步设置，所述步骤三中，在用含10%H2O2的甲醇阻断内源性过氧化物酶后，用PBS冲洗切片3次，每次持续5min。

[0014] 作为本发明的进一步设置，所述步骤三中，切片在4℃下与CYP11A1兔克隆抗体培育过夜后，用PBS洗涤切片3次，每次5分钟。

[0015] 作为本发明的进一步设置，所述步骤三中还包括睾丸间质细胞的统计，在20X物镜下对每个区域中的CYP11A1阳性细胞进行计数，在至少四只SD大鼠的睾丸切片上统计多个区域。

[0016] 这样设置的有益效果是：这样设置的检测方法，通过BPA的暴露重置了成年SD大鼠的生长激素（GH）和卵泡刺激素（FSH）内稳态，并导致血清睾酮水平显著升高，但不影响血清促黄体生成素水平。在qPCR和Western bolt的实验显示BPA显著上调睾丸间质细胞类固醇合成途径中的选择性基因/蛋白质，包括类固醇生成急性调节蛋白、细胞色素P450 11A1、细胞色素P450 17A和低密度脂蛋白受体，通过睾丸的RNA‑Seq分析表明BPA显著影响睾丸间质体细胞群体的基因，但不影响生殖细胞。此外，GO和KEGG分析表明与IGF‑1相关的PI3K/AKT信号被激活，IRS1、AKT1和CREB的磷酸化增加证实了这一点。通过这种测试方法，可以得出青春期短期低剂量的BPA暴露会影响成年雄性SD大鼠垂体（GH和FSH）和睾丸（睾酮）的稳态，并有效验证BPA对成年SD大鼠的作用机理。

附图说明

[0017] 图1A为本发明实施例中不同浓度双酚A对血清睾酮激素水平影响的示意图；

[0018] 图1B为本发明实施例中不同浓度双酚A对血清DHT影响的示意图；

[0019] 图1C为本发明实施例中不同浓度双酚A对血清LH影响的示意图；

[0020] 图1D为本发明实施例中不同浓度双酚A对血清FSH影响的示意图；

[0021] 图1E为本发明实施例中不同浓度双酚A对血清GH影响的示意图；

[0022] 图1F为本发明实施例中不同浓度双酚A对血清IGF‑1影响的示意图；

[0023] 图1G为本发明实施例中不同浓度双酚A对睾丸间质细胞数量影响的示意图；

[0024] 图2A为本发明实施例中不同浓度双酚A对Cyp11a1基因表达的影响示意图；

[0025] 图2B为本发明实施例中不同浓度双酚A对Cyp17a1基因表达的影响示意图；

[0026] 图2C为本发明实施例中不同浓度双酚A对Star基因表达的影响示意图；

[0027] 图2D为本发明实施例中不同浓度双酚A对Hsd3b1基因表达的影响示意图；

[0028] 图2E为本发明实施例中不同浓度双酚A对Hsd17b3基因表达的影响示意图；

[0029] 图2F为本发明实施例中不同浓度双酚A对Scarb1基因表达的影响示意图；

[0030] 图2G为本发明实施例中不同浓度双酚A对Lhcgr基因表达的影响示意图；

[0031] 图2H为本发明实施例中不同浓度双酚A对Akr1c14基因表达的影响示意图；

[0032] 图2I为本发明实施例中不同浓度双酚A对Srd5a1基因表达的影响示意图；

[0033] 图2J为本发明实施例中不同浓度双酚A对Sox9基因表达的影响示意图；

[0034] 图2K为本发明实施例中不同浓度双酚A对Dhh基因表达的影响示意图；

[0035] 图3A为本发明实施例中每种蛋白质在不同浓度双酚A下的代表性蛋白质印迹示意图；

[0036] 图3B为本发明实施例中CYP11A1在不同浓度双酚A下的量化示意图；

[0037] 图3C为本发明实施例中CYP17A1在不同浓度双酚A下的量化示意图；

[0038] 图3D为本发明实施例中STAR在不同浓度双酚A下的量化示意图；

[0039] 图3E为本发明实施例中SCARB1在不同浓度双酚A下的量化示意图；

[0040] 图3F为本发明实施例中LHCGR在不同浓度双酚A下的量化示意图；

[0041] 图3G为本发明实施例中HSD3B1在不同浓度双酚A下的量化示意图；

[0042] 图3H为本发明实施例中HSD17B3在不同浓度双酚A下的量化示意图；

[0043] 图3I为本发明实施例中INSL3在不同浓度双酚A下的量化示意图；

[0044] 图3J为本发明实施例中SOX9在不同浓度双酚A下的量化示意图；

[0045] 图3K为本发明实施例中DHH在不同浓度双酚A下的量化示意图；

[0046] 图3L为本发明实施例中SRD5A1在不同浓度双酚A下的量化示意图；

[0047] 图4A为本发明实施例中6个样本的PCA主成分分析图；

[0048] 图4B为本发明实施例中差异表达基因的热图；

[0049] 图4C为本发明实施例中差异表达基因的火山图；

[0050] 图4D为本发明实施例中585个DEGs与之前报道的睾丸细胞特异性标记基因的匹配列表图；

[0051] 图5A为本发明实施例中上调的差异表达基因的GO分析图；

[0052] 图5B为本发明实施例中上调的差异表达基因的KEGG分析图；

[0053] 图6为本发明实施例中的类固醇合成相关的信号通路图；

[0054] 图7A为本发明实施例中每种蛋白质在不同浓度双酚A下的代表性蛋白质印记表达图；

[0055] 图7B为本发明实施例中IRS1蛋白在不同浓度双酚A下的Western blot结果量化图；

[0056] 图7C为本发明实施例中AKT1蛋白在不同浓度双酚A下的Western blot结果量化图；

[0057] 图7D为本发明实施例中CREB蛋白在不同浓度双酚A下的Western blot结果量化图；

[0058] 图7E为本发明实施例中p‑IRS1蛋白在不同浓度双酚A下的Western blot结果量化图；

[0059] 图7F为本发明实施例中p‑AKT1蛋白在不同浓度双酚A下的Western blot结果量化图；

[0060] 图7G为本发明实施例中p‑CREB蛋白在不同浓度双酚A下的Western blot结果量化图；

[0061] 图7H为本发明实施例中p‑IRS1/IRS1在不同浓度双酚A下的Western blot结果量化图；

[0062] 图7I为本发明实施例中p‑AKT1/AKT1在不同浓度双酚A下的Western blot结果量化图；

[0063] 图7J为本发明实施例中p‑CREB/CREB在不同浓度双酚A下的Western blot结果量化图；

具体实施方式

[0064] 本发明通过双酚A暴露重置脑垂体稳态增加成年大鼠血清睾酮的测试方法的实施例：包括以下步骤，

[0065] 步骤一，将SD大鼠置于23±2℃的恒温、55±5%的恒湿环境中，进行12小时的光/暗循环，并允许SD大鼠自由获取食物和水，将BPA溶解在乙醇中，然后稀释在生育酚剥脱玉米油中，直至乙醇的最终浓度为0.5%，用于胃内给药；将SD大鼠随机分为三组，分别设置为对照组、低剂量组和高剂量组，每组8只，从出生后第21天开始，每天分别灌胃0、0.1或10 mg/kg BW剂量的BPA，持续14天；在出生第56天处死SD大鼠；收集干血，在4℃下以250g的转速离心15min，以收集血清，并在‑80℃下储存；将其中一个睾丸浸入液氮中冷冻，并在‑80℃下保存；将另一个睾丸固定在4%多聚甲醛中；

[0066] 步骤二，使用IMMULITE®2000试剂盒，通过免疫化学发光法测定血清T水平；采用酶联免疫吸附试验（ELISA）试剂盒测定血清黄体生成素（LH）、卵泡刺激素（FSH）、双氢睾酮（DHT）、生长激素（GH）和胰岛素样生长因子1（IGF‑1）的浓度；T、LH、FSH、DHT、GH和IGF‑1试剂盒的灵敏度分别为0.15 ng/ml、0.2 mIU/ml、0.1 mIU/ml、0.5 pg/ml、0.1 ng/ml和3 ng/ml；

[0067] 步骤三，通过定量PCR（qPCR）检测睾丸间质细胞类固醇合成基因和支持细胞标记基因的转录水平，使用RNeasy®微型试剂盒从睾丸中提取总RNA，使用NanoDrop ND‑1000测量RNA的数量和纯度，然后通过逆转录系统合成cDNA；在含有7.5ul LightCycler®480 SYBR®Green I Master、0.75ul反向引物、0.75ul正向引物、10ul总量为1ng的 RNA和4ul去酶水的容器中进行PCR反应；将核糖体蛋白S16（RPS16）作为qPCR结果的标准内参物；

[0068] 步骤四，用RIPA裂解缓冲液对睾丸组织进行均质和裂解获得总蛋白样品，使用增强型BCA蛋白质分析试剂盒测量蛋白质浓度，取60µg总蛋白样品，在含有10%十二烷基硫酸钠的聚丙烯酰胺凝胶中进行电泳分离，然后转移到聚偏氟乙烯膜上；用含5%BSA的TBST封闭2h后，在4℃下将聚偏氟乙烯膜与一抗一起培育过夜；然后，将聚偏氟乙烯膜与山羊抗小鼠IgG（H+L）HRP结合物或山羊抗兔IgG（H+L）在室温下培育2小时；通过Image J 软件进行蛋白质条带密度分析，将β‑肌动蛋白作为内部对照值；

[0069] 步骤五，选择CYP11A1作为识别睾丸间质细胞的生物标记物，先将睾丸固定，在乙醇和二甲苯中梯度脱水，将各SD大鼠睾丸切成4部分后包埋于石蜡块中，然后在每个睾丸中随机取相同大小的部分作为组织阵列包埋在石蜡中，切割获得横截面厚度为5微米的切片，并将切片粘附在载玻片上；完成脱蜡和复水后，用含10%H2O2的甲醇孵育10min以阻断内源性过氧化物酶，用PBS冲洗切片；切片在PH6.0的0.01M柠檬酸钠中加热，并在室温下冷却；用二抗宿主血清进行封闭后，切片在4℃下与CYP11A1兔克隆抗体孵育过夜，用PBS洗涤，然后在室温下与Dylight 488山羊抗兔IgG（H+L）孵育1小时；在用PBS洗涤后，向切片中添加含有DAPI的抗荧光淬灭剂，以显示胞质为绿色的CYP11A1阳性细胞和蓝色细胞核；

[0070] 步骤六，采用LC Bio Technology Co对对照组和低剂量组的睾丸RNA样本进行RNA‑seq分析；通过生物分析仪2100对完整性RIN值>7.0的RNA进行评估，并通过变性琼脂糖凝胶电泳进行确认；将片段化的RNA通过逆转录酶的作用下合成cDNA；序列由Illumina NovaseqTM 6000平台完成匹配，采用R‑package‑edgeR软件包筛选差异倍数>1.5，P值<0.05的差异表达基因（DEGs）；采用基因本体（GO）、京都基因和基因组百科全书（KEGG）分析DEGs功能；

[0071] 步骤七，将10种不同睾丸细胞类型的100个最特异标记基因与DEGs对照；计算并比较10种不同睾丸细胞类型的匹配基因百分比。

[0072] 作为本发明的进一步设置，所述步骤三中，在用含10%H2O2的甲醇阻断内源性过氧化物酶后，用PBS冲洗切片3次，每次持续5min。

[0073] 作为本发明的进一步设置，所述步骤三中，切片在4℃下与CYP11A1兔克隆抗体培育过夜后，用PBS洗涤切片3次，每次5分钟。

[0074] 作为本发明的进一步设置，所述步骤三中还包括睾丸间质细胞的统计，在20X物镜下对每个区域中的CYP11A1阳性细胞进行计数，在至少四只SD大鼠的睾丸切片上统计多个区域。

[0075] 为了研究双酚A对动物整体生长的影响，按35和56天的年龄记录体重,如下表S1。

[0076]

[0077] Mean ± SEM, n=6, *Significantly different from 0.1mg/kg dose group at P<0.05.

[0078] 在暴露于BPA的后第35天，体重没有发现差异，而在低剂量BPA组中，成年SD大鼠在第56天的体重略有轻微增加，增加了约8%。高剂量组体重明显低于低剂量组，P＜0.05。除此之外，在整个实验过程中，SD大鼠都是健康的，没有其他异常体征。此外，睾丸重量没有发现差异。

[0079] 为了评估在青春期暴露于双酚A对成年SD大鼠生殖功能的潜在影响，分析了主要生殖激素的血清水平，如图1所示。如图1A可以看出0.1 mg/kg剂量的BPA显著增加了T水平,但是对其直接代谢产物DHT水平没有影响,如图1B所示；进一步的，检查T水平的升高是否与促性腺激素的变化有关，实验方法中测量了血清LH和FSH水平，从图1C‑D可以高低两种剂量的BPA都能显著降低FSH，而任何剂量的BPA都不会显著影响LH，这表明血清T的增加不太可能是由于垂体促性腺激素的增加，与之相反的，LH的轻度和不显著降低是由于血清T增加的反馈抑制；除生殖激素外，如图1E所示，低剂量BPA也增加了血清GH水平，而不影响IGF‑1水平，见图1F。

[0080] 进一步的，CYP11A1是鉴定成年SD大鼠睾丸间质细胞的最佳标记物之一，通过CYP11A1染色可以定量ALC的数量。通过测试发现，双酚A处理对成年SD大鼠睾丸中CYP11A1阳性间质细胞数量没有显著影响。

[0081] 进一步的，通过试验操作，认为血清T水平升高表示睾丸间质细胞类固醇合成功能可能发生改变，为了证实在青春期双酚A的暴露是否确实影响成年SD大鼠睾丸间质细胞和支持细胞的功能，通过qPCR分析了睾丸间质细胞类固醇合成基因，包括Cyp11a1、Cyp17a1、Star、Hsd3b1、Hsd17b3、Lhcgr、Scarb1、Akr1c14和Srd5a1，以及支持细胞的标记基因Dhh和Sox9，如图2所示。BPA在0.1mg/kg剂量下特异性上调Cyp11a1、Cyp17a1和Star的表达水平，如图2A‑K所示，在低剂量下显著下调Sox9的表达水平，而不影响其他基因。通过测试可以得出，高水平血清T是由于一些关键的类固醇合成蛋白的特异性变化所致的。

[0082] 为了进一步研究BPA如何影响类固醇生成，通过Western blot分析了主要的类固醇合成蛋白，如图3所示，与mRNA水平的变化一致，0.1 mg/kg剂量的BPA也显著增加了CYP11A1、CYP17A1和STAR的蛋白质水平；其中，STAR是受影响最显著的蛋白质，此外，其他的类固醇生成蛋白与RNA的变化水平一致，通过测试还可以表明，两种剂量的BPA都显著降低了SOX9蛋白水平，而不影响DHH。进一步的，采用RNA‑Seq技术分析了对照组和低剂量BPA组的睾丸总RNA，以研究BPA暴露如何影响睾丸细胞的整体转录。两组每个样本均检测到约32k个独特转录本；主成分分析（PCA）显示，来自对照组和BPA组的样本分别聚集成2个可区分的组，如图4A，表明低剂量BPA显著影响整个睾丸细胞转录组，以P<0.05和比率>1.5倍为筛选标准，试验中鉴定出了585个差异表达基因（DEGs），其中，512个上调，73个下调，可以从热图（图4B）和火山图（图4C）中清楚地看到。S2和S3列出了前50个上调和下调DEGs。

[0083]

[0084] Supplemental Material Table S3: Top 50 Differentially down‑regulated genes by BPA

[0085]

[0086] 在该方法中，由于RNA‑seq是用整个睾丸的RNA进行的检测，因此很难知道哪些睾丸细胞负责检测到的DEGs，为了确定可能发生差异表达的特定细胞类型，通过现有技术，即单细胞RNA‑Seq研究（Green等人，2018年），将确定的10种睾丸细胞类型的前100个标记基因与当前研究检测到的DEGs进行对比，如图4D，大多数DEGs与间质室的体细胞群相匹配，几乎没有与生精小管内的细胞相匹配。与每种细胞类型的前100个标记基因相匹配的DEGs数量分别为：睾丸间质细胞26个、间充质细胞25个、血管内皮细胞23个、淋巴细胞12个、巨噬细胞16个和平滑肌/肌样细胞15个，很少或没有基因与支持细胞、精原细胞、精母细胞以及精子细胞匹配，他们的数量仅为支持细胞4个、精原细胞1个、精母细胞0个或精子细胞0个。从而进一步论证，BPA特异性地影响了间质室的睾丸体细胞群。

[0087] 进一步的，通过该方法，对上调的差异表达基因的GO分析显示，DEGs在细胞外空间、类固醇生物合成过程、IGF相关生物过程、细胞周期调节和内皮细胞形态等过程发生了富集，由图5A可见。其中IGF相关过程出现4次，表明其重要性。通过DEGs的KEGG分析可以看出类固醇生成、自身免疫性甲状腺疾病、皮质醇合成、胆固醇代谢、Jak/STAT信号通路和PI3K/Akt信号通路的富集，如图5B。

[0088] 总的来说，GO和KEGG分析都表明BPA主要影响细胞分泌、细胞膜成分、胆固醇代谢和类固醇生成。通过图6可以看出涉及类固醇生成的详细成分；其涉及的主要信号通路是IGF‑1/Jak‑STAT/PI3k‑Akt和LHCGR/G‑stimulating蛋白/腺苷酸环化酶（AC）。

[0089] 进一步的，为了证实PI3K/AKT信号通路确实被激活，该方法通过测试该通路关键成分的磷酸化来进一步分析该通路，如图7所示。低剂量BPA可显著上调胰岛素受体底物‑1（IRS1）、AKT1和cAMP反应元件（CRE）结合蛋白（CREB）3个关键成分的磷酸化，而不影响其总蛋白水平。通过这些结果可以体现，青春期暴露于BPA可显著重置成人PI3K/AKT信号通路的基础活性，这可能在成年睾丸间质细胞依赖BPA增加类固醇合成中起主要作用，从而解释BPA的作用机理。

[0090] 以上实例，只是本发明优选地具体实例的一种，本领域技术人员在本发明技术方案范围内进行的通常变化和替换都包含在本发明的保护范围内。