植物乳植杆菌（Lactiplantibacillusplantarum）属于乳植杆菌属，其广泛存在于食品（如蔬菜、乳肉制品等）和肠胃道等多种环境中。L.plantarum具备改善肠道菌群平衡、调节免疫代谢、抵抗病原体、降低胆固醇、缓解乳糖不耐症和加强营养成分吸收等益生功能，是益生菌开发的潜力贮备。高通量测序技术的快速发展导致对益生菌的研究也上升到基因组学层面，越来越多的菌株完成了全基因组测序，因此，利用比较基因组学对L.plantarum功能基因分析，对其功能开发与应用具有重要意义。
内蒙古农业大学的杨淑娟和云南农业大学食品科学技术学院的周金萍、林秋叶*等人对L.plantarumST及NCBI已公开的152株L.plantarum全基因组序列和1株模式菌株ATCC14197T基因组序列进行比较基因组学分析，从基因组水平全面分析该菌株的功能基因，同时利用API50CHL测定该菌株的碳水化合物利用情况，旨在为使其成为功能性食品的生产菌株或微生态制剂奠定理论基础。
1ANI分析与泛-核心基因构建
ANI是通过比对基因组的同源序列以鉴定菌株亲缘关系。在比较基因组学分析过程中通常认为ANI值大于95%为同一物种，ANI可用于判断基因组间的相似性及评估基因组间多态性程度。本研究对154株L.plantarum菌株进行ANI计算，并构建聚类热图（图1A）。结果表明，全部菌株之间ANI值均大于98%，表明研究所用菌株均为同一物种即L.plantarum。其中，L.plantarumST和模式菌株L.plantarumATCC14917T的ANI为99.05%。
154株L.plantarum的泛基因集包括了17378个基因，其中核心基因1262个，特有基因6427个。由图1B可知，随基因组数量的增加，泛基因个数依然表现出增加趋势，而核心基因个数逐渐趋于稳定。这一现象符合开放式基因组的特点，所以推测L.plantarum为开放式基因组。
2 系统发育树分析
为研究L.plantarumST与其余153株L.plantarum种内遗传进化关系，采用NJ法，以Roary软件识别得到的1262个核心基因构建系统发育树（图2），Bootstrap值为1000。154株L.plantarum中包括76株蔬菜与谷物分离源、14株乳制品分离源、10株肉制品分离源、22株肠道分离源、4株果蝇分离源及28株其他分离源。系统发育树结果显示，154株L.plantarum在进化过程中由于遗传多样性使其被分为两大分支，集中在1个分支说明亲缘关系更接近。ST处于第II分支中的单独一小分支，与肠道分离株BCC9546遗传距离较近。从系统发育树中分支中选取几株代表性、已有相关益生特性研究报道的菌株，作为对照，与L.plantarumST继续进行后续的比较基因组学分析。选取菌株为：JDM1、LIP-1、P8、ZJ316、ZS2058、WCSF1、ST-III、ATCC14197T、KP。
3RAST注释结果
10 株L.plantarum共注释到25 个功能大类。由图3可知，10株L.plantarum中参与碳水化合物代谢的基因数量最多（21.6%），其次是氨基酸及其衍生物合成与降解（14.62%）和蛋白质代谢（11.87%）相关基因，其中氮代谢（0.21%）相关功能编码基因仅在菌株JDM1、WCFS1、ZS2058中存在。L.plantarumST注释到24个功能类别。
注释结果表明这些L.plantarum参与碳水化合物代谢的基因数量最多。碳水化合物代谢能力与基因组中存在多种编码、运输糖代谢的基因及相应的调控蛋白有关。因为L.plantarum可以广泛利用多种糖源，所以可以应对更多复杂的环境。很多磷酸转移酶系统（PTS）转运系统存在于L.plantarum中，为机体输送糖类物质是该系统的作用。RAST注释结果显示WCFS1中共注释到57个PTS相关功能基因，ST共注释到49个。已有研究报道，L.plantarumWCFS1株包括完整的PTS酶II复合物以及运输系统，与WCFS1相比，本研究发现ST也有相同的糖类物质运送系统，从而确保菌株对糖的利用。
4KEGG注释结果
10株L.plantarum在KEGG数据库中分别在新陈代谢、细胞过程、遗传信息处理、环境信息处理、生物系统和人类疾病6大功能37个通路上得到功能注释，结果如图4所示。分析发现所有菌株的基因功能注释主要分布于代谢通路中，碳水化合物代谢是11个代谢通路中基因数量最多的，共2008个基因占代谢通路注释基因的23.2%。氨基酸代谢次之。环境信息处理中最多注释的是信号转导和膜运输，其中信号转导通路主要集中于ABC转运蛋白。多种糖类的代谢通路包含在碳水化合物代谢通路中，L.plantarumST中碳水化合物代谢多达200个基因注释量，推测L.plantarumST可能有较强的糖类分解及代谢功能，为机体提供能量。
免疫相关基因的通路信息及其通路ID如表1所示。可调控免疫和炎症的通路相关基因，包括过氧化物酶体增生物激活受体（PPAR）和维甲酸诱导基因I（RIG-I）通路。其中，STGL000154、STGL002738参与PPAR信号通路的调控，PPAR通路与细胞分化、炎性反应和能量代谢等紧密相连。在RIG-I样受体信号通路的调控中基因STGL000333参与，RIG-I样受体是固有免疫的模式识别受体，这个受体可以辨别非自身的病毒RNA，通过激活该受体信号通路最终达到促进细胞因子生成的目的，从而发挥出抗病毒的功效，这可能是L.plantarum具有抗病毒作用的因素。推测L.plantarumST能够帮助宿主抵抗病原菌的侵袭、进而调节肠道菌群的平衡。
5CAZy注释结果
为了从基因组层面探究L.plantarumST对碳水化合物的利用能力，利用CAZy数据库分析了10株L.plantarum的CAZy功能基因。5种CAZy类在这10株L.plantarum中被注释到，分别为碳水化合物碳结合结构域（CBMs）、糖苷水解酶类（GHs）、辅模块酶类（AAs）、糖苷转移酶类（GTs）、碳水化合物酯酶类（CEs）。10株L.plantarum共注释到47个CAZy基因家族，包括24个GHs家族、12个GTs家族、4个CEs家族、4个CBMs家族和3个AAs家族。
由图5可知，ST及另外9株菌注释到的GHs、GTs含量最为丰富，无显著差异。GHs可水解或重排糖苷键，将多糖水解为单糖，给予菌体能量。GTs催化糖基与其他糖类或非糖物质的基团进行转移，负责糖苷键的形成。其中，GT4在L.plantarum中拷贝数最多（平均约13个），其次是GT2（平均约11个）和GH1（平均约10个）。本研究10株L.plantarum注释到了24个GHs，含量丰富。ST的GH家族中含量最多的是GH1、GH13_31，还包括GH170、GH25、GH65、GH73、GH78、GH36等。在纤维素降解的过程中糖苷水解酶类中一些酶类有水解作用，利用关键限速酶β-葡萄糖苷酶能把纤维二糖在内的纤维寡糖分解为单糖，故推测L.plantarumST具有良好的降解纤维素潜力。由图可知注释到CBM32、CBM34、CBM48，其中CBM48主要来自糖苷水解酶，多为普鲁兰酶、异淀粉酶、分支酶等。CBM的多样性使其在底物特异识别、促进酶与底物结合及酶稳定性等方面均有重要作用。通过碳水化合物注释结果推测L.plantarumST有着较强的水解或重排糖苷键能力，同时可能具有良好的降解纤维素潜力。良好糖苷水解能力为L.plantarumST用于益生菌制品奠定了基础。
6 毒力和耐药基因检索分析
将L.plantarumST与其他9株L.plantarum的基因组经毒力因子、耐药基因数据库比对后发现均不存在毒力因子和耐药基因，故认为L.plantarumST是安全的，并且没有潜在的耐药性。
7 特有基因和益生基因分析
为进一步挖掘L.plantarumST功能基因同时解析其潜在的益生功能基因，继续对菌株的特有基因和益生基因进行分析。L.plantarumST存在41个特有基因，除去一些假定蛋白编码基因后，L.plantarumST有3个特有功能基因，表2为特有功能基因信息，分别为乙二醛还原酶基因（yvgN）、酪氨酸重组酶基因（xerC）、能量耦合因子转运蛋白基因（ecfT）。yvgN是一种加氢酶，对乙二醛进行加氢反应。ecfT是一种转运酶，负责摄入一些维生素及其他微量元素。
基于前期对L.plantarumST益生功能的研究，使用Prokka软件注释到L.plantarumST含有与益生特性相关的基因，具体信息见表3。群体感应AI-2信号分子合成的最关键蛋白酶是S-核糖同型半胱氨酸裂解酶（LuxS），AI-2在微生物种内和种间交流起到重要作用，在L.plantarumST中发现了该基因。目前发现在超过80种细菌中均存在LuxS基因，细菌许多生理现象离不开LuxS介导的群体感应系统，因为这种系统对细菌有重要的调节作用。L.plantarumST还含有谷胱甘肽合成（garB、gshAB）、核黄素合成（ribBA、ribE、ribF）相关基因，核黄素（VB2），作为辅酶有助于代谢是其主要功能。NAPPH还原酶类（nfrA），参与脂肪酸、核苷酸和脂类的合成等多种合成代谢反应，该物质起到递氢体的作用在生物体内的化学反应中。L.plantarumST还注释到了黏附因子（atpC）基因，黏附因子有助于细菌在宿主体内的定植因为可使菌株牢固黏附于宿主细胞表面。研究还发现，L.plantarumST存在提高宿主代谢能力（tagE）、乳酸和转运乳酸（pyk、ldh）、细胞壁特殊组分成分（ItaS）和生物活性肽（lmrA）产生相关的基因，可提高菌株益生功能。特有基因和益生基因分析发现，L.plantarumST含有与黏附能力、信号分子传递和抗氧化能力等相关的益生功能基因。因此，本研究从基因组水平揭示L.plantarumST具有益生特性相关基因，认为其是1株具有潜在益生功能的菌株。
8 碳水化合物代谢特性
为了评估L.plantarumST对不同碳源的利用能力，利用API 50 CHL碳水化合物鉴定试剂条对ST的碳水化合物代谢进行研究，挖掘出可被ST利用的碳源，为碳水化合物代谢中的功能基因分析提供表型数据支持，同时为该菌株在之后的应用提供借鉴。
如表4所示，API50CHL碳水化合物鉴定试剂条所涉及的49种碳水化合物中，可被ST利用的共有28种，划分为单糖类（9）、糖苷类和二糖类（12）、多糖类（1）、醇类（4）、盐类（2）物质五大类。单糖类包括：L-阿拉伯糖、D-核糖、D-半乳糖、D-葡萄糖、D-甘露糖、L-鼠李糖、甲基-α-D-吡喃葡萄糖苷、N-乙酰葡萄糖胺；糖苷类和二糖类包括：熊果苷、水杨苷、苦杏仁苷、七叶灵、D-纤维二糖、D-麦芽糖、D-乳糖、D-蜜二糖、D-蔗糖、D-海藻糖、D-龙胆糖、D-松二糖；多糖类包括：D-松三糖；醇类包括：丙三醇、肌醇、甘露醇、山梨醇；盐类包括：葡萄糖酸钾、5-酮基葡萄糖酸钾。糖发酵实验结果显示，L.plantarumST可代谢不同种类丰富的碳源，包括9种单糖、12种糖苷和二糖、1种多糖及4种醇类和2种盐类物质，再次说明了ST具有较强的碳水化合物代谢能力。该结果与本实验KEGG注释以及CAZy分析吻合，利用基因组学和表型分析相结合，发现L.plantarumST拥有较好的产CAZy能力，具备丰富的碳水化合物利用能力，可为开发相关酶制剂提供一定的参考价值。
结 语
对L.plantarumST结合NCBI已公开的152株L.plantarum全基因组序列和1株模式菌株ATCC14197T的基因组序列进行比较基因组学分析。ANI分析鉴定ST为L.plantarum，系统发育树分析发现L.plantarumST与肠道分离株BCC9546遗传距离较近，乳、肉制品分离株因来源不同而存在差异，分别集中在第II分支的上半部分和下半部分，同时果蝇分离株有明显聚集趋势。RAST和KEGG功能注释发现，L.plantarumST参与碳水化合物代谢的基因数量最多，包含大量的PTS相关基因，还发现PPAR和RIG-1免疫调控通路相关基因。CAZy注释发现，L.plantarumST有较多水解或重排糖苷键的基因，可能具有良好的纤维素降解能力。毒力和耐药注释比对后发现ST不存在毒力因子和耐药基因。特有基因分析发现，L.plantarumST携带特有的与能量转运功能有关功能基因ecfT；益生特性基因显示，L.plantarumST具有与群体感应信号分子AI-2（LuxS）、黏附分子（atpC）和谷胱甘肽（garB、gshAB）合成相关的功能基因。API50CHL碳水化合物代谢结果显示，L.plantarumST可代谢单糖类、糖苷和二糖类及多糖类等中的29种碳源，ST能够利用的碳源非常丰富。综上，本研究从基因组水平分析了L.plantarumST功能基因组特征，结合糖发酵表型结果，认为L.plantarumST是1株安全且具有潜在益生特性基因的益生菌，为后续L.plantarumST益生功能开发及其生产应用奠定了遗传学基础。
本文《基于比较基因组学解析植物乳植杆菌ST的功能基因组》来源于《食品科学》2023年44卷第14期86-93页，作者：杨淑娟，周金萍，李海燕，曹振辉，孙志宏，林秋叶。DOI:10.7506/spkx1002-6630-20220807-086。
