# Risk factors for HTLV-1 infection in Central Africa: A rural population-based survey in Gabon

Delia Doreen Djuicy, Augustin Mouinga-Ondémé, Olivier Cassar, Jill-Léa

Ramassamy, Antony Idam Mamimandjiami, Rodrigue Bikangui, Arnaud

Fontanet, Antoine Gessain

► To cite this version:

Delia Doreen Djuicy, Augustin Mouinga-Ondémé, Olivier Cassar, Jill-Léa Ramassamy, Antony Idam Mamimandjiami, et al. Risk factors for HTLV-1 infection in Central Africa: A rural population-based survey in Gabon. PLoS Neglected Tropical Diseases, 2018, 12 (10), pp.e0006832. 10.1371/journal.pntd.0006832. pasteur-02975273

HAL Id: pasteur-02975273

https://pasteur.hal.science/pasteur-02975273

Submitted on 22 Oct 2020

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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

RESEARCH ARTICLE

# Risk factors for HTLV-1 infection in Central Africa: A rural population-based survey in Gabon

Delia Doreen Djuicy¹,²,³,⁴, Augustin Mouinga-Ondémé³, Olivier Cassar¹,², Jill-Léa Ramassamy¹,², Antony Idam Mamimandjiami¹,²,³,⁴, Rodrigue Bikangui³, Arnaud Fontanet⁵,⁶, Antoine Gessain¹,²*

¹ Institut Pasteur, Unité d'Epidémiologie et Physiopathologie des Virus Oncogènes, Département de Virologie, Paris, France, ² CNRS, UMR 3569, Paris, France, ³ Centre International de Recherches Médicales de Franceville, Groupe des Rétrovirus Animaux, BP: Franceville, Gabon, ⁴ Ecole Doctorale Régionale d'Afrique centrale, Infectiologie Tropicale, BP: Franceville, Gabon, ⁵ Institut Pasteur, Unité de Recherche et d'Expertise Epidémiologie des Maladies Emergentes, Département d'Infection et Epidémiologie, Paris, France, ⁶ Conservatoire National des Arts et Métiers, Unité PACRI, Paris, France

* agessain@pasteur.fr

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

OPEN ACCESS

**Citation:** Djuicy DD, Mouinga-Ondémé A, Cassar O, Ramassamy J-L, Idam Mamimandjiami A, Bikangui R, et al. (2018) Risk factors for HTLV-1 infection in Central Africa: A rural population-based survey in Gabon. PLoS Negl Trop Dis 12(10): e0006832. https://doi.org/10.1371/journal.pntd.0006832

**Editor:** Edgar M. Carvalho, Hospital Universitário Professor Edgard Santos, BRAZIL

**Received:** June 22, 2018

**Accepted:** September 10, 2018

**Published:** October 12, 2018

**Copyright:** © 2018 Djuicy 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:** All relevant data are within the paper and its Supporting Information files.

**Funding:** DDD was supported by the "Agence Universitaire de la Francophonie" (AUF), the Institut Pasteur in Paris and the "Centre international de Recherches Médicales de Franceville" (CIRMF) in Gabon. This study received funding from the CNRS (UMR 3569) (AG); the Institut Pasteur, France (AG), and through the "Investissement d'Avenir" as

## Abstract

### Background

Human T-Lymphotropic Virus type 1 (HTLV-1) is a human oncoretrovirus that infects at least 5 to 10 million people worldwide and is associated with severe diseases. Africa appears as the largest HTLV-1 endemic area. However, the risk factors for the acquisition of HTLV-1 remain poorly understood in Central Africa.

### Methods

We conducted an epidemiological survey between 2013 and 2017, in rural areas of 6 provinces of Gabon, in a rainforest environment. Epidemiological data were obtained and blood samples were collected after informed consent. Plasma were screened for HTLV-1 antibodies by ELISA and the positive samples were then tested by Western blot (WB). Genomic DNA derived from buffy-coat was subjected to two semi-nested PCRs amplifying either HTLV-1 *env* gene or LTR region fragments.

### Results

We recruited 2,060 individuals over 15 years old, including 1,205 men and 855 women (mean age: 49 years). Of these, 299 were found to be ELISA HTLV-1/2 seropositive. According to WB criteria, 136 were HTLV-1 (6.6%), 25 HTLV-1/2 (1.2%) and 9 HTLV seroneffective (0.4%). PCR results showed that 146 individuals were positive for at least one PCR: 104 for the *env* gene and 131 for the LTR region. Based on both serological and molecular results, 179 individuals were considered infected with HTLV-1, leading to an overall prevalence of 8.7%. The distribution of HTLV-1 infection was heterogeneous across the country. Based on multivariable analyses, female gender, increasing age, ethnicity (Pygmy) and multiple hospitalizations (more than 5 times) were found to be independent risk factors for
part of a "Laboratoire d'Excellence" (LabEx) French
research program: Integrative Biology of Emerging
Infectious Diseases (ANR10-LBX-62 IBEID) (AG).
JLR was financially supported by the European
Union (FOOD/2016/379-660). The funders had no
role in study design, data collection and analysis,
decision to publish, or preparation of the
manuscript.

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

HTLV-1 infection. Furthermore, a non-human primate bite appeared to be marginally associ-
ated with a higher risk of HTLV-1 infection.

Conclusion

Based on state-of-the-art serological and molecular methods, we have demonstrated that rural adult populations in Gabon are highly endemic for HTLV-1. Our results regarding risk factors should lead to public health actions aiming to reduce HTLV-1 transmission.

Author summary

Human T-Lymphotropic Virus type 1 (HTLV-1) infects at least 5 to 10 million people worldwide and is associated with severe diseases. Africa, particularly the sub-Saharan region, is considered as the largest HTLV-1 endemic area. HTLV-1 can be acquired through three main routes: mother-to-child, sexual and through contaminated blood products. In addition, zoonotic transmission from non-human primates still occurs in specific rural African populations. Gabon is considered to be a very high HTLV-1 endemic area. However, risk factors for such a retroviral infection remain poorly known in Central Africa. In this study, using state-of-the-art serological and molecular tests, we aimed to get new insights on risk factors associated with HTLV-1 infection in rural Gabon. We report an overall prevalence of HTLV-1 of 8.7%, demonstrating that HTLV-1 is highly endemic in rural adult Gabonese populations. In the study population, we identified female gender, increasing age, ethnicity (Pygmy) and multiple hospitalizations (more than 5 times) as independent risk factors for HTLV-1 infection. Furthermore, a non-human primate bite appeared to be marginally associated with a higher risk of HTLV-1 infection. Our results regarding the risk factors for such a retroviral infection are important for public health actions aiming to reduce HTLV-1 transmission.

**Introduction**

Human T-lymphotropic virus type 1 (HTLV-1), the first human retrovirus discovered [1], is the etiological agent of several pathologies, mainly a very severe T-cell lymphoproliferation named Adult T-Cell Leukemia Lymphoma (ATLL) and a chronic disabling neuro-myelopathy, the Tropical Spastic Paraparesis/HTLV-1 Associated Myelopathy (TSP/HAM) [1-4]. HTLV-1 is not ubiquitously distributed worldwide. Indeed, it is mainly present in foci where viral prevalence can reach 2 to 40% in adults, depending on age, sex and geography. The most important HTLV-1 endemic areas are: the Southern part of the Japanese archipelago, several areas in South America and the Caribbean basin as well as some areas of Australo-Melanesia, Iran and large regions of sub-Saharan Africa. This human oncoretrovirus is estimated to infect at least 5 to 10 million people worldwide [5, 6]. While the great majority of HTLV-1 infected individuals remains asymptomatic throughout their life, ATLL and TSP/HAM occur in 2 to 7% of them [7]. The three main routes of HTLV-1 transmission are: mother-to-child through prolonged breastfeeding (mostly over 6 months) [8], sexual (mainly from male to female) [9] and by blood products contaminated with infected cells [10]. HTLV-1 originates from its simian counterpart STLV-1, which is highly prevalent in several Non-Human Primates (NHPs) species. Zoonotic transmission of STLV-1 still occurs mainly through severe NHP bites, at least in Central and West Africa [11, 12].
Sub-Saharan Africa is considered as the largest HTLV-1 endemic area accounting for at least half of the infected individuals worldwide (2.5 to 5 millions). However, the situation of HTLV-1 in Africa is not well known. Indeed, the majority of previous studies have been carried out either on very specific populations such as pregnant women, blood donors or hospitalized patient series, or in heterogeneous and relatively small groups of rural or urban inhabitants of a specific town or village or area [5, 13]. All these groups are far from being representative of the population of a given region or country. Furthermore, many of these prevalence studies only include serological analyses, without molecular detection of HTLV-1 proviral DNA. Lastly, there are still few epidemiological studies investigating HTLV-1 acquisition risk factors in the African continent [14–17]. Gabon, located in Central Africa, appears to be highly endemic for HTLV-1 and is an appropriate area for the study of HTLV-1 acquisition risk factors [18–21]. Indeed, to date, no public health action has been undertaken to prevent dissemination of this virus in Gabon, while the four modes of transmission are still prevalent in this country.

In this work, we performed a large epidemiological study among adult inhabitants of rural Gabon. We used state-of-the-art assays for both serological and molecular diagnoses. Our first objective was to provide better knowledge of the prevalence and the geographical distribution of HTLV-1 infection in Gabon. The second was to get new insights into the risk factors for the acquisition of HTLV-1 among adults living in this highly endemic area.

# Material and methods

## Study population

Between 2013 and 2017, we conducted epidemiological surveys involving inhabitants of rural areas living mainly in a rainforest environment, in 6 out the 9 provinces of Gabon. A systematic approach for the enrolment of adults (over 15 years old) was carried out in the populations from all reachable villages and settlements, scattered along side roads and tracks across the primary tropical forest. We included all volunteers. A standardized questionnaire was used to collect epidemiological data. Besides demographic information, we focused our questions on the potential risk factors for HTLV-1 acquisition including scarification, hospitalization, transfusion and contacts with wild animals and bush meat during hunting and/or butchering activities (especially with monkeys and apes). A 5 to 10 ml blood sample was collected on EDTA from all consenting individuals. Plasma and buffy coat were obtained 48 to 72 h after sampling and kept frozen at -80°C.

## Ethics

Ethical approval was obtained from the Comité National d'Ethique of Gabon (Permit number PROT 0011/2013/SG/CNE). Prior to field sampling, individual written informed consent was obtained from all participants after detailed information and explanations of the study were provided to the community. Written informed consent for children was obtained from their parents or recognized guardians.

## HTLV-1 serological tests

Plasma samples were screened for HTLV-1/2 antibodies by ELISA tests (HTLV-I/II ELISA 4.0, MP Biomedicals), and confirmatory Western blot assays (HTLV BLOT 2.4, MP Biomedicals) for all ELISA positive samples. Results were interpreted according to the manufacturer's instructions. Thus, the samples were classified as HTLV-1 seropositive, HTLV-2 seropositive, both HTLV-1/2 seroreactive, HTLV seroreactive, indeterminate and seronegative.
## HTLV-1 molecular tests

High molecular weight DNA was extracted from peripheral blood buffy coats (PB-BC) using the QIAamp blood minikit (Qiagen). DNA samples were determined as amplifiable after a polymerase chain reaction (PCR) positive result, performed on the human β-globin gene. One microgram of each DNA sample was then tested by two different semi-nested PCRs: according to protocols previously described [22]: a first one to amplify a 522-bp fragment of the *env* gene, and a second one to amplify a 450-bp fragment of the LTR-3' region of HTLV-1. A sample was considered as PCR positive when a band of the expected size was clearly detected on agar gel.

## Positive criteria for HTLV-1 infection

An individual was considered as infected with HTLV-1 if the WB profile was either HTLV-1, or HTLV-1/2 or HTLV. Furthermore, an individual with an indeterminate WB profile associated with at least one HTLV-1 positive PCR, was also considered as HTLV-1 infected.

## Statistical methods

Age was compared between Pygmy and Bantu ethnic groups using the Student *t* test. HTLV-1 infection prevalence was compared between provinces using the Chi-square test, and between age groups using the Chi-square test for trends. Logistic regression was used to identify factors associated with HTLV-1 infection and to estimate odds ratios (OR) and 95% confidence intervals (95%CI). ORs were presented as crude or age and gender adjusted (univariable analysis) and all variables adjusted (multivariable analysis). Variables retained in the final model were considered statistically significant if *p* values were < 0.05, and marginally significant if *p* values were < 0.07. Multiple imputations were performed to replace missing values for the “number of hospitalizations” variable, using the multivariable imputation by chained equations (MICE) method [23]. All analyses were performed using STATA 15.0 software (Stata Corporation, College Station, TX, USA).

## Results

We recruited 2,060 individuals (older than 15), including 1,205 men and 855 women (median age: 48, interquartile range: 33–64), originating from the 6 investigated provinces of Gabon (Fig 1A). A large majority of them (1,797) belonged to the Bantu ethnic group, and 261 belonged to the Pygmy group (Fig 1B). Of note, the Pygmy group included younger individuals than the Bantu one (mean age 37 years vs 51, p<0.001).

## HTLV-1 serological results

Among the 2,060 plasma samples tested by ELISA, 299 were found HTLV-1/2 seropositive. Using WB criteria provided by the manufacturer, 136 were HTLV-1 seropositive (6.6%), 25 HTLV-1/2 (1.2%), 10 HTLV-2 (0.5%), 9 HTLV seroreactive (0.4%) and 85 sero-indeterminate (4.1%). Moreover, 34 were seronegative. Based on strict WB results as defined by the manufacturer's assay (HTLV Blot 2.4, MP Biomedicals) and including only HTLV-1 and HTLV-1/2 profiles, the overall HTLV-1 seroprevalence was 7.8% (95%CI 6.7–9.1, 161/2,060) (Table 1).

## HTLV-1 molecular detection results

Molecular amplification was performed on DNA extracted from the PB-BC of 299 ELISA positive individuals. The results showed that 146 were positive for at least one PCR (104 for the *env* gene, 131 for the LTR region), and 89 for both PCRs (Table 1).
<figure><img src="image_4.png" /><figcaption>B</figcaption></figure>

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

Fig 1. Age distribution according to gender (A) and ethnic group (B).

https://doi.org/10.1371/journal.pntd.0006832.g001
<table><caption>Table 1. PCR results according to the HTLV Western blot profiles.</caption><thead><tr><th colspan="2">WB results</th><th colspan="5">PCR results</th></tr><tr><th>WB profiles</th><th>Number of<br>WB profile<br>(N)</th><th>env</th><th>LTR-3'</th><th>env<br>and LTR-3'</th><th>env<br>or LTR-3'<br>(n)</th><th>%<br>(n/N)</th></tr></thead><tbody><tr><td>HTLV-1</td><td>136</td><td>77</td><td>101</td><td>67</td><td>111</td><td>82<br>(111/136)</td></tr><tr><td>HTLV-2</td><td>10</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0<br>(0/10)</td></tr><tr><td>HTLV-1/2</td><td>25</td><td>15</td><td>21</td><td>15</td><td>21</td><td>84<br>(21/25)</td></tr><tr><td>HTLV</td><td>9</td><td>4</td><td>3</td><td>2</td><td>5</td><td>55.6<br>(5/9)</td></tr><tr><td>INDETERMINATE</td><td>85</td><td>8</td><td>6</td><td>5</td><td>9</td><td>10.6<br>(9/85)</td></tr><tr><td>NEGATIVE</td><td>34</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0<br>(0/34)</td></tr><tr><td>Total</td><td>299</td><td>104</td><td>131</td><td>89</td><td>146</td><td>49<br>(146/299)</td></tr></tbody></table>

%: proportion of at least one PCR positive result (env gene or LTR-3' region) according to WB profiles

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

## HTLV-1 prevalence in the studied population

Based on both serological and molecular results, 179 individuals were considered as HTLV-1 infected. This includes: i) the 170 individuals with a WB profile being either HTLV-1, HTLV-1/2 or HTLV, and ii) the 9 individuals with an indeterminate WB profile associated with at least one positive HTLV-1 PCR test (Table 1). Thus, 179 HTLV-1 infected persons were included in the epidemiological analysis (Fig 2A and 2B). This leads to an overall HTLV-1 infection prevalence of 8.7% (95%CI 7.5–10.0).

## Risk factors for HTLV-1 infection in the studied population

Univariable analysis revealed that women are significantly more likely to be HTLV-1 infected than men (age-adjusted OR 1.69; 95%CI 1.2–2.3, Table 2). Indeed, 97 of 855 women (11.35%) and 82 of 1,205 men (6.8%) were HTLV-1 infected. As illustrated (Fig 2A and Table 2), the age specific seroprevalence gradually increases in the study population (p <0.001), reaching 17.3% among women of the 61 to 70 age group. Furthermore, Pygmy individuals are at increased risk of being HTLV-1 seropositive compared to Bantu people, after adjustment for age and sex (Adjusted OR 1.67; 95%CI 1.04–2.68, Table 2 and Fig 2B).

Individuals included in this study originate from 6 of the 9 Gabonese provinces (Fig 3). HTLV-1 seroprevalence varies according to the geographic origin of infected individuals (Fig 3 and Table 2). Ogooué-Ivindo and Ogooué-Lolo provinces present a significantly higher HTLV-1 prevalence (respectively 14% and 11%; p = 0.027) compared to Ngounié, Woleu-Ntem and Nyanga provinces, where lower infection rates were recorded. Moreover, the Haut-Ogooué province shows an intermediate HTLV-1 seroprevalence (9.5%), which is not significantly different from the other provinces (Fig 3 and Table 2).

Indicators for HTLV-1 parenteral transmission were also assessed. According to crude OR, HTLV-1 infection is associated with behaviors leading to scarifications (crude OR 1.46; 95%CI 1.04–2.05), an history of blood transfusion (crude OR 1.61; 95%CI 1.06–2.44), and an increased number of hospitalization (more than five times, crude OR 2.96; 95%CI 1.56–5.62),
A

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

B

<figure><img src="image_7.png" /><figcaption>**Fig 2.** HTLV-1 seroprevalence according to age and sex (A), and ethnic group (B).</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006832.g002
<table><caption>Table 2. Univariable analysis of risk factors associated with HTLV-1 infection in the study population.</caption><thead><tr><th rowspan="2"></th><th colspan="2">HTLV-1 results</th><th rowspan="2">Crude OR (95% CI)</th><th rowspan="2">P-value</th><th rowspan="2">Adjusted OR<sup>a</sup> (95% CI)</th><th rowspan="2">P-value</th></tr><tr><th>Negative</th><th>Positive (%)</th></tr></thead><tbody><tr><td><strong>Gender</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>M</td><td>1,123</td><td>82 (6.80)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>F</td><td>758</td><td>97 (11.35)</td><td>1.75 (1.29-2.39)</td><td>0.0003***</td><td>1.69 (1.23-2.3)</td><td>0.001***</td></tr><tr><td><strong>Age</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>15-30</td><td>404</td><td>18 (4.27)</td><td>1</td><td>&lt;0.001***</td><td>1</td><td>&lt;0.001***</td></tr><tr><td>31-40</td><td>324</td><td>17 (4.99)</td><td>1,18 (0.60-2.32)</td><td></td><td>1.20 (0.61-2.37)</td><td></td></tr><tr><td>41-50</td><td>312</td><td>37 (10.60)</td><td>2.66 (1.49-4.76)</td><td></td><td>2.76 (1.54-4.95)</td><td></td></tr><tr><td>51-60</td><td>296</td><td>26 (8.07)</td><td>1.97 (1.06-3.66)</td><td></td><td>1.99 (1.07-3.7)</td><td></td></tr><tr><td>61-70</td><td>255</td><td>39 (13.27)</td><td>3.43 (1.92-6.13)</td><td></td><td>3.34 (1.87-5.97)</td><td></td></tr><tr><td>&gt;70</td><td>290</td><td>42 (12.65)</td><td>3.25 (1.83-5.76)</td><td></td><td>3.11 (1.75-5.53)</td><td></td></tr><tr><td>Missing</td><td>0</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Ethnic group</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bantus</td><td>1,644</td><td>153 (8.51)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Pygmies</td><td>235</td><td>26 (9.96)</td><td>1.19 (0.77-1.84)</td><td>0.4382</td><td>1.67 (1.04-2.68)</td><td>0.034*</td></tr><tr><td>Missing</td><td>2</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Province</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Nyanga</td><td>171</td><td>10 (5.52)</td><td>1</td><td>0.027**</td><td>1</td><td>0.013**</td></tr><tr><td>Ngounié</td><td>334</td><td>31 (8.49)</td><td>1.59 (0.76-3.31)</td><td></td><td>1.48 (0.7-3.11)</td><td></td></tr><tr><td>Haut-Ogooué</td><td>494</td><td>52 (9.52)</td><td>1.8 (0.89-3.62)</td><td></td><td>2.09 (1.03-4.24)</td><td></td></tr><tr><td>Ogooué-Lolo</td><td>228</td><td>28 (10.94)</td><td>2.1 (0.99-4.44)</td><td></td><td>2.24 (1.05-4.78)</td><td></td></tr><tr><td>Ogooué-Ivindo</td><td>129</td><td>21 (14.00)</td><td>2.78 (1.27-6.12)</td><td></td><td>3.54 (1.59-7.87)</td><td></td></tr><tr><td>Woleu-Ntem</td><td>525</td><td>37 (6.58)</td><td>1.21 (0.59-2.47)</td><td></td><td>1.58 (0.76-3.3)</td><td></td></tr><tr><td>Missing</td><td>0</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Hunting (among men)</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>273</td><td>14 (4.88)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>841</td><td>68 (7.48)</td><td>1.58 (0.87-2.85)</td><td>0.1284</td><td>1.34 (0.75-2.41)</td><td>0.321</td></tr><tr><td>Missing</td><td>9</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Hunting NHP (among men)</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>363</td><td>22 (5.71)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>689</td><td>58 (7.76)</td><td>1.39 (0.84-2.31)</td><td>0.20</td><td>1.21 (0.72-2.01)</td><td>0.472</td></tr><tr><td>Missing</td><td>71</td><td>2 (2.74)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Keeping NHP as pets</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>1,712</td><td>161 (8.60)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>134</td><td>15 (10.07)</td><td>1.19 (0.68-2.08)</td><td>0.54</td><td>1.46 (0.82-2.6)</td><td>0.193</td></tr><tr><td>Missing</td><td>35</td><td>3 (7.89)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Bitten by a NHP</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>1,824</td><td>171 (8.57)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>45</td><td>8 (15.09)</td><td>1.89 (0.87-4.08)</td><td>0.097</td><td>2.31 (1.05-5.11)</td><td>0.038*</td></tr><tr><td>Missing</td><td>12</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Butchering</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>264</td><td>26 (8.97)</td><td>1</td><td></td><td></td><td></td></tr><tr><td>Yes</td><td>1,298</td><td>124 (8.72)</td><td>0.99 (0.99-1.0)</td><td>0.795</td><td>1.00 (0.99-1.01)</td><td>0.674</td></tr><tr><td>Missing</td><td>319</td><td>29 (8.33)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td><strong>Hospitalization</strong></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>No</td><td>699</td><td>62 (8.15)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>1,175</td><td>116 (8.99)</td><td>1.11 (0.81-1.54)</td><td>0.5148</td><td>0.94 (0.68-1.31)</td><td>0.730</td></tr><tr><td>Missing</td><td>7</td><td>1 (12.50)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

(Continued)
<table><caption>Table 2. (Continued)</caption><thead><tr><th rowspan="2"></th><th colspan="2">HTLV-1 results</th><th rowspan="2">Crude OR (95% CI)</th><th rowspan="2">P-value</th><th rowspan="2">Adjusted OR<sup>a</sup> (95% CI)</th><th rowspan="2">P-value</th></tr><tr><th>Negative</th><th>Positive (%)</th></tr></thead><tbody><tr><td colspan="7"><strong>Number of hospitalization</strong></td></tr><tr><td>0</td><td>575</td><td>41 (6.66)</td><td>1</td><td>0.005**</td><td>1</td><td>0.109</td></tr><tr><td>1</td><td>537</td><td>43 (7.41)</td><td>1.12 (0.72-1.75)</td><td></td><td>1.02 (0.65-1.61)</td><td></td></tr><tr><td>2-5</td><td>547</td><td>64 (10.47)</td><td>1.64 (1.09-2.47)</td><td></td><td>1.30 (0.86-1.99)</td><td></td></tr><tr><td>&gt;5</td><td>71</td><td>15 (17.44)</td><td>2.96 (1.56-5.62)</td><td></td><td>2.27 (1.18-4.37)</td><td></td></tr><tr><td>Missing</td><td>151</td><td>16 (9.58)</td><td>1.49 (0.81-2.72)</td><td></td><td>1.36 (0.74-2.51)</td><td></td></tr><tr><td colspan="7"><strong>Transfusion</strong></td></tr><tr><td>No</td><td>1,658</td><td>148 (8.19)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>209</td><td>30 (12.55)</td><td>1.61 (1.06-2.44)</td><td>0.026*</td><td>1.35 (0.88-2.07)</td><td>0.165</td></tr><tr><td>Missing</td><td>14</td><td>1 (6.67)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="7"><strong>Circumcision (among men)</strong></td></tr><tr><td>No</td><td>41</td><td>3 (6.82)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>1,079</td><td>79 (6.82)</td><td>1.00 (0.30-3.30)</td><td>0.999</td><td>1.19 (0.35-3.95)</td><td>0.777</td></tr><tr><td>Missing</td><td>3</td><td>0 (0)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="7"><strong>Scarification</strong></td></tr><tr><td>No</td><td>682</td><td>50 (6.83)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>1,180</td><td>126 (9.65)</td><td>1.46 (1.04-2.05)</td><td>0.031*</td><td>1.23 (0.87-1.75)</td><td>0.242</td></tr><tr><td>Missing</td><td>19</td><td>3 (13.64)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td colspan="7"><strong>Tattoo</strong></td></tr><tr><td>No</td><td>1,553</td><td>154 (9.02)</td><td>1</td><td></td><td>1</td><td></td></tr><tr><td>Yes</td><td>306</td><td>21 (6.42)</td><td>0.69 (0.43-1.11)</td><td>0.127</td><td>0.96 (0.59-1.57)</td><td>0.880</td></tr><tr><td>Missing</td><td>22</td><td>4 (15.38)</td><td>-</td><td>-</td><td>-</td><td>-</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0006832.t002

but not with tattoo nor circumcision. However, when adjusted for age and sex, only individuals with 5 hospitalizations or more were at increased risk of being HTLV-1 infected (Table 2).

Since our study was conducted in primary tropical forest areas and focused on individuals potentially at risk of acquiring HTLV-1 through interspecies transmission, we also obtained information on contacts with NHPs. No association was observed with hunting and butchering activities or keeping NHPs as pets. However, a higher risk was observed among individuals who reported having been bitten by a NHP. Indeed, 8 of 53 individuals (15.1%) were HTLV-1 seropositive (adjusted OR 2.31; 95%CI 1.05–5.11, Table 2).

In the final multivariable analysis model (Table 3), we found an increased independent risk of HTLV-1 infection in women (OR 1.61; 95%CI 1.16–2.23), the elderly people (e.g. >70; OR 3.47; 95%CI 1.89–6.37), persons having a history of multiple hospitalizations (more than 5 times, OR 2.38; 95%CI 1.21–4.66), and individuals belonging to the Pygmy ethnic group (OR 1.93; 95%CI 1.18–3.14). In addition, a NHP bite appears to be marginally associated with a higher risk of HTLV-1 infection (OR 2.16; 95%CI 0.97–4.8).

Discussion

This study represents one of the largest epidemiological study conducted in Central Africa. We show that HTLV-1 is highly endemic among villagers living in the different provinces of

<footnote>*p<0.05</footnote>

<footnote>**p<0.01</footnote>

<footnote>***p<0.001</footnote>

<footnote>a Adjusted OR on age and sex</footnote>
<figure><img src="image_8.png" /><figcaption>**Fig 3.** Prevalence of HTLV-1 according to the geographical distribution of the six provinces studied in Gabon. The data given in the figure indicate, for each province, the HTLV-1 prevalence, the number of positive individuals among those tested, the mean age and the sex ratio. This map was modified from https://commons.wikimedia.org/wiki/Atlas_of_Gabon using Inkscape 0.92.2.</figcaption></figure>

https://doi.org/10.1371/journal.pntd.0006832.g003

the primary tropical forest of Gabon. These data confirm the high prevalence of HTLV-1 in Gabon and extend it to most rural areas, since we studied adult populations of both sexes, from 6 of the 9 rural provinces of this country. Previous studies mainly carried out in pregnant women or in populations from some regions of Gabon (mostly Franceville area and/or Eastern Gabon regions), have observed a high level of HTLV-1 infection [18–21, 24, 25]. However, such studies were mostly based on serological analyses with, for the oldest studies, Western
<table><caption>Table 3. Multivariable analysis of risk factors associated with HTLV-1 infection in the study population.</caption><thead><tr><th>Risk factors</th><th>OR (95% CI)</th><th>P-value</th></tr></thead><tbody><tr><td><strong>Gender</strong></td><td></td><td></td></tr><tr><td>Male</td><td>1</td><td rowspan="2">0.004**</td></tr><tr><td>Female</td><td>1.61 (1.16-2.23)</td></tr><tr><td><strong>Age group</strong></td><td></td><td></td></tr><tr><td>15-35</td><td>1</td><td rowspan="6">&lt;0.001***</td></tr><tr><td>31-40</td><td>1.28 (0.64-2.54)</td></tr><tr><td>41-50</td><td>3.03 (1.66-5.53)</td></tr><tr><td>51-60</td><td>2.25 (1.17-4.31)</td></tr><tr><td>61-70</td><td>3.75 (2.04-6.91)</td></tr><tr><td>&gt;70</td><td>3.47 (1.89-6.37)</td></tr><tr><td><strong>Number of hospitalization</strong></td><td></td><td></td></tr><tr><td>0</td><td>1</td><td rowspan="4">0.05</td></tr><tr><td>1</td><td>1.03 (0.66-1.6)</td></tr><tr><td>2-5</td><td>1.35 (0.87-2.1)</td></tr><tr><td>&gt;5</td><td>2.38 (1.21-4.66)</td></tr><tr><td><strong>Ethnic Group</strong></td><td></td><td></td></tr><tr><td>Bantu</td><td>1</td><td rowspan="2">0.008**</td></tr><tr><td>Pygmies</td><td>1.93 (1.18-3.14)</td></tr><tr><td><strong>Bitten by a NHP</strong></td><td></td><td></td></tr><tr><td>No NHP Bite</td><td>1</td><td rowspan="2">0.06</td></tr><tr><td>NHP Bite</td><td>2.16 (0.97-4.8)</td></tr></tbody></table>

https://doi.org/10.1371/journal.pntd.0006832.t003

blot interpretation criteria that are now considered as poorly specific. The results of our paper confirm the low specificity of the ELISA tests used for HTLV-1 serological screening. Furthermore, most of the indeterminate WB profiles (89.4% in our study) were PCR negative in DNA extracted from blood. These findings are frequently observed in central Africa where some of these false positives correspond to non-specific cross-reactivities, particularly between HTLV-1 Gag proteins and antigens derived from malaria [26]. Here, on the basis of state-of-the-art serological and molecular assays, we confirm a high level of retroviral infection. Surprisingly, this prevalence (8.7%) appears to be higher than in neighboring countries. Indeed, if we analyze the results obtained on similar populations (pregnant women or Bantu adults from rural areas of the rainforest), and using comparable assays, it is obvious that the overall level of HTLV-1 prevalence is at least 2 to 3 times higher in Gabon than in South Cameroun, in Congo or in Central African Republic [27–29]. The reasons for these differences remain to be elucidated. Of note, the existence of high prevalence foci close to low or very low endemic areas is a hallmark of HTLV-1 epidemiology. In South Japan, areas of low HTLV-1 prevalence can be very close to highly endemic areas, especially in rural regions. Furthermore, in a previous study performed in a geographically limited area of Eastern Gabon, the HTLV-1 seroprevalence rates greatly varies between villages, from 5 to 17.8% [21].

In our study, we observed a significant increase in HTLV-1 infection with age, as well as a higher prevalence in women. Such findings are well-known major characteristics of HTLV-1 epidemiology, regardless the endemic level and the studied areas. For instance, this has already been observed in rural populations of Guinea Bissau [15, 17], as well as in patients hospitalized in Cameroon [28] or in Benin [30]. Lastly, this increased prevalence in women and with age is

<footnote>**p<0.01</footnote>

<footnote>***p<0.001</footnote>
commonly found in various populations, such as inhabitants of Salvador de Bahia in Brazil, workers in Jamaica [31] or urban inhabitants of Southern Japan [32]. Such consistency indicates that this is a major specificity of this virus. It has been considered that cumulative sexual transmission over the years, predominantly occurring from male to female, may explain this specificity [33]. In some cases, a cohort effect can also be partly responsible for such an increase with age, as demonstrated in blood donors in Japan [34]. Surprisingly, a higher HTLV-1 seroprevalence in women was not observed in two previous studies performed in Gabon [19, 35]. These discrepancies might be related to the age of the studied population. Indeed, significant sex-specific differences in HTLV-1 prevalence are not found in young adult age groups [21, 32, 33, 36].

Despite the small number of Pygmy individuals included in this study as compared to Bantu (261 vs 1,797), the prevalence of HTLV-1 is significantly higher in Pygmies. Besides, in our multivariate analysis, being a Pygmy was associated with an increased independent risk of HTLV-1 infection (OR 1.93; 95%CI 1.18–3.34). Such findings were previously observed for HTLV-1 and also for HTLV-2 infection in Bakola's Pygmies living in the Ocean region, located in the Western part of South Cameroon [37]. The reasons why these retroviral infections have a higher endemic level in Pygmy populations remain unclear. It has been hypothesized that it may be linked to a founder effect with persistent high level of viral transmission. This may either be due to specific intrafamilial activities (for example prolonged breastfeeding or breastfeeding of several children by the same woman...), or to behavioral factors such as initiation rites, leading to blood contamination by non-sterile needles or sharp objects such as razor blades. Some of these hypotheses have also been proposed for other populations highly endemic for HTLV-1 or HTLV-2, such as Indigenous Australians living in Central Australia [38] or Amerindians from South American countries [39].

Concerning the distribution of HTLV-1 prevalence in the different provinces of Gabon, our data indicate a heterogeneous distribution ranging from 5.5% to 14%. Contrary to what was previously suggested, we have not observed a gradient of HTLV-1 prevalence following the North-South axis. Such discrepancy could be related to various factors including the non-homogeneous and therefore non-comparable structure of the populations previously studied (urban vs rural). In a previous study, the HTLV-1 age-adjusted prevalence rate was significantly higher in rural areas of Gabon than in urban areas [19]. Moreover, the gradient was especially suggested when looking at the population of pregnant women living in urban areas [24].

In this study, 15% of inhabitants bitten by a NHP were HTLV-1-infected compared to only 8.5% for people not bitten by such primates. In addition, NHP bite appears to be marginally associated with a high risk of HTLV-1 infection in our multivariable analysis. These results are reminiscent of those recently published that demonstrated that the bite by a monkey is a risk factor for HTLV-1 infection in a population of hunters living in tropical forest villages of Southern Cameroon [12]. They also echo the recent case-report of a young girl from a village in Eastern Gabon, which showed likely interspecies transmission of STLV-1 upon a severe bite by a *C. nictitans* [40]. Lastly, several publications also indicate that simian foamy viruses (belonging to another family of simian retroviruses) are highly endemic in Central African monkeys and apes, and can easily be transmitted to humans by severe bites from a NHP, leading to chronic infection [41].

In our study, one other risk factor associated with HTLV-1 is hospitalization. Indeed, HTLV-1 prevalence increases from 6% to 17% with the number of hospital stays. In the multivariable analysis, there is an increased independent risk of HTLV-1 infection in individuals with a history of multiple hospitalizations (more than 5 times). The reasons for the existence of such a risk are not clear. It could be linked to the nosocomial acquisition of this retroviral
infection in hospital, either by transfusion of cells infected with HTLV-1, or by contamination
through the use of infected syringes or non-sterile utensils. However, transfusion is not associ-
ated with HTLV-1 infection in our model. One may suppose that HTLV-1 is not actually
transmitted at the hospital. For instance, the virus could be acquired through severe bite,
which could secondarily lead to hospitalization. Similarly, a given illness could lead the patient
to first visit a traditional healer with risky practices (e.g. use of non sterile material), before
heading to more classical hospitalization.

We hoped that this study would alert Gabonese public health authorities and promote
appropriate measures to prevent HTLV-1 infection throughout the country, at least in rural
populations in which we reported high prevalence. Despite the fact that Gabon is both a highly
HTLV-1 endemic area and a high development index country, it has not yet introduced rou-
tine screening for HTLV-1, especially among blood donors. Importantly, a previous study in
Gabon has already demonstrated that transfusion is a major risk factor for HTLV-1 acquisition
in children [16]. In addition, several studies around the world have clearly demonstrated that
HAM/TSP is frequently linked to HTLV-1 acquisition of blood by transfusion of contami-
nated cells [42]. Routine testing of blood donors seems the most necessary and easiest measure
to implement quickly. In this context, we are currently conducting a survey, in close collabora-
tion with the Libreville blood bank, in order to get new insights in the prevalence of HTLV-1
among blood donors. Another control measure that could be implemented is the screening for
HTLV-1 infection in pregnant women. However, this raises the delicate and complex issue of
decreasing or stopping breastfeeding in HTLV-1 infected women in a context of resource-lim-
ited countries. Finally, public health information campaigns on the zoonotic risks associated
with hunting monkeys are also necessary.

**Supporting information**

S1 Checklist. STROBE checklist for cross-sectional study.

(DOC)

Acknowledgments

We acknowledge the “Centre international de Recherches Médicales de Franceville” in Gabon,
which provided logistic assistance for field missions.

We thank Celine Curis and Philippe V. Afonso for critical review of the manuscript.

**Author Contributions**

**Conceptualization:** Delia Doreen Djuicy, Augustin Mouinga-Ondémé, Antoine Gessain.

**Formal analysis:** Delia Doreen Djuicy, Olivier Cassar, Jill-Léa Ramassamy, Arnaud Fontanet.

**Funding acquisition:** Antoine Gessain.

**Investigation:** Delia Doreen Djuicy, Augustin Mouinga-Ondémé, Antony Idam Mamimandjiami, Rodrigue Bikangui, Antoine Gessain.

**Methodology:** Delia Doreen Djuicy, Augustin Mouinga-Ondémé, Olivier Cassar, Antoine Gessain.

**Project administration:** Antoine Gessain.

**Resources:** Antoine Gessain.

**Software:** Delia Doreen Djuicy.
**Supervision:** Augustin Mouinga-Ondémé, Olivier Cassar, Antoine Gessain.

**Validation:** Delia Doreen Djuicy.

**Writing – original draft:** Delia Doreen Djuicy, Olivier Cassar, Antoine Gessain.

**Writing – review & editing:** Delia Doreen Djuicy, Augustin Mouinga-Ondémé, Olivier Cassar, Jill-Léa Ramassamy, Arnaud Fontanet, Antoine Gessain.

## References

1. Poiesz BJ, Ruscetti FW, Gazdar AF, Bunn PA, Minna JD, Gallo RC. Detection and isolation of type C retrovirus particles from fresh and cultured lymphocytes of a patient with cutaneous T-cell lymphoma. Proc Natl Acad Sci U S A. 1980; 77(12):7415-9. Epub 1980/12/01. PMID: 6261256; PubMed Central PMCID: PMCPMC350514.

2. Gessain A, Barin F, Vernant JC, Gout O, Maurs L, Calender A, et al. Antibodies to human T-lymphotropic virus type-I in patients with tropical spastic paraparesis. Lancet. 1985; 2(8452):407-10. Epub 1985/08/24. PMID: 2863442.

3. Osame M, Izumo S, Igata A, Matsumoto M, Matsumoto T, Sonoda S, et al. Blood transfusion and HTLV-I associated myelopathy. Lancet. 1986; 2(8498):104-5. Epub 1986/07/12. PMID: 2873363.

4. Takatsuki T. Adult T-cell leukemia in Japan. In: Seno S TS, Irino S, editor. Topics in Hematology. Amsterdam: Excerpta Medica; 1977. p. 73-7.

5. Gessain A, Cassar O. Epidemiological Aspects and World Distribution of HTLV-1 Infection. Front Microbiol. 2012; 3:388. Epub 2012/11/20. https://doi.org/10.3389/fmicb.2012.00388 PMID: 23162541; PubMed Central PMCID: PMCPMC3498738.

6. Proietti FA, Carneiro-Proietti AB, Catalan-Soares BC, Murphy EL. Global epidemiology of HTLV-I infection and associated diseases. Oncogene. 2005; 24(39):6058-68. Epub 2005/09/13. https://doi.org/10.1038/sj.onc.1208968 PMID: 16155612.

7. Iwanaga M, Watanabe T, Yamaguchi K. Adult T-cell leukemia: a review of epidemiological evidence. Front Microbiol. 2012; 3:322. Epub 2012/09/14. https://doi.org/10.3389/fmicb.2012.00322 PMID: 22973265; PubMed Central PMCID: PMCPMC3437524.

8. Percher F, Jeannin P, Martin-Latil S, Gessain A, Afonso PV, Vidy-Roche A, et al. Mother-to-Child Transmission of HTLV-1 Epidemiological Aspects, Mechanisms and Determinants of Mother-to-Child Transmission. Viruses. 2016; 8(2). Epub 2016/02/06. https://doi.org/10.3390/v8020040 PMID: 26848683; PubMed Central PMCID: PMCPMC4776195.

9. Roucoux DF, Wang B, Smith D, Nass CC, Smith J, Hutching ST, et al. A prospective study of sexual transmission of human T lymphotropic virus (HTLV)-I and HTLV-II. J Infect Dis. 2005; 191(9):1490-7. Epub 2005/04/06. https://doi.org/10.1086/429410 PMID: 15809908.

10. Murphy EL. Infection with human T-lymphotropic virus types-1 and -2 (HTLV-1 and -2): Implications for blood transfusion safety. Transfus Clin Biol. 2016; 23(1):13-9. Epub 2016/01/19. https://doi.org/10.1016/j.tracli.2015.12.001 PMID: 26778839; PubMed Central PMCID: PMCPMC5042452.

11. Locatelli S, Peeters M. Cross-species transmission of simian retroviruses: how and why they could lead to the emergence of new diseases in the human population. AIDS (London, England). 2012; 26(6):659-73. Epub 2012/03/24. https://doi.org/10.1097/QAD.0b013e328350fb68 PMID: 22441170.

12. Filippone C, Betsem E, Tortevoye P, Cassar O, Bassot S, Froment A, et al. A Severe Bite From a Non-human Primate Is a Major Risk Factor for HTLV-1 Infection in Hunters From Central Africa. Clin Infect Dis. 2015; 60(11):1667-76. Epub 2015/02/28. https://doi.org/10.1093/cid/civ145 PMID: 25722199.

13. Fox JM, Mutalima N, Molyneux E, Carpenter LM, Taylor GP, Bland M, et al. Seroprevalence of HTLV-1 and HTLV-2 amongst mothers and children in Malawi within the context of a systematic review and meta-analysis of HTLV seroprevalence in Africa. Trop Med Int Health. 2016; 21(3):312-24. Epub 2015/12/25. https://doi.org/10.1111/tmi.12659 PMID: 26700941; PubMed Central PMCID: PMCPMC4819700.

14. Pepin J, Labbe AC, Mамадоу-Yaya F, Mbelesso P, Mbadingai S, Deslandes S, et al. Iatrogenic transmission of human T cell lymphotropic virus type 1 and hepatitis C virus through parenteral treatment and chemoprophylaxis of sleeping sickness in colonial Equatorial Africa. Clin Infect Dis. 2010; 51(7):777-84. Epub 2010/08/26. https://doi.org/10.1086/656232 PMID: 20735238.

15. van Tienen C, van der Loeff MF, Peterson I, Cotten M, Holmgren B, Andersson S, et al. HTLV-1 in rural Guinea-Bissau: prevalence, incidence and a continued association with HIV between 1990 and 2007. Retrovirology. 2010; 7:50. Epub 2010/06/08. https://doi.org/10.1186/1742-4690-7-50 PMID: 20525366; PubMed Central PMCID: PMCPMC2894744.
16. Delaporte E, Peeters M, Bardy JL, Ville Y, Placca L, Bedjabaga I, et al. Blood transfusion as a major risk factor for HTLV-I infection among hospitalized children in Gabon (Equatorial Africa). J Acquir Immune Defic Syndr. 1993; 6(4):424-8. Epub 1993/04/01. PMID: 8455148.

17. Larsen O, Andersson S, da Silva Z, Hedegaard K, Sandstrom A, Naucler A, et al. Prevalences of HTLV-1 infection and associated risk determinants in an urban population in Guinea-Bissau, West Africa. J Acquir Immune Defic Syndr. 2000; 25(2):157-63. Epub 2000/12/05. PMID: 11103046.

18. Bertherat E, Makuwa M, Renaut A, Nabias R, Georges-Courbot MC. HIV-1, HTLV-I, and HTLV-II in a semiurban population in East Gabon. J Acquir Immune Defic Syndr Hum Retrovirol. 1998; 19(4):430-2. Epub 1998/12/02. PMID: 9833756.

19. Delaporte E, Dupont A, Peeters M, Josse R, Merlin M, Schrijvers D, et al. Epidemiology of HTLV-I in Gabon (Western Equatorial Africa). Int J Cancer. 1988; 42(5):687-9. Epub 1988/11/15. PMID: 3182104.

20. Delaporte E, Monplaisir N, Louwagie J, Peeters M, Martin-Prevel Y, Louis JP, et al. Prevalence of HTLV-I and HTLV-II infection in Gabon, Africa: comparison of the serological and PCR results. Int J Cancer. 1991b; 49(3):373-6. Epub 1991/09/30. PMID: 1917135.

21. Le Hesran JY, Delaporte E, Gaudebout C, Trebuck A, Schrijvers D, Josse R, et al. Demographic factors associated with HTLV-1 infection in a Gabonese community. Int J Epidemiol. 1994; 23(4):812-7. Epub 1994/08/01. PMID: 8002196.

22. Mahieux R, Ibrahim F, Mauclere P, Herve V, Michel P, Tekaia F, et al. Molecular epidemiology of 58 new African human T-cell leukemia virus type 1 (HTLV-1) strains: identification of a new and distinct HTLV-1 molecular subtype in Central Africa and in Pygmies. J Virol. 1997; 71(2):1317-33. Epub 1997/02/01. PMID: 8995656; PubMed Central PMCID: PMCPMC191187.

23. van Buuren S, Boshuizen HC, Knook DL. Multiple imputation of missing blood pressure covariates in survival analysis. Stat Med. 1999; 18(6):681-94. Epub 1999/04/16. PMID: 10204197.

24. Etenna SL, Caron M, Besson G, Makuwa M, Gessain A, Mahe A, et al. New insights into prevalence, genetic diversity, and proviral load of human T-cell leukemia virus types 1 and 2 in pregnant women in Gabon in equatorial central Africa. J Clin Microbiol. 2008; 46(11):3607-14. Epub 2008/10/11. https://doi.org/10.1128/JCM.01249-08 PMID: 18845819; PubMed Central PMCID: PMCPMC2576568.

25. Pegha Moukandja I, Ngoungou EB, Lemamy GJ, Bisvigou U, Gessain A, Toure Ndouo FS, et al. Non-malarial infectious diseases of antenatal care in pregnant women in Franceville, Gabon. BMC Pregnancy Childbirth. 2017; 17(1):185. Epub 2017/06/14. https://doi.org/10.1186/s12884-017-1362-0 PMID: 28606185; PubMed Central PMCID: PMCPMC5469160.

26. Mahieux R, Horal P, Mauclere P, Mercereau-Puijalon O, Guillotte M, Meertens L, et al. Human T-cell lymphotropic virus type 1 gag indeterminate western blot patterns in Central Africa: relationship to Plasmodium falciparum infection. J Clin Microbiol. 2000; 38(11):4049-57. Epub 2000/11/04. PMID: 11060067; PubMed Central PMCID: PMCPMC87540.

27. Filippone C, Bassot S, Betsem E, Tortevoye P, Guillotte M, Mercereau-Puijalon O, et al. A new and frequent human T-cell leukemia virus indeterminate Western blot pattern: epidemiological determinants and PCR results in central African inhabitants. J Clin Microbiol. 2012; 50(5):1663-72. Epub 2012/03/10. https://doi.org/10.1128/JCM.06540-11 PMID: 22403426; PubMed Central PMCID: PMCPMC3347141.

28. Mauclere P, Le Hesran JY, Mahieux R, Salla R, Mfoupouendoun J, Abada ET, et al. Demographic, ethnic, and geographic differences between human T cell lymphotropic virus (HTLV) type I-seropositive carriers and persons with HTLV-I Gag-indeterminate Western blots in Central Africa. J Infect Dis. 1997; 176(2):505-9. Epub 1997/08/01. PMID: 9237719.

29. Tuppin P, Makuwa M, Guerma T, Bazabana MM, Loukaka JC, Jeannel D, et al. Low HTLV-I/II seroprevalence in pregnant women in Congo and a geographic cluster of an HTLV-like indeterminate western blot pattern. J Acquir Immune Defic Syndr Hum Retrovirol. 1996b; 11(1):105-7. Epub 1996/01/01. PMID: 8528728.

30. Dumas M, Houinato D, Verdier M, Zohoun T, Josse R, Bonis J, et al. Seroepidemiology of human T-cell lymphotropic virus type I/II in Benin (West Africa). AIDS Res Hum Retroviruses. 1991; 7(5):447-51. Epub 1991/05/01. https://doi.org/10.1089/aid.1991.7.447 PMID: 1873079.

31. Murphy EL, Figueroa JP, Gibbs WN, Holding-Cobham M, Cranston B, Malley K, et al. Human T-lymphotropic virus type I (HTLV-I) seroprevalence in Jamaica. I. Demographic determinants. Am J Epidemiol. 1991; 133(11):1114-24. Epub 1991/06/01. PMID: 2035515.

32. Tajima K, Kamura S, Ito S, Ito M, Nagatomo M, Kinoshita K, et al. Epidemiological features of HTLV-I carriers and incidence of ATL in an ATL-endemic island: a report of the community-based co-operative study in Tsushima, Japan. Int J Cancer. 1987; 40(6):741-6. Epub 1987/12/15. PMID: 2891624.

33. Mueller N, Okayama A, Stuver S, Tachibana N. Findings from the Miyazaki Cohort Study. J Acquir Immune Defic Syndr Hum Retrovirol. 1996; 13 Suppl 1:S2-7. Epub 1996/01/01. PMID: 8797696.
34. Oguma S, Imamura Y, Kusumoto Y, Nishimura Y, Yamaguchi K, Takatsuki K, et al. Stable human T-lymphotropic virus type I carrier rates for 7 years among a teenaged blood donor cohort of 1986 in Kumamoto, Japan. Leuk Res. 1995; 19(8):567-71. Epub 1995/08/01. PMID: 7658703.

35. Delaporte E, Peeters M, Durand JP, Dupont A, Schrijvers D, Bedjabaga L, et al. Seroepidemiological survey of HTLV-I infection among randomized populations of western central African countries. J Acquir Immune Defic Syndr. 1989a; 2(4):410-3. Epub 1989/01/01. PMID: 2754613.

36. Kajiyama W, Kashiwagi S, Nomura H, Ikematsu H, Hayashi J, Ikematsu W. Seroepidemiologic study of antibody to adult T-cell leukemia virus in Okinawa, Japan. Am J Epidemiol. 1986; 123(1):41-7. Epub 1986/01/01. PMID: 3000171.

37. Mauclere P, Afonso PV, Meertens L, Plancoulaine S, Calattini S, Froment A, et al. HTLV-2B strains, similar to those found in several Amerindian tribes, are endemic in central African Bakola Pygmies. J Infect Dis. 2011; 203(9):1316-23. Epub 2011/04/05. https://doi.org/10.1093/infdis/jir031 PMID: 21459818.

38. Einsiedel L, Woodman RJ, Flynn M, Wilson K, Cassar O, Gessain A. Human T-Lymphotropic Virus type 1 infection in an Indigenous Australian population: epidemiological insights from a hospital-based cohort study. BMC Public Health. 2016; 16:787. Epub 2016/08/17. https://doi.org/10.1186/s12889-016-3366-5 PMID: 27526923; PubMed Central PMCID: PMCPMC4986258.

39. Ferrer JF, Esteban E, Dube S, Basombrio MA, Segovia A, Peralta-Ramos M, et al. Endemic infection with human T cell leukemia/lymphoma virus type IIB in Argentinean and Paraguayan Indians: epidemiology and molecular characterization. J Infect Dis. 1996; 174(5):944-53. Epub 1996/11/01. PMID: 8896494.

40. Kazanji M, Mouinga-Ondeme A, Lekana-Douki-Etenna S, Caron M, Makuwa M, Mahieux R, et al. Origin of HTLV-1 in hunters of nonhuman primates in Central Africa. J Infect Dis. 2015; 211(3):361-5. Epub 2014/08/26. https://doi.org/10.1093/infdis/jiu464 PMID: 25147276.

41. Betsem E, Rua R, Tortevoye P, Froment A, Gessain A. Frequent and recent human acquisition of simian foamy viruses through apes' bites in central Africa. PLoS Pathog. 2011; 7(10):e1002306. Epub 2011/11/03. https://doi.org/10.1371/journal.ppat.1002306 PMID: 22046126; PubMed Central PMCID: PMCPMC3203161.

42. Gout O, Baulac M, Gessain A, Semah F, Saal F, Peries J, et al. Rapid development of myelopathy after HTLV-I infection acquired by transfusion during cardiac transplantation. N Engl J Med. 1990; 322 (6):383-8. Epub 1990/02/08. https://doi.org/10.1056/NEJM199002083220607 PMID: 2300089.