- 1University Hospital of Santiago de Compostela, Servizo Galego de Saúde (SERGAS), Santiago de Compostela, Spain
- 2Health Research Institute of Santiago de Compostela, Santiago de Compostela, Spain
Diffuse Large B-cell Lymphoma (DLBCL) is the most common type of aggressive lymphoma. Approximately 60% of fit patients achieve curation with immunochemotherapy, but the remaining patients relapse or have refractory disease, which predicts a short survival. Traditionally, risk stratification in DLBCL has been based on scores that combine clinical variables. Other methodologies have been developed based on the identification of novel molecular features, such as mutational profiles and gene expression signatures. Recently, we developed the LymForest-25 profile, which provides a personalized survival risk prediction based on the integration of transcriptomic and clinical features using an artificial intelligence system. In the present report, we studied the relationship between the molecular variables included in LymForest-25 in the context of the data released by the REMoDL-B trial, which evaluated the addition of bortezomib to the standard treatment (R-CHOP) in the upfront setting of DLBCL. For this, we retrained the machine learning model of survival on the group of patients treated with R-CHOP (N=469) and then made survival predictions for those patients treated with bortezomib plus R-CHOP (N=459). According to these results, the RB-CHOP scheme achieved a 30% reduction in the risk of progression or death for the 50% of DLBCL patients at higher molecular risk (p-value 0.03), potentially expanding the effectiveness of this treatment to a wider patient population as compared with other previously defined risk groups.
Introduction
Diffuse Large B-cell Lymphoma (DLBCL) is the most common type of aggressive lymphoma. Approximately 60% of patients achieve curation with the standard first line treatment, which is based on the combination of an anti-CD20 antibody (rituximab) with chemotherapy (cyclophosphamide, doxorubicin and vincristine) and prednisone (R-CHOP). The remaining patients have either refractory disease or relapse after achieving a remission, and this predicts an adverse prognosis (1). Traditionally, risk stratification has been based on scores that combine the value of different prognostic variables. Examples of these methods are the International Prognostic Index (IPI), the revised IPI (R-IPI), and the National Comprehensive Cancer Network IPI (NCCN-IPI) (2). Nevertheless, the accuracy of these scores is far from optimal, and other strategies are actively being explored based on novel molecular features. Earlier studies based on transcriptomic signatures revealed 3 prognostic groups based on their cell-of-origin (COO) status: activated B-cell–like (ABC), germinal-center B-cell–like (GCB) and unclassified (3). Furthermore, recent research reports proved that high risk lymphomas can also be identified as those which share a gene expression signature with either double & triple-hit DLBCLs or with Burkitt lymphomas (4). These lymphomas have been termed as molecular high-grade (MHG) by the academics (4). Finally, comprehensive classifications of DLBCL based on patterns of somatic mutations also exist, which are also associated with divergent clinical outcomes (5).
A few years ago, we presented a new prognostic tool in DLBCL based on a 102-gene expression profile (6). The data from this profile, when interpreted with machine learning tools, enabled the inference of personalized survival outcomes that were prognostically superior to those of the COO classification. Afterwards, we reproduced this profile in another cohort, and reduced the total number of genes in the signature to 19 variables which were prognostically independently of the IPI-related variables (7). The model was named LymForest-25. Finally, we validated the prognostic value of this signature in the UK population-based Haematological Malignancy Research Network database (8), confirming its superiority with respect to the COO and MHG classifications. Notably, the performance of the predictor continued to be high despite the exclusion of 2 genes which were not represented in the gene expression panel used in that study.
At the same time, a growing interest for improved treatments in DLBCL has emerged, and several trials have evaluated new upfront combinations during the last years. The ROBUST study was a randomized, phase III trial which explored the addition of lenalidomide to R-CHOP (R2-CHOP) vs standard R-CHOP, but failed to provide significant results (9). However, a tendency for an improved progression-free survival (PFS) with R2-CHOP was observed among patients with high risk disease (IPI ≥ 3). More recently, the POLARIX phase III trial evaluated a modified scheme of R-CHOP (pola-R-CHP), in which vincristine was replaced with polatuzumab vedotin, as compared with standard R-CHOP in patients with previously untreated intermediate-risk or high-risk DLBCL (10). A significant benefit in PFS was observed in the pola-R-CHP treatment branch, with a hazard ratio (HR) of 0.73. Notably, exploratory subgroup analysis evidenced that this benefit was more pronounced among patients with IPI ≥ 3 and in those with ABC phenotype. A different therapeutic strategy has been based on the incorporation of the proteasome inhibitor bortezomib into the R-CHOP scheme (RB-CHOP). Preclinical evidence indicated that bortezomib can exert antitumoral activity in B-cell lymphoma cell lines (11). This promoted clinical studies that ended up in the development of the REMoDL-B trial, a randomized phase III trial testing RB-CHOP vs R-CHOP in previously untreated DLBCL patients (12). The results of this trial indicated no evidence for a benefit of RB-CHOP over R-CHOP neither in PFS nor in overall survival (OS). However, exploratory post-hoc analysis evidenced a benefit for RB-CHOP in PFS in the MHG group, and a tendency towards a benefit in the ABC group (13).
Methods
In the present report, we aimed to reproduce the prognostic value of the LymForest gene expression profile in the publicly available data of the REMoDL-B trial (13), as well as to evaluate the possible predictive value of this signature. Briefly, normalized gene expression estimates were downloaded from the Gene Expression Omnibus (GEO), with ID GSE117556. This cohort contained data for 928 patients, out of which 469 were treated with R-CHOP and 459 were treated with RB-CHOP. Median follow-up was 29.37 months, and median overall survival was not reached. We created a random forest model of survival following previous specifications (7, 8), and this model was exclusively trained on the group of patients treated with R-CHOP. Out-of-bag metrics were derived for patients in this subgroup, and new predictions on patients treated with RB-CHOP were made based on the results of the training set. The values of the cumulative hazard function were used to calculate the c-indexes.
Results
Firstly, we decided to reproduce the machine learning predictions based on the expression of 17 out of 19 original genes. This was due to the fact that 2 genes (FAM208B and TRAV6) were not included in the Illumina HumanHT-12 WG-DASL V4.0 R2 expression beadchips. In the original UK population-based Haematological Malignancy Research Network database, the c-index of this signature was 0.612. In the case of the REMoDL-B trial cohort, the c-indexes were 0.619 and 0.640 for the R-CHOP and RB-CHOP treated patients, respectively. Then, we reasoned that the expression of one of the missing genes (the T-cell receptor alpha subunit variable region gene; TRAV6), could be substituted by the expression of the CD3 T-cell specific marker genes, namely CD3D, CD3G and CD3E. We observed that this strategy improved the c-index in the UK population-based Haematological Malignancy Research Network database (original c-index, 0.612; new c-index, 0.621). Hence, we performed the same modification in the REMoDL-B trial cohort, obtaining a c-index of 0.668 in the group of patients treated with R-CHOP, and a small reduction of the c-index to 0.631 in those treated with RB-CHOP.
In a second approach, we evaluated the possible predictive value of this signature in the REMoDL-B trial. With this aim, we extracted the 2-year survival probabilities from the machine learning predictions. We chose this threshold because most of the relapses and lymphoma-related deaths are known to occur during this period of time (14). Initially, we explored the possible utility of the 17-gene model by splitting the patients into 2 halves and 3 tertiles of risk (Table 1). No statistically significant difference in PFS was observed between patients treated with RB-CHOP and R-CHOP in either the high or the low 50% risk groups. However, a significant advantage of RB-CHOP for those patients assigned to the higher 33% of risk was identified (p-value 0.03, HR 0.66). Then, we reproduced the same procedure with the model enriched in T-cell markers (Table 1). In this case, we observed a significantly higher PFS with RB-CHOP for the 50% of patients at higher risk (p-value 0.03, HR 0.70; Figure 1A), whereas no significant differences were observed for those patients assigned to the lower 50% of risk (Figure 1B). This effect appeared to be even more pronounced among patients in the higher 33% of risk (p-value 0.03, HR 0.66) (Figures 1C, D).
Table 1 Results of the cox models testing for differential PFS outcomes of the different groups of patients analyzed in the text.
Figure 1 Kaplan-Meier curves representing the PFS of patients treated with RB-CHOP and R-CHOP according to their biological risk predicted by the CD3 markers & 17 gene expression signature. (A, B) Representation of RB-CHOP and R-CHOP curves for those patients in the higher 50% risk group (A) and in the lower 50%risk group (B). (C, D) Representation of RB-CHOP and R-CHOP curves for those patients in the higher 33% risk group (A) and in the lower 67%risk group (B).
Discussion
Our data indicates a role for bortezomib-containing upfront treatments in patients with DLBCL who have high-risk molecular features. According to these results, the RB-CHOP scheme achieved a 30% reduction in the risk of progression or death for the 50% of DLBCL patients at higher molecular risk, potentially expanding the effectiveness of this treatment to a wider patient population as compared with other previously defined risk groups. Furthermore, we confirmed that the inclusion of T-cell markers in the gene expression signature enriches the prognostic performance of the signature in patients treated with R-CHOP, although their importance appears to diminish in patients treated with RB-CHOP. In conclusion, the standardization and implementation of machine learning-guided molecular risk scores based on transcriptomic features should be performed in the context of clinical trials evaluating novel upfront combinations in the upfront treatment of DLBCL. Additionally, the LymForest molecular profile improves previous transcriptomic signatures for both prognostication and drug-response prediction in patients with DLBCL requiring systemic immunochemotherapy. This strategy could also be explored to enrich the results of other trials aiming to improve R-CHOP as upfront treatment in DLBCL, such as those based on polatuzumab (pola-R-CHP) (10) and those aiming to incorporate immunotherapy (bispecific antibodies or CAR-T cells) (15, 16). This is relevant because most of the new drugs in the frontline setting face a substantial difficulty to improve R-CHOP due to its high effectivity in the global population, and therefore the development of biomarkers to guide their use is of the utmost interest.
Machine learning has the potential to play an important role in the post hoc analysis of clinical trials by enabling more comprehensive and accurate analysis of trial data. AI-based methods can process large amounts of data and identify patterns, relationships, and insights that may not be immediately apparent through traditional statistical methods. Additionally, AI techniques can help to identify potential safety concerns, optimize dosing regimens, and identify subgroups of patients who may benefit the most from a particular treatment. For example, machine learning algorithms can be used to identify the best predictive biomarkers or clusters of patients for a particular treatment, which can inform future trial design and clinical decision making. Several studies have demonstrated the potential of AI in the post hoc analysis of clinical trials. For instance, a recent study by Yan et al. (2021) used a machine learning algorithm to predict clinical outcomes in patients with colorectal cancer treated with immunotherapy, achieving better performance than traditional statistical methods (17). Recently, newer approaches in the prediction of response to targeted drugs and drug combinations from patients treated in routine clinical practice have been presented. For example, Kong et al. (2022) presented an approach to predict the response to immune check-point inhibitors based on the construction of a network of genes and proteins that are known to be involved in the immune response (18). Using machine learning algorithms, they identified patterns in the network that were associated with drug response, proving that their approach can accurately predict responses in several different types of cancer. In another approach, our group explored new methods to predict risk in multiple myeloma (MM) by the integration of clinical and biochemical data with gene expression profiling. By applying machine learning algorithms, we created a 50-variable model that can predict OS with high concordance (19). The model included patient age, ISS stage, serum B2-microglobulin, first-line treatment, and the expression of 46 genes as covariates. Importantly, we found that patients treated with the best-predicted drug combination were significantly less likely to die than patients treated with other schemes, particularly those treated with a triplet combination including bortezomib, an immunomodulatory drug and dexamethasone.
Validated and transparent machine learning algorithms are essential in medical applications because they can provide accurate and reliable predictions, which can aid clinicians in making optimal decisions (20). Despite this great potential, it is important to recognize their limitations and potential biases. It is crucial to fully understand the strengths and weaknesses of each algorithm and to ensure that they are appropriately validated and transparent. This requires ongoing research and collaboration between machine learning experts and clinicians (21). In the particular context of clinical trials, post hoc analysis can be used to analyze data from clinical trials and determine if a drug is effective and safe (22). Machine learning can be used to identify patterns and relationships that can later be used to develop new drugs or optimize existing ones. However, post hoc studies also have limitations, including the possibility of data overfitting and the inability to control for confounding variables due to their retrospective nature (23). In our particular case, we retrained a previously described prognostic model in the group of patients treated with R-CHOP, because this population was the target of our predictor (6–8). Afterwards, we used these predictions to risk stratify those patients treated with RB-CHOP and compare outcomes. However, though the use of an external cohort for training the predictor could be an option, this should have the same structure (e.g., inclusion and exclusion criteria, baseline characteristics…) as the original clinical trial data. This highlights the crucial importance of external validity in clinical trials, particularly for the construction of new machine learning predictors, and the need to discuss these issues with regulatory agencies for drug approval based on such types of evidence (24). Surely, a prospective validation of the findings on a new clinical trial or in a post-authorization real world cohort would provide more reliable information. Eventually, the growing application of machine learning in clinical trials will make these post hoc analysis more relevant, and regulators should pursue the development of good clinical practices to ensure a reliable and fair application (25, 26). This includes using appropriate statistical methods, validating the model on multiple datasets, and being transparent about its possible limitations. Another issue of relevance for the application of this technology relies on the need for genomic standardization, which should be pursued in order to lead to reliable results for patient care. Standardization of genomic tests involves ensuring that the tests are performed in a consistent and reliable manner across different laboratories, using well-defined protocols, standardized testing platforms and quality control measures (27). This will help to ensure that the results of the tests are accurate and can be compared across different settings and over time.
In conclusion, we present an evaluation of LymForest-25 machine-learning-based gene expression profile to risk stratify patients and predict treatment responses in patients with DLBCL within the REMoDL-B trial. The results suggest a role for bortezomib-containing upfront treatments in molecular high-risk patients. The standardization and implementation of machine learning-guided molecular risk scores based on transcriptomic features should be pursued in the context of clinical trials evaluating novel upfront combinations in the upfront treatment of DLBCL.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/supplementary material.
Author contributions
AMO had the idea, designed the paper, analyzed the data and wrote the paper. JADA, MCL, APR, LBP, MSGP, MMPE, MFFR and JLBL reviewed the results, suggested modifications and approved the final publication. All authors contributed to the article and approved the submitted version.
Funding
The authors declare that no support was provided in the forms of grants and/or equipment and grants for the development of this study.
Acknowledgments
The authors wish to thank the Supercomputing Center of Galicia (CESGA) for providing computational resources; and the REMoDL-B trial organizers for publicly sharing the clinical and gene expression data.
Conflict of interest
AMO reports honoraria for lectures and participation in advisory boards from Janssen, Takeda, Abbey, Amgen, Novartis, Gilead and AstraZeneca; research grants from Roche, Pfizer and Celgene-BMS and funds for conference organization from Jassen, Takeda, Abbey, Amgen, Novartis, Gilead, Roche, Bristol-Myers-Squibb, Glaxo-Smith-Klyne, Incyte and Pfizer.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be constructed as a potential conflict of interest.
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.
References
1. Feugier P, Van Hoof A, Sebban C, Solal-Celigny P, Bouabdallah R, Fermé C, et al. Long-term results of the r-CHOP study in the treatment of elderly patients with diffuse large b-cell lymphoma: a study by the groupe d’Etude des lymphomes de l’Adulte. J Clin Oncol (2005) 23(18):4117–26. doi: 10.1200/JCO.2005.09.131
2. Ruppert AS, Dixon JG, Salles G, Wall A, Cunningham D, Poeschel V, et al. International prognostic indices in diffuse large b-cell lymphoma: a comparison of IPI, r-IPI, and NCCN-IPI. Blood (2020) 135(23):2041–8. doi: 10.1182/blood.2019002729
3. Scott DW, Wright GW, Williams PM, Lih CJ, Walsh W, Jaffe ES, et al. Determining cell-of-origin subtypes of diffuse large b-cell lymphoma using gene expression in formalin-fixed paraffin-embedded tissue. Blood (2014) 123(8):1214–7. doi: 10.1182/blood-2013-11-536433
4. Sha C, Barrans S, Cucco F, Bentley MA, Care MA, Cummin T, et al. Molecular high-grade b-cell lymphoma: defining a poor-risk group that requires different approaches to therapy [published correction appears in J clin oncol. 2019 Apr 20;37(12):1035]. J Clin Oncol (2019) 37(3):202–12. doi: 10.1200/JCO.18.01314
5. Schmitz R, Wright GW, Huang DW, Johnson CA, Phelan JD, Wang JQ, et al. Genetics and pathogenesis of diffuse Large b-cell lymphoma. N Engl J Med (2018) 378(15):1396–407. doi: 10.1056/NEJMoa1801445
6. Mosquera Orgueira A, Díaz Arias JÁ, Cid López M, Peleteiro Raíndo A, Antelo Rodríguez B, Aliste Santos C, et al. Improved personalized survival prediction of patients with diffuse large b-cell lymphoma using gene expression profiling. BMC Cancer. (2020) 20(1):1017. doi: 10.1186/s12885-020-07492-y
7. Mosquera-Orgueira A, Cid-Lopez M, Peleteiro-Raindo A, Díaz Arias JA, González Pérez MS, Antelo Rodríguez B, et al. LymForest-25: personally-tailored survival prediction of patients with diffuse large b-cell lymphoma using clinico-genomic prognostic models. Supplement: 16th Int Conf Malignant Lymphoma Virtual Edition (2021) 39(2):238. doi: 10.1002/hon.79_2880
8. Mosquera Orgueira A, Díaz Arías JÁ, Cid López M, Peleteiro Raíndo A, López García A, Abal García R, et al. Prognostic stratification of diffuse Large b-cell lymphoma using clinico-genomic models: validation and improvement of the LymForest-25 model. Hemasphere (2022) 6(4):e706. doi: 10.1097/HS9.0000000000000706
9. Nowakowski GS, Chiappella A, Gascoyne RD, Scott DW, Zhang Q, Jurczak W, et al. ROBUST: a phase III study of lenalidomide plus r-CHOP versus placebo plus r-CHOP in previously untreated patients with ABC-type diffuse Large b-cell lymphoma. J Clin Oncol (2021) 39(12):1317–28. doi: 10.1200/JCO.20.01366
10. Tilly H, Morschhauser F, Sehn LH, Friedberg JW, Trněný M, Sharman JP, et al. Polatuzumab vedotin in previously untreated diffuse Large b-cell lymphoma. N Engl J Med (2022) 386(4):351–63. doi: 10.1056/NEJMoa2115304
11. Strauss SJ, Higginbottom K, Jüliger S, Maharaj L, Allen P, Schenkein D, et al. The proteasome inhibitor bortezomib acts independently of p53 and induces cell death via apoptosis and mitotic catastrophe in b-cell lymphoma cell lines. Cancer Res (2007) 67(6):2783–90. doi: 10.1158/0008-5472.CAN-06-3254
12. Davies A, Cummin TE, Barrans S, Maishman T, Mamot C, Novak U, et al. Gene-expression profiling of bortezomib added to standard chemoimmunotherapy for diffuse large b-cell lymphoma (REMoDL-b): an open-label, randomised, phase 3 trial. Lancet Oncol (2019) 20(5):649–62. doi: 10.1016/S1470-2045(18)30935-5
13. Davies AJ, Stanton L, Caddy J, Barrans S, Wilding S, Saunders GN, et al. Five-year survival results from the remodl-b trial (ISRCTN 51837425) show improved outcomes in diffuse Large b-cell lymphoma molecular subgroups from the addition of bortezomib to r-CHOP chemoimmunotherapy. Blood (2022) 140(Supplement 1):1770–2. doi: 10.1182/blood-2022-159976
14. Jakobsen LH, Bøgsted M, Brown PN, Arboe B, Jørgensen J, Larsen TS, et al. Minimal loss of lifetime for patients with diffuse Large b-cell lymphoma in remission and event free 24 months after treatment: a Danish population-based study. J Clin Oncol (2017) 35(7):778–84. doi: 10.1200/JCO.2016.70.0765
15. Dickinson MJ, Carlo-Stella C, Morschhauser F, Bachy E, Corradini P, Iacoboni G, et al. Glofitamab for relapsed or5 refractory diffuse Large b-cell lymphoma. N Engl J Med (2022) 387(24):2220–31. doi: 10.1056/NEJMoa2206913
16. Neelapu SS, Dickinson M, Munoz J, Ulrickson ML, Thieblemont C, Oluwole OO, et al. Axicabtagene ciloleucel as first-line6 therapy in high-risk large b-cell lymphoma: the phase 2 ZUMA-12 trial. Nat Med (2022) 28(4):735–42. doi: 10.1038/s41591-022-01731-4
17. Yan L, Wu Y, Wu Y, Zhang Q, Wang X. Prediction of clinical outcomes for immunotherapy of colorectal cancer using a machine learning algorithm. Front Genet (2021) 12:631321.
18. Kong J, Ha D, Lee J, Kim I, Park M, Im SH, et al. Network-based machine learning approach to predict immunotherapy response in cancer patients. Nat Commun (2022) 13(1):3703. doi: 10.1038/s41467-022-31535-6
19. Mosquera Orgueira A, González Pérez MS, Díaz Arias JÁ, Antelo Rodríguez B, Alonso Vence N, Bendaña López , et al. Survival prediction and treatment optimization of multiple myeloma patients using machine-learning models based on clinical and gene expression data. Leukemia (2021) 35(10):2924–35. doi: 10.1038/s41375-021-01286-2
20. Goldstein BA, Navar AM, Pencina MJ, Ioannidis JP. Opportunities and challenges in developing risk prediction models with electronic health records data: a systematic review. J Am Med Inform Assoc (2017) 24(1):198–208. doi: 10.1093/jamia/ocw042
21. Jha S, Topol EJ. Adapting to artificial intelligence: radiologists and pathologists as information specialists. JAMA (2016) 316(22):2353–4. doi: 10.1001/jama.2016.17438
22. Srinivas TR, Ho B, Kang J, Kaplan B. Post hoc Analyses: after the facts. Transplantation (2015) 99(1):17–20. doi: 10.1097/TP.0000000000000581
23. Curran-Everett D, Milgrom H. Post-hoc Data analysis: benefits and limitations. Curr Opin Allergy Clin Immunol (2013) 13(3):223–4. doi: 10.1097/ACI.0b013e3283609831
24. Futoma J, Simons M, Panch T, Doshi-Velez F, Celi LA. The myth of generalisability in clinical research and machine learning in health care. Lancet Digit Health (2020) 2(9):e489–92. doi: 10.1016/S2589-7500(20)30186-2
25. Hollis S, Fletcher C, Lynn F, Urban HJ, Branson J, Burger HU, et al. Best practice for analysis of shared clinical trial data. BMC Med Res Methodol (2016) 16 Suppl 1(Suppl 1):76. doi: 10.1186/s12874-016-0170-y
26. Liu J, Chow SC. A proposal for Post Hoc subgroup analysis in support of regulatory submission. Ther Innov Regul Sci (2023) 57(2):196–208. doi: 10.1007/s43441-022-00459-0
27. Burke W, Ginsburg G, Terry S, Aronson N, Ashley EA, Billings PR, et al. Roundtable on translating genomic-based research for health; board on health sciences policy; institute of medicine; center for medical technology policy. genome-based diagnostics: demonstrating clinical utility in oncology: workshop summary. Washington (DC: National Academies Press (US (2013) 39(2):238.
Keywords: machie learning, DLBCL - diffuse large B cell lymphoma, bortezomib, R-CHOP, lymphoma, genomics, gene expression
Citation: Mosquera Orgueira A, Díaz Arías JÁ, Serrano Martín R, Portela Piñeiro V, Cid López M, Peleteiro Raíndo A, Bao Pérez L, González Pérez MS, Pérez Encinas MM, Fraga Rodríguez MF, Vallejo Llamas JC and Bello López JL (2023) A prognostic model based on gene expression parameters predicts a better response to bortezomib-containing immunochemotherapy in diffuse large B-cell lymphoma. Front. Oncol. 13:1157646. doi: 10.3389/fonc.2023.1157646
Received: 02 February 2023; Accepted: 17 April 2023;
Published: 28 April 2023.
Edited by:
Francesco Piazza, University of Padua, ItalyReviewed by:
Monica Balzarotti, Humanitas Research Hospital, ItalyLiang Wang, Capital Medical University, China
Copyright © 2023 Mosquera Orgueira, Díaz Arías, Serrano Martín, Portela Piñeiro, Cid López, Peleteiro Raíndo, Bao Pérez, González Pérez, Pérez Encinas, Fraga Rodríguez, Vallejo Llamas and Bello López. 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: Adrián Mosquera Orgueira, YWRyaWFuLm1vc3F1ZXJhLm9yZ2VpcmFAc2VyZ2FzLmVz