摘要：為了闡明甜椒響應(yīng)低溫脅迫的分子機(jī)制，對低溫脅迫時甜椒幼苗的 MDA 含量、SOD 活性以及光合作用相關(guān)參數(shù)進(jìn)行了測定，并結(jié)合轉(zhuǎn)錄組測序分析。結(jié)果表明：低溫脅迫3d 時甜椒 MDA 含量顯著升高，SOD 活性顯著降低，光合相關(guān)參數(shù)顯著下降；與低溫處理前相比，轉(zhuǎn)錄組測序分析結(jié)果顯示檢測出差異基因主要涉及酶系統(tǒng)、光合作用、信號傳導(dǎo)以及相關(guān)轉(zhuǎn)錄因子，其中SOD相關(guān)基因CaSOD呈示上調(diào)表達(dá)，光合作用關(guān)鍵基因 CaCP4與CaHCR的表達(dá)量則顯著下調(diào)。研究結(jié)果為深入解析甜椒幼苗響應(yīng)低溫脅迫的分子機(jī)制及甜椒抗低溫育種奠定了基礎(chǔ)。

關(guān)鍵詞：甜椒；低溫脅迫；轉(zhuǎn)錄組；差異表達(dá)基因

中圖分類號：S641.3文獻(xiàn)標(biāo)志碼：A文章編號：0253?2301（2024）04?0006?11

DOI：10.13651/j.cnki.fjnykj.2024.04.002

Physiological Response and Transcriptome Analysis of Sweet Pepper SeedlingsUnder Low Temperature Stress

MA Xin-yi1 ，LIN Yi-hui1 ，YU Song-jin1 ，NI Hui2 ，LIN Wei3 ，WEI Mi1 ，CHEN Han-xin1*

（1. Zhangzhou Institute of Agricultural Sciences， Zhangzhou， Fujian 363005， China;2. School of Marine Food andBiotechnology， Jimei University， Xiamen， Fujian 361021， China;3. Subtropical Agriculture Research Institute，F(xiàn)ujian Academy of Agricultural Sciences， Zhangzhou， Fujian 363005， China）

Abstract： In order to clarify the molecular mechanism of sweet pepper in response to low temperature stress， the MDA content， SOD activity of sweet pepper seedlings and the relative parameters of photosynthesis under the low temperature stress were determined， and the transcriptome was sequenced and analyzed. The results showed that the contentofMDAinsweetpepperincreasedsignificantly，theactivityofSODdecreasedsignificantly，andthe photosynthetic relative parameters decreased significantly under the low temperature stress for 3 days. Compared with those before the treatment of low temperature， the transcriptome sequencing analysis showed that the differentially expressedgenesweremainlyinvolvedintheenzymesystem，photosynthesis，signaltransductionandrelated transcription factors. Among them， the expression of SOD-related gene CaSOD was up-regulated， and the expression of photosynthesis key genes CaCP4 and CaHCR was significantly down-regulated. The results laid a foundation for the further analysis of the molecular mechanism of sweet pepper seedlings in response to the low temperature stress and the anti-low temperature breeding of sweet pepper.

Key words： Sweet pepper；Low temperature stress；Transcriptome；Differentially expressed genes

低溫脅迫作為植物常見的非生物脅迫之一，對其生長與發(fā)育有著較為嚴(yán)重的影響。植物在遭遇低溫脅迫時，葉片常常表現(xiàn)為皺縮萎蔫、光合效率降低等特征，植株則整體表現(xiàn)出生長受阻與產(chǎn)量下降，對農(nóng)業(yè)生產(chǎn)影響深遠(yuǎn)[1]。相關(guān)研究常以超氧化物歧化酶（SOD）、丙二醛（MDA）等生理生化指標(biāo)作為植物面對低溫脅迫的應(yīng)答。低溫下植物活性氧代謝系統(tǒng)受到影響，超氧陰離子自由基（O2?）等活性氧（ROS）含量增加，因此作為活性氧清除劑的 SOD 活性下降[2]。與此同時，植物在逆境脅迫下 MDA 含量會提高[3]。此外，對紫花苜蓿等植物的研究結(jié)果表明，低溫可導(dǎo)致1，5-二磷酸核酮糖羧化酶（RuBPcase）等光合反應(yīng)的相關(guān)酶降解或活性降低，從而影響植物的凈光合速率[1]。但是，有關(guān)低溫脅迫對植物生長影響的信號傳導(dǎo)途徑及基因表達(dá)調(diào)控機(jī)制還不明確。

近年來，高通量測序技術(shù)逐漸成為研究植物應(yīng)對低溫脅迫的分子調(diào)控機(jī)制的普遍技術(shù)途徑。相關(guān)研究采用組學(xué)技術(shù)揭示了低溫信號通過植物激素或者 Ca2+及 ROS 等第二信使向下傳遞。ROS 會激活包括 MAPK 級聯(lián)、鈣依賴性蛋白激酶（CDPKs）、鈣調(diào)神經(jīng)磷酸酶 B 型蛋白（CBLs）、與 CBL 相互作用的蛋白激酶（CIPK）等蛋白激酶（PKs）以及蛋白磷酸酶（PPs）。隨后，PKs 和 PPs傳遞的信息至下游并引發(fā)一系列包括轉(zhuǎn)錄因子（TF）在內(nèi)的磷酸化或去磷酸化級聯(lián)，最終完成植物在低溫逆境中信號傳導(dǎo)和與低溫脅迫有關(guān)基因的表達(dá)調(diào)控[4]。

甜椒 Capsicum annuum var.grossumL.又名青椒、菜椒、燈籠椒，屬于茄科 Solanaceae辣椒屬 Capsicum ，富含豐富的維生素 C 、多酚以及β-胡蘿卜素等[5]，是我國常見的設(shè)施栽培作物。甜椒屬于冷敏感作物，當(dāng)溫度低于7℃時，可導(dǎo)致葉表面點(diǎn)蝕、葉片發(fā)黃等生理現(xiàn)象，進(jìn)而嚴(yán)重影響其產(chǎn)量[6]。目前，對于甜椒低溫脅迫的研究多集中在生理生化等方面[7]。近年來已有部分針對甜椒在低溫逆境下轉(zhuǎn)錄組測序的研究。例如，Li 等[8]利用 RNA 測序技術(shù)研究甜椒幼苗在低溫脅迫應(yīng)答時發(fā)現(xiàn)鈣離子（Ca2+）信號傳導(dǎo)相關(guān)基因如 CBL 等下調(diào)，多數(shù)轉(zhuǎn)錄因子如 WRKY 等也下調(diào)。而 Shin 等[9]卻發(fā)現(xiàn)甜椒在低溫時次生代謝過程中的一些轉(zhuǎn)錄因子如 WRKY40、MKK9、MYB77等顯示上調(diào)，同時相關(guān)低溫脅迫響應(yīng)基因如 POD 等也顯著上調(diào)。Kong 等[10]測序發(fā)現(xiàn)甜椒果實(shí)中 AP2、MYB 、bHLH、WRKY 等轉(zhuǎn)錄因子也參與了低溫反應(yīng)。然而，甜椒整體響應(yīng)低溫脅迫機(jī)制以及關(guān)鍵基因尚不充分。因此，本研究以冷敏感甜椒品種博頓方椒 F1為試驗(yàn)材料，基于轉(zhuǎn)錄組測序技術(shù)結(jié)合生理生化指標(biāo)對低溫脅迫下的甜椒基因調(diào)節(jié)進(jìn)行深入研究，探究甜椒及其他植物低溫脅迫的分子響應(yīng)機(jī)制。

1材料與方法

1.1試驗(yàn)材料

試驗(yàn)甜椒品種博頓方椒 F1為冷敏感品種，購自漳州美利德生物工程有限公司。

1.2試劑與儀器

主要試劑：液氮購于漳州市新興氣體有限公司；植物丙二醛（MDA）測試盒購于南京建成生物工程研究所；超氧化物歧化酶（SOD）活性檢測試劑盒購于北京索萊寶科技有限公司

主要儀器：人工氣候箱（RON-500C ，寧波揚(yáng)輝儀器有限公司）；分光光度計(jì)（L5S ，上海儀電分析儀器有限公司）；便攜式光合作用儀（GFS-3000，德國 WALZ）。

1.3試驗(yàn)方法

低溫脅迫試驗(yàn)于2023年在漳州市農(nóng)業(yè)科學(xué)研究所進(jìn)行。種子經(jīng)過溫水浸泡催芽后，播種于50孔育苗盤內(nèi)。當(dāng)幼苗長至6～7片真葉時定植于小方盆中，并放置于人工氣候箱，以溫度28℃/18℃（晝/夜）、光照強(qiáng)度240μmol? m?2?s?1、濕度70%的條件下緩苗3 d 。之后進(jìn)行低溫處理，光照強(qiáng)度與濕度不變，溫度控制在15℃/5℃（晝/夜），處理7 d 。分別在處理后的0、3、5、7 d 時取功能葉（從上往下數(shù)第3～4葉）測定生理生化指標(biāo)以及轉(zhuǎn)錄組測序，每次隨機(jī)選取10株取樣，重復(fù)3次。所有樣品在采樣后均在液氨中快速冷凍，并于?80℃中保存。

1.3.1生理指標(biāo)測定使用植物丙二醛（MDA）測試盒和超氧化物歧化酶（SOD）活性檢測試劑盒測定甜椒葉片在低溫處理0、3、5、7 d 時 MDA 含量及 SOD 活性，方法按照試劑盒說明書進(jìn)行測定。每個處理重復(fù)3次。

1.3.2葉片光合氣體交換參數(shù)測定利用便攜式光合作用儀 GFS-3000對功能葉片的凈光合速率、蒸騰速率、氣孔導(dǎo)度、胞間 CO2濃度進(jìn)行測定。測定時間為08：00～11：30，測定光照度為800μmol·m?2·s?1，流速為500μmol·s?1， CO2濃度為400μmol·mol?1。每個處理重復(fù)3次。

1.3.3轉(zhuǎn)錄組測序委托深圳華大基因科技服務(wù)有限公司對0 d（ck）與3 d（cold）低溫處理的甜椒葉片，共6個樣品（各3個重復(fù)，分別為 ck1、ck2、 ck3、cold1、cold2、cold3）中總 RNA 進(jìn)行提取，用 Agilent 2100檢測提取結(jié)果，RIN（RNA integ- rity， RIN）≥8，結(jié)果 A 類；基于 DNBSEQ 平臺進(jìn)行轉(zhuǎn)錄組測序并對原始數(shù)據(jù)去冗余過濾。

1.3.4轉(zhuǎn)錄組差異表達(dá)基因qRT-PCR 驗(yàn)證為驗(yàn)證轉(zhuǎn)錄組測序的準(zhǔn)確性，本研究選擇了20個差異表達(dá)的低溫相關(guān)基因進(jìn)行熒光定量 PCR 驗(yàn)證。基于基因組數(shù)據(jù)庫獲得的基因全長信息，使用 primer 5.0設(shè)計(jì)熒光定量特異引物（表1），內(nèi)參基因?yàn)?ACTIN ，qRT-PCR 試劑盒（湖南艾科瑞生物工程有限公司）反應(yīng)體系包括10μL 2×SYBR GreenPro Tap HS premix ，1μL 上下游引物，1μL 稀釋后的 cDNA 模板（100 ng·μL?1），無菌水補(bǔ)充至20μL，反應(yīng)步驟：95℃變性30 s， 95℃退火5 s ，60℃延伸20 s ，40個循環(huán)。最后以2?△△CT 計(jì)算差異表達(dá)基因的相對表達(dá)量。

1.4數(shù)據(jù)分析

1.4.1生理及光合參數(shù)數(shù)據(jù)分析使用 R 語言對數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析，使用 Origin 2021及 GraphPad Prism 10.1.2繪制柱形圖。

1.4.2轉(zhuǎn)錄組數(shù)據(jù)分析使用過濾軟件 SOAP- nuke 過濾掉低質(zhì)量、接頭污染以及未知堿基 N 含量大于5%的 reads ，得到 clean reads 。不使用參考基因，使用 Trinity 對 clean reads 進(jìn)行組裝得到Unigene，用 BUSCO 對組裝的轉(zhuǎn)錄本進(jìn)行質(zhì)量評估。之后對Unigene進(jìn)行功能注釋、SSR 檢測，之后每個樣品在 All-Unigene的基礎(chǔ)上使用 Bowtie2將 clean reads 比對到基因組序列上，然后使用RSEM 計(jì)算各個樣品的基因表達(dá)水平計(jì)算表達(dá)量。使用 DESeq2分析基因差異表達(dá)顯著性，并使用 R 軟件中的phyper對差異表達(dá)基因（DEGs）做 GO 和 KEGG 富集分析，以Qvalue≤0.05的功能視為顯著富集。最后使用 R 軟件的 ggplot2繪制基因統(tǒng)計(jì)柱狀圖、火山圖、氣泡圖、GO 及 KEGG 富集圖。

2結(jié)果與分析

2.1低溫脅迫對甜椒葉片 MDA 含量和 SOD活性的影響

MDA 作為膜脂過氧化的最終產(chǎn)物，是衡量植物低溫逆境下細(xì)胞膜受損程度的常用理化指標(biāo)之一。MDA 含量測定發(fā)現(xiàn)（圖1a），在3 d 時甜椒幼苗葉片 MDA 含量與處理前相比顯著提高，在5 d 和7 d 時又降低，并趨于穩(wěn)定。孟雅寧等[11]發(fā)現(xiàn)甜椒低溫脅迫時 MDA 含量在7 d 后才開始下降，而劉凱歌等[12]則發(fā)現(xiàn)在10 d 。本研究發(fā)現(xiàn) MDA 含量在低溫脅迫后迅速上升，說明甜椒產(chǎn)生了低溫脅迫損傷，這與劉凱歌等研究結(jié)果相似。

植物在面對低溫等逆境脅迫時，葉片中的 ROS 含量會提高，攻擊其細(xì)胞膜，引起膜脂過氧化[13]。而 SOD 作為生物體內(nèi)的活性氧清除劑，在植物面臨逆境脅迫時處于重要地位，其活性隨著 ROS 含量的增高而下降[14]。本試驗(yàn) SOD 活性分析（圖1b）表明，葉片中 SOD 在低溫脅迫后先下降，在3 d 時甜椒幼苗葉片 SOD 活性與處理前相比顯著下降，在7 d 時 SOD 活性恢復(fù)至處理前水平。本研究顯示葉片中 SOD 在低溫脅迫后先下降，說明產(chǎn)生了低溫脅迫損傷；后期升高可能是 SOD 合成量增加。

2.2低溫脅迫對甜椒葉片光合作用的影響

對低溫肋迫下甜椒葉片的光合相關(guān)參數(shù)測定（圖2）發(fā)現(xiàn)，低溫脅迫下，甜椒葉片的凈光合速率、蒸騰速率、氣孔導(dǎo)度及胞間 CO2含量均發(fā)生顯著下降。其中，凈光合速率在3 d 時顯著降低，5 d 時小幅提高，但7 d 時又降至3d 時的水平；胞間 CO2含量在3 d 時顯著降低，5 d 時開始緩慢上升，但是7 d 時仍顯著低于處理前；蒸騰速率與氣孔導(dǎo)度在3 d 時顯著下降后一直保持在同一地位水平，并無明顯回升。莫億偉等[15]認(rèn)為低溫逆境下，光合速率主要受到相關(guān)酶活的影響，RuBPcase活性降低導(dǎo)致 ATP 供應(yīng)不足，暗反應(yīng)受到影響；果糖-1，6-二磷酸酯酶（FBPase）在低溫下也受到抑制，使電子傳遞速率降低；參與卡爾文循環(huán)的其他酶活也隨之降低，影響葉綠素含量及 CO2固定能力，最終影響光合速率。本試驗(yàn)中低溫脅迫對甜椒葉片的光合作用的影響較大，無法自行調(diào)節(jié)回到常溫狀態(tài)下的水平，這與相關(guān)研究報(bào)道的低溫脅迫損傷導(dǎo)致光合作用下降的現(xiàn)象具有相似性。

2.3甜椒低溫脅迫轉(zhuǎn)錄組結(jié)果分析

2.3.1低溫脅迫下甜椒轉(zhuǎn)錄組數(shù)據(jù)質(zhì)控及序列比對分析采用低溫處理0 d（ck）及3 d（cold）的甜椒葉片共6個樣本進(jìn)行轉(zhuǎn)錄組測序（表2），一共測得38.48 Gb 數(shù)據(jù)。各樣品過濾后堿基數(shù)均達(dá)到6.35 Gb 以上，Q20堿基百分比均在98.09%以上， Q30堿基百分比均在93.37%以上，過濾后 Reads 比例均在96.9%以上。由表3可知，對 clean reads 進(jìn)行組裝并去冗余后得到163349個Unigene，總長度、平均長度、N50以及 GC 含量分別為168775018 bp 、1033 bp 、1804 bp 和39.02%。以上指標(biāo)數(shù)據(jù)說明本次測序質(zhì)量良好，可進(jìn)行下一步分析。將 clean reads 對比到參考基因序列得到對比序列，總對比率達(dá)90.13%以上。

2.3.2低溫脅迫下甜椒轉(zhuǎn)錄組差異表達(dá)基因篩選

以 P＜0.05、Qvalue≤0.05、|Log2FC|≥1為篩選條件，進(jìn)一步對 ck 、cold 組間的 DEGs進(jìn)行統(tǒng)計(jì)分析。結(jié)果顯示，本次轉(zhuǎn)錄組測序 DEGs共有21862個，其中12156個 DEGs 上調(diào)（55.6%），9706個 DEGs 下調(diào)（44.4%）（圖3）。從本次篩選出的 DEGs 中在蛋白/酶系統(tǒng)、光合作用、信號傳導(dǎo)與轉(zhuǎn)錄因子等幾個方面進(jìn)一步篩選出低溫脅迫相關(guān)基因（表4）。在酶相關(guān) DEGs 中， SOD 差異表達(dá)基因CaSOD（ID ：DN1491）、過氧化物酶（POD）基因（ID： DN338、DN39066、DN9372）、椒谷胱甘肽轉(zhuǎn)移酶（GST）基因（ID： DN1851，ID ： DN340、DN7501等）、肌醇半乳糖苷合成酶（GolS）（ID：DN13841、 DN2302）均呈現(xiàn)上調(diào)表達(dá)上調(diào)；在光合作用相關(guān) DEGs 中，甜椒葉綠素 a/b 結(jié)合蛋白（CAB）4基因 CaCP4（ID ：DN25）與7-羥甲基葉綠素 a 還原酶（HCR）基因CaHCR（ID ： DN1414）呈現(xiàn)下調(diào)表達(dá)；信號傳導(dǎo)相關(guān) DEGs 中，鈣依賴性蛋白激酶（CDPK）基因（ID ：DN-10103、DN11238）、絲裂原活化蛋白激酶（MAPK）基因（ID：DN124、DN895、DN1474等）呈現(xiàn)上調(diào)表達(dá)；甜椒類鈣調(diào)磷酸酶亞基 B 蛋白質(zhì)（CBL）基因 CaCBL1（ID： DN3892）與 CaCBL3（ID ： DN2025）上調(diào)；轉(zhuǎn)錄因子作為植物脅迫的主要調(diào)節(jié)因子也在本試驗(yàn)中篩選出大量 DEGs ，其中甜椒 AP2基因 CaAP2（ID：DN11723）等上調(diào)，而 CaAP2 ANT（ID： DN790）下調(diào)；甜椒 MYB 基因（ID ： DN2522、DN1151等）上調(diào)；甜椒bHLH轉(zhuǎn)錄因子基因（ID： DN52、DN303、DN433等）上調(diào)；甜椒 WRKY 基因（ID ：DN1758、DN7795等）上調(diào)；甜椒 CAMTA 轉(zhuǎn)錄因子全部基因（ID：DN415、 DN5968、DN395）上調(diào)。

2.3.3低溫脅迫下甜椒轉(zhuǎn)錄組差異基因 GO 、KEGG 富集分析將篩選出的Unigene進(jìn)行 NR 注釋，統(tǒng)計(jì)不同的物種比例。由圖4可知，注釋上共計(jì)70.40%的基因?qū)儆谔鸾?辣椒（Capsicum annuum）。將在 NR 數(shù)據(jù)庫上比對的Unigene結(jié)果注釋到 GO 數(shù)據(jù)庫3個方面：生物過程（Biological Process）、細(xì)胞組成（CellularComponent）、分子功能（Molecular Function）的分類圖，見圖5。其中，在生物過程中，細(xì)胞過程（cellular process）與代謝過程（metabolic process）這兩個亞類占比最大；在細(xì)胞組成中，僅有細(xì)胞解剖實(shí)體（cellular anato- mical entity）與蛋白質(zhì)復(fù)合物（protein-containingcomplex）這兩個亞類；在分子功能中，結(jié)合（binding）與催化活性（catalytic activity）這兩個亞類占比最大。對比 DEGs 與 KEGG 數(shù)據(jù)庫，富集結(jié)果見圖6。結(jié)果顯示，本試驗(yàn) DEGs 主要富集在碳代謝（Carbon metabolism）途徑、植物激素信號轉(zhuǎn)導(dǎo)（Plant hormone signal transduction）、氨基酸合成（Biosynthesis of amino acids）途徑、輔助因子生物合成（Biosynthesis of cofactors）途徑等。

2.3.4低溫脅迫下甜椒轉(zhuǎn)錄組差異表達(dá)基因qRT- PCR 驗(yàn)證為了驗(yàn)證轉(zhuǎn)錄組數(shù)據(jù)的準(zhǔn)確性，對20個低溫相關(guān) DEGs進(jìn)行qRT-PCR 驗(yàn)證，結(jié)果表明這20個基因在甜椒葉片低溫脅迫中的變化趨勢與轉(zhuǎn)錄組的測序結(jié)果基本相符（圖7）。

3結(jié)論與討論

在甜椒生物酶系統(tǒng)中，SOD 及 POD 的 DEGs- CaSOD與CaPOD在低溫脅迫時上調(diào)，與以往試驗(yàn)中水稻非生理脅迫[16]與甜椒低溫脅迫[9]結(jié)果相符。這是因?yàn)橹参镌谟龅降蜏貢r會引起其氧化應(yīng)激產(chǎn)生過量的如 O2?、OH?等 ROS ，導(dǎo)致細(xì)胞損傷。而為防止其引起的氧化損傷，植物進(jìn)化出了一些酶或者非酶的系統(tǒng)，其中就包括 SOD 及 POD 等[9]。當(dāng) O2?過量產(chǎn)生時， SOD 能將其轉(zhuǎn)化為 H2O2， POD 能進(jìn)一步的使其轉(zhuǎn)化為水[17]。因此，本試驗(yàn)中上調(diào)表達(dá)的CaSOD與CaPOD證明了甜椒在低溫脅迫時已啟動相關(guān)保護(hù)機(jī)制。本研究中甜椒在低溫環(huán)境下光合作用受到較大影響，而轉(zhuǎn)錄組結(jié)果也證實(shí)了這一點(diǎn)。甜椒葉綠素 a/b 結(jié)合蛋白（CAB）作為植物葉綠體中的重要組成成分之一，具有吸收光能并將其能量傳遞至光反應(yīng) I（PSI）與光反應(yīng) II（PSII）中等重要作用[18]。本次測序檢測到甜椒中 CAB 的調(diào)控基因 CaCP4表達(dá)顯著下調(diào)，與以往研究顯示 CAB 在低溫脅迫下調(diào)結(jié)果相似[19?20]，說明甜椒葉片在低溫脅迫時光合作用的光能吸收與傳遞受到影響。與此同時，植物進(jìn)行光合作用時，葉綠素 a 能夠吸收光能并將其轉(zhuǎn)化成化學(xué)能，但是其本身不直接參與光合作用，需要被還原成一種更加活躍的酶，而 HCR 就是植物體內(nèi)主要還原葉綠素 a 的酶[21]。作為光合作用中關(guān)鍵酶 HCR 的調(diào)控基因，CaHCR在本次測序結(jié)果中顯著下調(diào)，說明甜椒葉片中光合反應(yīng)中的電子傳遞鏈?zhǔn)茏瑁M(jìn)而影響了光合作用中的化學(xué)反應(yīng)。

甜椒在應(yīng)對低溫脅迫時，其對低溫信號的傳導(dǎo)也是脅迫應(yīng)答的重要環(huán)節(jié)之一。此次轉(zhuǎn)錄組測序結(jié)果表明，甜椒低溫脅迫時多數(shù)信號傳導(dǎo)基因上調(diào)。其中，CDPK 與 CBL 作為鈣結(jié)合蛋白，在植物低溫脅迫下鈣傳遞通路中發(fā)揮重要作用[22]。大米中 OsCDPK7過表達(dá)可提高其在低溫等逆境下的抗性[23]，這與本試驗(yàn)中 CaCDPK7上調(diào)結(jié)果相似，說明甜椒中 CDPK 對于傳遞低溫信號有著重要的作用。除 CDPK 外，甜椒 CBL 調(diào)控基因 CaCBL1過表達(dá)。根據(jù)Conglin等[24]研究結(jié)果，CBL 相互作用蛋白激酶7（CIPK7）與 CBL1一同作用于低溫信號傳導(dǎo)，而 CBL1的功能缺失可導(dǎo)致植物對低溫的敏感以及抗寒能力的降低。本次實(shí)驗(yàn)發(fā)現(xiàn)甜椒中 CaCBL1上調(diào)，說明低溫脅迫時已啟動相關(guān)抗寒基因。

轉(zhuǎn)錄因子（TF）是植物應(yīng)對生物與非生物脅迫的重要調(diào)節(jié)因子[25]。在本試驗(yàn)中，甜椒轉(zhuǎn)錄因子bHLH、WRKY 與 CAMTA 等調(diào)控基因多數(shù)上調(diào)。其中，bHLH是 CBF1與 CBF2的調(diào)控因子，廣泛存在于各種真核生物中[4， 24]。研究發(fā)現(xiàn)，bHLH-92、bHLH49與 AabHLH35在植物非生物脅迫中起到了重要作用[26]。這與本試驗(yàn)篩選出的 CabHLH-92、CabHLH35、CabHLH49等在低溫環(huán)境中上調(diào)結(jié)果一致，證明了其在甜椒低溫脅迫的作用。而 WRKY 作為綠色植物中常見的 TF ，常參與到植物的非生物脅迫中[4]。以往研究中， AtWRKY15、 WRKY25被證實(shí)在擬南芥非生物脅迫中作用顯著[27?28]，本研究中多數(shù) WRKY 上調(diào)，與以往研究基本相似。除此之外，本試驗(yàn)還檢測出CaCA- MTA2、CaCAMTA3與 CaCAMTA4這3個 CAMTA 轉(zhuǎn)錄因子基因在低溫時過表達(dá)。在擬南芥中這3種 CAMTA 調(diào)控基因都對低溫做出迅速反應(yīng)， CAMTA3更是被證明能快速對 DREB1B 做出調(diào)控，而 DREB 也被認(rèn)為是低溫反應(yīng)轉(zhuǎn)錄級聯(lián)的第一步[29]，由此可以看出 CAMTA 在植物尤其是甜椒低溫脅迫反應(yīng)中的重要性。

綜上，本試驗(yàn)通過生理生化及轉(zhuǎn)錄組測序技術(shù)，比較分析甜椒幼苗在低溫脅迫下的理化指標(biāo)的變化及相關(guān)差異基因表達(dá)。理化試驗(yàn)從側(cè)面證明了低溫導(dǎo)致甜椒脂膜過氧化，活性氧含量增高，細(xì)胞膜受到損傷。轉(zhuǎn)錄組測序分析結(jié)果表明，低溫促進(jìn)了CaSOD的表達(dá)，消除了 ROS 等對細(xì)胞的毒性作用，提高了甜椒在低溫環(huán)境下的適應(yīng)性。轉(zhuǎn)錄組篩選出的光合作用關(guān)鍵基因 CaCP4與CaHCR顯著下調(diào)，說明光合作用的下降可能是由于光能吸收與能量傳遞受阻。除上述差異基因外，轉(zhuǎn)錄組測序還發(fā)現(xiàn)甜椒在低溫脅迫時其信號傳導(dǎo)與轉(zhuǎn)錄因子中均存在大量差異性表達(dá)基因，且多數(shù)顯著上調(diào)。本研究為甜椒低溫脅迫下生理代謝及防御機(jī)制提供了更加全面的見解，為今后鑒定其低溫脅迫過程中的關(guān)鍵基因和損傷機(jī)制提供理論基礎(chǔ)。

參考文獻(xiàn)：

[1]苗宇，趙堯堯，張麗佳，等.低溫脅迫后紫花苜蓿葉片光合作用的恢復(fù)特性研究［J］.中國農(nóng)業(yè)科技導(dǎo)報(bào)，2024，26（2）：80?89.

[2]李琬，項(xiàng)洪濤，何寧，等.初花期冷害對綠豆葉片生理的損傷及烯效唑的緩解效應(yīng)［J］.西北農(nóng)業(yè)學(xué)報(bào)，2023，32（10）：1566?1578.

[3]邱雨后，陳金玲，謝敏譞，等.堿脅迫對不同基因型培忠杉苗木生理生化特性的影響［J］.中南林業(yè)科技大學(xué)學(xué)報(bào)，2024（1）：97?108.

[4] KHAN S A ，LI M Z ，WANG S M ， et al. Revisiting the Role of Plant Transcription Factors in the Battle against Abiotic Stress［J］. International journal of molecular sciences，2018，19（6）：1634.

[5] Marín A ，F(xiàn)erreres F ，Tomás-Barberán F A ，et al. CharacterizationandQuantitationofAntioxidantConstituentsofSweetPepper （Capsicum annuum L.）［J］. J Agric Food Chem ，2004，52：3861?3869.

[6]王洪濤，艾希珍，秦競，等.幾個甜椒砧木品種對低溫弱光耐受性的差異［J］.西北農(nóng)業(yè)學(xué)報(bào)，2008（1）：233?237.

[7]張?zhí)煜?，曹明華，蔡坤秀，等.不同甜椒品種抗寒性研究初報(bào)［J］.福建農(nóng)業(yè)學(xué)報(bào)，2014，29（2）：160?164.

[8] LI J ，YANG P ，KANG J ，et al. Transcriptome Analysis of Pepper （Capsicumannuum） RevealedaRoleof 24-EpibrassinolideinResponse to Chilling［J］. Frontiersin Plant Science ，2016，7：1281.

[9] SHINSY ， PARKMH ， CHOIJW ， etal. Genenetwork underlying the response of harvested pepper to chilling stress［J］. Journal of Plant Physiology，2017，219：112?122.

[10] KONG X M ，ZHOU Q ，LUO F ，et al. Transcriptome analysis ofharvested bell peppers（Capsicum annuum L.） in response to coldstress［J］. PlantPhysiologyandBiochemistry ，2019，139：314?324.

[11]孟雅寧，嚴(yán)立斌，張紅肖，等.低溫弱光脅迫對甜椒幼苗生長及生理指標(biāo)的影響［J］.辣椒雜志，2019，17（1）：22?25，30.

[12]劉凱歌，龔繁榮，宋云鵬，等.低溫弱光對甜椒幼苗生長和生理生化的影響及其與品種耐性的關(guān)系［J］.北方園藝，2020（3）：8?14.

[13] SALEEM M H ，PARVEEN A ，PERVEEN S ，et al. AllewTTdjFuykmC/kXgwXLk00yWQiQZRBCBvUvKYaoATHQs=viation ofcadmium toxicity in pea （Pisum sativum L.） through Zn? Lyssupplementationanditseffectsongrowthandantioxidantdefense［J］. EnvironmentalScienceandPollutionResearch，2024，31：10594?10608.

[14] HOULDSWORTH A. Role of oxidative stress in neurodegenerativedisorders： A review of reactive oxygen species and prevention byantioxidants［J］. BrainCommunications ，2024，6（1）： fcad356.

[15]莫億偉，郭振飛，謝江輝.溫度脅迫對柱花草葉綠素?zé)晒鈪?shù)和光合速率的影響［J］.草業(yè)學(xué)報(bào)，2011，20（1）：96?101.

[16]李德江，王薇，盧鳳丹，等.鹽脅迫對不同倍體水稻生理及抗氧化酶相關(guān)基因的影響［J/OL］.吉林農(nóng)業(yè)大學(xué)學(xué)報(bào)，1?13[2024-04-17]. http：//kns.cnki.net/kcms/detail/22.1100.S.20220424.1137.005. html.

[17] GILLSS， TUTEJAN. Reactiveoxygenspeciesandantioxidant machineryinabioticstresstoleranceincropplants［J］. Plant PhysiolBiochem，2010，48：909?930.

[18] WEI X ， SU X ， CAO P ， et al. Structure of spinach photosystemII–LHCIIsupercomplexat 3.2 ? resolution［J］. Nature ，2016，534（7605）：69?74.

[19] XUE T ，WAN H ，CHEN J ，et al. Genome-wide identification and expression analysis of the chlorophyll a/b binding protein gene family in oilseed（Brassica napus L.） under salt stress conditions［J］. Plant Stress，2024，11：100339.

[20] LIXW ， ZHUYL ， CHENCY ， etal.Cloningand characterizationof twochlorophyll A/B binding proteingenesand analysis of their genefamily inCamelliasinensis［J］. Scientific Reports，2020，10（1）：4602.

[21] WANG X ，LIU L. Crystal structure and catalytic mechanism of 7- hydroxymethyl chlorophyll a reductase［J］. Journal of Biological Chemistry，2016，291（25）：13349?13359.

[22] RENL ， SUNJ ， CHENS ， etal. Atranscriptomicanalysisof Chrysanthemum nankingense provides insights into the basis of low temperature tolerance［J］. BMC Genomics，2014，15（1）：844.

[23] GAO H ，ZHU L ，LIU T ，et al. Identification of a Novel OsCYP2 AllelethatWasInvolvedinRiceResponsetoLowTemperature Stress［J］. Phyton（0031?9457），2023，92（6）：1743.

[24] FRANZ S E B ， LIESE A ， KURTH J ， et al. Calcium-dependent proteinkinaseCPK21 functionsinabioticstressresponsein Arabidopsis thaliana［J］. Mol Plant，2011，4：83?96.

[25] RITONGA F N ，NGATIA J N ，WANG Y ， et al. AP2/ERF ， animportant cold stress-related transcription factor family in plants：Areview［J］. Physiologyand Molecular Biology of Plants ，2021，27（9）：1953?1968.

[26] GHORBANIM ， KAHRIZID ， CHAGHAKABOODIZ ， etal. EvaluationoftheexpressionofinducedgenesinresponsetodehydrationstressofCamelina （Camelinasativa）calli： dehydration stress genes in Camelina［J］. Cellular and MolecularBiology，2023，68（12）：1?4.

[27] LIS ， FUQ ， CHENL ， etal. ArabidopsisthalianaWRKY25， WRKY26， andWRKY33coordinateinductionofplantthermotolerance［J］. Planta，2011，233（6）：1237?1252.

[28] VANDERAUWERA S ，VANDENBROUCKE K ， INZé A ， et al. AtWRKY15 perturbation abolishes the mitochondrial stress responsethat steers osmotic stress tolerance in Arabidopsis［J］. ProceedingsoftheNationalAcademyofSciences ，2012，109（49）：20113?20118.

[29] KIDOKOROS ， YONEDAK ， TAKASAKIH ， etal. DifferentCold-SignalingPathwaysFunctionintheResponsestoRapidandGradualDecreasesinTemperature［J］. ThePlantCell ，2017，29（4）：760?774.

（責(zé)任編輯：柯文輝）