레이블이 Nanopore인 게시물을 표시합니다. 모든 게시물 표시
레이블이 Nanopore인 게시물을 표시합니다. 모든 게시물 표시

금요일, 7월 24, 2026

작아서 안 보이면 키워서 보자: DNA를 50배 늘려서 읽는 SBX 기술

오늘은 DNA 서열 분석에 대한 재미있는 논문(Sequencing by Expansion (SBX) – a novel, high-throughput single-molecule sequencing technology)을 하나 가져와봤습니다. 우리 몸의 설계도인 DNA를 읽어내는 기술들이 계속 발전하고 있는데, 기존에는 아주 작은 구멍(나노포어)에 DNA를 통과시키면서 읽다 보니 신호가 너무 작고 휙 지나가버려서 에러가 나거나 겹치는 문제가 꽤 있었거든요.그래서 이번 연구진들은 원래 DNA를 억지로 읽으려 애쓰는 대신, 이걸 50배 정도로 길게 쫙 늘린 특수한 복사본(엑스판도머)을 만들어서 훑어보는 방식을 개발했다고 하네요. 제목은 Sequencing by Expansion 으로, 작아서 읽기 힘들면 그냥 길게 늘려서 읽으면 되지 않냐는.. "안 보이면 돋보기 쓰자" 마인드의 발상의 전환 되겠습니다.

기존처럼 DNA를 직접 측정하며 고생하지 않고, 기계가 읽기 편한 형태로 먼저 변환해서 정확도를 높였다는 점이 꽤 흥미롭네요. 아래는 여러 AI들이 이 논문에 대해 각자의 스타일로 정리해준 내용들이니 커피 한잔 하시면서 편하게 읽어보시면 좋지 않을까 합니다. :)


※ kimi결과가 살짝 아쉽네요 좀 바쁜가봅니다;; 

gemini
**[블로그 업로드용] 세상을 바꾸는 DNA 읽기 기술, '확장 시퀀싱(SBX)'을 아시나요?**
안녕하세요! 오늘은 생명과학계의 거대한 패러다임을 바꿀 수 있는 흥미진진하고 혁신적인 최신 연구 한 편을 소개해 드리려고 합니다. 바로 최근 바이오 분야에서 뜨겁게 주목받고 있는 **‘확장 시퀀싱(Sequencing by Expansion, 이하 SBX)’** 기술에 대한 이야기입니다. 박사 학위와 현장 실무 경험을 바탕으로, 중학생도 단번에 이해할 수 있도록 최대한 쉽고 친근하게 풀어내 보았으니 편하게 읽어보세요!
---
### 1. 연구 배경: 기존 DNA 읽기 기술의 숨은 ‘장벽’
우리 몸의 설계도라고 불리는 DNA는 A(아데닌), T(티민), G(구아닌), C(사이토신)라는 네 가지 화학 기호(염기)가 아주 길게 나열된 형태입니다. 이 순서를 정확하게 읽어내는 기술을 'DNA 시퀀싱'이라고 부르며, 현대 의학이나 암 진단, 유전병 연구의 핵심 기반이 됩니다.
최근 과학자들은 아주 미세한 나노 크기의 구멍(나노포어, Nanopore)에 DNA를 통과시키면서 발생하는 미세한 전류 변화를 감지해 염기를 읽는 방식을 개발했습니다. 마치 카세트테이프가 재생기를 통과할 때 소리가 나는 것처럼 말이죠. 하지만 치명적인 문제가 있었습니다. DNA 분자가 너무 작고 촘촘하다 보니 구멍을 지날 때 신호 구분이 잘 안 되어 전기적 잡음(노이즈)이 심했고, 결과적으로 읽기 오류가 자주 발생한다는 점이었습니다.
### 2. 연구 목적: "DNA를 늘려서 크게 읽으면 어떨까?"
이번 연구의 목적은 바로 이 고질적인 '신호 대 잡음비' 한계를 깨부수는 것이었습니다. 연구진은 복잡하게 뭉쳐 있고 해상도가 낮은 원래 DNA를 직접 읽으려 애쓰는 대신, "DNA 정보는 그대로 유지하되, 크기를 수십 배 뻥튀기하고 신호가 잘 나는 특수한 분자로 변환해서 읽자!"라는 획기적인 아이디어를 제안했습니다. 즉, 글씨가 너무 작아 읽기 힘들다면 '확대 복사'를 해서 읽기 쉽게 만들겠다는 전략입니다.
### 3. 연구 방법: 첨단 화학과 효소 공학이 만들어낸 3단계 마법
이 놀라운 변환 과정을 현실로 만들기 위해 연구팀은 다음과 같은 정밀한 실험 과정을 거쳤습니다.
* **1단계: 특수 벽돌(XNTP) 합성하기**
연구진은 인공적으로 조작된 네 가지 형태의 특수한 핵산 벽돌(XNTP)을 만들었습니다. 이 벽돌들에는 각각 나노 구멍을 통과할 때 아주 뚜렷하고 깨끗한 전기 신호를 발생시키는 ‘리포터(신호 발생기)’가 꼬리표처럼 달려 있습니다.
* **2단계: 인공 DNA 복사본(Xpandomer) 만들기**
자연계에 존재하는 일반적인 효소는 이 커다란 인공 벽돌들을 연결하지 못합니다. 그래서 연구진은 수많은 돌연변이 실험을 거쳐 이 특수 벽돌들을 척척 이어 붙일 수 있는 맞춤형 합성 효소(Xp Synthase)와 효소를 돕는 보조 물질(PEM)을 개발했습니다. 이를 통해 기존 DNA 템플릿을 바탕으로 한 인공 복사본 분자인 '엑스팬도머(Xpandomer)'를 성공적으로 합성해 냈습니다.
* **3단계: 산성 용액으로 50배 늘리기**
이렇게 만들어진 엑스팬도머를 약한 산성 조건에 노출하면, 결합 부위가 풀리면서 원래 DNA보다 **무려 50배 이상 길게 쭉 늘어나는 ‘확장’ 현상**이 일어납니다. 길이가 늘어나면서 리포터 꼬리표들도 일렬로 길게 간격을 두고 배치되어, 나노 구멍이 하나씩 신호를 명확하게 읽을 수 있는 최적의 상태가 됩니다.
### 4. 연구 결과: 완벽하게 분리된 네 가지 신호와 높은 정확도
연구팀은 실제 '스트레프토코쿠스(Streptococcus)'라는 박테리아의 222개 염기 길이 DNA를 사용해 이 기술을 검증했습니다.
결과는 대성공이었습니다. 50배 길어진 인공 분자가 나노 구멍을 통과할 때, A, T, G, C 각각의 리포터가 내는 전류 신호가 그래프 상에서 겹치지 않고 완전히 독립된 네 개의 뚜렷한 높낮이(피크)로 측정되었습니다. 이 덕분에 복잡한 계산 없이도 신호를 직관적으로 바코드 읽듯이 해석할 수 있게 되었습니다. 실험 데이터 분석 결과, 단 한 번만 읽었음에도 불구하고 **평균 98.3%라는 압도적인 높은 정확도**로 DNA 서열을 정확하게 맞추는 쾌거를 이루어냈습니다.
### 5. 고찰: 한계를 극복하기 위한 여정과 과제
성공적인 결과 뒤에는 정교한 시스템 엔지니어링의 노력이 숨어 있었습니다. 연구진은 인공 분자가 나노 구멍을 지나갈 때 간혹 구멍에 끼거나 멈추는 '막힘 현상(Blockage)'이 전체 측정 시간의 약 46%를 차지한다는 점을 확인했습니다.
이를 해결하기 위해 순간적으로 전압의 방향을 반대로 바꾸어 걸린 분자를 튕겨내 주는 '이젝트 웨이브폼(Eject waveform)' 기술을 적용하여 연속적인 측정이 가능하도록 시스템을 최적화했습니다. 또한, 합성 효소와 보조 물질, 그리고 신호 리포터가 완벽한 삼박자를 이루어야만 중간에 끊기지 않고 긴 DNA를 완전히 변환할 수 있다는 유기적 메커니즘을 밝혀냈습니다.
### 6. 의의와 시사점: 화학과 측정의 '완전한 분리'
이 연구가 가지는 가장 무서운 가치는 "화학적 변환 단계"와 "물리적 측정 단계"를 완전히 분리(Decoupling)해 놓았다는 점입니다.
기존 기술들은 DNA를 합성하면서 동시에 실시간으로 읽어야 했기 때문에, 화학 반응 조건과 기계의 측정 조건을 모두 만족시키느라 이도 저도 아닌 타협을 해야만 했습니다. 하지만 SBX 기술은 실험실 칩 위에서 안전하고 정밀하게 대량의 DNA를 긴 인공 분자로 완벽히 바꾸어 놓은 뒤, 나노포어 기계에서는 오직 신호를 읽는 데만 집중할 수 있게 설계되었습니다. 따라서 앞으로 기계 성능이 발전하거나 새로운 리포터 화학 물질이 개발되면, 서로 방해받지 않고 각각 독립적으로 기술을 무한히 업그레이드할 수 있는 확장성을 확보하게 되었습니다.
---
### 💡 이 연구가 왜 중요한가요?
과거에는 한 사람의 전체 DNA 정보를 읽는 데 수년의 시간과 수십억 원의 비용이 들었습니다. 지금은 비용이 많이 낮아졌지만, 여전히 대형 병원이나 전문 연구소의 값비싼 장비와 복잡한 분석 과정을 거쳐야만 합니다.
SBX 기술이 상용화되어 수백만 개의 나노 구멍이 뚫린 초소형 센서 칩과 결합하면, **마치 스마트폰만 한 조그만 기기에 장착해 몇 분 만에 저렴하고 완벽한 개인 유전자 지도를 완성**하는 시대가 열립니다. 스마트폰 센서로 혈당을 재듯, 아주 적은 비용으로 언제 어디서나 암 세포의 미세한 변이를 실시간으로 추적하거나 신종 바이러스 감염 여부를 초고속으로 정확하게 진단할 수 있게 되는 것이죠. 결국 이 연구는 미래 인류가 질병을 극복하고, 맞춤형 정밀 의료를 일상화하는 데 있어 가장 단단하고 획기적인 ' 디딤돌'을 놓았다는 점에서 엄청난 가치를 지니고 있습니다.

chatGPT
# DNA를 직접 읽지 않고 더 크게 펼쳐 읽는다? 차세대 유전체 분석기술 'SBX' 이야기
유전체 분석 기술은 지난 20여 년 동안 눈부시게 발전했습니다. 예전에는 한 사람의 DNA를 분석하는 데 수년의 시간과 수천억 원이 필요했지만, 지금은 며칠 안에 훨씬 저렴한 비용으로 분석할 수 있습니다. 이러한 발전 덕분에 암 연구, 희귀질환 진단, 감염병 분석, 신약 개발 등 다양한 분야에서 유전체 분석은 없어서는 안 될 핵심 기술이 되었습니다. 하지만 연구가 점점 대규모로 진행되면서 더 빠르고, 더 정확하며, 더 저렴한 분석 기술에 대한 요구도 커지고 있습니다. 이번에 소개하는 연구는 이러한 한계를 해결하기 위해 기존 방식과는 완전히 다른 개념의 새로운 염기서열 분석 기술인 **Sequencing by Expansion(SBX)**를 개발한 연구입니다.
기존의 나노포어(Nanopore) 시퀀싱은 DNA 분자를 아주 작은 구멍인 나노포어에 통과시키면서 전류의 변화를 측정해 염기서열을 읽습니다. 긴 DNA를 한 번에 읽을 수 있다는 큰 장점이 있지만, DNA를 직접 측정하기 때문에 네 가지 염기(A, T, G, C)가 만들어 내는 전기 신호의 차이가 매우 작습니다. 또한 DNA가 빠르게 이동하면서 신호가 서로 겹치는 문제가 발생해 정확도가 떨어질 수 있으며, 같은 염기가 반복되는 구간에서는 오류가 증가하는 한계도 있습니다. 연구진은 이러한 문제의 근본 원인이 DNA 자체를 직접 측정하는 방식에 있다고 판단했습니다.
이 연구의 목적은 DNA를 직접 읽는 대신, DNA의 정보를 훨씬 측정하기 쉬운 새로운 분자로 바꾼 뒤 이를 분석하는 새로운 기술을 개발하는 것이었습니다. 연구진은 이를 **Sequencing by Expansion(SBX)**라고 이름 붙였습니다. 핵심 아이디어는 DNA의 염기서열을 그대로 유지하면서 약 50배 이상 길게 확장된 인공 분자인 **Xpandomer(Xp)**를 만든 후, 이 분자를 나노포어에서 읽는 것입니다. 이렇게 하면 서로 다른 염기를 훨씬 뚜렷한 신호로 구분할 수 있어 정확도를 높일 수 있을 것으로 기대했습니다.
이를 위해 연구진은 여러 가지 새로운 기술을 동시에 개발했습니다. 먼저 일반적인 DNA 재료 대신 **XNTP**라는 특수한 염기를 새롭게 설계했습니다. 이 염기에는 각각의 DNA 염기를 구별할 수 있는 신호가 미리 포함되어 있으며, 이후 화학 처리를 거치면 긴 구조로 펼쳐질 수 있도록 만들어졌습니다. 또한 이러한 특수 염기를 안정적으로 연결할 수 있도록 기존 DNA 중합효소를 크게 개량한 **Xp Synthase**라는 새로운 효소도 개발했습니다. 여기에 DNA 합성을 더욱 길고 안정적으로 진행하도록 돕는 **PEM(Polymerase Enhancing Moiety)**이라는 보조 분자까지 추가하여 수백 개 이상의 염기를 연속적으로 합성할 수 있도록 만들었습니다.
SBX는 크게 두 단계로 진행됩니다. 첫 번째는 DNA를 Xpandomer로 변환하는 단계입니다. DNA를 고정된 기판 위에서 특수 효소와 XNTP를 이용해 복제한 뒤, 화학 처리를 통해 각각의 염기를 길게 펼쳐 줍니다. 이렇게 만들어진 Xpandomer는 원래 DNA보다 훨씬 길어지고, 각 염기가 서로 명확하게 구분되는 신호를 갖게 됩니다. 두 번째 단계에서는 완성된 Xpandomer를 나노포어에 통과시키면서 전류 변화를 측정합니다. 연구진은 각 염기가 일정한 위치에서 잠시 멈추도록 설계하여 하나씩 안정적으로 측정할 수 있도록 만들었습니다. 이 과정에서 네 가지 염기는 서로 다른 전류 신호를 나타내므로 컴퓨터는 이를 다시 원래의 DNA 염기서열로 변환할 수 있습니다. 기존처럼 DNA를 직접 읽는 것이 아니라, 읽기 쉬운 형태로 변환한 후 분석한다는 점이 SBX의 가장 큰 특징입니다.
연구진은 새롭게 개발한 기술을 검증하기 위해 222개의 염기로 이루어진 세균 DNA를 이용해 실제 시퀀싱 실험을 수행했습니다. 총 966개의 읽기(read)를 분석한 결과 약 26%는 전체 길이를 끝까지 읽는 데 성공했으며, 중앙값 기준 읽기 길이는 181염기로 나타났습니다. 무엇보다 전체 염기 판독 정확도는 약 **98.3%**를 기록했습니다. 또한 반복 염기나 머리핀(hairpin) 구조처럼 기존 나노포어 기술에서 오류가 자주 발생하는 구간도 안정적으로 분석할 수 있음을 확인했습니다. 이는 SBX가 단순한 아이디어 수준이 아니라 실제 염기서열 분석 기술로 충분한 가능성을 가지고 있음을 보여주는 중요한 결과입니다.
연구진은 이번 연구를 통해 DNA를 직접 측정해야 한다는 기존의 고정관념을 바꾸었습니다. DNA를 먼저 측정하기 쉬운 형태로 변환한 뒤 분석하는 방식은 신호를 더욱 선명하게 만들고, 정확도를 높이며, 화학 반응과 측정 과정을 각각 독립적으로 최적화할 수 있다는 장점이 있습니다. 특히 향후 수백만 개의 나노포어를 동시에 사용하는 대규모 시스템과 결합된다면 전장 유전체 분석, 단일세포 RNA 분석, 공간전사체 분석 등 대용량 데이터가 필요한 다양한 분야에서 기존보다 훨씬 빠르고 경제적인 분석이 가능할 것으로 기대하고 있습니다.
물론 이번 연구는 아직 bioRxiv에 공개된 프리프린트(preprint) 논문으로, 동료 심사를 모두 마친 연구는 아닙니다. 또한 실제 임상 현장이나 대규모 연구에서 기존 상용 플랫폼을 완전히 대체하기 위해서는 더 긴 DNA에서의 성능, 비용, 대량 생산성, 안정성 등에 대한 추가 검증이 필요합니다. 그럼에도 불구하고 DNA를 직접 읽는 기존 방식의 한계를 새로운 분자공학 기술로 해결하려는 접근은 매우 독창적이며, 차세대 시퀀싱 기술의 새로운 방향을 제시했다는 점에서 큰 의미를 갖습니다.
이 연구가 중요한 이유는 단순히 새로운 장비를 개발했기 때문이 아닙니다. 연구진은 DNA를 그대로 읽으려 하지 않고, **읽기 쉬운 형태로 먼저 바꾼 뒤 분석한다**는 새로운 개념을 제안했습니다. 이러한 발상의 전환은 앞으로 유전체 분석의 정확도와 속도를 동시에 높이고 비용까지 낮출 가능성을 보여줍니다. 만약 SBX 기술이 실제 상용화된다면 암의 조기 진단, 희귀질환의 유전자 검사, 감염병 감시, 정밀의료, 신약 개발 등 다양한 분야에서 더 많은 사람들에게 빠르고 정확한 유전체 분석 서비스를 제공할 수 있을 것으로 기대됩니다. 즉, 이번 연구는 단순한 기술 개발을 넘어 미래 유전체 분석 기술의 새로운 패러다임을 제시한 매우 의미 있는 연구라고 할 수 있습니다.

grok
**DNA 염기서열 읽는 새로운 기술, SBX 이야기**
DNA 서열 분석은 현대 생물학과 의학에서 정말 중요한 기술입니다. 우리 몸의 유전 정보를 읽어 질병 원인을 찾거나 새로운 치료법을 개발하는 데 쓰이죠. 하지만 기존 기술은 아직 비용이 비싸고 속도가 느려서, 더 많은 사람들이 쉽게 이용하기 어려웠습니다. 특히 나노포어(nanopore) 방식은 DNA를 직접 구멍으로 통과시키며 전기 신호로 읽는 방식인데, 신호가 약하고 오류가 많아 실용화가 힘들었습니다.
이 논문에서는 'Sequencing by Expansion(SBX)'라는 새로운 기술을 소개합니다. SBX는 DNA를 직접 읽지 않고, 먼저 DNA를 훨씬 길고 쉽게 읽을 수 있는 'Xpandomer(엑스판도머, Xp)'라는 특별한 분자로 바꾸는 아이디어입니다. 마치 작은 글씨를 큰 글씨로 확대해서 읽는 것처럼, DNA 정보를 확장시켜 더 정확하게 측정할 수 있게 만든 거예요.
연구팀은 2007년부터 이 기술을 개발해왔고, 로슈(Roche)에서 이어가고 있습니다. 목적은 저렴하고 빠르며 정확한 단일분자 서열 분석 기술을 만드는 것입니다. 이를 통해 유전체 연구, 암 진단, 개인 맞춤 의학 등이 더 널리 퍼지길 바랍니다.
방법은 두 단계로 나뉩니다. 먼저 합성 단계에서 DNA 템플릿에 특수한 효소(Xp Synthase)와 변형된 핵산(XNTPs)을 이용해 Xp를 만듭니다. 여기에 PEM이라는 보조 물질을 넣어 긴 사슬을 잘 만들도록 돕습니다. 이 과정은 고체 지지체(작은 칩) 위에서 진행되어 불순물을 쉽게 제거할 수 있어요. 만든 Xp는 산 처리로 확장되고, 나노포어로 이동합니다.
측정 단계에서는 알파-헤몰리신 단백질 구멍을 이용합니다. Xp가 구멍을 통과할 때 특정 구조(TCE)가 멈추게 해서 각 염기 신호를 1밀리초 동안 정확히 읽습니다. 전압 펄스를 주며 한 칸씩 이동시키는 방식으로, 네 가지 염기(C, T, A, G)를 뚜렷한 전류 차이로 구분합니다.
결과는 놀라웠습니다. 222개 염기로 된 세균 DNA를 테스트한 결과, 10분 동안 966개의 읽기를 얻었고 평균 정확도가 98.3%에 달했습니다. 전체 길이의 26%가 완벽하게 읽혔고, 중앙값 읽기 길이는 181개 염기였어요. 기존 나노포어 기술보다 오류가 적고 신호가 선명했습니다.
고찰에서 연구팀은 Xp 합성과 측정 과정을 분리한 점이 큰 장점이라고 강조합니다. 각 부분을 독립적으로 최적화할 수 있어서 미래 발전이 빠를 거예요. 아직 블록킹(막힘) 같은 기술적 과제가 있지만, 8백만 개 구멍을 가진 대형 칩으로 확대하면 초당 수억 개 염기를 읽을 수 있을 전망입니다.
이 연구의 의의는 DNA 서열 분석의 한계를 크게 뛰어넘었다는 데 있습니다. 기존 기술보다 빠르고 정확하며 비용을 낮출 잠재력이 크기 때문에, 연구실뿐 아니라 병원과 일상에서 유전자 검사가 보편화될 수 있어요. 특히 실시간 분석이 가능해져 응급 진단이나 현장 검사에 유리합니다.
마지막으로, 이 연구가 왜 중요한가요? 인간 게놈 프로젝트 이후로 우리는 유전 정보를 더 많이, 더 빠르게 읽어야 새로운 생물학적 발견과 의료 혁신을 이룰 수 있습니다. SBX는 그 꿈을 한 걸음 더 가까이 가져온 기술로, 암 조기 발견, 희귀질환 진단, 개인화 의학 시대를 앞당길 수 있습니다. 과학자들이 수년간 쌓아온 분자 설계와 효소 공학의 결실이 빛을 발하는 순간입니다. 앞으로 더 발전된 버전이 나오길 기대하며, 일반인으로서도 이런 혁신이 우리 삶을 어떻게 바꿀지 지켜보는 게 흥미로울 것 같아요.

NotebookLM
안녕하세요! 유전체 해독 기술의 최전선에서 연구와 실무를 병행해 온 전문가입니다. 오늘은 2025년 발표된 매우 혁신적인 논문을 하나 소개해 드리려고 합니다. 우리가 흔히 '생명의 설계도'라고 부르는 DNA를 읽어내는 방식을 완전히 뒤바꿀 **'확장 해독 기술(SBX: Sequencing by Expansion)'**에 대한 이야기입니다. 중학생도 "아, 이렇게 DNA를 늘려서 읽는구나!"라고 이해할 수 있게 핵심 위주로 쉽게 설명해 드릴게요.
**[블로그 포스팅] DNA를 50배로 늘려서 읽는다고? '확장 해독(SBX)'이 여는 초고속 유전자 분석 시대**
우리 몸의 설계도인 DNA를 읽어내는 기술은 현대 의학에서 아주 중요합니다. 하지만 지금까지의 기술은 DNA 조각이 너무 작고 움직임이 빨라서 정확하게 읽어내는 데 한계가 있었습니다. 오늘 소개할 연구는 이 작은 DNA 정보를 무려 50배나 길게 늘린 '복제본'을 만들어, 마치 돋보기를 대고 보듯 선명하게 읽어내는 놀라운 기술입니다.
**1. 연구 배경: 너무 작고 빠른 DNA, 읽기가 너무 힘들어요**
기존에 기대를 모았던 '나노포어(Nanopore)' 기술은 아주 작은 구멍에 DNA를 통과시키며 흐르는 전류의 변화로 염기(A, T, G, C)를 구분하는 방식이었습니다. 하지만 DNA 분자 자체가 너무 작고 구멍을 통과하는 속도가 일정하지 않아, 신호가 뭉개지거나(노이즈) 오류가 생기는 문제가 컸습니다. 비유하자면, 아주 빠르게 달리는 기차의 창문에 적힌 작은 글씨를 읽어야 하는 상황과 같았죠.
**2. 연구 목적: DNA를 읽기 쉬운 '확장판'으로 변신시키자!**
연구진은 이 문제를 해결하기 위해 **'확장 해독(Sequencing by Expansion, SBX)'**이라는 새로운 개념을 도입했습니다. 목표는 DNA 정보를 그대로 담고 있으면서도, 크기는 50배 더 길고 신호는 훨씬 뚜렷한 **'엑스판도머(Xpandomer, Xp)'**라는 가짜 폴리머를 만드는 것이었습니다. 이렇게 하면 기존 기술의 한계였던 신호 대 잡음비(Signal-to-noise) 문제를 근본적으로 해결할 수 있기 때문입니다.
**3. 연구 방법: 특수 효소와 '일시 정지' 버튼의 만남**
이 기술의 핵심은 크게 세 단계입니다.
첫째, **'Xp 신타아제(Xp Synthase)'**라는 특수 설계된 효소를 사용합니다. 이 효소는 원래의 DNA를 본떠서 커다란 부품(XNTP)들을 하나씩 조립해 '엑스판도머'를 만듭니다. 
둘째, 조립이 끝나면 산(Acid) 처리를 해서 이 엑스판도머를 **50배 길이로 쫙 펼칩니다.** 
셋째, 나노포어 구멍에 이 긴 분자를 통과시킵니다. 이때 **'전이 조절 장치(TCE)'**라는 일종의 '일시 정지' 버튼 역할을 하는 구조물이 각 염기 정보를 구멍 속에 1밀리초(ms) 동안 딱 멈춰 세워줍니다. 덕분에 컴퓨터는 각 염기가 무엇인지 아주 정확하게 읽어낼 수 있습니다.
**4. 주요 연구 결과: 98.3%의 높은 정확도와 엄청난 잠재력**
연구팀은 222개의 염기로 이루어진 미생물 DNA를 대상으로 이 기술을 시험했습니다. 결과는 매우 성공적이었습니다. 
첫째, **98.3%라는 높은 정확도**로 DNA 서열을 읽어냈습니다. 이는 기존의 나노포어 방식이 겪던 신호 혼란을 완벽하게 극복했음을 보여줍니다. 
둘째, 복잡하게 꼬인 구조나 똑같은 글자가 반복되는 구간도 엑스판도머로 길게 늘린 덕분에 실수 없이 판독해 냈습니다. 
셋째, 한 번에 처리하는 속도가 매우 빨라, 이론적으로는 **초당 5억 개의 염기**를 읽어낼 수 있는 대규모 장치로 확장 가능하다는 것을 증명했습니다.
**5. 고찰 및 시사점: '화학'과 '측정'의 완벽한 분리**
이 연구가 기존 기술과 차별화되는 점은 DNA를 복제하는 '화학적 과정'과 신호를 읽어내는 '측정 과정'을 완전히 분리했다는 것입니다. 이전에는 두 과정이 동시에 일어나서 하나를 개선하면 다른 하나가 나빠지는 문제가 있었지만, SBX 기술은 각각을 따로 최적화할 수 있습니다. 예를 들어, 더 정확한 효소를 개발하거나 더 정밀한 센서를 만드는 작업을 독립적으로 진행할 수 있어 기술 발전 속도가 훨씬 빨라질 것입니다.
**6. 의의와 시사점: 누구나 저렴하고 빠르게 유전자를 읽는 세상**
SBX 기술은 유전자 해독의 문턱을 크게 낮출 것입니다. 복잡한 슈퍼컴퓨터 분석이나 비싼 장비 없이도, 반도체 칩(CMOS) 기술과 결합해 아주 작고 빠른 해독기를 만들 수 있기 때문입니다. 이는 병원 현장에서 환자의 암 유전자를 즉석에서 확인하거나, 전염병의 정체를 단 몇 분 만에 밝혀내는 데 핵심적인 역할을 할 것입니다.
**💡 이 연구가 왜 중요한가요?**
이 연구는 **'생물학적 정보를 디지털 정보로 바꾸는 가장 선명한 번역기'**를 만들었기 때문에 중요합니다. 지금까지 우리는 희미한 신호를 인공지능으로 억지로 해석하느라 고생했지만, 이제는 DNA 자체를 '크고 선명하게' 만들어 버림으로써 근본적인 해결책을 찾았습니다. 이 기술이 상용화되면 개인 맞춤형 의료가 훨씬 저렴해지고 보편화될 것입니다. 결국 보이지 않는 생명의 암호를 누구나 쉽게 읽을 수 있는 '유전자 고속도로'를 닦았다는 점이 이 연구의 가장 큰 가치입니다.

Perplexity.ai
이 논문은 DNA를 그대로 읽는 대신, DNA 정보를 더 잘 읽히는 “확장된 분자”로 바꾼 뒤 나노포어로 분석하는 새 시퀀싱 기술, SBX(Sequencing by Expansion)를 소개합니다. 핵심은 원래 DNA를 직접 재는 방식이 가진 낮은 신호 문제를 줄여, 더 정확하고 빠르게 유전정보를 읽게 하려는 데 있습니다.
## 연구 배경
유전체 분석은 생명과학 연구와 진단에서 매우 중요해졌지만, 더 넓게 쓰이려면 아직도 비용과 정확도, 속도 문제가 남아 있습니다. 특히 나노포어 시퀀싱은 한 분자씩 읽을 수 있다는 장점이 있지만, DNA를 직접 읽는 과정에서 신호가 충분히 선명하지 않아 정확도가 떨어지는 한계가 있었습니다. 이 연구는 바로 그 한계를 해결하기 위해 출발했습니다.
## 연구 목적
이 연구의 목적은 DNA의 염기서열 정보를 더 잘 구분되는 형태로 바꾼 뒤, 나노포어에서 쉽게 읽을 수 있게 만드는 것이었습니다. 즉, “DNA를 잘 읽는 방법”이 아니라 “DNA를 읽기 쉬운 다른 분자로 바꿔서 읽는 방법”을 제시한 것입니다. 연구진은 이 기술이 고정밀·고속 시퀀싱에 도움이 될 수 있는지 보여주고자 했습니다.
## 연구 방법
연구진은 DNA를 바로 측정하지 않고, 먼저 DNA 정보를 Xpandomer(Xp)라는 확장된 분자로 바꾸는 과정을 만들었습니다. 이 과정에서 네 가지 염기에 맞는 특수한 XNTP, Xp Synthase 효소, 그리고 변환을 돕는 보조 물질(PEM)을 사용해 DNA 정보를 길고 잘 구분되는 분자로 확장했습니다. 그 다음 이 Xp를 나노포어에 통과시켜 전류 변화 신호를 읽고, 그 패턴으로 원래 DNA 염기서열을 해석했습니다.
## 연구 결과
이 기술은 DNA보다 50배 이상 길어진 Xpandomer를 만들어낼 수 있었고, 나노포어에서 각 염기를 구분할 수 있는 선명한 전류 신호를 보여주었습니다. 연구진은 Xp 생성과 측정이 서로 분리된 구조 덕분에, 각각의 과정을 따로 최적화할 수 있었다고 설명합니다. 즉, 읽기 어려운 DNA를 “더 읽기 쉬운 형태”로 바꿔 높은 신호대잡음비를 확보한 것이 가장 큰 성과입니다.
## 고찰
이 연구의 의미는 단순히 새 장비를 만든 것이 아니라, 시퀀싱의 약점을 분자 설계와 공정 설계로 해결하려 했다는 점입니다. 직접 DNA를 재는 대신, 정보를 한 번 다른 분자로 옮겨 읽는 방식이라서 신호가 훨씬 안정적이고 해석이 쉬워졌습니다. 다만 이 논문은 아직 기술 소개와 원리 검증 성격이 강하므로, 실제 임상이나 대규모 현장 적용에서는 속도, 비용, 재현성, 시스템 규모 확대가 더 검증되어야 합니다.
## 의의와 시사점
이 기술이 널리 발전하면 유전체 검사 비용을 낮추고, 더 많은 기관에서 빠르고 정확한 DNA 분석을 할 수 있는 길이 열릴 수 있습니다. 특히 질병 진단, 대규모 인구집단 분석, 단일세포 연구처럼 많은 데이터를 빠르게 읽어야 하는 분야에서 활용 가능성이 큽니다. 쉽게 말해, 이 연구는 “DNA를 읽는 기술”의 속도와 정확도를 한 단계 끌어올릴 수 있는 새 방향을 제시했다는 점에서 중요합니다.
이 연구가 중요한 이유는, 생명과학과 의료에서 가장 기본이 되는 정보인 DNA를 더 정확하고 효율적으로 읽게 해 주기 때문입니다.

Kimi
=== SBX 논문 핵심 데이터 ===
1. 연구 배경:
   - DNA 시퀀싱은 생물학/의학 연구의 핵심 도구
   - 기존 나노포어 시퀀싱(ONT)은 직접 DNA를 읽어 raw error rate >3%, 호모폴리머 분해능 낮음
   - 신호 대 잡음(signal-to-noise) 문제가 근본적 한계
2. SBX 개념:
   - DNA → Xpandomer(Xp)라는 확장된 서로게이트 폴리머로 변환
   - Xp는 원본 DNA보다 50배 이상 길어짐
   - 4가지 염기(C,T,A,G)마다 고유한 리포터(신호 패턴) 부착
3. 핵심 구성요소:
   - XNTP: 확장 가능한 뉴클레오타이드(4종류)
   - SSRT: 대칭 합성 리포터 테더
   - TCE: 이동 제어 요소(나노포어에서 1ms 정지)
   - PEM: 폴리머라제 증강 부위(400베이스 이상 합성 가능하게 함)
   - Xp Synthase: Dpo4 변형 효소(36개 아미노산 치환)
   - EO: 확장 올리고(리더+컨센트레이터+프라이머)
4. 방법:
   - 고체상(solid-phase) 합성: cyclic ole린 칩 사용
   - 산 처리로 P-N 결합 절단 → Xp 확장
   - alpha-hemolysin 나노포어, DPhPE 지질 이중막
   - 전압 파형: -70mV 측정(1ms) + -600mV 펄스(6μs) 반복
5. 결과 (222베이스 Streptococcus 템플릿):
   - 총 966개 읽기(≥10베이스)
   - 평균 정확도: 98.3%
   - Insertion: 0.56%, Deletion: 0.36%, Substitution: 0.76%
   - 26%가 전체 222베이스 풀랭스
   - 중간 읽기 길이: 181베이스
   - 처리량: 28,038베이스/분 (단일 나노포어)
   - 1kHz 펄스 → 1kb/s 판독 속도
6. 비교:
   - ONT MinION R10.4: hac 95.6%, sup 96.8% (필터링 후)
   - SBX: 98.3% (단순 4-범위 분류기, ≥10베이스 필터)
   - 8M 나노포어 어레이 확장 시 >500M베이스/초 예상
7. 의의:
   - 화학(합성)과 측정(나노포어) 분리 → 독립 최적화
   - 실시간(near real-time) 데이터 출력 가능
   - 고처리량 CMOS 어레이와 결합 가능
   - 인구규모 유전체, 단세포 RNA-seq, 공간 시퀀싱 등 응용 가능
System is currently busy. Please try again later.




출처: @ye._.vely618


토요일, 4월 13, 2024

ONT, cfDNA 그리고 Methylation

나노포어를 이용해서 cell-free DNA의 메틸레이션 패턴을 분석한 논문이 있어서 한번 가져와봤습니다.
이번에는 제목은 어떻게 쓰면 좋을지 한번 물어보았습니다.

제목 답변

해당 연구의 제목은 "Detecting cell-of-origin and cancer-specific methylation features of cell-free DNA from Nanopore sequencing"입니다. 이 제목은 세포유리 DNA의 메틸화 패턴을 분석하여 암세포에서 유래한 세포유리 DNA를 식별하는 방법을 개발한 연구라는 것을 나타내고 있습니다. 이 연구는 생명공학 분야에서 매우 중요한 연구 중 하나이며, 세포유리 DNA의 메틸화 패턴을 분석하는 기술은 암 진단과 치료에 큰 도움을 줄 수 있습니다.
"세포유리 DNA의 메틸화 패턴을 이용한 암 진단 연구"와 같은 제목을 사용할 수 있습니다. 이 글에서는 연구의 목적과 방법, 결과 등을 소개하고, 세포유리 DNA의 메틸화 패턴이 암 진단에 어떻게 활용될 수 있는 지에 대해 설명할 수 있습니다.

내용 요약

이 연구는 세포유리 DNA의 메틸화 패턴을 분석하여 암을 진단하고 모니터링하는 방법을 개발하는 것을 목적으로 합니다. 암은 우리 주변에서 매우 흔하게 발생하는 질병 중 하나이며, 암을 조기에 발견하고 치료하는 것은 매우 중요합니다. 세포유리 DNA의 메틸화 패턴을 분석하는 기술은 암을 조기에 발견하고 치료하는 데 큰 도움을 줄 수 있습니다.

세포유리 DNA는 혈액 내에 존재하는 DNA 조각으로, 암세포에서 유래한 것일 수 있습니다. 이 연구에서는 세포유리 DNA의 메틸화 패턴을 분석하여 암세포에서 유래한 세포유리 DNA를 식별하는 방법을 개발하였습니다. 이를 위해 Nanopore sequencing 기술을 사용하여 세포유리 DNA의 염기서열을 분석하고, DNA methylation 분석을 수행하였습니다.

연구 결과, 세포유리 DNA의 메틸화 패턴은 암세포에서 유래한 세포유리 DNA와 건강한 세포에서 유래한 세포유리 DNA에서 차이가 있음을 확인하였습니다. 이는 세포유리 DNA의 메틸화 패턴을 분석하여 암을 진단할 수 있는 가능성을 보여줍니다. 또한, 이 연구는 세포유리 DNA의 메틸화 패턴이 시간에 따라 변화한다는 것을 발견하였습니다. 이는 암의 진행 상황을 모니터링하는 데 유용한 정보를 제공할 수 있습니다.

이 연구는 세포유리 DNA의 메틸화 패턴을 분석하여 암을 진단하고 모니터링하는 데 새로운 가능성을 제시한 것으로 평가됩니다. 이 연구 결과를 바탕으로, 세포유리 DNA의 메틸화 패턴을 분석하는 기술이 발전하면 암 진단과 치료에 더욱 효과적인 방법이 개발될 것으로 기대됩니다.

이 연구는 생명공학 분야에서 매우 중요한 연구 중 하나입니다. 세포유리 DNA의 메틸화 패턴을 분석하는 기술은 암 진단과 치료에 큰 도움을 줄 수 있으며, 미래에는 더욱 발전하여 다양한 질병을 진단하고 치료하는 데 활용될 수 있을 것입니다. 또한 미래의 건강과 의료 산업에 큰 영향을 미칠 것이며, 우리의 삶을 더욱 건강하고 안전하게 만들어줄 것입니다.

요약해줘 하니깐, 그래 이 연구는 중요해, 중요하니깐 중요해만 반복적으로 답변하는 느낌은 어쩔 수 없는듯합니다. 디테일 한 연구 내용을 확인하려면 질문을 계속 해서 확인해야 할 듯 합니다.

다시 잠깐 더 살펴보자면 단순히 ONT로 실험한 것이 아니고 기존의 표준실험 방법인 WGBS(Whole Genome Bisulfite Sequencing)와 비교한 데이터로 ONT로 연구해도 cfDNA의 메틸레이션을 확인 할 수 있다라고 주장하고 있는데... 생각해보면 cfDNA는 long read 시퀀서의 대명사인 ONT를 사용할 매력이 반감되기는 하나 native한 methylation을 볼 수 있다는 또 다른 압도적인 장점이 존재하다 보니..
그러나 보통 cfDNA가 150bp 이내의 짧은 것들만 생각하나 di나 tri-nucleosome도 함께 확인 할 수 있기 때문에... di나 tir nucleosome은 300~400bp 이상의 길이가 있기 때문에, 꼭 ONT의 매력이 반감되지는 않을 듯 합니다.
근데 di-, tri-nucleosome이 왜 중요한지는 제가 공부를 아직 안 해서 잘 모르겠습니다.
나중에 공부하게 되면 다시 글을 하나 써보는 걸로 :)


수요일, 3월 13, 2024

Long-read 시대의 정렬 알고리즘

사실 아직 Long-read 시대라고 하기에는 전세계 적으로 생산되는 시퀀싱 데이터의 비율로 보면 꼬꼬마이긴 하나 그래도 비용이나 지금까지의 축적된 자료 활용이라는 몇몇 점들을 제외하면 짧은 애들시대는 조만간 저물지 않을까 합니다.

물론 짧은 애들 또한 계속 필요는 합니다. 그러나 Long read 이전과는 같지 않을 것이기에 Long read 시대를 준비해서 나쁘지는 않을겁니다.

그래서 "A survey of mapping algorithms in the long-reads era" 라는 논문 하나 준비해봤습니다.

DIO: 10.1186/s13059-023-02972-3

시대가 시대이고 영어 읽기 귀찮으니 구글에게 한번 맡겨보았습니다.


Q: 내 드라이브에서"s13059-023-02972-3.pdf" 파일을 읽고 자세히 한글로 요약해줘

서론:

DNA 염기서열 분석은 생명과학 분야의 핵심 기술이며, 롱리드 시퀀싱 기술의 발전은 이 분야에 큰 변화를 가져왔습니다. 롱리드 기술은 기존의 짧은 리드 시퀀싱 기술에 비해 훨씬 더 긴 염기서열 정보를 얻을 수 있지만, 이에 따라 데이터 분석 과정에서도 새로운 문제들이 발생했습니다. 그 중 하나가 롱리드 데이터를 참고 유전체에 매핑하는 과정입니다. 

본 논문에서는 롱리드 시대의 DNA 염기서열 분석에 사용되는 다양한 매핑 알고리즘에 대한 조사 연구 결과를 제시합니다.


1. 롱리드 기술 개요:

롱리드 기술은 10kb 이상의 긴 염기서열 정보를 얻을 수 있는 기술입니다.

롱리드 기술의 장점은 다음과 같습니다.

  • 유전체 구조 변형 (SV) 분석에 유리
  • 단일 분자 haplotyping 가능
  • 염기서열 조립 (assembly) 과정 개선

롱리드 기술의 단점은 다음과 같습니다.

  • 높은 오류율
  • 높은 컴퓨팅 비용


2. 매핑 알고리즘 종류:

롱리드 매핑 알고리즘은 크게 씨앗-확장 (seed-and-extend) 프레임워크와 씨앗-체인 (seed-and-chain) 프레임워크로 나눌 수 있습니다.

씨앗-확장 프레임워크는 짧은 리드 매핑 알고리즘에서 사용되는 기본적인 방법입니다.

씨앗-체인 프레임워크는 롱리드 데이터의 특징을 고려하여 개발된 새로운 방법입니다.


3. 주요 매핑 알고리즘 비교 분석:

본 논문에서는 다음과 같은 주요 매핑 알고리즘을 비교 분석했습니다.

  • Minimap2
  • NGMLR
  • GraphMap
  • LSC
  • BWA-MEM
  • LAST
  • BLASR
  • RapMap
  • SMALT

비교 분석 기준은 다음과 같습니다.

  • 정확도
  • 속도
  • 메모리 사용량
  • 특징


4. 매핑 알고리즘 선택 가이드:

연구 목적에 따라 적절한 매핑 알고리즘을 선택해야 합니다.

정확도가 가장 중요한 경우 Minimap2 또는 NGMLR을 사용하는 것이 좋습니다.

속도가 가장 중요한 경우 LSC 또는 BWA-MEM을 사용하는 것이 좋습니다.

메모리 사용량이 가장 중요한 경우 GraphMap 또는 LAST를 사용하는 것이 좋습니다.


5. 향후 연구 방향:

  • 롱리드 기술 발전에 따라 새로운 매핑 알고리즘 개발
  • 멀티 샘플 데이터 분석을 위한 효율적인 알고리즘 개발
  • 염기서열 변형 (variant) 분석을 위한 정확한 알고리즘 개발


결론:

본 논문은 롱리드 시대의 DNA 염기서열 분석에 사용되는 다양한 매핑 알고리즘에 대한 종합적인 개요를 제공합니다. 각 알고리즘의 특징, 장단점, 성능 비교 등을 분석하여 연구 목적에 맞는 최적의 알고리즘 선택에 도움을 줍니다.

오... 비약도 있고 생략도 많지만 (물론 클로드 3에 비해서는 아직 갈길이 멀듯합니다)그래도 pdf파일을 요약해줍니다. 근데.. 추가 정보에 롱리드 시퀀싱 기술의 url이..  신박한 링크를 넣어놨네요;; Orz


출처: @ye._.vely618

목요일, 3월 07, 2024

나노포어용 정렬 프로그램을 비교해보자

오늘은 최근에 논란의 중심에 서 있는 claude3를 사용해보았습니다.

Comparative analysis of alignment tools for application on Nanopore sequencing data

이 논문을 선택한 이유는 claude.ai에서 무료로 해주는 양의 논문이라서 선택하였습니다.

 

Q. 요약해줘

이 논문에서는 나노포어 시퀀싱 데이터에 대한 6가지 다양한 정렬 도구(BLASR, GraphMap, LAST, minimap2, NGMLR 및 BWA MEM)의 성능과 정렬 품질을 실험 데이터셋을 사용하여 비교 분석했습니다. 주요 내용은 다음과 같습니다:

3개의 서로 다른 실험 나노포어 시퀀싱 데이터셋(람다 파지 전장 유전체, 인간 MEFV 유전자 암플리콘, SARS-CoV-2 전장 유전체)을 사용하여 정렬 도구를 비교했습니다.

각 도구의 CPU 시간, 메모리 사용량, 일치율, 불일치율, 오류율 등의 성능 및 품질 지표를 측정하고 비교했습니다.

Minimap2가 모든 데이터셋에서 속도와 메모리 효율성 면에서 가장 우수한 성능을 보였습니다.

BLASR, GraphMap, NGMLR는 상대적으로 느린 속도를 보였고, NGMLR는 높은 메모리 사용량을 보였습니다.

일치율, 불일치율, 오류율 측면에서 LAST를 제외한 도구들이 대체로 유사한 수준의 높은 정렬 품질을 보였습니다.

데이터셋 특성과 알고리즘 접근법에 따라 도구 간 차이가 있었으며, 실험 설계와 추가 분석 과정을 고려하여 적절한 도구를 선택해야 합니다.


Q. 그럼, 나노포어 왜 전용 정렬 도구가 필요합니까?

논문에 따르면 나노포어 시퀀싱 기술의 가장 큰 단점은 다른 시퀀싱 기술에 비해 높은 에러율(error rate)이라고 설명하고 있습니다. 구체적인 문장은 다음과 같습니다:

"So far, the most important disadvantage of this technology is the increased error rate compared to other sequencing technologies."

이처럼 나노포어 시퀀싱 데이터는 기존 단분자 시퀀싱 데이터와는 다른 에러 프로파일을 가지고 있어, 이러한 높은 에러율을 보완하기 위해 전용의 정렬 도구가 필요하다고 설명하고 있습니다.


Q. 좀더 자세히 설명해줘

네, 더 자세히 설명드리겠습니다.

나노포어 시퀀싱은 전기화학적 신호를 읽어 DNA/RNA 서열을 분석하는 새로운 기술입니다. 이 기술의 장점은 긴 읽기(long reads)가 가능하고 빠르며 비용이 저렴하다는 점입니다. 하지만 단점으로 다른 기존 시퀀싱 기술에 비해 상대적으로 높은 에러율을 가지고 있습니다.

이런 높은 에러율은 나노포어 시퀀싱 데이터의 고유한 에러 프로파일에서 비롯됩니다. 삽입(insertion), 삭제(deletion), 치환(substitution) 등의 에러가 많이 발생합니다. 이러한 에러들을 효과적으로 다루지 않으면 후속 데이터 분석 단계에서 오류가 커질 수 있습니다.

따라서 나노포어 시퀀싱 데이터를 정확하게 분석하기 위해서는 이러한 높은 에러율과 특수한 에러 프로파일을 잘 보완할 수 있는 전용 정렬(alignment) 도구가 필요합니다. 기존의 단분자 시퀀싱 데이터 분석용 도구들은 나노포어 데이터의 특성을 제대로 반영하지 못해 정확성이 낮을 수 있기 때문입니다.

요약하면, 나노포어 시퀀싱의 높은 에러율과 특수한 에러 프로파일을 보완하여 정확한 데이터 분석이 가능하도록 최적화된 전용 정렬 도구가 필요한 것입니다.


음.. 일단 무료 요금에서 사용하는 양인 경우 일단 쓸만한데 대부분의 논문 양은 유료버전 사용하시라고 뜰껍니다. :)



출처: @ye._.vely618


화요일, 2월 13, 2024

생명정보를 끄집어 내는 나노포어!!

음.. 21년도에 옥스포드 나노포어에 관련하여 잘 설명한 글이 있었네요

"Beyond sequencing: machine learning algorithms extract biology hidden in Nanopore signal data"

DOI: 10.1016/j.tig.2021.09.001


그래서 당연히 이것도 gemini에게 블로그에 게재할 요약 글 하나 적어달라고 해봤습니다.


제목: 나노포어 시퀀싱, 숨겨진 생명 정보를 끄집어내는 미래 기술

나노포어에 대해 잘 모르는 분들을 위한 블로그 글을 쓰고 싶다고 하셨네요. 좋은 선택이에요! 이 글에서는 나노포어 시퀀싱 기술이 어떻게 작동하고, 어떤 흥미로운 활용 방식이 있는지 알아볼 예정입니다.


나노포어 시퀀싱의 기본 원리

먼저 나노포어 시퀀싱이 무엇인지 알아볼까요? 간단히 말하면 DNA나 RNA를 작은 구멍(나노포어)을 통과시켜 염기 서열 정보를 읽는 기술입니다. DNA/RNA 가닥이 나노포어를 통과하며 발생하는 전류 변화를 분석하여 염기 서열을 추측합니다. 하지만 이 전류 변화 패턴은 복잡하고 해석하기 어렵죠.


머신러닝과 나노포어의 만남

여기서 머신러닝이 중요한 역할을 합니다! 머신러닝 알고리즘은 대量の 나노포어 전류 데이터를 분석하여 이 패턴을 알파벳 코드(A, C, T, G)로 변환하는 방법을 배우게 됩니다. 이를 베이스 콜링(base calling)이라고 합니다. 

초기에는 단순한 통계적 방법과 규칙 기반 알고리즘이 사용되었지만, 최근에는 머신러닝 기반 알고리즘이 주목받고 있습니다. 이 알고리즘은 실험 데이터를 대량 학습하여 전류 패턴과 염기 사이의 복잡한 관계를 파악합니다. 대표적인 알고리즘으로는 히든 마르코프 모델 (Hidden Markov Model, HMM)과 신경망 (Neural Network)이 있습니다.

 - HMM: 각 염기마다 고유한 상태를 가정하고, 전류 패턴의 변화를 통해 상태 이동을 예측합니다. 이를 통해 전체 염기 서열을 추측합니다.

 - 신경망: 복잡한 전류 패턴을 직접 분석하여 염기 예측 확률을 계산합니다. 최근에는 딥러닝 기반의 신경망 모델들이 우수한 성능을 보여주고 있습니다.

이러한 알고리즘은 지속적으로 발전하고 있으며, 베이스 콜링 정확도는 점점 향상되고 있습니다. 하지만 아직 완벽하지 않으며, 특히 동일한 전류 패턴을 유발하는 modified base 구분은 어려움이 있습니다. 연구는 지속적으로 진행되고 있으며, 더 정확하고 신뢰할 수 있는 베이스 콜링 알고리즘 개발이 활발히 이루어지고 있습니다.

또한 베이스 콜링을 통해 염기 서열을 얻었지만, 이것만으로는 충분하지 않습니다. 마치 흩어진 퍼즐 조각처럼 이 염기 서열을 참고 게놈 전체에 위치를 정확하게 매핑해야 합니다. 이 과정을 시퀀스 매핑(sequence mapping)이라고 합니다.

시퀀스 매핑 알고리즘은 크게 두 가지 방법이 있습니다.

 - 참조 게놈 기반 매핑: 이미 알려진 참조 게놈을 사용하여 읽은 염기 서열과 유사한 부분을 찾아 매핑합니다. 빠르고 정확하지만 돌연변이가 있는 경우 정확도가 떨어질 수 있습니다.

 - De novo 매핑: 참조 게놈 없이 읽은 염기 서열 자체를 조합하여 전체 게놈을 재구성합니다. 돌연변이가 있어도 정확하게 매핑할 수 있지만, 계산량이 많고 어렵습니다.

최근에는 이 두 방법을 결합하여 정확도와 속도를 높이는 하이브리드 방식도 사용되고 있습니다. 또한 롱 리드 시퀀싱 기술인 나노포어는 다른 시퀀싱 기술에 비해 더 긴 염기 서열을 읽을 수 있어 매핑 과정의 정확도를 높이는 데 유리합니다. 머신러닝 알고리즘은 이 과정에도 활용되어, 정확한 매핑을 위해 지속적으로 발전하고 있습니다.


나노포어 활용의 무한 가능성

이제 흥미로운 부분, 나노포어를 활용하여 무엇을 할 수 있는지 알아볼까요? 기존 시퀀싱 기술에 비해 장점이 많은 나노포어는 다양한 분야에서 혁신을 가져올 것으로 기대됩니다.

 - DNA 메틸화 분석: DNA에 달린 메틸 그룹은 유전자 발현을 조절합니다. 나노포어는 이 메틸화 정보를 전류 패턴에서 직접 읽을 수 있어 유전자 발현 연구에 유용합니다.

 - 단백질-DNA 상호 작용 분석: DNA와 상호 작용하는 단백질을 나노포어를 통해 직접 확인할 수 있습니다. 이는 질병 유발 기작 연구 등에 활용될 수 있습니다.

 - RNA 구조 분석: RNA는 독특한 구조를 가지고 있으며, 나노포어는 이 구조 정보도 전류 패턴에서 추측할 수 있습니다. 이는 RNA 기능 연구에 중요한 정보를 제공합니다.

이 외에도 나노포어는 미생물 동정, 식물 연구, 개인 맞춤 의료 등 다양한 분야에서 활용 가능성이 높습니다.


결론

나노포어 시퀀싱은 빠르게 발전하는 기술이며, 머신러닝과의 결합은 더욱 강력한 생명 과학 연구 도구로서 활약할 것입니다. 이 글을 통해 나노포어에 대한 기본적인 이해를 얻고, 앞으로의 발전 가능성에 대해 기대를 갖게 되었기를 바랍니다!


참고

* 이 글은 생명공학 분야 기초 지식을 가진 독자를 대상으로 이해하기 쉽게 설명했습니다.

* 논문 원본은 영어이지만, 주요 내용을 한국어로 요약했습니다.

* 더 자세한 정보는 제공된 논문 링크를 참조하세요.

저도 나노포어 기술의 발전을 지켜볼 예정이니, 앞으로 더 흥미로운 소식을 공유하도록 하겠습니다!


응? 왜 니가 나노포어 기술의 발전을 지켜볼 예정이니? ㄴ(°0°)ㄱ 내가 블로그에 작성하는 글이라고 하니 사람이 쓰는 것처럼 소식을 공유하겠다는 멘트도 깨알같이 놓치지 않았네요.   

그럼 다음에도 어딘가에 제가 아직 보지 못한 논문들을 찾아서 AI를 활용하여 글을 또 써보도록 하겠습니다. 물론 AI가 글을 쓰겠지요 :)



출처: @ye._.vely618


일요일, 11월 12, 2023

Long Read Sequencing을 전적으로 믿으셔야 합니다.

항상 느끼는것이지만 사람들은 익숙한것에 익숙하다는...
물론 이 글을 쓰는 본인도 별반 차이 없다는게 현실
그럼에도 불구하고 세상에는 익숙한것을 거부하고 한 걸음 나아가는 분들이 있어서 발전한다는...

오늘은 유전체분야에서 다들 익숙한 짧은길이의 시퀀싱 플랫폼 대신 Long read 시퀀싱 플랫폼을 왜 써야하는지 보여주는 논문이 있어 들고와봤습니다.

The blooming of long-read sequencing reforms biomedical research

DOI: 10.1186/s13059-022-02604-2


Long read 시퀀싱의 대표주자인 PacBio와 Nanopore가 나온지도 꽤 된것 같으나, 아직 주류 시퀀싱 플랫폼으로서는 자리매김을 하지 못한... 조명받지 못하고 있지만.. 그 진가를 하나둘씩 알게되고 안쓸수 없지 않을까 합니다.

물론 이 Long read 플랫폼을 안 쓸수 없지만 PacBio/Nanopore 플랫폼을 이용해서 진단 키트를 개발하여 돈을 벌 수 있게 된다는 것은 아닙니다. 오해 없으시기 바랍니다. Long read 플랫폼에서 기존 NGS 시스템에서 해온 것 같이 하려고 하면 개인적인 생각으로는 그냥 하지 마세요 라고 말하고 싶네요. 그런 생각이라면..


여튼 Long read 플랫폼은 genome assembly과 transcriptome 분야에서 활약을 하고 있다고 하는데 이는 기존 NGS의 짧은 read로 인한 한계를 극복하였기 때문에 가능한 당연한 결과였을 듯 합니다.

genome assembly를 설명하면서 이런 저런 설명을 하였는데 사실 T2T genome이 세상에 나왔는데 무슨 설명이 더 필요할지... 물론 다배체면서, repeat 서열이 엄청 긴 식물과 같은 다양한 생물의 유전체 연구에서는 아직 할일이 많이 남아있을듯 합니다.

그리고 더불어 transcriptome, 전사체에서도 두각을 나타내고 있다고 합니다. 이전에 PcaBio사에서 나온 Iso-Seq이라는 플랫폼이 있었는데 기존 숏리드 플랫폼으로는 할 수 없었던 full-length 유전자 서열을 확인하여, 다양한 gene의 isoform 을 확인 할 수 있었는데, 이제는 이는 당연한것이고, read count를 활용하여 정확한 발현 측정까지도.. 가능하다고 합니다. 또한 이 논문에서 처음 보았는데 exitron이라는 것도 확인하였다고 하네요.

더불어 당연히 암연구에서도 long-read 플랫폼이 중요한 역할을 하고 있는데, 시퀀싱을 위해 dna나 rna가닥을 증폭시키지 않은 특징 때문에 가능 응용법이긴하죠. 특히 nanopre 플랫폼을 사용하여 급성 골수성 백혈병 (AML) 환자의 RNA변이와 fusion gene 구조를 확인하기도 하였고, 대장암 세포주에서 5mC 위치와 양을 정량하여 각각 AML와 대장암 연구에 도움이 되고 있다고 합니다.

이렇듯 Long Read Seq 플랫폼은 이제 쓰지 않을 이유가 없는 플랫폼으로 우리 곁에 생각보다 가깝게 다가왔으나 생각에는 Long이든 Short든 정확하게 내 병의 이유를 빠르고 "싸게" 알아낼 수 있으면 모두 행복한 것 아니겠습니까?

그럴 수 있는 방법이 어딘가에 있겠죠, 없으면, 언제나 그러했듯이 우리는 또 해답을 찾아내겠죠. :)




출처: @ye._.vely618


월요일, 2월 06, 2023

우리들의 Nanopore는 결핵균을 잘 시퀀싱하여 세계 공중 보건에 기여할 수 있을까?

Evaluation of Nanopore sequencing for Mycobacterium tuberculosis drug susceptibility testing and outbreak investigation: a genomic analysis

https://doi.org/10.1016/S2666-5247(22)00301-9 


그래서 한번 알아 보았습니다.

일루미나 일색인 시퀀싱 시장에 신성처럼 찾아온 (찾아왔지 떨어지지 않았습니다.) long-read sequencing의 이단아 ONT을 사용한 결핵균의 drug susceptibility testing (DST) 결과가 기존 일루미나 결과에 비해 부족하지 않고, 동등한 결과와 정보를 제공해 줄 수 있는지 확인해 보았습니다.


영국, 남아프리카등의 샘플을 MiSeq, HiSeq, NextSeq, MinION, GridION, PacBio 기기를 사용하여 시퀀싱 진행하였습니다. 물론 개별 시퀀서마다 동일한 샘플을 시퀀싱한 것은 아니고 영국샘플은 MiSeq, GridION에서, 남아프리카나 마다가스카르 샘플은 HiSeq, NextSeq, MinION에서 진행하였으나, 전체적으로 short와 long read 시퀀서간에 비교 할 수 있도록 디자인은 되었습니다.

일루미나는 UK에서 사용하는 COMPASS를 사용해서, Nanopore는 BCFtools를 사용하여 변이를 확인하였습니다.

그래서 일루미나와 Nanopore의 변이 결과를 확인해보니(본문의 figure1을 참고) 재현율이 다소 떨어지는 경향을 확인 할 수 있었습니다. 

Nanopore는 정밀도: 99.3%, 재현율: 90.2%

일루미나는 정밀도: 99.6%, 재현율: 91.9%

그리고 151개의 isolate(시퀀싱한 isolate는 208개이나 일루미나와 Nanopore를 동시에 비교 할 수 있는 데이터는 151개 였음)에서 66,537개의 저항성 변이가 call되었는데 일루미나와 Nanopore간에 4개의 차이밖에 보이지 않았다고 합니다. (이정도면 거의 동일한 수준 아닌가?)

그래서 결론은 우리 모두 예상하듯이 Nanopore 데이터로도 정밀한 변이 call이 가능함을 확인하였고, 그 결과 또한 일루미나와 동등한 수준이므로, 깨작 깨작 하지 말고 Nanopore를 이용해서 MTB 전체 서열을 한번에 분석해서 결과를 주면 좋지 않을까 기대하고 싶다. 되지 안하을까합니다.


본 저자는 ONT.L 주식을 (아직) 보유하고 있지 않습니다.










수요일, 7월 20, 2022

Long Read 조립은 누가누가 잘하나

Piroplasm를 나노포어를 사용하여 genome project를 진행했고 나노포어를 활용한 assembler들에 대한 성능 비교 논문 되겠습니다.


제목은 Systematic Comparison of the Performances of De Novo Genome Assemblers for Oxford Nanopore Technology Reads From Piroplasm

doi는 https://doi.org/10.3389/fcimb.2021.696669


piroplasm이 몬지는 모르겠으나 일단 그렇게 엄청나게 복잡하지는 않은 원생동물이나 사람이나 동물들에게 질병을 일으키는 녀석 되는 것 같습니다. funding중에 동물 전염병 및 인수공통전염병 관련 프로그램이 있는것으로 보아하니...


여튼 중요한건 nanopore로 읽어낸 서열을 사용하여 genome 조립할때 어떤 어떤 조립 프로그램이 제일 좋은지를 검토해본 것이니 OLC (Over-Layout-Consensus)나 전통적인 de-Brujin graph, string graph-based 방법 등등의 NECAT, Canu, wtdbg2, Miniasm, Smartdenovo, FlyeNextDenovo, Shasta와 같이 일반적으로 long-read에 사용하는 assembler들을 비교 테스트 하였다고 합니다.

대신 여기서는 assembly의 정확도와 함께 CPU 사용량, 메모리 사용량, 분석시간 사용방법 등등에 대해서도 함께 평가했다고 합니다. 참 바람직한 태도라고 봅니다. 모든 연구팀들이 그래픽카드 4개꼽히고 6T 메모리의 4U 서버를 가지고있는것은 아니니 말입니다.


실험 방법은 prioplasm free한 양 2마리(??)에게 prioplasm을 감염시켜 잘 배양(??)시킨 다음 Qiagen 사용 prep kit을 가지고 DNA 추출하고 PromethION으로 시퀀싱하였고 데이터 셋트에 따른 assembly 결과 평가를 위해 6가지 생산량 (약 15x, 30x, 50x, 70x, 100x, 120x)의 셋트를 만들었다고 합니다. 그리고 추가적(aka error correction)으로 (일루미나와 특허 소송에서 승리한) MGI로도 시퀀싱을 하였다고 합니다.


여튼 결과적으로

N50과 contig개수(적을수록 좋음)는 생산량과 밀접하고,
분석 시간은 생산량이 많으면 어떤 assembler를 사용하던 길어졌고,
polishing은 안하는것보다 하는것이 좋은것 같고 각 tool의 장단점은 Figure3에 방사형 그래프로 이쁘게 표현하였으니 한번 참고하시면 좋을것 같습니다.

그래서 Miniasm, Flye, wtdbg2는 그닥 좋은 선택지는 아닌것 같고 평균 커버리지가 30x 이상 확보된다면 NECAT, Canu, NextDenovo, Smartdenovo가 더 나은것 같다 정도 되겠습니다.

(사실 위의 tool들을 실행시키려면 평균 30x 이상은 있어야 작동을 합니다. 안그러면 작동안하던지 말도안되는 결과들을 뱉어내곤 합니다.)


그리고 시간이 충분했는지 각 assembler 결과들을 병합/후처리하는 작업을 하여 더 나은 assembly 결과를 보여주는지 테스트 했고 몇몇 조합에서 결과물이 향상된것을 확인했다는데... dramatically 좋은 결과는 보여주지 않은것 같았습니다. 

만약 병합/후처리하는 결과가 좋았다면 논문 결과가 single assembler 쓰지말고 ensemble방법을 추천드립니다라고 했었을테니 말이죠.. 



출처: @candyz_hyojung


일요일, 11월 07, 2021

HLA를 빠르고 정확하게 할 수 없을까

그래서 나노포어가 준비했습니다.

절대적으로 빠른 건지, 정확도의 수준은 정확히 모르겠으나 일단 샘플당 $80되겠습니다.
물론 이게 재료비값인지 소비자가인지는.....
아마 본인들이 했을때 가격인듯.... (일반 소비자가 이 가격에 가능할지는 모르겠네요)

Rapid, highly accurate and cost-effective open-source simultaneous complete HLA typing and phasing of class I and II alleles using nanopore sequencing

그럼 HLA(Human leukocyte antigen)검사는 무엇이냐? 쉽게 얘기하면 장기이식할때 장기이식할수 있는지 여부를 확인하는 검사되겠습니다. 이 HLA가 내거와 내거 아닌거를 구분하게 해주는 녀석되겠습니다. 

면역학 책 들춰보시면... 책을 덮습니다.

일반적으로 HLA는 A,B,C인 Class I, DRA1, DRB1/3/4/5, DQA1, DQB1, DPA1, DPB1인 Class II인데 우리나라는 A,B,C를 포함해서 일부만 본다고 합니다. 

여튼 나노포어는 긴 길이로 시퀀싱하는 장점과 빠른 시퀀싱속도를 바탕으로 시퀀싱 시간은 8시간내 (물론 라이브러리 제작시간은... 따로 좀 봐야할것 같습니다.) 분석은 샘플당 30분(근데 사용한 서버 수준이... 32core 256g 메모리 서버)소요 됐다고하는데...

1 샘플(aka 1검체)의 경우 나노포어 한개 사용해서 기존 일루미나나 써모 플랫폼보다 빠르고 싸게 결과를 도출할 수 있도록 해준다는데...

영어의 짧음인지...  가끔 나노포어가 더 안좋아 보이게 설명하는거 같기도합니다. 


결론은 나노포어가 더 길게 시퀀싱하기 떄문에 복잡한 HLA에서 더 정확하게 구별할 수 있고,
기존에 short 시퀀싱에 적합한 pcr primer말고 long range PCR에 적합하게 프라이머를 변경해서 진행하면 HLA의 복잡한 부분들을 더 잘 잡아낼 수 있을것이라고 털고있습니다.


여튼 개인적으로 기존 short 시퀀서 플랫폼 보다는 나노포어가 HLA에서는 좀더 나은 성능을 보여주지 않을까합니다. HLA의 복잡도 때문에 나노포어의 장점을 활용할 수 있는 PCR 프라이머로 재 디자인해서 진행하면 시퀀싱 시간/분석 시간/정확도에서 가장 적합하기 않을까하는 생각입니다.




@ye._.vely618


화요일, 9월 14, 2021

Long-read sequencing identified a causal structural variant in an exome-negative case and enabled preimplantation genetic diagnosis

아니.. 2018년도에 이런 논문을 블로그에 keep놓고 그냥 놀고 있었다니..
여하튼 눈에 띄어서 한편 사브작 하나 올려봅니다.
역시 일이 많을때 딴짓은 국룰이라는 명제를 증명하듯......


WGS(Whole Genome Sequencing)이 아무리 저렴해졌다고하나 내 genome에 이상이 있는지 진단을 위해서 시퀀싱하는 비용이 억대에서 천만원 정도로 떨어진거지 아직은 WGS가지고 일반인이 무엇을 하기에는.... 물론 Private Premium Platinum Gold Special Society 같은 수식어의 멤버 라면... 할려면 할수는 있겠다만..

여튼 현재 아직까지는 유전체 진단에서 가장 현실적이고 합리적인 방법은 WES이나.. 짧은 read들을 가지고는 structural variant를 탐지하기에는 무리가 좀 있으니 long read로 해야하지 않겠냐라고 얘기하는 2018년도 논문되겠습니다.

해당 논문에서 언급하는 환자는 G6PC 유전자에 돌연변이가 생겨 간비대, 신장비대를 초래하는 recessive inheritance 질환으로 긴 지역에 deletion이 있었으나 WES만으로 진단하였을 때는 확인하지 못했고 nanopore (long read seq)를 사용하여 확인하였을 때에는 정확히 긴 길이의 SV를 확인 할 수 있어다 라는게 주제되겠습니다. 지금까지 short read만으로 시퀀싱하여 진단되지 않거나 오진되는 환자들의 돌연변이를 정확하게 확인 할 수 있는 도구가 임상 진단에서 하나더 생겼고 그로 인해 정확한 진단을 할 수 있게 되었다는 자화자찬의 평을 내면서 논문을 마무리하였습니다.

사례를 통해 short read만으로는 진단에 한계가 있다는것을 확인하였고 long read 방법이 이제 진단에 슬슬 자리 잡아야 하지 않나 싶은 생각이 드는....

(근데 길게 뽑으려면 나노포어 실험이 그렇게 쉽지 않은듯하던데.... )





출처: @ye._.vely618
출처: @ye._.vely618


토요일, 1월 02, 2021

High Heterozygosity genome 어셈블리 할 때 해결사로 자처하고 나온 어셈블러

2021년 첫 포스팅 시작합니다. :)

오늘은 de novo assembly관련된 tools 소개 글입니다.

모델 생물 혹은 그외의 생명체에 대해서 genome을 알고 싶다면...
그냥 절래절래 하거나 아니면 필요하거나 알고 싶은 부분만 클로닝해서 슬쩍 슬쩍 알아내던 시절... 

돈없는 연구자들은 손가락이나 빨고 눈물이나 흘리던지 어딘가에 끼어서 연구를 진행했었어야 했으나...

이제는 바야흐로 2021년!! (사실 글의 초안을 작성하고 있었던 때는 2020년 11월;;;)

454 따위 역사속으로 사라지고 킹왕짱 long read인 PacBio와 나노포어가 활발히 사용되는 시대에 살고 있습니다.

(물론 de novo는 454와 일루미나로 이어 붙여야 제맛이지 하는 라떼들도 있지만...)

저는 박테리아에서 PacBio 써본 후 박테리아 연구에서 PacBio 사용하지 않고 일루미나 플랫폼을 사용한다고하면 일루미나 주식을 가지고 있나 생각하게 되었고, 그외에 genome을 de novo작업은 nanopore를 추천하고 있는데 굳이 나노포어를 사용하지 않겠다라고 한다면... 굳이? 라는 궁금증으로 가지게 되었다능.. ㅋㅋ
(당연히 무조건 저걸 써야하는건 아니쥬 ㅎㅎ 연구 목적에 따라 봐야 하는 결과물에 따라 플랫폼을 잘 선택하시면되겠습니다.)

※ 여기서 나노포어 base의 quality는 굳이 논할 이유는 없고 제대로된 데이터 만져보면 나노포어 쓰지 않을 이유가 없으실 겁니다.


이제는 PacBio나 나노포어를 굳이 사용하지 않을 이유가 없는 것이 예전이라면 PacBio나 나노포어를 지원해주는 프로그램이 많지 않아서 좀 꺼려질수 있었지만.. 지금은 반대로 너무 많아셔져서 어떤 tools을 사용해야 하나 할 정도니...

그중에서 어셈블리를 하면서 문제가 되는 부분이 High heterozygous region들이 있는데 (그 와 함께 저세상 텐션을 보여주는 polyploidy;; ) 이런 문제들을 해결에 주겠다는 해결사로 자처하고 나온 tools이 있어서 한번 끄적여 보았습니다.

Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies

물론 이전에도 heterozygous 문제를 해결하는 tools이 없지는 않았습니다.
이전 글에서도 잠시 소개했었던 HaploMerger2 도 있고 저는 잘 몰랐지만 Redundans라는 도구도 있었다고 합니다. 다만 이전에 나온 tools의 단점은 사용자가 수동으로 contig를 재 할당해야한다는 문제가 있다고 합니다.
(음.. 저도 이전에 한두번 HaploMerge2를 사용해었는데;;; 여기서 얘기하는 contig 재할당에 대한 얘기가 정확히 어떤 의미인지는 정확하게 모르겠네요;; 여튼... 좀 단계 단계를 수동으로 작업을 하기는 했었습니다만 여기서 얘기하는 "수동"이 이 얘기가 아닌거 같은데...)

여튼... purge의 분석 pipeline은 다음과 같은 단계들로 진행됩니다.
purge의 분석 Flow chart >Figure1<


음... 확실히 장점으로는 draft로 조립된 genome에서 중복되는 contig들을 제거해서 draft assembly 서열의 크기를 줄여서 실제 genome size에 가깝게 된다는 것 이긴 합니다.

그리고 이 tool을 사용할 시 참고할 점으로는 일루미나 데이터를 가지고 purge를 진행할때는 지양했으면 합니다.
제가 해봤을때에 nanopore-raw 서열가지고는 분석이 가능했는데 일루미나 데이터로 작업하였을 때 purge_haplotigs의 hist 명령어를 사용하여 cov의 input 값이 -l, -m, -h 값을 구하는 작업을 할때 -l, -m, -h 값을 특정 할 수 있는 문제가 좀 있었습니다. 너무 값들이 낮아서 어떤 값을 low, mid, high을 사용할지가 모호하더라구요..

여튼..

de novo  작업을 하시다가 생각보다 genome 크기가 큰 경우 내 genome 중간에 heterozygous한 지역이 있구나 생각하시고 이 tools한번 돌려보시면 좋을것 같습니다. :)



출처: @ye._.vely618





토요일, 12월 26, 2020

Long read는 Long read alignment로....

간만에 글 투척합니다.

오늘은 Long read align관련된 내용 투척 하도록 하겠습니다.

시퀀싱 기술이 좋아져서 PacBio와 함께 nanopore가 (최근 covid도 있었고) 함께 많이 이용되고 있는데 질좋은 long read를 생산하셨다면 bwa/bowtie와 같은 short read alignment보다는 long read에 최적화된 long read alignment인 minimap2같은 tool를 사용하시는게 정신건강에 좋을 것 같습니다.

최근 제가 별생각없이 long read서열을 human genome에 aling할때 bwa mem의 -x ont2d옵션을 사용하여 진행했었는데...

통수를 후려 갈기는... (원래 통수는 후려 갈기는 맛이 좋다능)


bwa mem -x ont2d를 사용하는 경우 원래 생산된 read(raw read 개수)보다 더 많은 read 들이 aligned되는 것 처럼 보인다는... (?? 생산된게 100개 read인데 align된 read는 200개라고?)

근데 minimap2를 사용하였을 때에는 reference에 align된 read의 개수와 생산된 read 개수(raw 리드 개수)가 유사한(aka 생산된 read 개수보다 적은) read 개수를 확인 할 수 있었습니다.

개발자도 같고(Heng Li), 같은 align하는 tool인데 왜 이렇게 차이가 나느냐 

음... bwa와 minimap2를 자세히 까보지는 않았지만(못하지만) bwa의 경우 mem -x ont2d를 사용한다고 하더라도 태생부터가 short read를 위해서 만들어진 alginment이고 nanopore의 경우 마음먹고 시퀀싱이 된다면 수십 kb의 길이가 나오는 관계로 bwa mem에서 seed를 기준으로  align할때 mismatch나 error에 대해서 관대하게 조건을 잡으면서 확장을 한다로 하더라도 수십 kb까지 확장하지 못하는 경우가 발생 하고 그 경우 split이 되어 다른 reference에 align되는 경우가 발생하는것으로 보였습니다.

대신 minimap2의 경우 long read를 고려해서 만들다 보니 확장이 비상식적으로 read가 길더라도 확장을 하지 split하지 않아서 bwa의 경우와 같이 생산된 read보다 많은 read가 align된 것 처럼 보이는 이상한 문제는 발생하지 않는 것 처럼 보였습니다.

결론은 내 데이터에 맞는 mapper를 사용해서 분석하자 되겠습니다. ㅠ.ㅜ


추신: 2018년도 이런 글(Minimap2 and the future of BWA)도 있었군요;; 



출처: @ye._.vely618


일요일, 4월 26, 2020

Benchmarking of long-read assemblers for prokaryote whole genome sequencing

나노포어는 현존하는 시퀀서중에 가장 긴 서열을 뽑아내는 시퀀서임에는 그 누구도 부인하지 못할것입니다. 근데.. 생산된 리드의 각 base의 phred score를 보자면.. 왜 갑자기 눈에서 물이나오는 이유는 왜때문일까요?
(그렇지만 저는 de-novo할때 보수적인 그룹이 아니라면 나노포어를 권장하는건 비밀..)

여하튼.. 현재 나노포어 어셈블리 용으로 이런저런 어셈블러가 판치고 있는 난세에 누가누가 좋은지 확인하는 작업을 해서 투고하신분이 나타나셨습니다.
제목도 정직합니다. 단, prokaryote대상입니다.
Benchmarking of long-read assemblers for prokaryote whole genome sequencing

prokaryote에서도 개판이면 굳이 사용할 이유가 있겠느냐? 주의 되겠습니다.
일단 가장 좋은것은 모르겠지만 최악은 걸러내야 해야 시간 낭비, 전기 낭비 하지 않지 않겠습니까?

현재(aka 당시에) 돌려볼 수 있는 7개 어셈블러 (Canu, Flye, Miniasm/Minipolish, NECAT, Raven, Redbean, Shasta)의 성능을 비교 평가 했습니다.
어셈블리의 정확성은 당연하고, prokarypte다 보니 circularisation도 중요하고, 계산시 사용되는 리소스와 분석 시간등을 평가했다고 합니다.

아름다운 figure는 상단에 링크된 논문에서 감상하시면 되고,
canu는 그나마 볼만한 서열들을 제공해줬고
flye는 canu다음으로 괜찮은 서열로 어셈블리 했다고 합니다.
redbean(wtdbg2) 과 shasta는 계산 리소스와 분석 시간에서는 효율적이었지만 결과는 그다지 효율적이지 않았고 하네요.

그래서 종합해서 논문에서 결론을 냈는데
모.. de-novo aseeembly 해보신분이라면 알고계시다 싶이.. 다들 장단점이 있었고, 원탑인 어셈블러는 없었지만 그 중에서 Flye, Miniasm / Minipolish와 raven이지 않나 싶다고 하네요

Flye는 믿을만한 서열을 제공했고(low depth에서도 나름..)
Miniasm / Minipolish는 circularisation이 좋았고
raven은 identity가 낮은 read set들에서 tolerant가 있었다고 합니다.

역시 최적의 어셈블리를 위한 정도는 당신이 사용 가능한 리소스를 동원해서 다양하게 돌려보고 비교한 게 킹왕짱이지 남의말 믿고 쓰면 너만 바보 되고
개발자님들인 이런 상황이니 개발좀 굽신굽신 :)



출처: @sana_twice.09

금요일, 3월 31, 2017

nanonet 사용법


Nanonet Git Hub

Nanopore MinION raw 데이터인 fast5에서 2D read 추출하려고 하는데
nanopore로 검색하면 나오는 poretools, nanopolish, nanook들을 사용했음에도
어느하나 2D 서열을 추출 못하고 있어서 전전긍긍하고 있을때

혜식처럼 나타나는 그분이 알려주신 바로 그 tool!!

nanonet되시겠습니다.

우선 CPU기반으로 사용하였습니다. nanonet은 GPU도 지원하고 있는데
GPU로 하였을때 CPU기반으로 돌렸을 때 확인되던 2D가 나오지 않았다능..

일단 어쩔수 없지만 cpu기반으로 수행

fastq로 저장을 해보자!
$nanonet2d --fastq --jobs 8 --chemistry r9.4 pass output.r94

fasta로 저장을 해보자!
$nanonet2d --jobs 8 --chemistry r9.4 pass output.r94

다른 tools들과 다른 점은 chemistry버전을 선택할 수 있는 것과 무엇인지 모르겠지만 간지나 보이는 trmplate와 complement의 model 선택 옵션
그리고 이를 뒷받침하듯이 nanonettrain명령어도 있고....
poretools같은 일반적인 tools이 안먹힐때는 nanonet도 한번 사용해 보시기 바랍니다. ;)

화요일, 3월 14, 2017

poretools 사용법

poretools doc사이트

poretools git url

Oxford Nanopore 데이터 핸들링을 위한 tools중 하나로
현재 다양한 nano로 시작하는 tools들이 있긴 하지만
우선 이 tools을 먼저 살펴보았습니다.

이유는 간단합니다. 제가 사용하고 있으니 작업 로그 삼아 글쓰는거죠 ㅎㅎ

설치는 간단합니다.
먼저 poretools에 필요한 hdf5, numpy, cython 라이브러리들을 설치해 주시고
다음과 같이 샤샤샥

git clone https://github.com/arq5x/poretools
cd poretools
python setup.py install
물론 이렇게 해서 한큐에 되는 경우는 그렇게 많지는 않죠..
물론 위는 root권한이 있어야 가능하고요
저같은 일반 유저들은 다음과 같이 해야겟죠(물론 이런 저런 라이브러리는 알아서 잘 설치를...)

git clone https://github.com/arq5x/poretools
cd poretools
python setup.py install --home=/path/to/install 
그럼 일단 큰 문제 없이 설치는 될 겁니다.

poretools는 nanopore read 추출/통계정도의 작업제외하고는 할 수 있는게...
그래서 여기서는 extract read관련한것만 간단히 적고자 합니다.

nanopore데이터는 hdf5로 PacBio와 format이 유사해서 직접 데이터를 확인 할 수 없습니다.
그래서 fasta파일이나 fastq파일로 추출 작업을 거쳐야 합니다.

다행히 추출 방법은 그리어렵지 않습니다(사실 어려울게 없죠).

poretools {fasta|fastq} /path/to/reads/pass/ > reads.{fa|fq}

이렇게 하면 fast5파일에서 fasta나 fastq정보를 추출해서 파일로 만들어 줍니다.
근데 nanopore보시다 보면 1D다 2D다 이상한 얘기를 하는것을 볼 수 있습니다.
1D, 2D는 sequence type입니다. 쉽게 얘기해서 1D는 한번 읽고 2D는 두번 읽은거

그래서 2D가 먼데?  요기

그러면 2D 서열은 어떻게 추출하는건가?

다음과 같이 각 type에 따라 저장 할 수 있습니다.
poretools {fasta|fastq} --type {all|fwd|rev|2D|fwd,rev|best} /path/to/reads/pass/ > reads.{fa|fq}

수요일, 2월 15, 2017

Oxford Nanopore 쉽게 이해하기


이제 옥스포드 나노포어가 나온지 좀 됐죠?

롱롱롱 read도 시퀀싱(아직 accuracy문제는 있긴 하지만..) 하고
시퀀싱 후에 듣보잡 서버에서만 쿵짝쿵짝했어야
수십기가 txt파일 받았었는데
이녀석은 내 데스크탑에서 usb꼽으면 빨강파랑초록 왔다갔다거리고
그래프도 막 나오고 조만간에는 아이퐁에서도
시퀀싱을 해볼 수 있다는 찌라시도 가끔나오다가 업자들의 입에서도 심심치 않게 나오고..

그리고 이거 돌아가는 꼬라지를 보니 그냥 시덥지 않은 우스개 소리가 아니라
이제 곧 현실이 될 것 같기도 한데 그림은 봐도 잘 모르겠고..

그래서 옥스포드 나노포어를 조금 쉽게 이해해보는 글을 하나 투척합니다.
이게 real은 아닙니다. 이해를 쉽게 돕고자 제가 이해한 바를 풀어 해친겁니다.




http://blogs.nature.com/naturejobs/2017/10/16/techblog-the-nanopore-toolbox/


최근 NGS시대를 맞이하면서 일루미나가 한번런에 몇샘플 할 수 있다 생산량이 500G다. 시퀀싱기계를 10EA 묶음 판매한다. 몇일이면 몇백명 exome 시퀀싱을 해서 이제 시퀀싱 단가 10만원이 실현됐다, 팩은 10kbp read 시퀀싱한다 많이들 떠들고 있습니다.

일단 현재 전세계에서 범용적으로 운영되고 있는 시퀀싱플랫폼은 대게 3개정도 입니다.
일루미나/PacBio/PGM 이 3대 시퀀싱 플랫폼이 모두 엄밀히 말하면 시퀀싱의 부산물을 가지고 시퀀싱(basecalling)하고 있는 것입니다. 일루미나와 PacBio은 형광, PGM은 H+

근데 가장 핫한 Nanopore의 경우 현존하는 시퀀싱 플랫폼과 달리 직접 서열을 확인하는 기기라고 말씀드릴 수 있습니다.
엄밀히 얘기하자면 직접 들여다 보는건 아니고 전기신호로 각 base의 그림자를 확인해서 실체를 알아낸다고 하는게 더 맞을 수 있습니다.


이게 무슨말이냐!!!

옛날 옛적에 8-90년대 국민학교시절때 드래곤볼인가 어떤 게임 카드 뒤에 아래와 같은 글씨

가 적혀 있던 것을  기억하시는분!!! 분명 있으실겁니다.
이것은 숫자인데 특정 카드 집에 카드를 꼽으면 해당 카드의 공격력(?)을 알 수 있게 해놓은 카드 입니다. 쉽게 숫자를 판독할 수 없게 해서 꼭 카드집에 꼽아야 해당 숫자를 확인 할 수 있게... 물론 시간이 흐르면 대충 저 기호만 봐도 공격력(?)을 알 수 있게 학습된다는... ㅋㅋ









쉽게 생각하면 나노포어가 이와 같이 시퀀싱을 해석한다고 보시면됩니다.
원래 숫자(A/C/G/T)가 있으나 우리는 그것을 직접 알 수 없지만 저런 기호(전기적신호)로는 저장을 할 수 있다!



그래서 저 기호를 해독 할 수 있는 카드집에 저 기호를 대보면 원래 숫자를 알 수 있게 되는것이죠



좀 이해가 되시나요?

근데 이 숫자들이 단지 숫자로만 되어있다면 간단한데
숫마뿐만 아니라 영문자도 섞여 있다면 이 카드집으로도 간단할까요?

위의 그림은 숫자와 영문자가 있는 기호라고 합니다.

이것을 지금까지 숫자만 판독하던 카드집에 넣고 확인을 하면 다음과 같습니다.


어랏!! 일단 읽혀지긴 합니다.
그럼 이것을 어떻게 읽으시겠습니까?

A125c947로 읽으실건가요?
AI2Sc947? 아니면 A12Sc947? 이런 문제가 여기서 뿐만아니라
나노포어에서도 생기게 됩니다.

나노포어가 기존 시퀀서들과 다른점이 여기에 있습니다.

기존 시퀀서들의 경우 A/G/C/T에 각각에 형광을 부착하던지
혹은 현재 합성되는 base가 무엇인지 이미 알고 있습니다. 시퀀싱, basecalling을 합니다.
(그럼에도 불구하고 플랫폼의 한계들로 에러가 포함되고 있습니다.)

그러나 나노포어의 경우 사전정보가 없습니다. 그냥 AGCT가닥을 pore에 집어넣어서
전기적 신호의 차이를 읽어서 시퀀싱을 하기때문에 이 전기적 신호에 대한 연습이 필요하게 됩니다.

위에서 나온 두번째와 네번째 기호에 대한 구분을 어떻게 할지 연습이 필요한것 처럼 말입니다.


위의 방법은 두가지가 있습니다.
측정한 자료를 마사지해서 지금써오던 단순한 카드집을 계속 사용한다.
아니면 새로운 카드집으로 변경한다.

첫번때는 정보의 손실과 기존 카드집에서보다 더큰 에러를 발생할 수 잇는것이죠 A라는 정보는 원래 없었으니 A와 비슷한 8로 치환, S자도 어차피 없는 정보니 S를 5로 오독하는것이죠

이러면 안되겠죠
그래서 지금까지 써왔던 단순한 카드집 대신
조금더 세련된 1과 I, 5와 S를 구분 할 수 있는
카드집이 필요한것이죠


나노포어도 계속 기존 정보를 마사지하는것이 아니라 새로운 카드집을 계속 개선해나가고 있다고 생각하시면 될것 같습니다.



조만간 영국에서 Brexit말고 나노포어의 얘기가 건너오길 기대하면서
낚시글은 이만 줄이도록 하겠습니다. :)


(2017년 2월 18일 글 중간 중간 수정 및 첨언을 하였습니다.)
(2021년 1월 1일 이미지를 나노포어 설명과 관련된 이미지를 새로 업로드하였습니다.)


금요일, 9월 14, 2012

파일의 포맷을 변환하는데 필요한 것들

내가 아니란 말이닷!!! ㅋㅋ

python에서 Biopython을 이용하여
간단하게 convert하는 샘플 코드를 제공하고 있으니
여러분들도 쉽게 만들수 있어요~ :)
Biopython에서 제공하는 Tutorial 


오늘 문의가 들어온 파일은 sff파일
Roche의 454 GS FLX? sequencing 결과파일로....
ABI와 함께 illumina한테 밀려서 뒷방으로 들어앉은 파일 포맷입니다.
그러나 아직도 쓰는 이유는 read 길이가 길기때문 :)

그렇습니다. PacBio도 Nanopore다 디립다 길게 sequencing해준다는
애들이 있습니다. 그런데 왜 옛날꺼 쓰냐?? PacBio는 base quality가 안습이고,
Nanopore는.... 언제 출시일지 전 잘 모르겠습니다. 업자가 아닌관계로 ㅎㅎ

그래서 위의 길게 sequencing 해준다는 시퀀서를 제외하고는 Roche의 454가 read 길이가 가장 길다고 할 수 있겠습니다. NGS중에선 말이죠

그런데 sff파일을 보려고 하면 문제가 생깁니다.
권모씨께서 문의를 한것이 그것때문인지는 모르겠지만 걍 일반인이
sff파일을 걍 직접 볼수가 없습니다. 왜냐 binary파일이니깐요(sff파일이 binary라고
알고 있는데  직접 다뤄본적이 없어서... ㅎㅎ )

그래서 사람이 볼수 있게 파일을 변환시켜줘야 한다는 겁니다.

convertSff.py
#!/usr/bin/python

import os, sys
from Bio import SeqIO

try:
inputSFF = sys.argv[1]
outputPREFIX = sys.argv[2]

except:
print "Usage: python convertSFF <input.sff> <output_name>"
print ""
exit(1)


SeqIO.convert(inputSFF,"sff","%s.fasta"%(outputPREFIX), "fasta")
SeqIO.convert(inputSFF,"sff","%s.quality"%(outputPREFIX), "qual")
SeqIO.convert(inputSFF,"sff","%s.fastq"%(outputPREFIX), "fastq")


권모씨의 요청으로 급조한 날림 convert python 코드 ㅋㅋ
이 스크립트를 수행하면 세개의 파일이 나오게 될것으로 예상됩니다. ㅎㅎ
안나오면 어쩔수없고... ㅎㅎ


아.. 그리고 사족으로 LT사의 SOLiD의 경우 우리가 알고 있는 서열과 달리
첫 염기 서열만 서열이고 그 다음부터는 A/G/T/C 알파벳이 아닌 숫자로 되어있는데..
이걸 굳이 변환해서 reference geneome에 mapped 작업하지 말라고 합니다.
Re-sequencing하는 경우라면 변환해서 mapping하지 말고 원래 원본 파일 그대로를
input으로 하는 align 프로그램을 사용해서 mapped한 다음에 그 다음 작업을
일반적으로 사용하는 samtools나 GATK같은 프로그램을 사용하라고 합니다.
(다들 알고있는거 한번더 상기 시켜드렸습니다. 혹시 아나요 SOLiD 포맷을 분석하게 될지.. ㅎㅎ)

분석시 raw 파일을 사용해야 하는 이유는 SOLiD만의 월등한 quality 효과를 볼수 있어서
그러지 않겠나하는....  믿거나 말거나 저 혼자만의 생각입니다.. ㅎㅎ
다만, 타사 제품과 다르게 복잡하게 숫자로 표현한건 아니겠죠...
나름의 숨은 뜻이.... 쿨럭.. (설마... 간지용;;;;; )

Re-sequencing이 아닌 denovo일 경우 모 어쩔수 없이 fastq파일로 변환을 해야 하지 않을까 합니다. assembly 프로그램을 작동시키려면 아무래도 SOLiD format보다는 fastq 포맷이
수월하니깐요.. :)

그럼....