王建婷 劉廣達 吳歸 安天琪 李娜

摘要 使用轉(zhuǎn)錄組學的方法研究水楊酸和茉莉酸甲酯處理微型月季的基因表達，發(fā)現(xiàn)2種處理的植物-病原菌相互關(guān)系通路均有數(shù)量較多的差異表達基因，并對該通路中差異表達基因連鎖的微衛(wèi)星序列進行了統(tǒng)計。該研究有助于進一步理解水楊酸和茉莉酸甲酯提升微型月季抗逆性的機理。篩選與抗逆基因連鎖的微衛(wèi)星序列可為今后的育種研究提供幫助。

關(guān)鍵詞 微型月季;微衛(wèi)星序列;水楊酸;茉莉酸甲酯

中圖分類號 S 943.2? 文獻標識碼 A? 文章編號 0517-6611（2021）21-0124-04

doi：10.3969/j.issn.0517-6611.2021.21.030

開放科學（資源服務(wù)）標識碼（OSID）：

Expression of Stress Resistance Genes and Linkage Analysis of Microsatellite Sequences in Miniature Rose

WANG Jian-ting， LIU Guang-da， WU Gui et al

（ School of Pharmacy， Baotou Medical College， Inner Mongolia University of Science Technology， Baotou， Inner Mongolia 014040）

Abstract In this study， we used transcriptomics methods to study the gene expression of miniature rose treated with salicylic acid and methyl jasmonate， it was found that there were a number of differentially expressed genes in the plant-pathogen interaction pathway after the two treatments， and the microsatellite sequences linked to the differentially expressed genes in the pathway were statistically analyzed. This study will help to further understand the mechanism of salicylic acid and methyl jasmonate improving the stress resistance of miniature rose. Screening microsatellite sequences linked to stress resistance genes will provide help for future breeding research.

Key words Miniature rose;Microsatellite sequences;Salicylic acid;Methyl jasmonate

基金項目? 自治區(qū)級大學生創(chuàng)新創(chuàng)業(yè)訓練計劃項目（201910127017）。

作者簡介 王建婷（1998—），女，四川宜賓人，從事藥用植物研究。*通信作者，講師，博士，從事藥用植物生物技術(shù)研究。

收稿日期 2021-05-27

微型月季（Rosa hybrid minima）為薔薇科植物，由于具有花色奇異、全年開放和適于盆栽等特點[1]，深受消費者喜愛，具有較高的市場價值。生物和非生物脅迫是微型月季生產(chǎn)中的重要問題。在栽培過程中，微型月季會受到多種病原微生物的生物脅迫[2]。在運輸過程中，微型月季易受到黑暗、缺水和震動等非生物脅迫[3]。提升微型月季的抗逆性對其商品化生產(chǎn)有著重要作用。水楊酸和茉莉酸甲酯是植物體內(nèi)普遍存在的抗性信號分子，可以對生物和非生物脅迫產(chǎn)生響應(yīng)，提升植物的抗性[4-6]。外施水楊酸和茉莉酸甲酯可以促進多種植物抗逆性[7]。

微衛(wèi)星序列又稱簡單重復序列（simple sequence repeats，SSR），是存在于大部分真核生物基因組中的短串重復序列[8]，SSR分子標記技術(shù)在植物育種和近緣植物親緣關(guān)系研究方面有著廣泛的應(yīng)用[9]。隨著轉(zhuǎn)錄組技術(shù)的發(fā)展，許多植物開展了針對轉(zhuǎn)錄組中SSR的研究[10-11]。相比于基因組中的SSR，轉(zhuǎn)錄組中的SSR存在于發(fā)生轉(zhuǎn)錄的基因序列中，這些SSR可能與功能基因直接相關(guān)，從而與性狀有更好的相關(guān)性[12-13]。

該研究以水楊酸和茉莉酸甲酯處理微型月季，利用轉(zhuǎn)錄組學的方法研究處理后差異表達基因富集程度較高的代謝通路，分析與抗逆性相關(guān)的代謝通路，探討外施水楊酸和茉莉酸甲酯對微型月季抗逆性作用機制的差異，尋找與主要抗逆性代謝通路中基因相連鎖的SSR序列，為提升微型月季抗逆性的育種工作提供科學依據(jù)。

1 材料與方法

1.1 材料

選擇盆栽微型月季（Rosa hybrid minima），采取葉面噴施的方法，水楊酸處理噴施100 mL 1 mmol/L水楊酸，茉莉酸甲酯處理噴施100 mL 1 mmol/L茉莉酸甲酯，對照噴施100 mL蒸餾水（CK）。處理24 h后，采集微型月季的花瓣。

1.2 方法

采集花瓣后液氮速凍，使用Invitrogen Trizol試劑提取總RNA，經(jīng)質(zhì)檢后送聯(lián)川生物進行文庫的構(gòu)建和高通量測序，測序采用Illumina Hiseq 4000平臺，測序讀長為雙端2×150 bp。

1.3 數(shù)據(jù)處理

測序產(chǎn)生的原始數(shù)據(jù)（raw data）去除接頭序列，過濾掉不合格序列后得到有效數(shù)據(jù)（clean data）。采取混合組裝的策略，歸一化得到Unigene，使用DIAMOND軟件[14]對組裝的Unigenes進行注釋。使用Salmon軟件[15]計算Unigenes的表達水平，使用R軟件包edge R篩選差異表達Unigenes[16]，使用perl腳本進行KEGG富集分析，設(shè)置差異表達的標準為│log2（fold change）│>1（P≤0.05）。富集因子（rich factor）為位于該KEGG通路的差異基因數(shù)/位于該KEGG通路的總基因數(shù)。對顯著差異的基因占某個通路所有基因數(shù)的比例和所有通路的總基因數(shù)目中的差異基因比例進行檢驗，得出富集分析中通路的P值。SSR的鑒定使用MISA軟件（http：∥pgrc.ikp-gatersleben.de/misa/）對Unigenes進行分析。

2 結(jié)果與分析

2.1 水楊酸對微型月季基因轉(zhuǎn)錄的影響

水楊酸處理組與CK相比，在差異表達基因富集程度最高的20條KEGG代謝通路中，異黃酮生物合成通路（isoflavonoid biosynthesis）的富集因子值最大，為0.3 差異表達Unigene基因數(shù)為8，通路的P值為3.9×10-4。植物-病原菌相互關(guān)系（plant-pathogen interaction）通路差異表達Unigene基因數(shù)最多，為13 富集因子為0.10，通路的P值為1.4×10-3（圖1）。

2.2 茉莉酸甲酯對微型月季基因轉(zhuǎn)錄的影響

茉莉酸甲酯處理組與對照相比，在差異表達基因富集程度最高的20條代謝通路中，異黃酮生物合成（isoflavonoid biosynthesis）通路的富集因子值最大，為0.3 差異表達Unigene基因數(shù)為8，通路的P值為2.2×10-6。植物-病原菌相互關(guān)系（plant-pathogen interaction）通路差異表達Unigene基因的數(shù)量最多，為8 富集因子為0.06，通路的P值為9.2×10-6（圖2）。

2.3 抗逆性相關(guān)代謝通路基因表達分析

將差異表達的Unigene在KEGG數(shù)據(jù)庫中進行注釋，水楊酸處理注釋到27個KO條目中，注釋到108個基因（表1）。茉莉酸甲酯處理注釋到24個KO條目中，注釋到64個基因（表2），其中有3個KO條目未注釋到基因。

2.4 植物-病原菌相互關(guān)系通路差異表達基因與微衛(wèi)星序列連鎖分析

與CK相比，水楊酸處理和茉莉酸甲酯處理后植物-病原菌相互關(guān)系通路中顯著差異表達的基因數(shù)均為最多，且與抗逆性有著緊密的聯(lián)系，有利于發(fā)掘抗逆性相關(guān)的SSR。設(shè)置連鎖的條件為KEGG注釋基因與SSR存在于同一Unigene中，進行統(tǒng)計并記錄基因名與具體SSR序列（表3）。水楊酸處理與CK相比，植物-病原菌相互關(guān)系通路中有134個Unigene數(shù)量差異顯著，注釋到108個基因，其中有26個基因與SSR序列連鎖。茉莉酸甲酯處理與CK相比，植物-病原菌相互關(guān)系通路中有82個Unigene數(shù)量差異顯著，注釋到64個基因，其中有15個基因與SSR序列連鎖。WRKY24、NQR、DREB3和PBL5 4個基因在水楊酸和茉莉酸甲酯處理下均與SSR存在連鎖關(guān)系。WRKY65、MYB108、PDPK2、WNK5、At2g26730、CRCK2、WRKY53、WRKY48和WRKY24基因與2個SSR連鎖。WRKY22基因與3個SSR連鎖。

3 討論

在植物應(yīng)對生物脅迫時，植物體內(nèi)水楊酸、茉莉酸和乙烯通過影響各種信號通路相互作用提升抗性;而在應(yīng)對非生物脅迫時，脫落酸調(diào)節(jié)的信號通路起到重要作用。在提升植物抗逆性過程中，這2類信號通路間有著復雜的相互影響，整體上處于拮抗關(guān)系[17]。異黃酮類物質(zhì)有助于提升植物的生物抗性，如大豆異黃酮可以抵抗植物病原微生物和食草性昆蟲，誘導根瘤菌結(jié)瘤以利于植物抵御逆境[18]。在該試驗中，水楊酸和茉莉酸甲酯處理后，異黃酮生物合成通路的富集因子值均為最大，異黃酮生物合成通路中的基因差異表達可能有助于特定種類異黃酮的合成，從而提升植物的抗性。水楊酸處理后，與CK相比，植物-病原菌相互關(guān)系、淀粉和蔗糖代謝和植物激素信號轉(zhuǎn)導通路差異表達基因數(shù)量較多。而茉莉酸甲酯處理后，僅植物-病原菌相互關(guān)系通路差異表達基因的數(shù)量較多，且2種處理后植物-病原菌相互關(guān)系通路的P值存在差異，說明2種誘導子的作用機制存在差異。

水楊酸和茉莉酸甲酯均可以通過植物-病原菌相互關(guān)系通路提升植物抗性，尤其是提升抵抗病原菌的能力[19]，也可以產(chǎn)生系統(tǒng)性獲得抗性（systemic acquired resistance，SAR）。

SAR可以長效提升植物抗性和緩解植物受到的非生物脅迫[ 17，20]。在該研究中，與CK相比，水楊酸處理組和茉莉酸甲酯處理組中植物-病原菌相互關(guān)系通路中顯著差異表達的基因數(shù)均為最多，說明該通路對微型月季的抗逆性有著重要作用。在差異表達基因與微衛(wèi)星序列連鎖分析中，有一定數(shù)量的SSR與植物-病原菌相互關(guān)系通路中的差異表達基因相連鎖，為開發(fā)抗逆性相關(guān)代謝通路的分子標記提供了豐富的信息位點。

參考文獻

[1] 馮歡，易姝利，謝佳恒，等.微型月季愈傷組織誘導及植株再生[J].植物學報，201 49（5）：595-602.

[2] 竇坦德，劉和風.出口盆栽微型月季的病蟲害控制及檢疫技術(shù)研究[J].檢驗檢疫科學，2006，16（S1）：55-58.

[3] MLLER R，SISLER E C，SEREK M.Stress induced ethylene production，ethylene binding，and the response to the ethylene action inhibitor 1-MCP in miniature roses[J].Scientia horticulturae，2000，83（1）：51-59.

[4] HALIM V A，VESS A，SCHEEL D，et al.The role of salicylic acid and jasmonic acid in pathogen defence[J].Plant biology，2006，8（3）：307-313.

[5] HORVTH E，SZALAI G，JANDA T.Induction of abiotic stress tolerance by salicylic acid signaling[J].Journal of plant growth regulation，2007，26（3）：290-300.

[6] 嚴俊鑫，遲德富，張永強，等.JA、SA對重瓣玫瑰生長發(fā)育和防御酶活性的影響[J].北京林業(yè)大學學報，201 35（3）：128-136.

[7] 包穎，魏琳燕，陳超.水楊酸和茉莉酸甲酯對鹽脅迫下月季品種月月粉生理特性的影響[J].云南農(nóng)業(yè)大學學報（自然科學），2020，35（6）：1040-1045.

[8] POWELL W，MACHRAY G C，PROVAN J.Polymorphism revealed by simple sequence repeats[J].Trends in plant science，1996，1（7）：215-222.

[9] 楊夢婷，黃洲，干建平，等.SSR分子標記的研究進展[J].杭州師范大學學報（自然科學版），2019，18（4）：429-436.

[10] 杜曉華，楊雅萍，朱小佩，等.三色堇轉(zhuǎn)錄組SSR分析及分子標記開發(fā)[J].園藝學報，2019，46（4）：797-806.

[11] 王榮香，徐子健，于平，等.海南粗榧轉(zhuǎn)錄組SSR特征分析[J].分子植物育種，2019，17（6）：1951-1957.

[12] EUJAYL I，SORRELLS M E，BAUM M，et al.Isolation of EST-derived microsatellite markers for genotyping the A and B genomes of wheat[J].Theoretical and applied genetics，200 104（2/3）：399-407.

[13] WANG H Y，WEI Y M，YAN Z H，et al.EST-SSR DNA polymorphism in durum wheat（Triticum durum L.）collections[J].Journal of applied genetics，2007，48（1）：35-42.

[14] BUCHFINK B，XIE C，HUSON D H.Fast and sensitive protein alignment using DIAMOND[J].Nature methods，2015，12（1）：59-60.

[15] PATRO R，DUGGAL G，LOVE M I，et al.Salmon provides fast and bias-aware quantification of transcript expression[J].Nature methods，2017，14（4）：417-419.

[16] ROBINSON M D，MCCARTHY D J，SMYTH G K.edgeR：a Bioconductor package for differential expression analysis of digital gene expression data[J].Bioinformatics，2010，26（1）：139-140.

[17] YASUDA M，ISHIKAWA A，JIKUMARU Y，et al.Antagonistic interaction between systemic acquired resistance and the abscisic acid-mediated abiotic stress response in Arabidopsis[J].Plant cell，2008，20（6）：1678-1692.

[18] 周三，周明，張碩，等.鹽生野大豆的異黃酮積累及其生態(tài)學意義[J].植物生態(tài)學報，2007，31（5）：930-936.

[19] 嚴俊鑫，遲德富，張永強，等.水楊酸誘導重瓣玫瑰對白粉病的抗性[J].東北林業(yè)大學學報，201 41（8）：95-101.

[20] BISWAS C，DEY P，KARMAKAR P G，et al.Next-generation sequencing and micro RNAs analysis reveal SA/MeJA1/ABA pathway genes mediated systemic acquired resistance（SAR）and its master regulation via production of phased，trans-acting siRNAs against stem rot pathogen Macrophomina phaseolina in a RIL population of jute（Corchorus capsularis）[J].Physiological and molecular plant pathology，201 87：76-85.