{'fraction_non_alphanumeric': 0.033925847317105524, 'fraction_numerical': 0.013497380115407574, 'mean_word_length': 5.023771474230923, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 1, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
NAIRscore as a biomarker for the quality of immune response to neoantigens is related with an increased overall survival in multiple myeloma

Neoantigen provides a promising breakthrough in tumor immunotherapy, although only a subset of patients responds well due to the quality of their immune response. However, few biomarkers have been reported to measure the quality of immune response to neoantigens and to predict prognosis of patients with multiple myeloma (MM). Here, we first developed a neoantigen-prediction pipeline starting from outcomes of somatic mutations and gene-expression profiles. Given the expression of some specific marker genes, the human leukocyte antigen (HLA)-I score and the cytolytic score were evaluated respectively to reflect HLA-I molecular expression and CD8+ T/natural killer (NK) cell abundance. According to the process of the immune response to neoantigens, we comprehensively took neoantigen load, cytolytic score, and HLA-I score to construct a neoantigen immune response score (NAIRscore), in which the HLA-I score presented a hazard ratio (HR) of less than 1, while the cytolytic score and neoantigen load presented a HR of greater than 1. Meanwhile, NAIRscore presented a competitive advantage to stratify MM samples. Especially, those exhibiting high NAIRscore correlated with an increased overall survival (OS), echoing the underlying molecular signatures of lower driver-gene mutations and down-regulated immune response. Notably, an online tool based on this study is provided to identify neoantigens and predict OS.


INTRODUCTION
Multiple myeloma (MM) is the second most frequent hematological malignancy as characterized by the clonal expansion of malignant plasma cells within bone marrow, affecting >30,000 individuals per year in the United States. 1,2 In recent decades, many explorations and advances in therapeutics have prolonged the survival for most patients with MM, but MM remains incurable and prone to relapse in the vast majority of cases. [3][4][5][6] The treatment of MM remains challenging, especially for relapsed/refractory patients. Due to clinical and pathological heterogeneity, response and survival in MM treat-ment vary to prognostic risk factors, such as disease stage at diagnosis and cytogenetic abnormalities. 7,8 Therefore, correlating prognostic factors with biological function could enable the development of novel and targeted therapeutic classes.
In recent years, cancer immunotherapy with neoantigen vaccines has been considered a promising strategy for stimulating the body's own immune system to eliminate tumor cells; however, only small subsets of patients within large cohorts favorably respond to the treatment due to the quality of their immune response. 9,10 Neoantigens are a class of short peptides with aberrant expression in tumor cells, which result from somatic mutations and are used as tumor-specific antigens. Therein, neoantigens can bind to human leukocyte antigen (HLA) molecules and be presented on the surface, thereby being recognized as non-self by T cell receptors (TCRs) to elicit a T cellmediated anti-tumor immune response. 11,12 Previously, our team developed a customized proteogenomic workflow for neoantigen prediction and selection, i.e., ProGeo-neo, 13 and we recently updated this tool (v.2.0) for expanding the source of candidate neoantigens. 14 Also, we issued a manually curated database of human tumor neoantigen peptides, i.e., dbPepNeo, 15 which helps screen immunogenic neoantigens from candidates. This database was updated recently to v.2.0. 16 Besides, we have developed additional tools for identification of immunogenic neoantigens and for the subsequent design of candidate peptide vaccines. [17][18][19] Several clinical trials have shown that personalized neoantigen vaccines have the ability to induce CD8+ and CD4+ T cell responses and improve the overall survival of patients with melanoma and glioblastoma, highlighting the considerable potential of neoantigen used in immunotherapy. 20,21 In addition, neoantigen load was correlated with somatic mutations, and both were reported to relate to prolonged survival with checkpoint inhibitors in melanoma and lung cancer. 22,23 High tumor mutational burden (H-TMB) has received FDA approval for selecting suitable patients for immunotherapy; 24 however, little has been reported with regard to low-TMB cancers, such as MM. Studies have shown that the quality of immune response to neoantigens may act as a better biomarker for immunotherapy. 10,25 Here, based on the 516 paired outcomes (i.e., Maf and FPKM files) of 478 patients with MM downloaded from the MMRF-CoMMpass project, we developed a neoantigen-prediction pipeline to identify neoantigens in each tumor. Then, to represent the quality of immune response to neoantigens, we took the neoantigen load, the cytolytic score, and the HLA-I score to construct a neoantigen immune response score (NAIRscore). A high NAIRscore was found to be associated with an increased overall survival. For ease of use by clinicians, based on this study, we have established a friendly web-based tool with URL at http://www.biostatistics.online/MMPrgognosis/.

Higher mutation burden and neoantigen load in recurrent patients compared with primary patients
In this study, based on the MMRF-CommPass project on the GDC data portal, we downloaded 1,092 simple nucleotide variation files (i.e., Maf) and 859 transcriptomic profiles (i.e., FPKM), respectively. These files were derived from whole-exon sequencing (WES) and RNA sequencing (RNA-seq) of CD138+ cells from patients with myeloma, respectively. As shown in Figure S1, after sample selection, a total of 516 matched Maf files and FPKM files derived from 478 patients with MM were retained, including 462 primary samples and 54 recurrent samples (Table 1). In particular, paired primary and recurrent samples from 32 patients with MM were simultaneously obtained.
We first carried out analysis of somatic mutations in 516 Maf files using R package maftools. 26 In 462 primary samples, a total of 27,953 nonsynonymous variations were detected, ranging from 5 to 1,864 variations per sample, with a median of 46 variations ( Figure S2A). Similarly, in 54 recurrent samples, a total of 3,815 nonsynonymous variations were identified, ranging from 4 to 207 variations per sample, with a median of 60 variations ( Figure S2B). Comparatively, more nonsynonymous mutations were found in recurrent samples (Figure 1A). And, in recurrent samples, significantly higher TMB was observed as well ( Figure 1B). When compared with the 33 TCGA tumors cohorts by function "tcgaCompare" of R package maftools, 26 MM expresses as a low-TMB tumor ( Figure S3).
As shown in Figure S1, a workflow for predicting neoantigens starting with the Maf file and the FPKM file was established through referring to the published neoantigen-prediction pipeline. 27 As the nonsynonymous variations in chromosomes may lead to the changes in amino acid sequences of mutated genes in various forms, including single amino acid variation (SAAV), insertion or deletion of amino acids (in-delAA), and frameshift mutation of amino acids (FSAA). Based on the mutated genes and their altered amino acids, the 8-11 amino acid short peptides containing those changed amino acid(s) were computationally cut from the reference amino acid sequence of mutated genes, and then their binding to HLA type was separately predicted for each patient. HLA-A$02:01 was applied to all patients with MM, as this HLA type is most prevalent among White people in the United States, according to the Allele Frequency Net Database (AFND). 28 Statistically, a total of 6,329 neoantigens were identified in 462 primary samples, ranging between 1 and 437 per sample, with a median of 10. In 54 recurrent samples, a total of 887 neoantigens were identified, ranging between 1 and 65 per sample, with a median of 12.
Here, we directly defined the number of neoantigens as neoantigen load in each sample. Comparatively, a significantly higher neoantigen load in recurrent samples was observed when compared with the primary samples ( Figure 1C). Also, a higher neoantigen load was found when comparing those recurrent samples with their paired primary samples ( Figure 1D). Furthermore, when we compared those samples diagnosed at different International Staging System (ISS) stages, no significance was observed on neoantigen load in primary and recurrent samples ( Figures S4A and S4B), suggesting that the production of neoantigens has no significant correlation with ISS stage.
In addition, in both primary ( Figure 1E) and recurrent ( Figure 1F) samples, a significant positive correlation was discovered between the nonsynonymous mutations and neoantigen load. Additional nonsynonymous mutations may be responsible for relapse in recurrent patients. At the same time, in recurrent samples, more nonsynonymous mutations resulted in a relatively higher TMB, but this may also lead to production of more neoantigens. Consistent with the top mutated genes ( Figures S2A and S2B), we found that the neoantigens' peptides were predominantly derived from immunoglobulin heavy-and light-chain genes IRF4 and TP53.
A high NAIRscore was related to an increased overall survival Neoantigens are considered a class of tumor-specific antigens that can bind to HLA moleculars and present on the surface of tumor cells, thereby being recognized specifically by T cells. 12 Since only HLA-I neoantigens were predicted in this study, we took HLA-I molecular expression and CD8+ T cell abundance into consideration. In light of the expression level of some specific marker genes according to a prior study, 29 the HLA-I score and the cytolytic score can be evaluated to represent HLA-I molecular expression and CD8+ T cell abundance in each sample, respectively. Marker genes B2M, HLA-A, HLA-B, and HLA-C were used to evaluate the HLA-I score, while marker genes GZMA, GZMH, PRF1, GNLY, and GZMM were used to calculate the cytolytic score.
Given the process of immune response to neoantigens, we used the cytolytic score, the HLA-I score, and neoantigen load to construct a NAIRscore to measure the quality of immune response to tumor neoantigens by performing a multivariable Cox proportional hazards regression model. As shown in Figure 2A, we observed that the median of NAIRscore significantly stratified samples into high-and low-NAIRscore groups, and the high-NAIRscore group was associated with an increased overall survival (OS) (hazard ratio [HR] = 0.37; log rank test, p = 0.0014). Due to the significant difference in neoantigen load between primary and recurrent samples above, the NAIRscore was separately applied for sample stratification. With one accord, high NAIRscores were observed to relate to an increased OS in both recurrent ( Figure 2B; HR = 0.07, p = 3e-04) and primary samples ( To further investigate the association between NAIRscore and OS in primary samples, we divided those primary samples into four quartiles based on the rank of the NAIRscore ( Figure 2D). Interestingly, those primary samples in the third and fourth quartiles can both be divided significantly into high-and low-NAIRscore groups by the median NAIRscore, respectively ( Figures 2E and 2F), and the high-NAIRscore group exhibited better OS. In contrast, those primary samples in the first and second quartiles cannot both be significantly stratified by the median of their NAIRscore ( Figures 2G and 2H). However, when those primary samples in the first and second quartiles were further differentiated between high-and low-NAIRscore groups by the optimal p value, respectively, we also observed that the high-NAIRscore group exhibited better outcomes ( Figures S5A and S4B). Indeed, the quality of immune response to tumor antigens can play a pivotal role in survival. 10 In addition, we also considered that the NAIRscore has no significant correlation with the ISS stage, as the comparisons presented in Figures S4C and S4D reveal.
NAIRscore as a biomarker has a competitive advantage Considering the observation of significant stratification in recurrent samples and in the third and fourth quartiles of those primary samples above, we merged these samples to re-construct NAIRscore by re-computing a multivariable Cox regression model. Expectedly, these samples were significantly stratified into high-and low-NAIRscore groups, and the high-NAIRscore group was found to be associated with increased OS ( Figure 3A). Consequently, in the subgroups of MM samples, i.e., ISS stage, gender, and age, NAIRscore can be used as an effective biomarker for sample stratification, and a high NAIRscore was associated with an increased OS ( Figures S6A-S6C). In addition, in this fitted Cox model, a competitive concordance index was measured, and the three features were all significant ( Figure 3B). Therein, an HLA-I score with an HR of less than 1 indicated that the higher expression of the HLA-I molecular expression, the more conducive to the patients' OS. In contrast, a cytolytic score and neoantigen load with an HR of greater than 1 indicated that the higher value of cytolytic score and neoantigen load, the worse the patients' OS. Meanwhile, the relative contribution of each variable in this model was estimated ( Figure 3C), suggesting that the three variables all critically contribute to construct NAIRscore.
Furthermore, by performing a 1,000Â bootstrap resampling to calculate the integrated area under the receiver operating characteristic (ROC) curve (iAUC), we further compared the predictive accuracy of NAIRscore with other factors. As shown in Figure 3D, the iAUCs of NAIRscore and ISS stage were comparable, and both were superior to the other features. Interestingly, neoantigen load, TMB, and nonsynonymous mutation were found to yield the same iAUC, suggesting that they had similar predictive accuracy, possibly due to their correlation.
Moreover, combining NAIRscore with three clinical features, i.e., ISS type, gender, and age, we further carried out a prognostic analysis and fitted the SCORE. Uniformly, when the median of SCORE were used to divide those samples into high-and low-SCORE groups, much smaller p values were observed ( Figure 3E). In the forest plot shown in Figure 3F, we observed a larger concordance index, indicating that the addition of clinical features improved the predictive power. Meanwhile, as shown in Figure 3G, NAIRscore was found to make a greater contribution to this model, followed by ISS stage, age, and gender. In addition, a nomogram was constructed for clinicians, illustrating that a quantitative probability can be evaluated in predicting 1-, 2-, and 4-year OS for an individual ( Figure 3H). The calibration curves of the nomogram in predicting 1-, 2-, and 4-year survival rates were also validated ( Figure 3I), suggesting that a good predictive robustness can be achieved using this nomogram. Also, the NAIRscore was observed to contribute the most risk points when compared with others, indicating that the NAIRscore was a more important feature.

Comparison of molecular characteristics between high-and low-NAIRscore groups
To explore the underlying molecular characteristics of different OSs between samples with high and low NAIRscores, we separately compared their somatic mutations and differentially expressed genes (DEGs). Primarily, driver genes in high-and low-NAIRscore groups were identified (Figures 4A and 4B). Meanwhile, genes NRAS, KRAS, and BRAF were considered as drivers in both high-and low-NAIRscore groups. In the low-NAIRscore group, however, we identified more driver genes, especially genes KRAS, TP53, and IRF4 with higher mutational frequency ( Figure 4C), which may explain in part the poor prognosis in the low-NAIRscore group. Compared with wild-type samples, samples with KRAS mutant ( Figure 4D) and samples with TP53 mutant (Figure 4E) all exhibited a decreased survival. Indeed, KRAS mutation was reported to relate to tumor-promoting inflammation and immune modulation, leading to immune escape in the tumor microenvironment. 30,31 The loss or mutation of TP53 resulted in immune evasion by affecting the activity of myeloid and T cells and modulating the immune response during cancer development. 32 The transcription factor IRF4 is required during immune response, including lymphocyte activation and the generation of immunoglobulin-secreting plasma cells. 33 Moreover, the DEGs between high and lowNAIRscores were investigated as well. As shown in Figure 4F, we identified 14 up-regulated genes and 171 down-regulated genes in total (Table S1). In particular, genes GZMA, GZMH, GNLY, and PRF1 marked to the cytolytic score were down-regulated in the high-NAIRscore group. Therefore, lower cytolytic scores were evaluated in the high-NAIRscore group, which corresponded to the cytolytic score with an HR greater than 1 ( Figure 3B). Further, these down-regulated DEGs in the high-NAIRscore group were applied for GO enrichment analysis (Table S2). Consistently, several significant biology processes (BPs) with regard to immune response were enriched ( Figure 4G), such as T cell activation, immuneresponse-activating cell-surface receptor signaling pathway, antigenreceptor-mediated signaling pathway, and regulation of T cell activation. The cell component (CC) and molecular function (MF) also were enriched, including TCR complex, antigen binding, TCR binding, and major histocompatibility complex (MHC) class I protein binding. More specifically, these down-regulated immune response induced by neoantigens may lead to better survival in the high-NAIRscore group.
A web-based tool to identify neoantigens and predict OS for patients with MM Based on the neoantigen-predictive pipeline proposed above and the prognostic model constructed by 516 patients with MM, we developed a web-based tool to identify neoantigens and predict the OS for each patient with MM (http://www.biostatistics.online/ MMPrgognosis/). In this tool, users need to provide two text files from the same patient with MM patient, i.e., a Maf file to reflect somatic mutations in tumor and an FPKM file to represent expression level of genes in tumor. After uploading the two files, those neoantigens in the tumor can be identified within several minutes, and users can directly download results for further exploration, such as personalized neoantigen vaccines. In the meantime, users also can obtain the predictive OS table and risk curve.

DISCUSSION
In this study, an effective neoantigen-prediction pipeline starting from Maf and FPKM profiles was proposed and integrated into a web-based tool, in which neoantigens resulted from somatic mutations can be identified within minutes. It can serve as an alternative when raw sequencing data is not available. Here, consistent with a prior report, 34 we found that larger mutational burden and neoantigens load in recurrent MM samples compared with the primary MM samples. This indicated that the larger mutational burden may be the cause of relapse in recurrent patients, as well as more neoantigen production.
According to the process of immune response to tumor antigens, we treated neoantigen load, the HLA-I score, and the cytolytic score as features to construct the NAIRscore, representing the quality of immune response to tumor antigens in samples. Meanwhile, we observed that neoantigen load and the cytolytic score have HRs of greater than 1, and the HLA-I score has an HR of less than 1. Consistently, neoantigen load has been reported as an HR of greater than 1, and high somatic mutation and neoantigen load were correlated with decreased progression-free survival in MM. 35 These suggested that a higher cytolytic score resulting from neoantigen load did not yield a better prognosis in MM, although a high neoantigen load and tumorinfiltrating CD8+ T cells have been reported to imply a better survival in some solid tumors, such as melanoma and neuroblastoma. 36,37 This might indicate a difference between low-and high-TMB tumors.
In addition, patients whose samples had a high NAIRscore were found to be significantly associated with a better OS, indicating that this is consistent with previous research that a high quality of immune response to tumor antigens was beneficial to prognosis. 10 In particular, we found that the NAIRscore was able to significantly stratify those MM samples with a high NAIRscore. That suggests that a comprehensive consideration of immune response processes may better serve as a biomarker, rather than any one individually, i.e., neoantigen load, tumor-infiltrating CD8+ T cells, and HLA molecular expression. Indeed, in the high-NAIRscore group, we observed fewer driver genes, such as KRAS, TP53, and IRF4, and a down-regulated BP of immune response induced by neoantigens, such as immuneresponse-activating cell-surface receptor signaling pathway and antigen-receptor-mediated signaling pathway. Of course, in the future, NAIRscore needs to be validated in more MM samples, including taking monoclonal gammopathy of undetermined significance (MGUS) and smoldering MM (SMM) into account. Especially, NAIRscore needs to be verified and compared in more tumors, including highand low-TMB tumors.
As we have known, H-TMB has received FDA approvals for selecting patients for promising immune checkpoint inhibitor (ICI) treatment. However, compared with TMB, neoantigen load was observed to be of equal advantage, while NAIRscore presented a greater advantage (Figure 3D). TMB refers to the number of somatic mutations per megabase of interrogated genomic sequence, while neoantigens denote those immunogenic short peptides originated from somatic mutations. To some extent, TMB and neoantigen load are all positively correlated with somatic mutations in tumor. Meanwhile, emerging research showed that the quality of the immune response is more important in tumors, 10 which may explain why NAIRscore exhibits a better advantage.
In summary, we constructed NAIRscore for the first time to evaluate the quality of immune response to tumor neoantigens, in which a high NAIRscore was found to be an effective biomarker to relate with increased OS. A friendly web-based tool was implemented to identify neoantigens and predict prognosis for patients with MM.

Sample collection and selection
In this study, open sequencing files of CD138+ cells from patients with MM were downloaded from the MMRF-CoMMpass project (dbGaP study accession: phs000748) through the GDC data portal (https://portal.gdc.cancer.gov/projects/MMRF-COMMPASS), including WES and RNA-seq ( Figure S1). We obtained 1,092 simple nucleotide variation files (Maf format) derived from WES and 859 transcriptomic profiles (FPKM) derived from RNA-seq. Here, we only retained the matched WES and RNA-seq samples with clinical information. Further, we selected those files sequenced from bone marrow samples from White patients with MM. Eventually, a total of 516 matched Maf files and FPKM files derived from 478 patients with MM were retained, including 462 primary samples and 54 recurrent samples. Therein, the paired primary and recurrent samples from 32 patients with MM were obtained simultaneously.

Sequencing and statistical analysis
As shown in GDC, the tumor cells (CD138+ cells sorted from bone marrow) and normal cells (whole blood cells from peripheral blood) from matched samples were evaluated by WES. Based on the workflow of aliquot ensemble somatic variant merging and masking and human reference genome (GRCh38), the simple nucleotide variations of each tumor sample were separately filtered and saved in Maf format. On the R (v.4.0.3) platform, the mutated classifications and mutated genes in each Maf file were filtered and analyzed by R package maftools (v.2.6.05), 26 the driver genes in different groups were implemented by the oncodrive function with mutations of greater than 8 and adjusted p values of less than 0.05.
In addition, the tumor cells were evaluated by RNA-seq, and the merged FPKM expression matrix was downloaded and used in this study. The FPKM matrix was applied for filtering of candidate neoantigens and evaluation of the cytolytic score and the HLA-I score. 29 The FPKM matrix was used to screen DEGs between groups by performing Wilcoxon rank-sum test and Benjamini-Hochberg (BH) adjustment, in which we set the cutoff at the absolute value of log2, a fold change of greater than 1, and an adjusted p value of less than 0.05.

Neoantigen-prediction pipeline
Referring to the published neoantigen-prediction pipeline, 27 as shown in Figure S1, we constructed a similar workflow for a neoantigen-prediction pipeline starting from Maf file and FPKM file. In this workflow, the mutated genes and their altered amino acids were first filtered from each Maf file. We downloaded the amino acid sequences of all human genes from UniprotKB. Then, based on those mutated genes and their altered amino acids, we can select those short peptides containing the mutated amino acids from their reference amino acid sequences. Only HLA-I neoantigens consisting of 8-11 amino acids were considered. Subsequently, the binding of these short peptides to HLA type was predicted by netMHCpan4.0, 38 and we screened those short peptides with %Rank of less than 2. Eventually, the expression level of genes mapped by those screening short peptides in the paired RNA-seq files was applied for selection, and we selected those peptides with FPKM of greater than 2.
Somatic mutations in chromosomes may lead to the changes in amino acid sequences of mutated genes in various forms, including SAAV, indelAA, or FSAA. For SAAV, we first found the mutant site in the amino acid sequence of the reference gene and replaced it with the mutant amino acid. Then, we can theoretically cut the sequence into 8-11 amino acid short peptides containing the mutated position. Similarly, for indelAA, we found the mutant site in the amino acid sequence of the reference gene and inserted or deleted the mutant amino acid(s). Also, those 8-11 amino acid short peptides containing the mutant amino acid(s) can be cut from the sequence. For FSAA, all mutant sites in the amino acid sequence of the reference gene were first identified, and the sequence containing these mutation sites was cut into 8-11 peptides.
Due to lack of sequence files (e.g., fastq), we cannot accurately predict HLA types for each patient with MM. In this study, as only White patients with MM were selected, according to the AFND, 28 we found that HLA-A$02:01 is the major HLA type among White people in the United States. Therefore, those short peptides were applied to predict the binding to HLA-A$02:01 by using netMHCpan4.0.
Evaluation of cytolytic score and HLA-I score A previous report 29 showed that the cytolytic score and the HLA-I score in hematologic tumors can be evaluated by the geometric mean of expression of several specific marker genes as follows: in which i denotes the i-th marker genes. Accordingly, the marker genes GZMA, GZMH, PRF1, GNLY, and GZMM were used to calculate the cytolytic score, which is positively related to the cytolytic activity of CD8+ T/natural killer (NK) cells. The marker genes B2M, HLA-A, HLA-B, and HLA-C was used to detect the HLA-I score, which is significantly associated with the presentation of HLA-I neoantigens.
NAIRscore construction and prognostic analysis In this study, we directly treated the number of neoantigens as neoantigen load in each sample. Taking neoantigen load, the cytolytic score, and the HLA-I score as features, a multivariable Cox proportional hazards regression model was constructed using the R package survival (v.3.2-13) and survminer (v.0.4.9), in which the coefficient (weight) of each feature can be generated in view of the importance of this feature. Meanwhile, NAIRscore can be constructed as follows: in which i denotes the i-th feature, and Feature and Coef represent its value and coefficient in the fitted model, respectively. The NAIRscore was designed to reflect the quality of immune response to neoantigens in each sample. We applied the median of NAIRscore to stratify samples into high-and low-NAIRscore groups.
Further, when combining NAIRscore and three clinical features, i.e., ISS stage, age, and gender, the scoring system was further upgraded as follows: in which i denotes the i-th feature. When the continuous features (i.e., NAIRscore, age) were applied, Feature and Coef represent its value and coefficient, respectively. When the feature was gender, "female" was treated as reference, and then "male" was applied by its Coef. When the feature was ISS stage, "ISS-I" was treated as reference, and then "ISS-II" and "ISS-III" were individually applied by their Coef.
The Kaplan-Meier survival curve was implemented by R package survminer (v.0.4.9), and the log rank test was used to compare the OS between high-and low-NAIRscore groups. The nomogram and calibration plots were visualized by using the R package rms (v.6.2-0).
To compare the predictive accuracy of the different variables or models, the iAUC with 100,00Â bootstrap resampling was performed and analyzed. Referring to two previous studies, 39

Establishment of online web tool
An easy-to-use web tool starting from the paired Maf and FPKM files was developed, which can allow users to identify neoantigens in tumor and predict a patient's OS. To implement an interaction interface and improve user experience, the front-end web framework was built by Bootstrap and Jquery, and the back-end was implemented by PHP.

Supplemental information
NAIRscore as a biomarker for the quality of immune response to neoantigens is related with an increased overall survival in multiple myeloma    . The individual median of NAIRscore in each prognostic model was used to divide samples into high-and low-NAIRscore groups, and then the HR, 95% CI, and P-value were determined by univariate Cox analysis. :   Table S1. The DEGs between high-and low-NAIRscore groups. Significance was determined by Wilcoxon test, and P values were adjusted by BH method.  Figure 4G.