KR20140121523A - 위암에 대한 예후 예측 시스템 - Google Patents

위암에 대한 예후 예측 시스템 Download PDF

Info

Publication number
KR20140121523A
KR20140121523A KR20130037335A KR20130037335A KR20140121523A KR 20140121523 A KR20140121523 A KR 20140121523A KR 20130037335 A KR20130037335 A KR 20130037335A KR 20130037335 A KR20130037335 A KR 20130037335A KR 20140121523 A KR20140121523 A KR 20140121523A
Authority
KR
South Korea
Prior art keywords
expression
value
gastric cancer
normlogtransvalue
rna
Prior art date
Application number
KR20130037335A
Other languages
English (en)
Other versions
KR101504818B1 (ko
Inventor
허용민
노성훈
서진석
정재호
박은성
Original Assignee
연세대학교 산학협력단
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 연세대학교 산학협력단 filed Critical 연세대학교 산학협력단
Priority to KR1020130037335A priority Critical patent/KR101504818B1/ko
Publication of KR20140121523A publication Critical patent/KR20140121523A/ko
Application granted granted Critical
Publication of KR101504818B1 publication Critical patent/KR101504818B1/ko

Links

Images

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6844Nucleic acid amplification reactions
    • C12Q1/6851Quantitative amplification
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/48Biological material, e.g. blood, urine; Haemocytometers
    • G01N33/50Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
    • G01N33/53Immunoassay; Biospecific binding assay; Materials therefor
    • G01N33/574Immunoassay; Biospecific binding assay; Materials therefor for cancer
    • G01N33/57407Specifically defined cancers
    • G01N33/57446Specifically defined cancers of stomach or intestine
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2537/00Reactions characterised by the reaction format or use of a specific feature
    • C12Q2537/10Reactions characterised by the reaction format or use of a specific feature the purpose or use of
    • C12Q2537/165Mathematical modelling, e.g. logarithm, ratio
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/118Prognosis of disease development

Landscapes

  • Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Organic Chemistry (AREA)
  • Molecular Biology (AREA)
  • Biochemistry (AREA)
  • Urology & Nephrology (AREA)
  • Biotechnology (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Microbiology (AREA)
  • Analytical Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Hematology (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Genetics & Genomics (AREA)
  • Hospice & Palliative Care (AREA)
  • Oncology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Cell Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Biophysics (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • General Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

본 발명은 위암의 예후 예측이 가능한 신규한 예후 예측 시스템에 관한 것으로, 보다 상세하게는, 유전자 집합 발현 비교 분석법을 통해 위암의 수술에 의한 절제 후 임상 결과를 예측할 수 있다.

Description

위암에 대한 예후 예측 시스템{Novel system for predicting prognosis of gastric cancer}
본 발명은 유전자 발현 비교 분석법을 통해 위암의 예후 예측이 가능한 신규한 예후 예측 시스템에 관한 것이다.
위선암(Gastric adeno-carcinoma)은 2000년 700,349명의 사망에서 두 번째 원인으로, 세계에서 가장 일반적으로 진단된 네 번째 암이다. 몇 가지 역학적 및 조직병리학적 특징들을 갖는 단일 이질적 질환으로 간주하고 있다. 위암 치료는 주로 환자를 수술 만으로 또는 수술과 화학요법으로 치료하여야 하는지를 결정하는 TNM(tumor, node, metastasis) 병기결정 같은 임상적 파라미터에 근거한다. 위암은 유방암과 대장암 등과 달리 TNM 병기 시스템에 따라서 1기에서 4기까지 명확하게 차이가 난다. 즉, 1기의 경우에는 5년 생존율이 90% 이상이며, 4기의 경우에는 20% 이하로 큰 차이를 보인다. 그러므로 TNM 병기 시스템의 예후 예측력이 매우 뛰어남을 알 수 있다 [참고문헌, 7th edition of the AJCC cancer staging Manual: stomach. Ann Surg Oncol 2010;17:3077-3079]. 상기 병기 시스템에 기반을 두어 위암은 흔히 조기 위암(Early Gastric Cancer), 국소진행형(Locally Advanced Gastric Cancer), 국소 침윤형(Locally Advanced Invasive Gastric Cancer) 및 전이 위암(Metastatic Gastric Cancer) 등으로 나눌 수 있다.
비록 수술이 실시 가능한 위암의 주요 치료이긴 하나, 진행성인 경우 재발률이 높다. 재발을 예방하고, 위암 환자들의 예후를 개선하기 위해 화학요법과 화학-방사선요법을 포함한 집학적 치료가 도입되었다. 그러나, 이들 치료 방법이 환자들에서 일반적인 임상 결과를 개선하기는 하나 종양의 임상병리학적 이질성과, 같은 병기에 있는 환자들의 다른 결과는 어쥬번트 화학요법의 임무를 예측하는데 한계가 있어 개별 환자들에 대한 최적 접근이 부족한 상태이다.
최근의 진전에도 불구하고 발병적으로 별개의 종양 유형에 대해 특이적 치료 섭생을 표적화하고, 궁극적으로 성과를 최대화시키기 위하여 종양 치료를 개인화하기 위한 암 치료의 도전과제들이 남아있다. 따라서, 각종 치료 선택사항들에 대한 환자 반응에 관한 예측적 정보를 동시에 제공하는 시험을 필요로 하고 있다.
1. 대한민국 공개특허 제2007-0022694호, 2007.02.27 2. 대한민국 공개특허 제2010-0072283호, 2010.06.30
1. Journal of clinical oncology, vol. 23(29), pp.7286-7295, 2005.10.10
본 발명의 목적은 TNM 병기와 상관없는 전체 위암 환자들의 임상 결과에 영향을 주는 유전자 발현 예후 지표 또는 위험도 점수에 기초한 새로운 예후 예측 시스템을 제공하는 것이다.
상기 목적을 달성하기 위하여, 본 발명은 대상으로부터 얻은 암세포를 포함하는 생물학적 샘플에서 CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 결정하는 단계; 및 상기 단계에서 결정된 RNA 전사체의 발현도에 기초하여 상기 생물학적 샘플의 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 계산하고, 상기 RS(%)에 따라 예후를 판단하는 단계를 포함하는, 위암으로 진단된 대상에서 예후를 예측하는 방법을 제공한다.
본 발명은 또한 대상으로부터 얻은 암세포를 포함하는 생물학적 샘플에서 CSK, CD8A 및 CHEK2의 RNA 전사체의 발현 수준을 측정하는 단계; 및 상기 전사체의 발현 증가는 긍정적인 임상 결과 가능성의 증가로 판단하는 단계를 포함하는, 위암으로 진단된 대상에서 예후를 예측하는 방법을 제공한다.
상기 예후 예측 방법은 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는 것일 수 있다.
본 발명은 또한 위암의 예후 예측을 실행하는 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체에 있어서, 환자로부터 얻은 핵산 시료에서 CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 결정하는 단계; 및 상기 단계에서 결정된 RNA의 발현도에 기초하여 상기 생물학적 샘플의 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 계산하고, 전체 생존율(Overall Survival, OS) 측면에서 상기 RS(%)의 설정치 범위가 75% 이상을 갖는 경우 고 위험군 환자, 25% 이상 내지 75% 미만인 경우 중간 위험군 환자, 25% 미만인 경우 저 위험군 환자로 분류하는 단계를 컴퓨터에 실행시키는 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체를 제공한다.
상기 기록 매체는 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는 것일 수 있다.
본 발명은 TNM 병기와 상관없는 전체 위암 환자군에 대한 전체 생존율 측면에서 예측 모델을 만든 후 통계적으로 유의한 생존에 영향을 미치는 RNA 전사체의 발현도를 결정하여 이로부터 위험도 평점 시스템을 만들어 예후 지표 값을 산출하는 방식으로 위암 수술에 의한 절제 후 임상 결과를 예측할 수 있다.
도 1은 연세대학교 세브란스 병원에서 2001년부터 2004년까지 총 520명을 대상으로 위암 병기에 따른 분포를 조사한 결과이다.
도 2는 레퍼런스 유전자에 대한 표준화 작업을 수행한 결과를 나타낸 것이다.
도 3은 RNA 전사체 발현도에 따른 예후가 좋은 군과 나쁜 군의 생존율 분석 결과를 나타낸 것이다.
도 4는 BRB-Array Tools를 이용한 PCA 분석에 따라 생존 위험 예측 결과를 나타낸 것이다.
도 5는 본 발명의 예후 모델의 정확도를 검증한 결과를 나타낸 것이다.
도 6은 위암 환자군을 대상으로 위험도 점수(RS(%))에 따른 생존 분석 결과를 나타낸 것이다.
이하 본 발명의 구성을 구체적으로 설명한다.
본 발명은 대상으로부터 얻은 암세포를 포함하는 생물학적 샘플에서 CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 결정하는 단계; 및 상기 단계에서 결정된 RNA 전사체의 발현도에 기초하여 상기 생물학적 샘플의 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 계산하고, 상기 RS(%)에 따라 예후를 판단하는 단계를 포함하는, 위암으로 진단된 대상에서 예후를 예측하는 방법에 관한 것이다.
본 발명에 따른 예후를 예측하는 방법은 위암 환자들의 예후와 관련을 보이는 유전자와 종양생물학적인 특성을 강하게 나타내는 유전자군을 선정하여 qPCR을 수행하여 콕스 회귀 분석에서 통계적 유의(p<0.001)를 갖는 유전자를 예후와 관련된 유전자 대상으로 선정하고, 상기 유전자들의 위험 계수(hazard ratio)를 유전자의 발현 값에 곱하여 하기 수학식 1 및 2에 따라 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 구한 후 전체 생존율(Overall Survival, OS) 측면에서 상기 RS(%)가 75% 이상을 갖는 샘플을 고 위험군, 25% 이상 내지 75% 미만인 경우 중간 위험군, 25% 미만인 경우 저 위험군으로 분류함으로써 위암으로 진단된 대상에서 예후를 예측할 수 있는 것을 특징으로 한다.
상기 RS 및 RS(%)는 하기 수학식 1 및 2에 따라 계산할 수 있다:
[수학식 1]
RS=HR1*normLogTransValue1 + HR2*normLogTransValue2 + ... + HRn* normLogTransValuen
[수학식 2]
RS(%)=100×(생물학적 샘플의 RS - 모집단의 RS 최소값)/(모집단의 RS 최고값 - 모집단의 RS 최소값)
상기 식에서,
HRn는 n번째 RNA 전사체의 위험 계수(hazard ratio)를 나타내고, 상기 HRn이 1 미만인 경우 -1/HRn로 변환하여 사용하며,
normLogTransValuen는 RNA 전사체의 발현과 관련된 값이고, 이 값은 해당 유전자의 전체 값을 대상으로 하여 중간값을 중심으로 스케일을 변화시킨 값이며,
상기 모집단은 TNM 병기와 상관없는 전체 위암을 갖는 일정 수의 집단을 의미하고, 일정 수는 RS 최고값과 최소값을 산출할 수 있는 정도의 임의의 정수이다.
상기 모집단의 수는 특별히 제한하지는 않으며, 일 구체예에서는 TNM 병기와 상관 없는 전체 위암 조직 520개를 모집단으로 사용하였다.
본 명세서에서, 용어, 위험 계수(Hazard Ratio: HR)란 암의 진행, 재발, 또는 요법 반응에 대한 기여도를 반영하는 계수를 의미한다. 위험 계수는 다양한 통계적 기법에 의하여 도출될 수 있다. 상기 위험 계수, HR 값은 다양한 통계적 모델에서 결정할 수 있으며, 예컨대, 다변량 콕스 비례 위험 회귀 분석에서 결정할 수 있다. 일 구체예에서, HR 값을 PI 수식에 사용함에 있어서, HR 값이 1보다 크거나 같을 경우 HR 값을 그대로 사용하고, HR 값이 1 보다 작을 경우 -1/HR 값을 사용할 수 있다.
또한, 상기 식에서, 용어, RNA 전사체의 발현 값이란 개별유전자, 즉 RNA 전사체의 발현과 관련된 값을 의미한다. 상기 값은 공지된 다양한 통계적 수단을 사용하여 결정할 수 있다. 예를 들어 발현 값은 콕스 회귀 분석에 의해 측정된 p 값을 log2 함수값으로 변형 후 사분위수 표준화(quantile normalization) 후의 값을 사용할 수 있다. 상기 수학식 1에서 사용된 발현 값은 해당 유전자의 전체 값을 대상으로 하여 중간값을 중심으로 스케일을 변화시킨 값을 사용하였다.
일 구체예에 따르면, RS는 다음과 같이 결정할 수 있다:
RS = - CSK×1.317 - CD8A×1.18 - CHEK2×1.26
상기 수학식 1에 따라 계산된 RS는 상기 수학식 2에 따라 RS(%)로 나타낼 수 있다. 위의 식에서 결정된 값을 모집단에서 해당하는 순위로 변환하여 전체 생존율(Overall Survival, OS) 측면에서 상기 RS(%)가 75% 이상을 갖는 샘플을 고 위험군, 25% 이상 내지 75% 미만인 경우 중간 위험군, 25% 미만인 경우 저 위험군으로 분류한다. 고 위험군은 나쁜 예후를, 저 위험군은 좋은 예후인 것으로 판단할 수 있다. 즉, RS(%) 값이 75% 이상을 갖는 샘플을 고 위험군은 3년 이상, 6년 이상, 10년 이상의 기간 동안 전체 생존율이 낮음을 의미하고, 25% 미만인 저 위험군은 3년 이상, 6년 이상, 10년 이상 동안 전체 생존율이 높음을 의미한다. 상기 용어, 좋은 예후는 임상 결과의 긍정적 임상 결과 가능성의 증가로 표현될 수 있고, 나쁜 예후는 임상 결과의 긍정적 임상 결과 가능성의 감소로 표현될 수 있다.
상기 방법은 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는데 유용할 수 있다.
본 발명은 또한 대상으로부터 얻은 암세포를 포함하는 생물학적 샘플에서 CSK, CD8A 및 CHEK2의 RNA 전사체의 발현 수준을 측정하는 단계; 및 상기 전사체의 발현 증가는 긍정적인 임상 결과 가능성의 증가로 판단하는 단계를 포함하는, 위암으로 진단된 대상에서 예후를 예측하는 방법을 제공한다.
상기 방법은 qPCR 기반 방법일 수 있다.
상기 발현 수준이 하나 이상의 RNA 전사체의 발현 수준에 대해 표준화되는 것일 수 있다.
상기 임상 결과가 전체 생존율(Overall Survival, OS) 측면에서 표현되는 것일 수 있다.
상기 방법은 RNA 전사체 전체의 발현 수준을 측정하고 발현 증가를 분석하여 긍정적인 임상 결과 가능성의 증가 또는 감소를 판단하여 예후를 예측할 수 있다.
상기 방법은 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는데 유용할 수 있다.
본 발명은 또한 위암의 예후 예측을 실행하는 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체에 있어서, 환자로부터 얻은 핵산 시료에서 CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 결정하는 단계; 및 상기 단계에서 결정된 RNA의 발현도에 기초하여 상기 시료의 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 계산하고, 전체 생존율(Overall Survival, OS) 측면에서 상기 RS(%)의 설정치 범위가 75% 이상을 갖는 경우 고 위험군 환자, 25% 이상 내지 75% 미만인 경우 중간 위험군 환자, 25% 미만인 경우 저 위험군 환자로 분류하는 단계를 컴퓨터에 실행시키는 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체를 제공한다.
상기 기록 매체는 TNM 병기와 상관 없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는데 유용할 수 있다.
상기 RS 및 RS(%)는 상기 수학식 1 및 2에 따라 계산할 수 있다.
상기 기록 매체는 전체 생존율(Overall Survival, OS) 측면에서 RS(%)의 설정치 범위가 75% 이상을 갖는 경우 고 위험군, 25% 이상 내지 75% 미만인 경우 중간 위험군, 25% 미만인 경우 저 위험군으로 판단할 수 있다. 즉, RS(%) 값이 75% 이상을 갖는 샘플을 고 위험군은 3년 이상, 6년 이상, 10년 이상의 기간 동안 전체 생존율이 낮다고 판단하고, 25% 미만인 저 위험군은 3년 이상, 6년 이상, 10년 이상 동안 전체 생존율이 높다고 판단한다. 상기 용어, 좋은 예후는 임상 결과의 긍정적 임상 결과 가능성의 증가로 표현될 수 있고, 나쁜 예후는 임상 결과의 긍정적 임상 결과 가능성의 감소로 표현될 수 있다.
달리 정의되지 않는다면, 본원에 사용된 기술 및 과학 용어들은 당업자 수준에서 일반적으로 이해하는 바와 같은 의미가 있다. 본 발명은 어떤 방식으로든 설명된 방법 및 재료로 제한되지 않는다. 본 발명의 목적상, 하기 용어들이 아래에서 정의된다.
용어, "폴리뉴클레오티드"는 일반적으로 임의의 폴리리보뉴클레오티드 또는 폴리데옥시리보뉴클레오티드를 말하고, 예컨대 변형 또는 비-변형 RNA 또는 DNA일 수 있다. 본 명세서에서, "폴리뉴클레오티드"는 구체적으로 cDNA를 포함한다.
용어, "올리고뉴클레오티드"는 비제한적으로 한-가닥 데옥시리보뉴클레오티드, 한- 또는 두-가닥 리보뉴클레오티드, RNA:DNA 하이브리드 및 두-가닥 DNA를 포함하는, 비교적 짧은 폴리뉴클레오티드를 말한다. 올리고뉴클레오티드, 예를 들면 한-가닥 DNA 프로브 올리고뉴클레오티드는 종종 예를 들면 상업적으로 입수가능한 자동화 올리고뉴클레오티드 합성기를 사용하는 화학적 방법에 의해 합성된다. 그러나, 올리고뉴클레오티드는 시험관내 재조합 DNA-매개 기술을 포함하는 각종 다른 방법들에 의해 및 세포 및 유기체 중에서의 DNA 발현에 의해 제조될 수 있다.
용어, "차등적으로 발현된 유전자" 또는 "차등적인 유전자 발현"은 정상 또는 대조용 대상체의 발현에 비하여 위암과 같은 암을 앓는 대상체 중에서 더 높거나 낮은 수준으로 활성화되는 유전자를 말한다. 또한, 동일한 질병의 다른 병기에서 더 높거나 낮은 수준으로 활성화되는 유전자를 포함한다. 차등적으로 발현되는 유전자는 핵산 수준 또는 단백질 수준에서 활성화 또는 억제되거나, 다른 스플라이싱을 받아 상이한 폴리펩타이드 산물을 야기하는 경우일 수 있다. 이러한 차이는 예를 들면 폴리펩타이드의 mRNA 수준, 표면 발현, 분비 또는 다른 분배에 있어서의 변화에 의해 입증될 수 있다. 본 발명의 목적상, "차등적인 유전자 발현"은 정상 및 질병에 걸린 대상체에서 또는 질병에 걸린 대상체의 다양한 병기에서 주어진 유전자의 발현 사이에 약 1.5배 이상, 약 4배 이상, 약 6배 이상, 약 10배 이상의 차이가 있을 때 존재하는 것으로 간주 된다.
유전자 전사체 또는 유전자 발현 생성물에 관한 용어 "표준화된"은 기준 유전자 세트의 전사체/생성물의 평균 수준에 대한 전사체 또는 유전자 발현 생성물의 수준을 말하는데, 여기서 레퍼런스 유전자들은 환자, 조직 또는 치료에 걸쳐 이들의 최소한의 변동에 기준하여 선택되거나 ("하우스키핑 유전자(housekeeping genes)"), 또는 레퍼런스 유전자는 시험된 유전자들 전체를 말한다. 후자의 경우, 일반적으로 "전체 표준화(global normalization)"로 언급되는데, 시험된 유전자들의 총 수가 비교적 큰, 바람직하게는 50 초과인 것이 중요하다. 구체적으로, RNA 전사체에 관한 용어 '표준화된'은 기준 유전자 세트의 전사 수준의 평균에 대한 전사 수준을 말한다.
용어, "발현 역치" 및 "정의된 발현 역치"는 혼용하여 사용하며, 이 수준 이상에서는 유전자 또는 유전자 생성물이 환자 반응에 대한 예측 마커로서 사용되는 해당 유전자 또는 유전자 생성물의 수준을 말한다. 역치는 대표적으로 임상적 연구로부터 실험적으로 정의된다. 발현 역치는 최대 민감성, 또는 최대 선택성(예를 들면 한 약물에 대한 반응자들만을 선택하도록), 또는 최소 오차로 선택될 수 있다.
용어, "유전자 증폭"은 특정 세포 또는 세포주에서 유전자 또는 유전자 단편의 다수개의 복사물이 형성되는 과정을 말한다. 복제된 영역 (증폭된 DNA의 신장)은 종종 "암플리콘"으로 언급된다. 종종, 생산된 mRNA의 양, 즉 유전자 발현도는 또한 특정 유전자의 만들어진 복제 수에 비례하여 증가한다.
본 명세서에서, "예후"는 본원에서 암에 의한 사망 또는 위암과 같은 신생물성 질환의 진행(재발, 전이성 확산 및 내약물성 포함)의 가능성의 예측을 말하는데 사용된다. 용어 "예측"은 본원에서 환자가 주요 종양의 수술 제거 후에 암 재발 없이 특정 기간 동안 살아남게 될 가능성을 말하는데 사용된다. 이러한 예측은 임의의 특정 환자에 대하여 가장 적절한 치료 기법을 선택함으로써 치료를 결정하는데 임상적으로 사용될 수 있다. 이러한 예측은 환자가 치료 섭생, 예를 들면 수술 시술에 대하여 유리하게 반응하기 쉬운지, 또는 수술 종료 후에 환자의 장기간 생존이 가능한지를 예측하는데 있어서 귀중한 수단이 된다.
본 명세서에서, "장기간" 생존은 본원에서 수술 또는 다른 치료 후 3년 이상, 보다 바람직하게는 6년 이상, 가장 바람직하게는 10년 이상 생존하는 것을 말하는데 사용된다.
달리 지시하지 않는 한, 분자 생물학(재조합 기술 포함), 미생물학, 세포 생물학 및 생화학의 종래 기술을 사용하여 본 발명을 수행할 수 있다.
1. 유전자 발현 프로파일 작성(Profiling)
유전자 발현 프로파일 작성 방법은 폴리뉴클레오티드의 혼성화 분석에 기초한 방법, 폴리뉴클레오티드의 서열화에 기초한 방법, 및 프로테오믹스 기재 방법을 포함한다. 예를 들어 mRNA 발현의 정량화를 위한 방법은 노던 블랏팅(northern blotting) 및 인 시츄 혼성화(in situ hybridization); RNAse 보호 검정시험; 및 PCR-기재 방법, 예를 들면 역 전사 폴리머라제 연쇄 반응(RT-PCR) 등을 포함한다. 또는, DNA 두 가닥, RNA 두 가닥, 및 DNA-RNA 하이브리드 두 가닥 또는 DNA-단백질 두 가닥을 포함하는 특정 두 가닥을 인식할 수 있는 항체들이 사용될 수 있다. 서열화-기재 유전자 발현 분석에서 대표적인 방법은 유전자 발현의 연속 분석(Serial Analysis of Gene Expression, SAGE) 및 대량적으로 평행한 시그너쳐 서열화(massively parallel signature sequencing, MPSS)에 의한 유전자 발현 분석을 포함한다.
2. PCR-기반 유전자 발현 프로파일 작성 방법
a. 역 전사효소 PCR(RT-PCR)
가장 민감하고 가장 유연한 정량적 PCR-기재 유전자 발현 프로파일작성 방법들 중 하나는 RT-PCR이고, 이것은 약물 치료와 함께 또는 약물 치료 없이 정상 조직 및 종양 조직에서의 상이한 샘플 집단에 있어서의 mRNA 수준을 비교하여 유전자 발현 패턴을 특성화하고, 밀접하게 관련된 mRNA들을 판별하고, RNA 구조를 분석하는데 사용될 수 있다.
제1단계는 표적 샘플로부터 mRNA의 단리이다. mRNA는 예를 들면 냉동되거나 보관된 파라핀-매립 및 고정된 (예를 들면, 포르말린-고정된) 조직 샘플로부터 추출될 수 있다.
mRNA 추출을 위한 일반적인 방법은 당업계에 공지되어 있으며, 파라핀에 매립된 조직으로부터의 RNA 추출 방법은 문헌 ([Rupp and Locker, Lab Invest. 56: A67 (1987)] 및 [De Andres et al ., BioTechniques 18: 42044 (1995)]) 등에 개시되어 있다. 특히, RNA 단리는 상업적 제조업체, 예를 들면 퀴아겐(Qiagen)으로부터의 정제 키트, 완충제 세트 및 프로테아제를 사용하여 제조업체의 설명서에 따라 수행할 수 있다. 다른 상업적으로 입수가능한 RNA 단리 키트는 마스터푸어(MasterPure)TM 완전 DNA 및 RNA 정제 키트(에피센트레(EPICENTRE)®, 위스콘신주 매디슨) 및 파라핀 블록(Paraffin Block) RNA 단리 키트(앰비온, 인크.(Ambion, Inc.))를 포함한다. 조직 샘플로부터의 전제 RNA는 RNA Stat-60 (Tel-Test)를 사용하여 단리할 수 있다. 종양으로부터 제조된 RNA는 예를 들면 염화세슘 밀도 구배 원심분리에 의해 단리될 수 있다.
RNA는 PCR을 위한 주형으로 사용될 수 없기 때문에, RT-PCR에 의한 유전자 발현 프로파일 작성의 제1단계는 RNA 주형의 cDNA로의 역 전사이고, 이후 그의 PCR 반응으로의 지수적 증폭이 이어진다. 역 전사효소는 주로 조류 골수아세포증 바이러스 역 전사효소(AMV-RT) 및 몰로니(Moloney) 쥐 백혈병 바이러스 역 전사효소(MMLV-RT)를 사용할 수 있다. 역 전사 단계는 대표적으로 발현 프로파일작성의 환경 및 목표에 따라, 특정 프라이머, 무작위 헥사머, 또는 올리고-dT 프라이머를 사용하여 초회항원자극된다. 예를 들면, 추출된 RNA는 진앰프(GeneAmp) RNA PCR 키트(퍼킨 엘머(Perkin Elmer), 미국 캘리포니아주)를 사용하여 제조업체의 설명서에 따라 역-전사될 수 있다. 유도된 cDNA는 이어서 후속되는 PCR 반응에서 주형으로서 사용될 수 있다.
비록 PCR 단계가 각종 열안정성 DNA-의존성 DNA 폴리머라제를 사용할 수 있지만, 이것은 전형적으로는 Taq DNA 폴리머라제를 사용하는데, Taq DNA 폴리머라제는 5'-3' 뉴클레아제 활성을 갖지만, 3'-5' 판독방지(proofreading) 엔도뉴클레아제 활성은 부족하다. 따라서, 타크맨® PCR은 전형적으로는 Taq 또는 Tth 폴리머라제의 그의 표적 암플리콘에 결합된 혼성화 프로브를 혼성화시키는 5'-뉴클레아제 활성을 이용하지만, 5' 뉴클레아제 활성을 갖는 임의의 효소가 사용될 수 있다. 2개의 올리고뉴클레오티드 프라이머들을 사용하여 PCR 반응의 대표적인 암플리콘을 생성시킨다. 제3 올리고뉴클레오티드 또는 프로브는 2개의 PCR 프라이머들 사이에 위치한 뉴클레오티드 서열을 검출하도록 설계된다. 프로브는 Taq DNA 폴리머라제 효소에 의해 비-연신성이고, 리포터 형광 염료 및 소광제(quencher) 형광 염료로 표지된다. 리포터 염료로부터 임의의 레이저-유도 방출은 2개의 염료가 이들이 프로브 상에 있을 때와 같이 함께 가깝게 위치할 때 소광제 염료에 의해 소거된다. 증폭 반응 동안, Taq DNA 폴리머라제 효소는 주형-의존적 방식으로 프로브를 절단한다. 생성되는 프로브 단편들은 용액 중에서 해리되고, 방출된 리포터 염료로부터의 신호에는 제2 형광단의 소거 효과가 없다. 리포터 염료의 한 분자가 합성된 새로운 분자 각각으로부터 방출되고, 소거되지 않은 리포터 염료의 검출은 데이터의 정량적 해석에 대한 기준을 제공한다.
타크맨® RT-PCR은 상업적으로 입수 가능한 장비, 예를 들면 ABI 프리즘(PRISM) 7700TM 시퀀스 검출 시스템(Sequence Detection System)TM(Perkin-Elmer-Applied Biosystems), 또는 라이트사이클러(Lightcycler)(Roche Molecular Biochemicals)를 사용하여 수행될 수 있다.
5'-뉴클레아제 검정시험 데이터는 초기에 Ct, 또는 역치 사이클로서 표현된다. 형광 값은 매 사이클 동안 기록되고, 증폭 반응으로 그 지점으로까지 증폭된 생성물의 양을 나타낸다. 형광 신호가 처음 통계학적으로 유의한 것으로 기록될 때의 지점이 역치 사이클(Ct)이다.
샘플 간의 변동 효과 및 오차를 최소화시키기 위하여, RT-PCR은 일반적으로 기준 RNA(이것은 이상적으로는 상이한 조직들 사이에서 일정 수준으로 발현됨)를 사용하여 수행되고, 실험 치료에 의해 영향을 받지 않는다. 유전자 발현 패턴을 표준화하는데 가장 자주 사용되는 RNA는 하우스키핑 유전자 글리세르알데히드-3-포스페이트-데히드로게나제 (GAPD) 및 β-액틴(ACTB)의 mRNA이다.
또한, 새로운 실시간 정량적 PCR은 이중-표지된 형광원성 프로브(즉, 타크맨® 프로브)를 통한 PCR 생성물 축적을 측정하는 기술로, 정량적 경쟁적 PCR(여기서는, 각 표적 서열에 대한 내부 경쟁자가 표준화에 사용됨) 및 샘플 내에 포함된 표준화 유전자 또는 RT-PCR에 대한 하우스키핑 유전자를 사용하는 정량적 비교용 PCR 모두와 상용 가능하다.
b. 매스어레이(MassARRAY) 시스템
시커놈, 인크.(Sequenom, Inc.)가 개발한 매스어레이-기반 유전자 발현 프로파일 작성 방법에서는 RNA의 단리 및 역 전사 후에 얻은 cDNA를 합성 DNA 분자(경쟁자)(이것은 단일 염기를 제외한 모든 위치에서 표적화 cDNA 구역과 일치함)로 스파이크하고(spiked), 내부 표준으로 사용한다. cDNA/경쟁자 혼합물을 PCR 증폭시키고, 후-PCR 새우 알칼리성 포스파타제(SAP) 효소 처리를 가하여 남아있는 뉴클레오티드의 데포스포릴화를 야기한다. 알칼리성 포스파타제의 불활성화 후에, 경쟁자 및 cDNA로부터의 PCR 생성물을 프라이머 신장시키고, 이것은 경쟁자- 및 cDNA-유래 PCR 생성물에 대한 별도의 질량 신호들을 생성시킨다. 정제 후, 이들 생성물들을 매트릭스-보조 레이저 탈착 이온화 흐름 시간 질량 분광측정법(MALDI-TOF MS) 분석을 이용한 분석에 필요한 성분들이 이미-부하 되어 있는 칩 어레이 상에 계량분배한다. 반응에 존재하는 cDNA를 이어서 생성된 질량 스펙트럼 내의 피크 면적 비를 분석하여 정량화한다.
c. 기타 PCR-기반 방법
기타 시차 디스플레이; 증폭된 단편 길이 다형성 (iAFLP); 비드어레이(BeadArray) TM 기술 (일루미나(Illumina), 캘리포니아주 샌 디에고); 유전자 발현에 대한 신속 검정시험에 상업적으로 입수 가능한 루미넥스(Luminex)100 LabMAP 시스템 및 다색-코딩된 미소구(Luminex Corp.)를 사용하는, 유전자 발현 검출용 비즈어레이(BeadsArray for Detection of Gene Expression, BADGE); 및 고 피복 발현프로파일작성 (HiCEP) 분석이 있다.
3. mRNA 단리, 정제 및 증폭의 일반적인 설명
파라핀에 매립된 조직을 사용하여 유전자 발현 프로파일을 작성하는 기술은 상술한 바와 같다. 최종적으로 얻은 데이터를 분석하여 관찰된 종양 샘플에서 확인된 특징적인 유전자 발현 패턴에 기초하여 환자에게 이용할 수 있는 최상의 치료 선택사항(들)을 판별해낸다.
본 발명의 중요한 면은 위암 조직에 의한 특정 유전자의 측정된 발현을 사용하여 예후 정보를 제공하는 것이다. 이러한 목적을 위해, 검정 시험된 RNA의 양, 사용된 RNA 품질에 있어서의 변동, 및 다른 인자, 예를 들면 기계 및 작업자 차이에 있어서의 차이에 대해 보정하는(표준화) 것이 필수적이다. 그러므로, 검정시험은 전형적으로 GAPD 및 ACTB와 같은 공지된 하우스키핑 유전자로부터 전사된 것들을 포함하는, 기준 RNA의 사용을 측정하여 혼입시킨다. 또는, 표준화는 검정시험된 유전자들 또는 이들의 많은 서브세트 전부의 평균 또는 중간 신호(Ct)를 기준으로 할 수 있다(전체 표준화 접근법). 하기 실시예에서는 중심 표준화 전략을 사용하였는데, 이것은 표준화를 위해 임상적 성과와의 상관성 부족에 기초하여 선택된 스크리닝된 유전자의 서브세트를 이용하였다.
용어, "트레이닝 세트(Training set)"라 함은 예후에 통계적으로 유의한 RNA 전사체를 추출한 대상 샘플을 의미한다.
용어, "검증 세트(Validation set)" 또는 "테스트 세트(Test set)"라 함은 상기 추출된 변수가 실제로 예후의 좋고 나쁨을 판단할 수 있는 정확도를 테스트하는 세트를 의미한다. 이러한 방법을 사용하는 이유는 특정 샘플 군에서만 효과적으로 예후 판단 능력이 있는 것이 아니라 독립적 샘플에서도 효능이 있는 것을 판단하기 위해서이다.
4. 재발에 대한 예후 지표 및 이들의 응용
위암 재발의 가능성에 대한 암 예후 방법을 구분 짓는 연산의 특징은 1) 재발 가능성을 측정하는데 사용된 독특한 시험 mRNAs 세트 (또는 대응하는 유전자 발현 생성물), 2) 발현 데이터를 식으로 합치는데 사용된 특정 가중치, 및 3) 환자들을 상이한 수준의 위험을 갖는 군, 예를 들면 저, 중간 및 고 위험 군으로 나누는데 사용된 역치를 포함한다. 이 연산은 수치적인 위험도 점수(RS) 및 RS(%)를 산출해낸다.
시험은 명시된 mRNA 또는 이들의 발현 생성물의 수준을 측정하기 위한 실험실 검정시험을 필요로 하지만, 신선한 조직이거나 또는 냉동된 조직, 또는 이미 반드시 환자들로부터 수집되어 보관되어 있는 고정되어 파라핀에 매립된 종양 생검 시험편을 매우 소량으로 이용할 수 있다. 따라서, 시험은 비침입성일 수 있다. 예를 들면 코어 생검 또는 미세침 흡인을 통해 수확된 종양 조직의 몇 가지 상이한 방법들과 상용성이기도 하다. 이 방법에 따르면, 암 위험도 점수(RS)는
(a) 상기 대상체로부터 얻은 암 세포를 포함하는 생물학적 샘플로 유전자 또는 단백질 발현 프로파일을 작성하고;
(b) 다수개의 개별 유전자의 발현도 [즉, mRNA 수준]를 정량화하여 각 유전자에 대한 발현 값을 정하고;
(c) 각각 암-관련 생물학적 함수에 의해 및(또는) 동시발현에 의해 연결된 유전자들에 대한 발현 값을 포함하는, 유전자 발현 값들의 서브세트를 생성시키고;
(d) 한 서브세트 내의 각 유전자의 발현도에 상기 서브세트 내에서의 그의 암 재발 반응에 대한 상대적 기여도를 반영하는 계수를 곱하고 곱한 값을 더하여 상기 서브세트에 대한 값을 산출하고;
(e) 각 서브세트의 그 값에 그의 암 재발 반응에 대한 기여도를 반영하는 계수를 곱하고;
(f) 상기 계수를 곱한 각 서브세트에 대한 값들의 합을 구하여 위험도 점수(RS) 및 RS(%)를 얻음으로써 결정되는데, 여기서, 암 재발과 선형 상관관계를 보이지 않는 각 서브세트의 기여도는 단지 소정의 역치 값 이상에서만 포함되고,
명시된 유전자의 증가된 발현이 암 재발 위험을 감소시키는 서브세트에는 음의 값을 부여하고, 명시된 유전자의 발현이 암 재발 위험을 증가시키는 서브세트에는 양의 값을 부여한다.
구체적인 실시태양에서, RS 및 RS(%)는
(a) CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 측정하고,
(b) 하기 수학식 1 및 2에 의해 RS 및 RS(%)를 계산함으로써 결정된다:
[수학식 1]
RS=HR1*normLogTransValue1 + HR2*normLogTransValue2 + ... + HRn* normLogTransValuen
[수학식 2]
RS(%)=100×(생물학적 샘플의 RS - 모집단의 RS 최소값)/(모집단의 RS 최고값 - 모집단의 RS 최소값)
상기 식에서,
HRn는 n번째 RNA 전사체의 위험 계수(hazard ratio)를 나타내고, 상기 HRn이 1 미만인 경우 -1/HRn로 변환하여 사용하며,
normLogTransValuen는 RNA 전사체의 발현과 관련된 값이고, 이 값은 해당 유전자의 전체 값을 대상으로 하여 중간값을 중심으로 스케일을 변화시킨 값이며,
상기 모집단은 TNM 병기와 상관없는 전체 위암을 갖는 일정 수의 집단을 의미하고, 일정 수는 RS 최고값과 최소값을 산출할 수 있는 정도의 임의의 정수이다.
여기서, RS(%) 값이 75% 이상이면 나쁜 예후이고, 25% 미만이면 좋은 예후인 것으로 판단한다.
이하, 본 발명을 실시예에 의해 상세히 설명한다. 단, 하기 실시예는 본 발명을 예시하는 것일 뿐, 본 발명의 내용이 하기 실시예에 한정되는 것은 아니다.
<제조예> 예후 예측 대상 선정 및 실험설계
예후 예측 대상 선정을 위해, 2001년 내지 2004년까지 연세대학교 세브란스 병원에서 일차 치료로 위 절제술을 받은 위암 환자들(YUSH, n=520)에서 종양 표본 및 임상 데이터를 얻었다. 1기에서 4기까지 대상으로 하였으나 주로 2기-3기에 집중되어 있었다(도 1). 모든 샘플들은 환자로부터 상세하게 기재된 동의서를 받은 후 수집하였고, 연구는 연세대학교 세브란스 병원의 조사윤리위원회의 승인을 받았다. 임상 데이터는 소급하여 얻었다. 전체생존기간은 수술부터 사망에 이르는 시간으로 정하였고, 데이터는 환자가 마지막 접촉 시 살아있을 때 검열을 받은 것으로 간주하였다.
표 1은 위암 2-3기에서 RNA 전사체를 이용하여 예후 지표를 산출하기 위한 대상 유전자군(44개)과 보정을 위한 레퍼런스 유전자(10개)를 나열한 것이다.
상기 유전자들의 mRNA 종들의 수준을 RT-PCR로 측정하였다. mRNA를 마스터퓨어TM RNA 정제 키트 (에피센터 테크놀로지스(Epicentre Technologies))를 사용하여 추출하고 리보그린(RiboGreen)® 형광 방법(Molecular probes)에 의해 정량화하였다. 정량적 유전자 발현의 분자 검정시험을 ROCHE 480을 사용하여 수행하였다. ROCHE 480은 써모사이클러, 레이저, 전하-커플링된 장치(CCD), 카메라 및 컴퓨터로 이루어진다. 시스템은 샘플을 써모사이클러 상에서 384-웰 포맷으로 증폭시켰다. 증폭 동안, 레이저-유도 형광 신호가 모든 384-웰에 대하여 실시간으로 수집되고 CCD에서 검출되었다. 시스템은 기기를 실행하기 위한 및 데이터를 분석하기 위한 소프트웨어를 포함한다.
각 환자에 대한 역치 사이클(CT) 값을 임상적 성과와의 상관성 부족에 기초하여 선택된, 특정 환자에 대한 스크리닝된 유전자의 한 서브세트의 평균에 기초하여 표준화하였다(중심 표준화 전략).
유전자 발현을 분석하기 위한 보정 수단으로 레퍼런스 유전자를 활용하는 것이 상례이므로, "Normfiner" 프로그램을 사용하여 10개의 레퍼런스 유전자 중에서 6개의 유전자를 활용하여 표준화 작업을 수행하였다. 표준화에 사용된 ATP5E, GPX1, PGK1, UBB, VDAC2, ACTB이다(도 2).
qRT-PCR에 사용된 프라이머는 미국 특허 제7526387호의 표 1 및 2에 기재된 것을 사용하였다.
Figure pat00001
예후 예측 모형을 만들기 위해, 클러스터 분석은 Cluster and Treeview(http://rana.lbl.gov/EigenSoftware.htm)에 따라 수행하였다. 자율 클러스터링 분석은 연속 분산 여과 후 수행하였고, 2개의 주요 클러스터로 이루어진 2개의 분류군의 예후적 차이는 log rank test과 Kaplan Meier Plot에 따라 시험하였다.
<실시예 1> 위암 환자들의 유전자 발현 프로파일 조사
도 3은 520개 샘플을 qPCR한 결과의 히트맵을 나타낸 것으로, 붉은색은 저 발현, 녹색 계열은 고발현을 나타낸 것이며, 520개 샘플은 2개의 예후가 좋은 군(CL2)과 나쁜 군(CL1)으로 나눌 수 있는 것을 볼 수 있었다. Kaplan Meier Plot과 log rank test에서 p-값이 0.0233 이 나온 것을 알 수 있었다.
<실시예 2> BRB-Array Tools를 이용한 생존 위험 예측(PCA 분석)
하기 표 2는 위험 군의 분류에 사용된 유전자이며, BRB ArrayTools Version 4.1(http://linus.nic.nih.gov./BRB-ArrayTools.html)를 사용하여 PCA(Principal Component Analysis)를 통하여 예후 지표(Prognostic Index)를 계산하였다. 즉 wi 값과 유전자의 중간값화 된 발현수치를 곱하여 아래의 산식으로 산출한 값을 두 계층으로 나눈 것이다.
i w i x i - 0.007356
도 4는 BRB-Array Tools를 이용한 PCA 분석에 따라 생존 위험 예측 결과를 나타낸 것이다.
상기 산식에 따른 이 모델의 정확도 측정을 위해 CCP(Compound Covariate Predictor), LDA(Diagonal Linear, Discriminant Analysis), 1-NN(1-Nearest Neighbor), 3-NN(3-Nearest Neighbors), NC(Nearest Centroid), SVM(Support Vector Machines), BCCP(Bayesian Compound Covariate Predictor)로 실시한 클래스 예측 결과는 아래와 같았다(표 3. 도 5).
p-값 % CV Support Weights (wi) 명칭 심볼
1 0.00173 100 0.119466 CSK CSK
2 0.0012782 100 0.112631 CD8A CD8A
3 0.014359 90 0.106045 CHEK2 CHEK2
상기 예후 지표를 사용하여 상위 50% 하위 50%의 두 계층으로 나눈 결과 도 5와 같이 각 예측 모델 별로 CCP, LDA, 1-NN, 3-NN, NC, SMV, BCCP에서 구분되는 것을 알 수 있다.
CCP LDA 1-NN 3-NN NC SVM BCCP
76 80 83 82 78 92 89
<실시예 3> BRB-Array Tools를 이용한 생존 위험 예측(HR 분석)
도 3의 히트맵과 같이 유전자 발현을 콕스 회귀 분석(Univariate Cox's proportional harzard test)을 실행하여 3개의 유전자를 선정하였다. 사용 툴은 BRB ArrayTools Version 4.1(http://linus.nic.nih.gov./BRB-ArrayTools.html)를 사용하였고, 이를 이용하여 위험 계수(Hazard ratio)와 관련된 알고리즘을 설정하였다.
qPCR에 의한 예후 위험도 평점 시스템은 표준화된 유전자 발현량을 RG(Relative Quantity)로 변환하고, RQ로 변화된 유전자 발현량(Cq) value를 콕스 회귀 분석에 따라 위험 계수(hazard ratio, HR)를 구하였다. 위험 계수를 -1/HR로 변환하고, 유전자 발현량을 곱하여 총합(RS, Risk Score)을 산출하였다(컷 오프: -23 이하는 -23으로 간주함). 위험도 점수의 백분율(percentile risk score)은 다음 식에 따라 계산하였다:
RS(%)=100×(샘플의 RS - 모집단의 RS 최소값)/(모집단의 RS 최고값 - 모집단의 RS 최소값)
75% 이상의 RS(%)을 갖는 샘플은 고 위험군으로 분류하고, 25% 이상 내지 75% 미만의 RS(%)을 갖는 샘플은 중간 위험군으로 분류하였다. 마지막으로 25% 미만의 RS(%)을 갖는 샘플은 저 위험군으로 분류하였다.
상기 표 4는 위험 군의 분류에 사용된 유전자 및 이의 위험 계수(Hazard Ratio)이다.
Parametric p-value Hazard Ratio UniqueID
0.0011978 0.759 CSK
0.0013091 0.847 CD8A
0.0153688 0.79 CHEK2
도 6에 나타난 바와 같이, 총 520개의 샘플은 저 위험군이 108명 중 24명이 사망하였고(전체 20% 차지함; 사망률 22.2%), 중간 위험군은 299명 중 100명이 사망하였으며(전체 57% 차지함; 사망률 33.5%), 고 위험군은 113명 중 48명이 사망(전체 21.6% 차지함; 사망률 42.5%)하였다.
결과적으로, 위암은 분자적 특징을 이용하여 병기 결정의 가능성을 보여주었다. 위암의 FFPE 샘플에 대한 qPCR를 결과로부터 선정된 유전자 리스트에서 얻은 RS(%)를 통해 위 절제술 후 임상 결과의 예후 예측이 가능함을 확인하였다.

Claims (11)

  1. 대상으로부터 얻은 암세포를 포함하는 생물학적 샘플에서
    CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 결정하는 단계; 및
    상기 단계에서 결정된 RNA 전사체의 발현도에 기초하여 상기 생물학적 샘플의 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 계산하고,
    상기 RS(%)에 따라 예후를 판단하는 단계를 포함하는, 위암으로 진단된 대상에서 예후를 예측하는 방법.
  2. 제1항에 있어서,
    RS 및 RS(%)은 하기 수학식 1 및 2에 따라 계산하는 방법:
    [수학식 1]
    RS=HR1*normLogTransValue1 + HR2*normLogTransValue2 + ... + HRn* normLogTransValuen

    [수학식 2]
    RS(%)=100×(생물학적 샘플의 RS - 모집단의 RS 최소값)/(모집단의 RS 최고값 - 모집단의 RS 최소값)

    상기 식에서,
    HRn는 n번째 RNA 전사체의 위험 계수(hazard ratio)를 나타내고, 상기 HRn이 1 미만인 경우 -1/HRn로 변환하여 사용하며,
    normLogTransValuen는 RNA 전사체의 발현과 관련된 값이고, 이 값은 해당 유전자의 전체 값을 대상으로 하여 중간값을 중심으로 스케일을 변화시킨 값이며,
    상기 모집단은 TNM 병기와 상관 없는 전체 위암을 갖는 일정 수의 집단을 의미하고, 일정 수는 RS 최고값과 최소값을 산출할 수 있는 정도의 임의의 정수이다.
  3. 제1항에 있어서,
    상기 방법이 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는 것인 방법.
  4. 제1항에 있어서,
    상기 방법은 전체 생존율(Overall Survival, OS) 측면에서 생물학적 샘플의 RS(%)가 75% 이상을 갖는 경우 고 위험군, 25% 이상 내지 75% 미만인 경우 중간 위험군, 25% 미만인 경우 저 위험군으로 판단하는 것인 방법.
  5. 대상으로부터 얻은 암세포를 포함하는 생물학적 샘플에서
    CSK, CD8A 및 CHEK2의 RNA 전사체의 발현 수준을 측정하는 단계; 및
    상기 전사체의 발현 증가는 긍정적인 임상 결과 가능성의 증가로 판단하는 단계를 포함하는, 위암으로 진단된 대상에서 예후를 예측하는 방법.
  6. 제5항에 있어서,
    상기 방법이 qPCR 기반 방법인 방법.
  7. 제5항에 있어서,
    상기 임상 결과가 전체 생존율(Overall Survival, OS) 측면에서 표현되는 것인 방법.
  8. 제5항에 있어서,
    상기 방법이 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는 것인 방법.
  9. 위암의 예후 예측을 실행하는 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체에 있어서, 환자로부터 얻은 핵산 시료에서
    CSK, CD8A 및 CHEK2의 RNA 전사체의 발현도를 결정하는 단계; 및
    상기 단계에서 결정된 RNA의 발현도에 기초하여 상기 시료의 위험도 점수(RS, Risk Score) 및 RS 백분율(RS(%))을 계산하고,
    전체 생존율(Overall Survival, OS) 측면에서 상기 RS(%)의 설정치 범위가 75% 이상을 갖는 경우 고 위험군 환자, 25% 이상 내지 75% 미만인 경우 중간 위험군 환자, 25% 미만인 경우 저 위험군 환자로 분류하는 단계를 컴퓨터에 실행시키는 프로그램을 기록한 컴퓨터로 판독가능한 기록 매체.
  10. 제9항에 있어서,
    상기 RS 및 RS(%)은 하기 수학식 1 및 2에 따라 계산하는 것인 기록 매체:
    [수학식 1]
    RS=HR1*normLogTransValue1 + HR2*normLogTransValue2 + ... + HRn* normLogTransValuen

    [수학식 2]
    RS(%)=100×(시료의 RS - 모집단의 RS 최소값)/(모집단의 RS 최고값 - 모집단의 RS 최소값)

    상기 식에서,
    HRn는 n번째 RNA 전사체의 위험 계수(hazard ratio)를 나타내고, 상기 HRn이 1 미만인 경우 -1/HRn로 변환하여 사용하며,
    normLogTransValuen는 RNA 전사체의 발현과 관련된 값이고, 이 값은 해당 유전자의 전체 값을 대상으로 하여 중간값을 중심으로 스케일을 변화시킨 값이며,
    상기 모집단은 TNM 병기와 상관없는 전체 위암을 갖는 일정 수의 집단을 의미하고, 일정 수는 RS 최고값과 최소값을 산출할 수 있는 정도의 임의의 정수이다.
  11. 제9항에 있어서,
    상기 기록 매체는 TNM 병기와 상관없는 전체 위암의 수술에 의한 절제 후 임상 결과를 예측하는 것인 기록 매체.
KR1020130037335A 2013-04-05 2013-04-05 위암에 대한 예후 예측 시스템 KR101504818B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020130037335A KR101504818B1 (ko) 2013-04-05 2013-04-05 위암에 대한 예후 예측 시스템

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020130037335A KR101504818B1 (ko) 2013-04-05 2013-04-05 위암에 대한 예후 예측 시스템

Publications (2)

Publication Number Publication Date
KR20140121523A true KR20140121523A (ko) 2014-10-16
KR101504818B1 KR101504818B1 (ko) 2015-03-24

Family

ID=51992981

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020130037335A KR101504818B1 (ko) 2013-04-05 2013-04-05 위암에 대한 예후 예측 시스템

Country Status (1)

Country Link
KR (1) KR101504818B1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018199589A1 (ko) * 2017-04-24 2018-11-01 (주) 노보믹스 위암의 생물학적 특성에 기반한 군 구분 및 예후 예측 시스템

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070128636A1 (en) 2005-12-05 2007-06-07 Baker Joffre B Predictors Of Patient Response To Treatment With EGFR Inhibitors
CA2700257A1 (en) 2007-10-02 2009-04-09 Hartmut Land Methods and compositions related to synergistic responses to oncogenic mutations
KR20120061010A (ko) * 2010-11-02 2012-06-12 서울대학교산학협력단 위암 진단 방법
WO2012081898A2 (ko) * 2010-12-13 2012-06-21 사회복지법인 삼성생명공익재단 위암의 예후 예측용 마커 및 이를 이용하는 위암의 예후 예측 방법

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018199589A1 (ko) * 2017-04-24 2018-11-01 (주) 노보믹스 위암의 생물학적 특성에 기반한 군 구분 및 예후 예측 시스템
US11365450B2 (en) 2017-04-24 2022-06-21 Novomics Co., Ltd. Group classification and prognosis prediction system based on biological characteristics of gastric cancer

Also Published As

Publication number Publication date
KR101504818B1 (ko) 2015-03-24

Similar Documents

Publication Publication Date Title
JP7042784B2 (ja) 遺伝子発現を用いた前立腺癌の予後を定量化する方法
JP7228896B2 (ja) 乳がん患者の予後の予測方法
CA2585571C (en) Predicting response to chemotherapy using gene expression markers
KR101530689B1 (ko) 직장결장암용 예후 예측
CA2585561C (en) Esr1, pgr, bcl2 and scube2 group score as indicators of breast cancer prognosis and prediction of treatment response
KR101504817B1 (ko) 국소 진행형 위암에 대한 예후 예측 시스템
KR101501826B1 (ko) 위암에 대한 예후 예측 모형의 제조방법
EP3556867A1 (en) Methods to predict clinical outcome of cancer
BRPI0706511A2 (pt) marcadores de expressão genética para prognóstico de cáncer colorretal
US10718030B2 (en) Methods for predicting effectiveness of chemotherapy for a breast cancer patient
US9890430B2 (en) Copy number aberration driven endocrine response gene signature
US20160222461A1 (en) Methods and kits for diagnosing the prognosis of cancer patients
WO2010101916A1 (en) Methods for predicting cancer response to egfr inhibitors
US20120004127A1 (en) Gene expression markers for colorectal cancer prognosis
WO2014089055A1 (en) Tivozanib response prediction
KR101504818B1 (ko) 위암에 대한 예후 예측 시스템

Legal Events

Date Code Title Description
A201 Request for examination
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
FPAY Annual fee payment

Payment date: 20190306

Year of fee payment: 5

FPAY Annual fee payment

Payment date: 20200303

Year of fee payment: 6