CC BY
8
Technology of Informative Fe for Immunosignature Analysis
DOI: 10.17691/stm2020.12.5.02 Received February 10, 2020
A.A. Koshechkin, Assistant, Department of Theoretical Foundations of Informati O.V. Romanovich, Associate Professor, Department of Theoretical Foundations of Leading Engineer, Institute of Applied Mathematics and Computer Science1;
D. Stamate, DSc, Senior Lecturer2;
S.A. Johnston, DSc, Center Director and Professor3;
A.V. Zamyatin, DSc, Head of the Department of Theoretical Foundations of Informatics1;
Director of the Institute of Applied Mathematics and Computer Science1
1National Research Tomsk State University, 36 Lenin Avenue, Tomsk, 634050, Russia;
2Data Science Department of Computing, Goldsmiths, University of London, New Cross, London, SE14 6NW, UK;
3Biodesign Center for Innovations in Medicine, Arizona State University, Tempe, AZ 85281, USA
The main difficulty in practical work with data obtained via immunosignature analysis is high dimensionality and the presence of a significant number of uninformative or false-informative features due to the specific character of the technology. To ensure practically relevant quality of data analysis and classification, it is necessary to take due account of this specific character.
The aim of the study is to create and test the technology for effective reduction of immunosignature data dimensionality, which provides practically relevant and high quality of classification with due regard for the properties of the data obtained.
Materials and Methods. The study involved the use of two normalized data sets obtained from the public biomedical repository and containing the results of immunosignature analysis.
The technology for selecting informative features was proposed within the framework of the study. It consisted of three successive steps: 1) breaking a multiclass task into a series of binary tasks using the "one vs all" strategy; 2) screening of false-informative features is performed for each binary comparison by comparing the values of the median of the sets "one" and "all"; 3) ranking of the remaining features according to their informative value and selection of the most informative ones for each binary comparison.
To assess the quality of the proposed technology for informative feature selection, we used the results obtained after application of classification based on the filtered data. Support vector method that proved itself in the problems of high-dimensional data classification was used as a classification model.
Results. Effectiveness of the proposed technology for informative feature selection was determined. This technology allows us to provide high quality of classification while significantly reducing the feature space. The number of features eliminated in the second step is approximately 50% for each data set under consideration, which greatly simplifies subsequent data analysis. After the third step, when the feature space is reduced to 15 features, the quality of classification by the macro-average F1-score metric is assessed as 98.9% for the GSE52581 dataset. For the GSE52581 dataset, with the feature space reduced to 266 features, the quality of classification by the macro-average F1-score metric is 91.3%.
Conclusion. The results of the work demonstrate the promising outlook of the proposed technology for informative feature selection as applied to the data of immunosignature analysis.
Key words: early diagnosis of diseases; immunosignature; feature selection in the sample; machine learning.
How to cite: Koshechkin A.A., Romanovich O.V., Stamate D., Johnston S.A., Zamyatin A.V. Technology of informative feature selection for immunosignature analysis. Sovremennye tehnologii v medicine 2020; 12(5): 19, https://doi.org/10.17691/stm2020.12.5.02
Introduction
In 2018, in Russia there were 624 thousand patients who were diagnosed with an oncological disease for the first time in their lives: of them, 30.6% had stage I, 25.8% — stage II, 18.2% — stage III, 20.3% — stage IV. In Russia, cancer mortality amounted to more than
293 thousand people in 2018. At the same time, no statistically significant changes in the absolute number of deaths due to malignant neoplasms have been observed over the past 5 years [1].
Cancer treatment efficacy directly depends on timely diagnosis. Early detection of cancer requires effective and patient-specific, easy-to-use, patient-friendly, and
Corresponding author: Alexander A. Koshechkin, e-mail: [email protected]
inexpensive diagnostic methods [2]. The technology of immunosignature analysis based on the idea of human antibody profiling is one of the most promising methods [3]. This technology is based on a microarray, which is a set of peptides with random amino acid sequences that provide a map of immune activity when interacting with human blood serum. There is a wide variety of peptide arrays containing 10 thousand to 330 thousand peptides.
At present, the applicability of various data mining and classification methods for the analysis and interpretation of data obtained via immunosignature analysis is actively studied. To build effective classification models, researchers need relevant and high quality data. The feature space is based on randomly created peptides, therefore not all features are likely to be informative; hence, their selection is one of the most important stages of data analysis. Discarding useless and redundant features not only improves model performance but also facilitates its interpretation [4]. In this regard, each article devoted to investigation of immunosignature data is bound to pay special attention to the stage of selecting informative features.
Work [5] dwells on the applicability of immunosignature analysis for detecting four different pancreatic diseases (cancer and precancerous condition, type 2 diabetes, and pancreatitis). At an early stage, these diseases have similar symptoms, which complicates the diagnosis. The authors used the Student's t-test to select the best features for further analysis. The average classification accuracy amounted to 92%. At the same time, each disease was found to have unique immunological characteristics.
The authors of study [6] demonstrate that the technology of immunosignature analysis has the potential to meet the requirements of a universal test for cancer diagnosis. An intellectual analysis of two data sets of the 6th and 15th grades was carried out. As a result, it was experimentally shown that immunosignature analysis makes it possible to separate different types of diseases with high precision. U-test was used to select informative features.
Work [7] investigates the possibilities of using the technology of immunosignature analysis using the example of a microarray with 330 thousand peptides for diagnosing breast cancer. The main idea of the study was to use the method of Projection to Latent Structures to identify effective data dimensionality. This was supposed to reduce the negative effect of model overfitting and improve object recognition quality. This approach goes against the main idea of immunosignature analysis aimed at finding possible antigens for various diseases. On applying the method of Projection to Latent Structures, the initial feature space is transformed into a new space of latent structures. In this regard, it becomes impossible to interpret the feature space in the context of "antigen-antibody" interaction. U-test for selecting the best features, also used by the
researchers from the previous article, is considered to be alternative to Projection to Latent Structures.
The use of statistical criteria for informative feature selection is an example of using the filtering methods. These methods are characterized by such problems as non-obviousness of the threshold for separating uninformative features and preservation of feature space redundancy. Analysis of data with redundant features generally requires a lot of memory and computing power, and can also cause such undesirable effect as classification model overfitting [8]. At the same time, the origin of data is not taken into account, which can result in non-obvious errors.
The aim of the study is to create and test the technology for effective reduction of immunosignature data dimensionality, which provides practically relevant and high quality of classification.
Materials and Methods
The study used two normalized data sets (GSE52580 and GSE52581) obtained from the public biomedical repository and containing the results of immunosignature analysis in patients with various oncological and infectious diseases as well as healthy individuals representing the control group [9, 10]. Previously, the datasets were transposed to correspond to the tidy data format [11]. The resulting materials are a set of data (a table) of peptide fluorescence intensity values, where peptide names are columns (features), and class labels are rows (samples).
The GSE52580 dataset has the following characteristics: the number of samples — 240; the number of features — 9787; the number of classes — 6.
The number of each class samples is the same in the GSE52580 dataset.
The GSE52581 dataset has the following characteristics: the number of samples — 1516; the number of features — 10,372; the number of classes — 15. The number of each class samples is different in this set (Table 1).
Technology description. The technology for informative feature selection was proposed within the framework of the study. It consisted of three successive steps:
1) breaking a multiclass task into a series of binary tasks using the "one vs all" strategy;
2) screening of false-informative features is performed for each binary comparison by comparing the values of the median of sets "one" and "all";
3) ranking of the remaining features according to their informative value and selection of the most informative ones for each binary comparison.
The first stage of the technology is the application of
/////////////////////^^^^
20 CTM | 2020 | V0l. 12 j N0.5 A.A. Koshechkin, O.V. Romanovich, D. Stamate, S.A. Johnston, A.V. Zamyatin
Ta bl e 1
Description of data set GSE52581
Class The number of samples
Healthy 249
Astrocytoma 166
Coccidioidomycosis 142
Breast cancer 141
Pancreatic cancer 136
Multiple myeloma 112
Lung cancer 107
Mixed oligoastrocytoma 97
Ovarian carcinoma 86
Pancreatitis 82
Recurrent breast cancer 61
Oligodendroglioma 48
Stage IV breast cancer 42
Glioblastoma multiforme 27
Ewing's sarcoma 20
the "one vs all" strategy. In order to select informative features, multiclass tasks are broken into several binary tasks using the "one vs all" or "one vs one" strategy with subsequent selection of the best features for each binary comparison [12]. This study, we use "one vs all" strategy. This decision stems from the special aspect of immunosignature analysis technology consisting in imitation of disease antigens using random amino acid sequence peptides. The purpose of the analysis at this stage is to find a peptide that plays the role of an antigen for a specific disease, i.e. antibodies produced against
the given disease bind to this peptide, while this does not occur in patients from groups with a different diagnosis.
The second stage of the technology is the median filter. Filtering methods are feature ranking methods that assess feature relevance by considering the intrinsic data properties [13]. After ranking, the features with informative value estimated as lower than the threshold value are removed. The resulting subset of features is used for further data analysis (for example, classification via machine learning methods). Filtering methods are easily scalable for high-dimensional data and have low computational complexity, but most of them are applicable only to binary tasks considering each feature separately and ignoring feature dependencies, which can negatively affect subsequent data analysis.
The first filter is a comparison of the median values for each feature for the set "one" and the set "all". The aim of this filter is to remove from the dataset all features for which the median fluorescence intensity value of set "one" is less than the median fluorescence intensity of set "all" as they are considered to be uninformative. The reason for this lack of information value is illustrated by the example of two features from the GSE52580 dataset (Figure 1).
If we look at these two features in isolation from the subject area, both of them will be informative for separating the class "esophageal cancer" from all the others. However, when distribution of fluorescence intensity of the class "esophageal cancer" is to the right of all other classes, this means that antibodies from the blood of patients diagnosed with esophageal cancer bound to the peptide in a greater number than in other classes presented (Figure 1 (a)). The case, when fluorescence intensity distribution is to the left of all other classes, implies that this peptide bound to antibodies of all classes, including "healthy" ones, but excluding "esophageal cancer" (Figure 1 (b)). This means that
Figure 1. Comparison of informative (peptide CSGTMNSEFQNTTRHVYIMS) (a) and false-informative (peptide CSGVFMLSHHQFHPSWYQPN) (b) features for the class "esophageal cancer"
there have been recorded antibodies developed for some other disease somehow uniting people of these groups, but unrelated to the diagnosis of esophageal cancer. All the earlier studies devoted to immunosignature analysis that we reviewed did not cover this aspect in any way and both variants of features were considered as informative. In turn, there are actually no informative features for objects of the class "healthy"; in fact, these are exceptions that did not fall into any other class.
The third stage of the technology is ranking and selection of features. The criterion of symmetric uncertainty (SU) is used as the second filter for estimating features. This criterion for estimating the correlation between both the features and the target variable is an improved version of the information gain criterion [14]. SU values are in the range [0; 1], where 0 means the complete absence of correlation and, consequently, irrelevance of the feature.
SU(X; Y)= 2;/G(X; , v ' H (X) + H (Y)
where /G(X; Y)=H(X)-X(X|Y) — the information gain for the feature and label of class Y; H(X) — the entropy of feature X; H(Y) — the entropy of feature Y.
The next step is to select a subset of informative features. For this purpose, the best features are selected for each binary comparison based on the estimates by the SU method.
Evaluation of informative feature selection efficiency. The main goal of immunosignature analysis is the diagnosis of diseases, which in terms of data analysis is a classification task. In this regard, to assess the quality of the proposed technology for informative feature selection, we used the results of classification obtained after its application.
There are many different classification methods that can be used to accomplish this task. Let us look at the machine learning technique that has already proven its high efficiency in previous studies [15], the support vector machine (SVM), which is based on the construction of a hyperplane maximally separating the classes [16]. Depending on the kernel settings, it is possible to build separating hyperplanes of various kinds. There is no general approach to automatic kernel selection, so this study evaluates efficiency of each of them.
In addition, it is necessary to standardize features in linear models (for example, SVM) due to the following circumstance. One of the most important assumptions when working with linear models whose parameters are estimated by the least squares method is that the residuals of the model are independent (i.e. not correlated) and have normal distribution with the mean value of 0 and some fixed standard deviation (for example, 1). Therefore, the features were standardized in this study according to the following formula:
Z=
X, - X
where Xi is an individual value for the feature; X — mean value for the feature; a — standard deviation for the feature.
There are many different classification quality assessment metrics that are suitable for the presented task. In this study, we use precision, recall, balanced accuracy, and F1-score [17] due to the objectives of the analysis: it is necessary to effectively separate the group of persons with a certain disease from all others, and the presented metrics make the separation possible. These will be computed for each "one vs all" binary comparison, where "one" is a positive class and "all" is a negative class.
Experimental research. Due to the small number of samples in the dataset, cross-validation analysis is required. Thus, the original dataset is divided into 5 approximately equal parts, observing the proportion of classes. At each iteration, 4 parts form a training set, while the 5th part forms a test set with a subsequent change. At each iteration, informative features are found based on training and assessed using the test set. The following is a description of one cross-validation iteration.
1. Evaluation of the informative value of features in the training set using the SU method and applying the "one vs all" strategy.
2. Conversion of false-informative feature estimates to 0 based on the filtering results and comparison of medians.
3. Selection of the best features for each binary comparison (N) until acceptable classification results are achieved.
4. SVM training of the training set, based on the obtained subset of features, and evaluating the classification quality using the test set.
5. Evaluation of the results obtained.
The work was done using the R programming language and libraries available at the CRAN and Bioconductor public repositories.
Results
First of all, let us look at the results of experiments on evaluating the efficiency of various SVM kernels and selecting informative features with the proposed technology for the GSE52580 dataset. Figure 2 shows the classification results for each SVM kernel according to the macro-average F1-score metric.
Information in Figure 2 demonstrates that the best results are shown by the linear, sigmoidal, and radialbasis kernels, though the linear kernel is much superior to its analogs in performance. In this regard, only the linear kernel will be used further.
Figure 3 shows the results of classification with SVM (linear kernel) using the F1-score metric, depending on the number of selected best features for each class.
The figure clearly shows that the best results are achieved when selecting the best three features for
a
/////////////////////^^^^
22 CTM j 2020 j vol. 12 j No.5 A.A. Koshechkin, O.V. Romanovich, D. Stamate, S.A. Johnston, A.V. Zamyatin
Figure 2. Classification results for various SVM kernels depending on the number of the best features (the GSE52580 set)
1.00
0.99
0.98
0.97
0.96
re or 0.95
c
CO 0.94
1
0.93
0.92
0.91
0.90
0.89
each class, except for "healthy" ones (there are no informative features for the "healthy" class objects, in fact, these are exceptions that did not fall into any other class), which means a total of 15 features for 6 classes.
Table 2 shows the results of classification on the test set using only the best three features for each disease. The values in the table are various quality assessment metrics, averaged when doing cross-validation. The ratings of the classification quality presented in the table indicate separability of classes, despite the significant reduction in the feature space.
Now, let us look at the results of experiments on evaluating the efficiency of various SVM kernels and selecting informative features using the proposed technology for the GSE52581 dataset. Figure 4 shows the results of classification using the macro-average F1-score metric for each SVM kernel.
The figure clearly shows the advantage of the linear kernel over others; therefore, only this kernel is used in the further analysis of this set.
Figure 5 shows the results of classification with the F1-score metric, depending on the number of selected best features for each class.
Table 3 shows the results of classification on the test set using only 19 best features for each class, except for the "healthy" (there are no informative features for the "healthy" class objects, in fact, these are exceptions that did not fall into any other class). The values in the table are various quality assessment metrics, averaged
3
4
5
Healthy
Esophageal carcinoma Breast cancer
6 7 8 9 10 Number of the best features
» Brain cancer
► Multiple myeloma
► Pancreatic cancer
Figure 3. Classification results for each class depending on the number of the best features (the GSE52580 set)
Table 2
The results of classification on the test sample based on the best three features for each disease
Metrics
Class Precision Recall F1-score Balanced accuracy
Brain cancer 111 1
Breast cancer 0.976 1 0.988 0.997
Esophageal cancer 111 1
Healthy 0.969 0.993 0.981 0.993
Multiple myeloma 1 0.975 0.987 0.987
Pancreatic cancer 0.968 0.944 0.956 0.969
2
0.90 0.88 0.86 0.84 0.82 0.80 0.78 0.76 0.74 0.72 0.70 0.68 0.66 0.64 0.62
1 2 3 4 5 6 7 8
9 10 11 12 13 14 15 16 17 18 19 20 Number of the best features j Linear kernel Polynomial kernel
Radial-basis function kernel ♦ Sigmoidal kernel
Figure 4. Classification results for various SVM kernels depending on the number of the best features (the GSE52581 set)
19
Astrocytoma Glioblastoma multiforme Breast cancer Healthy
Oligodendroglioma Stage IV breast cancer Ovarian carcinoma Pancreatitis Pancreatic cancer Coccidioidomycosis Recurrent breast cancer Ewing's sarcoma Multiple myeloma Lung cancer Mixed oligoastrocytoma
Figure 5. Classification results for each class depending on the number of the best features (the GSE52581 set)
Table 3
The results of classification on the test sample based on 19 best features for each disease
Class Metrics
Precision Recall F1-score Balanced accuracy
Astrocytoma 0.915 0.903 0.908 0.946
Breast cancer 0.971 0.979 0.975 0.988
Stage IV breast cancer 0.915 0.978 0.942 0.988
Glioblastoma multiforme 0.756 0.767 0.738 0.881
Healthy 0.893 0.968 0.928 0.972
Lung cancer 0.867 0.963 0.912 0.976
Mixed oligoastrocytoma 0.99 0.865 0.922 0.932
Multiple myeloma 0.962 0.757 0.84 0.877
Oligodendroglioma 0.938 0.85 0.887 0.924
Ovarian carcinoma 0.989 0.989 0.989 0.994
Pancreatic cancer 0.958 0.991 0.974 0.993
Pancreatitis 0.973 0.962 0.967 0.980
Recurrent breast cancer 0.786 0.771 0.768 0.881
Ewing's sarcoma 1 1 1 1
Coccidioidomycosis 0.964 0.956 0.960 0.976
///////////////^^^^^
24 CTM f 2020 f V0l. 12 f N0.5 A.A. Koshechkin, O.V. Romanovich, D. Stamate, S.A. Johnston, A.V. Zamyatin
when doing cross-validation. It is clearly seen that high classification quality has been achieved for the vast majority of classes, despite the significant reduction in the feature space.
The obtained classification results are consistent with previous studies in this field. However, the key aspect of the work is the absence of false-informative features in the final feature space, which was not given due regard earlier. This will produce a positive effect on the subsequent analysis and identification of antigens for various diseases.
Conclusion
The main difficulty in practical work with data obtained via immunosignature analysis is high dimensionality and the presence of a significant number of uninformative or false-informative features due to the specific character of the technology. To ensure practically relevant quality of data analysis and classification, it is necessary to take due account of this specific character. The proposed technology for informative feature selection provides high estimates of classification quality while significantly reducing the feature space.
The number of features eliminated in the second step is approximately 50% for each data set under study, which greatly simplifies subsequent data analysis. After the third step, when the feature space is reduced to 15 features, the quality of classification by the macro-average F1-score metric is assessed as 98.9% for the GSE52580 dataset. For the GSE52581 dataset, with the feature space reduced to 266 features, the quality of classification by the macro-average F1-score metric is 91.3%.
The results of the investigation demonstrate the promising outlook of the proposed technology for informative feature selection as applied to the data of immunosignature analysis.
Study funding. This study was not supported by any financial sources.
Conflicts of interest. The authors have no conflicts of interest to declare.
1. Zlokachestvennye novoobrazovaniya v Rossii v 2018 godu (zabolevaemost i smertnost') [Malignant neoplasms in Russia in 2018 (morbidity and mortality)]. Pod red. Kaprina A.D., Starinskogo V.V., Petrovoy G.V. [Kaprin A.D., Starinskiy V.V., Petrova G.V. (editors)]. Moscow: MNIOI im. P.A. Gertsena — filial FGBU "NMITs radiologii" Minzdrava Rossii; 2019; 250 p.
2. World Health Organization. Guide to cancer early diagnosis. World Health Organization; 2017. URL: https://apps.who.int/iris/bitstream/handle/10665/254500/ 9789241511940%20eng.pdf;jsessionid=F414948FB143i 3D7C21E675BA9C8?sequence=1.
V. ficant nosis 10(3):
networks.
3. Stafford P., Halperin R., Legutki J.B., Magee D.M., Galgiani J., Johnston S.A. Physical characterization of the "immunosignaturing effect". Mol Cell Proteomics 2012; 11(4): M111.011593, https://doi.org/10.1074/mcp.m111.011593.
4. Blum A.L., Langley P. Selection of relevant features and examples in machine learning. Artificial Intelligence 1997; 97(1-2): 245-271, https://doi.org/10.1016/s0004-3702 (97)00063-5.
5. Kukreja M., Johnston S.A., Stafford P. Immunosignaturing microarrays distinguish antibody profiles of related pancreatic diseases. J Proteomics Bioinform 2013; S6(1): 1-5, https://doi.org/10.4172/jpb.s6-001.
6. Stafford P., Cichacz Z., Woodbury N.W., Johnston S.A. Immunosignature system for diagnosis of cancer. Proc Natl Acad Sci US A 2014; 111(30): E3072-E3080, https://doi. org/10.1073/pnas.1409432111.
7. Anisimov D.S., Podlesnykh S.V., Kolosova E.A., Shcherbakov D.N., Petrova V.D., Dzhonston S.A., Lazarev A.F., Oskorbin N.M., Shapoval A.I., Ryazanov M.A. Projection to latent structures as a strategy for peptides microarray data analysis. Matematicheskaa biologia i bioinformatika 2017; 12(2): 435-445, https://doi.org/10.17537/2017.12.435.
8. Subramanian J., Simon R. Overfitting in prediction models — is it a problem only in high dimensions. Contemp Clin Trials 2013; 36(2): 636-641, https://doi.org/10.1016/j cct. 2013.06.011.
9. Stafford P., Zbigniew C., Johnston S. An immunosignature system for diagnosis of cancer [Cancer immunosignaturing — test 1]. National Center for Biotechnology Information Search database; 2013. URL: https://www.ncbi. nlm.nih.gov/geo/query/acc.cgi?acc=GSE52580.
10. Stafford P., Zbigniew C., Johnston S. An immunosignature system for diagnosis of cancer [Cancer immunosignaturing — test 2]. National Center for Biotechnology Information Search database; 2013. URL: https://www.ncbi. nlm.nih.gov/geo/query/acc.cgi?acc=GSE52581.
11. Wickham H. Tidy data. J Stat Softw 2014; 59(10), https://doi.org/10.18637/jss.v059.i10.
12. Izetta J., Verdes P.F., Granitto P.M. Improved multiclass feature selection via list combination. Expert Syst Appl 2017; 88: 205-215, https://doi.org/10.1016/j. eswa.2017.06.043.
13. Bommert A., Sun X., Bischl B., Rahnenfuhrer J.,
Lang M. Benchmark for filter methods for feature_____________
high-dimensional classification data. Computational Statistics & Data Analysis 2020; 143: 106839, https://doi.org csda.2019.106839.
14. Shannon C.E. A mathematical ■ Bell System Technical Journal 1 doi.org/10.1002/j.1538-7305.
15. Andryushchenko V.S
Statistical classification ................
reduction of the feature space dime of diseases. Sovremennye tehnolog 14-20, https://d
16. Cortes C., Vapnik Mach BFC
