A high HIV-1 strain variability in London, UK, revealed by full-genome analysis: Results from the ICONIC project

HIV subtypes and molecular transmission characteristics among elderly HIV-infected individuals aged 50 and above: a cross-sectional study in Huzhou City, Eastern China

OBJECTIVE

To analyse the HIV-1 subtypes and molecular transmission characteristics of HIV-infected older individuals aged 50 and above in Huzhou City, and provide a scientific basis for prevention and treatment strategies for them.

DESIGN

A cross-sectional study with clustered molecular transmission network cases was performed, and basic epidemiological information was retrieved from the Chinese Centres for Disease Prevention and Control (CDC) Information System.

SETTING AND PARTICIPANTS

A molecular epidemiological study was conducted in 899 newly diagnosed HIV-infected individuals from January 2019 and March 2023 in Huzhou city, Zhejiang province, Eastern China. Out of these, HIV sequences were successfully obtained from 673 individuals, including 274 who were older individuals aged 50 and above.

PRIMARY AND SECONDARY OUTCOMES

Reverse transcription-polymerase chain reaction (PCR) and nested PCR were used to amplify the polymerase gene of HIV-1, and gene sequencing was performed. We used univariate and multivariate logistic regression to describe the association of clustered molecular transmission network cases.

RESULTS

In total, 274 valid HIV sequences of older individuals were obtained, which revealed 14 subtypes. Circulating recombinant forms (CRF) 07_BC accounted for 55.8% and CRF01_AE accounted for 20.1% of the subtypes. Data of 150 older individuals were included in the molecular transmission network, and the proportion of elderly individuals in clustered cases is 52.26% (150/287). The results of multivariable logistic regression analysis showed that the older age group (60-82 years) and CRF07_BC subtype were associated with case clustering (transmission risk).

CONCLUSIONS

The key high-risk transmission network was mainly composed of the older age group (60-82 years) and CRF07_BC subtype. It is necessary to further strengthen AIDS health promotion and education for individuals aged 60 years and above, as well as for patients with the CRF07_BC subtype, to reduce HIV transmission and clustering risk.

Identification of a novel first-generation HIV-1 circulating recombinant form (CRF152_DG) among people living with HIV in Karachi, Pakistan

The objective of this study was to characterize a novel circulating recombinant form of human immunodeficiency virus type 1 (HIV-1) among people living with HIV in Karachi, Pakistan. We conducted near-full-length genome (NFLG) sequencing on eight samples exhibiting D/G recombination signals in the pol gene region. We successfully obtained NFLG sequences (790-9,614; with reference to the HXB2 genome) from four of the eight samples and then conducted phylogenetic and recombination analyses on them. The four NFLG sequences from our study and one DG unique recombinant form previously identified in the United Kingdom (GenBank accession: MF109700) formed a distinct monophyletic cluster with an Shimodaira-Hasegawa approximate likelihood ratio test node support value of 100%. Bootscan analyses of the five NFLG sequences of DG recombinants showed that all five NFLGs shared the same unique mosaic pattern of recombination breakpoints between D and G clades, with two D fragments in the pol and vif regions inserted into a G backbone. Subregion phylogenetic analyses confirmed these sequences to be a novel circulating recombinant form (CRF) composed of subtypes D and G. The DG recombinant sequences were eventually designated as CRF152_DG by the Los Alamos HIV Sequence Database staff.

IMPORTANCE

In Pakistan, the genetic diversity of human immunodeficiency virus type 1 (HIV-1) is becoming increasingly complex, compared to the early years of the epidemic that started after the detection of the first cases of HIV-1 in 1987 in Karachi. Based on the available molecular studies, two dominant HIV-1 clades, sub-subtype A1 and CRF02_AG, have been found to co-circulate with other clades, namely B, C, D, G, CRF01_AE, CRF35_A1D, and CRF56_cpx, in various urban areas of Pakistan. Several novel recombinant forms have also been detected. This first report of CRF152_DG highlights the complex nature of the HIV epidemic in Pakistan and emphasizes the importance of continual molecular surveillance (ideally based on whole-genome sequences) of HIV.

Identification of the Near Full-length Genome of a Novel HIV-1 CRF01_AE/CRF07_BC Recombinant with a Complex Genomic Structure Isolated in Hebei Province, China

BACKGROUND

During HIV genotypic drug resistance testing of patient samples in Baoding, Hebei Province, China, in 2022, a recombinant fragment was detected in the pol region of an HIV-1 strain.

OBJECTIVE

The objective of the study was to analyze the near full-length genome of a novel HIV-1 CRF01_AE/CRF07_BC recombinant with a complex genomic structure.

METHODS

Viral RNA was extracted from the blood of the infected individual and reverse transcribed to cDNA. Two overlapping segments of the HIV-1 genome were amplified using a nearendpoint dilution method and sequenced. Recombinant breakpoints were determined using RIP, jpHMM, and SimPlot 3.5.1 software. MEGA 6.0 software was used to construct a neighbor-joining phylogenetic tree.

RESULTS

We obtained the near full-length genome sequence (8680 bp) of a novel HIV-1 CRF01_AE/CRF07_BC recombinant. Recombination analysis showed that the genome comprised at least 12 overlapping segments, including six CRF07_BC and six CRF01_AE segments, with CRF07_BC as the backbone. The emergence of CRF01_AE/CRF07_BC recombinant strains indicated that HIV-1 co-infection is common. However, the increasing genetic complexity of the HIV-1 epidemic in China warrants continued investigation.

CONCLUSION

The increase in CRF01_AE/CRF07_BC recombinant viruses suggests that HIV-1 has a high genetic mutation rate in Hebei, China. This highlights the need for close monitoring of HIV-1 molecular epidemiologic changes to provide accurate, up-to-date information for effective disease control.

Identification of two near-identical novel HIV-1 unique recombinant forms (CRF01_AE/B) among men who have sex with men in baoding, hebei, China

Men who have sex with men (MSM) are the most frequent infection route of the human immunodeficiency virus (HIV) in Baoding, China, creating chances for the occurrence of unique recombinant forms (URFs) of the virus, i.e., recombination of different subtypes caused by co-circulation of multiple subtypes. In this report, two near-identical URFs (BDD002A and BDD069A) isolated from MSM in Baoding were identified. Phylogenetic tree analysis based on nearly full-length genomes (NFLGs) revealed that the two URFs formed a distinct monophyletic cluster with a bootstrap value of 100%. Recombinant breakpoints analysis identified that the NFLGs of BDD002A and BDD069A were both composed of CRF01_AE and subtype B, with six subtype B mosaic segments inserted into the CRF01_AE backbone. The CRF01_AE segments of the URFs clustered closely with the CRF01_AE reference sequences, and the B subregions clustered with the B reference sequences. The recombinant breakpoints of the two URFs were almost identical. These results suggest that effective interventions are urgently needed to prevent the formation of complex HIV-1 recombinant forms in Baoding, China.

Identification of a New HIV-1 Circulating Recombinant Form CRF112_01B Strain in Baoding City, Hebei Province, China

BACKGROUND

A large number of HIV-1 recombinants that originated from CRF01_AE and B strains are constantly emerging in men who have sex with men populations in China and deserve more attention and further monitoring.

OBJECTIVE

To analyze the near-full-length genome structure and recombination characteristics of a new HIV-1 strain (BD226AJ) detected in Baoding City and determine its subtype.

CASE REPRESENTATION

Viral RNA was extracted from a blood sample collected from an infected individual and reverse transcribed to cDNA. Two overlapping segments of the HIV-1 genome were amplified using a near-endpoint dilution method and sequenced. Recombinant breakpoints were determined using RIP, jpHMM, and SimPlot 3.5.1 software. MEGA v6.0 was used to construct a neighbor-joining phylogenetic tree to determine the homology relationships of this strain.

RESULTS AND DISCUSSION

We obtained 8830 nucleotides (nt) of the HIV-1 genome sequence by amplification and sequencing, and four recombinant fragments were identified by recombination analysis, namely CRF01_AE (HXB2, 823-4224 nt), subtype B (HXB2, 4225-5991 nt), CRF01_AE (HXB2, 5992-9295 nt), and subtype B (HXB2, 9296-9406 nt). The BLAST results showed that 96% of the sequence was similar to CRF112_01B. The jpHMM results confirmed that BD226AJ was the CRF112_01B strain.

CONCLUSION

Our results confirm the first epidemic of CRF112_01B in Hebei Province. This finding suggests that HIV-1 CRF112_01B may have been introduced into Hebei by men who have sex with men and indicates that the epidemic trend of this strain should be closely monitored.

Recombination in Coronaviruses, with a Focus on SARS-CoV-2

Recombination is a common evolutionary tool for RNA viruses, and coronaviruses are no exception. We review here the evidence for recombination in SARS-CoV-2 and reconcile nomenclature for recombinants, discuss their origin and fitness, and speculate how recombinants could make a difference in the future of the COVID-19 pandemics.

Revisiting the recombinant history of HIV-1 group M with dynamic network community detection

Recombination is a major mechanism through which HIV type 1 (HIV-1) maintains genetic diversity and interferes with viral eradication efforts. There is growing evidence demonstrating a recombinant origin of primate lentiviruses including HIV-1 group M (HIV-1/M). Inferring the extent of recombination across the entire HIV-1/M genome is of great importance as it provides deeper insights into the origin, dynamics, and evolution of the global pandemic. Here we propose an alternative method that can reconstruct the extent of genome-wide recombination in HIV-1, uncover reticulate patterns, and serve as a framework for HIV-1 classification. Our method provides an alternative approach for understanding the roles of virus recombination in the early evolutionary history of zoonosis for other emerging viruses.

The prevailing abundance of full-length HIV type 1 (HIV-1) genome sequences provides an opportunity to revisit the standard model of HIV-1 group M (HIV-1/M) diversity that clusters genomes into largely nonrecombinant subtypes, which is not consistent with recent evidence of deep recombinant histories for simian immunodeficiency virus (SIV) and other HIV-1 groups. Here we develop an unsupervised nonparametric clustering approach, which does not rely on predefined nonrecombinant genomes, by adapting a community detection method developed for dynamic social network analysis. We show that this method (dynamic stochastic block model [DSBM]) attains a significantly lower mean error rate in detecting recombinant breakpoints in simulated data (quasibinomial generalized linear model (GLM), P<8×10−8), compared to other reference-free recombination detection programs (genetic algorithm for recombination detection [GARD], recombination detection program 4 [RDP4], and RDP5). When this method was applied to a representative sample of n = 525 actual HIV-1 genomes, we determined k = 29 as the optimal number of DSBM clusters and used change-point detection to estimate that at least 95% of these genomes are recombinant. Further, we identified both known and undocumented recombination hotspots in the HIV-1 genome and evidence of intersubtype recombination in HIV-1 subtype reference genomes. We propose that clusters generated by DSBM can provide an informative framework for HIV-1 classification.

The Emergence and Transmission Dynamics of HIV-1 CRF07_BC in Mainland China

A total of 1155 partial pol gene sequences of human immunodeficiency virus (HIV)-1 CRF07_BC were sampled between 1997 and 2015, spanning 13 provinces in Mainland China and risk groups [heterosexual, injecting drug users (IDU), and men who have sex with men (MSM)] to investigate the evolution, adaptation, spatiotemporal and risk group dynamics, migration patterns, and protein structure of HIV-1 CRF07_BC. Due to the unequal distribution of sequences across time, location, and risk group in the complete dataset (‘full1155’), subsampling methods were used. Maximum-likelihood and Bayesian phylogenetic analysis as well as discrete trait analysis of geographical location and risk group were carried out. To study mutations of a cluster of HIV-1 CRF07_BC (CRF07-1), we performed a comparative analysis of this cluster to the other CRF07_BC sequences (‘backbone_295’) and mapped the mutations observed in the respective protein structure. Our findings showed that HIV-1 CRF07_BC most likely originated among IDU in Yunnan Province between October 1992 to July 1993 [95 per cent hightest posterior density (HPD): May 1989–August 1995] and that IDU in Yunnan Province and MSM in Guangdong Province likely served as the viral sources during the early and more recent spread in Mainland China. We also revealed that HIV-1 CRF07-1 has been spreading for roughly 20 years and continues to cause local transmission in Mainland China and worldwide. Overall, our study sheds light on the dynamics of HIV-1 CRF07_BC distribution patterns in Mainland China. Our research may also be useful in formulating public health policies aimed at controlling acquired immune deficiency syndrome in Mainland China and globally.

Assessing hepatitis C virus distribution among vulnerable populations in London using whole genome sequencing: results from the TB-REACH study

Background: Injecting drugs substantially increases the risk of hepatitis C virus (HCV) infection and is common in vulnerable population groups, such as the homeless and prisoners. Capturing accurate data on relative genotype distribution within these groups is essential to inform strategies to reduce HCV transmission. The aim of this study was to utilise a next-generation whole-genome sequencing method recently validated by Public Health England, in order to produce near complete HCV genomes. Methods: In total, 98 HCV positive patients were recruited from homeless hostels and drug treatment services through the National Health Services (NHS) Find and Treat (F&T) Service between May 2011 and June 2013 in London, UK. Samples were sequenced by Next-generation sequencing, with 88 complete HCV genomes constructed by a de novo assembly pipeline. They were analysed phylogenetically for an estimate of their genetic distance. Results: Of the 88 complete HCV genomes, 50/88 (56.8%) were genotype 1; 32/88 (36.4%) genotype 3; 4/88 (4.5%) genotype 2; and 1/88 (1.1%) for genotypes 4 and 6 each. Subtype 1a had the highest number of samples (51.1%), followed by subtype 3a (35.2%), 1b (5.7%), and 2b (3.4%). Samples collected from drug treatment services had the highest number of genotype 1 (69%); genotypes 4 and 6 were only found from samples collected in homeless shelters. Small clusters of highly related genomic sequences were observed both across and within the vulnerable groups sampled. Conclusions: Subsequent phylogenetic analysis provides a first indication that there are related HCV sequences amongst the three vulnerable population groups, reflecting their overlapping social behaviours. This study is the first presentation of whole genome HCV sequences from such vulnerable groups in London and paves the way for similar research in the future.

Identification of CRF89_BF, a new member of an HIV-1 circulating BF intersubtype recombinant form family widely spread in South America

Circulating recombinant forms (CRFs) contribute substantially to the HIV-1 pandemic. Among 105 CRFs described in the literature, 16 are BF intersubtype recombinants, most of South American origin, of which CRF12_BF is the most widely spread. A BF recombinant cluster identified in Bolivia was suggested to represent a new CRF_BF. Here we find that it belongs to a larger cluster incorporating 39 viruses collected in 7 countries from 3 continents, 22 of them in Spain, most from Bolivian or Peruvian individuals, and 12 in South America (Bolivia, Argentina, and Peru). This BF cluster comprises three major subclusters, two associated with Bolivian and one with Peruvian individuals. Near full-length genome sequence analyses of nine viruses, collected in Spain, Bolivia, and Peru, revealed coincident BF mosaic structures, with 13 breakpoints, 6 and 7 of which coincided with CRF12_BF and CRF17_BF, respectively. In a phylogenetic tree, they grouped in a clade closely related to these CRFs, and more distantly to CRF38_BF and CRF44_BF, all circulating in South America. These results allowed to identify a new HIV-1 CRF, designated CRF89_BF. Through phylodynamic analyses, CRF89_BF emergence was estimated in Bolivia around 1986. CRF89_BF is the fifth CRF member of the HIV-1 recombinant family related to CRF12_BF.

Genetic Diversity of HIV-1 in Krasnoyarsk Krai: Area with High Levels of HIV-1 Recombination in Russia

More than a quarter of HIV-infected individuals registered in Russia live in Siberia. Unlike Central Russia where HIV-1 subtype A6 is predominant, in most Siberian regions since 2012, a new HIV-1 CRF63_02A1 genetic variant has spread, with the share of this variant attaining 75-85% among newly identified HIV cases. Krasnoyarsk Krai is considered to be a high-risk territory according to morbidity rate and HIV infection incidence among the population. The current paper aims to study the molecular epidemiologic characteristics of HIV-1 spreading in Krasnoyarsk Krai. Phylogenetic and recombination analyses of pol (PR-RT, IN) and env regions of the virus were used for genotyping 159 HIV-1 isolated in Krasnoyarsk Krai. 57.2% of the isolates belonged to subtype A (A6) specific to Russia, 12.6% to CRF63_02A1, and 0.6% to CRF02_AGСА, and in 29.6% HIV-1 URFs were detected, including URF63/А (23.9%), URFА/В (4.4%), and URF02/А (1.3%). In 6 of 7, HIV-1 URFА/В identical recombination model was detected; the origin of 38 URF63/А was proven to be the result of individual recombination events. Since 2015, a share of the population with newly diagnosed HIV who were infected with HIV-1 URF reached an exceptionally high rate of 38.6%. As distinct from adjacent Siberian regions, the HIV-1 CRF63_02A1 prevalence rate in Krasnoyarsk Krai is within 16%; however, the increased contribution of new HIV-1 into the regional epidemic development was observed due to the recombination of viruses of subtypes А, В, and CRF63_02A1. The difference between the described molecular epidemiologic picture in Krasnoyarsk Krai and in adjacent areas is likely caused by differences in predominant routes of HIV transmission and by more recent HIV-1 CRF63_02A1 transmission in the PWID group, which had a high prevalence of HIV-1 subtype A by the time of the new virus transmission, resulting in increased possibility of coinfection with various HIV-1 genetic variants.

Factors influencing HIV-1 phylogenetic clustering

PURPOSE OF REVIEW

A major goal of public health in relation to HIV/AIDS is to prevent new transmissions in communities. Phylogenetic techniques have improved our understanding of the structure and dynamics of HIV transmissions. However, there is still no consensus about phylogenetic methodology, sampling coverage, gene target and/or minimum fragment size.

RECENT FINDINGS

Several studies use a combined methodology, which includes both a genetic or patristic distance cut-off and a branching support threshold to identify phylogenetic clusters. However, the choice about these thresholds remains an inherently subjective process, which affects the results of these studies. There is still a lack of consensus about the genomic region and the size of fragments that should be used, although there seems to be emerging a consensus that using longer segments, allied with the use of a realistic model of evolution and a codon alignment, increases the likelihood of inferring true transmission clusters. The pol gene is still the most used genomic region, but recent studies have suggested that whole genomes and/or sequences from nef and gp41 are also good targets for cluster reconstruction.

SUMMARY

The development and application of standard methodologies for phylogenetic clustering analysis will advance our understanding of factors associated with HIV transmission. This will lead to the design of more precise public health interventions.

Phylogenetic and Demographic Characterization of Directed HIV-1 Transmission Using Deep Sequences from High-Risk and General Population Cohorts/Groups in Uganda

Across sub-Saharan Africa, key populations with elevated HIV-1 incidence and/or prevalence have been identified, but their contribution to disease spread remains unclear. We performed viral deep-sequence phylogenetic analyses to quantify transmission dynamics between the general population (GP), fisherfolk communities (FF), and women at high risk of infection and their clients (WHR) in central and southwestern Uganda. Between August 2014 and August 2017, 6185 HIV-1 positive individuals were enrolled in 3 GP and 10 FF communities, 3 WHR enrollment sites. A total of 2531 antiretroviral therapy (ART) naïve participants with plasma viral load >1000 copies/mL were deep-sequenced. One hundred and twenty-three transmission networks were reconstructed, including 105 phylogenetically highly supported source–recipient pairs. Only one pair involved a WHR and male participant, suggesting that improved population sampling is needed to assess empirically the role of WHR to the transmission dynamics. More transmissions were observed from the GP communities to FF communities than vice versa, with an estimated flow ratio of 1.56 (95% CrI 0.68–3.72), indicating that fishing communities on Lake Victoria are not a net source of transmission flow to neighboring communities further inland. Men contributed disproportionally to HIV-1 transmission flow regardless of age, suggesting that prevention efforts need to better aid men to engage with and stay in care.

Near Full-length Genomic Sequencing and Molecular Analysis of HIV-Infected Individuals in a Network-based Intervention (TRIP) in Athens, Greece: Evidence that Transmissions Occur More Frequently from those with High HIV-RNA

Background:

TRIP (Transmission Reduction Intervention Project) was a network-based, contact tracing approach to locate and link to care, mostly people who inject drugs (PWID) with recent HIV infection.

Objective:

We investigated whether sequences from HIV-infected participants with high viral load cluster together more frequently than what is expected by chance.

Methods:

Paired end reads were generated for 104 samples using Illumina MiSeq next-generation se-quencing.

Results:

63 sequences belonged to previously identified local transmission networks of PWID (LTNs) of an HIV outbreak in Athens, Greece. For two HIV-RNA cut-offs (105 and 106 IU/mL), HIV transmissions were more likely between PWID with similar levels of HIV-RNA (p<0.001). 10 of the 14 sequences (71.4%) from PWID with HIV-RNA >106 IU/mL were clustered in 5 pairs. For 4 of these clusters (80%), there was in each one of them at least one sequence from a recently HIV-infected PWID.

Conclusion:

We showed that transmissions are more likely among PWID with high viremia.

Bigger and Better? Representativeness of the Influenza A Surveillance Using One Consolidated Clinical Microbiology Laboratory Data Set as Compared to the Belgian Sentinel Network of Laboratories

Infectious diseases remain a serious public health concern globally, while the need for reliable and representative surveillance systems remains as acute as ever. The public health surveillance of infectious diseases uses reported positive results from sentinel clinical laboratories or laboratory networks, to survey the presence of specific microbial agents known to constitute a threat to public health in a given population. This monitoring activity is commonly based on a representative fraction of the microbiology laboratories nationally reporting to a single central reference point. However, in recent years a number of clinical microbiology laboratories (CML) have undergone a process of consolidation involving a shift toward laboratory amalgamation and closer real-time informational linkage. This report aims to investigate whether such merging activities might have a potential impact on infectious diseases surveillance. Influenza data was used from Belgian public health surveillance 2014–2017, to evaluate whether national infection trends could be estimated equally as effectively from only just one centralized CML serving the wider Brussels area (LHUB-ULB). The overall comparison reveals that there is a close correlation and representativeness of the LHUB-ULB data to the national and international data for the same time periods, both on epidemiological and molecular grounds. Notably, the effectiveness of the LHUB-ULB surveillance remains partially subject to local regional variations. A subset of the Influenza samples had their whole genome sequenced so that the observed epidemiological trends could be correlated to molecular observations from the same period, as an added-value proposition. These results illustrate that the real-time integration of high-throughput whole genome sequencing platforms available in consolidated CMLs into the public health surveillance system is not only credible but also advantageous to use for future surveillance and prediction purposes. This can be most effective when implemented for automatic detection systems that might include multiple layers of information and timely implementation of control strategies.

PCR Amplification Strategies Towards Full-length HIV-1 Genome Sequencing

The advent of next-generation sequencing has enabled greater resolution of viral diversity and improved feasibility of full viral genome sequencing allowing routine HIV-1 full genome sequencing in both research and diagnostic settings. Regardless of the sequencing platform selected, successful PCR amplification of the HIV-1 genome is essential for sequencing template preparation. As such, full HIV-1 genome amplification is a crucial step in dictating the successful and reliable sequencing downstream. Here we reviewed existing PCR protocols leading to HIV-1 full genome sequencing. In addition to the discussion on basic considerations on relevant PCR design, the advantages as well as the pitfalls of the published protocols were reviewed.

Diverse Large HIV-1 Non-subtype B Clusters Are Spreading Among Men Who Have Sex With Men in Spain

In Western Europe, the HIV-1 epidemic among men who have sex with men (MSM) is dominated by subtype B. However, recently, other genetic forms have been reported to circulate in this population, as evidenced by their grouping in clusters predominantly comprising European individuals. Here we describe four large HIV-1 non-subtype B clusters spreading among MSM in Spain. Samples were collected in 9 regions. A pol fragment was amplified from plasma RNA or blood-extracted DNA. Phylogenetic analyses were performed via maximum likelihood, including database sequences of the same genetic forms as the identified clusters. Times and locations of the most recent common ancestors (MRCA) of clusters were estimated with a Bayesian method. Five large non-subtype B clusters associated with MSM were identified. The largest one, of F1 subtype, was reported previously. The other four were of CRF02_AG (CRF02_1; n = 115) and subtypes A1 (A1_1; n = 66), F1 (F1_3; n = 36), and C (C_7; n = 17). Most individuals belonging to them had been diagnosed of HIV-1 infection in the last 10 years. Each cluster comprised viruses from 3 to 8 Spanish regions and also comprised or was related to viruses from other countries: CRF02_1 comprised a Japanese subcluster and viruses from 8 other countries from Western Europe, Asia, and South America; A1_1 comprised viruses from Portugal, United Kingom, and United States, and was related to the A1 strain circulating in Greece, Albania and Cyprus; F1_3 was related to viruses from Romania; and C_7 comprised viruses from Portugal and was related to a virus from Mozambique. A subcluster within CRF02_1 was associated with heterosexual transmission. Near full-length genomes of each cluster were of uniform genetic form. Times of MRCAs of CRF02_1, A1_1, F1_3, and C_7 were estimated around 1986, 1989, 2013, and 1983, respectively. MRCA locations for CRF02_1 and A1_1 were uncertain (however initial expansions in Spain in Madrid and Vigo, respectively, were estimated) and were most probable in Bilbao, Spain, for F1_3 and Portugal for C_7. These results show that the HIV-1 epidemic among MSM in Spain is becoming increasingly diverse through the expansion of diverse non-subtype B clusters, comprising or related to viruses circulating in other countries.