꿀벌 6종 주요 병원체에 대한 초고속 다중 PCR 검출법의 개발

시료 수집과 꿀벌 성봉 및 유충의 핵산 순수 분리

꿀벌 유충 시료는 2016년 7월 경기대 양봉장에서 질병에 감염되지 않은 건강한 시료를 채집하여 사용하였다. 채집한 유충 시료는 -70°C에서 냉동보관 하였으며, 핵산의 분리는 Patho Gene-spin^TM DNA, RNA extraction kit (iNtRON Co., Ltd, Korea)를 사용하여 DNA를 순수분리 하였다. 분리한 핵산은 Biophotometer (Eppendorf Co., Ltd, Germany)를 사용하여 정량한 후, -70°C에 보관하며 실험에 사용하였다.

꿀벌 병원체의 특이 유전자와 꿀벌 housekeeping 유전자

사용한 꿀벌과 꿀벌 병원체 특이 유전자 재조합 DNA들은, 경기대학교 생명과학과 분자생물학 실험실에 확보된 것으로, 각 재조합 DNA들을 각 병원체의 유전자 정보들과 비교하여 각각 특이 염기서열이 존재함을 확인하였으며, 각각의 유전자 정보 및 참고문헌은 Table 1과 같다.

Table 1.

The recombinant DNAs of 6 honeybee pathogens and β-actin gene

6종 꿀벌 질병 진단을 위한 PCR-chip용 프라이머

6종의 꿀벌 병원체 진단을 위해 16S rRNA (P. larvae, M. plutonius,N. apis, N. ceranae), DNA gyrase subunit B (M. plutonius), 18S rRNA (A. apis와 A. flavus)의 특이 염기서열에서 ultra-rapid multi-PCR (UR-multi-PCR)용 프라이머 쌍을 선정하였으며, 꿀벌 핵산의 유무와 PCR 반응의 판별을 위한 대조군으로 사용하기 위하여 housekeeping 유전자인 β-actin 프라이머 쌍을 설계하였다(Table 2). 또한 본 프라이머 쌍들은 Table 1의 재조합 DNA의 특이 염기서열에 근거하여 설계하였으며, 프라이머의 합성은 Bionics사(Korea)에 의뢰 합성하였다.

Table 2.

Primer sets of 6 honeybee pathogens and β-actin primer set for UR-multi-PCR

6종 병원체 프라이머 쌍에 대한 목적 유전자의 검출 한계

각 꿀벌 병원체의 특이 유전자 재조합 DNA 주형(Table 1)을 10⁶ copies에서 10⁰ copies DNA까지 연속 희석하여 초기 주형량에 따른 검출 한계를 확인하였다. Multi-PCR은 GENECHECKER (GENESYSTEM Co., Ltd, Korea)에서 Rapi:Chip (GENESYSTEM Co., Ltd, Korea)을 사용하여 수행하였으며(왕 등, 2016a), 각 well은 5μl의 2×Rapi mix (GENESYSTEM Co., Ltd, Korea)와 1μl씩의 각각의 forward와 reverse 병원체 특이 프라이머 쌍(각 최종 농도 1μM)과 1μl의 연속 희석한 재조합 DNA로 총 10μl PCR 반응액으로 조성하였다. Multi-PCR 반응은 95°C 초기 변성 30초, 95°C 변성 4초, 52°C 중합 4초, 72°C 신장 4초로, 총 40회전 또는 50회전 수행하였다.

질병 진단을 위한 최적 꿀벌 핵산 사용 범위의 확립

진단에 적합한 꿀벌 애벌레의 핵산 사용량을 확립하기 위하여, 꿀벌 애벌레에서 추출한 핵산을 주형으로 β-actin 유전자의 증폭량을 확인하였다. PCR 반응은 5μl의 2×Rapi mix와 1μl씩의 β-actin 151-F(최종 농도 1μM)와 β-actin 151-R(최종 농도 1μM) 프라이머에 꿀벌 애벌레의 핵산(최고 1μg에서 최저 1ng까지 희석), 총 10μl로 조성하여 수행했으며, 이때의 positive control로써 1.48×10⁵ copies의 β-actin 유전자를 포함하는 재조합 DNA를 사용하여 β-actin PCR 산물의 합성 유무와 특이 Tm값을 비교하였다. Multi-PCR 반응은 95°C 초기변성 30초, 95°C 변성 4초, 52°C 중합 4초, 72°C 신장 4초로, 총 40회전 수행하였다.

꿀벌 genomic DNA 존재 시 6종 병원체 프라이머 쌍에 대한 목적 유전자의 검출 한계

12종 병원체(Black queen cell virus, Chronic bee paralysis virus, Deformed wing virus, Israeli acute paralysis virus, Sacbrood virus, Korean sacbrood virus, Paenibacillus larvae, Melissococcus plutonius, Ascosphaera apis, Aspergillus flavus, Nosema apis, Nosema ceranae)의 존재를 확인한 실시간 PCR의 결과에서 어떠한 병원체도 검출되지 않은 꿀벌 애벌레의 핵산 50ng과 6종 병원체의 각각의 재조합 DNA 주형을 10⁶ copies에서 10⁰ copies DNA까지 연속 희석 하여, 각 well당 5μl의 2×Rapi mix와 1μl씩의 각각의 forward 와 reverse 병원체 특이 프라이머 쌍과 총 10μl PCR 반응액으로 조성하였다. Multi-PCR 반응은 95°C 초기변성 30초, 95°C 변성 4초, 52°C 중합 4초, 72°C 신장 4초로, 총 40회전 수행하였다.

6종 질병 초고속 진단을 위한 최적 프라이머 농도

각 질병의 진단 민감도를 높이기 위한 최적 프라이머 농도를 확립하기 위하여, 각 병원체에 해당하는 프라이머는 최종 농도를 4μM, 3μM, 2μM, 1μM, 500nM, 250nM로 조절하여 multi-PCR을 진행하였다. 각 well당 총 10μl로 조성하였으며 꿀벌 애벌레의 핵산 50ng에 각 병원체의 특이 재조합 DNA를 1×10⁵ copies로 하여 5μl의 2×Rapi mix와 위의 각기 다른 프라이머 쌍의 농도로 multi-PCR 반응을 수행하였다. Multi-PCR 반응은 95°C 초기변성 15초, 95°C 변성 4초, 52°C 중합 4초, 72°C 신장 4초로, 총 40회전 수행하였다.

6종 질병 초고속 진단을 위한 최적 혼성 온도

UR multi-PCR에서 각 병원체 민감도를 높이기 위한 최적 혼성 온도의 확립을 위하여, 52~ 66°C 범위에서 측정되는 Ct값을 온도별로 비교 분석하였다. 각 well당 총 10μl로 조성하였으며 꿀벌 애벌레의 핵산 50ng에 각 병원체의 특이 재조합 DNA 10⁵ 분자와 각각의 병원체의 최적의 forward 와 reverse 특이 프라이머 농도를 적용하여 multi-PCR을 95°C 초기변성 15초, 95°C 변성 4초, 52~66°C 중합 4초, 72°C 신장 4초로, 총 40회전 조건으로 수행하였다.

꿀벌 6종 질병 초고속 다중 PCR의 최소 검출 시간 및 검출 한계 확립

6종 꿀벌 질병의 최적 및 최소 검출 시간을 확립하기 위하여, UR multi-PCR의 각 단계에서 초기변성은 95°C에서 15초로, PCR step 중 변성시간은 1초로 설정하였으며, 중합은 앞선 실험의 최적의 온도 조건에서 서로 다른 시간조건 3초, 2초 또는 1초로 하였고, 신장은 72°C, 3초로 설정하여 40회전으로 진행하였다. 사용된 PCR 조성은 총 반응액 10μl에 5μl의 2×Rapi mix와 각 병원체별 최적 농도의 프라이머를 첨가하였으며, 꿀벌 애벌레의 핵산 50ng과 재조합 DNA 10⁵~10³ copies를 사용하였으며, 이때 사용한 각각의 재조합 DNA 분자는 최소 검출 시간 내에 병원체의 검출 한계를 조사하기 위함이었다.

꿀벌 성체를 이용한 초고속 다중 PCR의 유효성 평가

UR multi-PCR을 통한 6종 꿀벌 병원체의 현장 적용실험을 평가하기 위하여, 경기대 양봉장에서 채집된 감염 시료로써, Aspergillus flavus와 Nosema ceranae에 감염된 성봉에서 각기 핵산을 추출하여 PCR을 진행하였다. Patho Gene-spin^TM DNA, RNA extraction kit를 이용하여 꿀벌 한 마리로부터 핵산을 추출해 내었으며, 각각 꿀벌 병원체 검사를 위한 multi-PCR은 well당 총 반응액 10μl로 하여 50ng의 꿀벌 핵산이 사용되었으며, 2×Rapi mix와 각 병원체별 최적화된 프라이머 농도로 조성하여, 최소 검출시간 조건에 따라 진행하였다.

6종 병원체 프라이머 쌍에 대한 목적 유전자의 검출 한계

초고속중합효소연쇄반응(UR-PCR)에서 6종 꿀벌 병원체 프라이머 쌍에 따른 목적 유전자의 검출 한계를 측정하기 위해, 각 병원체에 해당하는 재조합 DNA를 약 1×10⁶에서 1×10⁰ 분자까지 연속 희석하여 이를 주형으로 하였다. Paenibacillus larvae 16S rRNA 유전자를 포함하는 재조합 DNA 주형의 경우, AFB-DC-F1/R1 프라이머는 최소 1.23×10³ 분자까지 검출 가능함을 확인하였고(Fig. 1A), 검출산물의 융점분석에서 특이 Tm (melting Temperature) 값은 88.68±0.35으로 비특이 Tm값과 극명한 차이를 보였으며, 특이Tm값은 이들 검출산물의 유의성을 판별하는 기준이 되었다(Fig. 1B). 또한 다른 꿀벌 병원체의 특이 유전자(Mellisococcus plutonius 16S, Mellisococcus plutonius DNA gyrase, Ascosphaera apis 18S, Aspergillus flavus 18S, Nosema ceranae 16S와 Nosema apis 16S rRNA)와 꿀벌의 housekeeping 유전자(β-actin)가 포함되어 있는 재조합 DNA를 사용한 중합효소 연쇄반응에서 각 프라이머 쌍에 대한 검출 한계는 각각 1.24×10¹, 1.35×10², 1.17×10³, 1.06×10³, 1.03×10², 1.36×10³, 1.48×10¹ 분자였으며(Fig. 2A), 또한 목적 유전자 산물의 특정 Tm값을 통해 증폭의 유의성을 확인하였다(Fig. 2B).

Fig. 1.

Specificity of PCR product according to initial templates of target gene for Paenibacillus larvae. Target gene of Paenibacillus larvae was amplified with AFB-DC-F1/R1 primer pair. Panel A is the fluorescent graph of amplified target gene according to the serially diluted initial templates. Panel B is a peak graph of Tm value. The PCR condition was 95°C, 30 sec of pre-denaturation, 95°C, 4 sec of denaturation, 52°C, 4 sec of annealing, 72°C, 4 sec of polymerization in each cycle for 50 cycles. AFB-DC-F1/R1 primers can detect 1.23 ×103 copies of recombinant plasmid inserted with Paenibacillus larvae 16S rRNA.

Fig. 2.

Ct value and detection limit according to initial templates of target gene. Panel A showed the Ct value of each pathogens according to the initial templates and positive amplifications were classified with considering the Tm values of them (*gy: DNA gyrase). Panel B revealed the detection limit of each primer pairs and range of positive Tm values of each pathogen.

질병 진단을 위한 최적 꿀벌 핵산 사용 범위의 확립

병원체 감염 시료의 초기 총 핵산을 최소화하는 것은 시료로부터 병원체 특이 유전자를 검출해 내는데 있어 불리하며 또한 과도한 핵산의 존재 역시 중합효소 연쇄반응에 영향을 주기 때문에 진단에 적합한 꿀벌 핵산의 범위를 정하는 것은 중요하다. UR-PCR에 적합한 꿀벌 핵산 사용 범위의 확립을 위해, 게놈 DNA (genomic DNA, gDNA)는 꿀벌 애벌레 한 마리(평균 73mg)로부터 평균 1.59μg을 추출하였고, 추출한 핵산은 꿀벌의 housekeeping 유전자인 β-actin 유전자의 검출을 위해 각각 1,000ng, 500ng, 250ng, 100ng, 50ng, 25ng, 10ng, 1ng을 반응에 사용하였다. β-actin 유전자 증폭 형광 그래프 상에서 정규 곡선을 형성하며, 초기 형광 강도의 2배 이상의 최종 형광 값을 갖는 핵산 농도는 100ng, 50ng, 25ng과 10ng이었으며, 25ng의 gDNA를 주형으로 하였을 때 가장 높은 최종 형광 값을, 1,000ng과 500ng에서는 가장 낮은 형광 값을 나타내었다(Fig. 3A). 꿀벌 핵산 농도별 β-actin의 Tm값은 앞선 실험의 81.79±0.84 범위의 유의 값에 분포하였다(Fig. 3B). 또한 100~25ng의 주형에서 Ct값 24로써 가장 빠르게 측정되었으며, 이어서 10ng과 250ng에서 Ct값 25, 1,000ng과 500ng에서는 Ct값이 측정되지 않았다(Table 3). 이는 과량의 핵산(1,000ng과 500ng)을 사용할 경우 목적 유전자의 증폭 유무를 판단하기에 부적합함을 보여 주었다. 따라서 과량의 핵산에 의한 영향이 없으며, 병원체 유전자를 포함하는 최대의 핵산을 사용하기 위하여, 형광그래프 상에서 관측되는 정규 곡선, 최종 형광 값, 그리고 측정된 Ct값을 토대로, PCR 증폭에 최적인 꿀벌 핵산 사용량은 50ng으로 판정하였다.

Fig. 3.

The optimal condition of honey bee nucleic acid to UR-PCR. The amplification of β-actin gene was performed using 1,000-1ng gDNA of honeybee larvae. Panel A is the fluorescent graph of amplification according to initial gDNA and Panel B is peak graph of Tm value. The amplification of β-actin gene with 1,000ng and 500ng initial gDNA was done poorly. Only in 100-10ng gDNA as a template, the amplification curve was formed regularly and the optimal usage of gDNA was 50ng.

Table 3.

The Ct value and Tm value of β-actin gene depending on the concentration of initial templates

꿀벌 gDNA 존재 시 6종 병원체 프라이머 쌍에 대한 목적 유전자의 검출 한계

꿀벌의 gDNA내의 각 병원체 DNA가 함께 존재할 경우 검출 한계를 확인하고자, 각 병원체에 해당하는 재조합 DNA를 1×10⁸에서 1×10⁰ 분자까지 연속 희석하여 50ng의 꿀벌 gDNA와 혼합하여 이를 각각의 UR-PCR에 사용하였다. Paenibacillus larvae 16S rRNA 특이 프라이머쌍은 꿀벌 gDNA 존재 하에서 최소 1.23×10³ 분자까지 검출 가능하였으며(Fig. 4A), 이때 특이적인 증폭에 대하여 Tm은 88.95±0.35로 확인되었다(Fig. 4B). 이외에 Mellisococcus plutonius 16S, Mellisococcus plutonius DNA gyrase, Ascosphaera apis 18S, Aspergillus flavus 18S, Nosema ceranae 16S와 Nosema apis 16S rRNA의 검색에 대하여 각각 1.24×10³, 1.35×10², 1.17×10², 1.06×10³, 1.03×10², 1.36×10³의 분자까지 검출이 가능하였으며, 이때 각 Tm값의 범위는 87.16±0.35, 80.29±0.50, 87.31±0.42, 89.57±0.36, 81.80±0.65, 78.96±1.26으로 나타났으며, 이는 실전시료에서의 각 병원체별 유전자 증폭 유무를 판단하는 특이 지표로써 사용이 가능할 것이다 (Fig. 5B).

Fig. 4.

Specificity of PCR product according to initial templates of target gene for Paenibacillus larvae mixed with honey bee gDNA. The target gene of Paenibacillus larvae was amplified with AFB-DC-F1/R1 primer pair. Panel A is the fluorescent graph of amplified target gene of Paenibacillus larvae according to the serially diluted initial templates. Panel B is a peak graph of Tm value. The PCR condition was 95°C, 30 sec of pre-denaturation, 95°C, 4 sec of denaturation, 52°C, 4 sec of annealing, 72°C, 4 sec of polymerization in each cycle for 40 cycles. AFB-DC-F1/R1 primers (Final conc. 1uM) can detect 1.23×103 copies of recombinant plasmid inserted with Paenibacillus larvae gene with 50 ng genomic DNA of honeybee.

Fig. 5.

Ct value and detection limit according to initial templates of target gene mixed with honeybee gDNA. Panel A showed the Ct value of each pathogens according to the initial templates with gDNA and positive amplifications were classified with considering the Tm values of them (*gy: DNA gyrase). Panel B revealed the detection limit of each primer pairs and range of positive Tm values of each pathogen.

애벌레 한 마리당 회수한 gDNA는 평균 1.59μg에 해당하였고, 50ng의 genomic DNA 존재 하에 Paenibacillus larvae의 16S rRNA 재조합 유전자는 1.23×10³의 분자까지 검출할 수 있었으므로 이를 역으로 계산(1.59μg / 50ng) × detection limit of each pathogen)하였을 때, 1.59μg에 해당하는 애벌레 한 마리에 3.91×10⁴ copies 이상의 병원체가 있을 경우 검출이 가능한 것으로 확인하였다. 나머지 유전자 역시 애벌레 당 회수한 평균 gDNA 1.59μg를 기준으로, Mellisococcus plutonius는 3.94×10⁴, Mellisococcus plutonius (gyrase)는 4.29×10³, Ascosphaera apis는 3.72×10³, Aspergillus flavus는 3.37×10⁴, Nosema ceranae는 3.28×10³, Nesema apis는 각각 4.32×10⁴ copies의 병원체가 존재할 경우 검출 가능한 것으로 산출하였다.

6종 질병 초고속 진단을 위한 최적 프라이머 농도

초고속 PCR에 적합한 6종 질병(7종의 유전자)와 β-actin 유전자 프라이머의 최적 농도를 확립하기 위하여 최종 농도 4μM, 3μM, 2μM, 1μM, 500nM과 250nM의 범위에서 UR-PCR반응을 진행하였다. 3회 반복 시행하여 각 농도별 측정되는 Ct값과 최종 형광 값, melting graph에서 비특이적인 그래프가 형성되지 않는 조건에서 최적 프라이머 농도를 설정하고자 하였다. 그 결과 AFB-DC-F1/R1, EFB-DC-F1/R1, EFB-PCF1/R1, AA-DC-F1/R1, AF-DC-F1/R1은 최종 농도 2μM을 사용하였을 때 민감도가 가장 우수하였으며, No-DC-F2/R1와 NoA-F3(179)/No-R2는 최종 농도 4μM을 사용하는 것이 검출에 가장 유리하였고 또한 β-actin 151-F/R는 최종 농도 1μM에서도 우수한 증폭을 보였다(Table 4).

Table 4.

Ct and range of Tm value according to each concentration of 8 primer pairs

6종 질병 초고속 진단을 위한 최적 중합 온도

UR-PCR을 지향함에 있어 높은 중합 온도는, PCR 단계 간의 온도 상승속도와 하강속도(Ramping time)를 줄여, 전체 반응 시간을 단축시킬 수 있는 중요 요소이다. 따라서 6종 질병(7종의 유전자)와 β-actin 유전자의 증폭이 모두 가능한 높은 중합 온도의 검색이 요구되었다. 이에 중합 온도를 52°C에서 66°C까지 각기 달리하여 UR-PCR을 수행하였으며, 이들 유전자의 검출이 가능한 온도를 확인한 결과, P. larvae 16S, A. flavus 18S와 β-actin gene은 66°C의 중합 온도 에서도 검출이 가능하였으며, 그 다음으로 M. plutonius 16S, A. apis 18S는 64°C 그리고 M. plutonius DNA gyrase, N. apis 16S는 62°C, N. ceranae 16S는 60°C의 중합 온도까지 증폭 가능하였다. 그리하여 본 실험에서 목적하는 UR-PCR 지향을 위해, 각각의 다른 8개 유전자의 검출이 용이한 높은 중합 온도의 설정은 동일분자 수(10⁵ copies)에서의 Ct값 분포를 고려하여 60°C로 설정하였다(Fig. 6).

Fig. 6.

The optimal annealing temperature of 6 species pathogens and β-actin gene. For confirming the optimal annealing temperature of UR multi-PCR, the annealing temperature was controlled at 52°C to 66°C Paenibacillus larvae, Aspergillus flavus, β-actin gene can be detected up to 66°C annealing temperature and Mellisococcus plutonius, Ascosphaera apis gene can be detected up to 64°C Mellisococcus plutonius gyrase (gy: DNA gyrase), Nosema apis can be detected up to 62°C annealing temperature and Nosema ceranae can be detected up to 60°C annealing temperature. As a result, the optimal annealing temperature was 60°C (*). 105 copies of DNA were used as template and the concentration of each primer pairs was based on optimal condition.

꿀벌 6종 질병 초고속 다중 PCR의 최소 검출 시간 및 검출 한계 확립

UR multi-PCR 최소 검출 시간 확립을 위하여, 앞서 확립한 각 프라이머 쌍의 최적 농도와 최적 중합 온도 60°C를 기반으로 하였으며, 각 PCR 단계의 반응시간을 조절함으로써 병원체 최소 검출 시간을 확립 하고자 하였다. 중합 온도를 제외하고, 변성 90~95°C, 신장 70~72°C 사이에서 조절하여 40회전 반응하였으나, 변성과 신장의 온도는 기본 조건인 95°C 이하, 72°C 이하로 조절하였을 때 증폭에 불리함을 확인하였다. 또한 초기 변성은 15초 미만일 때 PCR 효율이 감소함을 확인하였다(자료 미제시). 변성과 중합 단계의 반응시간을 각 1초로 단축하여도 목적 유전자의 증폭에는 큰 변화가 없음이 확인되어, 초기변성 95°C, 15초, 변성 95°C, 1초, 중합 60°C 1초, 신장 72°C로 PCR 조건을 확립하였다. 최종적으로 신장 시간을 1~3초로 하여 꿀벌 애벌레의 핵산 50ng과 각 재조합 DNA 10⁵ copies를 주형으로 PCR을 수행한 결과, 6종의 병원체 유전자의 증폭은 신장 3초에서 검출에 가장 유리한 형광 값 및 Ct 값의 분포를 보였다 (Fig. 7A). 또한 각각의 병원체 유전자의 PCR 산물은 특이 Tm값을 나타냈으며(Fig. 7B), UR multi-PCR의 형광창과 아가로스 전기영동에서도 특이 산물이 합성됨을 확인하였다(Fig. 7C). 이 때 6종 질병(7종의 유전자)에 대하여 검출이 가능한 UR-multi PCR의 최소 검출시간(Ramping time 포함)은 12분 57초였다.

Fig. 7.

Minimum running time for detecting 6 honeybee diseases and β-actin gene. In order to confirm the minimum detection time which can amplify seven kinds of gene and β-actin gene, the polymerization time was controlled from 1 sec to 3 sec. Template DNAs were 105 copies of recombinant DNA. PCR condition was 95°C, 15 sec of pre-denaturation, 95°C, 1 sec of denaturation, 60°C, 1 sec of annealing, 72°C of polymerization in each cycle for 40 cycles. The time of polymerization which can amplify all 8 genes was 3 sec. Accordingly, the total running time of UR multi-PCR was 12 min 57 sec. Panel A is the fluorescent graph of amplification and panel B is peak graph of Tm value of each pathogen. Panel C is a fluorescence panel in window (upper) and agarose gel electrophoresis (below) of amplified DNA.

더 나아가 각각의 6종의 질병(7종의 유전자)에대한 UR multi-PCR의 최소 검출시간(Ramping time 포함)은 12분 57초에서 검출 한계는 10³ copies까지 검출 가능하였다(Table 5). 또한 더 빠른 시간 내에 병원체 검출 가능성을 타진하기 위해 다른 신장 시간 2초를 적용한 전체 소요 시간 12분 17초(Ramping time 포함)에서는 10³ copies의 N. ceranae 16S를 제외한 5종 질병(6종의 유전자)만을, 그리고 신장 시간 1초, 전체 소요 시간 11분 37초(Ramping time 포함)에서는 N. ceranae 16S와 M. plutoniuns 16S를 제외한 5종 질병(5종의 유전자)만을 검출할 수 있었다 (Table 6).

Table 5.

Detection limit of UR multi-PCR and Ct and Tm values

Table 6.

The minimum running time for UR multi-PCR detection

꿀벌 성봉을 이용한 초고속 다중 PCR의 유효성 평가

UR multi-PCR의 유효성을 평가하기 위하여, 질병에 걸린 꿀벌 시료로 부터 gDNA를 분리하였으며, 이를 주형으로 감염된 병원체의 검색을 수행하였다. 꿀벌질병시료로 부터 6종 꿀벌 병원체 중 Aspergillus flavus와 Nosema ceranae의 존재를 확인할 수 있었으며, 이는 해당 꿀벌 한 마리로부터 회수한 gDNA 50ng을 사용하여 UR multi-PCR을 진행하였다(Fig. 8).

Fig. 8.

UR multi-PCR for detecting Aspergillus flavus and Nosema ceranae from a honeybee sample. UR multi-PCR was performed with 50ng gDNA from an infected honeybee sample as initial templates. The PCR condition was 95°C, 15 sec of pre-denaturation, 95°C, 1 sec of denaturation, 60°C, 1 sec of annealing, 72°C, 3 sec of polymerization in each cycle for 40 cycles. Aspergillus flavus, Nosema ceranae and β-actin gene were amplified. Panel A is the fluorescent graph of amplification with gDNA of infected honeybee sample. Panel B is the peak graph of Tm value.

UR multi-PCR 결과에서, Aspergillus flavus와 Nosema ceranae 특이 유전자 증폭과 함께 비 특이적인 증폭이 관찰되고 있으나, 각 산물의 융점분석에서 특이 병원체 유전자의 증폭과 비특이적 유전자 증폭을 쉽게 구분할 수 있었다. 측정된 Aspergillus flavus와 Nosema ceranae의 Tm값은 각각 89.21°C, 82.07°C로, 각 병원체의 양성 증폭의 Tm값 범위(89.57±0.36°C과 81.80±0.65°C) 안에 해당함으로 해당 병원체의 증폭임을 판정할 수 있었다(Table 7).

Table 7.

Ct and Tm value of UR multi-PCR detection with honey sample

본 실험에서 확인한 두 가지 질병 이외 나머지 4가지 질병에 대해서 형광 그래프 및 융점 분석을 통한 질병 진단이 가능하였으며, 성봉과 같은 천연 시료 사용에 따른 비 특이적 증폭현상이 관찰되었지만, 모두 전형적 비특이적 증폭에 따른 온도변화, 즉 온도변화의 곡선이 약하게(dF/dT), 다수의 Tm값이 형성되며, 특이 유전자 증폭과는 쉽게 구분될 수 있었다(Fig. 8B).

본 PCR법을 양봉 현장에서 보다 쉽게 사용하기 위하여, 성봉 또는 유충으로부터 보다 간소화된 핵산 추출 방법의 개발이 필요할 것으로 보며, 이런 간편한 핵산 분리법과 본 초고속 PCR법이 연합될 수 있다면, 바로 현장에서 유전자 검사를 빠르게 수행할 수 있을 것이다. 더 나아가 본 연구에서 적용한 DNA기반 병원체뿐만 아니라 꿀벌의 RNA 바이러스에 대하여도 본 연구와 같은 신속한 검출법 개발로 연구가 확장될 것으로 기대하는 바이다.