Unravelling the complicated evolutionary and dissemination history of HIV-1M subtype A lineages

Female Sex Hormones Upregulate the Replication Activity of HIV-1 Sub-Subtype A6 and CRF02_AG but Not HIV-1 Subtype B

More than 50% of all people living with HIV worldwide are women. Globally, HIV/AIDS is the leading cause of death among women aged 15 to 44. The safe and effective methods of hormonal contraception are an essential component of preventive medical care in order to reduce maternal and infant mortality. However, there is limited knowledge regarding the effect of hormones on the rate of viral replication in HIV infection, especially non-B subtypes. The goal of the present work was to study in vitro how the female hormones β-estradiol and progesterone affect the replication of the HIV-1 subtypes A6, CRF02_AG, and B. The findings show that high doses of hormones enhanced the replication of HIV-1 sub-subtype A6 by an average of 1.75 times and the recombinant variant CRF02_AG by 1.4 times but did not affect the replication of HIV-1 subtype B. No difference was detected in the expression of CCR5 and CXCR4 co-receptors on the cell surface, either in the presence or absence of hormones. However, one of the reasons for the increased viral replication could be the modulated TLRs secretion, as it was found that high doses of estradiol and progesterone upregulated, to varying degrees, the expression of TLR2 and TLR9 genes in the PBMCs of female donors infected with HIV-1 sub-subtype A6.

Origin and evolution of HIV-1 subtype A6

Background

HIV outbreaks in the Former Soviet Union (FSU) countries were characterized by repeated transmission of the HIV variant AFSU, which is now classified as a distinct subtype A sub-subtype called A6. The current study used phylogenetic/phylodynamic and signature mutation analyses to determine likely evolutionary relationship between subtype A6 and other subtype A sub-subtypes.

Methods

For this study, an initial Maximum Likelihood phylogenetic analysis was performed using a total of 553 full-length, publicly available, reverse transcriptase sequences, from A1, A2, A3, A4, A5, and A6 sub-subtypes of subtype A. For phylogenetic clustering and signature mutation analysis, a total of 5961 and 3959 pol and env sequences, respectively, were used.

Results

Phylogenetic and signature mutation analysis showed that HIV-1 sub-subtype A6 likely originated from sub-subtype A1 of African origin. A6 and A1 pol and env genes shared several signature mutations that indicate genetic similarity between the two subtypes. For A6, tMRCA dated to 1975, 15 years later than that of A1.

Conclusion

The current study provides insights into the evolution and diversification of A6 in the backdrop of FSU countries and indicates that A6 in FSU countries evolved from A1 of African origin and is getting bridged outside the FSU region.

Genetic Characterization of a New HIV-1 Sub-Subtype A in Cabo Verde, Denominated A8

Previous molecular characterization of Human immunodeficiency virus (HIV-1) samples from Cabo Verde pointed out a vast HIV-1 pol diversity, with several subtypes and recombinant forms, being 5.2% classified as AU-pol. Thus, the aim of the present study was to improve the characterization of these AU sequences. The genomic DNA of seven HIV-1 AU pol-infected individuals were submitted to four overlapping nested-PCR fragments aiming to compose the full-length HIV-1 genome. The final classification was based on phylogenetic trees that were generated using the maximum likelihood and bootscan analysis. The genetic distances were calculated using Mega 7.0 software. Complete genome amplification was possible for two samples, and partial genomes were obtained for the other five. These two samples grouped together with a high support value, in a separate branch from the other sub-subtypes A and CRF26_A5U. No recombination was verified at bootscan, leading to the classification of a new sub-subtype A. The intragroup genetic distance from the new sub-subtype A at a complete genome was 5.2%, and the intergroup genetic varied from 8.1% to 19.0% in the analyzed fragments. Our study describes a new HIV-1 sub-subtype A and highlights the importance of continued molecular surveillance studies, mainly in countries with high HIV molecular diversity.

Elucidation of Early Evolution of HIV-1 Group M in the Congo Basin Using Computational Methods

The Congo Basin region is believed to be the site of the cross-species transmission event that yielded HIV-1 group M (HIV-1M). It is thus likely that the virus has been present and evolving in the region since that cross-species transmission. As HIV-1M was only discovered in the early 1980s, our directly observed record of the epidemic is largely limited to the past four decades. Nevertheless, by exploiting the genetic relatedness of contemporary HIV-1M sequences, phylogenetic methods provide a powerful framework for investigating simultaneously the evolutionary and epidemiologic history of the virus. Such an approach has been taken to find that the currently classified HIV-1 M subtypes and Circulating Recombinant Forms (CRFs) do not give a complete view of HIV-1 diversity. In addition, the currently identified major HIV-1M subtypes were likely genetically predisposed to becoming a major component of the present epidemic, even before the events that resulted in the global epidemic. Further efforts have identified statistically significant hot- and cold-spots of HIV-1M subtypes sequence inheritance in genomic regions of recombinant forms. In this review we provide ours and others recent findings on the emergence and spread of HIV-1M variants in the region, which have provided insights into the early evolution of this virus.

Genetic variability of the U5 and downstream sequence of major HIV-1 subtypes and circulating recombinant forms

The critical role of the regulatory elements at the 5′ end of the HIV-1 genome in controlling the life cycle of HIV-1 indicates that this region significantly influences virus fitness and its biological properties. In this study, we performed a detailed characterization of strain-specific variability of sequences from the U5 to upstream of the gag gene start codon of diverse HIV-1 strains by using next-generation sequencing (NGS) techniques. Overall, we found that this region of the HIV-1 genome displayed a low degree of intra-strain variability. On the other hand, inter-strain variability was found to be as high as that reported for gag and env genes (13–17%). We observed strain-specific single point and clustered mutations in the U5, PBS, and gag leader sequences (GLS), generating potential strain-specific transcription factor binding sites (TFBS). Using an infrared gel shift assay, we demonstrated the presence of potential TFBS such as E-box in CRF22_01A, and Stat 6 in subtypes A and G, as well as in their related CRFs. The strain-specific variation found in the sequence corresponding at the RNA level to functional domains of the 5ʹ UTR, could also potentially impact the secondary/tertiary structural rearrangement of this region. Thus, the variability observed in this 5′ end of the genomic region of divergent HIV-1 strains strongly suggests that functions of this region might be affected in a strain-specific manner. Our findings provide new insights into DNA–protein interactions that regulate HIV-1 replication and the influence of strain characterization on the biology of HIV-1 infection.

HIV-1 Transmission Patterns Within and Between Risk Groups in Coastal Kenya

HIV-1 transmission patterns within and between populations at different risk of HIV-1 acquisition in Kenya are not well understood. We investigated HIV-1 transmission networks in men who have sex with men (MSM), injecting drug users (IDU), female sex workers (FSW) and heterosexuals (HET) in coastal Kenya. We used maximum-likelihood and Bayesian phylogenetics to analyse new (N = 163) and previously published (N = 495) HIV-1 polymerase sequences collected during 2005-2019. Of the 658 sequences, 131 (20%) were from MSM, 58 (9%) IDU, 109 (17%) FSW, and 360 (55%) HET. Overall, 206 (31%) sequences formed 61 clusters. Most clusters (85%) consisted of sequences from the same risk group, suggesting frequent within-group transmission. The remaining clusters were mixed between HET/MSM (7%), HET/FSW (5%), and MSM/FSW (3%) sequences. One large IDU-exclusive cluster was found, indicating an independent sub-epidemic among this group. Phylodynamic analysis of this cluster revealed a steady increase in HIV-1 infections among IDU since the estimated origin of the cluster in 1987. Our results suggest mixing between high-risk groups and heterosexual populations and could be relevant for the development of targeted HIV-1 prevention programmes in coastal Kenya.

Phylogenetic measures of indel rate variation among the HIV-1 group M subtypes

The transmission fitness and pathogenesis of HIV-1 is disproportionately influenced by evolution in the five variable regions (V1–V5) of the surface envelope glycoprotein (gp120). Insertions and deletions (indels) are a significant source of evolutionary change in these regions. However, the rate and composition of indels has not yet been quantified through a large-scale comparative analysis of HIV-1 sequences. Here, we develop and report results from a phylogenetic method to estimate indel rates for the gp120 variable regions across five major subtypes and two circulating recombinant forms (CRFs) of HIV-1 group M. We processed over 26,000 published HIV-1 gp120 sequences, from which we extracted 6,605 sequences for phylogenetic analysis. We reconstructed time-scaled phylogenies by maximum likelihood and fit a binomial-Poisson model to the observed distribution of indels between closely related pairs of sequences in each tree (cherries). By focusing on cherries in each tree, we obtained phylogenetically independent indel reconstructions, and the shorter time scales in cherries reduced the bias due to purifying selection. Rate estimates ranged from 3.0×10−5 to 1.5×10−3 indels/nt/year and varied significantly among variable regions and subtypes. Indel rates were significantly lower in V3 relative to V1, and were also lower in HIV-1 subtype B relative to the 01_AE reference. We also found that V1, V2, and V4 tended to accumulate significantly longer indels. Furthermore, we observed that the nucleotide composition of indels was distinct from the flanking sequence, with higher frequencies of G and lower frequencies of T. Indels affected N-linked glycosylation sites more often in V1 and V2 than expected by chance, consistent with positive selection on glycosylation patterns within these regions. These results represent the first comprehensive measures of indel rates in HIV-1 gp120 across multiple subtypes and CRFs, and identifies novel and unexpected patterns for further research in the molecular evolution of HIV-1.

PANGEA-HIV 2: Phylogenetics And Networks for Generalised Epidemics in Africa

PURPOSE OF REVIEW

The HIV epidemic in sub-Saharan Africa is far from being under control and the ambitious UNAIDS targets are unlikely to be met by 2020 as declines in per-capita incidence being largely offset by demographic trends. There is an increasing number of proven and specific HIV prevention tools, but little consensus on how best to deploy them.

RECENT FINDINGS

Traditionally, phylogenetics has been used in HIV research to reconstruct the history of the epidemic and date zoonotic infections, whereas more recent publications focus on HIV diversity and drug resistance. However, it is also the most powerful method of source attribution available for the study of HIV transmission. The PANGEA (Phylogenetics And Networks for Generalized Epidemics in Africa) consortium has generated over 18 000 NGS HIV sequences from five countries in sub-Saharan Africa. Using phylogenetic methods, we will identify characteristics of individuals or groups, which are most likely to be at risk of infection or at risk of infecting others.

SUMMARY

Combining phylogenetics, phylodynamics and epidemiology will allow PANGEA to highlight where prevention efforts should be focussed to reduce the HIV epidemic most effectively. To maximise the public health benefit of the data, PANGEA offers accreditation to external researchers, allowing them to access the data and join the consortium. We also welcome submissions of other HIV sequences from sub-Saharan Africa to the database.

Transmission patterns of HIV-1 non-R5 strains in Poland

HIV-1 env sequencing enables predictions of viral coreceptor tropism and phylogenetic investigations of transmission events. The aim of the study was to estimate the contribution of non-R5 strains to the viral spread in Poland. Partial proviral env sequences were retrieved from baseline blood samples of patients with newly diagnosed HIV-1 infection between 2008-2014, including 46 patients with recent HIV-1 infection (RHI), and 246 individuals with long-term infection (LTHI). These sequences were subjected to the genotypic coreceptor tropism predictions and phylogenetic analyses to identify transmission clusters. Overall, 27 clusters with 57 sequences (19.5%) were detected, including 15 sequences (26.3%) from patients with RHI. The proportion of non-R5 strains among all study participants was 23.3% (68/292), and was comparable between patients with RHI and LTHI (11/46, 23.9% vs 57/246, 23.2%; p = 1.000). All 11 patients with non-R5 strains and RHI were men having sex with men (MSM). Among these patients, 4 had viral sequences grouped within phylogenetic cluster with another sequence of non-R5 strain obtained from patient with LTHI, indicating potential acquisition of non-R5 HIV-1 for at least 4/46 (8.7%) patients with RHI. We were unable to confirm the contribution of patients with RHI to the forward transmission of non-R5 strains, but a relatively high proportion of non-R5 strains among them deserves attention due to the limited susceptibility to CCR5 antagonists.

Geographically-stratified HIV-1 group M pol subtype and circulating recombinant form sequences

Accurate classification of HIV-1 group M lineages, henceforth referred to as subtyping, is essential for understanding global HIV-1 molecular epidemiology. Because most HIV-1 sequencing is done for genotypic resistance testing pol gene, we sought to develop a set of geographically-stratified pol sequences that represent HIV-1 group M sequence diversity. Representative pol sequences differ from representative complete genome sequences because not all CRFs have pol recombination points and because complete genome sequences may not faithfully reflect HIV-1 pol diversity. We developed a software pipeline that compiled 6,034 one-per-person complete HIV-1 pol sequences annotated by country and year belonging to 11 pure subtypes and 70 CRFs and selected a set of sequences whose average distance to the remaining sequences is minimized for each subtype/CRF and country to generate a Geographically-Stratified set of 716 Pol Subtype/CRF (GSPS) reference sequences. We provide extensive data on pol diversity within each subtype/CRF and country combination. The GSPS reference set will also be useful for HIV-1 pol subtyping.

Patterns of genomic site inheritance in HIV-1M inter-subtype recombinants delineate the most likely genomic sites of subtype-specific adaptation

Recombination between different HIV-1 group M (HIV-1M) subtypes is a major contributor to the ongoing genetic diversification of HIV-1M. However, it remains unclear whether the different genome regions of recombinants are randomly inherited from the different subtypes. To elucidate this, we analysed the distribution within 82 circulating and 201 unique recombinant forms (CRFs/URFs), of genome fragments derived from HIV-1M Subtypes A, B, C, D, F, and G and CRF01_AE. We found that viruses belonging to the analysed HIV-1M subtypes and CRF01_AE contributed certain genome fragments more frequently during recombination than other fragments. Furthermore, we identified statistically significant hot-spots of Subtype A sequence inheritance in genomic regions encoding portions of Gag and Nef, Subtype B in Pol, Tat and Env, Subtype C in Vif, Subtype D in Pol and Env, Subtype F in Gag, Subtype G in Vpu-Env and Nef, and CRF01_AE inheritance in Vpu and Env. The apparent non-randomness in the frequencies with which different subtypes have contributed specific genome regions to known HIV-1M recombinants is consistent with selection strongly impacting the survival of inter-subtype recombinants. We propose that hotspots of genomic region inheritance are likely to demarcate the locations of subtype-specific adaptive genetic variations.