謝勇 , 王友紹, 張維仕

1. 熱帶海洋環(huán)境國家重點(diǎn)實(shí)驗(yàn)室(中國科學(xué)院南海海洋研究所), 廣東 廣州 510301;

2. 中國科學(xué)院大亞灣海洋生物綜合試驗(yàn)站, 廣東 深圳 518121;

3. 中國科學(xué)院南海生態(tài)環(huán)境工程創(chuàng)新研究院, 廣東 廣州 510301;

4. 中國科學(xué)院大學(xué), 北京 100049;

5. 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(廣州), 廣東 廣州 511458;

6. 臨沂市白沙埠中學(xué), 山東 臨沂 276035

紅樹林作為四大高生產(chǎn)力的海洋生態(tài)系統(tǒng)之一,具有“四高”特性, 在改善近海生態(tài)環(huán)境、凈化水質(zhì)、防浪護(hù)岸和保護(hù)生物多樣性等方面具有不可替代的重要作用 (王友紹, 2019; Wang et al, 2021)。近年來,隨著工農(nóng)業(yè)的快速發(fā)展, 不規(guī)范的海水養(yǎng)殖等導(dǎo)致沿海地區(qū)污染嚴(yán)重, 其中重金屬污染較為常見(王友紹 等, 2019)。紅樹林沉積物中重金屬濃度的增加在全世界都有記錄(Tam et al, 1997)。多項(xiàng)研究證實(shí),紅樹林已成為重金屬污染物的重要儲(chǔ)存庫(Usman et al, 2013; Bourgeois et al, 2019)。大量重金屬污染物的輸入給紅樹植物的生長發(fā)育帶來了嚴(yán)峻的挑戰(zhàn)。紅海欖是一種廣泛分布于我國南方沿海的紅樹科植物, 經(jīng)常作為紅樹林恢復(fù)種進(jìn)行人工種植, 具有較強(qiáng)的抗重金屬脅迫能力。木質(zhì)素合成反應(yīng)是紅海欖耐受重金屬脅迫的一個(gè)重要機(jī)制。植物細(xì)胞木質(zhì)素的合成有利于增加細(xì)胞壁厚度, 阻止金屬離子進(jìn)入細(xì)胞, 減小重金屬脅迫對植物的損害(芮海云 等,2019)。肉桂酸-4-羥基化酶是植物中分布最廣的細(xì)胞色素P450 單加氧酶家族成員, 是苯丙烷途徑中第二步反應(yīng)的關(guān)鍵酶(Fahrendorf et al, 1993), 催化反式肉桂酸生成 4-香豆酸, 對木質(zhì)素的合成有促進(jìn)作用。迄今為止, 許多陸地植物的C4H基因被克隆分析, 如荸薺C4H基因(宋慕波 等, 2020)、海南龍血樹C4H基因(梁惠楨 等, 2018)、紫蘇C4H基因(宋西紅 等, 2015)、茶樹C4H基因(姚勝波 等, 2015)等, 但關(guān)于紅樹植物C4H基因序列及其抗逆功能研究很少 (宋暉 等, 2012)。

本研究利用同源克隆和RACE 技術(shù), 從紅海欖中克隆得到了肉桂酸-4-羥基化酶基因的cDNA 全長,并運(yùn)用生物信息學(xué)方法對該基因進(jìn)行預(yù)測分析, 以紅海欖18S rRNA 作為內(nèi)參基因, 檢測紅海欖葉片中C4H基因表達(dá)情況與重金屬處理時(shí)間的關(guān)系。本研究為進(jìn)一步了解C4H基因在紅海欖響應(yīng)重金屬脅迫過程中的作用, 以及利用基因工程手段提高紅樹植物耐受重金屬脅迫的能力提供了基礎(chǔ)資料。

1 材料與方法

1.1 試驗(yàn)材料

植物材料: 實(shí)驗(yàn)用的紅海欖幼苗采自廣東省湛江市紅樹林育苗基地, 挑選長勢良好的幼苗移至26℃人工氣候培養(yǎng)箱進(jìn)行培養(yǎng), 用人工配置的混合重金屬水溶液(Cu2+112.5mg·L–1、Pb2+22.5mg·L–1、Cd2+4.5mg·L–1、Hg2+4.5mg·L–1)進(jìn)行脅迫處理。

試劑盒: RNAprep Pure Plant Kit(Polysaccharides$ polyphenolics-rich)多糖多酚植物總RNA 提取試劑盒(天根生化科技有限公司, Cat: DP121221);PrimeScript? RT reagent Kit with gDNA Eraser(Perfect Real Time)反轉(zhuǎn)錄試劑盒(TAKARA, Code No.RR047A); CISTROME(Cat,#E0104S) 混合酶;StarPrep Gel Extraction Kit 膠回收試劑盒(GenStar 公司 ); SMARTer? RACE 5′/3′ Kit (Clontech,No.634858)試劑盒。

1.2 實(shí)驗(yàn)方法

1.2.1 總RNA 提取與反轉(zhuǎn)錄

紅海欖葉片總RNA 的提取采用試劑盒的方法進(jìn)行, 提取完成后利用NanoDrop Lite 型號紫外分光光度計(jì)檢測 RNA 樣品的濃度和純度(OD260nm/OD280nm), 要求RNA 樣品濃度高于30ng·μL–1, OD 比值在1.8～2.1 之間, 紅海欖總RNA 經(jīng)檢測合格后進(jìn)行反轉(zhuǎn)錄反應(yīng)生成cDNA, 反轉(zhuǎn)錄反應(yīng)采用試劑盒方法, 操作方法見說明書。

1.2.2 中間片段克隆

從NCBI 網(wǎng)站下載薔薇超目植物的C4H 氨基酸序列, 將序列導(dǎo)入Base-by-Base 尋找保守區(qū)域, 將所得結(jié)果導(dǎo)入 j-CODEHOP 設(shè)計(jì)兼并引物C4H-F 和C4H-R 進(jìn)行中間片段克隆(表1), 以紅海欖cDNA 為模板, 使用CISTROME(Cat,#E0104S)混合酶進(jìn)行中間片段克隆擴(kuò)增, PCR 程序設(shè)置為95℃預(yù)變性3min,94 ℃ 30s、55℃ 30s、72 ℃ 1min, 35個(gè)循環(huán), 72℃延伸10min, 反應(yīng)結(jié)束后對PCR 產(chǎn)物進(jìn)行1%瓊脂糖凝膠電泳得到C4H基因中間片段, 利用試劑盒方法進(jìn)行膠回收后送往華大基因公司測序。

表1 引物列表Tab.1 Primer list

1.2.3C4H基因cDNA 全長擴(kuò)增

以紅海欖C4H基因中間片段為模板, 利用Promer5.0 設(shè)計(jì)RACE 反應(yīng)中的特異性引物。3′端反應(yīng)體系的配制與RACE 試劑盒說明書一致, 擴(kuò)增程序如下, 94 ℃ 30s、62 ℃ 30s、72 ℃ 3min, 28個(gè)循環(huán)。5′端按照RACE 試劑盒方法得到條帶不清晰,故以RACE 試劑盒中反轉(zhuǎn)錄生成的第一鏈cDNA 為模板, C4H-NGSP-5′和UPM 為引物進(jìn)行PCR 擴(kuò)增,程序如下95℃預(yù)變性3min, 94 ℃ 30s、53 ℃ 30s、72 ℃ 1min, 35個(gè)循環(huán), 72℃延伸10min, PCR 產(chǎn)物進(jìn)行1%瓊脂糖凝膠電泳后膠回收送往華大基因公司測序, 得到3′和5′端序列后進(jìn)行拼接得到C4H基因cDNA 全長。

1.2.4C4H基因生物信息學(xué)分析

生物信息學(xué)分析所利用的軟件、網(wǎng)址等詳見表2。

表2 生物信息學(xué)分析軟件或網(wǎng)址Tab. 2 Bioinformatics analysis software or website

1.2.5 基因表達(dá)差異分析

為研究重金屬處理下RsC4H基因的表達(dá)變化情況, 在重金屬處理后的第0h、3h、6h、12h、1d、2d、3d、4d、5d、6d、7d、10d 分別剪取3 片新鮮的、無機(jī)械損傷的嫩葉進(jìn)行qRT-PCR 反應(yīng), 本實(shí)驗(yàn)采用非特異性 SYBRGreenI 染料法在杭州博日FQD-96AQPCR 儀中進(jìn)行。qRT-PCR 引物詳見表1,以紅海欖 18S rRNA 作為內(nèi)參基因, 實(shí)驗(yàn)結(jié)果用2–ΔΔCT表示。利用SPSS21 進(jìn)行單因素方差分析, 分析不同處理時(shí)間下紅海欖葉片RsC4H基因的表達(dá)量與初始表達(dá)量之間的差異, 標(biāo)注*的數(shù)據(jù)代表與對照組具有顯著性差異(P＜0.05)。

2 結(jié)果與分析

2.1 RsC4H 基因的克隆

將3′和5′端序列進(jìn)行拼接后得到紅海欖C4H基因的cDNA 全長, 將其命名為RsC4H基因。RsC4H基因全長為2006bp, BLAST 比對結(jié)果顯示其與銀白楊C4H基因(Populus alba, XM_035065468.1)同源性高達(dá)83.49%, 與山葡萄C4H基因(Vitis amurensis,MH045992.1)同源性達(dá)84.46%, 與橡膠樹C4H基因(Hevea brasiliensis, XM_021825151.1) 同源性達(dá)82.69%, 證明克隆得到的基因全長確實(shí)是紅海欖的C4H基因序列。RsC4H基因從271 號核苷酸到1842號核苷酸為開放閱讀框區(qū)域, 開放閱讀框長1572bp,共編碼523 個(gè)氨基酸(圖1)。

圖1 RsC4H 基因的ORF 與編碼氨基酸序列Fig. 1 ORF and coding amino acid sequence of RsC4H gene

2.2 RsC4H 蛋白理化性質(zhì)分析

RsC4H 蛋白的理化性質(zhì)結(jié)果顯示, 紅海欖C4H蛋白相對分子量為60.18kD, 預(yù)測等電點(diǎn)為8.96。其中包含異亮氨酸(Leu)60 個(gè)占整個(gè)氨基酸序列的11.5%, 纈氨酸(Val)42 個(gè)占8.0%, 帶負(fù)電荷的殘基總數(shù) Asp + Glu 共63 個(gè), 帶正電荷的殘基總 Arg +Lys 共 71 個(gè), 原子總數(shù) 8567 個(gè), 分子式為C2748H4317N743O742S17, 脂肪族指數(shù)為99.87, 親水性的平均值(GRAVY)為–0.168(＜0), 是親水性蛋白質(zhì),不穩(wěn)定指數(shù)為44.09(>40), 屬于不穩(wěn)定蛋白。將紅海欖C4H 氨基酸序列輸入ProtScale 在線分析紅海欖C4H 蛋白的疏水-親水性, 結(jié)果顯示RsC4H 蛋白在N 端有一個(gè)明顯的疏水峰, 在C 端有兩個(gè)比較明顯的疏水峰。在第19 位氨基酸的疏水性最強(qiáng), 分值達(dá)到3.611, 在第416 位氨基酸的親水性最強(qiáng), 分值是–2.711, 整體上親水區(qū)更多, 親水性的平均值為–0.168, 為親水性蛋白。

2.3 RsC4H 蛋白二級結(jié)構(gòu)與三級結(jié)構(gòu)預(yù)測

利用SOPMA 在線軟件對RsC4H 蛋白的二級結(jié)構(gòu)進(jìn)行預(yù)測, 結(jié)果顯示RsC4H 蛋白包含47.42%的α 螺旋(Alpha-helix), 15.11%的β 折疊(extended strand), 4.78%的β 轉(zhuǎn)角(beta-turn)和32.70%的無規(guī)則卷曲(random coil), 紅海欖C4H 蛋白的二級結(jié)構(gòu)以α 螺旋和無規(guī)則卷曲為主。利用SWISS-MODEL在線軟件進(jìn)行三維模型構(gòu)建, 以高粱C4H(6vby.1.A)三維模型為模板, 氨基酸序列同源性達(dá)76.53%, 構(gòu)建出紅海欖的C4H 蛋白的三維模型(圖2)。結(jié)果顯示RsC4H 蛋白含有大量的α-螺旋,中間由無規(guī)則卷曲進(jìn)行連接, 這與SOPMA 二級結(jié)構(gòu)的預(yù)測結(jié)果保持一致。

圖2 紅海欖C4H 蛋白的三維模型(SWISS-MODEL)Fig. 2 Three-dimensional model of RsC4H protein(SWISS-MODEL)

2.4 RsC4H 蛋白跨膜螺旋與信號肽預(yù)測

利用TMHMM 2.0 在線工具對RsC4H 蛋白進(jìn)行跨膜螺旋預(yù)測, 結(jié)果顯示 RsC4H 蛋白共有兩個(gè)跨膜螺旋結(jié)構(gòu), 在肽鏈的N 端和C 端各含一個(gè), 1—3號氨基酸在膜結(jié)構(gòu)外側(cè), 4—23 號氨基酸為跨膜螺旋,24—503 號氨基酸在膜結(jié)構(gòu)內(nèi)側(cè), 504—522 號氨基酸為跨膜螺旋, 523 號氨基酸在膜結(jié)構(gòu)外側(cè)。利用SignalP 4.0 Server 對RsC4H 蛋白信號肽進(jìn)行分析,結(jié)果顯示其Y值、S值很小, 均小于0.5, 表明該蛋白不存在信號肽酶切位點(diǎn), 紅海欖C4H 蛋白沒有信號肽。

2.5 RsC4H 蛋白保守結(jié)構(gòu)域查詢與亞細(xì)胞定位

利用NCBI 網(wǎng)站Conserved Domain Search 功能尋找保守結(jié)構(gòu)域, 結(jié)果顯示紅海欖C4H 蛋白含一個(gè)保守結(jié)構(gòu)域PLN02394, 屬于細(xì)胞色素P450 超家族。對RsC4H 蛋白進(jìn)行亞細(xì)胞定位分析, 預(yù)測結(jié)果顯示紅海欖的C4H 蛋白是一種膜結(jié)合蛋白, 主要存在于膜結(jié)構(gòu)上或內(nèi)質(zhì)網(wǎng)上。

2.6 RsC4H 蛋白系統(tǒng)發(fā)育樹構(gòu)建

由于目前對紅樹植物的C4H基因研究很少, 除秋茄C4H基因的部分cDNA 序列被克隆出來以外,基本找不到紅樹植物的C4H 氨基酸序列或基因序列,所以將紅海欖的C4H 蛋白與從NCBI BLAST 結(jié)果中選取的麻風(fēng)樹(Jatropha curcas, XP_01207817 6.1)、決明(Senna tora, KAF7816240.1)、筍瓜(Cucurbita maxima, XP_023007159.1)、南瓜(Cucurbita moschata, XP_022947682.1)、雷公藤(Tripterygium wilfordii, XP_038693505.1)、蓮花(Nelumbo nucifera, XP_010253046.1)、芍藥(Paeonia lactiflora, AGG55322.1)、橄欖(Canarium album,ACR10242.1)、喜樹(Camptotheca cuminata, ANR763 95.1)、可可樹(Theobroma cacao, EOY20175.1)、地黃(Rehmannia glutinosa, QOJ53928.1)美國黑樹莓(Rubus occidentalis, ACM17896.1)等12 個(gè)植物的C4H 蛋白, 用 MEGA6 中 Phylogeny 模塊下的Neighber－joining 方法, Bootstrap 值設(shè)置為1000,構(gòu)建了系統(tǒng)發(fā)育樹, 分析不同物種的C4H 蛋白之間的進(jìn)化關(guān)系(圖3) , 系統(tǒng)發(fā)育樹結(jié)果顯示紅海欖和筍瓜(CuCurbita maxima, XP_023007159.1)、南瓜(CuCurbita moschata, XP_022947682.1)的C4H 蛋白聚在一個(gè)小分支, 表明它們的親緣關(guān)系更近, 而南瓜、筍瓜與紅海欖都屬于薔薇超目, 系統(tǒng)發(fā)育樹的結(jié)果與它們在植物分類系統(tǒng)上的地位一致, 表明系統(tǒng)發(fā)育樹具有可信度。

圖3 C4H 蛋白的系統(tǒng)發(fā)育樹Fig. 3 Phylogenetic tree of C4H protein

2.7 RsC4H 基因表達(dá)差異分析

利用qRT-PCR 技術(shù)檢測紅海欖幼苗在遭受重金屬脅迫10 天內(nèi)葉片中C4H基因的表達(dá)情況, 結(jié)果如圖4 所示。RsC4H基因在重金屬脅迫下除在第3h表達(dá)量顯著高于初始值外, 其他時(shí)間內(nèi)表達(dá)量與初始值無顯著差異(P>0.05)。RsC4H基因在重金屬脅迫下的3h 內(nèi)大幅上調(diào), 在第3h 時(shí)表達(dá)量達(dá)到最大值, 顯著高于初始表達(dá)量(P＜0.05), 達(dá)到處理前的4.29 倍, 隨后迅速回落, 從第6h 到第6 天RsC4H基因表達(dá)量均略低于初始表達(dá)水平, 小幅度波動(dòng), 從第7 天之后逐漸穩(wěn)定, 在第10 天時(shí)表達(dá)量為初始表達(dá)量的1.66 倍左右, 但是差異并不顯著。整體上,重金屬脅迫下紅海欖C4H基因的表達(dá)水平有所提高, 這有利于木質(zhì)素的合成, 增加細(xì)胞壁厚度, 進(jìn)而阻止金屬離子進(jìn)入細(xì)胞, 減弱金屬離子對植株的損害。

圖4 重金屬脅迫下紅海欖C4H 基因的表達(dá)情況Fig. 4 The expression of RsC4H gene under heavy metal stress

3 討論

苯丙烷代謝途徑可合成多種植物次生代謝產(chǎn)物,包括黃酮類化合物、木質(zhì)素、香豆素等, 是植物重要的次生代謝途徑, 在植物生長發(fā)育、機(jī)械支持、防御反應(yīng)等方面發(fā)揮重要作用(Winkel-Shirley,2001)。苯丙烷代謝途徑由苯丙氨酸在苯丙氨酸解氨酶的作用下生成反式肉桂酸, 反式肉桂酸在肉桂酸-4-羥化酶催化下生成香豆酸, 香豆酸在4-香豆酸輔酶A 連接酶的作用下產(chǎn)生香豆酰輔酶A, 然后進(jìn)入黃酮、木質(zhì)素等化合物的下游合成途徑。當(dāng)植物中C4H基因的表達(dá)受到抑制時(shí), 植物木質(zhì)化程度降低,總木質(zhì)素含量以及木質(zhì)素單體的含量降低, 據(jù)此認(rèn)為C4H基因的表達(dá)直接影響植物木質(zhì)素的合成(Blee et al, 2001; Sykes et al, 2015) 。而植物木質(zhì)素含量的增加可以增加細(xì)胞壁的厚度, 阻止重金屬離子進(jìn)入細(xì)胞, 細(xì)胞壁的阻擋作用是植物耐受重金屬脅迫的第一道屏障, 因此C4H基因的表達(dá)對于植物耐受重金屬脅迫有重要作用。

本研究利用同源克隆技術(shù)和RACE 方法擴(kuò)增得到了紅海欖C4H基因的cDNA 全長, 其開放閱讀框長1572bp, 共編碼523 個(gè)氨基酸。紅海欖C4H 蛋白是不穩(wěn)定的親水性蛋白, 有一個(gè)P450 功能結(jié)構(gòu)域,在N 端和C 端分別含有一個(gè)跨膜螺旋, 無信號肽, 亞細(xì)胞定位顯示其是一種膜結(jié)合蛋白, 主要在膜結(jié)構(gòu)或內(nèi)質(zhì)網(wǎng)上發(fā)揮功能, 二級結(jié)構(gòu)主要為α 螺旋和無規(guī)則卷曲。通過構(gòu)建系統(tǒng)發(fā)育樹發(fā)現(xiàn), 紅海欖和筍瓜、南瓜的C4H 蛋白聚在同一個(gè)小分支, 它們的C4H 蛋白親緣關(guān)系更近。系統(tǒng)發(fā)育樹中的標(biāo)尺為0.01, 指遺傳變異度為1%, 表明發(fā)育樹中各物種的C4H 蛋白變化度很小, C4H 蛋白保守性很強(qiáng), 預(yù)測由于C4H 蛋白是植物響應(yīng)逆境脅迫的重要酶, 所以具有比較高的遺傳穩(wěn)定性, 物種間的差異較小。

肉桂酸-4-羥基化酶基因可以受多種脅迫條件誘導(dǎo)表達(dá), 比如寒冷、水楊酸、脫落酸處理等都會(huì)誘導(dǎo)植物C4H基因的表達(dá)(Cheng et al, 2018), 海南龍血樹在脫落酸、油菜素內(nèi)酯、茉莉酸、紫外線等多種脅迫處理時(shí),DcC4H基因的表達(dá)量均會(huì)上調(diào)(梁惠楨 等, 2018), 而肉桂酸-4-羥基化酶可以介導(dǎo)木質(zhì)素的合成(Umemoto et al, 2016)。C4H基因在大麻葉中的高度表達(dá)導(dǎo)致了木質(zhì)素含量的增加(Docimo et al, 2013), 蓮霧果實(shí)C4H基因的表達(dá)水平和木質(zhì)素含量呈正相關(guān)關(guān)系, 且兩者的相關(guān)系數(shù)高達(dá)0.972 (黃利娜 等, 2020), 以上研究均證實(shí)了C4H基因的表達(dá)對于木質(zhì)素的合成有促進(jìn)作用。本研究中,紅海欖在遭受重金屬脅迫時(shí),C4H基因的表達(dá)量先快速提高, 在第3h 時(shí)表達(dá)量達(dá)最大值, 這表明該基因?qū)χ亟饘倜{迫的響應(yīng)十分迅速, 在受到重金屬脅迫后快速提高了C4H基因的表達(dá)量, 促進(jìn)木質(zhì)素的生成, 阻止金屬離子進(jìn)入細(xì)胞, 隨后表達(dá)量出現(xiàn)快速回調(diào), 但隨著脅迫時(shí)間的延長,C4H的表達(dá)量也開始逐步穩(wěn)定, 雖然與對照組無顯著差異, 但有高于初始表達(dá)水平的趨勢。在重金屬脅迫下紅海欖C4H基因表達(dá)量整體上表現(xiàn)出上調(diào)趨勢, 據(jù)此認(rèn)為RsC4H基因參與了紅海欖抗重金屬脅迫的響應(yīng)過程,促進(jìn)了木質(zhì)素的生成, 阻止金屬離子進(jìn)入細(xì)胞, 對于紅海欖耐受重金屬脅迫有積極作用。

4 結(jié)論

本研究克隆得到了紅海欖C4H基因的 cDNA全長, 其開放閱讀框長1572bp, 共編碼 523 個(gè)氨基酸, 編碼蛋白的相對分子量為60.18kD, 包含2 個(gè)跨膜螺旋結(jié)構(gòu)和一個(gè)P450 家族功能結(jié)構(gòu)域, 是一個(gè)不穩(wěn)定的親水性蛋白, 無信號肽, 二級結(jié)構(gòu)以α 螺旋和不規(guī)則卷曲為主, 亞細(xì)胞定位顯示主要在膜結(jié)構(gòu)或內(nèi)質(zhì)網(wǎng)上發(fā)揮作用。紅海欖C4H 蛋白與其他物種的C4H 蛋白同源性很高, 預(yù)測它們的功能也相似,包括催化反式肉桂酸生成香豆酸、促進(jìn)木質(zhì)素的合成及參與多種脅迫應(yīng)急反應(yīng)。在重金屬脅迫下, 紅海欖C4H基因可以快速響應(yīng)重金屬脅迫, 在紅海欖機(jī)體抗重金屬脅迫響應(yīng)的初始階段就開始發(fā)揮作用, 促進(jìn)木質(zhì)素合成, 增加細(xì)胞壁厚度, 阻止金屬離子進(jìn)入細(xì)胞對植株造成損傷。紅海欖C4H基因的成功克隆, 豐富了紅海欖抗性基因庫的內(nèi)容, 為進(jìn)一步探究紅海欖耐受重金屬脅迫的分子機(jī)制提供了基礎(chǔ)資料。