Incidence of non-B subtypes of HIV-1 in Galicia, Spain: high frequency and diversity of HIV-1 among men who have sex with men

Transmission Clusters, Predominantly Associated With Men Who Have Sex With Men, Play a Main Role in the Propagation of HIV-1 in Northern Spain (2013–2018)

Viruses of HIV-1-infected individuals whose transmission is related group phylogenetically in transmission clusters (TCs). The study of the phylogenetic relations of these viruses and the factors associated with these individuals is essential to analyze the HIV-1 epidemic. In this study, we examine the role of TCs in the epidemiology of HIV-1 infection in Galicia and the Basque County, two regions of northern Spain. A total of 1,158 HIV-1-infected patients from both regions with new diagnoses (NDs) in 2013–2018 were included in the study. Partial HIV-1 pol sequences were analyzed phylogenetically by approximately maximum-likelihood with FastTree 2. In this analysis, 10,687 additional sequences from samples from HIV-1-infected individuals collected in Spain in 1999–2019 were also included to assign TC membership and to determine TCs’ sizes. TCs were defined as those which included viruses from ≥4 individuals, at least 50% of them Spaniards, and with ≥0.95 Shimodaira-Hasegawa-like node support in the phylogenetic tree. Factors associated to TCs were evaluated using odds ratios (OR) and their 95% CI. Fifty-one percent of NDs grouped in 162 TCs. Male patients (OR: 2.6; 95% CI: 1.5–4.7) and men having sex with men (MSM; OR: 2.1; 95% CI: 1.4–3.2) had higher odds of belonging to a TC compared to female and heterosexual patients, respectively. Individuals from Latin America (OR: 0.3; 95% CI: 0.2–0.4), North Africa (OR: 0.4; 95% CI: 0.2–1.0), and especially Sub-Saharan Africa (OR: 0.02; 95% CI: 0.003–0.2) were inversely associated to belonging to TCs compared to native Spaniards. Our results show that TCs are important components of the HIV-1 epidemics in the two Spanish regions studied, where transmission between MSM is predominant. The majority of migrants were infected with viruses not belonging to TCs that expand in Spain. Molecular epidemiology is essential to identify local peculiarities of HIV-1 propagation. The early detection of TCs and prevention of their expansion, implementing effective control measures, could reduce HIV-1 infections.

Cross-Continental Dispersal of Major HIV-1 CRF01_AE Clusters in China

Since the 1990s, several distinct clusters of human immunodeficiency virus-type 1 (HIV-1) CRF01_AE related to a large epidemic in China have been identified, but it is yet poorly understood whether its transmission has dispersed globally. We aimed to characterize and quantify the genetic relationship of HIV-1 CRF01_AEs circulating in China and other countries. Using representative sequences of Chinese clusters as queries, all relevant CRF01_AE pol sequences in two large databases (the Los Alamos HIV sequence database and the UK HIV Drug Resistance Database) were selected with the online basic local alignment search (BLAST) tool. Phylogenetic and phylogeographic analyses were then carried out to characterize possible linkage of CRF01_AE strains between China and the rest of the world. We identified that 269 strains isolated in other parts of the world were associated with five major Chinese CRF01_AE clusters. 80.7% were located within CN.01AE.HST/IDU-2, most of which were born in Southeast Asia. 17.8% were clustered with CN.01AE.MSM-4 and -5. Two distinct sub-clusters associated with Chinese men who have sex with men (MSM) emerged in HK-United Kingdom and Japan after 2000. Our analysis suggests that HIV-1 CRF01_AE strains related to viral transmission in China were initially brought to the United Kingdom or other countries during the 1990s by Asian immigrants or returning international tourists from Southeast Asia, and then after having circulated among MSM in China for several years, these Chinese strains dispersed outside again, possibly through MSM network. This study provided evidence of regional and global dispersal of Chinese CRF01_AE strains. It would also help understand the global landscape of HIV epidemic associated with CRF01_AE transmission and highlight the need for further international collaborative study in this field.

Expanding HIV-1 subtype B transmission networks among men who have sex with men in Poland

INTRODUCTION

Reconstruction of HIV transmission links allows to trace the spread and dynamics of infection and guide epidemiological interventions. The aim of this study was to characterize transmission networks among subtype B infected patients from Poland.

MATERIAL AND METHODS

Maximum likelihood phylogenenetic trees were inferred from 966 HIV-1 subtype B protease/reverse transcriptase sequences from patients followed up in nine Polish HIV centers. Monophyletic clusters were identified using 3% within-cluster distance and 0.9 bootstrap values. Interregional links for the clusters were investigated and time from infection to onward transmission estimated using Bayesian dated MCMC phylogeny.

RESULTS

Three hundred twenty one (33.2%) sequences formed 109 clusters, including ten clusters of ≥5 sequences (n = 81, 8.4%). Transmission networks were more common among MSM (234 sequences, 68.6%) compared to other infection routes (injection drug use: 28 (8.2%) and heterosexual transmissions: 59 (17.3%) cases, respectively [OR:3.5 (95%CI:2.6-4.6),p<0.001]. Frequency of clustering increased from 26.92% in 2009 to 50.6% in 2014 [OR:1.18 (95%CI:1.06-1.31),p = 0.0026; slope +2.8%/year] with median time to onward transmission within clusters of 1.38 (IQR:0.59-2.52) years. In multivariate models clustering was associated with both MSM transmission route [OR:2.24 (95%CI:1.38-3.65),p<0.001] and asymptomatic stage of HIV infection [OR:1.93 (95%CI:1.4-2.64),p<0.0001]. Additionally, interregional networks were linked to MSM transmissions [OR:4.7 (95%CI:2.55-8.96),p<0.001].

CONCLUSIONS

Reconstruction of the HIV-1 subtype B transmission patterns reveals increasing degree of clustering and existence of interregional networks among Polish MSM. Dated phylogeny confirms the association between onward transmission and recent infections. High transmission dynamics among Polish MSM emphasizes the necessity for active testing and early treatment in this group.

Identification of an HIV-1 BG Intersubtype Recombinant Form (CRF73_BG), Partially Related to CRF14_BG, Which Is Circulating in Portugal and Spain

HIV-1 exhibits a characteristically high genetic diversity, with the M group, responsible for the pandemic, being classified into nine subtypes, 72 circulating recombinant forms (CRFs) and numerous unique recombinant forms (URFs). Here we characterize the near full-length genome sequence of an HIV-1 BG intersubtype recombinant virus (X3208) collected in Galicia (Northwest Spain) which exhibits a mosaic structure coincident with that of a previously characterized BG recombinant virus (9601_01), collected in Germany and epidemiologically linked to Portugal, and different from currently defined CRFs. Similar recombination patterns were found in partial genome sequences from three other BG recombinant viruses, one newly derived, from a virus collected in Spain, and two retrieved from databases, collected in France and Portugal, respectively. Breakpoint coincidence and clustering in phylogenetic trees of these epidemiologically-unlinked viruses allow to define a new HIV-1 CRF (CRF73_BG). CRF73_BG shares one breakpoint in the envelope with CRF14_BG, which circulates in Portugal and Spain, and groups with it in a subtype B envelope fragment, but the greatest part of its genome does not appear to derive from CRF14_BG, although both CRFs share as parental strain the subtype G variant circulating in the Iberian Peninsula. Phylogenetic clustering of partial pol and env segments from viruses collected in Portugal and Spain with X3208 and 9691_01 indicates that CRF73_BG is circulating in both countries, with proportions of around 2–3% Portuguese database HIV-1 isolates clustering with CRF73_BG. The fact that an HIV-1 recombinant virus characterized ten years ago as a URF has been shown to represent a CRF suggests that the number of HIV-1 CRFs may be much greater than currently known.

Phylogeny and Phylogeography of a Recent HIV-1 Subtype F Outbreak among Men Who Have Sex with Men in Spain Deriving from a Cluster with a Wide Geographic Circulation in Western Europe

We recently reported the rapid expansion of an HIV-1 subtype F cluster among men who have sex with men (MSM) in the region of Galicia, Northwest Spain. Here we update this outbreak, analyze near full-length genomes, determine phylogenetic relationships, and estimate its origin. For this study, we used sequences of HIV-1 protease-reverse transcriptase and env V3 region, and for 17 samples, near full-length genome sequences were obtained. Phylogenetic analyses were performed via maximum likelihood. Locations and times of most recent common ancestors were estimated using Bayesian inference. Among samples analyzed by us, 100 HIV-1 F1 subsubtype infections of monophyletic origin were diagnosed in Spain, including 88 in Galicia and 12 in four other regions. Most viruses (n = 90) grouped in a subcluster (Galician subcluster), while 7 from Valladolid (Central Spain) grouped in another subcluster. At least 94 individuals were sexually-infected males and at least 71 were MSM. Seventeen near full-length genomes were uniformly of F1 subsubtype. Through similarity searches and phylogenetic analyses, we identified 18 viruses from four other Western European countries [Switzerland (n = 8), Belgium (n = 5), France (n = 3), and United Kingdom (n = 2)] and one from Brazil, from samples collected in 2005–2011, which branched within the subtype F cluster, outside of both Spanish subclusters, most of them corresponding to recently infected individuals. The most probable geographic origin and age of the Galician subcluster was Ferrol, Northwest Galicia, around 2007, while the Western European cluster probably emerged in Switzerland around 2002. In conclusion, a recently expanded HIV-1 subtype F cluster, the largest non-subtype B cluster reported in Western Europe, continues to spread among MSM in Spain; this cluster is part of a larger cluster with a wide geographic circulation in diverse Western European countries.

Automated analysis of phylogenetic clusters

BACKGROUND

As sequence data sets used for the investigation of pathogen transmission patterns increase in size, automated tools and standardized methods for cluster analysis have become necessary. We have developed an automated Cluster Picker which identifies monophyletic clades meeting user-input criteria for bootstrap support and maximum genetic distance within large phylogenetic trees. A second tool, the Cluster Matcher, automates the process of linking genetic data to epidemiological or clinical data, and matches clusters between runs of the Cluster Picker.

RESULTS

We explore the effect of different bootstrap and genetic distance thresholds on clusters identified in a data set of publicly available HIV sequences, and compare these results to those of a previously published tool for cluster identification. To demonstrate their utility, we then use the Cluster Picker and Cluster Matcher together to investigate how clusters in the data set changed over time. We find that clusters containing sequences from more than one UK location at the first time point (multiple origin) were significantly more likely to grow than those representing only a single location.

CONCLUSIONS

The Cluster Picker and Cluster Matcher can rapidly process phylogenetic trees containing tens of thousands of sequences. Together these tools will facilitate comparisons of pathogen transmission dynamics between studies and countries.

The global transmission network of HIV-1

Human immunodeficiency virus type 1 (HIV-1) is pandemic, but its contemporary global transmission network has not been characterized. A better understanding of the properties and dynamics of this network is essential for surveillance, prevention, and eventual eradication of HIV. Here, we apply a simple and computationally efficient network-based approach to all publicly available HIV polymerase sequences in the global database, revealing a contemporary picture of the spread of HIV-1 within and between countries. This approach automatically recovered well-characterized transmission clusters and extended other clusters thought to be contained within a single country across international borders. In addition, previously undescribed transmission clusters were discovered. Together, these clusters represent all known modes of HIV transmission. The extent of international linkage revealed by our comprehensive approach demonstrates the need to consider the global diversity of HIV, even when describing local epidemics. Finally, the speed of this method allows for near-real-time surveillance of the pandemic's progression.

HIV-1 Genetic Variability and Clinical Implications

Despite advances in antiretroviral therapy that have revolutionized HIV disease management, effective control of the HIV infection pandemic remains elusive. Beyond the classic non-B endemic areas, HIV-1 non-B subtype infections are sharply increasing in previous subtype B homogeneous areas such as Europe and North America. As already known, several studies have shown that, among non-B subtypes, subtypes C and D were found to be more aggressive in terms of disease progression. Luckily, the response to antiretrovirals against HIV-1 seems to be similar among different subtypes, but these results are mainly based on small or poorly designed studies. On the other hand, differences in rates of acquisition of resistance among non-B subtypes are already being observed. This different propensity, beyond the type of treatment regimens used, as well as access to viral load testing in non-B endemic areas seems to be due to HIV-1 clade specific peculiarities. Indeed, some non-B subtypes are proved to be more prone to develop resistance compared to B subtype. This phenomenon can be related to the presence of subtype-specific polymorphisms, different codon usage, and/or subtype-specific RNA templates. This review aims to provide a complete picture of HIV-1 genetic diversity and its implications for HIV-1 disease spread, effectiveness of therapies, and drug resistance development.

HIV prevalence and sexual risk behaviors associated with awareness of HIV status among men who have sex with men in Paris, France

A cross-sectional survey, using self-sampled finger-prick blood on blotting paper and anonymous behavioral self-administrated questionnaires was conducted in Paris in 2009 among MSM attending gay venues. Paired biological results and questionnaires were available for 886 participants. HIV seroprevalence was 17.7 % (95 % CI: 15.3-20.4). Four groups were identified according to their knowledge of their HIV biological status. Among the 157 found to be seropositive, 31 (19.7 %) were unaware of their status and reported high levels of sexual risk behaviors and frequent HIV testing in the previous 12 months. Among the 729 MSM diagnosed HIV-negative, 183 were no longer sure whether they were still HIV-negative, or had never been tested despite the fact that they engaged in at-risk sexual behaviors. This study provides the first estimate of HIV seroprevalence among MSM in Paris and underlines the specific need for combined prevention of HIV infection in this MSM population.

Reconstituting the epidemic history of HIV strain CRF01_AE among men who have sex with men (MSM) in Liaoning, northeastern China: implications for the expanding epidemic among MSM in China

HIV CRF01_AE accounted for 84% of the recent infections among men who have sex with men (MSM) in Liaoning Province of northeastern China. CRF01_AE strains were grouped into two distinct clusters (designated clusters 1 and 2) that were also detected in other regions in China. Phylodynamics study revealed that these two CRF01_AE strains were independently introduced into the population of MSM in China in the early and mid-1990s. Our study elucidated unique features of dynamics and interrelationships of MSM epidemics in China.

Molecular surveillance of HIV-1 in Madrid, Spain: a phylogeographic analysis

The molecular epidemiology of HIV-1 is constantly changing, mainly as a result of human migratory flows and the high adaptive ability of the virus. In recent years, Spain has become one of Europe's main destinations for immigrants and one of the western European countries with the highest rates of HIV-positive patients. Using a phylogeographic approach, we have analyzed the relationship between HIV-1 variants detected in immigrant and native populations of the urban area of Madrid. Our project was based on two coincidental facts. First, resistance tests were extended to naïve and newly diagnosed patients, and second, the Spanish government legislated the provision of legal status to many immigrants. This allowed us to obtain a large data set (n = 2,792) from 11 Madrid hospitals of viral pol sequences from the two populations, and with this unique material, we explored the impact of immigration in the epidemiological trends of HIV-1 variants circulating in the largest Spanish city. The prevalence of infections by non-B HIV-1 variants in the studied cohort was 9%, rising to 25% among native Spanish patients. Multiple transmission events involving different lineages and subsubtypes were observed in all the subtypes and recombinant forms studied. Our results also revealed strong social clustering among the most recent immigrant groups, such as Russians and Romanians, but not in those groups who have lived in Madrid for many years. Additionally, we document for the first time the presence of CRF47_BF and CRF38_BF in Europe, and a new BG recombinant form found in Spaniards and Africans is tentatively proposed. These results suggest that the HIV-1 epidemic will evolve toward a more complex epidemiological landscape.

Epidemiology of non-B clade forms of HIV-1 in men who have sex with men in the UK

OBJECTIVES

To describe the frequency and risk factors of non-B HIV-1 subtypes in men who have sex with men (MSM) in the UK.

DESIGN

Observational study.

METHODS

MSM diagnosed with HIV-1 infection from 1980-2007, with HIV genotype held in the UK HIV Drug Resistance Database were identified. Protease and reverse transcriptase sequences were collected and viral clade determined using the REGA algorithm. Associations between demographic variables and subtype were analysed using logistic regression.

RESULTS

The prevalence of non-B HIV-1 infection amongst MSM in the UK was 5.4% (437/8058). In the UK this increased with year of diagnosis from pre1996 to 2002, and has subsequently remained relatively stable at around 7-9% after 2002, with a recent increase in 2007 to 13%. Multivariate analysis showed that acquisition of non-B HIV-1 infection was independently associated with later year of HIV diagnosis (P < 0.001), black ethnicity (P < 0.001) and non-European country of birth (P = 0.01). Age was also associated with subtype with individuals aged 25-39 years being less likely to have non-B virus than those aged less than 25 years (P = 0.01). Restricting the analysis to white men born in the UK, the association between subtype and year of diagnosis remained statistically significant (P < 0.001), as did the association with age (P < 0.001).

DISCUSSION

The number of MSM in the UK infected with non-B clade HIV-1 is increasing, suggesting that the sociodemographic boundaries between HIV-1 viral subtypes globally are diminishing. Should viral subtypes be relevant to clinical disease progression or vaccine design, the changing pattern of distribution will need to be taken into account.

Short communication: Biological and genetic characterization of HIV type 1 subtype B and nonsubtype B transmitted viruses: usefulness for vaccine candidate assessment

Due to the extraordinary degree of genetic diversity of HIV-1 and the structural complexity of its envelope glycoproteins, designing an effective vaccine is difficult, requiring the development of viral reagents to assess vaccine-elicited neutralizing antibodies. The aim of this study was to improve on our previously developed panel of HIV-1 strains of different genetic forms, focusing on strains from acute and recently acquired infections as the most representative of the transmitted viruses. HIV-1 primary isolates were expanded in peripheral blood mononuclear cells. Viral stocks of 40 ml each were produced. Syncytium-inducing (SI) phenotype, coreceptor use, and TCID(50)/ml were determined. Near full-length HIV-1 genomes were amplified by RT-nested PCR in four overlapping segments. Phylogenetic analyses were performed with neighbor-joining trees and bootscanning. Forty-four HIV-1 strains were included in the panel. Twenty-four (54.1%) strains were from early infections (16 acute and 8 recent); of them, 21 (87%) were sexually transmitted. NSI/R5 phenotype was detected in 37 (84.1%) viruses and SI/R5,X4 in another 7 (15.9%). TCID(50)/ml ranged between 10(4) and 10(6.6). Twelve different genetic forms constituted this panel: subtypes A1, B, C, F1, and G; circulating recombinant forms CRF02_AG, CRF14_BG, and CRF24_BG; and unique recombinant forms CRF02_AG/A3, BF1, CRF12_BF/B, and DF1G. In conclusion, in this study, we report the development of a comprehensive and well-characterized panel of HIV-1 isolates for assessing neutralization in HIV vaccine research. This panel is available for distribution through the Programme EVA Centre for AIDS Reagents, National Institute for Biological Standard and Control (NIBSC).

Identification of a new HIV type 1 circulating BF intersubtype recombinant form (CRF47_BF) in Spain

We report the identification of a new HIV-1 circulating recombinant form (CRF47_BF) derived from subtypes B and F. It was initially identified in protease-reverse transcriptase sequences from nine individuals from three separate regions of Spain who acquired HIV-1 infection via sexual contact. All nine sequences formed a strongly supported phylogenetic cluster, branching apart from all known CRFs, and in bootscan analyses were BF mosaics with two coincident breakpoints. Two epidemiologically unlinked viruses were sequenced in near full-length genomes, which exhibited identical mosaic structures, with 16 intersubtype breakpoints in a genome predominantly of subtype B. Subtype F segments of the new CRF failed to cluster with any of the near full-length genome subtype F sequences available in public databases. Recent dates of HIV-1 diagnoses and short genetic distances suggest a recent origin of this CRF. This is the tenth reported CRF_BF, the first apparently having originated outside of South America.