Subtype-Dependent Co-receptor Tropism in Cuban HIV-1–Infected Patients: Implications for Maraviroc Treatment

Yenisleidys Martínez-Montesino Vivian Kourí-Cardellá Lissette Pérez-Santos Rui Han Yanet Pintos-Saavedra Yoan Alemán-Campos Yudira Soto-Brito Yoanna Baños-Morales Yaniris Caturla-Fernández About the authors

Abstract

INTRODUCTION Unlike most high-income countries where subtype B viruses predominate, the Cuban HIV-1 epidemic is characterized by a great diversity of subtypes and circulating recombinant forms. Some studies have shown that HIV variants exhibiting a preference for the CXCR4 co-receptor (X4-tropic) could have impacts on disease pathogenesis, with clinical implications for antiviral treatment plans. Determination of HIV co-receptor tropism is crucial for clinicians in deciding whether maraviroc is an appropriate antiviral.

OBJECTIVE Characterize V3 sequence variability and its relation to viral tropism across different subtypes circulating in Cuba and explore how this may affect treatment success with maraviroc.

METHODS We designed a cross-sectional study that included 72 plasma samples obtained at the Pedro Kourí Tropical Medicine Institute in Havana, Cuba. We sequenced the C2V3 env region and assessed subtype based both on env and pol sequences; tropism was predicted by Geno2pheno analysis. Additionally, 35 V3-loop Cuban sequences, obtained from a previous study, were incorporated into the analysis. Statistical associations among virological, clinical and epidemiological variables were assessed by a chi-square test.

RESULTS Tropism prediction for 72 variants revealed that CRF19_cpx was associated with dual-tropic R5X4 viruses (p = 0.034). Moreover, when 35 sequences from a former study were added, the association was significant not only for R5X4 (p = 0.019) but also for X4-tropic variants (p = 0.044). Alignment of 107 V3-loop sequences showed wide diversity among the different HIV-1 subtypes circulating in Cuba.

CONCLUSIONS In accordance with G2P, CRF19_cpx is a genetic variant with a high proportion of X4 and R5X4-tropic viruses. The results from the present study suggest that the Cuban recombinant could be a more pathogenic variant and that maraviroc may not be suitable for patients infected with CRF19_cpx.

Keywords
HIV; CCR5 receptor antagonists; maraviroc; HIV envelope glycoprotein gp120; Cuba

INTRODUCTION

During entry into the host cell, HIV-1 surface glycoprotein gp120 interacts with the CD4 receptor and, generally, with one of two chemokine co-receptors: CCR5 (chemokine [C-C] motif receptor 5) or CXCR4 (CXC chemokine receptor 4). This interaction defines viral tropism in R5-tropic or X4-tropic strains, respectively.[11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,22. Pérez-Álvarez L, Delgado E, Vega Y, Montero V, Cuevas T, Fernández-García A, et al. Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections. J Antimicrob Chemother [Internet]. 2014 Jan [cited 2020 Dec 2];69(1):246–53. Available at: https://pdfs.semanticscholar.org/3ada/c6787c3af51dfa8b2b51764550cfe6f57608.pdf?_ga=2.176869410.60519488.1608083556-726339794.1548872336
https://pdfs.semanticscholar.org/3ada/c6...
] Co-receptor selectivity is determined by genetic sequences within the HIV-1 gp120, particularly in a highly variable and structurally flexible region called the ‘V3 loop’, which is involved in co-receptor binding.[33. Soulie C, Morand-Joubert L, Cottalorda J, Charpentier C, Bellecave P, Le Guen L, et al. Performance of genotypic algorithms for predicting tropism for HIV-1 CRF01_AE recombinant. J Clin Virol [Internet]. 2018 Feb–Mar [cited 2020 Dec 2];99–100:57–60. Available at: https://linkinghub.elsevier.com/retrieve/pii/S1386-6532(17)30354-2
https://linkinghub.elsevier.com/retrieve...
,44. Panos G, Watson DC. Effect of HIV-1 subtype and tropism on treatment with chemokine coreceptor entry inhibitors; overview of viral entry inhibition. Crit Rev Microbiol [Internet]. 2015 [cited 2020 Dec 2];41(4):473–87. Available at: https://www.tandfonline.com/doi/full/10.3109/1040841X.2013.867829?scroll=top&needAccess=true
https://www.tandfonline.com/doi/full/10....
]

Viral tropism has been proposed as an influence on HIV-1 pathogenesis and replication, with implications for treatment. R5-tropic variants have shown predominance during early stages of the disease, are generally not syncytium-inducing and have low replication capacity.[55. Woollard SM, Kanmogne GD. Maraviroc: a review of its use in HIV infection and beyond. Drug Des Devel Ther [Internet]. 2015 Oct 1 [cited 2020 Dec 2];9:5447–68. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/pmid/26491256/
https://www.ncbi.nlm.nih.gov/pmc/article...
,66. Schuitemaker H, van 't Wout AB, Lusso P. Clinical significance of HIV-1 coreceptor usage. J Transl Med [Internet]. 2011 Jan 27 [cited 2020 Dec 15];9 Suppl 1 (Suppl 1):S5. Available at: https://translational-medicine.biomedcentral.com/articles/10.1186/1479-5876-9-S1-S5
https://translational-medicine.biomedcen...
] It is believed that during the course of infection, R5 viruses experience a switch in tropism and thus, X4 variants emerge, characterized by syncytium induction and higher replication capacity. In addition, a third group of viruses has been classified as dual-tropic or R5X4-tropic variants, referring to their ability to interact with both CCR5 and CXCR4 co-receptors.[22. Pérez-Álvarez L, Delgado E, Vega Y, Montero V, Cuevas T, Fernández-García A, et al. Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections. J Antimicrob Chemother [Internet]. 2014 Jan [cited 2020 Dec 2];69(1):246–53. Available at: https://pdfs.semanticscholar.org/3ada/c6787c3af51dfa8b2b51764550cfe6f57608.pdf?_ga=2.176869410.60519488.1608083556-726339794.1548872336
https://pdfs.semanticscholar.org/3ada/c6...
,77. Mild M, Esbjornsson J, Fenyo EM, Medstrand P. Frequent intrapatient recombination between human immunodeficiency virus type 1 R5 and X4 envelopes: implications for coreceptor switch. J Virol [Internet]. 2007 Apr [cited 2020 Dec 16];81(7):3369–76. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1866041/pdf/1295-06.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
]

IMPORTANCE

CRF19_cpx is one of the most prevalent HIV-1 subtypes circulating in Cuba. X4- and dual-tropic viruses prevail in CRF19_cpx-infected patients, which suggest this may be a more pathogenic variant and less susceptible to treatment with maraviroc.

Determining viral tropism is important when choosing antivirals like maraviroc, the first CCR5 inhibitor approved for treatment of R5-tropic HIV-1 infection in both treatment-naive and treatment-experienced adult patients. However, maraviroc does not inhibit entry of X4-tropic or dual-tropic HIV-1 viruses and a test that predicts HIV-1 co-receptor use is highly recommended before treatment prescription.[11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,88. Phuphuakrat A, Phawattanakul S, Pasomsub E, Kiertiburanakul S, Chantratita W, Sungkanuparph S. Coreceptor tropism determined by genotypic assay in HIV-1 circulating in Thailand, where CRF01_AE predominates. HIV Med [Internet]. 2014 May [cited 2020 Dec 15];15(5):269–75. Available at: https://onlinelibrary.wiley.com/doi/epdf/10.1111/hiv.12108
https://onlinelibrary.wiley.com/doi/epdf...
,99. Asin-Milan O, Chamberland A, Wei Y, Haidara A, Sylla M, Tremblay CL. Mutations in variable domains of the HIV-1 envelope gene can have a significant impact on maraviroc and vicriviroc resistance. AIDS Res Ther [Internet]. 2013 Jun 7 [cited 2020 Dec 20];10(1):15. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3700831/pdf/1742-6405-10-15.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
]

Two different methods are employed to predict co-receptor use. Phenotypic assays have greater accuracy; however, they are expensive, time-consuming, and require special facilities and trained personnel. On the other hand, genotypic assays rely on bioinformatic tools to predict co-receptor use based mainly on the sequence of the gp120's V3 loop.[1010. Recordon-Pinson P, Soulie C, Flandre P, Descamps D, Lazrek M, Charpentier C, et al. Evaluation of the genotypic prediction of HIV-1 coreceptor use versus a phenotypic assay and correlation with the virological response to maraviroc: the ANRS GenoTropism study. Antimicrob Agents Chemother [Internet]. 2010 Aug [cited 2020 Dec 3];54(8):3335–40. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2916345/pdf/0148-10.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,1111. Poveda E, Alcami J, Paredes R, Córdoba J, Gutiérrez F, Llibre JM, et al. Genotypic determination of HIV tropism - clinical and methodological recommendations to guide the therapeutic use of CCR5 antagonists. AIDS Rev [Internet]. 2010 [cited 2020 Dec18];12(3):135–48. Available at: https://www.researchgate.net/publication/230569071_Genotypic_Determination_of_HIV_TropismClinical_and_methodological_recommendations_to_guide_the_therapeutic_use_of_CCR5_antagonists/link/00b4952ea736ab0cf1000000/download
https://www.researchgate.net/publication...
] These latter methods are rapid and easier to use; however, their development has been based mostly on sequences from HIV-1 subtype B isolates, without considering the viral diversity of the epidemic within countries.[11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
22. Pérez-Álvarez L, Delgado E, Vega Y, Montero V, Cuevas T, Fernández-García A, et al. Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections. J Antimicrob Chemother [Internet]. 2014 Jan [cited 2020 Dec 2];69(1):246–53. Available at: https://pdfs.semanticscholar.org/3ada/c6787c3af51dfa8b2b51764550cfe6f57608.pdf?_ga=2.176869410.60519488.1608083556-726339794.1548872336
https://pdfs.semanticscholar.org/3ada/c6...
]

In Cuba, the HIV-1 epidemic is atypical when compared to the rest of the Latin American and Caribbean region. Cuba has a great diversity of HIV-1 subtypes including many recombinant forms (CRFs) such as BG recombinants (CRF20_BG, CRF23_BG and CRF24_BG) and complex CRFs, such as CRF18_cpx and CRF19_cpx.[1515. Kourí V, Khouri R, Alemán Y, Abrahantes Y, Vercauteren J, Pineda-Pena AC, et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine [Internet]. 2015 Jan 28 [cited 2020 Dec 16];2(3):244–54. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4484819/pdf/main.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,1616. Pérez L, Kouri V, Alemán Y, Abrahantes Y, Correa C, Aragones C, et al. Antiretroviral drug resistance in HIV-1 therapy-naive patients in Cuba. Infect Genet Evol [Internet]. 2013 Jun [cited 2020 Dec 15];16:144–50. Available at: https://www.sciencedirect.com/science/article/abs/pii/S1567134813000361
https://www.sciencedirect.com/science/ar...
]

A Cuban study identified that in patients harboring CRF19_cpx viruses, CXCR4-using variants prevailed and disease progression to AIDS occurred more rapidly than in patients harboring subtype B viruses.[1515. Kourí V, Khouri R, Alemán Y, Abrahantes Y, Vercauteren J, Pineda-Pena AC, et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine [Internet]. 2015 Jan 28 [cited 2020 Dec 16];2(3):244–54. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4484819/pdf/main.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] An association between HIV-1 subtype and co-receptor use has also been demonstrated in other studies.[22. Pérez-Álvarez L, Delgado E, Vega Y, Montero V, Cuevas T, Fernández-García A, et al. Predominance of CXCR4 tropism in HIV-1 CRF14_BG strains from newly diagnosed infections. J Antimicrob Chemother [Internet]. 2014 Jan [cited 2020 Dec 2];69(1):246–53. Available at: https://pdfs.semanticscholar.org/3ada/c6787c3af51dfa8b2b51764550cfe6f57608.pdf?_ga=2.176869410.60519488.1608083556-726339794.1548872336
https://pdfs.semanticscholar.org/3ada/c6...
,66. Schuitemaker H, van 't Wout AB, Lusso P. Clinical significance of HIV-1 coreceptor usage. J Transl Med [Internet]. 2011 Jan 27 [cited 2020 Dec 15];9 Suppl 1 (Suppl 1):S5. Available at: https://translational-medicine.biomedcentral.com/articles/10.1186/1479-5876-9-S1-S5
https://translational-medicine.biomedcen...
,1717. Fenyö EM, Esbjörnsson J, Medstrand P, Jansson M. Human immunodeficiency virus type 1 biological variation and coreceptor use: from concept to clinical significance. J Intern Med [Internet]. 2011 Dec [cited 2020 Dec 17];270(6):520–31. Available at: https://onlinelibrary.wiley.com/doi/epdf/10.1111/j.1365-2796.2011.02455.x
https://onlinelibrary.wiley.com/doi/epdf...
]

The objective of this study was to characterize V3 sequence variability across different Cuban HIV-1 subtypes and its relation to viral tropism in the context of the Cuban epidemic, in order to achieve a more complete understanding of HIV-1 pathogenesis and improved management of antiviral treatment for Cuban patients.

METHODS

Study Design We analyzed additional plasma samples from Cuban HIV-1–infected patients in order to extend our previous analysis on prevalent subtypes[1818. Kourí V, Alemán Y, Díaz D, Pérez L, Limia CM, Soto Y, et al. Co-receptor tropism determined by genotypic assay in HIV-1 circulating in Cuba. J AIDS Clin Res [Internet]. 2016 [cited 2020 Dec 2];7(7):592. Available at: https://www.hilarispublisher.com/open-access/coreceptor-tropism-determined-by-genotypic-assay-in-hiv1-circulating-in-cuba-2155-6113-1000592.pdf
https://www.hilarispublisher.com/open-ac...
] and include data on viral tropism. Samples were selected from HIV-1–infected patients who were tested for antiretroviral resistance at the Pedro Kourí Tropical Medicine Institute (IPK) in Havana, Cuba, from January 2015 through July 2016. Those patients infected with the following subtypes and recombinant forms (classified regarding pol sequence employed in resistance testing) were selected: subtype B, CRF_BG (for the purpose of this study CRF20_BG, CRF23_BG and CRF24_BG were grouped together as CRF_BG), CRF18_cpx and CRF19_cpx.[1515. Kourí V, Khouri R, Alemán Y, Abrahantes Y, Vercauteren J, Pineda-Pena AC, et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine [Internet]. 2015 Jan 28 [cited 2020 Dec 16];2(3):244–54. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4484819/pdf/main.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,1616. Pérez L, Kouri V, Alemán Y, Abrahantes Y, Correa C, Aragones C, et al. Antiretroviral drug resistance in HIV-1 therapy-naive patients in Cuba. Infect Genet Evol [Internet]. 2013 Jun [cited 2020 Dec 15];16:144–50. Available at: https://www.sciencedirect.com/science/article/abs/pii/S1567134813000361
https://www.sciencedirect.com/science/ar...
,1919. Kourí V, Alemán Y, Pérez L, Pérez J, Fonseca C, Correa C, et al. High frequency of antiviral drug resistance and non-B subtypes in HIV-1 patients failing antiviral therapy in Cuba. J Clin Virol [Internet]. 2012 Dec 5 [cited 2020 Dec 2];55(4):348–55. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4225368/pdf/JIAS-17-19754.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] After selection, a total of 72 patients were included and C2V3 sequences of the glycoprotein gp120 were obtained and employed to assess viral subtype regarding env region and predict co-receptor use.

An additional group of 35 V3-loop sequences from Cuban patients[1818. Kourí V, Alemán Y, Díaz D, Pérez L, Limia CM, Soto Y, et al. Co-receptor tropism determined by genotypic assay in HIV-1 circulating in Cuba. J AIDS Clin Res [Internet]. 2016 [cited 2020 Dec 2];7(7):592. Available at: https://www.hilarispublisher.com/open-access/coreceptor-tropism-determined-by-genotypic-assay-in-hiv1-circulating-in-cuba-2155-6113-1000592.pdf
https://www.hilarispublisher.com/open-ac...
] was included for analysis of association between viral subtype and tropism, to better characterize V3 sequence variability. We decided to include these sequences because they corresponded to a similar study performed by our group from January 2014 through January 2015.[1818. Kourí V, Alemán Y, Díaz D, Pérez L, Limia CM, Soto Y, et al. Co-receptor tropism determined by genotypic assay in HIV-1 circulating in Cuba. J AIDS Clin Res [Internet]. 2016 [cited 2020 Dec 2];7(7):592. Available at: https://www.hilarispublisher.com/open-access/coreceptor-tropism-determined-by-genotypic-assay-in-hiv1-circulating-in-cuba-2155-6113-1000592.pdf
https://www.hilarispublisher.com/open-ac...
] For sequence selection, the same criteria mentioned above were applied. We included a total of 14 sequences from subtype B, 13 from CRF_BG, 2 from CRF18_cpx and 6 from CRF19_cpx, classified according to pol sequence.

Epidemiological, clinical, virological and immunological information Epidemiological and clinical data were collected at sampling as well as from a screening of selected databases stored at IPK (SIDATRAT). Patient viral load and CD4+ cell counts were performed as previously described.[1818. Kourí V, Alemán Y, Díaz D, Pérez L, Limia CM, Soto Y, et al. Co-receptor tropism determined by genotypic assay in HIV-1 circulating in Cuba. J AIDS Clin Res [Internet]. 2016 [cited 2020 Dec 2];7(7):592. Available at: https://www.hilarispublisher.com/open-access/coreceptor-tropism-determined-by-genotypic-assay-in-hiv1-circulating-in-cuba-2155-6113-1000592.pdf
https://www.hilarispublisher.com/open-ac...
]

HIV sequencing and subtype assignment RNA extraction from plasma samples, C2V3 amplification and sequencing were obtained as previously described.[1818. Kourí V, Alemán Y, Díaz D, Pérez L, Limia CM, Soto Y, et al. Co-receptor tropism determined by genotypic assay in HIV-1 circulating in Cuba. J AIDS Clin Res [Internet]. 2016 [cited 2020 Dec 2];7(7):592. Available at: https://www.hilarispublisher.com/open-access/coreceptor-tropism-determined-by-genotypic-assay-in-hiv1-circulating-in-cuba-2155-6113-1000592.pdf
https://www.hilarispublisher.com/open-ac...
] Sequences were edited with Sequencher 4.9 (Gene Codes Coorporation, USA) employing the HXB2 strain (Genbank access number K03455.1, a. Available at: https://www.ncbi.nlm.nih.gov/genbank) as the reference sequence. Subtype was assessed via COMET version 2 (Luxembourg Institute of Health, Luxembourg) and REGA version 3 (Stanford University, USA).[2020. Pineda-Pena AC, Faria NR, Imbrechts S, Libin P, Abecasis AB, Deforche K, et al. Automated subtyping of HIV-1 genetic sequences for clinical and surveillance purposes: performance evaluation of the new REGA version 3 and seven other tools. Infect Genet Evol [Internet]. 2013 Oct [cited 2020 Dec 14];19:337–48. Available at: https://reader.elsevier.com/reader/sd/pii/S1567134813001810?token=224CF28032AD88B452AD006B04DDF3EEA3E86330D469FD9040A10BEC8BE99E5839D54348182D867AB69D2307C8D37009
https://reader.elsevier.com/reader/sd/pi...
,2121. de Oliveira T, Deforche K, Cassol S, Salminen M, Paraskevis D, Seebregts C, et al. An automated genotyping system for analysis of HIV-1 and other microbial sequences. Bioinformatics [Internet]. 2005 Oct 1 [cited 2020 Dec 16];21(19):3797–800. Available at: https://www.researchgate.net/publication/7682559_An_automated_genotyping_system_for_analysis_of_HIV-1_and_other_microbial_sequences
https://www.researchgate.net/publication...
] In addition, subtype assignment was confirmed by manual phylogenetic analysis employing PhyML/One Click software (Institut Pasteur, France).[2222. Lemoine F, Correia D, Lefort V, Doppelt-Azeroual O, Mareuil F, Cohen-Boulakia S, et al. NGPhylogeny.fr: new generation phylogenetic services for non-specialists. Nucleic Acids Res [Internet]. 2019 Apr 27 [cited 2020 Dec 3];47(W1):W260–W5. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6602494/pdf/gkz303.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] Sequences and HIV-1 subtype in the pol region were obtained from antiretroviral resistance testing routinely performed in our laboratory at IPK.[2323. Alemán Y, Vinken L, Kouri V, Pérez L, Álvarez A, Abrahantes Y, et al. Performance of an in-house human immunodeficiency virus type 1 genotyping system for assessment of drug resistance in Cuba. PLoS One [Internet]. 2015 Feb 11 [cited 2020 Dec 14];10(2):e0117176. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4324769/pdf/pone.0117176.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
]

Prediction of co-receptor use Tropism prediction based on V3-loop sequence was performed using the Geno2pheno (G2P) algorithm (Max Planck Institute, Germany). The false positive rate (FPR) that defines the probability of classifying an R5 virus (incorrectly) as an X4 virus was set at a cut-off value of 5% following German guidelines.[1515. Kourí V, Khouri R, Alemán Y, Abrahantes Y, Vercauteren J, Pineda-Pena AC, et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine [Internet]. 2015 Jan 28 [cited 2020 Dec 16];2(3):244–54. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4484819/pdf/main.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,2424. Lengauer T, Sander O, Sierra S, Thielen A, Kaiser R. Bioinformatics prediction of HIV coreceptor usage. Nat Biotechnol [Internet]. 2007 Dec [cited 2020 Dec 4];25(12):1407–10. Available at: https://www.nature.com/articles/nbt1371
https://www.nature.com/articles/nbt1371...
,2525. Jensen MA, van 't Wout AB. Predicting HIV-1 coreceptor usage with sequence analysis. AIDS Rev [Internet]. 2003 Apr–Jun [cited 2020 Dec 4];5(2):104–12. Available at: https://www.researchgate.net/publication/10648055_Predicting_HIV1_coreceptor_usage_with_sequence_analysis
https://www.researchgate.net/publication...
] Patients with FPR ≤5% were considered infected with X4 viruses; patients with FPR ≥20% were predicted to harbor R5 viruses; and patients with FPR values >5% and <20% were considered to have dual-tropic R5X4 virus infections. Additionally, V3 net charge and 11/25 rule were employed in predicting viral tropism.[2626. Seclén E, Garrido C, González M M, González-Lahoz J, de Mendoza C, Soriano V, et al. High sensitivity of specific genotypic tools for detection of X4 variants in antiretroviral-experienced patients suitable to be treated with CCR5 antagonists. J Antimicrob Chemother [Internet]. 2010 Jul [cited 2020 Dec 16];65(7):1486–92. Available at: https://www.researchgate.net/publication/43351431_High_sensitivity_of_specific_genotypic_tools_for_detection_of_X4_variants_in_antiretroviralexperienced_patients_suitable_to_be_treated_with_CCR5_antagonists
https://www.researchgate.net/publication...
]

Statistical analysis Comparison of FPR values among groups was performed using GraphPad Prism version 5 software (GraphPad, USA) and the non-parametric Kruskall Wallis Test and Dunn Test for post-hoc comparisons. Analysis of the association between co-receptor use, subtype, clinical, virological and immunological variables was performed via a chi-square test using SPSS v.22 (IBM, USA). For all comparisons, a p-value ≤0.05 was considered statistically significant.

EthicsStudy protocols were designed in accordance with the Helsinki Declaration and approved by the IPK Ethics Committee.

RESULTS

Study population characteristics The study universe consisted of 72 plasma samples corresponding to the same number of patients. Two patients lacked clinical and epidemiological information, and were therefore only included in the analysis of co-receptor tropism prediction and its association with HIV-1 subtype. For the remaining 70 patients, epidemiological, virological and immunological information demonstrated that the majority were male (74.3%), men who have sex with men (MSM, 68.6%), aged 31–45 years (41.4%) who had been diagnosed with HIV-1 >10 years ago (Table 1); 11.4% of patients had AIDS at the time of HIV-1 diagnosis and 64.3% had AIDS at the time of sampling. For 41.4% of patients, viral load at sampling was 10,000–100,000 copies/mL and 47.1% had a T CD4+ cell count ≤200 cells/mL. These patients had been assigned a previous HIV-1 subtype based on the pol sequence. Twenty patients were infected with subtype B, 24 with CRF19_cpx, 9 with CRF18_cpx and 19 with CRF_BG.

CharacteristicsN (%)Total number of patients*70 (100)Male52 (74.3)Sexual orientation: MSM48 (68.6)Years since HIV-1 diagnosis≤19 (12.9)>1–39 (12.9)>3–53 (4.3)>5–1020 (28.6)>1029 (41.4)Age (years)<152 (2.9)15–3017 (24.3)31–4529 (41.4)46–6120 (28.6)62–772 (2.9)AIDS at diagnosis8 (11.4)Evolution to AIDS after HIV-1 diagnosis (years)<119 (27.1)1–38 (11.4)>318 (25.7)No AIDS25 (35.7)Current AIDS45 (64.3)ARV treated68 (97.1)Viral load (copies/mL)<10003 (4.3)1,000–10,00013 (18.6)>10,000–100,00029 (41.4)>100,00025 (35.7)T CD4+ cell count (cells/mL)<20033 (47.1)200–35020 (28.6)>350–50011 (15.7)>5006 (8.6)

Env subtyping Samples were processed to obtain partial env sequences. Of these, 12.5% (9/72) samples had a C2V3 region too short for phylogenetic subtyping (Figure 1), but were still able to be subtyped using REGA version 3 and COMET version 2 bioinformatics platforms.

After env subtyping, 21 sequences were subtype B, 23 were CRF19_cpx, 7 were CRF18_cpx, and 21 were CRF_BG. Considering both pol and env sequences, an overall subtype was assigned to samples, resulting in 19 classified as subtype B, 21 as CRF19_cpx, 7 as CRF18_cpx, 18 as CRF_BG and 7 as unique recombinant forms (URF) since they had different subtypes in pol and env sequences.

Co-receptor use prediction Co-receptor use was analyzed based on G2P algorithm, 11/25 rule and V3 net charge. Co-receptor prediction by G2P showed that 59.7% (43/72) of analyzed samples were R5-tropic viruses, 16.7% (12/72) were X4-tropic viruses, 23.6% (17/72) were classified as dual tropic viruses and 40.3% (29/72) were predicted to be dual-tropic or X4 variants (Table 1). There was a 42.9% agreement in predicting X4-tropic viruses between G2P and the 11/25 rule for subtype B, while there was 50.0% agreement for non-B subtypes. Calculation of V3-loop net charge showed a 42.9% agreement with G2P for subtype B and 50.0% agreement for non-B subtypes.

Analysis of co-receptor use by G2P showed that CRF19_cpx viruses were more likely to be R5X4-tropic (p = 0.034; OR = 3.295; CI: 1.065–10.191) than other subtypes. This statistical difference was strengthened when we analyzed R5X4 and X4 together using variants (p = 0.015; OR = 3.526; CI: 1.253–9.921). No other statistical association was found for remaining subtypes. A comparison of FPR values across subtypes also showed that CRF19_cpx had the lowest mean value (Figure 2), even though the difference was only significant with regard to CRF18_cpx and CRF_BG.

Patient epidemiological, virological and immunological information was analyzed in relation to co-receptor use and viral env subtype (Table 3). Nevertheless, few statistical associations were found. X4-tropic viruses were associated with >10 years since diagnosis (p = 0.022), with a viral load at >10,000–100,000 copies/mL (p = 0.022) and with AIDS at the time of sampling (p = 0.045). R5-tropic viruses prevailed in non-AIDS patients (p = 0.002), in those with viral loads at 1,000–10,000 copies/mL (p = 0.008) and a T CD4+ cell count at 200–350 cells/mL (p =0 .031). R5X4 viruses corresponded to patients with a viral load >100,000 copies/mL (p = 0.004) and with T CD4+ cell counts <200 cells/mL (p = 0.005). CRF19_cpx prevailed in patients who were 15–30 years old (p = 0.028).

Our analysis incorporated an additional 35 HIV-1 V3-loop sequences from Cuban patients, obtained from January 2014 through January 2015 under the same criteria outlined above. Thus, this set of 107 sequences spanned a period from January 2014 through July 2016. Determination of co-receptor use in this group of viruses (Table 4) revealed an association of CRF19_cpx with X4-tropic viruses (p = 0.044), R5X4-tropic viruses (p = 0.019) or either R5X4 or X4 (p = 0.0004). Additionally, most CRF_BG were R5-tropic viruses (p = 0.032).

Genetic diversity of V3 loop After alignment of 107 V3-loop sequences, we identified some differences in aminoacidic patterns expressed among subtypes. Arginine (R) at position 13 was more prevalent in CRF_BG (p ≤0.0001) and CRF19_cpx (p ≤0.0001) than in subtype B (Supplementary Figure 3).

Valine (V) at positions 12 and 19 prevailed in CRF19_cpx compared to subtype B (p ≤0.0001). On the other hand, for non-B subtypes (CRF19_cpx, CRF18_cpx and CRF_BG), the sequence GPGQ in the tip of the V3 loop was more frequent than GPGR, characteristic of subtype B (p <0.0001). Some motifs found at the tip of the loop other than GPG were RPA/G, APG, GAG, DAG, GRG and GLG. All sequences with APG motifs were R5-tropic, while all CRF19_cpx sequences with substitution 34Y were either X4- or R5X4-tropic strains. However, these differences were not significant.

DISCUSSION

CRF19_cpx is currently one of the most prevalent subtypes circulating in Cuba. This recombinant is thought to be of African origin, resulting in a mosaic of subtypes A1, D and G; though the virus has spread most successfully in Cuba. An earlier study performed in our laboratory showed an association between infection with CRF19_cpx and rapid progression to AIDS.[1515. Kourí V, Khouri R, Alemán Y, Abrahantes Y, Vercauteren J, Pineda-Pena AC, et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine [Internet]. 2015 Jan 28 [cited 2020 Dec 16];2(3):244–54. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4484819/pdf/main.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] These results led us to hypothesize that this recombinant virus may be a more pathogenic form of HIV-1. Although it has not been demonstrated whether X4 viruses are a cause or rather a consequence of immune system deterioration, many studies reported co-receptor use preference based on the subtype.[66. Schuitemaker H, van 't Wout AB, Lusso P. Clinical significance of HIV-1 coreceptor usage. J Transl Med [Internet]. 2011 Jan 27 [cited 2020 Dec 15];9 Suppl 1 (Suppl 1):S5. Available at: https://translational-medicine.biomedcentral.com/articles/10.1186/1479-5876-9-S1-S5
https://translational-medicine.biomedcen...
,1717. Fenyö EM, Esbjörnsson J, Medstrand P, Jansson M. Human immunodeficiency virus type 1 biological variation and coreceptor use: from concept to clinical significance. J Intern Med [Internet]. 2011 Dec [cited 2020 Dec 17];270(6):520–31. Available at: https://onlinelibrary.wiley.com/doi/epdf/10.1111/j.1365-2796.2011.02455.x
https://onlinelibrary.wiley.com/doi/epdf...
] Based on G2P prediction, we found that CRF19_cpx had preferential tropism towards R5X4, which has been described as capable of using both CCR5 and CXCR4 co-receptors for entry into the host cell.[44. Panos G, Watson DC. Effect of HIV-1 subtype and tropism on treatment with chemokine coreceptor entry inhibitors; overview of viral entry inhibition. Crit Rev Microbiol [Internet]. 2015 [cited 2020 Dec 2];41(4):473–87. Available at: https://www.tandfonline.com/doi/full/10.3109/1040841X.2013.867829?scroll=top&needAccess=true
https://www.tandfonline.com/doi/full/10....
]

When we expanded our analysis to 107 samples, we found CRF19_cpx to be associated with dual-tropic forms of the virus and also ratified our previous observations about a link between these viruses and X4-tropic variants.

Since dual-tropic viruses have been considered an intermediate stage during the switch from R5 to X4 variants that may even result from intrapatient recombination once X4 viruses emerge in the host, CRF19_cpx association is expected not only with X4 but also with R5X4 viruses. In addition, understanding that some patients never experience this switch in tropism and continue harboring R5 viruses,[77. Mild M, Esbjornsson J, Fenyo EM, Medstrand P. Frequent intrapatient recombination between human immunodeficiency virus type 1 R5 and X4 envelopes: implications for coreceptor switch. J Virol [Internet]. 2007 Apr [cited 2020 Dec 16];81(7):3369–76. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1866041/pdf/1295-06.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,2727. Fogel GB, Lamers SL, Liu ES, Salemi M, McGrath MS. Identification of dual-tropic HIV-1 using evolved neural networks. Biosystems [Internet]. 2015 Nov [cited 2020 Dec 3];137:12–9. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4921197/pdf/nihms730074.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] we speculate that some HIV-1 subtypes or CRFs, such as CRF19_cpx, are likely to evolve towards dual-tropic and X4 variants during the course of infection. Whether this evolution is faster with respect to other subtypes or is actually related to HIV-1 pathogenesis has yet to be demonstrated. Furthermore, the cause of this subtype-dependent co-receptor use may be explained by intrinsic virological properties of each subtype (i.e. enhanced replicative capacity, antigenic diversity) or by the differential effects of viruses on host immune systems, which may eventually lead to expansion of the X4 virus population.[2727. Fogel GB, Lamers SL, Liu ES, Salemi M, McGrath MS. Identification of dual-tropic HIV-1 using evolved neural networks. Biosystems [Internet]. 2015 Nov [cited 2020 Dec 3];137:12–9. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4921197/pdf/nihms730074.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
]

Even though we measured no direct association between CRF19_cpx and disease progression stage, R5X4 viruses were related to a CD4+ cell count ≤200 cells/mL. Characteristically, X4 viruses have been correlated with advanced disease stage and low levels of CD4+ T cell count.[2828. Matume ND, Tebit DM, Gray LR, Hammarskjold ML, Rekosh D, Bessong PO. Next generation sequencing reveals a high frequency of CXCR4 utilizing viruses in HIV-1 chronically infected drug experienced individuals in South Africa. J Clin Virol [Internet]. 2018 Jun [cited 2020 Dec 17];103:81–7. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC7229640/pdf/nihms-1584881.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
3030. Gupta S, Neogi U, Srinivasa H, Banerjea AC, Shet A. HIV-1 coreceptor tropism in India: increasing proportion of X4-tropism in subtype C strains over two decades. J Acquir Immune Defic Syndr [Internet]. 2014 Apr 1 [cited 2020 Dec 17];65(4):397–404. Available at: https://journals.lww.com/jAIDS/Fulltext/2014/04010/HIV_1_Coreceptor_Tropism_in_India Increasing.3.aspx
https://journals.lww.com/jAIDS/Fulltext/...
] Nevertheless, it has been hypothesized that dual-tropic viruses could represent variants with the same potential as R5 viruses to evade immune response and also infect cells expressing CXCR4. The broadening in cell tropism of the viral population to include CXCR4-expressing cells would result in increased CD4+ T cell death and further immune impairment.[77. Mild M, Esbjornsson J, Fenyo EM, Medstrand P. Frequent intrapatient recombination between human immunodeficiency virus type 1 R5 and X4 envelopes: implications for coreceptor switch. J Virol [Internet]. 2007 Apr [cited 2020 Dec 16];81(7):3369–76. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1866041/pdf/1295-06.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
]

figure 1
Figure 1: Phylogenetic tree constructed by PhyML/One Click software[2222. Lemoine F, Correia D, Lefort V, Doppelt-Azeroual O, Mareuil F, Cohen-Boulakia S, et al. NGPhylogeny.fr: new generation phylogenetic services for non-specialists. Nucleic Acids Res [Internet]. 2019 Apr 27 [cited 2020 Dec 3];47(W1):W260–W5. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6602494/pdf/gkz303.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
]

Sixty-three sequences encompassing C2V3 env region of HIV-1 (in bold) were compared with reference sequences from Genbank. Bootstrap values (as a result of 500 replicates), above 50% are shown at the nodes. The tree was rooted in strain K.CD.97. EQTB11C.AJ249235.

Sample number (%)Subtype in env regionAll subtypesSubtype BCRF19 _cpxCRF18 _cpxCRF_BG (20-23-24)72 (100.0)21 (100.0)23 (100.0)7 (100.0)21 (100.0)R5 Samples number (%)43 (59.7)13 (61.9)9 (39.1)6 (85.7)15 (71.4)R5X4 Sample number (%)17 (23.6)2 (9.5)9 (39.1)1 (14.3)5 (23.8)p-value R5X4p = 0.034 OR = 3.295 CI: 1.065–0.191X4 Sample number (%)12 (16.7)6 (28.6)5 (21.7)0 (0.0)1 (4.8)R5X4/X4 Sample number (%)29 (40.3)8 (38.1)14 (60.9)1 (14.3)6 (28.6)p-value R5X4/X4*p = 0.015 OR = 3.526 CI: 1.253–9.921
Figure 2
Figure 2: Alignment of 107 V3-loop amino acid sequences obtained from Cuban HIV-1 patients.

FPR: false positive rate

Subtype B identity is indicated by dots.

CharacteristicsN (%)R5R5X4X4Subtype B envCRF18_cpx envCRF19_cpx envCRF_BG env70* (100.0)42 (100.0)17 (100.0)11 (100.0)21 (100.0)7 (100.0)22 (100.0)20 (100.0)Male sex52 (74.3)29 (69.0)13 (76.5)10 (90.9)15 (71.4)6 (85.7)14 (63.6)17 (85.0)Sexual orientation: MSM48 (68.6)25 (59.5)13 (76.5)10 (90.9)14 (66.7)6 (85.7)13 (59.1)15 (75.0)Years since HIV-1 diagnosis≤19 (12.9)7 (16.7)2 (11.8)0 (0.0)1 (4.8)0 (0.0)4 (18.2)4 (20.0)>1–39 (12.9)6 (14.3)2 (11.8)1 (9.1)2 (9.5)0 (0.0)4 (18.2)3 (15.0)>3–53 (4.3)3 (7.1)0 (0.0)0 (0.0)1 (4.8)0 (0.0)0 (0.0)2 (10.0)>5–1020 (28.6)12 (28.6)6 (35.3)2 (18.2)8 (38.1)2 (28.6)5 (22.7)5 (20.0)>1029 (41.4)14 (33.3)7 (41.2)8 (72.7)9 (42.9)5 (71.4)9 (40.9)6 (30.0)Age (years)<152 (2.9)2 (28.6)0 (0.0)0 (0.0)0 (0.0)1 (14.3)0 (0.0)1 (5.0)15–3017 (24.3)13 (31.0)4 (23.5)0 (0.0)3 (14.3)0 (0.0)9 (40.9)5 (25.0)31–4529 (41.4)18 (42.9)7 (41.2)4 (36.4)9 (42.9)4 (57.1)6 (27.3)10 (50.0)46–6120 (28.6)9 (21.4)6 (35.3)5 (45.5)9 (42.9)2 (28.6)5 (22.7)4 (20.0)62–772 (2.9)0 (0.0)0 (0.0)2 (18.2)0 (0.0)0 (0.0)2 (9.1)0 (0.0)AIDS at diagnosis8 (11.4)4 (9.5)2 (11.8)2 (18.2)4 (19.0)0 (0.0)2 (9.1)2 (10.0)ARV treated68 (97.1)40 (95.2)17 (100.0)11 (100.0)21 (100.0)7 (100.0)21 (95.5)19 (95.0)Evolution to AIDS after HIV-1 diagnosis (years)<119 (27.1)10 (23.8)5 (29.4)4 (36.4)7 (33.3)1 (14.3)6 (27.3)5 (25.0)1–38 (11.4)3 (7.1)3 (17.6)2 (18.2)2 (9.5)1 (14.3)1 (4.6)4 (20.0)>318 (25.7)8 (19.0)6 (35.3)4 (36.4)5 (23.8)3 (42.9)6 (27.3)4 (20.0)No AIDS25 (35.7)21 (50.0)3 (17.6)1 (9.1)7 (33.3)2 (28.6)9 (40.9)7 (35.0)Actual AIDS45 (64.3)21 (50.0)14 (82.4)10 (90.9)14 (66.7)5 (71.4)13 (59.1)13 (65.0)Viral load (copies/mL)<10003 (4.3)1 (2.4)1 (5.9)1 (9.1)1 (4.8)0 (0.0)1 (4.6)1 (5.0)1000–10,00013 (18.6)12 (28.6)0 (0.0)1 (9.1)3 (14.3)1 (14.3)3 (13.6)6 (30.0)>10,000–100,00029 (41.4)16 (38.1)5 (29.4)8 (72.7)10 (47.6)2 (28.6)11 (50)6 (30.0)>100,00025 (35.7)13 (31.0)11 (64.7)1 (9.1)7 (33.3)4 (57.1)7 (31.8)7 (35.0)T CD4+ cell count/mL<20033 (47.1)14 (33.3)13 (76.5)6 (54.5)8 (38.1)4 (57.1)10 (45.5)11 (55.0)200–35020 (28.6)16 (38.1)0 (0.0)4 (36.4)6 (28.6)2 (28.6)8 (36.4)4 (20.0)>350–50011 (15.7)8 (19.0)2 (11.8)1 (9.1)6 (28.6)1 (14.3)2 (9.1)2 (10.0)>5006 (8.6)4 (9.5)2 (11.8)0 (0.0)1 (4.8)0 (0.0)2 (9.1)3 (15.0)
Subtype in env regionSample number (%)All SubtypesSubtype BCRF19 _cpxCRF18_cpxCRF_BG (20-23-24)107 (100)35 (100.0)29 (100.0)9 (100.0)34 (100.0)R5 Sample number (%)66 (61.7)22 (62.9)10 (34.5)8 (88.9)26 (76.5)p-value R5p = 0.032 OR: 2.681 CI: 1.072–6.706R5X4 Sample number (%)24 (22.4)6 (17.1)11 (37.9)1 (11.1)6 (17.6)p-value R5X4p = 0.019 OR: 3.056 CI:1.173–7.962X4 Sample number (%)17 (15.9)7 (20.0)8 (27.6)0 (0.0)2 (5.9)p-value X4p = 0.044 OR: 2.921 CI:1.001–8.518R5X4/X4 Sample number (%)41 (38.3)13 (37.1)19 (65.5)1 (11.1)8 (23.5)p-value R5X4/X4p = 0.0004 OR:4.836 CI:1.945–12.025

To our knowledge, this is the first report on the association of CRF_BG viruses and R5 tropism. Previously, a link between CRF18_cpx viruses and R5 tropism had been reported.[1515. Kourí V, Khouri R, Alemán Y, Abrahantes Y, Vercauteren J, Pineda-Pena AC, et al. CRF19_cpx is an evolutionary fit HIV-1 variant strongly associated with rapid progression to AIDS in Cuba. EBioMedicine [Internet]. 2015 Jan 28 [cited 2020 Dec 16];2(3):244–54. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4484819/pdf/main.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] We found no significant relationship between those variables in the current study.

We also observed some changes in the aminoacidic sequence of the V3 loop among different subtypes. Particularly, some motifs other than GPG were found at the tip of the loop. In accordance with previous observations, GPGR motifs prevailed in subtype B, while the GPGQ motif was predominant among Cuban recombinants.[3131. Guo JL, Yan Y, Zhang JF, Ji JM, Ge ZJ, Ge R, et al. Co-receptor tropism and genetic characteristics of the V3 regions in variants of antiretroviral-naive HIV-1 infected subjects [Internet]. Epidemiol Infect. 2019 Jan [cited 2020 Dec 18];147:e181. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6518647/pdf/S0950268819000700a.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,3232. Morgado MG, Guimaraes ML, Gripp CB, Neves Junior I, Costa CI, dos Santos VG, et al. Polymorphism of the predictive antigenic sites on the V3 loop of Brazilian HIV-1 subtype B strains. HEC/FIOCRUZ AIDS Clinical Research Group. Mem Inst Oswaldo Cruz [Internet]. 1996 May–Jun [cited 2020 Dec 18];91(3):339–42. Available at: https://www.scielo.br/pdf/mioc/v91n3/15.pdf
https://www.scielo.br/pdf/mioc/v91n3/15....
] The GPGR motif is a very important target of neutralizing antibodies in subtype B viruses, so its difference with other subtypes should be considered for assessing monoclonal antibodies and peptide vaccine design.

In 1997, Milich et al. reported a high frequency of valine (V) at position 19 in syncytium-inducing viruses. According to these authors, the interaction between 19V and phenylalanine (F) at position 20 would contribute to a switch in tropism from non-syncytium–inducing to syncytium-inducing viruses.[3333. Milich L, Margolin BH, Swanstrom R. Patterns of amino acid variability in NSI-like and SI-like V3 sequences and a linked change in the CD4-binding domain of the HIV-1 Env protein. Virology [Internet]. 1997 Dec 8 [cited 2020 Dec 18];239(1):108–18. Available at: https://linkinghub.elsevier.com/retrieve/pii/S0042-6822(97)98821-8
https://linkinghub.elsevier.com/retrieve...
] This hypothesis could also explain the association between CRF19_cpx viruses and R5X4 and X4 tropism found in the current study. Considering that almost all sequences analyzed had V19, we could ask ourselves if this recombinant is capable of evolving to a change in tropism faster than other subtypes.

This study has some limitations. The number of sequences we have obtained so far is still small for ascertaining an association between co-receptor tropism and subtype; hence results need to be confirmed by future research including a greater number of patients. Additionally, currently-employed tools for HIV-1 tropism prediction, including G2P, were developed on the basis of a restricted set of subtypes, mainly subtypes B and C.[11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
] Therefore, their accuracy in predicting tropism for other subtypes and recombinant forms including the ones circulating in Cuba have to be assessed in the future based on their correlation with phenotypic methods. A phenotypic test was not performed due to restricted availability.

In Cuba, maraviroc could be a treatment alternative for patients experiencing failure with other HIV-1 drugs or be considered an option in combination with other regimens; however determining HIV-1 co-receptor use is mandatory before its prescription.[11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,88. Phuphuakrat A, Phawattanakul S, Pasomsub E, Kiertiburanakul S, Chantratita W, Sungkanuparph S. Coreceptor tropism determined by genotypic assay in HIV-1 circulating in Thailand, where CRF01_AE predominates. HIV Med [Internet]. 2014 May [cited 2020 Dec 15];15(5):269–75. Available at: https://onlinelibrary.wiley.com/doi/epdf/10.1111/hiv.12108
https://onlinelibrary.wiley.com/doi/epdf...
,3434. Heera J, Valluri SR, Craig C, Fang A, Thomas N, Meyer RD, et al. First prospective comparison of genotypic versus phenotypic tropism assays in predicting virologic responses to maraviroc in a phase 3 study. New Microbiol [Internet]. 2019 Apr [cited 2020 Dec 19];41(2):101–7. http://www.newmicrobiologica.org/PUB/allegati_pdf/2019/2/101.pdf
http://www.newmicrobiologica.org/PUB/all...
] A more thorough understanding of the possible associations between co-receptor tropism and subtypes circulating in Cuba, together with phenotypic assays and the possible creation of new bioinformatic algorithms more suitable for Cuban HIV-1 recombinants’ tropism prediction would facilitate proper decisions regarding treatment choice with maraviroc or other co-receptor inhibitors.

CONCLUSION

In the current study, 72 subtypes obtained from Cuban HIV-1 infected patients were sequenced and tropism was predicted employing a G2P bioinformatic tool.

Additionally, 35 other viruses that were sequenced in a former study were included for analysis. An association between CRF19_cpx and X4 and R5X4 tropism was found, which suggests that patients infected with this recombinant are probably less suitable to receive maraviroc than those infected with other HIV-1 subtypes. The study reinforces findings from previous studies on this subject and highlights the importance of HIV-1 diversity when considering pathogenesis and treatment options.[11. Riemenschneider M, Cashin KY, Budeus B, Sierra S, Shirvani-Dastgerdi E, Bayanolhagh S, et al. Genotypic Prediction of Co-receptor Tropism of HIV-1 Subtypes A and C. Sci Rep [Internet]. 2016 [cited 2020 Dec 3];6:24883. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4850382/pdf/srep24883.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,3535. Esbjörnsson J, Månsson F, Martínez-Arias W, Vincic E, Biague AJ, da Silva ZJ, et al. Frequent CXCR4 tropism of HIV-1 subtype A and CRF02_AG during late-stage disease--indication of an evolving epidemic in West Africa. Retrovirology [Internet]. 2010 Mar 22 [cited 2020 Dec 18];7:23. DOI: 10.1186/1742-4690-7-23. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC2855529/pdf/1742-4690-7-23.pdf
https://www.ncbi.nlm.nih.gov/pmc/article...
,3636. Church JD, Huang W, Mwatha A, Toma J, Stawiski E, Donnell D, et al. HIV-1 tropism and survival in vertically infected Ugandan infants. J Infect Dis [Internet]. 2008 May 15 [cited 2020 Dec 18];197(10):1382–8. Available at: https://pdfs.semanticscholar.org/be30/beae653aaaea6cd4305dd42ae01d7232364d.pdf?_ga=2.8524978.60519488.1608083556-726339794.1548872336
https://pdfs.semanticscholar.org/be30/be...
]

REFERENCES

Publication Dates

  • Publication in this collection
    04 Nov 2022
  • Date of issue
    Jul-Dec 2021

History

  • Received
    13 Jan 2021
  • Accepted
    21 July 2021
Medical Education Cooperation with Cuba Oakland - California - United States
E-mail: editors@medicc.org