張偉 張亞龍 臺(tái)德衛(wèi) 張效忠 汪婉琳

摘要 利用水稻苗期轉(zhuǎn)錄組數(shù)據(jù)，分析在高溫處理?xiàng)l件下不同時(shí)間段基因表達(dá)水平變化，篩選出可能參與高溫?zé)岷Φ幕?。LOC_Os04g04970、LOC_Os03g15960和LOC_Os03g51459，其在受到熱處理1 h的表達(dá)水平顯著上調(diào)，隨著熱處理的時(shí)間延長(zhǎng)，其表達(dá)水平逐漸降低，可能與熱處理相關(guān)。該研究可為進(jìn)一步揭示水稻苗期高溫?zé)岷Ψ肿由飳W(xué)研究提供參考。

關(guān)鍵詞 水稻;苗期;高溫脅迫;轉(zhuǎn)錄組

中圖分類號(hào) S 511? 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2021）20-0120-03

doi：10.3969/j.issn.0517-6611.2021.20.030

開放科學(xué)（資源服務(wù)）標(biāo)識(shí)碼（OSID）：

High Temperature RNA-SEQ Analysis of Rice Seedling

ZHANG Wei ZHANG Ya-long TAI De-wei1 et al （1.Rice Research Institute，Anhui Academy of Agricultural Science，Hefei，Anhui 230031;2.University of Washington，WA，U.S.A 98195）

Abstract This study used rice seedling transcriptome data to analyze the changes in gene expression levels at different time periods under high temperature treatment conditions， and screened out three genes （LOC_Os04g04970， LOC_Os03g15960 and LOC_Os03g51459） with significant changes in expression levels at each time period. Their expression level was significantly up-regulated after 1 hour of heat treatment， and gradually decreased with the extension of heat treatment， which may be related to heat treatment. This study will provide a reference for further revealing the molecular biology of high temperature heat damage in rice seedlings.

Key words Rice;Seeding stage;Heat stress;RNA-seq

基金項(xiàng)目 安徽省科技重大專項(xiàng)（h202003a06020005）;科技部國(guó)際合作司常規(guī)性科技援助項(xiàng)目（KY202001010）。

作者簡(jiǎn)介 張偉（1985—），男，安徽臨泉人，助理研究員，博士，從事分子生物學(xué)與生物信息學(xué)研究。

收稿日期 2020-12-17;修回日期 2021-01-19

近年來，隨著環(huán)境的變化，高溫天氣頻繁出現(xiàn)，高溫正影響我國(guó)長(zhǎng)江中下游的水稻生產(chǎn)[1]。當(dāng)日最高溫度高于水稻適宜生長(zhǎng)溫度范圍時(shí)，會(huì)影響水稻的生理生長(zhǎng)和生殖生長(zhǎng)，造成植物光合作用系統(tǒng)、細(xì)胞膜穩(wěn)定性、激素等方面發(fā)生變化，甚至導(dǎo)致植物死亡[2]。高溫通過影響植株生理生長(zhǎng)和生殖生長(zhǎng)，最終表現(xiàn)為產(chǎn)量下降，品質(zhì)降低[3-4]。

在水稻生殖生長(zhǎng)階段和成熟階段遇到高溫時(shí)危害較嚴(yán)重，水稻在整個(gè)生長(zhǎng)發(fā)育過程中對(duì)高溫敏感的時(shí)期首先是抽穗開花期[5-7]，花粉成熟、花藥開裂、花粉在柱頭上發(fā)芽和花粉管伸長(zhǎng)等均受到影響，在水稻減數(shù)分裂期如遇35 ℃以上高溫并持續(xù)3 d及以上，會(huì)導(dǎo)致水稻花器發(fā)育不全，花藥開裂率下降，花粉發(fā)育不良，花粉育性和活力下降等[8-9]，最終影響結(jié)實(shí)率。抽穗至成熟階段高溫會(huì)加快灌漿速率，縮短灌漿的持續(xù)期，籽粒光合產(chǎn)物不足，淀粉及其他有機(jī)物積累減少，使籽粒的充實(shí)度（千粒重）受到影響，米粒堊白增加，透明度變差[10]。高溫不但在孕穗期、開花期、灌漿期等生育時(shí)期影響水稻生產(chǎn)，在營(yíng)養(yǎng)生長(zhǎng)時(shí)期受到高溫時(shí)對(duì)植株的生長(zhǎng)也產(chǎn)生重要影響，如導(dǎo)致植株生長(zhǎng)速度下降，葉色變淡等[11]，尤其是在水稻幼苗時(shí)期，受到高溫脅迫時(shí)水稻葉片中蛋白種類和含量變化較大，植株高度、干物質(zhì)含量性狀也會(huì)發(fā)生變化[12-13]，對(duì)水稻植株產(chǎn)生過氧化傷害作用[14]，進(jìn)而影響水稻植株后期的生長(zhǎng)進(jìn)程。

當(dāng)前很多研究集中于水稻生殖生長(zhǎng)階段，對(duì)生理生長(zhǎng)階段研究較少，水稻苗期是水稻生理生長(zhǎng)重要階段之一，而關(guān)于水稻苗期受到高溫處理其基因表達(dá)方面的研究較少。隨

著生物信息學(xué)與分子生物學(xué)的發(fā)展，為水稻基因組、轉(zhuǎn)錄組、生理生化等方面的研究提供了技術(shù)保障，也為揭示水稻耐高溫研究注入新的動(dòng)力。該研究利用水稻苗期高溫處理的RNA-seq數(shù)據(jù)，使用Hisat2+Deseq2分析表達(dá)水平差異顯著的基因[15-16]，進(jìn)而預(yù)測(cè)和發(fā)掘水稻苗期與高溫處理相關(guān)基因，旨在為水稻高溫相關(guān)基因進(jìn)一步研究提供參考。

1 材料與方法

1.1 數(shù)據(jù)來源 從 The Rice Annotation? Project （RAP）數(shù)據(jù)庫(kù)下載 Os-Nipponbare-Reference- IRGSP-1.0 的水稻基因組序列作為參考基因組和基因組注釋文件（http：∥rice.plantbiology.msu.edu/downloads.shtml）。從NCBI（https：∥www.ncbi.nlm.nih.gov/sra）利用sratoolkit下載并拆分水稻的 SRA（PRJNA530826）數(shù)據(jù)。

1.2 試驗(yàn)方法 利用trimmomatic過濾低質(zhì)量的reads和接頭，用fastqc評(píng)估數(shù)據(jù)質(zhì)量是否可以用于后續(xù)分析。利用Hisat2將reads比對(duì)到基因組，接著利用stringtie將比對(duì)上的reads進(jìn)行轉(zhuǎn)錄組組裝，利用stringtie merge進(jìn)行合并，合并后對(duì)組裝的轉(zhuǎn)錄本進(jìn)行表達(dá)水平的計(jì)算。利用prepDE.py將計(jì)數(shù)轉(zhuǎn)換為Deseq2輸入文件。利用Deseq2進(jìn)行基因差異表達(dá)分析，篩選條件（閾值）：絕對(duì)log2FC> 并且Padj<0.0 用ggplots作圖。每個(gè)時(shí)段均有差異表達(dá)基因，利用Python code將每個(gè)時(shí)段均差異顯著的基因取出，進(jìn)行表達(dá)量計(jì)算分析，用pheatmap繪heatmap圖。

2 結(jié)果與分析

從NCBI數(shù)據(jù)庫(kù)下載水稻日本晴14 d苗期莖和葉受到熱處理的RNA seq原始數(shù)據(jù)（表1），處理時(shí)間分別為1、3、6、12、24、36和48 h，每個(gè)時(shí)間段3次重復(fù)。去除接頭和低質(zhì)量的reads后獲得高質(zhì)量的reads用于下一步分析。

利用Hisat2軟件建立日本晴基因組索引，利用Hisat2分別將reads比對(duì)到基因組上，運(yùn)行stringtie獲得含有count的gtf文件，利用prepDE.py轉(zhuǎn)換得到的數(shù)據(jù)用于下一步統(tǒng)計(jì)分析。

Deseq2分析差異表達(dá)基因，得到每個(gè)處理時(shí)段差異表達(dá)的基因，表明當(dāng)受到熱處理時(shí)，一系列相關(guān)基因表達(dá)水平發(fā)生變化，如圖1所示。進(jìn)一步分析發(fā)現(xiàn)，有些基因在所有處理時(shí)段均上調(diào)或者下調(diào)。利用Python code提取這些基因的counts，用于分析這些基因在受到熱處理時(shí)不同時(shí)期的表達(dá)水平（圖2、3），表明這些基因在受到熱處理后，其表達(dá)水平顯著下調(diào)或者上調(diào)。其中，LOC_Os04g04970、LOC_Os03g15960和LOC_Os03g51459 3個(gè)基因在受到熱處理1 h的表達(dá)水平顯著上調(diào)，隨著熱處理的時(shí)間延長(zhǎng)，其表達(dá)水平逐漸降低，可能與熱處理相關(guān)。LOC_Os03g15960編碼heat shock protein，另外2個(gè)功能鮮見報(bào)道，仍需要進(jìn)一步研究。LOC_Os02g35750和LOC_Os03g11910在處理3 h表達(dá)水平顯著上調(diào)，隨著處理時(shí)間的延長(zhǎng)，其表達(dá)水平逐漸降低。LOC_Os02g35750可能在葉發(fā)育早期與葉綠體生成有關(guān)[17]，然而，其表達(dá)水平受到熱處理影響，是否也參與水稻苗期高溫脅迫途徑仍需進(jìn)一步研究。LOC_Os03g11910是Hsp70家族成員，在保護(hù)植物免受非生物脅迫方面發(fā)揮著重要作用[18]，尚需揭示該基因的生物學(xué)功能。圖3顯示，這些基因經(jīng)過高溫處理后，其表達(dá)水平顯著下調(diào)，需要在進(jìn)一步的生物學(xué)研究過程中揭示它們的生物學(xué)功能。

3 討論

該研究分析表明，在水稻苗期高溫脅迫處理?xiàng)l件下，一些基因在各處理時(shí)段均上調(diào)，其中有些基因在處理1 h表達(dá)水平顯著上調(diào)，隨著處理時(shí)間的推移，其表達(dá)水平逐漸降低，它們與已報(bào)道的功能基因功能相符，這進(jìn)一步證明該研究的可靠性。采用取交集的方式，發(fā)掘每個(gè)時(shí)期均上調(diào)或下調(diào)基因，在方法上具有創(chuàng)新性，同時(shí)增加了發(fā)掘基因的可靠性，為苗期耐熱基因功能研究提供有價(jià)值的參考。該研究所用材料為日本晴，作為水稻參考材料適用于多方面的研究，然而若選用耐熱的極端材料進(jìn)行研究，將能更加充分發(fā)掘耐熱關(guān)鍵基因，將繼續(xù)開展這方面的研究。該研究利用多種處理時(shí)段的數(shù)據(jù)挖掘水稻苗期耐熱關(guān)鍵基因，為水稻苗期耐熱基因功能與機(jī)制分析提供基礎(chǔ)具有實(shí)際意義。

4 結(jié)論

該研究分析出3個(gè)基因（LOC_Os04g04970、LOC_Os03g15960和LOC_Os03g51459）很可能與水稻苗期受到熱脅迫相關(guān)。

參考文獻(xiàn)

[1] 黃義德，曹流儉，武立權(quán)，等.2003年安徽省中稻花期高溫?zé)岷Φ恼{(diào)查與分析[J].安徽農(nóng)業(yè)大學(xué)學(xué)報(bào)，200 31（4）：385-388.

[2] 鞠冠華，崔麗潔，張書，等.植物耐熱性的分子機(jī)制研究進(jìn)展[J].長(zhǎng)江蔬菜，2012（24）：5-11.

[3] 唐婷，陳艷艷，楊洋，等.水稻耐高溫性研究進(jìn)展[J]. 分子植物育種，2017，15（9）：3694-3700.

[4] 楊軍，章毅之，賀浩華，等.水稻高溫?zé)岷Φ难芯楷F(xiàn)狀與進(jìn)展[J].應(yīng)用生態(tài)學(xué)報(bào)，2020，31（8）：2817-2830.

[5] 陳慶全，余四斌，李春海，等.水稻抽穗開花期耐熱性 QTL 的定位分析[J].中國(guó)農(nóng)業(yè)科學(xué)，2008，41（2）：315-321.

[6] 李飛，卓壯，U.A.KAPILA SIRI UDAWELA ，等.水稻高溫?zé)岷Πl(fā)生機(jī)理與耐高溫遺傳基礎(chǔ)研究[J].植物遺傳資源學(xué)報(bào)，201 14（1）：97-103.

[7] 周建霞，張玉屏，朱德峰，等.高溫下水稻開花習(xí)性對(duì)受精率的影響[J].中國(guó)水稻科學(xué)，2014，28（3）：297-303.

[8] 楊永杰，符冠富，熊杰，等.高溫對(duì)水稻的影響及水稻耐熱性測(cè)評(píng)方法研究[J]. 中國(guó)稻米，201 18（1）：39-40.

[9]? 鄢圣敏，王開，楊川航，等.自貢市高溫?zé)岷?duì)水稻生長(zhǎng)的影響及對(duì)策探討[J].安徽農(nóng)業(yè)科學(xué)，2020，48（9）：34-36.

[10] 段驊，楊建昌.高溫對(duì)水稻的影響及其機(jī)制的研究進(jìn)展[J]. 中國(guó)水稻科學(xué)，201 26（4）：393-400.

[11] 李金軍，劉明放，陸金根，等.粳稻穗期耐高溫鑒定技術(shù)的研究[J]. 上海農(nóng)業(yè)學(xué)報(bào)，2003，19（4）：25-27.

[12] BARNABS B，JGER K，F(xiàn)EHR A.The effect of drought and heat stress on reproductive processes in cereals[J].Plant，cell & environment，2008，31（1）：11-38.

[13] HAN F，CHEN H，LI X J，et al.A comparative proteomic analysis of rice seedlings under various high-temperature stresses[J].Biochimica et biophysica acta （BBA）-Proteins and proteomics，2009，1794（11）：1625-1634.

[14] 郭培國(guó)，李榮華.夜間高溫脅迫對(duì)水稻葉片光合機(jī)構(gòu)的影響[J].植物學(xué)報(bào)，2000，42（7）：673-678.

[15] PERTEA M，KIM D，PERTEA G M，et al.Transcript-level expression analysis of RNA-seq experiments with HISAT，StringTie and Ballgown[J].Nature protocols，2016，11（9）：1650-1667.

[16] LOVE M I，HUBER W，ANDERS S.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome biology，2014，15（12）：1-21.

[17] WANG Y，REN Y L，ZHOU K N，et al.WHITE STRIPE LEAF4 encodes a novel P-type PPR protein required for chloroplast biogenesis during early leaf development[J].Frontiers in plant science，2017，8：1-17.

[18] WANG Y F，LIN S K，SONG Q，et al.Genome-wide identification of heat shock proteins （Hsps） and Hsp interactors in rice： Hsp70s as a case study[J].BMC Genomics，2014，15（1）：1-15.