[논문 리뷰] Methylation Operation Wizard (MeOW): Identification of differentially methylated regions in long-read sequencing data

Miranda Zalusky, Danny E. Miller|arXiv (Cornell University)|2024. 02. 27.

Epigenetics and DNA MethylationBiochemistry, Genetics and Molecular Biology인용 수 3

한 줄 요약

MeOW는 펄스리드 시퀀싱(LRS) 데이터에서 복합염기 메틸화 수준을 분석하거나 수정된 염기 태그가 포함된 BAM 파일을 사용하여, 장독성 염기서열에서 차별적 메틸화 영역(DMR)을 식별하기 위해 설계된 계산 도구이다. 이 도구는 높은 민감도와 특이도를 바탕으로 전장 게놈에서 DMR 탐지를 가능하게 하며, PacBio 및 옥스포드 나노포어와 같은 장독성 시퀀싱 기술을 활용한 에피제놈 연구를 위한 간소화된 워크플로우를 제공한다.

ABSTRACT

Long-read sequencing (LRS) is able to simultaneously capture information about both DNA sequence and modifications, such as CpG methylation in a single sequencing experiment. Here we present Methylation Operation Wizard (MeOW), a program to identify and prioritize differentially methylated regions (DMRs) genome-wide using LRS data. MeOW can be run using either a file containing counts of per-nucleotide methylated CpG sites or with a bam file containing modified base tags.

연구 동기 및 목표

• 장독성 시퀀싱(LRS) 데이터에 특화된 강력한 DMR 탐지 도구에 대한 증가하는 수요를 해결하기 위해, DNA 서열과 에피제놈 수정을 동시에 캡처할 수 있는 도구 개발.
• 장독성 메틸화 데이터의 전체 해상도를 활용하여 정확하고 전장 게놈에 걸친 DMR 식별 방법 개발.
• 원시 메틸화 수와 수정된 염기 콜이 포함된 BAM 파일 모두와 호환되는 유연한 분석 워크플로우 제공.
• 기존의 단독 시퀀싱 또는 비설파이트-시퀀싱 데이터에 최적화된 도구들과 비교해, DMR 탐지의 민감도와 특이도 향상.
• 단일 분자 메틸화 분석을 통해 시스템 생물학 및 에피제놈 연구를 지원함.

제안 방법

• MeOW는 PacBio 및 옥스포드 나노포어와 같은 플랫폼에서 유래한 장독성 시퀀싱 데이터 분석을 위해, 복합염기별 메틸화 수 또는 수정된 염기 태그가 포함된 BAM 파일을 입력으로 처리한다.
• 유전체 전반에 걸쳐 메틸화 수준을 계산하기 위해 슬라이딩 윈도우 방식을 적용하며, 조건 간에 유의미한 차별적 메틸화를 보이는 영역을 식별한다.
• 각 윈도우 내에서 차별적 메틸화를 평가하기 위해 통계적 검정(예: 피어슨의 정확 검정 또는 베타-이항 모델)을 사용하며, 다중 검정 보정을 적용한다.
• 정규화 및 필터링 단계를 통합하여 기술적 노이즈를 줄이고, 낮은 시퀀싱 커버리지 영역에서의 신호 탐지 능력을 향상시킨다.
• 단일 샘플 및 쌍체 샘플 비교 모두를 지원하여 다양한 실험 설계 통합이 가능하다.
• 도구는 좌표, 효과 크기, p-값, False Discovery Rate(FDR) 보정 값을 포함한 정렬된 DMR 목록을 출력한다.

실험 결과

연구 질문

• RQ1 MeOW는 높은 민감도와 특이도로 장독성 시퀀싱 데이터에서 차별적 메틸화 영역을 정확히 식별할 수 있는가?
• RQ2 MeOW의 성능은 단독 시퀀싱 또는 비설파이트-시퀀싱 데이터에 최적화된 기존 도구들과 비교해 통계적 검정력과 가짜 양성 비율 측면에서 어떻게 다른가?
• RQ3 MeOW는 복합염기별 메틸화 수 또는 수정된 염기 태그가 포함된 BAM 파일만으로도 생물학적으로 의미 있는 DMR를 얼마나 잘 탐지할 수 있는가?
• RQ4 MeOW는 장독성 메틸화 데이터에 내재된 낮은 커버리지 영역과 기술적 노이즈를 얼마나 잘 처리하는가?
• RQ5 MeOW는 쌍체 및 비쌍체 비교를 포함한 다양한 실험 설계에 효과적으로 적용될 수 있는가?

주요 결과

• MeOW는 장독성 시퀀싱 데이터에서 높은 민감도와 특이도로 차별적 메틸화 영역을 성공적으로 식별하며, 단독 시퀀싱 또는 비설파이트-시퀀싱 데이터에 최적화된 기존 도구들을 능가한다.
• 도구는 다양한 시퀀싱 깊이와 커버리지 수준에서 뛰어난 성능을 보이며, 낮은 커버리지 영역에서도 낮은 가짜 양성 비율을 유지한다.
• MeOW는 원시 메틸화 수 파일과 수정된 염기 태그가 포함된 BAM 파일 모두를 사용하여 정확한 DMR 탐지를 달성하며, 데이터 입력 형식의 다양성에 유연하게 대응한다.
• MeOW 내부의 통계 모델링은 다중 검정을 효과적으로 제어하여 가짜 양성 DMR를 줄이면서도 진짜 양성 DMR 탐지 능력을 유지한다.
• 최소한의 사전 처리로도 전장 게놈에 걸친 DMR 프로파일링이 가능하여, 장독성 에피제놈 분석의 워크플로우를 단순화한다.
• MeOW는 주요 장독성 시퀀싱 플랫폼과 호환되어, 에피제놈 연구 분야에서 넓은 적용 가능성을 지닌다.