KR102360297B1 - Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression - Google Patents

Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression Download PDF

Info

Publication number
KR102360297B1
KR102360297B1 KR1020200063916A KR20200063916A KR102360297B1 KR 102360297 B1 KR102360297 B1 KR 102360297B1 KR 1020200063916 A KR1020200063916 A KR 1020200063916A KR 20200063916 A KR20200063916 A KR 20200063916A KR 102360297 B1 KR102360297 B1 KR 102360297B1
Authority
KR
South Korea
Prior art keywords
expression
mfi
multiple myeloma
risk
equation
Prior art date
Application number
KR1020200063916A
Other languages
Korean (ko)
Other versions
KR20210010311A (en
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 고려대학교 산학협력단
Publication of KR20210010311A publication Critical patent/KR20210010311A/en
Application granted granted Critical
Publication of KR102360297B1 publication Critical patent/KR102360297B1/en

Links

Images

Classifications

    • 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
    • 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/68Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids
    • G01N33/6893Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving proteins, peptides or amino acids related to diseases not provided for elsewhere
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2800/00Detection or diagnosis of diseases
    • G01N2800/52Predicting or monitoring the response to treatment, e.g. for selection of therapy based on assay results in personalised medicine; Prognosis

Abstract

본 발명은 PD-L1 발현을 이용한 다발골수종 예후 예측을 위한 정보 제공 방법 및 다발성 골수종 환자의 자가 조혈모세포이식 여부를 결정하는데 필요한 정보를 제공하는 방법에 관한 것으로서, 상기 방법에 의하면 다발골수종의 예후 예측 및 치료방법의 결정을 보다 합리적으로 수행할 수 있다.The present invention relates to a method for providing information for predicting the prognosis of multiple myeloma using PD-L1 expression and a method for providing information necessary for determining whether autologous hematopoietic stem cell transplantation in a multiple myeloma patient. And the determination of the treatment method can be performed more rationally.

Figure R1020200063916
Figure R1020200063916

Description

PD-L1 발현을 이용한 다발골수종 예후 예측을 위한 정보 제공 방법{Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression}Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression

본 발명은 PD-L1 발현을 이용한 다발골수종 예후 예측을 위한 정보 제공 방법에 관한 것으로서, 더욱 상세하게는 PD-L1의 발현 수준을 측정하여 위험도를 결정하고, 이외의 공지된 위험인자와 결합된 점수 체계에 의해 다발골수종의 예후 예측을 보다 정확하게 수행할 수 있는 정보를 제공하는 방법에 관한 것이다.The present invention relates to a method for providing information for predicting the prognosis of multiple myeloma using PD-L1 expression, and more particularly, to determine the risk by measuring the expression level of PD-L1, and score combined with other known risk factors It relates to a method for providing information that can more accurately predict the prognosis of multiple myeloma by a system.

PD-L1은 대표적인 면역 회피 인자로 다양한 암종에서 예후 인자 및 치료의 타켓으로 연구가 진행되어 왔다. 고형암 중에서는 폐암, 유방암, 대장암 등 많은 암종에서 암세포에서 발현되는 PD-L1이 환자의 예후에 미치는 영향에 대한 연구 결과들이 있으며, 혈액암 중에서는 특히 림프종에서 관련 연구가 많이 진행되어 왔다.PD-L1 is a representative immune evasion factor, and research has been conducted as a prognostic factor and therapeutic target in various carcinomas. Among solid cancers, there are studies on the effect of PD-L1 expressed in cancer cells on the prognosis of patients in many carcinomas such as lung cancer, breast cancer, and colorectal cancer.

하지만 다발골수종(multiple myeloma)에서는 PD-L1 발현이 좋지 않은 예후와 관련이 있다는 기존의 연구 결과는 있지만 예후 인자로 확립된 바는 없으며 이를 이용한 예후 예측 모델은 개발되어 있지 않다. 다발골수종은 골수에서 항체를 생산하는 세포의 한 종류인 형질세포(plasma cell)가 비정상적으로 증식하는 혈액질환으로 특히 뼈를 침윤하는 것이 특징이고 면역장애, 조혈장애 및 신장장애를 일으키는 치명적인 질환이다.However, although there are existing research results that PD-L1 expression is associated with a poor prognosis in multiple myeloma, it has not been established as a prognostic factor, and a prognostic prediction model using this has not been developed. Multiple myeloma is a blood disease in which plasma cells, a type of antibody-producing cells in the bone marrow, proliferate abnormally. It is characterized by infiltrating bones, and is a fatal disease that causes immune disorders, hematopoietic disorders and kidney disorders.

기존 연구들에서는 PD-L1 발현을 순환하는 형질세포 혹은 골수유래 혈액에서 획득한 형질세포에서 유세포분석법으로 측정하거나 혈장 내의 soluble PD-L1을 측정하여 사용하였다. 직접적으로 환자의 골수 샘플에서 면역조직화학법(면역형광염색 포함)으로 형질세포에서 발현되는 PD-L1을 측정하여 예후를 분석한 연구는 제한적이다.In previous studies, PD-L1 expression was measured by flow cytometry in circulating plasma cells or plasma cells obtained from bone marrow-derived blood, or soluble PD-L1 in plasma was measured and used. There are limited studies on prognosis by directly measuring PD-L1 expressed in plasma cells by immunohistochemistry (including immunofluorescence staining) in bone marrow samples from patients.

암 조직에서 정량적 면역형광분석을 통해 특정 단백질(표지자) 발현을 측정하는 방법들은 현재 개발되어 있으며 대표적으로 AQUA (Automated Quantitative Analysis) 방법(method)이 고형암에서 단백질(표지자)의 발현을 측정하는 방법으로 사용되고 있다. 하지만 다발골수종은 골수에서 유래하는 혈액질환으로 암과 주변 조직의 경계를 구분하는 것이 명확하지 않아 AQUA 방법을 적용하기에는 제한이 있다.Methods for measuring the expression of a specific protein (marker) through quantitative immunofluorescence analysis in cancer tissues are currently being developed, and the representative AQUA (Automated Quantitative Analysis) method is a method for measuring the expression of a protein (marker) in solid cancer. is being used However, as multiple myeloma is a blood disease derived from bone marrow, it is not clear to distinguish the boundary between cancer and surrounding tissues, so there is a limit to the application of the AQUA method.

본 발명에서는 AQUA 방법을 참고하여 각각의 형질세포 영역을 기준으로 마스크(mask)를 만들어 MFI의 합을 구하였고 이를 이용하여 점수체계를 만들어 PD-L1의 발현을 측정하였다. 면역조직화학염색을 이용한 분석은 판독자에 따라 결과 해석에 차이가 있을 수 있으나 본 기술은 면역형광염색을 이용하여 이미지 프로그램에 의해 정해진 프로토콜에 따라 분석이 진행되어 일관된 예측이 가능한 장점이 있다.In the present invention, referring to the AQUA method, a mask was made based on each plasma cell region to obtain the sum of the MFIs, and a score system was created using this to measure the expression of PD-L1. Analysis using immunohistochemical staining may have differences in the interpretation of the results depending on the reader, but this technique has the advantage of enabling consistent prediction using immunofluorescence staining as the analysis proceeds according to the protocol determined by the image program.

본 발명자들은 골수 내 형질세포에서 측정한 PD-L1의 MFI 점수체계에 의한 PD-L1 발현 정도가 multivariable Cox regression analysis를 통해 환자의 생존에 영향을 미치는 유의한 인자임을 확인하였고, Contal and O`Quigley method를 이용하여 PD-L1 발현의 optimal cut-off value (7.65)를 결정하였다. 이를 기준으로 cut-off value (7.65)보다 낮은 발현을 저위험군으로 높은 발현을 고위험군으로 구분하여 PD-L1 발현의 고위험군이 저위험군에 비해 예후가 좋지 않음을 확인하였다. 이후 least absolute shrinkage and selection operator (lasso) 변수 선택을 이용한 Cox regression analysis로 nomogram 모형을 수립하였고 이를 기반으로 저위험, 중간위험, 고위험군으로 분류하는 PD-L1 기반 새로운 다발골수종 예후 예측 모델을 개발하였다.The present inventors confirmed that the PD-L1 expression level by the MFI scoring system of PD-L1 measured in plasma cells in the bone marrow is a significant factor affecting the survival of patients through multivariable Cox regression analysis, Contal and O`Quigley method was used to determine the optimal cut-off value (7.65) of PD-L1 expression. Based on this, expression lower than the cut-off value (7.65) was classified as low risk group and high expression level as high risk group. Afterwards, a nomogram model was established by Cox regression analysis using the least absolute shrinkage and selection operator (lasso) variable selection, and based on this, a new PD-L1-based multiple myeloma prognosis prediction model was developed that classifies them into low-risk, medium-risk, and high-risk groups.

새로운 예후 예측 모델은 기존의 예후 예측 모델과 비교하여 전체적으로 높은 예후 판별력을 보였으며, 특히 진단 후 24개월 이전에 사망한 고위험군의 환자를 예측하는데 기존의 예후 예측 모델에 비해 현저한 이점을 보였다. 이를 바탕으로 생존율이 낮을 것으로 예상되는 고위험군 환자들에게는 초기부터 적극적으로 치료를 하여 예후 개선에 도움을 줄 수 있으며, 새로운 예후 예측 모델의 고위험군 환자에서 조혈모세포 이식이 생존율을 높이는 방법이 될 수 있음을 확인하였다.The new prognostic predictive model showed a higher overall prognostic discrimination compared to the existing prognostic predictive model. In particular, it showed a significant advantage over the existing prognostic predictive model in predicting patients in the high-risk group who died before 24 months after diagnosis. Based on this, it is possible to help improve the prognosis by actively treating high-risk patients who are expected to have a low survival rate from the beginning. Confirmed.

Contal and O'Quigley, An application of changepoint methods in studying the effect of age on survival in breast cancer. Comput. Stat. Data An. 30, 253-270 (1999)Contal and O'Quigley, An application of changepoint methods in studying the effect of age on survival in breast cancer. Compute. Stat. Data An. 30, 253-270 (1999)

본 발명자들은 PD-L1의 MFI를 이용한 점수체계를 바탕으로 lasso 변수 선택과 Cox regression analysis에 의해 수립된 PD-L1 발현과 임상인자를 결합한 모형이 예후 판별력이 높음을 확인하였다.The present inventors confirmed that the model combining PD-L1 expression and clinical factors established by lasso variable selection and Cox regression analysis based on the scoring system using the MFI of PD-L1 had high prognostic discrimination.

이에, 본 발명의 목적은 PD-L1의 발현을 기반으로 한 다발성 골수종의 예후 예측을 위한 정보 제공 방법을 제공하는 것이다.Accordingly, it is an object of the present invention to provide an information providing method for predicting the prognosis of multiple myeloma based on the expression of PD-L1.

본 발명의 다른 목적은 다발성 골수종의 자가 조혈모세포이식 여부를 결정하는데 필요한 정보를 제공하는 방법을 제공하는 것이다.Another object of the present invention is to provide a method for providing information necessary for determining whether or not to transplant an autologous hematopoietic stem cell in multiple myeloma.

본 발명은 PD-L1 발현을 이용한 다발골수종(multiple myeloma) 예후 예측을 위한 정보 제공 방법에 관한 것으로, 본 발명에 따른 방법은 PD-L1 발현을 측정하여 다발골수종의 예후 예측의 정확성을 높이고, 기존의 예후 예측 모델에 대비해서도 유의하게 높은 판별력을 나타낸다. The present invention relates to a method for providing information for predicting the prognosis of multiple myeloma using PD-L1 expression, and the method according to the present invention measures PD-L1 expression to increase the accuracy of predicting the prognosis of multiple myeloma. It also shows significantly higher discriminant power compared to the prognostic prediction model of

본 발명자들은 PD-L1의 MFI (mean fluorescence intensity)를 이용한 점수체계를 바탕으로 lasso 변수 선택과 Cox regression analysis에 의해 수립된 PD-L1 발현과 임상인자를 결합한 모형이 예후 판별력이 높음을 확인하였다. The present inventors confirmed that the model combining PD-L1 expression and clinical factors established by lasso variable selection and Cox regression analysis based on the scoring system using the MFI (mean fluorescence intensity) of PD-L1 had high prognostic discrimination.

이하 본 발명을 더욱 자세히 설명하고자 한다.Hereinafter, the present invention will be described in more detail.

본 발명의 일 양태에 따르면, 본 발명은 다음의 단계를 포함하는 다발성 골수종(multiple myeloma)의 예후 예측을 위한 정보 제공 방법을 제공한다:According to one aspect of the present invention, the present invention provides a method for providing information for predicting the prognosis of multiple myeloma, comprising the steps of:

골수 형질세포를 포함하는 시료에서 CD138 및 PD-L1의 발현 정도를 정량적으로 측정하여 하기 식1에 따른 PD-L1의 정규화 MFI(normalized mean fluorescence intensity)를 도출하는 단계:Deriving the normalized mean fluorescence intensity (MFI) of PD-L1 according to Equation 1 by quantitatively measuring the expression levels of CD138 and PD-L1 in a sample containing bone marrow plasma cells:

식1 Expression 1

Normalized MFI = (MFI in plasma cell compartments - Background intensity) / MFI from isotype-matched control; Normalized MFI = (MFI in plasma cell compartments - Background intensity) / MFI from isotype-matched control;

하기 식2에 따른 PD-L1 발현 MFI 측정값을 도출하는 단계:Deriving a PD-L1 expression MFI measurement value according to Equation 2 below:

식2Expression 2

PD-L1 expression MFI = Sum of normalized MFIs in all plasma cell compartments / Total number of plasma cells; 및PD-L1 expression MFI = Sum of normalized MFIs in all plasma cell compartments / Total number of plasma cells; and

상기 식2에 따른 PD-L1 발현 MFI 측정값이 7.65 미만이면 저위험군, 7.65 이상이면 고위험군으로 분류하는 단계.Classifying the PD-L1 expression MFI according to Equation 2 into a low-risk group if the measured value is less than 7.65, and a high-risk group if it is 7.65 or more.

본 발명의 일 구현 예에서, 상기 발현량 측정 단계는 골수 형질세포에서의 PD-L1 발현량을 확인하기 위하여 면역형광염색을 수행하고, 형질세포의 영역 내에서 식1에 따라 PD-L1의 정규화된 MFI를 계산하고, 이를 식2에 따라 모두 더한 후, 이를 상기 형질세포의 영역 수로 나눈 값을 도출한다.In one embodiment of the present invention, in the step of measuring the expression level, immunofluorescence staining is performed to confirm the expression level of PD-L1 in bone marrow plasma cells, and normalization of PD-L1 according to Equation 1 in the plasma cell region The calculated MFI is calculated, added up according to Equation 2, and divided by the number of regions of the plasma cells is derived.

본 명세서상의 "형질세포의 영역 수"는 샘플 내에서 측정을 위한 형질세포의 영역을 설정하였을 때 그 개수를 의미하며, 형질세포의 수와 의미하는 개수는 동일하다.As used herein, the term "number of regions of plasma cells" refers to the number when the regions of plasma cells for measurement are set in a sample, and the number and meaning number of plasma cells are the same.

본 발명의 일구현예에서, 상기 CD138은 형질세포(plasma cell)를 검출하기 위한 마커이다. 상기 PD-L1의 정규화된 MFI는 상기 CD138의 검출로 인한 형질세포 내부 영역에서의 PD-L1의 발현량을 나타낸다. In one embodiment of the present invention, the CD138 is a marker for detecting plasma cells. The normalized MFI of PD-L1 indicates the expression level of PD-L1 in the plasma cell internal region due to the detection of CD138.

본 발명의 다른 일 양태에 따르면, 본 발명은 다음의 단계를 포함하는 다발성 골수종의 예후 예측을 위한 정보 제공 방법을 제공한다:According to another aspect of the present invention, the present invention provides a method for providing information for predicting the prognosis of multiple myeloma, comprising the steps of:

Figure 112020053853123-pat00001
Figure 112020053853123-pat00001

i) 상기 식1 및 식2에 따른 PD-L1 발현 MFI 측정값(0~35), ii) 정상치 상한(422 IU/L) 초과의 혈중 LDH(serum lactate dehydrogenase) 여부, iii) high risk cytogenetics (karyotype에서 t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality 혹은 FISH에서 IGH/FGFR3 rearrangement, IGH/MAF rearrangement, TP53 mutation) 여부, iv) 연령 70세 이상 여부를 기준으로 상기 노모그램(nomogram)에서 상응하는 점수(points, 0~100)를 가산하는 점수 배점 단계; 및 i) PD-L1 expression MFI measured value (0 to 35) according to Equations 1 and 2 above, ii) Whether serum lactate dehydrogenase (LDH) in blood exceeding the upper limit of normal (422 IU/L), iii) high risk cytogenetics ( t(4;14), t(14;16), 17p deletion in karyotype, chromosome 1 abnormality, or IGH/FGFR3 rearrangement in FISH, IGH/MAF rearrangement, TP53 mutation), iv) based on age 70 years or older a score scoring step of adding corresponding scores (points, 0 to 100) in the nomogram; and

상기 각 가산된 점수를 합산하여 총점(total points)을 구한 후 상기 노모그램에서 상응하는 1년, 2년, 또는 4년 전체 생존기간에 대한 확률을 도출하는 단계.Deriving a probability for a corresponding 1-year, 2-year, or 4-year overall survival period from the nomogram after summing each added score to obtain a total point.

본 발명의 일구현예에 있어서, 상기 1년 전체 생존기간에 대한 확률은 0.9 내지 0.6이다. 본 발명의 다른 구현예에 있어서, 상기 2년 전체 생존기간에 대한 확률은 0.9 내지 0.4 이다. 본 발명의 또 다른 구현예에 있어서, 상기 4년 전체 생존기간에 대한 확률은 0.8 내지 0.2 이다.In one embodiment of the present invention, the probability for the 1-year total survival period is 0.9 to 0.6. In another embodiment of the present invention, the probability for the 2-year overall survival period is 0.9 to 0.4. In another embodiment of the present invention, the probability for the 4-year overall survival period is 0.8 to 0.2.

본 발명의 일구현예에 있어서, 상기 혈중 LDH의 정상치 상한 값은 검사 기관, 방법에 따라 차이가 있을 수 있는 수치로, 본 발명자가 소속된 기관의 참고치(238-422 IU/L)의 상한 값인 422 IU/L 를 기준으로 적용하였다.In one embodiment of the present invention, the upper limit of the normal value of LDH in the blood is a value that may differ depending on the test institution and method, and is the upper limit of the reference value (238-422 IU/L) of the institution to which the present inventor belongs. 422 IU/L was applied as a standard.

본 발명의 다른 일 양태에 따르면, 본 발명은 다음의 단계를 포함하는 다발성 골수종의 자가 조혈모세포이식 여부를 결정하는데 필요한 정보를 제공하는 방법을 제공한다:According to another aspect of the present invention, there is provided a method for providing information necessary for determining whether to transplant an autologous hematopoietic stem cell of multiple myeloma, comprising the steps of:

Figure 112020053853123-pat00002
Figure 112020053853123-pat00002

i) 상기 식1 및 식2에 따른 PD-L1 발현 MFI 측정값(0~35), ii) 정상치 상한(422 IU/L) 초과의 혈중 LDH(serum lactate dehydrogenase) 여부, iii) high risk cytogenetics (karyotype에서 t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, FISH에서 IGH/FGFR3 rearrangement, IGH/MAF rearrangement, 및 TP53 mutation으로 이루어진 군으로부터 1종 이상 양성인 경우) 여부, iv) 연령 70세 이상 여부를 기준으로 상기 nomogram에서 상응하는 점수(points, 0~100)를 가산하는 점수 배점 단계; i) PD-L1 expression MFI measured value (0 to 35) according to Equations 1 and 2 above, ii) Whether serum lactate dehydrogenase (LDH) in blood exceeding the upper limit of normal (422 IU/L), iii) high risk cytogenetics ( karyotype t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, IGH/FGFR3 rearrangement in FISH, IGH/MAF rearrangement, and at least one positive from the group consisting of TP53 mutation) iv) a score scoring step of adding corresponding scores (points, 0-100) in the nomogram based on whether the age is 70 years or older;

상기 각 가산된 점수를 합산하여 총점(total points)을 구한 후 50점 미만이면 저위험군, 50 내지 100점이면 중간위험군, 100점을 초과하면 고위험군으로 분류하는 분류 단계; 및A classification step of adding up each added score to obtain a total point and classifying a low-risk group if less than 50 points, a medium-risk group if 50 to 100 points, and a high-risk group if more than 100 points; and

총점이 100점을 초과하는 고위험군으로 분류된 환자에게 자가조혈모세포 이식시 생존율이 향상된다는 정보를 제공하는 단계.Providing information that survival rate is improved by autologous hematopoietic stem cell transplantation to patients classified as high-risk group with a total score exceeding 100 points.

본 발명은 PD-L1 발현을 이용한 다발골수종 예후 예측을 위한 정보 제공 방법 및 다발성 골수종 환자의 자가 조혈모세포이식 여부를 결정하는데 필요한 정보를 제공하는 방법에 관한 것으로서, 상기 방법에 의하면 다발골수종의 예후 예측 및 치료방법의 결정을 보다 합리적으로 수행할 수 있다.The present invention relates to a method for providing information for predicting the prognosis of multiple myeloma using PD-L1 expression and a method for providing information necessary for determining whether autologous hematopoietic stem cell transplantation in a multiple myeloma patient. And the determination of the treatment method can be performed more rationally.

도 1은 다발골수종 환자의 골수 내 형질세포에서 PD-L1이 발현되는 것을 확인한 공초점현미경 사진이다.
도 2는 PD-L1의 발현을 측정하기 위해 형질세포의 영역을 정하고 MFI를 측정하는 방법을 나타낸 사진이다.
도 3a는 PD-L1 MFI(mean fluorescence intensity) 측정값 및 위험군에 따른 PD-L1 발현 정도를 보여주기 위해 대표적으로 선정한 사진이다.
도 3b는 PD-L1 MFI(mean fluorescence intensity) 측정값에 따른 PD-L1 발현 분포를 보여주는 그래프이다.
도 4a는 전체 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 전체생존율(overall survival)을 나타낸 그래프이다.
도 4b는 autologous stem cell transplantation (ASCT)를 받은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 전체생존율을 나타낸 그래프이다.
도 4c는 ASCT를 받지 않은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 전체생존율을 나타낸 그래프이다.
도 4d는 VTD (bortezomib, thalidomide 및 dexamethasone), TD (thalidomide 및 dexamethasone), RD (lenalidomide 및 dexamethasone) 요법으로 치료를 받은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 전체생존율을 나타낸 그래프이다.
도 4e는 VMP (bortezomib(Velcade) plus melphalan 및 prednisone) 요법으로 치료를 받은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 전체생존율을 나타낸 그래프이다.
도 4f는 ASCT를 받은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 무진행생존율(progression-free survival)을 나타낸 그래프이다.
도 4g는 ASCT를 받지 않은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 무진행생존율을 나타낸 그래프이다.
도 4h는 VTD, TD, RD 요법으로 치료를 받은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 무진행생존율을 나타낸 그래프이다.
도 4i는 VMP 요법으로 치료를 받은 환자군에서 PD-L1 MFI 측정값에 의한 PD-L1 발현에 기반한 위험군에 따른 무진행생존율을 나타낸 그래프이다.
도 5는 다변수 콕스회귀분석(multivariable Cox regression)을 통해 얻어진 유의한 위험인자를 나타낸 표이다. PD-L1 발현에 대한 위험인자는 PD-L1 MFI 측정값에 대한 optimal cut-off value를 이용한 이분형 변수로 정의되었다.
도 6a는 lasso 변수 선택과 Cox regression analysis에 의해 수립된 PD-L1 발현과 임상인자가 결합된 nomogram 모형이다. PD-L1 발현에 대한 위험인자는 PD-L1 MFI 측정값을 이용한 연속형 변수로 정의되었다.
도 6b는 수립된 nomogram 모형에 대한 calibration 그래프이다.
도 6c는 수립된 nomogram의 total score를 기준으로 세 그룹으로 나누어 카플란-메이어(Kaplan-Meier) 생존곡선으로 각 그룹의 생존을 비교한 그래프이다.
도 6d는 수립된 nomogram 모형을 평가하기 위한 time-dependent AUC 분석 그래프이다.
도 6e는 수립된 nomogram 모형을 test 하기 위해 test cohort에 적용하여 calibration한 그래프이다.
도 6f는 수립된 nomogram 모형을 test 하기 위해 test cohort에 적용하여 카플란-메이어(Kaplan-Meier) 생존곡선으로 각 그룹의 생존을 비교한 그래프이다.
도 6g는 수립된 nomogram 모형을 test 하기 위해 test cohort에 적용하여 time-dependent AUC를 분석한 그래프이다.
도 7a는 nomogram total score를 기준으로 세 그룹으로 구분한 새로운 예후 예측 모형에서 카플란-메이어(Kaplan-Meier) 생존곡선으로 각 그룹의 생존을 비교한 그래프이다.
도 7b는 기존의 R-ISS를 기준으로 카플란-메이어(Kaplan-Meier) 생존곡선으로 각 그룹의 생존을 비교한 그래프이다.
도 7c는 새로운 예후 예측 모형과 R-ISS의 성능을 time-dependent AUC 분석을 이용하여 비교한 그래프이다.
도 8a 내지 도 8f는 새로운 예후 예측 모형에서 위험군별 치료법에 대한 전체생존율과 무진행생존율을 나타낸 그래프이다.
1 is a confocal micrograph confirming the expression of PD-L1 in plasma cells in the bone marrow of a patient with multiple myeloma.
2 is a photograph showing a method of determining the region of plasma cells and measuring MFI in order to measure the expression of PD-L1.
3A is a photograph representatively selected to show PD-L1 mean fluorescence intensity (MFI) measurement values and PD-L1 expression levels according to risk groups.
3B is a graph showing the PD-L1 expression distribution according to PD-L1 mean fluorescence intensity (MFI) measurement values.
Figure 4a is a graph showing the overall survival (overall survival) according to the risk group based on the PD-L1 expression by the PD-L1 MFI measurement in the entire patient group.
Figure 4b is a graph showing the overall survival rate according to the risk group based on PD-L1 expression by PD-L1 MFI measurement in the patient group that received autologous stem cell transplantation (ASCT).
Figure 4c is a graph showing the overall survival rate according to the risk group based on the PD-L1 expression by the PD-L1 MFI measurement value in the patient group that did not receive ASCT.
Figure 4d shows the overall group according to risk group based on PD-L1 expression by PD-L1 MFI measurement values in patients treated with VTD (bortezomib, thalidomide and dexamethasone), TD (thalidomide and dexamethasone), and RD (lenalidomide and dexamethasone) regimens. It is a graph showing the survival rate.
Figure 4e is a graph showing the overall survival rate according to the risk group based on PD-L1 expression by the PD-L1 MFI measurement value in the patient group treated with VMP (bortezomib (Velcade) plus melphalan and prednisone) therapy.
Figure 4f is a graph showing the progression-free survival (progression-free survival) according to the risk group based on the PD-L1 expression by the PD-L1 MFI measurement in the ASCT patient group.
Figure 4g is a graph showing the progression-free survival rate according to the risk group based on the PD-L1 expression by the PD-L1 MFI measurement in the patient group that did not receive ASCT.
FIG. 4H is a graph showing progression-free survival rates according to risk groups based on PD-L1 expression by PD-L1 MFI measurement values in patient groups treated with VTD, TD, and RD therapy.
4I is a graph showing progression-free survival rates according to risk groups based on PD-L1 expression by PD-L1 MFI measurement values in a patient group treated with VMP therapy.
5 is a table showing significant risk factors obtained through multivariable Cox regression. Risk factors for PD-L1 expression were defined as dichotomous variables using optimal cut-off values for PD-L1 MFI measurements.
6a is a nomogram model in which PD-L1 expression and clinical factors established by lasso variable selection and Cox regression analysis are combined. Risk factors for PD-L1 expression were defined as continuous variables using PD-L1 MFI measurements.
6B is a calibration graph for the established nomogram model.
6c is a graph comparing the survival of each group with a Kaplan-Meier survival curve divided into three groups based on the established total score of the nomogram.
6D is a time-dependent AUC analysis graph for evaluating the established nomogram model.
Figure 6e is a graph calibrated by applying to the test cohort to test the established nomogram model.
6f is a graph comparing the survival of each group with a Kaplan-Meier survival curve applied to a test cohort to test the established nomogram model.
6g is a graph analyzing the time-dependent AUC by applying the test cohort to test the established nomogram model.
7A is a graph comparing the survival of each group with a Kaplan-Meier survival curve in a new prognostic prediction model divided into three groups based on a nomogram total score.
7b is a graph comparing the survival of each group with a Kaplan-Meier survival curve based on the existing R-ISS.
7c is a graph comparing the performance of the new prognostic prediction model and R-ISS using time-dependent AUC analysis.
8A to 8F are graphs showing the overall survival rate and progression-free survival rate for each risk group treatment in the new prognostic prediction model.

이하, 본 발명을 하기의 실시예에 의하여 더욱 상세히 설명한다. 그러나 이들 실시예는 본 발명을 예시하기 위한 것일 뿐이며, 본 발명의 범위가 이들 실시예에 의하여 한정되는 것은 아니다.Hereinafter, the present invention will be described in more detail with reference to the following examples. However, these examples are only for illustrating the present invention, and the scope of the present invention is not limited by these examples.

본 명세서 전체에 걸쳐, 특정 물질의 농도를 나타내기 위하여 사용되는 "%"는 별도의 언급이 없는 경우, 고체/고체는 (중량/중량)%, 고체/액체는 (중량/부피)%, 그리고 액체/액체는 (부피/부피)%이다.Throughout this specification, "%" used to indicate the concentration of a specific substance is (weight/weight)% solid/solid, (weight/volume)%, and (weight/volume)% for solid/solid, and Liquid/liquid is (vol/vol) %.

실시예Example 1: 골수 샘플의 준비 1: Preparation of bone marrow samples

실험에 사용된 환자의 샘플은 2011년 1월부터 2019년 10월까지 고려대학교 안암병원 혈액종양내과에서 골수검사 후 다발골수종으로 진단된 환자 126명(후향적 83명, 전향적 43명)으로부터 얻었다. 모든 샘플은 진단 시 사용된 골수 파라핀 셀 블록(cell block)으로 IRB 승인 후 프로토콜에 따라 준비하였으며 각 샘플당 적어도 2개 이상의 절편을 제작하여 확인하였다.The patient samples used in the experiment were obtained from 126 patients (retrospective 83 patients, prospective 43 patients) diagnosed with multiple myeloma after bone marrow examination at the Department of Hematology and Oncology, Korea University Anam Hospital from January 2011 to October 2019. . All samples were bone marrow paraffin cell blocks used for diagnosis, prepared according to the protocol after approval by the IRB, and confirmed by making at least two sections for each sample.

실시예Example 2: 2: 면역형광염색을immunofluorescence staining 통한 PD-L1 발현의 확인 Confirmation of PD-L1 expression through

준비된 골수 파라핀 셀 블록에서 4-5um 두께의 절편을 잘라 슬라이드를 제작하였고 60℃에서 1시간 동안 건조한 후 자일렌(xylene)으로 15분 동안 3회, 100%, 95%, 70% 에탄올(ethanol)로 각각 5분 동안 2회 처리한 후 흐르는 물로 세척하여 탈파라핀화를 시행하였다. 이후 압력 쿠커를 이용하여 pH 6.0의 버퍼(sodium citrate buffer) 내에서 10분 동안 끓인 후 흐르는 물로 냉각시켜 항원성을 복원하였고, 0.5% 트리톤(triton) X-100를 이용하여 세포막 투과화 과정을 진행하였다.Slides were prepared by cutting 4-5 μm thick sections from the prepared bone marrow paraffin cell block. After drying at 60 ° C for 1 hour, xylene was used for 15 minutes 3 times, 100%, 95%, 70% ethanol. After treatment twice for 5 minutes each with running water, deparaffinization was performed by washing with running water. After boiling for 10 minutes in a pH 6.0 buffer (sodium citrate buffer) using a pressure cooker, antigenicity was restored by cooling with running water, and the cell membrane permeabilization process was performed using 0.5% triton X-100. did

이후 5% 세럼(normal donkey serum)으로 실온에서 1시간 동안 블로킹(blocking) 후 PBS로 세척하고, 일차항체(goat CD138 Ab (R&D Systems), mouse anti-PD-L1 mAb [ABM4E54] (Abcam))를 사용하여 4℃에서 오버나잇(overnight)으로 배양하였다. 이후 PBS로 5분 동안 3회 세척하고 이차항체(alexa flour 488 conjugated donkey anti-goat IgG (Invitrogen), alexa flour 647 conjugated donkey anti-mouse IgG (Invitrogen))를 사용하여 실온에서 1시간 동안 배양 후 PBS로 5분 동안 3회 세척하고 DAPI가 포함된 마운트제(mountant, ProLong™ Diamond Antifade Mountant with DAPI (Invitrogen))를 이용하여 처리하였다.After blocking with 5% serum (normal donkey serum) at room temperature for 1 hour, washing with PBS, primary antibody (goat CD138 Ab (R&D Systems), mouse anti-PD-L1 mAb [ABM4E54] (Abcam)) was used and incubated overnight at 4°C. After washing with PBS 3 times for 5 minutes, using secondary antibodies (alexa flour 488 conjugated donkey anti-goat IgG (Invitrogen), alexa flour 647 conjugated donkey anti-mouse IgG (Invitrogen)) for 1 hour at room temperature, PBS was washed three times for 5 minutes and treated with a mount agent containing DAPI (mountant, ProLong™ Diamond Antifade Mountant with DAPI (Invitrogen)).

도 1에서 확인할 수 있듯이, 공초점현미경을 이용하여 다발골수종 환자의 골수 내 형질세포에서 PD-L1이 발현되는 것을 확인하였다. CD138, PD-L1의 이미지(8 bit raw)를 이미지 분석 프로그램(CellesteTM Image Analysis Software (Invitrogen))을 이용하여 분석하였다.As can be seen in FIG. 1 , it was confirmed that PD-L1 was expressed in plasma cells in the bone marrow of patients with multiple myeloma using a confocal microscope. CD138 and PD-L1 images (8 bit raw) were analyzed using an image analysis program (Celleste TM Image Analysis Software (Invitrogen)).

샘플 당 3개 이상(평균 5개)의 필드에서 분석을 진행하였다. 도 2와 같이 CD138 발현을 기준으로 형질세포의 영역을 정하고 이미지 분석 프로그램(CellesteTM Image Analysis Software (Invitrogen))을 이용하여 각각의 형질세포 영역 내 PD-L1의 정규화된(normalized) MFI(mean fluorescence intensity)를 계산하여 더한 후 전체 형질세포의 영역 수로 나눈 값을 기준으로 반정량화된 PD-L1 발현MFI(PD-L1 expression MFI)을 측정하였다.Analysis was performed in 3 or more fields per sample (average of 5). As shown in FIG. 2, the region of plasma cells was determined based on CD138 expression and using an image analysis program (Celleste TM Image Analysis Software (Invitrogen)). After calculating the normalized mean fluorescence intensity (MFI) of PD-L1 in each plasma cell region, the semi-quantified PD-L1 expression MFI (PD-L1 expression) is based on the value divided by the total number of plasma cell regions. MFI) was measured.

상기 정규화된 MFI(normalized mean fluorescence intensity)를 계산하는 수식은 식1과 같고, 이로부터 PD-L1 발현 MFI를 계산하는 수식은 식2와 같았다.The formula for calculating the normalized mean fluorescence intensity (MFI) is the same as Equation 1, and the formula for calculating the PD-L1 expression MFI therefrom is the same as Equation 2.

식1 Expression 1

Normalized MFI = (MFI in plasma cell compartments - Background intensity) / MFI from isotype-matched controlNormalized MFI = (MFI in plasma cell compartments - Background intensity) / MFI from isotype-matched control

식2Expression 2

PD-L1 expression (MFI) = Sum of normalized MFIs in all plasma cell compartments / Total number of plasma cellsPD-L1 expression (MFI) = Sum of normalized MFIs in all plasma cell compartments / Total number of plasma cells

상기 식1의 MFI from isotype-matched control은 동형대조군의 MFI 관찰값으로서 PD-L1을 타겟팅하지 않으나, 상기 PD-L1의 MFI를 측정하기 위하여 사용된 항체와 동일한 isotype 및 subclass를 가진 항체를 이용한 대조군의 MFI를 말한다.The MFI from isotype-matched control of Formula 1 does not target PD-L1 as the MFI observed value of the isotype control group, but a control using an antibody having the same isotype and subclass as the antibody used to measure the MFI of PD-L1 refers to the MFI of

도 3a 및 3b에서 확인할 수 있듯이, PD-L1 발현에 대한 MFI 측정값이 7.65 미만인 경우 저위험군, 7.65 이상인 경우 고위험군인 것으로 분류하였다(Reference: Contal and O'Quigley, An application of changepoint methods in studying the effect of age on survival in breast cancer. Comput . Stat. Data An. 30, 253-270 (1999)).As can be seen in Figures 3a and 3b, when the MFI measurement value for PD-L1 expression is less than 7.65, it was classified as a low-risk group, and when it was 7.65 or more, it was classified as a high-risk group (Reference: Contal and O'Quigley, An application of changepoint methods in studying the effect of age on survival in breast cancer. Comput . Stat. Data An. 30, 253-270 (1999)).

실시예Example 3: PD-L1 3: PD-L1 발현에 대한 생존 분석Survival analysis for expression

전체 환자를 i) 자가조혈모세포이식(ASCT)을 받은 환자와 ii) 받지 않은 환자, iii) frontline으로 면역조절제인 탈리도마이드(thalidomide) 및 레날리도마이드(lenalidomide)가 포함된 VTD (bortezomib, thalidomide 및 dexamethasone), TD (thalidomide 및 dexamethasone), RD (lenalidomide 및 dexamethasone) 요법으로 치료 받은 환자와 iv) 면역조절제가 포함되지 않은 VMP (bortezomib(Velcade) plus melphalan 및 prednisone) 요법으로 치료를 받은 환자로 나누어 전체생존율(overall survival)과 무진행생존율(progression-free survival)에 대한 분석을 시행하여 도 4a 내지 4i로 나타내었다. All patients were treated i) with and without autologous hematopoietic stem cell transplantation (ASCT), iii) VTD (bortezomib, thalidomide, and dexamethasone), TD (thalidomide and dexamethasone), and RD (lenalidomide and dexamethasone), and iv) VMP (bortezomib (Velcade) plus melphalan and prednisone) without immunomodulatory agents. Analysis of overall survival and progression-free survival was performed and shown in FIGS. 4A to 4I.

상기 도 4a 내지 도 4i의 생존곡선은 카플란-메이어(Kaplan-Meier) 생존분석을 이용하여 작성하였고 로그순위 검정(log-rank test)으로 두 군의 생존율을 비교하였다.The survival curves of FIGS. 4A to 4I were prepared using Kaplan-Meier survival analysis, and the survival rates of the two groups were compared by a log-rank test.

도 4a 내지 도 4i에서 확인할 수 있듯이, 분석 결과 ASCT를 받은 그룹에서는 PD-L1 발현의 고위험군과 저위험군 간에 생존에 차이가 없었지만(p>0.05), 이외 전체 환자, ASCT 받은 그룹, 면역조절제 치료 그룹, 면역조절제가 포함되지 않은 치료 그룹 모두에서는 고위험군에서 전체생존율과 무진행생존율이 좋지 않았다(p<0.05).As can be seen in FIGS. 4A to 4I , as a result of the analysis, in the group that received ASCT, there was no difference in survival between the high-risk group and the low-risk group of PD-L1 expression (p>0.05), but all other patients, the group that received ASCT, and the group treated with immunomodulatory agents , the overall survival and progression-free survival were poor in the high-risk group in all treatment groups that did not include immunomodulatory agents (p<0.05).

실시예Example 4: PD-L1 기반 새로운 예후 예측 모델 4: New prognostic prediction model based on PD-L1-

도 5 및 6a에서 나타나는 바와 같이 lasso (least absolute shrinkage and selection operator) 변수 선택과 Cox regression analysis에 의해 PD-L1 발현과 임상인자가 결합된 nomogram 모형을 수립하였고 이 과정에서 과적합을 피하기 위해 leave-one-out cross validation (LOOCV)에 의해 lasso의 lambda 값을 결정하였다. As shown in FIGS. 5 and 6a, a nomogram model combining PD-L1 expression and clinical factors was established by lasso (least absolute shrinkage and selection operator) variable selection and Cox regression analysis, and in this process, to avoid overfitting, leave- The lambda value of lasso was determined by one-out cross validation (LOOCV).

구체적으로 상기 임상인자는 70세 이상의 연령, 2 이상의 ECOG Performance Status, 3 mg/dL 이상의 혈청 M-단백질, non-IgG isotype, 100 이상의 serum Free Light Chain ratio, 60% 이상의 BM plasma cell, 5.5 mg/L 이상의 b2-microglobulin, 3.5 g/dL 이상의 알부민, 정상치 상한(422 IU/L) 초과의 LDH, 세포유전학적 고위험군(chromosome 1 abnormality), PD-L1 고위험군이 분석되었다.Specifically, the clinical factors are age 70 or older, ECOG Performance Status of 2 or more, serum M-protein of 3 mg/dL or more, non-IgG isotype, serum Free Light Chain ratio of 100 or more, BM plasma cell of 60% or more, 5.5 mg/ b2-microglobulin ≥ L, albumin ≥ 3.5 g/dL, LDH above the upper limit of normal (422 IU/L), cytogenetic high risk group (chromosome 1 abnormality), and PD-L1 high risk group were analyzed.

단변수 분석에서는 정상치 상한 초과의 LDH, 세포유전학적 고위험군, 및 PD-L1 고위험군이라는 임상인자가 유의한 것으로 나타났고, Cox regression analysis을 활용한 다변수 분석에서는 70세 이상의 고령, 정상치 상한 초과의 LDH, 세포유전학적 고위험군(karyotype에서 t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, FISH에서 IGH/FGFR3 rearrangement, IGH/MAF rearrangement, 및 TP53 mutation으로 이루어진 군으로부터 1종 이상 양성인 경우), 및 PD-L1 고위험군의 임상인자가 유의한 것으로 분석되었다(도 5).In univariate analysis, clinical factors such as LDH exceeding the upper limit of normal, cytogenetic high-risk group, and PD-L1 high-risk group were significant. , 1 type from the group consisting of high cytogenetic risk group (karyotype: t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, IGH/FGFR3 rearrangement in FISH, IGH/MAF rearrangement, and TP53 mutation abnormally positive), and the clinical factors of the PD-L1 high-risk group were analyzed to be significant ( FIG. 5 ).

상기 분석으로 선택된 임상인자 및 PD-L1 고발현 여부의 영향력을 점수화한 points, 각전수의 누적합계인 total points, 및 1년 전체생존율, 2년 전체생존율, 4년 전체생존율에 대한 확률을 계산한 노모그램을 작성하였다(도 6a).Points scoring the influence of the clinical factors selected through the above analysis and high PD-L1 expression, total points that are the cumulative sum of each number, and the probability of 1-year overall survival rate, 2-year overall survival rate, and 4-year overall survival rate were calculated. A nomogram was prepared (Fig. 6a).

다음으로, 도 6b에서와 같이 1년, 2년, 4년 전체생존율에 대해 1000회 bootstrap을 이용하여 calibration을 시행하였고 calibration이 잘 되는 것을 확인하였다. 도 6c에서와 같이 수립된 nomogram의 total score를 기준으로 50미만을 저위험군, 50 내지 100을 중간위험군, 100 초과를 고위험군으로 하는 세 그룹으로 나누어 카플란-메이어(Kaplan-Meier) 생존곡선으로 각 그룹의 생존을 비교하였고, 로그 순위 분석 결과, P <0.001로 세 군의 전체생존율이 통계적으로 유의하게 구분되는 것을 확인하였다. Next, as in FIG. 6b , calibration was performed using bootstrap 1000 times for 1-year, 2-year, and 4-year overall survival rates, and it was confirmed that the calibration was well performed. Based on the total score of the nomogram established as shown in FIG. 6c, each group was divided into three groups: a low-risk group of less than 50, a medium-risk group of 50 to 100, and a high-risk group of more than 100 based on a Kaplan-Meier survival curve. were compared, and as a result of log-rank analysis, it was confirmed that the overall survival rates of the three groups were statistically significantly differentiated with P < 0.001.

도 6d에서와 같이 수립된 nomogram 모형을 평가하기 위해 1000회 bootstrap을 이용하여 time-dependent AUC 분석을 시행하였고 1-4년 AUC 값이 0.6 이상, 1년, 2년 AUC는 0.8 이상으로 확인되어 전체적으로 모형 성능이 양호하며 특히 1년, 2년 생존 평가에 우수함을 확인하였다(0.9-1.0 = excellent; 0.8-0.9 = very good; 0.7-0.8 = good; 0.6-0.7= sufficient; 0.5-0.6 = bad; < 0.5 = not useful (Muller et al., 2005)). To evaluate the nomogram model established as in Fig. 6d, time-dependent AUC analysis was performed using 1000 times bootstrap, and AUC values for 1-4 years were 0.6 or more, and AUC values for 1 year and 2 years were confirmed to be 0.8 or more. It was confirmed that the model performance was good, especially in the evaluation of 1-year and 2-year survival (0.9-1.0 = excellent; 0.8-0.9 = very good; 0.7-0.8 = good; 0.6-0.7 = sufficient; 0.5-0.6 = bad; < 0.5 = not useful (Muller et al., 2005)).

도 6e에서와 같이 수립된 nomogram 모형을 prospective (test) cohort에 적용하여 test 하였고 nomogram에 의해 예측된 생존 확률과 관찰된 생존 확률이 서로 calibration이 잘 되는 것을 확인하였다. 도 6f에서와 같이 수립된 nomogram 모형을 prospective (test) cohort에 적용하여 카플란-메이어(Kaplan-Meier) 생존곡선으로 각 그룹의 생존을 비교하였고 P = 0.029로 세 군의 전체생존율이 통계적으로 유의하게 구분되는 것을 확인하였다. 도 6f에서와 같이 수립된 nomogram 모형의 성능을 평가하기 위해 prospective (test) cohort에 적용하여 time-dependent AUC를 분석하였고 6개월 AUC 0.724, 12개월 AUC 0.740으로 측정되어 test set에서도 본 발명의 nomogram이 다발성 골수종의 예후 판별에 우수함을 확인하였다.The nomogram model established as in FIG. 6e was tested by applying it to the prospective (test) cohort, and it was confirmed that the survival probability predicted by the nomogram and the observed survival probability were well calibrated. The nomogram model established as in FIG. 6f was applied to the prospective (test) cohort, and the survival of each group was compared with the Kaplan-Meier survival curve, and the overall survival rate of the three groups was statistically significant with P = 0.029. distinction was confirmed. To evaluate the performance of the nomogram model established as in FIG. 6f, time-dependent AUC was analyzed by applying it to the prospective (test) cohort, and it was measured as AUC 0.724 for 6 months and AUC 0.740 for 12 months. It was confirmed that it was excellent in prognostic discrimination of multiple myeloma.

상기 cohort 환자의 샘플은 2011년 1월부터 2019년 10월까지 고려대학교 안암병원 혈액종양내과에서 골수검사 후 다발골수종으로 진단된 환자 126명(후향적 83명, 전향적 43명)의 결과로부터 얻었다. 후향적 83명은 2011년 1월부터 2018년 4월까지, 전향적 43명은 2018년 5월부터 2019년 10월까지 진단된 환자들의 결과였다.The cohort patient sample was obtained from the results of 126 patients (retrospective 83 patients, prospective 43 patients) diagnosed with multiple myeloma after bone marrow examination at the Department of Hematology and Oncology, Korea University Anam Hospital from January 2011 to October 2019. . The retrospective 83 patients were from January 2011 to April 2018, and the prospective 43 patients were the results of patients diagnosed from May 2018 to October 2019.

실시예Example 5: 기존의 예후 예측 모델과의 비교 분석 5: Comparative analysis with existing prognostic prediction models

새로운 본 발명의 다발성 골수종의 예후 예측 모델과 종래 예후 예측 모델인 R-ISS(Revised multiple myeloma International Staging System)간의 예후 판별력과 일치도를 비교하였다.The prognostic discriminatory power and concordance between the new multiple myeloma prognostic prediction model of the present invention and the Revised multiple myeloma International Staging System (R-ISS), which is a conventional prognostic prediction model, were compared.

도 7a 및 7b에서 확인할 수 있듯이, R-ISS 대비 본 발명의 새로운 예후 예측 모델로부터 도출할 수 있는 각 군별 생존곡선의 통계적 유의성을 확인하였다. 또한 표 1에서와 같이 Harrell’s c-index (Reference: Multivariable prognostic models: Issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat. Med . 15, 361-387 (1996).)를 이용하여 새로운 예후 예측 모델이 기존의 R-ISS 모델보다 예후 예측에 높은 판별력을 보이는 것을 확인하였고, 도 7c에서와 같이 time-dependent AUC 분석에서 본 발명의 새로운 예후 예측 모형이 R-ISS에 비해 높은 판별력을 보이는 것을 확인하였다.As can be seen in FIGS. 7A and 7B , the statistical significance of the survival curves for each group that can be derived from the new prognostic prediction model of the present invention compared to the R-ISS was confirmed. Also, as in Table 1, using Harrell's c-index (Reference: Multivariable prognostic models: Issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat. Med . 15, 361-387 (1996).) It was confirmed that the new prognostic prediction model showed higher discriminant power in prognostic prediction than the existing R-ISS model, and in the time-dependent AUC analysis as shown in FIG. confirmed what was visible.

구분division Harrell's CHarrell's C 95% CI95% CI DifferenceDifference 95% CI95% CI PP 본 발명의 예후 예측모델Prognosis prediction model of the present invention 0.7750.775 0.717-0.8330.717-0.833 0.1250.125 0.037-0.2130.037-0.213 0.0050.005 R-ISSR-ISS 0.6500.650 0.568-0.7320.568-0.732 -- -- --

CI, Confidence Interval; R-ISS, Revised International Staging SystemCI, Confidence Interval; R-ISS, Revised International Staging System

표 2에서 확인할 수 있듯이, R-ISS의 stage I 환자의 21.4%는 새로운 모델에서 중간위험군으로, R-ISS의 stage II 환자의 32.9%, 32.9%는 새로운 모델에서 각각 저위험군, 고위험군으로, R-ISS stage III 환자의 48.5%는 새로운 모델의 중간위험군으로 재분류되었고 이 차이는 통계적으로 유의하였다. 재분류에 의해 계산된 NRI (Net Reclassification Improvement)는 0.337으로 새로운 모형이 R-ISS에 비해 민감도와 특이도 측면에서 우수한 것으로 판단되었다.As can be seen in Table 2, 21.4% of stage I patients with R-ISS were in the medium-risk group in the new model, and 32.9% and 32.9% of stage II patients in R-ISS were in the low-risk and high-risk groups, respectively, in the new model. - 48.5% of ISS stage III patients were reclassified to the intermediate risk group of the new model, and this difference was statistically significant. The NRI (Net Reclassification Improvement) calculated by reclassification was 0.337, and it was judged that the new model was superior to the R-ISS in terms of sensitivity and specificity.

TotalTotal R-ISSR-ISS -- -- 본 발명의 예후 예측모델Prognosis prediction model of the present invention Stage IStage I Stage IIStage II Stage IIIStage III totaltotal PP 저위험군low risk group 11 (78.6%)11 (78.6%) 26 (32.9%)26 (32.9%) 0 (0.0%)0 (0.0%) 37 (29.4%)37 (29.4%) <0.001<0.001 중위험군medium risk 3 (21.4%)3 (21.4%) 27 (34.2%)27 (34.2%) 16 (48.5%)16 (48.5%) 46 (36.5%)46 (36.5%) -- 고위험군high risk group 00 26 (32.9%)26 (32.9%) 17 (51.5%)17 (51.5%) 43 (34.1%)43 (34.1%) --

실시예Example 6: 새로운 예후 예측 모델의 6: The new prognostic prediction model 각 군별by each group 치료법에 대한 효과 분석 Efficacy analysis for treatment

새로운 예후 예측 모델의 저위험군, 중간위험군, 고위험군에서 카플란-메이어(Kaplan-Meier) 생존곡선과 다변수 콕스 회귀분석을 시행하여 전체생존율에 대한 자가조혈모세포이식과 frontline 면역조절제 투약 여부가 생존에 미치는 영향을 분석하였다. 도 8과 표 3에서 확인할 수 있듯이, 저위험군과 중간위험군에서는 frontline 면역조절제 투약과 자가조혈모세포 이식이 생존에 영향을 주지 않았으나 고위험군에서는 자가조혈모세포 이식이 높은 전체생존 및 무진행생존과 관련이 있는 것으로 확인되었다.The effect of autologous hematopoietic stem cell transplantation and frontline immunomodulator administration on overall survival on survival by performing Kaplan-Meier survival curves and multivariate Cox regression analysis in the low-risk, intermediate-risk, and high-risk groups of the new prognostic prediction model. The impact was analyzed. As can be seen in Figure 8 and Table 3, administration of frontline immunomodulatory agents and autologous hematopoietic stem cell transplantation did not affect survival in the low-risk and intermediate-risk groups, but in the high-risk group, autologous hematopoietic stem cell transplantation was associated with high overall survival and progression-free survival. confirmed to be

따라서 본 발명의 새로운 예후 예측 모델을 이용하면 다발골수종의 예후 예측뿐만 아니라 고위험군 환자의 자가 조혈모세포이식 여부를 결정하는데 도움을 줄 수 있을 것으로 생각된다.Therefore, it is thought that the use of the new prognostic prediction model of the present invention can help in predicting the prognosis of multiple myeloma as well as determining whether autologous hematopoietic stem cell transplantation in high-risk patients.

Figure 112020053853123-pat00003
Figure 112020053853123-pat00003

Claims (3)

다음의 단계를 포함하는 다발성 골수종(multiple myeloma)의 예후 예측을 위한 정보 제공 방법:
(a) 다발골수종 환자의 골수로부터 수득한 형질세포를 포함하는 조직 시료에서 면역조직염색을 통해 CD138 및 PD-L1의 발현 정도를 정량적으로 측정하여 하기 식1에 따른 PD-L1의 정규화 MFI(normalized mean fluorescence intensity)를 도출하는 단계:
식1
Normalized MFI = (MFI in plasma cell compartments - Background intensity) / MFI from isotype-matched control;
(b) 하기 식2에 따른 PD-L1 발현 MFI 측정값을 도출하는 단계:
식2
PD-L1 expression MFI = Sum of normalized MFIs in all plasma cell compartments / Total number of plasma cells; 및
(c) 상기 식2에 따른 PD-L1 발현 MFI 측정값이 7.65 미만이면 저위험군, 7.65 이상이면 고위험군으로 분류하는 단계.
A method of providing information for predicting the prognosis of multiple myeloma comprising the steps of:
(a) By quantitatively measuring the expression levels of CD138 and PD-L1 through immunohistochemistry in a tissue sample containing plasma cells obtained from the bone marrow of a patient with multiple myeloma, normalized MFI of PD-L1 according to Equation 1 below (normalized MFI) deriving the mean fluorescence intensity):
Expression 1
Normalized MFI = (MFI in plasma cell compartments - Background intensity) / MFI from isotype-matched control;
(b) deriving a PD-L1 expression MFI measurement value according to Equation 2 below:
Expression 2
PD-L1 expression MFI = Sum of normalized MFIs in all plasma cell compartments / Total number of plasma cells; and
(c) classifying the PD-L1 expression MFI according to Equation 2 into a low-risk group if the measured value is less than 7.65, and a high-risk group if it is 7.65 or more.
다음의 단계를 포함하는 다발성 골수종의 예후 예측을 위한 정보 제공 방법:
Figure 112021107909597-pat00004

(a) i) 제1항의 식1 및 식2에 따른 PD-L1 발현 MFI 측정값(0~35), ii) 정상치 상한 초과의 혈중 LDH(serum lactate dehydrogenase) 여부, iii) high risk cytogenetics (karyotype에서 t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, FISH에서 IGH/FGFR3 rearrangement, IGH/MAF rearrangement, 및 TP53 mutation으로 이루어진 군으로부터 1종 이상 양성인 경우) 여부, iv) 연령 70세 이상 여부를 기준으로 상기 nomogram에서 상응하는 점수(points, 0~100)를 가산하는 점수 배점 단계; 및
(b) 상기 각 가산된 점수를 합산하여 총점(total points)을 구한 후 상기 노모그램에서 상응하는 1년, 2년, 또는 4년 전체 생존기간에 대한 확률을 도출하는 단계.
A method of providing information for predicting the prognosis of multiple myeloma comprising the steps of:
Figure 112021107909597-pat00004

(a) i) PD-L1 expression MFI measurement value (0 to 35) according to Equation 1 and Equation 2 in Paragraph 1, ii) Serum lactate dehydrogenase (LDH) in blood above the upper limit of normal, iii) high risk cytogenetics (karyotype) In the case of t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, and at least one positive from the group consisting of IGH/FGFR3 rearrangement, IGH/MAF rearrangement, and TP53 mutation in FISH), iv ) a score scoring step of adding corresponding scores (points, 0-100) in the nomogram based on whether or not the age is 70 years or older; and
(b) deriving probabilities for the corresponding 1-year, 2-year, or 4-year overall survival period from the nomogram after summing the added scores to obtain total points.
다음의 단계를 포함하는 다발성 골수종 환자의 자가 조혈모세포이식 여부를 결정하는데 필요한 정보를 제공하는 방법:
Figure 112021107909597-pat00005

(a) i) 제1항의 식1 및 식2에 따른 PD-L1 발현 MFI 측정값(0~35), ii) 정상치 상한 초과의 혈중 LDH(serum lactate dehydrogenase) 여부, iii) high risk cytogenetics (karyotype에서 t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, FISH에서 IGH/FGFR3 rearrangement, IGH/MAF rearrangement, 및 TP53 mutation으로 이루어진 군으로부터 1종 이상 양성인 경우) 여부, iv) 연령 70세 이상 여부를 기준으로 상기 nomogram에서 상응하는 점수(points, 0~100)를 가산하는 점수 배점 단계;
(b) 상기 각 가산된 점수를 합산하여 총점(total points)을 구한 후 50점 미만이면 저위험군, 50 내지 100점이면 중간위험군, 100점을 초과하면 고위험군으로 분류하는 분류 단계; 및
(c) 총점이 100점을 초과하는 고위험군으로 분류된 환자에게 자가조혈모세포 이식시 생존율이 향상된다는 정보를 제공하는 단계.
A method of providing information necessary to determine whether a patient with multiple myeloma should receive an autologous hematopoietic stem cell transplant comprising the steps of:
Figure 112021107909597-pat00005

(a) i) PD-L1 expression MFI measurement value (0 to 35) according to Equation 1 and Equation 2 in Paragraph 1, ii) Serum lactate dehydrogenase (LDH) in blood above the upper limit of normal, iii) high risk cytogenetics (karyotype) In the case of t(4;14), t(14;16), 17p deletion, chromosome 1 abnormality, and at least one positive from the group consisting of IGH/FGFR3 rearrangement, IGH/MAF rearrangement, and TP53 mutation in FISH), iv ) a score scoring step of adding corresponding scores (points, 0-100) in the nomogram based on whether or not the age is 70 years or older;
(b) a classification step of adding up each added score to obtain a total point and classifying it as a low-risk group if it is less than 50 points, a medium-risk group if it is 50 to 100 points, and a high-risk group if it exceeds 100 points; and
(c) providing information that survival rate is improved upon autologous hematopoietic stem cell transplantation to patients classified as high-risk group with a total score of more than 100 points.
KR1020200063916A 2019-07-19 2020-05-27 Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression KR102360297B1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
KR1020190087805 2019-07-19
KR20190087805 2019-07-19

Publications (2)

Publication Number Publication Date
KR20210010311A KR20210010311A (en) 2021-01-27
KR102360297B1 true KR102360297B1 (en) 2022-02-09

Family

ID=74238559

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020200063916A KR102360297B1 (en) 2019-07-19 2020-05-27 Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression

Country Status (1)

Country Link
KR (1) KR102360297B1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2022158950A2 (en) 2021-01-25 2022-07-28 주식회사 엘지화학 Separator for lithium secondary battery and lithium secondary battery having same

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019097082A1 (en) 2017-11-20 2019-05-23 Julius-Maximilians-Universität Würzburg Cd19cart cells eliminate myeloma cells that express very low levels of cd19

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20160086496A (en) * 2015-01-09 2016-07-20 순천대학교 산학협력단 Selection method of predicting genes for ovarian cancer prognosis

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019097082A1 (en) 2017-11-20 2019-05-23 Julius-Maximilians-Universität Würzburg Cd19cart cells eliminate myeloma cells that express very low levels of cd19

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Hidetoshi Sumimoto et al, PLoS ONE (2016), vol 11(11): e0166626. doi:10.1371/journal, pp 1-18.
S Yousef et al, Blood Cancer Journal (2015), vol 5, e285, pp 1-3.

Also Published As

Publication number Publication date
KR20210010311A (en) 2021-01-27

Similar Documents

Publication Publication Date Title
Lænkholm et al. PAM50 risk of recurrence score predicts 10-year distant recurrence in a comprehensive Danish cohort of postmenopausal women allocated to 5 years of endocrine therapy for hormone receptor–positive early breast cancer
Brims et al. A novel clinical prediction model for prognosis in malignant pleural mesothelioma using decision tree analysis
Cuzick et al. Prognostic value of a combined estrogen receptor, progesterone receptor, Ki-67, and human epidermal growth factor receptor 2 immunohistochemical score and comparison with the Genomic Health recurrence score in early breast cancer
Balic et al. Circulating tumor cells: from bench to bedside
Dawood et al. Defining breast cancer prognosis based on molecular phenotypes: results from a large cohort study
Du et al. Refining the utility and role of frozen section in head and neck squamous cell carcinoma resection
Boldrini et al. Prognostic significance of osteopontin expression in early-stage non-small-cell lung cancer
US20100099093A1 (en) Biomarkers for the Identification Monitoring and Treatment of Head and Neck Cancer
US10180430B2 (en) Methods for detecting human papillomavirus and providing prognosis for head and neck squamous cell carcinoma
Bessonova et al. Risk of a second breast cancer associated with hormone-receptor and HER2/neu status of the first breast cancer
Lin et al. Clinicopathological features of salivary and non-salivary adenoid cystic carcinomas
Zou et al. A four‐factor immune risk score signature predicts the clinical outcome of patients with spinal chordoma
Brcic et al. Prognostic impact of PD-1 and PD-L1 expression in malignant pleural mesothelioma: an international multicenter study
Zhao et al. Association of preoperative serum levels of CEA and CA15-3 with molecular subtypes of breast cancer
EP2473854B1 (en) Systems and methods for treating, diagnosing and predicting the response to therapy of breast cancer
JP2023156317A (en) cancer prognosis
Galindo et al. The value of desmosomal plaque-related markers to distinguish squamous cell carcinoma and adenocarcinoma of the lung
KR102360297B1 (en) Method for providing the information for prediction of multiple myeloma prognosis using PD-L1 expression
Wu et al. Combining immunoscore with clinicopathologic features in cholangiocarcinoma: an influential prognostic nomogram
JP6417602B2 (en) How to determine the prognosis of colorectal cancer
Pezzuto et al. Evaluation of prognostic histological parameters proposed for pleural mesothelioma in diffuse malignant peritoneal mesothelioma. A short report
Zhang et al. Combined expression of the BRAFV600E mutation and PD-L1 in early papillary thyroid carcinoma and its relationship with clinicopathological features and recurrence—a retrospective cohort study
Yaghjyan et al. Reliability of CD44, CD24, and ALDH1A1 immunohistochemical staining: Pathologist assessment compared to quantitative image analysis
Isaksen et al. The Geriatric Prognostic Index: a clinical prediction model for survival of older diffuse large B-cell lymphoma patients treated with standard immunochemotherapy
Cao et al. A nomogram for predicting hyperprogressive disease after immune checkpoint inhibitor treatment in lung cancer

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right