[논문 리뷰] Modelling Breakage-Fusion-Bridge Cycles as a Stochastic Paper Folding Process

CD Greenman, SL Cooke외 3명|arXiv (Cornell University)|2012. 11. 10.

Algorithms and Data Compression참고 문헌 8인용 수 4

한 줄 요약

이 논문은 암 게놈에서의 Breakage-Fusion-Bridge (BFB) 사이클을 확률적 종이 접기 과정으로 모델링하며, 2차원 트리와 재귀적 접기 수열을 사용해 구조적 진화를 표현한다. 연구에서는 $n$ 사이클에 대해 $2^{n(n-1)/2}$개의 서로 다른 BFB 진화 경로가 존재한다고 유도하고, 시퀀싱 데이터로부터 접기 위치, 복제 수 증폭, 진화 가능성 등을 추론할 수 있는 확률적 프레임워크를 수립하여, 주요 암 샘플에서 복잡한 증폭 구역을 재구성하는 데 기여한다.

ABSTRACT

Breakage-Fusion-Bridge cycles in cancer arise when a broken segment of DNA is duplicated and an end from each copy joined together. This structure then `unfolds' into a new piece of palindromic DNA. This is one mechanism responsible for the localised amplicons observed in cancer genome data. The process has parallels with paper folding sequences that arise when a piece of paper is folded several times and then unfolded. Here we adapt such methods to study the breakage-fusion-bridge structures in detail. We firstly consider discrete representations of this space with 2-d trees to demonstrate that there are 2^(n(n-1)/2) qualitatively distinct evolutions involving n breakage-fusion-bridge cycles. Secondly we consider the stochastic nature of the fold positions, to determine evolution likelihoods, and also describe how amplicons become localised. Finally we highlight these methods by inferring the evolution of breakage-fusion-bridge cycles with data from primary tissue cancer samples.

연구 동기 및 목표

• 종이 접기 수열의 원리들을 활용하여 암 게놈에서의 Breakage-Fusion-Bridge (BFB) 사이클의 복잡한 구조적 진화를 모델링하기 위해.
• 2차원 트리를 사용하여 BFB 진화의 이산적 공간을 표현하고, 질적으로 다른 진화 경로를 수량화하기 위해.
• 접기 위치에 대한 확률적 성격을 통합하여 BFB 진화의 가능성을 계산하고 복제 수 증폭을 예측하기 위해.
• 실제 암 샘플의 복제 수 프로파일에서 유래한 접기 구조와 BFB 사이클의 진화적 역사를 추론하기 위해.
• genomic 복제 수 프로파일과 BFB에 의해 유도된 증폭 구역의 접기 패턴 및 절단점 분포 간의 수학적 프레임워크를 수립하기 위해.

제안 방법

• 2차원 트리를 사용하여 BFB 진화 경로의 이산적 공간을 표현하며, $n$ 사이클에 대해 $2^{n(n-1)/2}$개의 서로 다른 진화 경로를 보여준다.
• 각 접기가 DNA 절단 및 복구 사건에 해당하는 오리가미와 유사한 재귀적 접기 수열을 적용한다.
• BFB 과정을 마코프 길이 과정 $L_n$으로 모델링하며, $P(L_n|L_{n-1}) = \frac{1}{2L_{n-1}}$ for $0 \leq L_n \leq 2L_{n-1}$를 만족하여 마틴게일 성질 $\mathbb{E}[L_n|L_{n-1}] = L_{n-1}$을 확보한다.
• BFB 세그먼트 길이의 확률적 진화를 모델링하기 위해 연합 밀도함수 $P(L_1, \dots, L_n) = \frac{1}{2^n L \prod_{k=1}^{n-1} L_k}$를 도출한다.
• 각 사이클 동안 접기의 생존 확률을 순환적으로 계산하기 위해 함수 $W_k(x,y) = \int_x^{2y} \frac{W_{k-1}(x,z)}{z} dz$를 도입한다.
• 실험적 리드 깊이와 예측된 BFB 접기 패턴 및 증폭 수준 간의 연결을 위해 복제 수 프로파일과 접기 수 벡터 ($cn$, $fn$)를 사용한다.

실험 결과

연구 질문

• RQ1 BFB 사이클 $n$개에서 질적으로 다른 진화 경로는 총 몇 개인가?
• RQ2 임의의 절단점 선택에 기반해 주어진 BFB 접기 구조의 확률적 가능성은 무엇인가?
• RQ3 BFB에 의해 유도된 증폭 구역에서의 복제 수 증폭은 접기 수와 접기 순서와 어떻게 관련이 있는가?
• RQ4 실험적 복제 수 및 리드 깊이 데이터로부터 BFB 증폭 구역의 기저 접기 구조를 추론할 수 있는가?
• RQ5 BFB 사이클 전반에 걸쳐 최소 세그먼트 길이의 분포는 어떻게 되며, 이는 이전 접기의 생존에 어떤 영향을 미치는가?

주요 결과

• 2차원 트리 표현의 조합론적 분석을 통해 $n$ 개의 BFB 사이클에 대해 정확히 $2^{n(n-1)/2}$개의 질적으로 다른 진화 경로가 존재한다.
• 마틴게일 성질 $\mathbb{E}[L_n|L_{n-1}] = L_{n-1}$ 덕분에 초기 길이 $L$에서 시작하여도 BFB 세그먼트의 기대 길이가 사이클 전반에 걸쳐 일정하게 유지된다.
• $N$ 개의 BFB 중 $k^{\text{th}}$ BFB의 최소 길이가 $x$일 확률은 $M_{k,N}(x,L) = \frac{1}{2^k L} W_k(x,L) \cdot \left(1 - \sum_{i=1}^{N-k} \frac{1}{2^i} W_i(x,x) \right)$로 주어지며, 이는 전진 및 후행 생존 확률을 통합한 것이다.
• 함수 $W_k(x,y)$는 초기 조건 $W_1(x,y) = \log(2y) - \log(x)$를 만족하는 재귀적 적분 방정식 $W_k(x,y) = \int_x^{2y} \frac{W_{k-1}(x,z)}{z} dz$를 만족하며, 접기 생존 가능성의 반복적 계산을 가능하게 한다.
• 모델은 BFB 증폭 구역에서 최고 복제 수 영역(예: 영역 II)이 최대 18개의 복제를 달성할 수 있다고 예측하며, 이 중 16개는 BFB 복제로, 나머지 2개는 원래 염색체 복제로 기인한다. 이는 관측된 리드 깊이 프로파일과 일치한다.
• 모델은 실제 주요 암 샘플 데이터로부터 BFB 사이클의 접기 구조를 성공적으로 추론하였으며, 예측된 접기 위치가 실험적으로 확인된 절단점과 복제 수 변화와 일치한다.