삼중음성유방암에서 높게 나타난 MEAK7

공개 유방암 데이터를 통합 분석한 연구에서 MEAK7 유전자 발현이 삼중음성유방암(TNBC)에서 높고, 높은 발현이 짧은 무병생존기간과 연관된 것으로 나타났다. 다만 TNBC 환자의 전체생존이나 원격전이 없는 생존과는 통계적으로 유의한 연관성이 확인되지 않았다. Biology에 실린 이번 연구는 예후 표지자와 후속 기능 연구의 후보를 제시한 생물정보학 분석이다.

TNBC는 에스트로겐 수용체, 프로게스테론 수용체, HER2가 음성인 유방암이다. 연구진은 세포 성장·대사에 관여하는 mTORC1 신호의 대체 조절자로 알려진 MEAK7에 주목했다. TCGA·GTEx의 RNA 시퀀싱 자료와 bc-GenExMiner의 마이크로어레이 자료 등을 이용해 정상 조직과 종양, 유방암 아형 사이의 발현 차이와 생존 지표를 비교했다.

UALCAN 분석에서 MEAK7 발현 중앙값은 TNBC 12.3 TPM으로, 정상 인접 조직 7.8 TPM, 루미날 아형 5.9 TPM, HER2 양성 아형 9.7 TPM보다 높았다. TPM은 RNA 발현량을 비교하기 위한 정규화 단위다. 마이크로어레이 분석에서도 TNBC 및 기저세포양 아형에서 높은 발현이 관찰됐다. 다만 연구진은 TNBC와 비TNBC의 발현 차이를 중간 정도의 효과 크기로 평가했다.

무병생존과 전체생존에서 달랐던 결과

TNBC에서 MEAK7 발현이 중앙값을 초과한 환자군은 중앙값 이하인 군보다 무병생존(DFS)이 짧았다. 위험비는 1.27, 95% 신뢰구간은 1.03–1.57이었다. 이는 무병생존 분석에서 사건 발생의 상대적 위험이 높았다는 뜻이며, 생존기간이 27% 줄었다는 의미는 아니다. 반면 전체생존(OS)과 원격전이 없는 생존(DMFS)의 차이는 유의하지 않았다. 기저세포양 유방암에서도 DFS와의 연관성은 나타났지만 OS·DMFS에서는 유의한 차이가 없었다.

나이와 종양 병기를 보정한 분석에서도 MEAK7 발현과 전체생존의 연관성이 남았다. 그러나 이는 TNBC만이 아니라 TCGA-BRCA 유방암 환자 1,057명을 대상으로 한 결과다. 위험비는 1.21, 95% 신뢰구간은 1.02–1.45였다. 이를 TNBC에서 독립적 예후 가치가 확립됐다는 뜻으로 확대할 수 없다. 보정 변수에도 종양 등급, 림프절 상태, 치료 방식 등이 포함되지 않았다.

저메틸화와 RNA 상관관계가 제공한 단서

TNBC에서는 MEAK7의 발현 조절 영역인 프로모터의 메틸화 수준이 정상 인접 조직과 루미날 아형보다 낮았다. 연구진은 이를 높은 발현을 설명할 잠재적 조절 기전으로 해석했다. 하지만 저메틸화가 발현 증가를 직접 유발했다는 인과관계는 입증하지 않았다.

유방암 코호트에서 MEAK7 발현은 miR-149-3p·miR-135a-5p·LINC00993 발현과 음의 상관관계를, miR-135b-5p·HIF1A-AS2 발현과 양의 상관관계를 보였다. 이는 발현이 함께 변하는 패턴이며, 서로 직접 억제하거나 활성화한다는 증거는 아니다. HPA의 유방암 면역조직화학 자료에서는 MEAK7 단백질 염색이 낮음 또는 중간 수준으로 나타났다. RNA 발현 증가가 단백질 수준에서도 일관되게 나타나는지는 추가 확인이 필요하다.

CRISPR 설계에서 기능 검증으로

연구진은 MEAK7을 겨냥하는 CRISPR 가이드 RNA 20개를 계산적으로 설계했다. 일부 후보는 높은 예측 점수를 보였고 분석 범위에서 비표적 부위가 확인되지 않았지만, 다른 후보에서는 잠재적 비표적 부위가 나타났다. 실제 절단 효율이나 안전성, MEAK7 억제의 항암 효과를 확인한 결과는 아니다.

여러 분석 플랫폼이 TCGA 환자 자료를 일부 공유하므로, 반복된 결과를 모두 독립적 검증으로 볼 수는 없다. 다수의 생존 지표·하위집단 분석에 대한 정식 다중검정 보정이 없다는 한계도 있다. 제공된 본문에는 TNBC 생존 분석의 추적기간이 명시되지 않았다. 연구진은 이번 결과를 가설 생성 단계의 관찰로 규정하며 기능 실험과 전향적 검증이 필요하다고 밝혔다. 분석 방법과 결과는 논문 전문에서 확인할 수 있다.