칸디다 알비칸스 미토콘드리아 비교 유전체 분석 중립 진화 비코딩 영역 발견

초록

칸디다 알비칸스의 미토콘드리아 유전체를 3개 균주(참조주 SC5314와 임상주 L296, L757)에서 비교 분석하였다. 전체 33 928 bp 중 372개의 다형성을 확인했으며, 특히 tRNAGly‑COX1, NAD3‑COB, ssurRNA‑NAD4L 사이의 세 개 인터제닉 영역(IG1, IG2, IG3)이 중립적 변이를 많이 보였다. 라틴아메리카 18개 임상주와 ATCC 표준주에서 IG1‑IG3을 PCR 증폭·시퀀싱한 결과, 서열 및 길이 변이가 크게 나타났고, 아르헨티나 균주에만 특이적인 삽입이 관찰되었다. 중립 진화 특성과 높은 변이성을 가진 이 영역은 기존 다중위치 서열 타이핑(MLST)을 보완하는 새로운 분자 타이핑 마커로 활용될 수 있다.

상세 요약

본 연구는 칸디다 알비칸스(Candida albicans)의 미토콘드리아 DNA(MtDNA)가 핵 DNA에 비해 높은 변이율과 빠른 진화 속도를 가지고 있다는 전제 하에, 미토콘드리아 전반의 인트라스페시픽 변이를 체계적으로 탐색하고자 했다. 8배 정도의 시퀀싱 커버리지를 확보한 뒤, 임상주 L296과 L757의 미토콘드리아 전장을 조립하고, 공개된 참조주 SC5314와 정밀 정렬하였다. 전체 33 928 bp 중 372개의 다형점(Polymorphic sites)을 발견했으며, 이 중 230개는 코딩 영역, 142개는 비코딩 영역에 위치한다. 특히 코딩 영역 내 변이는 대부분 비동의성(Non‑synonymous) 변이와 동의성(Synonymous) 변이의 비율이 낮아, 기능적 제약이 존재함을 시사한다. 반면, 비코딩 영역, 특히 tRNAGly‑COX1, NAD3‑COB, ssurRNA‑NAD4L 사이에 위치한 IG1, IG2, IG3은 변이 밀도가 현저히 높고, 대부분 삽입·삭제(indel)와 단순 반복 서열 변이로 구성되어 중립적 진화(neutral evolution) 특성을 보였다.

이 세 인터제닉 영역을 표적으로 PCR을 수행하고, 라틴아메리카 전역에서 수집한 18개 임상주와 ATCC 표준주 2개를 추가 시퀀싱함으로써, 지역별 유전적 구조를 파악하였다. IG1, IG2, IG3은 각각 56 bp까지 길이 차이를 보였으며, 아르헨티나 균주에만 49 bp 규모의 삽입이 독점적으로 존재한다는 점이 눈에 띈다. 이러한 삽입은 특정 지리적 군집(cluster)을 형성하는 데 기여할 가능성이 있다. 계통수 분석에서는 IG1‑IG3을 각각 독립적인 마커로 사용했을 때, 세 개의 주요 클러스터가 도출되었으며, 이는 기존의 다중위치 서열 타이핑(MLST) 결과와 부분적으로 일치하지만, 미토콘드리아 특유의 고변이성을 이용해 보다 미세한 분류가 가능함을 보여준다.

연구진은 IG1‑IG3이 PCR 증폭이 용이하고, 전체 길이 시퀀싱만으로도 충분히 변이를 포착할 수 있다는 실용적 장점을 강조한다. 또한, 중립 진화 영역이므로 선택 압력에 의한 왜곡이 적어, 역학 조사와 감염원 추적에 신뢰할 수 있는 마커가 될 수 있다. 향후 대규모 샘플에 적용한다면, 지역별 전파 경로와 클론 다양성을 정밀하게 파악할 수 있을 것으로 기대된다.