바이러스 유전자 서열을 코돈끼리 서로 옆에 등장하는 관계망(그래프)으로 바꿔 변이를 구분하는 법
바이러스 유전자 서열을 코돈끼리 서로 옆에 등장하는 관계망(그래프)으로 바꿔 변이를 구분하는 법
코로나19 변이(베타, 감마, 델타, 오미크론)를 구분하던 기존 방법은 유전자 서열을 일렬로 늘어놓고 비교하는 방식이라 서열 안의 복잡한 상호관계를 놓쳤다. 연구팀은 GenEx라는 파이프라인을 만들어 유전자 서열을 코돈(염기 3개 단위) 단위로 쪼갠 뒤, 이 코돈들이 서로 이웃하거나 함께 나타나는 관계를 그래프(점과 선으로 이루어진 관계망)로 표현했다. 이렇게 만든 그래프에서 25개 넘는 구조적 특징을 뽑아 23개의 머신러닝 모델에 학습시킨 결과, 최고 98.75%의 정확도로 변이를 구분해냈다.
METAL MEDIA 해설 도표
바이러스 유전자 서열을 코돈끼리 서로 옆에 등장하는 관계망(그래프)으로 바꿔 변이를 구분하는 법
- 01유전자 서열을 코돈 단위로 자른 뒤, 코돈이 옆에 붙어 나오는 관계(LAPCG)와 여러 거리에서 함께 나타나는 관계(MSCG) 두 가지 방식으로 그래프를 만들었다.
- 02이 그래프에서 노드 개수, 최단경로 길이, 밀집도 같은 구조적 특징과 고유값 기반의 스펙트럴 특징 등 25개 이상을 뽑아냈다.
- 03특이값 분해(SVD)에서 나온 값을 제곱해서(σ²) 쓰는 방식으로, 변이 간 차이를 더 뚜렷하게 구분되게 만들었다.
- 0423개 머신러닝 모델로 검증한 결과 MSCG 방식이 가장 우수했고, 정확도 98.75%를 기록하면서도 기존 PMI 방식보다 1.88배 빠르게 처리됐다.
- 05각 변이마다 고유한 그래프 구조(예: 오미크론은 매칭 넘버가 다른 변이의 18배, 델타는 가장 안정적인 구조)를 확인해, 서열만 보는 계통분석과는 다른 관점의 진화 정보를 얻었다.
무엇을 했나
- 유전자 서열을 코돈 단위로 자른 뒤, 코돈이 옆에 붙어 나오는 관계(LAPCG)와 여러 거리에서 함께 나타나는 관계(MSCG) 두 가지 방식으로 그래프를 만들었다.
- 이 그래프에서 노드 개수, 최단경로 길이, 밀집도 같은 구조적 특징과 고유값 기반의 스펙트럴 특징 등 25개 이상을 뽑아냈다.
- 특이값 분해(SVD)에서 나온 값을 제곱해서(σ²) 쓰는 방식으로, 변이 간 차이를 더 뚜렷하게 구분되게 만들었다.
- 23개 머신러닝 모델로 검증한 결과 MSCG 방식이 가장 우수했고, 정확도 98.75%를 기록하면서도 기존 PMI 방식보다 1.88배 빠르게 처리됐다.
- 각 변이마다 고유한 그래프 구조(예: 오미크론은 매칭 넘버가 다른 변이의 18배, 델타는 가장 안정적인 구조)를 확인해, 서열만 보는 계통분석과는 다른 관점의 진화 정보를 얻었다.
| Aspect | Prior Work | Our Method |
|---|---|---|
| Biological unit | Nucleotide / segment | Codon |
| Representation | Sequence / static graph | Co-occurrence graph |
| Temporal modeling | No | Yes |
| Learning framework | CNN / Transformer / rules | ML on graph features |
| Explainability | Post-hoc | Intrinsic (graph-based) |
| Multi-task reuse | Limited | Yes |
| Method | Best Model(s) | Accuracy | Runtime (s) |
|---|---|---|---|
| PMI baseline | CatBoost | 96.25% | 272.3 |
| LAPCG | LightGBM, Grad. Boost | 96.25% | 196.2 |
| MSCG | MLP, Bagging Classifier | 98.75% | 144.6 |
![Fig. 3: Parallel-coordinate comparison of ten structural graph features across four SARS-CoV-2 variants of concern (Beta, Oct 2020; Gamma, Nov 2020; Delta, Apr 2021; Omicron, Nov 2021). Each polyline represents one variant; color encodes variant class. The ten axes correspond to: second eigenvalue (sev), average shortest path length (aspl), Wiener index (wi), top eigenvalue (tev), radius (r), graph energy (ge), diameter (d), average clustering coefficient (aclu), average closeness centrality (aclo), and transitivity (t). All values are min–max normalised to [0,1]. Omicron consistently occupies the highest band across spectral and centrality axes (ge, aclu, aclo), while Delta remains near the minimum on most features, highlighting a sharp structural divergence between these two variants.](https://media.metallab.ai/papers/2608.18238/f2.png)
| Method | Type | Task | Accuracy (reported) | Whole Genome | Reference |
|---|---|---|---|---|---|
| GENSCAN | Ab initio HMM | Gene annotation | ∼70 - 80% (protein-level) | ✓ | [2] |
| GeneMark-ES/ET | Ab initio (self-train) | Gene annotation | 35.7 – 75.8% (nucleotide) | ✓ | [10] |
| SNAP | Ab initio HMM | Gene annotation | 77 - 80% (nucleotide) | ✓ | [9] |
| GlimmerHMM | Ab initio HMM | Gene annotation | ∼9 – 43% (some datasets) | ✓ | [11] |
| MAKER2 | Annotation pipeline | Gene annotation | 68.60% (nucleotide) | ✓ | [8] |
| AUGUSTUS | Ab initio + evidence | Gene annotation | 82–92% (gene-level) | ✓ | [14] |
| BRAKER2 | Hybrid (RNA-seq + HMM) | Gene annotation | >AUGUSTUS (+2–3%) | ✓ | [7] |
| Helixer (DL) | Deep learning | Gene annotation | 86.8% (reported) | ✓ | [15] |
| k-mer + SVM | Traditional ML | Variant classification | ∼92.01% | ✓ | Various |
| CNN-LSTM | Deep learning | Variant classification | ∼95–97% | Partial | Various |
| GenEx LAPCG (Ours) | Graph + ML (O(n)) | Variant classification | 96.25% | ✓ | This work |
| GenEx MSCG (Ours) | Graph + ML (multi-scale) | Variant classification | 98.75% | ✓ | This work |
![Fig. 4: Hierarchical Clustering of SARS-CoV-2 Variants by Graph Features. The ten features in x-axis correspond to: second eigenvalue (sev), average shortest path length (aspl), Wiener index (wi), top eigenvalue (tev), radius (r), graph energy (ge), diameter (d), average clustering coefficient (aclu), average closeness centrality (aclo), and transitivity (t). All values are min–max normalised to [0,1]. Omicron consistently occupies the highest band across spectral and centrality axes (ge, aclu, aclo), while Delta remains near the minimum on most features, highlighting a sharp structural divergence between these two variants.](https://media.metallab.ai/papers/2608.18238/f3.png)
왜 중요한가
이 방법은 기준 서열과 정렬 작업 없이도 대규모 바이러스 유전자 데이터를 빠르게 분류할 수 있어, 새로운 변이가 등장했을 때 실시간에 가깝게 감시하는 데 활용될 수 있다. 코로나19에 국한되지 않고 다른 바이러스나 생물종의 유전자 분석에도 확장 가능한 틀을 제시한다는 점에서 의미가 있다.
이 논문의 용어
- 코돈 · DNA/RNA 염기 3개가 모여 하나의 아미노산 정보를 나타내는 단위
- PMI(점별 상호정보량) · 두 항목이 우연보다 얼마나 자주 함께 나타나는지를 수치화한 값
- 특이값 분해(SVD) · 행렬을 여러 성분으로 분해해 데이터의 핵심 구조를 뽑아내는 수학적 기법
- ANOVA · 여러 집단 간 평균 차이가 통계적으로 유의한지 검정하는 방법
- 매칭 넘버 · 그래프에서 서로 겹치지 않는 최대 개수의 연결(간선) 집합의 크기
최신 논문
- AI 코딩 에이전트에게 과학 소프트웨어 수리를 시켜보니, 절반도 제대로 못 고쳤다AI 코딩 에이전트에게 과학 소프트웨어 수리를 시켜보니, 절반도 제대로 못 고쳤다
- 논문 속 시연이 아니라 실제 서비스에 넣을 수 있는 희소 어텐션 만들기논문 속 시연이 아니라 실제 서비스에 넣을 수 있는 희소 어텐션 만들기
- 고객상담 AI 상담원이 규정을 '한 번의 행동'이 아니라 '전체 절차'로 지키게 만드는 방법고객상담 AI 상담원이 규정을 '한 번의 행동'이 아니라 '전체 절차'로 지키게 만드는 방법
- 로봇 팔에게 사람의 시연 없이 새 일 시키기, 말 잘하는 AI가 대신 가르친다로봇 팔에게 사람의 시연 없이 새 일 시키기, 말 잘하는 AI가 대신 가르친다
- 에이전트 학습용 환경을 새로 만드는 대신, 기존 환경에 '패치 부품'을 씌워 그 에이전트의 약점에 맞게 바꾸는 방법에이전트 학습용 환경을 새로 만드는 대신, 기존 환경에 '패치 부품'을 씌워 그 에이전트의 약점에 맞게 바꾸는 방법
- AI 모델을 '소유'하지 못한 조직은 안전 통제도 절반밖에 못 한다AI 모델을 '소유'하지 못한 조직은 안전 통제도 절반밖에 못 한다
- AI가 선생님 모델을 따라 배우다가, 정답에 다가가는 '좋은 생각'까지 억누르는 문제를 잡아낸다AI가 선생님 모델을 따라 배우다가, 정답에 다가가는 '좋은 생각'까지 억누르는 문제를 잡아낸다
- AI가 특정 사람 말투를 흉내내도록 시켜봤더니, 결국 AI 자신의 말투에서 못 벗어난다AI가 특정 사람 말투를 흉내내도록 시켜봤더니, 결국 AI 자신의 말투에서 못 벗어난다
METAL MEDIA 최신 기사
그림 출처: Arefin Amin et al., arXiv:2608.18238, CC BY 4.0