Next Article in Journal
Disease Course of Korean African Swine Fever Virus in Domestic Pigs Exposed Intraorally, Intranasally, Intramuscularly, and by Direct Contact with Infected Pigs
Next Article in Special Issue
Torque Teno Virus DNA Load in Blood as an Immune Status Biomarker in Adult Hematological Patients: The State of the Art and Future Prospects
Previous Article in Journal
Development of Glycyrrhizinic Acid-Based Lipid Nanoparticle (LNP-GA) as An Adjuvant That Improves the Immune Response to Porcine Epidemic Diarrhea Virus Spike Recombinant Protein
Previous Article in Special Issue
High-Density Lipoprotein Particles and Torque Teno Virus in Stable Outpatient Kidney Transplant Recipients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Torque Teno Virus (TTV) in Renal Transplant Recipients: Species Diversity and Variability

1
Virology Unit, Centro de Educación Médica e Investigaciones Clínicas (CEMIC) University Hospital, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Av. Galván 4102, Buenos Aires C1631FWO, Argentina
2
Laboratory of Virology, National Institute for Infectious Diseases Lazzaro Spallanzani—IRCCS, 00149 Rome, Italy
3
Department of Translational Research, University of Pisa, 56127 Pisa, Italy
4
Department of Nephrology, Centro de Educación Médica e Investigaciones Clínicas (CEMIC) University Hospital, Buenos Aires C1631FWO, Argentina
5
Genetic Epidemiology Laboratory, Centro de Educación Médica e Investigaciones Clínicas (CEMIC) University Hospital, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires C1631FWO, Argentina
*
Authors to whom correspondence should be addressed.
Viruses 2024, 16(3), 432; https://doi.org/10.3390/v16030432
Submission received: 24 January 2024 / Revised: 26 February 2024 / Accepted: 28 February 2024 / Published: 11 March 2024
(This article belongs to the Special Issue Advancing Research of Anelloviruses)

Abstract

:
Torque Teno Virus (TTV) is a nonpathogenic and ubiquitous ssDNA virus, a member of the Anelloviridae family. TTV has been postulated as a biomarker in transplant patients. This study aimed to determine the TTV species diversity and variability in renal transplant recipients and to associate species diversity with the corresponding TTV viral load. From 27 recipients, 30 plasma samples were selected. Viral load was determined using two real-time PCR assays, followed by RCA-NGS and ORF1 phylogenetic analysis. The TTV diversity was determined in all samples. Variability was determined in three patients with two sequential samples (pre- and post-transplantation). Most of the samples presented multiple TTV species, up to 15 different species were detected. In the pre-transplant samples (n = 12), the most prevalent species were TTV3 (75%) and TTV13 (75%), and the median number of species per sample was 5 (IQR: 4–7.5). TTV3 was also the most prevalent (56%) in the post-transplant samples (n = 18), and the median number of species was 2 (IQR: 1.8–5.5). No significant correlation between the number of species and viral load was found. The number and type of TTV species showed total variability over time. We report high TTV species diversity in Argentinian recipients, especially in pre-transplant period, with total intra-host variability. However, we found no significant correlation between this high diversity and TTV viral load.

1. Introduction

Torque Teno Virus (TTV) is a non-enveloped virus, with a small, single-stranded, negative-sense DNA genome. This virus is nonpathogenic in humans and highly ubiquitous [1]. TTV has been postulated as a surrogate marker of immunosuppression in transplant patients, especially among solid organ recipients [2].
TTV genome is approximately 3.8 Kb in length, organized into a coding region sequence (CDS) of around 2.6 Kb and a highly conserved non-coding region (UTR) of approximately 1.2 kb [1,3,4]. The CDS presents at least four partially overlapping open reading frames (ORF1–4) [5,6]. Like most circular DNA viruses, TTV DNA replication can be carried out by a rolling circle mechanism [7].
The TT designation of this virus was in honor of a Japanese patient with transfusion-acquired non-AG hepatitis, from whom the virus was first isolated in 1997 [8]. TTV is a member of the highly diverse Anelloviridae family. Three members of this family can infect humans: TTV (genus Alphatorquevirus), Torque Teno Mini Virus (TTMV, genus Betatorquevirus) and Torque Teno Midi Virus (TTMDV, genus Gammatorquevirus) [9,10].
According to the previous International Committee on Taxonomy of Viruses (ICTV) guidelines, the 29 TTV species described were divided into five genogroups: genogroup I (TTV 1–5), genogroup II (TTV 6–8), genogroup III (TTV 13–24), genogroup IV (TTV 25–29) and genogroup V (TTV 9–12). In the latest ICTV update, TTV species were re-categorized into 22 species, based on a sequence identity threshold of 69% for the ORF1 coding region [11].
TTV can be detected in peripheral blood regardless of age, sex or socioeconomic level, with 90% prevalence in healthy individuals. TTV transmission occurs mainly through the fecal–oral and respiratory routes. Transmission may also occur in utero, through the umbilical cord or amniotic fluid [12,13,14,15].
Several studies on blood donors demonstrated that TTV prevalence is linearly correlated with age. In a Romanian study including 701 blood donors with a median age of 38 years old, the TTV prevalence was 58%. In the largest cohort study published to date, including 1017 Italian blood donors with a median age of 44 years old, the detected TTV prevalence was 65%. A higher TTV prevalence (76%) was detected in a cohort of 313 Austrian blood donors with a median age of 53 years [16,17,18].
Despite being a non-pathogenic virus, currently not associated with any known pathologies, TTV appears to establish persistent infections in its host. Furthermore, it is postulated that TTV inhibits part of the host immune system [5,19].
TTV ORF2 encodes a protein (~200 aa) with a tyrosine phosphatase motif, capable of interfering with the translocation of NF-κB. NF-κB is a well-characterized signal transcription factor, crucial for lymphocyte development, inflammation and cell growth, among other processes. Therefore, by inhibiting the expression of these pro-inflammatory genes, TTV can evade the host’s immune system [20,21].
Another mechanism employed by TTV to enhance evasion and facilitate persistence in the host is the expression of microRNAs (miRNA, small, non-coding, 22-nucleotide-long RNAs). Potential miRNA sequences in the TTV genome have been discovered using computational methods (TTV miRNA t1a of genogroup 1, TTV miRNA t3b of genogroup 3 and TTV miRNA tth8 of genogroup 5). This miRNA-mediated immune evasion may contribute to the ubiquity and persistence of TTV in the population. The role and importance of TTV miRNA expression have not yet been fully studied [22].
Multiple TTV species can simultaneously be detected in humans, which may be acquired at a single time or sequentially over time, resulting in a mixture of TTV species and even different members of the Anelloviridae family. Only a few studies have been conducted in Latin America, reporting the most prevalent TTV species were from genogroup III, followed by species from genogroup I. This was also observed in Europe in healthy populations [23,24,25].
Advances in next-generation sequencing (NGS) have deepened our knowledge on the microbial world. The combination of molecular biology technics like rolling circle amplification (RCA) with NGS on different human samples (both fluid and tissues) demonstrated that anelloviruses are one of the major components of the human virome, with TTV being one of the most prevalent viruses detected [24,26,27].
The aims of this study were to determine the TTV species diversity and variability in renal transplant patients over time and to associate TTV species diversity with the corresponding TTV viral load.

2. Materials and Methods

Plasma samples from adult recipients (>18 years old) who underwent renal transplantation between November 2018 and April 2021 at CEMIC University Hospital (Buenos Aires, Argentina) were prospectively collected and evaluated for TTV viral load before and during the first year after transplantation.
The inclusion criteria included renal transplant patients who received thymoglobulin or basiliximab for induction therapy, followed by an immunosuppression maintenance scheme with tacrolimus, mycophenolic acid and steroids.
The exclusion criteria included renal transplant patients with primary renal graft failure or changes in the immunosuppression maintenance scheme. Patients who changed their ambulatory care center before the third month post-transplantation were also excluded.
Diversity was defined as the presence of more than one TTV species in the same sample, and variability was defined as the presence of different species (variation) over time within the same host.
To rule out graft rejection, all patients were evaluated using protocol renal biopsies between 3–6 months and one year after transplantation. Additional renal biopsies were performed in case the patient presented renal dysfunction or to confirm BK nephropathy. Cellular and humoral graft rejections were documented and scored according to the BANFF classification (2019 update). Borderline changes were also included as a graft rejection event since all cases of allograft rejection (cellular, humoral and borderline changes) led to clinical intervention with antirejection treatment.
All recipients were regularly screened for Cytomegalovirus (CMV), both after the completion of routine valganciclovir prophylaxis or if the patient, despite prophylaxis, had symptoms consistent with CMV infection. BK polyomaviruses (BKPyV) were tested by PCR every 2 months for the first six months after renal transplantation and every 3 months thereafter.

2.1. Plasma Samples

A total of 30 samples from 27 patients were selected for evaluating the TTV species diversity using next-generation sequencing (NGS).
The selection criteria were a TTV viral load higher than 3 Log10 copies/mL or viral load discrepancy between two PCR assays. Twelve were pre-transplant samples and eighteen were post-transplant samples. Of these 18 samples, 9 showed a discrepant result between the two detection methods used.
To evaluate the TTV species variability, 3 patients were evaluated before and after transplantation.

2.2. TTV Detection and Quantification

TTV was detected and quantified using two real-time PCR assays:
TTV Home Brew: A home brew PCR method was performed, using primers from a highly conserved UTR region, as previously reported [28]. Positive and negative controls were included in each run. Quantification was performed using a quantified plasmid.
Commercial PCR: The TTV R-GENE® kit (bioMérieux, Marcy-l’Etoile, France). TTV quantification was performed according to the manufacturer’s instructions. The sensitivity control, negative control, positive control and four quantification standards were included in each run. The reported limit of detection is 2.4 Log10 copies/mL.
TTV R-GENE® PCR can occasionally detect high levels of TTMV and TTMDV, other members of the Anelloviridae family.

2.3. TTV Genotyping

A whole-genome, amplicon-based, next-generation sequencing (NGS) technique was used to detect individual TTV species, as previously described [29].
Briefly, rolling circle amplification (RCA) technology was used to enrich the circular DNA of the anelloviruses. The RCA reaction mixture contained a DNA sample of less than 10 ng, 25 µM of exonuclease-resistant random primers (Thermo Fisher Scientific, Carlsbad, CA, USA), 4 mM of deoxynucleotides (dNTPs) (Solis BioDyne, Tartu, Estonia) and 10 U of φ29 DNA polymerase (Thermo Fisher Scientific, Carlsbad, CA, USA). Amplification was carried out at 30 °C for 18 h, followed by inactivation of the φ29 DNA polymerase at 65 °C for 10 min. The resulting product was quantified using a NanoDrop Lite instrument (Thermo Fisher Scientific) and diluted for use as a template in the universal anellovirus inverse (UAvI) PCR.
The UAvI PCR was then performed in order to obtain a full-genome amplicon. The reaction mixture of 50 μL contained 1× PrimeSTAR® GXL Buffer (Takara, Otsu, Japan), a 200 μM dNTP mixture (Takara, Otsu, Japan), 0.3 μM forward primer, 0.3 μM reverse primer, 1.25 U of PrimeSTAR® GXL DNA polymerase (Takara, Otsu, Japan) and 500 to 1000 ng of template. The primers used were 206INV_TTVFor and 205INV_TTVRev. The PCR was carried out for 35 cycles. The PCR product was purified and quantified. Libraries were generated using the Illumina DNA Prep kit, with an initial DNA input of 150 ng. The pool of libraries was quantified using Qubit, and the sizes of the fragments were estimated using the Bioanalyzer. The denatured libraries were then sequenced using an Illumina MiSeq V3 cartridge (600 cycles) with a paired-end 2 × 250 bp and 2 × 150 bp layout.
In the bioinformatic analysis, a pipeline that integrated both reference-based and de novo assembly approaches was used. Firstly, low-quality reads and adapters were trimmed using the fastp software v0.220. Human-origin reads were subsequently removed with the use of the Kraken 2 softwarev2.1.2, using a human genome database (Accession: GRCh38.p13).
The remaining reads were aligned against a comprehensive anelloviruses database. Mapped reads were used to reconstruct complete (or nearly complete) genomic consensus sequences. Non-human reads were then aligned using the BWA-MEM software v.0.7.17 to cross-validate the viral assembly. The obtained ORF1 gene sequences using the orfipy algorithm were then subjected to visual inspection for quality assurance.

2.4. Phylogenetic Analysis

The ORF1 phylogenetic analysis was performed using the sequences obtained from the selected plasma samples plus reference sequences (n = 190) from GenBank. These reference sequences were manually edited using BioEdit v7.0.5.3 [30] to be approximately as long as ORF1. All sequences were then aligned using ClustalW v1.81 [31].
The Maximum Likelihood method was used with GTR + F + I + G4 as the model of nucleotide substitution. The phylogenetic tree was made using the IQ-TREE software v.1.6.12, supported by UltraFast Boostrap and SH-aLRT branch testing (1000 replicates) [32,33].

2.5. Statistical Analysis

Continuous variables were expressed as means and standard deviation or medians and interquartile ranges (IQRs, 25th–75th percentile) according to their distribution. Categorical variables were expressed as frequencies and percentages. The Mann–Whitney method was used to compare the number of TTV species with the sample period (pre- and post-transplantation). The correlations between the number of TTV species and TTV viral load, risk of infection or risk of graft rejection were calculated using Pearson’s correlation coefficient (r). A value of p < 0.05 was considered statistically significant. The data were analyzed using the GraphPad Prism software (version 5; GraphPad software, La Jolla, CA, USA).
This study was approved by the CEMIC Ethics Committee (qualified by the Department of Health and Human Services, Buenos Aires, Argentina, HHS no. 124).

3. Results

A total of 107 renal transplant patients were enrolled between November 2018 and April 2021 at CEMIC University Hospital for an ongoing prospective study. Of them, 27 patients were selected, and 30 plasma samples were successfully typed by RCA-NGS sequencing.
The patients’ mean age was 49.5 ± 13.9 years old, and 14/27 (51.8%) were female. Most of them received a graft from a deceased donor (74.1%), and 96.3% received thymoglobulin as immunosuppression induction. All the patients were on steroids, tacrolimus and mycophenolic acid for maintenance immunosuppression (Table 1). All the transplant patients were monitored for HIV, HBV and HCV. None of the patients in our cohort were positive for these viruses. No patient lost the graft up to the first year after transplantation.
A total of 15 TTV species were detected in the 30 plasma samples (Figure 1).

3.1. TTV Diversity

The most prevalent species were TTV 3 (63.3%) and TTV 24 (53.3%), followed by TTV 9 (43.3%) and TTV 13 (43.3%) (Table 2).
The median number of TTV species per patient was 4 (IQR: 2–6.3). Up to 10 TTV species were detected in a single patient.
The difference in the total number of TTV species between the pre- and post-transplant periods was statistically significant (p = 0.021).
In the pre-transplant period (n = 12), the most prevalent species were TTV 3 (75%) and TTV 13 (75%), followed by TTV 18 (58.3%) and TTV 24 (58.3%). The median number of TTV species per sample was 5 (IQR: 4–7.5). In the samples from patients with a living donor, 3/12 (25%), the median number of TTV species was 4 (IQR: 4–5). In the samples from patients with a cadaveric donor, 9/12 (75%), the median number of TTV species was 6 (IQR: 4.5–8) (Table 3, Table S1).
In the post-transplant period (n = 18), the most prevalent species was TTV 3 (55.6%), followed by TTV 24 (44.4%) and TTV 5 (38.9%). The median number of TTV species per sample was 2 (IQR: 1.8–5.5). Particularly, in the period of greatest immunosuppression, at month 3 (n = 7), the median number of TTV species per sample was 3.5 (IQR: 2–7) (Table 3).
Different variants of the same species were detected in five samples, particularly species that were grouped in the latest taxonomic update. Specifically, two different variants of TTV 9 (previously named TTV 11 and TTV 12) were detected in three pre-transplant samples. In one of these samples, two variants of TTV 24 (formerly TTV 22 and TTV 24) were also detected. A different pre-transplant sample also had two variants of TTV 24 (TTV 22 and TTV 24). Finally, two variants of TTV 29 (formerly TTV 28 and TTV 29) were detected in a post-transplant sample (Table 3).

3.1.1. TTV Species Diversity and Viral Load

The TTV median viral load in pre-transplant samples was 4.5 Log10 copies/mL (IQR: 4.2–5.3) and 4.1 Log10 copies/mL (IQR: 3.9–4.8) using the Home Brew PCR and R-GENE® PCR, respectively. At month 3, the TTV median viral load was 8.4 Log10 copies/mL (IQR: 8.2–8.5) and 7.4 Log10 copies/mL (IQR: 6.9–7.7) using the Home Brew PCR and R-GENE® PCR, respectively (Table 3).
No associations between the number of TTV species and the TTV viral load pre- and post-transplantation were found using either assay.
In the pre-transplant samples, Pearson’s correlation coefficients were r = −0.1844 (p = 0.566) and r = −0.0074 (p = 0.982) using the Home Brew PCR and R-GENE® PCR, respectively (Figure 2A). In the post-transplant samples, the coefficients were r = 0.0135 (p = 0.958) and r = −0.2430 (p = 0.331) using the Home Brew PCR and R-GENE® PCR, respectively (Figure 2B).

3.1.2. TTV Species Diversity and Patient Characteristics

No associations between the number of TTV species and patient age or time on dialysis were found.
Pearson’s correlation coefficients were r = 0.062 (95% CI: −0.3051–0.4131, p = 0.745) between the number of TTV species and patient age (Figure S1) and r = 0.1637 (95% CI: −0.2089–0.4949, p = 0.387) between the number of TTV species and time on dialysis (Figure S2).

3.2. TTV Variability

The variability testing performed in three patients with pre- and post-transplant samples showed a variation in the number of TTV species. Specifically, the total number of TTV species went from 4 species to 1; from 5 species to 0 and from 2 species to 1 in each patient (Table 3).
Accordantly, total variability was also observed in the species detected throughout time. The post-transplant samples showed a different species in two patients and a TTV species that could not be sequenced in one patient.

3.3. TTV Species and Viral Infection or Graft Rejection

Among 27 renal transplant recipients, 4 had viral infection (2 CMV cases and 2 BKPyV cases), 3 had graft rejection, 2 had graft rejection and viral infection (CMV and BKPyV) and 18 had neither (Table 3).
No specific TTV species were associated with graft rejection or viral infection.
The number of TTV species detected did not correlate with the risk of graft rejection or viral infection. Pearson’s correlation coefficients were r = 0.112 (p = 0.858) between the number of TTV species and graft rejection and r = −0.308 (p = 0.526) between the number of TTV species and viral infection.

3.4. Other Anelloviruses Detection

Other anelloviruses (such as TTMV and TTMDV) were detected in the plasma samples.
TTMV was simultaneously detected with TTV in a sample obtained 6 months post-transplantation. The same patient had only TTV in the pre-transplant sample.
TTMDV was simultaneously detected with TTV in two patients at months 6 and 12 post-transplantation.
None of the three patients with either TTMV or TTMDV developed viral infection or graft rejection (Table 3). The presence of these anelloviruses was only detected using the TTV R-GENE® PCR.

4. Discussion

In the healthy human virome, one of the most frequent viruses is TTV, a member of the highly diverse Anelloviridae family. There are 22 TTV species (including more than 100 strains), divided into five genogroups (I, II, III, IV and V).
TTV is highly prevalent in the general population (90%) and can reach 100% in renal transplant patients. Studies performed in South America (specifically Brazil and Uruguay) in patients with HCV, HBV or HIV infection have reported that genogroup III species (TTV 13–24) were the most prevalent, followed by species from genogroup I (TTV 1–5) [24,25].
In renal transplant recipients from France and Italy, the most prevalent TTV species was TTV 3 (genogroup I). For French recipients, this was followed by former TTV 27 (now TTV 29, genogroup IV), TTV 24 and TTV 22 (both from genogroup III). In Italian recipients, TTV 13, TTV 15 and TTV 24 followed [23,34].
In our study, plasma samples (taken before and/or after renal transplantation) from a previously studied cohort [35] were selected to determine the TTV species diversity in this subpopulation and the variability over time in three of these patients. Genotyping was performed by whole-genome sequencing, using the RCA-NGS method and ORF1 phylogenetic analysis.
We were able to detect 15 TTV species, being TTV3 and TTV 24 the most prevalent.
In our study, high TTV species diversity was found. Specifically, up to 10 TTV species were detected in a single patient. In contrast, a single TTV species was detected in only three samples.
A high diversity of anelloviruses (mostly TTV) was previously described in blood and tissue samples taken from multiple sites within the same patient, including coinfections with closely related members of the Anelloviridae family (TTMV and TTMDV) [36,37]. In three samples from three different patients, TTMV or TTMDV were simultaneously detected along with TTV at different times post-transplantation.
TTV viral load, determined in plasma samples, has been associated with the immunosuppression state in renal transplant patients [23,35,38]. Therefore, TTV viral load can be used as a biomarker of immunosuppression in transplanted patients.
Given the abundance of TTV species described, a possible hypothesis was that a higher TTV viral load would indicate a higher number of TTV species present. In our study, no correlation between viral load and number of species was found. However, a significantly higher TTV species diversity was detected in the pre-transplant samples compared to post-transplant and in patients with grafts from a deceased donor compared to patients with a living donor.
This high diversity could be attributed to multiple factors, such as blood transfusions, the age of recipients, the earlier start of immunosuppression maintenance therapy in patients with a living donor or/and increased exposure to TTV infection during dialysis.
None of the patients in our cohort had received blood transfusions. Patients with a deceased donor started immunosuppression maintenance as late as 8 days post-transplantation; therefore, their immune system should have been able to control the TTV viral infection more efficiently. However, the observed higher TTV diversity could be attributed to the virus’ evasion mechanisms. Patients with a living donor initiated immunosuppression maintenance before transplantation, causing the initial lymphocyte depletion. This could explain the lower diversity, since there are fewer cells in which TTV can replicate.
Despite previous reports of TTV prevalence increasing with age and the possibility of higher exposure to different TTV species throughout life, no association between the number of TTV species and the age of the recipients was found in our study. Moreover, it was not possible to determine an association between the number of TTV species and the risk of graft rejection or viral infection.
For variability studies (evaluating the same patient at different time points), Kulifaj et al. described a high variability of TTV species when comparing pre- and post-transplantation samples [34]. In our study, all patients showed a variation in the number and type of TTV species detected using these experimental methodologies.
However, in this study, comparison between pre- and post-transplantation for the same individual was limited to three patients. Further studies, including a larger cohort with pre-and post-transplantation samples, are needed to better understand how TTV species change over time and their clinical relevance.
TTV replicates in mononuclear cells, mainly in lymphocytes [12,39]. Renal transplant immunity is largely diminished by both induction immunosuppression therapy (which causes a pronounced initial depletion) and immunosuppression maintenance throughout the post-transplant period. This could explain the observed variability of TTV species over time in the same host. Whether the detection of different TTV species is related to changes in the immune status remains to be clarified.
Since TTV is the most prevalent component in the healthy human virome and multiple species can be concomitantly detected within the same host, it is reasonable to consider that not all TTV species can equally affect or impact the human host. Therefore, it is important to further investigate which species are the most prevalent across the globe, their interaction with the host’s immune system and possible changes in the quantity and proportion of each TTV species over time in the same host.
The limitations of this study are the low number of samples genotyped and the single-center design. However, we were able to demonstrate that a combination of techniques like RCA, UAvI-PCR and NGS is a valuable tool for obtaining whole-genome TTV sequences. This combination was previously used by Spezia et al., 2023 [29] to successfully identify TTV species in patients with Kawasaki disease. These tools will contribute to a deeper understanding of TTV viral epidemiology due to the small number of studies currently available.
Research in this area is ongoing and has become increasingly established since TTV has been proposed as a biomarker of immune status. TTV determination may be an important factor in elucidating viral–host interactions in relation to immune responses.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/v16030432/s1. Figure S1: Association between the number of TTV species and patient age; Figure S2: Association between the number of TTV species and time on dialysis; Table S1: Number of TTV species in plasma samples from renal transplant recipients (pre- and post-transplantation).

Author Contributions

Conceptualization: N.S.R. and M.E.; methodology: N.S.R., P.G.S. and F.F.; validation: N.S.R. and P.G.S.; formal analysis: N.S.R. and P.G.S.; investigation: N.S.R., R.J. and M.E.; resources: N.S.R., G.L., F.M. and M.E.; data curation: N.S.R., writing—original draft preparation: N.S.R., P.G.S., R.J., E.H. and M.E.; writing—review and editing: N.S.R., P.G.S., R.J., F.F., N.B., G.G., E.H., F.A.P., G.L., M.P., F.M. and M.E.; visualization: N.S.R., P.G.S. and R.J.; supervision: N.S.R., F.M. and M.E.; project administration: N.S.R. and M.E.; funding acquisition: N.S.R., N.B., F.M. and M.E. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially supported by a CEMIC and Fundación “Norberto Quirno” grant (Young Investigator Award 2019, given to N.S.R.), a Fundación Allende grant (2019, given to N.B., MD), a Roemmers Foundation grant (2019, given to M.E., PhD), a grant from Ricerca Corrente Linea 1, Italian Ministry of Health; and the European Union’s Horizon 2020 research and innovation program (grant agreement 896932).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by CEMIC Institutional Ethics Committee (qualified by the Department of Health and Human Services, HHS no. 124).

Informed Consent Statement

Written informed consent was obtained from all the subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author (N.S.R., [email protected]).

Acknowledgments

We are very grateful to Philippe Bourgeois and Carole Janis for all the support and effort in providing the TTV R-GENE® kits (bioMérieux, Marcy-l’Etoile, France). From our hospital, we want to thank to Jorge Lopez Camelo for proofreading the manuscript, and the technicians Victoria Ceñal, Walter Tripelhorn and Karina Insua for their sampling assistance.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bendinelli, M.; Pistello, M.; Maggi, F.; Fornai, C.; Freer, G.; Vatteroni, M.L. Molecular properties, biology, and clinical implications of TT virus, a recently identified widespread infectious agent of humans. Clin. Microbiol. Rev. 2001, 14, 98–113. [Google Scholar] [CrossRef]
  2. Rezahosseini, O.; Drabe, C.H.; Sørensen, S.S.; Rasmussen, A.; Perch, M.; Ostrowski, S.R.; Nielsen, S.D. Torque-Teno virus viral load as a potential endogenous marker of immune function in solid organ transplantation. Transplant. Rev. J. 2019, 33, 137–144. [Google Scholar] [CrossRef] [PubMed]
  3. Focosi, D.; Antonelli, G.; Pistello, M.; Maggi, F. Torquetenovirus: The human virome from bench to bedside. Clin. Microbiol. Infect. 2016, 22, 589–593. [Google Scholar] [CrossRef] [PubMed]
  4. Wei, Y.; Chen, M.; Yang, X.; Zhang, L.; Rao, L.; Yuan, F.; Wang, Y.; Gong, J.; Li, L. Molecular characterization of human Torque Teno virus. Biomed. Rep. 2015, 3, 821–826. [Google Scholar] [CrossRef] [PubMed]
  5. Spandole, S.; Cimponeriu, D.; Berca, L.M.; Mihăescu, G. Human anelloviruses: An update of molecular, epidemiological and clinical aspects. Arch. Virol. 2015, 160, 893–908. [Google Scholar] [CrossRef] [PubMed]
  6. Kakkola, L.; Bondén, H.; Hedman, L.; Kivi, N.; Moisala, S.; Julin, J.; Ylä-Liedenpohja, J.; Miettinen, S.; Kantola, K.; Hedman, K.; et al. Expression of all six human Torque teno virus (TTV) proteins in bacteria and in insect cells, and analysis of their IgG responses. Virology 2008, 382, 182–189. [Google Scholar] [CrossRef] [PubMed]
  7. Rosario, K.; Duffy, S.; Breitbart, M. A field guide to eukaryotic circular single-stranded DNA viruses: Insights gained from metagenomics. Arch. Virol. 2012, 157, 1851–1871. [Google Scholar] [CrossRef] [PubMed]
  8. Nishizawa, T.; Okamoto, H.; Konishi, K.; Yoshizawa, H.; Miyakawa, Y.; Mayumi, M. A novel DNA virus (TTV) associated with elevated transaminase levels in posttransfusion hepatitis of unknown etiology. Biochem. Biophys. Res. Commun. 1997, 241, 92–97. [Google Scholar] [CrossRef]
  9. Takahashi, K.; Iwasa, Y.; Hijikata, M.; Mishiro, S. Identification of a new human DNA virus (TTV-like mini virus, TLMV) intermediately related to TT virus and chicken anemia virus. Arch. Virol. 2000, 026930, 979–993. [Google Scholar] [CrossRef]
  10. Ninomiya, M.; Nishizawa, T.; Takahashi, M.; Lorenzo, F.R.; Shimosegawa, T.; Okamoto, H. Communication Identification and genomic characterization of a novel human torque teno virus of 3.2 kb. J. Gen. Virol. 2007, 88, 1939–1944. [Google Scholar] [CrossRef]
  11. Varsani, A.; Opriessnig, T.; Celer, V.; Maggi, F.; Okamoto, H. Taxonomic update for mammalian anelloviruses (family Anelloviridae). Arch. Virol. 2021, 166, 2943–2953. [Google Scholar] [CrossRef]
  12. Maggi, F.; Fornai, C.; Zaccaro, L.; Morrica, A.; Vatteroni, M.L.; Isola, P.; Marchi, S.; Ricchiuti, A.; Pistello, M.; Bendinelli, M. TT Virus (TTV) Loads Associated with Different Peripheral Blood Cell Types and Evidence for TTV Replication in Activated Mononuclear Cells. J. Med. Virol. 2001, 194, 190–194. [Google Scholar] [CrossRef]
  13. Matsubara, H.; Michitaka, K.; Horiike, N.; Kihana, T.; Yano, M.; Mori, T.; Onji, M. Existence of TT virus DNA and TTV-like mini virus DNA in infant cord blood: Mother-to-neonatal transmission. Hepatol. Res. 2001, 21, 280–287. [Google Scholar] [CrossRef]
  14. Morrica, A.; Maggi, F.; Vatteroni, M.L.; Fornai, C.; Pistello, M.; Ciccorossi, P.; Grassi, E.; Gennazzani, A.; Bendinelli, M. TT Virus: Evidence for Transplacental Transmission. J. Infect. Dis. 2000, 181, 3–5. [Google Scholar] [CrossRef] [PubMed]
  15. Navarro, D.; Redondo, N.; Fernández-ruiz, M.; Aguado, J.M. Viruses, friends, and foes: The case of Torque Teno Virus and the net state of immunosuppression. Transpl. Infect. Dis. 2022, 24, e13778. [Google Scholar]
  16. Focosi, D.; Spezia, P.G.; Macera, L.; Salvadori, S.; Navarro, D.; Lanza, M.; Antonelli, G.; Pistello, M.; Maggi, F. Assessment of prevalence and load of torquetenovirus viraemia in a large cohort of healthy blood donors. Clin. Microbiol. Infect. 2020, 26, 1406–1410. [Google Scholar] [CrossRef] [PubMed]
  17. Spandole-Dinu, S.; Cimponeriu, D.G.; Crăciun, A.M.; Radu, I.; Nica, S.; Toma, M.; Alexiu, O.A.; Iorga, C.S.; Berca, L.M.; Nica, R. Prevalence of human anelloviruses in Romanian healthy subjects and patients with common pathologies. BMC Infect. Dis. 2018, 17, 334. [Google Scholar] [CrossRef] [PubMed]
  18. Haloschan, M.; Bettesch, R.; Görzer, I.; Weseslindtner, L.; Kundi, M.; Puchhammer-Stöckl, E. TTV DNA plasma load and its association with age, gender, and HCMV IgG serostatus in healthy adults. Age 2014, 36, 9716. [Google Scholar] [CrossRef] [PubMed]
  19. Maggi, F.; Bendinelli, M. Immunobiology of the Torque Teno Viruses and Other Anelloviruses. Curr. Top. Microbiol. Immunol. 2009, 331, 35–37. [Google Scholar]
  20. Peters, M.A.; Jackson, D.C.; Crabb, B.S.; Browning, G.F. Chicken Anemia Virus VP2 Is a Novel Dual Specificity Protein Phosphatase. J. Biol. Chem. 2002, 277, 39566–39573. [Google Scholar] [CrossRef]
  21. Zheng, H.; Ye, L.; Fang, X.; Li, B.; Wang, Y.; Xiang, X.; Kong, L.; Wang, W.; Zeng, Y.; Ye, L.; et al. Torque Teno Virus (SANBAN Isolate) ORF2 Protein Suppresses NF- kB Pathways via Interaction with IkB Kinases. J. Virol. 2007, 81, 11917–11924. [Google Scholar] [CrossRef]
  22. Kincaid, R.P.; Burke, J.M.; Cox, J.C.; de Villiers, E.M.; Sullivan, C.S. A Human Torque Teno Virus Encodes a MicroRNA That Inhibits Interferon Signaling. PLoS Pathog. 2013, 9, e1003818. [Google Scholar] [CrossRef]
  23. Görzer, I.; Haupenthal, F.; Maggi, F.; Gelas, F.; Kulifaj, D.; Brossault, L.; Puchhammer-Stöckl, E.; Bond, G. Validation of plasma Torque Teno viral load applying a CE-certified PCR for risk stratification of rejection and infection post kidney transplantation. J. Clin. Virol. 2023, 158, 105348. [Google Scholar] [CrossRef]
  24. Cancela, F.; Ramos, N.; Mirazo, S.; Mainardi, V.; Gerona, S.; Arbiza, J. Detection and molecular characterization of Torque Teno Virus (TTV) in Uruguay. Infect. Genet. Evol. 2016, 44, 501–506. [Google Scholar] [CrossRef]
  25. Rosa, A.S.; Araujo, O.C.; Savassi-Ribas, F.; Fernandes, C.A.; Coelho, H.S.; Niel, C.; Villela-Nogueira, C.A.; Araujo, N.M. Prevalence of occult hepatitis B virus infection and Torque teno virus infection and their association with hepatocellular carcinoma in chronic hepatitis C patients. Virus Res. 2017, 242, 166–172. [Google Scholar] [CrossRef]
  26. Laubscher, F.; Hartley, M.-A.; Kaiser, L.; Cordey, S. Genomic Diversity of Torque Teno Virus in Blood Samples from Febrile Paediatric Outpatients in Tanzania: A Descriptive Cohort Study. Viruses 2022, 14, 1612. [Google Scholar] [CrossRef]
  27. Abbas, A.A.; Young, J.C.; Clarke, E.L.; Diamond, J.M.; Imai, I.; Haas, A.R.; Cantu, E.; Lederer, D.J.; Meyer, K.; Milewski, R.K.; et al. Bidirectional transfer of Anelloviridae lineages between graft and host during lung transplantation. Am. J. Transplant. 2019, 19, 1086–1097. [Google Scholar] [CrossRef] [PubMed]
  28. Maggi, F.; Pifferi, M.; Fornai, C.; Andreoli, E.; Tempestini, E.; Vatteroni, M.; Presciuttini, S.; Marchi, S.; Pietrobelli, A.; Boner, A.; et al. TT Virus in the Nasal Secretions of Children with Acute Respiratory Diseases: Relations to Viremia and Disease Severity. J. Virol. 2003, 77, 2418–2425. [Google Scholar] [CrossRef] [PubMed]
  29. Spezia, P.G.; Matsudaira, K.; Filippini, F.; Miyamura, T.; Okada, K.; Nagao, Y.; Ishida, T.; Sano, T.; Pistello, M.; Maggi, F.; et al. Viral load of Torquetenovirus correlates with Sano’s score and levels of total bilirubin and aspartate aminotransferase in Kawasaki disease. Sci. Rep. 2023, 13, 18033. [Google Scholar] [CrossRef] [PubMed]
  30. Hall, T.A. BioEdit: A user friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999, 41, 95–98. [Google Scholar]
  31. Thompson, J.D.; Higgins, D.G.; Gibson, T.J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22, 4673–4680. [Google Scholar] [CrossRef]
  32. Hoang, D.T.; Chernomor, O.; Haeseler AVon Minh, B.Q.; Vinh, L.S. UFBoot2: Improving the Ultrafast Bootstrap Approximation. Mol. Biol. Evol. 2017, 35, 518–522. [Google Scholar] [CrossRef]
  33. Nguyen, L.; Schmidt, H.A.; Haeseler AVon Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol. 2014, 32, 268–274. [Google Scholar] [CrossRef]
  34. Kulifaj, D.; Tilloy, V.; Scaon, E.; Guerin, E.; Essig, M.; Pichon, N.; Hantz, S.; De Bernardi, A.; Joannes, M.; Barranger, C.; et al. Viral metagenomics analysis of kidney donors and recipients: Torque teno virus genotyping and prevalence. J. Med. Virol. 2020, 92, 3301–3311. [Google Scholar] [CrossRef] [PubMed]
  35. Reyes, N.S.; Laham, G.; Boccia, N.; García, G.; Jara, R.; Hermida, E.; Ricarte, C.; Diaz, C.; Soler Pujol, G.; Poletta, F.A.; et al. Prospective cohort study of Torque Teno Virus (TTV) viral load kinetics and the association with graft rejection in renal transplant patients. J. Clin. Virol. 2023, 165, 105501. [Google Scholar] [CrossRef] [PubMed]
  36. Abbas, A.A.; Diamond, J.M.; Chehoud, C.; Chang, B.; Kotzin, J.J.; Young, J.C.; Imai, I.; Haas, A.R.; Cantu, E.; Lederer, D.J.; et al. The Perioperative Lung Transplant Virome: Torque Teno Viruses are Elevated in Donor Lungs and Show Divergent Dynamics In Primary Graft Dysfunction. Am. J. Transplant. 2017, 17, 1313–1324. [Google Scholar] [CrossRef]
  37. Giménez, E.; Monzó, C.; Albert, E.; Fuentes-Trillo, A.; Seda, E.; Piñana, J.L.; Hernández Boluda, J.C.; Solano, C.; Chaves, J.; Navarro, D. Diversity and dynamic changes of anelloviruses in plasma following allogeneic hematopoietic stem cell transplantation. J. Med. Virol. 2020, 93, 5167–5172. [Google Scholar] [CrossRef]
  38. Doberer, K.; Schiemann, M.; Strassl, R.; Haupenthal, F.; Dermuth, F.; Görzer, I.; Eskandary, F.; Reindl-Schwaighofer, R.; Kikić, Ž.; Puchhammer-Stöckl, E.; et al. Torque teno virus for risk stratification of graft rejection and infection in kidney transplant recipients—A prospective observational trial. Am. J. Transplant. 2020, 20, 2081–2090. [Google Scholar] [CrossRef] [PubMed]
  39. Focosi, D.; Macera, L.; Boggi, U.; Nelli, L.C.; Maggi, F. Short-term kinetics of torque teno virus viraemia after induction immunosuppression confirm T lymphocytes as the main replication-competent cells. J. Gen. Virol. 2015, 96, 115–117. [Google Scholar] [CrossRef]
Figure 1. Phylogenetic tree of the TTV species detected in 30 samples from Argentine renal transplant recipients. Reference sequences (n = 190) from GenBank. Maximum Likelihood method was used with GTR + F + I + G4 as the model of nucleotide substitution. Phylogenetic tree was made using IQ-TREE software v1.6.12, supported by UltraFast Boostrap and SH-aLRT branch testing (1000 replicates). TTV species detected in samples from our study are colored.
Figure 1. Phylogenetic tree of the TTV species detected in 30 samples from Argentine renal transplant recipients. Reference sequences (n = 190) from GenBank. Maximum Likelihood method was used with GTR + F + I + G4 as the model of nucleotide substitution. Phylogenetic tree was made using IQ-TREE software v1.6.12, supported by UltraFast Boostrap and SH-aLRT branch testing (1000 replicates). TTV species detected in samples from our study are colored.
Viruses 16 00432 g001
Figure 2. Association between the number of TTV species and TTV viral load. (A) Pearson’s correlation between the number of TTV species and TTV viral load pre-transplantation. (B) Pearson’s correlation between the number of TTV species and TTV viral load post-transplantation.
Figure 2. Association between the number of TTV species and TTV viral load. (A) Pearson’s correlation between the number of TTV species and TTV viral load pre-transplantation. (B) Pearson’s correlation between the number of TTV species and TTV viral load post-transplantation.
Viruses 16 00432 g002
Table 1. Clinical and demographics characteristics of renal transplant patients (n = 27) with TTV sequencing.
Table 1. Clinical and demographics characteristics of renal transplant patients (n = 27) with TTV sequencing.
CharacteristicsRenal Transplant Patients (n = 27)
Age, years old [mean ± SD]49.4 ± 13.9
Gender (female) [n (%)]14 (51.8)
Number of kidney transplantation [n (%)]
   First22 (81.5)
   Second4 (14.8)
   Third1 (3.7)
Cause of end-stage renal disease [n (%)]
Diabetes mellitus 5 (18.5))
Glomerular disease7 (25.9)
Hypertensive kidney disease2 (7.4)
Polycystic kidney disease2 (7.4)
Other specified disorders of the kidney and ureter2 (7.4)
Unspecified chronic renal failure9 (33.3)
CMV serostatus [n (%)]
   D+/R+22 (81.5)
   D−/R+2 (7.4)
   D+/R−3 (11.1)
   D−/R−-
Pre-transplant renal replacement therapy [n (%)]
   Hemodialysis23 (85.2)
Time on dialysis, months [median (IQR)]30.6 (16.3–53.2)
Donor age, years old [mean ± SD]43.5 ± 20.2
Donor gender (male) [n (%)]17 (63.0)
Type of donor [n (%)]
   Deceased donor20 (74.1)
   Living donor7 (25.9)
Number of HLA mismatches [median (IQR)]5 (4–5)
Induction therapy [n (%)]
   Thymoglobulin26 (96.3)
   Basiliximab1 (3.7)
Initial maintenance immunosuppression scheme [n (%)]
   Tacrolimus, mycophenolic acid and steroids27 (100)
SD: standard deviation; IQR: interquartile range; D: donor R: recipient.
Table 2. TTV species prevalence in 30 plasma samples from 27 renal transplant recipients. Clinical and demographics characteristics of renal transplant patients (n = 27) with TTV sequencing.
Table 2. TTV species prevalence in 30 plasma samples from 27 renal transplant recipients. Clinical and demographics characteristics of renal transplant patients (n = 27) with TTV sequencing.
TTV Speciesn SamplesPrevalence (%)
TTV 31963.3
TTV 24 (TTV 22 *, TTV 23 * and TTV 24)1653.3
TTV 9 (TTV 9, TTV 11 * and TTV 12 *)1146.7
TTV 131343.3
TTV 181343.3
TTV 51136.7
TTV 1930.0
TTV 15 (TTV 15 and TTV 16 *)930.0
TTV 29 (TTV 27 *, TTV 28 * and TTV 29)930.0
TTV 21620.0
TTV 17310.0
TTV 19310.0
TTV 20310.0
TTV 613.3
TTV 1013.3
* Former classification.
Table 3. TTV viral load and TTV species in plasma samples from renal transplant patients (n = 27) using two PCR assays.
Table 3. TTV viral load and TTV species in plasma samples from renal transplant patients (n = 27) using two PCR assays.
n Patients ID Samples Time after Transplantation (month)TTV Viral Load
(Log10 copies/mL)
n TTV SpeciesAnelloviruses Viral Infection
(Virus, Days after Transplantation)
Graft Rejection
(Days after Transplantation)
Home Brew PCRR-GENE® PCRTTV SpeciesOther Anelloviruses
1M 1_50945Pre Tx5.45.063, 13, 15, 17, 18, 29-NONO
2M 2_48365Pre Tx4.33.981, 9 (11/12), 13, 17, 18, 19, 21, 24-NONO
3M 3_10226Pre Tx5.03.951, 3, 9 (11/12) *, 13, 24 (22/24) *-NONO
4M 4_50378Pre Tx4.74.063, 9 (12), 13, 15 (16), 18, 24-NONO
5M 5_46007Pre Tx4.44.0413, 15 (16), 18, 24-NONO
6M 6_53410Pre Tx4.22.943, 9 (12), 15 (16), 20-YES (CMV, +163)YES (+38)
7M 7_48016Pre Tx5.64.981, 3, 5, 9 (11/12) *, 10, 13, 19, 24-NONO
8M 8_10348Pre Tx4.64.253, 13, 18, 24, 29-YES (BKPyV, +101)NO
9M 9_49489Pre Tx4.34.1103, 5, 13, 15 (16), 17, 18, 20, 21, 24 (22,24), 29-NONO
10M 10_5117313.90.033, 9 (12), 13-NONO
11M 11_4974814.62.371, 5, 9 (12), 13, 15 (16), 24 (22), 29-NONO
12M 12_5380638.57.053, 5, 15, 18, 21-YES (CMV, +105)NO
13M 13_4996438.27.841, 3, 18, 21-NOYES (+220)
14M 14_5157738.27.773, 6, 9 (11), 18, 21, 24, 29 (28)-YES (BKPyV, +156)YES (+36)
15M 15_5034338.57.7215 (16), 29 (28/29) *-NONO
16M 16_4795538.56.723, 5-NONO
17M 17_4659238.27.273, 5, 9 (11), 13, 18, 24 (22), 29 (28)-YES (CMV, +221)NO
18M 18_4874863.76.123, 5TTMDVNONO
19M 19_4983368.56.913-NONO
20M 20_4966767.35.571, 3, 9, 18, 21, 24, 29 (27)-NONO
21M 21_1023268.47.623, 24-NOYES (+52)
22M 22_10175124.17.151, 5, 13, 18, 21TTMDVNONO
23M 23_11085126.95.0220, 24 (22)-NOYES (+86)
24M 24_10839128.67.225, 24-NONO
25M 25_46868Pre Tx4.03.441, 3, 9 (12), 13-NONO
M 26_51364128.47.2119-
26M 27_51795Pre Tx3.64.351, 5, 15 (16), 18, 29-NONO
M 28_1027463.08.90 TTMV
27M 29_46175Pre Tx6.45.523, 5-YES (BKPyV, +307)NO
M 30_4918398.17.2121-
Pre Tx: pre-transplant sample. TTMV: Torque Teno Mini Virus; TTMDV: Torque Teno Midi Virus. (): TTV species assigned denomination in the previous taxonomical classification. (/) *: variants of TTV species named according to the previous taxonomical classification.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Reyes, N.S.; Spezia, P.G.; Jara, R.; Filippini, F.; Boccia, N.; García, G.; Hermida, E.; Poletta, F.A.; Pistello, M.; Laham, G.; et al. Torque Teno Virus (TTV) in Renal Transplant Recipients: Species Diversity and Variability. Viruses 2024, 16, 432. https://doi.org/10.3390/v16030432

AMA Style

Reyes NS, Spezia PG, Jara R, Filippini F, Boccia N, García G, Hermida E, Poletta FA, Pistello M, Laham G, et al. Torque Teno Virus (TTV) in Renal Transplant Recipients: Species Diversity and Variability. Viruses. 2024; 16(3):432. https://doi.org/10.3390/v16030432

Chicago/Turabian Style

Reyes, Noelia Soledad, Pietro Giorgio Spezia, Raquel Jara, Fabio Filippini, Natalia Boccia, Gonzalo García, Eliana Hermida, Fernando Adrian Poletta, Mauro Pistello, Gustavo Laham, and et al. 2024. "Torque Teno Virus (TTV) in Renal Transplant Recipients: Species Diversity and Variability" Viruses 16, no. 3: 432. https://doi.org/10.3390/v16030432

APA Style

Reyes, N. S., Spezia, P. G., Jara, R., Filippini, F., Boccia, N., García, G., Hermida, E., Poletta, F. A., Pistello, M., Laham, G., Maggi, F., & Echavarria, M. (2024). Torque Teno Virus (TTV) in Renal Transplant Recipients: Species Diversity and Variability. Viruses, 16(3), 432. https://doi.org/10.3390/v16030432

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop