耐久肠球菌C11菌株的环境胁迫耐受性及其低温适应相关基因的基因组学鉴定

姚文婷，陈兰明*

1(上海海洋大学 食品学院，上海，201306) 2(农业部水产品贮藏保鲜质量安全风险评估实验室(上海)，上海，201306)

国内外大量研究表明，乳酸菌(lactic acid bacteria, LAB)具有调节宿主肠道微生态平衡、降低胆固醇、调节免疫功能及延缓衰老等生理功能[1]。耐久肠球菌(Enterococcusdurans)隶属于硬杆菌门(Firmicutes)，杆菌纲(Bacilli)，乳酸杆菌目(Lactobacillales)，肠球菌科(Enterococcaceae)，肠球菌属(Enterococcus)。该菌是动物肠道微生物的重要组分，且常见于发酵食品，如酸奶和奶酪中[2]。不同来源的耐久肠球菌的生物学特性和生理功能各不相同[3-4]。例如，黄坚等[5]从牦牛发酵酸奶中筛选到耐久肠球菌SWUN5857菌株，发现该菌株能够适应小鼠胃肠环境，有效提高小鼠体质量，促进小鼠胸腺和脾脏发育，增强小鼠免疫功能。LIU等[6]从内蒙古天然发酵奶油中分离到耐久肠球菌KLDS6.0930菌株，该菌则具有耐胆盐，吸收胆固醇等作用。

耐低温LAB不仅能够有效抑制腐败菌在低温下的生长和繁殖，而且，在食品发酵过程中因其低温适应性，有助于降低生产能耗，节约企业成本。因此，耐低温LAB的筛选对于食品加工与贮藏具有重要意义。然而，大多数LAB属于中温型微生物，一般最适生长温度为30～40 ℃；在低于10 ℃的温度下，LAB生长微弱，甚至不生长[7]。迄今为止，国内外有关嗜冷LAB的研究已有报道。例如，LI等[8]报道了3株在15 ℃生长的植物乳杆菌(Lactobacillusplantarum) LZ95、CY2和CY3菌株，其中，LZ95菌株在de Man-Rogosa-Sharpe(MRS)培养基中培养48 h后，OD620nm值增加了10倍。LI等[9]从藏牦牛瘤胃中分离到具有纤维素分解能力的粪肠球菌(Enterococcusfaecalis) JF85和Y83菌株，发现它们能够在15 ～55 ℃的温度下生长。

本研究基于前期研究基础，聚焦于1株分离自中国传统发酵食品且具有良好的体外抗氧化性的耐久肠球菌C11菌株(E.duransC11)[10]，分析该菌株在不同环境胁迫条件下的耐受性，并运用比较基因组学方法探讨E.durans的低温适应机制。

1 材料与方法

1.1 材料

1.1.1 菌株

耐久肠球菌C11菌株分离自我国传统发酵食品——乳饼[10]，由本实验室筛选、鉴定并保藏。

1.1.2 培养基、主要试剂及仪器

MRS培养基，北京陆桥技术有限公司；磷酸盐缓冲液(phosphate-buffered saline，PBS) (pH 7.4)、λ/HindIII DNA Marker (MD202)、6×Loading buffer，天根生化科技有限公司；Takara MiniBEST细菌基因组DNA提取试剂盒，上海皓嘉科技发展有限公司；胆盐和琼脂，国药集团化学试剂有限公司；胃蛋白酶(活性1∶3 000)、胰酶(活性1∶250)、微孔过滤膜(0.22 μm)、NaCl，上海生工生物工程有限公司；多功能酶标仪，BioTek公司；雷磁pHS-3C pH仪，上海雷磁仪器厂；Research®plus pipette移液器、台式高速低温冷冻离心机5417R、高速冷冻离心机5424、PCR仪，Eppendorf公司；凝胶成像系统Molecular Imager®Gel DocTMXR+ System，Bio-Rad公司；全自动生长曲线分析仪，Bioscreen公司。

1.2 实验方法

1.2.1E.duransC11菌株的培养

参考XU等[10]的方法，E.duransC11菌株以体积分数1%接种于MRS液体培养基中，置于37 ℃厌氧培养至对数生长期，连续传代培养2代，得到新鲜培养物。4 000 r/min离心10 min，用无菌PBS洗涤菌体3次(避免带入接种物中的酸性物质从而改变pH)。用适量的PBS悬浮菌体沉淀作为种子液备用。

1.2.2 温度、酸、渗透压的耐受性分析

按照1%体积分数的接种量将新鲜种子液接种到MRS液体培养基中，在不同温度(10、15、20、25、37、45 ℃)培养72 h，采用全自动生长曲线分析仪实时检测培养物在600 nm波长处的吸光度值(OD600 nm)，3次重复。

按照1%体积分数的接种量将新鲜种子液分别接种到不同起始pH值(2.0、3.0、4.0、5.0、6.0、7.0)、不同NaCl质量分数(0、5.0%、8.0%、10.0%、12.0%、15.0%)、以及不同胆盐质量分数(0、0.05%、0.10%、0.20%、0.30%)的MRS液体培养基中，于37 ℃培养24 h，实时检测培养物的OD600 nm值，3次重复。

1.2.3 人工胃液及肠液的耐受性分析

参照ZHANG等[11]的方法制备人工模拟胃液(pH 2.0)和人工模拟肠液(pH 6.8)。按照10%体积分数将新鲜种子液分别接种于经膜过滤除菌的1.0 mL人工胃液和人工肠液中，充分混匀，分别于37 ℃孵育180 min和240 min后取样，10倍梯度稀释，以标准平板法计数，计算细菌的存活率，实验3次重复。

1.2.4 细菌基因组DNA的制备

采用Takara MiniBEST细菌基因组DNA提取试剂盒，按照试剂盒说明书步骤提取E.duransC11菌株的基因组DNA。采用0.7%琼脂糖凝胶电泳检测DNA样品，凝胶成像系统拍照记录实验结果，多功能酶标仪测定DNA样品的浓度和纯度(A260/A280)。

1.2.5 基因组DNA序列的测定、质控和装配

运用Illumina HiSeq Xten测序平台，由北京诺禾致源生物信息科技有限公司测定E.duransC11全基因组序列。采用FASTQC软件[12]分析原始序列reads的质量。采用SOAP denovo软件(version 2.04)(http://soap.genomics.org.cn/)进行序列装配。

1.2.6 基因组注释

采用Glimmer (Version 3.0)[13]、tRNA_scan-SE (Version 1.3.1)[14]、RNAmmer (Version 1.2)软件[15]分别预测编码基因、tRNA、rRNA基因。采用Prophage Finder预测噬菌体序列(http://phast.wishartlab.com/)。采用CRISPRFinder软件[16]预测CRISPR (clustered regularly interspaced short palindromic repeats)位点。采用PGAP[17]、CD-Search[18]、SignalP 4.1 Serve[19]、TMHMM[20]软件分别预测假基因、Pfam域的基因、信号肽、跨膜结构域。通过比对毒力因子数据库VFDB(virulence factors database) (http://www.mgc.ac.cn/VFs/)、耐药基因数据库ARDB (antibiotic resistance genes database) (http://arpcard.mcmaster.ca) 分别预测毒力、耐药基因(Identity>70%, E<1e-10)。

1.2.7 比较基因组分析

截止到2018年12月20日，从美国国家生物技术信息中心(National Center for Biotechnology Information, NCBI)基因组数据库(http://www.ncbi.nlm.nih.gov/genome)中检索到17株E.durans菌株的基因组序列，其中3株菌为全基因组序列，14株菌为基因组草图，及E.duransC11共18株耐久肠球菌用于比较基因组分析(表1)。采用BLAST(basic local alignment search tool)软件比对分析菌株特异基因(阈值e-10)。采用CD-HIT[21]、Muscle软件[22]分析同源基因并对其进行多序列比对。采用PhyML软件[23]以邻接算法构建系统发育树，并且进行1 000次bootstrap校验，使用Evolview软件绘制系统发育树[24]。

2 结果与讨论

2.1 不同温度对E. durans C11菌株生长的影响

本研究测定了E.duransC11菌株在不同温度下(10～45 ℃)的生长曲线，结果如图1所示。在45 ℃高温条件下，E.duransC11在MRS培养基(pH 6.8)中几乎不能生长；与最适生长温度(37 ℃)相比，在25 ℃条件下，E.duransC11仍生长较快，约6 h进入对数期，21 h达到对数生长中期，30 h到达稳定期；在20 ℃和15 ℃条件下，该菌株约6 h进入对数期，27 h到达对数生长中期，42 h到达稳定期；在10 ℃条件下，菌株生长缓慢，约12 h进入对数期，40 h到达对数生长中期，54 h到达稳定期。上述结果表明，E.duransC11菌株在MRS培养基中能够适应10～25 ℃的低温条件。

注：“-”未知。

图1 温度对耐久肠球菌C11菌株生长的影响Fig.1 The effect of temperatures on the growth ofE. durans C11

2.2 不同起始pH值对E. durans C11菌株生长的影响

如图2所示，当MRS培养基的pH为2.0～4.0时，E.duransC11菌株生长缓慢甚至不生长，其OD600nm值无明显改变。当pH为5.0～6.0时，E.duransC11菌株能够生长，但与最适生长pH (6.8～7.0)相比，酸性培养条件(pH 5.0～6.0)显著抑制了E.duransC11的生长。

图2 不同起始pH值对耐久肠球菌C11生长的影响Fig.2 The effect of different initial pH on the growth ofE. durans C11

2.3 E. durans C11菌株对NaCl的耐受性

与没有添加NaCl的MRS培养基相比，添加5.0%的NaCl显著抑制了E.duransC11菌株的生长，约6 h进入对数期，16 h到达对数中期，24 h到达稳定期；当MRS培养基中添加的NaCl质量分数大于8.0%时(8.0%～15.0%)，E.duransC11菌株几乎不能生长(图3)。

图3 耐久肠球菌C11对不同质量分数NaCl的耐受性Fig.3 Tolerance of E. durans C11 to different concentrations of NaCl

2.4 E. durans C11菌株对胆盐的耐受性

如图4所示，当胆盐质量分数大于0.05%时，E.duransC11菌株不能生长，表明该菌株对胆盐(0.05%～0.3%)无耐受性。

图4 耐久肠球菌C11对不同质量分数胆盐的耐受性Fig.4 Tolerance of E. durans C11 to different concentrations of bile salt

2.5 E. durans C11菌株对人工胃肠液的耐受性

E.duransC11菌株对人工胃液(pH 2.0)、人工肠液(pH 6.8)的耐受性如表2所示。在人工胃液中处理180 min后，E.duransC11菌株的存活率为52.5%；而在人工肠液中处理240 min后，E.duransC11的存活率为85.1%，表明E.duransC11菌株对人工肠液的耐受性较强。

表2 耐久肠球菌C11对人工胃液及肠液的耐受性Table 2 Tolerance of E. durans C11 to artificialgastric and intestinal fluids

2.6 E. durans C11基因组草图特征

运用Illumina Hiseq二代测序技术，本研究测定获得了E.duransC11的全基因组草图，共获得3 126 666条reads，序列全长为2 988 204 bp，GC含量为37.68%。序列装配得到111个Scaffolds，测序深度为151×。通过基因组序列分析，预测到2 986 个基因，包括2 715个蛋白编码基因，77个RNA基因，194个假基因。E.duransC11基因组中比对到Pfam域、信号肽、跨膜螺旋结构的基因分别有2 278、137、713个。另外，鉴定到3个CRISPR重复序列，其中包含1个未知的CRISPR序列，位于scaffold 18的16 679～16 757 bp之间，长度为24 bp；2个已知的CRISPR序列，分别位于scaffold 30 的 3 942～4 175 bp和scaffold 53的16 419 ～16 618 bp之间，长度分别为24、36 bp，分别编码假设蛋白和复制蛋白(replication protein, RepA)。还鉴定到1个完整的原噬菌体基因簇，位于scaffold 3的580～16 832 bp之间，序列全长为16.2 kb，共编码25个蛋白，其中20个编码原噬菌体相关蛋白，5个编码假设蛋白(表3)。

表3 E. durans C11基因组测序数据Table 3 E. durans C11 genome statistics

注：“-”:未知

E.duransC11基因组中还鉴定到2个耐药基因(tetm、baca)，分别与四环素(tetracycline)、杆菌肽(bacitracin)抗性相关。此外，还发现了17个编码黏附素、产胞外酶、膜表面蛋白、荚膜多糖等相关基因(efaA、eno、gapA、tuf、bsh、GroEL、bopD、srt2/srt1/bee1/bee3、pilF、pilB、hasC)，可能与耐久肠球菌对宿主细胞的黏附、定殖、侵染等相关[27]。

2.7 比较基因组分析

2.7.1E.duransC11基因组预测蛋白质编码基因的COG功能分类

比较基因组分析揭示了18株E.durans菌株的基因组中48 112个蛋白质编码基因，它们在C11菌株基因组中的COG功能分类如表4所示。

表4 E. durans C11基因组预测蛋白质编码基因的COG功能分类Table 4 The number of genes classified into the 24 COGfunctional categories in E. durans C11 genome

续表4

分类E. durans C11基因个数占比/%功能描述F752.76核苷酸转运与代谢S54920.22未知功能P1224.49无机盐离子转运与代谢U230.85细胞内运输、分泌、小泡运输I481.77脂质转运与代谢R00.00胞外结构O622.28翻译后修饰、蛋白质转换、伴侣蛋白L2489.13复制、重组和修复A10.04RNA加工和修饰B00.00染色质结构和动力学Q210.77次级代谢产物的生物合成、转运、分解代谢T612.25信号转导机制K1846.78转录J1545.67翻译、核糖体结构和生物合成Y00.00核结构-1 21017.61未匹配到数据库

E.duransC11基因组中约2 237个蛋白质编码基因比对到COG数据库，分属24个功能分类(表4)。其中，编码未知功能(S)，复制、重组和修复(L)，碳水化合物转运与代谢(G)，转录(K)，翻译、核糖体结构与生物合成(J)基因的占比位居前5位，分别占所有蛋白编码基因的20.22%、9.13%、7.40%、6.78%、5.67%。

2.7.2E.durans基因组系统发育树

基于18株E.durans基因组的同源基因构建了系统发育树(图5)。由图5可见，18株E.durans的基因组聚为3簇，分别为cluster α、 β、γ。其中，cluster γ又聚为2个亚类，分别为Subcluster I和Subcluster II。在Subcluster II中，E.duransC11与E.duransNCTC8130、E.durans18S菌株具有相近的系统发育关系，其中，E.durans18S分离自奶酪，而E.duransNCTC8130的来源未知，它们的生物学特性均未有报道。

图5 18株耐久肠球菌基因组的系统发育树Fig.5 A phylogenetic tree of the 18 E. durans genomes

2.7.3E.duransC11菌株的特异基因

比较基因组分析揭示了E.duransC11基因组中33个菌株特异基因，其中，26个编码假设蛋白，7个分别编码螺旋转角螺旋域蛋白(helix-turn-helix domain-containing protein) (EIA52_00535)、环内酯自诱导肽(cyclic lactone autoinducer peptide) (EIA52_00555)、I型限制性内切酶R亚基(type I restriction endonuclease subunit R) (EIA52_06120)、DUF3037域结合蛋白(DUF3037 domain-containing protein) (EIA52_12720)、切除酶(excisionase) (EIA52_13810)、CPBP家族膜内金属蛋白酶(CPBP family intramembrane metalloprotease) (EIA52_13980, EIA52_13990)，它们与E.duransC11的低温适应性关系机制有待进一步的研究。

2.8 E. durans C11可能的低温适应机制

2.8.1 调节细胞膜脂组成

当环境温度下降时细胞膜的流动性降低，细胞难以发挥正常生理功能。然而，耐低温菌通过改变膜脂组成，如提高不饱和脂肪酸比例来调节细胞膜的流动性[28]。在大肠杆菌(Escherichiacoli)中，参与不饱和脂肪酸合成的2个关键酶为3-羟基癸脂酰ACP脱水异构酶(3-hydroxydecanoyl-ACP dehydratase/isomerase, FabA)和3-酮脂酰ACP合成酶 I(3-Ketoacyl-ACP synthase I, FabB)[29]。研究发现，革兰氏阳性菌粪肠球菌(Enterococcusfaecalis)基因组中含有fabZ1 (EIA52_09905)和fabF1 (EIA52_10800)基因，分别与大肠杆菌fabA和fabB序列高度同源[30]。在本研究中，从E.duransC11基因组中鉴定到fabZ1和fabF1基因，它们与E.faecalis中的fabZ1和fabF1有较高的氨基酸序列同源性(97%, 91%)，推测E.duransC11可能通过FabZ1-FabF1途径合成不饱和脂肪酸，以增加细胞膜的流动性，减少低温对细胞的损伤。

2.8.2 增加细胞内相容性溶质(compatible solutes)

据文献报道，相容性溶质(如甘氨酸、甜菜碱、甘油、海藻糖、甘露醇和山梨醇等)能够在细胞质内积累到很高浓度，并保护细胞免受低温、高温、干燥等应激损害[31]。在本研究中，通过比较基因组分析，在E.duransC11基因组中鉴定到编码甘氨酸/甜菜碱/胆碱转运(glycine betaine/carnitine transport)的基因簇(OpuABCD) (EIA52_01405, EIA52_12825, EIA52_12830, EIA52_12835, EIA52_12840)。研究发现，在15 ℃培养条件下，枯草芽孢杆菌(Bacillussubtilis)利用甜菜碱转运蛋白OpuA、OpuC和OpuD将甜菜碱摄入细胞，或利用前体胆碱合成甜菜碱，从而保护细胞免受低温损伤[32]。

2.8.3 改变蛋白质的氨基酸组成

KREIL等[33]研究发现，当温度降低时，蛋白质中谷氨酰胺(Gln)、丝氨酸(Ser)、苏氨酸(Thr)、组氨酸(His)的含量增加。另外，嗜冷酶通过疏水性脯氨酸(Pro)含量的减少来提高酶的柔韧性[34]，使酶的活性中心更容易与底物接近，使之在低耗能时发挥作用。本研究中E.duransC11在低温条件下表达蛋白质中Gln、Ser、Thr、His和Pro的含量变化情况有待蛋白质组学的进一步研究。

2.8.4 冷激蛋白和热激蛋白

冷激蛋白(cold shock proteins, Csps)参与细胞内多种代谢途径，如转录、翻译、蛋白质折叠，以及细胞膜流动性的调控等[35]。热激蛋白(heat shock proteins, Hsps)不仅在细菌的热应激,而且在冷应激中也发挥作用,原因在于许多Hsps是分子伴侣，在低温及其他应激条件下同样被诱导表达[36]。研究发现，当温度发生变化时，分子伴侣可以在原位调控蛋白质暴露的活性区域之间的相互作用，使蛋白折叠成为正确的构象[37]。在本研究中，比较基因组分析揭示了E.duransC11基因组4个分子伴侣：GrpE (EIA52_03480)、DnaJ (EIA52_03490)、DnaK (EIA52_03485)、hslO (EIA52_12115)；2个Hsps (EIA52_03735, EIA52_12085)；1个Hsps转录调节因子hrcA (EIA52_03475)；6个Csps (EIA52_03155, EIA52_04230, EIA52_04370, EIA52_09445, EIA52_13120, EIA52_11270)。此外，在E.duransC11基因组中还鉴定到1个DEAD-box ATP-dependent RNA解旋酶编码基因(CshB，EIA52_09480)，与DEAD-box RNA解旋酶CsdA(Psyc_1082)氨基酸序列同源性为32%。当温度下降时，CsdA促进mRNA的解旋，调控微生物基因表达[38]。

2.8.5 应激相关调节子

转录组学研究发现，在低温条件下，副溶血性弧菌(Vibrioparahaemolyticus) CHN25的应激相关调节子编码基因表达上调，例如LysR、GntR、MerR家族转录调节子、3’-5’-环磷酸腺苷(cAMP)受体蛋白(CRP)。在本研究中，在E.duransC11基因组中鉴定到CRP(EIA52_03380, EIA52_06660, EIA52_08790, EIA52_13485)、LYSR(EIA52_00465, EIA52_09350, EIA52_09710, EIA52_12995)、GntR(EIA52_00775, EIA52_00995, EIA52_03005, EIA52_04360, EIA52_05535, EIA52_11320)、MerR(EIA52_02670, EIA52_03215, EIA52_03870) 家族转录调节子等编码基因[39]。CRP是一种参与糖代谢的调节因子，在大肠杆菌的冷适应中发挥重要作用。据研究报道，细菌对于零下温度的适应涵盖了能量的合成、翻译、运输等方面[40]。例如，低温嗜冷杆菌(Psychrobactercryohalolentis) K5菌株能够在-4 ℃下诱导一些重要调控子的表达，如调节能量需求的F1/F0 ATP合成酶，在低温下排除有害物质的外膜外排系统蛋白(TolC)，加快翻译过程的延伸因子Ts(EF-Ts)等，这些基因在E.duransC11基因组中均被鉴定到(EIA52_00860, EIA52_06915, EIA52_06625)，可能与其低温适应性相关。

3 结论

本研究发现耐久肠球菌C11菌株能够在15～25 ℃温度条件下生长。该菌株的基因组草图已被测定，全长2 988 204 bp，GC含量37.68%，GenBank序列登录号为RQWF00000000。比较基因组学分析揭示了E.duransC11的33个菌株特异性基因，以及大量参与细胞膜脂组成、相容性溶质吸收或合成、环境胁迫应激调控等相关基因，可能与其低温适应相关。转录组学和蛋白质组学的进一步分析将有助于阐释耐久肠球菌低温适应的分子调控网络。