为研究L.plantarum ST与其余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.plantarum ST继续进行后续的比较基因组学分析。选取菌株为：JDM1、LIP-1、P8、ZJ316、ZS2058、WCSF1、ST-III、ATCC14197T、KP。

10 株L.plantarum共注释到25 个功能大类。由图3可知，10株L.plantarum中参与碳水化合物代谢的基因数量最多（21.6%），其次是氨基酸及其衍生物合成与降解（14.62%）和蛋白质代谢（11.87%）相关基因，其中氮代谢（0.21%）相关功能编码基因仅在菌株JDM1、WCFS1、ZS2058中存在。L.plantarum ST注释到24个功能类别。

注释结果表明这些L.plantarum参与碳水化合物代谢的基因数量最多。碳水化合物代谢能力与基因组中存在多种编码、运输糖代谢的基因及相应的调控蛋白有关。因为L.plantarum可以广泛利用多种糖源，所以可以应对更多复杂的环境。很多磷酸转移酶系统（PTS）转运系统存在于L.plantarum中，为机体输送糖类物质是该系统的作用。RAST注释结果显示WCFS1中共注释到57个PTS相关功能基因，ST共注释到49个。已有研究报道，L.plantarum WCFS1株包括完整的PTS酶II复合物以及运输系统，与WCFS1相比，本研究发现ST也有相同的糖类物质运送系统，从而确保菌株对糖的利用。

为了从基因组层面探究L.plantarum ST对碳水化合物的利用能力，利用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.plantarum ST具有良好的降解纤维素潜力。由图可知注释到CBM32、CBM34、CBM48，其中CBM48主要来自糖苷水解酶，多为普鲁兰酶、异淀粉酶、分支酶等。CBM的多样性使其在底物特异识别、促进酶与底物结合及酶稳定性等方面均有重要作用。通过碳水化合物注释结果推测L.plantarum ST有着较强的水解或重排糖苷键能力，同时可能具有良好的降解纤维素潜力。良好糖苷水解能力为L.plantarum ST用于益生菌制品奠定了基础。

将L.plantarumST与其他9株L.plantarum的基因组经毒力因子、耐药基因数据库比对后发现均不存在毒力因子和耐药基因，故认为L.plantarum ST是安全的，并且没有潜在的耐药性。

为进一步挖掘L.plantarum ST功能基因同时解析其潜在的益生功能基因，继续对菌株的特有基因和益生基因进行分析。L.plantarum ST存在41个特有基因，除去一些假定蛋白编码基因后，L.plantarum ST有3个特有功能基因，表2为特有功能基因信息，分别为乙二醛还原酶基因（yvgN）、酪氨酸重组酶基因（xerC）、能量耦合因子转运蛋白基因（ecfT）。yvgN是一种加氢酶，对乙二醛进行加氢反应。ecfT是一种转运酶，负责摄入一些维生素及其他微量元素。

基于前期对L.plantarum ST益生功能的研究，使用Prokka软件注释到L.plantarum ST含有与益生特性相关的基因，具体信息见表3。群体感应AI-2信号分子合成的最关键蛋白酶是S-核糖同型半胱氨酸裂解酶（LuxS），AI-2在微生物种内和种间交流起到重要作用，在L.plantarum ST中发现了该基因。目前发现在超过80种细菌中均存在LuxS基因，细菌许多生理现象离不开LuxS介导的群体感应系统，因为这种系统对细菌有重要的调节作用。L.plantarum ST还含有谷胱甘肽合成（garB、gshAB）、核黄素合成（ribBA、ribE、ribF）相关基因，核黄素（VB2），作为辅酶有助于代谢是其主要功能。NAPPH还原酶类（nfrA），参与脂肪酸、核苷酸和脂类的合成等多种合成代谢反应，该物质起到递氢体的作用在生物体内的化学反应中。L.plantarum ST还注释到了黏附因子（atpC）基因，黏附因子有助于细菌在宿主体内的定植因为可使菌株牢固黏附于宿主细胞表面。研究还发现，L.plantarum ST存在提高宿主代谢能力（tagE）、乳酸和转运乳酸（pyk、ldh）、细胞壁特殊组分成分（ItaS）和生物活性肽（lmrA）产生相关的基因，可提高菌株益生功能。特有基因和益生基因分析发现，L.plantarum ST含有与黏附能力、信号分子传递和抗氧化能力等相关的益生功能基因。因此，本研究从基因组水平揭示L.plantarum ST具有益生特性相关基因，认为其是1株具有潜在益生功能的菌株。

8 碳水化合物代谢特性

对L.plantarum ST结合NCBI已公开的152株L.plantarum全基因组序列和1株模式菌株ATCC14197T的基因组序列进行比较基因组学分析。ANI分析鉴定ST为L.plantarum，系统发育树分析发现L.plantarum ST与肠道分离株BCC9546遗传距离较近，乳、肉制品分离株因来源不同而存在差异，分别集中在第II分支的上半部分和下半部分，同时果蝇分离株有明显聚集趋势。RAST和KEGG功能注释发现，L.plantarum ST参与碳水化合物代谢的基因数量最多，包含大量的PTS相关基因，还发现PPAR和RIG-1免疫调控通路相关基因。CAZy注释发现，L.plantarum ST有较多水解或重排糖苷键的基因，可能具有良好的纤维素降解能力。毒力和耐药注释比对后发现ST不存在毒力因子和耐药基因。特有基因分析发现，L.plantarum ST携带特有的与能量转运功能有关功能基因ecfT；益生特性基因显示，L.plantarum ST具有与群体感应信号分子AI-2（LuxS）、黏附分子（atpC）和谷胱甘肽（garB、gshAB）合成相关的功能基因。API50CHL碳水化合物代谢结果显示，L.plantarum ST可代谢单糖类、糖苷和二糖类及多糖类等中的29种碳源，ST能够利用的碳源非常丰富。综上，本研究从基因组水平分析了L.plantarum ST功能基因组特征，结合糖发酵表型结果，认为L.plantarum ST是1株安全且具有潜在益生特性基因的益生菌，为后续L.plantarum ST益生功能开发及其生产应用奠定了遗传学基础。