基于比较基因组学解析植物乳植杆菌ST的功能基因组.pdf

86 2023,Vol.44,No.14 食品科学 生物工程基于比较基因组学解析植物乳植杆菌ST的 功能基因组杨淑娟1，周金萍2，李海燕2，曹振辉3，孙志宏1，林秋叶2,*（1.内蒙古农业大学 乳品生物技术与工程教育部重点实验室，农业农村部奶制品加工重点实验室，内蒙古自治区乳品生物技术与工程重点实验室，乳酸菌与发酵乳制品省部共建协同创新中心，内蒙古 呼和浩特 010018；2.云南农业大学食品科学技术学院，云南 昆明 650201；3.云南农业大学动物科学技术学院，云南 昆明 650201）摘 要：解析植物乳植杆菌（Lactiplantibacillus plantarum）ST的功能基因组特征，为其开发利用奠定基础。对L.plantarum ST与1 株模式菌株ATCC14197T基因组序列及NCBI中已公开的152 株L.plantarum全基因组序列进行比较基因组学分析。基于154 株L.plantarum的1 262 个核心基因构建了系统发育树，发现ST与肠道分离株BCC9546遗传距离最近。乳、肉制品分离株分别集中在第II分支的上半部分和下半部分，因来源不同而存在差异，且存在一定的聚集性。同时果蝇分离株具有明显聚集趋势。功能注释分析发现，ST基因组中包含磷酸转移酶系统糖类转运系统相关基因，还注释到了免疫调控通路相关基因，同时有较多水解或重排糖苷键的基因。ST不存在毒力因子和耐药基因。与其他L.plantarum相比，ST携带特有的与能量转运功能有关的功能基因ecfT。此外，ST具有与群体感应信号分子AI-2、黏附分子和谷胱甘肽合成相关的功能基因。API50CHL碳水化合物代谢结果显示，L.plantarum ST能利用丰富的碳源，可代谢包括单糖类、糖苷和二糖类及多糖类等中29 种碳源。本研究从基因组层面分析L.plantarum ST功能基因组特征，同时结合糖发酵表型实验结果，认为L.plantarum ST是1 株安全且具有潜在益生特性的益生菌，为其开发利用提供遗传学基础。关键词：植物乳植杆菌ST；比较基因组学；功能基因Comparative Genomic Analysis of Functional Genomics of Lactiplantibacillus plantarum ST YANG Shujuan1,ZHOU Jinping2,LI Haiyan2,CAO Zhenhui3,SUN Zhihong1,LIN Qiuye2,*(1.Key Laboratory of Dairy Biotechnology and Engineering,Ministry of Education,Key Laboratory of Dairy Products Processing,Ministry of Agriculture and Rural Affairs,Inner Mongolia Key Laboratory of Dairy Biotechnology and Engineering,Collaborative Innovative Center of Ministry of Education for Lactic Acid Bacteria and Fermented Dairy Products,Inner Mongolia Agricultural University,Hohhot 010018,China;2.College of Food Science and Technology,Yunnan Agricultural University,Kunming 650201,China;3.Faculty of Animal Science and Technology,Yunnan Agricultural University,Kunming 650201,China)Abstract:In this study,the functional genomic characteristics of Lactiplantibilus plantarum ST were analyzed.Comparative genomic analysis was performed on the whole genome sequences of L.plantarum ST and the type strain L.plantarum ATCC 14197T as well as those of 152 L.plantarum strains published in the NCBI GenBank.A phylogenetic tree was constructed based on 1 262 core genes of these 154 L.plantarum strains,and the genetic distance between ST and the intestinal isolate BCC9546 was found to be the closest.The isolates originated from milk and meat products were concentrated in the upper part and the lower part of the second branch,respectively,and they were different due to different sources and could be aggregated.Meanwhile,Drosophila isolates showed a significant aggregation trend.Functional annotation analysis suggested that the ST genome contained genes related to the phosphotransferase system(PTS)as a major carbohydrate transport system.Genes related to the immune regulatory pathway were also annotated.There were many 收稿日期：2022-08-07基金项目：国家自然科学基金地区科学基金项目（31760448）；云南省孙志宏专家工作站（202105AF150055）；“兴滇英才支持计划”青年人才项目（YNWR-QNBJ-2018-137）第一作者简介：杨淑娟（1997）（ORCID:0000-0002-8494-1128），女，博士研究生，研究方向为乳酸菌资源开发与利用。E-mail:*通信作者简介：林秋叶（1979）（ORCID:0000-0002-8287-0047），女，教授，博士，研究方向为民族发酵食品乳酸菌资源挖掘及利用。E-mail:生物工程 食品科学 2023,Vol.44,No.14 87genes related to the hydrolysis or rearrangement of glycosidic bonds,but no virulence factors or antibiotics resistance genes found in ST.Compared with the other L.plantarum isolates,ST carried the unique functional gene ecfT related to energy transport function.In addition,ST had functional genes involved in the synthesis of the quorum-sensing signaling molecule AI-2,glutathione,and adhesion molecules.The results of API 50 CHL carbohydrate metabolism showed that L.plantarum ST could utilize a wide range of carbon sources and metabolize 29 carbon sources including monosaccharides,glycosides,disaccharides and polysaccharides.This study provides a genetic basis for the development and utilization of L.plantarum ST as a safe probiotic strain.Keywords:Lactiplantibacillus plantarumST;comparativegenomics;functionalgenesDOI:10.7506/spkx1002-6630-20220807-086中图分类号：TS201.3 文献标志码：A 文章编号：1002-6630（2023）14-0086-08引文格式：杨淑娟,周金萍,李海燕,等.基于比较基因组学解析植物乳植杆菌ST的功能基因组J.食品科学,2023,44(14):86-93.DOI:10.7506/spkx1002-6630-20220807-086.http:/YANG Shujuan,ZHOU Jinping,LI Haiyan,et al.Comparative genomic analysis of functional genomics of Lactiplantibacillus plantarum STJ.Food Science,2023,44(14):86-93.(in Chinese with English abstract)DOI:10.7506/spkx1002-6630-20220807-086.http:/植物乳植杆菌（Lactiplantibacillus plantarum）属于乳植杆菌属，是乳酸菌的一种，其广泛存在于食品（如蔬菜、乳肉制品等）和肠胃道等多种环境中1，己被人类食用或使用上千年，是一类公认安全的微生物2。L.plantarum具备改善肠道菌群平衡、调节免疫代谢、抵抗病原体、降低胆固醇、缓解乳糖不耐症和加强营养成分吸收等益生功能，是益生菌开发的潜力贮备3-4。因此，L.plantarum作为潜在益生菌，在食品工业和医疗保健等领域得到广泛开发与应用。高通量测序技术的快速发展导致对益生菌的研究也上升到基因组学层面，越来越多的菌株完成了全基因组测序5-6，这种方法在基因组层面上提供了对菌株功能的新认知。比较基因组学能够深入探究菌株的基因组及功能基因7-8。MaoBingyong等9利用比较基因组学分析不同来源的L.plantarum基因组差异，发现L.plantarum特异基因数量与其所处的生态位有关。此外，还有研究发现L.plantarum ZJ316具有胆盐抗性和黏附于宿主肠壁的相关基因，可在胃肠道存活并定植，维持肠道菌群平衡10。L.plantarum DMDL9010存在降解亚硝酸盐的相关基因，为降低食品中的亚硝酸盐含量提供理论参考11。因此，利用比较基因组学对L.plantarum功能基因分析，对其功能开发与应用具有重要意义。本研究所用L.plantarum ST分离自中国云南省德昂酸茶，已完成对L.plantarum ST的全基因组测序和基因组组装工作12。前期研究发现L.plantarum ST对鼠伤寒沙门氏菌和大肠杆菌表现出较强的抗菌活性，在肠道中具有良好的胃肠液耐受能力、较强的黏附和自聚力，是1 株潜在的益生菌株13。本研究对L.plantarum ST及NCBI已公开的152 株L.plantarum全基因组序列和1 株模式菌株ATCC14197T基因组序列进行比较基因组学分析，从基因组水平全面分析该菌株的功能基因，同时利用API50CHL测定该菌株的碳水化合物利用情况，旨在为使其成为功能性食品的生产菌株或微生态制剂奠定理论基础。1 材料与方法1.1 材料与试剂所用菌株L.plantarum ST分离自中国云南德昂族传统食品德昂酸茶中，由云南农业大学提供，GenBank数据库生物保藏号为PRJNA792432。截至2021年10月10日，把NCBI（NationalCoalitionBuildingInstitute，https:/www.ncbi.nlm.nih.gov/）Refseq数据库中的152 株L.plantarum基因组完成图全部下载，同时下载模式菌株L.plantarum ATCC14917T基因组完成图。MRS培养基 广州环凯微生物科技有限公司；API50CHL碳水化合物鉴定试剂条 法国生物梅里埃生物公司。1.2 仪器与设备电热恒温培养箱 上海一恒科技有限公司；台式高速离心机 德国Eppendorf公司。1.3 方法1.3.1 比较基因组分析1.3.1.1 平均核苷酸一致性（averagenucleotideidentity，ANI）计算154 株L.plantarum的ANI值计算以Goris等14方法为参照。再利用TBtools14软件绘制ANI聚类热图。1.3.1.2 泛-核心基因集构建菌株基因组的基因预测采用Prokka15软件，之后利用Roary16软件识别核心基因集与泛基因集，编码蛋白氨88 2023,Vol.44,No.14 食品科学 生物工程基酸相似性大于95%是识别核心基因的原则8。泛-核心基因趋势图绘制使用软件R（v.4.0.4）进行。1.3.1.3 系统发育树构建核心基因序列是通过Roray软件分析后获得，再使用TreeBest软件（http:/ 基因功能注释利用子系统技术的快速标注（RapidAnnotationUsingSubsystemTechnology，RAST）（https:/rast.nmpdr.org/）数据库、京都基因与基因组百科全书（KyotoEncyclopediaofGenesandGenomes，KEGG）网站对具有代表性的、已有相关益生特性研究报道的L.plantarum全基因组功能注释。1.3.1.5 碳水化合物活性酶（carbohydrate-activeenzymes，CAZy）注释将L.plantarum基因组序列上传至在线注释dbCAN（http:/bcb.unl.edu/dbCAN2/）进行注释，下载其注释文件。统计L.plantarum基因组信息。1.3.2 毒力和耐药基因检索分析将L.plantarum氨基酸序列上传至毒力因子数据库（VirulenceFinder，https:/cge.cbs.dtu.dk/services/VirulenceFinder/）进行比对，标准为选择ID的阈值90%、最小长度60%。把L.plantarum氨基酸序列与耐药基因数据库（ComprehensiveAntibioticResearchDatabase，CARD）（http:/arpcard.mcmaster.ca）进行比对，耐药基因的确定为比对结果中E值l1010、相似度50%的蛋白序列。1.3.3 特有基因和益生基因分析利用Prokk软件对菌株基因组进行基因预测后，对菌株的特有基因及益生基因进行深入分析。1.3.4 碳水化合物代谢特性分析利用API 50CHL碳水化合物鉴定生化试剂条对L.plantarum ST碳水化合物代谢情况进行表型研究，并按照试剂盒说明书操作17。24、48 h各观察一次并记录结果，以48 h观察到结果作为最终结果。1.4 数据分析及绘图利用Perl脚本统计基因组信息。利用TBtools绘制ANI聚类热图。2 结果与分析2.1 ANI分析与泛-核心基因构建ANI是通过比对基因组的同源序列以鉴定菌株亲缘关系18。在比较基因组学分析过程中通常认为ANI值大于95%为同一物种，ANI可用于判断基因组间的相似性及评估基因组间多态性程度19。本研究对154 株L.plantarum菌株进行ANI计算，并构建聚类热图（图1A）。结果表明，全部菌株之间ANI值均大于98%，表明研究所用菌株均为同一物种即L.plantarum。其中，L.plantarum ST和模式菌株L.plantarumATCC14917T的ANI为99.05%.通过Prokka和Roary软件进行分析比较菌株间的基因差别，统计核、泛基因个数用大于90%的氨基酸一致性为标准。154 株L.plantarum的泛基因集包括了17 378 个基因，其中核心基因1 262 个，特有基因6 427 个。由图1B可知，随基因组数量的增加，泛基因个数依然表现出增加趋势，而核心基因个数逐渐趋于稳定。这一现象符合开放式基因组20的特点，所以推测L.plantarum为开放式基因组。100.0?99.598.597.597.096.595.596.098.099.0A0050200150100B15 00010 00020 0005 000?图 1 L.plantarum ANI值聚类热图（A）与泛基因集和核心 基因集曲线（B）Fig.1 Heatmap of ANI(A)and pan-core gene family curve(B)2.2 系统发育树分析菌株间遗传进化关系和群体结构可以通过系统发育树反映21。为研究L.plantarum ST与其余153 株L.plantarum种内遗传进化关系，采用NJ法，以Roary软件识别得到的1 262 个核心基因构建系统发育树（图2），Bootstrap值为1 000。154 株L.plantarum中包括76 株蔬菜与谷物分离源、14 株乳制品分离源、10 株肉制品分离源、22 株肠道分离源、4 株果蝇分离源及28 株其他分离源。系统发育树结果显示，154 株L.plantarum在进化过程中由于遗传多样性使其被分为两大分支，集中在1 个分支说明亲缘关系更接近。ST处于第II分支中的单独一小分支，与肠道分离株BCC9546遗传距离较近。乳制品分离株集中在第II分支生物工程 食品科学 2023,Vol.44,No.14 89的上半部分，肉制品分离株集中在下半部分，且存在一定的聚集性。可见，乳、肉制品分离株因来源不同而存在差异。此外，呈现明显聚集性有果蝇分离株，刘亚华等22探究发现果蝇分离株因为分离环境特别，它对环境的适应会获得有利于其在该环境下生存的环境特异性基因，从而导致其具有一定分离源特异性。从系统发育树中分支中选取几株代表性、已有相关益生特性研究报道的菌株，作为对照，与L.plantarum ST继续进行后续的比较基因组学分析。选取菌株为：JDM1、LIP-1、P8、ZJ316、ZS2058、WCSF1、ST-III、ATCC14197T、KP。0.01?III图 2 基于核心基因序列的系统发育树Fig.2 Phylogenetic tree based on core gene sequences2.3 RAST注释结果利用RAST软件对10 株具有良好益生特性的L.plantarum基因组序列进行功能注释，分析ST与其他菌株功能差异。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可以广泛利用多种糖源，所以可以应对更多复杂的环境23。很多磷酸转移酶系统（phosphotransferasesystem，PTS）转运系统存在于L.plantarum中，为机体输送糖类物质是该系统的作用24。RAST注释结果显示WCFS1中共注释到57 个PTS相关功能基因，ST共注释到49 个。已有研究报道，L.plantarum WCFS1株包括完整的PTS酶II复合物以及运输系统25，与WCFS1相比，本研究发现ST也有相同的糖类物质运送系统，从而确保菌株对糖的利用。9.10%5.19%3.27%21.60%14.62%0.37%0.71%0.64%1.94%1.41%0.54%0.21%3.01%5.51%0.36%1.48%0.37%11.87%7.96%3.40%0.46%3.11%1.05%1.26%0.55%AB14.57%0.29%0.67%21.38%9.68%2.97%0.48%1.34%0.58%0.48%3.64%8.34%12.94%0.38%1.44%0.38%3.07%5.56%0.00%0.58%1.53%1.92%0.77%3.36%3.64%?A.10 株L.plantarum；B.L.plantarum ST。图 3 L.plantarum RAST功能注释结果Fig.3 Function annotation of L.plantarum strains based on the RAST database2.4KEGG注释结果L.plantarum ST与另外9 株L.plantarum的氨基酸序列，上传至KEGG数据库进行比对得到注释结果。10 株L.plantarum在KEGG数据库中分别在新陈代谢、细胞过程、遗传信息处理、环境信息处理、生物系统和人类疾病6 大功能37 个通路上得到功能注释，结果如图4所示。分析发现所有菌株的基因功能注释主要分布于代谢通路中，碳水化合物代谢是11 个代谢通路中基因数量最多的，共2 008 个基因占代谢通路注释基因的23.2%。氨基酸代谢次之。环境信息处理中最多注释的是信号转导和膜运输，其中信号转导通路主要集中于ABC转运蛋白。多种糖类的代谢通路包含在碳水化合物代谢通路中，L.plantarum ST中碳水化合物代谢多达200 个基因注释量，推测L.plantarum ST可能有较强的糖类分解及代谢功能，为机体提供能量6。90 2023,Vol.44,No.14 食品科学 生物工程STJDM114917ST-3KPP8ZJ316ZS2058LPL-1WCSF1050100150200250?-?图 4 L.plantarum KEGG数据库注释结果Fig.4 KEGG annotation of L.plantarum免疫相关基因的通路信息及其通路ID如表1所示。可调控免疫和炎症的通路相关基因，包括过氧化物酶体增生物激活受体（peroxisomeproliferatoractivatedreceptor，PPAR）和维甲酸诱导基因I（retinoicacid-inducedgeneI，RIG-I）通路。其中，STGL000154、STGL002738参与PPAR信号通路的调控，PPAR通路与细胞分化、炎性反应和能量代谢等紧密相连。在RIG-I样受体信号通路的调控中基因STGL000333参与，RIG-I样受体是固有免疫的模式识别受体，这个受体可以辨别非自身的病毒RNA，通过激活该受体信号通路最终达到促进细胞因子生成的目的，从而发挥出抗病毒的功效26，这可能是L.plantarum具有抗病毒作用的因素。推测L.plantarum ST能够帮助宿主抵抗病原菌的侵袭、进而调节肠道菌群的平衡。表 1 L.plantarum ST基因组免疫调控通路相关基因Table 1 Genes related to genomic immune regulatory pathway in L.plantarum ST 通路ID描述基因编号ko04622RIG-I样受体信号通路STGL000333ko03320PPAR信号通路STGL000154；STGL0027382.5CAZy注释结果为了从基因组层面探究L.plantarum ST对碳水化合物的利用能力，利用CAZy数据库分析了10 株L.plantarum的CAZy功能基因。5 种CAZy类在这10 株L.plantarum中被注释到，分别为碳水化合物碳结合结构域（carbohydrate-bindingmodules，CBMs）、糖苷水解酶类（glycosidehydrolases，GHs）、辅模块酶类（auxiliaryactivities，AAs）、糖苷转移酶类（glycosyltransferases，GTs）、碳水化合物酯酶类（carbohydrateesterases，CEs）。10 株L.plantarum共注释到47 个CAZy基因家族，包括24 个GHs家族、12 个GTs家族、4 个CEs家族、4 个CBMs家族和3 个AAs家族。由图5可知，ST及另外9 株菌注释到的GHs、GTs含量最为丰富，无显著差异。GHs可水解或重排糖苷键，将多糖水解为单糖，给予菌体能量27。GTs催化糖基与其他糖类或非糖物质的基团进行转移，负责糖苷键的形成28。其中，GT4在L.plantarum中拷贝数最多（平均约13 个），其次是GT2（平均约11 个）和GH1（平均约10 个）。本研究10 株L.plantarum注释到了24 个GHs，含量丰富。ST的GH家族中含量最多的是GH1、GH13_31，还包括GH170、GH25、GH65、GH73、GH78、GH36等。在纤维素降解的过程中糖苷水解酶类中一些酶类有水解作用，利用关键限速酶-葡萄糖苷酶能把纤维二糖在内的纤维寡糖分解为单糖4，故推测L.plantarum ST具有良好的降解纤维素潜力。所有生物工程 食品科学 2023,Vol.44,No.14 91L.plantarum GT家族中GT2、GT4含量最多，主要编码蔗糖合酶、纤维素合酶等众多酶类，与蔗糖等二糖以及脂多糖、纤维素等合成有关。CBM家族能够促进催化活性结构域与碳水化合物的结合。由图可知注释到CBM32、CBM34、CBM48，其中CBM48主要来自糖苷水解酶，多为普鲁兰酶、异淀粉酶、分支酶等。CBM的多样性使其在底物特异识别、促进酶与底物结合及酶稳定性等方面均有重要作用29。通过碳水化合物注释结果推测L.plantarum ST有着较强的水解或重排糖苷键能力，同时可能具有良好的降解纤维素潜力。良好糖苷水解能力为L.plantarum ST用于益生菌制品奠定了基础。18?15129630GT4GT2GH1GH13_31GH25GH36GH78GH4GT51GT113GT32GH8GT14CBM67GT101GH126GH92GH125GH20CBM32AA4AA1AA10CBM34CBM48CE2CE4CE9GH13GH13_29GH13_39GH31GH32GH85GT26GT5GT35GH65GH170GH73P-8LPL-1JDM1STWCFS1ZS2058ZJ316KPST-IIISTCC14917GT28GH2GH109GH38GH42GT111GE1图 5 L.plantarum CAZy数据库注释结果Fig.5 Annotation results of L.plantarum based on CAZy database2.6 毒力和耐药基因检索分析安全性评价是判断一株菌是否为益生菌最为关键的一步，一般通过研究菌株的抗生素抗性、有毒代谢物等指标评估乳酸菌的安全性30。抗生素滥用严重导致细菌产生耐药性。将L.plantarum ST与其他9 株L.plantarum的基因组经毒力因子、耐药基因数据库比对后发现均不存在毒力因子和耐药基因，故认为L.plantarum ST是安全的，并且没有潜在的耐药性。2.7 特有基因和益生基因分析为进一步挖掘L.plantarum ST功能基因同时解析其潜在的益生功能基因，继续对菌株的特有基因和益生基因进行分析。L.plantarum ST存在41 个特有基因，除去一些假定蛋白编码基因后，L.plantarum ST有3 个特有功能基因，表2为特有功能基因信息，分别为乙二醛还原酶基因（yvgN）、酪氨酸重组酶基因（xerC）、能量耦合因子转运蛋白基因（ecfT）。yvgN是一种加氢酶，对乙二醛进行加氢反应。ecfT是一种转运酶，负责摄入一些维生素及其他微量元素31。表 2 154 株L.plantarum中L.plantarum ST特有功能基因Table 2 Unique genes of L.plantarum ST 菌株基因注释L.plantarum STyvgN乙二醛还原酶 xerC酪氨酸重组酶ecfT能量耦合因子转运体跨膜蛋白基于前期对L.plantarum ST益生功能的研究13，使用Prokka软件注释到L.plantarum ST含有与益生特性相关的基因，具体信息见表3。群体感应AI-2信号分子合成的最关键蛋白酶是S-核糖同型半胱氨酸裂解酶（LuxS），AI-2在微生物种内和种间交流起到重要作用，在L.plantarum ST中发现了该基因。目前发现在超过80 种细菌中均存在LuxS基因32-33，细菌许多生理现象离不开LuxS介导的群体感应系统，因为这种系统对细菌有重要的调节作用34。L.plantarum ST还含有谷胱甘肽合成（garB、gshAB）、核黄素合成（ribBA、ribE、ribF）相关基因，有研究表明35谷胱甘肽不仅有助于维持正常的免疫系统功能，同时有整合解毒与抗氧化作用。此外，菌株在胃肠道存活率的增加也与谷胱甘肽有关36，该物质有助于菌株益生功能的发挥。核黄素（VB2），作为辅酶有助于代谢是其主要功能。NAPPH还原酶类（nfrA），参与脂肪酸、核苷酸和脂类的合成等多种合成代谢反应，该物质起到递氢体的作用在生物体内的化学反应中37。L.plantarum ST还注释到了黏附因子（atpC）基因，黏附因子有助于细菌在宿主体内的定植因为可使菌株牢固黏附于宿主细胞表面。研究还发现，L.plantarum ST存在提高宿主代谢能力（tagE）、乳酸和转运乳酸（pyk、ldh）、细胞壁特殊组分成分（ItaS）和生物活性肽（lmrA）产生相关的基因，可提高菌株益生功能。特有基因和益生基因分析发现，L.plantarum ST含有与黏附能力、信号分子传递和抗氧化能力等相关的益生功能基因。因此，本研究从基因组水平揭示L.plantarum ST具有益生特性相关基因，认为其是1 株具有潜在益生功能的菌株。92 2023,Vol.44,No.14 食品科学 生物工程表 3 L.plantarum ST具有的益生特性相关基因Table 3 Genes related to probiotic characteristics of L.plantarum ST基因注释益生特性LuxSS-核糖基同型半胱氨酸裂解酶 群体感应nfrA黄素单核苷酸还原酶（NADPH）氧化还原酶garB谷胱甘肽酰胺还原酶谷胱甘肽生物合成gshAB谷胱甘肽生物合成双功能蛋白lmrA多药耐药ABC转运蛋白ATP结合和渗透酶蛋白生物活性肽的生产ItaS脂磷壁酸合酶细胞壁成分tagE甘油-磷酸酯-葡萄糖基转移酶改善宿主新陈代谢atpCATP合酶链黏附因子的产生ribBA核黄素生物合成蛋白核黄素合酶ribE核黄素合酶ribF双功能核黄素激酶pyk丙酮酸激酶糖酵解过程ldhL-乳酸脱氢酶2.8 碳水化合物代谢特性为了评估L.plantarum ST对不同碳源的利用能力，利用API50CHL碳水化合物鉴定试剂条对ST的碳水化合物代谢进行研究，挖掘出可被ST利用的碳源，为碳水化合物代谢中的功能基因分析提供表型数据支持，同时为该菌株在之后的应用提供借鉴。表 4 L.plantarum ST API 50 CHL分析Table 4 Analysis of API 50 CHL in L.plantarum ST 碳水化合物24 h48 h类别碳水化合物24 h48 h类别阴性对照七叶灵B丙三醇D水杨苷B赤藻糖醇DD-纤维二糖BD-阿拉伯糖AD-麦芽糖BL-阿拉伯糖AD-乳糖BD-核糖AD-蜜二糖BD-木糖AD-蔗糖BL-木糖AD-海藻糖BD-侧金盏花醇D菊糖C甲基-D-吡喃木糖苷AD-松三糖CD-半乳糖AD-棉子糖CD-葡萄糖A淀粉CD-果糖A果糖CD-甘露糖A木糖醇DL-山梨糖AD-龙胆糖BL-鼠李糖AD-松二糖B半乳糖醇DD-来苏糖A肌醇DD-塔格糖A甘露醇DD-岩藻糖A山梨醇DL-岩藻糖A甲基-D-吡喃甘露糖苷AD-阿拉伯醇D甲基-D-吡喃葡萄糖苷AL-阿拉伯醇DN-乙酰葡萄糖胺A葡萄糖酸钾E苦杏仁苷B2-酮基葡萄糖酸钾E熊果苷B5-酮基葡萄糖酸钾E注：.利用；.不利用。A.单糖类；B.糖苷类和二糖类；C.多糖类；D.醇类；E.盐类。如表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.plantarum ST可代谢不同种类丰富的碳源，包括9 种单糖、12 种糖苷和二糖、1 种多糖及4 种醇类和2 种盐类物质，再次说明了ST具有较强的碳水化合物代谢能力。该结果与本实验KEGG注释以及CAZy分析吻合，利用基因组学和表型分析相结合，发现L.plantarum ST拥有较好的产CAZy能力，具备丰富的碳水化合物利用能力，可为开发相关酶制剂提供一定的参考价值。3 结 论对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益生功能开发及其生产应用奠