Lesson 13 of 20
거르개가 지운 변이
논문 Results · Supplementary Figure 1 · EGFR 재분석
11강의 파이프라인을 다 돌린 뒤, 논문 저자들은 이상한 것을 발견했습니다.
EGFR 변이가 한 건도 없다
이것이 왜 이상한가. 앞선 연구들이 IP-SNSCC 에서 EGFR 엑손 20 삽입이 반복해서 나온다고 보고했기 때문입니다. 이 병에서 가장 자주 언급되는 표적 가능 변이입니다.
없을 리가 없는 것이 없었습니다.
그래서 EGFR 만 따로, COSMIC 조건을 빼고 다시 돌렸습니다. 결과입니다.
엑손 19 결손1명
엑손 20 삽입3명
합10명 중 4명
같은 자료입니다. 같은 리드입니다. 거르개 한 칸을 뺐을 뿐입니다.
이 재분석에는 대가도 있다
논문은 이것을 숨기지 않습니다. Discussion 에 이렇게 적혀 있습니다 —
EGFR 에 초점을 맞춘 재분석은 최초 파이프라인에서 EGFR 변이가
없다는 것 때문에 촉발되었으므로 편향을 들여올 수 있다.
이것이 사후 분석(post-hoc) 의 문제입니다. 결과를 본 뒤에 규칙을 바꾸면, 원하는 답이 나올 때까지 규칙을 바꿀 수 있게 됩니다.
정직한 처리는 이렇습니다 — 재분석 결과를 본문에 넣되, 왜 다시 돌렸는지 밝히고, 같은 완화를 다른 유전자에는 적용하지 않았음을 적는다. 논문은 그렇게 했습니다.
규칙을 바꿀 때마다 위양성 기회가 늘어납니다. 그래서 무엇을 왜 바꿨는지 적는 것이 유일한 방어입니다.
Python
Runs in your browser
그리고 이번에는 다른 것도 나왔다
기존 보고들은 대부분 인프레임(in-frame) 엑손 20 삽입만 보고했습니다. 이 논문은 인프레임과 프레임시프트 를 둘 다 찾았습니다.
인프레임 삽입
아미노산 몇 개가 끼어듭니다.
단백질은 그대로 만들어지고, 그 자리가 살짝 벌어져 계속 켜집니다.
→ EGFR 억제제가 노릴 수 있습니다.
단백질은 그대로 만들어지고, 그 자리가 살짝 벌어져 계속 켜집니다.
→ EGFR 억제제가 노릴 수 있습니다.
프레임시프트
읽는 틀이 어긋납니다.
그 뒤로는 엉뚱한 아미노산이 붙다가 대개 일찍 끝납니다.
→ 켜질 리가 없습니다. 약도 듣지 않습니다.
그 뒤로는 엉뚱한 아미노산이 붙다가 대개 일찍 끝납니다.
→ 켜질 리가 없습니다. 약도 듣지 않습니다.
논문의 지적 — 이 둘은 생물학적으로 다르며, EGFR 표적 치료의 반응에도 다르게 걸릴 것입니다. 'EGFR 변이가 있다' 는 한 줄로 묶어 버리면 이 차이가 사라집니다.
🤔 Check yourself
첫 분석에서 EGFR 변이가 하나도 나오지 않은 직접적 원인은?
Answer 논문의 설명 그대로입니다. 엑손 20 삽입은 heterogeneous in sequence and length 라서 COSMIC 에 incompletely represented 되어 있습니다. 필터는 잡음뿐 아니라 데이터베이스가 아직 모르는 것도 지웁니다.
🤔 Check yourself
EGFR 엑손 20 프레임시프트 변이를 인프레임 삽입과 따로 봐야 하는 이유는?
Answer 논문의 문장 — 'As frameshift EGFR alterations are unlikely to be activating, they may differ biologically from in-frame exon 20 insertions and potentially affect sensitivity to EGFR-directed therapies.' 변이를 유전자 이름으로만 묶으면 이 차이가 사라집니다.