This application claims benefit of priority from Taiwan Application No. 111102676, filed on Jan. 21, 2022. The entire disclosures of all the above applications are hereby incorporated by reference herein.
The present invention relates to a method and application for in vitro predicting survival time of individual with bladder cancer after surgery. More specifically, the present invention relates to a method, kit, and a computer-implemented method for predicting survival time of individual with bladder cancer after surgery from an individual's biological sample according to expression levels of specific target gene combination as criteria to molecular typing.
Bladder cancer is a common aggressive disease of the urinary system. Though bladder cancer can occur at any age, its risk rises with age, especially in the age group of 50 to 70 years old. Bladder cancer is three to four times more frequent in men than women. In Taiwan, bladder cancer is not among the top ten cancer deaths (men and women combined), but bladder cancer ranked 9th in male cancers in Taiwan in 2017. Thus, bladder cancer is one of the common urological malignancies.
The known risk factors of bladder cancer include smoking, race, genetics, chronic bladder inflammation in a long term, receiving anti-cancer medication, and the environment (such as those working in rubber, chemical dyes, printing, leather shoes, hair dye, paint, etc.).
According to the depth of tumor invasion, bladder cancer can be classified into non-muscle invasive bladder cancer (NMIBC) and muscle invasive bladder cancer (MIBC). Approximately 75% of newly diagnosed bladder cancers are NMIBC, and MIBC accounts for approximately 25% of bladder cancers. The survival rate of NMIBC patients is relatively high, but approximately 60% of the NMIBC patients develop recurrence even after surgery, and approximately 80% of the NMIBC patients develop recurrence within one year after surgery, while approximately 15% progress to MIBC. Once diagnosed with metastatic or advanced unresectable bladder cancer, approximately 70% of patients will die within 2 years. The 5-year survival rate of those patients is about 60% with radical surgery, but the 5-year survival rate of the patients with lymph node metastasis is only approximately 15%, and the 5-year survival rate of patients with other organ metastasis is approximately 5%.
Arriving at a diagnosis and a treatment of bladder cancer may involve intravenous pyelography (IVP) and cystoscope, both of which are often used to detect bladder cancer and to assist the surgery (e.g., transurethral resection of bladder tumor, TUR-BT). Chemotherapy and radiation therapy can be optionally combined after surgery. In general, bladder cancer can be treated successfully when detected at an early stage, and the average 5-year survival rate after surgery can reach approximately 60%. However, when the bladder is assessed by IVP using a cystoscope, there are problems such as dead zones of the lens and blurry images sometimes difficult to ascertain, leading in a high recurrence rate after surgery. Periodic follow-up observations have to be organized, and TUR-BT can be repeated if necessary during the follow-up period. However, before, during and after the TUR-BT operation, various pains, inconvenience and discomfort make the patients feel stressed, be less willing to undergo surgery, even ignore the importance of regular check-ups, resulting in possible deterioration and metastasis of cancers.
The bladder cancer is a highly heterogeneous disease in the molecular genetics and histopathology. Most NMIBC carcinomas are grown as papillary shapes, while MIBC carcinomas are associated with flat atypical hyperplasia and carcinoma in situ. Fibroblast growth factor receptor 3 (FGFR3) promoter mutations are frequent in NMIBC. Inactivation or mutation of p53 gene may confer MIBC with highly genetic instability.
Since the bladder cancer is a highly heterogeneous disease, the conventional histopathological classification cannot meet the clinical requirement. There is an urgent need to discover molecular typing markers of the bladder cancer for predicting the survival time of an individual with the bladder cancer after surgery, thereby assisting to choose a more suitable treatment strategy.
In an aspect, the invention provides a method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. In the method, expression levels of target gene combination of in vitro aggressive bladder cancer specimen of an individual are detected, in which the target gene combination can include at least two of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, a fragment, a homologue, a variant and a derivative thereof. The expression levels of the target gene combination are respectively compared with the reference expression levels of the target gene combination of a reference database, and converted to a risk score sum. When the risk score sum of an individual is equal to or greater than a threshold, the individual is classified into a high risk group, thereby increasing the predicting accuracy of an averaged survival time of an individual having aggressive bladder cancer after surgery.
In yet a further aspect, the invention also provides a computer-implemented method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. The computer-implemented method includes the steps of providing a biological specimen, detecting expression levels of the target gene combination of the biological specimen, comparing one of the expression levels with a corresponding reference expression level of the target gene combination of a reference database respectively, obtaining a difference and a risk score, and calculating a risk score sum of at least two of the target gene combination of the biological specimen, in which the step of comparing one of the expression levels with a corresponding reference expression level of the target gene combination of a reference database respectively, the step of obtaining a difference and a risk score, and the step of calculating a risk score sum of at least two of the target gene combination of the biological specimen can be executed within a computer system for implementing the aforementioned method.
In yet a further aspect, the invention also provides a system for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. The system includes a detecting module, a comparing module, a determining module and a controlling system, thereby increasing the predicting accuracy of an averaged survival time of an individual having aggressive bladder cancer after surgery.
According to the aforementioned aspect, the present invention provides a method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. In an embodiment, the method can include the step of establishing a reference database. Next, a biological sample is provided. And then, expression levels of the target gene combination of the biological specimen are detected, in which the target gene combination can include at least two of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, a fragment, a homologue, a variant and a derivative thereof. Later, one of the expression levels with a corresponding reference expression level of the target gene combination of a reference database respectively, for obtaining a difference and a risk score. A risk score of the one is given as 1 when an absolute value of the difference is equal to or more than a first threshold that is more than at least 5% of the reference expression level. Consequently, a risk score sum of at least two of the target gene combination of the biological specimen is calculated. When the risk score sum is equal to or greater than a second threshold that is 1 or 2, the individual having the biological specimen is classified to a high risk group, wherein the high risk group is defined by less 25 months of an averaged survival time after surgery of the individual with 1 or 2 of the second threshold, or by less 10 months of an averaged survival time after surgery of the individual with 3 of the second threshold.
According to another aspect, the invention provides a kit for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. The kit includes a reaction solution, nucleic acid probes and/or antibodies, in which the nucleic acid probes and/or antibodies react with a target gene combination of a biological specimen and generate expression levels, the biological specimen is originated from in vitro aggressive bladder cancer of an individual, and the target gene combination is selected from at least two of the group consisting of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, a fragment, a homologue, a variant and a derivative thereof.
According to yet another aspect, the invention provides a computer-implemented method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. In an embodiment, the computer-implemented method includes the steps of storing a reference data bank in a system, detecting expression levels of the target gene combination of a biological specimen, comparing one of the expression levels with a corresponding reference expression level of the target gene combination, calculating a risk score sum of at least two of the target gene combination, and identifying the individual from the risk score sum. A program executes the aforementioned steps of storing, comparing, calculating and identifying, for implementing the method on a computer system.
According to still another aspect, the invention provides a system for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. In an embodiment, the system includes a detecting module, a processing module, an identifying module and a controlling module.
In the aforementioned embodiment, the detecting module can include a detecting device, a reaction solution, nucleic acid probes and/or antibodies. The nucleic acid probes and/or antibodies react with a target gene combination of a biological specimen and generate expression levels, and the detecting device can detect the expression levels. The biological specimen can be originated from in vitro aggressive bladder cancer of an individual. The target gene combination can include but be not limited to at least two of the group consisting of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, a fragment, a homologue, a variant and a derivative thereof.
The processing module can be coupled to the detecting module, for receiving the expression levels from the detecting module and comparing one of the expression levels with a corresponding reference expression level of the target gene combination of a reference database respectively, thereby obtaining a difference and a risk score; and for calculating a risk score sum of at least two of the target gene combination of the biological specimen. The risk score of the one is given as 1 when an absolute value of the difference is equal to or more than a first threshold that is more than at least 5% of the reference expression level.
The identifying module can be coupled to the processing module, for identifying the individual from the risk score sum. When the risk score sum is equal to or greater than a second threshold that is 1 or 2, the individual having the biological specimen is classified to a high risk group, in which the high risk group is defined by less 25 months of an averaged survival time after surgery of the individual with 1 or 2 of the second threshold, or by less 10 months of an averaged survival time after surgery of the individual with 3 of the second threshold.
The controlling module can be coupled to the detecting module, the processing module and the identifying module. The system executes a program that includes instructions for activating the detecting module, the processing module and the identifying module.
In the aforementioned embodiments, the computer system further optionally includes a pretreating module coupling to the detecting module, for providing a nucleic acid sample and/or a protein sample of the biological specimen.
With application to the method and the kit for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample, which detect expression levels of target gene combination of in vitro aggressive bladder cancer specimen of an individual, respectively compare those expression levels with the reference expression levels of the target gene combination of a reference database, the results can be converted to a risk score sum, thereby increasing the predicting accuracy of an averaged survival time of an individual having aggressive bladder cancer after surgery, leading in applications on the computer-implemented method and the system for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample.
It is to be understood that both the foregoing general description and the following detailed description are by examples, and are intended to provide further explanation of the invention as claimed.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawings will be provided by Office upon request and payment of the necessary fee. The disclosure can be more fully understood by reading the following detailed description of the embodiment, with reference made to the accompanying drawings as follows.
Where a definition or use of a term in an incorporated reference is inconsistent or contrary to the definition of that term provided herein, the definition of that term provided herein applies and the definition of that term in the reference does not apply. Moreover, unless the context otherwise requires, singular terms can include the plural and vice versa.
As aforementioned, the present invention provides a method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample, in which expression levels of target gene combination of in vitro aggressive bladder cancer specimen of a patient are detected and compared with the reference expression levels of the target gene combination of a reference database, thereby predicting an averaged survival time of a patient having aggressive bladder cancer after surgery.
The terms “subject”, “individual” and “patient” are used interchangeably herein and refer generally to a subject with bladder cancer, and the subject can be a human or a non-human mammal. In an embodiment, the subjects with bladder cancer include NMIBC and MIBC individuals.
The term “survival time after surgery of the bladder cancer” generally refers to an averaged survival time after surgery of the individual with the bladder cancer, and it is also called as an averaged survival time after surgery of the bladder cancer. The “after surgery” typically refers to an individual with the bladder cancer after a radical surgery (e.g., radical cystectomy or partial cystectomy) or transurethral resection of bladder tumor (TUR-BT).
The term “biological sample” described herein can be originated from an in vitro aggressive bladder cancer from an individual. The reference expression level can be originated from an in vitro healthy bladder. The biological sample and the in vitro healthy bladder sample can include but be not limited to ex vivo organs, tissues, cells, body fluid, lymphatic liquid, urine, whole bloods, serum and/or cell culture supernatant, as well as nucleic acid extract (such as genomic DNA extract, mRNA extract, cDNA or cRNA obtained from the mRNA extract, etc.) or protein extract obtained from the aforementioned biological sample.
The term “target gene” or “target gene combination” described herein can include but be not limited to at least one of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, and the fragment, the homologue gene, the variant gene or the derivative gene of the aforementioned genes. In other examples, the target gene can be a biomarker.
The term “expression levels of target genes” can include expression levels of nucleic acids (DNA or RNA) and/or proteins. In some examples, one of the expression levels of target genes can be obtained as follows. A biological sample can be optionally pretreated to obtain nucleic acid extract and/or protein extract. Next, an expression level of one gene of the target gene combination (e.g., nucleic acid expression level and/or protein expression level) can be detected. In an example, the relative expression level of the gene can be normalized by a reference expression level (as a normalization data) of the target gene combination in a reference database corresponding to the gene from the nucleic acid extract and/or protein extract. In the aforementioned example, the reference database can include a plurality of expression levels of the target gene combination originated from at least one in vitro normal bladder sample.
The term “standardization” or “normalization” herein refers to standardizing or normalizing the data detected from the biological sample according to a standard the data from a normal (or healthy) sample, for further comparing and analyzing the data.
The term “risk score” herein refers to a difference obtained by comparing one of the aforementioned expression levels with the reference expression levels of the target gene combination in the reference database respectively, in which the absolute value of the difference that is equal to or greater than the reference expression level reaches a first threshold, a risk score of the one is given as 1. In this embodiment, there is no specific limitation to the first threshold, for example, at least 5%. In some embodiments, the expression levels of P4HB, SLC1A6 and ARID3A in the cancer lesion can be significantly greater than the ones in the normal tissue, so that the three biomarkers can be potentially applied to predict the averaged survival time after surgery of an individual with the bladder cancer in future. In other embodiments, the expression levels of PPT2 and ARMH4 in the surrounding tissue of the cancer lesion can be significantly higher than the ones in the normal tissue, so that the two biomarkers can be potentially as the biomarkers for predicting the precancerous lesion in future.
It should be clarified that, the difference, the first threshold or the risk score as aforementioned is merely illustrated but is not limited within the specific ranges or specific values. In other examples, the first threshold can be given to at least 6% to at least 10% or other values depending upon the actual requirements.
In an embodiment, a risk score sum of at least two of the target gene combination of the biological specimen can be calculated. When the risk score sum is equal to or more than a second threshold that is 1 or 2, the individual having the biological specimen is classified to a high risk group. In this embodiment, there is no specific requirement to the second threshold depending upon the definition of the risk score. For example, when the risk score is 1, the second threshold can be exemplified as 2; in some example, when the risk score is 0.5, the second threshold can be exemplified as 1; in certain examples, when the risk score is 2, the second threshold can be exemplified as 4. In the aforementioned embodiment, the “high risk group” is defined herein by less 25 months of an averaged survival time after surgery of the individual with 1 or 2 of the second threshold. In other embodiments, the “high risk group” is defined by less 10 months of an averaged survival time after surgery of the individual with 3 of the second threshold.
In an embodiment, a combination of other treatment can be optionally subjected to an individual after surgery. The term “treatment” can include but be not limited to chemotherapy and/or radiotherapy and/or surgery, for healing or improving the cancer or extending a lifespan of an individual. It should be supplemented that, the five target genes (or biomarkers) of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, as well as their upstream and downstream gene modulating pathways, can be as molecular typing markers, for designing personalized precision medicine, thereby elevating the diagnostic accuracy and treating effects.
It should be mentioned that, in the method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample, the specific target gene combinations are used in the aforementioned method according to specific criteria (i.e., the difference, the risk scores, the first threshold, the second threshold), so as to accurately predict the survival time of an individual with bladder cancer after surgery. If other genes excepted from the specific target gene combinations were used in the aforementioned method or the criteria changed, such result according to the modified method could not accurately predict the survival time of an individual with bladder cancer after surgery.
In some embodiments, a computer-implemented method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample can be applied to a computer program within a computer system for implementing the method, in which the computer-implemented method can include steps as follows. Firstly, a reference data bank that comprises reference expression levels of a target gene combination originated from at least one in vitro normal bladder specimen is stored in a system. Next, expression levels of the target gene combination of a biological specimen originated from in vitro aggressive bladder cancer of the individual are detected. And then, one of the expression levels of the biological specimen is compared with a corresponding reference expression level of the target gene combination of a reference database respectively, for obtaining a difference and a risk score. Subsequently, a risk score sum of at least two of the target gene combination of the biological specimen is calculated. A program executes the aforementioned steps of comparing one of the expression levels with a corresponding reference expression level of the target gene combination of a reference database respectively, obtaining a difference and a risk score, and calculating a risk score sum of at least two of the target gene combination of the biological specimen, for implementing the method on a computer system. In some examples, the program used for executing the computer-implemented method for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample can be also called as a biomarker screening software.
In other embodiments, the computer program can be applied to a system for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. The system includes a detecting module, a processing module, an identifying module and a controlling module. In this embodiment, the detecting module can include a detecting device, a reaction solution, nucleic acid probes and/or antibodies. The detecting module can be performed in combination with conventional detection reagents, a detecting equipment, a microarray chip and so on, which are common knowledge in the technical field of the present invention instead of reciting in detail. The nucleic acid probes and/or the antibodies react with the target gene combination of the biological sample for generating expression levels that are detected by the detecting device. The aforementioned biological specimen can be originated from in vitro aggressive bladder cancer of an individual, and the target gene combination is selected from at least two of the group consisting of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A, a fragment, a homologue, a variant and a derivative thereof. The processing module can be coupled to the detecting module, for receiving the expression levels from the detecting module and comparing one of the expression levels with a corresponding reference expression level of the target gene combination of a reference database respectively, thereby obtaining a difference and a risk score. The processing module also can calculate a risk score sum of at least two of the target gene combination of the biological specimen, and the risk score of the one is given as 1 when an absolute value of the difference is equal to or more than a first threshold that is more than at least 5% of the reference expression level. The identifying module can be coupled to the processing module, for identifying the individual from the risk score sum. When the risk score sum is equal to or greater than a second threshold that is 1 or 2, the individual having the biological specimen is classified to a high risk group, the high risk group is defined by less 25 months of an averaged survival time after surgery of the individual with 1 or 2 of the second threshold, or by less 10 months of an averaged survival time after surgery of the individual with 3 of the second threshold. The controlling module can be coupled to the detecting module, the processing module and the identifying module. The system executes a program that includes instructions for activating the detecting module, the processing module and the identifying module. In the above embodiment, the system optionally includes a pretreating module that is coupled to the detecting module, for providing a nucleic acid sample and/or a protein sample of the biological specimen.
Thereinafter, it will be understood that particular cancer subjects, biological sample, target gene combinations, detection methods, determination criteria, aspects, examples and embodiments described hereinafter are shown by way of illustration and not as limitations of the invention. The principal features of this invention can be employed in various embodiments without departing from the scope of the invention. Thus, one skilled in the art can easily ascertain the essential characteristics of the present invention and, without departing from the spirit and scope thereof, can make various changes and modifications of the invention to adapt it to various usages and conditions.
1.1 Use and Statistical Analysis of Data from The Cancer Genome Atlas (TCGA)
Reference was made to
Firstly, a bladder cancer dataset (TCGA, Cell 2017) having the most samples (n=413) in the bladder/urinary tract in the cBioPortal online database was taken in this Example, the patients with bladder cancer that had complete RNA expression and clinical data were adopted (n=408) among these, and the patients with bladder cancer that had events of overall survivals (n=405) were further chosen (n=367). And then, the gene lists (n=18883) excluded NA (not applicable/not available) value from 20,435 gene lists of RNA expressions of patients (n=18883). Next, the survival-related genes (n=1279) were calculated by univariate Cox regression analysis (p-value<0.05) from the first screening result (n=7204, the screening criteria: the expression rate>5%, Z-score>2), and then corrected by False Discovery Rate (FDR) (adjusted p value<0.05), so as to remove erroneous genes as much as possible (as shown in
The 145 genes in
For the purpose of realizing the predictability of the model, this Example proposed to find genes that were related to bladder cancer, overexpressed and reported in the journals. According to the first screening criteria (i.e., expression rate>5%, Z-score>2), two gene lists (n=5, also called as 5g, known oncogenes with abnormal overexpression, as shown in the vertical axis of
Reference was made to
In addition, 26 genes screened from by the computer program for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample of
The genes (n=26, as shown in
In this Example, a survival rate of a MIBC patient dataset (GSE13507) from GEO database and the genes (n=26) screened from
The data of the MIBC patient (GSE13507) of GEO database and the patient (cell 2017) with bladder cancer of TCGA database were analyzed, for assessing the overall risk trend and producing the heat map (as shown in
Reference was made to
Reference was made to
As shown in the result of
In addition, PPT2, ARMH4, P4HB, SLC1A6 and ARID3A were subjected to Kaplan-Meier survival analysis and the result was shown in
Reference was made to
As shown in previous research, SLC1A5 in the SLC family was associated with glutamate and cancer, but there was no correlation between them found in the analysis of SLC1A5 and the survival rate of patients with bladder cancer. However, this Example proved that genes of ARID3A, ARMH4, P4HB, PPT2, SLC1A6 and so on were correlated with the survival rate of patients with bladder cancer. Moreover, the aforementioned results confirmed that when at least two of the target gene combination of PPT2, ARMH4, P4HB, SLC1A6 and ARID3A were detected in the biological sample, the survival rate of the clinical prognosis of the individual corresponding to the biological sample was actually poor, and an averaged survival time after surgery of the individual would be less than 25 months.
Reference was made to
Reference was made to
Reference was made to
Reference was made to
Reference was made to
Reference was made to
In
In addition, reference was made to
Reference was made to
According to the immunohistochemical staining results of
In summary, specific cancers, specific biological samples, specific target gene combinations, specific analysis models or specific evaluating methods are exemplified for clarifying the method and the kit, the computer-implemented method and the system for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample. However, as is understood by a person skilled in the art, other cancers, other biological samples, other target gene combinations, other analysis models or other evaluating methods can be also adopted in the method and the kit, the computer-implemented method and the system for in vitro predicting survival time of an individual with bladder cancer after surgery from the individual's biological sample without departing the spirit and scope of the present invention rather than being limited as aforementioned. For example, the aforementioned target gene combinations can be added with other genes as molecular typing markers, for optimizing the method and the kit, the computer-implemented method and the system, thereby elevating the accuracy for predicting an averaged survival time of a patient having aggressive bladder cancer after surgery, and beneficially choosing the more appropriate strategies of treatments.
According to the embodiments of the present invention, the method, the kit, the computer-implemented method and the system for in vitro predicting survival time of an individual with bladder cancer after surgery from an individual's biological sample, which beneficially detect expression levels of target gene combination of in vitro aggressive bladder cancer specimen of an individual, respectively compare those expression levels with the reference expression levels of the target gene combination of a reference database, the results can be converted to a risk score sum, thereby increasing the predicting accuracy of an averaged survival time of an individual having aggressive bladder cancer after surgery.
Although the present invention has been described in considerable detail with reference to certain embodiments thereof, other embodiments are possible. Therefore, the spirit and scope of the appended claims should not be limited to the description of the embodiments contained herein.
Number | Date | Country | Kind |
---|---|---|---|
111102676 | Jan 2022 | TW | national |