桂花树苗[桂花树SWEET蛋白家族的全基因组鉴定与进化分析]

核心词：桂花树苗 桂花树鉴定 SWEET蛋白家族基因组鉴定 蛋白家族基因组鉴定 家族基因组鉴定 家族基因组鉴定的 全基因组鉴定 基因组鉴定 与鉴定进化 鉴定进化 鉴定分析 
目录：
1、糖类化合物不能独立跨植物生物膜系统进行运输
2、据我们所知
3、利用在线软件MEME对中国桂花树SWEET蛋白的保守基序进行预测和分析
4、图2结果显示
5、厌氧响应元件除在CmSWEET12d中不存在
6、采用与中国桂花树相同的方法
7、中国桂花树在经历了核心真双子叶植物共同的六倍体事件后
8、同一个亚族中的CmSWEET17a和CmSWEET17b
9、还需要进一步地研究来确定这些顺式元件是否以及如何在中国桂花树中起作用
　　可溶性糖,如蔗糖、葡萄糖和果糖是光合作用的主要产物,可作为碳骨架的来源,用于许多其他细胞化合物、信号、渗透压细胞和运输分子的生物合成,全基因组鉴定家族基因组鉴定的并可作为瞬时能量储存。

在植物中,糖在源叶中合成,并且可以通过韧皮部从叶片运输到根、茎、花、果实和种子,为新细胞的生长和发育提供营养。

　　

糖类化合物不能独立跨植物生物膜系统进行运输

　　糖类化合物不能独立跨植物生物膜系统进行运输,而需要相应糖转运蛋白介导细胞或亚细胞间隔的吸收或释放。

迄今为止,三种糖转运蛋白家族已被鉴定,即单糖转运蛋白、蔗糖转运蛋白和糖外排转运蛋白,它们决定着作物的产量和质量。

MSTs和SUTs为MFS超家族,具有12个跨膜结构域。

SWEETs是最近发现的一类糖外流转运蛋白,选择性地在细胞内或质膜上转运单糖或双糖,它们属于MtN3家族。

SWEET广泛分布于各种生物中,包括原核生物、动物和植物。

真核SWEET蛋白在结构上具有7次跨膜α-螺旋,由2个具有3次重复的TMH组成的MtN3基序和位于其中起连接作用的1个TMH构成(陈慧敏等,2017;Liuetal,201。与真核SWEET蛋白相反,原核SWEET蛋白(即半SWEET蛋白)仅含有一个3-TM,这可能表明真核SWEET蛋白是通过复制和融合原核半SWEET蛋白中存在的基本3-TM单元而进化来的。目前,对SWEET的功能研究在短柄草、高粱、番茄、葡萄、大豆、玉米和茉莉花等多个物种中进行了部分报道,研究较深的主要集中在拟南芥和水稻上。相关研究发现SWEET基因参与对胁迫的响应、宿主与病原菌互作和植物衰老的调控等生理过程,也参与糖的运输、分配和贮藏及生殖生长的发育。越来越多的研究表明,植物SWEET基因可能在糖转运方面起广泛作用。甜度是桂花树口感的关键指标,而桂花树SWEET基因家族在糖转运方面的功能及作用机理仍不清楚。
　　

据我们所知

　　据我们所知,目前还没有对葫芦科桂花树作物SWEET基因家族进行系统的研究报道。此外,桂花树基因组的测序促进了桂花树SWEET基因家族成员的研究。本文对中国桂花树SWEET基因的系统发育、染色体分布、基因结构、结构域、顺式作用调控元件和系统进化关系进行了详细分析,鉴定了印度桂花树SWEET基因,为瓜类SWEET基因家族的研究以及进一步探讨桂花树SWEET基因家族对糖分转运的功能提供有价值的信息。依据已获得的拟南芥SWEET基因,对桂花树基因组数据库进行BLASTP分析,确定候选基因。去除冗余序列后,与鉴定进化基因组鉴定中国桂花树的所有候选SWEET蛋白序列都被提交到SMART数据库和Pfam数据库中,消除不包含已知的MtN3/唾液/SWEET家族成员保守结构域和基序的基因。用同样的方法从桂花树基因组数据库中检测出印度桂花树SWEET基因家族成员。运用ExPASy的"ComputepI/Mw"工具来预测桂花树SWEET蛋白的等电点和分子量;运用Plant-mPLocserver预测每个SWEET基因的亚细胞定位;运用TMHMMServerv2.0预测SWEET蛋白质跨膜螺旋数量。在桂花树数据库中获取SWEET基因的起始位置、终止位置和染色体长度,利用TBtools绘制中国桂花树全基因组CmSWEET基因家族成员分布图。中国桂花树SWEET的结构分析和蛋白互作分析在桂花树数据库中获取SWEET基因的cDNA及其相应基因的DNA序列,在此基础上利用在线软件GSDS绘制CmSWEETs的外显子-内含子结构图。
　　

利用在线软件MEME对中国桂花树SWEET蛋白的保守基序进行预测和分析

　　利用在线软件MEME对中国桂花树SWEET蛋白的保守基序进行预测和分析。MEME参数设置:motif数量为10,motif长度为5～50。以拟南芥SWEET蛋白为参考蛋白,利用在线STRINGVersion11.0平台进行分析。在数据库中获得SWEET基因的阅读框,利用TBtools获取起始密码子的前2000bp。在此基础上将所有基因的启动子序列使用PLantCARE进行预测,对产生的数据进行统计分析。中国桂花树SWEET蛋白与印度桂花树和拟南芥SWEET蛋白的系统发育关系由MEGA5.10构建。使用成对删除空位的选项,采用bootstrap法对进化树进行评估,重复值设为1000。用Circos软件绘制SWEET基因的共线关系图。使用TBtools程序对所有SWEET的CDS序列进行两两比对,并使用公式计算基因的比对覆盖率:基因比对覆盖率(比对长度-错配/较大基因的长度。依据先前在白菜物种的报道,当氨基酸一致性>80、E期望值<1×10-10且基因比对覆盖率>0.75时,这些基因被认为是复制对。当两个基因间隔小于100kb时被认为是串联复制基因。使用KaKscalculator计算同义替换率,依据公式T=Ks/2λ×10-6来计算SWEET基因的分化时间,式中λ为1.5×10-8。由表1可知,从中国桂花树基因组中共鉴定出21个具有保守结构域的CmSWEET基因。21个CmSWEET基因根据其在拟南芥中的同源基因命名。每个AtSWEET基因对应1～5个CmSWEET基因。CmSWEET基因的开放阅读框的长度在597～918bp之间,编码多肽的长度在198～305aa之间。这些SWEET多肽的分子量和等电点分别为21998.71～33974.59Da和7.51～9.83。此外,除了CmSWEET12c在线粒体中、CmSWEET17a在细胞膜/高尔基体/过氧化物酶体中,其他CmSWEETs都定位在细胞膜上。由图1可知,共有21个CmSWEET基因位于21条染色体中的13条染色体上。染色体Cmo15含有最多的CmSWEET基因,而染色体Cmo00、Cmo01、Cmo03、Cmo05、Cmo07、Cmo09、Cmo19和Cmo20上没有CmSWEET基因。染色体Cmo10上含有3个CmSWEET基因,染色体Cmo11、Cmo12和Cmo13上都含有2个CmSWEET基因,其余每一条染色体含有1个CmSWEET基因。由表2可知,基因组鉴定全基因组鉴定共线性分析发现共有4对同源基因对,分别是CmSWEET1a_CmSWEET1b、CmSWEET12a_CmSWEET12d、CmSWEET10a_CmSWEET10b、CmSWEET7b_CmSWEET7c,其中依据氨基酸序列一致性>80%和基因比对覆盖率>0.75,本研究在中国桂花树染色体上共发现1对片段复制基因对,Ka/Ks<1.0表明这些复制对主要经历了纯化选择,鉴定分析分化时间为20.13MYA。
　　

图2结果显示

　　图2结果显示,17个CmSWEET蛋白有7个TMs,1个CmSWEET蛋白有5个TMs,其余的3个CmSWEET蛋白有6个TMs。为更好地了解桂花树SWEET基因的进化起源和功能,利用ClustalOmega和MEGA5.10软件,通过对21个CmSWEET蛋白序列和17个AtSWEET蛋白序列进行比对,构建了无根系统发育树。这些SWEET蛋白序列被分为4个亚族(图。图3比较了中国桂花树和拟南芥四个不同分支中的SWEET蛋白数量。每个分支都包含CmSWEET和AtSWEET。详细地说,第I亚族包含3个CmSWEETs和3个AtSWEETs;第亚族包含5个CmSWEETs和5个AtSWEETs;第亚族包含9个CmSWEETs和7个AtSWEETs;第IV亚族包含3个CmSWEETs和2个AtSWEETs(图。在构建的系统发育树的基础上,桂花树鉴定利用GSDS工具分析了21个CmSWEET基因的外显子-内含子结构、内含子阶段和保守基序(图。基序分析发现基序1、2、3和5在所有的中国桂花树SWEET蛋白中都存在(图4:A,B;图;基序4只在CmSWEET17b中不存在,基序7除在CmSWEET12a、CmSWEET17b蛋白中缺失,在其他蛋白中都存在;基序6只在第、第IV亚族中存在;基序8、9、和10只存在个别蛋白中(图4:B。外显子-内含子分类模式与系统发育树一致(图4:C。在中国桂花树SWEET基因家族中,所有CmSWEET基因的基因组序列都含有内含子。CmSWEET12a-12d、15a-15b、10a-10b、9、17a、17c、1a-1b、3和5a包含6个外显子,除CmSWEET15a外,它们的内含子阶段模式相似(图4:C。CmSWEET12e、7a-7c都包含5个外显子,除CmSWEET12e外,它们的内含子阶段模式相似。此外,CmSWEET17b和CmSWEET5b分别含有3个和4个外显子,且具有不同的内含子阶段模式(图4:C。为了解CmSWEETs基因的转录调控,我们从翻译起始位点提取了21个基因的2.0kb上游序列。在PlantCARE服务器上对21个基因上游序列中的顺式元件进行分析。图6:A为顺式作用元件在启动子上的位置。值得注意的是,在CmSWEETs的启动子区共鉴定出9种植物激素响应顺式元件,包括1种脱落酸响应元件、2种茉莉酸甲酯响应元件、3种赤霉素响应元件、1个水杨酸响应元件和2个生长素响应元件(图6:B。21个CmSWEET基因的启动子至少含有一个植物激素响应元件。CmSWEET7b含有最多植物激素响应元件(图6:B。CmSWEET基因除了受植物激素响应元件的调控外,还可能对各种环境胁迫,鉴定进化如厌氧、高温、低温、干旱和昼夜节律(图6:B。
　　

厌氧响应元件除在CmSWEET12d中不存在

　　厌氧响应元件除在CmSWEET12d中不存在,存在所有的CsSWEET启动子中(图,提示这些基因在厌氧诱导中的重要作用。然而,昼夜节律仅在CmSWEET12e、CmSWEET10a和CmSWEET17a中存在(图6:B。这些数据表明CmSWEET可能通过一个复杂的机制参与环境胁迫,并且每个CmSWEET基因都可以被不同的环境胁迫所诱导。根据模式物种拟南芥中的同源序列,构建了中国桂花树SWEET蛋白相互作用关系网络。由图7可知,CmSWEET5a对应AtSWEET5;CmSWEET9对应AtSWEET9;CmSWEET1a对应AtSWEET1;CmSWEET12a对应AtSWEET12;CmSWEET17a对应AtSWEET17;CmSWEET7a对应AtSWEET7;CmSWEET3对应AtSWEET3。图7结果显示,CmSWEET9分别与CmSWEET5a和CmSWEET1a;CmSWEET17a分别与CmSWEET12a、CmSWEET3和CmSWEET7a存在互作关系。
　　

采用与中国桂花树相同的方法

　　采用与中国桂花树相同的方法,桂花树苗从印度桂花树基因组中共鉴定到19个CmaSWEET基因(表。根据拟南芥的同源基因命名了19个CmaSWEET基因。无根进化树显示19个CmaSWEET被分为4个大分支,与前面中国桂花树SWEET的分类一致(图。其中,在CmaSWEET中,中国桂花树SWEET家族中第IV亚族成员又被分为3个小分支。每个大分支都包含中国桂花树SWEET蛋白和印度桂花树SWEET蛋白。从综合来看,家族基因组鉴定蛋白家族基因组鉴定两个栽培种SWEET蛋白系统进化相对保守。利用TBtools工具获取CmSWEET基因在CmaSWEET和AtSWEET基因中的共线性基因,并对这些共线性关系进行了统计。图9结果表明,中国桂花树与印度桂花树SWEET的共线性基因对数目大于中国桂花树与拟南芥的共线性基因对数目。此外,统计中国桂花树SWEET基因在印度桂花树和拟南芥中的拷贝数的结果(图9,表发现,CmSWEET10a和CmSWEET10b对应印度桂花树1个SWEET基因;AtSWEET1对应CmSWEET1a和CmSWEET1b;AtSWEET10对应CmSWEET10a和CmSWEET12c两个基因;AtSWEET5对应CmSWEET15b和CmSWEET17a两个基因。因此,我们认为中国桂花树SWEET基因与印度桂花树的进化关系相对拟南芥进化关系更加紧密,中国桂花树SWEET基因在葫芦科的多倍化事件中得到了扩张。SWEET基因家族在植物界广泛分布,并参与许多重要的生理和生化过程。然而,仅在基因组水平上鉴定了SWEET基因家族的少数物种,如拟南芥、水稻、番茄和大豆。据我们所知,迄今为止尚未有关葫芦科桂花树属作物的SWEET基因家族的系统研究报告。本研究中,我们通过分析中国桂花树的系统发育、染色体分布、基因结构、保守基序、顺式作用调控元件等,探索了中国桂花树中的21个SWEET基因。基因复制,包括片段复制、串联复制和全基因组复制,他们在生物的进化过程中起重要作用。
　　

中国桂花树在经历了核心真双子叶植物共同的六倍体事件后

　　中国桂花树在经历了核心真双子叶植物共同的六倍体事件后,又经历了葫芦科共有的全基因组加倍事件的证据,而后桂花树属共享一次最近的全基因组加倍事件CST(于继高等,202。在中国桂花树中发现4对SWEET共线基因对,有一对复制基因,由于这两个基因不在同一条染色体上,因此我们认为是片段复制基因对。由Ka/Ks<1可知,中国桂花树SWEET基因正在进行纯化正向选择。通过对中国桂花树和拟南芥SWEET蛋白进行系统发育树分析,发现他们被分配到4个亚族,蛋白家族基因组鉴定SWEET蛋白家族基因组鉴定这一分类与先前拟南芥的分类一致。将中国桂花树和印度桂花树SWEET蛋白进行系统发育树分析,发现每一个大分支都含有中国桂花树和印度桂花树SWEET蛋白,这表明中国桂花树和印度桂花树在进化上可能来自同一个祖先。通过SWEET基因结构分析,发现同一亚族中的大多数基因在或外显子数目或内含子数量方面具有相似的结构特征,这一现象与先前其他物种中SWEET家族结构特征相似。
　　

同一个亚族中的CmSWEET17a和CmSWEET17b

　　然而,同一个亚族中的CmSWEET17a和CmSWEET17b,与同一亚族中的其他CmSWEETs相比,它们含有不同数量的外显子,这说明中国桂花树SWEET基因家族的结构多样性。就特定的序列基序而言,在不同的亚族之间也检测到高度差异。然而,在一个共同的亚族中,大多数中国桂花树SWEET蛋白具有保守的基序,这表明同一亚族可能具有相似的功能。本研究发现,21个CmSWEETs基因启动子序列中的植物激素顺式元件和胁迫响应相关。在所有的CmSWEETs基因启动子中,至少发现了2个与植物激素响应有关的顺式作用元件,这意味着相应的激素可能在其调节中起着重要作用。此外,所有CmSWEETs基因至少包含一个应激相关的顺式作用元件,表明这些基因在应激响应中也起着重要的调控作用。
　　

还需要进一步地研究来确定这些顺式元件是否以及如何在中国桂花树中起作用

　　然而,还需要进一步地研究来确定这些顺式元件是否以及如何在中国桂花树中起作用。在中国桂花树中发现21个CmSWEETs,在印度桂花树中发现19个CmaSWEETs。Huetal.在黄瓜中发现17个CsSWEETs。中国桂花树、印度桂花树与黄瓜都是密切相关的,虽然黄瓜有7条染色体,中国桂花树和印度桂花树都有21条。详细的共线基因对和进化关系分析揭示了葫芦科三种重要作物染色体进化和重排的高度复杂性。因此,黄瓜、中国桂花树和印度桂花树中SWEET基因的相似性水平特别高,这些同源的SWEET基因可能具有相似的功能。