- 1Department of Ophthalmology, Eye & Ear, Nose, and Throat (ENT) Hospital of Fudan University, Shanghai, China
- 2Shanghai Key Laboratory of Visual Impairment and Restoration, Fudan University, Shanghai, China
- 3Center of Clinical Epidemiology and Evidence-based Medicine, Fudan University, Shanghai, China
- 4Department of Pathology, Eye & Ear, Nose, and Throat (ENT) Hospital of Fudan University, Shanghai, China
Purpose: To establish an easy and widely applicable prognostic prediction model for uveal melanoma (UM) based on a Chinese population.
Patients and Methods: A total of 295 consecutive cases treated at the Eye & ENT Hospital of Fudan University were included as the primary cohort, and 256 cases were included in the validation cohorts from two external Caucasian databases. Clinicopathological data were collected retrospectively, and nomogram models were formulated based on multivariable analysis. The concordance index (C-index), AUC (area under the Receiver Operating Characteristic, ROC curve), and Brier score were calculated and compared.
Results: Based on the training cohort, a nomogram model was established with five relevant variables: age, tumor size, ciliary body involvement, non-spindle cell type and extra-scleral extension. The C-index was 0.737, the 3- and 5-year AUCs were 0.767 and 0.742, and the Brier scores for 3- and 5-year survival were 0.082 and 0.129, respectively, which showed superior prediction compared to that of the Tumor, Node and Metastasis staging system. The model also displayed good discrimination and calibration in the external validation cohorts. By risk stratification, patients could be divided into low- and high-risk groups, and the overall survival curves displayed significant differences in the training and validation cohorts.
Conclusion: Our nomogram model was simple and accurate at predicting the overall survival of patients with UM. It was established based on Asian patients and proved suitable for Caucasian patients; thus, it has a wide range of potential applications, especially for patients living in less medically developed countries and regions.
Introduction
Uveal melanoma (UM) is the most common intraocular primary malignancy in adults involving the choroid, iris and ciliary body, and has a high tendency to metastasize, resulting in high mortality (1–3). Approximately half of the patients develop metastatic disease during a long follow-up (1, 2). Currently, there is no proven effective treatment for patients with metastasis. Once metastasis is detected, the patient survival rate rapidly decreases to approximately 15% at one year, and the median survival time is estimated to be 4 to 15 months (4). Some clinical trials are being conducted with patients with metastasis, which may provide some benefits for patients (5, 6). Therefore, early screening of high-risk patients and early detection of metastasis are very important. Moreover, establishing risk stratification and making different follow-up plans accordingly can also save many medical resources, reduce unnecessary invasive examinations and testing expenses, and thus result in rational allocation of medical resources.
Worse prognostic indicators for UM include 1) clinical factors, such as older age at presentation, male sex, larger basal tumor diameter, high tumor thickness, ciliary body involvement, and extraocular spread (7, 8); 2) histopathological factors, such as epithelioid melanoma cytomorphology, certain extravascular matrix patterns such as closed loops (9), a high mitotic count, and high tumor microvascular density (8); and 3) genetic abnormalities, including chromosomal changes (chromosome 3 loss, 8q gain, 1p loss and 6p gain) (8), class 2 gene expression profiles (10), and other biomarkers (BAP1, SF3B1, EIF1AX, and PRAME) (11, 12). Most of the prognostic prediction models have been established based on these risk indicators, such as the TNM staging system, the gene expression profile (10, 13), the Liverpool Uveal Melanoma Prognosticator Online (14), the Cancer Genome Atlas Classification (15) and other models (16–19). Most of the present prognostic models are relatively complicated and less feasible. The high expenses of these tests may also be limitations for some patients, especially those living in developing regions and countries.
To our knowledge, no prognostic models have been analyzed and established based on Asian patients. Asian patients, compared to the Caucasian population, have some unique features, such as a lower incidence probability, earlier onset age, and larger tumor diameter (20, 21). The two largest clinical prognosis studies in China showed that Chinese patients have a younger onset age and no sex difference and tend to have a better prognosis (22, 23). Studies from other Asian countries (Japan (24), India (25), Korea (26) and Singapore (27) also showed similar results.
Based on the long-term follow-up results of the patients in our hospital, we established a novel and simple nomogram prediction model to predict the survival of UM patients. Two public databases from Western countries were analyzed for external verification.
Methods
Patients and Clinicopathological Data of the Training Cohort
We retrospectively and consecutively reviewed all patients with UM who were treated between January 2002 and October 2017 at the Eye & ENT Hospital of Fudan University in Shanghai, one of the three largest tertiary eye centers in China. All patients were diagnosed by experienced ocular surgeons. All patients received enucleation or proton radiotherapy during this time period. Informed consent was obtained from the subjects after explanation of the nature and possible consequences of the study. Institutional Review Board (IRB)/Ethics Committee approval was obtained, and all data for the study were collected and analyzed in accordance with the tenets set forth in the Declaration of Helsinki.
Patients’ clinical data included age at diagnosis, sex, laterality, best corrected visual acuity (BCVA), tumor location, largest basal tumor diameter and thickness determined by type B ultrasonography and survival status. The tumor classification of each patient was recorded according to the seventh edition of the Tumor, Node and Metastasis (TNM) staging system. All tumor specimens stained with hematoxylin and eosin (H&E) were examined by one experienced pathologist to identify the cell type, tumor location and extraocular spread (extra-scleral extension) pathologically. Based on the modifications of Callender’s classification of UM (28), cell types were classified as (1) spindle, epithelioid and mixed; and (2) spindle, and non-spindle types (epithelioid and mixed). In addition, overall survival (OS) time was recorded.
Data Collection From the Validation Cohort
The external validation data came from The Cancer Genome Atlas (TCGA) database and the Surveillance, Epidemiology, and End Results (SEER) database. From the SEER database, we selected all patients diagnosed with “ciliary body melanoma”, “choroid melanoma” or “iris melanoma”, together with the 7th edition of the TNM classification.
Statistical Analysis
Quantitative variables are expressed as the means ± standard deviations. Frequency counts and percentages of participants within each category were calculated for categorical data. Kaplan–Meier curves were used to analyze the survival data. Univariate and multivariate Cox proportional hazards regression was performed to assess factors related to patient survival. Hazard ratios (HRs) and their corresponding 95% confidence intervals (CIs) were derived from Cox models. The proportionality assumption of the Cox model was assessed by testing the correlation of the Schoenfeld residuals and the logarithm of time using the cox-zph test in R.
We searched for predictors of OS that were repeatedly reported in studies or systematic reviews and can be easily ascertained in different settings with various levels of clinical experience. As described herein, 14 variables shown in Table 1 were entered into the selection process. We used a backward stepwise elimination approach with a total of five elimination steps to simplify the model on the basis of the Akaike information criterion. The criterion estimates the fit of each statistical model, penalizes overfitting, and provides a means to select relevant variables that improve the model even if they do not reach the threshold for significance (P<0.05) (29). The continuous variables in the primary regression model were tested for linearity. The global test of any interaction was not statistically significant, so the full prediction model included only main effects. On the basis of the results of the multivariable analysis, a nomogram was formulated by the R package.
Table 1 Demographics and clinicopathologic characteristics of the patients in the training cohort and the validation cohorts.
We assessed the predictive accuracy of the prognostic nomogram in the internal testing set and two external validation cohorts with discrimination and calibration. Discrimination for predicting outcome was evaluated by calculating the Harrell concordance index (C-index) (30), together with the area under the time-dependent ROC curve (AUC). The 95% CIs were generated with bootstrapping methods to account for residual uncertainty. Calibration of the nomogram for 3-, 5-, and 10-year OS was assessed with calibration plots and Brier scores. A 95% perfect calibration is implied by a 45° diagonal line, whereas relevant deviation above or below this line reflects underprediction or overprediction. We determined the clinical significance of using nomograms by performing decision curve analysis (DCA) (31).
We did not impute missing information because there were no missing values in the derivation cohort. To classify individuals into high-risk and low-risk categories, individual predicted risks were converted into binary categories using a cutoff value. The optimal cutoff value was determined using the maximally selected rank statistics from the “survminer” R package. Decision curve analysis was performed using the “rmda” R package. The significance level was set at P<0.05 (two-sided probability). All analyses were performed using R version 4.0.1 (http://www.rproject.org). The results are reported in compliance with the Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD) criteria (Supplementary Table S1).
Results
Screening Process and Clinicopathologic Characteristics of Patients
In total, 318 patients were treated in our hospital in this period. All the patients were Asian. Patients without complete follow-up data (4.4%) and patients treated by proton radiotherapy lacking pathological data (2.8%) were excluded. Thus, 295 patients, all treated by enucleation, with a mean follow-up time of 70.4 ± 47.5 months (range, 1.5 to 207.7 months) were included in the training cohort. The patients included and excluded were comparable (Supplementary Table S2).
In the TCGA database, excluding three patients without TNM staging or pathological data or with contradictory data, 77 patients with a mean follow-up time of 26.2 ± 17.5 months (range, 0.1 to 86.7 months) were included. In the SEER database, patients with incomplete or ambiguous TNM stage or patients without exact pathological data were excluded, and a total of 179 patients with a mean follow-up time of 44.0 ± 21.6 months (range, 1 to 94 months) were analyzed (Figure 1).
Figure 1 Flow chart of the inclusion and exclusion criteria for the training and validation cohorts. (A) Training cohort; (B) TCGA validation cohort; (C) SEER validation cohort. *Patients with the 7th edition of the TNM classification were diagnosed between 2010 and 2015 in the SEER database.
The clinicopathologic characteristics of the patients in the training cohort and the validation cohorts are listed and compared in Table 1. In the three cohorts, our patients showed a younger onset age, lower male tendencies, less non-spindle cell type, less ciliary body involvement and less extra-scleral extension. The corresponding OS curves are shown in Figure 2.
Independent Prognostic Factors in the Training Cohort
By Cox multivariate regression analysis in the training cohort, among all the clinicopathologic variables, the five most relevant variables were screened out as age, largest basal diameter, ciliary body involvement, non-spindle cell type and extra-scleral extension, and accordingly, prediction Model I was established (Table 2).
In the TNM staging system, the largest basal diameter and tumor height were combined as the tumor size. In the TCGA and SEER databases, tumor size was recorded. Thus, to facilitate external validation and maintain consistency, we developed Model II in which tumor size was included instead of the largest basal diameter. Although the tumor size categories did not reach the threshold for significance (P<0.05), the HR markedly increased with tumor size, indicating that it could also serve as a predictor in the prediction model (Table 2, Model II). We compared the discrimination and calibration between Model II and Model I and compared them with the traditional TNM staging system.
Prognostic Nomogram for OS
The nomograms for Model I and Model II incorporating the prognostic factors were established. In the nomograms, each subtype within the variable was assigned a score on the point scale. Summing up the scores and then locating them on the total points scale, we could easily draw a straight line down to the estimated probability of survival at each time point for each patient (Figure 3).
Discrimination, Calibration and Validation of the Nomogram Model
We calculated two validation measures of accuracy: discrimination and calibration. In the training cohort, the C-index and 3- and 5-year AUC for the established Model I and Model II were significantly higher than those of the TNM staging system (Table 3, Figure 4). For the calibration, Models I and II showed the same good results, and both were better than the TNM staging system (Table 3).
Figure 4 ROC curves for the training cohort. (A) The 3-, 5-, and 10-year ROC curves of the training cohort; (B) ROC curves of Nomogram Model I, Model II and the TNM staging system for the training cohort.
Considering that Model I and Model II were comparable in the training cohort and that Model II was more suitable for external verification, we chose Model II to conduct subsequent analyses. In the validation cohorts, Model II showed a higher C-index and 3- and 5-year AUCs than the TNM staging system (Table 3). Its 3- and 5-year Brier scores were lower than those of the TNM staging system, indicating a better calibration of Model II (Table 3). The calibration plots of prediction closely approximated the 45° line, which presented good consistency in the primary cohort for 3-, 5-, and 10-year OS and excellent agreement in the validation cohorts between the nomogram prediction and actual observation for 3-year OS (Figure 5). Thus, we chose Model II as the final prediction model.
Figure 5 Calibration curves of Model II for the training and validation cohorts. Nomogram-predicted OS is plotted on the x-axis; actual OS is plotted on the y-axis. A plot along the 45° line indicates a perfect calibration model in which the predicted OS is identical to the actual outcomes. (A–C) The 3-, 5-, and 10-year calibration curves for the training cohort; (D) The 3-year calibration curve for the TCGA validation cohort; (E) The 3-year calibration curve for the SEER validation cohort.
DCA was adopted to assess the nomogram’s clinical significance. This method offers insight into clinical consequences on the basis of threshold probability, from which the net benefit could be derived (32). The lines of Model II in Figure 6 were all far from both extreme curves and from the TNM line in the training cohort and the validation cohorts, which suggested that the Model II nomogram had superior prediction compared to that of the TNM staging system.
Figure 6 Decision curve analysis (DCA) curves with net benefit score on the vertical axis and high risk thresholds on the horizontal axis for the training and validation cohorts. The net benefit is determined by calculating the difference between the expected benefit and the expected harm associated with each prediction model. The gray line denotes the assumption that all patients had outcome events (death) during follow-up. The dark black line represents the assumption that no patients had outcome events (death) during follow-up. Other curves represent different prediction models. The curve with the highest benefit score at that threshold is the best choice (17). (A) Training cohort; (B) TCGA validation cohort; (C) SEER validation cohort.
Performance of the Nomogram in Stratifying Patient Risk
Based on the maximally selected rank statistics in the training cohort, the optimal cutoff value was determined to be 170. Accordingly, our patients could be divided into two groups: low-risk and high-risk groups. Their OS curves displayed the most significant difference, with a P value <0.0001 (Figure 7A). In the validation cohorts, risk stratification also showed a significant distinction between the two groups (Figures 7B, C).
Figure 7 Risk group stratification of the training and validation cohorts by the cutoff value of 170. (A) Training cohort; (B) TCGA validation cohort; (C) SEER validation cohort.
Discussion
The TNM staging system is used worldwide for UM and integrates the basal tumor diameter, tumor height, tumor location, extraocular extension, lymph glands and distant metastasis and could thus be used as a prediction model. The Cancer Genome Atlas Class groups tumors as A, B, C or D based on their chromosome status and can also be used to predict metastasis and death (15, 33). The class 1 and class 2 gene expression profiles are genetic assessment tools that are widely applied in Western countries, especially in the United States. Based on genetic analysis, the gene expression profile can distinguish patient prognosis with a significant difference in that nearly all metastatic deaths fall into the class 2 category, and it has been validated on multiple independent datasets (13, 34). Another promising estimating model is the Liverpool Uveal Melanoma Prognosticator Online, which combines pathological, clinical and genetic data to provide individual patients with an estimated survival time as well as the risk of metastasis (14). This model has been validated with several international cohorts and can serve as a valuable tool for predicting all-cause mortality (35, 36). In addition, there are several other prognostic models, such as the artificial neural network by Kaiserman (16), the nomogram-SEER model by Zeng (17), the parsimonious model by Damato (18), and the Markov multistate model by Eleuteri (19). The main characteristics of these models are listed in Table 4. Although gene and chromosome testing increase the accuracy of predictions, the high expense and intricate testing technology may limit their wide application in developing regions and countries. Moreover, from Table 4, we found that these models were developed mostly with the Caucasian population.
Our prediction models screened out only five clinicopathologic variables and did not include gene testing or chromosome typing. However, this model still achieves rather high prediction accuracy and showed even better discrimination and calibration than those of the TNM staging system in our patients. This result suggested that these clinicopathologic features, even without the genetic phenotypes, highly represent the inherent characteristics of the disease, further indicating that the gene phenotypes may have been outwardly expressed by these clinicopathologic features in our patients. Thus, the prediction accuracy was not sacrificed when evaluating without genetic testing. Our model was very easy to operate and has low requirements for testing equipment and a low testing cost; hence, it is highly feasible for use by patients from less medically developed countries and regions. Moreover, we conducted validation in the TCGA and SEER databases, and the model showed good discrimination and calibration in these two databases, which indicated that our nomogram model is also suitable for Western patients.
Comparing the baseline of our patients with the Western databases, the clinicopathology and prognosis of our patients were different: our patients had a younger onset age, less male tendencies, less non-spindle cell type, less ciliary body involvement and less extra-scleral extension, and our patients showed a better survival trend. It is inferred that the better prognosis may be related to the inherent genetic features. A study sequenced the whole genome of Chinese patients and found that the GNA11 mutation rate was only 21.5% in those patients, which was much lower than the reported rate of 32.6% in a Caucasian population; instead of BAP1, SF3B1, and EIF1AX, their patients with HIF1A and FOXO1 mutations exhibited worse OS (37). This finding suggests that there may be inherent genetic differences between the two groups of patients.
Based on this model, the risk stratification displayed an obvious distinction between the low- and high-risk groups in both the training and validation cohorts. This result suggests that the model could provide an excellent differential ability in prognosis prediction. Thus, high-risk patients could be easily screened out and closely followed up to detect and treat metastasis earlier. Furthermore, we found that the curves in the training cohort were very flat and smooth, especially the one indicating the low-risk group, while the corresponding curves in the TCGA and SEER cohorts showed a decreasing tendency. This result indicated that this model reflected the prognosis of Chinese patients well, while there may be other independent variables affecting the prognosis of Caucasian patients. Therefore, it may also serve as proof that there are physiological differences between Asian and Caucasian patients.
There are some limitations of this study: 1. Although we used two external databases, the follow-up time of these two databases was limited, so we can only verify the 5-year survival rate instead of the long-term prediction effect of the 10-year survival rate. 2. We only selected patients with complete information, especially in the SEER database, which could have potentially introduced selection bias. 3. Due to the characteristics of the retrospective study, we didn’t compare the efficiency of our model with other UM prediction models. In the future, some prospective studies will be carried out to further verify the efficacy of the present model, and its comparisons with other UM prediction models will be designed as well.
Model II has the following strengths: 1. This model was established based on Asian patients with UM, and also proved applicable to the Caucasian population, thus indicating its wide scope of application. 2. The model requires only a few variables, which are easy to obtain clinically. It could not only provide the estimated probability of survival for each patient, but also differentiate the low- and high-risk patients (with <170 or ≥170 cores) and thus make different follow-up plans accordingly. For example, we could reduce unnecessary invasive examinations and testing expenses for low-risk patients, and give close follow-ups or preventive treatments, such as clinical drug trials, to patients with high risk. In this way, the model is simple to operate and will contribute to rational allocation of medical resources.
In conclusion, we developed a simple and easy-to-promote nomogram model to predict the OS of patients with UM based on our Chinese population, and it was more accurate than the TNM staging system. Our prediction model also showed high discrimination and good calibration when verified by two Caucasian databases, suggesting that the model is suitable for a wide population.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics Statement
The studies involving human participants were reviewed and approved by Ethics Committee of Eye & ENT Hospital of Fudan University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
Conception and design (HY and JQ); Data collection (HY, BX, YB, KX, JGu., RZ, HR, YY, and JQ); Data analysis and interpretation (JGa.); Manuscript preparation (HY, BX, and JQ); Obtained funding (KX). All authors contributed to the article and approved the submitted version.
Funding
Project supported by the Shanghai Committee of Science and Technology, China (Grand No. 20Y11911200). The funding had no role in the design or conduct of this research.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.879394/full#supplementary-material
References
1. Singh AD, Turell ME, Topham AK. Uveal Melanoma: Trends in Incidence, Treatment, and Survival. Ophthalmology (2011) 118:1881–5. doi: 10.1016/j.ophtha.2011.01.040
2. Isager P, Ehlers N, Overgaard J. Prognostic Factors for Survival After Enucleation for Choroidal and Ciliary Body Melanomas. Acta Ophthalmol Scan (2004) 82:517–25. doi: 10.1111/j.1600-0420.2004.00330.x
3. Russo A, Avitabile T, Reibaldi M, Bonfiglio V, Pignatelli F, Fallico M, et al. Iris Melanoma: Management and Prognosis. Appl Sci (2020) 10:8766. doi: 10.3390/app10248766
4. Carvajal RD, Schwartz GK, Tezel T, Marr B, Francis JH, Nathan PD. Metastatic Disease From Uveal Melanoma: Treatment Options and Future Prospects. Br J Ophthalmol (2017) 101:38–44. doi: 10.1136/bjophthalmol-2016-309034
5. Middleton MR, McAlpine C, Woodcock VK, Corrie P, Infante JR, Steven NM, et al. Tebentafusp, A TCR/Anti-CD3 Bispecific Fusion Protein Targeting Gp100, Potently Activated Antitumor Immune Responses in Patients With Metastatic Melanoma. Clin Cancer Res (2020) 26:5869–78. doi: 10.1158/1078-0432.CCR-20-1247
6. Bol KF, van den Bosch T, Schreibelt G, Mensink HW, Keunen JE, Kilic E, et al. Adjuvant Dendritic Cell Vaccination in High-Risk Uveal Melanoma. Ophthalmology (2016) 123:2265–7. doi: 10.1016/j.ophtha.2016.06.027
7. Damato B, Eleuteri A, Taktak AFG, Coupland SE. Estimating Prognosis for Survival After Treatment of Choroidal Melanoma. Prog Retin Eye Res (2011) 30:285–95. doi: 10.1016/j.preteyeres.2011.05.003
8. Broggi G, Russo A, Reibaldi M, Russo D, Varricchio S, Bonfiglio V, et al. Histopathology and Genetic Biomarkers of Choroidal Melanoma. Appl Sci (2020) 10:8081. doi: 10.3390/app10228081
9. Hendrix MJ, Seftor EA, Hess AR, Seftor RE. Vasculogenic Mimicry and Tumour-Cell Plasticity: Lessons From Melanoma. Nat Rev Cancer (2003) 3:411–21. doi: 10.1038/nrc1092
10. Onken MD, Worley LA, Char DH, Augsburger JJ, Correa ZM, Nudleman E, et al. Collaborative Ocular Oncology Group Report Number 1: Prospective Validation of a Multi-Gene Prognostic Assay in Uveal Melanoma. Ophthalmology (2012) 119:1596–603. doi: 10.1016/j.ophtha.2012.02.017
11. Field MG, Durante MA, Anbunathan H, Cai LZ, Decatur CL, Bowcock AM, et al. Punctuated Evolution of Canonical Genomic Aberrations in Uveal Melanoma. Nat Commun (2018) 9(1):116. doi: 10.1038/s41467-017-02428-w
12. Field MG, Durante MA, Decatur CL, Tarlan B, Oelschlager KM, Stone JF, et al. Epigenetic Reprogramming and Aberrant Expression of PRAME are Associated With Increased Metastatic Risk in Class 1 and Class 2 Uveal Melanomas. Oncotarget (2016) 7:59209–19. doi: 10.18632/oncotarget.10962
13. Onken MD, Worley LA, Ehlers JP, Harbour JW. Gene Expression Profiling in Uveal Melanoma Reveals Two Molecular Classes and Predicts Metastatic Death. Cancer Res (2004) 64:7205–9. doi: 10.1158/0008-5472.CAN-04-1750
14. Eleuteri A, Damato B, Coupland SE, Taktak AFG. Enhancing Survival Prognostication in Patients With Choroidal Melanoma by Integrating Pathologic, Clinical and Genetic Predictors of Metastasis. Int J BioMed Eng Technol (2012) 8:18–35. doi: 10.1504/IJBET.2012.045355
15. Vichitvejpaisal P, Dalvin LA, Mazloumi M, Ewens KG, Ganguly A, Shields CL. Genetic Analysis of Uveal Melanoma in 658 Patients Using the Cancer Genome Atlas Classification of Uveal Melanoma as A, B, C and D. Ophthalmology (2019) 126:1445–53. doi: 10.1016/j.ophtha.2019.04.027
16. Kaiserman I, Rosner M, Pe’er J. Forecasting the Prognosis of Choroidal Melanoma With an Artificial Neural Network. Ophthalmology (2005) 112:1608. doi: 10.1016/j.ophtha.2005.04.008
17. Zeng Q, Yao Y, Zhao M. Development and Validation of a Nomogram to Predict Cancer-Specific Survival of Uveal Melanoma. BMC Ophthalmol (2021) 21:230. doi: 10.1186/s12886-021-01968-6
18. Damato B, Eleuteri A, Hussain R, Kalirai H, Thornton S, Taktak A, et al. Parsimonious Models for Predicting Mortality From Choroidal Melanoma. Invest Ophthalmol Vis Sci (2020) 61:35. doi: 10.1167/iovs.61.4.35
19. Eleuteri A, Taktak AFG, Coupland SE, Heimann H, Kalirai H, Damato B. Prognostication of Metastatic Death in Uveal Melanoma Patients: A Markov Multi-State Model. Comput Biol Med (2018) 102:151–6. doi: 10.1016/j.compbiomed.2018.09.024
20. Manchegowda P, Singh AD, Shields C, Kaliki S, Shah P, Gopal L, et al. Uveal Melanoma in Asians: A Review. Ocul Oncol Pathol (2021) 7:159–67. doi: 10.1159/000512738
21. Shields CL, Kaliki S, Cohen MN, Shields PW, Furuta M, Shields JA. Prognosis of Uveal Melanoma Based on Race in 8100 Patients: The 2015 Doyne Lecture. Eye (Lond) (2015) 29:1027–35. doi: 10.1038/eye.2015.51
22. Yue H, Qian J, Yuan Y, Zhang R, Bi Y, Meng F, et al. Clinicopathological Characteristics and Prognosis for Survival After Enucleation of Uveal Melanoma in Chinese Patients: Long-Term Follow-Up. Curr Eye Res (2017) 42:759–65. doi: 10.1080/02713683.2016.1245422
23. Liu YM, Li Y, Wei WB, Xu X, Jonas JB. Clinical Characteristics of 582 Patients With Uveal Melanoma in China. PloS One (2015) 10:e0144562. doi: 10.1371/journal.pone.0144562
24. Toyama S, Tsuji H, Mizoguchi N, Nomiya T, Kamada T, Tokumaru S, et al. Long-Term Results of Carbon Ion Radiation Therapy for Locally Advanced or Unfavorably Located Choroidal Melanoma: Usefulness of CT-Based 2-Port Orthogonal Therapy for Reducing the Incidence of Neovascular Glaucoma. Int J Radiat Oncol Biol Phys (2013) 86:270–6. doi: 10.1016/j.ijrobp.2012.12.022
25. Kashyap S, Venkatesh P, Sen S, Khanduja S, Shrey D, Tinwala S, et al. Clinicopathologic Characteristics of Choroidal Melanoma in a North Indian Population: Analysis of 10-Year Data. Int Ophthalmol (2014) 34:235–9. doi: 10.1007/s10792-013-9821-8
26. Lee CS, Lee J, Choi JJ, Yang WI, Yoon JS, Lee SY, et al. Cytogenetics and Prognosis for Uveal Melanoma in Korean Patients. Acta Ophthalmol (2011) 89:e310–4. doi: 10.1111/j.1755-3768.2010.02099.x
27. Wong W, Sundar G, Chee C, Zhao PS, Rajagopalan R, Gopal L. Clinical Spectrum, Treatment and Outcomes of Uveal Melanoma in a Tertiary Centre. Singapore Med J (2019) 60:474–8. doi: 10.11622/smedj.2019054
28. McLean IW, Foster WD, Zimmerman LE, Gamel JW. Modifications of Callender’s Classification of Uveal Melanoma at the Armed Forces Institute of Pathology. Am J Ophthalmol (1983) 96:502–9. doi: 10.1016/s0002-9394(14)77914-0
29. Bozdogan H. Model Selection and Akaike’s Information Criterion (AIC): The General Theory and its Analytical Extensions. Psychometrika (1987) 52:345–70. doi: 10.1007/BF02294361
30. Pencina MJ, D’Agostino RB. Evaluating Discrimination of Risk Prediction Models: The C Statistic. JAMA (2015) 314:1063–4. doi: 10.1001/jama.2015.11082
31. Fitzgerald M, Saville BR, Lewis RJ. Decision Curve Analysis. JAMA (2015) 313:409–10. doi: 10.1001/jama.2015.37
32. Huang YQ, Liang CH, He L, Tian J, Liang CS, Chen X, et al. Development and Validation of a Radiomics Nomogram for Preoperative Prediction of Lymph Node Metastasis in Colorectal Cancer. J Clin Oncol (2016) 34:2157–64. doi: 10.1200/JCO.2015.65.9128
33. Jager MJ, Brouwer NJ, Esmaeli B. The Cancer Genome Atlas Project: An Integrated Molecular View of Uveal Melanoma. Ophthalmology (2018) 125:1139–42. doi: 10.1016/j.ophtha.2018.03.011
34. Onken MD, Worley LA, Tuscan MD, Harbour JW. An Accurate, Clinically Feasible Multi-Gene Expression Assay for Predicting Metastasis in Uveal Melanoma. J Mol Diagn (2010) 12:461–8. doi: 10.2353/jmoldx.2010.090220
35. Cunha Rola A, Taktak A, Eleuteri A, Kalirai H, Heimann H, Hussain R, et al. Multicenter External Validation of the Liverpool Uveal Melanoma Prognosticator Online: An OOG Collaborative Study. Cancers (Basel) (2020) 12:477. doi: 10.3390/cancers12020477
36. DeParis SW, Taktak A, Eleuteri A, Enanoria W, Heimann H, Coupland SE, et al. External Validation of the Liverpool Uveal Melanoma Prognosticator Online. Invest Ophthalmol Vis Sci (2016) 57:6116–22. doi: 10.1167/iovs.16-19654
Keywords: uveal melanoma (UM), prognostic model, nomogram, survival probabilities, external validation
Citation: Yue H, Xu B, Gao J, Bi Y, Xue K, Guo J, Zhang R, Ren H, Yuan Y and Qian J (2022) A Novel and Easy-to-Promote Prognostic Model for Patients With Uveal Melanoma. Front. Oncol. 12:879394. doi: 10.3389/fonc.2022.879394
Received: 23 February 2022; Accepted: 09 May 2022;
Published: 02 June 2022.
Edited by:
Wenyin Shi, Thomas Jefferson University, United StatesReviewed by:
Leonard Girnita, Karolinska Institutet (KI), SwedenCorrado Spatola, University of Catania, Italy
Copyright © 2022 Yue, Xu, Gao, Bi, Xue, Guo, Zhang, Ren, Yuan and Qian. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jiang Qian, cWlhbmppYW5nQGZ1ZGFuLmVkdS5jbg==
†These authors have contributed equally to this work and share first authorship