Next Article in Journal
Identification, Virulence, and Molecular Characterization of a Recombinant Isolate of Grass Carp Reovirus Genotype I
Previous Article in Journal
To Include or Occlude: Rational Engineering of HCV Vaccines for Humoral Immunity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years

1
Department of Microbiology, Asan Medical Center, University of Ulsan College of Medicine, Seoul 05505, Korea
2
HIV Databases, Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87544, USA
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(5), 806; https://doi.org/10.3390/v13050806
Submission received: 5 April 2021 / Revised: 25 April 2021 / Accepted: 29 April 2021 / Published: 30 April 2021
(This article belongs to the Section Animal Viruses)

Abstract

:
We aimed to investigate whether the sequence length of HIV-1 increases over time. We performed a longitudinal analysis of full-length coding region sequences (FLs) during an HIV-1 outbreak among patients with hemophilia and local controls infected with the Korean subclade B of HIV-1 (KSB). Genes were amplified by overlapping RT-PCR or nested PCR and subjected to direct sequencing. Overall, 141 FLs were sequentially determined over 30 years in 62 KSB-infected patients. Phylogenetic analysis indicated that within KSB, two FLs from plasma donors O and P comprised two clusters, together with 8 and 12 patients with hemophilia, respectively. Signature pattern analysis of the KSB of HIV-1 revealed 91 signature nucleotide residues (1.1%). In total, 48 and 43 signature nucleotides originated from clusters O and P, respectively. Six positions contained 100% specific nucleotide(s) in clusters O and P. In-depth FL analysis for over 30 years indicated that the KSB FL significantly increased over time before combination antiretroviral therapy (cART) and decreased with cART. This increase occurred due to the significant increase in env and nef genes, originating in the variable regions of both genes. The increase in sequence length of HIV-1 over time suggests an evolutionary direction.

1. Introduction

We previously conducted a nationwide genetic analysis of HIV-1 using sera from individuals in the early stages of HIV-1 infection (before 1994) to identify the cause of the HIV-1 outbreak among patients with hemophilia in Korea in 1990–1994. These molecular epidemiological studies revealed that the viruses from two cash-paid plasma donors were incompletely inactivated in the process of manufacturing clotting factor IX and were identified as the agents of infection among 20 HIV-1-infected patients with hemophilia [1,2,3,4,5,6]. The viruses in 8 and 12 patients with hemophilia infected with the Korean subclade of HIV-1 subtype B (KSB) originated from plasma donors O and P, respectively. In these studies [1,2,3,4,5,6,7], we conducted an in-depth genetic analysis, but a few genes are yet to be explored. The sequence length significantly affects the extent of clustering in a phylogenetic tree [8,9]. KSB is a distinct, monophyletic clade within the HIV-1 subtype B and is presumed to have originated from strains in the USA through a founder effect [8,10,11,12,13]. The most recent common ancestor virus is estimated to have been active around 1984 [14]; however, the earliest infection case was diagnosed in 1988 [2,3,4,5,6].
In this study, we identified signature pattern residues in the full-length coding region sequence (FL) of KSB. We also performed a phylogenetic analysis at the FL level in 64 patients, including 20 patients with hemophilia B (HPs) [1,2,3,4]. We confirmed a previously postulated epidemiological link between the viruses that infected 20 HPs and 2 plasma donors and the viruses that infected local control patients [2,3,4,5,6]. Additionally, their longitudinal sequence analyses conducted over approximately 30 years indicate that the sequence length in FL significantly increases over time before the administration of combination antiretroviral therapy (cART). It has been reported that the progression in subtype D-infection is four-fold faster than that in subtype-B infections [15]. The sequence length of subtype D is the shortest among group M.
To date, our study is one of the most comprehensive and longest longitudinal studies conducted to examine the evolution of the HIV-1 subtype B originated from a single source of HIV-1 [1,2,3,4,5,6], and its results provide novel insights into the pathogenesis of HIV-1 infection over time.

2. Materials and Methods

2.1. Ethical Statement

This study was approved by the institutional review board of the Asan Medical Center (Code 2012-0390, 4 June 2012). All subjects provided their informed consent for inclusion before participating in the study. The study was conducted in accordance with the Declaration of Helsinki.

2.2. Patients and Samples

Four HIV-1 infected plasma donors were diagnosed during primary infection in 1990–1992. Their plasma was used to manufacture the domestic clotting factor IX. Viruses from donors O and P were incompletely inactivated. The details of the procedure have been described previously [1,2,3,4,5,6,16,17]. Briefly, 20 patients with hemophilia (HPs 1–20) were diagnosed with HIV-1 infection between 1990 and 1994. HP 21 was infected with HIV-1 via imported factor 9 before 1987 and was diagnosed in 1987. In this study, FLs were sequenced from 62 KSB-infected patients, including 3 plasma donors (O, P, and R) and 20 patients with hemophilia (designated as HP 1–20), 1 infected with CRF02-AG, 1 infected with subtype D, and 1 infected with subtype B (Table S1). In this study, the FL sequence denotes the length from the start codon of Gag to Nef with the terminal stop codon.

2.3. RNA/DNA Preparation and FL Gene Amplification

Blood samples were collected from 20 HPs at 6-month intervals for CD4+ T cell measurements. Sera were used before the year 2000, and peripheral blood mononuclear cells (PBMC) were used after 2000 for PCR amplification.
Total RNA was extracted from 300 μL serum samples using a QIAamp UltraSens Viral RNA kit (Qiagen, Hilden, Germany), and 2 μL aliquots of RNA were reverse transcribed by mixing with 1 μL of oligo(dT), 1 μL of dNTPs, and 6 μL of DEPC treated water, followed by incubating the mixture at 65 ℃ for 5 min and then on ice for 1 min [2,3,4,5,6]. Next, 4 μL of 5× buffer, 2 μL of 0.1M DTT, 3 μL of DEPC-treated water (Ambion Inc., Foster City, CA, USA), and 1 μL of Superscript III reverse transcriptase (Invitrogen, Carlsbad, CA, USA) were added to each sample. The samples were incubated for 50 min at 50 ℃. The reaction was terminated by incubation for 5 min at 85 ℃ and then on ice for 1 min.
DNA was extracted from 400 μL of PBMC samples using a QIAamp DNA Mini kit (Qiagen, Hilden, Germany), and 5 μL aliquots of DNA were used for the nested PCR.
The vif, vpr, tat/rev, and vpu regions (1.2-kb) were amplified via nested PCR using TaKaRa r-taq (Takara Bio Inc., Shiga, Japan). The first and second PCRs were performed in 20 μL and 50 μL reaction mixtures, respectively. The outer primer pairs used were 545 and KMK2, whereas the inner primer pairs used were 548F and LA106 (Table 1). After an initial denaturation at 95 ℃ for 10 min, 35 cycles were run at 95 ℃ for 30 s, at 52 ℃ for 30 s, and at 72 ℃ for 2 mins and 30 s, followed by a final extension step at 72 ℃ for 10 mins. The second PCR was performed with 1 μL of the first PCR product; the cycling conditions were set as follows: 95 ℃ for 30 s, 57 ℃ for 30 s, and 72 ℃ for 1 min and 30 s, as well as a final extension at 72 ℃ for 10 mins. The procedure used for the amplification of the remaining genes is described elsewhere [2,4,5,6]. A maximum of five PCR, including that using a negative control, were performed per sample at a given time. To avoid selection bias, all positive PCR amplicons were sequenced and used for FL. Amplification was performed via five overlapping PCR steps (Table 1) [16,17], followed by direct sequencing using an Applied Biosystems 3730XL DNA Analyzer (Foster City, CA, USA) [5].

2.4. Phylogenetic Tree Analysis

In total, 70 FL sequences were obtained from 20 patients with hemophilia. The sequences from 42 local controls and 1 subtype B infected patient were aligned against the HIV-1 subtype reference dataset from the HIV Sequence Database (http://hiv-web.lanl.gov/content/hiv-db/Subtype_REF/align.html) (accessed on 20 July 2020) Phylogenetic trees were constructed using IQ TREE with 1000 bootstrap replicates [18].

2.5. Viral Signature Pattern Analysis (VESPA)

The VESPA program (http://www.hiv.lanl.gov/content/sequence/VESPA/vespa.html) (accessed on 20 July 2020). was used to identify the sites within each sequence group that was distinct from those in other groups [19].

2.6. Statistical Analysis

Data are presented as means ± standard deviation. Statistical significance was determined using Student’s two-tailed t-tests, paired t-test, chi-square tests, Fisher’s exact test, and Pearson’s correlation coefficient using MedCalc. Results were deemed statistically significant when the p-value was <0.05.

2.7. Nucleotide Sequence Data

The GenBank accession numbers for the sequences in this study are listed as follows: AF224507, AY839827, DQ054367, DQ295192-96, DQ837381, JQ316126-37, JQ341411, JQ429433, KF561435-43, KJ140245-66, KU869610, MK577478-81, MK871374, MG461319-22, MN237642-46, MN043576-607, MT101871, MT224125-27, MT559044-66, MT582420-24, MT679550-53, MW405263-343, and MW881608-49. Subtype B sequences were randomly retrieved from the Los Alamos National Laboratory (LANL) HIV Database.

3. Results

3.1. Origin of the KSB of Subtype B

A major strength of this study is that the sequences of the earliest KSB-infected patients were included. We found that patient BGO, who was diagnosed in July 1988, was the first patient infected with KSB (MT559045). Hence, we obtained the FLs in 36, 23, and 3 patients diagnosed in 1988–1991, 1992–1993, and after 1993, respectively (Table S1 and Figure 1). These 36 patients who were diagnosed in 1988–1991 correspond to 84% of all 43 KSB-infected patients diagnosed in this period [4].

3.2. Molecular Epidemiologic Data on the FL HIV-1 Gene

In 20 HPs, 71 FLs were obtained at 60 time points over 287 ± 99 months (about 24 years) from the outbreak, in January 1990. We obtained 169 FLs from 65 patients. Of these, the sequences were obtained from 2 or more samples collected on different dates from 30 patients, including 19 patients with hemophilia (Table S1). Phylogenetic analysis revealed that the earliest 62 FLs obtained from 62 patients (20 HPs and 42 local controls) belonged to KSB, whereas 2 FLs obtained from 2 patients belonged to subtypes B and D (Figure 1). The 62 KSB FLs were subdivided into several clusters, including two large clusters (“O,” which comprised 9 sequences and “P,” which comprised 13 sequences) that included 20 HPs and plasma donors O and P. The bootstrap values of the nodes for clusters O, P, and KSB were all 100%, as determined using 1000 bootstrap replicates (Figure 1).

3.3. Korean Signature Pattern Amino Acid Residues

We previously reported the signature pattern of amino acids at residues 12 and 26 in the Gag and Env proteins [5,6]. Additionally, 8 and 9 signature pattern amino acids were determined in the Vif and Nef proteins, respectively [20,21]. We found 31 novel Korean signature nucleotides in the pol gene based on our previous report [2] compared with 31 sequences from 15 subtype B infected Korean patients (Table S2). Of those, 11 were nonsynonymous substitutions, and 20 were synonymous substitutions, compared to those in subtype B.
Overall, in the FL sequences that were over 8609–8618 bp long, the signature pattern analysis indicated 91 signature nucleotides (16, 21, 6, 2, 2, 2, 35, and 7 in gag, pol, vif, vpr, tat/rev, vpu, env, and nef genes, respectively; 1.1%) that distinguished 20 HPs and 42 local controls within KSB (p < 0.05). In total, 48 and 43 signature nucleotides originated from clusters O and P, respectively. Of those, only six positions in gag, pol, vif, and env genes contained 100% specific nucleotide(s) positions in clusters O and P [2,3,4,5,6], compared to 0% in local controls (Table S3).

3.4. Sequence Identities of HPs Compared to Plasma Donors O and P

The earliest FLs from the donors O and P were 8609 bp and 8618 bp long, respectively. The sequence identity between the earliest sequences determined in October 1991 and the last sequences determined in January 2002 (8627 bp) from donor O was 96.5%. The sequence similarity in cluster O between the earliest sequences from donor O and the sequence of each HP averaged 97.7% ± 0.8%. In cluster P, the sequence similarity between the earliest sequences from donor P and each HP was 97.4% ± 1.4%.
We determined the correlation between the sampling intervals after the outbreak (January 1990) and the number of nucleotide differences observed, relative to those of the corresponding plasma donor. Four patients who were first sampled in 2002 exhibited the lowest sequence identity (Figure 2A). The sequence identity significantly dropped over time since the outbreak. The lowest sequence identity was 89.8% at 153 months in October 2002 in HP-20. The overall correlation coefficient, γ, was estimated to be −0.85 in 19 HPs (p < 0.001) (Figure 2A).
When we analyzed 40 FLs before the administration of cART, the γ value was 0.77 (p < 0.001) (Figure 2B), whereas at the time of administration of cART, the γ value was−0.17, compared with that observed just before the administration of cART (Figure 2C). In the same context, the correlation between the sampling year and sequence length before the administration of cART was also significant in 20 HPs (r = 0.79; p < 0.001) (Figure 2D).
In addition, we found this phenomenon in 42 local control patients without hemophilia (r = 0.38; p < 0.01) (Figure 2E), and the correlation was higher in 12 local controls with sequences determined at 2 time points before cART (r = 0.61; p = 0.001).

3.5. Sequence Length of HIV-1 KSB Significantly Increases over Time before cART Administration

Among the earliest 20 FLs in 20 HPs, we excluded 4 HPs because the sequences were obtained after 10 years since the outbreak in January 1990. The earliest sequences obtained over 32 ± 5 months from the outbreak revealed 8614 ± 17 nucleotides (n = 16) with an increase of 4.1 ± 16 nucleotides per year compared to the corresponding sequences of the donors. The second sequences obtained more than 5 years after the outbreak were determined from samples collected from 17 HPs before cART (n = 23). This revealed 8650 ± 22 nucleotides with a significant increase of 39 ± 22 nucleotides over 143 ± 41 months from the outbreak (p < 0.0001); however, the first sequences of donors O and P obtained in 1991 and 1993 were used.
We found that the sequence length increased significantly over time before cART (Figure 2B). The sequence length by gene was analyzed to confirm which gene demonstrated an increase in sequence length. There was a significant increase in the sequence length of env (p < 0.001) and nef (p < 0.05) genes (Table 2). This increase originated in the variable regions of both genes (data not shown).
In contrast, we analyzed the changes in sequence length within each HP. The sequence length increased by 39 ± 30 nucleotides over 110 ± 43 months between the first and last sequences in 16 HPs, with sequences determined at ≥2 time points before cART (p = 0.0001) (Table 2). This corresponds to an annual increase of 4.4 nucleotides and a requirement of 47 and 32 years to reach the mean length of SIVsm and HIV-2, respectively.
The sequence length in the env gene increased by 31.7 ± 25 nucleotides during the same period (p = 0.0001).
In contrast, the sequence length decreased by 8.0 ± 23 nucleotides over 181 ± 47 months at the administration of cART in 16 HPs, compared to the sequence length determined just before the administration of cART (Figure S1).
To confirm whether this phenomenon occurs in patients without hemophilia, we analyzed the sequences in 12 KSB-infected patients at ≥2 time points among 42 local controls. There was also a significant increase in sequence length between the FLs over 104 ± 43 months (33 ± 29 nucleotides) (p < 0.01). The increase in the env sequence length (24 ± 25 nucleotides) over 104 ± 43 months (p < 0.05) significantly contributed to the 33 ± 29 nucleotide increase, whereas the increase in sequence length of the nef gene was not significant.
Combining 20 HPs with 42 local controls, we observed significant correlations between 100 FLs and the sampling year (p < 0.001) (Figure 2F), and between the sequence length of env and the sampling year (r = 0.45, p < 0.001); no significant correlation was observed with cART.
Additionally, we found a significant correlation even in individual patients with long-term slow progression (Figure 3).

3.6. The Sequence Length in HIV-1 Subtype B Significantly Increases over Time

To confirm whether the sequence length increases in subtypes other than KSB, we randomly selected 64 FLs (subtype B) from the LANL HIV Database from 1983 to 1998 before cART. The correlation between the sampling year and the sequence length was significant for FLs (r = 0.43, p < 0.001), env (r = 0.26, p < 0.05), and nef (r = 0.30, p < 0.05) genes (Figure 4A–C).

3.7. Clinical Significance of the Increase in Sequence Length in HIV-1 Subtype B

Several studies have reported that the V1-V2 envelope loops and the env sequence length increase by approximately 1% per year in the early phases of typical infections [24,25,26]. There had been several reports on the elongation of the V2 region in long-term nonprogressors (LTNPs) [27,28,29,30]. We analyzed the correlation between the number of amino acids in the V2 region and the duration since the diagnosis of KSB and the subtype B infection. Consistent with the increase in sequence length over time, there were significant correlations between the number of amino acids in the V2 region and the duration of the infection (n = 213, r = 0.60, p < 0.001) as well as between CD4+ T cell count and the number of amino acids in the V2 region in 75 patients (n = 213, r = −0.17, p < 0.05). Additionally, there were significant inverse correlations between the FLs and the CD4+ T cell counts > 100/μL in all 65 patients (r = −0.30, p = 0.001) and between the sequence length of the env gene and the CD4+ T cell counts > 100/μL before cART (r = −0.34, p < 0.001) (Table S1).

4. Discussion

This study provides evidence suggesting that the sequence length increases over time, as demonstrated by extensive sequence data collected from 20 patients with a well-known history of hemophilia and KSB-infected local controls. Owing to well-documented primary HIV-1 infection [1,2,3,4,5,6,7] and extensive sequence data obtained for over 30 years from 20 HPs with a common source of HIV-1, we were able to confirm that the FL significantly increases over time; the strength of the correlation was stronger in 20 HPs (Figure 2D) than that in local controls (Figure 2E). Otherwise stated, the more homogeneous the cohort, the higher ts the strength of the correlation.
To our knowledge, this is the first report on the association between the sequence length in FLs and the duration of infection. At any given time, viral populations will be dominated by those strains that are the best fit at that time [27]. However, several reports have focused on V1 and V2 elongation in an elite controller [28,29,30,31] and on the consistent usage of the CCR5 coreceptor [32]. The sequence length of the nef gene in HIV-1, a virulence gene, is significantly shortened compared to that of other genes SIVsm [33] (Figure S2), implying pathogenicity [34].
It is important to reemphasize that viruses adopt a symbiotic pathway rather than killing the host or evolving toward attenuated pathogenicity [27]. In fact, the replicative capacity of HIV-1 in the 2000s was significantly lower than that of the virus in the 1980s [35]. Consequently, the survival period was longer in people infected in the 2000s; however, various factors might be involved. The order of replicative fitness by subtype is as follows: D > A > C [36]. Based on this, the increase in sequence length could be considered as a strategy or evolutionary direction for the virus to adapt under the pressure exerted by the immune system of the host and to coexist with the host.
The increase in sequence length implies that the longer the coding sequence (CDS) length, the lower is the density of ribosomes, resulting in less efficient protein synthesis [37] and fewer virus copies. This may be why an increase in sequence length, such as V1 or V2 elongations, occurs in the elite controller or LTNP [27,28,29], as observed in this study. In the same context, the prognosis of patients with SIVsm and HIV-2 exhibiting longer CDSs is better than that of patients with HIV-1. The increase in CDSs in HIV-1 over time is likely to be related to an attenuated pathogenicity and a different evolutionary direction.
This study had a few limitations. First, sera sampled from different individuals at various time points were used, and the time lag between the primary infection and sampling was particularly long (over 10 years in four HPs). Second, owing to the lack of sample availability, before 2000 and after 2000, sequences were obtained from sera and PBMCs, respectively. However, to our knowledge, the sequence length of HIV-1 is not affected by the sample type, although the sequence identity is affected by G-to-A hypermutations due to ApoBec3 proteins [38]. In addition, we confirmed that there was no difference in sequence length of the nef gene in the samples used for RT-PCR and nested-PCR, which were obtained on the same day. Therefore, we do not think that the two distinct methodologies have impacted the data or the conclusions.
In this study, the sequence length of HIV-1 increased by 4.4 nucleotides per year before cART. When viruses were transmitted to another patient, among quasispecies, the possibility of infection by HIV-1 with a shorter sequence might be higher than that by HIV-1 with long sequences. Thus, at the population level, the accumulation effect of the increase in sequence length might be slower than that at the individual level because of the bottleneck effect. Our novel data suggest that the increase in CDSs over time might indicate an evolutionary direction.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/v13050806/s1, Figure S1. Comparison of the increase in the length of the full-length HIV-1 sequence before cART and upon cART administration in 16 patients with hemophilia; Figure S2. Comparison of the length of the Nef proteins of HIV-1, SIVcpz, HIV-2, and SIVsm; Table S1. Characteristics of the 62 patients determined near full-length sequences of KSB HIV-1 and 3 non-KSB patients; Table S2. Signature pattern residues in the Pol proteins of 96 Korean subclade of HIV-1 subtype B; Table S3. Comparison of the frequency of the signature pattern nucleotides in 62 KSB-infected patients.

Author Contributions

Y.-K.C. designed the experiments, Y.-K.C. and J.-E.K. performed the experiments, Y.-K.C., and B.T.F. analyzed the data and wrote the paper. All authors—reviewed and edited. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant from the Korean Society of Ginseng (2020).

Institutional Review Board Statement

The institutional review board of the Asan Medical Center approved the study protocol (Code 2012-0390, 4 June 2012). The study was conducted in accordance with the Declaration of Helsinki.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

All data are available from the corresponding author on request.

Acknowledgments

We thank the patients with hemophilia for their cooperation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cho, Y.K.; Foley, B.T.; Sung, H.; Kim, Y.B.; Kim, J.H. Molecular epidemiologic study of a human immunodeficiency virus 1 outbreak in hemophiliacs B infected through clotting factor 9 after 1990. Vox Sang. 2007, 92, 113–120. [Google Scholar] [CrossRef]
  2. Cho, Y.K.; Jung, Y.S.; Foley, B.T. Phylogenetic analysis of full-length pol gene from Korean hemophiliacs and plasma donors infected with Korean subclade B of HIV-1. AIDS Res. Hum. Retrovir. 2011, 27, 613–621. [Google Scholar] [CrossRef] [PubMed]
  3. Cho, Y.K.; Jung, Y.; Lee, J.S.; Foley, B.T. Molecular evidence of HIV-1 transmission in 20 Korean individuals with hemophilia; phylogenetic analysis of vif gene. Haemophilia 2012, 18, 291–299. [Google Scholar] [CrossRef]
  4. Cho, Y.K.; Kim, J.E.; Foley, B.T. Phylogenetic analysis of the earliest nef gene from hemophiliacs and local controls in Korea. BioRes. Open Access 2012, 1, 41–49. [Google Scholar]
  5. Cho, Y.K.; Kim, J.E.; Jeong, D.; Foley, B.T. Signature pattern analysis for the full-length env gene of the earliest Korean subclade B of HIV-1: Outbreak among Korean hemophiliacs. Virus Genes 2017, 53, 789–796. [Google Scholar] [CrossRef] [PubMed]
  6. Cho, Y.K.; Kim, J.E.; Foley, B.T. Genetic analysis of the full-length gag gene from the earliest Korean subclade B of HIV-1: An outbreak among Korean hemophiliacs. Viruses 2019, 11, 545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Cho, Y.K.; Kim, J.E.; WOO, J.H. Korean Red Ginseng increases defective pol gene in peripheral blood mononuclear cells of HIV-1-infected patients; inhibition of its detection during ginseng-based combination therapy. J. Ginseng Res. 2019, 43, 684–691. [Google Scholar] [CrossRef] [PubMed]
  8. Leitner, T.; Escanilla, D.; Franzén, C.; Uhlén, M.; Albert, J. Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis. Proc. Natl. Acad. Sci. USA 1996, 93, 10864–10869. [Google Scholar] [CrossRef] [Green Version]
  9. Novitsky, V.; Moyo, S.; Lei, Q.; DeGruttola, V.; Essex, M. Importance of viral sequence length and number of variable and informative sites in analysis of HIV clustering. AIDS Res. Hum. Retrovir. 2015, 31, 531–542. [Google Scholar] [CrossRef] [Green Version]
  10. Daniels, R.S.; Kang, C.; Patel, D.; Xiang, Z.; Douglas, N.W.; Zheng, N.N.; Cho, H.W.; Lee, J.S. An HIV type 1 Subtype B founder effect in Korea: gp160 signature patterns infer circulation of CTL-escape strains at the population level. AIDS Res. Hum. Retrovir. 2003, 19, 631–641. [Google Scholar] [CrossRef]
  11. Korber, B.; Myers, G. Signature pattern analysis: A method for assessing viral sequence relatedness. AIDS Res. Hum. Retrovir. 1992, 8, 1549–1560. [Google Scholar] [CrossRef] [PubMed]
  12. Korber, B.T.; Foley, B.; Gaschen, B.; Kuiken, C. Epidemiological and Immunological Implications of the Global Variability of HIV-1; Pataleo, G., Walker, B.D., Eds.; Humana Press: Totowa, NJ, USA, 2001; pp. 1–32. ISBN 978-1-61737-128-8. [Google Scholar]
  13. Junqueira, D.M.; Almeida, S.E. HIV-1 subtype B: Traces of a pandemic. Virology 2016, 495, 173–184. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, M.S.; Jang, S.Y.; Park, C.S.; Lee, K.M.; Lee, D.H.; Lee, C.H. Timing and evolution of the most recent common ancestor of the Korean clade HIV subtype B based on nef and vif sequences. J. Microbiol. 2009, 47, 85–90. [Google Scholar] [CrossRef]
  15. Easterbrook, P.J.; Smith, M.; Mullen, J.; O’Shea, S.; Chrystie, I.; de Ruiter, A.; Tatt, L.D.; Geretti, A.M.; Zuckerman, M. Impact of HIV-1 viral subtype on disease progression and response to antiretroviral therapy. J. Int. AIDS Soc. 2010, 13, 4. [Google Scholar] [CrossRef] [Green Version]
  16. Cho, Y.K.; Kim, J.E.; Foley, B.T. Phylogenetic analysis of near full-length HIV-1 genomic sequences from 21 Korean individuals. AIDS Res. Hum. Retrovir. 2013, 29, 738–743. [Google Scholar] [CrossRef] [Green Version]
  17. Cho, Y.K.; Sung, H.; Bae, I.G.; Oh, H.B.; Kim, N.J.; Woo, J.H.; Kim, Y.B. Full sequence of HIV type 1 Korean subtype B in an AIDS case with atypical seroconversion: TAAAA at TATA box. AIDS Res. Hum. Retrovir. 2005, 21, 961–964. [Google Scholar] [CrossRef] [PubMed]
  18. Nguyen, L.T.; Schmidt, H.A.; Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Bio. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  19. Ou, C.Y.; Ciesielski, C.A.; Myers, G.; Bandea, C.I.; Luo, C.C.; Korber, B.T.; Mullins, J.I.; Schochetman, G.; Berkelman, R.L.; Economou, A.N.; et al. Molecular epidemiology of HIV transmission in a dental practice. Science 1992, 256, 1165–1171. [Google Scholar] [CrossRef] [Green Version]
  20. Park, C.S.; Kim, M.S.; Lee, S.D.; Kim, S.S.; Lee, K.M.; Lee, C.H. Molecular phylogenetic analysis of HIV-1 vif gene from Korean isolates. J. Microbiol. 2006, 44, 655–659. [Google Scholar]
  21. Park, C.S.; Lee, D.H.; Lee, K.M.; Lee, C.H. Characterization and signature pattern analysis of Korean clade HIV-1 using nef gene sequences. J. Microbiol. 2008, 46, 88–94. [Google Scholar] [CrossRef]
  22. Cho, Y.K.; Kim, J.E.; Woo, J.H. Genetic defects in the nef gene are associated with Korean Red Ginseng intake: Monitoring of nef 452 sequence polymorphisms over 20 years. J. Ginseng Res. 2017, 41, 144–150. [Google Scholar] [CrossRef] [Green Version]
  23. Cho, Y.K.; Kim, J.E.; Lee, S.U.; Foley, B.T.; Choi, B.S. Impact of HIV-1 subtype and Korean Red Ginseng on AIDS progression: Comparison of subtype B and subtype D. J. Ginseng Res. 2019, 43, 312–319. [Google Scholar] [CrossRef] [PubMed]
  24. Shankarappa, R.; Margolick, J.B.; Gange, S.J.; Rodrigo, A.G.; Upchurch, D.; Farzadegan, H.; Gupta, P.; Rinaldo, C.R.; Learn, G.H.; He, X.; et al. Consistent viral evolutionary changes associated with the progression of human immunodeficiency virus type 1 infection. J. Virol. 1999, 73, 10489–10502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Sagar, M.; Wu, X.; Lee, S.; Overbaugh, J. Human immunodeficiency virus type 1 V1–V2 envelope loop sequences expand and add glycosylation sites over the course of infection, and these modifications affect antibody neutralization sensitivity. J. Virol. 2006, 80, 9586–9598. [Google Scholar] [CrossRef] [Green Version]
  26. Bunnik, E.M.; Euler, Z.; Welkers, M.R.; Boeser-Nunnink, B.D.; Grijsen, M.L.; Prins, J.M.; Schuitemaker, H. Adaptation of HIV-1 envelope gp120 to humoral immunity at a population level. Nat. Med. 2010, 16, 995–997. [Google Scholar] [CrossRef]
  27. Malim, M.H.; Emerman, M. HIV-1 sequence variation, drift, shift, and attenuation. Cell 2001, 104, 469–472. [Google Scholar] [CrossRef] [Green Version]
  28. Silver, Z.A.; Dickinson, G.M.; Seaman, M.S.; Desrosiers, R.C. A highly unusual V1 region of Env in an elite controller of HIV infection. J. Virol. 2019, 93, e00094-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Shioda, T.; Oka, S.; Xin, X.; Liu, H.; Harukuni, R.; Kurotani, A.; Fukushima, M.; Shiino, T.; Takebe, Y.; Lwamoto, A.; et al. In vivo sequence variability of human immunodeficiency virus type 1 envelope gp120: Association of V2 extension with slow disease progression. J. Virol. 1997, 71, 4871–4881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Wang, B.; Spira, T.J.; Owen, S.; Lal, R.B.; Saksena, N.K. HIV-1 strains from a cohort of American subjects reveal the presence of a V2 region extension unique to slow progressors and non-progressors. AIDS 2000, 14, 213–223. [Google Scholar] [CrossRef]
  31. Masciotra, S.; Owen, S.M.; Rudolph, D.; Yang, C.; Wang, B.; Saksena, N.; Spira, T.; Dhawan, S.; Lal, R.B. Temporal relationship between V1V2 variation, macrophage replication, and coreceptor adaptation during HIV-1 disease progression. AIDS 2000, 16, 1887–1898. [Google Scholar] [CrossRef]
  32. Daniels, R.S.; Wilson, P.; Patel, D.; Longhurst, H.; Pattersom, S. Analysis of full-length HIV type 1 env genes indicates differences between the virus infecting T cells and dendritic cells in peripheral blood of infected patients. AIDS Res. Hum. Retrovir. 2004, 20, 409–413. [Google Scholar] [CrossRef]
  33. Hirao, K.; Andrews, S.; Kuroki, K.; Kusaka, H.; Tadokoro, T.; Kita, S.; Ose, T.; Rowland-Jones, S.L.; Maenaka, K. Structure of HIV-2 Nef reveals features distinct from HIV-1 involved in immune regulation. iScience 2020, 23, 100758. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Schindler, M.; Münch, J.; Kutsch, O.; Li, H.; Santiago, M.L.; Bibollet-Ruche, F.; Müller-Trutwin, M.C.; Novembre, F.J.; Peeters, M.; Courgnaud, V.; et al. Nef-mediated suppression of T cell activation was lost in a lentiviral lineage that gave rise to HIV-1. Cell 2006, 125, 1055–1067. [Google Scholar] [CrossRef] [Green Version]
  35. Ariën, K.K.; Troyer, R.M.; Gali, Y.; Colebunders, R.L.; Arts, E.J.; Vanham, G. Replicative fitness of historical and recent HIV-1 isolates suggests HIV-1 attenuation over time. AIDS 2005, 19, 1555–1564. [Google Scholar] [CrossRef] [Green Version]
  36. Venner, C.M.; Nankya, I.; Kyeyune, F.; Demers, K.; Kwok, C.; Chen, P.L.; Rwambuya, S.; Munjoma, M.; Chipato, T.; Byamugisha, J.; et al. Infecting HIV-1 subtype predicts disease progression in women of sub-Saharan Africa. EBioMedicine 2016, 13, 305–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Fernandes, L.D.; de Moura, A.P.S.; Ciandrini, L. Gene length as a regulator for ribosome recruitment and protein synthesis: Theoretical insights. Sci. Rep. 2017, 7, 17409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Okada, A.; Iwatani, Y. APOBEC3G-Mediated G-to-A hypermutation of the HIV-1 Genome: The missing link in antiviral molecular mechanisms. Front. Microbiol. 2016, 7, 2027. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Phylogenetic tree analysis of the earliest near full-length sequences (about 8615 bp from the gag to the nef genes) of 64 Korean patients infected with HIV-1 performed using IQ TREE with 1000 bootstrap replicates: 20 patients with hemophilia (HPs); 3 plasma donors (O, P, and R); 39 local controls infected with the Korean subclade of HIV-1 subtype B (KSB) and 2 non-KSB-infected patients. The upper 103 sequences belonged to KSB, and 2 sequences (05CSR3 and 04KBH8) belonged to subtypes B and D, respectively. In total, 9 and 13 patients, including donors O (Cluster O: donor O, 1–4, 6, 8, 10, and 18, as designated by the taxa in red) and P (Cluster P, as designated by the taxa in red), were strongly clustered within the KSB-infected local controls. The two digits before the ID of each patient and the one or two digits after the ID of each patient denote the year and month of sampling, respectively. The bootstrap values of the nodes for clusters O, P, and KSB (designated by bold in red) were 100%, as determined by 1000 bootstrap replicates.
Figure 1. Phylogenetic tree analysis of the earliest near full-length sequences (about 8615 bp from the gag to the nef genes) of 64 Korean patients infected with HIV-1 performed using IQ TREE with 1000 bootstrap replicates: 20 patients with hemophilia (HPs); 3 plasma donors (O, P, and R); 39 local controls infected with the Korean subclade of HIV-1 subtype B (KSB) and 2 non-KSB-infected patients. The upper 103 sequences belonged to KSB, and 2 sequences (05CSR3 and 04KBH8) belonged to subtypes B and D, respectively. In total, 9 and 13 patients, including donors O (Cluster O: donor O, 1–4, 6, 8, 10, and 18, as designated by the taxa in red) and P (Cluster P, as designated by the taxa in red), were strongly clustered within the KSB-infected local controls. The two digits before the ID of each patient and the one or two digits after the ID of each patient denote the year and month of sampling, respectively. The bootstrap values of the nodes for clusters O, P, and KSB (designated by bold in red) were 100%, as determined by 1000 bootstrap replicates.
Viruses 13 00806 g001
Figure 2. The full-length (FL) sequence length increases over time and the sequence identity decreases over time in KSB-infected patients. (AC): Sequence identity and FL sequences according to the interval since the outbreak in January 1990 in 20 HPs. (A) Correlation between the intervals from the outbreak to sampling and the sequence identity of each earliest FL sequence before cART in the 20 patients with hemophilia, compared with the corresponding earliest sequence of the plasma donors. (B) The correlation coefficient (CC), γ, was 0.77 for the 39 FLs before cART (P < 0.001). (C) The CC, γ, was −0.17 upon cART. (DF), Significant increase in the FL sequence length by sampling year before cART, in patients with hemophilia (D), 42 KSB-infected local controls (n = 60) (E), and all the 62 KSB-infected patients (F).
Figure 2. The full-length (FL) sequence length increases over time and the sequence identity decreases over time in KSB-infected patients. (AC): Sequence identity and FL sequences according to the interval since the outbreak in January 1990 in 20 HPs. (A) Correlation between the intervals from the outbreak to sampling and the sequence identity of each earliest FL sequence before cART in the 20 patients with hemophilia, compared with the corresponding earliest sequence of the plasma donors. (B) The correlation coefficient (CC), γ, was 0.77 for the 39 FLs before cART (P < 0.001). (C) The CC, γ, was −0.17 upon cART. (DF), Significant increase in the FL sequence length by sampling year before cART, in patients with hemophilia (D), 42 KSB-infected local controls (n = 60) (E), and all the 62 KSB-infected patients (F).
Viruses 13 00806 g002
Figure 3. The length of FLs increased over time in 3 long-term nonprogressors (LTNPs), whereas it decreased in subtype D-infected patients. The correlation coefficient between the sequence length and the duration of infection was significant in 3 LTNPs over 25 years (up to 294, 307, and 286 months). Patients infected with CRF02_AG (A), KSB (B), and subtype B (C) were diagnosed with HIV-1 infection in 1988, 1992, and 1987, respectively. Interestingly, the sequence length in patients infected with subtype B significantly decreased from 8697 ± 26 (n = 6), just before 286 months (March 2011), to 8623 ± 24 (n = 6; marked as black stars), after 286 months, when the plasma RNA copy number significantly increased (p < 0.001) [22]. The sequence length significantly decreased after 307 months. In contrast, it was inversely significant in patient 93-01, who was diagnosed in December 1992 (D). The first sexual contact and diagnosis of pulmonary tuberculosis of the patient were observed in 1988 and 1989, respectively. He was treated with Korean Red Ginseng (12,720 g) from April 1993 to August 2004 [23]. Thus, despite the most rapidly progressing subtype D infection, he remained healthy for at least 12 years [23]. The correlation analysis did not include the data marked as stars and as a white circle on cART.
Figure 3. The length of FLs increased over time in 3 long-term nonprogressors (LTNPs), whereas it decreased in subtype D-infected patients. The correlation coefficient between the sequence length and the duration of infection was significant in 3 LTNPs over 25 years (up to 294, 307, and 286 months). Patients infected with CRF02_AG (A), KSB (B), and subtype B (C) were diagnosed with HIV-1 infection in 1988, 1992, and 1987, respectively. Interestingly, the sequence length in patients infected with subtype B significantly decreased from 8697 ± 26 (n = 6), just before 286 months (March 2011), to 8623 ± 24 (n = 6; marked as black stars), after 286 months, when the plasma RNA copy number significantly increased (p < 0.001) [22]. The sequence length significantly decreased after 307 months. In contrast, it was inversely significant in patient 93-01, who was diagnosed in December 1992 (D). The first sexual contact and diagnosis of pulmonary tuberculosis of the patient were observed in 1988 and 1989, respectively. He was treated with Korean Red Ginseng (12,720 g) from April 1993 to August 2004 [23]. Thus, despite the most rapidly progressing subtype D infection, he remained healthy for at least 12 years [23]. The correlation analysis did not include the data marked as stars and as a white circle on cART.
Viruses 13 00806 g003
Figure 4. The sequence length increases over time in subtype-B-infected patients. A total of 64 FLs (subtype B) were randomly selected from Los Alamos National Laboratory Database over 16 years from 1983–1998 (evenly 4 per year). They were obtained before cART. The correlation between the sampling year and the sequence length was significant for FLs (A), env (B), and nef (C) genes. Of the 64 sequences, 44 originated from the United States, and if these alone are analyzed, the correlation is even higher (in order, r = 0.58, p < 0.001; r = 0.40, p < 0.01; and r = 0.46, p < 0.01 for FLs, env, and nef genes, respectively).
Figure 4. The sequence length increases over time in subtype-B-infected patients. A total of 64 FLs (subtype B) were randomly selected from Los Alamos National Laboratory Database over 16 years from 1983–1998 (evenly 4 per year). They were obtained before cART. The correlation between the sampling year and the sequence length was significant for FLs (A), env (B), and nef (C) genes. Of the 64 sequences, 44 originated from the United States, and if these alone are analyzed, the correlation is even higher (in order, r = 0.58, p < 0.001; r = 0.40, p < 0.01; and r = 0.46, p < 0.01 for FLs, env, and nef genes, respectively).
Viruses 13 00806 g004
Table 1. Primer sequences used for the polymerase chain reaction.
Table 1. Primer sequences used for the polymerase chain reaction.
Nested PCRPrimerSequences (5′-3′)
For the full-length gag gene [6]
First PCR503k5′-CCKTCTGTTGTGTGACTCTGGTAA-3′forward
5245′-CATTGTTTAACTTTTGGGCCATCC -3′reverse
Second PCR504F5′-TCTCTAGCAGTGGCGCCCGAAC-3′forward
5055′-GAGACATGGGTGCGAGAGCGT-3′forward
5225′-ACTGTCCTACTTTGATAAAACCTC-3′reverse
For the full-length pol gene [2,7]
First PCRHXB25′-GTGGGAGAAATCTATAAAAGATGG-3′forward
OBP25′-GAGACTCCCTGACCCAGATG-3′reverse
OBP2k5′-GAGACTCCCTGACCCAGATG-3′reverse
5505′-CCTAGTGGGATGTGTACTTCTGAA-3′reverse
Second PCRPO15′-AAAATTGCAGGGCCCCTAGGA-3′forward
PR3-15′-GAAGCAGGAGCCGATAGACA-3′forward
OBP45′-CAATCATCACCTGCCATCTG-3′reverse
P25′-AGGAAGGACACCAAATGAAAG-3′forward
P165′-GGATKAGTGCTTTCATAGTGA-3′reverse
For vif, vpr, tat, rev, and vpu genes
First PCR5455′-GCAGTACAAATGGCAGTATTCATC-3′forward
KMK25′-ATGGGAATTGGTTCAAAGGA-3′reverse
Second PCR548F5′-AGTGACATAAAAGTAGTRCCAAGAA-3′forward
LA1065′-TTCACTCTCATTGCCACT-3′reverse
For the full-length env gene [5]
First PCROWE15′-TCATCAAGTTTCTCTATCAAAGCA-3′forward
OWE25′-TCTGACTGGAAAGCCCACTT-3′reverse
Second PCROWE35′-GCAATATTAGCAATAGTTGTGTGG-3′forward
OWE45′-ATACTGCTCCCACCCCTTCT-3′reverse
For nef gene [4]
First PCRNef5′55′-AGGATTGTGGAACTTCTGGGAC-3′forward
LTR35′-AGGCTCAGATCTGGTCTAAC-3′reverse
Second PCRNef35′-ATGGGTGGCAAGTGGTCAAA-3′forward
N105′-CGTCCAGAATTCGGAAAGTCCCCAGCGGAAAGT-3′reverse
Table 2. Changes of the sequence length of HIV-1 over 110 ± 43 months by genes in patients with hemophilia.
Table 2. Changes of the sequence length of HIV-1 over 110 ± 43 months by genes in patients with hemophilia.
Sequences with IntervalsGagPolVifVprTatRevVpuEnvNefFL
First sequences
(n = 16)
1503 ± 03012 ± 057929121576247 ± 22563 ± 19622.5 ± 28620 ± 20
Second sequences
(n = 16)
1504 ± 33011 ± 357929121576247 ± 22595 ± 22627.2 ± 88659 ± 23
p-valueNSNSNSNSNSNSNS< 0.001< 0.05= 0.0001
There was an interval of approximately 110 ± 43 months between the sequences obtained at the two time points. NS = not significant. p-value calculated using a paired t-test. Data are expressed by mean ± standard deviation.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cho, Y.-K.; Kim, J.-E.; Foley, B.T. Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years. Viruses 2021, 13, 806. https://doi.org/10.3390/v13050806

AMA Style

Cho Y-K, Kim J-E, Foley BT. Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years. Viruses. 2021; 13(5):806. https://doi.org/10.3390/v13050806

Chicago/Turabian Style

Cho, Young-Keol, Jung-Eun Kim, and Brian T. Foley. 2021. "Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years" Viruses 13, no. 5: 806. https://doi.org/10.3390/v13050806

APA Style

Cho, Y. -K., Kim, J. -E., & Foley, B. T. (2021). Sequence Length of HIV-1 Subtype B Increases over Time: Analysis of a Cohort of Patients with Hemophilia over 30 Years. Viruses, 13(5), 806. https://doi.org/10.3390/v13050806

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

Article Metrics

Back to TopTop