- Research article
- Open Access
- Open Peer Review
This article has Open Peer Review reports available.
Possibility of multivariate function composed of plasma amino acid profiles as a novel screening index for non-small cell lung cancer: a case control study
© Maeda et al; licensee BioMed Central Ltd. 2010
Received: 7 January 2010
Accepted: 22 December 2010
Published: 22 December 2010
The amino-acid balance in cancer patients often differs from that in healthy individuals, because of metabolic changes. This study investigated the use of plasma amino-acid profiles as a novel marker for screening non-small-cell lung cancer (NSCLC) patients.
The amino-acid concentrations in venous blood samples from pre-treatment NSCLC patients (n = 141), and age-matched, gender-matched, and smoking status-matched controls (n = 423), were measured using liquid chromatography and mass spectrometry. The resultant study data set was subjected to multiple logistic regression analysis to identify amino acids related with NSCLC and construct the criteria for discriminating NSCLC patients from controls. A test data set derived from 162 patients and 3,917 controls was used to validate the stability of the constructed criteria.
The plasma amino-acid profiles significantly differed between the NSCLC patients and the controls. The obtained model (including alanine, valine, isoleucine, histidine, tryptophan and ornithine concentrations) performed well, with an area under the curve of the receiver-operator characteristic curve (ROC_AUC) of >0.8, and allowed NSCLC patients and controls to be discriminated regardless of disease stage or histological type.
This study shows that plasma amino acid profiling will be a potential screening tool for NSCLC.
Recently, computer-aided systems for data mining, for example by multivariate analysis, are now readily available and have shown promising results when applied to metabolic profiling for diagnostic purposes [1, 2]. Currently, several applications of metabolome analysis based on machine learning for human cancer diagnosis using peripheral blood or urine were demonstrated [3–10].
Among metabolites, the amino-acid balance in patients with various diseases often differs from that maintained in healthy individuals, as a result of metabolic changes. Amino acids are considered to be central compounds within metabolic networks. The blood serves as the medium linking the metabolic processes in the different organs of the human body. Human amino-acid metabolism in the blood has been monitored clinically for >30 years. Fischer's ratio, which is defined as the balance between branched-chain amino acids (BCAAs) and aromatic amino acids, has been used as an indicator of both the progression of liver fibrosis and the effectiveness of drug treatment . Specific abnormalities in amino-acid concentrations, as assessed using multivariate analysis, have also been reported in animal models of diabetes, in human liver fibrosis and in other pathologies [12–14].
The metabolism in cancer cells is known to be significantly altered compared with that in normal cells, and these changes are also reflected in the plasma amino-acid profiles of patients with various types of cancer. For example, a significant reduction in gluconeogenic amino acids (GAAs) and a significant increase in free tryptophan have been reported in lung cancer patients . Kubota et al. used plasma amino-acid profiles to discriminate between patients with breast cancer, gastrointestinal cancer, and head and neck cancers, and healthy controls . Therefore, detecting metabolic changes from amino-acid profiles could potentially be useful in cancer diagnosis.
Post-genomic technologies also offer possibilities for exploiting amino-acid profiling. Recently, novel methods for analyzing amino acids have been established using high-performance liquid chromatography (HPLC)-electrospray ionization (ESI)-mass spectrometry (MS) [17–19]. This will help to make amino-acid measurements easier and reduce both the time and the cost of analysis.
Therefore, one potentially useful metabolomics tool is the "AminoIndex", which could be a simple and versatile method for monitoring various pathological conditions . Here we investigated the possibility of "AminoIndex" as a novel diagnostic method for the screening of non-small-cell lung cancer (NSCLC).
All of the patients in the study had been diagnosed histologically with NSCLC at the Osaka Medical Centre for Cancer and Cardiovascular Diseases, Japan, between January 2006 and October 2008. While hospitalized, their informed consent for inclusion was obtained. Data from the first 141 patients enrolled between January 2006 and September 2007 were used as the study data set. A further 4,340 subjects without apparent cancers, who were undergoing comprehensive medical examinations at the Mitsui Memorial Hospital, Japan, in 2008, were recruited as control subjects. Of these, 423 were age-matched, gender-matched, and smoking status-matched with the patients in the study data set group. Data from the remaining patients and control subjects were used as the test data set. Data from an additional 15 SCLC patients, who were hospitalized at the Osaka Medical Centre for Cancer and Cardiovascular Diseases, Japan, between January 2006 and October 2008, were also used. Blood samples were collected from the controls and the NSCLC patients before any medical treatment. The study was conducted in accordance with the Declaration of Helsinki, and the protocol was approved by the ethics committees of the Osaka Medical Centre for Cancer and Cardiovascular Diseases and Mitsui Memorial Hospital. All subjects gave their informed consent for inclusion before they participated in the study.
Blood samples (5 ml) were collected from forearm veins, after overnight fasting, in tubes containing ethylenediaminetetraacetic acid (EDTA; Termo, Tokyo, Japan), and were immediately placed on ice. Plasma was prepared by centrifugation at 3,000 rpm and 4°C for 15 min, and then stored at -80°C until analysis. After plasma collection, all samples were stored and processed at the Life Science Institute of Ajinomoto Co., Inc. (Kawasaki, Japan). To reduce any bias introduced prior to analysis, samples were analyzed in random order. The plasma samples were deproteinized using acetonitrile at a final concentration of 80% before measurement. The amino-acid concentrations in the plasma were measured by HPLC-ESI-MS, followed by precolumn derivatization [17–19]. The analytical methods were described in detail previously . The concentrations of amino acids in the plasma were expressed as μM.
Statistical analysis of plasma amino-acid profile
The mean amino-acid concentrations ± standard deviations (SDs) were calculated. Differences between the plasma amino-acid concentrations in NSCLC patients and controls were assessed using the Mann-Whitney U-test and receiver-operator characteristic (ROC) curve. The area under the curve (AUC) for each ROC curve (the ROC_AUC) was calculated for each amino acid.
where z i,j was transformed concentration of of the i-th sample of the j-th amino acid, x i,j was concentration of the i-th sample of the j-th amino acid, n was sample size, and was the average concentration of j-th amino acid.
Machine learning and validation
First, an unconditional multiple logistic regression analysis with variable selection was used to construct a criterion for distinguishing NSCLC patients from controls using the study data set with the raw plasma concentrations of 21 amino acids as explanatory variables. The candidate variables of most appropriate logistic regression model, which had the minimum Akaike's information criterion (AIC) value, were selected from among all of the possible combinations in which the number of variables was below seven. A leave-one-out cross-validation (LOOCV) was performed to correct potential over-optimization for all models in parallel. Briefly, one sample was omitted from the study data set, and the logistic regression model was calculated for the remaining samples, to estimate coefficients for each amino acid. The logistic regression function values for the left-out sample were calculated based on the model. This process was repeated until every sample in the study data set had been left out once, and the function values generated were then used for AIC calculation. Finally, a case-control study was utilized for our study, and so a conditional logistic regression analysis, conditioned on the matching factors (i.e., gender, age, and smoking status), was performed in order to evaluate the association between the combination of amino acids obtained above and NSCLC. The discriminant score, which was defined as a logit of the conditional logistic regression function value, was constructed as a criterion. The degree of discriminancy of this score between NSCLC patients and controls was evaluated through the ROC curve. A distinct test data set, which had not been used in the model generation, was also used to confirm the stability of the obtained model, and to calculate the ROC_AUC values for the discriminant scores.
To assess the effects of cancer stage and histological type, both the study data set and the test data set was stratified according to the analysis parameters. To assess the effects of cancer stage and histological type on the discriminant scores of NSCLC patients, a subgroup analysis was performed using the ROC curve, in each data set. A two-sided P value of less than 0.05 was considered to indicate statistical significance.
All statistical analyses were performed using MATLAB (The Mathworks, Natick, MA), LogXact (Cytel, Cambridge, MA), and GraphPad Prism (GraphPad Software, La Jolla, CA).
Characteristics of patients and control subjects
Characteristics of study participants
Study data set
Test data set
(Male, Female, Unknown)
Squamous cell carcinoma
Changes in amino-acid concentrations in NSCLC patients
Plasma amino-acid concentration
Plasma concentration, μM
Patients (n = 141)
Controls (n = 423)
PCA of plasma amino-acid profile of study data set
Classifier for discriminating NSCLC patients
Furthermore, to verify the robustness of the resulting model, a ROC curve was generated using the split test data set, which had not been used to construct the model. A ROC_AUC of the ROC for the discriminant score was 0.812 in the test data set (Figure 1), again demonstrating that the obtained model performed well.
Subgroup analysis of the discriminant scores
From the point of view of cancer screening, attention might be paid to whether or not the obtained model also provides sufficient discriminating power to extract effectively patients with early-stage cancer and for all histological types. Thus, to investigate the consistency of the results based on the discriminant scores among different subpopulations defined by cancer stage and histological type, a subgroup analysis was performed using both study data set and the test data set. The discriminant scores of the SCLC patients were also calculated to verify whether the obtained model could discriminate them from the controls.
Lung cancer has been the leading cause of cancer death since 1998 and >60,000 patients have died since 2005 in Japan. The 5-year survival rate for patients undergoing surgery is only 61%, and an accurate screening method for lung cancer would be an important advance . In Japan, chest X-rays and sputum cytology are used for screening lung cancer. Although chest X-rays are useful for detecting peripheral lung cancer, two-thirds of patients diagnosed in this way have associated metastases, and this method is not sufficient to detect the early stages of the disease . In addition, highly skilled staffs are required to achieve sufficient accuracy. Sputum cytology might be useful for detecting upper respiratory-tract carcinoma, but this method has been reported to be inadequate for detecting peripheral lung cancer and lung cancer in asymptomatic non-smokers . Recently, low-dose helical computed tomography (CT) was reported to be capable of detecting small, early lung cancers in high-risk populations; however, it is not known whether using this method would affect the mortality rate due to lung cancer or whether it would be cost-effective .
In comparison to those methods, the "AminoIndex" would be easier to use, as it involves a relatively simple plasma assay, imposes a lower physical burden on patients and does not require advanced technical skills to perform . The current study demonstrated that plasma amino acid profiles were associated with NSCLC. The ROC_AUCs were 0.817 for the study data set under the conditional logistic regression analysis conditioned on the matching factors (Figure 1). Okamoto et al. recently reported that plasma amino-acid profiles might be used to screen colorectal and breast cancer . Despite the smaller sample size, they reported ROC_AUCs of 0.860 (with study data) and 0.910 (with test data) for colorectal cancer patients, and 0.906 (with study data) and 0.865 (with test data) for breast cancer patients . Our current study achieved similar discrimination power using data set with a larger sample size under controlling for potential confounders, thereby demonstrating the robustness of the model.
Many reports have shown that the metabolism, including that of amino acids, is notably altered in cancer cells [4, 24–26], and that the plasma amino-acid profiles are also changed [15, 16, 27–30]. Cascino et al. described significant increases in levels of Trp, Glu and Orn in lung cancer patients . Proenza and colleagues also reported an increased level of Orn in patients with lung cancer . Naini et al. reported reduced levels of plasma Arg in lung cancer patients .
Changes in the amino-acid balance and an increase in gluconeogenesis have been well documented, especially in cachexic patients with advanced cancer [32, 33]. In the current study, the obtained model identified patients at all stages of lung cancer and without cachexia equally well, suggesting that the method did not rely on detecting metabolic abnormalities associated with malnutrition, which might be present in advanced cancer patients (Figure 2A, B). Hirayama et al. reported no significant correlation between the levels of metabolites, including several amino acids, and the patients' tumour stage . And it was also reported that amino acids were frequently identified compounds among whole metabolites in blood in relation to cancer [3, 8]. The metabolism of specific amino acids is known to be associated with specific organs, such as muscle, liver or kidney, changes in the levels of amino acids are affected by their metabolism in, and excretion from, multiple organs of the body. Although it remains unclear how the metabolic changes occurring in tumour cells affect the systemic, plasma amino-acid profile, these results show that the metabolic changes caused by cancer development are at least partially responsible for the changes in plasma amino-acid profile seen even in lung cancer patients with early stage cancer. So, profiling the plasma free amino acids is similar to monitoring metabolic networks in multiple organs and it might better allow us to detect particular conditions in specific organs.
Since this study was designed as a case-control study, the obtained model could not be directly applied to further observation or prediction even though the robustness of the model was preliminarily demonstrated. Therefore model construction and validation using cohort with larger samples will be necessary to clarify its utility. Nonetheless, we believe that this screening technique could be a straightforward diagnostic method for the management of lung cancer.
The current study demonstrated that the plasma amino-acid profile of NSCLC patients differed from that of healthy subjects. And we showed that the multivariate classifier might be effective for discriminating lung cancer patients. Although further prospective validation will be necessary in the future, this method might be an effective and convenient screening tool for lung cancer patients.
We thank Dr. Jiro Okami, Dr. Kazuyuki Oda, Dr. Kiyonobu Ueno, Dr. Kazumi Nishino, Dr. Yuuki Akazawa and Dr. Junji Uchida for collecting the blood samples and background data on the patients. We thank Dr. Hiroshi Miyano, Mr. Kazutaka Shimbo, Mr. Hiroo Yoshida, Ms. Michiko Amao and Ms. Mina Nakamura for amino-acid analyses. We thank Dr. Katsuhisa Horimoto, and Dr. Mitsuo Takahashi for help with the statistical analysis. We also thank Ms. Tomoko Kasakura for help with data collection.
- Chadeau-Hyam M, Ebbels TM, Brown IJ, Chan Q, Stamler J, Huang CC, Daviglus ML, Ueshima H, Zhao L, Holmes E, et al: Metabolic profiling and the metabolome-wide association study: significance level for biomarker identification. J Proteome Res. 2010, 9 (9): 4620-4627. 10.1021/pr1003449.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee K, Hwang D, Yokoyama T, Stephanopoulos G, Stephanopoulos GN, Yarmush ML: Identification of optimal classification functions for biological sample and state discrimination from metabolic profiling data. Bioinformatics. 2004, 20 (6): 959-969. 10.1093/bioinformatics/bth015.View ArticlePubMedGoogle Scholar
- Urayama S, Zou W, Brooks K, Tolstikov V: Comprehensive mass spectrometry based metabolic profiling of blood plasma reveals potent discriminatory classifiers of pancreatic cancer. Rapid Commun Mass Spectrom. 2010, 24 (5): 613-620. 10.1002/rcm.4420.View ArticlePubMedGoogle Scholar
- Rocha CM, Barros AS, Gil AM, Goodfellow BJ, Humpfer E, Spraul M, Carreira IM, Melo JB, Bernardo J, Gomes A, et al: Metabolic profiling of human lung cancer tissue by 1 H high resolution magic angle spinning (HRMAS) NMR spectroscopy. J Proteome Res. 2010, 9 (1): 319-332. 10.1021/pr9006574.View ArticlePubMedGoogle Scholar
- Pasikanti KK, Esuvaranathan K, Ho PC, Mahendran R, Kamaraj R, Wu QH, Chiong E, Chan EC: Noninvasive urinary metabonomic diagnosis of human bladder cancer. J Proteome Res. 2010, 9 (6): 2988-2995. 10.1021/pr901173v.View ArticlePubMedGoogle Scholar
- Nolen B, Velikokhatnaya L, Marrangoni A, De Geest K, Lomakin A, Bast RC, Lokshin A: Serum biomarker panels for the discrimination of benign from malignant cases in patients with an adnexal mass. Gynecol Oncol. 2010, 117 (3): 440-445. 10.1016/j.ygyno.2010.02.005.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim Y, Koo I, Jung BH, Chung BC, Lee D: Multivariate classification of urine metabolome profiles for breast cancer diagnosis. BMC Bioinformatics. 2010, 11 (Suppl 2): S4-10.1186/1471-2105-11-S2-S4.View ArticlePubMedPubMed CentralGoogle Scholar
- Qiu Y, Cai G, Su M, Chen T, Zheng X, Xu Y, Ni Y, Zhao A, Xu LX, Cai S, et al: Serum metabolite profiling of human colorectal cancer using GC-TOFMS and UPLC-QTOFMS. J Proteome Res. 2009, 8 (10): 4844-4850. 10.1021/pr9004162.View ArticlePubMedGoogle Scholar
- Couzin J: Biomarkers. Metabolite in urine may point to high-risk prostate cancer. Science. 2009, 323 (5916): 865-10.1126/science.323.5916.865a.View ArticlePubMedGoogle Scholar
- Abate-Shen C, Shen MM: Diagnostics: The prostate-cancer metabolome. Nature. 2009, 457 (7231): 799-800. 10.1038/457799a.View ArticlePubMedGoogle Scholar
- Fischer JE, Funovics JM, Aguirre A, James JH, Keane JM, Wesdorp RI, Yoshimura N, Westman T: The role of plasma amino acids in hepatic encephalopathy. Surgery. 1975, 78 (3): 276-290.PubMedGoogle Scholar
- Kimura T, Noguchi Y, Shikata N, Takahashi M: Plasma amino acid analysis for diagnosis and amino acid-based metabolic networks. Curr Opin Clin Nutr Metab Care. 2009, 12 (1): 49-53. 10.1097/MCO.0b013e3283169242.View ArticlePubMedGoogle Scholar
- Noguchi Y, Zhang QW, Sugimoto T, Furuhata Y, Sakai R, Mori M, Takahashi M, Kimura T: Network analysis of plasma and tissue amino acids and the generation of an amino index for potential diagnostic use. Am J Clin Nutr. 2006, 83 (2): 513S-519S.PubMedGoogle Scholar
- Zhang Q, Takahashi M, Noguchi Y, Sugimoto T, Kimura T, Okumura A, Ishikawa T, Kakumu S: Plasma amino acid profiles applied for diagnosis of advanced liver fibrosis in patients with chronic hepatitis C infection. Hepatol Res. 2006, 34 (3): 170-177. 10.1016/j.hepres.2005.12.006.View ArticlePubMedGoogle Scholar
- Cascino A, Muscaritoli M, Cangiano C, Conversano L, Laviano A, Ariemma S, Meguid MM, Rossi Fanelli F: Plasma amino acid imbalance in patients with lung and breast cancer. Anticancer Res. 1995, 15 (2): 507-510.PubMedGoogle Scholar
- Kubota A, Meguid MM, Hitch DC: Amino acid profiles correlate diagnostically with organ site in three kinds of malignant tumors. Cancer. 1992, 69 (9): 2343-2348. 10.1002/1097-0142(19920501)69:9<2343::AID-CNCR2820690924>3.0.CO;2-S.View ArticlePubMedGoogle Scholar
- Shimbo K, Kubo S, Harada Y, Oonuki T, Yokokura T, Yoshida H, Amao M, Nakamura M, Kageyama N, Yamazaki J, et al: Automated precolumn derivatization system for analyzing physiological amino acids by liquid chromatography/mass spectrometry. Biomed Chromatogr. 2009, 24 (7): 683-691. 10.1002/bmc.1346.View ArticleGoogle Scholar
- Shimbo K, Oonuki T, Yahashi A, Hirayama K, Miyano H: Precolumn derivatization reagents for high-speed analysis of amines and amino acids in biological fluid using liquid chromatography/electrospray ionization tandem mass spectrometry. Rapid Commun Mass Spectrom. 2009, 23 (10): 1483-1492. 10.1002/rcm.4026.View ArticlePubMedGoogle Scholar
- Shimbo K, Yahashi A, Hirayama K, Nakazawa M, Miyano H: Multifunctional and highly sensitive precolumn reagents for amino acids in liquid chromatography/tandem mass spectrometry. Anal Chem. 2009, 81 (13): 5172-5179. 10.1021/ac900470w.View ArticlePubMedGoogle Scholar
- Asamura H, Goya T, Koshiishi Y, Sohara Y, Eguchi K, Mori M, Nakanishi Y, Tsuchiya R, Shimokata K, Inoue H, et al: A Japanese Lung Cancer Registry study: prognosis of 13,010 resected lung cancers. J Thorac Oncol. 2008, 3 (1): 46-52. 10.1097/JTO.0b013e31815e8577.View ArticlePubMedGoogle Scholar
- Sumitani M, Takifuji N, Nanjyo S, Imahashi Y, Kiyota H, Takeda K, Yamamoto R, Tada H: Clinical relevance of sputum cytology and chest X-ray in patients with suspected lung tumors. Intern Med. 2008, 47 (13): 1199-1205. 10.2169/internalmedicine.47.0777.View ArticlePubMedGoogle Scholar
- Jett JR: Limitations of screening for lung cancer with low-dose spiral computed tomography. Clin Cancer Res. 2005, 11 (13 Pt 2): 4988s-4992s. 10.1158/1078-0432.CCR-05-9000.View ArticlePubMedGoogle Scholar
- Okamoto N, Miyagi Y, Chiba A, Akaike M, Shiozawa M, Imaizumi A, Yamamoto H, Ando T, Yamakado M, Tochikubo O: Diagnostic modeling with differences in plasma amino acid profiles between non-cachectic colorectal/breast cancer patients and healthy individuals. Metabolic changes in patients with cancer lead to alterations in their amino-acid balances Thus, amino-acid profiles may be useful as biomarkers of cancers This study was conducted to analyze amino-acid profiles in plasma by multivar. 2009, 1 (1): 1-8.Google Scholar
- Hirayama A, Kami K, Sugimoto M, Sugawara M, Toki N, Onozuka H, Kinoshita T, Saito N, Ochiai A, Tomita M, et al: Quantitative metabolome profiling of colon and stomach cancer microenvironment by capillary electrophoresis time-of-flight mass spectrometry. Cancer Res. 2009, 69 (11): 4918-4925. 10.1158/0008-5472.CAN-08-4806.View ArticlePubMedGoogle Scholar
- Rodriguez PC, Quiceno DG, Zabaleta J, Ortiz B, Zea AH, Piazuelo MB, Delgado A, Correa P, Brayer J, Sotomayor EM, et al: Arginase I production in the tumor microenvironment by mature myeloid cells inhibits T-cell receptor expression and antigen-specific T-cell responses. Cancer Res. 2004, 64 (16): 5839-5849. 10.1158/0008-5472.CAN-04-0465.View ArticlePubMedGoogle Scholar
- Yamaguchi K, Saito H, Oro S, Tatebe S, Ikeguchi M, Tsujitani S: Expression of inducible nitric oxide synthase is significantly correlated with expression of vascular endothelial growth factor and dendritic cell infiltration in patients with advanced gastric carcinoma. Oncology. 2005, 68 (4-6): 471-478. 10.1159/000086990.View ArticlePubMedGoogle Scholar
- Laviano A, Cascino A, Muscaritoli M, Fanfarillo F, Rossi Fanelli F: Tumor-induced changes in host metabolism: a possible role for free tryptophan as a marker of neoplastic disease. Adv Exp Med Biol. 2003, 527: 363-366.View ArticlePubMedGoogle Scholar
- Norton JA, Gorschboth CM, Wesley RA, Burt ME, Brennan MF: Fasting plasma amino acid levels in cancer patients. Cancer. 1985, 56 (5): 1181-1186. 10.1002/1097-0142(19850901)56:5<1181::AID-CNCR2820560535>3.0.CO;2-8.View ArticlePubMedGoogle Scholar
- Proenza AM, Oliver J, Palou A, Roca P: Breast and lung cancer are associated with a decrease in blood cell amino acid content. J Nutr Biochem. 2003, 14 (3): 133-138. 10.1016/S0955-2863(02)00225-5.View ArticlePubMedGoogle Scholar
- Vissers YL, Dejong CH, Luiking YC, Fearon KC, von Meyenfeldt MF, Deutz NE: Plasma arginine concentrations are reduced in cancer patients: evidence for arginine deficiency?. Am J Clin Nutr. 2005, 81 (5): 1142-1146.PubMedGoogle Scholar
- Naini AB, Dickerson JW, Brown MM: Preoperative and postoperative levels of plasma protein and amino acid in esophageal and lung cancer patients. Cancer. 1988, 62 (2): 355-360. 10.1002/1097-0142(19880715)62:2<355::AID-CNCR2820620221>3.0.CO;2-E.View ArticlePubMedGoogle Scholar
- Heber D, Byerly LO, Chlebowski RT: Metabolic abnormalities in the cancer patient. Cancer. 1985, 55 (1 Suppl): 225-229. 10.1002/1097-0142(19850101)55:1+<225::AID-CNCR2820551304>3.0.CO;2-7.View ArticlePubMedGoogle Scholar
- Kern KA, Norton JA: Cancer cachexia. JPEN J Parenter Enteral Nutr. 1988, 12 (3): 286-298. 10.1177/0148607188012003286.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/10/690/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<url>http://creativecommons.org/licenses/by/2.0</url>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.