KR102710301B1 - 전장유전체 증폭을 위한 치쿤군야 바이러스 범용 프라이머 세트 및 이를 이용한 진단 키트 - Google Patents

Abstract

본 발명은 치쿤군야 바이러스(Chikungunya virus) 감염 여부를 판별하기 위한 전장 유전체 서열 획득 방법 및 이의 용도에 관한 것이다. 본 발명의 방법, 키트, 조성물은 치쿤군야 바이러스의 유전체 서열을 증폭할 수 있는 범용 프라이머 세트를 포함하고, 이에 따라 시료 내 목적 치쿤군야 바이러스의 존재 유무 및 관련 질환을 간편하고 빠르게 확인할 수 있다. 또한 전장 유전체 분석을 통해 치쿤군야 바이러스의 변이여부를 판별하여 치쿤군야 바이러스 동정에도 효과적으로 사용될 수 있다.

Description

전장유전체 증폭을 위한 치쿤군야 바이러스 범용 프라이머 세트 및 이를 이용한 진단 키트 {Chikungunya virus universal primer sets for whole genome amplification method and diagnosis kit}

본 발명은 치쿤군야 바이러스(Chikungunya virus)의 전장유전체를 증폭할 수 있는 프라이머 세트를 활용하여 미량의 시료에서 바이러스 유전체를 증폭하고 염기서열을 분석하여 전장 유전체 정보를 얻어내는 기술에 관한 것이다.

치쿤군야 바이러스(Chikungunya virus, CHIKV)는 Alphavirus 속과 Togaviridae family에 속하는 positive-sense single-stranded RNA 바이러스이다^1,2. 이 바이러스는 3개의 비번역 영역(UTR), 4개의 비구조 단백질(NSP) 코딩 서열 및 5개의 구조 단백질 코딩 서열을 포함하는 약 12kb 길이의 게놈을 가지고 있다^2,3. 1952-1953 탄자니아에서 대규모 열병으로 처음 인식되었다^4,5. CHIKV는 주요 곤충 매개체 중 하나인 이집트숲모기(Aedes aegypti)에 의해 아프리카에 광범위하게 퍼졌다^6,7. 동부 중앙 남부 아프리카(ECSA)와 일부 서아프리카(WA) 계통이 확인되었다. 동남아시아의 발병은 2004년 이전에 아시아 계통으로 자주 관찰되었다. 2004년 인도양 섬에서 ECSA 계통의 E1-A226V 변이체(E1 외피 당단백질의 위치 226에 있는 단일 아미노산 돌연변이)의 기초는 두 번째 주요 매개체인 Aedes albopictus의 높은 감수성을 유도했다⁹. 인도에서 동남아시아로 거대하고 반복적으로 발생했으며 2005-2007년 동안 아프리카와 일부 유럽에서도 발생했다^9-11. 이 변이체는 ECSA의 하위 계통으로 첫 번째 ECSA 유전자형-인도양 계보(IOL)로 클러스터링되었다^9,12,13. 2013년에 세인트 마틴의 카리브 섬은 아시아계 CHIKV 사례를 보고했으며 나중에 남미 국가와 미국으로 확대되었다^14-16. 재출현은 오늘 현재까지 지속적으로 제기되었다^17,18. 이전에 사례가 없었던 새로운 벡터 또는 새로운 국가를 암시할 가능성이 있다^19-21.

CHIKV 재출현의 위험에서 백신 및 치료 개발, 역학 및 진화 유전체학의 측면에서 전체 게놈 시퀀싱(WGS)의 중요성이 증가했다²². Whole Genome Sequencing은 완전한 염기서열 분석을 위한 염기서열분석이다. 그것은 샷건 접근 방식으로 전체 게놈 무작위 시퀀싱의 적용에서 발생한다²³. 무작위 크기의 게놈 시퀀싱을 통해 각 게놈 단편('읽기'라고 함)은 시퀀스에 반복적인 영역을 가진다. 읽기(read)는 이러한 반복 영역을 식별하고 서로를 조합하여 하나의 연속 시퀀스를 만든다. 결과적으로 전체 게놈이 시퀀싱된다^23,24. WGS는 이제 NGS 플랫폼에서 수행할 수 있다²⁵.

이전에 대규모 병렬 시퀀싱으로 알려진 NGS(next generation sequencing)는 Sanger 시퀀싱 이후에 등장한 대용량 시퀀싱 기술의 총칭이다. NGS는 실행 시간, 비용 및 출력 데이터를 개선하기 위해 Sanger 시퀀싱의 다양한 접근 방식을 기반으로 개발되었다. short-read 또는 long-read를 분석하기 위한 여러 플랫폼이 있다^26-28. Illumina는 short-read NGS28을 나타내는 플랫폼 중 하나이다. Illumina 플랫폼은 2개의 가역 어댑터의 브리지 증폭으로 페어 엔드 시퀀싱을 수행한다. 서로 다른 어댑터 서열이 단일 가닥 DNA 단편의 각 끝과 연결되어 플로우 셀 표면에 부착된다. 다른 하나는 하나의 어댑터를 부착한 후 표면의 상보적인 올리고에 결합한다. 그런 다음 두 개의 어댑터에 의해 유도된 브리지 모양에 따라 새로운 가닥이 합성된다. 각 형광 연결 염기는 합성 주기 동안 다른 신호를 방출한다. 방출된 신호는 DNA 서열을 결정한다²⁹. ONT(Oxford Nanopore Technology)는 PacBio 플랫폼을 탑재한 Long-Read NGS의 대표적인 플랫폼이다. ONT는 최대 수십만 kb의 읽기 길이를 처리할 수 있다²⁸. 단일 단백질 나노포어를 이용하여 나노포어 시퀀싱을 수행한다. 나노포어는 플로우 셀의 지질 이중층에 위치하며 단일 가닥 DNA 단편이 통과할 때 염기의 전기 신호를 읽는다. 각 염기는 다른 잔류 이온 전류를 제공하며 그 순서는 시퀀스를 나타낸다³⁰. 이러한 NGS 플랫폼의 등장은 WGS뿐만 아니라 RNA 시퀀싱, 메타게노믹스, 질병 전파 모니터링(추적) 및 생물 정보학에서도 우수한 연구를 가능하게 했다^31-36.

본 발명은 빠른 차세대 시퀀싱(NGS)을 위한 CHIKV 범용 프라이머 개발을 목표로 한다. 한국에서 분리되어 이전에 부분 또는 구조 CDS에 대해 연구된 균주 NCCP43132(Strain KNIH/2009/77)를 사용하였다. 첫째, 프라이머는 ECSA(IOL 포함) 및 아시아 계통의 전체 서열을 포함하도록 가능한 한 설계되었다. 다음으로, 우리는 전체 게놈 서열에 대해 이중 가닥 cDNA를 사용하여 illumine shotgun sequencing을 수행했다. 앰플리콘의 빠른 시퀀싱은 Oxford Nanopore Technique에 의해 수행되었으며 시퀀스는 CLC 워크벤치 소프트웨어를 사용하여 분석되었다. 프라이머의 효율성을 확인하기 위해 다른 국가의 여러 균주와 파키스탄의 임상 샘플을 동일한 실험을 수행했다.

한국등록특허 제2030245호는 치쿤군야 바이러스 검출용 올리고뉴클레오티드 세트 및 이의 용도에 관해 개시하고 있고, 한국공개특허 제2011-0118176호는 치쿤군야의 검출을 위한 프로브 및 프라이머에 관해 개시하고 있으며, 미국공개특허 제2019/0367998호는 치쿤군야 바이러스 검출 방법이 개시되어 있다.

하지만, 본 발명의 치쿤군야 바이러스의 전장 유전체 검출을 위한 범용 프라이머 세트에 대해서는 아직까지 개시된 바가 없다.

한국등록특허 제2030245호 한국공개특허 제2011-0118176호 미국공개특허 제2019/0367998호

본 발명은 상기와 같은 요구에 의해 도출된 것으로서, 본 발명자들은 치쿤군야 바이러스에서 증폭이 가능한 범용 프라이머를 활용하여 미량의 유전체로부터 PCR 등을 통하여 증폭된 유전체를 염기서열분석법으로 빠르고 간편하게 바이러스의 전장 유전체 염기서열을 확보할 수 있음을 확인하여 본 발명을 완성하였다.

상기 과제를 해결하기 위하여, 본 발명은 하기의 단계를 포함하는 분리된 시료 내 치쿤군야 바이러스의 전장 유전체 서열 확인 방법을 제공한다.

상기 서열 확인 방법은 (a) 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 제조하는 단계;

(b) 시료에서 RNA를 분리하는 단계;

(c) 분리된 RNA를 사용하여 역전사 반응을 수행하는 단계;

(d) 역전사 반응이 수행된 시료에 프라이머 세트를 첨가하여 중합효소 연쇄반응을 실시하는 단계; 및

(e) 증폭된 산물의 서열을 분석하는 단계를 포함한다.

또한, 상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 것이다.

본 발명은 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 포함하는 치쿤군야 바이러스 진단 키트로서,

상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 치쿤군야 바이러스 진단 키트를 제공한다.

본 발명은 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 포함하는 알파 코로나 바이러스 유발 질환 진단용 조성물로서,

상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 것인 치쿤군야 바이러스 유발 질환 진단용 조성물을 제공한다.

본 발명의 일 예에서, 상기 시료는 혈액, 혈청, 객담, 소변 또는 생체 조직으로부터 분리된 것이다.

본 발명의 다른 예에서, 상기 치쿤군야 바이러스는 KNIH/2009/77을 포함하는 것이다.

본 발명의 또 다른 예에서, 상기 치쿤군야 바이러스 유발 질환은 돌발성 발열, 다발성 관절통 및 발진으로 이루어진 군에서 선택된 하나 이상인 것이다.

본 발명의 특징 및 이점을 요약하면 다음과 같다.

본 발명은 치쿤군야 바이러스 전장 유전체의 서열 획득 방법 및 이의 용도에 관한 것이다.

본 발명의 방법, 키트 및 조성물은 치쿤군야 바이러스의 전장유전체를 증폭시킬 수 있는 정방향 및 역방향 프라이머 세트를 이용하고, 유전체 증폭을 통하기 때문에 미량의 바이러스 시료 및 미배양 임상시료 등에서 치쿤군야 바이러스의 존재 유무 및 변종유무를 간편하고 빠르게 확인할 수 있다.

따라서, 본 발명의 방법, 키트 및 조성물은 치쿤군야 바이러스의 전장유전체 서열을 정확하고 신속하게 얻을 수 있게 해주며, 이를 통해 치쿤군야 바9..이러스로 유발된 질환의 진단에 효과적으로 적용될 수 있다.

도 1은 본 발명의 일 구현 예에 따른 계통발생 분석결과를 나타낸 것이다. (1A) 치쿤군야 바이러스의 계통수는 3개의 계통으로 나뉜다; 서부 아프리카(WA), 동부 중남부 아프리카(ECSA), 아시아. (1B) 검은색 사각형은 KNIH/2009/77 균주와 다른 균주 사이의 계통발생적 거리를 보기 위해 (1A)의 빨간색 화살표를 기반으로 한 작은 나무 확대이다.

이하, 본 발명의 바람직한 구현예에 대하여 상세히 설명한다. 또한, 하기의 설명에서는 구체적인 구성요소 등과 같은 많은 특정 사항들이 도시되어 있는데, 이는 본 발명의 보다 전반적인 이해를 돕기 위해서 제공된 것일 뿐 이러한 특정 사항들 없이도 본 발명이 실시될 수 있음은 이 기술분야에서 통상의 지식을 가진 자에게는 자명하다 할 것이다. 그리고, 본 발명을 설명함에 있어서, 관련된 공지 기능 혹은 구성에 대한 구체적인 설명이 본 발명의 요지를 불필요하게 흐릴 수 있다고 판단되는 경우 그 상세한 설명을 생략한다.

본 발명의 목적을 달성하기 위하여, 본 발명은 하기의 단계를 포함하는 분리된 시료 내 치쿤군야 바이러스의 전장 유전체 서열 확인 방법을 제공한다.

상기 서열 확인 방법은 (a) 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 제조하는 단계;

(b) 시료에서 RNA를 분리하는 단계;

(c) 분리된 RNA를 사용하여 역전사 반응을 수행하는 단계;

(d) 역전사 반응이 수행된 시료에 프라이머 세트를 첨가하여 중합효소 연쇄반응을 실시하는 단계; 및

(e) 증폭된 산물의 서열을 분석하는 단계를 포함한다.

또한, 상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 것이다.

상기 프라이머의 농도는 각각 0.1 내지 1.0 pM, 0.1 내지 0.8 pM, 0.1 내지 0.6 pM, 0.2 내지 1.0 pM, 0.2 내지 0.8 pM, 0.2 내지 0.6 pM, 0.3 내지 1.0 pM, 0.3 내지 0.8 pM, 0.3 내지 0.6 pM, 0.4 내지 1.0 pM 또는 0.4 내지 0.8 pM, 예를 들어, 0.4 내지 0.6 pM일 수 있으나, 이에 한정되는 것은 아니다.

상기 중합효소 연쇄반응을 실시하는 단계는 30 내지 50 사이클, 30 내지 46 사이클, 30 내지 44 사이클, 33 내지 50 사이클, 33 내지 46 사이클, 33 내지 44 사이클, 36 내지 50 사이클, 36 내지 46 사이클, 36 내지 44 사이클, 38 내지 50 사이클 또는 38 내지 46 사이클, 예를 들어, 38 내지 44 사이클 조건으로 수행되는 것일 수 있으나, 이에 한정되는 것은 아니다.

본 발명의 일 예에서, 상기 시료는 혈액, 혈청, 객담, 소변 또는 생체 조직으로부터 분리된 것이다.

본 발명의 다른 예에서, 상기 치쿤군야 바이러스는 KNIH/2009/77을 포함하는 것이다.

상기 증폭된 산물의 서열을 분석하는 단계는 이 분야에 알려진 시퀀싱 방법이면 종류에 관계없이 사용할 수 있고, 바람직하게는 Sanger법 또는 차세대염기서열분석법(NGS; next generation sequencing)을 사용할 수 있으나, 이에 한정되는 것은 아니다.

염색체의 일부가 결핍 또는 중복되어 나타나는 DNA복제수 변이(CNVs)를 포함한 염색체 이상을 확인하기 위해 핵형분석, 형광동소보합법, 염색체 마이크로어레이, NGS기반의 스크리닝 검사와 같이 다양한 검사가 이루어지고 있다 (Capalbo A, et al. 2017, Hum Reprod. Vol. 32(3), pp. 492-498). 핵형분석은 다른 검사들에 비해 5Mb 정도의 낮은 해상도를 보이며 그보다 작은 크기의 염색체 결실/중복은 검출이 불가능하다. 5Mb 미만의 작은 크기의 염색체 결실 및 중복을 미세결실/중복이라고 하며, 단일유전자에 의한 질환 중 미세결실/중복에 의한 비율이 전체 변이의 15%에 해당한다(Vissers LE, et al. 2005, Hum Mol Genet. Vol. 15;14 Spec No. 2:R215-23.).

이러한 미세결실/중복을 검출해내기 위해서 특정 염기서열에 상보적인 탐침자를 활용한 형광동소보합법(FISH)과 염색체 마이크로어레이 검사가 이루어지고 있다. 형광동소보합법은 확인하려는 염기서열에 상보적인 탐침자에 형광라벨을 붙여 염색체 내에 특정 염기서열의 여부를 확인하는 검사법이다. 100kb-1Mb의 해상도를 보이기 때문에 미세결실/중복의 검출이 가능하지만 탐침자 서열에 상보적인 부분만 확인이 가능하기 때문에 기존에 알려진 변이에 대해서만 검출이 가능하다는 단점이 있다.

현재 염색체 미세결실/중복을 확인하는 가장 일반적인 검사법으로 마이크로어레이를 기반으로 하는 비교유전체 혼성화법(aCGH)이 활용되고 있다(Russo CD, et al. 2014, Cancer Discov. Vol. 4(1), pp. 19-21). 마이크로어레이를 통해 검출 가능한 CNV의 크기는 탐침자의 밀도에 의해 결정되며 대략 50kb 크기의 CNV까지 검출이 가능하다. (Watson CT, et al. 2014) 하지만 전좌 또는 역위와 같이 염색체 재배열에 의한 염색체 이상은 검출이 불가능하다.

생명체를 구성하는 유전자의 염기서열 순서를 확정하는 기술은 1세대 생어(Sanger) 방식과 차세대 염기서열 분석법인 NGS(Next Generation Sequencing)로 나뉜다. 1977년경에 개발된 생어 서열분석(Sanger sequencing)이라 불리는 1세대 방법은 유전자 증폭 반응(Polymerase Chain Reaction; PCR) 과정 중 ddNTPs(di-deoxynucleotidetriphosphates)가 DNA 가닥(strand)의 합성을 중단시키는 연쇄-정지(chain-termination) 원리를 응용하여 증폭된 DNA 절편(fragment)을 모아 염기서열을 확인하는 방법이다. 이 방법은 정확도는 높지만, 염기서열 확정에 시간이 오래 걸리며, 높은 비용이 발생한다는 문제점을 가지고 있다.

이와 같은 문제점을 극복하기 위하여 차세대 염기 서열 분석법이 개발되었다. 수많은 유전체의 유전정보를 분절화한 후, 증폭하여 빅데이타를 얻은 다음 생물정보학을 이용하여 원하는 생물체의 전장 유전체 염기서열(Whole genome sequences; WGS)을 확보하는 것이다. 차세대염기서열분석법(NGS)은 염색체를 작은 조각으로 나누고 각 조각의 유전정보를 병렬적으로 분석하는 염기서열분석법이다. NGS는 유전자분석 기술이 발전하면서 상대적으로 검사의 소요시간과 비용이 적고 단일염기 다형성(SNP), 삽입-결실(INDELs)까지 검출 가능한 높은 해상도 때문에 신생아의 유전성 질환 선별검사로 활용되고 있다. 그러나 염색체를 작게 나누어 분석하는 NGS의 원리적 특성상 큰 규모의 염색체의 구조적 변이나 CNVs을 검출하는데 기술적 한계가 있다(Yohe S, Thyagarajan B. 2017, Arch Pathol Lab Med. Vol. 141(11), pp. 1544-1557).

하지만 NGS는 탐침자를 기반으로 하는 마이크로어레이에서 검출할 수 없는 염색체 재배열에 의한 염색체 이상과 기존에 알려지지 않은 새로운 CNV의 검출이 가능하다(Talkowski ME, et al. 2011, Am J Hum Genet. Vol. 88(4), pp. 469-81). 또한 염색체를 작게 조각 내어 염기서열을 분석하는 특성으로 마이크로어레이 보다 더 높은 coverage와 해상도를 보이고 염색체 이상이 시작되는 구획점(breakpoint) 검출이 가능하다는 장점이 있다(Zhao M, et al. 2013, BMC Bioinformatics. Vol. 14, Suppl 11:S1).

본 발명의 목적을 달성하기 위하여, 본 발명은 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 포함하는 치쿤군야 바이러스 진단 키트로서,

상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 치쿤군야 바이러스 진단 키트를 제공한다.

또한, 본 발명은 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 포함하는 알파 코로나 바이러스 유발 질환 진단용 조성물로서,

상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 것인 치쿤군야 바이러스 유발 질환 진단용 조성물을 제공한다.

이때, 본 발명의 또 다른 예에서, 상기 치쿤군야 바이러스 유발 질환은 돌발성 발열, 다발성 관절통 및 발진으로 이루어진 군에서 선택된 하나 이상인 것일 수 있다.

본 발명의 이점 및 특징, 그리고 그것들을 달성하는 방법은 상세하게 후술되어있는 실시예들을 참조하면 명확해질 것이다. 그러나 본 발명은 이하에서 개시되는 실시예들에 한정되는 것이 아니라 서로 다른 다양한 형태로 구현될 것이며, 단지 본 실시예들은 본 발명의 개시가 완전하도록 하고, 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에게 발명의 범주를 완전하게 알려주기 위해 제공되는 것이며, 본 발명은 청구항의 범주에 의해 정의될 뿐이다

실시예 1. RNA standard and cDNA synthesis

균주 NCCP 43132의 바이러스 RNA(KNIH/2009/77)는 대한민국 국립병원체문화원(NCCP)에서 제공했다. 표 1은 이 연구에 사용된 게놈 RNA 및 임상 샘플에 대한 세부 정보를 제공한다. 임상 샘플은 파키스탄에서 cDNA 상태로 전달되었다. 게놈 RNA의 초기 농도는 제조사 매뉴얼에 따라 QuantiFluor RNA System(Promega, USA)과 Quantus™ Fluorometer(Promega, USA)를 이용하여 정량하였다.

Amplicon order	Name	5'-Sequence-3'	Temperature	Length	Site *	Size
1	CHIK9F	GTGAGACACACGTAGCCTAC (서열번호 1)	60.5	20	9-28	3029
	CHIK3018R	CACCTCCCACTCCTTAATAG (서열번호 2)	58.4	20	3018-3037
2	CHIK2900F_1	TCAGAGCACGTCAACGTACT (서열번호 3)	58.4	20	2900-2919	3173
	CHIK6053F_1	ACAGTTGGRTAGTTTCTAGC (서열번호 4)	54.3	20	6053-6072
3	CHIK5813F_1	AAGAAACTCCAGGAGARTGC (서열번호 5)	56.4	20	5813-5832	3113
	CHIK8906F	TGGGTGCGTACATGARTGAC (서열번호 6)	58.4	20	8906-8925
4	CHIK8725F	CAYGATTGGACCAAGCTGCG (서열번호 7)	60.5	20	8725-8744	3055
	CHIK11775F	TCCCTGTGGGTTCGGAGAAT (서열번호 8)	60.5	20	11775-11794

* Genomic position are based on the genomic sequence of strain S27-African prototype (NC_004162)

단일 가닥 cDNA(ss cDNA)는 게놈 RNA에서 합성되었다. 20μl의 반응 혼합물은 4 μl의 5X LunaScript RT supermix(New England Biolabs, USA), 2 μl의 RNA 및 뉴클레아제가 없는 물로 구성된다. 다음 조건에서 반응되었다; 2분 동안 25℃(프라이머 어닐링); 55℃에서 10분 동안(cDNA 합성); 1분 동안 95℃(열 비활성화).

실시예 2. Primer design and Screen

CHIKV의 게놈 서열은 Virus Pathogen Resource(ViPR; www.viprbrc.org) 데이터베이스에서 검색되었다. ViPR을 사용하여 CHIKV의 총 767개 게놈 서열을 정렬하고 CHIKV의 보존 영역을 결정했다. 프라이머는 보존된 서열로 설계되었다. 프라이머는 이전 연구와 유사한 기준으로 설계되었다³⁷. 모든 올리고머는 Macrogen(한국)에서 합성되었다.

PCR 반응은 TaKaRa Ex Taq® 키트(Takara, Japan)를 사용하여 96-Well Thermal Cycler(Veriti Instruments, INC., USA)로 수행되었다. 50μl의 Reaction은 1μl의 각 프라이머(10μM), 5μl의 Ex Taq Buffer(10X), 3μl의 dNTP(각 2.5mM), 2μl의 템플릿, 0.25μl의 TaKaRa Ex Taq Polymerase 및 멸균된 증류수를 포함했다. PCR 조건은 다음과 같다; 95℃에서 1분간(초기 변성); 95℃에서 15초(변성), 60℃에서 15초(어닐링), 72℃에서 3분(연장)의 30 사이클; 3분 동안 72℃(최종 연장). 생성된 앰플리콘은 100 전압에서 40분 동안 전기영동(ADVANCE CO., LTD., Japan)으로 확인하였다. 전체 게놈 시퀀싱을 위해 가장 개별적인 앰플리콘을 생성하는 프라이머 쌍을 선택했다. 선택한 프라이머 시퀀스는 표 1에 나와 있다. 선택한 쌍의 앰플리콘은 제조업체의 프로토콜에 따라 DNA Clean & Concetrator-25(ZYMOresearch, USA)를 사용하여 정제되었다.

실시예 3. Quantification of amplicons

앰플리콘의 양은 실시간 정량적 PCR(qPCR) 및 디지털 액적 PCR(ddPCR) 분석으로 2단계로 측정되었다. 분석에 사용된 프라이머와 프로브는 표 2^38-40에 나열되어 있다. 분석 믹스는 사용 준비를 위해 하나의 튜브에서 각 프라이머와 프로브를 혼합하여 준비했다. 각 프라이머와 프로브는 qPCR의 경우 0.5μM, ddPCR의 경우 0.4μM로 혼합되었다. 각 분석 믹스의 다음 부피는 단일 반응에 사용되어 반응에서 0.2μM의 최종 농도를 만든다. 3μl의 nsp2-1, C, E1 분석 및 5μl의 nsp2-2 분석.

Target amplicon	Target gene	Name	5'-Sequence-3'	Site*	Source
1	nsp2-1	CHIKV_nsp2_F1	CATCTGCACYCAAGTGTACCA	2578-2598	Jesse J.Waggoner et al.
		CHIKV_nsp2_R1	GCGCATTTTGCCTTCGTAATG	2654-2674
		CHIKV_nsp2_P1	GCGGTGTACACTGCCTGTGACYGC	2614-2637
2	nsp2-2	CHIKV_nsp2_F2ª	GAGCATAYGGTTACGCAGATAG	3855-3876	Barbara W. Johnson et al.
		CHIKV_nsp2_R2	TACTGGTGATACATGGTGGTTTC + TGCTGGTGACACATGGTGGTTTC	3934-3956
		CHIKV_nsp2_P2ª	ACGAGTAATCTGCGTACTGGGACGTA + ACGAGTCATCTGCGTAYTGGGACGCA	3886-3911
3	C	CHIKV_C_F	CAACTTGCCCAGCTGATCTC	7684-7703
		CHIKV_C_R	TTCTTATTCTTCCGATTCYTGCG	7750-7772
		CHIKV_C_P	ATAAACTGACAATGCGCGYGGTACC	7712-7736
4	E1	CHIKV_E1_F	AAGCTYCGCGTCCTTTACCAAG	10387-10408	Boris Pastorino et al.
		CHIKV_E1_R	CCAAATTGTCCYGGTCTTCCT	10575-10595
		CHIKV_E1_P	CCAATGTCYTCMGCCTGGACACCTTT	10486-10511

All Taq-man probes were labelled with 5' end FAM as reporter dye and 3' end BHQ1 as quencher dye.

* Genomic position are based on the genomic sequence of S27-African prototype (NC_004162)

ª Oligonucleotides have one base change form original sequences for wide coverage spectrums against CHIKV. Changed base highlight by extrabold.

+ Two Oligonucleotides locate in same site and can cover the different genome spectrum.

qPCR의 반응 혼합물은 Maxima Probe/ROX qPCR Master Mix 2X(Thermo Scientific™, USA)를 사용하여 제조되었다. 각 반응 혼합물의 25μl는 12.5μl 2X Master Mix, 2μl 템플릿, 분석 믹스 및 뉴클레아제가 없는 물로 구성되었다. PCR 조건은 다음과 같다. 95℃에서 10분간 (초기 변성); 95℃ 15초(변성) 및 60℃ 60초(어닐링/신장)의 40 사이클러. 형광은 LightCycler® 96 기기(Roche, Switzerland)의 확장 단계에서 FAM의 동적 범위를 기록하고 LightCycler® 96 소프트웨어(버전 1.1)로 분석했다.

ddPCR의 반응 혼합물은 부피가 20 μl이고 프로브용 ddPCR 슈퍼믹스(BioRad Laboratories, USA) 10 μl, 주형 2μl, assay mix 및 nuclease-free water로 구성되었다. 모든 유효한 사본 번호는 제조업체의 지침에 따라 동적 범위 FAM(470nM/514nM)에서 선택되었다. ddPCR 실험은 QX200 시스템(BioRad Laboratories, Hercules, CA, USA)을 사용하여 수행되었다. 원시 데이터는 처음에 QuantaSoft 소프트웨어(BioRad Laboratories, Hercules, CA, USA)를 사용하여 분석되었다.

앰플리콘의 복제 수는 이전에 설명한 대로 ddPCR 결과로 선형 회귀 분석을 사용하여 qPCR 결과에서 계산되었다⁴¹.

실시예 4. Whole Genome Shotgun Sequencing

1) Illumina Library preparation and sequencing

RNA sequencing에는 2가지 방법이 있다. 하나는 이중 가닥 cDNA(ds cDNA)를 만드는 것이고 다른 하나는 RNA library prep kit를 사용하는 것이다. 본 발명에서는 단일 가닥 cDNA에서 두 번째 가닥 cDNA를 합성하여 ds cDNA를 만들었다. 단일 가닥 cDNA는 random hexamer와 oligo-dt primer를 모두 포함하는 LunaScript RT supermix kit(NEB, England)를 사용하여 합성되었다. 혼합물은 LunaScript RT SuperMix(5X) 4μl, 템플릿 4μl 및 최대 20μl의 뉴클레아제가 없는 물 부피로 구성되었다. 두 번째 가닥 cDNA 합성은 제조사의 프로토콜에 따라 Maxima™ H Minus Double-Stranded cDNA Synthesis kit(Thermo Fisher Scientific, USA)를 사용하여 수행되었다. 그런 다음 100U의 RNase I을 100㎕의 ds cDNA 반응물에 첨가하고 실온에서 5분 동안 인큐베이션하여 잔류 RNA를 제거하였다. DNA Clean & Concentrator Kits(Zymo Research, USA)로 정제된 ds cDNA는 QuantiFluor ds DNA System(Promega, USA) 및 Quantus™ Fluorometer를 사용하여 제조업체 설명서에 따라 정량화했다.

라이브러리 준비는 제조사 매뉴얼에 따라 300-350bp PCR 단편을 만드는 Nextera DNA Flex Library Prep Kit(Illumina, USA)를 사용하여 수행되었다. 실제 어댑터 연결 DNA의 정량화를 위해 프로토콜에 따라 LightCycler® 96 기기(Roche, Switzerland)에서 KAPA 라이브러리 정량화 키트(Roche, Switzerland)를 사용하여 라이브러리 정량화를 확인했다. 라이브러리 샘플 정량화를 위한 정량화 대조군으로 10nM PhiX Control v3(Illumina, USA)의 10배 희석 시리즈를 사용했다. 샘플 농도는 LightCycler® 96 Software(version 1.1)의 Absolute Quantitative Analysis 프로그램에 의해 Cq 값을 기준으로 분석되었다. 100pM의 최종 로딩 샘플은 시퀀싱 실행 중에 각 기본 신호를 고르게 돕는 5% PhiX 컨트롤을 추가하여 준비했다. 20μl의 최종 로딩 샘플을 로딩한 다음 기기의 지침에 따라 시퀀싱 프로그램을 실행하였다.

2) Sequence analysis

CLC Genomic Workbench Software(20.0.4 버전, Qiagen, Denmark)를 판독의 트리밍 및 매핑에 사용했다. Illumina의 두 쌍 끝 판독(정방향-역방향)은 품질이 낮은 판독값을 제거하기 위해 다음 매개변수 조건으로 트리밍되었다. Phred 품질 점수 30점; 최대 모호성 수 중 1개; 호모폴리머(단말에서 20개 염기 중 9개 염기에 걸쳐 동일한 염기를 의미하는 호모폴리머); 5' 및 3' 터미널에서 각각 15 및 10개의 베이스를 제거했다. 길이가 60 미만인 읽기는 삭제한다. 트리밍된 읽기는 NCBI에서 다운로드한 FJ445484(11,790bp)에 매핑되었다. 매핑 조건은 다음과 같다. 경기 점수 1점; 불일치 점수 2점; 3 삽입 비용; 3 삭제 비용; 자동 감지 페어링 거리; 무작위로 매핑한다(정렬 점수 및 선형 간격 비용은 CLC 워크벤치 소프트웨어에서 제공하는 기본값임). CHIKV의 완전한 서열은 계통발생적 거리를 확인하기 위해 MEGA 소프트웨어(버전 7.0.26, MEGA, USA)에 의해 컨센서스 서열로 분석되었다(표 3). 균주 SG650(AF079456), Onyong-nyong 바이러스(ONNV)는 트리의 외부 그룹으로 사용되었다. Muscle 알고리즘을 이용하여 서열을 정렬하고, Neighbor-joining 알고리즘을 이용하여 계통수를 구성하였다. 합의 시퀀스는 NCBI(GenBank: - )에서 사용할 수 있다.

Strain	Geographic origin	Sampling year	GenBank ID
Ross low-psg	Tanzania	1953	HM045811
TH35	Thailand	1958	HM045810
LSFS	DRC *	1960	HM045809
Angola M2022	Angola	1962	HM045823
Gibbs 63-263	India	1963	HM045813
IbH35	Nigeria	1964	HM045786
3412-78	Thailand	1978	HM045808
HB78	CAR +	1978	HM045822
UgAg4155	Uganda	1982	HM045812
Hu/85/NR/001	Philippines	1985	HM045800
RSU1	Indonesia	1985	HM045797
CHIKV/Homosapiens/USA/HIMSTSSA287/1995	USA	1995	KY575574
HD 180760	Senegal	2005	HM045817
MY/06/37348	Malaysia	2006	FN295483
SL-CK1	Sri Lanka	2007	HM045801
SGEHICHT077808	Singapore	2008	FJ445484
MY/08/065	Malaysia	2008	FN295485
SD08Pan	China	2008	GU199351
Bali 2011	Indonesia	2011	MF773561
PNG 2013	Papua New Guinea	2013	MF773569
Caribbean 2014	Caribbean	2014	MF773560
Philippines 2014	Philippines	2014	MF773563
Pakistan-01/2016	Pakistan	2016	MF774613
hk02	Hong Kong	2016	MF499120
Bangladesh 2017	Bangladesh	2017	MF773566
19CHKYN01LY	China	2019	MN402883
KNIH/2009/77	South Korea	2009

* DRC ; Democratic Republic of the Congo, + CAR ; Central African Republic

실시예 5. MinION nanopore Sequencing

1) Library preparation and Sequencing Run

라이브러리는 제조업체의 프로토콜에 따라 Native Barcoding Expansion 1-12(EXP-NBD104, Oxford Nanopore Technologies, England)와 함께 Ligation kit(SQK-LSK109, Oxford Nanopore Technologies, England)의 1D 앰플리콘/cDNA를 사용하여 수행되었다. 라이브러리 준비 후 1D 플로우 셀(FLO-MIN 106, Oxford Nanopore Technologies, England)을 MinION 기기(Oxford Nanopore Technologies, 영국)에 놓고 컴퓨터의 MinKNOW 소프트웨어에 연결했다. Flow cell Check 및 Loading Sample도 프로토콜에 따라 수행되었다. 시퀀싱은 base-calling으로 24시간 동안 실행되었고 결과적으로 FAST5 파일을 만들었다.

2) Sequence Analysis

분석을 위해 FAST5 파일을 FASTAQ 형식으로 전송하였다. 그런 다음 MinION GUI 소프트웨어에 의해 낮은 품질의 읽기가 트리밍되었다. 트리밍된 판독 데이터는 분석을 위해 CLC Genomic Workbench Software로 가져왔다. 두 번째 트리밍은 다음과 같은 조건으로 읽혀지도록 처리되었다. 최대 모호성 수 중 1개; 호모폴리머 제거; 5' 및 3' 터미널에서 각각 15 및 10개의 베이스를 제거했다. 길이가 60 미만인 판독은 폐기하였다. 트리밍된 판독은 샷건 시퀀싱에 의해 전체 게놈 서열에서 가져온 참조 서열에 매핑되었다. CHIKV의 완전한 서열은 MEGA 소프트웨어에 의해 컨센서스 서열로부터 분석되었다.

실험예 1. Whole Genome Shotgun Sequencing

한국의 모든 치쿤군야 바이러스 사례는 여행 관련 환자를 확인한다. 균주 NCCP 43132(KNIH/2009/77)는 여행 관련 환자로부터 처음 분리되어 국립 병원균 수집(NCCP)에 기탁되었다⁴². 이 균주의 부분 코딩 서열(CDS) 또는 구조 단백질 CDS 영역만이 이전에 연구되었습니다^42-44. 완전한 게놈 기반 추가 연구를 위한 기반을 제공하기 위해 먼저 이중 가닥 cDNA를 사용하여 이 균주의 전체 게놈 샷건 시퀀싱을 분석했다. pair-end 읽기 데이터는 Illumina 시퀀싱에 의해 만들어졌고 CLC 워크벤치에 입력되었다. 읽은 데이터는 품질 관리 보고서를 기반으로 트리밍되었다. 트리밍된 읽기는 SGEHICHD122508(FJ445502)에 매핑되어 이 균주를 분리한 첫번째 연구에서 가장 유사한 읽기를 보여주었다. 가장 상동성인 게놈을 찾기 위해 BLAST 프로그램(NCBI)을 사용하여 컨센서스 서열을 확인하였다. BLAST 결과는 서열이 ECSA 유전자형-인도양 계보(IOL)45의 균주 SGEHICHT077808(FJ445484)과 가장 높은 유사성(쿼리 커버의 100%, 동일성 백분율의 99.94%)을 갖는 것으로 나타났다. SGEHICHT077808(FJ445484)에 다시 매핑하여 11,790bp 길이 크기의 결과 컨센서스 시퀀스를 얻었다. 표 4 및 표 5는 분석 소프트웨어에서 제공하는 트리밍 및 매핑 보고서를 보여준다. 추가 검사를 통해 시퀀스의 끝인 11763-11790 사이트에서 발견된 100 임계값 미만의 범위를 읽는다. 합의 서열은 MEGA7의 근육 정렬에 의해 여러 계통의 완전한 서열과 정렬되었다. Consensus sequence에서 7개의 차이점이 확인되었고 Sanger sequencing을 통해 올바른 것으로 확인되었습니다. 계통 발생 분석은 IOL에 클러스터링된 균주 NCCP 43132의 계통 발생 거리를 확인했습니다(도 1). 균주 NCCP 43132는 IOL^45,46에 클러스터링된 싱가포르, 말레이시아, 파푸아뉴기니의 균주와 유사하다.

(3A)
Trimmed Read		283,946
Not Mapped Read		146,548
Mapped Read	Read Count	137,398
	Number of bases	14,453,299
	Average Length	105.19
	GC contents%	50.08
Mapping Coverage	Minimum	2
	Median	1,225
	Average	1,211.46
	Maximum	1,845
Consensus sequence		11,790

(3A) Chart indicate the result of trimming, mapping of shotgun sequencing based on the number of reads.

(3B)
site	FJ445484	Consensus seq.	note
588	G	A	G/A
2875	T	C	C/T
3165	A	G	A/G
3366	C	T	C
8541	C	T	C
9312	C	T	C/T
9494	T	C	T/G/C

(3B) Seven bases and its site presented difference between reference and consensus sequences. Genomic site are provided based on the reference sequence CHKV strain S27-African prototype (NC_004162). Note show the bases at same site in all whole genome database of CHIKV.

실험예 2. MinION nanopore Sequencing

Chikungunya 바이러스에는 3개 또는 4개의 주요 혈통이 있다. 서아프리카(WA), 동부 중앙 남부 아프리카(ECSA), 아시아 및 인도양 계보(IOL). ECSA, Asian, IOL 혈통이 유행병에 크게 기여했지만 2004년 이전에는 WA 혈통이 Afica 대륙에 집중되었다. 다른 혈통과 확연히 다르다. 본 발명에서 프라이머는 WA 계통 이외의 다른 사람들을 대상으로 설계되었다. 디자인된 프라이머로 얻을 앰플리콘의 크기와 오버랩 크기는 각각 약 3Kb와 100bp를 기준으로 먼저 선정하였고, 최종적으로 젤 전기영동을 통해 비교적 깨끗한 단일 밴드를 나타내는 프라이머를 선정하였다. Amplicon 1의 경우 30여 개의 프라이머 조합에 걸쳐 강한 non-target 크기의 밴드가 연속적으로 형성되었다. 그래서 추가 스크리닝을 중단하고 프라이머 세트를 선택했다.

실험예 2. MinION nanopore Sequencing

Nanopore 시퀀싱 후 OxfordNamopore의 Guppy GPU basecalling 소프트웨어를 사용하여 15Gb 원시 데이터(fast5 형식)를 계산하고 fastq 파일로 변환했다. 변환하는 동안 base-calling, 어댑터 트리밍 및 Q-Score 7 이상으로 필터링하여 시퀀스 품질을 수정했다. CLC Workbench 소프트웨어를 이용하여 생성된 fastq 파일에 대해 추가적인 트리밍을 수행한 후 Illumina shot gun sequencing을 통해 얻은 전체 genome sequence(11,790bp)를 reference로 사용하였다. 획득한 읽기 커버리지 길이 그래프는 1Kb 간격으로 인위적으로 조정되었다. 그 결과 118,012개의 매핑된 읽기, 3,027.56bp의 평균 읽기 길이, 30,629.57bp 길이의 평균 커버리지 길이를 얻었다. Consensus sequence를 추출하여 11,784bp 길이의 전체 genome sequence를 얻었으며, MEGA7을 이용한 reference와 비교하여 돌연변이도 확인하였다.

실험예 3. 고찰

치쿤군야 바이러스는 2004년 이전에 3개의 큰 계보를 가지고 있다. 서아프리카(WA), 동부 중앙 남부 아프리카(ECSA) 및 아시아(아시아 도시라고도 함, AUL). 인도양 섬에서 발견된 바이러스는 ECSA47의 하위 그룹으로 인도양 섬(IOL) 하위 계통으로 분류되었다. 2013-2014년 카리브해와 남미의 치쿤군야 바이러스는 아시아 계통으로 분류되었으며 ECSA 계통의 새로운 미국 하위 그룹이 2016년에 제시되었다⁴⁸. 오늘날 3개 또는 4개의 계통이 CHIKV의 주요 계통으로 사용되었다; WA, ECSA, 아시아 및 IOL⁴⁶.

<참고문헌>

1. Fenner, F., Pereira, H. G., Porterfield, J. S., Joklik, W. K. & Downie, A. W. Family and Generic Names for Viruses Approved by the International Committee on Taxonomy of Viruses, June 1974. INT 3, 193-198 (1974).

2. Galan-Huerta, K. A., Rivas-Estilla, A. M., Fernndez-Salas, I., Farfan-Ale, J. A. & Ramos-Jim

nez, J. Chikungunya virus: A general overview. Medicina Universitaria 17, 175-183 (2015).

3. Rupp, J. C., Sokoloski, K. J., Gebhart, N. N. & Hardy, R. W. Alphavirus RNA synthesis and non-structural protein functions. J Gen Virol 96, 2483-2500 (2015).

4. Robinson, M. C. An epidemic of virus disease in Southern Province, Tanganyika territory, in 1952-1953. I. Clinical Features. Transactions of The Royal Society of Tropical Medicine and Hygiene 49, 28-32 (1955).

5. Ross, R. W. The Newala epidemic: III. The virus: isolation, pathogenic properties and relationship to the epidemic. J. Hyg. 54, 177-191 (1956).

6. Leta, S. et al. Global risk mapping for major diseases transmitted by Aedes aegypti and Aedes albopictus. International Journal of Infectious Diseases 67, 25-35 (2018).

7. Iwamura, T., Guzman-Holst, A. & Murray, K. A. Accelerating invasion potential of disease vector Aedes aegypti under climate change. Nature Communications 11, 2130 (2020).

8. Powers, A. M. & Logue, C. H. Y. 2007. Changing patterns of chikungunya virus: re-emergence of a zoonotic arbovirus. Journal of General Virology 88, 2363-2377.

9. Coffey, L. L., Failloux, A.-B. & Weaver, S. C. Chikungunya Virus-Vector Interactions. Viruses 6, 4628-4663 (2014).

10. Zeller, H., Van Bortel, W. & Sudre, B. Chikungunya: Its History in Africa and Asia and Its Spread to New Regions in 2013-2014. The Journal of Infectious Diseases 214, S436-S440 (2016).

11. Rezza, G. Chikungunya is back in Italy: 2007-2017. Journal of Travel Medicine 25, (2018).

12. Schuffenecker, I. et al. Genome Microevolution of Chikungunya Viruses Causing the Indian Ocean Outbreak. PLoS Med 3, (2006).

13. Kariuki Njenga, M. et al. Tracking epidemic Chikungunya virus into the Indian Ocean from East Africa. Journal of General Virology 89, 2754-2760 (2008).

14. Khan, K. et al. Assessing the Origin of and Potential for International Spread of Chikungunya Virus from the Caribbean. PLoS Curr 6, (2014).

15. Fischer, M. & Staples, J. E. Chikungunya Virus Spreads in the Americas ― Caribbean and South America, 2013-2014. MMWR Morb Mortal Wkly Rep 63, 500-501 (2014).

16. Johansson, M. A. Chikungunya on the move. Trends in Parasitology 31, 43-45 (2015).

17. Fritz, M. et al. Re-emergence of chikungunya in the Republic of the Congo in 2019 associated with a possible vector-host switch. International Journal of Infectious Diseases 84, 99-101 (2019).

18. Dudouet, P. et al. Chikungunya resurgence in the Maldives and risk for importation via tourists to Europe in 2019-2020: A GeoSentinel case series. Travel Medicine and Infectious Disease 36, 101814 (2020).

19. Appassakij, H., Silpapojakul, K., Promwong, C. & Rujirojindakul, P. The Potential Impact of Chikungunya Virus Outbreaks on Blood Transfusion. Transfusion Medicine Reviews 34, 23-28 (2020).

20. Bohers, C. et al. The recently introduced Aedes albopictus in Tunisia has the potential to transmit chikungunya, dengue and Zika viruses. PLOS Neglected Tropical Diseases 14, e0008475 (2020).

21. Kakarla, S. G. et al. Temperature dependent transmission potential model for chikungunya in India. Science of The Total Environment 647, 66-74 (2019).

22. Houldcroft, C. J., Beale, M. A. & Breuer, J. Clinical and biological insights from viral genome sequencing. Nature Reviews Microbiology 15, 183-192 (2017).

23. Fleischmann, R. D. et al. Whole-genome random sequencing and assembly of Haemophilus influenzae Rd. Science 269, 496-512 (1995).

24. Motahari, A. S., Bresler, G. & Tse, D. N. C. Information Theory of DNA Shotgun Sequencing. IEEE Trans. Inform. Theory 59, 6273-6289 (2013).

25. Phan, H. T. T. et al. Illumina short-read and MinION long-read WGS to characterize the molecular epidemiology of an NDM-1 Serratia marcescens outbreak in Romania. Journal of Antimicrobial Chemotherapy 73, 672-679 (2018).

26. Shendure, J. & Ji, H. Next-generation DNA sequencing. Nature Biotechnology 26, 1135-1145 (2008).

27. Voelkerding, K. V., Dames, S. A. & Durtschi, J. D. Next-Generation Sequencing: From Basic Research to Diagnostics. Clinical Chemistry 55, 641-658 (2009).

28. Kumar, K., Cowley, M. & Davis, R. Next-Generation Sequencing and Emerging Technologies. Seminars in Thrombosis and Hemostasis 45, (2019).

29. Zhou, X. et al. The next-generation sequencing technology and application. Protein Cell 1, 520-536 (2010).

30. Bayley, H. Nanopore Sequencing: From Imagination to Reality. Clin Chem 61, 25-31 (2015).

31. Stapleford, K. A. et al. Whole-Genome Sequencing Analysis from the Chikungunya Virus Caribbean Outbreak Reveals Novel Evolutionary Genomic Elements. PLOS Neglected Tropical Diseases 10, e0004402 (2016).

32. Wilson, J. A. C. et al. RNA-Seq analysis of chikungunya virus infection and identification of granzyme A as a major promoter of arthritic inflammation. PLOS Pathogens 13, e1006155 (2017).

33. Kafetzopoulou, L. E. et al. Assessment of metagenomic Nanopore and Illumina sequencing for recovering whole genome sequences of chikungunya and dengue viruses directly from clinical samples. Eurosurveillance 23, 1800228 (2018).

34. Suwannakarn, K., Theamboonlers, A. & Poovorawan, Y. Molecular genome tracking of East, Central and South African genotype of Chikungunya virus in South-east Asia between 2006 and 2009. Asian Pacific Journal of Tropical Medicine 4, 535-540 (2011).

35. Chen, R. et al. Comprehensive Genome Scale Phylogenetic Study Provides New Insights on the Global Expansion of Chikungunya Virus. J Virol 90, 10600-10611 (2016).

36. Volk, S. M. et al. Genome-Scale Phylogenetic Analyses of Chikungunya Virus Reveal Independent Emergences of Recent Epidemics and Various Evolutionary Rates. Journal of Virology 84, 6497-6504 (2010).

37. Kim, K. W. et al. A Half-Day Genome Sequencing Protocol for Middle East Respiratory Syndrome Coronavirus. Front. Microbiol. 12, (2021).

38. Jesse J. Waggoner et al. Clinical evaluation of a single-reaction real-time RT-PCR for pan-dengue and chikungunya virus detection. Journal of Clinical Virology 78, 57-61 (2016).

39. Johnson, B. W., Russell, B. J. & Goodman, C. H. Laboratory Diagnosis of Chikungunya Virus Infections and Commercial Sources for Diagnostic Assays. The Journal of Infectious Diseases 214, S471-S474 (2016).

40. Pastorino, B. et al. Development of a TaqMan® RT-PCR assay without RNA extraction step for the detection and quantification of African Chikungunya viruses. Journal of Virological Methods 124, 65-71 (2005).

41. Park, C. et al. Comparison of Digital PCR and Quantitative PCR with Various SARS-CoV-2 Primer-Probe Sets. 31, 358-367 (2021).

42. Cha, G. W. et al. Travel-Associated Chikungunya Cases in South Korea during 2009-2010. Osong Public Health Res Perspect 4, 170-175 (2013).

43. Hwang, J.-H. & Lee, C.-S. The First Imported Case Infected with Chikungunya Virus in Korea. Infection & Chemotherapy 47, 55-59 (2015).

44. Chung, W.-C., Hwang, K. Y., Kang, S.-J., Kim, J.-O. & Song, M. J. Development of a neutralization assay based on the pseudotyped chikungunya virus of a Korean isolate. J Microbiol. 58, 46-53 (2020).

45. de Bernardi Schneider, A. et al. Updated Phylogeny of Chikungunya Virus Reveals Lineage-Specific RNA Architecture. https://www.preprints.org/manuscript/201907.0161/v1 (2019) doi:10.20944/preprints201907.0161.v1.

46. Xavier, J. et al. Chikungunya virus ECSA lineage reintroduction in the northeasternmost region of Brazil. 11.

47. Parola, P. et al. Novel Chikungunya Virus Variant in Travelers Returning from Indian Ocean Islands. Emerg Infect Dis 12, 1493-1499 (2006).

48. White, S. K. et al. A new 'American' subgroup of African-lineage Chikungunya virus detected in and isolated from mosquitoes collected in Haiti, 2016. PLOS ONE 13, e0196857 (2018).

<110> KRISS <120> Chikungunya virus universal primer sets for whole genome amplification method and diagnosis kit <130> M21-0000 <160> 8 <170> KoPatentIn 3.0 <210> 1 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK9F <400> 1 gtgagacaca cgtagcctac 20 <210> 2 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK3018R <400> 2 cacctcccac tccttaatag 20 <210> 3 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK2900F_1 <400> 3 tcagagcacg tcaacgtact 20 <210> 4 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK6053F_1 <400> 4 acagttggrt agtttctagc 20 <210> 5 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK5813F_1 <400> 5 aagaaactcc aggagartgc 20 <210> 6 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK8906F <400> 6 tgggtgcgta catgartgac 20 <210> 7 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK8725F <400> 7 caygattgga ccaagctgcg 20 <210> 8 <211> 20 <212> DNA <213> Artificial Sequence <220> <223> CHIK11775F <400> 8 tccctgtggg ttcggagaat 20

Claims (10)

1. 하기의 단계를 포함하는 분리된 시료 내 치쿤군야 바이러스의 전장 유전체 서열 확인 방법
   (a) 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 제조하는 단계;
   (b) 시료에서 RNA를 분리하는 단계;
   (c) 분리된 RNA를 사용하여 역전사 반응을 수행하는 단계;
   (d) 역전사 반응이 수행된 시료에 프라이머 세트를 첨가하여 중합효소 연쇄반응을 실시하는 단계; 및
   (e) 증폭된 산물의 서열을 분석하는 단계
   상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 것임.
2. 제1항에 있어서, 상기 시료는 혈액, 혈청, 객담, 소변 또는 생체 조직으로부터 분리된 것인, 방법.
3. 제1항에 있어서, 상기 치쿤군야 바이러스는 KNIH/2009/77을 포함하는 것인, 방법.
4. 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 포함하는 치쿤군야 바이러스 진단 키트로서,
   상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 치쿤군야 바이러스 진단 키트.
5. 제4항에 있어서, 진단 시료는 혈액, 혈청, 객담, 소변 또는 생체 조직으로부터 분리된 것인, 키트.
6. 제4항에 있어서, 상기 치쿤군야 바이러스는 KNIH/2009/77을 포함하는 것인, 키트.
7. 치쿤군야 바이러스 전장 유전체에 상보적인 하나 이상의 정방향 및 역방향 프라이머 세트를 포함하는 치쿤군야 바이러스 유발 질환 진단용 조성물로서,
   상기 프라이머 세트는 서열번호 1 및 서열번호 2; 서열번호 3 및 서열번호 4; 서열번호 5 및 서열번호 6; 및 서열번호 7 및 서열번호 8의 프라이머 세트를 포함하는 것인 치쿤군야 바이러스 유발 질환 진단용 조성물.
8. 제7항에 있어서, 상기 치쿤군야 바이러스 유발 질환은 돌발성 발열, 다발성 관절통 및 발진으로 이루어진 군에서 선택된 하나 이상인 것인, 조성물.
9. 제7항에 있어서, 진단 시료는 혈액, 혈청, 객담, 소변 또는 생체 조직으로부터 분리된 것인, 조성물.
10. 제7항에 있어서, 상기 치쿤군야 바이러스는 KNIH/2009/77을 포함하는 것인, 조성물.