建立雙重?zé)晒舛縋CR方法檢測蜱傳疏螺旋體

楊小娜，張 琳，侯學(xué)霞，郝 琴

蜱是傳播多種人獸共患病病原體的重要媒介，寄主范圍廣泛，包括哺乳動物、鳥類、兩棲動物和其他脊椎動物。大多數(shù)蜱在幼蜱、若蜱、成蜱3個生長階段中通過吸食宿主血液來維持生長和發(fā)育[1]。硬蜱攜帶的病原體?？梢鹑R姆病、斑點熱、無形體病、埃立克體病、森林腦炎、新疆出血熱、Q熱、新型布尼亞病毒感染等疾病發(fā)生[2]。由于蜱攜帶的病原體較多，一些新的蜱傳病原體不斷發(fā)現(xiàn)，蜱傳疾病在全球的分布范圍不斷擴大，對家畜和人類健康造成重大危害[3]。

蜱傳疏螺旋體根據(jù)引起人類疾病的不同通常分為伯氏疏螺旋體和回歸熱螺旋體兩個分類群。伯氏疏螺旋體(Borreliaburgdorferi，BB)是引起萊姆病的病原體，萊姆病在世界70多個國家均有報道，是北美洲和歐亞大陸國家最常見的蜱傳疾病之一[4]。萊姆病可引起游走性紅斑、關(guān)節(jié)炎和其他神經(jīng)系統(tǒng)癥狀[5]。目前，萊姆病螺旋體可分為20多個基因型，其中至少有5個能引起人類疾病，包括B.burgdorferisensustricto(B.b.s.s)、B.garinii、B.afzelii、B.bavariensis和B.spielmani[6]。Borreliamiyamotoi(BM)于1994年在日本首次分離并被確認為一種新的螺旋體。2011年，俄羅斯首次出現(xiàn)了人類感染病例[7]，隨后在美國、歐洲、日本也發(fā)現(xiàn)了病例[8-10]。從系統(tǒng)進化上講，它屬于回歸熱螺旋體，與大多數(shù)回歸熱螺旋體不同的是它不引起典型的回歸熱，并通過硬蜱傳播，可導(dǎo)致嚴重的疾病，如腦膜腦炎[11]。

在亞洲，全溝硬蜱不僅是傳播B.burgdorferi的重要媒介，而且已被證實能夠傳播B.miyamotoi，白足鼠等嚙齒動物、鳥類等也是兩種螺旋體重要的宿主[5]。至今很多研究報道了B.burgdorferi與B.miyamotoi同時感染硬蜱的病原學(xué)證據(jù)，描述了二者在同一種群的流行狀況，也發(fā)現(xiàn)了二者混合感染的患者[12]。兩種病原具有相同的傳播媒介(硬蜱)和宿主動物，同時人類B.miyamotoi感染有一些與萊姆病難以區(qū)分的臨床表現(xiàn)，如皮膚紅斑、神經(jīng)系統(tǒng)受損癥狀等。在我國，尚未有報道利用熒光定量PCR方法用于檢測B.miyamotoi，因此，研發(fā)一種高效、快速和敏感的檢測方法用于檢測和區(qū)分這兩種病原體是非常必要的。在本研究中我們設(shè)計了一種能同時檢測這兩種病原體的雙重?zé)晒舛縋CR方法，具有良好的靈敏度和特異度，并用這種方法檢測了模擬樣本和新疆采集的200只硬蜱，對該方法進行了初步的應(yīng)用評價。

1 材料與方法

1.1 材料 對照菌株DNA：鉤端螺旋體、土拉弗朗西斯菌、巴爾通體、布魯氏菌、大腸埃希菌均由中國疾病預(yù)防控制中心傳染病預(yù)防控制所相關(guān)科室提供。試劑與儀器：TaKaRa 熒光定量taq酶(TaKaRa Taq Version 2.0 plus dye, RR390)；Roche LightCycler 480實時熒光定量PCR儀；天根DNA提取試劑(DNeasy Blood&Tissue Kit)。

1.2 方 法

1.2.1 引物探針合成 分別根據(jù)B.burgdorferirecA基因、B.miyamotoi甘油磷酸二酯酶生物合成基因(glpQ)設(shè)計引物和探針[13-14]。兩組引物探針見表1，利用Primer-BLAST(https://www.ncbi.nlm.nih.gov/tools/primer-blast/index.cgi)確定其對目的基因的特異性。

表1 多重PCR引物和探針Tab.1 PCR primers and probes

1.2.2 標準品制備B.burgdorferirecA基因擴增片段為GTTCTGCAACATTAACACCTAAAGCTTTTGCATAAACAGGATCAAGAGCA-TGCTCAGCATCAATAAAAGCAGCTATCCCA-CCT(83 bp)，B.miyamotoiglpQ基因擴增片段為ATGCACGACCCAGAAATTGACACAACCACA-AATGTTGCACAATTATTTCCCAATCGAGCT-AGAGAAAACGGACGATATTACGCCACTGA-CTT(78 bp)，擴增片段交由上海生工公司合成質(zhì)粒。將4 μg質(zhì)粒干粉加入1 mL去離子水，得到質(zhì)粒濃度4 μg/mL(4 ng/μL)。每個堿基的平均分子量660，pUC57質(zhì)粒含2 710個堿基，單質(zhì)粒分子量為660×(基因片段堿基數(shù)+2 710)，每μL拷貝數(shù)=(質(zhì)粒濃度/單質(zhì)粒分子量)×10-9×阿伏伽德羅常數(shù)。計算得到兩種質(zhì)?？截悢?shù)均為1.3×109拷貝/μL。

1.2.3 單重?zé)晒舛縋CR

1.2.3.1 靈敏度 對熒光定量PCR反應(yīng)的引物和探針濃度進行優(yōu)化后，確定最佳反應(yīng)體系。首先分別對recA、glpQ基因質(zhì)粒進行單基因靶點的檢測，將標準質(zhì)粒梯度稀釋為109～10-1拷貝/μL，每個濃度設(shè)置4個重復(fù)試驗，確保引物探針的可用性。

1.2.3.2 特異度 利用建立的單重?zé)晒舛縋CR體系，分別檢測鉤端螺旋體、布魯氏菌、土拉菌、巴爾通體、大腸埃希菌DNA。

1.2.4 雙重?zé)晒舛縋CR實驗

1.2.4.2 靈敏度 利用109～10-1拷貝/μL的陽性質(zhì)粒標準品，建立20 μL反應(yīng)體系，每個反應(yīng)中含recA基因上下游引物和探針(10 μmol/L)分別1 μL，glpQ基因上下游引物和探針(10 μmol/L)分別0.8 μL，去離子水0.6 μL，兩種陽性質(zhì)粒各2 μL。每個濃度設(shè)置4個重復(fù)實驗，4個實驗均為陽性結(jié)果則判定為陽性，確定該方法的靈敏度。

1.2.4.3 與單重?zé)晒舛縋CR方法差異性比較 根據(jù)1.2.4.2靈敏度檢測結(jié)果，將雙重與單重實驗的Ct值(熒光信號到達設(shè)定域值時所經(jīng)歷的循環(huán)數(shù))進行比較，利用配對設(shè)計資料的Wilcoxon signed -rank test 符號秩和檢驗，判斷單重與雙重檢測方法的結(jié)果有無明顯差別。

1.2.4.4 標準曲線 利用1.2.4.2的實驗數(shù)據(jù)，根據(jù)樣本濃度對數(shù)和Ct值繪制標準曲線。

1.2.4.5 特異度 利用建立的雙重?zé)晒舛縋CR體系，檢測鉤端螺旋體、布魯氏菌、土拉菌、巴爾通體、大腸埃希菌DNA。

1.2.4.6 濃度變化對檢測結(jié)果的影響 為探索使用雙重?zé)晒舛縋CR方法檢測病原體過程中，其中一種病原體在另一種病原體濃度變化背景下檢測結(jié)果的差異。設(shè)置第1組實驗：recA質(zhì)粒濃度固定為104拷貝/μL，glpQ質(zhì)粒濃度106～102拷貝/μL；設(shè)置第2組實驗：glpQ質(zhì)粒濃度固定為104拷貝/μL，recA質(zhì)粒濃度106～102拷貝/μL。如果在某一種病原體濃度變化時，另一種病原體的Ct值沒有顯著變化，說明在一種病原體濃度變化為背景時，對另外一種病原體的檢測沒有顯著影響。

1.2.4.7 探索其他病原體存在對檢測兩種病原體結(jié)果的影響 利用鉤端螺旋體、布魯氏菌、土拉菌、巴爾通體、大腸埃希菌DNA，設(shè)計如下實驗：

(1)只有recA和glpQ兩種陽性質(zhì)粒。

(2)recA和glpQ兩種陽性質(zhì)粒加上述5種病原體DNA，作為實驗組探索其他病原體存在對兩種病原體檢測的影響。

(3)只有上述5種病原體DNA，不含有recA和glpQ兩種陽性質(zhì)粒。

(4)不含有任何病原體DNA，用去離子水補齊體系。

1.2.5 樣品檢測

1.2.5.1 模擬樣本檢測 將梯度稀釋的標準模板(107～103拷貝/μL)分別取1 μL加入9 μL檢測陰性(巢式PCR和熒光定量PCR反應(yīng)均為陰性)的蜱DNA中，此時濃度變?yōu)?06～102拷貝/μL，分別進行雙重?zé)晒舛縋CR反應(yīng)。通過批內(nèi)及批間變異系數(shù)(CV)評價反應(yīng)的重現(xiàn)性。批內(nèi)重現(xiàn)性：同一批次內(nèi)，106～102拷貝/μL的模擬樣本分別設(shè)置3個平行樣本，分析各個稀釋度反應(yīng)Ct值的變異系數(shù)。批間重現(xiàn)性：以不同批次間(不同時間進行3次實驗)標準模板各個稀釋度反應(yīng)Ct值的變異系數(shù)評價批間重現(xiàn)性。

1.2.5.2 實際樣本檢測 利用本研究設(shè)計的雙重?zé)晒舛縋CR方法檢測采集于新疆200只蜱蟲B.burgdorferi和B.miyamotoi的攜帶情況。每次實驗均設(shè)計陽性對照和陰性對照進行質(zhì)控。

2 結(jié) 果

2.1 單重?zé)晒舛縋CR

2.1.1 靈敏度 用優(yōu)化好的反應(yīng)體系做單重?zé)晒舛縋CR反應(yīng)，取recA基因上下游引物探針濃度為500 nmol/L,glpQ基因上下游引物探針濃度為400 nmol/L。兩種基因能檢測到的最低濃度均為101拷貝/μL，一個反應(yīng)最少能檢測20個基因拷貝。

2.1.2 特異度 利用單重?zé)晒舛縋CR體系，檢測鉤端螺旋體、布魯氏菌、土拉菌、巴爾通體、大腸埃希菌DNA。檢測結(jié)果均為陰性，說明特異性良好。

2.2 雙重?zé)晒舛縋CR

2.2.1 靈敏度 雙重?zé)晒舛縋CR方法對不同濃度梯度的標準品有很好的檢測效果，重復(fù)樣本的Ct值變化很小，重復(fù)性較好(圖1)。recA和glpQ兩種陽性質(zhì)粒能檢測到的最低濃度均為102拷貝/μL，一個反應(yīng)最少能檢測200個基因拷貝，與單重?zé)晒舛縋CR方法相比靈敏度下降了一個數(shù)量級。

AB A：B. burgdorferi recA基因擴增曲線；B： B. miyamotoi glpQ基因擴增曲線圖1 雙重?zé)晒舛縋CR靈敏度實驗Fig.1 Sensitivity experiments of duplex real-time PCR

2.2.2 與單重?zé)晒舛縋CR方法差異性比較 利用雙重和單重方法的Ct平均值，進行配對設(shè)計資料的Wilcoxon signed -rank test符號秩和檢驗。經(jīng)檢驗，50.05，可認為兩種方法測定結(jié)果無差別(表2)。

表2 不同濃度梯度標準品測量值和期望值的比較Tab.2 Comparison of measured values and expected values

2.2.3 標準曲線B.burgdorferirecA基因和B.miyamotoiglpQ基因的濃度對數(shù)分別與檢測Ct值呈線性相關(guān)，R2分別為0.993 7、0.996 2，接近于1，擬合效果好(圖2)。

ABA：B. burgdorferi recA基因標準曲線； B：B. miyamotoi glpQ基因標準曲線圖2 雙重?zé)晒舛縋CR標準曲線Fig.2 Duplex real-time PCR standard curve

2.2.4 特異度 利用雙重?zé)晒舛縋CR體系，檢測鉤端螺旋體、布魯氏菌、土拉菌、巴爾通體、大腸埃希菌DNA。檢測結(jié)果均為陰性，說明特異性良好。

2.2.5 濃度變化影響 第1組實驗：glpQ質(zhì)粒濃度變化時(106～102拷貝/μL)，recA質(zhì)粒(104拷貝/μL)檢測Ct值無明顯差別；第2組實驗：recA質(zhì)粒濃度變化時(106～102拷貝/μL)，glpQ質(zhì)粒(104拷貝/μL)檢測Ct值無明顯差別(圖3)。說明對于這兩種疏螺旋體，其中一種病原基因濃度變化對另一種病原體檢測的結(jié)果無明顯影響。

AB圖3 一種病原濃度變化對另一種病原體檢測結(jié)果(Ct值)的影響Fig.3 Effect of concentration on detection results (Ct values)

2.2.6 其他病原體存在對檢測B.burgdorferi和B.miyamotoi結(jié)果的影響 陰性對照(3)、(4)檢測結(jié)果均為陰性，(1)、(2)BB和BM檢測均為陽性，Ct值差別不大，但是(2)的最大熒光強度比(1)低(圖4)。說明其他多種病原體DNA的存在會降低測量的熒光強度，但對于病原體的檢出沒有影響。

ABA：B. burgdorferi recA基因擴增曲線；B：B. miyamotoi glpQ基因擴增曲線圖4 其它病原體干擾實驗的基因擴增曲線Fig.4 Gene amplification curves when other pathogens interfere with the experiment

2.3 樣本檢測

2.3.1 模擬樣本檢測 通過對各濃度梯度(106～102拷貝/μL)模擬蜱樣本進行檢測，均產(chǎn)生陽性結(jié)果。同一批次內(nèi)模擬蜱樣本各濃度Ct值的變異系數(shù)為0.111～0.406(BB)、0.675～1.301(BM)，各批次間(不同時間進行的3次試驗)各濃度Ct值的變異系數(shù)為0.412～1.238(BB)、0.285～1.149(BM)，說明反應(yīng)具有較好的重現(xiàn)性和穩(wěn)定性(圖5)。

圖5 模擬蜱樣本批內(nèi)和批間變異系數(shù)Fig.5 Intra-assay and inter-assay coefficient of variation of simulated tick samples

2.3.2 實際樣本檢測 在新疆地區(qū)采集到的200只蜱中，B.burgdorferi的recA基因檢測有11個陽性(Ct值<37)，B.miyamotoi的glpQ基因檢測均為陰性。

3 討 論

B.miyamotoi與B.burgdorferi具有相同的傳播媒介和感染宿主，其感染具有相似的臨床癥狀和體征，能導(dǎo)致人類嚴重疾病。我國已在黑龍江省牡丹江地區(qū)和山西省祁縣地區(qū)發(fā)現(xiàn)了蜱感染B.miyamotoi的病原學(xué)依據(jù)，其它地區(qū)尚未有相關(guān)報道[12,15]。在我國，尚無鑒定B.miyamotoi的熒光定量PCR方法的建立。為了能同時檢測和鑒別兩種疏螺旋體，提高鑒定效率，本研究利用B.burgdorferi的recA基因和B.miyamotoi的glpQ基因，建立了一種雙重?zé)晒舛縋CR方法。

本研究首先利用兩種病原陽性質(zhì)粒分別進行了單重?zé)晒舛縋CR實驗，確保每對引物和探針都能擴增出對應(yīng)的病原體，同時確定反應(yīng)進行的最佳引物和探針濃度。為了增強實驗結(jié)果分析的可靠性，在單重實驗的基礎(chǔ)上建立雙重?zé)晒舛縋CR方法，實驗結(jié)果證明兩種病原體的引物探針沒有交叉反應(yīng)，對于不同濃度標準品的檢測具有良好的區(qū)分度。雙重檢測實驗靈敏度和特異度良好，標準曲線顯示Ct值與兩種病原標準品濃度的對數(shù)呈良好的線性相關(guān)關(guān)系。通過配對設(shè)計資料的Wilcoxon signed -rank test 符號秩和檢驗判定單重和雙重方法測得的Ct值結(jié)果無明顯差異。病原體濃度變化和其他病原體DNA的存在對檢測Ct值沒有明顯的影響。我們利用建立的新方法檢測模擬蜱樣本和實際蜱樣本，證明了雙重?zé)晒舛縋CR方法的可行性和實用性。

在多重檢測體系中，多對引物和探針及熒光基團之間會相互干擾，引起實驗結(jié)果靈敏度低、特異性不強的情況，我們使用了兩種方法來避免該問題。一是設(shè)計探針時使每種病原體探針的報告基團熒光發(fā)射最大波長處于不同的光譜范圍，本實驗選用的6-FAM(465-510)、VIC(533-580)兩種熒光染料，其波長處于不同光譜范圍。二是在建立雙重?zé)晒舛縋CR方法前先做了顏色補償實驗，以確保熒光定量PCR儀的檢測通道能區(qū)分每種病原體的探針，顏色補償結(jié)果用于后續(xù)結(jié)果分析。

本研究建立的雙重?zé)晒舛縋CR方法能快速檢測蜱樣本中B.burgdorferi和B.miyamotoi，但是能否用于臨床標本的檢測或者人類疾病的診斷分析還需要進一步的研究。該方法快速、靈敏、特異，為我們提供了一個有價值的檢測工具，可用于病媒檢測和風(fēng)險評估，為疾病控制及生物安全提供了技術(shù)支持，也為臨床檢測提供了思路。