方法
检测orfA基因拷贝数鉴定布鲁氏菌种型的实时荧光定量PCR方法研究
中华地方病学杂志, 2023,42(5) : 409-413. DOI: 10.3760/cma.j.cn231583-20220315-00079
摘要
目的

建立实时荧光定量PCR方法检测布鲁氏菌基因组内IS711转座酶基因(orfA基因)的拷贝数,并应用于布鲁氏菌种与生物型的鉴定。

方法

构建基于Taqman探针法实时荧光定量PCR技术的orfA基因拷贝数的检测体系。设计bcsp31基因和orfA基因的引物和探针,应用实时荧光定量PCR方法同时检测同一菌株相同DNA浓度下的bcsp31基因和orfA基因含量,获得2个基因的循环数(CT值),再根据bcsp31基因和orfA基因CT值的差异,换算出待检测布鲁氏菌菌株基因组内orfA基因的拷贝数。同时,将布鲁氏菌16M菌株DNA进行2倍递减稀释,验证检测体系的稳定性。

结果

当16M菌株DNA浓度相差2倍时,实时荧光定量PCR方法同时检测bcsp31基因和orfA基因,测得的CT值差值均值均为1.00,95%置信区间为0.95~1.05,标准差为0.17,变异系数为0.17。应用该检测体系检测30株布鲁氏菌菌株DNA中orfA基因的拷贝数,发现羊种生物1~3型菌株分别有6、9和7个拷贝数;猪种生物2型菌株有10个拷贝数,与其他4个猪种生物型(1、3~5)菌株拷贝数不同;绵羊附睾种菌株有37个拷贝数;8个牛种生物型(1~7、9)菌株拷贝数稳定在5~6个。

结论

成功建立了一种检测布鲁氏菌orfA基因拷贝数的实时荧光定量PCR方法,该方法能够鉴定部分布鲁氏菌种与生物型菌株。

引用本文: 田国忠. 检测orfA基因拷贝数鉴定布鲁氏菌种型的实时荧光定量PCR方法研究 [J] . 中华地方病学杂志, 2023, 42(5) : 409-413. DOI: 10.3760/cma.j.cn231583-20220315-00079.
参考文献导出:   Endnote    NoteExpress    RefWorks    NoteFirst    医学文献王
扫  描  看  全  文

正文
作者信息
基金 0  关键词  0
English Abstract
评论
阅读 0  评论  0
相关资源
引用 | 论文 | 视频

版权归中华医学会所有。

未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。

除非特别声明,本刊刊出的所有文章不代表中华医学会和本刊编委会的观点。

布鲁氏菌科包括苍白杆菌属和布鲁氏菌属等[1,2]。依据布鲁氏菌菌株的表型特征和系统发育关系,目前确认的布鲁氏菌属包括12个种:牛种布鲁氏菌(Brucella abortus)、羊种布鲁氏菌(Brucella melitensis)、猪种布鲁氏菌(Brucella suis)、犬种布鲁氏菌(Brucella canis)、绵羊附睾种布鲁氏菌(Brucella ovis)、沙林鼠种布鲁氏菌(Brucella neotomae)、Brucella papionisBrucella cetiBrucella pinnipedialisBrucella microtiBrucella inopinataBrucella vulpis[1,2]。布鲁氏菌生物学鉴定试验包括菌株生长CO2需求试验、硫化氢产生试验、噬菌体裂解试验、硫堇和碱性复红染料抑菌试验,以及单项特异性血清(A、M和R血清)凝集试验等[3]。生物学鉴定试验操作复杂、耗时,并且需要在三级生物实验室进行。而PCR技术以其安全、快速、灵敏度高和特异性强等特点,广泛应用于布鲁氏菌属、种与生物型的鉴定,其中以bcsp31基因为目标基因设计的布鲁氏菌属PCR鉴定方法最为经典[4]。布鲁氏菌种与生物型的鉴定主要有AMOS-PCR方法[5]、多重PCR方法[6]和多位点串联重复序列分析方法(MLVA)[7]。IS711插入序列是布鲁氏菌特有的核苷酸序列,Southern印迹杂交技术是检测IS711插入序列拷贝数最好的方法[8],但是其需要的DNA量较大,且费时费力,无法满足现代检测体系高通量的要求。本研究基于Taqman探针法实时荧光定量PCR技术,设计一种快速、特异、简便和低费用的IS711转座酶基因(orfA基因)拷贝数检测体系,并应用于布鲁氏菌的种与生物型鉴定。

1 材料与方法
1.1 实验菌株和菌株DNA:

实验菌株包括30株布鲁氏菌菌株(其中有27株布鲁氏菌参考菌株、3株临床分离菌株)和8株非布鲁氏菌菌株。参考菌株包括6个种19个生物型的19株菌(牛种8个生物型,猪种5个生物型,羊种3个生物型,以及绵羊附睾种、犬种和沙林鼠种各1个生物型),6株疫苗株和2株粗糙型菌株;临床分离菌株包括1株猪种生物3型菌株(10062)和2株羊种生物2型菌株(79054和17021)。其他非布鲁氏菌菌株包括钩端螺旋体、血清型O∶9小肠结肠炎耶尔森菌、霍乱弧菌、血清型O∶157大肠埃希菌、结核分枝杆菌、嗜麦芽窄食单胞菌、炭疽芽孢杆菌和伤寒沙门菌。30株布鲁氏菌菌株DNA和8株非布鲁氏菌菌株DNA来源于中国疾病预防控制中心传染病预防控制所,- 80 ℃冰箱保存。

1.2 主要仪器与试剂:

DNA微量分光光度测定仪(Spectrophotometer NanoDrop 1000,美国NanoDrop公司),实时荧光定量PCR仪(LightCycler 4800Ⅱ,美国Roche公司)。荧光定量PCR反应酶复合物Premix Ex TaqTM(Probe qPCR,日本TaKaRa公司)。

1.3 实时荧光定量PCR引物和探针:

bcsp31基因的引物和探针序列参照文献[9],正向引物B-F:5'-GCTCGGTTGCCAATATCAATGC-3',反向引物B-R:5'-GGGTAAAGCGTCGCCAGAAG-3',探针B-Probe:5'-AAATCTTCCACCTTGCCCTTGCCATCA-3′,其中5'端用荧光基团FAM标记,3'端用淬灭基团BHQ-1标记。orfA基因的引物和探针,依据orfA基因核苷酸序列,采用DNAstar(版本V8.1.3)软件进行设计,正向引物O-F:5'-AACAGCATGCAGCTTGGTCGTC-3',反向引物O-R:5'-ACCATATCGAAAGTCCACGCAG AT-3',探针O-Probe:5'-TCCACCGCGCGAGCGACC GATGC-3',其中5'端用荧光基团HEX标记,3'端用淬灭基团BHQ-1标记。引物和探针均由生工生物工程(上海)股份有限公司合成。

1.4 实时荧光定量PCR反应体系和反应条件:

将-80 ℃冰箱保存的30株布鲁氏菌菌株DNA和8株非布鲁氏菌菌株DNA取出,采用DNA微量分光光度测定仪测定DNA浓度,使用双蒸水将菌株DNA浓度均稀释为0.1 ng/μl,以供本实验使用。实时荧光定量PCR反应体系:酶复合物Premix Ex TaqTM 10.0 μl,正向引物、反向引物和探针(10 μmol/L)各0.4 μl,DNA模板2.0 μl,补充双蒸水至20.0 μl。实时荧光定量PCR反应条件:95 ℃预变性2 min;95 ℃变性5 s,60 ℃退火20 s,40个循环,在退火阶段检测荧光信号以及循环数(CT值),平行测定4次。bcsp31基因和orfA基因实时荧光定量PCR方法检测的反应体系和反应条件相同。使用实时荧光定量PCR仪进行数据分析,计算CT值的平均值及标准差和标准误。

1.5 方法验证:

①特异性试验:以非布鲁氏菌菌株DNA为模板,应用实时荧光定量PCR方法检测orfA基因和bcsp31基因。②稳定性试验:测定布鲁氏菌16M菌株DNA浓度,并进行2倍递减稀释,以稀释后的菌株DNA为模板,应用实时荧光定量PCR方法检测orfA基因和bcsp31基因。③浓度差值、CT值与拷贝数换算:操作同稳定性试验,4次平行测定后,将测得的CT值进行DNA浓度、拷贝数关系换算,获得拷贝数与CT值的换算公式,计算orfA基因拷贝数。同时自NCBI网站查找已公布布鲁氏菌全基因组内染色体Ⅰ和Ⅱ合计的orfA基因拷贝数(全基因组测序),与本研究结果进行比较。

1.6 统计分析:

采用Excel 2017软件进行数据的处理与分析。

2 结果
2.1 实时荧光定量PCR方法特异性试验:

基于orfA基因和bcsp31基因,应用实时荧光定量PCR方法同时检测钩端螺旋体、血清型O∶9小肠结肠炎耶尔森菌、霍乱弧菌、血清型O ∶ 157大肠埃希菌、结核分枝杆菌、嗜麦芽窄食单胞菌、炭疽芽孢杆菌和伤寒沙门菌共8种细菌DNA,均未检出特异荧光曲线(无CT值)。

2.2 实时荧光定量PCR方法稳定性试验:

布鲁氏菌16M菌株DNA的初始浓度为56.2 ng/μl,当DNA浓度相差2倍时,测得的CT值差值均值均为1.00,95%置信区间为0.95~1.05,标准差为0.17,变异系数为0.17。并且,在检测的30株布鲁氏菌菌株中,bcsp31基因CT值的总体标准差为0.04,标准误为0.02。

2.3 菌株DNA浓度差值与实时荧光定量PCR方法测得CT值及基因拷贝数的关系:

布鲁氏菌16M菌株DNA 2倍递减稀释后,经过4次平行测定,获得了拷贝数与CT值的换算公式:n = ROUNDDOWN(POWER(2,(CTbcsp31 - CTorfA)/1)),其中,n为orfA基因拷贝数,ROUNDDOWN为取整数,POWER为乘幂,CTbcsp31 - CTorfA为实时荧光定量PCR方法分别检测同一菌株同一浓度DNA的bcsp31基因和orfA基因CT值的差值,1为系数。

2.4 布鲁氏菌菌株orfA基因的拷贝数检测结果:

应用实时荧光定量PCR方法同时检测同一菌株同一浓度DNA的两个基因(bcsp31和orfA),每一菌株进行4次平行测定,获得CT值均值,按照公式n = ROUNDDOWN(POWER(2,(CTbcsp31 - CTorfA)/1))进行orfA基因的拷贝数换算,结果见表1。不同种与生物型的布鲁氏菌具有独特的orfA基因拷贝数,其中具有种与生物型区别的是羊种生物2型菌株,具有9个拷贝数,可与羊种生物1型(6个拷贝数)和3型菌株(7个拷贝数)鉴别;猪种生物2型菌株具有10个拷贝数,能够与其他4个猪种生物型菌株鉴别;绵羊附睾种菌株具有37个拷贝数;8个牛种生物型菌株拷贝数稳定在5~6个。同时,实时荧光定量PCR方法检测结果与全基因组测序结果比较,稳定性较好。见表2

点击查看表格
表1

实时荧光定量PCR方法检测布鲁氏菌bcsp31基因和orfA基因的循环数(CT值)与拷贝数关系表

表1

实时荧光定量PCR方法检测布鲁氏菌bcsp31基因和orfA基因的循环数(CT值)与拷贝数关系表

菌株编号菌株名称种与生物型bcsp31基因orfA基因CT值差值corfA基因拷贝数d
CT值均值aCT值标准差bCT值均值aCT值标准差b
1544A牛种生物1型26.630.0324.060.092.576
286/8/59牛种生物2型26.550.0324.170.022.385
3Tulya牛种生物3型26.670.0924.410.052.265
4292牛种生物4型27.130.0624.530.052.606
5B3196牛种生物5型27.170.0324.790.032.385
6870牛种生物6型26.400.0723.710.022.696
763/75牛种生物7型26.310.0323.600.022.716
8C68牛种生物9型26.400.0723.980.062.425
916M羊种生物1型27.180.0324.520.022.666
1063/9羊种生物2型27.560.0324.430.043.139
11Ether羊种生物3型27.670.0524.790.082.887
121330S猪种生物1型26.790.0524.150.052.646
13Thomson猪种生物2型27.140.0523.840.053.3010
14686猪种生物3型27.070.0724.370.072.706
1540猪种生物4型26.310.0324.130.072.185
16513猪种生物5型25.130.0122.420.072.716
17RM6/66犬种27.010.0524.030.062.988
1863/290绵羊附睾种27.140.0521.940.055.2037
195K33沙林鼠种26.690.0424.100.042.596
20M5羊种疫苗株25.780.0322.920.012.867
21M28羊种疫苗株26.920.0224.070.032.857
22REV1羊种疫苗株27.150.0524.740.052.415
23A19牛种疫苗株25.840.0223.320.022.526
24104M牛种疫苗株26.050.0423.520.012.536
25S2猪种疫苗株25.840.0122.960.022.887
2645/20牛种粗糙型27.380.0224.710.022.676
27B115羊种粗糙型25.880.0423.270.012.616
2810062猪种生物3型27.440.0424.950.062.496
2979054羊种生物2型25.910.0522.890.053.028
3017021羊种生物2型27.060.0624.070.052.998

注:a4次平行测定获得的CT值均值;b4次平行测定获得的CT值标准差;c实时荧光定量PCR方法检测bcsp31基因和orfA基因CT值的差值;d根据拷贝数与CT值的换算公式计算得出;28~30号菌株为临床分离菌株

点击查看表格
表2

实时荧光定量PCR方法检测实验菌株orfA基因拷贝数与全基因组测序菌株orfA基因拷贝数的比较

表2

实时荧光定量PCR方法检测实验菌株orfA基因拷贝数与全基因组测序菌株orfA基因拷贝数的比较

菌株名称实验菌株对照菌株
种与生物型CT值aorfA基因拷贝数borfA基因拷贝数c染色体Ⅰ染色体Ⅱ
544A牛种生物1型2.5766(A13334)CP003176.1CP003177.1
86/8/59牛种生物2型2.3856NZ_CP007765.1CP007764.1
Tulya牛种生物3型2.265---
292牛种生物4型2.606---
B3196牛种生物5型2.385---
870牛种生物6型2.6965NZ_CP007709.1NZ_CP007710.1
63/75牛种生物7型2.716---
C68牛种生物9型2.4255NZ_CP007705.1NZ_CP007706.1
16M羊种生物1型2.6666CP007763.1CP007762.1
63/9羊种生物2型3.13910NZ_CP007789.1NZ_CP007788.1
Ether羊种生物3型2.8876CP007760.1CP007761.1
1330S猪种生物1型2.6466NC_004310.3CP002998.1
Thomson猪种生物2型3.301013(23445)CP000911.1CP000912.1
686猪种生物3型2.7064NZ_CP007719.1NZ_CP007718.1
40猪种生物4型2.185---
513猪种生物5型2.716---
RM6/66犬种2.9885(52141)CP003174.1CP003175.1
63/290绵羊附睾种5.203738(25840)CP000708.1CP000709.1
5K33沙林鼠种2.596---

注:a实时荧光定量PCR方法检测bcsp31基因和orfA基因循环数(CT值)的差值;borfA基因实时荧光定量PCR方法检测获得的拷贝数;c已公布布鲁氏菌全基因组内染色体Ⅰ和Ⅱ合计的orfA基因拷贝数;括号内为全基因组测序菌株;"-"为无数据

3 讨论

随着病原菌不断被发现与研究,布鲁氏菌的宿主、种与生物型逐渐被更新,其分类与鉴别越来越多样化。国外报道了可以鉴别牛种、羊种、猪种、犬种、绵羊附睾种、沙林鼠种、Brucella ceti(鲸鱼和海豚)、Brucella pinnipedialis(鳍足类)8个种[10],以及鉴定猪种5个生物型的多重PCR方法[11];目前最常使用的MLVA技术[7]和多位点序列分型(MLST)技术[12]同样具有布鲁氏菌种与生物型鉴别能力,但是这些方法难以鉴别犬种与某些猪种生物型菌株,不能区别羊种生物2型和3型菌株。IS711插入序列在布鲁氏菌基因组内具有多个拷贝数[13],在分析布鲁氏菌不同种与生物型菌株全基因组DNA时发现,与IS711插入序列密切相关的orfA基因,同样存在着多个拷贝数,与IS711插入序列的拷贝数具有对应的关系。

bcsp31基因在布鲁氏菌菌株中具有属特异性和单拷贝数的特征,而orfA基因具有多个拷贝数。利用这两个基因(bcsp31和orfA)经实时荧光定量PCR方法测得CT值的差值,能够换算出orfA基因的拷贝数,该方法能够区分3个羊种生物型菌株:羊种生物1型菌株有6个拷贝数,羊种生物2型菌株有9个拷贝数,羊种生物3型菌株有7个拷贝数。对于绵羊附睾种菌株,临床分离的菌株很少,目前没有快捷的检测方法对其进行鉴定,实时荧光定量PCR方法发现其orfA基因有37个拷贝数,与其他菌种具有明显差别。猪种生物2型菌株有10个拷贝数,与其他猪种生物型具有明显的区别。

对布鲁氏菌疫苗株而言,牛种疫苗株A19和104M与牛种生物1型菌株orfA基因的拷贝数一致,而猪种疫苗株S2有7个拷贝数,与猪种生物1型菌株(6个拷贝数)有差别;羊种疫苗株M5和M28皆有7个拷贝数,羊种疫苗株REV1有5个拷贝数,而羊种生物1型菌株有6个拷贝数,三者之间可以进行鉴别。

本实验共测试了3株临床分离菌株,其中猪种生物3型菌株1株,其orfA基因的拷贝数与参考菌株(686)完全一致。另外2株羊种生物2型菌株测得的orfA基因有8个拷贝数,而实验室保存的羊种生物2型参考菌株(63/9)有9个拷贝数,全基因组测序羊种生物2型菌株(63/9)有10个拷贝数,三者间的差异是否为传代过程中菌株的基因发生了变化,有待于进一步研究。

由于实验室没有保存Brucella cetiBrucella pinnipedialisBrucella microtiBrucella inopinata等其他布鲁氏菌菌种,因此本实验未能测试这些菌种orfA基因的拷贝数,未来将进行相关的研究,以完善所建立的方法。

利益冲突
利益冲突

所有作者声明无利益冲突

参考文献
[1]
LeclercqSO, CloeckaertA, ZygmuntMS. Taxonomic organization of the family Brucellaceae based on a phylogenomic approach[J]. Front Microbiol, 2019, 10: 3083. DOI: 10.3389/fmicb.2019.03083.
[2]
Soler-LlorénsPF, QuanceCR, LawhonSD, et al. A Brucella spp. isolate from a Pac-man frog (Ceratophrys ornata) reveals characteristics departing from classical Brucellae[J]. Front Cell Infect Microbiol, 2016, 6: 116. DOI: 10.3389/fcimb.2016.00116.
[3]
刘秉阳孙玺尚德秋,.布鲁氏菌病学[M].北京:人民卫生出版社, 1989: 6-70.
LiuBY, SunX, ShangDQ, et al. Brucellosis[M]. Beijing: People's Medical Publishing House, 1989: 6-70.
[4]
BailyGG, KrahnJB, DrasarBS, et al. Detection of Brucella melitensis and Brucella abortus by DNA amplification[J]. J Trop Med Hyg, 1992, 95(4): 271-275.
[5]
BrickerBJ, HallingSM. Differentiation of Brucella abortus bv. 1, 2, and 4, Brucella melitensis, Brucella ovis, and Brucella suis bv. 1 by PCR[J]. J Clin Microbiol, 1994, 32(11): 2660-2666. DOI: 10.1128/jcm.32.11.2660-2666.1994.
[6]
李振军侯雪新孙渭歌,.布鲁氏菌种及种内部分生物型快速PCR鉴定方法的建立及应用研究[J].中华微生物学和免疫学杂志, 2013, 33(2): 133-137. DOI: 10.3760/cma.j.issn.0254-5101.2013.02.013.
LiZJ, HouXX, SunWG, et al. Development and application of one-step polymerase chain reaction (PCR) for rapid identification of Brucella species and some biovars[J]. Chin J Microbiol Immunol, 2013, 33(2): 133-137. DOI: 10.3760/cma.j.issn.0254-5101.2013.02.013.
[7]
聂守民罗波艳孙养信,. 2020年陕西省布鲁氏菌MLVA分型研究及流行病学分析[J].中华地方病学杂志, 2022, 41(3): 180-185. DOI: 10.3760/cma.j.cn231583-20210511-00154.
NieSM, LuoBY, SunYX, et al. MLVA typing and epidemiological analysis of Brucella in Shaanxi Province, 2020[J]. Chin J Endemiol, 2022, 41(3): 180-185. DOI: 10.3760/cma.j.cn231583-20210511-00154.
[8]
田国忠郑敏崔步云,.犬种布鲁菌病病原学和遗传特征分析[J].中华地方病学杂志, 2016, 35(10): 707-712. DOI: 10.3760/cma.j.issn.2095-4255.2016.10.002.
TianGZ, ZhengM, CuiBY, et al. Analysis of pathogens and genetic characteristics of Brucella canis[J]. Chin J Endemiol, 2016, 35(10): 707-712. DOI: 10.3760/cma.j.issn.2095-4255.2016.10.002.
[9]
孙宗祥姜勇波张丽娜,.荧光定量PCR方法检测人群血样中的布鲁杆菌DNA[J].中华地方病学杂志, 2015, 34(11): 808-812. DOI: 10.3760/cma.j.issn.2095-4255.2015.11.007.
SunZX, JiangYB, ZhangLN, et al. Detection of Brucella spp.DNA in patients serum samples by real-time PCR[J]. Chin J Endemiol, 2015, 34(11): 808-812. DOI: 10.3760/cma.j.issn.2095-4255.2015.11.007.
[10]
López-GoñiI, García-YoldiD, MarínCM, et al. Evaluation of a multiplex PCR assay (Bruce-ladder) for molecular typing of all Brucella species, including the vaccine trains[J]. J Clin Microbiol, 2008, 46(10): 3484-3487. DOI: 10.1128/JCM.00837-08.
[11]
López-GoñiI, García-YoldiD, MarínCM, et al. New Bruce-ladder multiplex PCR assay for the biovar typing of Brucella suis and the discrimination of Brucella suis and Brucella canis[J]. Vet Microbiol, 2011, 154(1-2): 152-155. DOI: 10.1016/j.vetmic.2011.06.035.
[12]
李文博赵士程高丹,. 2014-2020年四川省人间布鲁氏菌病分离菌株多位点序列分型研究[J].中华地方病学杂志, 2023, 42(1): 24-29. DOI: 10.3760/cma.j.cn231583-20211214-00405.
LiWB, ZhaoSC, GaoD, et al. Multi-locus sequence typing of human Brucella isolated in Sichuan Province from 2014 to 2020[J]. Chin J Endemiol, 2023, 42(1): 24-29. DOI: 10.3760/cma.j.cn231583-20211214-00405.
[13]
HinićV, BrodardI, ThomannA, et al. IS711-based real-time PCR assay as a tool for detection of Brucella spp. in wild boars and comparison with bacterial isolation and serology[J]. BMC Vet Res, 2009, 5: 22. DOI: 10.1186/1746-6148-5-22.
 
 
展开/关闭提纲
查看图表详情
回到顶部
放大字体
缩小字体
标签
关键词