畜牧兽医学报  2023, Vol. 54 Issue (7): 2858-2871. DOI: 10.11843/j.issn.0366-6964.2023.07.018    PDF    
基于转录组学研究双酚A对猪睾丸支持细胞炎症和氨基酸代谢通路的影响
胡婷, 张永红, 侯晓林, 姚华, 崔德凤, 潘早早, 张凌宇, 张家希, 吴琼     
北京农学院动物科学技术学院,北京 102206
摘要:旨在研究猪睾丸支持细胞(ST)暴露于环境雌激素双酚A(BPA)不同时间后的基因表达谱变化。本研究将猪ST细胞设为C、B两组,每组设置3个重复,C组为空白对照组,B组细胞暴露于50 μmol·L-1浓度的BPA。BPA暴露6、24 h后,收集C6、B6、C24、B24组细胞总RNA样品,构建测序文库,通过双端150 bp测序方法进行高通量转录组测序。对组装数据进行功能注释、差异基因分析以及GO和KEGG富集分析。通过Real-time PCR方法对关键差异基因表达进行验证。BPA暴露24 h后差异基因数量(6 928个)较BPA暴露6 h差异基因数量(3 940个)明显增加。编码分泌型磷蛋白1的SPP1基因在BPA暴露6 h后表达增加,但在BPA暴露24 h后表达降低。管腔结合蛋白编码基因BIP、泛素B编码基因UBB、泛素C编码基因UBC、鸟氨酸脱羧酶1编码基因ODC1表达量在BPA暴露6和24 h后均显著上调。富集分析结果表明,BPA暴露6 h后,TNF信号通路富集差异倍数最高,p53、IL-17以及MAPK信号通路等也均显著富集,暗示ST细胞表现出明显的促炎反应。在BPA暴露24 h后,氨基酸及糖代谢通路富集差异倍数最高,包括精氨酸和脯氨酸、组氨酸以及糖酵解代谢通路。Real-time PCR验证结果显示,BPA暴露6 h后炎性反应通路相关的PTGS2等以及BPA暴露24 h后氨基酸代谢通路相关ARG1等差异基因表达显著上调。综上可知,BPA对猪睾丸支持细胞的损伤机制呈现明显的时间效应,ST细胞在BPA暴露初期表现出以TNF等信号通路激活为特征的炎性反应,在BPA暴露后期氨基酸及糖代谢显著激活。
关键词双酚A    猪睾丸支持细胞    转录组测序    富集分析    
The Effects of Bisphenol A on Inflammation and Amino Acid Metabolism Pathways in Porcine Testis Sertoli Cells Based on Transcriptome Analysis
HU Ting, ZHANG Yonghong, HOU Xiaolin, YAO Hua, CUI Defeng, PAN Zaozao, ZHANG Lingyu, ZHANG Jiaxi, WU Qiong     
Animal Science and Technology College, Beijing University of Agriculture, Beijing 102206, China
Abstract: This experiment was conducted to compare the gene expression profile changes of porcine testicular sertoli cells (ST) exposed to environmental estrogen bisphenol A (BPA) for different times. The porcine ST cells were divided into group C and group B with 3 replicates per group. Group C was a blank control group, and group B cells were exposed to BPA at a concentration of 50 μmol·L-1. At 6 and 24 h after BPA exposure, total RNA samples from C6, B6, C24 and B24 group cells were collected, and a sequencing library was constructed. High-throughput transcriptome sequencing was performed using a paired-end 150 bp sequencing method. The assembly data were used for functional annotation, differentially expressed gene (DEG) analysis, and GO and KEGG enrichment analysis. The expression of some key DEGs were validated through real-time PCR method. The number of DEGs (6 928) significantly increased at 24 h after BPA exposure, compared with 3 940 DEGs at 6 h. The expression of SPP1 gene encoding secretory phosphoprotein 1 increased at 6 h after BPA exposure, but decreased at 24 h after BPA exposure. The expression of lumen binding protein-encoding gene BIP, ubiquitin B-encoding gene UBB, ubiquitin C-encoding gene UBC, and ornithine decarboxylase 1-encoding gene ODC1 were significantly upregulated both at 6 and 24 h after BPA exposure. Enrichment analysis showed that the TNF signal pathway was the most significantly enriched at 6 h, and p53, IL-17 and MAPK signal pathways were also significantly enriched, indicating that ST cells exhibited proinflammatory response. At 24 h after BPA exposure, DEGs in pig ST cells were enriched in amino acids and sugar metabolism pathways, including arginine and proline, histidine and glycolysis metabolic pathways. Real-time PCR results showed that PTGS2 and other DEGs related to the inflammatory response pathway at 6 h after BPA exposure and ARG1 and other DEGs related to the amino acid metabolic pathways at 24 h after BPA exposure were significantly up-regulated. In conclusion, the damage mechanism of BPA on porcine sertoli cells shows obvious time effect. Porcine ST cells exhibit inflammatory responses characterized by activation of TNF and other signaling pathways at the early stage of BPA exposure, and the amino acid and glucose metabolism pathways are significantly activated in ST cells at the late stage of BPA exposure.
Key words: bisphenol A    porcine testis sertoli cells    transcriptome sequencing    enrichment analysis    

猪的繁殖性能关乎生猪出栏量以及生猪供需稳定。现代化猪场中广泛采用人工授精方式,公猪出现繁殖障碍造成的经济损失远远超过母猪[1]。公猪处于不良环境时,精子生理功能受损,造成公猪繁殖力降低。引起公猪繁殖障碍的因素包括日粮中的能量、微量元素、蛋白质、矿物元素的水平、饲养管理不当、环境因素、遗传因素、传染病等[2]。公猪繁殖与精子发生、成熟、获能、顶体反应、受精等过程密不可分,精子的发生与睾丸生殖细胞相关[3]。睾丸支持细胞(sertoli cell,ST细胞)是睾丸生殖细胞的一种,呈不规则的高锥体形,细胞基部附着在基膜上,顶部伸至曲精小管腔面[4]。ST细胞参与构成血睾屏障,创造了相对稳定的生精环境[5],ST细胞能够支持、营养、促进生精细胞的产生[6],分泌雄激素结合蛋白,维持较高的雄激素水平[7],分泌睾丸液,有助于精子运送[8]。此外,ST细胞能够吞噬消化精子形成过程脱落的残余胞质[9]

环境雌激素是一种环境内分泌的干扰化学物质,具有与生物激素相似的性质[10],可干扰体内正常内分泌物质的合成、释放、运输、结合和代谢,从而激活或抑制内分泌系统的功能[11]。在众多环境雌激素中,双酚A(BPA)结构与雌激素相似,具有拟雌激素和抗雄激素的作用因而具有重要的生殖发育毒性[12-15],广泛存在于塑料制品及涂料中,主要通过消化道、呼吸道和皮肤接触进入机体[16]。多项研究结果显示,大多数受检者的体液和组织中均可检测到BPA[17]。试验动物和人类研究证实,BPA对间质细胞发育和功能具有内分泌干扰作用[18]。BPA破坏间质细胞的发育和功能,引起相关的生殖疾病,如睾丸发育不全综合征,青春期延迟和不育[19]。成年大鼠皮下注射0.5% BPA,睾丸重量显著减轻,精子发生降低[20]。小鼠在孕期暴露于BPA,其雄性后代每日产生的精子数量减少,生殖器官明显变小[21]。BPA诱导活性氧产生,干扰ST细胞类固醇的生成[22],诱导ST细胞凋亡,生殖细胞成熟中断,导致睾丸发育缺陷[23]。然而,BPA影响猪睾丸支持细胞功能并造成公猪生殖障碍的机制目前还不清楚。

为探究环境雌激素双酚A影响公猪繁殖能力的作用机制,本研究分析了不同时间BPA暴露下猪睾丸支持细胞转录组学变化,进一步通过富集分析明确BPA暴露下猪睾丸支持细胞信号通路的变化。

1 材料与方法 1.1 试验材料与设备

猪睾丸支持细胞(ST细胞)由东北农业大学徐世文教授馈赠,冻存于本实验室。胎牛血清、青/链霉素双抗、DMEM培养基、胰酶、DMSO溶液、Tri-zol试剂购自美国赛默飞世尔科技公司,BPA粉末购自罗恩公司。ST细胞在含10%胎牛血清、100 μg·mL-1链霉素和100 IU·mL-1青霉素的DMEM完全培养基中培养,置于5% CO2的37 ℃恒温培养箱中。测序平台使用美国Ilumina公司的Novaseq 6000。

1.2 BPA工作液配制

BPA粉末溶解于DMSO溶液中配制为1 mol·L-1原液,然后倍比稀释为10-3 mol·L-1 BPA溶液,-20 ℃保存备用。细胞暴露前,取500 μL的10-3mol·L-1 BPA溶液溶于10 mL的DMEM完全培养基,制备成终浓度为50 μmol·L-1的BPA工作液,0.22 μm滤膜过滤备用。

1.3 BPA暴露处理

ST细胞培养于6孔板,浓度为2.5×106个·孔-1,试验分为两组:空白对照组(C)和BPA暴露组(B)。贴壁后更换培养基为无血清DMEM培养基+1%双抗培养基,饥饿处理12 h,BPA处理组细胞更换含50 μmol·L-1的BPA的DMEM完全培养基,空白对照组细胞更换为含与处理组细胞等量0.1% DMSO的DMEM完全培养基。每种处理条件设置3个重复孔,BPA暴露6、24 h后收取细胞样品进行后续RNA提取。C组以及B组细胞重复孔按照不同培养时间命名为C6-1、C6-2、C6-3,C24-1、C24-2、C24-3,B6-1、B6-2、B6-3,B24-1、B24-2、B24-3

1.4 RNA提取、文库构建、库检及上机测序

每个细胞处理孔加入1 mL Trizol进行RNA提取,利用NanoDrop 2000与Qubit对提取的细胞RNA浓度及纯度进行检测。质检合格的RNA制备文库,对核糖体RNA进行酶切,随机分成250~300 bp片段。利用Novaseq 6000测序平台对RNA文库进行双端150 bp高通量转录组测序,将测序图像信号经CASAVA碱基识别转换成文字信号,以Fastq格式储存为原始数据。

1.5 数据比对组装与功能注释

利用FastQC去除原始Fastq数据中低质量和污染的reads后,获得clean reads通过Tophat v2.1.1与猪参考基因组Sus scrofa11.1进行比对,利用StringTiev1.3.3b进行转录组组装。对获得的基因分别基于NR、Swiss-Prot、Pfam、COG、GO和KEGG数据库进行BLAST功能注释。

1.6 基因表达与差异基因分析

利用表达定量软件RSEM v1.3.1对基因表达水平进行定量分析,定量指标为TPM。使用DESeq2 v1.24.0鉴定空白对照组与BPA暴露组细胞之间的差异表达基因(DEGs),使用Benjamini和Hochberg方法(BH)调整P值控制错误发现率,两组之间表达倍数变化(FC)≥ 2或≤0.5,且adjusted P<0.05的基因名被认定为DEGs。

1.7 功能富集分析

基于使用Goatools v0.6.5的Fisher精确测试进行GO富集分析,用以确定DEGs的潜在作用。使用KOBAS v2.1.1进行KEGG富集分析,用以评估显著富集的代谢或信号转导途径。经BH调整的P<0.05的GO功能或KEGG通路被认为显著富集。

1.8 荧光定量PCR

应用Trizol法对处理后的细胞进行总RNA提取后,应用核酸浓度检测仪进行RNA质量检测,使用Promega RT reagent Kit(Perfect Real Time)进行RNA反转录,用7500 FAST PCR仪进行qPCR。mRNA引物序列见表 1

表 1 mRNA引物序列 Table 1 mRNA primers sequence
2 结果 2.1 测序质控和样本间相关性分析

本次转录组测序12个样品共获得87.63 Gb高质量读取条带用于后续组装和分析,每个样品数据量平均为5.94 Gb。所有样本的Q20和Q30平均值分别为97.4%和93%,GC含量平均为53.7%,表明测序获得的读取条带质量高。基于Spearman相关性分析对所有12个样本的TPM分布构建相关矩阵,结果表明,不同BPA暴露时间下的细胞与空白对照组细胞样品明显分离,组内重复样品能够聚集在一起,表明各组测序结果一致性高,测序数据的再现性高(图 1)。

不同颜色的方块代表两个样本的相关性 The squares with different colors represent the correlation between the two samples 图 1 样本间相关性热图 Fig. 1 Correlation heatmap between samples
2.2 参考基因功能注释分析

从所有样本中共发现33 004个基因,包括31 908个参考基因(96.7%)和1 096个未注释的新基因(3.3%)。在31 908个参考基因中,共22 279(69.8%)个基因能够通过COG(18 137)、GO(18 528)、KEGG(13 218)、NR(21 166)、Swiss Prot(19 103)和Pfam(15 947)数据库进行注释(图 2A)。

A.基于不同数据库功能注释的Upset图;B.COG功能注释 A.Upset diagram based on different database function annotations; B. COG function annotation 图 2 参考基因功能注释图 Fig. 2 Function annotation diagram of referenced genes

随后将参考基因比对到COG数据库中,结果显示共获得18 137个COG功能注释信息,划分为23个功能分类(图 2B)。与细胞内运输、分泌和囊泡运输相关的基因分布最多(U),有3 431个基因。涉及翻译后修饰、蛋白质周转、分子伴侣蛋的功能基因(O)数量次之,有1 887个基因。基因功能尚未可知的基因(S类)数量最多,达到12 198个,这说明BPA暴露后ST细胞应答反应仍需深入挖掘。

GO注释结果显示,共有18 528个基因能够被注释为58个GO分类(图 3)。在生物学过程类别中,细胞过程、生物调节相关基因最多,分别为13 998和11 790个。在细胞成分类别中,注释到细胞膜结构和细胞组成成分相关基因最多,分别为15 988和9 927个。在分子功能类别中,与结合和催化活性相关基因最多,分别为12 354、8 852个。

图 3 GO注释分析 Fig. 3 GO annotation analysis

KEGG注释结果表明(图 4),注释到环境信息处理的信号转导的参考基因数量最多(2 036个),注释到免疫系统相关的参考基因有1 130个,与内分泌系统相关基因有1 024个。

图 4 KEGG注释分析 Fig. 4 KEGG annotation analysis
2.3 新转录本预测与样本间Venn分析

对新转录本进行功能注释(图 5A),与内含子链完全匹配的转录本数量最多,达63 041个。至少有一个剪接与参考转录本共享的潜在新转录本片段(j)数量为16 753个,未知的转录本(u)数量为1 328个,与参考转录本外显子所处链的反义链有交集的转录本(x)数量为727个,完全落入参考转录本内含子区的转录本(i)有611个,与参考转录本的外显子有一定的交集的转录本(o)数量为458个。

A.新转录本类型分布图;B. 组间参考基因Venn图 A. The novel transcript types distribution; B. Venn analysis of reference genes between groups 图 5 新转录本预测和Venn分析 Fig. 5 Prediction of novel transcripts and Venn analysis

对参考基因进行Venn分析获取不同处理组间共有或特有的基因(图 5B)。BPA暴露6、12 h的细胞与空白对照组细胞共有基因10 601个,BPA暴露6 h后的B6组细胞特有基因79个,而BPA暴露24 h后的B24组细胞特有基因增加到291个。

2.4 差异基因分析

与空白对照组细胞相比,BPA暴露6 h导致3 940个基因表达显著差异(图 6A),其中2 261个基因表达上调,1 679个基因表达下调。BPA暴露6 h后,基因表达上调差异倍数前10的基因为ACTG1、SPP1、TSP1、BIPTFPAI-1、UPAODC1、UBCUBB,基因表达下调差异倍数前10基因为TUBA1B、HSPG2、SOD3、S100A6、EMILIN2、TUBB2B、P311、LGALS1、AHCYFAT1。

BPA暴露6 (A)和24 h(B)的差异基因散点图 Scatter plots of DEGs from cells exposed to BPA for 6 (A) and 24 h (B) 图 6 差异基因散点图 Fig. 6 Scatter plots of DEGs

BPA暴露24 h后,差异基因数量达到6 928个,包括3 222个上调基因和3 706个下调基因(图 6B)。上调差异倍数前10基因分别为FTLBIPHSPA8、UBBUBCHsp90、ODC1、CLUANXA1和RBP4,下调差异倍数前10的基因分别为MT-CO1、FN1、GAPDHACTG1、SPP1、ENO1、APLP2、COIIILDHATMSB10。从结果可以看出,SPP1基因在BPA暴露6 h后表达增加,但在BPA暴露24 h后表达降低。在BPA暴露6和24 h后,BIPUBBUBCODC1四个基因表达均明显升高。

2.5 差异基因GO富集分析

对BPA暴露6、24 h的差异基因进行GO富集分析。如图 7A所示,展示了富集差异倍数前20的信号通路。ST细胞BPA暴露6 h后,与泛素依赖性蛋白质分解代谢过程调节相关的差异基因最多,达到46个,其次是T细胞分化调节(42个)。BPA暴露6 h的差异基因GO富集弦图结果显示(图 7B),参与泛素依赖性蛋白质分解代谢过程调节过程上调差异倍数前10的差异基因有IL-23、p35、p21、SOX9、TRIB3、SOCS1、PI3K、NLRP3、TBX20和NF-κB, 下调差异倍数前10的基因有SLA-DQB1、CDKN2C、SMAD6、SASH3、MYBDDA3、CDC20、PLK1、CCNB2和NR2F2。

BPA处理6 h后GO富集分析气泡图(A)和富集弦图(B);BPA处理24 h后GO富集分析气泡图(C)和富集弦图(D) GO enrichment analysis bubble map (A) and chordograms (B) after BPA exposure for 6 h; GO enrichment analysis bubble map (C) and chordograms (D) after BPA exposure for 24 h 图 7 GO富集分析 Fig. 7 GO enrichment analysis

ST细胞BPA暴露24 h后差异基因GO富集分析显示(图 7C),富集通路与BPA暴露后6 h明显不同,糖酵解过程富集差异倍数最高(adjusted P-value=3.78×10-4),富集基因23个。细胞外基质组织、有机酸分解代谢过程、羧酸分解代谢过程富集差异基因数量最多,分别为101、94、94个。GO富集弦图显示(图 7D),在参与糖酵解过程的基因中,除PFKFB2基因上调外,其他富集的22个基因均下调,下调差异倍数前10的差异基因有GPIPGK1、PGAM1、HKDC1、TPI1、PFKLHK1、ENO1、GAPDHALDOC

2.6 差异基因KEGG富集分析

图 8A所示,在KEGG富集分析中展示了富集差异倍数前20的信号通路。BPA暴露6 h后,TNF信号通路(adjusted P-value=5.90×10-15)、癌症相关microRNA(adjusted P-value=1.40×10-11)、p53信号通路(adjusted P-value=6.94×10-9)、IL-17信号通路(adjusted P-value=7.52×10-8)、MAPK信号通路(adjusted P-value=2.88×10-7)等炎症信号通路富集差异倍数较高。KEGG富集弦图显示(图 8B),参与TNF信号通路上调差异基因包括PTGS2、LIFCREB5、NF-κBCCL2、IL6、BIRC2、CREB3L3、CXCL10。BPA暴露24 h后,富集差异倍数前3的通路为精氨酸和脯氨酸代谢、组氨酸代谢、糖酵解途径等代谢通路(图 8C),而MAPK信号通路等炎症信号通路在BPA暴露后24 h富集并不明显。KEGG富集弦图显示(图 8D),参与精氨酸和脯氨酸代谢通路上调差异基因包括ARG1、SMOXSAT1、ODC1,下调差异基因包括AGMATDAOAOC1、GATMGAMTL3HYPDH

BPA处理6 h后KEGG富集分析气泡图(A)和富集弦图(B);BPA处理24h后KEGG富集分析气泡图(C)和富集弦图(D) KEGG enrichment analysis bubble map (A) and chordograms (B) after BPA exposure for 6 h; KEGG enrichment analysis bubble map (C) and chordograms (D) after BPA exposure for 24 h 图 8 KEGG富集分析 Fig. 8 KEGG enrichment analysis
2.7 关键差异基因表达验证

用Real-time PCR对转录组测序结果进行验证。BPA暴露后6 h,TNF等炎症信号通路相关差异基因PTGS2、LIFNF-κBIL-6、CCL2、TRAF3表达显著升高(图 9)。BPA暴露后24 h,精氨酸和脯氨酸等代谢通路相关差异基因ARG1、SMOXSAT1、ODC1、FTLPFKFB2表达显著上调(图 9B)。

Real-time PCR分析BPA处理6(A)、24 h(B)关键差异基因表达 Real-time PCR analysis of DEGs after BPA exposure for 6 (A) and 24 h (B) 图 9 差异基因表达Real-time PCR分析 Fig. 9 Real-time PCR analysis of DEGs
3 讨论

BPA属低毒性化学物,动物试验发现BPA有模拟雌激素的效果[24],低剂量也能导致动物雌性早熟、精子数下降、前列腺增长等。此外,双酚A具有一定的胚胎毒性和致畸性,可明显增加动物卵巢癌、前列腺癌、白血病等癌症的发生。生殖系统是BPA作用的主要靶器官,BPA对生殖系统的损伤作用体现在抑制间质细胞中睾丸类固醇生成酶和类固醇生成,增加炎症细胞数量,导致睾丸发育缺陷,破坏生殖细胞的成熟,改变雌激素相关基因表达,精子生成受限,降低精子质量及运动能力。本研究使用转录组测序分析了猪睾丸支持细胞暴露于BPA不同时间后的基因表达谱变化。与小鼠试验模型相比,猪与人体在生理结构上更加接近,本研究结果也为探究BPA对人生殖系统的作用提供了参考。

3.1 BPA对炎性反应通路的作用

本研究发现,猪睾丸支持细胞BPA暴露6 h后差异基因功能富集在TNF、MAPK、IL-17等炎性信号通路。TNF是细胞凋亡以及炎症和免疫的主要介质,TNFα受体分为2种(TNFRⅠ和TNFRⅡ),存在于多种细胞表面,TNF-α与TNFR的结合激活了多种信号通路,通常会引起细胞凋亡、炎症和肿瘤的发生等。BPA暴露诱导PK-15细胞TNF-α表达升高,利用qRT-PCR技术检测其他炎性因子IL-6、IL-8、IL-1β等均表达增高,证实BPA引起猪肾细胞的炎性反应[25]

MAPK信号通路是TNF/TNFR1信号下游通路。MAPK是信号从细胞表面传导到细胞核内部的重要传递者,MAPK通路有MAPK激酶激酶、MAPK激酶和MAPK三种重要的关键酶,3种酶共同调节细胞的生长分化、炎症反应、对环境的应激适应等重要的细胞生理过程。研究证实,双酚类化合物能够激活人卵巢皮质细胞的ERβ/GPR30-MAPK信号通路,导致细胞炎症的发生[26]。本研究发现,BPA诱导TNF通路中TFNR1以及TRAF2/5、TAK1显著上调,进而激活MKK4/7、MKK3/6和JNK高表达。

IL-17是一种促炎因子,能够诱导上皮细胞、内皮细胞、成纤维细胞分泌合成IL-6、IL-8等。JNK/NF-κB信号通路是典型的细胞炎症反应通路。巨噬细胞敲除IL-17A基因后,BPA激活JNK/NF-κB以及诱导脂肪组织炎症的作用明显减弱[27]。研究表明,血清双酚A水平与IL-23和IL-17呈正相关,诱导全身炎症反应[28]。本研究发现,BPA诱导IL-17基因上调,IL-17通路下游的NF-κB和炎症因子IL-23表达均明显增加。

研究发现,IL-17与TNF-α共同作用可上调皮肤角质形成细胞系HaCaT细胞表达IL-17R和p-p38 MAPK,进一步促进IL-6、IL-8、MIP-3α等炎症因子的分泌[29]。四氯化碳通过TNF/MAPK/NF-κB信号通路协同作用诱导小鼠细胞炎症,促进TNF-α、IL-6等促炎症因子上调,加剧了NF-κB、ERK、MAPK蛋白的磷酸化,导致小鼠肝纤维化的发生[30]。BPA和BPS通过MAPK/JUN通路诱导绵羊卵巢颗粒细胞的炎症反应,降低了类固醇的分泌,影响生殖系统的功能[31]。硫氢化钠通过抑制ERK/JNK/MAPK通路缓解BPA导致的鼠肺泡上皮细胞炎症反应的发生[32]。这些研究表明,不同促炎细胞通路之间互相联系影响。本研究结果表明,BPA暴露通过促进TNF、p53和MAPK通路协同激活诱导细胞炎症。

3.2 BPA对细胞增殖的作用及致癌机制

作为环境雌激素的一种类型,已经有大量的研究证实BPA具有促进细胞增殖的作用及致癌效应。本研究发现,猪睾丸支持细胞暴露于BPA后,p53信号通路显著富集。p53是一种抑癌蛋白,p53基因是细胞周期中的负调节因子,在周期调控、细胞分化、细胞凋亡、DNA修复等重要的生物学过程中发挥着重要作用,p53蛋白的过表达与肿瘤的转移复发等相关[33]。研究发现,BPA通过ERK、EGFR、AR和ERβ等信号途径,导致p53磷酸化,进而激活EGFR/ERK/p53信号通路,诱导原发性和前列腺癌细胞的细胞周期阻滞[34],证实生殖细胞暴露于BPA后,p53表达异常激活了p53相关的信号通路导致生殖细胞增殖,诱发生殖系统癌症。p53在细胞凋亡过程发挥了重要的作用,BPA暴露于小鼠精原细胞GC-2后,p53和Bcl-2基因上调,细胞凋亡率升高[35]。羊生殖细胞暴露于BPA后,p53通路激活,凋亡蛋白Caspase-3、p53和Bcl-2表达量显著上调,细胞发生凋亡,生殖细胞功能降低[36]。BPA暴露也上调细胞凋亡和p53通路相关基因表达[37]

3.3 BPA对糖酵解途径的作用

KEGG富集分析显示,BPA暴露24 h后,猪ST细胞糖酵解和氨基酸代谢显著富集。研究表明,高浓度的BPA暴露会增加基础葡萄糖摄取,促进糖酵解,但较低水平的双酚A没有此类效果[38]。亚毒性浓度的BPA暴露不仅促进癌细胞系OVCAR-3细胞的增殖,同时促进糖酵解代谢,增加细胞内ATP、乳酸和丙酮酸产生[39]。短期或长期暴露双酚类化合物也可以增加中性粒细胞糖酵解功能[40]

3.4 BPA对氨基酸代谢的作用

氨基酸代谢是BPA暴露24 h后富集差异倍数最高的信号通路。大量研究证明,BPA暴露是导致氨基酸代谢紊乱的原因之一。代谢组学研究表明,50 μg·kg-1的BPA暴露可诱导SD大鼠血浆中19种代谢物表达上调,32种代谢物表达下调,其中柠檬酸循环、丁酸代谢、丙氨酸、天冬氨酸和谷氨酸代谢以及甘油磷脂代谢显著紊乱[41]。妊娠期母鼠BPA暴露可以导致子代公鼠谷氨酸转运体GLAST以及GLT-1蛋白的上调表达,谷氨酸代谢酶GS、GLS、GDH下调表达[42]

4 结论

本研究利用转录组分析技术探究了猪睾丸支持细胞暴露于BPA不同时间后的基因表达谱。BPA暴露6 h后,猪睾丸支持细胞表现出明显的促炎反应,TNF信号通路富集差异倍数最高。与BPA暴露6 h不同的是,猪睾丸支持细胞在BPA暴露24 h后氨基酸及糖代谢通路富集,包括精氨酸和脯氨酸、组氨酸、糖酵解途径等代谢通路。本研究初筛了BPA损伤猪睾丸支持细胞功能的信号通路及关键基因,为进一步探究BPA诱导公猪生殖障碍的作用机制奠定了基础。

参考文献
[1]
田秀娥, 王永军, 胡建宏, 等. 环境温度与公猪繁殖力[J]. 家畜生态学报, 2021, 42(3): 86-89, 93.
TIAN X E, WANG Y J, HU J H, et al. Ambient temperature and boar fertility[J]. Journal of Domestic Animal Ecology, 2021, 42(3): 86-89, 93. (in Chinese)
[2]
闫益波, 焦福林, 李文刚. 热应激对种公猪繁殖性能影响的研究进展[J]. 养猪, 2020(4): 41-44.
YAN Y B, JIAO F L, LI W G. Effect of heat stress on reproductive performance of breeding boars[J]. Swine Production, 2020(4): 41-44. DOI:10.13257/j.cnki.21-1104/s.2020.04.018 (in Chinese)
[3]
YANG H, MA J Y, WAN Z, et al. Characterization of sheep spermatogenesis through single-cell RNA sequencing[J]. FASEB J, 2021, 35(2): e21187.
[4]
LIU Y, SUN B X, ZHANG S X, et al. Glycine alleviates fluoride-induced oxidative stress, apoptosis and senescence in a porcine testicular Sertoli cell line[J]. Reprod Domest Anim, 2021, 56(6): 884-896. DOI:10.1111/rda.13930
[5]
COUTURE R, MARTIN L J. The transcription factors SF-1 and SOX8 cooperate to upregulate Cx43 expression in mouse TM4 sertoli cells[J]. Biochem Biophys Rep, 2020, 24: 100828.
[6]
REHMAN Z U, WORKU T, DAVIS J S, et al. Role and mechanism of AMH in the regulation of Sertoli cells in mice[J]. J Steroid Biochem Mol Biol, 2017, 174: 133-140. DOI:10.1016/j.jsbmb.2017.08.011
[7]
张倩, 周可柔, 胡倩倩, 等. 大鼠睾丸支持细胞的分离培养及鉴定[J]. 安徽科技学院学报, 2018, 32(2): 1-6.
ZHANG Q, ZHOU K R, HU Q Q, et al. Isolation and identification of Sertoli cells in rats[J]. Journal of Anhui Science and Technology University, 2018, 32(2): 1-6. (in Chinese)
[8]
伊敏娜, 才文道力玛, 陶力, 等. 睾丸支持细胞的功能、分离纯化及鉴定研究进展[J]. 中国畜牧兽医, 2021, 48(8): 2947-2956.
YIN M N, CAI W D L M, TAO L, et al. Research progress in the function, isolation, purification and identification of testicular Sertoli cells[J]. China Animal Husbandry & Veterinary Medicine, 2021, 48(8): 2947-2956. (in Chinese)
[9]
申毅锋, 朱坤, 李文元, 等. 基于睾丸支持细胞结构功能探讨中药改善生精功能的动物实验系统评价[J]. 中国实验方剂学杂志, 2022, 28(21): 201-211.
SHEN Y F, ZHU K, LI W Y, et al. Systematic review of Chinese medicine in improving sperm quality in animal experiments based on structure and function of Sertoli cells[J]. Chinese Journal of Experimental Traditional Medical Formulae, 2022, 28(21): 201-211. (in Chinese)
[10]
ULLAH H, JAHAN S, UI AIN Q, et al. Effect of bisphenol S exposure on male reproductive system of rats: A histological and biochemical study[J]. Chemosphere, 2016, 152: 383-391. DOI:10.1016/j.chemosphere.2016.02.125
[11]
李红梅, 熊忆茗, 徐海明, 等. 双酚A类似物的雄性生殖毒性研究进展[J]. 中国环境科学, 2021, 41(6): 2939-2945.
LI H M, XIONG Y M, XU H M, et al. Research progress in male reproductive toxicity of bisphenol A analogues[J]. China Environmental Science, 2021, 41(6): 2939-2945. (in Chinese)
[12]
梁大廷, 戴旭超, 孙发. 双酚A暴露对男性精子生成的影响研究进展[J]. 中华男科学杂志, 2020, 26(11): 1030-1034.
LIANG D T, DAI X C, SUN F. Impacts of bisphenol A exposure on spermatogenesis in men: advances in studies[J]. National Journal of Andrology, 2020, 26(11): 1030-1034. (in Chinese)
[13]
QIU W H, ZHAO Y L, YANG M, et al. Actions of bisphenol A and bisphenol S on the reproductive neuroendocrine system during early development in zebrafish[J]. Endocrinology, 2016, 157(2): 636-647. DOI:10.1210/en.2015-1785
[14]
ULLAH H, AMBREEN A, AHSAN N, et al. Bisphenol S induces oxidative stress and DNA damage in rat spermatozoa in vitro and disrupts daily sperm production in vivo[J]. Toxicol Environ Chem, 2017, 99(5-6): 953-965. DOI:10.1080/02772248.2016.1269333
[15]
ULLAH A, PIRZADA M, JAHAN S, et al. Bisphenol A and its analogs bisphenol B, bisphenol F, and bisphenol S: comparative in vitro and in vivo studies on the sperms and testicular tissues of rats[J]. Chemosphere, 2018, 209: 508-516. DOI:10.1016/j.chemosphere.2018.06.089
[16]
YAMADA A, ARAKAKI R, SAITO M, et al. Dual role of Fas/FasL-mediated signal in peripheral immune tolerance[J]. Front Immunol, 2017, 8: 403.
[17]
WANG M, SU P. The role of the Fas/FasL signaling pathway in environmental toxicant-induced testicular cell apoptosis: An update[J]. Syst Biol Reprod Med, 2018, 64(2): 93-102. DOI:10.1080/19396368.2017.1422046
[18]
KATO T, MIZUNO K, NISHIO H, et al. Disorganization of claudin-11 and dysfunction of the blood-testis barrier during puberty in a cryptorchid rat model[J]. Andrology, 2020, 8(5): 1398-1408. DOI:10.1111/andr.12788
[19]
DANCE A, KASTELIC J, THUNDATHIL J. A combination of insulin-like growth factor I (IGF-I) and FSH promotes proliferation of prepubertal bovine Sertoli cells isolated and cultured in vitro[J]. Reprod Fertil Dev, 2017, 29(8): 1635-1641.
[20]
BAHELKA I, STUPKA R, AČG ÍTEK J, et al. The impact of bisphenols on reproductive system and on offspring in pigs-a review 2011-2020[J]. Chemosphere, 2021, 263: 128203.
[21]
KHALAF A A, AHMED W, MOSELHY W A, et al. Protective effects of selenium and Nano-selenium on bisphenol-induced reproductive toxicity in male rats[J]. Hum Exp Toxicol, 2019, 38(4): 398-408.
[22]
ZHANG Y L, HAN L, YANG H, et al. Bisphenol A affects cell viability involved in autophagy and apoptosis in goat testis sertoli cell[J]. Environ Toxicol Pharmacol, 2017, 55: 137-147.
[23]
LIAO C Y, KANNAN K. A survey of bisphenol A and other bisphenol analogues in foodstuffs from nine cities in China[J]. Food Addit Contam Part A, 2014, 31(2): 319-329.
[24]
PAWLICKI P, DULIBAN M, TUZ R, et al. Do G-protein coupled estrogen receptor and bisphenol A analogs influence on Leydig cell epigenetic regulation in immature boar testis ex vivo?[J]. Anim Reprod Sci, 2019, 207: 21-35.
[25]
周羽, 刘媛, 汪序忠, 等. 双酚A对猪肾细胞炎性因子表达及凋亡的影响[J]. 中国畜牧杂志, 2021, 57(12): 223-227, 233.
ZHOU Y, LIU Y, WANG X Z, et al. Effect of bisphenol A on the expression of inflammatory factors and apoptosis in porcine kidney cells[J]. Chinese Journal of Animal Science, 2021, 57(12): 223-227, 233. (in Chinese)
[26]
杨芷蕊, 李晓蒙, 熊围, 等. 双酚S经ERβ/GPR30-MAPK信号通路介导对人正常卵巢上皮细胞的毒性作用[C]//2021(第五届)毒性测试替代方法与转化毒理学(国际)学术研讨会会议论文集. 沈阳: 中国毒理学会, 2021: 136-137.
YANG Z R, LI X M, XIONG W, et al. The toxic effect of bisphenol S on human normal ovarian epithelial cells is mediated by ERβ/GPR30-MAPK signaling pathway[C]//Proceedings of the 2021(5th) International Symposium on Alternative Methods of Toxicity Testing and Translational Toxicology. Shenyang: Chinese Society of Toxicology, 2021: 136-137. (in Chinese)
[27]
周义. IL-17A调控脂肪组织巨噬细胞极化介导双酚A诱导的脂肪组织炎症及胰岛素抵抗[D]. 合肥: 安徽医科大学, 2022.
ZHOU Y. IL-17A regulates adipose tissue macrophage polarization to mediate BPA-induced adipose tissue inflammation and insulin resistance[D]. Hefei: Anhui Medical University, 2022. (in Chinese)
[28]
LINARES R, FERNÁNDEZ M F, GUTIÉRREZ A, et al. Endocrine disruption in Crohn 's disease: Bisphenol A enhances systemic inflammatory response in patients with gut barrier translocation of dysbiotic microbiota products[J]. FASEB J, 2021, 35(7): e21697.
[29]
沈奎亚, 刘梅, 吴新安. IL-17与TNF-α对HaCaT细胞IL-17R、p-p38 MAPK和炎症因子表达的影响[J]. 安徽医药, 2016, 20(7): 1292-1295.
SHEN K Y, LIU M, WU X A. Co-effect of IL-17 and TNF-α on the expressions of IL-17R, p-p38 MAPK and inflammatory factors in cultured keratinocytes[J]. Anhui Medical and Pharmaceutical Journal, 2016, 20(7): 1292-1295. (in Chinese)
[30]
CAI Q, WANG Z Q, ZHANG R, et al. Huangjia Ruangan Granule inhibits inflammation in a rat model with Liver fibrosis by regulating TNF/MAPK and NF-κB signaling pathways[J]. Evid Based Complement Alternat Med, 2022, 2022: 8105306.
[31]
TÉTEAU O, CARVALHO A V, PAPILLIER P, et al. Bisphenol A and bisphenol S both disrupt ovine granulosa cell steroidogenesis but through different molecular pathways[J]. J Ovarian Res, 2023, 16(1): 30.
[32]
ABO-ZAID O A R, MOAWED F S M, HASSAN H A, et al. Hydrogen sulfide attenuates lung injury instigated by Bisphenol-A via suppressing inflammation and oxidative stress[J]. BMC Pharmacol Toxicol, 2022, 23(1): 98.
[33]
ZHAO C, YONG T, ZHANG Y B, et al. Breast cancer proliferation and deterioration-associated metabolic heterogeneity changes induced by exposure of bisphenol S, a widespread replacement of bisphenol A[J]. J Hazard Mater, 2021, 414: 125391.
[34]
BILANCIO A, BONTEMPO P, DI DONATO M, et al. Bisphenol A induces cell cycle arrest in primary and prostate cancer cells through EGFR/ERK/p53 signaling pathway activation[J]. Oncotarget, 2017, 8(70): 115620-115631.
[35]
KASONGO A A, LEROUX M, AMROUCHE-MEKKIOUI I, et al. BPA exposure in L6 myotubes increased basal glucose metabolism in an estrogen receptor-dependent manner but induced insulin resistance[J]. Food Chem Toxicol, 2022, 170: 113505.
[36]
SHI X Y, WANG Z, LIU L Y, et al. Low concentrations of bisphenol A promote human ovarian cancer cell proliferation and glycolysis-based metabolism through the estrogen receptor-α pathway[J]. Chemosphere, 2017, 185: 361-367.
[37]
谷莉莉. 雌激素受体及P53/P21信号通路在孕哺期双酚A暴露致子代大鼠生长板衰老中的作用[D]. 沈阳: 中国医科大学, 2022.
GU L L. Effect of estrogen receptor and P53/P21 signaling pathway in growth plate senescence in offspring induced by maternal exposure to bisphenol A during pregnancy and lactation[D]. Shenyang: China Medical University, 2022. (in Chinese)
[38]
ZHANG H, ZHENG Y, XIAO Y X, et al. Autophagy attenuates placental apoptosis, oxidative stress and fetal growth restriction in pregnant ewes[J]. Environ Int, 2023, 173: 107806.
[39]
袁文博, 刘文斌, 陈洪强, 等. 双酚A致雄性生殖损害毒性通路的富集分析及关键凋亡信号通路的实验研究[J]. 第三军医大学学报, 2021, 43(1): 38-44.
YUAN W B, LIU W B, CHEN H Q, et al. Enrichment analysis and experimental study of BPA-induced male reproductive damage toxicity pathways and key apoptosis signaling pathways[J]. Journal of Third Military Medical University, 2021, 43(1): 38-44. (in Chinese)
[40]
PEILLEX C, KEREVER A, LACHHAB A, et al. Bisphenol A, bisphenol S and their glucuronidated metabolites modulate glycolysis and functional responses of human neutrophils[J]. Environ Res, 2021, 196: 110336.
[41]
MAO L L, FANG S H, ZHAO M R, et al. Effects of Bisphenol A and Bisphenol S exposure at low doses on the metabolome of adolescent male sprague-dawley rats[J]. Chem Res Toxicol, 2021, 34(6): 1578-1587.
[42]
刘娣. 孕期双酚A暴露对子代大鼠海马组织中谷氨酸代谢影响的机制研究[D]. 沈阳: 沈阳医学院, 2021.
LIU D. Effects of prenatal bisphenol A exposure on glutamate metabolism in the hippocampus of rats offspring and the related mechanism[D]. Shenyang: Shenyang Medical College, 2021. (in Chinese)

(编辑   郭云雁)