水稻響應熱脅迫核心基因的篩選與鑒定

張斌 楊昕霞 袁志輝

摘要： 基于熱脅迫條件下的水稻轉錄組數據，通過HTSeq和DESeq軟件篩選差異表達基因;對其進行功能富集和構建蛋白質互作網絡，鑒定最重要模塊中的核心基因。結果顯示：水稻熱脅迫0 h、1 h、6 h和12 h條件下，共有278個差異表達基因，生物學過程主要富集在刺激反應和蛋白質折疊上;鑒定出含有12個基因的重要模塊.基因表達模式顯示重要模塊中7個核心基因受熱脅迫誘導上調，推測這7個核心基因在水稻響應熱脅迫過程中發(fā)揮關鍵作用。

關鍵詞： 水稻;高溫脅迫;核心基因;轉錄組

中圖分類號： S511.01?? 文獻標識碼： A?? 文章編號： 1000-4440（2021）04-0817-06

Screening and identification of core genes responding to heat stress in rice

ZHANG Bin， YANG Xin-xia， YUAN Zhi-hui

（Hunan University of Science and Engineering， Yongzhou 425199， China）

Abstract： In this study， rice transcriptome data under heat stress were used to screen differentially expressed genes by HTSeq and DESeq software， then functional enrichment analysis was carried out and protein interaction network was constructed to identify the core genes in most important module. The results showed that there were 278 differentially expressed genes under heat stress for 0 h， 1 h， 6 h and 12 h， and the biological process was mainly concentrated in stimulation response and protein folding. The important module containing 12 genes was identified. The gene expression patterns showed that seven core genes in important modules were up-regulated under heat stress， suggesting that these seven core genes played a key role in the response of rice to heat stress.

Key words： rice;high temperature stress;core gene;transcriptome

水稻作為全球重要的糧食作物，其生長過程中極易受到各種生物脅迫和非生物脅迫的影響，其中高溫脅迫嚴重影響水稻的產量和稻米的品質。因此，研究者在水稻熱脅迫的危害及其應答熱脅迫的生理和分子機制方面做了大量的工作。在成熟期，高溫導致水稻籽粒灌漿提前結束[1];高溫還可以引起稻米品質下降[2]。植物可以啟動熱激響應基因表達從而抵御脅迫傷害[3]。擬南芥鈣調蛋白AtCaM3激活鈣離子/鈣調蛋白結合蛋白激酶，鈣離子內流激活鈣離子依賴的蛋白激酶，進而磷酸化熱脅迫反應的關鍵調控因子，從而啟動下游基因表達[4]。熱激蛋白在植物響應熱脅迫過程中發(fā)揮重要作用[5-7]。隨著研究的深入，人們發(fā)現植物響應熱脅迫的分子機制極其復雜。Takehara等發(fā)現OsSUS3基因導入日本晴后提高了該基因的表達水平，增強了抗熱性[8]。Moon等報道異源表達OsHSP1基因可提高擬南芥的耐熱性[9]。但是，到目前為止，水稻響應熱脅迫的生理和分子調控機理仍未被闡述清楚。因此，本研究利用轉錄組學手段探究水稻響應熱脅迫過程中的關鍵基因，并進一步闡釋其功能途徑。

1 材料與方法

1.1 材料

熱脅迫處理的水稻轉錄組原始數據從NCBI數據庫（https：//www.ncbi.nlm.nih.gov/bioproject/PRJNA530826/）下載。本研究進行試驗的水稻品種為日本晴 （Nipponbare）。

1.2 差異表達基因篩選

利用Hisat2軟件建立索引以及比對到參考基因組，通過Stringtie軟件進行轉錄本組裝、合并和定量，基因差異分析則通過DESeq2軟件進行。以|log2foldChange|>2且P值<0.01為條件分別篩選高溫脅迫1 h、6 h和12 h時間點的差異表達基因，再通過Ggplot2繪制所有差異基因的火山圖和韋恩圖。

1.3 蛋白質互作網絡構建及重要模塊分析

利用在線數據庫（http：//string-db.org），以參數combined score>0.4，構建蛋白質互作網絡， 并通過Cytoscape進行可視化。參數設置：scores>5，degree cut-off=2，node scorecut-off=0.2，Max depth=100和k-score=2。此外，利用MCODE（分子復合物檢測）軟件篩選最重要模塊。

1.4 核心基因篩選及表達分析

設置馬修斯相關系數（MCC）算法，利用Cytoscape軟件的cytoHubba插件鑒定相互作用最緊密的10個基因，與最重要模塊基因取交集，篩選核心基因。取log2TPM值（TPM為每千堿基記錄本），分析水稻在高溫脅迫0 h、1 h、6 h、12 h、24 h和48 h時間點的基因表達模式。

1.5 差異表達基因富集分析

差異表達基因使用在線軟件（http：//systemsbiology.cau.edu.cn/agriGOv2/）進行GO富集分析，取FDR（錯誤發(fā)現率）值≤0.05的通路，利用Ggplot2插件進行繪圖。差異表達基因使用在線軟件（http：//kobas.cbi.pku.edu.cn/kobas3/genelist/）進行KEGG富集分析，取P值≤0.05的通路，利用Ggplot2插件繪圖。

1.6 qRT-PCR驗證

種子發(fā)芽后培養(yǎng)14 d（28 ℃，14 h白天/10 h黑夜），45 ℃熱激1 h取樣。利用TRIzol法提取幼苗總RNA，添加脫氧核糖核酸酶I（Dnase I）去除DNA污染。使用反轉錄試劑盒SuperScript III First-Strand Synthesis SuperMix（Thermo）合成cDNA。利用ABI7500定量PCR儀進行qRT-PCR檢測，利用2-△△Ct法計算基因相對表達水平。反應程序：94 ℃預變性32 min;94 ℃變性10 s，60 ℃退火10 s，72 ℃延伸10 s，40個循環(huán)。總體系20 μl：10.0 μl SYBR Green Mix（Thermo），8.0 μl RNAase-free水，1.0 μl引物，1.0 μl模板。Os10g0510000基因為內參，引物序列見表1。

2 結果與分析

2.1 水稻熱脅迫下差異表達基因篩選

熱脅迫1 h，檢測到2 642個差異表達基因，1 624個上調，1 018個下調（圖1A）;熱脅迫6 h，檢測到1 704個差異表達基因，812個上調，892個下調（圖1B）;熱脅迫12 h，檢測到1 773個差異表達基因，899個上調，874個下調（圖1C）。共有278個重疊基因（圖1D）。

2.2 水稻熱脅迫下差異表達基因富集分析

278個基因富集到GO項上，生物學過程主要富集在刺激反應和蛋白質折疊;分子功能主要富集在結合功能;細胞組分無顯著富集（圖2A）。KEGG分析結果顯示278個差異基因主要富集在內質網蛋白質加工和次生代謝產物合成通路上（圖2B）。

2.3 蛋白質互作網絡構建與重要模塊中關鍵核心基因篩選

通過蛋白質互作網絡（圖3A）獲得了一個含有12個基因表達全部上調的最重要模塊（圖3B）。進一步通過Cytoscape軟件的cytoHubba插件鑒定出作用最緊密的10個基因（圖3C），其中7個基因存在最重要模塊上，確定為核心基因（圖3D）。

2.4 水稻響應熱脅迫核心基因表達模式分析

分析核心基因的表達模式，發(fā)現在熱脅迫條件下不同時間點（1 h、6 h、12 h、24 h和48 h）基因表達都表現上調（圖4）。隨熱脅迫時間的延長，基因表達水平有一定波動。Os01g0840100、Os07g0641700、Os12g0277500、Os03g0366000和Os02g0612000 5個基因的表達在熱脅迫1 h時達到最大值，呈現先上升，后下降的趨勢;Os08g0338700和Os09g0491772響應熱脅迫時間比較長，在48 h內基因表達表現出上升的趨勢。

2.5 qRT-PCR方法驗證轉錄組測序結果

利用qRT-PCR方法對轉錄組數據進行驗證。結果（圖5）顯示，野生型水稻日本晴幼苗經熱脅迫處理1 h后，篩選出的7個核心基因Os12g0277500、Os08g0338700、Os07g0641700、Os09g0491772、 Os01g0840100、Os02g0612000和Os03g0366000的表達水平都上調，結果與轉錄組測序數據基本一致，表明轉錄組測序結果可靠，也進一步證明了通過轉錄組數據篩選鑒定水稻響應熱脅迫過程中核心基因的可靠性。

3 討論

生物具有感知溫度變化的系統(tǒng)[10]。植物響應熱脅迫涉及受體激活、活性氧產生、熱激蛋白和熱激轉錄因子等多種信號途徑[11]，過程復雜是其分子機制沒有研究清楚的原因之一。因此，迫切要求判斷植物響應熱脅迫的潛在核心基因。轉錄組技術能夠分析熱脅迫中植物基因表達的改變，是一種確定核心基因有效方法。

本研究分析了水稻熱脅迫下4個轉錄組數據集，共鑒定出278個差異表達基因，GO分析結果顯示這些基因主要富集在對刺激的反應和蛋白質折疊的生物過程上，KEGG分析結果顯示這些基因主要富集在內質網蛋白質加工的代謝途徑上。蛋白質折疊在生物適應高溫脅迫的過程中起著重要的作用[12]。熱激蛋白（HSP）分為HSP100、HSP90、HSP70、HSP60和小HSP家族，具有阻止變性蛋白質聚集和協(xié)助變性蛋白質重新折疊的功能[13]。HSP101基因功能缺失的擬南芥變得對熱敏感[14]。因此， HSP家族成員在植物響應熱脅迫過程中發(fā)揮重要作用。葉綠體和線粒體HSP70s作為轉運子的一部分可以幫助蛋白質前體轉移[15-17]。研究結果表明擬南芥通過細胞周期轉錄因子介導葉綠體AtHSP70-4基因的表達適應溫度變化[18]。HSP70-4和E3泛素連接酶通過泛素-26S蛋白酶系統(tǒng)介導質體靶向蛋白質前體的降解，減輕細胞損傷[19]。線粒體mtHsc70-1基因在擬南芥細胞色素c氧化酶（COX）依賴性呼吸系統(tǒng)中發(fā)揮作用，基因敲除后表現出嚴重的生長缺陷[20]。Os01g0840100和Os09g0491772同屬HSP70家族，在水稻熱脅迫響應過程中的具體功能未見報道，本研究中其受熱脅迫誘導，可能發(fā)揮調控作用。質體Cpn60是線粒體Hsp60的同源物，參與核酮糖-1，5-二磷酸羧化酶的組裝，在此過程中需要蛋白質GroES/Cpn10進行底物封裝[21-23]。Cpn60家族在擬南芥生長發(fā)育過程中發(fā)揮著非常重要的作用，短日照條件下，擬南芥AtCpn60β1基因缺失導致幼苗死亡，突變體對熱脅迫敏感[24] 。OsCpn60α1（Os12g0277500）對核酮糖-1，5-二磷酸羧化酶大亞基的折疊至關重要，突變體表現出淡綠色和幼苗致死的表型，高溫對OsCpn60α1的轉錄有較強的誘導作用[25]。同源基因OsCpn60β1對水稻葉綠體發(fā)育至關重要[26]。這與我們的分析結果一致，說明OsCpn60在熱應激過程中發(fā)揮了重要的作用。擬南芥AT Cpn10基因與大腸桿菌groES基因功能互補，在不同器官中mRNA都有表達，受熱脅迫誘導[27]。Os07g0641700和Os03g0366000同屬Cpn10基因，表達水平在0～24 h持續(xù)上調，說明在熱脅迫過程中發(fā)揮了一定的作用。DnaK-DnaJ-GrpE系統(tǒng)是蛋白質穩(wěn)態(tài)的重要組成部分[28]。擬南芥GrpE表達受熱脅迫誘導，與大腸桿菌grpE基因功能互補[29]。Os08g0338700和Os02g0612000同屬GrpE基因，此基因的功能亦未見報道，本研究中受熱脅迫的誘導，推測其在熱脅迫響應中發(fā)揮調控作用。水稻響應熱脅迫與蛋白質折疊息息相關，7個核心基因同屬HSP家族，蛋白質互作網絡分析結果表明它們聯系非常緊密，同源蛋白質參與蛋白質的運輸、組裝和折疊，因此，我們認為這些基因在熱脅迫過程中協(xié)同發(fā)揮關鍵性作用。

總之，本研究旨在確定水稻熱脅迫響應過程中的核心基因。共鑒定出278個差異基因和7個核心基因，7個基因可能在水稻響應熱脅迫過程中發(fā)揮重要作用。但是，需要進一步研究來闡明這些基因在水稻響應熱脅迫過程中的生物學功能。

參考文獻：

[1] KIM J， SHON J， LEE C K， et al. Relation-ship between grain filling duration and leaf senescence of temperate rice under high temperature[J]. Field Crops Research， 2011， 122： 207-213.

[2] ZHONG L J， CHENG F M， WEN X， et al. The deterioration of eating and cooking quality caused by high temperature during grain filling in early-season indica rice cultivars[J]. Journal of Agronomy & Crop Science， 2005， 191： 218-225.

[3] LIM C J， YANG K A， HONG J K， et al. Gene expression profiles during heat acclimation in Arabidopsis thaliana suspension-culture cells[J]. J Plant Res， 2006， 119： 373-383.

[4] SUZUKI N， SEJIMA H， TAM R， et al. Identification of the MBF1 heat-response regulon of Arabidopsis thaliana[J]. Plant J， 2011， 66（5）： 844-851.

[5] GUAN J C， YEH C H， LIN Y P， et al. A 9 bp cis-element in the promoters of class I small heat shock protein genes on chromosome 3 in rice mediates L-azetidine-2-carboxylic acid and heat shock responses[J]. Journal of Experimental Botany， 2010， 61（15）： 4249-4261.

[6] JAGADISH S V， MUTHURAJAN R， OANE R， et al. Physiological and proteomic approaches to address heat tolerance during anthesis in rice（Oryza sativa L.）[J]. Journal of Expermental Botany， 2010， 61： 143-156.

[7] JUNG K H， GHO H J， NGUYEN M X， et al. Genome-wide expression analysis of HSP70 family genes in rice and identification of a cytosolic HSP70 gene highly induced under heat stress[J]. Functional & Integrative Genomics， 2013， 13（3）： 391-402.

[8] TAKEHARA K， MURATA K， YAMAGUCHI T， et al. Thermo-responsive allele of sucrose synthase 3 （Sus3） provides high-temperature tolerance during the ripening stage in rice （Oryza sativa L.）[J]. Breeding Science， 2018， 68（3）： 336-342.

[9] MOON J C， HAM D J， HWANG S G， et al. Molecular characterization of a heat inducible rice gene， OsHSP1， and implications for rice thermotolerance[J]. Genes Genom， 2014， 36： 151-161.

[10]MITTLER R， FINKA A， GOLOUBINOFF P. How do plants feel the heat？[J]. Trends Biochem Sci， 2012， 37（3）： 118-125.

[11]陳思婷，郭房慶.植物耐熱性及熱激信號轉導機制研究進展[J].中國科學：生命科學，2013，43（12）： 1072-1081.

[12]MORO F， MUGA A. Thermal adaptation of the yeast mitochondrial Hsp70 system is regulated by the reversible unfolding of its nucleotide exchange factor[J]. J Mol Biol，2006， 358（5）： 1367-1377.

[13]WANG W， VINOCUR B， SHOSEYOV O， et al. Role of plant heat-shock proteins and molecular chaperones in the abiotic stress response[J]. Trends Plant Sci， 2004， 9（5）： 244-252.

[14]HONG S W， VIERLING E. Mutants of Arabidopsis thaliana defective in the acquisition of tolerance to high temperature stress[J]. Proc Natl Acad Sci USA， 2000， 97（8）： 4392-4397.

[15]JACKSON C D， KEEGSTRA K. Arabidopsis genes encoding components of the chloroplastic protein import apparatus[J]. Plant Physiol， 2001， 125：1567-1576.

[16]SU P H， LI H M. Stromal Hsp70 is important for protein translocation into pea and Arabidopsis chloroplasts[J].Plant Cell， 2010， 22：1516-1531.

[17]ZHANG X P， GLASER E. Interaction of plant mitochondrial and chloroplast signal peptides with the Hsp70 molecular chaperone[J]. Trends Plant Sci， 2002， 7：14-21.

[18]ZHOU S， SUN H， ZHENG B， et al. Cell cycle transcription factor E2F2 mediates non-stress temperature response of AtHSP70-4 in Arabidopsis[J]. Biochem Biophys Res Commun， 2014， 455（3-4）： 139-146.

[19]LEE S， LEE D W， LEE Y， et al. Heat shock protein cognate 70-4 and an E3 ubiquitin ligase， CHIP， mediate plastid-destined precursor degradation through the ubiquitin-26S proteasome system in Arabidopsis[J]. Plant Cell， 2009， 21（12）： 3984-4001.

[20]WEI S S， NIU W T， ZHAI X T， et al. Arabidopsis mtHSC70-1 plays important roles in the establishment of COX-dependent respiration and redox homeostasis[J]. J Exp Bot， 2019， 70（20）： 5575-5590.

[21]BARRACLOUGH R， ELLIS R J. Protein synthesis in chloroplasts IX assembly of newly-synthesized large subunits into ribulose bisphosphate carboxylase in isolated intact pea chloroplasts[J]. Biochim Biophys Acta， 1980， 608（1）：19-31.

[22]HEMMINGSEN S M， ELLIS R J. Purification and properties of ribulosebisphosphate carboxylase large subunit binding protein[J]. Plant Physiology， 1986， 80（1）：269-276.

[23]HEMMINGSEN S M， WOOLFORD C， VAN D， et al. Homologous plant and bacterial proteins chaperone oligomeric protein assembly[J]. Nature， 1988， 333（6171）：330-334.

[24]ATSUSHI I， HIDEAKI T， MASATO N， et al. Deletion of a chaperonin 60 beta gene leads to cell death in the Arabidopsis lesion initiation 1 mutant.[J]. Plant and Cell Physiology， 2003（3）：255-261.

[25]KIM S R， YANG J I， AN G. OsCpn60α1， encoding the plastid chaperonin 60α subunit， is essential for folding of rbcL[J]. Moleculer Cells， 2013， 35（5） ：402-409.

[26]WU Q， ZHANG C， CHEN Y， et al. OsCpn60β1 is essential for chloroplast development in rice （Oryza sativa L.）[J]. International Journal of Molecular Sciences， 2020， 21（11） ：4023.

[27]KOUMOTO Y， TSUGEKI R， SHIMADA T， et al. Isolation and characterization of a cDNA encoding mitochondrial chaperonin 10 from Arabidopsis thaliana by functional complementation of an Escherichia coli groES mutant[J]. Plant J，1996， 10（6）： 1119-1125.

[28]BEN-ZVI A P， GOLOUBINOFF P. Review： mechanisms of disaggregation and refolding of stable protein aggregates by molecular chaperones[J]. J Struct Biol， 2001， 135（2）： 84-93.

[29]HU C， LIN S Y， CHARNG C， et al. Recent gene duplication and subfunctionalization produced a mitochondrial GrpE， the nucleotide exchange factor of the Hsp70 complex， specialized in thermotolerance to chronic heat stress in Arabidopsis[J]. Plant Physiology， 2012， 158（2）：747-758.

（責任編輯：張震林）

收稿日期：2020-01-25

基金項目：湖南省自然科學基金項目（2020JJ4030）;湖南省創(chuàng)新平臺與人才計劃項目（2020NK4222）

作者簡介：張 斌（1981-），男，湖南永州人，博士，講師，主要從事植物發(fā)育生物學研究，（E-mail）zhangbin27104@163.com

通訊作者：袁志輝，（E-mail）zhh_yuan@126.com