Skip to main content

BRIEF RESEARCH REPORT article

Front. Genet., 22 November 2022
Sec. Immunogenetics

Impact of polymorphisms in genes orchestrating innate immune responses on replication kinetics of Torque teno virus after kidney transplantation

Natalia Redondo,
Natalia Redondo1,2*Isabel Rodríguez-Goncer,Isabel Rodríguez-Goncer1,2Patricia ParraPatricia Parra1Eliseo AlbertEliseo Albert3Estela Gimnez,Estela Giménez2,3Tamara Ruiz-MerloTamara Ruiz-Merlo1Francisco Lpez-Medrano,,Francisco López-Medrano1,2,4Rafael San Juan,,Rafael San Juan1,2,4Esther GonzlezEsther González5ngel SevillanoÁngel Sevillano5Amado Andrs,Amado Andrés4,5David Navarro,,David Navarro2,3,6Jos María Aguado,,José María Aguado1,2,4Mario Fernndez-Ruiz,,Mario Fernández-Ruiz1,2,4
  • 1Unit of Infectious Diseases, Hospital Universitario ‘12 de Octubre’, Instituto de Investigación Sanitaria Hospital ‘12 de Octubre’ (imas12), Madrid, Spain
  • 2Centro de Investigación Biomédica en Red de Enfermedades Infecciosas (CIBERINFEC), Instituto de Salud Carlos III, Madrid, Spain
  • 3Department of Microbiology, Instituto de Investigación Sanitaria INCLIVA, Hospital Clínico Universitario, Valencia, Spain
  • 4Department of Medicine, School of Medicine, Universidad Complutense, Madrid, Spain
  • 5Department of Nephrology, Instituto de Investigación Sanitaria Hospital “12 de Octubre” (imas12), Hospital Universitario “12 de Octubre”, Madrid, Spain
  • 6Department of Microbiology, School of Medicine, University of Valencia, Valencia, Spain

Background: Torque teno virus (TTV) DNAemia has been proposed as a surrogate marker of immunosuppression after kidney transplantation (KT), under the assumption that the control of viral replication is mainly exerted by T-cell-mediated immunity. However, Tthe impact on post-transplant TTV kinetics of single genetic polymorphisms (SNPs) in genes orchestrating innate responses remains unknown. We aimed to characterize the potential association between 14 of these SNPs and TTV DNA levels in a single-center cohort of KT recipients.

Methods: Plasma TTV DNAemia was quantified by real-time PCR in 221 KT recipients before transplantation (baseline) and regularly through the first 12 post-transplant months. We performed genotyping of the following SNPs: CTLA4 (rs5742909, rs231775), TLR3 (rs3775291), TLR9 (rs5743836, rs352139), CD209 (rs735240, rs4804803), IFNL3 (rs12979860, rs8099917), TNF (rs1800629), IL10 (rs1878672, rs1800872), IL12B (rs3212227) and IL17A (rs2275913).

Results: The presence of the minor G allele of CD209 (rs4804803) in the homozygous state was associated with undetectable TTV DNAemia at the pre-transplant assessment (adjusted odds ratio: 36.96; 95% confidence interval: 4.72–289.67; p-value = 0.001). After applying correction for multiple comparisons, no significant differences across SNP genotypes were observed for any of the variables of post-transplant TTV DNAemia analyzed (mean and peak values, areas under the curve during discrete periods, or absolute increments from baseline to day 15 and months 1, 3, 6 and 12 after transplantation).

Conclusion: The minor G allele of CD209 (rs4804803) seems to exert a recessive protective effect against TTV infection in non-immunocompromised patients. However, no associations were observed between the SNPs analyzed and post-transplant kinetics of TTV DNAemia. These negative results would suggest that post-transplant TTV replication is mainly influenced by immunosuppressive therapy rather than by underlying genetic predisposition, reinforcing its clinical application as a biomarker of adaptive immunity.

Introduction

The study of the human virome in health and disease has gained growing attention over recent years (Webb et al., 2020; Dodi et al., 2021). Viruses belonging to Anelloviridae family are the most abundant eukaryotic viruses in the virome and may be detected in a variety of samples, such as blood, plasma, urine or saliva (Kaczorowska and van der Hoek, 2020; Arze et al., 2021). Anelloviruses are non-enveloped viruses with small circular replication-associated protein-encoding single-stranded DNA genomes (Biagini, 2009; Kaczorowska and van der Hoek, 2020), which lack attributable pathogenic roles (“orphan viruses”) (Focosi et al., 2016; Rezahosseini et al., 2019). Once primary infection occurs at early stages of life, anelloviruses remain in different body compartments and fluids—including peripheral blood mononuclear cells, feces, semen, throat swabs, umbilical cord blood, lungs, kidneys or cerebrospinal fluid—under the control of the immune system, resulting in a prevalence as high as 90% in the adult population (Redondo et al., 2022a). The precise underlying mechanisms on how this immune control is carry out largely remain to be determined, although a major role has been proposed for the cellular arm. Belonging to the Alphatorquevirus genus and discovered in 1997 (Nishizawa et al., 1997), Torque teno virus (TTV) has been proven by us and others to serve as a convenient surrogate marker of the overall status of immunosuppression after solid organ (SOT) and allogeneic hematopoietic stem cell transplantation (HSCT) (Fernandez-Ruiz et al., 2019; Rezahosseini et al., 2019; Mouton et al., 2020; Redondo et al., 2022a; Jaksch et al., 2022).

The innate immunity acts as a frontline defense against viruses through an orchestrated response, that is, triggered upon recognition of viral motifs by pathogen recognition receptors (PRRs) present in macrophages and dendritic cells (Takeuchi and Akira, 2010). The rationale for the use of TTV DNAemia as a biomarker of immune competence after SOT lies on the assumption that the viral kinetics is mainly dictated by the T-cell-mediated immunity (Redondo et al., 2022a; Jaksch et al., 2022). Indeed, various studies have shown a direct correlation between TTV DNA loads and calcineurin inhibitors trough levels (Gorzer et al., 2014; Jaksch et al., 2018). The role played by the innate immune arm in the setting of ongoing immunosuppression remains largely unknown, as is the potential impact of polymorphisms in genes coding for PRRs (such as toll-like receptors [TLRs]), interleukins (IL) or interferons (IFNs) (Prasetyo et al., 2015; Ramzi et al., 2019; Ramzi et al., 2021). Evidence of an individual genetic susceptibility to TTV regardless of the amount of immunosuppressive therapy would question the reliability of viral replication as clinical biomarker in the SOT population.

We aimed to investigate the association between 14 single genetic polymorphisms (SNPs) in different genes mainly involved in the orchestration of innate immune responses (Table 1) and TTV DNA levels at baseline and various points during the first post-transplant year in a well characterized cohort of kidney transplant (KT) recipients (Fernandez-Ruiz et al., 2019). The selection of these SNPs was dictated by previous research showing a potential impact on the susceptibility to viral infections. In the case of TLR3 (rs3775291), various pieces of evidence have shown an effect on the incidence of infection by cytomegalovirus (CMV) or BK polyomavirus (BKPyV), two relevant viral pathogens in the KT scenario, but also tick-borne encephalitis, chikungunya or hepatitis B virus (HBV) (Kindberg et al., 2011; Geng et al., 2016; Studzinska et al., 2017; Fischer et al., 2018; Bucardo et al., 2021; Redondo et al., 2022b). We have previously reported that certain SNPs in TLR9 (rs5743836, rs352139) modulate the risk of CMV infection in two independent cohorts of KT recipients (Fernandez-Ruiz et al., 2015; Redondo et al., 2022b). Regarding SNPs located in the CD209 gene, rs735240 appears to increase the incidence of CMV infection in seropositive KT recipients not receiving antiviral prophylaxis (Fernandez-Ruiz et al., 2015), whereas rs4804803 has been correlated with an increased susceptibility to dengue virus (Vargas-Castillo et al., 2018) and, more recently, BKPyV (Redondo et al., 2022c). We analyzed the SNPs located in IFNL4 (rs12979860, rs8099917) due to its well-established relevance in other viral infections, including CMV (Fernandez-Ruiz et al., 2015) and hepatitis C virus (HCV) (Ge et al., 2009; Thomas et al., 2009). Finally, we aimed to validate in the SOT population the associations reported by other authors between CTLA4 (rs5742909, rs231775), TNF (rs1800629), IL10 (rs1800872, rs1878672), IL12B (rs3212227) and IL17 (rs2275913) SNPs and the kinetics of TTV DNAemia after HSCT (Ramzi et al., 2019; Ramzi et al., 2021).

TABLE 1
www.frontiersin.org

TABLE 1. Candidate SNPs selected for the present study.

Material and methods

Study population and setting

The present research was performed as a post hoc retrospective analysis of a previous study that investigated the role of TTV DNA levels to predict the occurrence of serious and opportunistic infection and de novo malignancy in a cohort of KT recipients recruited at the University Hospital “12 de Octubre” (a 1,300-bed tertiary care center in Madrid with an active KT program since 1990) between November 2014 and December 2016 (Fernandez-Ruiz et al., 2019). As detailed elsewhere, adult patients with end-stage renal disease (ESRD) undergoing KT during the study period and providing informed consent were eligible for inclusion. Exclusion criteria included double organ transplantation and primary graft non-function. By applying these criteria, 221 KT recipients were eventually included. The study was performed in accordance with the ethical standards laid down in the Declarations of Helsinki and Istanbul. The local Clinical Research Ethics Committee approved the study protocol.

Study design

Participants were enrolled at the time of KT and followed-up for at least 12 months, unless graft loss (retransplantation or return to dialysis) or death occurred earlier. Plasma TTV DNA load was quantified at baseline (i.e., within 6 h prior to the transplant procedure), day 7, and months 1, 3, 6 and 12 by a polymerase chain reaction (PCR)-based quantitative nucleic acid amplification test. Immunosuppression and prophylaxis regimens are detailed as Supplementary Material.

Single genetic polymorphisms genotyping

Whole blood specimens that have been stored at −70°C were retrieved for SNP genotyping. DNA was extracted with the KingFisher Duo Prime system using the MagMax DNA Multi-Sample Ultra 2.0 kit (Thermo Fisher Scientific, Waltham, MA) following the manufacturer’s instructions. CTLA4 (rs5742909, rs231775), TLR3 (rs3775291), TLR9 (rs5743836, rs352139), CD209 (rs735240, rs4804803), IFNL3 (rs12979860, rs8099917), TNF (rs1800629), IL10 (rs1878672, rs1800872), IL12B (rs3212227) and IL17A (rs2275913) genotyping was performed by Taqman technology (Thermo Fisher Scientific) in a QuantStudio 3 real-time PCR system (Applied Biosystems, Foster City, CA). SNP and allele calling was made by means of the TaqMan Genotyper Software version 1.0 (Applied Biosystems) and the QuantStudio Design and Analysis Software version 1.5.1 (ThermoFisher Scientific).

Plasma torque teno virus DNA load quantification

TTV DNA extraction and quantification was performed as previously described (Fernandez-Ruiz et al., 2019). Briefly, DNA was extracted from 200 μL of plasma with the NucliSENSR easyMAGR automated system (bioMérieux, Marcy-l’Étoile, France), following the manufacturer’s instructions. DNA loads were quantified by means of a real-time PCR assay targeting a highly conserved segment of the 5′untranslated region of the viral genome (TTV R-gene kit, ARGENE range, bioMérieux). PCR amplification and amplicon detection was performed on an ABI Prism 7500 system (PE Biosystems, Foster City, CA). The viral load (in copy numbers per mL) was determined using a standard curve with known copy numbers and log10-transformed for statistical analyses. The lower limit of detection (LLoD) was 167 copies/mL [95% confidence interval (CI): 92-581] or 2.2 log10 copies/mL (95% CI: 2.0–2.8), with DNA quantitation in the linear range from 2.1 × 102 to 2.1 × 107 copies/ml. Specimens with undetectable DNA loads were assigned a value of 0.01 (−2.0 log10) copies/mL for analysis purposes. All samples from each patient were simultaneously assayed in singlets.

Statistical analysis

Quantitative data were reported as the mean ± standard deviation (SD) or the median with interquartile range (IQR). Qualitative variables were given as absolute and relative frequencies. Normality of the distributions was tested with the Kolgomorov-Smirnov test. Deviation from the Hardy-Weinberg equilibrium for each SNP was evaluated by the χ2 test with one degree of freedom. Comparisons of TTV kinetics at different points across SNP genotypes were performed by the χ2 test or the Fisher’s exact test for qualitative variables (i.e. detectable or undetectable [below the LLoD] DNAemia), or by the T-Student or U-Mann-Whitney tests for continuous variables (i.e. plasma DNA levels). In addition, other viral kinetic parameters were compared across SNPs: peak plasma TTV DNA levels and areas under the curve (AUCs) for TTV DNAemia through discrete time periods (1, 3, and 6 months after transplantation), and increments (Δ) in DNA levels from baseline to day 15 and months 1, 3, 6 and 12. Additional pairwise comparisons were conducted between different SNP genotype groups, either individually or in combination. The independent impact of selected SNPs on the probability of having undetectable TTV DNAemia was confirmed by logistic regression, with associations given as odds ratios (ORs) and 95% CIs. All the significance tests were two-tailed and considered as significant at a p-value < 0.05. To control for p-value inflation due to multiple comparisons, the Bonferroni method (corrected α value = nominal α value/total number of comparisons) was applied. Statistical analysis was performed using SPSS version 21 (Statistical Package for Social Sciences, Chicago, IL).

Results

We included 221 KT recipients, whose demographics, clinical characteristics and patient and graft outcomes are detailed in Table 2. Samples from all the patients were successfully genotyped for the 14 SNPs considered. The median number of assessments for plasma TTV DNA per patient was 5 (IQR: 4–5). The majority of recipients had detectable TTV DNAemia (i.e. above the LLoD) at every time point, ranging from 96.3% (180/187) at baseline to 99.4% (176/177) at post-transplant month 6. The genotypic frequencies of candidate SNPs are shown in Supplementary Table S1. The observed genotype frequency distributions did not deviate from those expected according to the Hardy-Weinberg equilibrium except for TLR9 (rs5743836) and IFNL3 (rs12979860).

TABLE 2
www.frontiersin.org

TABLE 2. Demographic and clinical characteristics of the study cohort (n = 221).

First, the effect of studied polymorphisms on plasma TTV DNAemia at discrete time points was investigated. In particular, we explored the impact of the minor alleles in each SNP in both dominant (heterozygous and homozygous) and recessive (homozygous only) models. Across the 14 SNPs considered, we did not find significant differences in TTV DNA levels at any of the monitoring points (Table 3).

TABLE 3
www.frontiersin.org

TABLE 3. Plasma TTV DNA levels at different post-transplant time points according to candidate SNPs.

Next, we analyzed if there was any association between candidate SNPs and the presence of undectectable plasma TTV DNAemia at baseline (before the initiation of immunosuppressive therapy). Seven (3.7%) patients had pre-transplant TTV DNA levels below the LLoD. We observed that carriers of the minor C allele of the IL10 (rs1878672) SNP in the homozygous state (CC) were more likely to have undetectable baseline TTV DNAemia compared to recipients bearing the reference G allele (GG/GC) [12.5% (3/24) versus 2.5% (4/163), respectively; nominal p-value = 0.046]. There were also significant differences within the TLR3 (rs3775291) SNP, since all the 7 patients with undetectable TTV DNAemia harbored the minor T allele either in the heterozygous or the homozygous state [7.1% (7/99) versus 0.0% (0/88) for CT/TT and CC carriers; nominal p-value = 0.015]. Finally, the minor allele of CD209 (rs4804803) in the homozygous state was also associated with undetectable TTV DNAemia before transplantation [37.5% (5/8) versus 2.2% (4/179) for GG and AA/AG carriers; nominal p-value = 0.0017]. Nevertheless, it should be noted that only the latter association was below the Bonferroni-corrected p-value threshold for statistical significance (which was settled at 0.00178) (Table 4). We further assessed whether the impact of the CD209 (rs4804803) SNP remained significant after adjusting for recipient demographics and pre-transplant clinical characteristics also associated with undetectable TTV DNAemia at baseline (Supplementary Table S2). In a logistic regression model that included recipient age and previous renal replacement therapy as covariates, the presence of the minor G allele of CD209 (rs4804803) in the homozygous state was still significantly associated with pre-transplant TTV DNA levels below the LLoD (adjusted OR: 36.96; 95% CI: 4.72–289.67; p-value = 0.001).

TABLE 4
www.frontiersin.org

TABLE 4. Association between undetectable TTV DNAemia at the baseline (pre-transplant) assessment and candidate SNPs in dominant (heterozygous and homozygous) and recessive (homozygous only) models.

In order to better characterize the genetic determinants of post-transplant TTV viral kinetics, we compared peak TTV DNA levels through different time intervals according to candidate SNPs. The only apparent correlation was observed for TNF (rs1800629), with carriers of the minor allele either in the heterozygous or homozygous state showing lower peak levels during the first 3 post-transplant months (4.2 ± 1.5 versus 5.0 ± 1.8 log10 copies/mL for GA/AA and GG carriers; nominal p-value = 0.008) (Supplementary Table S3). This comparison, however, did not attain the Bonferroni-corrected significance level (settled at 0.00059). Accordingly, recipients bearing the minor A allele of this SNP also showed a non-significant trend—by applying the Bonferroni correction—towards a lower AUC for plasma TTV DNAemia through month 6 (5.8 ± 1.7 versus 6.8 ± 1.7 log10 copies/mL for GA/AA and GG carriers; nominal p-value = 0.007) (Supplementary Table S4). Finally, no significant differences at the Bonferroni-adjusted α level were found in increments (Δ) in TTV DNA levels from baseline to day 15 and months 1, 3, 6 and 12 after transplantation either (Supplementary Table S5).

Discussion

There is increasing evidence on the usefulness of TTV as a surrogate marker of the immune status in a variety of clinical scenarios (Martin-Lopez et al., 2020; Honorato et al., 2021; Studenic et al., 2021), in particular SOT (Fernandez-Ruiz et al., 2019; Redondo et al., 2022a; Eldar-Yedidia et al., 2022; Jaksch et al., 2022), under the rationale that the T-cell-mediated immunity plays an instrumental role in controlling viral replication. Nevertheless, the relative contribution of the innate system—and its genetic determinants—has not been characterized so far. To our knowledge only three previous works have analyzed the impact of genetic polymorphisms on TTV replication in HSCT recipients (with two studies from the same group) and people living with human immunodeficiency virus (HIV) (Prasetyo et al., 2015; Ramzi et al., 2019; Ramzi et al., 2021). Ramzi et al. (2019); Ramzi et al. (2021) found a correlation between SNPs in IL10, CTLA4 and TNF genes and TTV infection in allogeneic HSCT recipients. In detail, the heterozygote genotypes of IL10 rs1800872 (−592C/A) and CTLA4 rs231775 (+49 A/G) were associated with a higher prevalence of TTV DNAemia, whereas the A allele of TNF rs1800629 (−308G/A) had a protective effect. On the other hand, Prasetyo et al. reported a correlation between the APOBEC3B deletion polymorphism status and TTV, HBV and HCV infection among HIV patients (Prasetyo et al., 2015; Ramzi et al., 2019; Ramzi et al., 2021). The present investigation is the first to evaluate to what extent the kinetics of TTV DNA levels following KT are influenced by SNPs in genes coding for PRRs (TLR3, TLR9 and CD209), ILs and cytokines (IL-12B, IL-17, IL-10, TNF), IFN-λ3 (IL-28B) and the costimulatory receptor CTLA-4). These candidate SNPs were chosen on the basis of prior studies performed in the HSCT population (Ramzi et al., 2019; Ramzi et al., 2021) or due to their well-established involvement in other viral infections in SOT recipients (Redondo et al., 2022d).

We found no clear association between any of the SNP genotypes considered and various parameters reflecting viral kinetics after transplantation, such as mean and peak DNA levels or AUCs for plasma TTV DNAemia during discrete periods, or absolute increments from baseline. The significant differences observed for the minor A allele of TNF (rs1800629) in terms of lower peak DNA levels and AUC through months 3 and 6 were not consistent across the entire post-transplant monitoring period and did not survive correction for multiple comparisons. Interestingly, Ramzi et al. (2010) Reported that the A allele of the TNF (rs1800629) SNP was associated with undetectable TTV DNAemia in a single-center cohort of HSCT recipients (OR: 0.46; 95% CI: 0.22–0.96; p-value = 0.025), although the timing for monitoring was unclear and no correction for multiple testing was performed. In line with our results, the same group observed no apparent impact of genotypes of CTLA4 (rs5742909) on the incidence of TTV infection after HSCT (Ramzi et al., 2019; Ramzi et al., 2021).

In addition to the longitudinal post-transplant monitoring of TTV replication, we have specifically investigated the associations between candidate SNPs and the presence of TTV infection at the baseline assessment, before immunosuppressive therapy was initiated. The cross-sectional comparison at this time point would reveal the potential role of genetic predisposition to TTV among ESRD patients in the absence of iatrogenic immunosuppression. In contrast to the negative results observed for the post-transplant period, we found that the minor G allele of CD209 (rs4804803) in the homozygous state exerted a protective effect even after the Bonferroni correction, and that this association with undetectable TTV DNAemia at baseline remained significant after adjusting for clinical covariates. Although caution must be exercised due to the low number of patients with pre-transplant TTV DNA levels below the LLoD, this finding is in accordance with a recent study by our group showing a protective effect against BK polyomavirus viremia linked to the G allele of CD209 (rs4804803) after KT (Redondo et al., 2022c). In addition, the minor G allele has been also associated with a lower susceptibility to tuberculosis (Vannberg et al., 2008) and severe dengue (Sakuntabhai et al., 2005). The CD209 gene codes for DC-SIGN, a transmembrane PRR belonging to the CLR family. It has been described that the presence of the G allele negatively affects gene transcription, thus downregulating the synthesis of DC-SIGN in dendritic cells (Sakuntabhai et al., 2005). In addition, DC-SIGN acts as the cell receptor for many viruses through its high affinity binding of mannose-containing carbohydrates expressed by viral glycoproteins (Lin et al., 2003). In view of the non-enveloped structure of anelloviruses, the mechanistic explanation for the association found between the CD209 (rs4804803) SNP and baseline TTV DNAemia remains to be determined and demands further investigation in healthy subjects (e.g., blood donors).

Despite its large sample size, frequent TTV DNA monitoring and comprehensive set of SNPs screened, some limitations to our study should be acknowledged. As previously described (Redondo et al., 2022a; Jaksch et al., 2022), the vast majority of recipients had TTV replication early after transplantation. Thus, associations between genetic polymorphisms and undetectable TTV DNAemia (below the LLoD of the PCR assay) were only analyzed at baseline. Since the relatively high number of SNPs imposed stringent thresholds for statistical significance, false-negative results due to insufficient statistical power cannot be excluded, particularly for those SNPs—such as CTLA4 (rs5742909) or TNF (rs1800629)—with very low absolute numbers of patients bearing the corresponding minor alleles.

In conclusion, the G allele of CD209 (rs4804803) in the homozygous state would play a protective role against TTV in non-immunocompromised patients listed for KT, whereas no significant associations have been found during the post-transplant period for any of the studied SNPs. Thus, the present results support the conception that variations in plasma TTV DNA levels after KT are mainly driven by the effect of immunosuppressive therapy rather than by underlying genetic predisposition, reinforcing its clinical usefulness as a surrogate marker of immunosuppression. Post-transplant TTV replication kinetics seems to be mainly under the control of the adaptive immune responses, with no meaningful effect of SNPs in genes orchestrating innate arm.

Data availability statement

The data that support the findings of this study are available from the corresponding author upon reasonable request. SNP genotyping data are registered in the BioProject database under the ID PRJNA898147.

Ethics statement

The studies involving human participants were reviewed and approved by the Clinical Research Ethics Committee Hospital 12 de Octubre (Study protocol number 14/030). The patients/participants provided their written informed consent to participate in this study.

Author contributions

NR and MF-R designed the study, performed statistical analyses and wrote the manuscript; EA, PP, and EG performed laboratory analyses; TR-M collected patient samples; IR-G, FL-M, RS, EG, NP, and AA participated in patient recruitment and performed data collection; AA, DN, and JMA critically reviewed the manuscript and provided significant input and feedback. All authors read and approved the final manuscript.

Funding

This study has been funded by Instituto de Salud Carlos III (ISCIII), Spanish Ministry of Science and Innovation, through the projects PIE13/00045, PI15/01953, and PI19/01300—co-funded by European Regional Development Fund/European Social Fund “A way to make Europe”/“Investing in your future”. IR-G holds a research training contract “Río Hortega” (CM19/00163) and MF-R holds a research contract “Miguel Servet” (CP18/00073), both from the ISCIII and also co-funded by the European Union.

Acknowledgments

The authors gratefully acknowledge all patients recruited in the institutional cohort of kidney transplant recipients for their participation.

Conflict of interest

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.

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/fgene.2022.1069890/full#supplementary-material

Abbreviations

AUC, area the under curve; BKPyV, BK polyomavirus; CI, confidence interval; CpG, cytosine-phosphate-guanine; CMV, cytomegalovirus; D, donor; dsRNA, double-stranded RNA; ESRD, end-stage renal disease; HBV, hepatitis B virus; HCV, hepatitis C virus; HIV, human immunodeficiency virus; HSCT, hematopoietic stem cell transplantation; IQR, interquartile range; KT, kidney transplantation; LLoD, lower limit of detection; OR, odds ratio; PAMP, pathogen-associated molecular pattern; PCR, polymerase chain reaction; R, recipient; PRR, pattern recognition receptor; SD, standard deviation; SNP, single-nucleotide polymorphism; SOT, solid organ transplantation; TLR, toll-like receptor; TTV, Torque teno virus.

References

Alagarasu, K., Kaushal, H., Shinde, P., Kakade, M., Chaudhary, U., Padbidri, V., et al. (2021). TNFA and IL10 polymorphisms and IL-6 and IL-10 levels influence disease severity in influenza A(H1N1)pdm09 virus infected patients. Genes. 12 (12), 1914. doi:10.3390/genes12121914

PubMed Abstract | CrossRef Full Text | Google Scholar

Ali, E. S. G., Bassyouni, R. H., Abdelaleem, O. O., Hassan, E. A., and Gaber, S. N. (2022). Association between SNPs of Cytotoxic T-lymphocyte-associated protein 4 (CTLA-4), programmed cell death 1 (PD-1)and the susceptibility to chronic Hepatitis C infection in virus C-infected patients. Virus Res. 310, 198684. doi:10.1016/j.virusres.2022.198684

PubMed Abstract | CrossRef Full Text | Google Scholar

Arze, C. A., Springer, S., Dudas, G., Patel, S., Bhattacharyya, A., Swaminathan, H., et al. (2021). Global genome analysis reveals a vast and dynamic anellovirus landscape within the human virome. Cell. Host Microbe 29 (8), 1305–1315.e6. doi:10.1016/j.chom.2021.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Biagini, P. (2009). Classification of TTV and related viruses (anelloviruses). Curr. Top. Microbiol. Immunol. 331, 21–33. doi:10.1007/978-3-540-70972-5_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bravo, D., Solano, C., Gimenez, E., Remigia, M. J., Corrales, I., Amat, P., et al. (2014). Effect of the IL28B Rs12979860 C/T polymorphism on the incidence and features of active cytomegalovirus infection in allogeneic stem cell transplant patients. J. Med. Virol. 86 (5), 838–844. doi:10.1002/jmv.23865

PubMed Abstract | CrossRef Full Text | Google Scholar

Bucardo, F., Reyes, Y., Morales, M., Briceno, R., Gonzalez, F., Lundkvist, A., et al. (2021). Association of genetic polymorphisms in DC-SIGN, toll-like receptor 3, and tumor necrosis factor alpha genes and the lewis-negative phenotype with chikungunya infection and disease in Nicaragua. J. Infect. Dis. 223 (2), 278–286. doi:10.1093/infdis/jiaa364

PubMed Abstract | CrossRef Full Text | Google Scholar

Chaaithanya, I. K., Muruganandam, N., Surya, P., Anwesh, M., Alagarasu, K., and Vijayachari, P. (2016). Association of oligoadenylate synthetase gene cluster and DC-SIGN (CD209) gene polymorphisms with clinical symptoms in chikungunya virus infection. DNA Cell. Biol. 35 (1), 44–50. doi:10.1089/dna.2015.2819

PubMed Abstract | CrossRef Full Text | Google Scholar

Chihab, H., Zaidane, I., Elhabazi, A., Jadid, F. Z., El Fihri, R., Elmessaoudi-Idrissi, M., et al. (2019). Toll-like receptor 9 polymorphisms and Hepatitis B virus clearance in Moroccan chronic carriers. Gene 687, 212–218. doi:10.1016/j.gene.2018.11.041

PubMed Abstract | CrossRef Full Text | Google Scholar

Corrales, I., Gimenez, E., Solano, C., Amat, P., de la Camara, R., Nieto, J., et al. (2015). Incidence and dynamics of active cytomegalovirus infection in allogeneic stem cell transplant patients according to single nucleotide polymorphisms in donor and recipient CCR5, MCP-1, IL-10, and TLR9 genes. J. Med. Virol. 87 (2), 248–255. doi:10.1002/jmv.24050

PubMed Abstract | CrossRef Full Text | Google Scholar

Czupryna, P., Parczewski, M., Grygorczuk, S., Pancewicz, S., Zajkowska, J., Dunaj, J., et al. (2017). Analysis of the relationship between single nucleotide polymorphism of the CD209, IL-10, IL-28 and CCR5 D32 genes with the human predisposition to developing tick-borne encephalitis. Postepy Hig. Med. Dosw. 71 (1), 788–796. doi:10.5604/01.3001.0010.3856

PubMed Abstract | CrossRef Full Text | Google Scholar

Dodi, G., Attanasi, M., Di Filippo, P., Di Pillo, S., and Chiarelli, F. (2021). Virome in the lungs: The role of anelloviruses in childhood respiratory diseases. Microorganisms 9 (7), 1357. doi:10.3390/microorganisms9071357

PubMed Abstract | CrossRef Full Text | Google Scholar

Egli, A., Levin, A., Santer, D. M., Joyce, M., O'Shea, D., Thomas, B. S., et al. (2014). Immunomodulatory Function of Interleukin 28B during primary infection with cytomegalovirus. J. Infect. Dis. 210 (5), 717–727. doi:10.1093/infdis/jiu144

PubMed Abstract | CrossRef Full Text | Google Scholar

Eldar-Yedidia, Y., Ben-Shalom, E., Hillel, M., Belostotsky, R., Megged, O., Freier-Dror, Y., et al. (2022). Association of post-transplantation anellovirus viral load with kidney transplant rejection in children. Pediatr. Nephrol. 37, 1905–1914. doi:10.1007/s00467-021-05336-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Eskandari-Nasab, E., Moghadampour, M., Tahmasebi, A., and Asadi-Saghandi, A. (2018). Interleukin-17 A and F gene polymorphisms affect the risk of tuberculosis: An updated meta-analysis. Indian J. Tuberc. 65 (3), 200–207. doi:10.1016/j.ijtb.2017.08.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandez-Ruiz, M., Albert, E., Gimenez, E., Ruiz-Merlo, T., Parra, P., Lopez-Medrano, F., et al. (2019). Monitoring of alphatorquevirus DNA levels for the prediction of immunosuppression-related complications after kidney transplantation. Am. J. Transpl. 19 (4), 1139–1149. doi:10.1111/ajt.15145

CrossRef Full Text | Google Scholar

Fernandez-Ruiz, M., Corrales, I., Arias, M., Campistol, J. M., Gimenez, E., Crespo, J., et al. (2015). Association between individual and combined SNPs in genes related to innate immunity and incidence of CMV infection in seropositive kidney transplant recipients. Am. J. Transpl. 15 (5), 1323–1335. doi:10.1111/ajt.13107

CrossRef Full Text | Google Scholar

Fischer, J., Koukoulioti, E., Schott, E., Fulop, B., Heyne, R., Berg, T., et al. (2018). Polymorphisms in the Toll-like receptor 3 (TLR3) gene are associated with the natural course of Hepatitis B virus infection in Caucasian population. Sci. Rep. 8 (1), 12737. doi:10.1038/s41598-018-31065-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Focosi, D., Antonelli, G., Pistello, M., and Maggi, F. (2016). Torquetenovirus: The human virome from bench to bedside. Clin. Microbiol. Infect. 22 (7), 589–593. doi:10.1016/j.cmi.2016.04.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Ge, D., Fellay, J., Thompson, A. J., Simon, J. S., Shianna, K. V., Urban, T. J., et al. (2009). Genetic variation in IL28B predicts hepatitis C treatment-induced viral clearance. Nature 461 (7262), 399–401. doi:10.1038/nature08309

PubMed Abstract | CrossRef Full Text | Google Scholar

Geng, P. L., Song, L. X., An, H., Huang, J. Y., Li, S., and Zeng, X. T. (2016). Toll-like receptor 3 is associated with the risk of HCV infection and HBV-related diseases. Med. Baltim. 95 (21), e2302. doi:10.1097/MD.0000000000002302

PubMed Abstract | CrossRef Full Text | Google Scholar

Goncalves de Albuquerque, S. D. C., da Costa Oliveira, C. N., Vaitkevicius-Antao, V., Silva, A. C., Luna, C. F., de Lorena, V. M. B., et al. (2019). Study of association of the rs2275913 IL-17A single nucleotide polymorphism and susceptibility to cutaneous leishmaniasis caused by Leishmania braziliensis. Cytokine 123, 154784. doi:10.1016/j.cyto.2019.154784

PubMed Abstract | CrossRef Full Text | Google Scholar

Gorzer, I., Haloschan, M., Jaksch, P., Klepetko, W., and Puchhammer-Stockl, E. (2014). Plasma DNA levels of Torque teno virus and immunosuppression after lung transplantation. J. Heart Lung Transpl. 33 (3), 320–323. doi:10.1016/j.healun.2013.12.007

CrossRef Full Text | Google Scholar

Gupta, K., Kaur, G., Pathak, T., and Banerjee, I. (2022). Systematic review and meta-analysis of human genetic variants contributing to COVID-19 susceptibility and severity. Gene 844, 146790. doi:10.1016/j.gene.2022.146790

PubMed Abstract | CrossRef Full Text | Google Scholar

Honorato, L., Witkin, S. S., Mendes-Correa, M. C., Conde Toscano, A. L. C., Linhares, I. M., de Paula, A. V., et al. (2021). The torque teno virus titer in saliva reflects the level of circulating CD4(+) T lymphocytes and HIV in individuals undergoing antiretroviral maintenance therapy. Front. Med. 8, 809312. doi:10.3389/fmed.2021.809312

PubMed Abstract | CrossRef Full Text | Google Scholar

Jablonska, A., Studzinska, M., Szenborn, L., Wisniewska-Ligier, M., Karlikowska-Skwarnik, M., Gesicki, T., et al. (2020). TLR4 896A/G and TLR9 1174G/A polymorphisms are associated with the risk of infectious mononucleosis. Sci. Rep. 10 (1), 13154. doi:10.1038/s41598-020-70129-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaksch, P., Gorzer, I., Puchhammer-Stockl, E., and Bond, G. (2022). Integrated immunologic monitoring in solid organ transplantation: The road towards torque teno virus-guided immunosuppression. Transplantation 106 (10), 1940–1951. doi:10.1097/TP.0000000000004153

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaksch, P., Kundi, M., Gorzer, I., Murakozy, G., Lambers, C., Benazzo, A., et al. (2018). Torque teno virus as a novel biomarker targeting the efficacy of immunosuppression after lung transplantation. J. Infect. Dis. 218 (12), 1922–1928. doi:10.1093/infdis/jiy452

PubMed Abstract | CrossRef Full Text | Google Scholar

Jedlinska-Pijanowska, D., Kasztelewicz, B., Dobrzanska, A., Dzierzanowska-Fangrat, K., Jaworski, M., and Czech-Kowalska, J. (2021). Association between single nucleotide polymorphisms and viral load in congenital cytomegalovirus infection. J. Mother Child. 24 (4), 9–17. doi:10.34763/jmotherandchild.20202404.d-20-00014

PubMed Abstract | CrossRef Full Text | Google Scholar

Kaczorowska, J., and van der Hoek, L. (2020). Human anelloviruses: Diverse, omnipresent and commensal members of the virome. FEMS Microbiol. Rev. 44 (3), 305–313. doi:10.1093/femsre/fuaa007

PubMed Abstract | CrossRef Full Text | Google Scholar

Kindberg, E., Vene, S., Mickiene, A., Lundkvist, A., Lindquist, L., and Svensson, L. (2011). A functional Toll-like receptor 3 gene (TLR3) may be a risk factor for tick-borne encephalitis virus (TBEV) infection. J. Infect. Dis. 203 (4), 523–528. doi:10.1093/infdis/jiq082

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y., Yang, L., Sha, K., Liu, T., and Zhang, L. (2016). Correlation of interferon-lambda 4 ss469415590 with the hepatitis C virus treatment response and its comparison with interleukin 28b polymorphisms in predicting a sustained virological response: A meta-analysis. Int. J. Infect. Dis. 53, 52–58. doi:10.1016/j.ijid.2016.10.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, G., Simmons, G., Pohlmann, S., Baribaud, F., Ni, H., Leslie, G. J., et al. (2003). Differential N-linked glycosylation of human immunodeficiency virus and Ebola virus envelope glycoproteins modulates interactions with DC-SIGN and DC-SIGNR. J. Virol. 77 (2), 1337–1346. doi:10.1128/jvi.77.2.1337-1346.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin-Lopez, M., Albert, E., Fernandez-Ruiz, M., Gonzalez-Alvaro, I., Rodriguez, E., Aguado, J. M., et al. (2020). Torque teno virus viremia in patients with chronic arthritis: Influence of biologic therapies. Semin. Arthritis Rheum. 50 (1), 166–171. doi:10.1016/j.semarthrit.2019.06.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Mezger, M., Steffens, M., Semmler, C., Arlt, E. M., Zimmer, M., Kristjanson, G. I., et al. (2008). Investigation of promoter variations in dendritic cell-specific ICAM3-grabbing non-integrin (DC-SIGN) (CD209) and their relevance for human cytomegalovirus reactivation and disease after allogeneic stem-cell transplantation. Clin. Microbiol. Infect. 14 (3), 228–234. doi:10.1111/j.1469-0691.2007.01902.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Miri, H. H., Fazeli, P., Ali-Hassanzadeh, M., Bemani, P., Kabelitz, D., and Kalantar, K. (2021). Correlation between IL-28 polymorphism and spontaneous clearance in HCV patients: Systematic review and meta-analysis. Arch. Virol. 166 (9), 2469–2478. doi:10.1007/s00705-021-05141-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Misra, M. K., Pandey, S. K., Kapoor, R., Sharma, R. K., and Agrawal, S. (2015). Cytotoxic T-lymphocyte antigen 4 gene polymorphism influences the incidence of symptomatic human cytomegalovirus infection after renal transplantation. Pharmacogenet. Genomics 25 (1), 19–29. doi:10.1097/FPC.0000000000000102

PubMed Abstract | CrossRef Full Text | Google Scholar

Mouton, W., Conrad, A., Bal, A., Boccard, M., Malcus, C., Ducastelle-Lepretre, S., et al. (2020). Torque teno virus viral load as a marker of immune function in allogeneic haematopoietic stem cell transplantation recipients. Viruses 12 (11), E1292. doi:10.3390/v12111292

PubMed Abstract | CrossRef Full Text | Google Scholar

Nishizawa, T., Okamoto, H., Konishi, K., Yoshizawa, H., Miyakawa, Y., and Mayumi, M. (1997). A novel DNA virus (TTV) associated with elevated transaminase levels in posttransfusion hepatitis of unknown etiology. Biochem. Biophys. Res. Commun. 241 (1), 92–97. doi:10.1006/bbrc.1997.7765

PubMed Abstract | CrossRef Full Text | Google Scholar

Prasetyo, A. A., Sariyatun, R., Reviono,, , Sari, Y., Hudiyono,, , Haryati, S., et al. (2015). The APOBEC3B deletion polymorphism is associated with prevalence of Hepatitis B virus, hepatitis C virus, Torque Teno virus, and Toxoplasma gondii co-infection among HIV-infected individuals. J. Clin. Virol. 70, 67–71. doi:10.1016/j.jcv.2015.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramzi, M., Arandi, N., Zarei, T., Saadi, M. I., Yaghobi, R., Moghadam, M., et al. (2019). Genetic variation of TNF-alpha and IL-10, IL-12, IL-17 genes and association with torque teno virus infection post hematopoietic stem cell transplantation. Acta Virol. 63 (2), 186–194. doi:10.4149/av_2019_210

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramzi, M., Iravani Saadi, M., Zarei, T., Yaghobi, R., and Arandi, N. (2021). Association between cytotoxic T-lymphocyte antigen 4 gene polymorphisms and torque teno virus infection after hematopoietic stem cell transplantation. Exp. Clin. Transpl. 19 (3), 259–263. doi:10.6002/ect.2017.0105

PubMed Abstract | CrossRef Full Text | Google Scholar

Redondo, N., Navarro, D., Aguado, J. M., and Fernandez-Ruiz, M. (2022). Human genetic polymorphisms and risk of viral infection after solid organ transplantation. Transpl. Rev. 36 (1), 100669. doi:10.1016/j.trre.2021.100669

PubMed Abstract | CrossRef Full Text | Google Scholar

Redondo, N., Navarro, D., Aguado, J. M., and Fernandez-Ruiz, M. (2022). Viruses, friends, and foes: The case of Torque Teno Virus and the net state of immunosuppression. Transpl. Infect. Dis. 24 (2), e13778. doi:10.1111/tid.13778

PubMed Abstract | CrossRef Full Text | Google Scholar

Redondo, N., Rodriguez-Goncer, I., Parra, P., Lopez-Medrano, F., Gonzalez, E., Hernandez, A., et al. (2022). Genetic polymorphisms in TLR3, IL10 and CD209 influence the risk of BK polyomavirus infection after kidney transplantation. Sci. Rep. 12 (1), 11338. doi:10.1038/s41598-022-15406-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Redondo, N., Rodriguez-Goncer, I., Parra, P., Ruiz-Merlo, T., Lopez-Medrano, F., Gonzalez, E., et al. (2022). Influence of single-nucleotide polymorphisms in TLR3 (rs3775291) and TLR9 (rs352139) on the risk of CMV infection in kidney transplant recipients. Front. Immunol. 13, 929995. doi:10.3389/fimmu.2022.929995

PubMed Abstract | CrossRef Full Text | Google Scholar

Rezahosseini, O., Drabe, C. H., Sorensen, S. S., Rasmussen, A., Perch, M., Ostrowski, S. R., et al. (2019). Torque-Teno virus viral load as a potential endogenous marker of immune function in solid organ transplantation. Transpl. Rev. 33 (3), 137–144. doi:10.1016/j.trre.2019.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Sakuntabhai, A., Turbpaiboon, C., Casademont, I., Chuansumrit, A., Lowhnoo, T., Kajaste-Rudnitski, A., et al. (2005). A variant in the CD209 promoter is associated with severity of dengue disease. Nat. Genet. 37 (5), 507–513. doi:10.1038/ng1550

PubMed Abstract | CrossRef Full Text | Google Scholar

Singh, A. K., Prakash, S., Garg, R. K., Jain, P., Kumar, R., and Jain, A. (2021). Study of single nucleotide polymorphisms in endosomal toll-like receptors-3, 7, and 9 genes in patients with dengue: A case-control study. Cureus 13 (5), e14883. doi:10.7759/cureus.14883

PubMed Abstract | CrossRef Full Text | Google Scholar

Studenic, P., Bond, G., Kerschbaumer, A., Becede, M., Pavelka, K., Karateev, D., et al. (2021). Torque Teno Virus quantification for monitoring of immunomodulation with biologic compounds in the treatment of rheumatoid arthritis. Rheumatol. Oxf. 61, 2815–2825. doi:10.1093/rheumatology/keab839

PubMed Abstract | CrossRef Full Text | Google Scholar

Studzinska, M., Jablonska, A., Wisniewska-Ligier, M., Nowakowska, D., Gaj, Z., Lesnikowski, Z. J., et al. (2017). Association of TLR3 L412F polymorphism with cytomegalovirus infection in children. PLoS One 12 (1), e0169420. doi:10.1371/journal.pone.0169420

PubMed Abstract | CrossRef Full Text | Google Scholar

Takeuchi, O., and Akira, S. (2010). Pattern recognition receptors and inflammation. Cell. 140 (6), 805–820. doi:10.1016/j.cell.2010.01.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Thomas, D. L., Thio, C. L., Martin, M. P., Qi, Y., Ge, D., O'Huigin, C., et al. (2009). Genetic variation in IL28B and spontaneous clearance of hepatitis C virus. Nature 461 (7265), 798–801. doi:10.1038/nature08463

PubMed Abstract | CrossRef Full Text | Google Scholar

Vannberg, F. O., Chapman, S. J., Khor, C. C., Tosh, K., Floyd, S., Jackson-Sillah, D., et al. (2008). CD209 genetic polymorphism and tuberculosis disease. PLoS One 3 (1), e1388. doi:10.1371/journal.pone.0001388

PubMed Abstract | CrossRef Full Text | Google Scholar

Vargas-Castillo, A. B., Ruiz-Tovar, K., Vivanco-Cid, H., Quiroz-Cruz, S., Escobar-Gutierrez, A., Cerna-Cortes, J. F., et al. (2018). Association of single-nucleotide polymorphisms in immune-related genes with development of dengue hemorrhagic fever in a Mexican population. Viral Immunol. 31 (3), 249–255. doi:10.1089/vim.2017.0069

PubMed Abstract | CrossRef Full Text | Google Scholar

Varshney, D., Singh, S., Sinha, E., Mohanty, K. K., Kumar, S., Kumar Barik, S., et al. (2022). Systematic review and meta-analysis of human Toll-like receptors genetic polymorphisms for susceptibility to tuberculosis infection. Cytokine 152, 155791. doi:10.1016/j.cyto.2021.155791

PubMed Abstract | CrossRef Full Text | Google Scholar

Webb, B., Rakibuzzaman, A., and Ramamoorthy, S. (2020). Torque teno viruses in health and disease. Virus Res. 285, 198013. doi:10.1016/j.virusres.2020.198013

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, S., Zhang, X., Zhang, Y. B., Tian, X., Liu, A., Cui, C., et al. (2020). Association of TLR3 (rs3775291) and IL-10 (rs1800871) gene polymorphisms with susceptibility to hepatitis B infection: A meta-analysis. Epidemiol. Infect. 148, e228. doi:10.1017/S0950268820002101

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: single-nucleotide polymorphisms, Torque teno virus, TTV replication kinetics, kidney transplantation, SNP

Citation: Redondo N, Rodríguez-Goncer I, Parra P, Albert E, Giménez E, Ruiz-Merlo T, López-Medrano F, San Juan R, González E, Sevillano Á, Andrés A, Navarro D, Aguado JM and Fernández-Ruiz M (2022) Impact of polymorphisms in genes orchestrating innate immune responses on replication kinetics of Torque teno virus after kidney transplantation. Front. Genet. 13:1069890. doi: 10.3389/fgene.2022.1069890

Received: 14 October 2022; Accepted: 09 November 2022;
Published: 22 November 2022.

Edited by:

Katarzyna Bogunia-Kubik, Polish Academy of Sciences, Poland

Reviewed by:

Daniele Focosi, Pisana University Hospital, Italy
Valéria de Lima Kaminski, Federal University of São Paulo, Brazil

Copyright © 2022 Redondo, Rodríguez-Goncer, Parra, Albert, Giménez, Ruiz-Merlo, López-Medrano, San Juan, González, Sevillano, Andrés, Navarro, Aguado and Fernández-Ruiz. 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: Natalia Redondo, bmF0YWxpYS5yZWRvbmRvLmltYXMxMkBoMTJvLmVz

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.