Integrated Genomics and Transcriptomics Provide Insights into Salt Stress Response in Bacillus subtilis
一株源自雷竹笋加工废弃物枯草芽孢杆菌AC P81盐胁迫响应机制的整合基因组与转录组解析
来源：Microorganisms, Volume 12, Article ID 285, 2024, DOI:10.3390/microorganisms12020285
《微生物》第12卷，文章编号285，2024年，DOI：10.3390/microorganisms12020285
摘要
土壤盐渍化是制约农业生产关键难题，耐盐微生物是盐碱地生物改良优质资源。本研究从雷竹笋加工废料分离得到中度耐盐枯草芽孢杆菌ACP81，该菌株最高耐受20% NaCl，6% NaCl为亚致死胁迫浓度。6%盐胁迫下菌株纤维素酶、β-淀粉酶活性显著上升，中性蛋白酶活性下降。基因组测序显示菌株全长4215636 bp，包含大量碳水水解酶编码基因；对照（0%NaCl）与6%NaCl组转录组对比共筛选932个差异基因，盐胁迫上调鞭毛组装、相容性溶质合成、磷酸转移酶PTS、糖代谢通路基因，celB、celC、malL等纤维素与淀粉酶基因显著上调，淀粉蔗糖代谢通路富集；RT-qPCR验证转录组测序结果。盐胁迫诱导菌株合成葡萄糖等相容性溶质抵御渗透损伤，同时重塑细胞膜结构、提升运动能力，为盐碱地改良与工业水解酶生产提供优良菌种资源。
关键词
枯草芽孢杆菌；盐胁迫；基因组；转录组；相容性溶质；纤维素酶；β-淀粉酶；鞭毛组装；磷酸转移酶系统；耐盐微生物
研究目的
1 从竹笋加工废料分离、鉴定耐盐枯草芽孢杆菌ACP81，测定其耐盐生长阈值与亚致死盐浓度。
2 解析6% NaCl亚致死胁迫下菌株形态、水解酶、抗氧化酶活性变化规律。
3 完成菌株全基因组组装注释，结合两组转录组数据筛选盐胁迫响应差异基因，富集关键代谢通路。
4 阐明纤维素、淀粉酶基因上调介导的葡萄糖积累耐盐分子机制。
5 综合基因组、转录组、生理表型数据，明确ACP81完整耐盐调控网络，评估菌株应用潜力。
研究思路
1 菌株分离鉴定：竹笋废料梯度稀释涂布功能筛选平板，依据水解透明圈筛选目标菌株，16S rRNA测序分类鉴定。
2 耐盐生长测定：设置0、2、6、11、16、20% NaCl梯度，芬兰Bioscreen仪器连续记录生长曲线，确定亚致死浓度。
3 生理生化检测：0%与6%NaCl培养菌体，测定纤维素酶、淀粉酶、中性蛋白酶、SOD/CAT/POD抗氧化酶活性，TEM观测细胞形态。
4 基因组测序：Illumina平台完成全基因组测序，注释蛋白编码、CAZyme碳水酶基因。
5 转录组测序：无盐、6%盐胁迫对数期菌体建库测序，DESeq2筛选差异基因，GO、KEGG富集分析。
6 关键基因验证：挑选celB、malL、betB等耐盐核心基因开展RT-qPCR，验证转录组结果。
7 整合所有组学与表型数据，构建枯草芽孢杆菌盐胁迫适应调控模型。
研究亮点
1 首次从竹笋加工废弃物分离得到高耐盐枯草芽孢杆菌ACP81，耐受20% NaCl，兼具纤维素、淀粉酶分泌能力，兼具土壤改良与工业酶双重应用价值。
2 亚致死盐胁迫下菌株纤维素、β-淀粉酶显著激活，葡萄糖作为相容性溶质抵御渗透胁迫，提出新型耐糖代谢调控机制。
3 多组学联合完整解析耐盐通路：鞭毛组装、细胞壁D-丙氨酸修饰、相容性甜菜碱合成、PTS糖转运共同协同应对高盐环境。
4 抗氧化系统差异化调控：CAT、POD上调清除过氧化氢，SOD下调减少超氧阴离子，形成独特氧化应激平衡模式。
5 完整提供该菌株基因组与盐胁迫转录组数据集，为耐芽孢杆菌基因工程改造提供大量候选靶点。
可延伸的方向
1 敲除/过表达malL、celB基因，验证淀粉酶、纤维素酶对菌株耐盐能力的直接调控作用。
2 改造betB、gbsB相容性溶质合成基因，进一步提升菌株高盐耐受上限。
3 田间盆栽试验验证ACP81对盐碱土壤改良、作物促生实际效果。
4 发酵工艺优化，高盐条件下诱导菌株大量生产工业纤维素酶、β-淀粉酶。
5 对比不同来源枯草芽孢杆菌转录组，挖掘通用耐盐核心保守基因。
测量的数据及研究意义
1 Bioscreen梯度盐浓度生长动力学OD曲线（图1a），直观反映盐浓度越高菌株生长抑制越强，确定6% NaCl为亚致死胁迫浓度，明确2% NaCl为菌株最优生长盐环境，为后续胁迫实验设置提供标准条件。
2 透射电镜细胞形态观测数据（图1b），无盐条件菌体细长完整，6%盐胁迫下菌体变圆、细胞壁破损，证明高盐破坏细胞膜结构，菌株通过上调细胞壁修饰基因修复屏障。
3 各类水解酶活性定量数据（表1），6% NaCl组纤维素酶、β-淀粉酶活性显著高于对照组，中性蛋白酶降低，证实盐胁迫特异性激活多糖水解通路，葡萄糖积累是重要渗透保护策略。
4 抗氧化酶与丙二醛、自由基含量数据（表2），盐胁迫CAT、POD上升，SOD、氧自由基含量下降，揭示菌株差异化抗氧化平衡机制，减轻盐诱导氧化损伤。
5 全基因组圈图、基础注释数据（图2），菌株携带144个碳水水解相关CAZyme基因，是其多糖降解遗传基础，提供完整基因背景，为后续基因编辑提供序列支撑。
6 转录组差异基因韦恩、PCA、火山图（图3、图4），对照组、盐胁迫组转录谱区分明显，共932个显著差异基因，证明盐胁迫触发全局性转录重编程，大量代谢通路协同重塑。
7 GO、KEGG富集分析图（图5、图6），盐胁迫显著富集鞭毛、相容性溶质、PTS糖转运通路，锁定四类核心耐盐功能模块，解析多层次适应机制。
8 RT-qPCR基因相对表达数据（图7），betB、gbsB、malL、celB等基因盐胁迫上调数十倍，体外qPCR验证转录组测序结果可靠，核心耐糖、渗透调控基因真实响应盐胁迫。
结论
1 从竹笋加工废料分离获得枯草芽孢杆菌ACP81，最高耐受20% NaCl，6% NaCl为亚致死胁迫浓度，2% NaCl为最适生长盐度。
2 6%盐胁迫下菌株细胞形态改变，纤维素酶、β-淀粉酶活性显著提升，葡萄糖作为相容性溶质缓解渗透压力；CAT、POD抗氧化酶上调，平衡盐诱导氧化损伤。
3 菌株基因组含有丰富碳水水解基因，盐胁迫诱导celB、celC、malL多糖代谢基因高表达，同时上调鞭毛组装、细胞壁D-丙氨酸修饰、甜菜碱合成、PTS糖转运通路基因。
4 betB、gbsB是相容性溶质甜菜碱合成关键基因，上调后提升细胞渗透耐受能力；鞭毛基因高表达增强菌体运动，利于营养摄取。
5 ACP81兼具高耐盐与多糖水解能力，可应用于盐碱地生物改良、工业水解酶发酵生产，基因组与转录组数据为耐盐芽孢杆菌机制研究提供资源。
使用芬兰 Bioscreen仪器测量数据的研究意义
1 本实验采用芬兰Bioscreen全自动微生物生长分析仪，96孔板同步设置6组NaCl浓度梯度，37℃恒温持续振荡，全程自动定时采集OD600数值，无需人工分批取样稀释测浊，大幅减少人为操作带来的读数偏差与实验人力成本。
2 仪器可连续监测数十小时完整生长周期，清晰区分延滞期、对数期、稳定期生长差异，精准量化不同盐浓度对增殖的抑制程度，仅单点终点OD无法区分生长动态细微差别。
3 全部培养体系、接种初始菌量、孵育条件完全标准化，各组数据统计学平行度高，可精准筛选最优生长盐度与亚致死胁迫浓度，为后续转录组、酶活实验提供标准化胁迫处理条件。
4 Bioscreen生长曲线作为基础生理表型证据，与酶活、电镜、转录组、qPCR分子数据形成完整证据链，保证盐胁迫响应机制推导严谨可靠。
5 高通量多孔并行检测模式一次性完成多盐梯度平行实验，无需分多批次摇瓶培养，大幅缩短菌株耐盐特性初筛周期。