Abstract
Background:
Recently, a novel programmed cell death mechanism, Cuproptosis, has been discovered and found to play an important role in the development and progression of diverse tumors. In the present study, we comprehensively investigated the core gene of this mechanism, GLS, in breast cancer.
Materials and methods:
Bulk RNA sequencing data were curated from the TCGA repository to investigate the aberrant expression of GLS over diverse cancer types. Then, we examined its efficacy as a diagnostic biomarker in breast cancer by Area Under Curve (AUC) of the Receiver Operative Characteristic (ROC) curve. Furthermore, by applying siRNA technique, we knocked down the GLS expression level in cancerous cell lines, measuring the corresponding effects on cell proliferation and metastasis. Afterward, we explored the potential implications of GLS expression in the tumor immune microenvironment quantitatively by using several R packages and algorithms, including ESTIMATE, CIBERSORT, etc.
Results:
Pan-cancer analysis suggested that GLS was aberrantly over-expressed in many cancer types, with breast cancer being typical. More in-depth analyses revealed the expression of GLS exerted a high ROC-AUC value in breast cancer diagnosis. Through the knock-down of GLS expression, it was found that GLS expression was strongly relevant to the growth and metastasis of tumor. Furthermore, it was also found to be correlated with the immune tumor microenvironment.
Conclusion:
We highlighted that GLS expression might be applicable as a diagnostic biomarker in breast cancer and possess significant implications in the growth and metastasis of tumor and the immune tumor microenvironment, sharing new insights into ontological and personalized medicine.
Introduction
Breast cancer is the most common malignant tumor in women and one of the main causes of disease burden worldwide (1). Despite the significant advances in diagnosis and prognosis, treating patients with advanced stages remains challenging (2–5). To improve the prognosis, early detection, and timely treatment are of great significance. Therefore, revealing new and highly sensitive molecular biomarkers remains critical. Cuproptosis, a newly discovered cell death mechanism, serves an important role in the modification of specific mitochondrial metabolic enzymes, through which it triggers a series of intracellular signaling events and ultimately leads to cellular caseation (6). Earlier studies have also shown that the aberration of copper levels in the serum and tissues of patients with different solid tumors indicates poor clinical outcomes in a fair number of cases (7–10). Indeed, copper as a basic metal ion in most aerobic organisms serves as a structural and catalytic cofactor for many essential biological enzymes (11). Furthermore, in some studies, certain copper chelators were applied in anticancer therapies, such as tetrathiomolybylate, which has shown a clear survival benefit for high-risk triple-negative breast cancer (12, 13). Therefore, it is of ration that utilizing Cuproptosis-related genes to assist in the early screening of breast cancer. GLS, one of the 10 Cuproptosis-related genes, is an encoding gene of glutaminase which is the enzyme catalyzing the hydrolysis of L-β-glutamine into L-glutamate and ammonia. As glutamine metabolism plays an important role in various tumors, GLS has also become a promising therapeutic target in cancer treatment (14, 15). However, back to date, the utilization of GLS in breast cancer diagnosis remains in need. From this end, it is of interest to explore the role of GLS from this aspect. In the present study, by integrating bioinformatics analytics and immunohistochemical (IHC) staining of real-world patient samples from the HPA database, as well as the siRNA technique, we proposed the idea that GLS may serve as a valuable diagnostic biomarker for breast cancer. Figure 1 demonstrates the general design of the present study.
Figure 1
Materials and methods
Data curation
The sequencing data of 10,534 TCGA pan-cancer samples and 15,776 GTEx normal samples were downloaded from the UCSC Xena database. In particular, for breast cancer, raw transcriptomic count data and the normalized mRNA profile fragments per kilobase of exon per million reads (FPKM) of 115 TNBC patients and 113 normal tissues were retrieved.
Logistics modeling and diagnostic nomogram construction
The logistics algorithm was used to model the data so that the predictive ability of GLS expression for breast cancer could be further improved. Subsequently, we established a nomogram for predicting the diagnostic probability of breast cancer using the expression of GLS.
Estimation of the tumor microenvironment condition
The tumor microenvironment condition was assessed quantitatively by calculating the levels of stromal and immune cell infiltration using the expression profiles obtained from the TCGA dataset. This was done by the R package “ESTIMATE” in which the stromal score, immune score, and ESTIMATE score were calculated (16).
Screening of immune cell infiltration
For the analysis of immune cell infiltration, we used the CIBERSORT algorithm which is a strong and robust immunoinformatic method to evaluate the immune cell population based on bulk RNA transcriptome. This was achieved by the R package “immunedeconv” (17).
Single-cell analysis
The single-cell transcriptomic profile of the GSE114727 dataset (18) was recorded in the Gene Expression Omnibus (GEO) database. The visualization was done by the online toolkit TISCH (http://tisch.comp-genomics.org/) (19).
Cell culturing and transfection
The breast cancer cell line, MCF7, was cultured in Huazhong University (April 2023). The medium was specially designed (i.e., DMEM + 10% FBS + 20ng/ml EGF + 0.5ug/mL Hydrocortisone + 10ug/mL Insulin + 1% NEAA + 1% P/S) and regularly changed when the medium colour turned yellowish. The culturing flasks were saved consistently at 37°C in a humidified incubator containing 5% CO2.
The vectors, si-NC and si-GLS were transfected into MCF7 cell line by Liposome 2000 transfection agent (Invitgen, USA) according to the manufacturer’s protocol. Sequences of si-NC was: UUCUCCGAACGUGUCACGUTT ACGUGACACGUUCGGAGAATT, and sequences of si-GLS was: GAUGGUGUCAUGCUAGACAAATT UUUGUCUAGCAUGACACCAUCTT.
CCK8 assay
The siRNA-transfected MCF7 cells were fully digested, inoculated into 96-well culture plates at 2 x 103 cells per well and incubated at 37°C. The absorptivity was measured every day using the CCK-8 kit, for a total of seven days of incubation. The above procedures were carried out at 450 nm on the advice of the reagent supplier and four secondary wells were retained for the experimental sessions.
Western blotting
We treated the MCF7 cell line with cell lysis buffer at 100°C for a total of 10 min with cells that had been previously rinsed with cold PBS to extract total protein. Prepared proteins were added to a 10% SDS/PAGE grid system and equal amounts of proteins were electrolysed at constant pressure (i.e., 200V) and then transferred to a PVDF membrane. After sealing the membrane with 5% skimmed milk powder for 1h, we incubated it overnight at 4°C with primary antibodies, followed by 1.5h incubation with secondary antibodies. Then, the PVDF membrane strips were washed and visualized, the relative protein abundance was assessed using Image Lab analysis software and ImageJ.
Statistical analysis
R Studio software (v 4.1.0) was used for bioinformatics analytics and GraphPad Prism 5.0 (GraphPad Software, Inc., San Diego, CA) was used for semi-quantification of the real-time PCR and western blot results. The software was also used to generate the relevant charts. The Chi-square test or the Wilcoxon signed-rank test was conducted to calculate the P-values. The correlation analyses were done under the Spearman method. The area under the curve (AUC) and the receiver operating characteristic (ROC) curve were analyzed to determine the predictive accuracy. P-value < 0.05 was deemed statistically significant.
Results
Pan-cancer analysis of GLS expression across various cancer types
Pan-cancer analysis is a widely used and efficient approach to identify the diagnostic value of a specific gene in regard of human tumors (20, 21). Therefore, we first conducted a pan-cancer analysis based on the transcriptomic data from the TCGA database, the result of which revealed that the expression of the GLS gene was significantly down-regulated in diverse cancer types, including but not limited to GBM, KICH, KIRC, and KIRP, etc. (Figure 2A), particularly in the breast cancer (abbreviation: BRCA or BC, Figure 2B). As the TCGA database lacks normal controls, we further combined the transcriptomic data from the GTEx repository for comparison. As a result, the previous observations were confirmed (Figures 2C, D).
Figure 2
Evaluation of the reliability of GLS expression as a diagnostic biomarker in breast cancer
We used the TCGA dataset alone (Figure 3A) and combined it with the GTEx database data (Figure 3B) to plot the receiver operating curves (ROC curves), through which we found that the ROC-AUC values were high (in the case of solely used TCGA dataset, ROC-AUC value = 0.68, and in the case of TCGA-GTEx combination, ROC-AUC value = 0.75), suggesting that GLS may serve as an important diagnostic biomarker breast cancer. Then, we used the logistics algorithm to optimize a model aiming for better diagnostic prediction and constructed a GLS-based nomogram (Figure 3C), with a C-index = 0.680, hindering it possessed a very good clinical diagnostic ability, which once again supported our idea.
Figure 3
By knocking down the GLS expression via the siRNA techinique, we found that the down expression of GLS promoted the vitality of tumor cells (Figure 3D).
As we know, the emt pathway has a very important role in cancer deterioration. From the very general aspect, the up regulation of emt pathway is regarded as one of the potential bases of tumor metastasis. On the other hand, E-adherin and N-cadherin, the 2 key proteins in the emt pathway are thought as the most active players that anchor the tumor tissues to their place of origin. Namely, when the 2 proteins are down regulated, the emt pathway will be activated, and so more likely to have tumor metastasis occurred. Our observation showed that after intervening with si-GLS, the expression levels of both E-cadherin and N-cadherin were significantly decreased, indicating a promoting effect on the pathway, and so the metastasis of the tumor (Figures 3E, F).
The expression of GLS was tightly associated with multiple clinical characteristics of breast cancer
A genetic signature coined as PAM50 composed of 50 genes is commonly used clinically to divide breast cancer into several different molecular subtypes, and even now is still being a critical indicator in the treatment of breast cancer (22). Here we used GLS instead of this signature (Figure 4A) and found that the classification can be performed more efficiently, although the classification is somehow rougher (Luminal type A, Luminal type B, and Her2 type were indistinguishable, thus hereby we termed it as mixed type). The ROC analyses of the comparison between Basal type and mixed type (Figure 4B), between Basal type and Normal type (Figure 4C), and between Normal type and mixed type (Figure 4D) were also conducted. As a result, we found that such GLS-based distinguishments were of high accuracy with high ROC-AUC values, especially in the case of Basal type and mixed type, the ROC-AUC value was up to 0.854.
Figure 4
PR-positive cancer cells may require progesterone to grow. Consequently, these cells may stop growing or die when treated with substances that block the binding and action of progesterone. Therefore, patients with positive PR in the short term have better treatment prospects. In this study, we found that GLS expression was significantly reduced in PR-positive patients than in PR-negative patients (Figure 4E). This conjecture was then confirmed by the subsequent ROC analysis results (AUC-value = 0.681, Figure 4F).
Similar to PR status, ER-positive cancer cells may require estrogen to grow. These cells may stop growing or die when treated with substances that block the binding and action of estrogen. Therefore, patients who are ER-positive in the short term have better treatment prospects. We found that GLS expression was also significantly less in ER-positive patients than in PR-negative patients (Figure 4G). This conjecture was also confirmed by the subsequent ROC analysis results (AUC-value = 0.757, Figure 4H).
Taking them altogether, based on the PR/ER status, it seems that the low expression of GLS may exert adverse effects on the treatment in the short term. In addition, its expression level in breast cancer is lower than that in normal tissues, as such, finding ways to increase its expression level may be a new adjuvant treatment strategy.
Exploration of the immune implications of GLS in breast cancer
We first assessed the association of GLS with the tumor microenvironment by using the ESTIMATE algorithm and found that it was significantly associated with the stromal cells and immune cells (Figure 5A). This conclusion was further strengthened by the group comparison between high- and low-GLS expression (Figure 5B). Then, we performed ssGSEA analysis on a variety of immune cells and the expression of GLS. Subsequently, it was found that GLS extensively affected the infiltration of almost all immune cell types (Figure 5C). This conclusion was also further strengthened by the group comparison of the expression level of GLS (Figure 5D), as well as its visualization at single-cell transcriptome (Figure 5E), which further highlighted the important impact of GLS on the immune microenvironment of breast cancer.
Figure 5
Functional analysis
Through the ssGSEA algorithm, we investigated the functionalities of GLS in breast cancer. It was found that except for DNA repair (Figure 6H), its expression was positively associated with the rest pathways (Figures 6A-G, I).
Figure 6
In addition, we supplemented an analysis of drug sensitivity targeting the GLS gene based on the GDSC repository (Supplementary Figure S1), from that we found KIN001−102, Phenformin, WZ3105, Vinorelbine, and SB52334 were positively related to GLS expression. They might have certain value to be tested in the future.
Discussion
Breast cancer is a life-threatening disease that poses a significant challenge to the world and has been becoming a pivotal cause of global women’s cancer death (1). Although the survival outcomes have significantly improved over the past decades, patients with metastatic breast cancer still show poor prognoses (2–5). As such, early detection and genetic screening of breast cancer are of great interest to public health.
On the other hand, Cuproptosis, a novel cell death mechanism induced by copper ion imbalance, has been shown to play a significant role in the development and deterioration of various cancers, indicating that it could serve as a potential biological target for diagnosing or treating these illnesses (6–13). Meanwhile, GLS, the gene encoding glutaminase, is an important Cuproptosis-related gene. The dependence of a plethora of malignancies on glutaminase has prompted the proposition of glutaminase inhibitors as a viable therapeutic strategy for cancer (14, 15).
To date, although previous studies have shown that imbalanced copper metabolism can significantly impact the tumor microenvironment, adjusting the immunotherapy responsiveness, and several researchers even established GLS-involved signatures associated with human cancer cell proliferation and metastasis, there is still a lack of focus on the correlation between GLS and breast cancer. From this end, the present study explored the possibility of using GLS as a diagnostic biomarker for breast cancer and further clarified the underlying mechanisms by the bioinformatics method combined with clinical investigations on real-world samples.
We first conducted a pan-cancer analysis, reviewing its expressions in different cancers. As a result, its aberrant expressions were observed in diverse cancer types, especially breast cancer. In breast cancer, it was found that GLS was abnormally expressed which raised our interest in the possibility of utilizing it as a diagnostic biomarker in the future. In fact, the ROC-AUC values in the TCGA cohort alone and combination with GTEx data were surprisingly satisfying. Therefore, we further modeled a diagnostic nomogram aiming to assist physicians in making clinical decisions in their daily practice. More important, we surprisingly found that the expression of GLS could be used for a more efficient classification of molecular subtypes in breast cancer, although certain more exact subtypes were not applicable.
Given that nowadays the detailed interplays of the tumor immune microenvironment (TIME) have been unraveled drastically, its major component, namely the vast diversity of immune cell populations such as regulatory T cells, cytotoxic T cells, helper T cells, B cells, monocytes, etc. are deemed to play a critical role regarding various diseases, including cancers (23–26). Therefore, it was of huge interest to investigate the TIME in the present study. To the best of our knowledge, this is the very first time to explore the relationship between GLS expression and the TIME in breast cancer. We first utilized the ESTIMATE algorithm to examine the association of GLS expression with the TIME. Through quantitatively measuring the stromal and immune landscapes by the ssGSEA method, we found that the expression of GLS possessed statistically significant impacts in both regards, which was further supported by the comparison of high and low GLS expression groups. More importantly, at the single-cell level, we found its expression was enriched in various immune cells including CD8+ and CD4+ T cells, B cells, and more. These findings suggested that GLS might extensively involve in the shaping of the tumor microenvironment in breast cancer.
In conclusion, in the present study, we performed in-depth bioinformatics analytics of GLS expression in breast cancer, based on which we established a practically useful diagnostic monogram, which to a large extent, cross-supporting the previous research, although form a different aspect (27). Collectively speaking, the idea of using GLS as a biomarker for breast cancer is convincing. Moreover, we provided significant mechanistic insights into GLS’s impacts on the tumor microenvironment, focusing on the immune aspect.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Author contributions
Conceptualization: XH, DZ, QW, ZL, and GW. Data curation: XH, DZ, MW, SZ, LW, and YT. Formal analysis: XH and DZ. Visualization: XH, YL, SZ, QW and ZL. Writing – original draft: XH, DZ, MW, LW, SZ, YT, QW, ZL, YT, and YL. Writing - review and edition: YL and GW. All authors contributed to the article and approved the submitted version.
Acknowledgments
We want to express our deep gratitude to the public databases, including the TCGA and HPA repositories, and more, for providing open-accessible and high-quality research resources. We also sincerely thank the Fundamental Research Funds for the Central Universities of China for their financial support of our present study.
Conflict of interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in the present study.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2023.1220038/full#supplementary-material
Supplementary Figure S1Chemosensitivity test targeting GLS gene expression.
References
1
SungHFerlayJSiegelRLLaversanneMSoerjomataramIJemalAet al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA: Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660
2
LiLDunsGJDessieWCaoZJiXLuoX. Recent advances in peptide-based therapeutic strategies for breast cancer treatment. Front Pharmacol (2023) 14:1052301. doi: 10.3389/fphar.2023.1052301
3
DenJSistiA. Recent advances in breast cancer diagnosis, treatment, psychology, management, and reconstruction. Medicina (Kaunas Lithuania) (2023) 59(2):212. doi: 10.3390/medicina59020212
4
SwainSMShastryMHamiltonE. Targeting HER2-positive breast cancer: advances and future directions. Nat Rev Drug Discovery (2023) 22(2):101–26. doi: 10.1038/s41573-022-00579-0
5
GeJZuoWChenYShaoZYuK. The advance of adjuvant treatment for triple-negative breast cancer. Cancer Biol Med (2021) 19(2):187–201. doi: 10.20892/j.issn.2095-3941.2020.0752
6
TsvetkovPCoySPetrovaBDreishpoonMVermaAAbdusamadMet al. Copper induces cell death by targeting lipoylated TCA cycle proteins. Sci (New York NY) (2022) 375(6586):1254–61. doi: 10.1126/science.abf0529
7
KahlsonMADixonSJ. Copper-induced cell death. Sci (New York NY) (2022) 375(6586):1231–2. doi: 10.1126/science.abo3959
8
CobinePABradyDC. Cuproptosis: Cellular and molecular mechanisms underlying copper-induced cell death. Mol Cell (2022) 82(10):1786–7. doi: 10.1016/j.molcel.2022.05.001
9
JiangYHuoZQiXZuoTWuZ. Copper-induced tumor cell death mechanisms and antitumor theragnostic applications of copper complexes. Nanomed (London England) (2022) 17(5):303–24. doi: 10.2217/nnm-2021-0374
10
DuanWJHeRR. Cuproptosis: copper-induced regulated cell death. Sci China Life Sci (2022) 65(8):1680–2. doi: 10.1007/s11427-022-2106-6
11
ChenLMinJWangF. Copper homeostasis and cuproptosis in health and disease. Signal Transduct Target Ther (2022) 7(1):378. doi: 10.1038/s41392-022-01229-y
12
EmamiFAliomraniMTangestaninejadSKazemianHMoradiMRostamiM. Copper-curcumin-bipyridine dicarboxylate complexes as anticancer candidates. Chem Biodivers (2022) 19(10):e202200202. doi: 10.1002/cbdv.202200202
13
IshidaSMcCormickFSmith-McCuneKHanahanD. Enhancing tumor-specific uptake of the anticancer drug cisplatin with a copper chelator. Cancer Cell (2010) 17(6):574–83. doi: 10.1016/j.ccr.2010.04.011
14
ChenYTanLGaoJLinCWuFLiYet al. Targeting glutaminase 1 (GLS1) by small molecules for anticancer therapeutics. Eur J Medicinal Chem (2023) 252:115306. doi: 10.1016/j.ejmech.2023.115306
15
JinHWangSZaalEAWangCWuHBosmaAet al. A powerful drug combination strategy targeting glutamine addiction for the treatment of human liver cancer. eLife (2020) 9:e56749. doi: 10.7554/eLife.56749
16
YoshiharaKShahmoradgoliMMartínezEVegesnaRKimHTorres-GarciaWet al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nat Commun (2013) 4:2612. doi: 10.1038/ncomms3612
17
SturmGFinotelloFPetitprezFZhangJDBaumbachJFridmanWHet al. Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinf (Oxford England) (2019) 35(14):i436–45. doi: 10.1093/bioinformatics/btz363
18
AziziECarrAJPlitasGCornishAEKonopackiCPrabhakaranSet al. Single-cell map of diverse immune phenotypes in the breast tumor microenvironment. Cell (2018) 174(5):1293–1308.e36. doi: 10.1016/j.cell.2018.05.060
19
HanYWangYDongXSunDLiuZYueJet al. TISCH2: expanded datasets and new tools for single-cell transcriptome analyses of the tumor microenvironment. Nucleic Acids Res (2023) 51(D1):D1425–31. doi: 10.1093/nar/gkac959
20
XieJZhangJTianWZouYTangYZhengSet al. The pan-cancer multi-omics landscape of FOXO family relevant to clinical outcome and drug resistance. Int J Mol Sci (2022) 23(24):15647. doi: 10.3390/ijms232415647
21
PanSSunSLiuBHouY. Pan-cancer Landscape of the RUNX Protein Family Reveals their Potential as Carcinogenic Biomarkers and the Mechanisms Underlying their Action. J Trans Internal Med (2022) 10(2):156–74. doi: 10.2478/jtim-2022-0013
22
XieJZouYGaoTXieLTanDXieX. Therapeutic landscape of human epidermal growth factor receptor 2-positive breast cancer. Cancer Control J Moffitt Cancer Center (2022) 29:10732748221099230. doi: 10.1177/10732748221099230
23
AlokASeokKWesolowJ. A case of abdominal pain and diarrhea post immunotherapy: hypophysitis associated with immune checkpoint inhibitors. J Trans Internal Med (2022) 10(2):178–80. doi: 10.2478/jtim-2022-0030
24
TangYTianWXieJZouYWangZLiNet al. Prognosis and dissection of immunosuppressive microenvironment in breast cancer based on fatty acid metabolism-related signature. Front Immunol (2022) 13:843515. doi: 10.3389/fimmu.2022.843515
25
LuoYTianWLuXZhangCXieJDengXet al. Prognosis stratification in breast cancer and characterization of immunosuppressive microenvironment through a pyrimidine metabolism-related signature. Front Immunol (2022) 13:1056680. doi: 10.3389/fimmu.2022.1056680
26
ZouYYeFKongYHuXDengXXieJet al. The single-cell landscape of intratumoral heterogeneity and the immunosuppressive microenvironment in liver and brain metastases of breast cancer. Adv Sci (Weinheim Baden-Wurttemberg Germany) (2023) 10(5):e2203699. doi: 10.1002/advs.202203699
27
BudcziesJPfitznerBMGyörffyBWinzerKJRadkeCDietelMet al. Glutamate enrichment as new diagnostic opportunity in breast cancer. Int J Cancer (2015) 136(7):1619–28. doi: 10.1002/ijc.29152
Summary
Keywords
breast cancer, bioinformatics, cuproptosis, biomarker, EMT pathway
Citation
Zhang D, Wang M, Huang X, Wang L, Liu Y, Zhou S, Tang Y, Wang Q, Li Z and Wang G (2023) GLS as a diagnostic biomarker in breast cancer: in-silico, in-situ, and in-vitro insights. Front. Oncol. 13:1220038. doi: 10.3389/fonc.2023.1220038
Received
12 May 2023
Accepted
31 July 2023
Published
18 August 2023
Volume
13 - 2023
Edited by
Ayaz Shahid, Western University of Health Sciences, United States
Reviewed by
Runtian Wang, Nanjing Medical University, China; Pengcheng Wang, Capital Medical University, China
Updates
Copyright
© 2023 Zhang, Wang, Huang, Wang, Liu, Zhou, Tang, Wang, Li and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Geng Wang, 12190816@qq.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.