Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Scientific Reports volume 13, Article number: 7659 (2023 ) Cite this article Hbv S Antigen
Around ten million people are infected with HTLV-1 worldwide, and 1–4% develop HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP), characterized by an important degeneration of the spinal cord, which can lead to death. Distinct HLA alleles have been associated with either HAM/TSP susceptibility or protection. However, these HLA alleles set may change according to the population studied. Brazil is the second country in the number of HTLV-1-infected people and there are few reports addressing the HLA influence on HTLV-1 infection as well as on disease outcome. The objective of this study was to evaluate the influence of HLA alleles as a risk factor for HAM/TSP and the proviral load (PVL) levels, clinical progression, and death outcomes in an admixed Brazilian population. The HLA-A, -B, -C, and -DRB1 were genotyped in 375 unrelated HTLV-1-infected individuals divided into asymptomatic carriers (AC) (n = 165) and HAM/TSP (n = 210) in a longitudinal cohort from 8 to 22 years of follow-up. Because locus B deviated from Hardy–Weinberg Equilibrium for the study groups, the results represented for HLA-B alleles were inconclusive. The alleles HLA-A*68 and -C*07 were related to HAM/TSP risk in multivariate analysis. The alleles HLA-A*33, and -A*36 were associated with protection against disease progression in HAM/TSP patients, while -C*12, -C*14, and -DRB1*08 were associated with increased risk of death. In the AC group, the presence of, -C*06 and -DRB1*15 alleles influenced an increased PVL, in an adjusted linear regression model, while -A*30, -A*34, -C*06, -C*17 and -DRB1*09 alleles were associated with increased PVL in HAM/TSP group compared to HAM/TSP individuals not carrying these alleles. All these alleles were also related to increased PVL associated with clinical progression outcome. Increased PVL associated with the death outcome was linked to the presence of HLA-A*30. PVL has been associated with HLA, and several alleles were related in AC and HAM/TSP patients with or without interacting with clinical progression outcomes. Understanding the prognostic value of HLA in HAM/TSP pathogenesis can provide important biomarkers tools to improve clinical management and contribute to the discovery of new therapeutic interventions.
The human T cell lymphotropic virus (HTLV-1) is a retrovirus associated with two primary diseases: adult T cell leukemia/lymphoma (ATLL)1 and the HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP)2,3. Approximately 5 to 10 million people are infected with HTLV-1 worldwide4. Japan and Brazil are countries with the highest estimated total numbers of HTLV-1-infected people5. The HTLV-1 prevalence in Brazil ranges from 146.3 to 390.2/100,000 inhabitants4. Although most infected individuals remain asymptomatic, around 0.25–3.8% develop HAM/TSP over time6. This neurological disease is characterized as a slowly progressing spinal cord disease. The leg muscles gradually become weak. The limbs are stiff, the movements are slow, and walking becomes more difficult. Muscle spasms in the legs are common. Other HTLV-1-associated diseases, such as neurogenic bladder7 and ocular manifestations8, can also be present among HTLV-1-positive individuals. These clinical outcomes are associated with increased HTLV-1 proviral loads (PVL)9; however, it is still unclear what determines the high PVL in these patients. The feature that only a small percentage of HTLV-1 carriers will develop some symptoms has not yet been fully elucidated.
The neurological injury in HAM/TSP may be a consequence of an inflammatory reaction triggered by the recognition of infected cells by cytotoxic T lymphocytes, followed by the release of cytokines and damage to the central nervous system. Factors related to the HTLV-1/host interaction may be involved in the risk of developing HAM/TSP9.
In the search for a biomarker that can identify who is more prone to develop the disease, depending on the population's geographic distribution, some studies have shown that HLA alleles can influence the course of the infection.
The HLA alleles A*31, B*07, B*54, C*07, DRB1*01, and DQB1*05 were described to be related to a higher risk of developing HAM/TSP in Japan10,11,12,13,14. DRB1*01, in the absence of A*02 and C*08, was related to a higher risk of developing HAM/TSP in the Iranian population15,16. In addition, B*07 and DRB1*01 were related to an increased risk of developing HAM/TSP in Spain, where the majority of the individuals studied were from Latin America17. In South America, the allele B*35 was associated with a higher risk of developing HAM/TSP in Argentina18. On the other hand, the HLA-A*02, -C*08, -B*40, -DRB1*15, and -DQB1*06 were related to protection in Japan11,12,13. However, in Iran, Caribbean, Peruvian, and Spain populations, the same HLA alleles were unrelated to HAM/TSP susceptibility or protection15,16,17,19,20,21,22. Only two studies were performed in Brazil, involving a very small sample size of patients with HAM/TSP (n = 9) or selecting alleles previously associated with protection and disease risk in other geographic populations23,24. Unfortunately, these studies did not assess the complete HLA alleles or compare the clinical outcome. In addition, the Brazilian population is highly admixed and a probable association of ethnicity should be considered. Our objective was to describe the genetic polymorphism in the HLA-A, -B, -C, and -DRB1 loci in a subset of the Brazilian population infected with HTLV-1 in a longitudinal follow-up cohort at a Reference Center for Infectious Diseases in Rio de Janeiro, Brazil.
Table 1 contains demographic and clinical data. Three hundred and seventy-five individuals infected with HTLV-1 were randomly selected from the HTLV-1 cohort, corresponding to 41.7% of all HTLV-1 infected individuals followed as outpatients at the INI/Fiocruz. The enrolled participants lived predominantly in the metropolitan area of Rio de Janeiro City. They were divided into HAM/TSP (n = 210) and AC (n = 165) subgroups. The median age of the HAM/TSP group was 61 (59.33–62.55, 95% CI) years, similar to the AC group (61 (54.92–59.82 years, p = 0.14). The median clinical follow-up time was 13 (12.24–13.96, 95% CI), and 17 (14.29–16.57) years for HAM/TSP and AC, respectively. The distribution of gender and ethnicity (a self-identified skin color, described as white or non-white (black and mixed)) was similar in both groups. Also, the HAM/TSP group presented a higher frequency of other clinical manifestations such as the neurogenic bladder (79.5%) and ophthalmic alterations (10.5% of the patients) compared to the AC group, p < 0.01. PVL was higher in the HAM/TSP group compared to the AC group (p < 0.01).
During the clinical follow-up, 24 (11.4%) patients died from secondary causes related to HTLV-1 infection, such as recurrent urinary infections, bedsores, and sepsis. The most common causes of death among HAM/TSP in our cohort were urinary sepsis and deep venous thrombosis associated with pulmonary thromboembolism. All of these individuals were HAM/TSP patients. The clinical progression was investigated during the follow-up, and 148 (70.5%) HAM/TSP and only seven (4.2%) AC patients, showed a deteriorating outcome (p < 0.01). The clinical follow-up in the AC progressors (v.g., from AC to HAM/TSP) subgroup ranged from three to 26 years (median 22 years, 9.98–25.73, 95% CI).
Figure 1 shows the clinical progression throughout 1–27 years (median 14 years, 13.53–15.34, 95% CI) of follow-up in the HAM/TSP progressors subgroup. The clinical progression occurred more rapidly according to the severity of the patient’s condition. Patients with the worse disability, such as those with bilateral walking support (HAM/TSP Bilateral-Progressors), progressed more quickly than those with unilateral walking support (HAM/TSP Unilateral-Progressors) or those able to walk without support (HAM/TSP walking). A Log-Rank Mantel-Cox test showed a statistically significant difference between subgroups (p < 0.0001). The difference in clinical evolution between those individuals with HAM/TSP with unilateral support and bilateral support was significant (p = 0.0462); however, there was no difference in clinical development between individuals who needed unilateral support and those with an independent gait.
Longitudinal clinical progression in HAM/TSP patients related to walking condition ability. A log-Rank Mantel-Cox test evaluated the statistical difference between HAM/TSP Progressors subgroups. The Y-axis represents the percentage of patients without disease progression through the time (in years). The difference between the HAM/TSP-Progresssor with Unilateral support versus HAM/TSP-Progresssor with Bilateral support was statistically significant (p = 0.0462). The comparison of clinical progression between the HAM/TSP-Progressor without walking support and the HAM/TSP-Progresssor with Bilateral support group was also statistically significant (p < 0.0001). No difference was observed when comparing HAM/TSP-Progressor without walking support versus HAM/TSP-Progresssor with Unilateral support patients.
The frequencies of the HLA -A, -B, -C, and -DRB1 alleles in the study population were compatible with those found at the Brazilian Bone Marrow Volunteer Donor Registry (REDOME)25, updated to January 2023, except for the HLA-A*68 allele. From 5,585,608 HLA-typed Brazilians HLA-A*68 represented 6.20% and in the Rio de Janeiro State 6.41% among 254,519 donors, in contrast with 10.37% observed in our study population. Supplemental Table S1 described all allelic frequencies found in our study population.
The B locus deviated from the Hardy–Weinberg Equilibrium (HWE) for both HAM/TSP and AC groups (p < 0.0000 and p < 0.00088, respectively, Supplemental Tables S2 and S3), and for this reason, the results were inconclusive. The remaining loci are in HWE for both groups. After univariate and multivariate analysis, alleles HLA-A*68 and -C*07 showed a higher frequency in the HAM/TSP group than in the AC group, even after being adjusted by age, gender, ethnicity, and time of clinical follow-up (Table 2, OR 2.03 (1.19–3.46, 95% CI), p = 0.01, and OR 1.61 (1.04–2.48, 95% CI), p = 0.03, respectively), and were associated to HAM/TSP risk. As HLA-C*07 in the absence of -A*02 has been previously described as a risk factor for HAM/TSP24, we tested its association in the absence of the -A*02 allele. The contingency table analysis showed that the presence or absence of the -A*02 did not influence the -C*07 allele distribution in both groups.
Due to the small number of affected patients with neurogenic bladder and ophthalmological manifestations was not possible to correlate these data with HLA genotyping.
Predicted three loci haplotypes are represented in Supplemental Tables S4 and S5, for AC and HAM/TSP groups, respectively. In total, 209 and 196 HLA- A ~ C ~ B haplotypes were observed in AC and HAM/TSP patients, respectively. The most common (> 1%) haplotypes were restricted to a minimum of four (A01 ~ C07 ~ B08; A01 ~ C04 ~ B35; A02 ~ C05 ~ B44, A68 ~ C04 ~ B44; A02 ~ C14 ~ B51; A68 ~ C07 ~ B07, A68 ~ C07 ~ B07) copies in AC group, with maximum 11 copies for the A02 ~ C04 ~ B35) haplotype. The most common haplotype in the AC group was observed in only one copy in the HAM/TSP group. Otherwise, the most common haplotype in HAM/TSP (A01 ~ C07 ~ B08) was expressed in five copies in this group and four in the AC group. No differences were observed between both groups.
Logistic regression was performed to analyze whether HLA-A, -C, and -DRB1 alleles were associated with different clinical progression outcomes in HAM/TSP patients. Only seven AC patients progressed to HAM/TSP during the clinical follow-up (17.86 ± 8.51 years), while 148 patients in HAM/TSP group showed clinical progression. Therefore, we only analyzed the association between HLA and disease progression in the HAM/TSP group (Table 3). The alleles -A*33 (OR 0.28 (0.09–0.91 95% CI), p = 0.03) and -A*36 (OR 0.12 (0.02–09 95% CI), p = 0.04) were associated with protection against disease progression in the adjusted model. However, these alleles have a low frequency in our population and these results should be tested in other populations. Supplemental Table S6 shows the frequency of all HLA-A, -B, -C, and -DRB1 loci in HAM/TSP patients related to the disease progression used in this study. Logistic regression was also performed to verify HLA alleles' association with death in the HAM/TSP group. The alleles -C*12 (OR 6.25 (1.71–22.8, CI 95%), p = 0.01), -C*14 (OR 8.85 (1.68–46.67, CI95%), p = 0.01) and -DRB1*08 (OR 9.1 (1.72–48.04, CI 95%), p = 0.01) were associated with death in HAM/TSP patients in the adjusted model (Table 4). Unfortunately, these alleles have a low frequency in our studied population and these results should be tested in other populations to prove their role in HAM/TSP death outcome. Supplemental Table S7 shows the frequency of all HLA-A, -B, -C and -DRB1 loci in HAM/TSP patients related to death used in this study.
Several HLA alleles have been associated with PVL in Japan11,26 but not in other populations. We analyzed whether HLA alleles in our population alter the PVL. We performed an independent case and control analysis because PVL was higher in the HAM/TSP group compared to the AC group (Table 1). Table 5 describes, through linear regression analysis, the role of the HLA alleles in the PVL in AC patients. Individuals in this group, who carry the alleles C*06 (β 3.904 (1.016; 6.792, 95% CI), p = 0.008) and -DRB1*15 (β 3.381 (0.741; 6.021, 95% CI), p = 0.012) had an increased PVL compared to those who do not carry those alleles (β, 95% CI as shown in Table 5), after adjustment. In the univariate analysis, carriers of the alleles -A*30 and -A*66 also presented higher PVL (Table 5). However, after adjusting the multivariate analysis, there was no difference between carriers and not-carriers in the PVL related to these alleles.
In the HAM/TSP group, individuals who carry the A*30 (β 3.033 (0.252;5.8130, 95% CI), p = 0.033), -A*34 (β 9.172 (2.755;15.588, 95% CI), p = 0.005), -C*06 (β 4.016 (0.883;7.15, 95% CI), p = 0.012), -C*17 (β 7.048 (0.696;13.399, 95% CI), p = 0.03) and -DRB1*09 (β 6.874 (2.341;11.408, 95% CI), p = 0.003) alleles had an increased PVL compared to those who do not carry those alleles (β, 95% IC as showed in Table 6). Indeed, these alleles were also related to an increased PVL and a worse clinical progression (Table 6). On the other hand, only the allele -A*30 (β 24.634 (11.282; 37.985, 95% CI) was associated with a higher PVL and death (Table 6). Due to the low frequency of -DR*09 carries, it was not possible to calculate the β value in the multivariate model. None of the alleles -A*34, -C*06, -C*17, and -DRB1*09 were associated with both death outcome and the PVL levels in HAM/TSP patients.
Brazil is a country with continental proportions and endemic to HTLV-1. Its population is composed of miscegenation among individuals of Indigenous, African, Asia, and European origin. Unfortunately, the Brazilian population is underrepresented in public genomic reference databases. Previous data comparing whole exome sequencing and SNP array from the Brazilian population and the 1000 Genomes datasets showed that there are significant differences in allele frequencies for both common and rare variants and this result indicates that none of the admixed American populations present in the 1000 Genomes dataset can be used as a surrogate for studies of the Brazilian population27. This is a limitation of our case–control study. Some studies have associated HLA genetic polymorphism with protection or susceptibility to HAM/TSP development. This study evaluated the influence of HLA class I and II gene polymorphism in HTLV-1 infected individuals in a longitudinal follow-up and its association with PVL levels, clinical progression, and death outcomes. AC and HAM/TSP patients were accompanied for two to a maximum of 27 years, and clinical progression was related to HAM/TSP patients. The lifetime risk of HAM/TSP is different among ethnic groups26. Some studies report that it can take up to 21 years from the onset of infection to using a wheelchair in HAM/TSP28. Our study could not specify the median time patients became wheelchair users because the walking condition was only assessed in the first and the last clinical consultation. However, HAM/TSP patients that presented a clinical progression were in the cohort for at least eight years. Regarding mortality in HAM/TSP patients, it was pointed out that it would be a consequence of complications of the disease itself and co-infections with HCV/HIV29. Although there are some case reports of patients with HAM/TSP developing ATL30, no large-scale prospective cohort studies have been conducted on patients with HAM/TSP, nor has the risk of ATL transformation been evaluated in these patients. None of the patients included in the present study died of ATL.
Unfortunately, both groups from our HTLV-1 study population were in HW disequilibrium for locus HLA-B, in opposition to the studied 2,624 renal transplant candidates from São Paulo State, Brazil31. However, recently a review article discussed the role of HLA-B in 1322 Brazilians compared to different admixed populations worldwide32. These authors found that HLA-B gene diversity was higher in admixed and European populations while was lower in African populations. Studying gene diversity in admixed Brazilian population may increase the dataset in others genes banks around the world. As an example, 23 new alleles of the HLA-B gene were exclusively identified in the Brazilian population32. The high variability in the HLA-B gene must be studied in light of the role of their alleles in infectious diseases, so common in developing countries. Further studies should address the association of the HLA-B locus with HAM/TSP risk in a large admixed population to increase the likelihood of identifying the disease relationship.
In our cohort, PVL and neurogenic bladder were significantly higher in HAM/TSP than in the AC group, reinforcing what has already been published33,34. The close relationship between HTLV-1 infection and uveitis has been identified concerning the ophthalmologic alterations, indicating that HTLV-1 is also highly associated with intraocular inflammatory disorder35,36. Some HLA alleles have been related to uveitis. The HLA-A*29 is strongly associated with Birdshot uveitis in individuals of Western-European ancestry, while HLA-B*27 is associated with uveitis, particularly in ankylosing spondylitis37,38. In our cohort, ophthalmic manifestations were more frequent in the HAM/TSP group, but it was not possible to correlate it with any HLA due to the small number of affected patients.
Only two studies with small sample sizes were carried out in Brazil to investigate the HLA association with HTLV-1 infection23,24. Borducchi et al.23 analyzed a small number of HTLV-1-infected patients (n = 71), including AC, ATLL, and HAM/TSP. No statistical difference was determined among these groups because of the small sample size. This study showed a tendency to a lower HLA-A2 frequency in HAM/TSP white patients and an association of HLA-DR11 with HAM/TSP only in the Mestizo patients, in disagreement with our study even after adjustment by ethnicity. Unfortunately, that report did not mention the number of HAM/TSP patients and did not compare their results using multivariable analysis. In another study24, only alleles previously associated with protection and disease risk were tested in a small cohort of 84 AC and 9 HAM/TSP patients. They identified that the alleles HLA-C*07 in the absence of -A*02 was associated with an increased risk of disease, and alleles C*08 and B*07 had no association with the development of HAM/TSP. In our study population, the presence or absence of A*02 did not influence the C*07 allele, and C*08 also showed no association with the disease.
In our study, the most frequent allele observed at locus A was A*02 (19.02%, Supplemental Table S1). Still, unlike in the Japanese studies11,12, neither A*02 nor the HLA-C*08 alleles were associated with the protection for HAM/TSP development.
Our results showed that the alleles -A*68 and -C*07, had a higher frequency among HAM/TSP patients than in the AC individuals, leading to an increased risk of HAM/TSP. To our knowledge, no previous studies indicated the association of -A*68 with HAM/TSP. However, we should take into account that the -A*68 allele had a higher frequency in our study cohort when compared to the frequency of this allele in the REDOME25. Although differences in the frequency could be found among some ethnic groups (greater in black (7.88%) and native Amerindians (7.32%) compared to white (5.66%) ones), the difference observed in patients with HTLV-1 was much more expressive (10.37%), and even greater in HAM/TSP patients (13.29%), suggesting a possible increased risk of this allele for HTLV-1 infection and disease progression in the Brazilian population. Indeed, in an HIV-positive African population, the allele -A*68 was associated with the rapid transmission of HIV in serum discordant heterosexual partners, with 21.3% of frequency39 and was associated with a protective role to SARS-COV-2 severe infection and death in a Mexican population, where this allele has a high frequency (12.6% in patients with COVID-19 and 28.6% in the healthy control group)40. HLA-A*68 is more frequently identified in South American indigenous populations, such as Argentina, Paraguay, and Chile41,42,43. Besides REDOME data there is a lack of reports describing the frequency of the HLA-A*68 allele in Brazil or in infectious diseases. More studies addressing the role of HLA-A*68 in the Brazilian population in the context of HTLV-1 infection might confirm his association with the augmented risk for infection and disease progression as observed in our study. The -C*07 allele was cited in the Japanese population as a risk factor for HAM/TSP13, as well as in our study. Interestedly, in an Argentine population, -C*07 was associated with an increased risk of ATLL but not with HAM/TSP18. This allele was underrepresented in AIDS Brazilian patients conferring protection against cytomegalovirus retinitis44.
This is the first study that analyzed whether HLA alleles could influence clinical progression in a longitudinal cohort. Patients were clinically followed for 20 years, and their walking clinical progression condition was determined. In our cohort -A*33 and -A*36 were associated with clinical progression and to our knowledge these alleles were not described in other populations. The role of these alleles should be clarified in other populations concerning their physiopathological properties.
There is a clear association between HLA alleles and PVL levels17,18,19. We identified new alleles that were associated with increased PVL in the AC group, carriers of the alleles -C*06, and -DRB1*15, and in the HAM/TSP group, carriers of the alleles -A*30, -A*34, -C*06, -C*17 and -DRB1*09. Although the presence of these alleles influenced an increased PVL, among these alleles, none of them was previously reported. Only B*40 was cited in a Japanese population study13, where its distribution was associated with a protective effect. In our cohort HLA-B alleles were in disequilibrium of HW in both cases and control groups.
Further, we analyzed if the presence of the HLA alleles were responsible for the increase in PVL regarding the clinical progression outcome. In the AC group, carriers of the alleles: -C*06, and -DRB1*15 had an increased PVL. In the HAM/TSP group, the alleles -A*30, -A*34, -C*06, -C*17, and -DRB1*09 also influenced PVL and risk for clinical progression. Intriguingly, the alleles that influence PVL and clinical progression in the AC group are not the same as in the HAM/TSP group. Among these alleles, the only correlation with other studies we found in the -DRB1*15 allele, which in the Japanese study had a protective effect13 and is cited with a tendency to a higher frequency in HAM/TSP patients than in HCs in a French Afro-Caribbean population22, but without association with PVL. It is the first time these associations are described in the Brazilian population.
Moreover, we can highlight that the presence of the HLA-C*06 allele is associated with PVL, clinical progression, and death outcomes for the first time in our population. Studies carried out in ethnically different populations have not demonstrated equanimity in their results regarding risk or protection for the development of HAM/TSP. The host's genetic background must have an important role in this issue.
HTLV-1 has remarkable genetic stability45 and two previous reports from Brazilian cities (São Paulo46 and Salvador47) showed no association with sequenced viral genotype and clinical outcomes. In Brazil, Genotype a (Cosmopolitan) and subgroup A (transcontinental) is the most frequently (> 95%) identified HTLV-1 genotype48. None of our enrolled patients was Japanese or their descendant, in which Japanese subgroup B could be identified.
A major limitation of our study is related to the Brazilian genetic background. Brazilian population, as well as all Latin American countries, was composed through centuries of European colonization and the African slave trade. The continuous admixture events among Europeans, Native Americans, and West-Central and Southeast Africans occurred differently in each region of the country, producing a stratified population, even in the current days. Despite our study population showed similar HLA frequencies compared to REDOME (including more than 5 million Brazilians from several regions), it may not reflect the entire genetic diversity of the Brazilian population. Even though, the Brazilian population has a large proportion of European white population49,50. Furthermore, Brazil is the large population in Latin America and the fifth-largest in the world, and the Brazilian population is underrepresented in public genomic reference databases27. In these databases, 80% of all genetic data is related to the European ancestry population which accounts for only 16% of the world´s population. This gap prevents the discovery of new rare allele variants. Few studies reported the ancestry background of the general Brazilian population51 and studies associating ancestry analysis and HLA frequencies could better relate the role of HLA polymorphisms in the HTLV-1 infection and its medical importance. Furthermore, the study was conducted with clinical data collected retrospectively in a convenience sample. It was impossible to match patients with similar and long clinical follow-up times and with similar disease stages. Another limitation of this study was the reduction of the HLA typing to the level of allelic group. Although we have used in part of the typing a medium-level resolution that could result in better discrimination of the HLA alleles of risk and susceptibility, we were also limited by the HLA frequencies that tend to be small and implying having a larger cohort to have statistical differences. Nonetheless, we are not aware that allelic level should be the final goal for HLA studies, as in many diseases that initially were associated with an antigen (allele group level) that with further studies have pointed out that only some alleles of the group were in association with the disease, as for spondyloarthropathies with HLA-B*2752 or rheumatic arthritis HLA-DRB*0453. Our results show an association between HLA frequencies and HAM/TSP risk, PVL, disease progression, and death. Because of the low number of participants presenting some HLA alleles, the results observed in this work have to be confirmed in other populations. Nevertheless, in light of our acknowledgment, this is the study involving the highest number of participants associating HLA alleles and HTLV-1 infection in the Brazilian population. Future studies, designed with an even greater number of patients, with a long and similar clinical follow-up period, and which include patients in the early stages of the disease, may increase the strength of comparisons between some HLAs regarding their role in HTLV-1 infection in the Brazilian population, reinforcing our findings.
The HLA system is vital in response to infections. In the case of HTLV-1 infection, where a minority of individuals carrying the virus will develop the neurological disease, the genetic determinants could be seen as an interplay between different immunological events. However, more studies should be conducted to monitor the HLA profile in AC patients at risk of developing the disease in longitudinal follow-up.
In this study, we investigated the influence of HLA polymorphisms on the susceptibility for the development of HAM/TSP and clinical progression in a longitudinal cohort of an admixed Brazilian population. We found that HLA-A*68 and -C*07 carriers presented an augmented risk for HAM/TSP development in HTLV-1-infected individuals. In addition, -C*12, -C*14, and -DRB1*08 were associated with an augmented risk of death in HAM/TSP individuals. In contrast, the alleles HLA-A*33, and -A*36 were related to protection against disease progression in these individuals. We also showed that asymptomatic -C*06 and -DRB1*15 carriers have increased proviral load, while -A*30, -A*34, -C*06, -C*17, and -DRB1*09 alleles were associated with this characteristic on HAM/TSP patients. We also demonstrated that patients with high PVL and carrying these alleles presented disease clinical progression, while the presence of HLA-A*30 in these patients was associated with death. This is the first study to evaluate the influence of HLA alleles on the risk for HAM/TSP development as well as in clinical progression and death, in an admixed Brazilian population. Brazil has one of the highest rates of HTLV-1 infection, and studies addressing genetic factors associated with HAM/TSP risk could prevent disease worsening through earlier medical care.
The Laboratory for Clinical Research in Neuroinfections at the Evandro Chagas National Institute of Infectious Diseases (Lapclin-Neuro), FIOCRUZ (INI/Fiocruz), Rio de Janeiro, RJ, Brazil, is a reference center for diagnostic and clinical follow-up of asymptomatic carriers and patients with neurological manifestations of HTLV-1 infection. Since 1999 more than 1000 AC and neurological patients, the majority of them with HAM/TSP, have been investigated and followed up in a longitudinal cohort by a team of neurologists. According to their clinical visit, HTLV-1 positive and genetically unrelated patients were randomly and consecutively included from the longitudinal study cohort. The Research Ethics Committee of the INI/FIOCRUZ (CAAE-0012.0.009.000-08) approved this study, and all subjects provided written informed consent. Volunteers were consecutively followed from 1999 to 2020 in a longitudinal cohort study from eight to 22 years. A sex, age, and ethnicity-matched unrelated case–control study was then designed to identify HLA's influence on the disease disability status, PVL, disease progression, and death outcomes. All clinical data were collected retrospectively, using a convenience sample provided by the cohort follow-up. A comprehensive medical history relating to neurological illnesses was taken at the time of inclusion in the cohort and annually for the AC patients and biannually for the HAM/TSP patients.
The clinical assessment and classification of the neurological cases followed the World Health Organization guidelines and diagnostic criteria for HAM/TSP54. Patients’ inclusion criteria were defined by positive serology for HTLV-1 in the blood (for AC) and blood and cerebrospinal fluid (CSF) for neurological patients (HAM/TSP and other neurological findings). The screening for HTLV-1/2 infection was performed throughout the serological test using a Recombinant ELISA HTLV1/2 kit from Wiener (Rosario, Argentina), and confirmation of HTLV-1 infection through western blot analysis (Cambridge Biotech, Worcester, MA, USA). Patients with HTLV-2 infection were not enrolled in the study. Individuals were divided into two groups: individuals with HAM/TSP, and other isolated neurological findings such as isolated neurogenic bladder, and HTLV-1 asymptomatic carriers (AC). The HAM/TSP participants were classified according to their disability status using the IPEC Clinical Disability Scale (IPEC-CDS) as described elsewhere55. The IPEC disability scale is a numerical scale ranging from 0 (no disability) to 29 (worse neurological status). This scale is a composite of motor, pain, sensory, and sphincter subscores55.
Clinical progression was determined by comparing patients' gait performance according to the IPEC-CDS motor subscore at the first clinical visit, when patients were included in the cohort, with the same score obtained at the last clinical visit. Any worsening in the IPEC-CDS motor score was considered disease progression, regardless of the initial score assigned to patients at the time of study enrollment. We also did not consider the degree of progression, meaning that patients who progressed to two or more levels on the motor scale were not separated from those who progressed to only one level.
A survival curve was generated to analyze the longitudinal clinical progression in HAM/TSP patients, using each of their initial clinical classifications regarding the ability to walk (walk without assistance, with unilateral or bilateral support, and the need for a wheelchair). A change in the motor scale of worsening prognosis was considered a sign of disease progression.
Peripheral blood (5 mL) was collected in EDTA tubes from all individuals, and DNA was extracted using the Puregene commercial kit (Gentra Systems Inc., Minneapolis, MN, USA), following the manufacturer’s instructions.
First, the HTLV-1 cohort was HLA genotyped through an in-house PCR-SSP (polymerase chain reaction-sequence-specific primer) from 1999 to 2007 (n = 273), and throughout PCR-SSO (polymerase chain reaction with specific oligonucleotide probes) technique (n = 102).
The HLA class I loci (-A, -B, and -C) genotyping was performed in an in-house PCR-SSP according to the methodology established by the Organ Transplant Laboratory from Oxford University56. Primers were donated by the Department of Immunology at the St Mary’s School of Medicine, Imperial College, London, UK, and the PCR-SSO technique, using DNA LABType® trading system (One Lambda Inc., Canoga Park, CA, USA) was performed, according to manufacturer`s instructions, that provides SSO probes for the sequence-specific oligonucleotide linked to fluorescent microspheres for the identification of HLA alleles in genomic DNA samples.
The HTLV-1 PVL DNA was measured by real-time PCR assay (SmartCycle II; Cepheid) using the TaqMan system (Applied Biosystems) to amplify a 159 bp fragment of tax gene with SK43 and SK44 primers. A standard curve was generated using the β-globin gene as a reference. DNA from the TARL-2 cell line, which contains a single copy of the provirus HTLV-1, was used to establish the standard curve for tax gene quantification. The proviral load (PVL) in 104 cells was calculated by the following equation: [(copy number of tax gene)/ (copy number of β-globin gene/2)] × 10,00033. The mean of proviral load measures during the follow-up was used to determine the proviral load of each patient.
A descriptive analysis of clinical and demographic characteristics was performed. Outcomes (clinical status—HAM/TSP or AC, disease progression, and death) and alleles (HLA-A, -B, -C, and -DR) were associated using absolute and relative frequencies or mean and standard deviation. Hardy–Weinberg equilibrium was calculated by Exact test using a Markov chain to define the equilibrium in the allelic distribution in HLA-A, -B, -C, and -DR locus between HAM/TSP and AC patients. Prism 7 was used to determine a survival curve to analyze the walking ability's clinical progression. The Mann–Whitney U test evaluated quantitative variables, whereas the univariate Odds Ratio (OR) was used for categorical variables. We used a multivariate generalized linear Binomial model to estimate OR and Ci = Is. Linear models were tested to estimate the average difference in PVL among alleles and outcomes. A p-value < 0.05 was considered significant. All analysis was performed using R software version 3.6.3.
All research was performed following relevant guidelines/regulations described by the policies of the Nature Portfolio journals, and the Declaration of Helsinki. This study was approved by the Research Ethics Committee of the INI/FIOCRUZ (CAAE-0012.0.009.000-08) and all participants provided written informed consent.
All data generated or analyzed during this study are included in this published article [and its supplemental information files].
Uchiyama, T., Yodoi, J., Sagawa, K., Takatsuki, K. & Uchino, H. Adult T-cell leukemia: clinical and hematologic features of 16 cases. Blood 50, 481–492 (1977).
Article CAS PubMed Google Scholar
Gessain, A. et al. Antibodies to human T-lymphotropic virus type-I in patients with tropical spastic paraparesis. Lancet 2, 407–410 (1985).
Article CAS PubMed Google Scholar
Osame, M. et al. HTLV-I associated myelopathy, a new clinical entity. Lancet 1, 1031–1032 (1986).
Article CAS PubMed Google Scholar
Gessain, A. & Cassar, O. Epidemiological aspects and world distribution of HTLV-1 infection. Front. Microbiol. 3, 388 (2012).
Article PubMed PubMed Central Google Scholar
Bangham, C. R. M., Araujo, A., Yamano, Y. & Taylor, G. P. HTLV-1-associated myelopathy/tropical spastic paraparesis. Nat. Rev. Dis. Primers. 1, 15012 (2015).
Yamano, Y. & Sato, T. Clinical pathophysiology of human T-lymphotropic virus-type 1-associated myelopathy/tropical spastic paraparesis. Front. Microbiol. 3, 389 (2012).
Article CAS PubMed PubMed Central Google Scholar
Castro, N. M. et al. Urinary symptoms associated with human T-cell lymphotropic virus type I infection: Evidence of urinary manifestations in large group of HTLV-I carriers. Urology 69, 813–818 (2007).
Rathsam-Pinheiro, R. H. et al. Revisiting Keratoconjunctivitis sicca associated with Human T-Cell Lymphotropic Virus Type 1: prevalence, clinical aspects and proviral load. Braz. J. Infect. Dis. 23, 95–101 (2019).
Article PubMed PubMed Central Google Scholar
Ijichi, S. et al. An autoaggressive process against bystander tissues in HTLV-I-infected individuals: a possible pathomechanism of HAM/TSP. Med. Hypotheses. 41, 542–547 (1993).
Article CAS PubMed Google Scholar
Nishimura, Y. et al. A possible association between HLA and HTLV-I-associated myelopathy (HAM) in Japanese. Tissue Antigens 37, 230–231 (1991).
Article CAS PubMed Google Scholar
Jeffery, K. J. et al. HLA alleles determine human T-lymphotropic virus-I (HTLV-I) proviral load and the risk of HTLV-I-associated myelopathy. Proc. Natl. Acad. Sci. U S A. 96, 3848–3853 (1999).
Article CAS PubMed PubMed Central ADS Google Scholar
Jeffery, K. J. et al. The influence of HLA class I alleles and heterozygosity on the outcome of human T cell lymphotropic virus type I infection. J. Immunol. 165, 7278–7284 (2000).
Article CAS PubMed Google Scholar
Penova, M. et al. Genome wide association study of HTLV-1-associated myelopathy/tropical spastic paraparesis in the Japanese population. Proc. Natl. Acad. Sci. U S A. 118, e2004199118 (2021).
Article CAS PubMed PubMed Central Google Scholar
Kitze, B. et al. Human CD4+ T lymphocytes recognize a highly conserved epitope of human T lymphotropic virus type 1 (HTLV-1) env gp21 restricted by HLA DRB1*0101. Clin. Exp. Immunol. 111, 278–285 (1998).
Article CAS PubMed PubMed Central Google Scholar
Sabouri, A. H. et al. Differences in viral and host genetic risk factors for development of human T-cell lymphotropic virus type 1 (HTLV-1)-associated myelopathy/tropical spastic paraparesis between Iranian and Japanese HTLV-1-infected individuals. J. Gen. Virol. 86, 773–781 (2005).
Article CAS PubMed Google Scholar
Rafatpanah, H. et al. Association between HLA-DRB1*01 and HLA-Cw*08 and outcome following HTLV-I infection. Iran J. Immunol. 4, 94–100 (2007).
Treviño, A. et al. Association between HLA alleles and HAM/TSP in individuals infected with HTLV-1. J. Neurol. 260, 2551–2555 (2013).
Benencio, P. et al. HLA-B*35 as a new marker for susceptibility to human T-cell lymphotropic virus type 1 (HTLV-1) Associated Myelopathy/Tropical Spastic Paraparesis (HAM/TSP) in patients living in Argentina. Retrovirology 17, 29 (2020).
Article CAS PubMed PubMed Central Google Scholar
Taghaddosi, M. et al. Association between HLA Class I Alleles and Proviral Load in HTLV-I associated myelopathy/tropical spastic paraperesis (HAM/TSP) patients in Iranian population. Iran. J. Basic Med. Sci. 16, 264–267 (2013).
PubMed PubMed Central Google Scholar
Tarokhian, H. et al. The effect of HTLV-1 virulence factors (HBZ, Tax, proviral load), HLA class I and plasma neopterin on manifestation of HTLV-1 associated myelopathy tropical spastic paraparesis. Virus Res. 228, 1–6 (2017).
Article CAS PubMed Google Scholar
Talledo, M. et al. Evaluation of host genetic and viral factors as surrogate markers for HTLV-1-associated myelopathy/tropical spastic paraparesis in Peruvian HTLV-1-infected patients. J. Med. Virol. 82, 460–466 (2010).
Article CAS PubMed Google Scholar
Deschamps, R. et al. Absence of consistent association between human leukocyte antigen-I and -II alleles and human T-lymphotropic virus type 1 (HTLV-1)-associated myelopathy/tropical spastic paraparesis risk in an HTLV-1 French Afro-Caribbean population. Int. J. Infect. Dis. 14, e986-990 (2010).
Borducchi, D. M. et al. Human leucocyte antigen and human T-cell lymphotropic virus type 1 associated diseases in Brazil. Br. J. Haematol. 123, 954–955 (2003).
Catalan-Soares, B. C. et al. HLA class I alleles in HTLV-1-associated myelopathy and asymptomatic carriers from the Brazilian cohort GIPH. Med. Microbiol. Immunol. 198, 1–3 (2009).
Article CAS PubMed Google Scholar
Genomic Profile of REDOME/REREME [Internet].Brazilian Immunogenetics Network.[cited 2022 Feb 17].Available from: http://imunogenetica.org/resultados/perfil-genomico-do-redome-rereme/.
Saito, M. Association between HTLV-1 genotypes and risk of HAM/TSP. Front. Microbiol. 10, 1101 (2019).
Article PubMed PubMed Central Google Scholar
Rocha, C. S., Secolin, R., Rodrigues, M. R., Carvalho, B. S. & Lopes-Cendes, I. The Brazilian Initiative on Precision Medicine (BIPMed): Fostering genomic data-sharing of underrepresented populations. NPJ Genom. Med. 5, 42 (2020).
Article PubMed PubMed Central Google Scholar
Araujo, A. Q. C. Neurologic complications of HTLV-1: A review. Rev. Bras. Neurol. 1, 5–10 (2019).
Marcusso, R. M. N. et al. Dichotomy in fatal outcomes in a large cohort of people living with HTLV-1 in São Paulo Brazil. Pathogens 9, E25 (2019).
Martin, F., Fedina, A., Youshya, S. & Taylor, G. P. A 15-year prospective longitudinal study of disease progression in patients with HTLV-1 associated myelopathy in the UK. J. Neurol. Neurosurg. Psychiatry. 81, 1336–1340 (2010).
Ravazzi-Gauch, C., Bajay, M. M., Caldas, H. C. & Abbud-Filho, M. HLA-A, -B, and -DRB1 allele and haplotype diversity in a cohort of Brazilian renal transplant candidates. Hum. Immunol. 77, 464–469 (2016).
Article CAS PubMed Google Scholar
Silva, N.D.S.B., Souza, A.S., Andrade, H.S., Pereira, R.N., Castro, C.F.B., & Vince, N., et al. Immunogenetics of HLA-B: SNP, allele, and haplotype diversity in populations from different continents and ancestry backgrounds. HLA (2023).
Silva, M. T. T. et al. Human T lymphotropic virus type 1 (HTLV-1) proviral load in asymptomatic carriers, HTLV-1-associated myelopathy/tropical spastic paraparesis, and other neurological abnormalities associated with HTLV-1 infection. Clin. Infect. Dis. 44, 689–692 (2007).
Silva, M. T. et al. Isolated bladder dysfunction in human T lymphotropic virus type 1 infection. Clin. Infect. Dis. 48, e34-36 (2009).
Mochizuki, M. et al. Uveitis associated with human T lymphotropic virus type I: seroepidemiologic, clinical, and virologic studies. J. Infect. Dis. 166, 943–944 (1992).
Article CAS PubMed Google Scholar
Kamoi, K. HTLV-1 in ophthalmology. Front. Microbiol. 11, 388 (2020).
Article PubMed PubMed Central Google Scholar
Kuiper, J. J. W. & Venema, W. J. HLA-A29 and Birdshot Uveitis: Further Down The Rabbit Hole. Front. Immunol. 11, 599558 (2020).
Article CAS PubMed PubMed Central Google Scholar
İnanç, M., Şimşek, M. & Çakar Özdal, M. P. Etiological and Clinical Characteristics of HLA-B27-associated Uveitis in a Tertiary Referral Center. Turk. J. Ophthalmol. 49, 10–14 (2019).
Article PubMed PubMed Central Google Scholar
Song, W. et al. Disparate associations of HLA class I markers with HIV-1 acquisition and control of viremia in an African population. PLoS ONE 6, e23469 (2011).
Article CAS PubMed PubMed Central ADS Google Scholar
Hernández-Doño, S. et al. Protective HLA alleles against severe COVID-19: HLA-A*68 as an ancestral protection allele in Tapachula-Chiapas Mexico. Clin. Immunol. 238, 108990 (2022).
Article PubMed PubMed Central Google Scholar
Rey, D. et al. HLA genetic profile of Mapuche (Araucanian) Amerindians from Chile. Mol. Biol. Rep. 40, 4257–4267 (2013).
Article CAS PubMed Google Scholar
Layrisse, Z. et al. Extended HLA haplotypes in a Carib Amerindian population: the Yucpa of the Perija Range. Hum. Immunol. 62, 992–1000 (2001).
Article CAS PubMed Google Scholar
Single, R. M. et al. Demographic history and selection at HLA loci in Native Americans. PLoS ONE 15, e0241282 (2020).
Article CAS PubMed PubMed Central Google Scholar
Biberg-Salum, T. G. et al. HLA-C alleles and cytomegalovirus retinitis in Brazilian patients with AIDS. J. Ophthalmol. 1, 3830104 (2018).
Gessain, A., Gallo, R. C. & Franchini, G. Low degree of human T-cell leukemia/lymphoma virus type I genetic drift in vivo as a means of monitoring viral transmission and movement of ancient human populations. J. Virol. 66, 2288–2295 (1992).
Article CAS PubMed PubMed Central Google Scholar
Pessôa, R. et al. Molecular characterization of human T-cell lymphotropic virus type 1 full and partial genomes by Illumina massively parallel sequencing technology. PLoS ONE 9, e93374 (2014).
Article PubMed PubMed Central ADS Google Scholar
Araújo, T. H. A. et al. Complete genome sequence of human T-cell lymphotropic type 1 from patients with different clinical profiles, including infective dermatitis. Infect. Genet. Evol. 79, 104166 (2020).
Afonso, P. V., Cassar, O. & Gessain, A. Molecular epidemiology, genetic variability and evolution of HTLV-1 with special emphasis on African genotypes. Retrovirology 16, 39 (2019).
Article CAS PubMed PubMed Central Google Scholar
Kehdy, F. S. G. et al. Origin and dynamics of admixture in Brazilians and its effect on the pattern of deleterious mutations. Proc. Natl. Acad. Sci. U S A. 112, 8696–8701 (2015).
Article CAS PubMed PubMed Central ADS Google Scholar
Secolin, R. et al. Distribution of local ancestry and evidence of adaptation in admixed populations. Sci. Rep. 9, 13900 (2019).
Article PubMed PubMed Central ADS Google Scholar
De Oliveira, T. C., Secolin, R. & Lopes-Cendes, I. A review of ancestrality and admixture in Latin America and the caribbean focusing on native American and African descendant populations. Front Genet. 14, 1091269 (2023).
Article PubMed PubMed Central Google Scholar
Fahed, H., Mauro, D., Ciccia, F. & Ziade, N. R. What does human leukocyte antigen B27 have to do with spondyloarthritis?. Rheum Dis Clin North Am. 46, 225–239 (2020).
Becart, S., Whittington, K. B., Prislovsky, A., Rao, N. L. & Rosloniec, E. F. The role of posttranslational modifications in generating neo-epitopes that bind to rheumatoid arthritis-associated HLA-DR alleles and promote autoimmune T cell responses. PLoS ONE 16, e0245541 (2021).
Article CAS PubMed PubMed Central Google Scholar
Osame, M. Review of WHO Kagoshima meeting and diagnostic guidelines for HAM/TSP. Human Retrovirology: HTLV. Raven; p. 191–7 (1990).
Lima, M. A., Bica, R. B. S. & Araújo, A. Q. C. Gender influence on the progression of HTLV-I associated myelopathy/tropical spastic paraparesis. J. Neurol. Neurosurg. Psychiatry. 76, 294–296 (2005).
Article CAS PubMed PubMed Central Google Scholar
Bunce, M. et al. Phototyping: Comprehensive DNA typing for HLA-A, B, C, DRB1, DRB3, DRB4, DRB5 & DQB1 by PCR with 144 primer mixes utilizing sequence-specific primers (PCR-SSP). Tissue Antigens 46, 355–367 (1995).
Article CAS PubMed Google Scholar
We want to thank Dr. Charles Bangham and Dr. Asna Siddiqui for supplying us with reagents and technical assistance to implement the SSP-PCR technique in our laboratory.
This work was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico [Grant No. 477734/07-3] and INI/Fiocruz. The funding agency had no role in the study and collection, analysis, and interpretation of data and in writing the manuscript.
Laboratory of Immunology and Immunogenetics in Infectious Diseases, National Institute of Infectology Evandro Chagas-INI/FIOCRUZ, Avenida Brasil, 4365, Manguinhos, Rio de Janeiro, RJ, 21040-900, Brazil
Doris Schor, Eric Henrique Roma & Maria Glória Bonecini-Almeida
Histocompatibility and Cryopreservation Laboratory, State University of Rio de Janeiro (UERJ), Rio de Janeiro, RJ, 20950-000, Brazil
Luís Cristóvão Porto & Anna Paula Villela
Clinical Research Platform, National Institute of Infectology Evandro Chagas-INI/FIOCRUZ, Rio de Janeiro, RJ, 21040-900, Brazil
Clinical Research Laboratory in Neuroinfections, National Institute of Infectology Evandro Chagas-INI/FIOCRUZ, Rio de Janeiro, RJ, 21040-900, Brazil
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
D.S., L.C.P. and A.P.V. performed the HLA typing experiments; A.Q.C.A. recruited the patients and collected the clinical data; D.S., L.C.P., M.G.B.A. and E.H.R. designed the study; J.C.A. performed the statistical analysis; D.S., L.C.P., E.H.R. and M.G.B.A. analyzed the data and interpreted the results; D.S., L.C.P., E.H.R., J.C.A., A.Q.C.A. and M.G.B.A. contributed to the writing of the manuscript. All authors contributed critically to the manuscript conception. All authors read and approved the final manuscript.
Correspondence to Eric Henrique Roma.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Schor, D., Porto, L.C., Roma, E.H. et al. Putative role of HLA polymorphism among a Brazilian HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) population. Sci Rep 13, 7659 (2023). https://doi.org/10.1038/s41598-023-34757-w
DOI: https://doi.org/10.1038/s41598-023-34757-w
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Scientific Reports (Sci Rep) ISSN 2045-2322 (online)
Hcv Core-Ns3-Ns5 Fusion Antigen-Bio Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.