🧬 biology

RNA folding below the nanopore impairs direct RNA sequencing fidelity and modification calling

본 연구는 직접 RNA 시퀀싱 중 나노포어 하단에서 접히는 RNA 이차 구조가 체계적인 이온 전류 왜곡과 위양성 변형 호출을 유발한다는 점을 식별하였으며, 이를 통해 시퀀싱 충실도를 향상시키기 위한 구조 인지 모델의 결정적인 필요성을 확립한다.

원저자: Jana Wiedermannová, Marek Schwarz, Josef Pánek, Saran Natarajan, Veronika Kočárková

게시일 2026-07-13

📖 3 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Jana Wiedermannová, Marek Schwarz, Josef Pánek, Saran Natarajan, Veronika Kočárková

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ⚕️ 이것은 동료 심사를 거치지 않은 프리프린트의 AI 생성 설명입니다. 의학적 조언이 아닙니다. 이 내용을 바탕으로 건강 관련 결정을 내리지 마세요. 전체 면책 조항 읽기

당신은 나노포어(nanopore)라고 불리는 작고 첨단 기술이 집약된 터널을 통과해 빠르게 지나가는 긴 리본 모양의 RNA에 적힌 비밀 메시지를 읽으려고 노력 중이라고 상상해 보세요. 이것이 Oxford Nanopore Direct RNA Sequencing (DRS)가 작동하는 방식입니다. 이 방식은 리본이 터널을 통 melalui 지나갈 때 발생하는 전기적 "웅성거림(hum)"을 듣고 글자(A, U, C, G)를 파악하며, 심지어 그 글자에 붙어 있는 작은 화학적 스티커(변형/modifications)까지 찾아냅니다.

하지만 반전이 있습니다. 이 논문은 RNA 리본이 단순히 직선 형태가 아니라는 점을 밝혀냈습니다. 리본이 터널의 반대편(trans compartment)으로 빠져나오자마자, 마치 머리핀처럼 스스로 말려 올라가 매듭을 지을 수 있다는 것입니다. 연구진은 이러한 매듭이 기계를 혼란스럽게 만들어 많은 가짜 알람(false alarms)을 유발한다는 사실을 발견했습니다.

"헤어핀(Hairpin)" 문제
나노포어를 지하철 개찰구라고 생각해 보세요. RNA 리본은 승객입니다. 보통 승객은 부드럽게 통과합니다. 하지만 승객이 개찰구를 통과하자마자 갑자기 신발 끈을 꽉 조여 매기 위해 멈춰 서서 매듭(헤어핀 구조)을 만든다면 어떻게 될까요? 이 동작은 리본을 잡아당깁니다. 이 당김은 승객의 움직임과 개찰구 센서가 그들을 읽는 방식을 변화시킵니다.

저자들은 RNA 헤어핀이 포어(pore) 바로 아래에서 형성될 때, 매듭이 생기기 약 15글자 전부터 "혼란의 구역(zone of confusion)"이 발생한다는 것을 발견했습니다. 이 구역에서 기계는 다음과 같은 실수를 저지르기 시작합니다:

실제로는 존재하지 않는 곳에 화학적 스티커(변형)가 있다고 착각합니다.
글자를 잘못 읽습니다(서열 오류).
글자의 품질 점수(Quality Score)를 낮춥니다.

"유령" 변형(Ghost Modifications)
논문은 이것이 단순한 우연이 아님을 보여줍니다. 연구진은 실제 화학적 스티커가 전혀 없는 것이 보장된 "in vitro"(실험실 제작) RNA로 테스트를 진행했습니다. 그럼에도 불구하고 기계는 헤어핀 바로 앞에서 스티커가 발견되었다고 보고했습니다. 이는 마치 금속 탐지기가 보물이 묻혀 있지 않은 금속 기둥 근처에서 크게 울리는 것과 같습니다. 그 "보물"(변형)은 매듭이 리본을 잡아당겨 만들어낸 유령이었습니다.

배제된 요인들
저자들은 이 문제가 RNA가 가닥의 끝부분 근처에 있거나 단순히 "지저도한" 서열 때문에 발생하는 것이 아님을 증명하기 위해 매우 신중하게 검증했습니다.

RNA 길이를 짧게 만들거나 길게 만들었음에도 불구하고, 문제는 리본의 끝이 아니라 매듭이 있는 위치에 정확히 머물러 있었습니다.
RNA가 스스로 묶이지 못하도록 몇 개의 글자를 변경하여 매듭을 끊어냈습니다. 그러자 매듭이 사라졌고, 유령 변형과 오류도 함께 사라졌습니다.
그들은 다른 소프트웨어 모델(Dorado 버전 5.0.0, 5.2.0, 5.3.0, 6.0.0)들도 테스트했습니다. 최신 모델들이 약간의 도움은 주었지만, 문제를 완전히 해결하지는 못했습니다. 유령 신호는 여전히 존재하며, 이는 문제가 단순히 소프트웨어 버그가 아니라 물리적인 현상(매듭이 리본을 잡아당기는 것)임을 시사합니다.

증거
연구팀은 단순히 추측한 것이 아니라 측정했습니다.

전류를 관찰한 결과, 헤어핀 직전에 전류가 이상하게 흔들리는 것을 발견했습니다.
"체류 시간(dwell time)"(RNA가 포어 안에 머무는 시간)을 측정했을 때도 변화가 있음을 확인했습니다.
박테리아(Bacillus subtilis), 인간 세포, E. coli에서 테스트했으나, 동일한 "매듭에 의한 혼란"이 모든 곳에서 발생했습니다.
실험실 제작 RNA의 경우, 낮은 품질의 지점 근처에서 나타나는 가짜 변형이 우연히 발생할 확률보다 약 7.5배 더 많이 나타난다는 것을 발견했습니다.

왜 중요한가
이 논문은 만약 과학자들이 질병을 이해하기 위해 RNA 상의 실제 화학적 스티커를 찾고 있다면, 매우 주의해야 한다고 제안합니다. 만약 강한 RNA 매듭 바로 앞에서 변형이 발견된다면, 그것은 매듭이 리본을 잡아당겨서 생긴 "유령"일 수 있기 때문입니다.

저자들은 이를 해결하기 위한 세 가지 방법을 제시합니다:

매듭 표시하기: 컴퓨터가 예측된 매듭의 15글자 전 범위 내에서 발견된 "변형"에 대해 더욱 회의적으로 판단하도록 설정합니다.
AI를 더 잘 훈련시키기: 소프트웨어 모델이 이러한 매듭을 인식하여 속지 않도록 가르칩니다.
매듭 방지하기: 용액에 화학 물질을 추가하여 RNA가 직선을 유지하게 함으로써, 포어 아래에서 헤어핀 형태로 묶이는 것을 막습니다.

요약하자면, 이 논문은 나노포어 아래에서의 RNA 폴딩(folding)이 기계를 속이는 숨겨진 주범임을 보여줍니다. 이것은 센서를 속이는 물리적인 줄다리기이며, 우리가 이러한 매듭을 고려하지 않는 한, 우리의 RNA 변형 지도는 가짜 랜드마크를 포함하게 될 것입니다.

기술 요약: 나노포어 하단의 RNA 폴딩이 Direct RNA 시퀀싱의 충실도 및 변형 호출을 저해함

문제 정의
Oxford Nanopore Technologies(ONT)의 Direct RNA Sequencing(DRS)은 전체 길이의 천연 RNA 분자를 시퀀싱하고 이온 전류 신호로부터 화학적 변형을 직접 추론할 수 있는 독보적인 능력을 제공한다. 그러나 현재 DRS는 short-read cDNA 시퀀싱에 비해 per-read 정확도가 낮으며, 인델(indel)과 위양성(false-positive) 변형 호출이 빈번하게 발생하는 문제가 있다. 기존 연구들은 이러한 오류를 낮은 신호 대 잡음비(SNR), 뉴클레오타이드 간의 화학적 유사성, 그리고 호모폴리머(homopolymeric) 또는 특정 서열 맥락에서의 전위(translocation) 가변성 때문으로 돌려왔으나, 체계적인 왜곡을 일으키는 중요한 원인은 식별되지 않은 채 남아 있었다. 특히, 시퀀싱 중 나노포어 하단(trans compartment)에서 발생하는 RNA 이차 구조의 재폴딩(refolding) 현상이 베이스콜링(basecalling) 오류와 가짜 변형 신호를 유발하는 동인으로서 체계적으로 규명된 바 없었다.

방법론
저자들은 in vivo 전사체 분석과 통제된 in vitro 실험을 결합하여 RNA 폴딩이 DRS 신호에 미치는 영향을 분리하여 조사하였다.

생물종 및 샘플: 본 연구는 Bacillus subtilis(야생형 총 RNA 및 변형이 없는 in vitro 전사(IVT) 대조군), E. coli, 그리고 인간 세포주(HEK293T) 데이터셋을 활용하였다.
실험 설계: 인과관계를 확립하기 위해, 연구진은 B. subtilis 유전자에 기반한 특정 IVT 전사체를 설계하였다(cspD, yĬH, iolW, yxkC, kdgT). 연구진은 고유한 전사 종결 헤어핀(hairpin) 또는 내부 이차 구조가 다음과 같이 변화된 변이체들을 생성하였다:
- 파괴(Disrupted): 헤어핀의 원위부(3' arm)를 절단하거나("half"), 염기 쌍 형성을 방지하기 위해 점 돌연변이를 도입함("mut").
- 이동(Relocated): 3' 말단을 연장하여("long") 읽기 끝(read end)에 대한 헤어핀의 위치를 이동시킴.
- 재배치(Repositioned): 헤어핀 서열을 고품질 영역에 융합함("del" 및 "del mut").
시퀀싱 및 분석: 데이터는 PromethION 장치를 사용하여 R10.4.1 포어에서 시퀀싱되었다. 베이스콜링은 변형된 염기 추론 기능이 활성화된 Dorado(버전 5.0.0부터 6.0.0까지)를 사용하여 수행되었다. 신호 정렬 및 변형 호출을 위해 modkit, m6Anet, SWARM, Uncalled4와 같은 도구들이 사용되었다.
신호 분석: 저자들은 예측된 이차 구조(최소 자유 에너지, $\Delta G$ )와 관련하여 원시 이온 전류 수준, 체류 시간(dwell times, event lengths), 그리고 베이스 품질 점수(QS)를 분석하였다.

주요 결과

헤어핀 관련 아티팩트: in vivo 및 변형이 없는 IVT 샘물 모두에서, 고신뢰도 변형 호출(HCM) 및 베이스콜링 오류(미스매치, 인델)가 안정적인 RNA 헤어핀(고유 전사 종결자 및 강력한 내부 구조 포함)의 바로 상류(~12 nt)에서 일관되게 재현 가능한 창(window) 내에 농축되어 나타났다.
파괴를 통한 인과관계 입증: 절단 또는 점 돌연변이를 통해 이러한 헤어핀 형성을 파괴했을 때, 위양성 변형 호출의 축적이 사라지고 베이스 품질 점수가 배경 수준으로 회복되었다. 이는 다운스트림의 폴딩이 서열 맥락 자체가 아니라 아티팩트의 원인임을 확인시켜 주었다.
신호 섭동(Perturbation): 신호 수준 분석 결과, 다운스트림 헤어핀의 존재가 헤어핀 상류의 국소적인 이온 전류 편차와 체류 시간(dwell time)의 분산을 유발함을 발견하였다. 이러한 섭동은 나노포어 출구에서의 변경된 전위 역학 또는 입체적 간섭(steric interference)과 일치한다.
일반성: 이 현상은 다양한 생물종(B. subtilis, E. coli, 인간)과 다양한 변형 호출기(Dorado, SWARM)에서 관찰되었으며, 이는 소프트웨어 특정 버그가 아닌 일반적인 물리적 메커니즘임을 나타낸다.
UC-Rich 맥락: UC-rich 서열이 베이스콜링 정확도를 낮춘다는 것은 알려져 있으나, 본 연구는 다운스트림 헤어핀이 UC-rich 모티프가 존재하는 맥락에서 아티팩트를 증폭시키는 독립적인 동인임을 밝혔다. 그러나 강력한 헤어핀은 UC-rich 모티프가 없는 경우에도 신호 왜곡을 일으킨다.
모델의 한계: 여러 Dorado 모델 버전(v5.0.0 to v6.0.0)을 분석한 결과, 일부 업데이트가 특정 아티팩트를 줄이기는 했으나, 구조와 연관된 근본적인 왜곡은 지속되었으며 새로운 모델에서도 변형이 없는 RNA에서 상당한 위양성 신호가 여전히 발생하였다.

의의 및 주장
본 논문은 trans 구획에서의 RNA 재폴딩이 DRS에서 이전에 인식되지 않았던 일반적인 신호 왜곡의 원인임을 밝혀냈다. 저자들은 이 메커니즘이 특히 안정적인 이차 구조의 상류에서 시퀀싱 충실도를 체계적으로 저해하고 위양성 변형 호출을 생성한다고 주장한다.

연구는 강력한 헤어핀의 상류에 위치한 후보 변형 부위는 직교 검증(예: 질량 분석법)이나 엄격한 음성 대조군(negative control)에 의해 뒷받침되지 않는 한 주의해서 해석되어야 함을 강조한다. 저자들은 세 가지 완화 전략을 제안한다:

변형 호출 알고리즘에 구조 인식 품질 관리(structure-aware QC)를 구현하여, 예측된 강력한 헤어핀의 약 ~15 nt 상류에 있는 베이스를 플래그 처리하거나 가중치를 낮춤.
정의된 헤어핀을 포함하는 합성 RNA 구조체 또는 구조 인식 데이터가 포함된 비변형 IVT 전사체를 사용하여 향로의 베이스콜링 모델을 훈련시킴.
플로우 셀(flow cell)의 trans 구획에서 헤어핀 형성을 물리적으로 억제함(예: 화학적 첨가제 또는 변성제 사용), 단 이는 정밀한 최적화가 필요함.

궁극적으로, 본 연구는 DRS의 맥락 의존적 오류에 대한 기계론적 설명을 제공하며, 에피전스크립토믹(epitranscriptomic) 연구 및 임상 적용의 신뢰성을 향상시키기 위한 구조 인식 데이터 해석의 필수적인 토대를 마련하였다.

유사한 논문