生物炭添加下枯草芽孢桿菌應對鎘脅迫的分子機制

丁靜，陳偉光，陳宇婷，何岸飛，姜晶，盛光遙

（蘇州科技大學環境科學與工程學院，江蘇 蘇州 215009）

鎘具有很高的毒性，不僅會破壞生物機能，影響生態系統的正常功能[1]，還會通過食物鏈在生物以及人體內富集，給人類健康帶來巨大風險[2]，因此對鎘污染水體進行修復至關重要。生物炭是由富碳生物質（如農副產品）在缺氧條件下通過熱加工（如熱解）而獲得的物質[3]，其孔隙和表面官能團可與鎘發生物理化學作用[4]，是一種新興、低成本的吸附劑。然而，生物炭對鎘的固定效果較差，在廢水長期處理過程中可能發生解吸，造成二次污染[5]?？莶菅挎邨U菌是一種廣泛分布于不同環境的革蘭氏陽性桿狀細菌，是原位生態水體中的代表菌種[6]。它被認為是一種良好的生物吸附劑，低成本及環境友好的特點更使其在鎘廢水修復領域具有不可替代的優勢[7-9]。然而，利用微生物吸附劑處理含鎘廢水存在難以分離和回收、微生物活性在生長初期易受鎘毒性影響等缺點[10]，從而限制了微生物吸附劑的應用。

有研究指出，生物炭是固定微生物的良好載體，不僅可以為微生物提供有利的生存環境，還能為微生物生長提供碳源、能源和礦物營養[11]。而微生物可以協助生物炭完成對污染物的長期固定，同時還可能將生物炭用作厭氧呼吸中的末端電子受體[12]。Yuan等[13]的研究也證實了生物炭可作為微生物的電子供體。因此，近年來有學者采用生物炭聯合枯草芽孢桿菌處理含鎘廢水[14-15]，該方法可以彌補單一生物炭和微生物處理時的缺陷，同時生物炭和枯草芽孢桿菌之間的相互作用可能會進一步加強對含鎘廢水的處理效果。然而，這些研究主要關注體系中鎘的吸附去除效果，并結合一系列的外在表征手段，從化學層面探索鎘、枯草芽孢桿菌和生物炭之間的相互作用過程，缺乏從微生物角度闡述其詳細的內在分子機制[16]。明確枯草芽孢桿菌在生物炭添加及重金屬脅迫環境下的生理響應和分子作用機制，對優化生物炭-枯草芽孢桿菌復合技術處理重金屬廢水至關重要。

微生物對外在環境的生理適應主要在轉錄水平上協調[17]，因此研究微生物在不同環境條件下基因轉錄表達的變化，可以理解環境條件對微生物的內在影響機制[18-20]。Yu 等[21]通過轉錄組分析系統研究了溶解氧對枯草芽孢桿菌遺傳調控和代謝的影響。轉錄組數據分析顯示，低氧供應對代謝有3 個主要影響：增強碳代謝（葡萄糖代謝、丙酮酸代謝和碳溢出）、抑制氮源降解（谷氨酸族氨基酸和黃嘌呤）和嘌呤合成。此外，參與能量、細胞類型分化、蛋白質合成的基因表達也受氧氣供應的影響。Oomes 等[22]通過轉錄組學研究了鈣離子對枯草芽孢桿菌孢子形成的影響，結果顯示，孢子形成、鞭毛形成和生物膜基質形成通路都受到了顯著影響，孢壁多糖生物合成基因以及生物膜形成相關基因的表達也被影響。此外，對于枯草芽孢桿菌暴露于多種營養和環境條件下平鋪陣列轉錄組的研究，為研究枯草芽孢桿菌與環境污染物的相互作用機制提供了強大的理論基礎[23]。

因此，本實驗利用轉錄組學測序，研究枯草芽孢桿菌在不同體系（正常、鎘脅迫、生物炭、生物炭-鎘脅迫）中的轉錄調控分子機理，討論生物炭添加對枯草芽孢桿菌抵抗鎘脅迫壓力及鎘去除能力的影響機制。研究結果可為未來生物炭-微生物聯用處理重金屬廢水提供理論支撐。

1 材料與方法

1.1 生物炭制備

生物炭原料為玉米秸稈，產自江蘇東海。將收集的玉米秸稈用粉碎機打碎（2～5 cm）后裝入砂鍋，用錫紙包裹，置于馬弗爐中，在500 ℃條件下高溫熱解3 h。待自然降溫6 h 后取出，在研缽中研磨，過100 目篩。洗凈烘干，最后冷凍干燥24 h后收集備用。

1.2 細菌來源及培養

本實驗采用的枯草芽孢桿菌購自中國普通微生物菌種保藏管理中心（CGMCC 編號1.4255）。培養基包括氯化鈉5 g·L-1、蛋白胨10 g·L-1、牛肉膏3 g·L-1，pH 7.0。在30 ℃、170 r·min-1下培養。

1.3 實驗體系搭建

實驗設置正常培養（BS）、生物炭添加（BS_C）、鎘脅迫（BS_Cd）、生物炭-鎘脅迫（BS_C_Cd）4 組體系。將對數期保藏的枯草芽孢桿菌按7%（V/V）的比例先后兩次接種于新鮮培養基中，依次完成活化（2 h）和擴大培養（8 h）?？紤]到枯草芽孢桿菌在生物炭添加和鎘脅迫環境中適宜的暴露時間并不一致，并且在BS_C_Cd體系中，需要將枯草芽孢桿菌和生物炭混合培養一段時間后再進入鎘脅迫環境，因此本實驗設置了兩組培養體系。將擴大培養的枯草芽孢桿菌等量（200 mL）分裝到兩組無菌錐形瓶中，選擇一組為BS體系，在另一組中加入生物炭（0.2 g），形成BS_C 體系。兩組體系均設置3 個生物學重復，并將其置于30 ℃和170 r·min-1的條件下繼續培養。培養3.5 h后，對兩組體系進行取樣，用于轉錄組測序。取樣結束后，迅速在兩組體系中加入滅菌后的CdCl2溶液，控制每組體系中的鎘濃度為50 mg·L-1，形成BS_Cd 體系和BS_C_Cd 體系，反應時間控制為20 min，用于分析鎘脅迫對枯草芽孢桿菌的即時壓力效應。在此期間，分別在0、10、20 min 時對兩組體系進行鎘濃度測定。反應結束時（20 min），取樣用于轉錄組測序。

1.4 鎘濃度測定

采用原子吸收光譜儀PinAAcle900T（PerkinElmer，Waltham，MA，USA）測定溶液中的鎘濃度。標準鎘溶液（1000 μg·mL-1）購自上海阿拉丁生化科技股份有限公司。樣品測定前使用0.45μm 濾膜過濾，并加入定量硝酸（1%）稀釋，以控制濃度在標準曲線范圍內，同時也抑制溶液中的鎘水解。

1.5 樣品收集及RNA分離純化

BS 體系和BS_Cd 體系：樣品采集后，利用離心機在4 000 r·min-1下離心3 min 收集菌體，用液氮快速冷凍后放入-80 ℃超低溫冰箱中保存。

BS_C體系和BS_C_Cd體系：樣品采集后，分兩次離心收集菌體。第一次以2 500 r·min-1的較低轉速離心，除去樣品中的生物炭。第二次以4 000 r·min-1轉速離心收集菌體，用液氮快速冷凍，放入-80 ℃超低溫冰箱中保存。

采用標準方法進行樣品總RNA 提取，通過瓊脂糖凝膠電泳分析RNA 完整性及是否存在DNA 污染，利用NanoPhotometer spectrophotometer（Implen，Munich，德國）檢測RNA的純度，通過Agilent 2100 bioanalyzer（Agilent，Palo Alto，CA，美國）精確檢測RNA 的完整性。

1.6 轉錄組學測序及分析

從提取的總RNA 中去除rRNA，獲得mRNA。隨后加入碎片化緩沖液，將得到的mRNA隨機打斷成短片段，按照鏈特異性方式建庫[24]。文庫構建完成后，使用Invitrogen Qubit 2.0 Fluorometer（ThermoFisher，Waltham，MA，美國）對其進行初步定量，然后稀釋文庫至1.5 ng·μL-1。隨后使用Agilent 2100 bioanalyzer（Agilent，Palo Alto，CA，美國）對文庫中的插入片段大小進行檢測。插入片段大小符合預期后，利用qRTPCR對文庫的有效濃度進行準確定量（文庫有效濃度高于2 nmol·L-1），并進行Illumina測序。

測序完成后，對獲得的轉錄組數據進行原始數據過濾，去除帶接頭的序列、含N（N 表示無法確定堿基信息）序列和低質量序列（Qphred≤20 的堿基數占整個序列長度50%以上），并進行測序錯誤率和C 含量分布檢查，最終獲得用于后續分析使用的高質量序列。用Bowtie2軟件對過濾后的序列進行基因組定位分析[25]。用Rockhopper軟件將測序結果根據參考基因組進行組裝，并與已注釋的基因模型進行比較，發現新的轉錄本區域[26]。通過Blastx與Nr庫比對，對新預測的轉錄本區域進行注釋。采用ClusterProfiler軟件對差異基因集進行GO功能富集和KEGG通路富集分析，找出與差異基因顯著性相關的生物學功能或通路[27]。

2 結果與討論

2.1 含鎘體系中的鎘去除情況

圖1 顯示，兩個含鎘體系（BS_Cd 和BS_C_Cd）均在0～10 min 階段即表現出明顯的鎘吸附，之后在10～20 min 階段吸附放緩。BS_C_Cd 體系中生物炭和細菌對鎘的復合吸附量遠大于BS_Cd 體系中單一細菌的吸附量，這可能會在一定程度上緩解體系中枯草芽孢桿菌的鎘脅迫毒性。因此本實驗將鎘脅迫的反應時間設定為20 min，盡量讓枯草芽孢桿菌在鎘脅迫環境中具有一定的暴露時間，同時控制兩個體系中的鎘濃度差，從而獲得較科學的鎘脅迫/生物炭添加對枯草芽孢桿菌生理代謝的影響機制。

圖1 BS_Cd和BS_C_Cd體系中的鎘去除情況Figure 1 Cadmium removal in BS_Cd and BS_C_Cd systems

2.2 枯草芽孢桿菌在不同體系中的基因表達差異

由于本實驗中設置了兩組培養體系（每個體系又分為兩個階段），為了更好地認識枯草芽孢桿菌受生物炭刺激活性后，其對鎘脅迫壓力的生理響應，本研究對BS_C_Cd 和BS_Cd 的轉錄組結果進行了對比分析。從圖2 中可以發現，生物炭和鎘的加入均會導致枯草芽孢桿菌基因表達差異。BS_CdvsBS 組共有1 932 個差異表達基因，其中上調基因967 個，下調基因965 個；BS_CvsBS 組有312 個差異表達基因，其中上調基因133 個，下調基因179 個；BS_C_CdvsBS_C組有1 722 個差異表達基因，其中上調基因834 個，下調基因888 個。除此之外，在對BS_C_Cd 和BS_Cd 樣品進行對比分析時發現，其存在1 541 個差異表達基因，其中上調基因691 個，下調基因850 個，差異基因數目遠高于BS_CvsBS 組，說明在鎘脅迫壓力下，生物炭對枯草芽孢桿菌生理代謝的影響要遠大于正常生長條件。BS_C_CdvsBS_C 組差異基因數目少于BS_CdvsBS 組，這進一步證明枯草芽孢桿菌受到生物炭的積極作用，緩解了外界鎘脅迫對其造成的壓力和影響。這可能是因為生物炭粗糙的表面孔道、較大的比表面積和孔隙容量為枯草芽孢桿菌提供了良好的生存環境，幫助枯草芽孢桿菌在一定程度上規避了鎘毒性[28]。同時在生物炭-鎘脅迫體系中，生物炭的存在為枯草芽孢桿菌的生長提供了營養[29]。Tao 等[30]以棉稈為原料制備生物炭作為枯草芽孢桿菌的高效接種載體，也發現與未添加生物炭的體系相比，細菌在生物炭體系中生長較快，細菌的最高濃度也較高?；鹕綀D直觀展示了每個組合的差異基因分布情況（圖2b～圖2f）。相比BS_CvsBS 組，其余有鎘添加的實驗組都表現出較高的基因表達差異及顯著性水平。這進一步說明當暴露在鎘脅迫壓力下時，枯草芽孢桿菌的生理代謝會受到明顯抑制。然而如前所述，可能由于生物炭能幫助枯草芽孢桿菌規避鎘毒性，并可為細菌生長提供營養，因此在細菌應對鎘脅迫的生理響應上產生明顯影響，從而導致其基因表達差異顯著。

圖2 各實驗/對照組間的差異表達基因數量（a）及其差異表達基因火山圖分布（P＜0.05，|log2差異倍數|＞0）（b）～（f）Figure 2 Numbers of differentially expressed genes（a）and volcano plot of the differentially expressed genes in the experimental/control groups（P＜0.05，|log2 Fold Change|＞0）（b）～（f）

2.3 枯草芽孢桿菌在不同體系中的GO功能富集

2.3.1 鎘脅迫/生物炭添加前后枯草芽孢桿菌重要功能條目表達變化

GO 功能富集揭示了各實驗/對照組重要功能條目相關基因的上/下調情況。BS_CdvsBS組的差異基因GO 功能富集分析結果如圖3（a）所示，BS_CdvsBS以BS_Cd 為實驗組、BS 為對照組。圖3（a）顯示，在鎘脅迫下，上調的差異表達基因主要歸類為蛋白質復合物（Protein-containing complex）、離子轉運（Ion transport）、細胞器（Organelle）、陽離子跨膜轉運蛋白活性（Cation transmembrane transporter activity）、運動（Locomotion）。這表明在鎘壓力下，枯草芽孢桿菌通過調控陽離子運輸和轉運基因的表達來加快溶液中鎘的去除。此外，差異基因在蛋白質復合物方面也出現了正向調控，特別是CotA基因表達上調。CotA 是枯草芽孢桿菌的芽孢外壁蛋白，可以保護細菌免受紫外線輻射和過氧化物的損害[31]，因此該生物功能的正向調控可能有助于細菌隔絕溶液中鎘的直接毒害。下調的差異基因主要歸類為氧化還原酶活性（Oxidoreductase activity）、肽酶活性（Peptidase activity，acting on L-amino acid peptides）、維生素生物合成過程（Vitamin biosynthetic process）。整體來看，GO 富集分析在細胞組分分類中的上調基因數量均多于下調基因，同時在生物學過程分類中的離子轉運和運動功能條目、分子功能中的陽離子跨膜轉運蛋白活性功能條目中，上調基因數量也多于下調基因。而在生物學過程分類中的維生素合成功能條目、分子功能分類中的氧化還原酶和水解酶活性功能條目中，下調基因數量要多于上調基因。此外，值得一提的是運動功能條目，其差異基因表達全部為上調，這意味著鎘壓力明顯增強了細胞運動功能。盡管枯草芽孢桿菌面對鎘毒性會產生上述自衛反應，但鎘脅迫對枯草芽孢桿菌的侵害仍然是直接且不可避免的，在鎘毒性下，枯草芽孢桿菌的部分生理代謝功能仍然受到了明顯抑制。

圖3 BS_Cd vs BS組（a）和BS_C vs BS組（b）的差異基因GO功能富集分析Figure 3 GO function enrichment analysis of differentially expressed genes in BS_Cd vs BS（a）and BS_C vs BS（b）groups

BS_CvsBS以BS_C為處理組、BS為對照組。圖3（b）顯示，添加生物炭后，GO 富集分析在細胞組分分類中上調基因數量均多于下調基因。上調的差異表達基因主要歸類為細胞（Cell）、胞內（Intracellular）、水解酶活性（Hydrolase activity，hydrolyzing O-glycosyl compounds）、催化復合體（Catalytic complex）、核糖核蛋白復合物的生物合成（Ribonucleoprotein complex biogenesis）。下調的差異基因主要歸類為各種生物大分子合成及代謝調控，包括RNA 合成和DNA 的轉錄過程，以及各種轉移酶活性。從增強的基因功能條目來看，生物炭的添加對枯草芽孢桿菌的細胞組分相關功能產生了積極影響，但對其生物合成和各種轉移酶活性產生了負面影響。

2.3.2 生物炭-鎘脅迫體系中枯草芽孢桿菌重要功能條目表達變化

BS_C_CdvsBS_C 是以BS_C_Cd 為處理組、BS_C為對照組。圖4（a）顯示，上調的差異基因包括細胞組分的各功能條目和分子功能下的各種底物/離子轉運功能條目，特別是轉運活性（Transporter activity）。此外，在生物學過程分類下的細胞運動（Cell motility）、趨化性（Chemotaxis）、纖毛或鞭毛依賴的細胞運動（Cilium or flagellum-dependent cell motility）、運動（Locomotion）、細胞或亞細胞組分的移動（Movement of cell or subcellular component）和趨避性（Taxis）功能條目中，差異基因表達全部為上調。下調的差異基因主要為維生素的生物合成和代謝過程，以及內肽酶活性。這表明相對于BS_C 體系，BS_C_Cd 體系中的枯草芽孢桿菌受到鎘脅迫時，其細胞內部的一些生物合成和代謝過程會受到抑制，但是細胞運動趨勢明顯增強（或許會向生物炭表面和內部運動），從而使其有效規避鎘毒性壓力。并且，細菌的離子轉運功能也得到明顯增強，這可能提高了枯草芽孢桿菌的鎘轉移能力。

圖4 BS_C_Cd vs BS_C組（a）和BS_C_Cd vs BS_Cd組（b）的差異基因GO功能富集分析Figure 4 GO function enrichment analysis of differentially expressed genes in BS_C_Cd vs BS_C（a）and BS_C_Cd vs BS_Cd（b）groups

如前所述，除上述3 組樣品間的差異基因分析外，本研究對BS_C_Cd 和BS_Cd 的轉錄組結果也進行了對比分析，用以參考性地認識生物炭刺激枯草芽孢桿菌活性后，枯草芽孢桿菌對鎘脅迫壓力的生理響應是否與正常狀態下的鎘應對措施有所不同。圖4（b）顯示，GO 富集分析在大部分功能條目中的上調基因數量都多于下調基因。上調的差異基因主要有小分子結合（Small molecule binding）、核苷酸結合（Nucleotide binding）、核苷磷酸結合（Nucleoside phosphate binding）、裂解酶活性（Lyase activity）、氧化還原酶活性（Oxidoreductase activity，acting on CH-OH group of donors）、膜的內在組分（Intrinsic component of membrane）、膜的必需組分（Integral component of membrane）、細胞（Cell）、細胞組分（Cell part）、細胞質（Cytoplasm）、胞內（Intracellular）、輔因子生物合成過程（Cofactor biosynthetic process）、輔因子代謝過程（Cofactor metabolic process）、單羧酸生物合成過程（Monocarboxylic acid biosynthetic process）、單羧酸代謝過程（Monocarboxylic acid metabolic process）、有機磷生物合成過程（Organophosphate biosynthetic process）、有機磷代謝過程（Organophosphate metabolic process）、含磷化合物代謝過程（Phosphate-containing compound metabolic process）、磷代謝過程（Phosphorus metabolic process）、四吡咯生物合成過程（Tetrapyrrole biosynthetic process）、四吡咯代謝過程（Tetrapyrrole metabolic process）。下調的差異基因有細胞外圍（Cell periphery）、質膜（Plasma membrane）和激酶活性（Kinase activity）?？梢钥吹?，相比于單一含鎘體系，生物炭的存在增強了枯草芽孢桿菌體內各種生物學過程、細胞組分和分子功能基因的表達，這與前文生物炭為枯草芽孢桿菌提供營養的說法相一致。張秋等[32]探究了生物炭對鎘污染土壤團聚體酶活性的影響，結果顯示，生物炭顯著提升了鎘污染土壤團聚體的酶活性，尤其是在碳循環酶中的蛋白酶和氧化還原酶中的過氧化氫酶處理上表現明顯。本研究認為，這可能是因為生物炭改良和提高了團聚體的結構和性能，形成了對微生物更有利的生存環境，同時對土壤中重金屬有效態的鈍化作用也促進了作物根系的生長，但對酶活性變化機理并未進行深入討論。生物炭也促進了細菌在鎘環境下對含磷化合物的合成和代謝。有研究指出[33]，微生物分泌的有機酸會在生物炭表面形成弱酸性的微環境，促進生物炭中磷的釋放，使重金屬與其形成穩定的磷化合物進而被長期固定在生物炭表面。此外值得注意的是，陽離子跨膜轉運蛋白相關基因czcD和重金屬外流泵蛋白相關基因cadA在BS_CdvsBS 組的表達均表現為下調，而在BS_C_Cd 和BS_Cd 的對比結果中則呈現為上調。czcD基因已被證明是金屬抗性基因[34]。并且研究發現，鎘的抗性基因系統主要有czc 和cad 兩種，cadA基因存在于包括葡萄球菌屬、微球菌屬及鹽芽孢桿菌屬等多種細菌基因組中，菌株通過該基因所編碼的酶將鎘泵出菌體外，從而起到解毒作用[35]。因此，czcD和cadA基因在BS_C_Cd 和BS_Cd 組中的表達上調，意味著生物炭可能激活了枯草芽孢桿菌體內鎘抗性基因的表達，從而增強了其對鎘脅迫的抵抗能力。并且膜相關功能條目的差異基因以上調為主。已有研究通過掃描電鏡觀察發現，在生物炭-鎘脅迫體系中確實形成了生物膜[15]。因此，本實驗中膜相關功能條目差異基因的上調，也進一步驗證了以上的實驗結果，預示著體系中早期生物膜發育跡象。

2.4 枯草芽孢桿菌在不同體系中的KEGG通路富集

2.4.1 鎘脅迫/生物炭添加前后枯草芽孢桿菌代謝通路表達變化

KEGG 富集分析結果表明（表1），體系中加入鎘后，枯草芽孢桿菌的鞭毛組裝和細菌趨化性通路都明顯增強。這兩條通路的增強是枯草芽孢桿菌在鎘脅迫下產生應激自衛反應的直接表現。鞭毛是細菌的運動器官，同時也在細菌的表面黏附上起到一定作用[36]。鎘的存在激活了枯草芽孢桿菌鞭毛運動的開關FliM，鞭毛組裝通路的增強說明在鎘環境中，枯草芽孢桿菌通過加快自身運動來躲避鎘毒性[37]。此外，氧化磷酸化通路也被增強。氧化磷酸化是細菌生成ATP 的重要途徑之一。相關研究發現，氧化磷酸化參與了細胞的抗毒害活動，通過促進細胞產能和提高能量轉化效率來解毒，證實了能量衰竭可能是造成毒害的原因[38]。而核黃素代謝和硫傳遞系統相關通路都受到了明顯抑制，說明細菌在鎘毒性壓力下部分物質代謝能力減弱，進一步體現了鎘毒性對枯草芽孢桿菌代謝活性的抑制。BS_CvsBS處理組的KEGG通路富集分析中，只有淀粉和蔗糖代謝通路表現出顯著富集，說明生物炭的添加顯著增強了枯草芽孢桿菌對能源物質的代謝功能。

表1 各處理/對照組KEGG顯著富集通路中的差異基因表達情況Table 1 Expression of the differential genes in KEGG enrichment pathways in the experimental/control groups

2.4.2 生物炭-鎘脅迫體系中枯草芽孢桿菌代謝通路表達變化

BS_C_CdvsBS_C 處理組的差異基因顯著富集通路情況與BS_CdvsBS 組相似，部分代謝相關通路（核黃素代謝、硫代謝、不同環境中的微生物代謝）的基因表達呈現下調，但與鞭毛組裝及細菌趨化性通路相關的基因表達則呈現上調。此外，核糖體和纈氨酸、亮氨酸及異亮氨酸生物合成通路的差異基因表達也表現出明顯上調?？梢婃k添加確實對細菌的代謝產生了部分抑制，但是會明顯增強細胞躲避壓力和毒性的運動行為。

而對BS_C_Cd 和BS_Cd 的差異基因進行顯著富集代謝通路分析時發現（表2），枯草芽孢桿菌物質和能源代謝相關通路的基因表達均呈現上調，包括氨基糖和核苷酸糖代謝、不同環境中的微生物代謝、碳代謝、糖酵解和糖異生以及磷酸戊糖途徑。這說明生物炭的存在可以緩解鎘脅迫對枯草芽孢桿菌的毒性效應，且可能為細菌提供了能源，繼而刺激了枯草芽孢桿菌的生命代謝活動。這也與周鳳等[39]的研究結果類似。

表2 BS_C_Cd和BS_Cd差異基因的KEGG顯著富集通路情況Table 2 KEGG enrichment pathway of differentially expressed genes between BS_C_Cd and BS_Cd

差異基因在氨基糖和核苷酸糖代謝、糖酵解和糖異生、磷酸戊糖途徑通路的顯著富集，也揭示了生物膜的早期發育現象。有研究指出，生物膜發育過程中，大多數被測的代謝途徑都發生了顯著的動態變化，包括TCA 循環、糖酵解、核苷酸和氨基酸生物合成等[40]。并且差異基因在細胞組分顯著富集的功能條目也與膜組成有關。而在BS_CdvsBS 組中，則沒有表現出任何生物膜形成跡象。已有研究發現，在生物炭協同體系中，枯草芽孢桿菌可以附著在生物炭表面及孔道內生長，形成復合體；并且菌體分泌出更多的胞外聚合物，使得細菌間粘連現象更明顯，這可能是形成生物膜的重要標志[15]。生物膜附著在固體載體上，對機械應力和生物量具有一定的保護作用，是一種有潛力的工業廢水處理方法[41]。因此生物膜的形成可以對細菌產生保護作用，這對于提高枯草芽孢桿菌對含鎘廢水的處理能力具有重要意義。

此外，在不顯著富集通路中，與谷胱甘肽代謝相關的差異基因表達呈現上調，但在BS_CdvsBS 組中，相關的差異基因表達則呈現下調。谷胱甘肽具有抗氧化和解毒功能，可防止活性氧物質如自由基、過氧化物等對細胞造成損害[42]。并且有研究指出，谷胱甘肽在金屬結合肽的合成中至關重要，金屬結合肽可通過形成細胞內穩定的金屬絡合物來螯合重金屬[43-44]。因此，該差異基因上調說明生物炭的存在提高了谷胱甘肽代謝活性，在提高細菌抗氧化和解毒能力的同時，可能也對鎘的去除起到了積極作用。嘌呤代謝、蛋白質輸出、細菌分泌系統這幾個通路的基因表達呈現下調，這可能是由于在BS_Cd 體系中，枯草芽孢桿菌直接暴露于鎘脅迫壓力下，需要更多刺激自身的免疫和保護功能。而在BS_C_Cd 體系中，受到生物炭的外在物理性保護和內在代謝促進，枯草芽孢桿菌受到的鎘脅迫壓力相對較小，因此對這些通路的基因表達刺激較小。

3 結論

（1）生物炭添加和鎘脅迫均會導致枯草芽孢桿菌基因表達差異，但生物炭可以緩解鎘脅迫對枯草芽孢桿菌的壓力和影響。

（2）在鎘脅迫壓力下，枯草芽孢桿菌的生理代謝功能受到了明顯抑制，但細菌通過上調細胞運動及鎘轉運相關基因的表達來規避鎘毒性，其鞭毛組裝和細菌趨化性通路都明顯增強。

（3）生物炭添加對枯草芽孢桿菌的細胞組分相關功能產生了積極影響，并且增強了細菌對能源物質的代謝功能。

（4）生物炭-鎘脅迫體系中，生物炭增強了枯草芽孢桿菌在鎘脅迫壓力下各種生物過程、細胞組分和分子功能基因的表達，包括細菌對含磷化合物的合成和代謝，以及細菌體內鎘抗性基因的表達。