葫蘆科作物抗氧化基因的全基因組鑒定及分析

錢玉磊 閆晉強(qiáng) 劉文睿 吳智明 謝大森 江彪

摘要：抗氧化（oxidation resistance，OXR）基因是含有TLDc結(jié)構(gòu)域的一類基因，能有效清除活性氧，在抗氧化過程中發(fā)揮著重要作用。基于葫蘆科基因組數(shù)據(jù)庫，從6個(gè)葫蘆科作物中共鑒定出37個(gè)OXR成員，其中冬瓜、西瓜、甜瓜、葫蘆和黃瓜均有6個(gè)成員，南瓜有7個(gè)成員。通過系統(tǒng)發(fā)育以及保守結(jié)構(gòu)域分析，將這些OXR成員分為 5個(gè)亞族：OXR1、OXR2、OXR3、OXR4和OXR5。OXR基因在染色體上呈不均勻分布，同一亞族OXR成員所編碼的氨基酸數(shù)、分子量、外顯子數(shù)、內(nèi)含子數(shù)以及保守基序的排列模式都極為相似，表明OXR基因具有高度保守性。順式作用元件分析結(jié)果顯示，OXR基因的啟動(dòng)子區(qū)域含有多種順式作用元件以及MYB轉(zhuǎn)錄因子的結(jié)合位點(diǎn)?；驈?fù)制分析發(fā)現(xiàn)，葫蘆科OXR2成員通過片段復(fù)制或全基因組復(fù)制發(fā)生了一定擴(kuò)張，這些復(fù)制基因?qū)艿郊兓x擇。共線性分析結(jié)果顯示，在親緣關(guān)系較近的葫蘆科作物之間，OXR基因所在的共線性區(qū)塊較大。不同組織的表達(dá)譜分析結(jié)果顯示，葫蘆科OXR4成員有相似的表達(dá)模式，在葉片中表達(dá)量均較高；OXR2亞族內(nèi)復(fù)制基因出現(xiàn)了非對(duì)稱表達(dá)的現(xiàn)象。

關(guān)鍵詞：葫蘆科；OXR基因；保守結(jié)構(gòu)域；順式作用元件；生物信息學(xué)；表達(dá)譜；親緣關(guān)系

中圖分類號(hào)：S188；S642.01 文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1002-1302（2023）08-0051-10

基金項(xiàng)目：國家自然科學(xué)基金（編號(hào)：31972403、32102406）；廣東省農(nóng)業(yè)科學(xué)院中青年學(xué)科帶頭人項(xiàng)目（編號(hào)：R2020PY-JG003）；廣東省農(nóng)業(yè)科學(xué)院優(yōu)秀博士項(xiàng)目（編號(hào)：R2021YJ-YB2004）。

作者簡介：錢玉磊（1997—），男，河南鄭州人，碩士研究生，研究方向?yàn)槭卟朔肿由飳W(xué)。E-mail：qylei2022@163.com。

通信作者：江 彪，博士，研究員，研究方向?yàn)槭卟诉z傳育種與分子生物學(xué)。E-mail：jiangbiao@gdaas.cn。

植物受到脅迫刺激時(shí)，細(xì)胞內(nèi)的電子傳遞和氧化還原的動(dòng)態(tài)平衡被破壞，導(dǎo)致細(xì)胞內(nèi)的氧由氧化態(tài)轉(zhuǎn)變?yōu)檫€原態(tài)，從而產(chǎn)生過量的活性氧，造成蛋白質(zhì)、膜脂、DNA以及其他細(xì)胞組分嚴(yán)重?fù)p傷［1］。然而，活性氧也有對(duì)植物有利的一面，其可作為一種信號(hào)分子，協(xié)調(diào)細(xì)胞內(nèi)許多重要的信號(hào)傳導(dǎo)，如蛋白質(zhì)的修飾、基因表達(dá)和激素調(diào)節(jié)、協(xié)調(diào)適應(yīng)或傳遞防御反應(yīng)等［2-4］。植物體內(nèi)具有復(fù)雜的活性氧清除系統(tǒng)，主要包括保護(hù)酶系和抗氧化物質(zhì)，使體內(nèi)活性氧產(chǎn)生與降解處于動(dòng)態(tài)平衡，不對(duì)自身造成損傷［5］。

抗氧化（oxidation resistance，簡稱OXR）基因所編碼的蛋白質(zhì)含有保守結(jié)構(gòu)域TLDc［TBC （Tre2/Bub2/Cdc16），LysM （lysine motif），domain catalytic］，可以影響氧自由基清除基因的表達(dá)，從而降低細(xì)胞中活性氧水平［6］。近年來，在擬南芥基因組中已鑒定出6個(gè)具有TLDc保守結(jié)構(gòu)域的OXR基因家族成員，其中AtOXR2是一種線粒體蛋白，過表達(dá)AtOXR2可以提高擬南芥的光合作用效率，增強(qiáng)對(duì)氧化脅迫的耐受能力［7］。在蔬菜作物中，張夢云等基于十字花科基因組數(shù)據(jù)庫，對(duì)7種十字花科植物基因組中的抗氧化基因進(jìn)行了鑒定和生物信息學(xué)分析［8］。

葫蘆科是世界上最重要的植物科之一，包含許多常見的蔬菜和瓜果，具有重要的經(jīng)濟(jì)價(jià)值。隨著基因組測序技術(shù)的發(fā)展，黃瓜（Cucumis sativus）［9］、甜瓜（C. melo）［10］、西瓜（Citrullus lanatus）［11］、葫蘆（Lagenaria siceraria）［12］、南瓜（Cucurbita moschata）［13］、冬瓜（Benincasa hispida）［14］等葫蘆科作物先后完成了全基因組測序，為系統(tǒng)研究葫蘆科作物的抗氧化基因家族奠定基礎(chǔ)。本研究通過系統(tǒng)鑒定6個(gè)葫蘆科作物抗氧化基因家族，并開展生物信息學(xué)分析和不同組織的表達(dá)譜分析，為深入闡明葫蘆科作物OXR家族的基因功能奠定理論基礎(chǔ)。

1 材料與方法

1.1 葫蘆科OXR成員的鑒定

從葫蘆科基因組數(shù)據(jù)庫（http：//cucurbitgenomics.org/）下載黃瓜、甜瓜、西瓜、葫蘆、南瓜和冬瓜的基因組信息文件，Pfam網(wǎng)站（https：//pfam.xfam.org/）下載TLDc保守結(jié)構(gòu)域的隱馬爾可夫模型（PF07534），使用hmmsearch指令，搜索閾值<0.000 01，在葫蘆科6種作物的蛋白序列文件中進(jìn)行搜索。以已報(bào)道的6個(gè)AtOXR基因家族成員的蛋白質(zhì)序列作為查詢序列［7］，與葫蘆科作物的所有蛋白質(zhì)序列進(jìn)行BLASTP（檢索閾值<0.000 01）。將上述2種方法所鑒定到的候選葫蘆科OXR成員的蛋白序列提交至NCBI、Pfam和SMART數(shù)據(jù)庫，查詢結(jié)構(gòu)域特征，最后下載結(jié)果文件，并人為鑒定TLDc結(jié)構(gòu)域是否存在，最終確定葫蘆科OXR成員。

1.2 系統(tǒng)發(fā)育分析、保守結(jié)構(gòu)域分析和物理化學(xué)性質(zhì)分析

利用MEGA11中的ClustalW（默認(rèn)設(shè)置）比對(duì)鑒定出的37個(gè)葫蘆科OXR成員和6個(gè)AtOXR成員的蛋白質(zhì)序列，構(gòu)建Neighbor-Joining樹，使用 1 000 個(gè)重復(fù)進(jìn)行 bootstrap分析。根據(jù)擬南芥的分類方法將葫蘆科OXR成員分類［7］，利用Evolview v3將分類結(jié)果及每個(gè)OXR成員的蛋白序列和TLDc 結(jié)構(gòu)域所在位置可視化［15］，利用在線工具ProtParam（https：//web.expasy.org/protparam/）預(yù)測OXR成員的分子量（MW）和等電點(diǎn)（PI）。

1.3 葫蘆科OXR成員的染色體定位

根據(jù)相應(yīng)的基因組注釋文件，提取37個(gè)OXR成員的基因位置信息，然后使用MG2C（http：//mg2c.iask.in/mg2c_v2.1/）將其定位到相應(yīng)的染色體上。

1.4 保守基序和基因結(jié)構(gòu)分析

從基因組注釋文件中獲得OXR基因的外顯子和內(nèi)含子位置信息，將葫蘆科OXR成員和擬南芥OXR成員的蛋白質(zhì)序列提交MEME在線網(wǎng)站（https：//meme-suite.org/meme/doc/meme.html），最大基序（motif）檢索數(shù)設(shè)置為15，利用CFVisual（https：//github.com/ChenHuilong1223/CFVisual）將結(jié)果可視化。

1.5 啟動(dòng)子分析

提取每個(gè)OXR成員起始密碼子前2 500 bp序列，提交PlantCare網(wǎng)站（http：//bioinformatics.psb.ugent.be/webtools/plantcare/html/），獲得順式作用元件信息，然后以 R 4.1.3的tidyverse和ggplot2函數(shù)繪制啟動(dòng)子熱圖。

1.6 復(fù)制基因分析及選擇壓力分析

使用MCScanX［16］中的detect_collinearity_within_gene_families.pl分析6個(gè)葫蘆科物種內(nèi)的OXR復(fù)制基因?qū)σ约胺荗XR復(fù)制基因?qū)?，然后提取所有?fù)制基因?qū)υ诨蚪M中的位置信息，并利用Circos軟件［17］將復(fù)制基因?qū)梢暬?/p>

利用Muscle軟件分別比對(duì)復(fù)制基因?qū)Φ陌被嵝蛄校缓笥肞araAT［18］將氨基酸序列比對(duì)引導(dǎo)對(duì)應(yīng)的編碼序列比對(duì)，進(jìn)一步使用KaKs_Kaculator 2.0［19］計(jì)算復(fù)制基因?qū)Φ腒a（非同義替換率）/Ks（同義替換率），通過Ka/Ks的值來判斷基因在進(jìn)化過程中受到的選擇壓力。

1.7 共線性分析

利用葫蘆科作物的基因組文件及基因組注釋文件，使用JCVI［20］進(jìn)行6個(gè)葫蘆科物種間的OXR基因共線性分析。

1.8 葫蘆科OXR成員在不同器官組織的表達(dá)譜分析

從葫蘆科基因組數(shù)據(jù)庫下載冬瓜、南瓜（PRJNA385310）、甜瓜（PRJNA383830）、黃瓜（PRJNA80169）和葫蘆（PRJNA387615）的不同組織（根、莖、葉、花）的表達(dá)譜數(shù)據(jù)，分別以O(shè)XR1在根中的FPKM（fragments per kilobase of exon model per million mapped fragments）值為基準(zhǔn)，計(jì)算OXR成員在根、莖、葉、花的相對(duì)表達(dá)量，使用TBtools中的HeatMap［21］繪制熱圖，將結(jié)果可視化。

2 結(jié)果與分析

2.1 葫蘆科OXR基因家族成員的鑒定

根據(jù)TLDc的隱馬爾可夫模型（PF07534），使用hmmsearch指令在葫蘆科6種作物全基因組中共搜索到37個(gè)候選OXR成員基因。同時(shí)，利用同源BLAST方法鑒定到38個(gè)候選OXR成員基因。隨后，將2種方法鑒定到的候選成員的蛋白質(zhì)序列提交數(shù)據(jù)庫進(jìn)行確認(rèn)，結(jié)果表明，hmmsearch指令鑒定到的成員都含有TLDc保守結(jié)構(gòu)域，同源BLAST檢測出的38個(gè)成員包含hmmsearch指令鑒定到的37個(gè)成員，另外1個(gè)候選OXR成員不含有TLDc保守結(jié)構(gòu)域，因此將其去除。最終，在葫蘆科6種作物中共鑒定出37個(gè)OXR成員，其中冬瓜、黃瓜、西瓜、甜瓜和葫蘆均有6個(gè)成員，南瓜有7個(gè)成員。

2.2 系統(tǒng)發(fā)育、保守結(jié)構(gòu)域和物理化學(xué)性質(zhì)分析

以葫蘆科37個(gè)OXR成員和擬南芥6個(gè)OXR成員的蛋白質(zhì)序列構(gòu)建系統(tǒng)發(fā)育樹，可將葫蘆科OXR基因家族分為5個(gè)OXR亞家族，分別為OXR1亞族、OXR2亞族、OXR3亞族（AtOXR6歸入OXR3亞族）、OXR4亞族、OXR5亞族。葫蘆科6個(gè)作物在OXR1亞族、OXR3亞族、OXR4亞族和OXR5亞族均僅有1個(gè)成員。在OXR2亞族中，擬南芥僅有1個(gè)成員，而冬瓜、黃瓜、西瓜、甜瓜和葫蘆有2個(gè)成員，南瓜有3個(gè)成員。因此，相較于擬南芥，葫蘆科OXR2成員在進(jìn)化史中有一定的擴(kuò)張。

根據(jù)葫蘆科OXR成員所在亞家族，按照“物種拉丁名縮寫+亞家族+小寫字母”的形式命名。同一亞族OXR成員的TLDc結(jié)構(gòu)域的長度和位置相似，不同亞族的OXR成員的TLDc結(jié)構(gòu)域的位置和長度有差別（圖1），如OXR4亞族的TLDc結(jié)構(gòu)域的長度短于其他OXR族，且位于蛋白質(zhì)的中部，并不靠近C端。

理化性質(zhì)分析結(jié)果（表1）表明，葫蘆科作物OXR家族成員的氨基酸數(shù)量在232～1 059之間，等電點(diǎn)在4.65～7.07之間。分子量最小的是LsiOXR4，僅有25.64 ku，其開放閱讀框長度為 699 bp，編碼232個(gè)氨基酸；而分子量最大的CmaOXR4，有118.62 ku，其開放閱讀框長度達(dá)到 3 180 bp，編碼1 059個(gè)氨基酸。通過對(duì)比葫蘆科和擬南芥的OXR基因家族成員的氨基酸數(shù)發(fā)現(xiàn)，同一亞家族內(nèi)，除南瓜CmaOXR1和CmaOXR4外，其他葫蘆科OXR成員和擬南芥OXR成員的氨基酸數(shù)量相似。

2.3 染色體定位

將37個(gè)葫蘆科OXR基因家族成員定位到相應(yīng)的染色體上，結(jié)果顯示OXR基因在染色體上呈不均勻分布。南瓜OXR基因分布于6條染色體上，葫蘆OXR基因分布于5條染色體上，冬瓜和甜瓜分布于4條染色體上，黃瓜和西瓜OXR基因位于3條染色體上。不同染色體上分布的OXR基因數(shù)量不同，最多有4個(gè)（黃瓜3號(hào)染色體）（圖2）。

2.4 基因結(jié)構(gòu)和保守基序分析

為了研究葫蘆科作物OXR基因家族成員的進(jìn)化關(guān)系，對(duì)37個(gè)葫蘆科OXR成員和6個(gè)擬南芥OXR成員的外顯子和內(nèi)含子的數(shù)量及排列模式進(jìn)行了分析。結(jié)果表明，位于不同亞家族的所有OXR成員的外顯子和內(nèi)含子的數(shù)量以及排列模式差異顯著，而同一亞族的葫蘆科OXR成員的內(nèi)含子-外顯子排列模式基本一致（圖3）。

根據(jù)MEME網(wǎng)站分析，在葫蘆科OXR成員和擬南芥OXR成員中共檢測到15個(gè)保守基序，每個(gè)OXR成員大致含有5～9個(gè)不等的保守基序，同一個(gè)亞族內(nèi)的OXR成員所含的保守基序的數(shù)量和保守基序的分布模式大致相同，推測同亞族的OXR可能發(fā)揮相似的功能。不同亞族的OXR成員的保守基序的分布模式有所不同，可能是造成其功能差異的主要原因。

2.5 啟動(dòng)子分析

提取葫蘆科OXR基因起始密碼子前2 500 bp序列，提交至Plantcare網(wǎng)站，預(yù)測其順式作用元件（圖4）。共發(fā)現(xiàn)49種順式作用元件，主要分為4類：調(diào)控生長發(fā)育、光反應(yīng)、激素響應(yīng)和逆境脅迫相關(guān)。其中調(diào)控生長發(fā)育的元件包括對(duì)柵欄葉肉細(xì)胞分化、胚乳表達(dá)、晝夜節(jié)律、分生組織表達(dá)等調(diào)控元件，而在BhiOXR2a和CmeOXR2a中還發(fā)現(xiàn)有參與類黃酮生物合成的MYB結(jié)合位點(diǎn)。光反應(yīng)元件包含24個(gè)不同的元件，說明葫蘆科OXR基因家族可能廣泛參與了光調(diào)控途徑。多個(gè)葫蘆科OXR基因中有許多與脫落酸、赤霉素、水楊酸和生長素等植物激素響應(yīng)相關(guān)的順式作用元件。另外，OXR3亞族成員都具有參與脫落酸響應(yīng)的ABRE順式作用元件，數(shù)量1～4個(gè)不等，推測葫蘆科OXR3可能參與脫落酸的響應(yīng)。逆境脅迫相關(guān)元件主要包含多種非生物脅迫響應(yīng)元件（干旱、低溫、缺氧和創(chuàng)傷響應(yīng)元件）以及參與干旱誘導(dǎo)的MYB結(jié)合位點(diǎn)。

2.6 復(fù)制分析及基因選擇壓力分析

6個(gè)葫蘆科作物物種內(nèi)的共線性分析發(fā)現(xiàn)，在葫蘆科OXR基因家族中，只有OXR2亞族發(fā)生了復(fù)制（圖5）。其中，南瓜OXR2復(fù)制基因?qū)^多，而冬瓜、黃瓜、甜瓜、西瓜和葫蘆由于近期沒有發(fā)生全基因組復(fù)制事件，其復(fù)制基因?qū)^少。

通過葫蘆科OXR2復(fù)制基因?qū)Φ腒a/Ks分析，發(fā)現(xiàn)Ka/Ks的值均小于1，表明這些復(fù)制基因?qū)κ艿郊兓x擇（表2）。在同義突變率Ks分析中，CmaOXR2b-CmaOXR2c復(fù)制基因?qū)Φ腒s值為0.364 422，南瓜近期發(fā)生特異的全基因組復(fù)制事件（specific whole-genome duplication，sWGD）的Ks值為0.32［13］，二者相差不大。此外，CmaOXR2b-CmaOXR2c共線對(duì)所在的共線區(qū)塊所含的基因?qū)?shù)（183對(duì)）遠(yuǎn)多于CmaOXR2a-CmaOXR2c共線對(duì)所在的共線塊所含的基因?qū)?shù)（72對(duì)），證實(shí)CmaOXR2b-CmaOXR2c是由南瓜的近期發(fā)生sWGD而產(chǎn)生的。

2.7 共線性分析

為研究OXR基因家族在葫蘆科內(nèi)的進(jìn)化歷史，根據(jù)之前報(bào)道的物種分化樹將葫蘆科作物兩兩開展共線性分析。結(jié)果發(fā)現(xiàn) 親緣關(guān)系越近 OXR所在的共線性區(qū)塊越大，如黃瓜和甜瓜、西瓜和葫蘆。ClaOXR1和ClaOXR3位于西瓜的10號(hào)染色體上，LsiOXR1和LsiOXR3位于葫蘆的3號(hào)染色體上，都在同一個(gè)較大的共線塊（OXR1+OXR3 block）中。此外，葫蘆科不同作物同一亞族成員間存在一定的共線性。不同作物OXR2a之間、不同作物OXR2b之間存在共線性，不同作物OXR2a和OXR2b間不存在共線性，說明在葫蘆科祖先物種中可能發(fā)生了古老的復(fù)制事件，產(chǎn)生了OXR2a和OXR2b（圖6）。

2.8 葫蘆科OXR基因家族的組織表達(dá)譜分析

基于葫蘆科基因組數(shù)據(jù)庫中已發(fā)表的表達(dá)譜數(shù)據(jù)，分析OXR基因在不同組織中的表達(dá)模式。結(jié)果（圖7）顯示，不同亞族基因的表達(dá)模式有所差異，如冬瓜BhiOXR3在根、莖、葉、花中表達(dá)量無明顯差異，而BhiOXR4在葉片中表達(dá)量明顯高于其他組織。此外，不同作物相同亞族基因表達(dá)也存在一定的差異，冬瓜、葫蘆、甜瓜、黃瓜OXR4基因在葉片中表達(dá)量明顯高于其他組織，而南瓜OXR4基因則無明顯的葉片高豐度表達(dá)。值得注意的是，復(fù)制基因?qū)Γ˙hiOXR2a-BhiOXR2b、LsiOXR2a-LsiOXR2b、CsaOXR2a-CsaOXR2b、CmaOXR2b-CmaOXR2c）在不同組織中的表達(dá)存在較大差異，出現(xiàn)表達(dá)分化現(xiàn)象。

3 討論與結(jié)論

抗氧化（OXR）基因有助于保護(hù)真核生物免受活性氧的傷害［22］。目前，OXR基因家族已在擬南芥、小麥、向日葵和某些十字花科植物的基因組中被鑒定出來，并深入研究了部分基因的功能［7-8，23-25］。

然而，關(guān)于葫蘆科OXR家族基因的研究還未見報(bào)道。本研究基于葫蘆科基因組數(shù)據(jù)庫，從6個(gè)葫蘆科作物中共鑒定出37個(gè)OXR基因，其中冬瓜、黃瓜、西瓜、甜瓜和葫蘆均有6個(gè)OXR基因，而南瓜有7個(gè)OXR基因。

根據(jù)氨基酸序列和保守結(jié)構(gòu)域的相似性，葫蘆科OXR基因家族可分為5個(gè)亞家族，這一結(jié)果同擬南芥、向日葵和十字花科植物的分類結(jié)果相似［7-8，25］。同一亞家族的擬南芥 OXR基因和葫蘆科OXR基因編碼的氨基酸個(gè)數(shù)、開放閱讀框和分子量的大小都相似，其內(nèi)含子-外顯子的數(shù)量及保守基序相似，說明同一亞家族內(nèi)的OXR基因可能具有相似的生物學(xué)功能。

OXR基因的啟動(dòng)子含有多種順式作用元件，表明OXR基因可能參與多種植物代謝調(diào)控途徑。過表達(dá)擬南芥OXR基因，可通過影響水楊酸代謝而清除多余的活性氧，從而能夠提高對(duì)氧化脅迫的耐受性［7，24］。葫蘆科OXR基因的啟動(dòng)子區(qū)域包含調(diào)控生長發(fā)育、光反應(yīng)、激素響應(yīng)和非生物脅迫響應(yīng)相關(guān)元件，這一結(jié)果與之前報(bào)道的十字花科植物OXR基因和小麥OXR基因的順式作用元件［8，23］相似，說明葫蘆科OXR基因可能在抵御生物和非生物脅迫時(shí)發(fā)揮重要作用。

基因復(fù)制是植物基因組進(jìn)化的主要驅(qū)動(dòng)力，為遺傳革新、表型多樣性以及物種分化提供了大量的原始遺傳物質(zhì)［26-27］。葫蘆科OXR2亞族成員通過片段復(fù)制或全基因組復(fù)制發(fā)生了擴(kuò)張，并且OXR2復(fù)制基因?qū)Φ腒a/Ks 比值均小于 1。Ka/Ks<1 表示基因受到純化選擇［28］，說明葫蘆科OXR2亞族在進(jìn)化上相對(duì)保守，有助于維持OXR2亞族成員功能的穩(wěn)定性。通過探尋OXR基因家族在葫蘆科內(nèi)的進(jìn)化歷史，發(fā)現(xiàn) 6個(gè)葫蘆科作物的OXR2a和OXR2b相互之間都具有共線性關(guān)系，推測OXR2亞族的成員在這6個(gè)葫蘆科作物的祖先物種中通過古老的復(fù)制事件，產(chǎn)生了OXR2a和OXR2b等2個(gè)同源基因?qū)?，伴隨著葫蘆科物種的分化，OXR2a和OXR2b也被保留在了6個(gè)葫蘆科物種的基因組中。南瓜基因組中有3個(gè)OXR2基因，其中CmaOXR2b和CmaOXR2c是同源基因，這一發(fā)現(xiàn)進(jìn)一步證實(shí)南瓜近期發(fā)生了sWGD事件［13］。

為了研究OXR基因的表達(dá)特征，本研究從葫蘆科基因組數(shù)據(jù)庫下載冬瓜、南瓜、甜瓜、黃瓜和葫蘆不同組織的表達(dá)譜數(shù)據(jù)。表達(dá)譜分析發(fā)現(xiàn)不同亞族OXR基因的表達(dá)量有差異，表明它們在不同組織器官中的抗氧化作用機(jī)制不盡相同。葫蘆科OXR4亞族的基因的表達(dá)模式較為相似，暗示OXR4亞族的基因可能具有相似的抗氧化功能。喬鑫通過對(duì)梨等141種植物基因組中的不同類型的復(fù)制基因的鑒定和表達(dá)分析，發(fā)現(xiàn)復(fù)制基因廣泛存在非對(duì)稱表達(dá)現(xiàn)象（2個(gè)同源基因拷貝的表達(dá)水平己經(jīng)偏離了1 ∶1的對(duì)稱性表達(dá)模式）［27］，然而葫蘆科OXR2亞族內(nèi)復(fù)制基因存在顯著的表達(dá)差異，也出現(xiàn)了非對(duì)稱表達(dá)的現(xiàn)象。這可能是由于基因復(fù)制發(fā)生后，2個(gè)同源基因拷貝在啟動(dòng)子區(qū)域的分化導(dǎo)致它們之間的表達(dá)分化［29］，但其表達(dá)產(chǎn)生分化的機(jī)制還需進(jìn)一步解析。

綜上所述，本研究利用多種生物信息學(xué)分析網(wǎng)站和軟件，系統(tǒng)鑒定了6種葫蘆科作物的OXR基因，并綜合分析了葫蘆科OXR基因的系統(tǒng)發(fā)育關(guān)系、物理化學(xué)性質(zhì)、染色體定位、基因結(jié)構(gòu)、保守基序、順式作用元件、在葫蘆科作物中的進(jìn)化歷史和在不同組織中表達(dá)模式等。鑒定和分析結(jié)果為進(jìn)一步研究葫蘆科作物OXR基因家族在響應(yīng)非生物脅迫方面的生物學(xué)功能奠定基礎(chǔ)。

參考文獻(xiàn)：

［1］李 格，孟小慶，蔡 敬，等. 活性氧在植物非生物脅迫響應(yīng)中功能的研究進(jìn)展［J］. 植物生理學(xué)報(bào)，2018，54（6）：951-959.

［2］Czarnocka W，Karpiński S.Friend or foe？Reactive oxygen species production，scavenging and signaling in plant response to environmental stresses［J］. Free Radical Biology and Medicine，2018，122：4-20.

［3］Foyer C H，Ruban A V，Graham N. Viewing oxidative stress through the lens of oxidative signalling rather than damage［J］. The Biochemical Journal，2017，474（6）：877-883.

［4］Mittler R. ROS are good［J］. Trends in Plant Science，2017，22（1）：11-19.

［5］趙 瀅，王振興，許培磊，等. 山葡萄雙豐和左優(yōu)紅葉綠素?zé)晒馓匦约盎钚匝醮x與低溫傷害的關(guān)系［J］. 園藝學(xué)報(bào)，2018，45（4）：650-658.

［6］Finelli M J，Oliver P L.TLDc proteins：new players in the oxidative stress response and neurological disease［J］. Mammalian Genome，2017，28（9）：395-406.

［7］Colombatti F，Mencia R，Garcia L，et al. The mitochondrial oxidation resistance protein AtOXR2 increases plant biomass and tolerance to oxidative stress［J］. Journal of Experimental Botany，2019，70（12）：3177-3195.

［8］張夢云，袁凌云，朱世東，等. 十字花科植物抗氧化基因（OXR）的鑒定與生物信息學(xué)分析［J/OL］. 分子植物育種，2021：1-11.（2021-08-26）［2022-12-20］. http：//kns.cnki.net/kcms/detail/46.1068.S.20210825.1923.021.html.

［9］Huang S W，Li R Q，Zhang Z H，et al. The genome of the cucumber，Cucumis sativus L.［J］. Nature Genetics，2009，41（12）：1275-1281.

［10］Garcia-Mas J，Benjak A，Sanseverino W，et al. The genome of melon （Cucumis melo L.）［J］. Proceedings of the National Academy of Sciences of the United States of America，2012，109（29）：11872-11877.

［11］Guo S G，Zhang J G，Sun H H，et al. The draft genome of watermelon （Citrullus lanatus） and resequencing of 20 diverse accessions［J］. Nature Genetics，2013，45（1）：51-58.

［12］Wu S，Shamimuzzaman M，Sun H，et al. The bottle gourd genome provides insights into Cucurbitaceae evolution and facilitates mapping of a Papaya ring-spot virus resistance locus［J］. The Plant Journal，2017，92（5）：963-975.

［13］Sun H H，Wu S，Zhang G Y，et al. Karyotype stability and unbiased fractionation in the paleo-allotetraploid Cucurbita genomes［J］. Molecular Plant，2017，10（10）：1293-1306.

［14］Xie D S，Xu Y C，Wang J P，et al. The wax gourd genomes offer insights into the genetic diversity and ancestral cucurbit karyotype［J］. Nature Communications，2019，10：5158.

［15］Subramanian B，Gao S H，Lercher M J，et al. Evolview v3：a webserver for visualization，annotation，and management of phylogenetic trees［J］. Nucleic Acids Research，2019，47（1）：270-275.

［16］Wang Y P，Tang H B，de Barry J D，et al. MCScanX：a toolkit for detection and evolutionary analysis of gene synteny and collinearity［J］. Nucleic Acids Research，2012，40（7）：49.

［17］Krzywinski M，Schein J，Birol I，et al. Circos：an information aesthetic for comparative genomics［J］. Genome Research，2009，19（9）：1639-1645.[HJ2mm]

［18］Zhang Z，Xiao J F，Wu J Y，et al. ParaAT：a parallel tool for constructing multiple protein-coding DNA alignments［J］. Biochemical and Biophysical Research Communications，2012，419（4）：779-781.

［19］Wang D，Zhang Y，Zhang Z et al. KaKs_Calculator 2.0：a toolkit incorporating gamma-series methods and sliding window strategies［J］. Genomics，Proteomics & Bioinformatics，2010，8（1）：77-80.

［20］Tang H，Bowers J E，Wang X，et al. Synteny and collinearity in plant genomes［J］. Science，2008，320（5875）：486-488.

［21］Chen C J，Chen H，Zhang Y，et al. TBtools：an integrative toolkit developed for interactive analyses of big biological data［J］. Molecular Plant，2020，13（8）：1194-1202.

［22］Blaise M，Alsarraf H M A B，Wong Jaslyn E M M，et al. Crystal structure of the TLDc domain of oxidation resistance protein 2 from zebrafish［J］. Proteins，2012，80（6）：1694-1698.

［23］李 寒，夏鵬亮，唐艷紅，等. 小麥OXR基因家族全基因組生物信息學(xué)及表達(dá)分析［J/OL］. 分子植物育種，2022：1-10.（2022-04-11）［2022-05-20］. http：//kns.cnki.net/kcms/detail/46.1068.S.20220411.0854.004.html.

［24］Mencia R，Céccoli G，F(xiàn)abro G，et al. OXR2 increases plant defense against a hemibiotrophic pathogen via the salicylic acid pathway［J］. Plant Physiology，2020，184（2）：1112-1127.

［25］Torti P，Raineri J，Mencia R，et al. The sunflower TLDc-containing protein HaOXR2 confers tolerance to oxidative stress and waterlogging when expressed in maize plants［J］. Plant Science，2020，300：110626.

［26］Flagel L E，Wendel J F. Gene duplication and evolutionary novelty in plants［J］. The New Phytologist，2009，183（3）：557-564.

［27］喬 鑫. 梨等植物基因組中重復(fù)基因的鑒定與進(jìn)化分析 ［D］. 南京：南京農(nóng)業(yè)大學(xué)，2018：72-76.

［28］周至銘，楊佳寶，張 程，等. 向日葵LACS家族鑒定及響應(yīng)非生物脅迫表達(dá)分析［J］. 園藝學(xué)報(bào)，2022，49（2）：352-364.

［29］Hahn M W. Distinguishing among evolutionary models for the maintenance of gene duplicates［J］. Journal of Heredity，2009，100（5）：605-617.