本发明属于生物工程领域，具体涉及一株新的可广谱裂解蜡样芽胞杆菌的噬菌体及其在防治病原蜡样芽胞杆菌中的应用。
背景技术：
：：蜡样芽胞杆菌是一种广泛存在于环境中的食源性条件致病微生物（drobniewskifa:bacillus-cereusandrelatedspecies.clinmicrobiolrev1993,6(4):324-338）。据研究报道，蜡样芽胞杆菌由于能够产生肠毒素和呕吐毒素而导致人体食物中毒，并引起腹泻和呕吐等症状（bottoneej:bacilluscereus,avolatilehumanpathogen.clinmicrobiolrev2010,23(2):382-398）。此外，蜡样芽胞杆菌还可在其生长期产生诸如磷脂酶、蛋白酶和溶血素等其他毒素而导致多种非肠道感染的人类疾病（drobniewskifa:bacillucereusandrelatedspecies.clinmicrobiolrev1993,6(4):324-338）。由于抗生素滥用问题的存在及抗生素抗性基因的广泛传播，越来越多的耐药蜡样芽胞杆菌菌株被发现，这使得寻求抗生素之外的新的蜡样芽胞杆菌防控手段变得极为迫切（simmr,vorosa,ekmanjvetal:bc4707isamajorfacilitatorsuperfamilymultidrugresistancetransportproteinfrombacilluscereusimplicatedinfluoroquinolonetolerance.plosone2012,7(5)）。由于噬菌体可特异性裂解病原微生物的性质，被认为是一种新型的抗菌物质（coopercj,mirzaeimk,nilssonas:adaptingdrugapprovalpathwaysforbacteriophage-basedtherapeutics.frontmicrobiol2016,7）。相较于传统抗生素而言，噬菌体是地球上丰度最高、最多样的生命体，它可实现与宿主共进化以克服细菌的噬菌体抗性产生的问题。然而，细菌仍有10-7的概率通过形成突变株而产生对噬菌体的抗性。通过发掘更多的噬菌体，以构建“噬菌体鸡尾酒”可更有效的降低细菌产生噬菌体抗性的概率。截止到目前为止，尽管许多可裂解蜡样芽胞杆菌的噬菌体已被发现并研究（gillisa,mahillonj:phagespreyingonbacillusanthracis,bacilluscereus,andbacillusthuringiensis:past,presentandfuture.viruses-basel2014,6(7):2623-2672;kongm,ryus:bacteriophagepbc1anditsendolysinasanantimicrobialagentagainstbacilluscereus.applenvironmicrob2015,81(7):2274-2283），但新的蜡样芽胞杆菌噬菌体发现的速度正在放慢，这是由于大部分新分离的噬菌体表现出与前期已分离噬菌体的高度基因组相似性(breitbartm,rohwerf:hereavirus,thereavirus,everywherethesamevirus.trendsinmicrobiology2005,13(6):278-284;yuany,gaom,pengq,wud,liup,wuy:genomicanalysisofaphageandprophagefromabacillusthuringiensisstrain.jgenvirol2014,95(pt3):751-761.)。在遗传信息上相似性较低的噬菌体可能具有更多样的感染机制；同时，利用不同遗传背景的噬菌体构建“噬菌体鸡尾酒”可更有效减少噬菌体抗性的出现，并且可拓宽噬菌体制剂的裂解谱(baij,kimyt,ryus,leejh:biocontrolandrapiddetectionoffood-bornepathogensusingbacteriophagesandendolysins.frontmicrobiol2016,7.)。技术实现要素：本发明着眼于上述问题，从医院污水中分离到一株新的芽胞杆菌噬菌体vb_bcem-hse3并研究了其生物学特性和遗传学特性。该噬菌体不仅可裂解人类食源性病原菌蜡样芽胞杆菌菌株6113，还可裂解其他蜡样芽胞杆菌菌群细菌，表现出了极其广泛的裂解活性、较高的ph及温度耐受性。基因组学研究发现该噬菌体与其他已测序噬菌体未显示出任何的相似性，因此在遗传学上该株噬菌体属于一株全新的噬菌体。本发明中噬菌体vb_bcem-hse3的全基因组序列已由发明人提交到genbank数据库，genbank序列号为mf418016。为了实现上述的目的，本发明的大概方案如下：利用临床分离的人类病原菌蜡样芽胞杆菌菌株6113（cctccm2017799，genbank序列号njno00000000）作为指示菌，利用双层琼脂平板法从污水环境中分离纯化得到一株噬菌体vb_bcem-hse3(方法参考yuanyh,pengq,wudd,kouz,wuy,liupm,gaomy:effectsofactin-likeproteinsencodedbytwobacilluspumilusphagesonunstablelysogeny,revealedbygenomicanalysis.applenvironmicrob2015,81(1):339-350)。本发明的噬菌体vb_bcem-hse3已经在中国典型培养物保藏中心保藏，地址：中国湖北武汉，武汉大学，保藏日期：2017年12月15日，保藏编号：cctccm2017798，分类命名为bacilluscereusphagevb_bcem-hse3。其特征在于该噬菌体为肌尾噬菌体，电镜观察结果显示该噬菌体由一个直径118nm的正二十面体头部和一个长度171nm的尾部构成，噬菌体尾管的宽度6.5nm和尾鞘的宽度为27nm（图1）。对该噬菌体的生物学特性进行了进一步的研究，采用双层琼脂平板法测定了噬菌体vb_bcem-hse3的裂解谱，结果显示该噬菌体可侵染包括蜡样芽胞杆菌、炭疽芽孢杆菌、苏云金芽胞杆菌等几个菌种的细菌菌株。该噬菌体会在蜡样芽孢杆菌6113的琼脂培养基表面形成透明的0.5mm大小的噬菌斑（图2）。同时探究了该噬菌体侵染其宿主菌蜡样芽胞杆菌6113的一步生长曲线（图3），在噬菌体感染宿主菌的初始期60min为潜伏期，而到180min后进入裂解稳定期，噬菌体的终浓度可达到1012pfu/ml，该噬菌体侵染细菌6113的释放量约为95pfu/细胞。一步生长曲线测定结果显示该噬菌体侵染宿主细菌后，可显著抑制宿主菌的生长，并且导致宿主菌浓度在150min左右发生显著下降。该噬菌体也显示出了较高的ph耐受性和温度耐受性。在低于37°c以下的不同温度处理30min后，vb_bcem-hse3极为稳定；在50°c的条件下处理30min，该噬菌体只能保留其初始裂解活性的1.2%（图4）。当在ph范围3.0到11.0的不同缓冲液中处理该噬菌体后，该噬菌体均具有较高的裂解活性，其中在ph5.0到ph7.0之间该噬菌体的裂解活性不会发生显著的改变（图5）。为了了解噬菌体vb_bcem-hse3的遗传背景，采用illuminahiseq2500测序仪对该噬菌体进行了基因组测序，测序结果显示该噬菌体基因组大小为124.002kb，gc含量为36.21%，包含144个编码序列(cdss)，编码序列的平均长度为791.01bp，两个编码序列平均间隔距离70.6bp。在该噬菌体编码的144个cds中，共有92个cds在ncbi的nr数据库中有同源序列(e-value£10-3)。其中，48个cds只显示出与蜡样芽胞杆菌噬菌体phagesp-15的序列相似性，其他44个cds显示出与其他蜡样芽胞杆菌噬菌体的相似性。除了这92个cds外，其余52个cds为噬菌体vb_bcem-hse3的特异基因，这些特异基因可能赋予了噬菌体vb_bcem-hse3某些独特的性状。在该噬菌体编码的144个基因中，共有55个cds有预测的功能，其中18个cds编码噬菌体结构蛋白，29个cds编码噬菌体复制和核酸代谢相关蛋白，2个cds编码噬菌体宿主裂解系统蛋白，其他cds编码基因表达调控相关蛋白，这些cds呈模块化存在于基因组中（图6）。将预测的基因与virulencefactor(vfdb)数据库进行比对，结果发现该噬菌体基因组上不存在毒力因子相关基因，这意味着在该噬菌体实际的应用中有着较高的安全性。通过将噬菌体的基因组与ncbinr数据库进行比对，结果在数据库中未发现任何与噬菌体vb_bcem-hse3具有较高相似性的基因组序列，仅芽胞杆菌噬菌体sp-15与vb_bcem-hse3的部分区域显示出较高的相似性。因此在分类学上，该株噬菌体vb_bcem-hse3为一株全新的噬菌体。附图说明图1为噬菌体vb_bcem-hse3的透射电镜显微形态观察。图2为噬菌体vb_bcem-hse3裂解蜡样芽孢杆菌6113的噬菌斑形态。图3为噬菌体vb_bcem-hse3的一步生长曲线。图4为噬菌体vb_bcem-hse3的温度稳定性研究。图5为噬菌体vb_bcem-hse3的ph稳定性研究。图6为噬菌体vb_bcem-hse3的基因组组成。具体实施方式为了便于对本发明的进一步理解，下面提供的实施例对其做了更详细的说明。但是这些实施例仅供更好的理解发明而并非用来限定本发明的范围或实施原则，本发明的实施方式不限于以下内容。实施例1：噬菌体分离和制备蜡样芽孢杆菌6113分离自湖北武汉梨园医院病人呕吐物，利用16srdna对该菌进行生物学鉴定。以该菌作为噬菌体分离指示菌，取梨园医院污水样品于12,000×g离心10min，弃固体沉淀。将上清加入luria-bertani(lb)肉汤培养基于30°c震摇过夜以富集样品中的噬菌体。随后，混悬液于12,000×g离心10min收集上清，使用0.22-mm微孔滤膜(millipore,massachusetts)将上清过滤除菌。按5%的比例将滤液加入指数生长期的蜡样芽孢杆菌6113液体培养基中培养8h。接着，将上述溶液离心收集上清，再以0.22-mm微孔滤膜过滤，噬菌体存在于滤液中，收集滤液于4°c保存备用。噬菌体的分离、纯化、繁殖、裂解量测定均是通过双层琼脂覆盖法来完成的，参照yuanyh等方法（yuanyh,pengq,wudd,kouz,wuy,liupm,gaomy:effectsofactin-likeproteinsencodedbytwobacilluspumilusphagesonunstablelysogeny,revealedbygenomicanalysis.applenvironmicrob2015,81(1):339-350），所纯化的噬菌体命名为bacilluscereusphagevb_bcem-hse3。实施例2：噬菌体生物学特性研究1)噬菌体形态观察首先将纯化的噬菌体vb_bcem-hse3悬液滴加于镀有碳膜的铜网上，接着以2%的磷钨酸钾溶液(ph7.2)进行染色，自然干燥后利用透射电镜(h-7000fa;hitachi,tokyo,japan)在加速电压100kv的条件下观察（方法参考yuanyh,pengq,wudd,kouz,wuy,liupm,gaomy:effectsofactin-likeproteinsencodedbytwobacilluspumilusphagesonunstablelysogeny,revealedbygenomicanalysis.applenvironmicrob2015,81(1):339-350）。噬菌体裂解谱测定用于宿主谱测定所用的菌株见表1，测试菌株包括蜡样芽胞杆菌、炭疽芽孢杆菌、苏云金芽胞杆菌、枯草芽孢杆菌、短小芽孢杆菌、大肠杆菌、铜绿假单胞菌、金黄色葡萄球菌、假结核耶尔森菌、鲍曼不动杆菌，其中芽孢杆菌培养条件为条件下于lb肉汤培养基培养。宿主谱的测定方法采用改良的双层琼脂覆盖法（方法参考kongm,ryus:bacteriophagepbc1anditsendolysinasanantimicrobialagentagainstbacilluscereus.applenvironmicrob2015,81(7):2274-2283）。方法如下：指数生长期的测试菌株加入45°c的lb半固体培养基混匀，迅速倾倒平板，稍干后，将108pfu/ml浓度的纯化噬菌体vb_bcem-hse3点接于平板上，在30°c培养过夜，观察其噬菌斑的形成情况。噬菌体生长曲线测定噬菌体vb_bcem-hse3以1.0的感染复数(moi)加入指数生长期的蜡样芽孢杆菌6113(od600=0.8)中吸附5min。吸附完毕后于10,000×g离心1min除去未吸附的噬菌体。所得沉淀重悬于50mllb肉汤培养基并继续培养，每隔30min取样检测噬菌体滴度。记录单独细菌培养物和加入噬菌体后的细菌培养物的生长曲线。检测噬菌体的释放量，方法参照catalaomj等描述的方法（catalaomj,gilf,moniz-pereiraj,pimentelm:themycobacteriophagems6encodesachaperone-likeproteininvolvedintheendolysindeliverytothepeptidoglycan.molmicrobiol2010,77(3):672-686），所有实验均重复三次。噬菌体稳定性测定噬菌体稳定性测定实验包括温度稳定性和ph稳定性。噬菌体悬液分别在4°c、28°c、37°c、50°c、60°c和70°c不同温度下处理30min，随后缓慢冷却，检测每份样品噬菌体的滴度，计算不同温度处理下其存活率。将噬菌体悬液分别加入不同ph条件的缓冲液于4°c处理30min，包括ph1缓冲液（50mmna2hpo4-cl,300mmnacl,10%甘油）、ph3缓冲液（50mmna2hpo4-citricacid,300mmnacl,10%甘油）、ph5缓冲液（50mmna2hpo4-citricacid,300mmnacl,10%甘油）、ph7缓冲液（50mmtris-cl,300mmnacl,10%甘油）、ph9缓冲液（50mmtris-cl,300mmnacl,10%甘油）、ph11缓冲液（50mmna2co3-nahco3,300mmnacl,10%甘油），处理完毕后，立即检测每份样品噬菌体的滴度，计算不同酸碱条件处理下其存活率。噬菌体存活率（viability）的计算方法为处理后存活的噬菌体数与原始噬菌体计数的比值，每组实验重复三次。实施例3：噬菌体基因组的抽提、纯化及测序提取上述制备噬菌体vb_bcem-hse3的dna，方法参照santos等的提取方法并加以改进（yuanyh,pengq,wudd,kouz,wuy,liupm,gaomy:effectsofactin-likeproteinsencodedbytwobacilluspumilusphagesonunstablelysogeny,revealedbygenomicanalysis.applenvironmicrob2015,81(1):339-350）。方法如下：每毫升噬菌体悬液分别加入dnase（终浓度6units/ml）和rnase（终浓度20µg/ml）于37°c处理30min，接着向处理后的噬菌体悬液加入20µlzncl2溶液（终浓度2m）于37°c孵育5min。随后将上述处理液于10,000×g离心1min，除去上清，将沉淀重悬于500µltes缓冲液(0.1mtris/hcl,ph8.0;0.1medta;0.3%sds)中并于60°c孵育15min。然后，向上述溶液中加入20µl蛋白酶k(20mg/ml)并于37°c处理90min，而后加入60µl醋酸钾溶液（3m，ph5.2），混匀后冰浴15min。向上述混悬液加入等体积的苯酚/氯仿/异戊醇（25:24:1，v/v/v）抽提两次，再以氯仿/异戊醇（24:1，v/v）抽提1次，吸取水相并加入等体积的异戊醇混匀，离心得沉淀，将沉淀用70%乙醇洗涤两次后溶解于10µl的无菌水中，完成噬菌体vb_bcem-hse3基因组dna样品的抽提，抽提所得dna经采用qubit2.0和0.8%琼脂糖凝胶电泳进行质检。质检合格的dna可进行后续测序实验。首先，对基因组dna进行片段化，末端修复、3’末端加a、连接接头、富集等步骤，完成测序样本文库构建。所建文库使用qubit®2.0fluorometer检测浓度，agilent2100检测文库的大小。按照cbotuserguide所示相应流程，在illuminahiseq测序仪配套的cbot上完成cluster生成和第一向测序引物杂交。按照illuminauserguide准备测序试剂，将携有cluster的flowcell上机。选用paired-end程序，进行双端测序。测序过程由illumina提供的datacollectionsoftware进行控制，并进行实时数据分析。噬菌体基因组分析所得测序数据除去总体质量较低、含有接头序列、末端质量偏低等不合格的reads，利用veletv1.2.07拼接软件对预处理后的illumina数据拼接，并使用primewalking设计引物补齐序列缺口得到完整序列。运用genemarks软件对完整序列进行基因预测（besemerj,lomsadzea,borodovskym:genemarks:aself-trainingmethodforpredictionofgenestartsinmicrobialgenomes.implicationsforfindingsequencemotifsinregulatoryregions.nucleicacidsres2001,29(12):2607-2618），将预测的基因的蛋白序列分别与non-redundant(nr)和cdd数据库进行blastp比对(blast2.2.28+)(marchler-bauera,panchenkoar,shoemakerba,thiessenpa,geerly,bryantsh:cdd:adatabaseofconserveddomainalignmentswithlinkstodomainthree-dimensionalstructure.nucleicacidsres2002,30(1):281-283.)，使用hhpred(sodingj:proteinhomologydetectionbyhmm-hmmcomparison.bioinformatics2005,21(7):951-960)和pfam(puntam,coggillpc,eberhardtry,mistryj,tatej,boursnellc,pangn,forslundk,cericg,clementsjetal:thepfamproteinfamiliesdatabase.nucleicacidsres2012,40(d1):d290-d301)软件进行蛋白质的功能域和结构组成的分析。采用cgview软件将噬菌体vb_bcem-hse3的基因注释结果进行可视化作图分析。通过gepard软件将该噬菌体与其相似性最高的噬菌体分别进行基因组和蛋白组点杂交比对研究。表1噬菌体vb_bcem-hse3的裂解谱aatcc,americantypeculturecollection;bgsc,bacillusgeneticstorecenter。当前第1页12当前第1页12