# Southern African Origin of HTLV-1 in Romania

Olivier Cassar, Ambroise Marçais, Olivier Hermine, Emilie Deruelle, Giovanni Begliomini, Alexandru Bardas, Horia Bumbea, Andrei Colita, Daniel Coriu, Viola Maria Popov, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Olivier Cassar, Ambroise Marçais, Olivier Hermine, Emilie Deruelle, Giovanni Begliomini, et al. Southern African Origin of HTLV-1 in Romania. PLoS Neglected Tropical Diseases, 2024, 18 (8), pp.e0012337. 10.1371/journal.pntd.0012337. pasteur-04712352

HAL Id: pasteur-04712352

https://pasteur.hal.science/pasteur-04712352v1

Submitted on 27 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

# Southern African Origin of HTLV-1 in Romania

Olivier Cassar<sup>1*</sup>, Ambroise Marçais<sup>2*</sup>, Olivier Hermine<sup>2</sup>, Emilie Deruelle<sup>2</sup>, Giovanni Begliomini<sup>1</sup>, Alexandru Bardas<sup>3</sup>, Horia Bumbea<sup>4</sup>, Andrei Colita<sup>5</sup>, Daniel Coriu<sup>3</sup>, Viola Maria Popov<sup>6</sup>, Alina Tanase<sup>3</sup>, Philippe Vicente Afonso<sup>1</sup>, Antoine Gessain<sup>1*</sup>

1 Institut Pasteur, Université Paris Cité, CNRS UMR 3569, Unité d'Épidémiologie et Physiopathologie des Virus Oncogènes, Département de Virologie, Paris, France, 2 Hôpital Necker-Enfants Malades, Service d'Hématologie, Paris, France, 3 Fundeni Clinical Institute, University of Medicine and Pharmacy Carol Davila, Bucharest, Romania, 4 Emergency University Hospital, University of Medicine and Pharmacy Carol Davila, Bucharest, Romania, 5 Coltea Hospital, University of Medicine and Pharmacy Carol Davila, Department of Hematology, Bucharest, Romania, 6 Colentina Clinical Hospital, Hematology 1 Department, Bucharest, Romania

* These authors contributed equally to this work.

* olivier.cassar@pasteur.fr (OC); antoine.gessain@pasteur.fr (AG)

<figure><img src="image_3.png" /></figure>

## Abstract

In Europe, most HTLV-1-infected individuals originate from highly endemic regions such as West Indies, sub-Saharan Africa, and South America. The only genuine endemic region for HTLV-1 in Europe is Romania where ATL series have been reported among Romanian patients. Our objective is to better understand the origin of this endemic focus based on a study of the genetic diversity of HTLV-1 in Romanians. DNA was obtained from PBMCs/buffy coats of 11 unrelated HTLV-1-infected individuals of Romanian origin. They include 9 ATL cases and 2 asymptomatic carriers. LTR sequences were obtained for all specimens. Complete genomic HTLV-1 sequences were obtained using four PCR series on 10 specimens. Phylogenetic trees were generated from multiple alignments using HTLV-1 prototypic sequences and the new generated sequences. Most of the complete LTR sequences (756-bp) showed low nucleotide diversity, ranging from 0% to 0.8% difference, and were closely related (less than 0.8% divergence) to the only previously characterized Romanian strain, RKI2. One strain, ROU7, diverged slightly (1.5% on average) from the others. Phylogenetic analyses both on partial LTR and the complete genome demonstrate that the 11 sequences belong to the HTLV-1a cosmopolitan genotype and 10 of them belong to the previously denominated a-TC Mozambique-Southern Africa A subgroup. In this study, we demonstrated that the HTLV-1 present in Romania most probably originated in Southern Africa. As most Romanian HTLV-1 strains are very closely related, we can assume that HTLV-1 has been introduced into the Romanian population recently. Further studies are ongoing to decipher the routes of arrival and dissemination of these HTLV-1 strains, and to date the emergence of this endemic focus in Central Europe.

<figure><img src="image_4.png" /></figure>

OPEN ACCESS

**Citation:** Cassar O, Marçais A, Hermine O, Deruelle E, Begliomini G, Bardas A, et al. (2024) Southern African Origin of HTLV-1 in Romania. PLoS Negl Trop Dis 18(8): e0012337. https://doi.org/10.1371/journal.pntd.0012337

**Editor:** Adly M. M. Abd-Alla, International Atomic Energy Agency, AUSTRIA

**Received:** January 12, 2024

**Accepted:** July 3, 2024

**Published:** August 6, 2024

**Peer Review History:** PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pntd.0012337

**Copyright:** © 2024 Cassar et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** The data that support the findings of this study are available from GenBank (Accession numbers: OR852771 to OR852781 and OR825456 to OR825465).

## Author summary

In Europe, HTLV-1 infections mainly originate in regions such as the West Indies, sub-Saharan Africa and South America. However, Romania has been identified as a true
**Funding:** o AG received funding from the French Government's Investissement d'Avenir Program, Laboratoire d'Excellence "Integrative Biology of Emerging Infectious Diseases" (Grant no. ANR-10-LABX-62-IBEID), the Centre National de la Recherche Scientifique (CNRS-UMR 3569) and from the Institut Pasteur, Paris, France. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. None of the authors received a salary from the funders.

**Competing interests:** The authors have declared that no competing interests exist.

endemic area, due to reported cases of ATL among Romanian patients. To understand this situation, we analyzed the genetic diversity of HTLV-1 in Romanians and collected DNA from 11 infected individuals, most of them ATL cases. The results showed a close genetic relationship between most of the sequences, closely matching the previously identified Romanian RKI2 strain. However, one strain, ROU7, showed a slight divergence. Phylogenetic analysis positioned these sequences in the cosmopolitan HTLV-1a genotype, belonging mainly to the a-TC Mozambique-Southern Africa subgroup. This suggests a probable origin of HTLV-1 in Romania from South Africa, perhaps recently introduced. Further studies aim to elucidate the transmission routes and emergence of this HTLV-1 endemicity in Central Europe.

## Introduction

Human T-cell lymphotropic virus type 1 (HTLV-1) is a human oncoretrovirus infecting at least 5 to 10 million individuals worldwide. It has major public health implications, mainly due to its causal association with severe diseases such as adult T-cell leukemia/lymphoma (ATL) and HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) [1]. Based on genetic variability, seven main HTLV-1 molecular genotypes (a to g) have been described in certain geographic areas and ethnic groups [2,3]. Although the virus is highly endemic in Southern Japan, sub-Saharan Africa, the Caribbean and Australo-Melanesia, its prevalence varies significantly from one population to another [2]. Moreover, HTLV-1 demonstrates remarkable genetic stability, with the observed polymorphism among viral strains correlating with the geographic origin of infected individuals. This limited genetic drift offers a molecular avenue for tracking viral transmission and the migrations of historical infected populations. In Europe, such as France, the United Kingdom or Spain, HTLV-1 is rare, except in people who have immigrated from highly endemic regions such as sub-Saharan Africa, the West Indies and South America, respectively [1,2].

Over the past 30 years, several studies indicate that Romania, a Central European country is a HTLV-1 endemic area in Europe [4–10]. This is based, on the one hand, on studies of blood donors (mandatory screening since 1999) with a seroprevalence of up to 0.053 in first time blood donors (ten times higher than in France and the UK) [11] and, on the other hand, on reports of ATL cases, mainly sporadic ones or small series [5–8], but also recently a large series of 56 patients [9]. However, no specific studies have been carried out to determine the origin of this virus. Therefore, to get new insights into this origin, we decided to study the genetic diversity of HTLV-1 present in Romanians.

## Material and methods

### Ethics statement

Formal written consent was obtained from participants in this study, which was approved by the Human Protection Committee Ile de France II (Comité de Protection des Personnes—CPP IdF II) registered as an Institutional Review Board (IRB registration number: 000001072) and the French Data Protection Authority (Commission Nationale de l'Informatique et des Libertés—CNIL registration number: 1692254).
<table><caption>Table 1. Characteristics of HTLV-1 patients from Romania, including 9 cases with Adult T-cell leukemia (ATL).</caption><thead><tr><th>ID</th><th>Birth country</th><th>Age at diagnosis (y)</th><th>Gender</th><th>ATL clinical status</th><th>HIV status</th><th>PVL (%)&</th></tr></thead><tbody><tr><td>ROU1*</td><td>Romania</td><td>27</td><td>F</td><td>Chronic</td><td>Negative</td><td>40</td></tr><tr><td>ROU2</td><td>Romania</td><td>32</td><td>F</td><td>Acute</td><td>Negative</td><td>&gt;50</td></tr><tr><td>ROU3</td><td>Romania</td><td>39</td><td>F</td><td>Acute</td><td>Negative</td><td>&gt;50</td></tr><tr><td>ROU4</td><td>Romania</td><td>26</td><td>M</td><td>Acute</td><td>Negative</td><td>31</td></tr><tr><td>ROU5</td><td>Romania</td><td>42</td><td>F</td><td>Chronic</td><td>Unknown</td><td>51</td></tr><tr><td>ROU6</td><td>Romania</td><td>48</td><td>M</td><td>Chronic</td><td>Negative</td><td>40</td></tr><tr><td>ROU7*</td><td>Romania</td><td>37</td><td>F</td><td>Carrier</td><td>Negative</td><td>&lt;1</td></tr><tr><td>ROU8*</td><td>Romania</td><td>30</td><td>F</td><td>Carrier</td><td>Negative</td><td>26</td></tr><tr><td>ROU9</td><td>Romania</td><td>53</td><td>F</td><td>Acute</td><td>Unknown</td><td>NA</td></tr><tr><td>PH523</td><td>Romania</td><td>77</td><td>F</td><td>Lymphoma</td><td>Unknown</td><td>NA</td></tr><tr><td>PH630*</td><td>Romania</td><td>17</td><td>F</td><td>Lymphoma</td><td>Unknown</td><td>NA</td></tr></tbody></table>

ATL: adult T-cell leukemia/lymphoma; NA: Not available.

&Proviral load (PVL) expressed as a percentage of infected PBMCs.

#HTLV-1 seropositive mother.

†Blood transfusion history

https://doi.org/10.1371/journal.pntd.0012337.t001

## Studied individuals and data collection

We included eleven HTLV-1 unrelated infected individuals living in Romania, mostly suffering from ATL (Table 1). They were referred for consultation and sometimes therapeutic treatment, mainly to the hematology department of Necker-Enfants Malades hospital (APHP) in Paris. They were all born in Romania, most of them living in Bucarest and selfidentify as Caucasian.

## HIV and HBV screening and blood transfusion history

All but 4 patients (ROU5, ROU9, PH523 and PH630) were tested for HIV. HBV infection status was known for only 5 individuals (ROU1, ROU2, ROU3, ROU7 and ROU9). Two individuals have reported transfusion events that took place in Romania during infancy (ROU1 and ROU8).

## HTLV-1 amplification and sequencing

High molecular weight DNA samples obtained from peripheral blood buffy coats were subjected to four series of PCR (F1-F4) using LTR-*gag* (F1), *pro-pol* (F2), *pol-env* (F3), and Px-LTR (F4) primers, which were designed from highly conserved regions for the major HTLV-1 subtypes. The LTR-*gag* primers are the following: Enh280: 5'-TGACGACAACCCCTCACC TCAA-3' and R2380: 5'-GTCCGGAAAGGGAGGCGTATTAG-3' corresponding to nucleotides 258 to 279 and 2,377 to 2,399 respectively of the prototype ATK-1 sequence (Genbank: J02029). The *pro-pol* primers are the following: F2279: 5'-GGAGCAGACATGACAGTCC TTCC-3' (nt 2,254 to 2,276) and R5005: 5'-GGCGGCTATTAAGACCAGGAAGC-3' (nt 5,002 to 5,024). The *pol-env* primers are the following: F4583: 5'-CAGGAGCCATCTCAGCTAC CC-3' (nt 4,560 to 4,580) and Env22: 5'-GGCGAGGTG GAGTCCTTGGAGGC-3' (nt 60752 to 6,774). Finally, the Px-LTR primers are the following: F4bisfor: 5'-CCCTCCTTTACCC ATCGTTAG-3' (nt 5,941 to 5,961) and F4bisRev: 5'-CGCAGAACAGAAAACGAAAC-3' (nt 8,867 to 8,886). The size of the different generated amplicons are the following: LTR-*gag*, 2,136-bp; *pro-pol*, 2,769-bp; *pol-env*, 2,078-bp and Px-LTR, 2,905-bp.
A PIKO thermocycler (Ozyme, Saint Quentin-en-Yvelines, France) was used with the fol-
lowing amplification conditions for the first 3 fragments (F1-F3): 98°C, 30 s; 40× (98°C, 5 s;
67°C, 5 s; 72°C, 1 mn); 72°C, 1 mn and then for the last fragment (F4): 98°C, 30 s; 40× (98°C,
5 s; 61°C, 5 s; 72°C, 1 mn); 72°C, 1 mn. Reaction tubes were prepared in a dedicated room out-
side the laboratory, with a final volume of 50 µl (DNA matrix, 250 ng; dNTP mix (Roche,
Basel, Switzerland), 40 µM; 5×Phire reaction buffer which contains 1.5 mM MgCl2 at final
reaction concentration (Life Technologies, Courtaboeuf, France), 5 µl; Phire hot start DNA
polymerase (Life Technologies, Courtaboeuf, France), 2 U and 0.5 µM of each oligonucleotide
primer (Eurofins Genomics, Germany). Five µl of amplified DNA was size fractionated by
1.5% agarose gel electrophoresis, and the PCR products (45 µl) were sent for purification and
sequencing reactions to the Eurofins Genomics Company (Eurofins Genomics, Cologne,
Germany).

The 4 fragments obtained were then sequenced using 18 pairs of primers [12] and the Clus-
talW algorithm (MacVector 18.6.1 software, Oxford Molecular) was implemented to align for-
ward and reverse sequences of each segment to derive a consensus sequence of the full LTR
(756-bp) region, and colinearized *gag-pol-env-tax* (7,567-bp) genes.

Phylogenetic analyses

For the LTR and concatenated gag-pol-env-tax genes analyses, sequences were aligned using
the DAMBE program (version 4.2.13). The final alignment was submitted to the Modeltest
program (version 3.6) and the best model was selected according to the Akaike information
criterion. This was then applied to phylogenetic analyses using the PAUP program (version
4.0b10) to infer trees according to Neighbor-Joining (NJ) method. The tree architecture was
later validated by the Maximum Likelihood method, using SeaView (version 5); the robustness
of the clades was estimated by approximate likelihood-ratio test (aLRT).

LTR analyses

Phylogenetic comparison was performed on 485-nucleotide-long LTR alignment of isolates,
including the 11 sequences generated in this study (in bold red) and 70 reference strains. The
analysis was only carried out on 485-bp because the majority of sequences available in the LTR
are not complete and this size, although reduced compared to the complete LTR (756-bp),
made it possible to include a significant number of sequences from the different subgroups
known in the transcontinental genotype. Two Transcontinental Japanese LTR sequences were
used as outgroup (a-JPN). The phylogenetic tree was derived by the neighbor-joining method
using the GTR model (gamma = 0.2884). Horizontal branch lengths are drawn to scale, with
the bar indicating 0.005 nucleotide replacement per site. Numbers on each node indicate the
support probability (approximate likelihood-ratio test, aLRT, calculated using the SeaView 5
software) for each cluster (Fig 1). Next to each sequence, three letters symbolize the country of
origin of the infected individual (mostly IOC country codes): MOZ-Mozambic, RSA-Republic
of South Africa, SWZ-Eswatini, IRN-Iran, KWT-Kuwait, SEN-Senegal, DRC- Democratic
Republic of the Congo, UGA-Uganda, CIV-Côte d'Ivoire, ROU-Romania, TGO-Togo,
AGO-Angola, NIG-Nigeria, JPN-Japan. The topology of the phylogenetic tree was validated
by Maximum Likelihood method. The groups of interest are colored as follows: grey, blue, pur-
ple, orange, green, yellow belong to Transcontinental (TC) HTLV-1a genotype, Mozambic
Southern Africa (B), Mozambic-Middle East, Mozambic, Mozambic Southern Africa (A) and
Japanese (a-JPN), respectively. The GenBank accession nos OR852771 to OR852781 corre-
spond to the 11 new LTR sequences obtained from Romanian ATL patients.
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Phylogenetic tree generated with neighbor-joining (NJ) method on a 485-bp fragment of the LTR region for 81 HTLV-1 available sequences including the 11 generated in this work (in bold red type). The numbers at some nodes of the tree correspond to the probability of each monophyletic groupe (a-LRT). The branch lengths are drawn to scale with bar indicating 0.005-nucleotide replacement by site. The ATK-1 and ATL-25 strains were used as outgroup.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0012337.g001

## Gag-pol-env-tax genes analyses

Phylogenetic analysis was derived from Neighbor-Joining method using the GTR model (gamma = 0.4815) and was confirmed using the Maximum Likelihood method. Phylogenetic tree was obtained from concatenated fragments of the 6,366-nucleotide-long *gag-pol-env-tax* genes, generated from 46 available complete HTLV-1 genomes comprising the sequences of genotypes a, b and c and the ten sequences generated in this work (in red bold type). Five complete Australo-Melanesian HTLV-1c sequences were used as outgroup. The topology of the tree was identical when the different genes were analyzed separately. Horizontal branch lengths are drawn to scale, with the bar indicating 0.01 nucleotide replacement per site. Numbers on each node indicate the aLRT for the supported cluster (Fig 2). The GenBank accession nos OR825456 to OR825465 correspond to the 10 new full HTLV-1 sequences obtained from Romanian ATL patients.

## Results

The 11 HTLV-1 infected individuals from Romania included 9 women (mean age 39.3, range 17–77) and 2 men (mean age 37 years) (Table 1). Of the eleven patients, four had acute ATL, three had chronic ATL and two had ATL lymphoma. In addition, two were asymptomatic carriers. The proviral load—expressed as a percentage of infected peripheral blood mononuclear cells (PBMCs)—ranged from less than 1% to over 50%. We were able to test for the presence of HIV antibodies in 6 people, all of whom were negative (ROU1, ROU2, ROU3, ROU4, ROU7 and ROU8) (Table 1). Regarding hepatitis B virus status, only 5 individuals were tested: one
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Phylogenetic tree generated with neighbor-joining (NJ) method on a 6,366-bp fragment of the HTLV-1 *Gag-Pol-Env-Tax* concatenated genes for 56 available sequences including the 10 generated in this work (in bold red type). The numbers at some nodes of the tree correspond to the probability of each monophyletic group (a-LRT). The branch lengths are drawn to scale with bar indicating 0.01-nucleotide replacement by site. Five HTLV-1c Australo-Melanesian strains were used as outgroup.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0012337.g002

was vaccinated (ROU1), one had antibodies indicating past infection (ROU2), two were negative (ROU3 and ROU7) and the last one was HbS carrier, indicating chronic infection (ROU9). Two patients had a history of blood transfusion performed in Romania at an early age: ROU1 was transfused during the neonatal period, and ROU8 had multiple transfusions during childhood for acute leukemia (5 years old).

## Long Terminal Repeat (LTR) region analysis

The complete LTR sequence was obtained for the 11 samples from PCR fragments F1 and F4. Three sequences were identical (ROU2, ROU5 and ROU8). Most sequences displayed a low nucleotide diversity, i.e. less than 0.8% difference, and were closely related (less than 0.8% divergence) to the only previously characterized Romanian strain, RKI2 [13]. One strain, ROU7, diverged slightly (1.5% on average) from the others.

Phylogenetic analysis was performed on 485-bp-long LTR alignment of isolates, including the 11 sequences generated in this study (in bold red) and 70 reference strains which represent the different subtypes within HTLV-1 transcontinental (a-TC) strains (Fig 1). These subgroups were defined and named according to previous studies based on analysis of the LTR region of strains from very different geographical origins [3,14]. The analysis shows that the 11 new Romanian strains belong to the transcontinental HTLV-1a subgroup (a-TC). In addition, all but one (ROU7) belong to the a-TC Mozambique-Southern Africa A subgroup, in which Romanian RKI2 strain is also present (Fig 1). ROU7 is positioned outside this group.

The 11 new characterized sequences from Romania do not belong to the TC-Japan or TC-Middle East subgroups.
## Colinearized gag-pol-env-tax genes fragment analysis

Due to limited DNA availability for sample ROU8, we were only able to generate the complete HTLV-1 proviral sequence for only 10 out of these 11 samples. Alignment of a 7523 nt long segment of the 10 genomes-the complete genome without LTRs-shows that, as with LTR, 9 strains were very close to each other: 2 are identical (ROU2 and ROU5) and the maximal nucleotide divergence is 0.2%. One strain (ROU7) was more distant with a divergence of around 1%.

Analysis of the sequences showed that the open reading frames (ORFs) were conserved, with the exception of strain PH630, which presents a stop codon in the *pol* gene.

Phylogenetic tree was obtained from concatenated fragments of the 6366-nucleotide-long *Gag-Pol-Env-Tax* genes, generated from 46 available complete HTLV-1 genomes comprising the sequences of genotypes a, b and c and the ten new sequences generated in this work (in red bold type) (Fig 2). This phylogenetic analysis confirms that the 10 new Romanian HTLV-1 strains all belong to the HTLV-1 a-TC subgroup. Notably, 9 of the 10 Romanian strains (sample ROU7 being the exception) form a specific "Romania" clade, within the highly supported a-TC Mozambic-Southern African A subgroup (Fig 2). The topology of the tree was identical when the different genes were analyzed separately.

## Discussion

In this molecular epidemiological study, we demonstrated that the HTLV-1 present in Romanians living in Romania most probably originated in Southern Africa. As most Romanian HTLV-1 strains are very closely related, and sometimes identical, we can assume that HTLV-1 has been introduced into the Romanian population fairly recently with a limited number of strains, which leads to a founder viral effect. Furthermore, analysis of the different HTLV-1 genes separately generated a similar topology of the phylogenetic tree, consistent with the fact that mutations are found scattered throughout the genome. This suggests that these viruses were not generated by recombination [15].

Concerning the presence of a stop codon in the *pol* gene of the proviral strain of sample PH630, it has been suggested that ATL cells are clonally selected notably by acquiring a number of genetic alterations in viral genes (*gag*, *pol*, *tax*) that favor escape from the host immune system [16].

ATL develops preferentially in people infected with HTLV-1 in childhood, and rarely occurs in those infected in adulthood [17]. It is therefore likely that the ATL patients included in this study were infected at an early age. Interestingly, two patients had HTLV-1-seropositive mothers (PH630 and ROU7), suggesting possible mother-to-child transmission. The transfusion of two patients (ROU1 and ROU8) at an early age in Romania may explain their HTLV-1 positivity.

On the basis of the data available at present, we cannot demonstrate, or even suggest, how this virus entered Romania and then spread there. It is tempting to think that there is some similarity between the situation regarding the spread of HTLV-1 and that which existed a few decades ago in Romania for the epidemics of HIV and hepatitis B, particularly in children [18,19]. However, only solid, in depth epidemiological studies specifically seeking for risk factors for HTLV-1 acquisition in the Romanian population and the families of infected individuals, combined with historical data, will be needed to try to better understand the origin of the current situation. The fact that cases of ATL appeared clinically several decades ago in patients who were then on average 40 years old (this report, 5, 9) suggests occurences of HTLV-1 infection in Romania dating back to nearly 60 years ago. It should be noted that this epidemiological work will certainly be difficult to carry out due to the retrospective aspect of these studies,
which therefore concern facts that are already old, associated with public health decisions taken at the time. The recent discovery that Moldavia, a small border country to the east of Romania, also appears to be an area of high HTLV-1 endemicity, based on blood donors studies, should also encourage local public health authorities to carry out this type of epidemiological study [20].

Above all, it is essential to pursue surveillance and research efforts to limit the spread of this virus within the Romanian population.

## Acknowledgments

We acknowledge the assistance of Drs Frenoy N. from Paul Brousse hospital (Villejuif, France) and Baumelou E. from Foch hospital (Suresnes, France), where some clinical specimens were collected.

## Author Contributions

**Conceptualization:** Olivier Cassar, Ambroise Marçais, Olivier Hermine, Philippe Vicente Afonso, Antoine Gessain.

**Data curation:** Olivier Hermine, Emilie Deruelle.

**Formal analysis:** Olivier Cassar, Giovanni Begliomini, Philippe Vicente Afonso, Antoine Gessain.

**Funding acquisition:** Antoine Gessain.

**Investigation:** Olivier Cassar, Ambroise Marçais, Olivier Hermine, Emilie Deruelle, Alexandru Bardas, Horia Bumbea, Andrei Colita, Daniel Coriu, Viola Maria Popov, Alina Tanase, Antoine Gessain.

**Methodology:** Olivier Cassar, Olivier Hermine, Philippe Vicente Afonso, Antoine Gessain.

**Project administration:** Olivier Hermine.

**Resources:** Ambroise Marçais, Emilie Deruelle.

**Supervision:** Olivier Hermine, Antoine Gessain.

**Validation:** Olivier Cassar, Olivier Hermine, Antoine Gessain.

**Writing – original draft:** Olivier Cassar, Philippe Vicente Afonso, Antoine Gessain.

**Writing – review & editing:** Olivier Cassar, Emilie Deruelle, Philippe Vicente Afonso, Antoine Gessain.

## References

1. Schierhout G, McGregor S, Gessain A, Einsiedel L, Martinello M, Kaldor J. Association between HTLV-1 infection and adverse health outcomes: a systematic review and meta-analysis of epidemiological studies. Lancet Infect Dis. 2020 Jan; 20(1):133–43. https://doi.org/10.1016/S1473-3099(19)30402-5 PMID: 31648940

2. Gessain A, Cassar O. Epidemiological Aspects and World Distribution of HTLV-1 Infection. Front Microbiol. 2012; 3:388. https://doi.org/10.3389/fmicb.2012.00388 PMID: 23162541

3. Afonso PV, Cassar O, Gessain A. Molecular epidemiology, genetic variability and evolution of HTLV-1 with special emphasis on African genotypes. Retrovirology. 2019; 16:39. https://doi.org/10.1186/s12977-019-0504-z PMID: 31842895

4. Paun L, Ispas O, Del Mistro A, Chieco-Bianchi L. HTLV-I in Romania. Eur J Haematol. 1994 Feb; 52 (2):117–8. https://doi.org/10.1111/j.1600-0609.1994.tb01297.x PMID: 8119382
5. Veelken H, Kohler G, Schneider J, Dierbach H, Mertelsmann R, Schaefer HE, et al. HTLV-I-associated adult T cell leukemia/lymphoma in two patients from Bucharest, Romania. Leukemia. 1996 Aug; 10 (8):1366-9. PMID: 8709646

6. Lin BT, Musset M, Szekely AM, Alexandre J, Fraitag S, Bodemer C, et al. Human T-cell lymphotropic virus-1-positive T-cell leukemia/lymphoma in a child. Report of a case and review of the literature. Arch Pathol Lab Med. 1997 Dec; 121(12):1282-6. PMID: 9431320

7. Tanase AD, Colita A, Craciun OG, Lipan L, Varady Z, Stefan L, et al. Allogeneic Stem Cell Transplantation for Adult T-Cell Leukemia/Lymphoma-Romanian Experience. J Clin Med. 2020 Jul 28; 9(8). https://doi.org/10.3390/jcm9082417 PMID: 32731502

8. Shtalrid M, Shvidel L, Korenfeld R, Duek A, Landau Z, Berrebi A. HTLV-1 associated adult T-cell leukemia/lymphoma in Israel: report of two patients of Romanian origin. Haematologica. 2005 Mar; 90(3): ECR13. PMID: 15753054

9. Bumbea H, Marçais A, Coriu D, Tanase AD, Colita A, Bardas A, et al. Results from Treatment of a Large Cohort of ATL Patients from a Country with High HTLV-1 Prevalence. Blood. 2019;134.

10. Stienlauf S, Yahalom V, Schwartz E, Shinar E, Segal G, Sidi Y. Epidemiology of human T-cell lymphotropic virus type 1 infection in blood donors, Israel. Emerg Infect Dis. 2009 Jul; 15(7):1116-8. https://doi.org/10.3201/eid1507.080796 PMID: 19624934

11. Laperche S, Worms B, Pillonel J. Blood safety strategies for human T-cell lymphotropic virus in Europe. Vox Sang. 2009 Feb; 96(2):104-10. https://doi.org/10.1111/j.1423-0410.2008.01136.x PMID: 19076337

12. Cassar O, Einsiedel L, Afonso PV, Gessain A. Human T-cell lymphotropic virus type 1 subtype C molecular variants among Indigenous Australians: New insights into the molecular epidemiology of HTLV-1 in Australo-Melanesia. PLoS Negl Trop Dis. 2013. https://doi.org/10.1371/journal.pntd.0002418 PMID: 24086779

13. Ellerbrok H, Fleischer C, Salemi M, Reinhardt P, Ludwig WD, Vandamme AM, et al. Sequence analysis of the first HTLV-I infection in Germany without relations to endemic areas. AIDS Res Hum Retroviruses. 1998 Sep 1; 14(13):1199-203. https://doi.org/10.1089/aid.1998.14.1199 PMID: 9737591

14. Vicente AC, Gudo ES, Iniguez AM, Otsuki K, Bhatt N, Abreu CM, et al. Genetic characterization of human T-cell lymphotropic virus type 1 in Mozambique: transcontinental lineages drive the HTLV-1 endemic. PLoS Negl Trop Dis. 2011; 5(4):e1038. https://doi.org/10.1371/journal.pntd.0001038 PMID: 21532745

15. Cassar O, Desrames A, Marcais A, Gout O, Taylor GP, Hermine O, et al. Multiple recombinant events in human T-cell Leukemia virus Type 1: complete sequences of recombinant African strains. Emerg Microbes Infect. 2020 Dec; 9(1):913-23. https://doi.org/10.1080/22221751.2020.1752117 PMID: 32249692

16. Kataoka K, Nagata Y, Kitanaka A, Shiraishi Y, Shimamura T, Yasunaga J, et al. Integrated molecular analysis of adult T cell leukemia/lymphoma. Nature genetics. 2015; 47(11), 1304-1315. https://doi.org/10.1038/ng.3415 PMID: 26437031

17. Watanabe T. Adult T-cell leukemia: molecular basis for clonal expansion and transformation of HTLV-1-infected T cells. Blood. 2017 Mar 2; 129(9): 1071-1081. https://doi.org/10.1182/blood-2016-09-692574 PMID: 28115366

18. Haukenes G, Brinchmann-Hansen K, & Macovei O. Prevalence of hepatitis B and C and HIV antibodies in children in a Romanian orphanage. APMIS. 1992; 100(7-12), 757-761.

19. Hersh BS, Popovici F, Jezek Z, Satten GA, Apetrei RC, Beldescu N, et al. Risk factors for HIV infection among abandoned Romanian children. Aids. 1993; 7(12), 1617-1624. https://doi.org/10.1097/00002030-199312000-00012 PMID: 8286071

20. Rosadas C., Harvala H., Davison K., & Taylor G. P. HTLV-1 screening of blood donations: We are systematically missing opportunities. Br J Haematol. 2023; 202(6), 1220-1223. https://doi.org/10.1111/bjh.18988 PMID: 37487701