Frontiers journals are at the top of citation and impact metrics

Original Research ARTICLE

Front. Med., 22 August 2016 | https://doi.org/10.3389/fmed.2016.00034

HLA-DQ B1*0201 and A1*0102 Alleles Are Not Responsible for Antituberculosis Drug-Induced Hepatotoxicity Risk in Spanish Population

imageVirginia Leiro-Fernández1*, imageDiana Valverde2, imageRafael Vázquez-Gallardo1, imageLucía Constenla-Caramés3, imageVíctor del Campo-Pérez4 and imageAlberto Fernández-Villar1
  • 1NeumoVigol+i Research Group, Pulmonary Department, EOXI Vigo, Biomedical Research Institute Vigo (IBIV), Alvaro Cunqueiro Hospital, Vigo, Spain
  • 2Biochemical, Genetics and Immunology Department, Biology Faculty, Vigo University, Vigo, Spain
  • 3Clinic Analysis Department, EOXI Vigo, Alvaro Cunqueiro Hospital, Vigo, Spain
  • 4Preventive Medicine Department, EOXI Vigo, Alvaro Cunqueiro Hospital, Vigo, Spain

Aims: To evaluate the role of human leukocyte antigen (HLA) class II DQB1*0201 and DQA1*0102 in the risk of antituberculosis drug (ATD)-induced hepatotoxicity (ATDH) in a cohort of tuberculosis patients of Caucasian origin from Spain.

Methods: Matched case-control study including active tuberculosis (TB) patients from Spain (Caucasian) treated with first-line ATD (Isoniazid, Rifampin, and Pyrazinamide). Presence or absence of HLA class II DQB1*0201 and DQA1*0102 alleles were compared between cases and controls.

Results: We included 110 TB patients, 55 ATDH cases, and 55 sex-matched controls. The analysis of the presence of HLA-DQB1*0201 and HLA-DQA*0102 did not show significative differences between both groups [presence of HLA-DQB1*0201 53.6% of the cases vs. 45.4% of the controls, OR: 1.63 95% CI (0.62–4.52) p = 0.38; presence of HLA-DQA*0102 7.5% of cases vs. 20% of controls, OR: 0.36 95% CI (0.08–1.23) p = 0.12]. After multivariate logistic regression analysis including in the model, other potential risk factors of hepatotoxicity HLA class II DQB1*0201 and DQA1*0102 alleles were not found significantly associated with the risk of development ATDH. We could not demonstrate an association between HLA-DQA1*0102 and HLA-DQB1*0201 with the risk of ATDH in this Caucasian population of Spanish origin.

Introduction

Tuberculosis (TB) continues to be a global health problem despite the availability of effective therapies. Hepatotoxicity is one of the most important side effects of first-line antituberculosis drugs (ATD) that could lead to increase in morbidity disease because of treatment interruptions and negative impact in treatment compliance including increasing costs (1, 2). The first-line ATD-induced hepatotoxicity (ATDH) is isoniazid (INH), rifampicin (RIF), and pyrazinamide (PZA). Many risk factors for ATDH have been reported, and its incidence depends on the number of risk factors, varying from 5.8 to 24.2% (15). Genetic susceptibility to ATDH has been the subject of study during last years focusing in genes encoding phase I and II metabolic enzymes, and the risk has been significantly associated with certain genetic polymorphisms (3, 611). Human leukocyte antigen (HLA) class I and II drug-induced hepatotoxicity associations have been identified with other drugs and also in healthy individuals with alanine transaminase (ALT) elevations (7). Each HLA associations points to a role for an immune mechanism. Even that the majority of recently published studies on ATDH have been concerned with polymorphisms affecting drug-metabolizing enzymes, there are evidence that some HLA class II alleles are related to the risk of ATDH in the Indian and Chinese population (12, 13). The purpose of this study is to evaluate the risk of ATDH related to HLA class II DQB1*0201 and DQA1*0102 in a cohort of TB patients of Caucasian origin from Spain.

Materials and Methods

Study Design

This individually matched case-control study was nested in a cohort of 1,200 patients from Caucasian origin with active TB that had been treated between January 1998 and December 2008 at the TB Units in Complexo Hospitalario Universitario de Vigo, Pontevedra, and A Coruña (Galicia, Spain). Some characteristics of this cohort were reported in previous publications (5, 8, 9). We included all cases of ATDH that fulfilled the inclusion criteria. The control group consisted of sex- and drug-matched patients that showed no evidence of ATDH, selected consecutively from the same cohort. The study participants should follow the inclusion criteria: (a) age between 15 and 75 years old, (b) microbiological demonstration of active TB, (c) treatment with regimens that included INH, RIF, and PZA at the usual drug dosages (INH 5 mg/kg/day – maximum 300 mg/day, RIF 10 mg/kg/day – maximum 600 mg/day, and PZA 25–30 mg/kg/day – maximum 2,500 mg/day), and (d) adequate compliance. Exclusion criteria were: (a) increased baseline serum transaminases [aspartate aminotransferase (AST) and/or alanine aminotransferase (ALT)], normal values ≤40 IU/L, (b) positive serological testing for the human immunodeficiency virus, hepatitis B virus, or hepatitis C virus, (c) regular alcohol intake or concomitant use of hepatotoxic drugs, (d) history of chronic liver disease, (e) pregnancy, or (f) no or poor adherence to treatment.

Diagnosis of Antituberculosis Drug-Induced Hepatotoxicity

Antituberculosis drug-induced hepatotoxicity was defined as an increase in serum transaminase (either aspartate aminotransferase or alanine aminotransferase) to values higher than three times the upper limit of normal (ULN) (i.e., >120 IU/L) at any time during the treatment period (1). The degree of severity of hepatotoxicity was gaged by the peak level of serum transaminases, and classified according to the Toxicity Classification Standards (10).

Standard Follow-up of Patients on Antituberculosis Drug Therapy

In this cohort, patients on antituberculosis drug therapy received a routine follow-up of clinical assessments every 2 weeks during the first month, and thereafter, monthly until completion of therapy. Routine liver function tests [blood cell counts and biochemical assays for serum AST, ALT, gamma glutamyl transpeptidase (GGT), and bilirubin (BIL)] were performed 15–30 days after the beginning of treatment, at the end of 2 and 4 months, and at 2-month intervals thereafter, when treatments lasted for more than 6 months. Clinical and laboratory check-ups were performed more frequently when hepatitis symptoms or abnormal serum transaminase levels were observed. All patients were warned verbally and in writing of possible hepatotoxicity symptoms and were provided with an emergency contact telephone number. Adherence to treatment was verified by detection of INH metabolites in urine specimens and by tablet counts at every visit.

Genotyping Protocol

Blood was collected by venipuncture, and genomic DNA was isolated with Flexi Gene DNA Kit (QIAGEN) following the instructions of manufacturer (Supplementary Material). The HLA class II region (second exon) was amplified by polymerase chain reaction (PCR). Specifically DQB1*0201 exon 2 genotype was detected by PCR as previously described with minor modifications (14, 15).

Statistical Analyses

Qualitative variables were expressed as percentages and absolute frequencies, and quantitative variables as the median and interquartile range (IQR) or median and SD. Qualitative variables were compared using the McNemar test. Quantitative variables were analyzed using the Wilcoxon test. A logistic regression model was developed in order to evaluate the risk of ATD-induced hepatitis in the presence or absence of HLA DQB1*0201 and DQA1*0201 alleles. Variables were adjusted for potential confounders, including age (in years), sex, nutritional status evaluated by the body mass index [BMI = weight/(height in meters)2], and baseline transaminase values. The patient’s sample size was determined from allele frequencies observed for healthy individuals with a genetic model analyzing the frequency for the disease gene with a RR value = 2 (α = 0.05). The statistical power for variant alleles calculated for the sample size using healthy control population as a reference is 80% for unilateral association.1 A p value less than 0.05 was considered statistically significant.

Ethical Issues

This investigation followed the tenets of the Helsinki declaration. It was approved by the Ethical Committee of Galicia (CEIC). Written informed consent was obtained from all participants. All participants were adults aged more than 18 years.

Results

We included 110 TB patients, 55 with ATDH, and 55 controls with no evidence of this complication. The baseline characteristics are showed in Table 1. Hepatotoxicity was mild in 34% [median peak (IQR) of ALT, AST, GGT, and BIL of 154 IU/L (133–181), 94 IU/L (72–117), 14 IU/L (30–59), and 0.4 mg/dL (0.3–0.65), respectively], moderate in 42.6% [median peak (IQR) of ALT, AST, GGT and BIL of 260 IU/L (220–337), 130 IU/L (108–198), 87 IU/L (46–120), and 0.45 mg/dL (0.3–0.53), respectively] and severe in 23% [median peak (IQR) of ALT, AST, GGT, and BIL of 953.5 IU/L (699–1126), 647 IU/L (239–826), 73 IU/L (61.2–139), and 0.45 mg/dL (0.3–1.01), respectively]. The median length of time from the beginning of treatment until the development of hepatotoxicity was 45 days (IQR, 22–69 days).

TABLE 1
www.frontiersin.org

Table 1. Baseline characteristics of cases and controls.

HLA Analysis

HLA DQA*0102 analysis was completed in 53 ATDH patients and in 53 controls. The frequency of the presence of HLA DQA*0102 was 7.5% in ATDH patients and 20% in non-ATDH patients with no significant differences between cases and controls (OR: 0.36 95% CI 0.08–1.23, p = 0.12) (Table 2). The analysis of the presence of HLA DQB1*0201 allele, completed in 110 patients, did not show differences between both groups of ATDH cases and non-ATDH controls [53.6 vs. 45.4%, respectively, OR: 1.63 95% CI (0.62–4.52) p = 0.38] (Table 2). When we performed the analysis by hepatotoxicity degree, we did not show differences in HLA DQA*0102 and HLA DQB1*0201 presence in the different subgroups (mild, moderate, and severe hepatotoxicity). After multivariate logistic regression analysis including in the model, other potential risk factors of hepatotoxicity as advancing age, BMI, and pre-treatment transaminase values alterations in HLA class II DQB1*0201 and DQA1*0102 alleles were not found significantly associated with the risk of development ATDH.

TABLE 2
www.frontiersin.org

Table 2. HLA II DQB1*0201 and DQA1*0102 analysis in patients with TB and with or without hepatotoxicity.

Discussion

This study performed in Caucasian origin TB patients showed no significant association between HLA DQA1*0102 absence and the risk of hepatotoxicity when the first-line standard therapy is used (INH, RIF, and PZA). However, there is a non-significant predominance of its absence in cases (92.5 vs. 80%). We also did not find an increased ATDH risk in HLA DQB1*0201 carriers. Most of the genetic risk in ATDH is focused in genes encoding drug-metabolizing enzymes that could explain interindividual variation in response and toxicity with ATD (611). However, other pathways of adverse drug reactions like hypersensitivity reactions should be deeply studied in ATDH. Different hypothesis may explain immune-mediated reactions, and the major histocompatibility complex (MHC) class I or II molecules are always involved due to its crucial role in immune recognition (11). Interestingly, HLA variants have been linked to autoimmune hepatitis (7, 11, 12). Association of HLA variants has been described with different drugs as amoxicilin-clavulanate (16, 17), flucoxacilin (11), ximelagatran (11), lapatinib (18), lumiracoxib (19, 20), and neviparine (21). Sharma et al. showed a fourfold increased risk of ATDH in patients without HLADQA1*0102 haplotype suggesting to play a role in pathogenic development of ATDH (12).

Our investigation did not show association in an ethnically different group of TB patients, with a Caucasian origin from Spain. In a case-control study from China, exploring HLADQB1 alleles found that only DQB1*05/*05 positive genotype patients were in risk of ATDH after being adjusted by other protective drugs and weight. Sharma et al. also found, relevantly, the presence of HLADQB1*0201 in ATDH risk. In our study, these finding were not replicated too. An explanation could be the different ethnical origin of Sharma et al. study population. Differences in HLA allele associations have been recently demonstrated between ethnical groups (7). There are few data to support an association between HLA alleles and ATDH risk; however, the results from the studies are heterogeneous probably due to small sample sizes and different ethnical origin (7).

Singer et al. have been described for the selective cyclooxygenase-2 inhibitor lumiracoxib an increased risk of liver disease secondary to lumiracoxib treatment linked to the presence of HLAQA1*0102 allele (OR = 5.0, 95% CI 3.6–7.0) (19). Probably, differences could be explained by the different hepatotoxicity mechanism of non-steroidal drugs, and the interactions between several factors inherent to subjacent disease and the patients.

One of the most important issues is to analyze the convenience of determining HLADQA1*0102 to predict possible hepatotoxic reaction to ATD that could lead to closer clinical and analytic monitoring of TB patients. However, after our results, we cannot make specific recommendations. The purpose of our study was not to analyze its relevance but rather a possible association between HLADQA1*0102 and the risk of ATDH in absence of other previously known ATDH risk factors. For that reason, the criteria of hepatotoxicity definition is low (three times de ULN AST/ALT) and not clinically relevant in most of the cases (1). Variations in the definitions of hepatotoxicity may have given rise to different results among studies. Another potential limitation of the study was that the sample size was small. It was due to the very strict inclusion criteria that exclude TB patients with any previously known ATDH risk factor. We are aware that the failure to find differences in could be due to small sample so that the null hypothesis could be rejected when the alternative hypothesis is true (type II error). The analysis should, therefore, be replicated in a larger study with effective sample size, i.e., containing proband number that achieves adequate statistical power. Finally, we must consider the possibility of frequent bias of genetic case control. This was a nested case-control study design. However, it utilizes the exposure and the confounder data originally collected prior to the onset of the disease, thus reducing potential recall bias and temporal ambiguity (minor concerns in genetic case-control studies). Furthermore, nested studies include cases and controls drawn from the same cohort, thus decreasing the likelihood of selection bias.

Conclusion

We can conclude that even the absence of HLADQA1*0102 could be a risk factor associated with ATDH risk, we cannot demonstrate this association in a cohort of TB patients from Caucasian origin Spanish. We also could not demonstrate an association between HLADQB1*0201 and ATDH. More genetics investigations should be performed in conjunction with the rest of the genetic and clinical ATDH risk factors to determine its routine practical application in TB patients subsidiaries of treatment with first-line ATD.

Author Contributions

VL-F had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Study concept and design: VL-F, DV, AF-V, and VC-P. Acquisition of data: VL-F, RV-G, and AF-V Analysis and interpretation of data: LC-C, DV, VL-F, and VC-P. Drafting of the manuscript: VL-F. Critical revision of the manuscript for important intellectual content: DV, VC-P, and AF-V. Statistical analysis: VL-F and VC-P. All authors declare that they have participated actively in the study, all the authors have read and approved the submitted manuscript, the study complies with current ethical considerations, the manuscript reports unpublished work that is not currently under consideration elsewhere and will not be submitted to another journal until a final decision has been made, and there are no conflicts of interest, unless those identified specifically by the authors in the letter.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

The authors thank N. Chouciño, M. Otero, and A. Gorís, Complexo Hospitalario Universitario de Vigo, L. Anibarro and A. Pena, Complexo Hospitalario de Pontevedra, and D. Díaz-Cabanela, Complexo Hospitalario Universitario A Coruña for their cooperation.

Funding

This research was supported by a predoctoral grant from Sociedad Española de Neumología (SEPAR), grant PIGDIT05SAN21PR from Xunta de Galicia, grant PI052461 from Fondos de Investigación Sanitaria (FIS), and funding from EU Seventh Framework Program [FP7/REGPOT-2012-2013-1; Code: N° 316265, Biomedical Capacities (BIOCAPS)].

Supplementary Material

The Supplementary Material for this article can be found online at https://www.frontiersin.org/article/10.3389/fmed.2016.00034

Abbreviations

ALT, alanine transaminase; AST, aspartate transaminase; ATD, antituberculosis drugs; ATDH antituberculosis drugs-induced hepatotoxixty; BIL, bilirubin; BMI, body mass index; CEIC, Ethical Committee of Galicia; GGT, gamma glutamyl transpeptidase; HLA, human leukocyte antigen; INH, isoniazid; IQR, interquartile range; PCR, polymerase chain reaction; PZA, pyrazinamide; RIF, rifampicin; SBT, sequence based typing; SD; TB, tuberculosis; IU, international units; ULN, upper limit of normal.

Footnote

  1. ^Available from: http://www.allelefrequencies.net/. Consulted May 2016.

References

1. Saukkonen JJ, Cohn DL, Jasmer RM, Schenker S, Jereb JA, Nolan CM, et al. Hepatotoxicity of Antituberculosis Therapy Subcommittee. An official ATS statement: hepatotoxicity of antituberculosis therapy. Am J Respir Crit Care Med (2006) 174:935–52. doi: 10.1164/rccm.200510-1666ST

CrossRef Full Text | Google Scholar

2. Yee D, Valiquette C, Pelletier M, Parisien I, Rocher I, Menzies D. Incidence of serious side effects from first-line antituberculosis drugs among patients treated for active tuberculosis. Am J Respir Crit Care Med (2003) 167:1472–7. doi:10.1164/rccm.200206-626OC

CrossRef Full Text | Google Scholar

3. Cai Y, Yi J, Zhou C, Shen X. Pharmacogenetic study of drug-metabolising enzyme polymorphisms on the risk of anti-tuberculosis drug-induced liver injury: a meta-analysis. PLoS One (2012) 7:e47769. doi:10.1371/journal.pone.0047769

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Fernández-Villar A, Sopeña B, Vázquez R, Ulloa F, Fluiters E, Mosteiro M, et al. Isoniazid hepatotoxicity among drug users: the role of hepatitis C. Clin Infect Dis (2003) 3:293–8. doi:10.1086/345906

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Fernández-Villar A, Sopeña B, Fernández-Villar J, Vázquez-Gallardo R, Ulloa F, Leiro V, et al. The influence of risk factors on the severity of anti-tuberculosis drug-induced hepatotoxicity. Int J Tuberc Lung Dis (2004) 8:1499–505.

PubMed Abstract | Google Scholar

6. Sun F, Chen Y, Xiang Y, Zhan S. Drug-metabolising enzyme polymorphisms and predisposition to anti-tuberculosis drug-induced liver injury: a meta-analysis. Int J Tuberc Lung Dis (2008) 12:994–1002.

PubMed Abstract | Google Scholar

7. Metushi I, Uetrecht J, Phillips E. Mechanism of isoniazid-induced hepatotoxicity: then and now. Br J Clin Pharmacol (2016) 81:1030–6. doi:10.1111/bcp.12885

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Leiro V, Fernández-Villar A, Valverde D, Constenla L, Vázquez R, Piñeiro L, et al. Influence of glutathione S-transferase M1 and T1 homozygous null mutations on the risk of antituberculosis drug-induced hepatotoxicity in a Caucasian population. Liver Int (2008) 28:835–9. doi:10.1111/j.1478-3231.2008.01700.x

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Leiro-Fernandez V, Valverde D, Vázquez-Gallardo R, Botana-Rial M, Constenla L, Agúndez JA, et al. N-acetyltransferase 2 polymorphisms and risk of anti-tuberculosis drug-induced hepatotoxicity in Caucasians. Int J Tuberc Lung Dis (2011) 15:1403–8. doi:10.5588/ijtld.10.0648

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Tostmann A, Boeree MJ, Aarnoutse RE, de Lange WC, van der Ven AJ, Dekhuijzen R. Antituberculosis drug-induced hepatotoxicity: concise up-to-date review. J Gastroenterol Hepatol (2008) 23:192–202. doi:10.1111/j.1440-1746.2007.05207.x

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Russmann S, Jetter A, Kullak-Ublick GA. Pharmacogenetics of drug-induced liver injury. Hepatology (2010) 52:748–61. doi:10.1002/hep.23720

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Sharma SK, Balamurugan A, Saha PK, Pandey RM, Mehra NK. Evaluation of clinical and immunogenetic risk factors for the development of hepatotoxicity during antituberculosis treatment. Am J Respir Crit Care Med (2002) 166:916–9. doi:10.1164/rccm.2108091

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Chen R, Zhang Y, Tang S, Lv X, Wu S, Sun F, et al. The association between HLA-DQB1 polymorphism and antituberculosis drug-induced liver injury: a Case-Control Study. J Clin Pharm Ther (2015) 40:110–5. doi:10.1111/jcpt.12211

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Sacchetti L, Tinto N, Calcagno G, Improta P, Salvatore F. Multiplex PCR typing of the three most frequent HLA alleles in celiac disease. Clin Chim Acta (2001) 310:205–7. doi:10.1016/S0009-8981(01)00575-7

CrossRef Full Text | Google Scholar

15. Robinson J, Waller MJ, Fail SC, McWilliam H, Lopez R, Parham P, et al. The IMGT/HLA database. Nucleic Acids Res (2009) 37:D1013–7. doi:10.1093/nar/gkn662

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Stephens C, López-Nevot MÁ, Ruiz-Cabello F, Ulzurrun E, Soriano G, Romero-Gómez M, et al. HLA alleles influence the clinical signature of amoxicillin-clavulanate hepatotoxicity. PLoS One (2013) 8:e68111. doi:10.1371/journal

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Lucena MI, Molokhia M, Shen Y, Urban TJ, Aithal GP, Andrade RJ, et al. Susceptibility to amoxicillin-clavulanate-induced liver injury is influenced by multiple HLA class I and II alleles. Gastroenterology (2011) 141:338–47. doi:10.1053/j.gastro.2011.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Spraggs CF, Budde LR, Briley LP, Bing N, Cox CJ, King KS, et al. HLA-DQA1*02:01 is a major risk factor for lapatinib-induced hepatotoxicity in women with advanced breast cancer. J Clin Oncol (2011) 29:667–73. doi:10.1200/JCO.2010.31.3197

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Singer JB, Lewitzky S, Leroy E, Yang F, Zhao X, Klickstein L, et al. A genome-wide study identifies HLA alleles associated with lumiracoxib-related liver injury. Nat Genet (2010) 42:711–4. doi:10.1038/ng.632

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Chitturi S, Farrell GC. Identifying who is at risk of drug-induced liver injury: is human leukocyte antigen specificity the key? Hepatology (2011) 53:358–62. doi:10.1002/hep.24094

CrossRef Full Text | Google Scholar

21. Phillips E, Bartlett JA, Sanne I, Lederman MM, Hinkle J, Rousseau F, et al. Associations between HLA-DRB1*0102, HLA-B*5801, and hepatotoxicity during initiation of nevirapine-containing regimens in South Africa. J Acquir Immune Defic Syndr (2013) 62:e55–7. doi:10.1097/QAI.0b013e31827ca50f

CrossRef Full Text | Google Scholar

Keywords: DQB1*0201, DQA1*0102, hepatitis, pharmacogenetics susceptibility, first-line antituberculosis drugs

Citation: Leiro-Fernández V, Valverde D, Vázquez-Gallardo R, Constenla-Caramés L, del Campo-Pérez V and Fernández-Villar A (2016) HLA-DQ B1*0201 and A1*0102 Alleles Are Not Responsible for Antituberculosis Drug-Induced Hepatotoxicity Risk in Spanish Population. Front. Med. 3:34. doi: 10.3389/fmed.2016.00034

Received: 16 February 2016; Accepted: 29 July 2016;
Published: 22 August 2016

Edited by:

Anne Hilgendorff, Ludwig Maximilian University of Munich, Germany

Reviewed by:

Martin Petrek, Palacký University, Czech Republic
Christophe Guignabert, French Institute of Health and Medical Research (INSERM), France

Copyright: © 2016 Leiro-Fernández, Valverde, Vázquez-Gallardo, Constenla-Caramés, del Campo-Pérez and Fernández-Villar. 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) or licensor 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: Virginia Leiro-Fernández, virginia.leiro.fernandez@sergas.es