Recent Evolutionary History of Human Immunodeficiency Virus Type 1 Subtype B: Reconstruction of Epidemic Onset Based on Sequence Distances to the Common Ancestor

We obtained and studied HIV-1 sequences with a known sampling year from three outbreaks of the HIV-1 epidemic: 141 env V3 (270 nt) sampled between 1984 and 1992 and 117 pol prot/RT (804 nt) sequences sampled between 1986 and 1999 from Dutch homosexual men and injecting drug users (IDUs), as well as 77 env V3 sequences sampled between 1983 and 1994 in the United States. Since retrospective serological and/or epidemiological data on these populations are available, providing estimates of the dates of the onset of the HIV-1 epidemics, we had the opportunity to test different phylogenetic models for their accuracy in deriving the recent evolutionary history of HIV-1 subtype B and the onset date of the HIV-1 epidemic. We observed that, in any given year, individual sequences vary widely in their distances to the common ancestor, and sequences close to the ancestors were found decades after the onset of the epidemic. Nevertheless, the mean evolutionary distances of virus strains to ancestors were increasing significantly during the course of the studied epidemics, which indicates that the molecular clock is operational in the recent evolution of HIV-1. When the relationship between the sampling years of sequences and their nucleotide distances to the common ancestor was extrapolated to the past, analysis of pol sequences provided accurate estimates of the onset years of the epidemics, whereas analysis of V3 sequences by the maximum-likelihood or neighbor-joining methods led to an overestimation of the age of the epidemics. Separate analysis of nonsynonymous and synonymous distances revealed that this overestimation results from nonsynonymous substitutions, whose numbers were not increasing significantly in all three virus populations over the observation period. In contrast, analysis of synonymous env V3 distances provided accurate estimates of the onset years for the outbreaks we studied. Received: 26 October 2001 / Accepted: 8 November 2001     

Molecular Evolution of the HIV-1 Thai Epidemic between the Time of RV144 Immunogen Selection to the Execution of the Vaccine Efficacy Trial

The RV144 HIV-1 vaccine trial (Thailand, 2003 to 2009), using immunogens genetically matched to the regional epidemic, demonstrated the first evidence of efficacy for an HIV-1 vaccine. Here we studied the molecular evolution of the HIV-1 epidemic from the time of immunogen selection to the execution of the efficacy trial. We studied HIV-1 genetic diversity among 390 volunteers who were deferred from enrollment in RV144 due to preexisting HIV-1 infection using a multiregion hybridization assay, full-genome sequencing, and phylogenetic analyses. The subtype distribution was 91.7% CRF01_AE, 3.5% subtype B, 4.3% B/CRF01_AE recombinants, and 0.5% dual infections. CRF01_AE strains were 31% more diverse than the ones from the 1990s Thai epidemic. Sixty-nine percent of subtype B strains clustered with the cosmopolitan Western B strains. Ninety-three percent of B/CRF01_AE recombinants were unique; recombination breakpoint analysis showed that these strains were highly embedded within the larger network that integrates recombinants from East/Southeast Asia. Compared to Thai sequences from the early 1990s, the distance to the RV144 immunogens increased 52% to 68% for CRF01_AE Env immunogens and 12% to 29% for subtype B immunogens. Forty-three percent to 48% of CRF01_AE sequences differed from the sequence of the vaccine insert in Env variable region 2 positions 169 and 181, which were implicated in vaccine sieve effects in RV144. In conclusion, compared to the molecular picture at the early stages of vaccine development, our results show an overall increase in the genetic complexity of viruses in the Thai epidemic and in the distance to vaccine immunogens, which should be considered at the time of the analysis of the trial results.     

Reviewing the history of HIV-1: spread of subtype B in the Americas

The dispersal of HIV-1 subtype B (HIV-1B) is a reflection of the movement of human populations in response to social, political, and geographical issues. The initial dissemination of HIV-1B outside Africa seems to have included the passive involvement of human populations from the Caribbean in spreading the virus to the United States. However, the exact pathways taken during the establishment of the pandemic in the Americas remain unclear. Here, we propose a geographical scenario for the dissemination of HIV-1B in the Americas, based on phylogenetic and genetic statistical analyses of 313 available sequences of the pol gene from 27 countries. Maximum likelihood and bayesian inference methods were used to explore the phylogenetic relationships between HIV-1B sequences, and molecular variance estimates were analyzed to infer the genetic structure of the viral population. We found that the initial dissemination and subsequent spread of subtype B in the Americas occurred via a single introduction event in the Caribbean around 1964 (1950-1967). Phylogenetic trees present evidence of several primary outbreaks in countries in South America, directly seeded by the Caribbean epidemic. Cuba is an exception insofar as its epidemic seems to have been introduced from South America. One clade comprising isolates from different countries emerged in the most-derived branches, reflecting the intense circulation of the virus throughout the American continents. Statistical analysis supports the genetic compartmentalization of the virus among the Americas, with a close relationship between the South American and Caribbean epidemics. These findings reflect the complex establishment of the HIV-1B pandemic and contribute to our understanding between the migration process of human populations and virus diffusion.     

中国HIV-1 B’亚型毒株不同基因区与全长基因组系统进化比较

本研究旨在探索不同基因区的使用对HIV-1 B’亚型毒株系统进化分析结果的影响。首先利用既往研究中已发表的共计47条来自泰国,缅甸和中国多个地区不同传播途径的B’毒株近全长基因组序列,将其按基因区分为不同的数据集 (gag, pol, vif, vpr, vpu, env, nef),并分别进行系统进化分析研究。比较不同基因区系统进化分析的结果发现,B’亚型毒株 pol基因在分析的基因区中,具有最低的复杂度和进化速率,可以较好的区分B’TH和B’YN毒株,重复近全长基因组序列的分析结果;尽管env基因则具有最高的复杂度和进化速率,但无法获得类似结果。本研究比较了不同基因区对HIV-1 B’亚型毒株系统进化分析结果的影响,对进一步开展HIV分子流行病学调查,分析我国B’毒株在我国的传播奠定了基础。     

Simple epidemiological dynamics explain phylogenetic clustering of HIV from patients with recent infection

Phylogenies of highly genetically variable viruses such as HIV-1 are potentially informative of epidemiological dynamics. Several studies have demonstrated the presence of clusters of highly related HIV-1 sequences, particularly among recently HIV-infected individuals, which have been used to argue for a high transmission rate during acute infection. Using a large set of HIV-1 subtype B pol sequences collected from men who have sex with men, we demonstrate that virus from recent infections tend to be phylogenetically clustered at a greater rate than virus from patients with chronic infection ('excess clustering') and also tend to cluster with other recent HIV infections rather than chronic, established infections ('excess co-clustering'), consistent with previous reports. To determine the role that a higher infectivity during acute infection may play in excess clustering and co-clustering, we developed a simple model of HIV infection that incorporates an early period of intensified transmission, and explicitly considers the dynamics of phylogenetic clusters alongside the dynamics of acute and chronic infected cases. We explored the potential for clustering statistics to be used for inference of acute stage transmission rates and found that no single statistic explains very much variance in parameters controlling acute stage transmission rates. We demonstrate that high transmission rates during the acute stage is not the main cause of excess clustering of virus from patients with early/acute infection compared to chronic infection, which may simply reflect the shorter time since transmission in acute infection. Higher transmission during acute infection can result in excess co-clustering of sequences, while the extent of clustering observed is most sensitive to the fraction of infections sampled.     

Evolutionary History of HIV-1 Subtype B and CRF01_AE Transmission Clusters among Men Who Have Sex with Men (MSM) in Kuala Lumpur,Malaysia

HIV-1 epidemics among men who have sex with men (MSM) continue to expand in developed and developing countries. Although HIV infection in MSM is amongst the highest of the key affected populations in many countries in Southeast Asia, comprehensive molecular epidemiological study of HIV-1 among MSM remains inadequate in the region including in Malaysia. Here, we reported the phylodynamic profiles of HIV-1 genotypes circulating among MSM population in Kuala Lumpur, Malaysia. A total of n = 459 newly-diagnosed treatment-naïve consenting subjects were recruited between March 2006 and August 2012, of whom 87 (18.9%) were self-reported MSM. Transmitted drug resistance mutations were absent in these isolates. Cumulatively, phylogenetic reconstructions of the pro-rt gene (HXB2∶2253–3275) showed that HIV-1 subtype B and CRF01_AE were predominant and contributed to approximately 80% of the total HIV-1 infection among MSM. In addition to numerous unique transmission lineages within these genotypes, twelve monophyletic transmission clusters of different sizes (2–7 MSM sequences, supported by posterior probability value of 1) were identified in Malaysia. Bayesian coalescent analysis estimated that the divergence times for these clusters were mainly dated between 1995 and 2005 with four major transmission clusters radiating at least 12 years ago suggesting that active spread of multiple sub-epidemic clusters occurred during this period. The changes in effective population size of subtype B showed an exponential growth within 5 years between 1988 and 1993, while CRF01_AE lineage exhibited similar expansion between 1993 and 2003. Our study provides the first insight of the phylodynamic profile of HIV-1 subtype B and CRF01_AE circulating among MSM population in Kuala Lumpur, Malaysia, unravelling the importance of understanding transmission behaviours as well as evolutionary history of HIV-1 in assessing the risk of outbreak or epidemic expansion.     

Phylogeographic analysis of HIV-1 subtype C dissemination in Southern Brazil

The HIV-1 subtype C has spread efficiently in the southern states of Brazil (Rio Grande do Sul, Santa Catarina and Paraná). Phylogeographic studies indicate that the subtype C epidemic in southern Brazil was initiated by the introduction of a single founder virus population at some time point between 1960 and 1980, but little is known about the spatial dynamics of viral spread. A total of 135 Brazilian HIV-1 subtype C pol sequences collected from 1992 to 2009 at the three southern state capitals (Porto Alegre, Florianópolis and Curitiba) were analyzed. Maximum-likelihood and Bayesian methods were used to explore the degree of phylogenetic mixing of subtype C sequences from different cities and to reconstruct the geographical pattern of viral spread in this country region. Phylogeographic analyses supported the monophyletic origin of the HIV-1 subtype C clade circulating in southern Brazil and placed the root of that clade in Curitiba (Paraná state). This analysis further suggested that Florianópolis (Santa Catarina state) is an important staging post in the subtype C dissemination displaying high viral migration rates from and to the other cities, while viral flux between Curitiba and Porto Alegre (Rio Grande do Sul state) is very low. We found a positive correlation (r(2) = 0.64) between routine travel and viral migration rates among localities. Despite the intense viral movement, phylogenetic intermixing of subtype C sequences from different Brazilian cities is lower than expected by chance. Notably, a high proportion (67%) of subtype C sequences from Porto Alegre branched within a single local monophyletic sub-cluster. These results suggest that the HIV-1 subtype C epidemic in southern Brazil has been shaped by both frequent viral migration among states and in situ dissemination of local clades.     

Evolution of the human immunodeficiency virus type 1 subtype-specific V3 domain is confined to a sequence space with a fixed distance to the subtype consensus.

Human immunodeficiency virus type 1 (HIV-1) strains can be separated into genetic subtypes based on phylogenetic analysis of the envelope gene. Once it had been shown that population-wide intrasubtype genetic variation of HIV-1 strains increases in the course of the AIDS epidemic, it remained uncertain whether HIV-1 subtypes are phenotypic entities spreading as distinct virus populations. To examine this, we applied Eigen's concepts of sequence geometry and fitness topography to the analysis of intrasubtype evolution of the gp120 V3 domain of HIV-1 subtypes A, B, C, and D in the course of the global AIDS epidemic. We observed that despite the high evolution rate of HIV-1, the nonsynonymous distances to the subtype consensus of sequences obtained early in the epidemic are similar to those obtained more than 10 years later, in contrast to the synonymous distances, which increased steadily over time. For HIV-1 subtype B, we observed that the evolution rate of the individual sequences is independent of their distance from the subtype B consensus, but for the individual sequences most distant from the consensus evolution away from the consensus is constrained. As a result, individual HIV-1 genomes fluctuate within a sequence space with fixed distance to the subtype consensus. Our findings suggest that the evolution of the V3 domain of HIV-1 subtypes A, B, C, and D is confined to an area in sequence space within a fixed distance to the consensus of a respective subtype. This in turn indicates that each HIV-1 subtype is a distinct viral quasispecies that is well adapted to the present environment, able to maintain its identity in the V3 region over time, and unlikely to merge during progression of the AIDS epidemic.     

Human Immunodeficiency Virus Type 1 (HIV-1) Subtype B Epidemic in Panama Is Mainly Driven by Dissemination of Country-Specific Clades

The Human immunodeficiency virus type-1 (HIV-1) subtype B is the most predominant clade in Central America; but information about the evolutionary history of this virus in this geographic region is scarce. In this study, we reconstructed the spatiotemporal and population dynamics of the HIV-1 subtype B epidemic in Panama. A total of 761 HIV-1 subtype B pol sequences obtained in Panama between 2004 and 2013 were combined with subtype B pol sequences from the Americas and Europe. Maximum Likelihood phylogenetic analyses revealed that HIV-1 subtype B infections in Panama derived from the dissemination of multiple founder viruses. Most Panamanian subtype B viruses (94.5%) belong to the pandemic viral strain proposed as originated in the US, whereas others (5.5%) were intermixed among non-pandemic Caribbean strains. The bulk (76.6%) of subtype B sequences from Panama grouped within 12 country-specific clades that were not detected in other Central American countries. Bayesian coalescent-based analyses suggest that most Panamanian clades probably originated between the early 1970s and the early 1980s. The root location of major Panamanian clades was traced to the most densely populated districts of Panama province. Major Panamanian clades appear to have experienced one or two periods of exponential growth of variable duration between the 1970s and the 2000s, with median growth rates from 0.2 to 0.4 year^{− 1}. Thus, the HIV-1 subtype B epidemic in Panama is driven by the expansion of local viral strains that were introduced from the Caribbean and other American countries at an early stage of the AIDS pandemic.     

Phylodynamics of HIV-1 subtype B among the men-having-sex-with-men (MSM) population in Hong Kong

The men-having-sex-with-men (MSM) population has become one of the major risk groups for HIV-1 infection in the Asia Pacific countries. Hong Kong is located in the centre of Asia and the transmission history of HIV-1 subtype B transmission among MSM remained unclear. The aim of this study was to investigate the transmission dynamics of HIV-1 subtype B virus in the Hong Kong MSM population. Samples of 125 HIV-1 subtype B infected MSM patients were recruited in this study. Through this study, the subtype B epidemic in the Hong Kong MSM population was identified spreading mainly among local Chinese who caught infection locally. On the other hand, HIV-1 subtype B infected Caucasian MSM caught infection mainly outside Hong Kong. The Bayesian phylogenetic analysis also indicated that 3 separate subtype B epidemics with divergence dates in the 1990s had occurred. The first and latest epidemics were comparatively small-scaled; spreading among the local Chinese MSM while sauna-visiting was found to be the major sex partner sourcing reservoir for the first subtype B epidemic. However, the second epidemic was spread in a large-scale among local Chinese MSM with a number of them having sourced their sex partners through the internet. The epidemic virus was estimated to have a divergence date in 1987 and the infected population in Hong Kong had a logistic growth throughout the past 20 years. Our study elucidated the evolutionary and demographic history of HIV-1 subtype B virus in Hong Kong MSM population. The understanding of transmission and growth model of the subtype B epidemic provides more information on the HIV-1 transmission among MSM population in other Asia Pacific high-income countries.     

HIV-1 Genetic Diversity and Its Impact on Baseline CD4+T Cells and Viral Loads among Recently Infected Men Who Have Sex with Men in Shanghai,China

The HIV-1 epidemic among men who have sex with men (MSM) has been spreading throughout China. Shanghai, a central gathering place for MSM, is facing a continuously increasing incidence of HIV-1 infection. In order to better understand the dynamics of HIV-1 diversity and its influence on patient’s immune status at baseline on diagnosis, 1265 newly HIV-1-infected MSM collected from January 2009 to December 2013 in Shanghai were retrospectively analyzed for genetic subtyping, CD4+T cell counts, and viral loads. HIV-1 phylogenetic analysis revealed a broad viral diversity including CRF01_AE (62.13%), CRF07_BC (24.51%), subtype B (8.06%), CRF55_01B (3.24%), CER67_01B (0.95%), CRF68_01B (0.4%), CRF08_BC (0.08%) and CRF59_01B (0.08%). Twenty-four unique recombination forms (URFs) (1.98%) were identified as well. Bayesian inference analysis indicated that the introduction of CRF01_AE strain (1997) was earlier than CRF07_BC strain (2001) into MSM population in Shanghai based on the time of the most recent common ancestor (tMRCA). Three epidemic clusters and five sub-clusters were found in CRF01_AE. Significantly lower CD4+T cell count was found in individuals infected with CRF01_AE than in those infected with CRF07_BC infection (P<0.01), whereas viral load was significantly higher those infected with CRF01_AE than with CRF07_BC (P<0.01). In addition, the patients with >45 years of age were found to have lower CD4+T cell counts and higher viral loads than the patients with <25 years of age (P<0.05). This study reveals the presence of HIV-1 subtype diversity in Shanghai and its remarkable influence on clinical outcome. A real-time surveillance of HIV-1 viral diversity and phylodynamics of epidemic cluster, patient’s baseline CD4+T cell count and viral load would be of great value to monitoring of disease progression, intervention for transmission, improvement of antiretroviral therapy strategy and design of vaccines.     

HIV-1 genetic diversity in Western Brittany,France

We describe human immunodeficiency virus type 1 (HIV-1) diversity in Western Brittany, France, and trace the dissemination of HIV-1 non-B subtype infection. The strategy for HIV-1 subtyping used involved subtype specific enzyme immunoassays, heteroduplex mobility assays and phylogenetic analysis of the sequences of env encoding the V3 loop region. Samples were obtained from 567 patients: 465 (82%) were of subtype B and 66 (11.6%) were not (20 were subtype A, 11 subtype C, four subtype D, seven subtype F, five subtype G and 19 others with circulating recombinant forms: 4CRF01_AE, 11CRF02_AG, 1H, 3CRF11_cpx). These findings are consistent with other studies of French populations. There is an epidemiological correlation between subtype B and homosexual or heterosexual contamination in France and between non-B subtype and heterosexual contamination in Africa.     

Characterization of HIV Transmission in South-East Austria

To gain deeper insight into the epidemiology of HIV-1 transmission in South-East Austria we performed a retrospective analysis of 259 HIV-1 partial pol sequences obtained from unique individuals newly diagnosed with HIV infection in South-East Austria from 2008 through 2014. After quality filtering, putative transmission linkages were inferred when two sequences were ≤1.5% genetically different. Multiple linkages were resolved into putative transmission clusters. Further phylogenetic analyses were performed using BEAST v1.8.1. Finally, we investigated putative links between the 259 sequences from South-East Austria and all publicly available HIV polymerase sequences in the Los Alamos National Laboratory HIV sequence database. We found that 45.6% (118/259) of the sampled sequences were genetically linked with at least one other sequence from South-East Austria forming putative transmission clusters. Clustering individuals were more likely to be men who have sex with men (MSM; p<0.001), infected with subtype B (p<0.001) or subtype F (p = 0.02). Among clustered males who reported only heterosexual (HSX) sex as an HIV risk, 47% clustered closely with MSM (either as pairs or within larger MSM clusters). One hundred and seven of the 259 sequences (41.3%) from South-East Austria had at least one putative inferred linkage with sequences from a total of 69 other countries. In conclusion, analysis of HIV-1 sequences from newly diagnosed individuals residing in South-East Austria revealed a high degree of national and international clustering mainly within MSM. Interestingly, we found that a high number of heterosexual males clustered within MSM networks, suggesting either linkage between risk groups or misrepresentation of sexual risk behaviors by subjects.     

The molecular epidemiological aspects of hiv-infection spreading in Perm region

Retrospective analysis of HIV-infection spreading in Perm region in conjunction with the genetic characterization of viral subtypes circulated on this territory from 1988 (when 1st case of infection was detected) until 2005 was performed. Analysis of epidemic process allowed to determine three periods of its development basing on both epidemic intensity and nature of circulating HIV-1 subtypes. During 1988 - 1996 (first period), when viral population was heterogenous (simultaneous circulation of three HIV-1 subtypes) with multiple routes of transmission, the epidemic process was characterized by low intensity. High incidence of HIV-infection among injection drug users and high homogeneity of circulated HIV-1 variants (98% of isolated variants belonged to HIV-1 subtype A with low level of genetic variability) were characteristics of the second period lasted from 1997 to 2001. Decrease in HIV-infection incidence in 2002-2005 was accompanied by the increase of HIV-1 transmission through heterosexual contacts and continuation of subtype A predominance between isolates. However increase in heterogeneity of viral population during this period, which manifested as increase of env and pol genes polymorphism, was detected.     

Detailed Molecular Epidemiologic Characterization of HIV-1 Infection in Bulgaria Reveals Broad Diversity and Evolving Phylodynamics

Limited information is available to describe the molecular epidemiology of HIV-1 in Bulgaria. To better understand the genetic diversity and the epidemiologic dynamics of HIV-1 we analyzed 125 new polymerase (pol) sequences from Bulgarians diagnosed through 2009 and 77 pol sequences available from our previous study from persons infected prior to 2007. Epidemiologic and demographic information was obtained from each participant and phylogenetic analysis was used to infer HIV-1 evolutionary histories. 120 (59.5%) persons were infected with one of five different HIV-1 subtypes (A1, B, C, F1 and H) and 63 (31.2%) persons were infected with one of six different circulating recombinant forms (CRFs; 01_AE, 02_AG, 04_cpx, 05_DF, 14_BG, and 36_cpx). We also for the first time identified infection with two different clusters of unique A-like and F-like sub-subtype variants in 12 persons (5.9%) and seven unique recombinant forms (3.5%), including a novel J/C recombinant. While subtype B was the major genotype identified and was more prevalent in MSM and increased between 2000–2005, most non-B subtypes were present in persons ≥45 years old. CRF01_AE was the most common non-B subtype and was higher in women and IDUs relative to other risk groups combined. Our results show that HIV-1 infection in Bulgaria reflects the shifting distribution of genotypes coincident with the changing epidemiology of the HIV-1 epidemic among different risk groups. Our data support increased public health interventions targeting IDUs and MSM. Furthermore, the substantial and increasing HIV-1 genetic heterogeneity, combined with fluctuating infection dynamics, highlights the importance of sustained and expanded surveillance to prevent and control HIV-1 infection in Bulgaria.     

Molecular and Phylogeographic Analysis of Human Immuno-deficiency Virus Type 1 Strains Infecting Treatment-naive Patients from Kigali, Rwanda

This study aimed at describing the genetic subtype distribution of HIV-1 strains circulating in Kigali and their epidemiological link with the HIV-1 strains from the five countries surrounding Rwanda. One hundred and thirty eight pol (RT and PR) sequences from 116 chronically- and 22 recently-infected antiretroviral therapy (ART)-naïve patients from Kigali were generated and subjected to HIV drug resistance (HIV-DR), phylogenetic and recombinant analyses in connection with 366 reference pol sequences from Rwanda, Burundi, Kenya, Democratic Republic of Congo, Tanzania and Uganda (Los Alamos database). Among the Rwandan samples, subtype A1 predominated (71.7%), followed by A1/C recombinants (18.1%), subtype C (5.8%), subtype D (2.9%), one A1/D recombinant (0.7%) and one unknown subtype (0.7%). Thirteen unique and three multiple A1/C recombinant forms were identified. No evidence for direct transmission events was found within the Rwandan strains. Molecular characteristics of HIV-1 were similar between chronically and recently-infected individuals and were not significantly associated with demographic or social factors. Our report suggests that the HIV-1 epidemic in Kigali is characterized by the emergence of A1/C recombinants and is not phylogenetically connected with the HIV-1 epidemic in the five neighboring countries. The relatively low level of transmitted HIV-DR mutations (2.9%) reported here indicates the good performance of the ART programme in Rwanda. However, the importance of promoting couples'' counseling, testing and disclosure during HIV prevention strategies is highlighted.     

河南郑州市男男同性恋HIV-1感染者gag区基因亚型分析

为了解河南郑州市男男同性恋人群HIV-1流行株的亚型分布情况,运用反转录及套式PCR方法从40例已经被确认为HIV-1阳性的MSM的全血样本中扩增gag全长基因并进行序列测定,应用BioEdit软件对序列进行校对编辑,利用MEGA3.1软件构件系统进化树,同时使用美国Los Alamos国家实验室HIV核酸序列库和美国NCBI提供的在线分析工具进行比对分析,确定基因亚型。结果成功获得24条gag基因序列,亚型分析结果显示,共存在B、CRF01-AE和CRF07-BC三种亚型,其中B亚型8例（33.33﹪）,CRF01-AE亚型10例（41.67﹪）,CRF07-BC亚型6例（25﹪）。小样本量的流行病学调查显示,河南郑州市同性恋人群中主要存在B,CRF01-AE和CRF07-BC三种亚型,CRF01-AE已成为河南郑州市MSM人群中HIV主要流行亚型,河南地区重组型毒株逐渐占优势,流行情况更加复杂。