利用高通量测序对沙田柚自交和异交花柱进行转录组测序,通过差异分析得到了沙田柚谷胱甘肽转移酶基因( CmGST)序列。 CmGST全长为872 bp (GenBank登录号为KY612456),开放阅读框为657 bp,编码219个氨基酸,推测蛋白质分子量为25.00 kDa,理论等电点pI为5.48。利用生物信息分析软件对CmGST基因进行同源性比对和系统进化分析表,结果显示该基因编码的氨基酸具有谷胱甘肽S转移酶典型的C端与 N端双结构域,属于GST Tau家族蛋白;与克莱门柚( Citrus maxima)的相似度约为98%。这些分析结果可为今后深入研究该蛋白的结构特征和功能提供参考。 In this paper, the transcriptome of the self-pollinated style and cross-pollinated style of Citrus maxima cv. Shatian Yu were sequences by high-throughput sequencing technology. The glutathione S-transferase gene ( CmGST) sequence of Citrus maxima cv. Shatian Yu was obtained through differential analysis method. The fu11 length of CmGST gene was 872 bp (GenBank accession number: KY612456), and the ORF (Open Reading Frame) length was 657 bp, encoding a protein polypeptide of 219 amino acids with a predicted molecular weight of 25.00 kDa and pI of 5.48. Bioinformatics analysis showed that the CmGST protein contained the Tau GST-specific N-terminal domain (G site) and the C-terminal domain (H site). The homology analysis of amino acid sequence indicated that the glutathione S-transferase protein shared high homology with glutathione S- transferase protein of Citrus clementina (98%). This work provides a systemic sequence analysis of the glutathione S-transferase protein of Citrus maxima cv. Shatian Yu. It will provide the useful reference for further investigation of its structure and function.
罗琼,肖君,李少梅,龙小霞,刘华英,秦新民*
广西师范大学生命科学学院,广西 桂林
收稿日期:2017年3月7日;录用日期:2017年3月25日;发布日期:2017年3月28日
利用高通量测序对沙田柚自交和异交花柱进行转录组测序,通过差异分析得到了沙田柚谷胱甘肽转移酶基因(CmGST)序列。CmGST全长为872 bp (GenBank登录号为KY612456),开放阅读框为657 bp,编码219个氨基酸,推测蛋白质分子量为25.00 kDa,理论等电点pI为5.48。利用生物信息分析软件对CmGST基因进行同源性比对和系统进化分析表,结果显示该基因编码的氨基酸具有谷胱甘肽S转移酶典型的C端与 N端双结构域,属于GST Tau家族蛋白;与克莱门柚(Citrus clementina)的相似度约为98%。这些分析结果可为今后深入研究该蛋白的结构特征和功能提供参考。
关键词 :沙田柚,谷胱甘肽转移酶基因,生物信息学
Copyright © 2017 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
植物自交不亲和性(Self-incompatibility, SI)是指能产生具有正常功能且同期成熟的雌雄配子的雌雄同体植物,在自花授粉或相同基因型异花授粉时无法受精的现象 [
SI是目前所知的最主要的一种显花植物控制受精的机制,是由雌蕊细胞和花粉之间相互作用产生的生理反应。这两种独特细胞间的接触是受空间和时间的授粉行为控制的,这远不同于其他细胞间的信号转导 [
沙田柚为无患子目、芸香科、柑橘属植物,为严格的配子体自交不亲和植物。关于沙田柚自交不亲和机理的研究取得一定得进展,薛妙男等确定沙田柚为配子体自交不亲和类型,花粉管在花柱1/2处生长受到抑制;通过免疫胶体金技术,确定自交1-3天花柱通道细胞S-蛋白产生部位及分布位置 [
为了进一步探讨沙田柚自交不亲和的机理,我们对自交和异交花柱进行了转录组测序,通过自交与异交花柱差异基因的比对,获得了一个沙田柚谷胱甘肽转移酶基因,并对该基因编码的蛋白质进行了生物化学特征分析,旨在为沙田柚自交不亲和分子机理的深入研究提供帮助。
沙田柚实验材料采自广西灵川县潮田乡大山口村果园十年生结果树。在盛花期对沙田柚样树进行人工自交(沙田柚 × 沙田柚)授粉和异交授粉(酸柚♂ × 沙田柚♀),分别收集自交和异交1~3 d的授粉花柱以及当天未开花的花柱,立即放入液氮中保存,并放入−80℃超低温冰箱备用。
总RNA的提取按改良Trizol法进行 [
使用DNAman、Finder、TMPRED、SWISS-MODEL、NetPhos2.0、SignaIP4.1等软件进行序列的开放阅读框(ORF)、编码氨基酸、编码蛋白的理化分析。将测定的基因序列与GenBank中的13种植物的同源序列进行比对,基于氨基酸序列谷胱甘肽S转移酶基因的系统进化树用DNAman软件进行构树。
谷胱甘肽S转移酶基因(Unigene16804_All)全长872 bp (GenBank登录号为KY612456),包含一个657 bp的开放阅读框(ORF)编码219个氨基酸(图1)。
通过在线软件protparam分析(http://www.expasy.org/tools/protparam.html),该基因基因编码的蛋白质分子质量为25.00 kDa;理论等电点为5.48;共由3517个原子组成,分子式为C1149H1753N283O322S10;不稳定指数为33.42 (<40),为稳定蛋白。
用DNAMAN预测该蛋白质的疏水性:该基因编码的蛋白质疏水性最大值为3.04,亲水性最小值为−3.11。综合表明该蛋白质疏水性平均值为−1.93,为亲水蛋白(图2)。
跨膜预测的结果显示该基因编码的蛋白质存在着2个跨膜区。跨膜方向由内向外有149~167位氨基酸;跨膜方向为由外向内的有6~22和143~165位氨基酸(图3)。
该蛋白质有两个保守结构域,一是GST_C_superfamily (谷胱甘肽转移酶超家族),位于该蛋白质的第85~215位,二是Thioredoxin_like superfamily (硫氧还蛋白超家族),位于该蛋白质的第1~75位(图4)。
该基因编码的蛋白质可能的磷酸化位点共有16个,其中Ser磷酸化位点11个,分别位于肽链的32、67、82、93、115、149、174、190、194、196、198位;Thr磷酸化位点1个,位于肽链第169位;Tyr磷酸化位点共有5个,分别位于肽链30、90、101、107和206位。表明该蛋白的磷酸化以Ser磷酸化为主,兼有Thr和Tyr磷酸化(图5)。
图1. Unigene16804_All序列和推测氨基酸序列
图2. Unigene16804_All蛋白疏水性分析
图3. Unigene16804_All跨膜区预测
图4. Unigene16804_All蛋白保守结构域
图5. Unigene16804_All蛋白的磷酸化位点分析
利用SignaIP4.1软件对沙田柚谷胱甘肽转移酶基因编码的氨基酸序列进行分析,结果表明该蛋白质没有明显的信号肽。
通过在线分析软件Predict Protein (http://www.predictprotein.org/)对谷胱甘肽转移酶蛋白二级结构进行预测,结果显示该蛋白含有螺旋、折叠和无规则卷曲,所占比例分别是54.34%、5.48%、40.18%。
利用SWISS-Model对沙田柚谷胱甘肽转移酶蛋白三级结构进行预测,结果表明:该蛋白三级结构含有14个α-螺旋,7个β-折叠,其间由无规则卷曲连接(图6)。
从GenBank数据库中下载13种植物谷胱甘肽转移酶基因编码的氨基酸序列,序列进行比对结果表明沙田柚谷胱甘肽转移酶基因编码的氨基酸与芸香科的克莱门柚(Citrus clementina, XP_006432070.1)谷胱甘肽转移酶蛋白的同源性为98%。利用DNAMan构建系统发育树,结果表明沙田柚谷胱甘肽转移酶基因编码的蛋白质与芸香科的克莱门柚(Citrus clementina)有很近的亲缘关系,属于同一个进化分支(图7)。
氨基酸序列同源性分析表明,本文所克隆的沙田柚基因(Unigene16804_All)与克莱门柚谷胱甘肽转移酶的同源性均为98%,通过对该基因编码蛋白的功能结构域分析,发现该蛋白具有谷胱甘肽转移酶保守结构域,表明Unigene16804_All属于谷胱甘肽转移酶超家族基因。
图6. Unigene16804_All蛋白的三级结构
图7. 基于氨基酸序列的谷胱甘肽转移酶蛋白系统发育树
植物的谷胱甘肽转移酶(GST)是一个多基因家族,植物的GST分为φ、τ、ζ、θ、λ和脱氢抗坏血酸还原酶(DHARs) 6类 [
本文克隆的谷胱甘肽转移酶基因在沙田柚未授粉,以及自交花柱和异交花柱中的表达存在较为明显的差异:未授粉花柱中该基因的表达量(FPKM)为0.106,自花授粉1,2,3 d花柱中的表达量(FPKM)分别为3.32、5.98、2.45,异花授粉1,2,3 d花柱中基因的表达量(FPKM)则分别为2.26、3.51和14.16。从上述结果可以看出无论自交授粉还是异交授粉,第1 d花柱中谷胱甘肽转移酶基因的表达都表现为迅速升高,在自交授粉状态下,基因的表达第2 d稍许升高,第3 d则下降。而异交授粉过程中,该基因的表达则表现为持续升高。
沙田柚谷胱甘肽转移酶基因的这种表达模式与其参与花粉-雌蕊的相互作用和影响花粉管的极性生长的功能相符。薛妙男等 [
国家自然科学基金(31360477);广西教育厅项目(2013YB036)。
罗琼,肖君,李少梅,龙小霞,刘华英,秦新民. 沙田柚谷胱甘肽转移酶基因(CmGST)的生物信息学及表达分析 Bioinformatics and Expression Analysis of GST Gene Encoding Glutathione S-Transferase in Citrus maxima cv. Shatian Yu[J]. 植物学研究, 2017, 06(02): 47-54. http://dx.doi.org/10.12677/BR.2017.62008