一种快速、高通量检测幽门螺旋杆菌耐药性的方法与流程

：
1.本发明属于微生物技术领域以及医药技术领域，具体涉及一种快速、高通量检测幽门螺旋杆菌耐药性的方法。


背景技术：

2.幽门螺杆菌(helicobacter pylori，hp)是一种革兰氏阴性、微需氧螺旋状细菌，与人类胃肠道疾病有关。hp感染可引起慢性胃炎，并导致多种并发症，如消化不良、消化性溃疡和胃癌。据统计，全球约90％的胃癌发病与hp感染密切相关。因此，hp感染严重威胁人类的健康、造成庞大的经济与社会损失，预防和治疗幽门螺旋杆菌感染具有巨大的应用研究价值。
3.抗生素治疗是目前清除hp的主要方案，但由于胃内存在大量胃酸及消化酶，可灭活绝大部分抗生素，因此现有的多数抗生素不能用于hp的治疗。迄今为止，仅5种抗生素在全球范围内获批用于hp的治疗：阿莫西林、克拉霉素、左氧氟沙星、甲硝唑及四环素。然而，大量研究数据指出，我国流行的hp对上述五种抗生素的耐药率在逐年爬升。至2017年，我国hp流行株对不同抗生素的耐药情况为：克拉霉素28.9％、甲硝唑63.8％、左氧氟沙星28.0％、阿莫西林3.1％、四环素3.9％。
4.hp耐药性检测的金标准为美国临床和实验室标准协会(clinical and laboratory standards institute，clsi)提出的琼脂稀释法。该方法对实验人员专业水平、实验设备、实验耗材均有较高要求，难以在临床应用中开展。因此，寻求简便、准确率高的hp耐药性评判新标准有重要的应用意义。
5.目前认为hp耐受抗生素主要由该菌基因组中特定基因发生随机突变导致。hp基因发生突变后，抗生素难以结合到其作用靶点或药物难以在菌体内由前体药物转化为作用药物，均导致hp对抗生素耐药。因此，研究者通过大数据分析手段统计分析hp基因突变与耐药表型的规律，继而建立了基于hp特定基因突变预测其耐药表型的模型，为临床医生提供了一套简便、快速的hp药敏预测方案。
6.近年来，随着新一代测序技术的发展，微生物学家们逐渐发现hp是一个在种内具有高度变异性的物种，不同大洲人群感染的hp在基因组水平上差异极大。


技术实现要素：

7.本发明的第一个目的是提供一种准确地预测hp中国流行株对抗生素的耐药情况的快速、高通量检测幽门螺旋杆菌耐药性的方法。
8.本发明人团队在前期研究中发现，既往文献报道的其他国家发生的hp耐药基因突变并不具有预测中国流行hp株药敏表型的意义(图2)。
9.因此，本研究结合我国hp流行株的基因组特征，建立起hp中国流行株耐药基因的分子表征库，并设计出hp中国流行株的耐药基因组分子表征关联模型，针对我国流行hp药敏情况，对于三种易发生耐受的抗生素，设计出适用于我国hp分离株耐药性的快速检测方
案。本方案对比现有方法，有利于协助临床医生对我国hp感染的治疗做出更为精准的判断。
10.本发明的快速、高通量检测幽门螺旋杆菌耐药性的方法，包括以下步骤：
11.分离获取幽门螺旋杆菌样品，对幽门螺旋杆菌进行如下检测：
12.检测克拉霉素耐药基因23s rrna的突变位点：109g
→
a、2056a
→
g、2095t
→
c、2834c
→
t和蛋白infb的突变位点794glu
→
asp，如发生任一突变则幽门螺旋杆菌具有克拉霉素耐药性；
13.或检测左氧氟沙星耐药基因的蛋白gyra的突变位点：87asn
→
lys、91asp
→
asn/gly/tyr和蛋白gyrb的突变位点：684glu
→
asp，如发生任一突变则幽门螺旋杆菌具有左氧氟沙星耐药性；
14.或检测甲硝唑耐药基因的蛋白rdxa的突变位点：16arg
→
cys/his和蛋白rpsu的突变位点：56val
→
ile和蛋白sodb的突变位点：54asp
→
ala，如发生任一突变则幽门螺旋杆菌具有甲硝唑耐药性。
15.本发明的第二个目的是提供一种快速、高通量检测幽门螺旋杆菌耐药性的试剂盒，其包括：
16.检测克拉霉素耐药基因23s rrna的突变位点：109g
→
a、2056a
→
g、2095t
→
c、2834c
→
t和蛋白infb的突变位点794glu
→
asp的试剂；
17.或检测左氧氟沙星耐药基因的蛋白gyra的突变位点：87asn
→
lys、91asp
→
asn/gly/tyr和蛋白gyrb的突变位点：684glu
→
asp的试剂；
18.或检测甲硝唑耐药基因的蛋白rdxa的突变位点：16arg
→
cys/his和蛋白rpsu的突变位点：56val
→
ile和蛋白sodb的突变位点：54asp
→
ala的试剂。
19.本发明与现有技术相比，其有益效果为：
20.本发明经研究发现现有的检测方案不能准确地预测hp中国流行株对抗生素的耐药情况，因此发明人建立了首个基于hp中国流行株耐药基因的分子表征库信息的基因
‑
表征关联模型，提出基于该模型的新型快速、高通量hp耐药检测方案。与现有技术相比，本发明能更精准地协助临床医生快速预测我国hp的耐药表型。
附图说明
21.图1是hp中国流行株耐药基因的分子表征库；hp中国流行株耐药基因的分子表征库部分内容：54株中国分离hp株对克拉霉素、左氧氟沙星、甲硝唑的药敏表型及耐药基因分子特征。54株中国分离hp株的进化关系如左侧进化树所示。长形方框代表菌株对相应抗生素的药敏表型，代表敏感、代表耐药；正方形代表菌株在特定基因位点的情况，代表与atcc26695一致，无突变发生，代表存在突变。深灰色背景的突变位点为国外文献报道与耐药基因或蛋白突变位点，白色背景的突变位点为本发明公布的与中国流行株耐药相关的基因或蛋白突变位点。
22.图2是不同耐药基因、蛋白突变对hp中国流行株的三种抗生素mic的影响；图a、b、c分别显示不同耐药基因或蛋白突变对hp中国流行株mic的影响。白色为无突变hp分离株的mic值，黑色为耐药基因突变hp分离株的mic值。可见发生基因突变的hp株mic值高于无突变株。图中*所标注的突变位点为本专利新发现的与hp中国流行株耐药相关的基因或蛋白突
变位点。
具体实施方式
23.以下实施例是对本发明的进一步说明，而不是对本发明的限制。
24.下述实施例中涉及的培养基及试剂如下：
25.培养基添加抗生素混悬液的制备：万古霉素10mg/ml、两性霉素b 5mg/ml、甲氧苄氨嘧啶20mg/ml、头孢磺啶10mg/ml，上述抗生素均采用灭菌水溶解。
26.hp选择性培养平板(g/l)：酵母浸粉3.0g/l、酪胨12.0g/l、动物组织消化物5.0g/l、牛肉浸粉3.0g/l、玉米淀粉1.0g/l、氯化钠5.0g/l、琼脂13.5g/l、无菌绵羊红细胞或胎牛血清70ml/l、抗生素混悬液10ml/l。hp增殖培养平板为不添加抗生素混悬液的hp选择性培养平板。
27.药敏检测用抗生素悬液的制备：(1)克拉霉素：称取100mg克拉霉素，溶解于0.2ml甲醇或冰醋酸中配成贮存液，使用时采用ph 6.5、浓度为0.1mol/l的磷酸盐缓冲液稀释至所需浓度；(2)左氧氟沙星：称取100mg左氧氟沙星，加0.2ml水后缓慢滴加浓度为0.1mol/l氢氧化钠直至溶解配成贮存液，使用时使用无菌水稀释至所需浓度；(3)甲硝唑：称取100mg甲硝唑，溶解于0.2ml二甲基亚砜中配成贮存液，使用时采用无菌水稀释至所需浓度；(4)阿莫西林：称取10mg阿莫西林，溶解于1ml ph 6.0、浓度为0.1mol/l的磷酸盐缓冲液中配成贮存液，使用时采用ph 6.0、浓度为0.1mol/l的磷酸盐缓冲液稀释至所需浓度；(5)四环素：称取抗生素粉末直接加无菌水配成所需浓度。
28.实施例1幽门螺杆菌的分离、培养与鉴定
29.(1)hp的分离及培养
30.①
hp的分离：在无菌环境下，取慢性胃炎患者胃粘膜组织，涂布于hp选择性培养平板上，后放置于37℃微需养环境中(5％o2、10％co2、85％n2)培养5
‑
7天。挑取平板上形态典型的菌落至新的hp选择性培养平板上进行划线纯化，纯化2次后获得纯化的hp分离株。
31.②
hp的培养：从含hp的选择性培养平板上挑取单菌落，接种于hp增殖培养平板上，37℃微需养环境中培养72
‑
96小时，加入生理盐水洗脱获得hp分离株菌液。
32.(2)hp的鉴定
33.采用基于16s rrna基因序列的方法对分离菌株进行分子鉴定：采用细菌dna提取试剂盒(mabio，china)进行提取hp的基因组dna，后采用2
×
pcr mix(dongshengbio，china)进行pcr扩增。pcr扩增引物采用16s rrna基因通用引物，上游引物序列为27f：5
’‑
aga gtt tga tcc tgg ctc ag
‑3’
；下游引物序列为1492r：5
’‑
ctac ggc tac ctt gtt acg a
‑3’
。pcr反应条件为：预变性95℃5min；95℃30s，56℃30s，72℃1min 30s共30个循环，72℃退火延伸10min。对pcr产物进行切胶回收，后进行sanger测序。获得的16s rrna基因序列后比对ncbi数据库(https://blast.ncbi.nlm.nih.gov)，结果提示其与hp同源性最高，且identity和coverage均大于99％，可确定分离菌株为hp。本研究采用的54株中国hp分离株16s rrna基因比对结果显示与hp标准株atcc26695高度同源，identity和coverage均大于99％。
34.实施例2幽门螺旋杆菌的药敏检测
35.采用琼脂稀释法对hp进行抗生素药敏检测，采用hp atcc 43504作为质控菌株。根
据美国临床和实验室标准协会(the clinical&laboratory standards institute，clsi)制定的方法和标准进行测定。
36.琼脂稀释法检测hp对不同抗生素的敏感性：将阿莫西林、克拉霉素、左氧氟沙星、甲硝唑及四环素五种药物的粉末根据clsi推荐的方式溶解于相应溶质中，配置成抗生素贮存液，后根据所需检测的药敏浓度，分别加入至hp选择性培养平板中，制成含相应浓度抗生素的hp药敏检测平板。将hp菌悬液浓度调整为1
×
107cfu/ml，并接种至含相应浓度抗生素的检测平板上，在37℃微需养环境中培养72小时，观察含相应浓度抗生素的检测平板上有无hp菌落形成。以不能形成菌落的最低浓度药物平板中的药物浓度为该种抗生素对所检hp的最小抑菌浓度(the minimum inhibitory concentrations，mic)。每个药物浓度进行平行操作3次，得出最小抑菌浓度的平均值。hp对抗生素耐受(resistant，r)的评判标准为：阿莫西林≥0.125mg/l、克拉霉素≥0.25mg/l、左氧氟沙星≥1mg/l、甲硝唑≥8mg/l、四环素≥1mg/l。
37.实施例3基于全基因组测序的幽门螺旋杆菌中国流行株耐药表征及分子特征库构建
38.3.1.hp的全基因组测序及基因注释
39.采用illumina nextseq 550测序平台对hp进行全基因组测序。细菌基因组dna提取方法同前。二代测序采用amt rapid dna
‑
seq kit for illumina(cistro,china)建库，采用high output v2.5 kit(illumina,usa)试剂盒测序。下机数据采用trimmomatic(v0.39)软件进行质控，后采用spades(v3.13.1)软件进行组装。组装后hp基因组采用quast(v5.0.2)软件进行组装质控评估，对于质控合格的基因组，采用prokka(v1.13)软件进行注释。
40.3.2.hp耐药基因的点突变检测
41.从hp注释后的核苷酸序列文件及氨基酸序列文件中提取待分析的耐药基因的核苷酸序列及耐药蛋白的氨基酸序列。选择待分析的目标基因或蛋白序列导入clc workbench(qiagen，germany)软件，采用classical sequence analysis中的alignment and tree分析模块，比对已知的hp抗生素敏感标准菌株atcc 26695的参考序列，进行基因/蛋白序列特征图谱的分析比对。与atcc 26695相应的基因核苷酸序列相比，若出现基因突变，则在数据库中录入相应基因位点及其核苷酸的突变类型。与atcc 26695的相应蛋白氨基酸序列相比，若出现氨基酸突变，则在数据库中录入相应突变的氨基酸位点及氨基酸的突变类型。
42.3.3.hp中国流行株耐药表型
‑
分子特征关联模型的构建与应用
43.综合hp中国流行株的耐药表型谱、耐药基因突变谱及耐药蛋白突变谱信息，整理获得hp中国流行株耐药基因的分子表征库(图1，表1和表2)。采用pearson相关分析检测耐药基因核苷酸或氨基酸序列突变与耐药表型的相关性，采用p<0.05认为两者存在相关性。由附图2分析可知，多数由其它地区报道的hp耐药分子特征与中国hp流行株的耐药表型并不存在相关性。而采用发明人所整理的hp中国流行株耐药基因的分子表征库，可建立hp中国流行株的耐药基因组分子表征关联模型，采用模型则可根据其它中国流行株的耐药基因突变预测其耐药表型，分析发现存在耐药突变的hp菌株对抗生素的mic高于野生株(图2)。
44.实施例4高通量、快速检测幽门螺旋杆菌耐药性
45.在患者胃黏膜分离获得hp单菌落后，参照前述方法进行全基因组hp全基因组测序，并提取23s rrna基因的核苷酸序列及infb、gyrb、rdxa、rpsu及sodb共5个蛋白的氨基酸序列。对比附表1、附表2检测相应的核苷酸位点及氨基酸位点，读取hp分离株对克拉霉素、左氧氟沙星及甲硝唑三种抗生素为耐药或是敏感。对比现有方案，采用本发明能检出采用旧方案不能检出的hp耐药情况，即本发明提供的检测方案对于耐药的hp中国流行株检出有更高的敏感性(附表3)。
46.表1 hp耐药基因核苷酸突变位点与对应抗生素药物参照表
[0047][0048]
表2 hp耐药蛋白基因氨基酸突变位点与对应抗生素药物参照表
[0049][0050][0051]
表3现有方案与本发明方案对54株hp中国分离株耐药情况的判读准确率对比
[0052][0053]
从上可知，本发明是基于hp中国流行株的基因及耐药谱特征建立的预测模型，本发明主要针对中国株普遍耐受的抗生素药敏表型进行检测，更契合中国流行株的基因组特征情况，具有高通量、高效、准确预测中国hp流行株耐药情况的作用。
[0054]
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。