微生物学报  2023, Vol. 63 Issue (7): 2880-2898   DOI: 10.13343/j.cnki.wsxb.20220838.
http://dx.doi.org/10.13343/j.cnki.wsxb.20220838
中国科学院微生物研究所，中国微生物学会

文章信息

伍亚龙, 杨恺, 史梅莓, 吕鹏军, 王勇, 张其圣. 2023
WU Yalong, YANG Kai, SHI Meimei, LYU Pengjun, WANG Yong, ZHANG Qisheng.
全基因组测序揭示两株泡菜源植物乳杆菌基因型差异和潜在益生特性
Whole-genome sequencing reveals genotypic differences and potential probiotic properties of two Paocai-derived Lactiplantibacillus plantarum strains
微生物学报, 63(7): 2880-2898
Acta Microbiologica Sinica, 63(7): 2880-2898

文章历史

收稿日期：2022-11-08
网络出版日期：2023-02-21
 Abstract              PDF               Figures               Tables
引用本文
伍亚龙, 杨恺, 史梅莓, 吕鹏军, 王勇, 张其圣. 全基因组测序揭示两株泡菜源植物乳杆菌基因型差异和潜在益生特性[J]. 微生物学报, 2023, 63(7): 2880-2898.
Wu Yalong, Yang Kai, Shi Meimei, Lyu Pengjun, Wang Yong, Zhang Qisheng. Whole-genome sequencing reveals genotypic differences and potential probiotic properties of two Paocai-derived Lactiplantibacillus plantarum strains[J]. Acta Microbiologica Sinica, 2023, 63(7): 2880-2898.
全基因组测序揭示两株泡菜源植物乳杆菌基因型差异和潜在益生特性
伍亚龙2,3 , 杨恺1 , 史梅莓1 , 吕鹏军1 , 王勇2,3 , 张其圣1,3     
1. 四川省食品发酵工业研究设计院有限公司, 四川 成都 611130;
2. 四川益动源生物科技有限公司, 四川 眉山 620000;
3. 四川东坡中国泡菜产业技术研究院, 四川 眉山 620000
收稿日期：2022-11-08；接受日期：2023-02-07；网络出版日期：2023-02-21
基金项目：四川省科技创新苗子工程培育项目(2022110)
摘要：[目的] 为了探究植物乳杆菌(Lactiplantibacillus plantarum)基因型差异和潜在益生特性，采用全基因组测序技术对其进行测序并解析基因组序列及生物特性。[方法] 本研究基于HiSeq和PacBio测序平台，对团队前期从四川多代泡菜中分离获得、体外益生特性评价良好的潜在益生菌菌株L. plantarum Eden-Star PC06和L. plantarum Eden-Star PC108的全基因组进行测序。利用相关生物信息学软件对原始数据进行组装及其后续的功能注释、分子进化、菌株安全性、次级代谢产物合成基因簇以及益生特性相关基因进行分析。[结果] 通过基因组装得到了2株植物乳杆菌的全基因组信息，L. plantarum Eden-Star PC06和Eden-Star PC108基因组大小分别为3 163 902 bp和3 205 054 bp；GC含量分别为44.68%和44.67%；分别包含3 161个和3 197个DNA编码序列；功能基因数据库比对结果显示2株菌在碳水化合物利用、氨基酸利用和糖基转移酶等基因上得到大量注释；通过比对数据库，在2株植物乳杆菌全基因组上发现了4个与肠液耐受相关的胆盐水解酶基因、完整的植物乳杆菌细菌素合成相关基因簇和抵御多种胁迫的益生相关基因。[结论] 本研究通过全基因组测序在基因水平上探究了L. plantarum Eden-Star PC06和Eden-Star PC108基因型差异和益生特性基因，证明L. plantarum Eden-Star PC06和Eden-Star PC108是2株有应用前景的益生菌菌株，以期为筛选优良益生菌菌株和评价其益生特性提供遗传学基础。
关键词：植物乳杆菌    全基因组测序    比较基因组    功能基因    胆盐水解酶    细菌素    
Whole-genome sequencing reveals genotypic differences and potential probiotic properties of two Paocai-derived Lactiplantibacillus plantarum strains
WU Yalong2,3 , YANG Kai1 , SHI Meimei1 , LYU Pengjun1 , WANG Yong2,3 , ZHANG Qisheng1,3     
1. Sichuan Food and Fermentation Industry Research & Design Institute Co., Ltd., Chengdu 611130, Sichuan, China;
2. Sichuan Eden Biology Technology Co., Ltd., Meishan 620000, Sichuan, China;
3. Sichuan Dongpo Chinese Paocai Industrial Technology Research Institute, Meishan 620000, Sichuan, China
Received: 8 November 2022; Accepted: 7 February 2023; Published online: 21 February 2023
*Corresponding author: ZHANG Qisheng, E-mail: bigbeastone@163.com.
Foundation item: This work was supported by the Sichuan Science and Technology Innovation Seedling Project Cultivation Project (2022110)
Abstract: [Objective] To explore the genotypic differences and potential probiotic properties of Lactiplantibacillus plantarum by using whole genome sequencing technology. [Methods] Based on the HiSeq and PacBio, this study sequenced the potential probiotic strains L. plantarum Eden-Star PC06 and L. plantarum Eden-Star PC108 with good evaluation of in vitro probiotic properties isolated from Sichuan multi-generation Paocai in the previous study. The relevant bioinformatics software was used to assemble the raw data and analyze the subsequent functional annotation, molecular evolution, strain safety, secondary metabolite synthesis gene clusters, and probiotic-related genes. [Results] The whole-genome information of two L. plantarum strains was obtained through gene assembly, and the genome sizes of L. plantarum Eden-Star PC06 and Eden-Star PC108 were 3 163 902 bp and 3 205 054 bp, respectively. The GC content of Eden-Star PC06 and Eden-Star PC108 was 44.68% and 44.67%, respectively, and they contained 3 161 and 3 197 DNA coding sequences, respectively. The comparison results of functional gene database showed that the two strains were extensively annotated on genes such as carbohydrate utilization, amino acid utilization, and glycosyltransferase. We found four genes related to intestinal fluid tolerance, intact gene clusters related to bacteriocin synthesis of L. plantarum, and probiotic genes resistant to multiple stresses on the whole genome of two strains of L. plantarum in the comparison database. [Conclusion] This study explores the genotypic differences and potential probiotic properties of two L. plantarum strains Eden-Star PC06 and Eden-Star PC108 at the gene level by whole-genome sequencing, and it has proved that L. plantarum Eden-Star PC06 and Eden-Star PC108 are potential probiotic strains with promising applications. This study is expected to provide a genetic basis for screening out excellent probiotic strains and evaluating their probiotic properties.
Keywords: Lactiplantibacillus plantarum    whole genome sequencing    comparative genomics    functional gene    bile salt hydrolase    bacteriocin    

植物乳杆菌(Lactiplantibacillus plantarum)是一种兼性异型乳酸菌，广泛存在于环境、人体胃肠道和发酵蔬菜食品中。近年来，大量研究表明，植物乳杆菌是一种多功能益生菌，具有抗氧化、降解胆固醇、维持宿主肠道菌群平衡和提高人体免疫等功能[1-2]。Choi等[3]从泡菜中分离出一株植物乳杆菌EM，经研究发现具有降低血胆固醇的作用。Wang等[4]的研究发现，植物乳杆菌LPL-1能合成具有抗菌活性的植物乳杆菌素，可以有效代替抗菌素抑制致病菌的生长。由于植物乳杆菌的益生特性，使其在食品和医用等领域都有广泛应用。

近年来高通量测序技术的不断发展进步，为研究人员开拓了一个全新的研究方向。2004年第1株分离自人体口腔的植物乳杆菌WCFS1完成全基因组测序后，截至2022年8月，已有超过100株植物乳杆菌完成全基因组测序[5]。而现有的生理生化方法对植物乳杆菌潜在益生特性的研究已十分有限，通过全基因组测序探究基因组与潜在益生特性之间关系，是一条筛选优质益生菌菌株的有效途径[6]。有研究指出，植物乳杆菌EM是一株具有抗菌活性、降胆固醇和耐酸胆汁的益生菌，通过全基因组分析其包含耐酸耐胆盐、耐温耐氧化等基因发现，其同时含有与抗菌活性相关的抗菌素基因簇以及与降低血清胆固醇相关的胆盐水解酶基因[7]。

本团队前期从四川多代泡菜中分离出114株植物乳杆菌，并对其生物特性和安全性进行评价，最终得到13株具有潜在益生特性的植物乳杆菌，从中选取2株植物乳杆菌Eden-Star PC06和Eden-Star PC108进行全基因组测序，其中Eden-Star PC06耐酸耐胆盐能力较为突出，而Eden-Star PC108拥有较强的自聚集和共聚集能力。本研究中完成图采用二代与三代结合即Illumina HiSeq+PacBio的测序方式，每个样品同时提供不低于基因组100×的PacBio测序数据和100×Illumina测序数据，保证更完整更精确的组装，完成图可以避免小质粒( < 15 kb)信息的丢失，保证获得包含质粒的完整基因组。利用比较基因组和功能基因注释等方法，从基因组水平上揭示菌株间的差异，并进一步挖掘其益生相关功能基因，以判断菌株作为益生菌菌株的潜力，同时为选优良益生菌菌株和评价其益生特性提供遗传学基础。

1 材料与方法 1.1 实验菌株信息

本研究测定的2株L. plantarum Eden-Star PC06 (下文中简称为PC06)和Eden-Star PC108 (下文中简称为PC108)由四川益动源生物科技有限公司保藏，均分离自四川多代泡菜，已在中国普通微生物菌种保藏管理中心(China General Microbial Culture Preservation and Management Center, CGMCC)注册专利，保藏号分别为No. 25060和No. 25061。项目已上传国家微生物科学数据中心(National Microbiology Data Center, NMDC)，项目编号(NMDC accession)为NMDC10018235，基因组编号(accession number)分别为NMDC60045920和NMDC60045921。

本研究选取12株完成全基因组完成图且有研究报道具有益生特性的植物乳杆菌作为参考菌株(表 1)，参与比较基因组的菌株信息均来自NCBI (National Center for Biotechnology Information, https://www.ncbi.nlm.nih.gov/) GenBank数据库。

表 1. 植物乳杆菌参考基因组信息 Table 1. Reference genome information of Lactiplantibacillus plantarum
Strains Whole genome size (bp) Origin GenBank ID
L. plantarum WCFS1 3 308 273 Human saliva NC_004567
L. plantarum JDM1 3 197 759 Human digestive tract NC_012984.1
L. plantarum NCU116 3 354 689 Pickles NZ_CP016071
L. plantarum ST-III 3 254 376 Pickles NC_014554.1
L. plantarum P9 3 016 749 Fermented sour porridge NZ_CP060172.1
L. plantarum Zhang-LL 2 952 218 Fermented rice products NZ_CP011769.1
L. plantarum Lp-6 3 101 304 Chinese koumiss NZ_CP060182.1
L. plantarum KLDS1.0391 2 886 607 Traditional fermented dairy products CP019348.1
L. plantarum ZJ316 3 203 964 Infant feces NC_020229.1
L. plantarum 10CH 3 311 056 Cheese NZ_CP023728.1
L. plantarum CAUH2 3 254 946 Sichuan pickles NZ_CP015126.1
L. plantarum EM 3 184 808 Korean kimchi NZ_CP037429.1
表选项

1.2 全基因组测序和组装

全基因组由上海美吉生物医药科技有限公司基于Illumina HiSeq和PacBio Sequel平台进行测序。细菌基因组扫描图是利用短序列组装软件SOAPdenovo2 (http://soap.genomics.org.cn/)对二代测序后的优化序列进行多个Kmer参数的拼接，得到最优的contigs组装结果，然后把reads比对到contig上，根据reads的配对端(paired-end)和重叠(overlap)关系，对组装结果进行局部组装和优化，形成scaffolds。细菌基因组完成图利用组装软件Unicycler[8]进行三代序列组装，组装过程中会借助Pilonjin软件进行序列校正，如果最终组装序列两端存在一定长度以上的overlap，则将序列成环并截去其中一端overlap序列[9]。最终可得到完整的染色体及质粒序列。

1.3 全基因组特征预测

利用Glimmer (http://ccb.jhu.edu/software/glimmer/index.shtml)、GeneMarkS、Prodigal软件对基因组中的编码序列(coding sequence, CDS)进行预测。扫描图组装结果默认使用Glimmer进行预测，完成图组装结果默认使用Glimmer预测染色体基因组，使用GeneMarkS预测质粒基因组[9]。利用tRNAscan-SE v2.0软件(http://trna.ucsc.edu/software/)对基因组中包含的tRNA进行预测，可以获得每个样本基因组中tRNA的核苷酸序列信息，反密码子信息及二级结构信息。利用Barrnap软件(https://github.com/tseemann/barrnap)对基因组中包含的rRNA进行预测，获得每个样本基因组中所有rRNA的种类、位置、序列信息[10]。最后采用CGView绘制基因组圈图。

1.4 构建系统发育树

对本研究中2株植物乳杆菌和12株参考菌株进行系统发育分析，利用MEGA X软件基于16S rRNA基因序列构建邻接(neighbor-joining, NJ)树，参数bootstrap值设为1 000。

1.5 平均核苷酸一致性(ANI)分析

平均核苷酸一致性(average nucleotide identity, ANI)是基于两两基因组之间所有直系同源核苷酸序列比较的一个平均值，主要用于在全基因组水平评估物种间的亲缘关系。本研究利用软件PYANIV0.2.10 (https://github.com/widdowquinn/pyani.git)计算L. plantarum PC06和L. plantarum PC108与参考菌株之间的相似性并绘制聚类热图。

1.6 基因组功能注释与分析

使用Diamond和blast2go软件分别通过COG (clusters of orthologous groups)、KEGG (Kyoto encyclopedia of genes and genomes)和CAZy (carbohydrate-active enzymes)数据库进行功能基因注释。使用VFDB (virulence factors database)数据库核心数据集(setA)和完整数据集(setB)中14个基础毒力因子类别对细菌所有基因组进行毒力因子预测，采用严格参数(相似性 > 80%，覆盖率 > 60%，E值≤1e-5)和一般参数(相似性 > 30%，覆盖率 > 70%，E值≤1e-5)[11]。使用CARD (comprehensive antibiotic resistance database, http://arpcard.mcmaster.ca)耐药基因标识符(RGI)软件(CARD比对条件设置为loose)[11]从细菌基因组序列中预测抗微生物耐药性(AMR)基因，提供由抗生素耐药性本体(ARO)分类标签和抗生素耐药基因(ARG)注释以初步评估菌株的安全性。此外，利用BLAST软件与参考菌株L. plantarum WCFS1进行对比，挖掘胃肠道耐受性相关基因；利用BAGEL4 (http://bagel.molgenrug.nl/index.php)在线分析软件对基因组序列进行产植物乳杆菌素相关基因簇的挖掘，综合分析益生相关基因以初步探究其益生潜力。

2 结果与分析 2.1 植物乳杆菌PC06和PC108全基因组概况

对测序样品L. plantarum PC06和PC108进行数据评估和序列组装，得到基因组组装结果(表 2)和基因组基本特征(表 3)。

表 2. 植物乳杆菌PC06与PC108基因组装结果详情 Table 2. Details of Lactiplantibacillus plantarum PC06 and PC108 gene assembly results
Characteristics L. plantarum PC06 L. plantarum PC108
Total reads No. 410 272 508 057
Total bases (bp) 4 648 952 898 5 873 062 770
Genome size (bp) 3 343 066 3 344 278
Largest (bp) 294 068 228 819
Average length (bp) 11 331.39 11 559.85
Clean Q20 (%) 96.59 96.70
Clean Q30 (%) 90.64 90.89
G+C (%) 44.33 44.50
表选项

表 3. 植物乳杆菌PC06和PC108基因组基本特征 Table 3. Basic characteristics of the Lactiplantibacillus plantarum PC06 and PC108 genomes
Genomic characteristics PC06 PC108
Number of CDS 3 161 3 197
Number of tRNA 73 69
Number of rRNA 16 16
Number of sRNA 48 49
表选项

测序结果显示，L. plantarum PC06全基因组包含1个片段长度为3 163 902 bp的环状染色体，其GC含量为44.68%，含6个完整质粒，其序列长度和GC含量分别为67 746 bp和38.30%、61 375 bp和39.23%、40 384 bp和36.68%、6 143 bp和36.82%、3 516 bp和37.32%、2 042 bp和37.95%。L. plantarum PC06预测有3 161个编码基因。

L. plantarum PC108全基因组包含1个片段长度为3 205 054 bp的环状染色体，其GC含量为44.67%，含5个完整质粒，其序列长度和GC含量分别为49 882 bp和41.79%、35 865 bp和39.79%、34 459 bp和41.37%、9 568 bp和36.95%、9 450 bp和35.95%。L. plantarum PC108预测有3 197个编码基因。

针对L. plantarum PC06和PC108组装的基因组序列，利用GCView对染色体(chromosome, Chr)进行绘制，包含全部基因、COG注释基因、prophage、CRISPR、ncRNA、GC含量以及GC skew结果。其图谱见图 1。

图 1 Lactiplantibacillus plantarum PC06与PC108基因组图谱 Figure 1 Genome mapping of Lactiplantibacillus plantarum PC06 and PC108. A: L. plantarum PC06. B: L. plantarum PC108. The information corresponding to the outer to inner circle is genome size identification, positive strand, gene information on negative strand, GC content, GC-Skew, and genome size identification.
图选项

2.2 系统发育树构建

为了探究L. plantarum PC06和PC108两株菌与其他具有益生特性的12株L. plantarum之间的关系，基于16S rRNA基因构建包含PC06、PC108以及12株参考菌株的系统发育树(图 2)。发现PC06与EM和CAUH2具有较近的亲缘关系，EM和CAUH2的来源均为泡菜。PC108与JDM1和ZJ316亲缘关系较近，而JDM1与ZJ316均来源于人体，从系统发育树中也可看出PC06与PC108亲缘关系相距较远，这意味着PC108可能往更适应人体环境的方向进化。

图 2 16S rRNA基因系统发育树分析 Figure 2 Phylogenetic tree analysis of 16S rRNA gene. Above the node is the self-spreading value; Below the branch is the length of the evolutionary branch, which represents the degree of change in the evolutionary branch, the shorter the difference means the smaller the difference and the closer the evolutionary distance.
图选项

2.3 平均核苷酸一致性(ANI)分析

ANI是基于两两基因组之间所有直系同源核苷酸/氨基酸序列比较的一个平均值，主要用于在全基因组水平评估物种间的亲缘关系，反映基因组之间进化距离关系，一般认为ANI≥95%为同种[12]。

本研究对14株植物乳杆菌全基因组序列进行ANI分析，并绘制了聚类热图(图 3)。结果显示全部菌株之间ANI均大于98.8%，可从基因组层面上确认PC06和PC108为植物乳杆菌，且不同植物乳杆菌基因组间存在差异。其中PC06和PC108之间ANI为99.05%；PC06和EM的ANI最高为99.84%；PC108和ZJ316的ANI最高为99.56%。

图 3 14株植物乳杆菌的ANI分析结果 Figure 3 Results of ANI analysis of 14 strains of Lactiplantibacillus plantarum.
图选项

2.4 基因组功能注释与分析 2.4.1 COG数据库注释

COG数据库可将蛋白功能划分为26类，每一类都是由直系同源序列构成。通过比对蛋白序列，将其注释到某一COG分类中，用以推测该序列的功能[13]。L. plantarum WCFS1是目前研究最深入的植物乳杆菌之一，有较高的参考价值，同时其具有较强的胆盐耐受能力，能在人体肠道存活和定殖[14]，具有较高的益生潜力。因此，以L. plantarum WCFS1作为参考菌株能一定程度上加深对L. plantarum PC06和PC108的了解。

L. plantarum PC06和PC108在COG数据库中分别得到了2 480个和2 508个功能基因注释，占CDS的78.46%和78.45%。注释结果显示(表 4)，L. plantarum PC06碳水化合物转运和代谢(G)相关基因得到了最多注释量，达到294个；而L. plantarum PC108转录(K)得到的基因注释量最多，有262个；同时一般功能预测(R)和氨基酸转运及代谢(E)相关功能蛋白也得到了较高的基因注释量，意味着这2株菌具有较强的利用碳水化合物和氨基酸的能力。除此之外，这2株菌还在细胞壁、生物膜和包膜生物发生中分别获得172个和169个基因注释量，可以辅证这2株菌可能具有较强的生物膜形成能力，并且防御机制的注释量分别为73个和87个，表明它们可能具有一定的抵抗外界环境潜力，但基因是否能表达需后续研究验证。

表 4. 三株植物乳杆菌COG分类详细统计表 Table 4. Detailed statistics of the COG classification of the three Lactiplantibacillus plantarum strains
COG code COG annotation PC06 PC108 WCFS1
A RNA processing and modification 0 0 0
B Chromatin structure and dynamics 0 0 0
C Energy production and conversion 122 115 114
D Cell cycle control, cell division, chromosome partitioning 46 45 26
E Amino acid transport and metabolism 223 227 267
F Nucleotide transport and metabolism 114 114 87
G Carbohydrate transport and metabolism 294 266 313
H Coenzyme transport and metabolism 134 140 91
I Lipid transport and metabolism 105 110 75
J Translation, ribosomal structure and biogenesis 213 218 155
K Transcription 268 268 305
L Replication, recombination and repair 132 131 157
M Cell wall/Membrane/Envelope biogenesis 172 169 156
N Cell motility 21 16 11
O Posttranslational modification, protein turnover, chaperones 105 105 62
P Inorganic ion transport and metabolism 126 133 169
Q Secondary metabolites biosynthesis, transport and catabolism 27 20 42
R General function prediction only 226 225 421
S Function unknown 128 124 224
T Signal transduction mechanisms 134 131 100
U Intracellular trafficking, secretion, and vesicular transport 22 20 26
V Defense mechanisms 73 87 60
W Extracellular structures 4 4 0
X Mobilome: prophages, transposons 78 117 0
Z Cytoskeleton 3 3 0
表选项

L. plantarum PC06和PC108在碳水化合物转运和代谢(G)、氨基酸转运及代谢(E)和转录(K)上得到了和L. plantarum WCFS1相近的基因注释量，意味着可能在这些方面有着和WCFS1相似的能力。与L. plantarum WCFS1相比，L. plantarum PC06和PC108关于前噬菌体和转座子(X)的功能蛋白得到了较多注释，这意味着PC06和PC108有获取外源基因片段的能力，而这能帮助细菌适应环境的变化或提高自身的生存竞争力，从而占据优势生态位。

2.4.2 KEGG数据库注释

KEGG数据库整合了多方面的数据，可用于系统分析基因表达产物的功能及代谢途径[15]。绘制的KEGG数据库注释结果见图 4，L. plantarumPC06和PC108分别有2 867个和2 822个基因在KEGG数据库中得到功能注释。在KEGG注释中PC06和PC108有相似的表现，大部分基因被注释到各个代谢通路，主要分布于全局和概述图谱，其次碳水化合物代谢以及氨基酸代谢，碳水化合物代谢通路包括各种糖类的代谢通路。环境信息处理中得到最多注释的是膜运输，主要集中于ABC转运蛋白、磷酸转移酶系统以及细菌分泌系统。值得一提的是，L. plantarum PC06和PC108在人类疾病方面得到注释的基因较少，主要是抗菌药物抗性。

图 4 Lactiplantibacillus plantarum PC06与PC108 KEGG功能注释图 Figure 4 Annotated diagram of the KEGG function of Lactiplantibacillus plantarum PC06 and PC108. A: L. plantarum PC06. B: L. plantarum PC108.
图选项

2.4.3 CAZy数据库注释

CAZy是关于合成或分解复杂碳水化合物和糖复合物酶类的专业数据库[16]。碳水化合物在很多生物学功能中具有重要地位，通过研究碳水化合物相关酶可以得到大量有意义的生物学信息。在CAZy注释结果中(图 5)，L. plantarum PC06和PC108在糖苷水解酶(glycoside hydrolases, GHs)上得到了最高的注释量，分别为43.27%和44.76%，其次是糖基转移酶(glycosyl transferases, GTS)。糖苷水解酶作用于含有2个及以上1, 4-α-d-葡萄糖基的多糖中的1, 4-α-d-葡萄糖苷键，多糖水解过程释放大量能量，进而支撑细菌的各种活动。糖基转移酶利用活化的供体形成糖苷键，将糖转移到特定的受体如蛋白质、脂质或其他聚糖上，形成新的聚合物从而参与各种生理过程[17]。碳水化合物酯酶(carbohydrate esterases, CEs)也得到了部分注释，能催化各种碳水化合物底物的去酯化，而辅助活性酶(auxiliary activity enzymes)包括一大类作用于碳水化合物的氧化还原活性酶[18]。

图 5 Lactiplantibacillus plantarum PC06与PC108 CAZy注释结果 Figure 5 Result of Lactiplantibacillus plantarum PC06 and PC108 CAZy annotations. A: L. plantarum PC06. B: L. plantarum PC108.
图选项

2.4.4 VFDB数据库注释

VFDB是全面的一体化的病原菌毒力因子信息管理网络数据库。来源于微生物，并对微生物自身侵染和引起特定宿主疾病具有促进作用的物质称为毒力因子，主要包括细菌毒素、调节细菌黏附作用的细胞表面蛋白、对细菌本身具有保护作用的蛋白、细胞表面碳水化合物以及具有细菌致病性的水解酶等[19]。

使用严格标准在VFDB数据库中比对后未发现毒力因子。使用一般标准时，L. plantarum PC06和PC108分别得到99个和97个有明确分类的基因注释量(图 6)，其中主要涉及铁离子摄取系统、细菌黏附与定殖和抗吞噬。铁是绝大多数细菌生存所必需的营养物质，2株菌均在铁离子摄取系统中得到了大量的基因注释，意味着PC06和PC108可能有着较强的铁离子吸收利用能力以提升其生存能力。L. plantarum PC06和PC108在抗吞噬作用上也得到了较多的注释量，抗吞噬的毒力因子能够帮助细菌对抗宿主免疫防御系统从而免于清除[20]。同时与定殖相关的基因也得到了较多注释，例如在L. plantarum PC06和PC108中找到了3种与黏附性有关的基因(VF0349、VF0327和VF0014)，而这些基因在细菌中广泛存在，同源基因比对结果发现植物乳杆菌WCFS1、NL42等菌株均有以上基因[14, 21]，而这些基因使菌株拥有与其他细胞黏附或让细胞相互黏附的能力，并非具有毒力效应。值得一提的是，在PC108中存在2个VF0014基因，位于不同的基因片段上，这可能会导致PC108的自聚性强于PC06，与团队之前的研究结果相吻合[22]。事实上这些基因在毒力因子数据库中被鉴定为毒力因子，因为它们也参与敌对/宿主环境中的致病细菌适应、存活或附着。然而，如果无其他发病机制，这些基因可以被认为是对细菌有益的，因为它们增加了细菌适应性，并且在需要细胞存活能力的情况下是可被接受的(例如增殖培养益生菌等)。

图 6 Lactiplantibacillus plantarum PC06与PC108 VFDB注释结果 Figure 6 Results of Lactiplantibacillus plantarum PC06 and PC108 VFDB annotations. A–J: Secretion system, toxin, invasion, adherence, serum resistance, stress protein, antiphagocytosis, regulation, magnesium uptake system, iron uptake system.
图选项

此外，在这些基因中检测出了毒力因子溶血素Ⅲ (VFG038904)和β-溶血素(CVF171)，事实上该基因在几种商业益生菌中也有报道，例如批准的公认安全(generally regarded as safe, GRAS)益生菌菌株植物乳杆菌299 V，中国广泛使用的商业益生菌植物乳杆菌JDM1，泡菜益生菌植物乳杆菌ST-III等，但团队前期研究证明L. plantarum PC06和PC108不具有溶血性[22]。

2.4.5 CARD数据库注释

CARD广泛包含来自各种生物体、基因组、质粒的抗生素抗性相关的参考基因，可用于指导环境、人体、动物菌群耐药组和抗生素抗性机制研究[23]。

L. plantarum PC06与PC108从CARD数据库中共计得到了154个和156个抗生素抗性基因注释。结果显示(图 7)，L. plantarum PC06与PC108对大环内酯类抗生素、四环素类抗生素和氟喹诺酮类抗生素具有较多的基因得到注释，说明它们对这3种抗生素可能具有较高的抗性，与团队之前的研究结果相同[22]。

图 7 Lactiplantibacillus plantarum PC06与PC108 CARD注释结果 Figure 7 Results of Lactiplantibacillus plantarum PC06 and PC108 CARD annotations. A: L. plantarum PC06. B: L. plantarum PC108.
图选项

通过分析L. plantarum PC06与PC108耐药基因所在位置发现，2株菌的质粒上均存在少量耐药基因(表 5)。其中PC06存在6种不同的耐药基因，分别分布在质粒A、质粒B和质粒C上，而PC108仅在质粒B和质粒C上存在3种耐药基因。从ARO描述来看，2株菌质粒上携带的耐药基因都是非植物乳杆菌特有的，可能来自环境中其他细菌。研究表明，大多数益生菌对抗生素的耐药性是固有抗性，是非转移性的，但有质粒的益生菌会存在与结合质粒相关的耐药性发生转移的可能性[24]。因此，由于L. plantarum PC06与PC108的质粒上存在少量耐药基因，理论上存在发生耐药基因转移的可能性，但确切关于菌株耐药性是否有转移风险以及其可能性大小的讨论需进一步研究论证。

表 5. Lactiplantibacillus plantarum PC06与PC108中位于质粒上的耐药基因统计 Table 5. Statistics of drug resistance genes located on plasmids in Lactiplantibacillus plantarum PC06 and PC108
Name Location ARO name ARO description Drug class Resistance mechanism
PC06 PlasmidA vanHM vanHM is a vanH variant in the vanM gene cluster Glycopeptide antibiotic Antibiotic target alteration
PC06 PlasmidA mepR MepR is an upstream repressor of MepA in Staphylococcus aureus. It is part of the mepRAB operon Glycylcycline; Tetracycline antibiotic Antibiotic efflux
PC06 PlasmidB Acinetobacter baumannii AbaF Expression of abaF in E. coli resulted in increased resistance to fosfomycin Fosfomycin Antibiotic efflux
PC06 PlasmidB cmrA cmrA is a transposon-encoded chloramphenicol exporter that is found in Rhodococcus rhodochrous Phenicol antibiotic Antibiotic efflux
PC06 PlasmidC vanHM vanHM is a vanH variant in the vanM gene cluster Glycopeptide antibiotic Antibiotic target alteration
PC06 PlasmidC macB MacB is an ATP-binding cassette (ABC) transporter that exports macrolides with 14- or 15- membered lactones. It forms an antibiotic efflux complex with MacA and TolC. macB corresponds to 1 locus in Pseudomonas aeruginosa PAO1 and 1 locus in Pseudomonas aeruginosa LESB58 Macrolide antibiotic Antibiotic efflux
PC108 PlasmidB lfrA lfrA is involved in the active efflux of quinolones and is found in Mycobacteroides abscessus Fluoroquinolone antibiotic antibiotic efflux
PC108 PlasmidB Acinetobacter baumannii AbaF Expression of abaF in E. coli resulted in increased resistance to fosfomycin Fosfomycin antibiotic efflux
PC108 PlasmidC tetA(58) TetA (58) is a tetracycline efflux pump described in Paenibacillus sp. LC231 Tetracycline antibiotic antibiotic efflux
表选项

2.5 胆盐水解酶相关基因

能在肠道内生存是益生菌在人体内发挥作用的重要指标之一。益生菌在通过胃后会到达十二指肠，在那里将会遇到各种压力条件，包括结合胆汁盐的存在。胆汁盐不仅能分散和吸收脂肪，还能作为表面活性剂破坏细胞膜的完整性，同时还能产生自由基降低细胞内的pH[25]，而胆汁盐能被胆盐水解酶(Bsh)解偶联并在结肠中重新吸收[26]。

L. plantarum WCFS1被证实有较好的胃肠道存活能力，有报道指出其在体外胃肠道存活模型中暴露1 h后，活细胞减少了3个数量级[27]。研究表明L. plantarum WCFS1中有4个被注释为胆盐水解酶的基因(Bsh1、Bsh2、Bsh3和Bsh4)[26]，通过基因目的片段查找发现L. plantarum PC06与PC108同样存在这4个胆盐水解酶基因，并且通过KEGG数据库注释到2个通路类别(表 6)。说明这2株菌可能具有与L. plantarum WCFS1相似的胃肠道存活能力，同时有报道指出胆盐水解酶与胆固醇清除能力有关[28]，可通过减少肠道胆固醇重吸收来降低甘油三酯和低密度脂蛋白胆固醇，有助于控制心血管疾病患者的高胆固醇[29]，但这2株菌的胆固醇降解能力和对于心血管疾病的帮助需进一步验证。

表 6. 两株植物乳杆菌胆盐水解酶(Bsh)相关基因和KEGG注释 Table 6. Bile salt hydrolase (Bsh)-related genes and KEGG description of two Lactiplantibacillus plantarum strains
Strains Gene Identify (%) EC KEGG path description
PC06 bsh1 98.6 3.5.1.24 Secondary bile acid biosynthesis; Primary bile acid biosynthesis
PC06 bsh2 98.6 3.5.1.11 Penicillin and cephalosporin biosynthesis
PC06 bsh3 98.6 3.5.1.11 Penicillin and cephalosporin biosynthesis
PC06 bsh4 98.6 3.5.1.11 Penicillin and cephalosporin biosynthesis
PC108 bsh1 98.6 3.5.1.24 Secondary bile acid biosynthesis; Primary bile acid biosynthesis
PC108 bsh2 98.6 3.5.1.11 Penicillin and cephalosporin biosynthesis
PC108 bsh3 98.6 3.5.1.11 Penicillin and cephalosporin biosynthesis
PC108 bsh4 98.6 3.5.1.11 Penicillin and cephalosporin biosynthesis
表选项

2.6 两株植物乳杆菌在细菌素合成基因簇比较结果

为了适应环境和生存，植物乳杆菌会对其他共生或竞争的微生物作出反应，这种反应称为群体感应，取决于细胞密度的基因表达，这对细菌的存活有着重要意义。细菌素是一种由细菌核糖体合成的具有抗菌活性的多肽，其中，植物乳杆菌素是由植物乳杆菌产生的细菌素，相关抗菌蛋白由pln基因簇编码[30]。

以L. plantarum WCFS1为例，其细菌素得到了充分的研究。该菌细菌素合成基因簇包含5个操纵子的pln基因座，plnABCD编码称为细菌素A (plnA)，操纵子plnEFI和plnJKLR用它们各自的免疫蛋白编码细菌素EF和JK，细菌素随后由操纵子plnGHSTUVWXY编码的ABC转运蛋白和辅助蛋白(PlnGH)转运和分泌[31-32]。

与L. plantarum WCFS1相比，L. plantarum PC06在plnJ和plnN之间插入了1个Na+/H+逆向转运蛋白(与L. plantarum ZJ316相似度为95.83%)，在plnJKLR操纵子中多了1个IS3家族转座酶和1个转座酶InsN，这可能是其他基因片段的转座子带来的，而COG数据库中得到较多转座子相关注释也能进一步佐证，但这对细菌素JK合成是否有影响尚未可知，同时还发现了1个包含67个氨基酸的未知功能假设蛋白，其功能有待进一步研究。L. plantarum PC108的细菌素合成簇则与L. plantarum WCFS1有较大不同，操纵子plnJKLR是相对保守的操纵子，但PC108缺失了plnJ和plnO (糖基转移酶GlyS)的相关基因片段，这可能会影响细菌素JK的合成。同时，PC108的操纵子plnJKLR上多出了1个包含61个氨基酸的核心肽，通过UniRef90比对发现这是1个具有双苷氨酸前导的IIb类细菌素lactobin A。有研究报道指出lactobin A首次被发现是由Lactobacillus amylovorus LMG P-13139产生，其抑菌谱较窄，对Lactobacillus delbrueckii subsp.最有效[33]。除此之外，PC108的编码蛋白plnG有一部分发生了氨基酸突变，突变区域与正常plnG的相似度为97.682%，由其中基因编码蛋白功能可知，此操纵子功能与细菌素的转运有关，这可能会影响细菌素的转运功能(图 8)。

图 8 三株植物乳杆菌细菌素合成基因簇 Figure 8 Cluster of bacteriocin synthesis genes in three strains of Lactobacillus plantarum.
图选项

2.7 两株植物乳杆菌益生特性相关基因

L. plantarum PC06和PC108的潜在益生特性在团队前期的研究中得到证实[22]，同时在全基因组分析数据中得到了支持。本研究中检测与酸耐受性相关的基因F0F1 ATP合成酶与胆盐抗性相关的胆盐水解酶(BSH)的基因(表 7)。益生菌在工业化生产过程中可能会遭受高/低温胁迫，暴露与高温环境中可诱导保守的热休克蛋白(HSPs)的表达，例如GroES、GroEL、GrpE、DnaK和DnaJ。在多种微生物中都发现了冷休克蛋白相关基因，并且这些基因与细菌在低温下存活的适应过程相关[34]，在L. plantarum PC06和PC108染色体上发现了CSP家族的冷休克蛋白基因。此外，L. plantarum PC06和PC108全基因组中还检测出过氧化氢酶相关基因、硫基过氧化酶相关基因和谷胱甘肽过氧化酶相关基因，这些酶可以防止氧化应激。

表 7. 植物乳杆菌PC06和PC108重要的益生相关基因信息 Table 7. Important genes of probiotic-related in Lactiplantibacillus plantarum PC06 and PC108
Strains Gene name Product EC Stress response
PC06/PC108 AtpA ATP synthase subunit alpha 7.1.2.2 Acid stress
PC06/PC108 AtpD ATP synthase subunit delta 7.1.2.2 Acid stress
PC06/PC108 AtpC ATP synthase epsilon chain 3.6.3.14 Acid stress
PC06/PC108 AtpG ATP synthase gamma chain 3.6.3.14 Acid stress
PC06/PC108 AtpH ATP synthase subunit delta 3.6.3.14 Acid stress
PC06/PC108 AtpB ATP synthase subunit a – Acid stress
PC06/PC108 AtpF ATP synthase subunit b 3.6.3.14 Acid stress
PC06/PC108 AtpE ATP synthase subunit c 3.6.3.14 Acid stress
PC06/PC108 Bsh1 Choloylglycine hydrolase 3.5.1.24 Bile tolerance
PC06/PC108 Bsh2 Choloylglycine hydrolase 3.5.1.11 Bile tolerance
PC06/PC108 Bsh3 Choloylglycine hydrolase 3.5.1.11 Bile tolerance
PC06/PC108 Bsh4 Choloylglycine hydrolase 3.5.1.11 Bile tolerance
PC06/PC108 GroES Heat shock protein 60 family co-chaperone – Temperature stress
PC06/PC108 GroEL Heat shock protein 60 family chaperone – Temperature stress
PC06/PC108 DnaJ Chaperone protein – Temperature stress
PC06/PC108 DnaK Chaperone protein – Temperature stress
PC06/PC108 CSP Cold shock protein of CSP family – Temperature stress
PC06/PC108 KatE Catalase KatE 1.11.1.6 Oxidative stress
PC06/PC108 Gpx Glutathione peroxidase 1.11.1.9 Oxidative stress
PC06/PC108 Tpx Thiol peroxidase, Tpx-type 1.11.1.15 Oxidative stress
PC06 YfeX DYP peroxidase 1.11.1.- Oxidative stress
PC108 YfeX Presumptive dye decolorization peroxidase (DyP) 1.11.1.- Oxidative stress
PC06/PC108 ArsC Arsenate reductase (thiorandoxin) 1.20.4.4 Oxidative stress
PC06/PC108 Trx Thiorededin family proteins – Oxidative stress
– indicates no EC number.
表选项

3 讨论与结论

四川泡菜是一个天然的微生物资源库，有大量有价值的微生物值得研究。本研究收集了四川多代泡菜母水样本，从中分离筛选出13株具有潜在益生特性的植物乳杆菌，其中包括L. plantarum PC06和PC108。在前期筛选阶段发现这2株菌表现出了较好的耐酸耐胆盐、胃肠液耐受性和自聚集共聚集能力，并且具有较好的安全性[22]。在此基础上对这2株菌进行的全基因组测序和基因组功能解析。通过系统发育树和ANI分析结果发现，L. plantarum PC06和PC108基因组存在差异，ANI为99.08%，在系统发育树中分支距离较远。

使用CAZy、KEGG和COG数据库注释结果表明，L. plantarum PC06和PC108得到了相近的基因注释数量，并且2株菌在碳水化合物利用、氨基酸利用和糖基转移酶等功能上得到了大量的注释，这意味着其可能具有较强的碳水化合物利用能力，为菌株在各种环境中的生存能力提供了基础；VFDB和CARD数据库比对结果表明，大部分得到注释的基因并非能够直接导致人体致病的相关基因，其中仅有的2种与溶血性相关的基因也在前期研究中证实并未表达，而PC06与PC108的质粒上带有极少量的耐药基因，理论上存在转移的可能性，需后续研究严重。除此之外，L. plantarum PC06和PC108均携带了4个与胆盐水解酶(Bsh)相关的基因片段，比对相似度均为98.6%，这意味着这2株菌具有耐受胃肠液的潜力。另外，我们还发现这2株菌均有完整的细菌素合成基因簇，其中L. plantarum PC06与L. plantarum WCFS1的细菌素合成基因比较相似，而L. plantarum PC108的缺失了合成plnJK的部分基因片段，但有一个IIb类细菌素lactobin A的相关基因嵌入到了操作子plnJKLR中，这意味着L. plantarum PC108产生的细菌素可能与PC06和WCFS1有差异，具体功能有待进一步研究。最后，我们汇总了L. plantarum PC06和PC108中与各种胁迫相关的基因，为其在人体内发挥益生特性提供了基础。

综上所述，本研究通过分离筛选与全基因组测序得到了2株泡菜源植物乳杆菌的全基因组信息，并且通过与L. plantarum WCFS1的比较分析，发现了与胆盐水解酶和与具有抗菌活性的细菌素合成基因簇，说明这2株菌能够作为具有潜在益生特性的益生菌菌株，具体益生功能还需进行进一步研究，如转录组学和蛋白质组学研究，将使我们能够阐明所鉴定的基因在细菌存活、宿主定殖和抗菌活性等的生理学重要性。本研究从基因层面揭示了L. plantarum PC06和PC108的潜在益生机理，为今后筛选优良益生菌菌株和基因层面上益生特性评价提供遗传学基础。

References
[1] la FATA G, WEBER P. Probiotics and the gut immune system: indirect regulation[J]. Probiotics and Antimicrobial Proteins, 2018, 10(1): 11-21 DOI:10.1007/s12602-017-9322-6.
[2] SIEZEN RJ, van HYLCKAMA VLIEG JET. Genomic diversity and versatility of Lactobacillus plantarum, a natural metabolic engineer[J]. Microbial Cell Factories, 2011, 10(Suppl 1): S3 DOI:10.1186/1475-2859-10-S1-S3.
[3] CHOI EA, CHANG HC. Cholesterol-lowering effects of a putative probiotic strain Lactobacillus plantarum EM isolated from kimchi[J]. LWT-Food Science and Technology, 2015, 62(1): 210-217 DOI:10.1016/j.lwt.2015.01.019.
[4] WANG Y, SHANG N, QIN YX, ZHANG Y, ZHANG JL, LI PL. The complete genome sequence of Lactobacillus plantarum LPL-1, a novel antibacterial probiotic producing class IIa bacteriocin[J]. Journal of Biotechnology, 2018, 266: 84-88 DOI:10.1016/j.jbiotec.2017.12.006.
[5] 刘超楠, 蔡熙姮, 张新宇, 刘晗璐. 植物乳杆菌基因组学研究进展[J]. 食品科技, 2020, 45(11): 1-7. DOI:10.13684/j.cnki.spkj.2020.11.002
LIU CN, CAI XH, ZHANG XY, LIU HL. Progress in genomics of Lactobacillus plantarum[J]. Food Science and Technology, 2020, 45(11): 1-7 (in Chinese).
[6] EVANOVICH E, de SOUZA MENDONÇA MATTOS PJ, GUERREIRO JF. Comparative genomic analysis of Lactobacillus plantarum: an overview[J]. International Journal of Genomics, 2019, 2019: 1-11.
[7] KIM E, CHANG HC, KIM HY. Complete genome sequence of Lactobacillus plantarum EM, A putative probiotic strain with the cholesterol-lowering effect and antimicrobial activity[J]. Current Microbiology, 2020, 77(8): 1871-1882 DOI:10.1007/s00284-020-02000-8.
[8] WICK RR, JUDD LM, GORRIE CL, HOLT KE. Unicycler: resolving bacterial genome assemblies from short and long sequencing reads[J]. PLoS Computational Biology, 2017, 13(6): e1005595 DOI:10.1371/journal.pcbi.1005595.
[9] GRISSA I, VERGNAUD G, POURCEL C. CRISPRFinder: a web tool to identify clustered regularly interspaced short palindromic repeats[J]. Nucleic Acids Research, 2007, 35(suppl_2): W52-W57.
[10] BUCHFINK B, XIE C, HUSON DH. Fast and sensitive protein alignment using DIAMOND[J]. Nature Methods, 2015, 12(1): 59-60 DOI:10.1038/nmeth.3176.
[11] CHOKESAJJAWATEE N, SANTIYANONT P, CHANTARASAKHA K, KOCHARIN K, THAMMARONGTHAM C, LERTAMPAIPORN S, VORAPREEDA T, SRISUK T, WONGSURAWAT T, JENJAROENPUN P, NOOKAEW I, VISESSANGUAN W. Safety assessment of a nham starter culture Lactobacillus plantarum BCC9546 via whole-genome analysis[J]. Scientific Reports, 2020, 10: 10241 DOI:10.1038/s41598-020-66857-2.
[12] SONG YQ. Study on population genetics and functional genomics of Lactobacillus delbrueckii subsp. bulgaricus[D]. Hohhot: Doctoral Dissertation of Inner Mongolia Agricultural University, 2018 (in Chinese).
宋宇琴. 德氏乳杆菌保加利亚亚种的群体遗传学和功能基因组学研究[D]. 呼和浩特: 内蒙古农业大学博士学位论文, 2018.
[13] TATUSOV RL, GALPERIN MY, NATALE DA, KOONIN EV. The COG database: a tool for genome-scale analysis of protein functions and evolution[J]. Nucleic Acids Research, 2000, 28(1): 33-36 DOI:10.1093/nar/28.1.33.
[14] van den NIEUWBOER M, van HEMERT S, CLAASSEN E, de VOS WM. Lactobacillus plantarum WCFS1 and its host interaction: a dozen years after the genome[J]. Microbial Biotechnology, 2016, 9(4): 452-465 DOI:10.1111/1751-7915.12368.
[15] OGATA H, GOTO S, SATO K, FUJIBUCHI W, BONO H, KANEHISA M. KEGG: Kyoto encyclopedia of genes and genomes[J]. Nucleic Acids Research, 1999, 27(1): 29-34 DOI:10.1093/nar/27.1.29.
[16] NGUYEN STC, FREUND HL, KASANJIAN J, BERLEMONT R. Function, distribution, and annotation of characterized cellulases, xylanases, and chitinases from CAZy[J]. Applied Microbiology and Biotechnology, 2018, 102(4): 1629-1637 DOI:10.1007/s00253-018-8778-y.
[17] GILFIX BM. Congenital disorders of glycosylation and the challenge of rare diseases[J]. Human Mutation, 2019, 40(8): 1010-1012.
[18] WARDMAN JF, BAINS RK, RAHFELD P, WITHERS SG. Carbohydrate-active enzymes (CAZymes) in the gut microbiome[J]. Nature Reviews Microbiology, 2022, 20(9): 542-556 DOI:10.1038/s41579-022-00712-1.
[19] CHEN LH, YANG J, YU J, YAO ZJ, SUN LL, SHEN Y, JIN Q. VFDB: a reference database for bacterial virulence factors[J]. Nucleic Acids Research, 2005, 33(suppl_1): D325-D328.
[20] 孙杰, 沈敏, 吴恬菲, 伦永志. 人肠源植物乳杆菌PUM1785全基因组测序分析与功能挖掘[J]. 微生物学杂志, 2021, 41(5): 69-79.
SUN J, SHEN M, WU TF, LUN YZ. Whole genome sequencing analysis and function mining of Lactobacillus plantarum PUM1785 from human intestine[J]. Journal of Microbiology, 2021, 41(5): 69-79 (in Chinese).
[21] ZHANG B, ZUO FL, YU R, ZENG Z, MA HQ, CHEN SW. Comparative genome-based identification of a cell wall-anchored protein from Lactobacillus plantarum increases adhesion of Lactococcus lactis to human epithelial cells[J]. Scientific Reports, 2015, 5: 14109 DOI:10.1038/srep14109.
[22] 史梅莓, 伍亚龙, 杨恺, 吕鹏军, 汪冬冬, 唐垚, 王勇, 张其圣. 四川泡菜中潜在益生性植物乳杆菌的筛选及安全性评价[J]. 食品工业科技, 2022, 43(22): 165-172.
SHI MM, WU YL, YANG K, LÜ PJ, WANG DD, TANG Y, WANG Y, ZHANG QS. Screening and safety evaluation of potential probiotic Lactobacillus plantarum in Sichuan Paocai[J]. Science and Technology of Food Industry, 2022, 43(22): 165-172 (in Chinese).
[23] ALCOCK BP, RAPHENYA AR, LAU TTY, TSANG KK, BOUCHARD M, EDALATMAND A, HUYNH W, NGUYEN AL V, CHENG AA, LIU SH, MIN SY, MIROSHNICHENKO A, TRAN HK, WERFALLI RE, NASIR JA, OLONI M, SPEICHER DJ, FLORESCU A, SINGH B, FALTYN M, et al. CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database[J]. Nucleic Acids Research, 2020, 48(D1): D517-D525.
[24] DEGIROLAMO C, RAINALDI S, BOVENGA F, MURZILLI S, MOSCHETTA A. Microbiota modification with probiotics induces hepatic bile acid synthesis via downregulation of the fxr-Fgf15 axis in mice[J]. Cell Reports, 2014, 7(1): 12-18 DOI:10.1016/j.celrep.2014.02.032.
[25] BRON PA, MARCO M, HOFFER SM, van MULLEKOM E, de VOS WM, KLEEREBEZEM M. Genetic characterization of the bile salt response in Lactobacillus plantarum and analysis of responsive promoters in vitro and in situ in the gastrointestinal tract[J]. Journal of Bacteriology, 2004, 186(23): 7829-7835 DOI:10.1128/JB.186.23.7829-7835.2004.
[26] LAMBERT JM, BONGERS RS, de VOS WM, KLEEREBEZEM M. Functional analysis of four bile salt hydrolase and penicillin acylase family members in Lactobacillus plantarum WCFS1[J]. Applied and Environmental Microbiology, 2008, 74(15): 4719-4726 DOI:10.1128/AEM.00137-08.
[27] van BOKHORST-VAN de VEEN H, LEE IC, MARCO ML, WELS M, BRON PA, KLEEREBEZEM M. Modulation of Lactobacillus plantarum gastrointestinal robustness by fermentation conditions enables identification of bacterial robustness markers[J]. PLoS One, 2012, 7(7): e39053 DOI:10.1371/journal.pone.0039053.
[28] KUMAR R, GROVER S, KAUSHIK JK, BATISH VK. IS30-related transposon mediated insertional inactivation of bile salt hydrolase (bsh1) gene of Lactobacillus plantarum strain Lp20[J]. Microbiological Research, 2014, 169(7/8): 553-560.
[29] COSTABILE A, BUTTARAZZI I, KOLIDA S, QUERCIA S, BALDINI J, SWANN JR, BRIGIDI P, GIBSON GR. An in vivo assessment of the cholesterol-lowering efficacy of Lactobacillus plantarum ECGC 13110402 in normal to mildly hypercholesterolaemic adults[J]. PLoS One, 2017, 12(12): e0187964 DOI:10.1371/journal.pone.0187964.
[30] SÁENZ Y, ROJO-BEZARES B, NAVARRO L, DÍEZ L, SOMALO S, ZARAZAGA M, RUIZ-LARREA F, TORRES C. Genetic diversity of the pln locus among oenological Lactobacillus plantarum strains[J]. International Journal of Food Microbiology, 2009, 134(3): 176-183 DOI:10.1016/j.ijfoodmicro.2009.06.004.
[31] STURME MHJ, FRANCKE C, SIEZEN RJ, de VOS WM, KLEEREBEZEM M. Making sense of quorum sensing in lactobacilli: a special focus on Lactobacillus plantarum WCFS1[J]. Microbiology (Reading, England), 2007, 153(Pt 12): 3939-3947.
[32] DIEP DB, STRAUME D, KJOS M, TORRES C, NES IF. An overview of the mosaic bacteriocin pln loci from Lactobacillus plantarum[J]. Peptides, 2009, 30(8): 1562-1574 DOI:10.1016/j.peptides.2009.05.014.
[33] de VUYST L, AVONTS L, NEYSENS P, HOSTE B, VANCANNEYT M, SWINGS J, CALLEWAERT R. The lactobin A and amylovorin L471 encoding genes are identical, and their distribution seems to be restricted to the species Lactobacillus amylovorus that is of interest for cereal fermentations[J]. International Journal of Food Microbiology, 2004, 90(1): 93-106 DOI:10.1016/S0168-1605(03)00298-8.
[34] KIM WS, KHUNAJAKR N, REN J, DUNN NW. Conservation of the major cold shock protein in lactic acid bacteria[J]. Current Microbiology, 1998, 37(5): 333-336 DOI:10.1007/s002849900387.
全基因组测序揭示两株泡菜源植物乳杆菌基因型差异和潜在益生特性
伍亚龙 , 杨恺 , 史梅莓 , 吕鹏军 , 王勇 , 张其圣