Effects of host restriction factors and the HTLV-1 subtype on susceptibility to HTLV-1-associated myelopathy/tropical spastic paraparesis
© The Author(s) 2017
Received: 18 July 2016
Accepted: 10 April 2017
Published: 19 April 2017
Although human T-lymphotropic virus type 1 (HTLV-1) infection is a prerequisite for the development of HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP), specific provirus mutations in HAM/TSP have not yet been reported. In this study, we examined whether HAM/TSP patients had the disease-specific genomic variants of HTLV-1 by analyzing entire sequences of HTLV-1 proviruses in these patients, including familial cases. In addition, we investigated the genetic variants of host restriction factors conferring antiretroviral activity to determine which mutations may be related to resistance or susceptibility to HAM/TSP.
The subjects included 30 patients with familial HAM/TSP (f-HAM/TSP), 92 patients with sporadic HAM/TSP (s-HAM/TSP), and 89 asymptomatic HTLV-1 carriers (ACs). In all 211 samples, 37 samples (18%) were classified into transcontinental subtype and 174 samples (82%) were classified as Japanese subtype. Among three groups, the percentage of transcontinental subtype in f-HAM/TSP, s-HAM/TSP and ACs was 33, 23 and 7%, respectively. The frequency of transcontinental subtype was significantly higher in both f-HAM/TSP (p < 0.001) and s-HAM/TSP (p < 0.001) than in ACs. Fifty mutations in HTLV-1 sequences were significantly more frequent in HAM/TSP patients than in ACs, however, they were common only in transcontinental subtype. Among these mutations, ten common mutations causing amino acid changes in the HTLV-1 sequences were specific to the transcontinental subtype. We examined host restriction factors, and detected a rare variant in TRIM5α in HAM/TSP patients. The patients with TRIM5α 136Q showed lower proviral loads (PVLs) than those with 136R (354 vs. 654 copies/104 PBMC, p = 0.003). The patients with the 304L variant of TRIM5α had significantly higher PVLs than those with 304H (1669 vs. 595 copies/104 PBMC, p = 0.025). We could not find any HAM/TSP-specific mutations of host restriction factors.
Transcontinental subtype is susceptible to HAM/TSP, especially in familial cases. Ten common mutations causing amino acid changes in the HTLV-1 gene were specific to the transcontinental subtype. TRIM5α polymorphisms were associated with PVLs, indicating that TRIM5α could be implicated in HTLV-1 replication.
Human T-lymphotropic virus type 1 (HTLV-1) is the etiologic agent of both adult T cell leukemia (ATL) [1, 2] and an inflammatory neurologic disease, HTLV-1-associated myelopathy/tropical spastic paraparesis (HAM/TSP) [3, 4]. HAM/TSP is clinically characterized by chronic progressive spastic paraparesis, urinary incontinence and mild sensory disturbance. The precise pathophysiological mechanism underlying this disease remains unclear. However, the immune response of the HTLV-1-specific CD8-positive cytotoxic T lymphocytes against HTLV-1-infected CD4-positive lymphocytes migrating into the central nervous systems is a pivotal factor for development of HAM/TSP .
HTLV-1 is a complex retrovirus. Its genome contains three typical genes (gag, pol, and env) flanked by long terminal repeat (LTR) sequences at both the 5′ and 3′ ends. In addition, a region known as pX encodes regulatory proteins (Tax, Rex and HBZ) and accessory proteins (p12, p13, and p30) . HTLV-1 is segregated into seven major genetic subtypes (a–g) based on the nucleotide diversity of its LTR region . The cosmopolitan subtype is further divided into five sub-subtypes: (A) Transcontinental, (B) Japanese, (C) West African, (D) North African, and (E) Peruvian Black . While HTLV-1 proviral loads (PVLs) is an important risk factor for HAM/TSP , most studies of HTLV-1 genotype have reported no correlation between nucleotide substitutions and the risk of HAM/TSP [10, 11]. The recent analysis of complete HTLV-1 sequence could not show the specific mutation in clinical states from Brazil . Only Furukawa et al.  reported the association between the variation of HTLV-1 tax gene and the risk of HAM/TSP.
Several host genetic factors, including human leukocyte antigen (HLA) and non-HLA gene polymorphisms affect the occurrence of HAM/TSP . We reported familial clusters of HAM/TSP . These findings indicate that genetic backgrounds may affect individual susceptibility to HAM/TSP. Host restriction factors are important as an initial or early line of defense against infection. In human immunodeficiency virus (HIV) infection, several host restriction factors such as APOBEC3 family, TRIM, tetherin, and SAMHD1 have been well investigated with regard to the biological mechanism of the antiviral response and mutations or polymorphisms of these genes [15, 16]. In HTLV-1 infection, APOBEC3G and SAHMD1 have been studied in terms of antiretroviral activity. APOBEC3G belongs to a family of cytidine deaminases and acts as a potent host restriction factor of retroviral replication through G-to-A hypermutations . In ATL, APOBEC3G generates nonsense mutation in HTLV-1 proviral genomes . SAMHD1 limits HTLV-1 infection of monocyte via STING-Mediated apoptosis . However, no clinical report has been published regarding the genetic variations of host restriction factors in HTLV-1-associated diseases.
DNA sequencing has been remarkably advanced by emerging high-throughput technology. An increased number of disease susceptibility genes in sporadic disease and disease-causing genes in Mendelian diseases have been identified in human genomes . In viral genomes, complete sequencing has been readily obtained and variants have been identified [21, 22].
In this study, we examined whether HAM/TSP patients had the disease-specific genomic variants of HTLV-1 by analyzing entire sequences of proviruses in the patients, including familial cases. In addition, host restriction genes were sequenced to investigate an association with the development of HAM/TSP.
Peripheral blood samples were obtained from 30 patients with familial HAM/TSP (f-HAM/TSP), 92 patients with sporadic HAM/TSP (s-HAM/TSP), and 89 asymptomatic carriers (ACs) after obtaining informed consent. Peripheral blood mononuclear cells (PBMCs) were separated from the blood by Ficoll gradient centrifugation. The study population was enrolled in our database from 1987 to 2014. Patients with HAM/TSP were diagnosed according to the World Health Organization diagnostic criteria . f-HAM/TSP cases were identified as patients with multiple family members suffering from HAM/TSP ; we were able to conduct paired analysis in 18 cases from nine families. This study was approved by the Institutional Review Board of Kagoshima University. All participants provided written informed consent.
Sequencing the complete provirus genomes
DNA was extracted from PBMCs using the DNeasy Blood and Tissue kit (Qiagen, GmbH, Hilden, Germany) according to the manufacturer’s instructions. HTLV-1 proviral loads (PVLs) in the PBMCs were assayed using quantitative PCR with the ABI PRISM 7700TM sequence detection system as reported previously . The complete provirus genome was amplified by nested PCR of the extracted DNA as reported previously . Initially, we amplified two overlapping fragments of the HTLV-1 gene from 100 ng of genomic DNA. PCR products were then subjected to nested PCR. The final PCR products were purified with the MinElute PCR purification kit (Qiagen, GmbH, Hilden, Germany). All PCR products were ligated with the specified indexes, and simultaneously screened on a MiSeq sequencing system in accordance with the manufacturer’s protocol.
The raw sequencing reads were analyzed using CLC genomics Workbench v7.5 (CLC Bio, Aarhus, Denmark). The reads were trimmed using the quality score limit of 0.05 and maximum limit of 2 ambiguous nucleotides. Initially, the reads were mapped to the HTLV-1 complete genome (GenBank accession no. AB513134), allowing up to 2 mismatches . All bases of the alignment were evaluated and single nucleotide polymorphisms were called using the Fixed Ploidy Variant Detection tools of the CLC Genomics Workbench. All comparison alignments were performed and phylogenetic tree was constructed by neighbor-joining method with 1000 bootstrap replications using CLC Genomics Workbench . Five complete HTLV-1 sequences, two transcontinental subtypes (Rk13 and BOI) and three Japanese subtypes (ATK, ATL-YS, and AB513134), were used to construct the phylogenetic tree. Next, to analyze the transcontinental subtype sequence, we constructed the consensus sequence from 37 transcontinental subtypes detected in this study. This consensus sequence was designated “transcontinental subtype reference” and was used as the reference to realign the sequence belonging to transcontinental subtypes and identify possible mutations.
Whole exome sequence
Three micrograms of genomic DNA was processed with the Agilent SureSelect Human All Exon v4 + UTRs capture kit (Agilent, Santa Clara, CA, USA) according to the manufacture’s instructions. The captured DNA was sequenced on an Ilumina Hiseq 2000 (San Diego, CA). Sequences were aligned to the human reference genome (NCBI37/hg19) using the Burrows-Wheeler Aligner (BWA) . Variant calling was performed using SAMtools . Variants were annotated using in-house scripts, which provided the variants list. Previously known variants were annotated from the 1000 Genomes and dbSNP137. We obtained the predicted functional scores of all non-synonymous variants with five prediction algorithms, including SIFT, PolyPhen2, LRT, PhyliP and MutationTaster, using scores obtained from the NSFP database . We obtained mutations or polymorphisms of the candidate gene such as APOBEC3D/F/G/H, TRIM5α, tetherin, and SAMHD1 from whole exome sequence. Rare variants defined minor allele frequency (MAF) ≦1%.
Data were analyzed using SPSS-21 software (SPSS, Chicago, Illinois). Statistical analyses were performed using parametric (t test) and non-parametric tests (Mann–Whitney U test) for continuous variables and χ2 (Pearson χ2 test/Fisher’s exact test) for categorical variables. Differences were considered significant when p < 0.05. Linkage disequilibrium (LD) was analyzed by Haploview software .
Nucleotide sequence accession numbers
The nucleotide sequences of all transcontinental subtypes used for constructing the consensus sequence in this study were submitted to the DDBJ (DNA Data Bank of Japan) database under accession numbers LC192500–192536.
Characteristics of study population
Clinical findings, HTLV-1 subtype, proviral loads, and HTLV-1 sequence changes in patients of HAM/TSP and ACs
f-HAM/TSP (n = 30)
s-HAM/TSP (n = 92)
ACs (n = 89)
Age (years old)
54.7 ± 13.3
58.4 ± 13.2
52.7 ± 14.4
Male (N, %)
Transcontinental/Japanese subtype (n)
20,992 ± 33,571 (N = 29)
30,707 ± 43,985 (N = 80)
7667 ± 17,209 (N = 80)
HTLV-1 proviral loads (copies/104 PBMCs)
819 ± 743
813 ± 936
366 ± 657
No. of SNPs (reference AB513134)a
34.2 ± 18.7
29.6 ± 19.2
22.3 ± 12.0
No. of SNPs (each reference)b
21.4 ± 6.8
20.3 ± 7.5
19.5 ± 5.3
No. of G-to-A mutations (reference AB513134)c
10.2 ± 4.6
8.9 ± 4.4
7.4 ± 3.2
No. of G-to-A mutations (each reference)d
6.6 ± 2.5
6.3 ± 2.4
6.7 ± 2.2
Characterization of HTLV-1 subtype
HTLV-1 sequence changes in HAM/TSP and ACs
Abortive genetic changes of HTLV-1 viral genes in HAM/TSP and ACs
Comparison of HTLV-1 sequences between HAM/TSP patients and ACs
Ten common mutations causing amino acid changes detected in transcontinental subtype in the coding regions of HTLV-1
HTLV-1 proviral gene
3963A > T
5237T > C
5593C > A
6903C > T
6908T > G
7247A > G
7273G > A
7962C > T
8211G > A
7247A > G
Paired analysis of HTLV-1 sequences in familial HAM/TSP patients
HTLV-1 provirus sequences in familial HAM/TSP by paired analysis
Number of variants within family
Host restriction factors in HAM/TSP
Rare variants and polymorphisms of host restriction factors identified in HAM/TSP and ACs
Number of G-to-A mutations of HTLV-1 sequence
8.8 ± 4.1
6.5 ± 2.3
6.8 ± 2.2
6.4 ± 2.4
5.2 ± 1.7
6.6 ± 2.4
6.4 ± 2.3
6.6 ± 2.4
Characterization of TRM5α SNPs and proviral loads in HAM/TSP and ACs
AA change in TRIM5α
Minor allele frequency
HTLV-1 proviral loads (copies/104 PBMC)
623 ± 894
606 ± 802
514 ± 344
612 ± 829
637 ± 790
608 ± 832
354 ± 401
654 ± 870
1669 ± 1411
595 ± 809
612 ± 828
In this study, we demonstrated that the transcontinental subtype was significantly more frequent in HAM/TSP patients, especially in f-HAM/TSP patients, than in ACs. We constructed a phylogenic tree based on the complete sequence data and divided it into subtypes. There are two distinct subtypes in Japan; the most frequently observed (nearly 80%) HTLV-1 subtypes belong to the Japanese subtype, while the minor transcontinental subtype (20%) seems to cluster in the southern islands of Kyushu and Ryukyu . Furukawa et al.  reported the association between transcontinental subtype and the risk of HAM/TSP. However, almost previous studies demonstrated that the nucleotide substitutions in some fragments of HTLV-1 genome are specific for the geographic origin of the patients rather than for the type of associated pathologies [31, 32]. In endemic area of HTLV-1, including the area examined in previous reports, one of the HTLV-1 subtypes was predominantly widespread. Kagoshima is a unique area where both HTLV-1 subgroups exist. Therefore, it is an appropriate region in which to evaluate whether the HTLV-1 subtype could influence disease development. Additionally, most previous reports except for a report by Furukawa et al. examined relatively few patients. We examined 211 samples and showed that the transcontinental subtype was more frequent in HAM/TSP patients, especially in familial cases, than in ACs. Comparison of the sequences of transcontinental subtype with those of Japanese subtype showed two transcontinental subtype-specific substitutions that changed amino acid in the tax gene . A recent study showed similar results with 42 HAM/TSP patients and that the HBZ mRNA expression was significantly higher in HAM/TSP patients with Japanese subtype than in those with transcontinental subtype , while we analyzed the complete sequence of the transcontinental subtype in Japan and detected the transcontinental subtype-specific substitutions. Analyzing viral genes other than tax and HBZ between HAM/TSP patients and ACs might be useful for clarifying the differences in development with HAM/TSP between two HTLV-1 subtypes.
HTLV-1 has remarkably low genetic variability. In this study, the number of single nucleotide polymorphisms (SNPs) was relatively high in HAM/TSP patients compared to ACs when we analyzed all HTLV-1 sequences using the same reference of the Japanese subtype, which was the most prevalent in Japan. However, when we segregated the HTLV-1 subtypes and analyzed them using the references adequate for each subtype, the number of SNPs was not different, and the overall divergence was less than one percent. We could not find any specific mutations in HAM/TSP patients. It is unclear whether HAM/TSP-specific mutations could be detected by using a larger number of samples because the mutations in the HTLV-1 sequence are relatively rare. We examined a complete consensus sequence in this study, but it may be beneficial to analyze the quasispecies and integration sites of the HTLV-1 provirus. Recently, quasispecies-derived mechanisms have been shown to mediate the adaptability for persistence to escape from the host immune responses and/or to acquire resistance to antiviral agents, especially in HIV, hepatitis C virus (HCV) and hepatitis B virus (HBV) . The frequency of mutations was small in the region with sufficient coverage in our study. The integration site has been studied in ATL. Hot spots of integration could not be identified, but the association of nearby host genes has been related to hematologic malignancies . Identifying the integration site would also be useful in studying HAM/TSP. Furthermore, G-to-A mutations were common and caused the most nonsense mutations. These results suggested that APOBEC3G might exert an effect on viral mutations . We examined the host APOCEC3G gene by whole exome sequencing and analyzed the association between polymorphisms and/or mutations in the host and the number of SNPs and/or G-to-A mutations in the HTLV-1 sequence, but no significant difference was detected. The frequencies of nonsense, insertion and deletion mutations of HTLV-1 viral genes in the HAM/TSP and AC groups were 6.6 and 5.6%, and which were fewer than previously reported for ATL cases (46.7%) . The lower frequencies of insertion and deletion mutations could be due to MiSeq using short reads.
The examination of the HTLV-1 sequences intra-family, even among HTLV-1 carriers, has scarcely reported. Examination of the sequences of gp46 showed no difference between carriers and their partners experienced of seroconversion . In this study, we found few differences within a family, despite examining complete sequences. Although there were only two parent–child pairs, pediatric cases had an earlier age of onset. We recently reported that the mean age of onset of 40 familial HAM/TSP cases was 41 years. Compared to this value, not only children but also parents had an earlier age of onset. Although we could not accurately identify the source of the infection, a carrier child whose parent developed HAM/TSP needs to be carefully followed up. A previous report for PVL between familial carriers showed variations. In this study, PVLs, anti-HTLV-1 antibody titer, and motor disability score were relatively similar in HAM/TSP cases.
We examined the host restriction factors known in HIV infection, such as APOBEC3, TRIM5α and SAMHD1, in HAM/TSP patients and carriers. APOBEC3 proteins are members of the cytidine deaminase family, which facilitates deamination during reverse transcription, resulting in nucleotide mutations. APOBEC3G deaminates cytosine residues in single-stranded DNA during reverse transcription, causing G-to-A hypermutations and inhibiting the viral replication. These editing enzymes likely evolved to inhibit retrovirus replication, including that of HIV [15, 16]. In HTLV-1, a nucleocapsid domain inhibits the incorporation of APOBEC3G into the virions and escape from the antiviral effect . While this activity is not foolproof, APOBEC3G induce G-to-A mutations and consequently causes nonsense mutations . An analysis of these restriction mechanisms has been conducted in vitro and in vivo; however, a population-based study has never been attempted to establish a relationship among host gene mutations, clinical states, and viral properties. This study is the first report of the host restriction factors in HTLV-1-associated diseases.
First, we examined APOBEC3G, and the patients having rare variants tended to show more G-to-A mutations than those without rare variants. However, no significance differences were observed in the clinical states, PVLs or viral mutations. Other APOBEC3 family members showed no significant difference. Next, we examined TRM5α, and patients having two polymorphisms in TRIM5α showed a significant difference in PVLs. TRIM5α is characterized by a RING finger zip binding domain, one or two B-box domains and a coiled-coil domain; it exhibits RING domain-dependent E3 ubiquitin ligase activity [38, 39]. TRIM5α recognizes the entry of retroviral capsids into the cytoplasm and blocks viral infection . Several reports have described the relationship between TRIM5α genetic-variants and HIV disease progression. In this study, the 136Q polymorphism showed lower PVLs than 136R. The 136Q allele was more frequent in HIV-1 resistant individuals in US-based natural history and Pumwani cohort studies [41, 42]. Codon 136 is located within the coiled-coil domain. A switch from arginine to glutamine could be implicated in TRIM5α oligomerization and increased antiviral activity. The rare variant of 304L showed higher PVLs than 304H and was exclusively harbored in HAM/TSP patients. H304L is a noncoding variant located in the 3′UTR. The influence of this variant remains unknown and may be an accidental coincidence. Collectively, TRIM5α may influence host resistance to HTLV-1 infection.
It is important to examine the host-virus interaction to understand infectious disorders. One of the approaches to elucidate such interaction is to identify the polymorphisms and/or mutations in human and viral genomes. Emerging new technologies for nucleotide sequencing allow us to obtain vast numbers of sequencing data. It is difficult to detect the mutations associated with the disease and the pathophysiology in these massive data. Many new mutations were identified to combine the various analyzing methods in ATL . Case–control studies such as genome-wide association studies, which require large numbers of cases and controls, are useful strategies by which to detect susceptible genes in common diseases. However, obtaining a large sample size of rare diseases, including HAM/TSP, and analyzing the causative genes is difficult. Optimal analyses such as advanced bioinformatics analyses need to be performed. Analyzing families with many affected individuals may be more powerful and fruitful in finding rare causal variants [44, 45]. Multiple rare variants in COQ2 have been reported to be associated with multiple system atrophy by analyzing familial clustering . It is important for us to conduct the studies in different ethnic populations to obtain a larger sample size and to examine the cases focusing on familial clustering that suggest genetic influence.
In conclusion, the transcontinental subtype is more frequent in patients with HAM/TSP, especially in familial cases. Ten common mutations with amino acid changes in the transcontinental subtype compared with Japanese subtypes were detected and implicated in the development of HAM/TSP. A R136Q polymorphism in TRIM5α results in lower PVLs. This result indicates that TIRM5α may be implicated in HTLV-1 infection.
SN designed the experiments, analyzed the data and wrote the manuscript. YT and DK performed PCR study. EM, TM, RK provided the information of the patients and supervised the project. SI and HT developed the concept and supervised the project. All authors read and approved the final manuscript.
We thank Yuji Okamoto MD, Junhui Yuan MD and Ms. Yuka Komai for their excellent technical assistances.
The authors declare that they have no competing interests.
Availability of data and materials
All data analysed during this study are included in this published article and supplementary information.
Ethics approval and consent to participate
This study was approved by the Institute Review Boards of Kagoshima University Graduate School of Medical and Dental Sciences. All participants provided their written informed consent. This report doesn’t contain any individual data.
This study was supported by the Research program for conquering intractable disease from Japan Agency for Medical Research and development, AMED and JSPS KAKENHI Grant Numbers 25293205 and 24133701.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hinuma Y, Nagata K, Hanaoka M, Nakai M, Matsumoto T, Kinoshita KI, et al. Adult T-cell leukemia: antigen in an ATL cell line and detection of antibodies to the antigen in human sera. Proc Natl Acad Sci USA. 1981;78:6476–80.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 USA. 1980;77:7415–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Osame M, Usuku K, Izumo S, Ijichi N, Amitani H, Igata A, et al. HTLV-I associated myelopathy, a new clinical entity. Lancet. 1986;1:1031–2.View ArticlePubMedGoogle Scholar
- 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:407–10.View ArticlePubMedGoogle Scholar
- Matsuura E, Kubota R, Tanaka Y, Takashima H, Izumo S. Visualization of HTLV-1-specific cytotoxic T lymphocytes in the spinal cords of patients with HTLV-1-associated myelopathy/tropical spastic paraparesis. J Neuropathol Exp Neurol. 2015;74:2–14.View ArticlePubMedGoogle Scholar
- Matsuura E, Yamano Y, Jacobson S. Neuroimmunity of HTLV-I infection. J Neuroimmune Pharmacol. 2010;5:310–25.View ArticlePubMedPubMed CentralGoogle Scholar
- Verdonck K, Gonzalez E, Van Dooren S, Vandamme AM, Vanham G, Gotuzzo E. Human T-lymphotropic virus 1: recent knowledge about an ancient infection. Lancet Infect Dis. 2007;7:266–81.View ArticlePubMedGoogle Scholar
- Proietti FA, Carneiro-Proietti AB, Catalan-Soares BC, Murphy EL. Global epidemiology of HTLV-I infection and associated diseases. Oncogene. 2005;24:6058–68.View ArticlePubMedGoogle Scholar
- Nagai M, Usuku K, Matsumoto W, Kodama D, Takenouchi N, Moritoyo T, et al. Analysis of HTLV-I proviral load in 202 HAM/TSP patients and 243 asymptomatic HTLV-I carriers: high proviral load strongly predisposes to HAM/TSP. J Neurovirol. 1998;4:586–93.View ArticlePubMedGoogle Scholar
- Mahieux R, de The G, Gessain A. The tax mutation at nucleotide 7959 of human T-cell leukemia virus type 1 (HTLV-1) is not associated with tropical spastic paraparesis/HTLV-1-associated myelopathy but is linked to the cosmopolitan molecular genotype. J Virol. 1995;69:5925–7.PubMedPubMed CentralGoogle Scholar
- Pessoa R, Watanabe JT, Nukui Y, Pereira J, Kasseb J, de Oliveira AC, et al. Molecular characterization of human T-cell lymphotropic virus type 1 full and partial genomes by Illumina massively parallel sequencing technology. PLoS ONE. 2014;9:e93374.View ArticlePubMedPubMed CentralGoogle Scholar
- Furukawa Y, Yamashita M, Usuku K, Izumo S, Nakagawa M, Osame M. Phylogenetic subgroups of human T cell lymphotropic virus (HTLV) type I in the tax gene and their association with different risks for HTLV-I-associated myelopathy/tropical spastic paraparesis. J Infect Dis. 2000;182:1343–9.View ArticlePubMedGoogle Scholar
- Saito M. Immunogenetics and the pathological mechanisms of human T-cell leukemia virus type 1-(HTLV-1-) associated myelopathy/tropical spastic paraparesis (HAM/TSP). Interdiscip Perspect Infect Dis. 2010;2010:478461.View ArticlePubMedPubMed CentralGoogle Scholar
- Nozuma S, Matsuura E, Matsuzaki T, Watanabe O, Kubota R, Izumo S, et al. Familial clusters of HTLV-1-associated myelopathy/tropical spastic paraparesis. PLoS ONE. 2014;9:e86144.View ArticlePubMedPubMed CentralGoogle Scholar
- Santa-Marta M, de Brito PM, Godinho-Santos A, Goncalves J. Host factors and HIV-1 replication: clinical evidence and potential therapeutic approaches. Front Immunol. 2013;4:343.View ArticlePubMedPubMed CentralGoogle Scholar
- Malim MH, Bieniasz PD. HIV restriction factors and mechanisms of evasion. Cold Spring Harb Perspect Med. 2012;2:a006940.View ArticlePubMedPubMed CentralGoogle Scholar
- Sheehy AM, Gaddis NC, Choi JD, Malim MH. Isolation of a human gene that inhibits HIV-1 infection and is suppressed by the viral Vif protein. Nature. 2002;418:646–50.View ArticlePubMedGoogle Scholar
- Fan J, Ma G, Nosaka K, Tanabe J, Satou Y, Koito A, et al. APOBEC3G generates nonsense mutations in human T-cell leukemia virus type 1 proviral genomes in vivo. J Virol. 2010;84:7278–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Sze A, Belgnaoui SM, Olagnier D, Lin R, Hiscott J, van Grevenynghe J. Host restriction factor SAMHD1 limits human T cell leukemia virus type 1 infection of monocytes via STING-mediated apoptosis. Cell Host Microbe. 2013;14:422–34.View ArticlePubMedGoogle Scholar
- Tsuji S. The neurogenomics view of neurological diseases. JAMA Neurol. 2013;70:689–94.View ArticlePubMedGoogle Scholar
- Dudley DM, Chin EN, Bimber BN, Sanabani SS, Tarosso LF, Costa PR, et al. Low-cost ultra-wide genotyping using Roche/454 pyrosequencing for surveillance of HIV drug resistance. PLoS ONE. 2012;7:e36494.View ArticlePubMedPubMed CentralGoogle Scholar
- Willerth SM, Pedro HA, Pachter L, Humeau LM, Arkin AP, Schaffer DV. Development of a low bias method for characterizing viral populations using next generation sequencing technology. PLoS ONE. 2010;5:e13564.View ArticlePubMedPubMed CentralGoogle Scholar
- Osame M, Janssen R, Kubota H, Nishitani H, Igata A, Nagataki S, et al. Nationwide survey of HTLV-I-associated myelopathy in Japan: association with blood transfusion. Ann Neurol. 1990;28:50–6.View ArticlePubMedGoogle Scholar
- Jalan N, Kumar D, Andrade MO, Yu F, Jones JB, Graham JH, et al. Comparative genomic and transcriptome analyses of pathotypes of Xanthomonas citri subsp. citri provide insights into mechanisms of bacterial virulence and host range. BMC Genom. 2013;14:551.View ArticleGoogle Scholar
- Lu QB, Zhang XA, Wo Y, Xu HM, Li XJ, Wang XJ, et al. Circulation of Coxsackievirus A10 and A6 in hand-foot-mouth disease in China, 2009–2011. PLoS ONE. 2012;7:e52073.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Durbin R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009;25:1754–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu X, Jian X, Boerwinkle E. dbNSFP: a lightweight database of human nonsynonymous SNPs and their functional predictions. Hum Mutat. 2011;32:894–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–5.View ArticlePubMedGoogle Scholar
- Sabouri AH, Saito M, Usuku K, Bajestan SN, Mahmoudi M, Forughipour M, 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. 2005;86:773–81.View ArticlePubMedGoogle Scholar
- Komurian F, Pelloquin F, de The G. In vivo genomic variability of human T-cell leukemia virus type I depends more upon geography than upon pathologies. J Virol. 1991;65:3770–8.PubMedPubMed CentralGoogle Scholar
- Segurado AA, Biasutti C, Zeigler R, Rodrigues C, Damas CD, Jorge ML, et al. Identification of human T-lymphotropic virus type I (HTLV-I) subtypes using restricted fragment length polymorphism in a cohort of asymptomatic carriers and patients with HTLV-I-associated myelopathy/tropical spastic paraparesis from Sao Paulo, Brazil. Mem Inst Oswaldo Cruz. 2002;97:329–33.View ArticlePubMedGoogle Scholar
- Yasuma K, Matsuzaki T, Yamano Y, Takashima H, Matsuoka M, Saito M. HTLV-1 subgroups associated with the risk of HAM/TSP are related to viral and host gene expression in peripheral blood mononuclear cells, independent of the transactivation functions of the viral factors. J Neurovirol. 2016;22:416–30.View ArticlePubMedGoogle Scholar
- Domingo E, Sheldon J, Perales C. Viral quasispecies evolution. Microbiol Mol Biol Rev. 2012;76:159–216.View ArticlePubMedPubMed CentralGoogle Scholar
- Cook LB, Melamed A, Niederer H, Valganon M, Laydon D, Foroni L, et al. The role of HTLV-1 clonality, proviral structure, and genomic integration site in adult T-cell leukemia/lymphoma. Blood. 2014;123:3925–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Iga M, Okayama A, Stuver S, Matsuoka M, Mueller N, Aoki M, et al. Genetic evidence of transmission of human T cell lymphotropic virus type 1 between spouses. J Infect Dis. 2002;185:691–5.View ArticlePubMedGoogle Scholar
- Derse D, Hill SA, Princler G, Lloyd P, Heidecker G. Resistance of human T cell leukemia virus type 1 to APOBEC3G restriction is mediated by elements in nucleocapsid. Proc Natl Acad Sci USA. 2007;104:2915–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Reymond A, Meroni G, Fantozzi A, Merla G, Cairo S, Luzi L, et al. The tripartite motif family identifies cell compartments. EMBO J. 2001;20:2140–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Meroni G, Diez-Roux G. TRIM/RBCC, a novel class of ‘single protein RING finger’ E3 ubiquitin ligases. Bioessays. 2005;27:1147–57.View ArticlePubMedGoogle Scholar
- Stremlau M, Owens CM, Perron MJ, Kiessling M, Autissier P, Sodroski J. The cytoplasmic body component TRIM5alpha restricts HIV-1 infection in Old World monkeys. Nature. 2004;427:848–53.View ArticlePubMedGoogle Scholar
- Javanbakht H, An P, Gold B, Petersen DC, O’Huigin C, Nelson GW, et al. Effects of human TRIM5alpha polymorphisms on antiretroviral function and susceptibility to human immunodeficiency virus infection. Virology. 2006;354:15–27.View ArticlePubMedGoogle Scholar
- Price H, Lacap P, Tuff J, Wachihi C, Kimani J, Ball TB, et al. A TRIM5alpha exon 2 polymorphism is associated with protection from HIV-1 infection in the Pumwani sex worker cohort. AIDS. 2010;24:1813–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Kataoka K, Nagata Y, Kitanaka A, Shiraishi Y, Shimamura T, Yasunaga J, et al. Integrated molecular analysis of adult T cell leukemia/lymphoma. Nat Genet. 2015;47:1304–15.View ArticlePubMedGoogle Scholar
- Ramagopalan SV, Dyment DA, Cader MZ, Morrison KM, Disanto G, Morahan JM, et al. Rare variants in the CYP27B1 gene are associated with multiple sclerosis. Ann Neurol. 2011;70:881–6.View ArticlePubMedGoogle Scholar
- Dyment DA, Cader MZ, Chao MJ, Lincoln MR, Morrison KM, Disanto G, et al. Exome sequencing identifies a novel multiple sclerosis susceptibility variant in the TYK2 gene. Neurology. 2012;79:406–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Multiple-System Atrophy Research Collaboration. Mutations in COQ2 in familial and sporadic multiple-system atrophy. N Engl J Med. 2013;369:233–44.Google Scholar