The Effects and Toxicity of Different Pyrene Concentrations on Escherichia coli Using Transcriptomic Analysis
基于转录组分析探究不同芘浓度对大肠杆菌的作用与毒性
来源：Microorganisms, Volume 12, Article ID 326, 2024, DOI:10.3390/microorganisms12020326
《微生物》，第12卷，文章编号326，2024年，DOI：10.3390/microorganisms12020326
摘要
芘是典型多环芳烃环境污染物，具备致畸、致癌毒性，微生物耐受芘的分子机制尚不清晰。本研究设置300、600、1000 mg/L三组芘胁迫浓度，利用转录组测序解析大肠杆菌全局基因转录变化。结果显示随芘浓度升高，差异表达基因数量同步上升；芘毒性主要扰动糖代谢、膜转运、硫酸盐还原、氧化还原酶、多重外排泵相关基因表达。研究筛选大量差异sRNA并构建sRNA与靶基因调控网络，锁定芘应激核心调控基因。该研究完整揭示大肠杆菌应对芘污染的多层次转录应答通路，为多环芳烃微生物降解、环境毒理相关研究提供转录组数据支撑。
关键词
芘；大肠杆菌；多环芳烃；转录组；差异表达基因；sRNA；膜转运；氧化应激；硫酸盐还原；外排泵
研究目的
1 明确梯度芘暴露下大肠杆菌生长抑制规律，量化不同浓度芘的毒性强弱。
2 解析300、600、1000 mg/L芘胁迫下大肠杆菌全局mRNA转录差异，筛选关键应激代谢通路。
3 鉴定芘响应差异sRNA，构建sRNA-靶基因调控网络，挖掘核心调控分子。
4 阐明糖代谢、硫酸盐还原、膜转运、多重外排泵通路在芘耐受中的分子功能。
5 建立多环芳烃胁迫细菌转录组研究范式，为微生物修复多环芳烃污染提供理论基础。
研究思路
1 梯度培养实验：设置0-10 g/L芘浓度梯度，使用芬兰Bioscreen全自动微生物生长仪连续测定大肠杆菌OD600生长曲线，确定亚致死转录组取样浓度（300、600、1000 mg/L）。
2 转录组建库测序：对照组与三组芘胁迫组分别提取总RNA，去除rRNA后构建cDNA文库，Illumina平台高通量测序，质控后比对参考基因组计算基因FPKM表达量。
3 差异基因筛选：以|log2FC|≥1、p≤0.05为阈值筛选差异基因，GO、KEGG功能富集分析各通路表达变化。
4 qRT-PCR验证：选取5个代表性差异基因，荧光定量PCR验证转录组测序结果可靠性。
5 sRNA分析：鉴定全部细菌sRNA，筛选胁迫差异小RNA，富集其靶基因功能并构建互作调控网络图。
6 关键通路深度解析：重点分析碳水代谢、CysB硫酸盐还原系统、氧化还原酶、膜转运与RND多重外排泵通路差异基因表达特征。
7 整合生长表型、转录组、sRNA调控数据，完整构建大肠杆菌芘胁迫应答分子模型。
研究亮点
1 首次系统覆盖低、中、高三种亚致死芘浓度，动态揭示细菌转录应答随污染物浓度梯度变化的规律。
2 同步解析mRNA与调控sRNA双重转录组，搭建芘应激sRNA-靶基因调控网络，挖掘CysB等核心耐受通路。
3 证实碳水代谢为细菌抵御芘胁迫核心供能通路，硫酸盐还原、氧化还原酶参与芘氧化解毒过程。
4 明确ABC、PTS转运蛋白与RND多重外排泵协同作用，通过摄取营养、外排芘缓解细胞毒性。
5 采用Bioscreen高通量生长定量精准划分芘毒性阈值，为转录组取样浓度设置提供可靠表型依据。
可延伸的方向
1 敲除CysB、sRNA0625等核心调控基因，体内验证其芘耐受功能。
2 改造大肠杆菌外排与氧化还原酶基因，构建高效芘降解工程菌株用于土壤水体修复。
3 拓展多种多环芳烃（苯并芘等）对比转录组，挖掘通用微生物抗PAH保守通路。
4 开展动物肠道模型，探究芘暴露对肠道大肠杆菌稳态与宿主毒理的联动影响。
5 基于差异启动子构建芘生物传感工程菌，实现环境芘快速定量检测。
测量的数据及研究意义
1 Bioscreen测得不同芘浓度大肠杆菌生长曲线OD时序数据（图1），清晰划分毒性区间：1-800 mg/L芘几乎不抑制生长，900 mg/L至4 g/L生长逐步下降，5 g/L以上完全停滞，精准筛选300/600/1000 mg/L三个亚致死适宜转录组检测浓度，避免细胞大量死亡导致测序样本失效。
2 转录组差异基因火山图、GO/KEGG富集柱状图（图2、图3），随芘浓度升高上调、下调差异基因数量同步增多，1000 mg/L胁迫转录扰动最强，碳水代谢、膜转运、硫酸盐还原、氧化应激通路显著富集，锁定芘应激核心功能模块。
3 qRT-PCR基因定量数据，5个目标基因表达趋势与转录组FPKM结果高度吻合，证明测序数据真实可靠。
4 差异sRNA及靶基因富集、互作网络图（图4），筛选得到180个大肠杆菌sRNA，其中sRNA0625为全局核心调控分子，关联硫酸盐代谢关键基因，证实sRNA通过调控糖、转运通路参与细菌芘适应性应答。
5 CysB硫酸盐还原系统基因表达倍数柱状图（图5），该通路全部基因在芘胁迫下显著上调，600 mg/L组诱导幅度最高，半胱氨酸合成提升细胞还原能力，参与芘氧化解毒，阐明抗氧化耐受机制。
6 碳水代谢、膜转运、外排通路差异基因统计表，糖酵解、丁酸代谢大量基因上调，为应激细胞提供充足能量，ABC、PTS转运系统摄取碳源，RND多重外排泵将胞内芘排出，双重降低毒性。
结论
1 芘对大肠杆菌存在浓度依赖性毒性：低浓度无明显生长抑制，浓度超过900 mg/L增殖能力持续下降，5 g/L以上完全阻断生长。
2 芘胁迫诱导大肠杆菌大规模转录重塑，上调基因以碳水代谢、硫酸盐还原、氧化还原酶、膜转运及多重外排基因为主，高浓度胁迫下转录扰动程度显著增强。
3 sRNA0625是芘应激核心调控小RNA，通过靶向硫酸盐、碳水代谢关键基因调控细菌抗毒生理过程。
4 CysB介导的硫酸盐还原通路可提升胞内还原力，氧化还原酶参与芘氧化分解；ABC、RND转运系统协同完成营养摄取与芘外排，共同构成细菌芘耐受分子网络。
5 本研究提供梯度芘胁迫下细菌mRNA与sRNA完整转录组资源，为多环芳烃微生物毒性、生物修复机制研究提供数据与理论支撑。
使用芬兰 Bioscreen仪器测量数据的研究意义
1 本实验采用芬兰FP-1100型Bioscreen全自动微生物生长分析仪，100孔专用培养板同步批量检测多梯度芘浓度组别，37℃恒温持续振荡，每小时自动采集OD600，全程无需人工定时取样涂布，大幅降低人为操作带来的读数误差与实验工作量。
2 仪器可连续记录长达48小时完整生长周期，清晰区分延滞期、对数期、稳定期的生长差异，精准量化不同芘浓度对增殖抑制程度，仅终点单点OD无法区分梯度毒性细微差别。
3 标准化统一培养基、芘溶剂丙酮添加量、接种初始OD、振荡条件，所有实验组对照基准完全一致，数据统计学可比性强，精准划分亚致死浓度区间，为转录组实验筛选最优取样梯度。
4 生长曲线宏观表型与后续转录组分子数据形成完整证据链，先明确毒性表型差异，再从转录层面解析内在分子调控机理，让组学实验设计具备可靠生理依据。
5 高通量多孔检测模式一次性完成0–10 g/L大范围浓度梯度平行测定，快速筛选毒性阈值，省去多次摇瓶梯度预实验，大幅缩短前期条件摸索周期。