番茄?擬南芥 PREs 及水稻 ILIs 基因生物信息學分析

郭鵬宇 楊志杰

摘要 為進一步了解常見模式植物番茄、擬南芥 PREs 以及水稻 ILIs 共 18 個基因的生物信息學數據，為基因功能的研究奠定理論基礎，結合 NCBI 等數據庫，運用 MG2C 等工具，對上述基因結構、蛋白理化性質等生物信息學數據作出預測與分析。除 OsILI6 外，其余基因均只含1個內含子，且 CDS 序列均較短。蛋白理化性質分析表明這些蛋白質穩定性較低，二級結構分析表明 α 螺旋與無規則卷曲構成蛋白質的主體部分。三維結構模擬表明這些蛋白質以二聚化的形式發揮功能，在結構上相對保守，分析數據可為后續基因功能研究提供支持。

關鍵詞 番茄;擬南芥;水稻;PREs;生物信息學分析

中圖分類號 Q 812? 文獻標識碼 A

文章編號 0517-6611（2021）18-0099-06

doi：10.3969/j.issn.0517-6611.2021.18.025

開放科學（資源服務）標識碼（OSID）：

Bioinformatics Analysis of Tomato，Arabidopsis PREs and Rice ILIs Genes

GUO Peng-yu，YANG Zhi-jie （Bioengineering College，Chongqing University，Chongqing 400044）

Abstract In order to further understand the bioinformatics data of 18 genes in common model plants tomato，Arabidopsis PREs and rice ILIs，and lay a theoretical foundation for the study of gene function，combined with NCBI and other databases，using MG2C and other tools，bioinformatics data such as the above gene structure and protein physical and chemical properties were predicted and analyzed .Except for OsILI6，all other genes contain only one intron，and the CDS sequence was relatively short.The analysis of physicochemical properties showed that these proteins had low stability and the analysis of the secondary structure showed that α helix and random coils constitute the main part of the protein.Three-dimensional structural simulations indicated that these proteins function in the form of dimerization.These proteins were relatively conservative in structure，and the analytical data could provide support for subsequent gene function studies.

Key words Tomato;Arabidopsis thaliana;Rice;PREs;Bioinformatics analysis

作者簡介 郭鵬宇（1995—），男，山西呂梁人，碩士研究生，研究方向：生物化學與分子生物學。

收稿日期 2021-02-04

番茄是管狀花目的一年生或多年生草本植物。據考證，番茄最早起源于南美的秘魯、厄瓜多爾，于 16 世紀末或 17世紀初的明萬歷年間傳入我國[1]。1753年，Linnaeus 根據雄蕊數目和雌雄蕊的著生習性將番茄劃歸于茄科茄屬中。目前，番茄已成為研究果實發育的重要模式植物[2]。作為一種兼具營養價值和商業價值的園藝作物，近年來世界范圍內的番茄產量不斷提升[3]。 擬南芥是十字花科的一年生細弱草本植物，分布廣泛，植株雖小卻結子較多。早在 2000 年，其基因組測序就已全部完成，其遺傳背景較為簡單，且作為一種自花授粉植物，基因高度純合，是公認的基因組學研究的模式植物[4]。

水稻屬禾本科稻屬作物，是世界公認的三大糧食作物之一。2002 年，我國宣布完成水稻基因組精細圖的繪制。憑借其較為簡單的遺傳背景及遺傳轉化的潛力，水稻成為谷類作物育種及基因功能研究的模式植物之一[5]。研究發現，擬南芥中存在一類能夠編碼對赤霉素合成抑制劑多效唑產生拮抗作用蛋白質的基因，稱為多效唑抗性基因（PACLOBUTRAZOL RESISTANCE GENE），簡稱 PREs，該類蛋白質為非典型的bHLH 轉錄因子。擬南芥 PRE 基因家族共有 6 個成員，即 AtPRE1-AtPRE6，在調控擬南芥赤霉素[6]、油菜素內酯[7]、生長素響應[8]以及調節細胞長度[9]、介導光信號[10]響應等方面發揮重要作用。番茄中共有 5 個PRE基因，即SlPRE1～SlPRE5[11]。Slstyle2.1，即 SlPRE1 參與調控花柱長度，并且可以促進栽培品種自花授粉的進化[12]。SlPRE2 的超表達會影響番茄植株形態的變化[11]，而SlPRE2-RNAi 會影響胎座以及果皮細胞的大小，改變植株對外源赤霉素的響應[13]。水稻共有7個PRE基因，即OsILI1～OsILI7。AtPRE1 同源基因 ILI1（INCREASED LAMINA INCLINATION1） 的過表達會導致水稻葉片傾斜角度增加，通過與水稻 IBH1 的相互作用參與油菜素內酯信號轉導來調控細胞的伸長[9]。ILI4 的過量表達增加了水稻籽粒的大小以及對油菜素內酯合成抑制劑的抗性[14]。該研究旨在利用生物信息學的方法對擬南芥、番茄 PRE家族、水稻 ILI 家族基因結構、染色體定位、蛋白質理化性質分析、進化關系等方面進行分析與預測，結合已有的功能研究報，為尚未進行功能研究的基因提供理論基礎。

1 材料與方法

1.1 鑒定與染色體定位

從茄科基因組數據庫（https：∥solgenomics.net/）、擬南芥基因組數據庫（https：∥www.arabidopsis.org/）及水稻基因組數據庫（http：∥rice.plantbiology.msu.edu/）下載得到全部蛋白質序列、cDNA 序列、CDS（coding-sequence）序列以及各基因在基因組上的序列。運行在線軟件MG2C（http：∥mg2c.iask.in/mg2c_v2.1/）繪制染色體定位圖譜。

1.2 基因結構分析

從各物種基因組數據庫下載各基因的 CDS 序列以及該基因在基因組上 DNA 序列的 FASTA 格式文件。運行在線工具 Genes Structure Display Server 2.0（http：∥gsds.gao-lab.org/index.php）繪制各基因編碼區及非編碼區序列分布圖。

1.3 蛋白質保守基序分析

根據各物種蛋白質序列，運行在線軟件 MEME（https：∥meme-suite.org/meme/），分析各物種蛋白質序列中的保守基序，motif 個數設置為 10，其他參數為默認值。

1.4 蛋白保守結構域序列比對和系統進化樹構建

利用 DNAMAN 6.0 軟件完成所有蛋白質序列的比對分析，利用軟件 MEGA 6.0 完成系統進化樹的構建。

1.5 蛋白理化性質與亞細胞定位預測構建

運行在線軟件 ExPasy （https：∥www.expasy.org/）分析各編碼蛋白質的理化性質，運行在線數據庫 GenScript（https：∥www.genscript.com/psort.html？src=leftbar），對各蛋白質進行亞細胞定位預測。

1.6 蛋白二級結構分析

運行在線軟件 SOPMA （https：∥npsa-prabi.ibcp.fr/cgi-bin/npsa_automat.pl？page=npsa_sopma.html）分析各蛋白質的二級結構組成。

1.7 蛋白三維結構模擬

根據同源建模法，運行 Swiss-Model 軟件（https：∥swissmodel.expasy.org/interactive）對各蛋白質的三維結構進行模擬。

1.8 啟動子區域順式作用元件分析

從各物種基因組數據庫下載得到各基因上游 3 000 bp 序列的 FASTA 文件，利用在線網站 Plant CARE （http：∥bioinformatics.psb.ugent.be/webtools/plantcare/html/）進行順式作用元件分析，利用 TBtools 1.072 繪制順式作用元件分布圖。

1.9 蛋白互作網絡預測分析

根據 STRING 11.0 （https：∥string-db.org/cgi/input？sessionId=bABx2NvfkwQD&input_page_active_form=single_sequence）在線數據庫進行各蛋白質互作網絡的預測。

2 結果與分析

2.1 鑒定與染色體定位

根據已有的文獻報道，結合 SGN （茄科基因組數據庫Solanaceae Genomics Network）、TAIR （擬南芥數據庫The Arabidopsis Information Resource）、RGAP （水稻基因組數據庫Rice Genome Annotation Project） 以及 NCBI （National Center for Biotechnology Information Search database）數據庫檢索分析，分別找到番茄 5 個 PREs、擬南芥 6個 PREs、水稻 7 個 ILIs 基因共 18 個相關基因。在各物種基因組數據庫中查找得到各基因在基因組染色體上的位置（圖1），結果表明，番茄的 5 個PRE 基因分布于第 2、4、5、6 染色體上，擬南芥的 PRE 基因分布于第 1、3、5 染色體上，而水稻的 ILIs 基因分布于第 2、3、4、6、10、11 染色體上。

2.2 基因結構

對番茄、擬南芥 PREs、水稻 ILIs 基因進行基因組序列與編碼序列分析（圖2），結果表明，除 OsILI6 的編碼序列含 3 個外顯子，其余基因外顯子數量均為 2;此外，OsILI3 與 OsILI4 只含有外顯子與內含子，而其他基因均包含 5′與 3′非編碼序列，且個別基因非編碼序列較長，如 SlPRE3 的 3′端以及 OsILI6 的 5′端。這表明該家族基因結構相對簡單。

2.3 蛋白質保守基序

對番茄、擬南芥 PREs 與水稻 ILIs 家族的蛋白質序列進行保守基序的預測分析（圖3），結果表明，這些蛋白質包含 3～4 個保守基序;所有序列均包含 Motif 1，而 OsILI6 無 Motif 2，SlPRE3、AtPRE3、AtPRE4 無 Motif 3，而 SlPRE4、OsILI4、OsILI5、OsILI7 則無 Motif 4。保守基序中字母越高，表明對應氨基酸殘基在該位點出現的頻率越大，相對保守程度也越高。綜合來看，4 個 Motif 的保守程度相對較高。這些保守基序的分析有利于深入分析蛋白質結構。

49卷18期??? 郭鵬宇等 番茄·擬南芥 PREs 及水稻 ILIs 基因生物信息學分析

2.4 蛋白保守結構域序列比對和系統進化樹分析

通過對番茄、擬南芥 PREs、水稻 ILIs基因家族蛋白質進行氨基酸組成的分析（圖 4A），結果表明，這些蛋白質存在有相對的氨基酸序列組成的保守結構域，即 HLH 結構域，而無典型的 Bas-ic 結構域，而在進化過程中，這些基因的功能高度保守與上述保守結構域的構成密切相關。

為了進一步明確這些基因之間的進化關系，利用氨基酸序列對上述 18 個基因構建系統進化樹（圖4B）? 。目前番茄 SlPRE1 與 SlPRE2 已有功能研究，且與 AtPRE1 及 AtPRE5 親緣關系較近。總體來看，番茄與擬南芥親緣關系較近，而二者與水稻的親緣關系相對較遠。

2.5 蛋白理化性質與亞細胞定位預測分析

利用 ExPasy 在線軟件對上述 18 個基因的編碼蛋白質進行理化性質分析（表 1）。結果表明，18 個蛋白質的氨基酸殘基數為86～130，其中 SlPRE4 最長，OsILI6 最短;蛋白質分子量為9.746～14.215;理論等電點為4.43～9.61，其中 SlPRE1、AtPRE2 最大，而 OsILI6 最小;脂溶指數為81.08～106.26，其中 OsILI6 最小，OsILI7 最大;不穩定系數為60.29～93.81;蛋白質的疏水性在 -0.766～-0.218，表明這些蛋白質均屬親水性蛋白;亞細胞定位預測結果表明這些蛋白均定位于細胞核內。

2.6 蛋白二級結構分析

運行 SOPMA 在線軟件對 18 個基因所編碼的蛋白質進行二級結構解析（表 2）。結果表明，在上述 18 個蛋白質的二級結構中，α 螺旋和無規則卷曲占主體地位;部分蛋白質無 β轉角及延伸鏈（如SlPRE4、AtPRE1、AtPRE2、AtPRE3、AtPRE4、AtPRE6、OsILI1、OsILI7）;除α螺旋與無規則卷曲外，個別蛋白質只含 β 轉角與延伸鏈之中的一種結構，如SlPRE2、SlPRE5、AtPRE5、OsILI5 無延伸鏈，而OsILI2、OsILI3、OsILI4 則無 β 轉角。總體來看，α 螺旋和無規則卷曲是這些蛋白質的主要組成成分。

2.7 蛋白三維結構模擬 為了進一步了解上述 18 個蛋白質的三維空間結構，根據同源建模法運行 Swiss-Model 在線軟件進行同源建模（圖 5），結果表明，這些蛋白質均以二聚化的形式完成建模。

2.8 基因啟動子區域順式作用元件分析

通過 SGN數據庫、TAIR 數據庫、RGAP 數據庫獲得各基因 ATG 上游 3 000 bp 的序列，運行在線數據庫 Plant CARE 網站對序列進行順式作用元件分析。結果顯示（圖6），這些基因上游存在大量的光響應元件 （G-box、GT1-motif、ACE），除此之外，激素響應元件（GARE-motif、ABRE、CGTCA-motif、AuxRR-core、P-box）、逆境脅迫響應元件（WUN-motif、TC-rich repeats、LTR、MBS、DRE）、分生組織表達元件（CAT-box）、胚乳表達元件（GCN4_motif）及種子特異性調控元件（RY-element）也存在與上游啟動子區域。SlPRE3、AtPRE4～AtPRE6、OsILI2～OsILI7 啟動子區域均含有干旱響應元件;除 SlPRE3、SlPRE4、AtPRE1、AtPRE2、OsILI1、OsILI2、OsILI7 外，其余基因啟動子區域均含有低溫響應元件，部分基因啟動子區域含有類黃酮物質合成基因調控元件（SlPRE3、AtPRE4），推測這2個基因有可能參與類黃酮物質的合成。

2.9 基因蛋白互作網絡預測分析

為進一步了解上述基因可能的互作蛋白，該部分參考 STRING 在線數據庫，對這些基因可能的互作網絡進行了預測。結果表明，除 OsILI2 與 OsILI3 外，其余蛋白質均有預測的蛋白質網絡互作圖。這表明大部分基因可以通過蛋白質的相互作用參與生長過程的調控（圖7）。

3 討論與結論

截至目前，已有功能研究報道的 PREs 以及水稻 ILIs 均屬于非典型的 bHLH 轉錄因子[8，14-15]。這些基因通過赤霉素、油菜素內酯[16]以及生長素參與植物營養生長與生殖生長[17]的調控。目前，擬南芥 AtPRE1～AtPRE6、SlPRE1～SlPRE2 以及 OsILI1、OsILI4 均有功能研究報道，而其余基因暫無功能研究。該研究通過生物信息學的方法，結合已有的功能研究報道，對番茄、擬南芥 PREs 以及水稻 ILIs 基因進行生物信息學分析，為這些基因功能的進一步研究奠定理論基礎。

該研究對3個物種共 18 個蛋白質的理化性質進行了初步分析，發現這些蛋白質均屬于疏水性蛋白質，且穩定性較差;此外，α 螺旋和無規則卷曲是這些蛋白質主要的二級結構組分，僅有較少的 β 轉角與延伸鏈存在其中，從三維結構模擬可以更直觀地發現這一點。蛋白質的三維空間模擬結果表明，這些蛋白質可以通過二聚化的形式參與調控植株生長。與其他蛋白質的互作網絡預測表明，這些蛋白通過與其他蛋白質的相互作用發揮調控作用。蛋白質保守基序分析與蛋白質多序列比對及進化關系分析表明這些蛋白質在調控植物生長發育方面存在功能相似性。

通過對3個物種共 18 個基因進行生物信息學分析，從染色體定位、基因結構、蛋白質理化性質以及二級、三級結構分析等方面做出分析預測，為進一步的功能研究提供理論基礎與數據支撐。

參考文獻

[1]

趙凌俠，李景富.番茄起源、傳播及分類的回顧[J].作物品種資源，1999（3）：29-31.

[2] QUINET M，ANGOSTO T，YUSTE-LISBONA F J，et al.Tomato fruit development and metabolism[J].Frontiers in plant science，2019，10：1-23.

[3] CUI J，JIANG N，ZHOU X X，et al.Tomato MYB49 enhances resistance to Phytophthora infestans and tolerance to water deficit and salt stress[J].Planta，2018，248（6）：1487-1503.

[4] MEINKE D W，CHERRY J M，DEAN C，et al.Arabidopsis thaliana：A model plant for genome analysis[J].Science，1998，282（5389）：662-682.

[5] LI X X，DUAN X P，JIANG H X，et al.Genome-wide analysis of basic/helix-loop-helix transcription factor family in rice and Arabidopsis[J].Plant physiology，2006，141（4）：1167-1184.

[6] LEE S，LEE S，YANG K Y，et al.Overexpression of PRE1 and its homologous genes activates gibberellin-dependent responses in Arabidopsis thaliana[J].Plant & cell physiology，2006，47（5）：591-600.

[7] BAI M Y，SHANG J X，OH E，et al.Brassinosteroid，gibberellin and phytochrome impinge on a common transcription module in Arabidopsis[J].Nature cell biology，2012，14（8）：810-817.

[8] CASTELAIN M，LE HIR R，BELLINI C.The non-DNA-binding bHLH transcription factor PRE3/bHLH135/ATBS1/TMO7 is involved in the regulation of light signaling pathway in Arabidopsis[J].Physiologia plantarum，2012，145（3）：450-460.

[9] ZHANG L Y，BAI M Y，WU J，et al.Antagonistic HLH/bHLH transcription factors mediate brassinosteroid regulation of cell elongation and plant development in rice and Arabidopsis[J].Plant cell，2009，21（12）：3767-3780.

[10] HYUN Y，LEE I.KIDARI，encoding a non-DNA Binding bHLH protein，represses light signal transduction in Arabidopsis thaliana[J].Plant molecular biology，2006，61（1/2）：283-296.

[11] ZHU Z G，CHEN G P，GUO X H，et al.Overexpression of SlPRE2，an atypical bHLH transcription factor，affects plant morphology and fruit pigment accumulation in tomato[J].Scientific reports，2017，7：1-11.

[12] CHEN K Y，CONG B，WING R，et al.Changes in regulation of a transcription factor lead to autogamy in cultivated tomatoes[J].Science，2007，318（5850）：643-645.

[13] ZHU Z G，LIANG H L，CHEN G P，et al.The bHLH transcription factor SlPRE2 regulates tomato fruit development and modulates plant response to gibberellin[J].Plant cell reports，2019，38（9）：1053-1064.

[14] TANAKA A，NAKAGAWA H，TOMITA C，et al.BRASSINOSTEROID UPREGULATED1，encoding a helix-loop-helix protein，is a novel gene involved in brassinosteroid signaling and controls bending of the lamina joint in rice[J].Plant physiology，2009，151（2）：669-680.

[15] IKEDA M，FUJIWARA S，MITSUDA N，et al.A triantagonistic basic helix-loop-helix system regulates cell elongation in Arabidopsis[J].Plant cell，2012，24（11）：4483-4497.

[16] WANG H，ZHU Y Y，FUJIOKA S，et al.Regulation of Arabidopsis brassinosteroid signaling by atypical basic helix-loop-helix proteins[J].The plant cell 2009，21（12）：3781-3791.

[17] MARA C D，HUANG T B，IRISH V F.The Arabidopsis floral homeotic proteins APETALA3 and PISTILLATA negatively regulate the BANQUO genes implicated in light signaling[J].Plant cell，2010，22（3）：690-702.