 Research
 Open Access
 Published:
Drug response prediction model using a hierarchical structural component modeling method
BMC Bioinformaticsvolume 19, Article number: 288 (2018)
Abstract
Background
Componentbased structural equation modeling methods are now widely used in science, business, education, and other fields. This method uses unobservable variables, i.e., “latent” variables, and structural equation model relationships between observable variables. Here, we applied this structural equation modeling method to biologically structured data. To identify candidate drugresponse biomarkers, we first used proteomic peptidelevel data, as measured by multiple reaction monitoring mass spectrometry (MRMMS), for liver cancer patients. MRMMS is a highly sensitive and selective method for proteomic targeted quantitation of peptide abundances in complex biological samples.
Results
We developed a componentbased drug response prediction model, having the advantage that it first combines collapsed peptidelevel data into proteinlevel information, facilitating subsequent biological interpretation. Our model also uses an alternating least squares algorithm, to efficiently estimate both coefficients of peptides and proteins. This approach also considers correlations between variables, without constraint, by a multiple testing problem. Using estimated peptide and protein coefficients, we selected significant protein biomarkers by permutation testing, resulting in our model for predicting liver cancer response to the tyrosine kinase inhibitor sorafenib.
Conclusions
Using data from a cohort of liver cancer patients, we then “finetuned” our model to successfully predict drug responses, as demonstrated by a high area under the curve (AUC) score. Such drug response prediction models may eventually find clinical translation in identifying individual patients likely to respond to specific therapies.
Background
Liver cancer (hepatic cancer), is predominantly found in the tissue parenchyma, and is thus known as hepatocellular carcinoma (HCC), the most common form of liver cancer in adults. HCC can exert different growth patterns from one tumor to the next [1, 2]. In Eastern Asia, HCC is the thirdmost common form of cancer, and the secondleading cause of cancer death, with a worldwide total of 600,000 deaths each year [3, 4]. However, as many treatment methods have been developed for treating HCC, overall, these have shown little benefit in improving patients’ prognosis [5]. More efficient treatment of HCC may lie in “personalized medicine,” i.e., tailoring therapies for individual patients [6]. Such ability to classify HCC patients, with therapies optimized for specific stage and growth patterns, would reduce time and cost, and likely prolong survival.
Toward that objective, accurate prediction models are essential. Historically, methods of building cancer prediction models were based on the classification methods of linear/logistic regression, support vector machine, or random forest [7,8,9]. While these models are effective for prediction, they do not consider any structural or hidden biological data, making it difficult to derive more meaningful biological interpretation.
Here, we built a drug response prediction model, by identifying candidate protein biomarkers, via multiple reaction monitoringmass spectrometry (MRMMS) technology. MRMMS is a targeted proteomics technology that is highly selective and sensitive for quantitating targeted proteins or peptides in biological samples [10, 11]. MRMMS can also measure several hundred protein targets per sample, simultaneously, generating consistent, precise, and reproducible datasets [12]. Consequently, MRMMS holds high potential for biomarker discovery. Unlike other protein data, MRMMS data is hierarchically structured.
Following mass spectrometry, our MRMMS data consisted of 231 peptides, from 124 proteins, with each protein containing ≥1 peptide. While classical methods for prediction model building only select the best peptides as variables, to optimize prediction performance, these do not consider any biological relationship between peptides and proteins.
In this study, we built a drug response prediction model, using a componentbased structural equation modeling method, based on the biological structure of MRMMS data (e.g., peptide to protein). Structural equation modeling (SEM) is used to analyze the structural relationship between unobserved (latent) variables and observed variables. SEM can be classified as factor based SEM and component based SEM. Confirmatory factor analysis (CFA) and partial least squares path modeling (PLSPM) analysis are the most popular methods of factor based SEM and component based SEM, respectively [13]. Our proposed model is based on generalized structured component analysis (GSCA) [14], resembling our earlier derivation of a pathwaybased approach. That analysis, using hierarchical components of collapsed rare variants (PHARAOH), uses a hierarchical structure of pathways and genes [15].
Using latent variables, we can collapse multiple peptides into a structured form of proteins that they comprise of, providing more feasible biological explanations of the results. In addition to hierarchical structure, we further showed that HisCoM effectively cover proteinlevel analysis, taking all peptides into account simultaneously. Moreover, for real biological data analysis, using MRMMS, we discovered possible protein biomarkers associated with patients’ response to the multiple tyrosine kinase inhibitor sorafenib (Nexavar®) [16]. Sorafenib is known as effective and safe drug for recovering liver cancer (hepatocellular carcinoma) patients not only AsianPacific region but also in the North American region [17, 18]. Using these protein biomarkers, we then evaluated the performance of our drug response prediction model. Additionally, we compared the performance of our prediction model, using area under the curve (AUC) scores, to performances by generalized linear models of logistic regression, without ridge parameters, and logistic regression, with ridge parameters. Furthermore, through extensive simulation studies, we compared the performance of our proposed method with other logistic regression methods. For hierarchical structuring, in this case, for proteins with multiple peptides, our HisCoM was shown to perform better than logistic regression, as assessed by AUC scores.
Methods
Preparing samples and materials
Hepatocellular carcinoma (HCC) patient serum samples (n = 115) were collected at Seoul National University Hospital, from 2013 to 2015 [19]. Upon diagnosis of liver cancer, patients were placed on a regimen using the tyrosine kinase inhibitor Sorafenib (Nexavar®, Bayer, Inc., Whippany, NJ, USA). Patients’ tumor sizes were first examined immediately following HCC diagnosis, at the start of hospital admission. Six weeks after first diagnosis (sufficient time to see a response), patients’ tumors were again measured, by contrastenhanced computed tomography or magnetic resonance imaging, and staged according to the standardized Modified Response Evaluation Criteria in Solid Tumors (mRECIST) [20]. After the second examination, patients were divided into two groups, based on positive and negative drug responses. The positive drug response group consisted of patients with complete response (CR), partial response (PR), or stable disease (SD), according to mRECIST [20]. CR and PR responses were diagnosed when the tumor size was reduced after 6 weeks. Also, SD was diagnosed when the size of the tumor remained stable, from the first to second visit. On the other hand, the negative drug response group consisted of patients with progressive disease (PD), wherein the size of their tumors increased, from first diagnosis to 6 weeks later. The study protocol was approved by the Institutional Review Board of Seoul National University Hospital (IRB No. 0506–150005), and written, informed consent was obtained from each patient or legally authorized representative.
Among all 115 patients (101 men and 14 women), 40 patients (37 men and 3 women) were grouped into the positive drug response group, and 75 (64 men and 11 women) were grouped into the negative drug response group. From each patient’s serum, data for 231 peptides was generated by multiple reaction monitoring mass spectrometry (MRMMS), a highly sensitive and selective method for targeted quantitation of peptide abundances, in complex biological samples [21]. Here, the 231 peptides can represent 124 proteins. Since the MRMMS technique measures the quantity of targeted peptides in patients’ serum, we used the log2transformed ratio of light peptide intensity to heavy peptide intensity. Light peptide intensity represented the quantity of peptides from specific patient’s blood, as measured by MRMMS, while heavy peptide intensity represented the quantity of artificially built, same sequences as the light peptides, but using heavier isotope elements, also as measured by MRMMS. The software Skyline was used to measure the intensity of light and heavy peptides by MRMMS [22]. Demographic information, such as age and sex, were also available. The range of age varied from 34 to 84, with 101 male and 14 female samples.
Constructing the drug response model
The overall schematic procedure is shown in Fig. 1. At the beginning, we selected protein level biomarkers by HisCoM method with 1000 permutation test, for possible prediction of sorafenib response, using MRMMS data. Second, we constructed prediction models, via a componentbased structural equationmodeling method. Finally, we evaluated the constructed drug response prediction models’ performances, by AUC scores.
An example of our proposed drug response prediction model is shown in Fig. 2. This model combines collapsed peptidelevel MRMMS data into proteinlevel information, and efficiently estimates both peptide and protein coefficients. In this example, two proteins were involved (K = 2), and each protein consisted of two or three peptides (T_{k} = 2 and 3). Weight (w) and path coefficients (β) were estimated using alternating least squares [23].
Here, suppose that there are K proteins, and the k^{th} protein contains T_{k} peptides, for k = 1,…,K. To estimate parameters, the following penalized log likelihood function was maximized. The y_{j} represent the drug response group, based on mRECIST: y_{j} = 0 for a negative response, and y_{j} = 1 for a positive response. Let y_{j} is distributed independently with a mean of μ_{j} = E[y_{j}] and η_{j} is defined as η_{j} = g(μ_{j}) by a logit link function g. Then, we can derive a penalized log likelihood function with dispersion parameter δ and canonical parameter γ_{i} as following:
Here, λ_{prot} and λ_{pep} are the ridge parameters for proteins and peptides represent as “tuning” parameters, respectively: one for the peptides within a protein and the other for the proteins themselves.
Let w_{k} = [w_{k1},⋯,\( {w}_{k{T}_k} \)], β = [β_{0,} β_{1},⋯, β_{K}], F = [f_{1},⋯, f_{N}], and f_{j} = [1, f_{j1},⋯, f_{jK}] where \( {f}_{jk}=\sum \limits_{t=1}^{T_k}{x}_{jk i}{w}_{ki} \). We define x_{jki}as the quantity of i^{th} peptide of the k^{th} protein in sample j. The w_{ki} as a weight coefficient of i^{th} peptide of the k^{th} protein. Also, the β_{k} as a path coefficient of k^{th} protein. Maximizing the eq. (1), via iteratively reweighted least squares, is identical to minimizing the following penalized least squares eq. (2):
with respect to w_{k} and β [15, 24]. Here, V is an N by N diagonal matrix with elements \( {v}_j=\frac{{\left({\partial}_{\mu_j}/{\partial}_{\eta_j}\right)}^2}{\tau_j} \). τ_{j} is the variance function evaluated at μ_{j}. The z is the adjusted response variable and an N × 1 vector with elements \( {z}_j={\eta}_j+\frac{\left({y}_j{\mu}_j\right)}{v_j} \) [25].
After estimating the w_{ki} and β_{k} coefficients, we constructed a drug response prediction model for π_{j} = P[y_{j} = 1] = μ_{j}, as follows, after standardization of x_{jki} The coefficients of age and sex are also estimated by maximizing the loglikelihood function simultaneously while penalizing the coefficients of peptides. In our final prediction model, the beta coefficients of age and sex are fixed across the individuals.
j: individual samples (j = 1, ⋯,115)
k: proteins (k = 1, ⋯,124)
i: peptides (i = 1, ⋯,231)
When our final drug response prediction model was constructed, we evaluated its performance by area under the receiver operating characteristic curve (AUC) score, based on the training, validation, and test sets. The approach for separating the training set from the validation and test sets, is depicted in Fig. 3. First, we then randomly selected 39 out of 115 samples (35 men and 4 women) as a test set, excluded from the modelling process, while assessing the remaining 76 samples. This test set will be used to measure the performance of our final drug response prediction model. The ratio of positive responses to negative responses was sustained (14 positive responses and 25 negative responses). The range of age in test set was distributed from 41 to 84. The remaining 76 samples were randomly (without replacement) divided into training and validation set. For a fair comparison, the ratio of positive responses to negative responses was retained (13 positive responses and 25 negative responses). The concept of a sample separation process was based on a previously developed intraductal papillary mucinous neoplasm (IPMN) patient prediction model [26].
Lastly, our drug response prediction model was compared to the generalized linear model with a binary response (GLM), and the generalized linear regression with a binary response via ridge parameter (GLMwR) methods. All the analyzes were calculated and computed via software R (Version R 3.2.3) [27].
Simulation design
For the simulation study, we designed two models: the first model composed of two significant proteins and the second model with both significant and nonsignificant protein in the presence of a hierarchical structure of MRMMS data (e.g., peptide to protein). Let the first simulation model contain JCHAIN and RBP4 with parameters estimated by HisCoM. Note that JCHAIN was a significant protein (pvalue: 0.0142), with 3 peptides, and RBP4 was also a significant protein (pvalue: 0.0031), with 2 peptides. The simulation model is given by
For the Simulation model 2, we assume the true model contains RBP4 and APOA1, with parameters estimated by HisCoM. Note that RBP4 was a significant protein (pvalue: 0.0031), with 2 peptides, and APOA1 was a nonsignificant protein (pvalue: 0.4794), with 7 peptides. The second simulation model is given by
In this case, x_{j, 1} represents the j^{th} individual’s peptide data (x_{j,1}, x_{j,2}, ⋯, x_{j,9}). From the estimated βs and ωs, derived from the data, we estimated π_{1}, π_{2},⋯, π_{115}. Then, the responses were generated from the Bernoulli distribution B(1, π_{j}), for j = 1, 2, ⋯, 115. We then constructed HisCoM, GLM, and GLMwR drug response prediction models, using MRMMS peptide data (x_{j,1}, x_{j,2},⋯, x_{j, 9}), to generate response variables. For each simulation model, we measured the AUC score. Using the same estimated values of π_{1}, π_{2},⋯, π_{115}, we repeated the whole process 1000 times, and obtained 1000 AUC scores for each of the HisCoM, GLM, and GLMwR models. We then calculated the mean of the 1000 AUC scores, based on those models.
Results
Biomarker discovery for the drug response prediction model
To evaluate our model, at the beginning, we randomly selected 39 out of 115 samples, as a separate, test set, to evaluate the overall performance of the final drug response prediction model. We performed crossvalidation analysis using remaining 76 samples. The dataset was randomly divided into training/validation sets (38 samples for each set). From the training data set, the significant proteins were selected based on pvalues. Then the prediction model was build, and its AUC score was computed from the validation set. We repeated this crossvalidation 100 times. Through 100 crossvalidation, we evaluated whether the significant proteins were selected repeatedly by HisCoM. Also, using the estimated path coefficients by training set, we evaluated the performance of the prediction model from the validation set.
The significances of the protein path coefficients were then determined, using a 1000 permutation test, for each replicate. The permutation test was performed by shuffling drug response variables randomly across subjects while retaining the ratio of positive response to negative response and then estimating the path coefficients. Using path coefficients estimated by 1000 permutation test, we construct distribution of each protein’s path coefficient. Through comparing the path coefficient value of the original data with those from the permuted data, pvalues were computed for each protein, with the significant (P < 0.05) ones being selected for further analysis.
Figure 4 shows the null distribution of path coefficients (β_{k}), derived from 1000 permutation tests. Figure 4a show the case of nonsignificant protein, while the Fig. 4b does the case of significant protein. The red dots indicate the estimated path coefficients from the data. In our analysis, the path coefficients (β_{k}) of six proteins were significant.
During the processing of the training/validation sets, with 76 samples, we repeated this process 100 times, to check the consistency of possibly significant proteins. Since our method uses two ridge parameters, we defined the same tuning parameter values as 10 for peptides and proteins, for computational efficiency. As a result, we selected the top 6 significant proteins (APOC4, CD163, CD5L, JCHAIN, SERPING1, and RBP4), which were repeatedly selected as significant by the process of 100 replications. We noted that these six proteins were previously identified as possible proteomic biomarkers, for hepatocellular carcinoma [28,29,30].
We then repeated the process once more, with only those 6 proteins, as MRMMS data, for more accurate estimation of pvalues and path coefficients, for the drug response prediction model. We next calculated pvalues and path coefficients. In Table 1, pvalues are shown for the six selected proteins.
Using the selected 6 proteins, we constructed a drug response prediction model, with estimated w and β values. We also constructed different prediction models, limiting the number of proteins. In this case, we constructed models using 1 of the 6 proteins, 2 of the 6, 3 of the 6, and all six. For all these models, age and sex were considered as covariates (see eqs. 3 and 4, below). All the analyzes were calculated and computed via software R (Version R 3.2.3) [27].
Model evaluation by AUC results
With the selected proteins, we first constructed a sorafenib drug response prediction model given in eq. (4), using HisCoM. The performance of the drug response prediction models was measured by AUC scores. In this case, the numbers k and i varied, depending on the number of proteins in the model.
Table 2 shows the AUC score of our single protein prediction model, compared to a corresponding the generalized linear model with a binary response (GLM), and the generalized linear model with a binary response via ridge parameters (GLMwR). The performance of the single protein prediction models showed similar AUC scores, across all three different statistical methods, while the AUC scores, for each individual protein, varied from 0.60 to 0.90. Figure 5 provides a visual comparison of AUC scores, by different single protein statistical models.
The prediction model, using 2 of the 6 proteins, had higher AUC scores, compared to the single protein models, across all three different statistical methods. Table 3 shows the AUC scores for the models with 2 of the 6 proteins. The AUC scores across each statistical model varied from 0.73 to 0.95, higher than those for the single protein prediction models. Figure 6 shows a visual comparison, of AUC scores, by different twoprotein statistical models. The HisCoM’s AUC score was similar to those of GLM and GLMwR, but had higher performance or lower performance, depending on the combination of proteins. The best performing protein combination, across each HisCoM, GLM, and GLMwR statistical model, was the combination of SERPING1 and JCHAIN.
Similarly, most of the prediction models, with 3 of the 6 proteins, scored higher than 0.9 AUC, using all three modeling methods. Table 4 shows the AUC scores for each model using exhaustive combinations of the three proteins, varying from 0.86 to 0.95. Figure 7 shows a visual comparison of AUC scores, by different 3protein statistical models.
Using all 6 proteins, we also constructed HisCoM drug response prediction models, with estimated ω and β as the covariates age and sex, respectively. In Fig. 8, the AUC score for our HisCoM model was 0.96, using the validation set. At first, we tried to compare our prediction model to the generalized linear model with a binary response. However, the latter had a convergence problem, due to high correlation among peptides, as shown in Fig. 9. To resolve this problem, we fit the logistic regression model with a ridge penalty, using the “GLMNET” R Package. The result is shown in Fig. 8, and the AUC score for the generalized linear model with a binary response via ridge parameter (GLMwR), was 0.949, using the same validation set. As a result, our HisCoM had a slightly better AUC score (0.96), compared to that of GLMwR (0.949).
Simulation results
The performance of Simulation model 1 (JCHAIN + RBP4) results, the mean AUC scores of HisCoM, GLM and GLMwR are shown in Table 5. The mean AUC score of HisCoM was 0.8362. Figure 10 shows the range of 1000 AUC scores, as depicted by box plots, with respect to each statistical method. It shows that the HisCoM performed better than others. The performance of Simulation model 2 (RBP4 + APOA1) results, the mean AUC scores of HisCoM, GLM and GLMwR, are shown in Table 6. These results show that HisCoM had the highest mean AUC score, compared to the two other statistical methods. The mean AUC score of model 2 by HisCoM was 0.7270, while the means of the other statistical methods were less than 0.7. Also, Fig. 11 shows the range of 1000 AUC scores, as depicted by box plots, with respect to each statistical method model.
In summary, both simulation model results show that HisCoM was the best performing model, compared to the other statistical methods, when there exists a hierarchical structure of MRMMS data (e.g., peptide to protein).
Discussion
In this study, we developed a prediction model for tumor response to the multiple tyrosine kinase inhibitor sorafenib (Nexavar®), for liver cancer patients [16], using a componentbased structural equation modeling method. We used HisCoM to construct the model, for Korean hepatocellular carcinoma (HCC) patients, using MRMMS proteomic data, including some demographic variables. HisCoM fit the whole data set at once. In this case, we measured 231 peptides’ weights, and 124 protein’s path coefficients, to the drug response variable, all at once. The positive or negative drug response variables were defined by tumor responses according to mRECIST [20]. Thus, this model can be used for largescale, structured data, with marker selection (as well as model building), simultaneously. The second, and most important advantage of HisCoM, is that it generates latent variables, which are not directly observed, while collapsing other (observed) variables. For example, our HisCoM combines several collapsed peptides’ MRMMS data, into several proteins, as latent variables. Unlike other classical methods, such as linear/logistic regression, support vector machine, and random forest, our HisCoM approach considers peptidetoprotein computational structure, and peptidetoprotein biological structure. In the analysis, we found 6 possible protein biomarkers that significantly associate with sorafenib drug response. On the other hand, other classical prediction modeling methods do not consider structure of biological information. Using peptidelevel data, we found significant proteins, as possible biomarkers, for building a sorafenib response prediction model for liver cancer patients. The overall work flow, with our statistical analysis, using a HisCoM schema, can be accurately applied not only to other cancers, but also to most any largescale structured data.
Conclusions
From possible biomarker selection, to AUC performance test scores, through a modelbuilding process, we compared the performance of our model, constructed using a HisCoM method, to other classical statistical methods such as generalized linear models, using logistic regression (alone) or logistic regression with ridge parameters. For possible drug response biomarkers, 6 significant proteins were statistically selected, using pvalues, as computed by permutation tests: APOC4 (pvalue: 0.0061), CD163 (pvalue: 0.0112), CD5L (pvalue: 0.0031), JCHAIN (pvalue: 0.0102), SERPING1 (pvalue: 0.0142), and RBP4 (pvalue: 0.0031). All six of these proteins were previously reported as possible biomarkers for hepatocellular cancer (HCC) [31,32,33]. Of these, CD5L is the bestknown HCC biomarker [28]. For the single protein model, using HisCoM, the AUC scores varied from 0.60 to 0.90, depending on the specific protein. For modeling combinations of 2 of the 6 proteins, by HisCoM, the AUC scores varied from 0.73 to 0.95, showing increased performance, compared to singleprotein prediction models. On the other hand, AUC scores varied from 0.86 to 0.95, for the 3protein model, by HisCoM. Finally, using all six of the abovementioned proteins in the model, we successfully constructed a drug response prediction model using 1, 2, 3, or all sixprotein models. By comparing our model’s performance with the generalized linear model with a binary response via ridge penalization, the performance of our sixprotein HisCoM prediction model was an AUC score of 0.96, slightly better than the generalized linear model with a binary response via ridge parameter, for the 6protein panel, with an AUC score of 0.949 (Fig. 8). Thus, both the HisCoM and GLMwR methods had high AUC scores. Overall, we conclude that our model was marginally superior to the classical model types.
For future research, we can apply this overall prediction modelbuilding approach, using HisCoM, to other cancer data especially derived from MRMMS platform. Since these potential biomarkers were identified in patients’ serum, these could be obtained by a minimally invasive procedure (e.g., as compared to biopsies, lumbar puncture, etc.). Such models could ultimately assist physicians in discerning which therapies might be effective, for individual patients.
Abbreviations
 AUC:

Area under the curve
 CR:

Complete response
 GLM:

Generalized linear model
 GLMwR:

Generalized linear model with ridge parameters
 GSCA:

Generalized structured component analysis
 HCC:

Hepato cellular carcinoma
 HisCoM:

Hierarchical structural Component Models
 mRECIST:

Modified Response Evaluation Criteria in Solid Tumors
 MRMMS:

Multiple reaction monitoring mass spectrometry
 PD:

Progress disease
 PHARAOH:

Pathwaybased approach using HierArchical components of collapsed RAre variants Of Highthroughput sequencing data
 PR:

Partial response
 SD:

Stable disease
References
 1.
Asnacios A, Fartoux L, Romano O, Tesmoingt C, Louafi S, Mansoubakht T, Artru P, Poynard T, Rosmorduc O, Hebbar M. Gemcitabine plus oxaliplatin (GEMOX) combined with cetuximab in patients with progressive advanced stage hepatocellular carcinoma. Cancer. 2008;112(12):2733–9.
 2.
DeVita VT Jr, Lawrence TS. DeVita, Hellman, and Rosenberg’s Cancer: principles & practice of oncology, vol. 10e; 2009. p. 696–714.
 3.
Ferenci P, Fried M, Labrecque D, Bruix J, Sherman M, Omata M, Heathcote J, Piratsivuth T, Kew M, Otegbayo JA. Hepatocellular carcinoma (HCC): a global perspective. J Clin Gastroenterol. 2010;44(4):239–45.
 4.
carcinoma APWPoPoH. Prevention of hepatocellular carcinoma in the Asia–Pacific region: consensus statements. J Gastroenterol Hepatol. 2009;25:657–63.
 5.
Lin S, Hoffmann K, Schemmer P. Treatment of hepatocellular carcinoma: a systematic review. Liver Cancer. 2012;1(3–4):144–58.
 6.
Villanueva A, Toffanin S, Llovet JM. Linking molecular classification of hepatocellular carcinoma and personalized medicine: preliminary steps. Curr Opin Oncol. 2008;20(4):444.
 7.
Visser H, le Cessie S, Vos K, Breedveld FC, Hazes JM. How to diagnose rheumatoid arthritis early: a prediction model for persistent (erosive) arthritis. Arthritis Rheum. 2002;46(2):357–65.
 8.
Spitz MR, Etzel CJ, Dong Q, Amos CI, Wei Q, Wu X, Hong WK. An expanded risk prediction model for lung cancer. Cancer Prev Res. 2008;1(4):250–4.
 9.
Huang CL, Liao HC, Chen MC. Prediction model building and feature selection with support vector machines in breast cancer diagnosis. Expert Syst Appl. 2008;34(1):578–87.
 10.
Gillette MA, Carr SA. Quantitative analysis of peptides and proteins in biomedicine by targeted mass spectrometry. Nat Methods. 2013;10(1):28–34.
 11.
Whiteaker JR, Lin C, Kennedy J, Hou L, Trute M, Sokal I, Yan P, Schoenherr RM, Zhao L, Voytovich UJ. A targeted proteomicsbased pipeline for verification of biomarkers in plasma. Nat Biotechnol. 2011;29(7):625–34.
 12.
Kennedy JJ, Abbatiello SE, Kim K, Yan P, Whiteaker JR, Lin C, Kim JS, Zhang Y, Wang X, Ivey RG. Demonstrating the feasibility of largescale development of standardized assays to quantify human proteins. Nat Methods. 2014;11(2):149–55.
 13.
Kline RB. Principles and practice of structural equation modeling. New York: Guilford publications; 2015. p. 724.
 14.
Hwang H, Takane Y. Generalized structured component analysis. Psychometrika. 2004;69(1):81–99.
 15.
Lee S, Choi S, Kim YJ, Kim BJ, Hwang H, Park T. Pathwaybased approach using hierarchical components of collapsed rare variants. Bioinformatics. 2016;32(17):i586–94.
 16.
Llovet JM, Ricci S, Mazzaferro V, Hilgard P, Gane E, Blanc JF, de Oliveira AC, Santoro A, Raoul JL, Forner A. Sorafenib in advanced hepatocellular carcinoma. N Engl J Med. 2008;359(4):378–90.
 17.
Cheng AL, Kang YK, Chen Z, Tsao CJ, Qin S, Kim JS, Luo R, Feng J, Ye S, Yang TS. Efficacy and safety of sorafenib in patients in the AsiaPacific region with advanced hepatocellular carcinoma: a phase III randomised, doubleblind, placebocontrolled trial. Lancet Oncol. 2009;10(1):25–34.
 18.
Society AC: Cancer facts & figures 2017. American Cancer Society journal, CA: A Cancer Journal for Clinicians. 2017;1718.
 19.
Kim H, Yu SJ, Yeo I, Cho YY, Lee DH, Cho Y, Cho EJ, Lee JH, Kim YJ, Lee S. Prediction of Response to Sorafenib in Hepatocellular Carcinoma: A Marker Panel by Multiple Reaction MonitoringMass Spectrometry. Mol Cell Proteomics. 2017;mcp. M116:066704.
 20.
Lencioni R, Llovet JM. Modified RECIST (mRECIST) assessment for hepatocellular carcinoma. In: Seminars in liver disease: © Thieme Medical Publishers. Stuttgart, Germany; 2010. p. 052–60.
 21.
Chambers AG, Percy AJ, Simon R, Borchers CH. MRM for the verification of cancer biomarker proteins: recent applications to human plasma and serum. Expert Rev Proteomics. 2014;11(2):137–48.
 22.
MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, Kern R, Tabb DL, Liebler DC, MacCoss MJ. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010;26(7):966–8.
 23.
Hwang H. Regularized generalized structured component analysis. Psychometrika. 2009;74(3):517–30.
 24.
Green PJ. Iteratively reweighted least squares for maximum likelihood estimation, and some robust and resistant alternatives. J R Stat Soc Ser B Methodol. 1984:149–92.
 25.
McCullagh P. Generalized linear models. Eur J Oper Res. 1984;16(3):285–92.
 26.
Jang JY, Park T, Lee S, Kim Y, Lee SY, Kim SW, Kim SC, Song KB, Yamamoto M, Hatori T. Proposed nomogram predicting the individual risk of malignancy in the patients with branch duct type Intraductal papillary mucinous neoplasms of the pancreas. Ann Surg. 2016;266(6):10628.
 27.
Team RC. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2014.
 28.
Gray J, Chattopadhyay D, Beale GS, Patman GL, Miele L, King BP, Stewart S, Hudson M, Day CP, Manas DM. A proteomic strategy to identify novel serum biomarkers for liver cirrhosis and hepatocellular cancer in individuals with fatty liver disease. BMC Cancer. 2009;9(1):271.
 29.
Braconi C, Meng F, Swenson E, Khrapenko L, Huang N, Patel T. Candidate therapeutic agents for hepatocellular cancer can be identified from phenotypeassociated gene expression signatures. Cancer. 2009;115(16):3738–48.
 30.
Kong LQ, Zhu XD, Xu HX, Zhang JB, Lu L, Wang WQ, Zhang QB, Wu WZ, Wang L, Fan J. The clinical significance of the CD163+ and CD68+ macrophages in patients with hepatocellular carcinoma. PLoS One. 2013;8(3):e59771.
 31.
Tseng GC, Cheng C, Yu YP, Nelson J, Michalopoulos G, Luo JH. Investigating multicancer biomarkers and their crosspredictability in the expression profiles of multiple cancer types. Biomark Insights. 2009;4:57.
 32.
Hutcheson J, Bourgo RJ, Balaji U, Ertel A, Witkiewicz AK, Knudsen ES. Retinoblastoma protein potentiates the innate immune response in hepatocytes: significance for hepatocellular carcinoma. Hepatology. 2014;60(4):1231–40.
 33.
Kinoshita M, Miyata M. Underexpression of mRNA in human hepatocellular carcinoma focusing on eight loci. Hepatology. 2002;36(2):433–8.
 34.
Cai Z, Daescu O, Li M. Bioinformatics research and applications: 13th international symposium, ISBRA 2017. Honolulu: Springer; 2017. May 29–June 2, 2017, proceedings, vol. 10330
Acknowledgements
The abridged abstract of this work was previously published in the Proceedings of the 13th International Symposium on Bioinformatics Research and Applications (ISBRA 2017), Lecture Notes in Computer Science: Bioinformatics Research and Applications [34].
The authors thank Dr. Curt Balch, for extensive English editing.
Funding
This work was supported by the BioSynergy Research Project (2013M3A9C4078158) of the Ministry of Science, ICT and Future Planning through the National Research Foundation and by grants from the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HI15C2165, HI16C2037).
Publication of this article was funded by a grant from the Korea Health Technology R&D Project through the Korea Health Industry Development Institute (KHIDI), funded by the Ministry of Health & Welfare, Republic of Korea (grant number: HI16C2037).
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
About this supplement
This article has been published as part of BMC Bioinformatics Volume 19 Supplement 9, 2018: Selected articles from the 13th International Symposium on Bioinformatics Research and Applications (ISBRA 2017): bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume19supplement9.
Author information
Affiliations
Contributions
TP and SK conceived and designed the research; SK and TP developed the method and were involved in drafting of manuscript; SK and SC implemented the software; JHY recruited liver cancer patients; YK generated MRMMS data from the patients; SL designed the simulation studies. All the authors read and approved the final manuscript.
Corresponding author
Correspondence to Taesung Park.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Published
DOI
Keywords
 Biomarkers
 Componentbased structural equation modeling
 Drug response
 Liver cancer
 Multiple reaction monitoring mass spectrometry (MRMMS)
 Prediction model
 Sorafenib