益生菌是一类达到一定数量时能够维持人体肠道菌群平衡，保护宿主健康的活性微生物[1-3]。乳酸菌在益生菌中扮演着重要角色，它们能够利用碳水化合物产生乳酸[4]，有助于蛋白质和钙的吸收，从而预防疾病[5]。1857年，法国微生物学家LOUIS PASTEUR首次发现乳酸菌，这标志着乳酸菌的出现[6]。之后越来越多的科学家开始专注研究乳酸菌，以期能够为人类带来更好的健康效果。产马乳酒乳杆菌（Lactobacillus kefiranofaciens）属于乳杆菌科（Lactobacteriaceae），乳杆菌属（Lactobacillus），呈长杆状的革兰氏阳性菌[7]，包含两个亚种，其中产马乳酒乳杆菌马乳酒亚种于2021年被国家卫健委列为国家新食品原料[8]，该亚种兼性厌氧且具有高产胞外多糖的生物学特性及抗氧化、抗肿瘤、抗过敏、改善肠道健康、免疫调节等益生功能，具有潜在应用价值[9]。此前，本实验室对西藏开菲尔粒的微生态组成进行了多年研究，并从中分离出多株产马乳酒乳杆菌，前期已经发现其中产马乳酒乳杆菌马乳酒亚种（Lactobacillus kefiranofaciens subsp.kefiranofaciens）ZW3可以黏附在小鼠肠道上并产生有益菌群[10]，通过脑肠轴调节色氨酸代谢缓解小鼠抑郁[11]。前期研究还发现，Lactobacillus kefiranofaciens subsp.kefiranofaciens ZW18有很好的免疫调节能力等益生特性。随着检测技术的发展，基因组测序为全面系统地认识乳酸菌提供了方便[12]，并为预测乳酸菌的功能提供帮助[13-14]。
为了深入研究该种属菌的特性及差异性，本研究对产马乳酒乳杆菌马乳酒亚种（Lactobacillus kefiranofaciens subsp.kefiranofaciens）ZW18的全基因组进行测序和分析，与亲缘菌株进行比较，并且发掘益生作用的相关基因，为其深入研究和进一步应用奠定基础。
1 材料与方法
1.1 材料与试剂
1.1.1 菌株产马乳酒乳杆菌马乳酒亚种（Lactobacillus kefiranofaciens subsp.kefiranofaciens）ZW18：保藏于天津科技大学食品科学与工程学院发酵食品与微生物资源开发研究室，美国国立生物技术信息中心（national venter for biotechnology information，NCBI）登录号为CP123735。
产马乳酒乳杆菌高加索酸奶粒亚种（Lactobacillus kefiranofaciens subsp.kefirgranum）、产马乳酒乳杆菌（Lactobacillus kefiranofaciens）、瑞士乳杆菌（Lactobacillus helveticus）、卷曲乳杆菌（Lactobacillus crispatus）、嗜酸乳杆菌（Lactobacillus acidophilus）等57株乳杆菌的全基因组序列均下载于NCBI的GenBank数据库。
1.1.2 培养基
MRS液体培养基：牛肉膏10.0 g，葡萄糖20.0 g，蛋白胨10.0 g，酵母膏5.0 g，氯化钠5.0 g，柠檬酸氢二胺2.0 g，磷酸氢二钾2.0 g，硫酸镁0.2 g，硫酸锰0.05 g，吐温-80 1.0 g，蒸馏水定容至1 000 mL，pH 6.2～6.4。115 ℃高压蒸汽灭菌15 min。MRS固体培养基在MRS液体培养基中加入2%琼脂粉。
1.1.3 试剂
液氮：天津市军粮城常福气体有限公司；其他试剂均为国产分析纯。
1.2 仪器与设备
Anaerobox Ⅳ厌氧培养箱：美国GeneScience公司；3-18K冷冻离心机：德国SIGMA公司；Nanopore测序仪：牛津纳米孔科技有限公司。
1.3 方法
1.3.1 菌株ZW18的活化与培养
挑取-80 ℃冰箱冻存管的菌株ZW18于MRS固体培养基平板上划线，37 ℃厌氧培养活化，至菌落长出；挑取单菌落于MRS液体培养基中，37 ℃厌氧培养，以2%（V/V）接种量接种于新鲜MRS液体培养基中培养至对数生长期，室温下14 000 r/min离心1 min，菌体沉淀于液氮中速冻后转移至-80 ℃储存以供后续测序。
1.3.2 全基因组测序
将样品交由北京百迈客生物科技有限公司进行高通量测序。细菌基因组完成图基于Illumina二代测序以及Oxford Nanopore三代测序平台实现完整组装，获得全面的基因组信息。
1.3.3 生物信息分析
通过软件Prodigal v2.6.3[15]进行基因预测；采用软件Infernal v1.1.3[16]预测基因组中的三类核糖体核酸（ribosomal ribonucleic acid，rRNA）；采用tRNAscan-SE v2.0[17]高精确预测基因组中的转运核糖核酸（transfer RNA，tRNA）；采用CRISPR Recognition Tool对基因组进行CRISPR预测；采用软件BPGA构建系统进化树；利用在线注释软件Ortho Venn 2（https：//orthovenn2.bioinfotoolkits.net/home）分析核心基因组、独特基因。使用基本局部比对搜索工具（basic local comparison search tool，BLAST）将预测得到的编码基因蛋白序列与公共蛋白数据库进行比对，如基因本体（gene ontology，GO）、基因的进化谱系：非监督直系群数据库（evolutionary genealogy of genes：non-supervised orthologous groups database，eggNOG）、京都基因与基因组百科全书（kyoto encyclopedia of genes and genomes，KEGG）等，获得编码基因的功能信息；采用在线注释网站RAST（https：//rast.nmpdr.org/）进行全基因组功能基因注释。使用Circos v0.66[18]完成全基因组精细圈图的绘制。
2 结果与分析
2.1 菌株ZW18基因组基本特征
为了更加全面直观的描述菌株ZW18基本特征，通过软件Circos v0.66制作基因组圈图，结果见图1。
图1 产马乳酒乳杆菌马乳酒亚种ZW18基因组圈图
Fig.1 Genome diagram of Lactobacillus kefiranofaciens subsp.kefiranofaciens
由图1可知，经全基因组测序和组装，菌株ZW18基因组共有一个完整闭合的染色体和两个质粒，全长为2 438 418 bp，鸟嘌呤（guanine，G）胞嘧啶（cytosine，C）含量为37.27%。菌株ZW18中编码tRNA基因64个，rRNA基因12个，共预测到11个潜在CRISPR序列，其中9个分布在染色体中，2个在质粒2中。菌株ZW18在GO、eggNOG和KEGG数据库注释结果统计见图2。由图2A可知，通过GO数据库共注释到的基因个数为1 843个，其中与细胞膜构建、参与催化反应及代谢相关基因占比最多；由图2B可知，通过eggNOG数据库注释可得关于细胞复制、重组和修复基因数量最多，其次为碳水化合物的利用与一般功能预测相关基因，分别为273个、175个、162个；由图2C可知，通过KEGG数据库注释，菌株ZW18中ABC转运蛋白数量最多，能够把代谢物、药物和各种毒素从细胞中转运出来，维持细胞稳态。
图2 GO数据库（A）、eggNOG数据库（B）、KEGG数据库（C）基因注释统计结果
Fig.2 Statistical results of gene annotation in GO database (A),eggNOG database (B), and KEGG database (C)
2.2 全基因组比较分析
2.2.1 系统发育树分析
为了确定近缘菌株，根据核心基因序列，利用基因组比较平台BPGA对菌株ZW18及其他57株乳杆菌进行系统发育树的构建，结果见图3。由图3可知，菌株ZW18与产马乳酒乳杆菌马乳酒亚种ZW3、产马乳酒乳杆菌（Lactobacillus kefiranofaciens）ATCC 43761、产马乳酒乳杆菌高加索酸奶粒亚种（Lactobacillus kefiranofaciens subsp. kefirgranum）DSM 10550及产马乳酒乳杆菌（Lactobacillus kefiranofaciens）1207遗传距离最近。除此之外与其较近的还有瑞士乳杆菌（Lactobacillus helveticus）、卷曲乳杆菌（Lactobacillus crispatus）与嗜酸乳杆菌（Lactobacillus acidophilus）。
图3 基于核心基因序列构建菌株ZW18的系统发育树
Fig.3 Phylogenetic tree of strain ZW18 constructed based on the core gene sequences
2.2.2 核心基因组与独特基因分析
利用在线注释软件OrthoVenn2选取与菌株ZW18遗传距离较近的产马乳酒乳杆菌马乳酒亚种ZW3、产马乳酒乳杆菌ATCC 43761、产马乳酒乳杆菌高加索酸奶粒亚种DSM 10550及产马乳酒乳杆菌1207进行全基因组序列比对，结果见图4。核心基因组指在所有个体中都存在的基因且非常保守。由图4可知，这5株产马乳酒乳杆菌共有的核心基因家族为1 553个。菌株ZW18的独特基因数量最多，有30个，说明菌株ZW18更具有多样性，有更多研究价值。
图4 核心基因组和独特基因Venn图
Fig.4 Venn diagram of core genome and specific genes
2.3 益生功能相关基因分析
2.3.1 与应激反应相关基因分析
在菌株ZW18中共注释出19种与应激反应相关的基因，具体结果见表1。
表1 菌株ZW18中与应激反应相关主要基因分析结果
Table 1 Analysis results of major genes related to stress response in strain ZW18
应激反应 基因编号 基因注释烷基过氧化氢还原酶（EC:4.1.1.44）谷胱甘肽还原酶（EC:1.8.1.7）铁蛋白类抗氧化蛋白非特异性DNA结合蛋白Dps铁氧化物酶（EC 1.16.3.1）NADH过氧化物酶氧化应激GE000333 GE000493 GE001957 GE000607 GE000607 GE001008 GE000495 GE000511 GE001004 GE001708 GE000065 GE000122 GE000972 GE001082 GE001099 GE001135 GE001249 GE002505转录调节因子，Crp/Fnr家族氧化还原敏感转录调节因子（富含AT的DNA结合蛋白）
续表
应激反应 基因编号 基因注释GE000340 GE000949 GE001366 GE002113烟酸磷酸核糖转移酶（EC 2.4.2.11）冷休克热休克GE001660 GE002011 GE002516 GE000145 GE000281 GE000727 GE000725 GE000724 GE000726 GE001071 GE000722烟酰胺酶（EC 3.5.1.19）SIR2家族NAD+依赖性蛋白脱乙酰酶NAD依赖性甘油醛-3-磷酸脱氢酶（EC 1.2.1.12）冷休克蛋白CspA热诱导转录抑制剂HrcA伴侣蛋白DnaK伴侣蛋白DnaJ伴侣蛋白GrpE tmRNA结合蛋白SmpB延伸因子LepA
由表1可知，菌株ZW18中与氧化应激相关的基因有12种。烷基过氧化氢还原酶（alkyl hydroperoxide reductase）是分布在原核生物和真核生物中的一大家族硫醇特异性抗氧化蛋白的关键成分[19]。非特异性脱氧核糖核酸（deoxyribonucleic acid，DNA）结合蛋白Dps属于铁蛋白超家族，通过氧化和储存铁可以有效去除H2O2和Fe2+的毒性，保护DNA、蛋白质和膜脂免受自由基损害[20]。烟酰胺腺嘌呤二核苷酸（nicotinamide adenine dinucleotide，NADH）过氧化物酶能够使菌株抵抗过氧化氢等的毒性，保护细胞生长[21]。Sir2蛋白在维持细胞抗胁迫能力、基因组稳定性等方面起重要作用[22]。除抵抗氧化应激外，菌株ZW18中还含有6种热休克、1种冷休克相关蛋白基因，帮助菌体抵御外界环境不良刺激。
2.3.2 与胃肠道耐受相关基因分析
人体胃肠道（gastrointestinal tract，GIT）具有许多防御机制来阻止病原体定植，包括胃酸和肠液，在保护机体、抵抗病原菌侵染方面起重要作用，同时也是考验益生菌能否在肠道中起作用的关键[23]。因此，对菌株ZW18基因组中与胃肠应激耐受性相关基因进行挖掘，结果见表2。由表2可知，菌株ZW18基因组中含有大量有助于细菌耐酸、耐胆盐的基因，表明该菌株有能力对胃肠道中的应激做出反应。F（1）F（0）-三磷酸腺苷（adenosine triphosphate，ATP）合成酶（F（1）F（0）-ATP synthetase）是重要的耐酸相关酶，活性随环境pH降低而增加。当pH下降时，F（1）F（0）-ATP合成酶会利用水解ATP的方式将H+从细胞内泵出以保持pH的平衡，同时保证细胞内正常代谢，从而在较低的pH环境下继续生长并产生有益的物质[24]。Na+/H+逆向转运蛋白在维持细胞pH稳态和Na+动态平衡中起重要作用，利用该Na+外排系统能够将多余Na+排至胞外从而维持细胞内外渗透压平衡[25]。在菌株ZW18中也发现了两个胆酰甘氨酸水解酶基因和两个胆盐转运蛋白基因。胆酰甘氨酸水解酶是一种胆盐水解酶，具有去交联胆盐的能力，被认为具有解除胆盐毒性的作用，延长菌体在肠道中的生存时间[26]。
表2 菌株ZW18中与胃肠应激耐受性相关基因分析结果
Table 2 Analysis results of genes related to gastrointestinal stress tolerance in strain ZW18
基因编号 基因注释 基因功能GE000203 GE000204 GE000205 GE002131 GE000450 GE000451 GE000367 GE000368 GE000365 GE000366 F（1）F（0）-ATP合成酶Na+/H+逆向转运蛋白天冬氨酸/谷氨酸渗透酶AcaP鸟氨酸脱羧酶胆酰甘氨酸水解酶酸耐受胆盐耐受胆盐转运蛋白
2.3.3 与黏附和聚集相关蛋白质基因分析
益生菌对肠上皮细胞的黏附能力是其发挥益生作用的首要标准，黏附能力强的菌株更容易在肠道中停留，同时竞争性地排除有害病原体[27]。菌株ZW18基因组中与黏附和聚集相关的多种蛋白基因预测结果见表3。由表3可知，ZW18基因组中预测出7种与黏附和聚集相关的蛋白质基因。烯醇化酶在多种细菌表面上普遍存在，能够有效促进细菌与宿主之间的结合，增强细菌黏附力[28]；纤维连接蛋白结合蛋白是一种重要的多结构域糖蛋白，具有多种黏附特性[29]；S-层蛋白存在于乳杆菌中，能够增强菌株黏附力，抵御致病菌入侵宿主[30]；LPxTG细胞壁锚定结构域蛋白可以增加菌株黏附率，在一定程度上可以引起免疫细胞的响应，进而调节肠道炎症[31]。
表3 菌株ZW18中与黏附和聚集相关蛋白相关基因的分析结果
Table 3 Analysis results of genes related to adhesion and aggregation proteins in strain ZW18
基因编号 基因注释GE000355 GE001758 GE000640 GE001015 GE000372 GE002007 GE000644 GE001048 GE000743烯醇化酶纤维连接蛋白结合蛋白S-层蛋白粘蛋白结合蛋白脂蛋白信号肽酶LPxTG细胞壁锚定结构域蛋白EF-Tu表面黏附蛋白
2.3.4 与免疫相关特性基因分析
益生菌由于其黏附特性和竞争性抑制病原菌的能力使得肠道内含有大量益生菌群，对营养物质的摄入、疾病治疗等具有调节作用[32]。前期研究发现，菌株ZW18中具有多种黏附蛋白、细菌素和抵抗应激反应的基因。这些益生基因能够帮助菌株定植于肠上皮细胞，抑制病原菌，对人体肠道进行免疫调节[33]。除此之外，菌株ZW18基因组中预测出胞外蛋白基因（GE001482、GE002271）。据文献报道，这些基因可以辅助肠上皮细胞与免疫细胞相互作用，促进细胞因子的分泌[34]。前期研究也发现，菌株ZW18具有免疫辅助治疗小鼠黑色素瘤的效果[35]。
2.3.5 与细菌素相关特性基因分析
细菌素是细菌在代谢过程中分泌的一种小分子抗菌肽，对其自身往往无抑制作用，却能够抑制其他相近种属菌株生长[36]。细菌素无毒且抑菌作用强，是胃肠道中抑制病原菌的关键物质，可作为抗菌材料，具有较好的开发利用价值。因此，对菌株ZW18基因组中与细菌素相关特性基因进行分析，结果见表4。由表4可知，菌株ZW18中共预测出5种与细菌素相关的基因，包括羊毛硫细菌素Class-Ⅳ、肠菌素A、瑞士乳杆菌素J同源物及细菌素免疫蛋白，其中，羊毛硫细菌素Class-Ⅳ数量最多，共69个。
表4 菌株ZW18中与细菌素相关特性基因的分析结果
Table 4 Analysis results of genes related to bacteriocin in strain ZW18
基因编号 基因注释GE002377～GE002394（染色体）GE002467～GE002478（质粒1）GE002395～GE002433（质粒2）羊毛硫细菌素Class ⅣGE001027 GE001386 GE001861 GE000967 GE000046肠菌素A瑞士乳杆菌素J同源物细菌素免疫蛋白
2.4 菌株ZW18合成营养物质相关基因分析
益生菌的营养合成（维生素和必需氨基酸）正成为新型发酵食品的一个研究热点。维生素是人体生理反应中重要的微量营养素，维生素缺乏会增加患传染病、过敏性疾病和炎症性疾病的风险[37]。对菌株ZW18基因组中与维生素合成的相关基因进行预测，结果见表5。由表5可知，在菌株ZW18基因组中预测到34种与维生素合成相关的基因，分别与硫胺素、核黄素、吡哆醇、生物素和叶酸合成有关。
表5 菌株ZW18中与合成维生素相关基因分析结果
Table 5 Analysis results of genes related to synthetic vitamins in strain ZW18
维生素种类 基因编号 基因注释硫胺素磷酸甲基嘧啶激酶（EC 2.7.4.7）硫胺素焦磷酸激酶（EC 2.7.6.2）1-脱氧-D-木酮糖5-磷酸合酶（EC 2.2.1.7）3-酮酰辅酶A硫醇酶（EC 2.3.1.16）核黄素 生物素ECF转运蛋白的底物特异性成分BioY生物素操纵子阻遏剂生物素-蛋白连接酶（EC 6.3.4.15）生物素ECF转运蛋白能量模块的ATP酶组分BioM代谢FMN腺苷酸转移酶（EC 2.7.7.2）假设蛋白YebC核黄素ECF转运蛋白的底物特异性成分RibU多抗菌挤出蛋白（Na（+）/药物抗转运蛋白），MDR外排泵吡哆醇GE001794 GE000786 GE001176 GE000069 GE001140 GE001375 GE001368 GE001373 GE001375 GE001266 GE000728 GE001841 GE000074 GE001207 GE001336 GE001600 GE002131 GE000045 GE000149 GE000627 GE001264 GE001544 GE000137 GE000728 GE000449 GE000885 5-氨基-6-（5'-磷酸比基氨基）尿嘧啶磷酸酶核黄素激酶（EC 2.7.1.26）吡哆胺5'-磷酸氧化酶（EC 1.4.3.5）生物素叶酸GE000466 GE000145 GE001176 GE001085 GE001250 GE001251 GE000581 GE000258 GE000359 GE001146 GE000927 GE000360 GE001446 GE000333 GE000498 GE000499 GE002043吡哆醛激酶（EC 2.7.1.35）NAD依赖性甘油醛-3-磷酸脱氢酶（EC 1.2.1.12）1-脱氧-D-木酮糖5-磷酸合酶（EC 2.2.1.7）假设的NagD样磷酸酶细胞分裂蛋白FtsH（EC 3.4.24.-）tRNA（Ile）-赖氨酸合成酶（EC 6.3.4.19）次黄嘌呤-鸟嘌呤磷酸核糖转移酶（EC 2.4.2.8）叶基聚谷氨酸合酶（EC 6.3.2.17）胸苷酸合酶（EC 2.1.1.45）叶酸ECF转运蛋白的底物特异性成分FolT 5-甲酰基四氢叶酸环连接酶（EC 6.3.3.2）二氢叶酸还原酶（EC 1.5.1.3）对氨基苯甲酸合酶，酰胺转移酶组分（EC 2.6.1.85）醇脱氢酶（EC 1.1.1.1）
续表
维生素种类 基因编号 基因注释GE000645 GE001649 GE001266 GE000441 GE001639 GE001640 GE001794 GE001640 GE001362甲酸盐--四氢叶酸连接酶（EC 6.3.4.3）丝氨酸羟甲基转移酶（EC 2.1.2.1）亚甲基四氢叶酸脱氢酶（NADP+）（EC 1.5.1.5）磷酸核糖甘氨酰胺甲酰转移酶（EC 2.1.2.2）磷酸甲基嘧啶激酶（EC 2.7.4.7）磷酸核糖氨基咪唑甲酰胺甲酰转移酶（EC 2.1.2.3）5,10-亚甲基四氢叶酸还原酶（EC 1.5.1.20）
人体需要从食物中摄取的8种必需氨基酸为异亮氨酸、亮氨酸、赖氨酸、蛋氨酸、苯丙氨酸、苏氨酸、缬氨酸、色氨酸[38]，而半胱氨酸、酪氨酸、组氨酸和精氨酸是对婴儿和正在成长中的儿童很重要的4种氨基酸[39]。对菌株ZW18基因组中与合成氨基酸相关的基因进行预测，结果见表6。由表6可知，从菌株ZW18基因组中共预测到24种与氨基酸合成相关的基因，分别与精氨酸、蛋氨酸、苏氨酸、赖氨酸生物合成有关。
表6 菌株ZW18中与合成氨基酸相关基因分析结果
Table 6 Analysis results of genes related to synthetic amino acids in strain ZW18
描述 氨基酸种类 基因编号 基因注释赖氨酸天冬氨酸激酶（EC:2.7.2.4）天冬氨酸-半醛脱氢酶（EC:1.2.1.11）二氨基庚二酸脱羧酶（EC:4.1.1.20）N-乙酰基二氨基庚二酸脱乙酰酶2,3,4,5-四氢吡啶-2,6-二羧酸N-乙酰转移酶（EC:2.3.1.89）N-乙酰基-L,L-二氨基庚二酸脱乙酰酶（EC 3.5.1.47）4-羟基-四氢二吡啶甲酸还原酶（EC:1.17.1.8）高丝氨酸脱氢酶（EC 1.1.1.3）同型半胱氨酸S-甲基转移酶（EC 2.1.1.10）成人所需氨基酸蛋氨酸ABC转运蛋白底物结合蛋白蛋氨酸GE000317 GE000700 GE000324 GE000318 GE000323 GE000319 GE001932 GE000323 GE000321 GE000322 GE000698 GE000461 GE000022 GE000031 GE000036 GE001171 GE001453 GE001469 GE001493 GE001550 GE001551 GE001807 GE001953 GE002019 GE000363 GE001553 GE001065 GE000697 GE000715 GE001361 GE000284蛋氨酸ABC转运蛋白ATP结合蛋白S-腺苷甲硫氨酸合成酶（EC 2.5.1.6）高丝氨酸激酶（EC 2.7.1.39）高丝氨酸O-琥珀酰转移酶（EC 2.3.1.46）5-甲基四氢蝶酰三谷氨酸-同型半胱氨酸甲基转移酶（EC 2.1.1.14）S-腺苷同型半胱氨酸核苷酶（EC 3.2.2.9）
续表
描述 氨基酸种类 基因编号 基因注释GE000331 GE000994 GE000995 GE001360 GE001894蛋氨酸ABC转运蛋白渗透蛋白S-核糖基同型半胱氨酸裂解酶（EC 4.4.1.21）婴儿和正在成长中的儿童所需氨基酸苏氨酸半胱氨酸精氨酸GE001362 GE000699 GE000716 GE001902 GE000710 GE001469 5,10-亚甲基四氢叶酸还原酶（EC 1.5.1.20）苏氨酸合成酶（EC:4.2.3.1）半胱氨酸合成酶（EC:2.5.1.47）核糖体蛋白-丝氨酸乙酰转移酶（EC:2.3.1.-）精氨酸和脯氨酸代谢（EC:4.1.1.17）精氨酸生物合成（EC:2.7.2.2）
3 结论
通过对产马乳酒乳杆菌马乳酒亚种ZW18进行全基因组测序可得其基因组中共有一个完整闭合的染色体和两个质粒，全长为2 438 418 bp，GC含量为37.27%。通过数据库注释可知其生长稳定且对碳水化合物利用基因较多。菌株ZW18有1 553个核心基因，30个独特基因，是产马乳酒乳杆菌中具有多样性的菌株。菌株ZW18中含有很多益生相关基因，在抵抗恶劣环境胁迫、抗氧化、耐受胃肠道、调节免疫、抗菌等方面有突出作用。除此之外，该菌株能够合成人体所需维生素，如硫胺素、核黄素等，也可以合成成人及婴幼儿生长所需氨基酸，如赖氨酸、蛋氨酸等，具有潜在研究价值。综上，从产马乳酒乳杆菌马乳酒亚种ZW18的全基因组分析及功能特性可以看出，菌株ZW18是一株益生功能强且具有潜在应用价值的菌株，为其今后的进一步深入研究与应用提供了基础数据。
[1]毛友辉，朱薇，邓放明，等.益生菌的研究进展及其在食品中的应用[J].现代生物医学进展，2011，11（20）：3978-3980，3977.
[2]孙雪姣，王一然，丁瑞雪，等.益生菌火龙果酸奶的研发[J].食品研究与开发，2018，39（19）：99-103.
[3] TORRES-MARAVILLA E, LENOIR M, MAYORGA-REYES L, et al.Identification of novel anti-inflammatory probiotic strains isolated from pulque[J].Appl Microbiol Biotechnol,2016,100:385-396.
[4]周莉，平洋，谭静，等.益生菌中乳酸菌概况及检测技术的研究进展[J].中国调味品，2019，44（10）：190-194.
[5]王怡明.乳酸菌及其生物工程研究新进展[J].科技风，2022（10）：152-154.
[6]弥超，高维东，王玉，等.益生菌的功能研究及其在功能性食品中的应用[C]//中国营养学会，中国疾病预防控制中心营养与健康所，农业农村部食物与营养发展研究所，中国科学院上海营养与健康研究所，华中科技大学公共卫生学院.中国营养学会第十五届全国营养科学大会论文汇编，2022：534.
[7]段晓微.马乳酒样乳杆菌ZW3高密度培养及抗过敏功能的初步研究[D].天津：天津科技大学，2019.
[8]国家卫生健康委.关于蝉花子实体（人工培植）等15种“三新食品”的公告2020年第9号[J].中国食品卫生杂志，2021，33（1）：68.
[9]王艳萍，严婷，王亚琦，等.马乳酒样乳杆菌马乳酒样亚种及其研究进展[J].微生物学报，2022，62（2）：414-420.
[10]XING Z Q,TANG W,YANG Y,et al.Colonization and gut flora modulation of Lactobacillus kefiranofaciens ZW3 in the intestinal tract of mice[J].Probiotics Antimicro,2018,10:374-382.
[11]SUN Y,GENG W T,PAN Y J,et al.Supplementation with Lactobacillus kefiranofaciens ZW3 from Tibetan Kefir improves depression-like behavior in stressed mice by modulating the gut microbiota[J].Food Funct,2019,10(2):925-937.
[12]余中节，赵洁，宋宇琴，等.全基因组关联分析在乳酸菌研究中的应用[J].中国食品学报，2019，19（10）：307-314.
[13] PASOLLI E, DE FILIPPIS F, MAURIELLO I E, et al.Large-scale genome-wide analysis links lactic acid bacteria from food with the gut microbiome[J].Nat Commun,2020,11(1):2610.
[14]张彦位，杨玲，路江浩，等.乳酸菌全基因组测序的应用进展[J].食品工业科技，2022，43（15）：444-450.
[15] HYATT D, CHEN G L, LOCASCIO P F, et al.Prodigal: prokaryotic gene recognition and translation initiation site identification[J].BMC Bioinformatics,2010,11(1):119.
[16] NAWROCKI E P, EDDY S R.Infernal 1.1: 100-fold faster RNA homology searches[J].Bioinformatics,2013,29(22):2933-2935.
[17] CHAN P P, LOWE T M.tRNAscan-SE: Searching for tRNA genes in genomic sequences[J].Methods Mol Biol,2019,1962:1-14.
[18] KRZYWINSKI M, SCHEIN J, BIROL I, et al.Circos: An information aesthetic for comparative genomics[J].Genome Res,2009,19(9):1639-1645.
[19]李海燕，胡梦蝶，刘洋，等.谷氨酸棒杆菌中烷基过氧化物还原酶（CgAhp）抵御环境胁迫的作用机制[J].曲阜师范大学学报（自然科学版），2022，48（4）：101-108.
[20]方丽.大肠杆菌DNA结合蛋白克隆表达纯化及体外自组装功能研究[D].淄博：山东理工大学，2014.
[21]段希宇，叶陵，刘成国，等.乳酸菌的抗氧化作用机制[J].微生物学杂志，2017，37（3）：111-115.
[22]朱桐，薛原野，李欣桐，等.去乙酰化酶SIRT7在细胞衰老中的研究进展[J].安徽医学，2021，42（11）：1314-1317.
[23]熊涛，刘妍妍，黄涛，等.副干酪乳杆菌NCU622耐酸耐胆盐及其黏附性能[J].食品科学，2015，36（5）：93-98.
[24]张阳玲.酸胁迫下乳酸乳球菌细胞壁相关基因调控机制研究[D].天津：天津大学，2020.
[25]徐宁，程海娇，刘清岱，等.细菌Na+/H+逆向转运蛋白的研究进展[J].微生物学通报，2015，42（10）：2002-2011.
[26]杨玲双，谢新强，李滢，等.屎肠球菌132胆盐水解酶基因的克隆表达与酶学特性[J].现代食品科技，2022，38（8）：35-43.
[27]周铄浛，杨柳，安彬，等.乳酸菌体外黏附性方法的评价及相关性研究[J].食品科技，2020，45（8）：20-28.
[28]王思淼，魏晓，赵向娜，等.长双歧杆菌NCC2705纯化蛋白烯醇化酶、Tuf蛋白的黏附作用研究[J].军事医学，2014，38（5）：333-336.
[29]武琪，张玉娟，刘桐，等.牛支原体黏附素及黏附素结合蛋白研究进展[J].畜牧兽医学报，2023（8）：3206-3216.
[30]石晓璐，张英春，张兰威，等.乳酸杆菌S-层蛋白的特性及功能性研究进展[J].微生物学通报，2017，44（5）：1206-1213.
[31]劳立峰，吴振，潘道东.罗伊氏乳杆菌SH23中LPxTG表面蛋白的原核表达及其抗炎特性研究[C]//中国食品科学技术学会.第十六届益生菌与健康国际研讨会摘要集，2021：18.
[32]吕秀莉，岳莹雪，平丽筠，等.益生菌黏附机制及其拮抗肠道致病菌研究进展[J].食品科学，2023，44（9）：313-320.
[33]唐曼玉，王晚晴，强敬雯，等.益生菌与肠道菌群、免疫调节的相互作用与机制研究进展[J].食品工业科技，2022，43（16）：486-493.
[34]任效东.益生菌分泌的胞外蛋白与肠道粘膜免疫关系的研究进展[J].农技服务，2016，33（3）：109.
[35]ZHAO J Q,WANG Y P,WANG J R,et al.Lactobacillus kefiranofaciens ZW18 from Kefir enhances the anti-tumor effect of anti-programmed cell death 1(PD-1)immunotherapy by modulating the gut microbiota[J].Food Funct,2022,13(19):10023-10033.
[36]滕坤玲，钟瑾.益生菌产生的细菌素及其功能机制[J].微生物学报，2022，62（3）：858-868.
[37]孔维溧，芦鑫荣，侯琳琳，等.维生素与免疫系统健康[J].四川大学学报（医学版），2023，54（1）：7-13.
[38]盛勤芳.人体内八种必需氨基酸的初步探究及其意义[J].科技视界，2014（28）：239.
[39] SAULNIER D M, SANTOS F, ROOS S, et al.Exploring metabolic pathway reconstruction and genome-wide expression profiling in Lactobacillus reuteri to define functional probiotic features[J].PloS one,2011,6(4):e18783.