바이설파이트 시퀀싱 (Bisulfite sequencing)
쉽게 풀면
바이설파이트 시퀀싱은 DNA에서 메틸화라는 화학적 표지가 붙어 있는 위치를 정확하게 찾아내는 방법입니다. 특수한 화학 처리를 하면 메틸화되지 않은 시토신은 다른 염기로 바뀌지만, 메틸화된 시토신은 그대로 남습니다. 이 처리 전후의 서열을 비교하면 원래 DNA에서 어느 시토신이 메틸화되어 있었는지 염기 하나 단위까지 정밀하게 알아낼 수 있습니다.
왜 중요한가
DNA 메틸화는 유전자 발현을 조절하는 대표적인 후성유전학적 표지로, 암을 비롯한 여러 질환에서 정상 조직과는 다른 메틸화 패턴이 관찰된다. 바이설파이트 시퀀싱은 이러한 메틸화 위치를 염기 단위로 정밀하게 파악할 수 있어, 후성유전체 연구뿐 아니라 발달생물학, 노화 연구, 암 바이오마커 발굴 등 폭넓은 분야에서 표준 기법으로 자리 잡았다. 최근에는 질병 진단이나 조직 기원 추정을 위한 메틸화 기반 바이오마커 개발에도 널리 활용되고 있다.
어떤 절차인가
- DNA 추출·단편화: 고품질 DNA를 쓰고, 필요하면 단편화하거나 메틸화된 어댑터를 붙입니다(어댑터의 시토신이 변환 대상이 되지 않도록 메틸화 어댑터를 쓰는 라이브러리가 흔합니다).
- 바이설파이트 처리: 아황산수소나트륨이 단일가닥 DNA에서 메틸화되지 않은 시토신(C)을 탈아미노화해 우라실(U)로 바꿉니다. 5-메틸시토신(5mC)은 이 반응이 훨씬 느려 C로 남습니다(Frommer et al., 1992).
- 증폭·시퀀싱: PCR 과정에서 U는 T로 읽힙니다. 따라서 시퀀싱 결과에서 기준 서열의 C가 그대로 C로 읽히면 메틸화, T로 읽히면 비메틸화로 해석합니다.
- 정렬과 계산: 변환으로 서열이 3글자(A, G, T 위주)에 가까워지므로 전용 정렬 도구(예: Bismark, Krueger와 Andrews, 2011)를 써서 위치마다 C 읽기와 T 읽기 수를 셉니다.
언제 쓰고 무엇을 얻나
- 전장유전체 바이설파이트 시퀀싱(WGBS): 유전체 거의 전체의 CpG(와 CHG, CHH)를 염기 단위로 얻습니다(Lister et al., 2009). 비용과 시퀀싱 양이 큽니다.
- 축소 대표(RRBS): MspI로 자른 CpG 밀집 조각만 모아 적은 비용으로 CpG 섬과 프로모터를 훑습니다(Meissner et al., 2005). 유전체의 나머지는 보지 못합니다.
- 한계: 5mC와 5-하이드록시메틸시토신(5hmC)을 구분하지 못합니다(둘 다 C로 읽힘). 이를 가르려면 oxBS-seq나 TAB-seq가 쓰입니다. 바이설파이트 처리는 DNA를 많이 손상하므로 효소로 변환하는 EM-seq 같은 방법이 대안으로 제시되었습니다(Vaisvila et al., 2021).
- 얻는 값: 위치별 메틸화 비율 = C 읽기 ÷ (C 읽기 + T 읽기). 그룹 비교는 차등메틸화영역(DMR) 분석으로 합니다.
가상의 예: 한 CpG 위치에서 C로 읽힌 읽기 38개, T로 읽힌 읽기 12개이면 메틸화 비율은 38 ÷ 50 = 0.76(76%)입니다. Wilson 방법의 95% 신뢰구간은 약 63~86%로, 읽기 수가 50개에서는 폭이 상당히 넓습니다. 또한 라이브러리에 섞어 넣은 비메틸화 대조 DNA(예: 람다 파지 DNA)에서 5,000개 중 15개가 C로 읽혔다면 변환 실패율은 15 ÷ 5,000 = 0.3%, 변환율은 99.7%입니다. 이 비율만큼 비메틸화 위치가 메틸화로 오인되므로 관측 비율 p에서 m = (p − e) ÷ (1 − e) = (0.76 − 0.003) ÷ (1 − 0.003) ≈ 0.7593로 보정하기도 합니다(변환 실패가 모든 위치에서 같다고 가정한 단순 보정).
논문에는 이렇게 적는다
- 입력 DNA 양, 변환 키트 이름, 라이브러리 제작 방법, 시퀀싱 장비와 읽기 길이.
- 변환율(비메틸화 대조 DNA, 또는 메틸화 가능성이 없는 서열에서 구한 값)과 그 값이 얼마인지.
- 정렬 도구와 기준 유전체 버전, 중복 읽기 제거 여부, 위치별 최소 읽기 깊이(예: 5× 또는 10× 이상만 사용하는 등 연구가 정한 기준), 메틸화 비율을 구한 방법.
- 통계 검정 방법(DMR 호출 도구, 다중 검정 보정).
- 예시 문장: “게놈 DNA 500 ng을 바이설파이트 처리한 뒤 라이브러리를 만들고 Illumina 장비로 150 bp 쌍 말단 시퀀싱하였다. 읽기는 Bismark로 인간 기준 유전체(hg38)에 정렬하였고, 비메틸화 람다 DNA 대조의 변환율은 99.5%였다.” (가상의 예문)
논문에서는 이렇게 쓰입니다
실제 논문에서
이 용어를 직접 쓴 국내 논문 문장을 아직 수집하지 못했습니다. 수집되는 대로 출처와 함께 싣겠습니다.
예시 문장
아래는 용어의 쓰임을 보여 주려고 만든 예시로, 실제 논문에서 가져온 문장이 아닙니다.
암 조직에서 특정 유전자 조절 부위의 메틸화 변화를 정량적으로 분석한 예문이다.
발달생물학 분야에서 세포 분화에 따른 메틸화 재프로그래밍을 관찰한 예문이다.
미생물생태학 연구에서 환경 요인과 메틸화 패턴 사이의 관계를 살펴본 예문이다.
조금 더 깊게 보면
바이설파이트 시퀀싱에는 유전체 전체를 다루는 전장유전체 바이설파이트 시퀀싱(WGBS)과, CpG가 밀집한 CpG 섬 위주로만 효율적으로 훑는 축소 대표 바이설파이트 시퀀싱(RRBS)이 있으며, 연구 목적과 비용에 따라 선택된다.
분석 결과는 흔히 개별 시토신의 메틸화 수준을 0에서 1 사이 비율로 나타내는 메틸화 비율(methylation ratio)이나, 특정 구간의 메틸화 차이를 비교하는 차등메틸화영역(DMR) 분석으로 요약된다.
다만 바이설파이트 처리는 5-메틸시토신과 5-하이드록시메틸시토신을 구분하지 못한다는 한계가 있어, 이를 구별하려면 옥시바이설파이트 시퀀싱 등 보완적인 기법이 함께 사용되기도 한다.
주의할 점
바이설파이트 처리 과정에서 DNA가 손상될 수 있어 충분한 양의 고품질 시료를 확보하는 것이 중요하다.
헷갈리기 쉬운 용어
| 용어 | 한 줄 설명 |
|---|---|
| 바이설파이트 시퀀싱 | 화학 변환 후 시퀀싱해 염기 단위로 5mC 위치를 읽음 |
| DNA 메틸화 프로파일링 | 메틸화 양상을 측정하는 방법 전반(바이설파이트 시퀀싱, 마이크로어레이, 농축 기반 방법 포함) |
| RRBS와 WGBS | RRBS는 CpG 밀집 구간만(저비용), WGBS는 거의 유전체 전체(고비용) |
| oxBS-seq·TAB-seq·EM-seq | 5hmC를 구분하거나 DNA 손상을 줄이는 변형·대안 방법 |
| DNA 메틸화 | 측정 대상인 화학적 표지 자체 |
참고 문헌
- Frommer, M., et al. (1992). A genomic sequencing protocol that yields a positive display of 5-methylcytosine residues in individual DNA strands. Proceedings of the National Academy of Sciences, 89.
- Meissner, A., et al. (2005). Reduced representation bisulfite sequencing for comparative high-resolution DNA methylation analysis. Nucleic Acids Research, 33.
- Lister, R., et al. (2009). Human DNA methylomes at base resolution show widespread epigenomic differences. Nature, 462.
- Krueger, F., & Andrews, S. R. (2011). Bismark: a flexible aligner and methylation caller for Bisulfite-Seq applications. Bioinformatics, 27.
- Vaisvila, R., et al. (2021). Enzymatic methyl sequencing detects DNA methylation at single-base resolution from picograms of DNA. Genome Research, 31.
표의 기준값은 위 문헌에 근거한 관례이며, 분야·학술지마다 다른 기준을 쓸 수 있습니다.