변이단위 및 유전자단위 분석을 통한 약물유전체 분석 방법론 연구

Abstract

서론: 차세대 시퀀싱 기술의 발전으로 인간 유전체 분석은 임상환경에서의 예후, 예측 및 진단에 도움을 주고 있다. 지금까지 단일 변이의 질병 연관성 분석은 특정 질병군에서 큰 성공을 거두었지만, 드문 변이에 대한 분석의 한계와 높은 유의성을 보인 변이에 대한 해석의 어려움이 수반되는 문제점이 있기 때문에 유전자 단위 혹은 변이 집단적인 분석 방법의 필요성이 대두되었다. 또한 기술적 한계로 인한 시퀀싱 정보의 품질 문제는 임상환경에서 요구하는 품질 수준에 미치지 못한다. 이에 시퀀싱 결과에 대한 재검증을 요구하고 있으며 유의한 결과에 대한 재현이 어렵다는 것도 지적하고 있다. 따라서 임상환경에서 요구하는 수준의 시퀀싱 정보를 가공하고 판별하는 방법뿐만 아니라 단일 변이 분석의 한계를 극복하는 유전자 단위의 복합적인 분석 방법을 활용하는 것이 필수적이다.

방법: 거짓 양성 변이 판별을 위하여 시퀀싱 된 리드를 효과적으로 파악할 수 있는 도구를 개발하여 675개의 매우 드문 변이에 대한 리드 정렬상태를 확인하였다. 거짓 양성으로 판별된 변이로부터 정렬 상태의 특징을 찾아 계통오류 프로파일을 정의하였다. 또한 유전자 점수 기법의 활용성을 파악하기 위해서 개인화된 유전자 점수 (Personalized Gene Score

PGS), Residual Variation Intolerance Score (RVIS) 그리고 Gene Damage Index (GDI) 를 질병과의 연관성이 보고된 유전자에 대하여 각 점수별 우선순위에 대한 민감도와 특이도를 측정하여 성능을 비교하였다.

결과: 매우 드문 변이 675개 중에서, 리드 정렬 상태를 파악하여 519개의 거짓 양성 변이를 판별하였다. 519개의 거짓 양성 변이를 야기한 리드 정렬 상태를 (1) Simplicity region, (2) SNV cluster, (3) Periperal sequence read 그리고 (4) Base inversion 계통오류로 분류 및 정의하였다. 최신의 유전자 점수의 성능을 다양한 질병의 범주에서 평가한 결과, 세 가지 유전자 점수 중, PGS는 다른 유전자 점수에 비해서 Pharmacogenomics과 관계된 유전자의 우선 순위에서 최고의 성능을 나타내었다 (AUCPGS=0.68, AUCRVIS=0.43, AUCGDI=0.62). 약물 관련 유전자 리스트 중에서 ADME core 유전자에 대해서도 역시 PGS가 압도적인 성능을 보여주었다 (AUCPGS=0.72, AUCRVIS=0.52, AUCGDI=0.65). 반면에 희귀질환 관련 유전자의 경우, OMIM Haploinsufficiency 유전자에 대해서는 RVIS가 높은 성능을 보여주었다 (AUCPGS=0.46, AUCRVIS=0.72, AUCGDI=0.47).

결론: 시퀀싱 결과를 재사용하는 방법으로 거짓 양성 변이를 최소화할 수 있도록 계통오류 네 가지를 정의하였고 이를 확인할 수 있는 최적의 도구를 개발하였다. 임상환경에서 요구하고 있는 고품질의 시퀀싱 결과를 위한 초석을 마련하였으며 계통오류를 감안한 변이 추출 알고리즘의 발전 또한 기대해볼 수 있다. 또한 세 가지 유전자 점수를 다양한 질병의 관점에서 비교한 결과, PGS는 약물 유발 표현형과 연관된 유전자, RVIS는 멘델의 법칙을 따르는 희귀질환 유전자에 대한 우선순위를 최적으로 제시해주고 있다. 따라서 분석의 목적에 맞는 유전자 점수를 선택하여 적용하는 것이 합리적이다.