沙田柚F—box蛋白基因的克隆及序列分析

摘要：采用生物信息學方法，對沙田柚[Citrus maxima （Burm.） Merr. cv. Shatian Yu.]F-box蛋白基因編碼的蛋白質(zhì)從序列特征、理化性質(zhì)、跨膜結構域、高級結構及功能域等方面進行了預測和分析。結果表明，該基因全長為1 427 bp（GenBank登錄號為KR363148），開放閱讀框（ORF）全長為1 101 bp，共編碼366個氨基酸，編碼蛋白質(zhì)的分子質(zhì)量43.09 ku，理論等電點5.15。沙田柚F-box蛋白基因編碼蛋白含有一個植物F-box蛋白家族的保守結構域，為親水性非分泌不穩(wěn)定蛋白，不跨膜運動，共有30個可能的磷酸化位點。氨基酸序列分析表明，其編碼的氨基酸與甜橙（Citrus sinensis，KDO71185）和克萊門柚（Citrus clementina，XP_006425495）F-box蛋白的同源性分別為99%、98%。系統(tǒng)進化樹表明，沙田柚F-box蛋白基因與與甜橙（Citrus sinensis，KDO71185）和克萊門柚（Citrus clementina，XP_006425495）親緣關系很近，屬于同一進化分支。

關鍵詞：沙田柚[Citrus maxima （Burm.） Merr. cv. Shatian Yu.]；F-box蛋白基因；序列分析

中圖分類號：S576 文獻標識碼：A 文章編號：0439-8114（2016）14-3723-04

DOI：10.14088/j.cnki.issn0439-8114.2016.14.046

Abstract： The characters of the F-box protein of Citrus maxima（Burm.） Merr. cv. Shatian Yu. were analyzed and predicted by bioinformatics method，including the composition of amino acid sequence，physicochemical parameters，hydrophobicity，transmembrane domain，secondary structure and function of protein etc. The results showed that the F-box gene was 1 427 bp（GenBank accession number： KR363148） in length with an open reading frame（ORF） of 1 101 bp，encoding 366 amino acids with deduced molecular weight of 43.09 ku，and theoretical pI value of 5.15，and contained an F-box domain in N end. Bioinformatics analysis showed that the F-box protein was a hydrophilic and unstable protein， has no transmembrane domain，and there are 30 phosphorylation sites within the polypeptide chain. The homology analysis of amino acid sequence indicated that the F-box protein shared high homology with that of Citrus sinensis，KDO71185（99%） and Citrus clementina， XP_006425495（98%）. Phylogenetic analysis revealed the F-box protein gene showed closer kinship with that of Citrus sinensis and Citrus clementina，indicating that they belong to the same evolutionary branch.

Key words： Citrus maxima （Burm.） Merr. cv. Shatian Yu.； F-box protein gene； sequence analysis

自交不親和是植物界的一種普遍現(xiàn)象，根據(jù)其遺傳背景可分為孢子體自交不親和（Sporophytic self-incompatibility，SSI）和配子體自交不親和（Gemetophytic self-incompatibility，GSI）兩種類型。植物自交不親和性在被子植物的有性生殖中起到避免近親繁殖的重要作用[1]，闡明自交不親和的分子機理在培育無子品種等方面具有重要的意義。近年來，薔薇科（Rosaseae）、茄科（Solanaseae）、玄參科（Scrophulariaceae）等自交不親和的研究取得了較大進展，已確定S-核酸酶（S-RNase）為花柱S-決定子，S-核酸酶（S-RNase）是一類具有核酸酶活性的堿性蛋白[2，3]。在花粉決定因子的研究方面，Lai等[4]采用S位點區(qū)域染色體步移法在自交不親和的玄參科金魚草（Antirrhinum）S-locus位點首次獲得一個在花粉組織特異表達的F-box基因AhSLF-S2（Antirrhinum S-locus F-box S2）。Sijacic等[5]在茄科矮牽牛（Petunia inflata）中發(fā)現(xiàn)的的S-locus F-box基因PiSLF2。薔薇科植物李屬梅（Prunus mume）以及扁桃（P. dulcis）中的自交不親和S型（S-haplotype）特異性的S-locus F-box基因Pm-SLF和PdSFB的克隆，李屬甜櫻桃（Prunus avium）不同型的S-locus中6個SFB基因，PaSFB1-6的發(fā)現(xiàn)，表明介導GSI反應的花粉S-決定基因可能是這些基因[6～8]。

沙田柚[Citrus maxima（Burm.） Merr. cv. Shatian Yu.]屬于蕓香科配子體高度自交不親和果樹。薛妙男等[9]確定了沙田柚為配子體自交不親和，花粉管生長抑制部位在花柱的1/2，分析了沙田柚花柱蛋白以及花柱在不同發(fā)育時期的幾種同工酶特征[10，11]；秦新民等[12]對沙田柚花粉管特異蛋白進行了分離和鑒定；薛妙男等[13]確定了花柱通道細胞中的特異蛋白；秦新民等[14，15]研究了花粉管中特異蛋白產(chǎn)生部分和分布位置。為了進一步探討沙田柚自交不親和的機理，對自交和異交花柱進行了轉(zhuǎn)錄組測序，通過自交與異交花柱差異基因的比對，獲得了沙田柚F-box蛋白基因，并對該基因編碼的蛋白質(zhì)進行了生物化學特征分析，旨在為沙田柚自交不親和分子機理的深入研究提供幫助。

1 材料與方法

1.1 材料

沙田柚試驗材料采自廣西靈川縣潮田鄉(xiāng)大山口村果園十年生結果樹。在盛花期對沙田柚樣樹進行人工自交（沙田柚×沙田柚）授粉和異交授粉（酸柚×沙田柚），分別收集自交和異交1～3 d的授粉花柱以及當天未開花的花柱，立即放入液氮中保存，并放入-80 ℃超低溫冰箱備用。

1.2 方法

1.2.1 RNA的提取、建庫和測序 按改良Trizol法對總RNA進行提取[16]，由深圳華大基因科技服務有限公司對檢測合格的RNA進行建庫和測序。用帶有Oligo（dT）的磁珠富集mRNA，并隨機將mRNA打斷成片段，然后以打斷的RNA為模板，由隨機引物合成第1條cDNA鏈，加入緩沖液、RNase H、DNA polymeraseⅠ以及dNTPs合成第2條cDNA鏈。使用試劑盒純化cDNA，用EB緩沖液洗脫后進行末端修復和加上poly（A）以及連接測序接頭，使用瓊脂糖凝膠電泳對cDNA文庫的片段大小進行檢測，再進行PCR擴增，用Illumina HiSeqTM 2000對制備好的文庫進行測序。

1.2.2 序列分析和系統(tǒng)樹構建 使用DNAman、ORF Finder、TMPRED、SWISS-MODEL、NetPhos 2.0等軟件對序列以及蛋白質(zhì)理化性質(zhì)進行分析，使用DNAman軟件對氨基酸序列F-box蛋白構建系統(tǒng)進化樹。

2 結果與分析

2.1 基因的生物信息學分析

沙田柚F-box蛋白基因（CL8619.Contig2_All）的基因組序列全長為1 427 bp（GenBank登錄號為KR363148），通過NCBI網(wǎng)站上的NCBI ORF Finder和生物信息學軟件DNAman進行分析，該序列包含一個1 101 bp的開放閱讀框（ORF），編碼366個氨基酸（圖1）。

2.2 編碼蛋白質(zhì)的分析和疏水性預測

通過DNAMan軟件分析，該蛋白質(zhì)的分子質(zhì)量43.09 ku，理論pI 5.15。從圖2中可以看出，所編碼的肽鏈疏水性最大值約為2.689，最小值約為-2.544，平均值為-0.233，屬于親水性蛋白。運用在線分析軟件（http：//www.expasy.org/tools/protparam.html）分析可知，該蛋白質(zhì)帶負電荷的氨基酸（Asp+Glu）總數(shù)為51，帶正電荷的氨基酸（Arg+Lys）總數(shù)為40，分子式為C1 964H2 966N492O561S20，不穩(wěn)定系數(shù)為46.99，屬于不穩(wěn)定性蛋白。

2.3 功能結構域分析

功能結構域分析結果（圖3）表明，將編碼的氨基酸序列用NCBI的Conserved Domain Architecture Retrieval Tool分析，發(fā)現(xiàn)該蛋白質(zhì)具有一個與F-box super-family蛋白相同的保守結構域。

2.4 跨膜預測

運用跨膜數(shù)據(jù)庫TMbase（http：//ch.embnet.org/software/TMPRED_form.html）預測F-box蛋白的跨膜區(qū)域。從圖4可以看出，該蛋白不存在跨膜結構蛋白，為膜內(nèi)蛋白。

2.5 蛋白質(zhì)磷酸化預測

通過在線蛋白磷酸化位點分析軟件NetPhos 2.0 Server，對F-box基因所編碼的蛋白進行預測。結果（圖5）表明，該多肽鏈中分值在0.5以上的可能的磷酸化位點共有30個。其中，絲氨酸（Ser）可能的磷酸化位點共有14個，分別位于肽鏈的8、13、28、38、130、196、206、213、240、242、286、288、298和328位；蘇氨酸（Thr）可能的磷酸化位點共有4個，分別位于肽鏈的7、62、148和311位；酪氨酸（Tyr）可能的磷酸化位點共有12個，分別位于肽鏈的9、40、50、67、71、96、134、184、218、233、302和339位。

2.6 蛋白質(zhì)二級結構及三級結構的預測

通過在線分析軟件Predict Protein（http：//www.predictprotein.org/）對F-box蛋白二級結構進行預測，結果表明，該蛋白含有α-螺旋為8.47%，β-折疊 31.15%，其他結構為60.38%。

利用SWISS-MODEL在線軟件對F-box蛋白三級結構進行預測，結果如圖6所示，F(xiàn)-box蛋白的三級結構含有3個α-螺旋、9個β折疊，其間由無規(guī)則卷曲連接。

2.7 不同物種間F-box蛋白同源性比較及進化分析

將沙田柚F-box蛋白的氨基酸序列與從GenBank數(shù)據(jù)庫中下載其他10種植物F-box基因編碼的氨基酸序列進行比對（圖7），結果表明，沙田柚F-box蛋白基因編碼的氨基酸與甜橙（Citrus sinensis，KDO71185）和克萊門柚（Citrus clementina，XP_006425495）F-box蛋白的同源性分別為99%和98%。利用DNAMan構建系統(tǒng)發(fā)育樹，結果表明沙田柚F-box蛋白基因編碼的蛋白質(zhì)與蕓香科的甜橙和克萊門柚有很近的親緣關系，屬于同一個進化分支（圖7）。

2.8 F-box蛋白基因在沙田柚自交和異交花柱中的表達

轉(zhuǎn)錄組測序結果分析表明，該基因在自交花柱和異交花柱中的轉(zhuǎn)錄水平存在明顯差異（圖8）。在未授粉花柱中其表達量為0.73，在自交1 d花柱（Z1）中迅速上升（1.61），2 d（Z2）仍維持較高水平（1.40），3 d（Z3）則迅速下降（0.54）。而在異交1 d的花柱（Y1）中，該基因的表達迅速下降（0.25），2 d（Y2）出現(xiàn)回升（0.73），3 d（Y3）則快速下降（0.13）。

利用RPKM（Reads Per Kb per Million reads）[17]對該基因在自交和異交花柱中的表達水平估算，設定FDRS≤0.001且差異倍數(shù)≥2倍的基因視作顯著差異表達基因，對以RPKM值取2的對數(shù)值（log2（RPKM））對自交1 d/異交1 d、自交2 d/異交2 d和自交3 d/異交3 d花柱中F-box基因的表達水平進行統(tǒng)計分析，其log2（RPKM ratio）分別為2.69、0.96和2.05。

3 討論

F-box蛋白是一類廣泛存在于真核生物，含有F-box結構域的蛋白家族[18]，在細胞周期調(diào)控、轉(zhuǎn)錄調(diào)控、細胞凋亡和信號轉(zhuǎn)導、激素響應等過程中發(fā)揮著重要的功能[19-21]。目前，人們已從多種植物中克隆了F-box蛋白基因，但尚未見關于沙田柚F-box蛋白基因克隆的報道。本研究克隆的CL8619.Contig2_All序列與甜橙和克萊門柚F-box蛋白基因編碼的氨基酸序列同源性分別為99%和98%，同時在N-端含有一個保守的F-box結構域；說明CL8619.Contig2_All序列是一個F-box蛋白基因。

F-box基因在植物自交不親和中也可能有重要作用。Lai等[4]首先在金魚草中克隆了一個有花粉組織特異性表達的F-box蛋白基因（AhSLF-S2）。隨后，在茄科矮牽牛中獲得了S-locus F-box基因PiSLF2[5]，在薔薇科梅和扁桃中分別鑒定出了特異性的S-locus F-box基因PmSLF和PdSFB[6-8]，相關研究表明，這些F-box基因可能是介導植物配子體自交不親和反應，與其自交不親和性密切相關[6-8，22]。沙田柚為配子體自交不親和果樹，本研究克隆的F-box在沙田柚自交花柱和異交花柱中的表達存在明顯差異：在自交1 d花柱中其表達量由0.73迅速上升至1.61，2 d仍維持較高水平（1.40），3 d則迅速下降到0.54。而在異交1 d的花柱中，該基因的表達迅速下降0.25，2 d出現(xiàn)回升（0.73），3 d則快速下降至0.13。自交1 d/異交1 d，自交3 d/異交3 d的log2（RPKM ratio）分別為2.69和2.05，達到差異表達的水平。但該F-box蛋白基因在沙田柚自交不親和反應中的功能尚有待進一步證實。

參考文獻：

[1] TAKAYAMA S， ISOGAI A. Self-incompatibility in plants[J]. Annual Review of Plant Biology，2005，56：467-489.

[2] SASSA H，HIRANO H，IKEHASHI H.Self-incompatibility-related RNase in styles of Japanese pear（Pyrus serotina Rehd.）[J]. Plant Cell Phsiology，1992，33：811-814.

[3] ANDERSON M A，CORNISH E C，MAU S-L，et al. Cloning of cDNA for a stylar glycoprotein associated with expression of self-incompatibility in Nicotiana alata[J].Nature，1986，321：38-44.

[4] LAI Z，MA W S，HAN B，et al.An F-box gene linked to the self-incompatibility（S） locus of Antirrhinum is expressed specifically in pollen and tapetum[J]. Plant Molecular Biology，2002，50：29-42.

[5] SIJACIC P，WANG X，SKIRPAN A L，et al. Identification of the pollen determinant of S-RNase-mediated self-incompatibility[J].Nature，2004，429：302-305.

[6] ENTANI T，IWANO M，SHIBA H，et al. Comparative analysis of the self-incompatibility（S-） locus region of Prunus mume： Identification of a pollen-expressed F-box gene with allelic diversity[J].Genes to Cells，2003，8：203-213.

[7] YAMANE H，IKEDA K，USHIJIMA K，et al. A pollen-expressed gene for a novel protein with an F-box motif that is very tightly linked to a gene for S-RNase in two species of cherry， Prunus cerasus and P. avium[J].Plant Cell Physiclogy，2003， 44：764-769.

[8] IKEDA K，IGIC B，USHIJIMA K，et al. Primary structural features of the S haplotype-specific F-box protein， SFB， in Prunus[J]. Sexual Plant Reproduction，2004，16：235-243.

[9] 薛妙男，陳騰土，楊繼華.沙田柚自交和異交親和性觀察[J].園藝學報，1995，22（2）：127-132

[10] 楊繼華，薛妙男，陸景峰.沙田柚自交、異交花柱蛋白質(zhì)的比較分析[J].植物學通報，1996，13（增刊）：45.

[11] 李紅艷，楊繼華，薛妙男.沙田柚花柱三種同工酶分析[J].廣西師范大學學報（自然科學版），1999，17（4）：80-84.

[12] 秦新民，李惠敏，薛妙男，等.沙田柚自交，異交花粉管蛋白的雙向電泳分析[J].廣西植物，2004，24（6）：566-569.

[13] 薛妙男，李義平，張杏輝，等.沙田柚自交花柱中識別蛋白的免疫金定位[J].園藝學報，2001，28（1）：59-61.

[14] 秦新民，莫花濃，萬 珊，等.沙田柚花粉管特異蛋白的免疫細胞化學研究[J].廣西師范大學學報，2008，26（4）：113-115.

[15] 秦新民，莫花濃，石菁萍，等.沙田柚花粉管S1-RNase免疫膠體金定位研究[J].廣西農(nóng)業(yè)科學，2009，40（5）：483-485.

[16] BARIOLA P A，HOWARD C J，TAYLOR C B，et al. The Arabidopsis thaliana ribonuclease gene RNA1 is tightly controlled in response to phosphate limitation[J].Plant J，1994，6：673-685.

[17] MORTAZAVI A，WILLIAMS B A，MCCUE K，et al. Mapping and quantifying mammalian transcriptomes by RNA-Seq[J].Nature Methods，2008，5（7）：621-628.

[18] KIPREOS E T，PAGANO M. The F-box protein family[J].Genome Biol，2000，1（5）：1-7.

[19] NAKAYAMAK K I，NAKAYAMAK K. Ubiquitin ligases：Cell cycle control and cancer[J].Nat Rev Cancer，2006，6（5）：369-381.

[20] MONTE E，AMADOR V，RUSSO E，et al. A U-Box GA signaling component with a role proteasome degradation[J].Journal Plant Growth Regulation，2003，22：152-162.

[21] YU H C，WU J，XU N F，et al. Roles of F-box proteins in plant hormone responses[J].Acts Biochim Biophys Sin，2007， 39（12）：915-922.

[22] QIAO H，WANG H Y，ZHAO L，et al.The F-box protein AhSLF-S2 physically interacts with S-RNases that may be inhibited by the ubiquitin/26S proteasome pathway of protein degradation during compatible pollination in antirrhinum[J].Plant Cell，2004，16：582-595.