嗜热链球菌KLDS3.1012胞外多糖合成途径的基因组学及表型特征分析

李柏良1,赵 莉1,王成凤1,靳 妲1,丁秀云1,2,刘 飞1,*,霍贵成1

(1.东北农业大学 乳品科学教育部重点实验室,黑龙江 哈尔滨 150030;2.广州基迪奥生物科技有限公司,广东 广州 510000)

摘 要: 在分子水平上挖掘嗜热链球菌KLDS3.1012合成胞外多糖的途径并且通过表型特征进行验证。首先,基于Illumina HiSeq与Illumina MiSeq测序平台对菌株KLDS3.1012进行基因组测序并构建基因组图谱;随后,从糖代谢、糖核苷酸合成、胞外多糖基因簇等方面进行生物信息学分析;最后,利用API 50CH测定菌株KLDS3.1012的糖发酵情况及采用高效离子交换色谱法测定胞外多糖的单糖组成。生物信息学分析结果表明菌株KLDS3.1012具有半乳糖、葡萄糖、果糖、甘露糖、乳糖和蔗糖转运系统和4 种糖核苷酸(UDP-葡萄糖、dTDP-鼠李糖、UDP-半乳糖和UDP-N-乙酰葡糖胺)合成的相关基因及1 个胞外多糖合成基因簇。表型特征分析结果表明菌株KLDS3.1012可以利用以上6 种碳源并能形成由鼠李糖、半乳糖和葡萄糖组成的胞外多糖。本研究为分析该菌株合成胞外多糖的遗传基础与胞外多糖结构的关系提供理论依据,并对将其开发为发酵剂具有一定指导意义。

关键词: 嗜热链球菌;基因组测序;生物信息分析;糖代谢;胞外多糖

嗜热链球菌是非致病性乳酸菌,常与德氏乳杆菌保加利亚亚种组合使用,作为酸奶的发酵剂[1-2]。部分嗜热链球菌可以产生胞外多糖(exopolysaccharides,EPS),EPS分泌到培养基中,其不仅可以改善发酵乳制品的流变学和感官特性[3-4],还具有诸多的益生特性,如抗氧化性和抗肿瘤活性、降低血压和血糖、改善肠道菌群、抑制致病菌和调节免疫系统等[5-9]。因此,能够产EPS的嗜热链球菌菌株在发酵方面具有更大的价值。

EPS的合成过程分为前体物质糖核苷酸的合成及EPS基因簇合成EPS两个阶段。单糖、双糖作为微生物主要利用的碳源,经转运系统转运至胞内,糖类物质的转运方式一般有3 种:磷酸烯醇式丙酮酸依赖磷酸转移酶系统(phosphotransferase system,PTS)、ABC型糖转运系统以及糖渗透[10]。其中PTS为主要的转运方式,由负责糖类物质磷酸化的可溶性磷酸载体蛋白HPr、能量耦合蛋白酶EI以及膜结合的具有糖特异性的通透酶EII组成。通常EII由3 个功能性的亚单位EIIA、EIIB与EIIC构成,转运甘露糖的PTS则由EIIA、EIIB、EIIC与EIID组成。根据序列及结构的特征,EII可以分为乳糖家族、葡萄糖家族、甘露糖家族及甘露醇家族。在糖转运过程中,如果某一PTS中的EII因缺失亚基而不完整,是不可以通过另一个非同家族的PTS的EII替代其功能。转运到细胞内的糖经过一系列的酶促反应,最后形成糖核苷酸,作为合成EPS的活性前体[11]

参与指导EPS重复单元合成、链长、聚合与输出的基因常成簇存在,称之为EPS基因簇。异型多糖重复单元的形成通过引导糖基转移酶转移相应的糖核苷酸到异戊二烯醇磷酸酯(脂载体)上,在其他的糖基转移酶作用下,高度特异性地转移相应的糖核苷酸形成重复单元的主、侧链。已经形成的重复单元经翻转酶从胞内输出至细胞外膜,随后经聚合酶聚合成完整的EPS,输出至细胞膜外[12-13]

嗜热链球菌KLDS3.1012是从中国新疆传统发酵酸奶中分离出来的,由东北农业大学乳品科学教育部重点实验室保藏。前期实验研究发现该菌株具产黏性能强及所产酸奶品质优良等特点,适合于制备高品质直投式发酵剂[14]。为更加深入分析该菌株合成ESP的遗传基础,本研究首先通过二代测序技术对该菌株进行基因组测序,基于生物信息学分析该菌株的糖代谢能力及EPS基因簇,同时利用API 50CH测定该菌株的糖发酵情况,并采用高效离子交换色谱法测定EPS的单糖组成。为研究该菌株EPS遗传基础与EPS结构之间关系及后续更加合理地应用该菌株提供了理论依据。

1 材料与方法

1.1 材料与试剂

嗜热链球菌KLDS3.1012由东北农业大学乳品科学教育部重点实验室工业微生物菌种保藏中心(KLDSDICC)提供,且通过16S rRNA基因测序鉴定。

M17肉汤培养基 青岛高科园海博生物技术有限公司;细菌基因组提取试剂盒 北京天根生物技术有限公司;API 50CH 法国梅里埃公司;DEAE-Sepharose Fast Flow和Sepharose CL-6B填料 美国GE医疗生命科学公司;无水乙醇 天津市天力化学试剂有限公司;三氯乙酸、三氟乙酸 天津市大茂化学试剂厂;重蒸苯酚北京博奥拓达科技有限公司;浓硫酸 哈尔滨理工化学试剂有限公司;其他试剂均为国产分析纯。

1.2 仪器与设备

LDZF-50KB-II立式蒸汽灭菌器 上海申安医疗器械厂;CJ-2D超净工作台 天津泰斯特仪器有限公司;DHP-927型电热恒温培养箱 上海一恒科技有限公司;透析袋 北京索莱宝生物科技有限公司;蛋白纯化仪美国GE医疗生命科学公司。

1.3 方法

1.3.1 菌株的活化及基因组提取

将甘油保藏的菌株KLDS3.1012以2%的体积分数接种于M17液体培养基,42 ℃培养24 h,连续转接两次,第3次培养12 h后用于基因组提取。按照细菌基因组提取试剂盒说明书的步骤提取菌株KLDS3.1012的基因组。

1.3.2 全基因组测序及组装

利用Illumina HiSeq(2×100 bp配对末端文库)和Illumina MiSeq(2×250 bp配对末端文库)平台组合测序技术对菌株KLDS3.1012进行基因组测序。通过进行质量控制和数据过滤,利用SOAPdenovo 2.0对所有有效数据(Clean Data)进行组装[15]

1.3.3 基因组注释

采用NCBI原核生物基因组注释流程(prokaryotic genome annotation pipeline,PGAP,https://www.ncbi.nlm.nih.gov/books/NBK174280/)进行基因预测及注释。基因组中蛋白序列通过WebMGA网站进行COG(Cluster of Orthologous Group)注释[16]

1.3.4 生物信息学分析

利用CGView Server构建基因组图谱[17]。糖代谢途径参照KEGG通路数据库进行分析[18]。ISfinder用于鉴定基因组序列中的转座酶[19]。使用ACT(Artemis Comparison Tool)软件对菌株KLDS3.1012与嗜热链球菌ND03的EPS基因簇进行分析[20]

1.3.5 表型特征的测定

采用API CHL培养基和API 50CH测试条按照试剂盒说明书测定菌株KLDS3.1012的碳水化合物利用情况。按照邵丽[21]所述的实验步骤,提取菌株KLDS3.1012所产的EPS。参考Ren Wei等[22]的方法利用DEAESepharose Fast Flow和Sepharose CL-6B对EPS进行纯化。参照邵丽[21]的方法,使用高效离子交换色谱(high performance anion exchange chromatography,HPAEC)法进行单糖组成的测定。

1.4 数据处理

利用HPAEC测定各标准品的保留时间和峰面积,以各浓度为横坐标,峰面积为纵坐标,绘制标准曲线。对样品中各单糖的峰面积进行积分,根据标准曲线计算相应的浓度。

2 结果与分析

2.1 基因组的基本特征及注释

图1 嗜热链球菌KLDS3.1012的基因组圈图
Fig. 1 Circular genome map ofS. thermophilus KLDS3.1012

由外至内依次为,第1圈表示正链的ncRNA分布,包括tRNA,rRNA及sRNA三种类型;第2圈表示正链的COG注释基因分布,以不同的颜色区分;第3圈表示负链的COG注释基因分布;第4圈表示负链的ncRNA分布;第5圈(黑色)表示GC含量,以平均GC为基准线,向外突出的表示高于均值,向内突出的表示低于均值;第6圈为GC skew值,紫色表示小于0,绿色表示大于0。

菌株KLDS3.1012的基因组序列已提交到GenBank数据库,登录号为:NZ_LHSK01000000。通过CGView Server构建其基因组圈图(图1)并利用PGAP流程进行注释,一共预测出1 992 个基因,其中包括1 705 个蛋白编码基因,213 个假基因,15 个rRNA,55 个tRNA和4 个ncRNA。

菌株KLDS3.1012的亲缘关系报告(表1)显示,菌株KLDS3.1012与嗜热链球菌ND03在基因组序列最为相似(99.572 2%),亲缘关系最近,而与嗜热链球菌M17PTZA496的基因组序列差异较大,对称一致性仅为83.842 8%。

表1 嗜热链球菌亲缘关系报告
Table 1 Genome neighbor report of S. thermophilus strains

含有缺口的一致性/%菌株 大小/Mb 对称一致性/%ND03 1.831 95 99.572 2 99.970 3 APC151 1.839 13 99.385 1 99.974 6 MN-ZLW-002 1.848 52 95.613 1 99.586 8 MN-BM-A01 1.876 52 94.905 0 99.574 4 TH1436 1.782 13 93.710 1 99.448 7 TH982 1.789 31 93.559 6 99.391 9 SMQ-301 1.861 79 92.786 3 99.350 4 LMD-9 1.864 18 92.738 5 99.321 9 ASCC 1275 1.845 50 92.704 6 99.421 8 MN-BM-A02 1.850 43 92.584 6 99.422 0 TH985 1.844 43 92.528 7 99.305 1 KLDS SM 1.856 79 92.436 8 99.432 9 ND07 1.869 51 92.120 3 99.431 9 ST3 1.865 06 92.009 0 99.277 5 ACA-DC 2 1.731 84 91.747 0 99.095 9 LMG 18311 1.796 85 91.295 5 99.043 9 S9 1.787 44 90.634 4 99.079 9 CS8 1.791 66 90.553 4 99.027 0 CNRZ1066 1.796 23 90.435 2 99.021 3 TH1435 1.753 46 90.378 7 99.318 1 1F8CT 1.746 94 90.237 3 98.936 4 KLDS 3.1003 1.899 96 89.936 0 99.007 3 MTH17CL396 1.828 22 89.788 0 99.108 9 JIM 8232 1.929 90 89.189 9 99.251 7 B59671 1.821 17 88.903 6 98.992 4 TH1477 1.889 56 87.512 2 99.008 2 M17PTZA496 2.065 46 83.842 8 98.696 9

2.2 糖代谢能力分析

2.2.1 糖转运系统

菌株KLDS3.1012不仅含有HPr(AKL23_RS05975)和EI(AKL23_RS05970),还含有10 个糖特异性通透酶(EII)的基因,其中包括蔗糖通透酶EII的基因、PTS果糖通透酶EII、PTS甘露糖通透酶EII和PTS葡萄糖通透酶EII(表2)。然而,PTS葡萄糖通透酶EII是由假基因编码的,而另外2 个EII基因也不完整,因此不能通过PTS转运葡萄糖。在菌株KLDS3.1012的基因组中鉴定到乳糖/半乳糖渗透酶(AKL23_RS06545)和葡萄糖渗透酶(AKL23_RS04150)。由于编码ABC型糖转运系统的基因不完整,因此ABC型糖转运系统不能将糖转运到菌株KLDS3.1012的细胞质中。综上所述,从基因组分析,菌株KLDS3.1012可以利用蔗糖、果糖、甘露糖、葡萄糖、半乳糖和乳糖。

表2 嗜热链球菌KLDS3.1012的PTS相关基因
Table 2 Putative genes related to phosphotransferase system (PTS) in S. thermophilus KLDS3.1012

编码蛋白 基因座PTS葡萄糖转运体亚基EIIA AKL23_RS00985 PTS葡萄糖转运体亚基EIIABC AKL23_RS09850(假基因)PTS葡萄糖转运体亚基EIIA AKL23_RS00990 PTS甘露糖转运体附属蛋白ManO AKL23_RS01625 PTS甘露糖转运体亚基EIIAB AKL23_RS01640 PTS甘露糖转运体亚基EIIC AKL23_RS01635 PTS甘露糖转运体亚基EIID AKL23_RS01630 PTS果糖转运体亚基EIIABC AKL23_RS02075 PTS果糖转运体亚基EIIA AKL23_RS06370 PTS蔗糖转运体亚基EIIABC AKL23_RS10585

2.2.2 糖核苷酸的合成途径

被转运到细胞质中的糖会通过各种代谢途径生成糖核苷酸。如图2所示,菌株KLDS3.1012可以通过β-半乳糖苷酶(AKL23_RS06540)水解乳糖生成葡萄糖和半乳糖,一方面,半乳糖通过UDP-半乳糖的途径(AKL23_RS06565和AKL23_RS06560)转化生成UDP-半乳糖和Leloir途径转化生成葡萄糖-1-磷酸。另一方面,葡萄糖被葡萄糖激酶(AKL23_RS03575)磷酸化生成葡萄糖-6-磷酸,再由磷酸变位酶(AKL23_RS03880)进一步形成葡萄糖-1-磷酸。葡萄-1-磷酸作为重要的中间产物,它可以通过UDP-葡萄糖焦磷酸化酶(AKL23_RS08610)转化成UDP-葡萄糖或通过一系列酶促反应(AKL23_RS05870、AKL23_RS05860、AKL23_RS05865、AKL23_RS03055和AKL23_RS06935)生成dTDP-鼠李糖。编码UDP-葡萄糖4-差向异构酶的基因具有3 个拷贝(AKL23_RS05370、AKL23_RS06390和AKL23_RS06555),这可能与UDP-葡萄糖和UDP-半乳糖之间的相互转换相关。甘露糖、果糖和蔗糖转运到细胞质中后分别被磷酸化形成甘露糖-6-磷酸、果糖-6-磷酸和蔗糖-6-磷酸,三者之间可以相互转化(AKL23_RS08125、AKL23_RS08110和AKL23_RS10365)。随后,果糖-6-磷酸由AKL23_RS04205、AKL23_RS05900和AKL23_RS02810编码的一系列酶转化成UDP-N-乙酰葡糖胺。总之,从基因组分析,菌株KLDS3.1012能够形成4 种糖核苷酸,即UDP-葡萄糖、dTDP-鼠李糖、UDP-半乳糖和UDP-乙酰葡糖胺。这些糖核苷酸可以形成EPS的重复单元,再通过EPS基因簇形成EPS。

图2 嗜热链球菌KLDS3.1012糖核苷酸的合成
Fig. 2 Synthesis of sugar nucleotides inS. thermophilus KLDS 3.1012

2.3 EPS基因簇的生物信息学分析

如图3所示,菌株KLDS3.1012存在1 个EPS合成基因簇(图4),由35 个基因组成。在该基因簇中的5'末端存在1 个嘌呤核苷磷酸化酶的编码基因deoD(AKL23_RS05125),该基因存在于大部分EPS基因簇附近[23-24]。基因AKL23_RS05120、AKL23_RS05115、AKL23_RS05110和AKL23_RS05105分别负责编码EPS的合成调节、链长决定和输出相关的保守基因。AKL23_RS05100编码的起始糖基转移酶,可以将UDP-半乳糖转移到脂质载体上。负责编码其他糖基转移酶的基因分别是AKL23_RS05095(假基因)、AKL23_RS05070、AKL23_RS05025(UDP-N-乙酰葡糖胺转移酶),AKL23_RS05020、AKL23_RS05015和AKL23_RS05000,这些糖基转移酶编码基因对决定EPS的单糖组成起到关键的作用。AKL23_RS05055编码的翻转酶具有将重复单位输出的功能。AKL23_RS004985编码的orf14.9与生长功能相关,该基因也存在于其他嗜热链球菌的EPS基因簇中[25]。AKL23_RS004980和AKL23_RS004960分别编码磷酸酶和通透酶。此外,该基因簇中还存在3 个磷酸甘油酸变位酶的编码基因(AKL23_RS004975、AKL23_RS004970和AKL23_RS004965)。

图3 嗜热链球菌KLDS3.1012的EPS基因簇
Fig. 3 Exopolysaccharides gene cluster in S. thermophilus KLDS3.1012

从5′(左)到3′(右),编码不同产物的基因用不同的图案表示。基因位点(AKL23_RS)与长度(bp)标记于基因上下。

前面的基因组亲缘关系研究表明,菌株KLDS3.1012与嗜热链球菌ND03的基因组序列亲缘关系最近,高达99.57%。通过ACT工具对2 株菌的EPS基因簇序列进行比对分析,由图4可知,菌株KLDS3.1012 EPS基因簇与嗜热链球菌ND03的EPS基因簇具有很高的同源性,仅存在7 个倒置区域,而且中间的空白区域是由于菌株KLDS3.1012的gap造成的。

图4 嗜热链球菌KLDS3.1012 EPS基因簇的比对分析
Fig. 4 Alignment of exopolysaccharides gene cluster in S. thermophilus KLDS3.1012

红色连接表明序列相似度从高(亮红)到低(粉色);蓝色连接表明序列发生倒转。

2.4 碳水化合物代谢能力

API 50CH碳水化合物代谢结果(图5)表明,菌株KLDS3.1012可以代谢半乳糖、葡萄糖、果糖、甘露糖、乳糖和蔗糖。虽然与其他嗜热链球菌碳水化合物代谢能力存在差异,但是菌株KLDS3.1012在碳水化合物利用方面,遗传分析与表型结果是相一致的。

图5 嗜热链球菌KLDS3.1012的API 50CH发酵结果
Fig. 5 Results of API 50CH test strip onS. thermophilus KLDS3.1012

试纸条呈黄色代表菌株可以利用该种糖,试纸条呈紫色代表菌株不可以利用该种糖,其中第10、11、12、13、29、31管分别代表半乳糖、葡萄糖、果糖、甘露糖、乳糖、蔗糖。

2.5 EPS的单糖组成

菌株KLDS3.1012所产的EPS通过DEAE-Sepharose Fast Flow和Sepharose CL-6B进行两步纯化,随后经三氟乙酸水解,最后用高效离子色谱测定EPS的单糖组成。由图6可知,菌株KLDS3.1012所产的EPS由鼠李糖、半乳糖和葡萄糖以2.6∶1.7∶1的物质的量比组成。在单糖组成中没有发现UDP-N-乙酰葡糖胺,然而,生物信息学分析表明菌株KLDS3.1012可以形成UDP-N-乙酰葡糖胺并且拥有UDP-N-乙酰葡糖胺转移酶,这可能是由于相关基因表达水平的不足造成的,还需要进一步的研究验证这一点。

图6 嗜热链球菌KLDS3.1012的EPS单糖离子色谱图
Fig. 6 HPAEC profile of monosaccharide composition of EPS from S. thermophilusKLDS3.1012

3 讨 论

嗜热链球菌作为重要的发酵剂,从功能上看,嗜热链球菌在治疗乳糖不耐症、抗氧化、改善肠道免疫反应和缓解一些癌症方面起着至关重要的作用[26]。从代谢能力上看,嗜热链球菌不仅能够代谢碳水化合物形成乳酸,酸化乳制品,还能形成糖核苷酸产生EPS[27-28]。嗜热链球菌所产的EPS具有结构复杂和功能多样的特点,因此作为近年来的持续热点问题。随着基因组测序技术的发展,研究人员可以利用生物信息学的手段充分挖掘嗜热链球菌所合成的EPS的相关分子元件,为揭示嗜热链球菌合成EPS的机理及调控提供基础。Bolotin等[29]研究了嗜热链球菌CNRZ1066和LMG13811的糖代谢相关遗传基础;Sun Zhihong等[30]完成了ND03的全基因组测序并对其独特的EPS基因簇进行分析;Wu Qinglong等[23]通过嗜热链球菌ASCC 1275的EPS基因簇分析,发现该菌株高产EPS的相关基因;Bai Ying等[31]报道了嗜热链球菌MN-BM-A01的EPS产量和相关基因的分析;Vendramin等[32]对8 株已完成基因组测序的嗜热链球菌epsA基因的启动子进行了比较基因组学分析。本研究不仅从生物信息学层面分析了菌株KLDS3.1012的糖代谢能力及EPS基因簇,同时还利用API 50CH测定该菌株的糖发酵情况及采用高效离子交换色谱法测定EPS的单糖组成,为研究EPS的基因与结构奠定了基础。

虽然嗜热链球菌在糖代谢水平上相对保守,然而不同的嗜热链球菌菌株的糖代谢情况并不一致,嗜热链球菌基本都可以利用葡萄糖、乳糖和果糖[33],而半乳糖、甘露糖、蔗糖、麦芽糖、蜜二糖和棉子糖的利用却有菌株特异性[34-35],如嗜热链球菌AR333可以利用葡萄糖、乳糖和蔗糖[22],嗜热链球菌MTCC 5461可以利用17 种碳水化合物[36],而本研究中的菌株KLDS3.1012可以利用半乳糖、葡萄糖、果糖、甘露糖、乳糖和蔗糖。能够利用半乳糖是工业发酵菌株的优良性状,可以避免半乳糖的过量积累而影响发酵乳制品的质量。此外,利用半乳糖可以提高EPS的产量,因为半乳糖可以形成糖核苷酸。嗜热链球菌所产EPS的单糖组成及比例也都不一样,AR333和GST-6所产的EPS都由半乳糖和葡萄糖组成,但是物质的量比差异较大[5,22],DGCC 7785所产的EPS由半乳糖、葡萄糖和N-乙酰半乳糖胺组成[37],而本研究中的KLDS3.1012所产的EPS由鼠李糖、半乳糖和葡萄糖,与嗜热链球菌DGCC 7698所产的EPS具有相似性[37]

糖核苷酸作为合成EPS的活性前体,是合成EPS的关键底物[11]。在分析菌株KLDS3.1012的糖核苷酸合成途径中发现,在其基因组中,催化UDP-葡萄糖与UDP-半乳糖相互转变的UDP-葡萄糖4-差向异构酶的编码基因高达3 个拷贝,利于二者的灵活转变,以满足形成重复单元过程中对不同糖核苷酸的需求。嗜热链球菌的EPS基因簇中参与指导EPS重复单元合成、链长、聚合与输出的基因相对保守,而糖基转移酶差别甚大,其具有多态性,目前,在51 株嗜热链球菌中存在118 个糖基转移酶,不同的糖基转移酶可以转移不同的糖核苷酸形成EPS的重复单元[38],因此糖基转移酶对EPS的结构起到决定性作用。本研究中菌株KLDS3.1012的糖基转移酶中只鉴定出半乳糖核苷酸转移酶和UDP-N-乙酰葡糖胺转移酶,此外,该菌株的EPS基因簇与嗜热链球菌ND03的EPS基因簇具有很高的同源性。

4 结 论

从基因组水平分析,菌株KLDS3.1012具有半乳糖、葡萄糖、果糖、甘露糖、乳糖和蔗糖转运系统和4 种糖核苷酸(UDP-葡萄糖、dTDP-鼠李糖、UDP-半乳糖和UDP-N-乙酰葡糖胺)合成的相关基因及1 个胞外多糖合成基因簇。表型测定结果表明菌株KLDS3.1012可以利用以上6 种碳源并能形成由鼠李糖、半乳糖和葡萄糖组成的EPS。

参考文献:

[1] LIU M, SIEZEN R J, NAUTA A. In Silico prediction of horizontal gene transfer events in Lactobacillus bulgaricus and Streptococcus thermophilus reveals protocooperation in yogurt manufacturing[J].Applied and Environmental Microbiology, 2009, 75(12): 4120-4129.DOI:10.1128/AEM.02898-08.

[2] SETTACHAIMONGKON S, NOUT M R, FERNANDES E C A, et al.Influence of different proteolytic strains of Streptococcus thermophilus in co-culture with Lactobacillus delbrueckii subsp.bulgaricus on the metabolite profile of set-yoghurt[J]. International Journal of Food Microbiology, 2014, 177: 29-36. DOI:10.1016/j.ijfoodmicro.2014.02.008.

[3] PURWANDARI U, SHAH N P, VASILJEVIC T. Effects of exopolysaccharide-producing strains of Streptococcus thermophilus on technological and rheological properties of set-type yoghurt[J].International Dairy Journal, 2007, 17(11): 1344-1352. DOI:10.1016/j.idairyj.2007.01.018.

[4] MENDE S, ROHM H, JAROS D. Influence of exopolysaccharides on the structure, texture, stability and sensory properties of yoghurt and related products[J]. International Dairy Journal, 2016, 52: 57-71.DOI:10.1016/j.idairyj.2015.08.002.

[5] ZHANG J, CAO Y, WANG J, et al. Physicochemical characteristics and bioactivities of the exopolysaccharide and its sulphated polymer from Streptococcus thermophilus GST-6[J]. Carbohydrate Polymers,2016, 146: 368-375. DOI:10.1016/j.carbpol.2016.03.063.

[6] PATTEN D A, LEIVERS S, CHADHA M J, et al. The structure and immunomodulatory activity on intestinal epithelial cells of the EPSs isolated from Lactobacillus helveticus sp. Rosyjski and Lactobacillus acidophilus sp. 5e2[J]. Carbohydrate Research, 2014, 384: 119-127.DOI:10.1016/j.carres.2013.12.008.

[7] MAEDA H, ZHU X, OMURA K, et al. Effects of an exopolysaccharide (kefiran) on lipids, blood pressure, blood glucose, and constipation[J]. Biofactors, 2004, 22(1/2/3/4): 197-200.DOI:10.1002/biof.5520220141.

[8] ZHANG J, ZHAO X, JIANG Y, et al. Antioxidant status and gut microbiota change in an aging mouse model as influenced by exopolysaccharide produced by Lactobacillus plantarum YW11 isolated from Tibetan kefir[J]. Journal of Dairy Science, 2017.DOI:10.3168/jds.2016-12480.

[9] DEEPAK V, RAMACHANDRAN S, BALAHMAR R M, et al.In vitro evaluation of anticancer properties of exopolysaccharides from Lactobacillus acidophilus in colon cancer cell lines[J]. In Vitro Cellular & Developmental Biology-Animal, 2016, 52(2): 163-173.DOI:10.1007/s11626-015-9970-3.

[10] LAWS A, GU Y, MARSHALL V. Biosynthesis, characterisation, and design of bacterial exopolysaccharides from lactic acid bacteria[J].Biotechnology Advances, 2001, 19(8): 597-625. DOI:10.1016/S0734-9750(01)00084-2.

[11] 郑华军. 保加利亚乳酸杆菌工业生产菌株2038的基因组学分析[D].上海: 复旦大学, 2010.

[12] MADHURI K V, PRABHAKAR K V. Microbial exopolysaccharides:biosynthesis and potential applications[J]. Oriental Journal of Chemistry, 2014, 30(3): 1401-1410. DOI:10.13005/ojc/300362.

[13] 丹彤, 王俊国, 张和平. 乳酸菌胞外多糖的结构、生物合成及其应用[J].食品科学, 2013, 34(7): 335-339. DOI:10.3969/j.issn.1001-2230.2013.03.007.

[14] 田辉. 嗜热链球菌高密度培养与直投式发酵剂开发[D]. 哈尔滨:东北农业大学, 2012.

[15] LUO R, LIU B, XIE Y, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler[J]. Gigascience, 2012,1(1): 18. DOI:10.1186/2047-217X-1-18.

[16] WU S, ZHU Z, FU L, et al. WebMGA: a customizable web server for fast metagenomic sequence analysis[J]. BMC Genomics, 2011, 12(1):444. DOI:10.1186/1471-2164-12-444.

[17] GRANT J R, STOTHARD P. The CGView Server: a comparative genomics tool for circular genomes[J]. Nucleic Acids Research, 2008,36(Suppl 2): W181-W184. DOI:10.1093/nar/gkn179.

[18] KANEHISA M, GOTO S. KEGG: kyoto encyclopedia of genes and genomes[J]. Nucleic Acids Research, 2000, 28(1): 27-30.DOI:10.1093/nar/28.1.27.

[19] SIGUIER P, PÉROCHON J, LESTRADE L, et al. ISfinder: the reference centre for bacterial insertion sequences[J]. Nucleic Acids Research, 2006, 34(Suppl 1): D32-D36. DOI:10.1093/nar/gkj014.

[20] CARVER T J, RUTHERFORD K M, BERRIMAN M, et al. ACT: the Artemis comparison tool[J]. Bioinformatics, 2005, 21(16): 3422-3423.DOI:10.1093/bioinformatics/bti553.

[21] 邵丽. 产胞外多糖乳杆菌的筛选及其多糖的分离, 结构和生物活性研究[D]. 无锡: 江南大学, 2015.

[22] REN W, XIA Y, WANG G, et al. Bioactive exopolysaccharides from a S. thermophilus strain: screening, purification and characterization[J].International Journal of Biological Macromolecules, 2016, 86: 402-407. DOI:10.1016/j.ijbiomac.2016.01.085.

[23] WU Q L, TUN H M, LEUNG F C, et al. Genomic insights into high exopolysaccharide-producing dairy starter bacterium Streptococcus thermophilus ASCC 1275[J]. Scientific Reports, 2014, 4: 4974.DOI:10.1038/srep04974.

[24] BROADBENT J R, MCMAHON D J, WELKER D L, et al.Biochemistry, genetics, and applications of exopolysaccharide production in Streptococcus thermophilus: a review[J]. Journal of Dairy Science, 2003, 86(2): 407-423. DOI:10.3168/jds.S0022-0302(03)73619-4.

[25] TYVAERT G, MOREL C, JOLY J, et al. The constant gene orf14.9,which belongs to the variable EPS (exopolysaccharide) cluster, is involved in the cell growth of Streptococcus thermophilus[J]. Canadian Journal of Microbiology, 2006, 52(9): 908-912. DOI:10.1139/w06-047.

[26] IYER R, TOMAR S K, UMA MAHESWARI T, et al. Streptococcus thermophilus strains: multifunctional lactic acid bacteria[J].International Dairy Journal, 2010, 20(3): 133-141. DOI:10.1016/j.idairyj.2009.10.005.

[27] GALIA W, PERRIN C, GENAY M, et al. Variability and molecular typing of Streptococcus thermophilus strains displaying different proteolytic and acidifying properties[J]. International Dairy Journal,2009, 19(2): 89-95. DOI:10.1016/j.idairyj.2008.08.004.

[28] DE VUYST L, WECKX S, RAVYTS F, et al. New insights into the exopolysaccharide production of Streptococcus thermophilus[J].International Dairy Journal, 2011, 21(9): 586-591. DOI:10.1016/j.idairyj.2011.03.016.

[29] BOLOTIN A, QUINQUIS B, RENAULT P, et al. Complete sequence and comparative genome analysis of the dairy bacterium Streptococcus thermophilus[J]. Nature Biotechnology, 2004, 22(12): 1554-1558.DOI:10.1038/nbt1034.

[30] SUN Z H, CHEN X, WANG J, et al. Complete genome sequence of Streptococcus thermophilus strain ND03[J]. Journal of Bacteriology,2011, 193(3): 793-794. DOI:10.1128/JB.01374-10.

[31] BAI Y, SUN E, SHI Y, et al. Complete genome sequence of Streptococcus thermophilus MN-BM-A01, a strain with high exopolysaccharides production[J]. Journal of Biotechnology, 2016,224: 45-46. DOI:10.1016/j.jbiotec.2016.03.003.

[32] VENDRAMIN V, TREU L, CAMPANARO S, et al. Genome comparison and physiological characterization of eight Streptococcus thermophilus strains isolated from Italian dairy products[J]. Food Microbiology, 2017, 63: 47-57. DOI:10.1016/j.fm.2016.11.002.

[33] CUI Y, XU T, QU X, et al. New insights into various production characteristics of Streptococcus thermophilus strains[J]. International Journal of Molecular Sciences, 2016, 17(10): 1701. DOI:10.3390/ijms17101701.

[34] ERKUS O, OKUKLU B, YENIDUNYA A F, et al. High genetic and phenotypic variability of Streptococcus thermophilus strains isolated from artisanal Yuruk yoghurts[J]. LWT-Food Science and Technology,2014, 58(2): 348-354. DOI:10.1016/j.lwt.2013.03.007.

[35] MORANDI S, BRASCA M. Safety aspects, genetic diversity and technological characterisation of wild-type Streptococcus thermophilus strains isolated from north Italian traditional cheeses[J]. Food Control,2012, 23(1): 203-209. DOI:10.1016/j.foodcont.2011.07.011.

[36] NATHANI N M, PATEL A K, SENAN S, et al. Genomic analysis of dairy starter culture Streptococcus thermophilus MTCC 5461[J].Journal of Microbiology and Biotechnology, 2013, 23(4): 459-466.DOI:10.3389/fmicb.2012.00435.

[37] PACHEKREPAPOL U, LUCEY J A, GONG Y, et al.Characterization of the chemical structures and physical properties of exopolysaccharides produced by various Streptococcus thermophilus strains[J]. Journal of Dairy Science, 2017, 100(5): 3424-3435.DOI:10.3168/jds.2016-12125.

[38] CUI Y, JIANG X, HAO M, et al. New advances in exopolysaccharides production of Streptococcus thermophilus[J]. Archives of Microbiology,2017, 199(6): 799-809. DOI:10.1007/s00203-017-1366-1.

Genomic and Phenotypic Characterization of the Biosynthesis Pathway of Exopolysaccharides in Streptococcus thermophilus KLDS3.1012

LI Bailiang1, ZHAO Li1, WANG Chengfeng1, JIN Da1, DING Xiuyun1,2, LIU Fei1,*, HUO Guicheng1
(1. Key Laboratory of Dairy Science, Ministry of Education, Northeast Agricultural University, Harbin 150030, China;2. Guangzhou Genedenovo Biotechnology Co. Ltd., Guangzhou 510000, China)

Abstract: The purpose of this study was to explore the biosynthesis pathway of exopolysaccharides (EPS) in Streptococcus thermophilus KLDS3.1012 at the molecular level and to verify it by the phenotypic traits. Firstly, genomic sequencing was performed using Illumina HiSeq and MiSeq sequencing platforms and a circular genomic map was skeletonized; subsequently, glycometabolism, sugar nucleotides synthesis and EPS gene clusters were analyzed by bioinformatics; fi nally, carbohydrate utilization was measured by API 50CH and the monosaccharide composition of the EPS produced by S. thermophilus KLDS3.1012 was analyzed using high performance anion exchange chromatography (HPAEC).In silico bioinformatic analysis unveiled that the genome of S. thermophilus KLDS3.1012 was related to sucrose, fructose,mannose, glucose, galactose and lactose uptake system, and biosynthesis of four sugar nucleotides (UDP-gluctose, dTDP-rhamnose, UDP-galactose and UDP-N-acetylglucosamine) as well as an EPS gene cluster. The phenotypic traits revealed that S. thermophilus KLDS3.1012 could utilize these six sugars and the EPS produced it was composed of glucose, rhamnose and galactose. These results provide a theoretical basis for better understanding the relationship between the genetic elements and the structure of EPS in S. thermophilus KLDS3.1012 and for facilitating its potential applications as starter culture in the dairy industry.

Keywords: Streptococcus thermophilus; genomic sequencing; bioinformatic analysis; glycometabolism; exopolysaccharides DOI:10.7506/spkx1002-6630-20180122-290

收稿日期:2018-01-22

基金项目:“十三五”国家重点研发计划重点专项(2017YFD0400303);国家自然科学基金青年科学基金项目(31401512)

第一作者简介:李柏良(1989—)(ORCID: 0000-0001-8282-6401),男,博士研究生,研究方向为食品科学。E-mail: 15846092362@163.com

*通信作者简介:刘飞(1980—)(ORCID: 0000-0001-8827-4409),男,副研究员,博士,研究方向为乳酸菌资源开发。E-mail: david.as@163.com

中图分类号:R151.1

文献标志码:A

文章编号:1002-6630(2019)06-0136-07

引文格式:李柏良, 赵莉, 王成凤, 等. 嗜热链球菌KLDS3.1012胞外多糖合成途径的基因组学及表型特征分析[J]. 食品科学, 2019,40(6): 136-142. DOI:10.7506/spkx1002-6630-20180122-290. http://www.spkx.net.cn

LI Bailiang, ZHAO Li, WANG Chengfeng, et al. Genomic and phenotypic characterization of the biosynthesis pathway of exopolysaccharides in Streptococcus thermophilus KLDS3.1012[J]. Food Science, 2019, 40(6): 136-142. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-20180122-290. http://www.spkx.net.cn