Original Research ARTICLE
Dynamics of HIV-1 Molecular Networks Reveal Effective Control of Large Transmission Clusters in an Area Affected by an Epidemic of Multiple HIV Subtypes
- 1NHC Key Laboratory of AIDS Immunology (China Medical University), National Clinical Research Center for Laboratory Medicine, The First Affiliated Hospital of China Medical University, Shenyang, China
- 2Units of Medical Laboratory, Chinese Academy of Medical Sciences, Shenyang, China
- 3Key Laboratory of AIDS Immunology, Chinese Academy of Medical Sciences, Shenyang, China
- 4Key Laboratory of AIDS Immunology of Liaoning Province, Shenyang, China
- 5Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou, China
This study reconstructed molecular networks of human immunodeficiency virus (HIV) transmission history in an area affected by an epidemic of multiple HIV-1 subtypes and assessed the efficacy of strengthened early antiretroviral therapy (ART) and regular interventions in preventing HIV spread. We collected demographic and clinical data of 2221 treatment-naïve HIV-1–infected patients in a long-term cohort in Shenyang, Northeast China, between 2008 and 2016. HIV pol gene sequencing was performed and molecular networks of CRF01_AE, CRF07_BC, and subtype B were inferred using HIV-TRACE with separate optimized genetic distance threshold. We identified 168 clusters containing ≥ 2 cases among CRF01_AE-, CRF07_BC-, and subtype B-infected cases, including 13 large clusters (≥ 10 cases). Individuals in large clusters were characterized by younger age, homosexual behavior, more recent infection, higher CD4 counts, and delayed/no ART (P < 0.001). The dynamics of large clusters were estimated by proportional detection rate (PDR), cluster growth predictor, and effective reproductive number (Re). Most large clusters showed decreased or stable during the study period, indicating that expansion was slowing. The proportion of newly diagnosed cases in large clusters declined from 30 to 8% between 2008 and 2016, coinciding with an increase in early ART within 6 months after diagnosis from 24 to 79%, supporting the effectiveness of strengthened early ART and continuous regular interventions. In conclusion, molecular network analyses can thus be useful for evaluating the efficacy of interventions in epidemics with a complex HIV profile.
The rapid evolution of human immunodeficiency virus (HIV) leaves measurable footprints in the viral genome that can be used for epidemic surveillance by phylogenetic analysis (Hassan et al., 2017). In recent years, a simplified genetic distance-based method has increasingly been used to infer HIV-1 networks in the population, in which a molecular cluster represents a group of individuals infected with genetically similar HIV strains (Smith et al., 2009; US-CDC, 2018; Han et al., 2020). The expansion of molecular clusters represents recent and ongoing HIV transmission, and key subpopulations associated with such clusters are targets for prioritized interventions such as partner tracing and HIV testing through partner services (Green et al., 2017) for diagnosis of unknown HIV-positive cases (Green et al., 2017). Immediate antiretroviral therapy (ART) is recommended for individuals diagnosed with HIV infection (Hoenigl et al., 2016b), along with pre-exposure prophylaxis (PrEP) for socially linked individuals who test negative (Kasaie et al., 2017).
HIV molecular networks can be used not only to guide targeted intervention in key subpopulations (Aldous et al., 2012; Lubelchek et al., 2015; Castley et al., 2017; Chaillon et al., 2017; Zhang et al., 2017; Stecher et al., 2018; Volz et al., 2018; Ragonnet-Cronin et al., 2019; Board et al., 2020), but also to reconstruct the history of HIV spread between populations (Kouyos et al., 2010; Brenner et al., 2017; Patino-Galindo et al., 2017; Pineda-Pena et al., 2018; Delgado et al., 2019; Fabeni et al., 2019). Several local and national studies conducted in recent years have monitored the dynamics of HIV molecular clusters and evaluated their expansion speed (Chan et al., 2015; Mehta et al., 2018; Ragonnet-Cronin et al., 2018; Wertheim et al., 2018; Jovanovic et al., 2019; Dennis et al., 2020). However, most of these studies focused on areas with only subtype B or C HIV-1 epidemics, and few local studies have reconstructed HIV transmission history in an area with multiple HIV subtypes, which may be more complicated because of the variable transmission dynamics and evolution rates of different HIV-1 strains in the population.
China is among the countries with the highest numbers of HIV-1 subtypes in the world outside of West and Central Africa (Li et al., 2016; Hemelaar et al., 2019). According to the latest National HIV Molecular Epidemiological Survey in China, as many as 18 known HIV-1 subtypes and circulating recombinant forms (CRFs) have been detected, with four main subtypes accounting for 90% of total infections: CRF07_BC (41.9%), CRF01_AE (33.2%), CRF08_BC (10.9%), and subtype B (4.0%) (Li et al., 2016). Shenyang, a city in northeastern China, records about 1000 new cases of HIV infection annually, which represents a mid-level HIV incidence in the country (Li et al., 2016). However, Shenyang has experienced multiple–HIV-subtype epidemics, initially in the heterosexual population and later in the homosexual population (Han et al., 2010; Han et al., 2013). Chinese guidelines for ART initiation in HIV-infected patients have been updated several times: in 2002, when China’s free ART program was introduced, it was recommended for the World Health Organization stage III or IV, symptomatic disease, extrapulmonary tuberculosis, or CD4 + T cell counts < 200 cells/μl; this was amended to CD4 + T cell counts ≤ 350 cells/μl in 2008; CD4 + T cell counts ≤ 500 cells/μl in 2014; and finally, to a recommendation of immediate treatment for all diagnosed cases of HIV infection in 2016 (China-CDC, 2005; AIDS-Professional-Group, 2011; AIDS-Professional-Group, 2015; China-CDC, 2018; Xu et al., 2020). In 2008, a large-scale prospective cohort of men who have sex with men (MSM) was established in Shenyang that included thousands of HIV-1–negative MSM who were regularly followed up and screened for HIV infection status through serologic and pooled nucleic acid testing (Xu et al., 2010). The impact of the abovementioned strengthened treatment and intervention policies on the local HIV epidemic has yet to be systematically evaluated.
In this study, we reconstructed the HIV-1 molecular networks of three major subtypes of HIV based on partial pol gene sequences of HIV-1 who were newly diagnosed with HIV infection between 2008 and 2016 in Shenyang. Demographic and clinical data were also analyzed to characterize the cases associated with larger clusters. We then evaluated the expansion dynamics of each large cluster (≥ 10 cases) and assessed the effects of strengthened early ART and regular interventions on the local HIV epidemic.
Materials and Methods
The study enrolled 2221 individuals with newly diagnosed HIV infection at the First Affiliated Hospital of China Medical University from 2008 to 2016. This hospital is the largest general hospital in Shenyang city, and admits nearly half of all HIV infection cases. All individuals who were diagnosed at or were referred to the hospital for treatment between 2008 and 2016 were included in the study. Blood samples were collected at the time of diagnosis or before ART. Demographic data including sex, age, ethnic group, occupation, education, marital status, HIV risk behaviors, date of diagnosis, date of ART initiation, and resident city; and clinical data including viral load and CD4 + T cell count were collected. The study was approved by the ethics committee of the First Affiliated Hospital of China Medical University. All study participants signed informed consent forms.
HIV-1 Limiting Antigen (Lag) Avidity Enzyme Immunoassay
Recent HIV infection (RHI) was distinguished from chronic HIV infection (CHI) using the LAg-Avidity EIA kit (Maxim Biomedical, Rockville, MD, United States) according to the manufacturer’s instructions. The normalized optical density (OD) of each sample was calculated as OD of the sample divided by that of the calibrator. RHI was defined as OD ≤ 2.0 in the screening test, and OD ≤ 1.5 in the confirmatory test (Kouyos et al., 2010).
HIV-1 Sequences and Cluster Identification
RNA extraction and partial pol gene amplification and sequencing were performed as previously described (Zhao et al., 2011). The sequences were aligned using the online HIVAlign program1 and manually edited. The retention length was 1015 bp (HXB2: 2253–3267). HIV-1 subtypes were determined by phylogenetic analysis after constructing an approximate maximum likelihood tree using Fast Tree 3.0 (Price et al., 2010), in which subtype N was used as the outgroup, the nucleotide substitution model was GTR + G + I, and support values of the nodes were calculated with a Shimodaira Hasegawa-like test (Liu et al., 2011). HIV-1 molecular networks were constructed using HIV Transmission Cluster Engine (HIV-TRACE) (Kosakovsky Pond et al., 2018)2 according to a previously described protocol (Wertheim et al., 2014; Oster et al., 2015; Whiteside et al., 2015; Wertheim et al., 2016). Briefly, all sequences were aligned with a reference HIV-1 pol sequence and the Tamura–Nei 93 pairwise distance was calculated for each pair of sequences. To obtain a high-resolution molecular network, we optimized the genetic distance threshold of three major subtypes to identify the largest number of molecular clusters (Wertheim et al., 2017). Pairwise distances of 0.5%, 0.5%, and 0.7% were used as the optimized genetic thresholds for CRF01_AE, CRF07_BC, and subtype B, respectively (Supplementary Figure S3 and Supplementary Table S1). All codons associated with antiretroviral drug resistance were included in this study. Previous studies have shown that the transmission of drug resistance in Shenyang occurs at a low rate (Zhao et al., 2011; Zhao, 2015). We removed codons associated with antiretroviral drug resistance and found that the results were unchanged (data not shown).
Relationship Between Large Clusters and Others
Clusters with ≥ 10 cases and between two and nine cases were defined as large and small/medium clusters, respectively. The definition of a large cluster is similar to that used in previous studies (Hughes et al., 2009; Leigh Brown et al., 2011; Patino-Galindo et al., 2017; Lorenzin et al., 2019; Rhee et al., 2019; Dennis et al., 2020). We evaluated demographic features of the study population including the time of diagnosis, sex, risk group, age, ethnic group, resident city, marital status, education, and occupation. HIV risk behaviors were categorized as MSM, heterosexual (hetero), injection drug user (IDU), and other/unknown. The following clinical data were analyzed: CD4 + lymphocyte count, HIV-1 RNA viral load, RHI or CHI, and the time between HIV infection diagnosis and ART initiation. Because the standards of ART initiation were updated in 2008 and 2014, the study period was divided into three 3-year phases: 2008–2010, 2011–2013, and 2014–2016. Early and delay ART were defined as initiated ART within 6 months and above 2 years after diagnosis.
Proportional Detection Rate, Cluster Growth Predictor, and Effective Reproductive Number (Re)
To describe the dynamics of a given cluster, three parameters–i.e., PDR, cluster growth predictor, and Re–were calculated as follows. PDR for a given year (j) was calculated as the cumulative number of cases in the cluster sampled up to and including year j, divided by the cumulative number of cases up to and during the last sampling year (i), per observation time between years j and i. PDR ≥ 2 (i.e., a 2-fold increase in size in 1 year) was considered as a significant change (Dennis et al., 2020).
Cluster growth predictor was calculated as previously described (Wertheim et al., 2018) as the number of newly diagnosed individuals in a given year divided by the square root of cluster size at the end of that year. A declining curve indicated that a given cluster had a very low probability of causing an outbreak.
Re of each large cluster (≥ 10 cases) was estimated with the birth–death skyline serial model in BEAST v2.4.2 (Jovanovic et al., 2019; Vasylyeva et al., 2019; Vinken et al., 2019; Dennis et al., 2020). Re represents the average number of secondary infections caused by a typical infected individual when only part of the population is susceptible. The value is often used to describe temporal changes of an epidemic in a population, with Re > 1 and Re < 1 indicating the growth or decline of the epidemic, respectively.
Results were analyzed with standard statistical tests. Categoric data were compared with the chi-squared test or Fisher’s exact test using SPSS v20.0 (SPSS Inc, Chicago, IL, United States). P < 0.05 indicated a statistically significant difference.
The Subtype Profile of Study Population and Molecular Network Characteristics
A total of 2221 individuals who were newly diagnosed with HIV infection between January 2008 and December 2016 were included in the study. Of these cases, 75% (1669/2221) were CRF01_AE, 12% (264/2221) were CRF07_BC, 7% (154/2221) were subtype B, and 6% (134/2221) were other CRFs/unique recombinant forms (URFs) (Supplementary Figures S1, S2). After CRF01_AE, CRF07_BC, and subtype B, the top 5 CRFs were CRF59_01B (n = 16), CRF55_01B (n = 12), CRF08_BC (n = 10), CRF67_01B (n = 9), and CRF33_01B (n = 6). The prevalence of both CRF01_AE and subtype B gradually declined between 2008 and 2016 (from 81.4 to 67.4% and from 10.5 to 6.2%, respectively), while that of CRF07_BC and other CRFs/URFs increased (from 4.7 to 18.1% and from 3.5 to 8.3%, respectively) (Supplementary Figure S2).
Molecular networks were constructed for 2087 sequences of CRF01_AE, CRF07_BC, and subtype B using HIV-TRACE. A total of 788 (37.8%) sequences (81.9% CRF01_AE, 10.3% CRF07_BC, and 7.9% subtype B) were linked to at least one other sequence and formed 168 transmission clusters, including 138 of CRF01_AE, 16 of CRF07_BC, and 14 of subtype B, with cluster size ranging from 2 to 107 sequences. Of the 788 clustered sequences, 89.0%, 92.6%, and 88.7% were MSM nodes for CRF01_AE, CRF07_BC, and subtype B, respectively. Of the 168 transmission clusters, 66.1% (111/168) comprised only MSM nodes; 32.7% (55/168) contained hetero nodes, and the percentage of hetero nodes in a hetero-related cluster ranged from 2.6 to 100%. However, 14 hetero-dominated clusters (i.e., in which hetero nodes accounted for more than half of those in the cluster) were all small (cluster size of 2 or 3). Only two clusters of CRF01_AE were IDU-related. Of the 168 clusters, there were 13 large clusters each comprising at least 10 patients, including nine of CRF01_AE clusters, one of CRF07_BC, and three of subtype B (Figure 1).
Figure 1. Distribution of risk group clusters. (A–C) Risk group clusters in transmission networks of CRF01_AE (A), CRF07_BC (B), and subtype B (C). Blue, orange, and black columns represents clusters containing 100% MSM and hetero- and IDU-related clusters, respectively.
Population Characteristics of Large Clusters and Other Groups
Of the 2087 individuals newly diagnosed with HIV infection, 444 (21.3%) and 344 (16.5%) belonged to small/medium (2–9 cases) and large (≥ 10 cases) clusters, respectively, whereas 1299 (62.2%) were non-clustered. We evaluated factors associated with clustering in these individuals and found that those in large clusters had distinct characteristics from individuals in the other two groups (Table 1), and were more likely to be male (85.8 vs 77.7% in small/medium clusters and 71.6% in non-clustered individuals; P < 0.001) and younger (29.4% of individuals < 25 years of age vs. 20.7% and 18.6%; P < 0.001); have RHI status (33.4% vs 30.9% and 22.8%; P < 0.001); report MSM contact as their main risk behavior (94.5% vs 85.4% and 84.2%; P < 0.001); and have a high CD4 + cell count (18% with ≥ 500 cells/μl vs 14% and 11.6%; P < 0.001). No significant differences in viral load were observed between groups. Other factors that increased the probability of clustering were Han ethnicity, single marital status, and residence in cities in Liaoning other than Shenyang (P < 0.001).
Table 1. Comparison of characteristics of HIV-1 infected individuals by cluster inclusion and cluster size.
Progressive Decline in the Proportion of Large Clusters Over Time
Individuals in large clusters (≥ 10 cases) tended to be diagnosed earlier (2008–2010) than those who were not in a cluster (29.4% vs 13.3%, P < 0.001) (Table 1). Further we evaluated the contribution of large clusters to the local HIV epidemic over time (Figure 2). The proportion of individuals in large clusters gradually declined from 30% in 2008 to 8% in 2016. Importantly, the proportion of individuals with RHI in large clusters also decreased from 66% in 2009 to 20% in 2016 (Figure 3), with the number of RHI cases decreasing from 29 to 6 during that period. In contrast, the percentage of non-clustered individuals increased steadily from 52% in 2008 to 73% in 2016, while no substantial changes were observed in small/medium clusters (18% in 2008 and 19% in 2016).
Figure 2. Composition of newly diagnosed HIV-1–infected cases. Blue, orange, and gray columns represent cases belonging to large clusters (≥ 10 cases), small/medium clusters (2–9 cases), and cases not in a cluster, respectively.
Figure 3. Expanding dynamics of the three groups including 13 local large clusters (≥ 10 cases) from 2008 to 2016. Years are shown along the x-axis and the number of new diagnoses for each cluster or all large clusters (≥ 10 cases) is shown along the y-axis. Orange and blue columns indicate recent and chronic HIV infection, respectively.
Expansion History and Dynamics of Large Clusters in the Period of 2008–2016
To clarify the dynamics of the 13 large clusters (≥ 10 cases), we analyzed the expansion history of them during the period of 2008–2016 (Figure 3). These clusters were roughly divided into historical (2008–2010), middle-phase (2011–2013), and recently active (2014–2016) according to their period of most rapid expansion (defined as > 45% of the final cluster size reached by the end of 2016). For example, the cluster AE-1 within which 48% cases were diagnosed between 2008 and 2010 belong to “Historical Group”. PDR, cluster growth predictor, and Re were retrospectively calculated to evaluate the expansion speed of each large cluster. Although the shape of the curves varied, the three parameters confirmed that transmission of most large clusters declined during the study period (Figures 4A,B).
Figure 4. PDR trends, cluster growth predictor trends, and Re of 13 large clusters (≥ 10 cases). (A) PDR trends (blue line), cluster growth predictor trends (orange line) (A), and Re number (secondary infections per infected person) (B) during the period of 2008–2016. Re was inferred from the birth–death skyline plot. Thick solid lines represent the estimated mean, and 95% highest posterior density credible regions are shown as gray areas.
For the historical group (clusters AE-1, AE-3, and B-2), both newly diagnosed (48%, 51%, and 58%) and RHI (48%, 75%, and 100%) cases were concentrated between 2008 and 2010 (Figures 3, 4). Similarly, both PDR and cluster growth predictor declined, with slight fluctuations. The peak of the Re curves for AE-1, AE-3, and B-2 occurred in 2008, 2009, and 2010, respectively. The Re of all of three clusters declined and remained at 1 for > 5 years, implying that the clusters were historical and have receded in recent years.
The middle-phase group included clusters AE-5, B-1, AE-6, and AE-8, which had more cases diagnosed between 2011 and 2013 than at any other time (59%, 50%, 47%, and 50%, respectively); three of the clusters (B-1, AE-5, and AE-8) had even higher proportions of RHI cases diagnosed during this period (75%, 71.4%, and 42.9%, respectively) (Figures 3, 4). These data were consistent with the trends observed for PDR and cluster growth predictor and demonstrated that all clusters underwent rapid expansion in 2011–2013 and declined thereafter. The Re of clusters B-1 and AE-6 were relatively high before 2010, while AE-5 and AE-8 were expanding during 2012–2014 and 2010–2014, respectively. However, all four clusters showed a stable Re of 1 in recent years.
The recently active group (clusters 07BC-1, AE-2, AE-4, AE-7, AE-9, and B-3) had more newly diagnosed cases (> 45%) in 2014–2016 than at any other time (Figures 3, 4). AE-9 and B-3 were recently emerged clusters that appeared in 2013 and 2014, respectively. The PDR of 5 of the 6 recently active clusters decreased or remained at a constant low level; the exception was cluster B-3, for which PDR increased from 1.25 in 2015 to 2 in 2016. Similarly, the cluster growth curves predictor of AE-2, AE-4, and AE-9 fluctuated but declined toward the end of the study period, while the curves for 07BC-1, AE-7, and B-3 showed an upward trend. According to the birth–death model, cluster 07BC-1 and AE-7 declined and remained at 1 in 2016, as well as cluster AE-2, AE-4 and AE-9. On the contrary, the Re of B-3 increased between 2014 and 2016 and reached 2.676 in 2016. Of the cluster with upward curves (Supplementary Table S2), cases in B-3 tended to be younger (mean age, 26.6 years) and had a higher proportion of local MSM (100%).
Reduced Ongoing Expansion of Large Clusters Coincident With Earlier Initiation of ART
Since the establishment of this long-term cohorts in Shenyang in 2008, the standards of ART for HIV-infected patients in China have incrementally improved. Our data showed that the proportion of cases who initiated ART within 6 months after diagnosis increased from 24% in 2008–2010 to 54% in 2011–2013 and 79% in 2014–2016 (Figure 5). In contrast, the proportion of cases who initiated ART > 2 years after diagnosis decreased from 29% in 2008–2010 to 15% in 2011–2013 and 1% in 2014–2016. Meanwhile, the proportion of patients starting ART between 0.5 and 2 years post diagnosis also decreased (from 17% in 2008–2010 to 6% in 2014–2016), as did the proportion without medical care (from 30% in 2008–2010 to 13% in 2014–2016). It is worth noting that large clusters had a higher percentage of patients who delayed ART or did not receive treatment (38.1%) compared to non-clustered cases (24.7%) and small/medium clusters (27.5%) (P < 0.0001; Table 1), and were mainly concentrated in the period of 2008–2010 (data not shown).
Figure 5. Distribution and composition of the time lag between HIV infection diagnosis and ART initiation. Blue, orange, gray, and yellow columns represent ART initiation within 6 months of diagnosis, ART initiation between 0.5 and 2 years post diagnosis, ART initiation > 2 years post diagnosis, and ART naïve/lost, respectively.
In this study we retrospectively reconstructed the molecular networks of three main HIV-1 subtypes among patients who were newly diagnosed with HIV-1 infection in Shenyang in 2008–2016. Our results show that the expansion of large clusters was progressively controlled, coinciding with and supporting the effectiveness of strengthened early ART and continuous regular interventions.
Threshold Selection in an Area With Multiple HIV-1 Subtypes
The threshold is a key factor in molecular network construction (Hassan et al., 2017). Studies on HIV evolution-which have mainly focused on HIV-1 subtype B- have supported threshold selection and molecular network-guided applications; 1.5% was selected as the optimal threshold based on a rate of evolution of 1% every 10 years for the pol gene (Smith et al., 2009; Hightower et al., 2013). Nearly all molecular network studies on non-B HIV have used the subtype B threshold (Bon et al., 2010; Parczewski et al., 2012; Rose et al., 2017; Fabeni et al., 2019), however, the threshold for subtype B are also appropriate for non-B viral strains have not been fully explored (Han et al., 2020).
Unlike areas with “a single subtype” of HIV- e.g., Western and central Europe and North America (83.3% subtype B); and Southern Africa (98.8% subtype C) (Hemelaar et al., 2019)-there were three subtypes accounting for up to 94% of cases in Shenyang. The prevalence of CRF01_AE and subtype B decreased by around 20% (from 91.9 to 73.6%) between 2008 and 2016, whereas that of CRF07_BC and other CRFs or URFs increased. The few network studies that have been carried out in areas with multiple HIV subtypes and have typically used one single genetic threshold (Hoenigl et al., 2016a; Stecher et al., 2018). Rates of HIV-1 evolution vary across subtypes (Patino-Galindo and Gonzalez-Candelas, 2017; Bbosa et al., 2019a), and no previous studies have focused on HIV evolution or threshold selection for our local epidemic strains. We selected an optimal threshold according to the principle outlined in a previous study (Wertheim et al., 2017) in order to identify the maximum number of clusters in the genetic network. Above the threshold, clusters began to coalesce and the network lost resolution (Supplementary Figure S3). This principle has been used in several recent studies (Chaillon et al., 2019; Ragonnet-Cronin et al., 2019; Zai et al., 2020). We used separate optimal thresholds for CRF01_AE, CRF07_BC, and subtype B, rather than one threshold for multiple subtypes.
Control of Local Large Clusters Has Coincided With Strengthened Early ART
We used HIV molecular networks to identify closely related transmission events; a cluster was formed if the pairwise genetic distance of any two sequences was less than the optimal threshold. A large cluster was treated as a large-scale spreading event in which already highly connected individuals made proportionally more contacts over time (Ragonnet-Cronin et al., 2016). In our study, more individuals with delayed ART or without ART in large clusters were diagnosed in 2008–2010 than in the other two phases of the study, and were likely a source of infection that contributed to the increase in the number of new diagnoses in subsequent years. However, the transmission of most large clusters showed a declining trend at the end of the study, as evidenced by the reductions in the number of new cases in each cluster and in the three parameters of cluster growth (Figures 3, 4). These results suggest that the contribution of large clusters to the local HIV epidemic decreased over time and that these clusters may not be the main driving force of future HIV epidemics.
Moreover, transmissibility is highest in the early stage of HIV infection and plays an important role in ongoing transmission (Brenner et al., 2007; Powers et al., 2011; Marzel et al., 2016). Therefore, clusters with more RHIs are thought to be more active and should be prioritized for intervention. In our study, the rate of RHIs detected by HIV-1 LAg-Avidity EIA or determined from seroconversion records decreased over time in 12 of the 13 large clusters (with cluster B-3 being the exception) (Figure 3), providing further evidence of HIV large transmission clusters declined.
A series of major policies and regulations have been promulgated by the Chinese government to control the spread of HIV. Since the initiation of the National Free Antiretroviral Treatment Program (Zhang et al., 2007; Zhang et al., 2009; Cao et al., 2020) and “Four Frees and One Care” (Sun et al., 2010) program in 2002 and 2003, respectively, HIV testing and access to care have markedly improved, and the free ART program has been rapidly scaled up (Wu et al., 2017; Cao et al., 2020). In 2008, a large prospective cohort of HIV-1–seronegative high-risk individuals were established in Shenyang that received continuous education on HIV and free HIV testing and counseling, and were regularly followed up. Nonetheless, the most significant change during the study period was the improvement of treatment standards and the expansion of treatment coverage, with the criterion for ART initiation changing from CD4 + T cell counts ≤ 350 cells/μl in 2008 to ≤ 500 cells/μl in 2014, with immediate ART now recommended for all patients after a diagnosis of HIV infection. According to our data, early ART (initiated within 6 months after diagnosis) increased from 24 to 79% between 2008 and 2016 (Figure 5), which coincided with decreases in both the ongoing transmission of large clusters and the number of newly diagnosed RHI cases in large clusters, indicating the control of HIV transmission and the effectiveness of preventative treatments (Figures 4, 5). Similarly, a study from Belgium demonstrated that the expansion of an outbreak cluster (subtype F1 outbreak among MSM) was controlled by 2012 (Re < 1), coinciding with a decrease in rates of delayed ART initiation as a result of implementation of the immediate ART initiation policy (Vinken et al., 2019).
Significance of Large Clusters and Other Groups Contributing to the Local Epidemic
In a previous study by the U.S. National HIV Surveillance System, the transmission rates of 11 large clusters were 11 times higher than the national average and were targeted for intervention (France et al., 2018; Oster et al., 2018a,b). Another study conducted in Canada showed that large cluster cases were progressively increasing, contributing to ≥ 40% of ongoing local HIV-1 transmission events (Brenner et al., 2017). Our results also showed that individuals in large clusters were younger, single, had high CD4 counts, were recent infected, were MSM and had delayed/no ART compared to individuals in other groups (Table 1). Thus, individuals in large clusters may have higher risk of further transmission in later years compared to those in small/medium clusters and the non-clustered group (Oster et al., 2018a). We therefore focused on the expansion of large clusters and speculate that their declining contributions over time indicate the progressive control of the local epidemic of main HIV-1 strains in the studied area.
During the study period, the number of newly diagnosed cases in Shenyang was still increasing. HIV infection has a long asymptomatic period, and an infected case can be diagnosed at any time during that period. Thus, the number of newly diagnosed cases was largely affected by testing and did not accurately reflect the severity of the HIV epidemic. RHIs were detected with the HIV-1Lag avidity EIA, and our results showed that the incidence of RHIs was stable during the study period, while that of CHI increased markedly in 2014 and declined slightly thereafter, which we think was due to the start of medical referral from other sites to the study center for treatment. Indeed, the sampling depth of our study increase from 39% in 2013 to 54% in 2014. Moreover, a large number of patients with CHIs who had actually been infected in the past were diagnosed during this period by continuous testing. Given these conditions, new strategies are needed to evaluate the dynamics of HIV epidemics and the effectiveness of interventions.
Non-clustering refers to an absence of links to any cases or clusters (Bbosa et al., 2019b). The non-clustered cases may be: a) the referred patients from other surveillance sites; b) inflow of cases who infected HIV at other cities; c) CHIs with long-term within-host evolution history; and d) recombinant of intra-lineages strains leading the genetic distance gap from other pure sequences. In fact, non-clustering is one of the most overlooked topics in the current literatures; how to manage the growing non-clustered population is an important question that should be addressed.
Perspectives and Recommendations
Our molecular cluster analysis provided evidence that large clusters in Shenyang have been gradually controlled, while highlighting the need to identify individuals who should be prioritized for intervention based on their association with high-risk clusters. To achieve this goal, we recommend the followings. Firstly, both the percentage of RHI cases in recent years and potential for future growth of a cluster should be considered. Cluster B-3 in our study was a rapidly expanding cluster of subtype B that showed increasing PDR and cluster growth predictor curves and Re > 1. It was also a newly emerged cluster with 2 RHIs, which met the criteria of an active cluster for prioritized intervention according to China Center for Disease Control and Prevention guidelines (China-CDC, 2019) Cluster B-3 was shown to have undergone an expansion between 2014 and 2016, and is expected to grow further in coming years. Secondly, interventions should aim to control viral replication in rapidly expanding clusters, particularly those with cases in the early stages of infection. Thirdly, partner services as well as increased testing and intervention options should be provided to persons associated with high-risk clusters. Lastly, additional specific interventions should be considered depending on the characteristics of the high-risk cluster (e.g., cluster B-3), and individuals with the same characteristics should be closely monitored in terms of HIV status and referred for PrEP even if they are HIV-1–seronegative.
In this study, complementary epidemiologic and phylodynamic analyses were used to evaluate the growth tendency of large clusters and identify priorities for intervention. Both PDR and cluster growth predictor are simplified algorithms (Wertheim et al., 2018; Dennis et al., 2020) based on epidemiologic data (number of cases diagnosed each year), and can be easily determined. The curves for both parameters have similar shapes but the PDR curve is more stable, whereas that of cluster growth predictor fluctuates and reveals small changes, especially in larger clusters. For example, in our study a large number of new cases were referred from other surveillance sites to the study center for treatment since 2014; therefore, a peak was observed in the cluster growth predictor curve of many clusters in 2014 (orange lines in Figure 4), whereas no corresponding change was observed in the PDR curves (blue lines in Figure 4). Re, the third phylodynamic parameter that is based on sequence diversity, reflects the efficiency with which an infectious agent is transmitted and is frequently used to model infection dynamics (Novitsky et al., 2015; Ragonnet-Cronin et al., 2018; Jovanovic et al., 2019; Lorenzin et al., 2019; Vasylyeva et al., 2019; Vinken et al., 2019; Zai et al., 2020). We found that Re was accurate and reliable. However, unlike the PDR independent of data volume, Re was unsuitable for small clusters. Thus, the three parameters can be applied in different ways to real-time monitoring of molecular networks for the identification of rapidly expanding clusters: PDR is more suitable for smaller clusters while cluster growth predictor is better for larger ones, and Re can be used for final verification of mid-sized or large clusters.
The assessment of HIV molecular clusters is influenced by sampling depth. In our study, the average sampling rate of 45% may have limited the ability to detect transmission clusters. Additionally, recombination events between diverse HIV strains made it challenging to distinguish clusters, as recombinants are typically removed from analytic datasets (Grabowski et al., 2018). In this study, we analyzed only three major subtypes of HIV-1 (CRF01_AE, CRF07_BC, and subtype B); other CRFs and URFs (6%) were not included because the optimal thresholds could not be determined given the diverse origins and rates of evolution of the different strains.
In summary, the results of this study show that large HIV transmission clusters declined in Shenyang between 2008 and 2016, coinciding with the implementation of early ART and continuous regular interventions and confirming the effectiveness of these strategies. We also demonstrated that molecular network analyses can be used to evaluate the efficacy of interventions in areas of epidemic with a complex HIV profile, which can guide the implementation of targeted interventions.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/genbank/, MT368043-369927; https://www.ncbi.nlm.nih.gov/genbank/, MT336755-336776.
HS and XH conceived and designed the study. BZ, ML, ZW, and YQ performed experimental work, BZ, and HD for the data collection. ML performed molecular and phylodynamic analyses. ML and XH wrote the first draft. All authors read and approved the final manuscript.
This work was supported by the Mega-projects of National Science Research for the 13th 5-Year Plan (2017ZX10201101).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank for Ping Zhong and MA careful review of the manuscript. We also thank Ping Zhong, MA, and BZ for helpful discussions.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.604993/full#supplementary-material
Supplementary Figure 1 | Maximum-likelihood phylogenetic tree of HIV-1 pol sequences from 2221 infected individuals in Shenyang, China. The maximum likelihood (ML) tree was reconstructed using the partial pol sequences under GTR + I + G sites substitution model. The ML tree is rooted with subtype N sequence as outgroup. The support values of the nodes was estimated with the approximate likelihood-ratio (aLRT) SH-like test. CRF01_AE sequences are in red, CRF07_BC sequences are in blue and subtype B are in green. Other CRFs, URFs and their references are in black.
Supplementary Figure 2 | Distribution of HIV-1 subtypes in newly diagnosed infected individuals between January 2008 and December 2016. Columns represent the number of newly diagnosed HIV-infected individuals annually; blue represents CRF01_AE, orange represents CRF07_BC, gray represents subtype B, and yellow represents other CRFs or URFs. The pie chart shows the proportion of each subtype across all individuals.
Supplementary Figure 3 | Number of nodes and transmission clusters as a function of Tamura–Nei 93 distance thresholds. (A) CRF01_AE. (B) CRF07_BC. (C) Subtype B. Blue and orange lines represent the number of clusters and nodes, respectively
Supplementary Table 1 | Characteristics of Molecular Networks at Optimal Genetic Threshold of CRF01_AE, CRF07_BC, and subtype B.
Supplementary Table 2 | Demographic Characteristics of the 13 Large Transmission Clusters.
AIDS-Professional-Group (2011). AIDS professional group, society of infectious diseases chinese medical association: guidelines for diagnosis and treatment of HIV/AIDS in China. Zhonghua Chuan Ran Bing Za Zhi. 29, 629–640.
AIDS-Professional-Group (2015). AIDS professional group, society of infectious diseases, chinese medical association: third edition of the guidelines for diagnosis and treatment of HIV/AIDS. Chin. J. Clin. Infect. Dis. 8, 385–401.
Aldous, J. L., Pond, S. K., Poon, A., Jain, S., Qin, H., Kahn, J. S., et al. (2012). Characterizing HIV transmission networks across the United States. Clin. Infect. Dis. 55, 1135–1143. doi: 10.1093/cid/cis612
Bbosa, N., Ssemwanga, D., Nsubuga, R. N., Salazar-Gonzalez, J. F., Salazar, M. G., Nanyonjo, M., et al. (2019b). Phylogeography of HIV-1 suggests that Ugandan fishing communities are a sink for, not a source of, virus from general populations. Sci. Rep. 9:1051. doi: 10.1038/s41598-018-37458-x
Board, A. R., Oster, A. M., Song, R., Gant, Z., Linley, L., Watson, M., et al. (2020). Geographic distribution of HIV transmission networks in the United States. J. Acquir. Immune Defic. Syndr. 85, e32–e40. doi: 10.1097/QAI.0000000000002448
Bon, I., Ciccozzi, M., Zehender, G., Biagetti, C., Verrucchi, G., Lai, A., et al. (2010). HIV-1 subtype C transmission network: the phylogenetic reconstruction strongly supports the epidemiological data. J. Clin. Virol. 48, 212–214. doi: 10.1016/j.jcv.2010.03.021
Brenner, B. G., Ibanescu, R. I., Hardy, I., Stephens, D., Otis, J., Moodie, E., et al. (2017). Large cluster outbreaks sustain the HIV epidemic among MSM in Quebec. AIDS 31, 707–717. doi: 10.1097/QAD.0000000000001383
Brenner, B. G., Roger, M., Routy, J. P., Moisi, D., Ntemgwa, M., Matte, C., et al. (2007). High rates of forward transmission events after acute/early HIV-1 infection. J. Infect. Dis. 195, 951–959. doi: 10.1086/512088
Cao, W., Hsieh, E., and Li, T. (2020). Optimizing treatment for adults with HIV/AIDS in China: successes over two decades and remaining challenges. Curr. HIV/AIDS Rep. 17, 26–34. doi: 10.1007/s11904-019-00478-x
Castley, A., Sawleshwarkar, S., Varma, R., Herring, B., Thapa, K., Dwyer, D., et al. (2017). A national study of the molecular epidemiology of HIV-1 in Australia 2005-2012. PLoS One 12:e0170601. doi: 10.1371/journal.pone.0170601
Chaillon, A., Delaugerre, C., Brenner, B., Armero, A., Capitant, C., Nere, M. L., et al. (2019). In-depth sampling of high-risk populations to characterize HIV transmission epidemics among Young MSM Using PrEP in france and quebec. Open Forum Infect. Dis. 6:ofz080. doi: 10.1093/ofid/ofz080
Chaillon, A., Essat, A., Frange, P., Smith, D. M., Delaugerre, C., Barin, F., et al. (2017). Spatiotemporal dynamics of HIV-1 transmission in France (1999-2014) and impact of targeted prevention strategies. Retrovirology 14:15. doi: 10.1186/s12977-017-0339-4
Chan, P. A., Hogan, J. W., Huang, A., DeLong, A., Salemi, M., Mayer, K. H., et al. (2015). Phylogenetic investigation of a statewide HIV-1 epidemic reveals ongoing and active transmission networks among men who have sex with men. J. Acquir. Immune Defic. Syndr. 70, 428–435. doi: 10.1097/QAI.0000000000000786
China-CDC (2019). Guidelines for HIV Transmission Network Monitoring and Intervention (Trial version). Available online at: http://www.chinacdc.cn/jkzt/crb/zl/azb/ (accessed December 20, 2019).
China-CDC (2018). AIDS and hepatitis C professional group, society of infectious diseases chinese medical association, China CDC: [Chinese guidelines for diagnosis and treatment of HIV/AIDS]. Zhonghua Nei Ke Za Zhi. 57, 867–884.
Delgado, E., Benito, S., Montero, V., Cuevas, M. T., Fernandez-Garcia, A., Sanchez-Martinez, M., et al. (2019). Diverse large HIV-1 Non-subtype B clusters are spreading among men who have sex with men in Spain. Front. Microbiol. 10:655. doi: 10.3389/fmicb.2019.00655
Dennis, A. M., Hue, S., Billock, R., Levintow, S., Sebastian, J., Miller, W. C., et al. (2020). Human immunodeficiency Virus Type 1 phylodynamics to detect and characterize active transmission Clusters in North Carolina. J. Infect. Dis. 221, 1321–1330. doi: 10.1093/infdis/jiz176
Fabeni, L., Alteri, C., Berno, G., Scutari, R., Orchi, N., De Carli, G., et al. (2019). Characterisation of HIV-1 molecular transmission clusters among newly diagnosed individuals infected with non-B subtypes in Italy. Sex Transm. Infect. 95, 619–625. doi: 10.1136/sextrans-2019-054017
France, A., Panneer, N., Ocfemia, C., Saduvala, N., Campbell, E., Switzer, W., et al. (2018). “Rapidly growing HIV transmission clusters in the Unites States, 2013–2016,” in 2018 Conference on Retroviruses and Opportunistic Infections, San Francisco, CA.
Green, N., Hoenigl, M., Chaillon, A., Anderson, C. M., Kosakovsky Pond, S. L., Smith, D. M., et al. (2017). Partner services in adults with acute and early HIV infection. AIDS 31, 287–293. doi: 10.1097/QAD.0000000000001308
Han, X., Dai, D., Zhao, B., Liu, J., Ding, H., Zhang, M., et al. (2010). Genetic and epidemiologic characterization of HIV-1 infection in liaoning province, China. J. Acquir. Immune Defic. Syndr. 53(Suppl. 1), S27–S33. doi: 10.1097/QAI.0b013e3181c7d5bf
Han, X. X., An, M. H., Zhang, M., Zhao, B., Wu, H., Liang, S., et al. (2013). Identification of 3 distinct HIV-1 founding strains responsible for expanding epidemic among men who have sex with men in 9 Chinese Cities. J. Acquir. Immune Defic. Syndr. 64, 16–24. doi: 10.1097/qai.0b013e3182932210
Hemelaar, J., Elangovan, R., Yun, J., Dickson-Tetteh, L., Fleminger, I., Kirtley, S., et al. (2019). Global and regional molecular epidemiology of HIV-1, 1990–2015: a systematic review, global survey, and trend analysis. Lancet Infect. Dis. 19, 143–155. doi: 10.1016/s1473-3099(18)30647-9
Hightower, G. K., May, S. J., Perez-Santiago, J., Pacold, M. E., Wagner, G. A., Little, S. J., et al. (2013). HIV-1 clade B pol evolution following primary infection. PLoS One 8:e68188. doi: 10.1371/journal.pone.0068188
Hoenigl, M., Chaillon, A., Kessler, H. H., Haas, B., Stelzl, E., Weninger, K., et al. (2016a). Characterization of HIV Transmission in South-East Austria. PLoS One 11:e0151478. doi: 10.1371/journal.pone.0151478
Hoenigl, M., Chaillon, A., Moore, D. J., Morris, S. R., Mehta, S. R., Gianella, S., et al. (2016b). Rapid HIV Viral load suppression in those initiating antiretroviral therapy at first visit after HIV Diagnosis. Sci. Rep. 6:32947. doi: 10.1038/srep32947
Hughes, G. J., Fearnhill, E., Dunn, D., Lycett, S. J., Rambaut, A., Leigh Brown, A. J., et al. (2009). Molecular phylodynamics of the heterosexual HIV epidemic in the United Kingdom. PLoS Pathog. 5:e1000590. doi: 10.1371/journal.ppat.1000590
Jovanovic, L., Siljic, M., Cirkovic, V., Salemovic, D., Pesic-Pavlovic, I., Todorovic, M., et al. (2019). Exploring evolutionary and transmission dynamics of HIV epidemic in serbia: bridging socio-demographic with phylogenetic approach. Front. Microbiol, 10:287. doi: 10.3389/fmicb.2019.00287
Kasaie, P., Pennington, J., Shah, M. S., Berry, S. A., German, D., Flynn, C. P., et al. (2017). The impact of preexposure prophylaxis among men who have sex with men: an individual-based model. J. Acquir. Immune Defic. Syndr. 75, 175–183. doi: 10.1097/QAI.0000000000001354
Kosakovsky Pond, S. L., Weaver, S., Leigh Brown, A. J., and Wertheim, J. O. (2018). HIV-TRACE (TRAnsmission Cluster Engine): a tool for large scale molecular epidemiology of HIV-1 and other rapidly evolving pathogens. Mol. Biol. Evol. 35, 1812–1819. doi: 10.1093/molbev/msy016
Kouyos, R. D., von Wyl, V., Yerly, S., Boni, J., Taffe, P., Shah, C., et al. (2010). Molecular epidemiology reveals long-term changes in HIV type 1 subtype B transmission in Switzerland. J. Infect. Dis. 201, 1488–1497. doi: 10.1086/651951
Leigh Brown, A. J., Lycett, S. J., Weinert, L., Hughes, G. J., Fearnhill, E., Dunn, D. T., et al. (2011). Transmission network parameters estimated from HIV sequences for a nationwide epidemic. J. Infect. Dis. 204, 1463–1469. doi: 10.1093/infdis/jir550
Li, X., Li, W., Zhong, P., Fang, K., Zhu, K., Musa, T. H., et al. (2016). Nationwide trends in molecular epidemiology of HIV-1 in China. AIDS Res. Hum. Retroviruses 32, 851–859. doi: 10.1089/AID.2016.0029
Lorenzin, G., Gargiulo, F., Caruso, A., Caccuri, F., Foca, E., Celotti, A., et al. (2019). Prevalence of Non-B HIV-1 Subtypes in North Italy and analysis of transmission clusters based on sequence data analysis. Microorganisms 8:36. doi: 10.3390/microorganisms8010036
Lubelchek, R. J., Hoehnen, S. C., Hotton, A. L., Kincaid, S. L., Barker, D. E., and French, A. L. (2015). Transmission clustering among newly diagnosed HIV patients in Chicago, 2008 to 2011: using phylogenetics to expand knowledge of regional HIV transmission patterns. J. Acquir. Immune Defic. Syndr. 68, 46–54. doi: 10.1097/QAI.0000000000000404
Marzel, A., Shilaih, M., Yang, W. L., Boni, J., Yerly, S., Klimkait, T., et al. (2016). HIV-1 transmission during recent infection and during treatment interruptions as major drivers of new infections in the Swiss HIV cohort study. Clin. Infect. Dis. 62, 115–122. doi: 10.1093/cid/civ732
Mehta, S. R., Chaillon, A., Gaines, T. L., Gonzalez-Zuniga, P. E., Stockman, J. K., Almanza-Reyes, H., et al. (2018). Impact of public safety policies on human immunodeficiency virus transmission Dynamics in Tijuana Mexico Clin. Infect. Dis. 66, 758–764. doi: 10.1093/cid/cix884
Novitsky, V., Kuhnert, D., Moyo, S., Widenfelt, E., Okui, L., and Essex, M. (2015). Phylodynamic analysis of HIV sub-epidemics in Mochudi. Botswana. Epidemics 13, 44–55. doi: 10.1016/j.epidem.2015.07.002
Oster, A. M., France, A. M., Panneer, N., Banez Ocfemia, M. C., Campbell, E., Dasgupta, S., et al. (2018b). Identifying clusters of recent and rapid HIV transmission through analysis of molecular surveillance data. J. Acquir. Immune Defic. Syndr. 79, 543–550. doi: 10.1097/QAI.0000000000001856
Oster, A. M., Wertheim, J. O., Hernandez, A. L., Ocfemia, M. C., Saduvala, N., and Hall, H. I. (2015). Using molecular HIV surveillance data to understand transmission between subpopulations in the United States. J. Acquir. Immune Defic. Syndr. 70, 444–451. doi: 10.1097/QAI.0000000000000809
Parczewski, M., Leszczyszyn-Pynka, M., Bander, D., Urbanska, A., and Boron-Kaczmarska, A. (2012). HIV-1 subtype D infections among Caucasians from Northwestern Poland–phylogenetic and clinical analysis. PLoS One 7:e31674. doi: 10.1371/journal.pone.0031674
Patino-Galindo, J. A., Torres-Puente, M., Bracho, M. A., Alastrue, I., Juan, A., Navarro, D., et al. (2017). The molecular epidemiology of HIV-1 in the Comunidad Valenciana (Spain): analysis of transmission clusters. Sci. Rep. 7:11584. doi: 10.1038/s41598-017-10286-1
Pineda-Pena, A. C., Theys, K., Stylianou, D. C., Demetriades, I., Program, S. E., Abecasis, A. B., et al. (2018). HIV-1 infection in cyprus, the eastern mediterranean european frontier: a densely sampled transmission dynamics analysis from 1986 to 2012. Sci. Rep. 8:1702. doi: 10.1038/s41598-017-19080-5
Powers, K., Ghani, A., Miller, W., Hoffman, I., Pettifor, A., Kamanga, G., et al. (2011). The role of acute and early HIV infection in the spread of HIV and implications for transmission prevention strategies in Lilongwe. Malawi 378, 256–268. doi: 10.1016/s0140-6736(11)60842-8
Ragonnet-Cronin, M., Hu, Y. W., Morris, S. R., Sheng, Z. J., Poortinga, K., and Wertheim, J. O. (2019). HIV transmission networks among transgender women in Los Angeles County, CA, USA: a phylogenetic analysis of surveillance data. Lancet Hiv 6, E164–E172. doi: 10.1016/S2352-3018(18)30359-X
Ragonnet-Cronin, M., Jackson, C., Bradley-Stewart, A., Aitken, C., McAuley, A., Palmateer, N., et al. (2018). Recent and rapid transmission of HIV Among people who inject drugs in scotland revealed through phylogenetic analysis. J. Infect. Dis. 217, 1875–1882. doi: 10.1093/infdis/jiy130
Ragonnet-Cronin, M., Lycett, S. J., Hodcroft, E. B., Hue, S., Fearnhill, E., Brown, A. E., et al. (2016). Transmission of Non-B HIV subtypes in the United Kingdom is increasingly driven by large non-heterosexual transmission clusters. J. Infect. Dis. 213, 1410–1418. doi: 10.1093/infdis/jiv758
Rhee, S., Magalis, B. R., Hurley, L., Silverberg, M. J., Marcus, J. L., Slome, S., et al. (2019). National and international dimensions of human immunodeficiency Virus-1 sequence clusters in a northern california clinical cohort. Open Forum Infect. Dis. 6:ofz135.
Rose, R., Lamers, S. L., Dollar, J. J., Grabowski, M. K., Hodcroft, E. B., Ragonnet-Cronin, M., et al. (2017). Identifying transmission clusters with cluster picker and HIV-TRACE. AIDS Res. Hum. Retroviruses 33, 211–218. doi: 10.1089/AID.2016.0205
Smith, D. M., May, S. J., Tweeten, S., Drumright, L., Pacold, M. E., Kosakovsky Pond, S. L., et al. (2009). A public health model for the molecular surveillance of HIV transmission in San Diego. California. AIDS 23, 225–232. doi: 10.1097/QAD.0b013e32831d2a81
Stecher, M., Chaillon, A., Eberle, J., Behrens, G. M. N., Eis-Hubinger, A. M., Lehmann, C., et al. (2018). Molecular epidemiology of the HIV epidemic in three german metropolitan regions - cologne/bonn. Munich and Hannover, 1999-2016. Sci. Rep. 8:6799. doi: 10.1038/s41598-018-25004-8
Sun, X., Lu, F., Wu, Z., Poundstone, K., Zeng, G., Xu, P., et al. (2010). Evolution of information-driven HIV/AIDS policies in China. Int. J. Epidemiol. 39(Suppl. 2), ii4–ii13. doi: 10.1093/ije/dyq217
US-CDC (2018). National Center for HIV/AIDS, Viral Hepatitis, STD, and TB Prevention; Division of HIV/AIDS:Prevention Detecting and Responding to HIV Transmission Clusters: A Guide for Health Department. Avaliable at: https://www.cdc.gov/hiv/pdf/funding/announcements/ps18-1802/CDC-HIV-PS18-1802-AttachmentEDetecting-Investigating-and-Responding-to-HIV-TransmissionClusters.pdf (accessed October 17, 2018).
Vasylyeva, T. I., du Plessis, L., Pineda-Pena, A. C., Kuhnert, D., Lemey, P., Vandamme, A. M., et al. (2019). Tracing the impact of Public Health Interventions on HIV-1 transmission in portugal using molecular epidemiology. J. Infect. Dis. 220, 233–243. doi: 10.1093/infdis/jiz085
Vinken, L., Fransen, K., Cuypers, L., Alexiev, I., Balotta, C., Debaisieux, L., et al. (2019). Earlier initiation of antiretroviral treatment coincides with an initial control of the HIV-1 sub-subtype F1 outbreak among men-having-sex-with-men in flanders. Belgium. Front. Microbiol. 10:613. doi: 10.3389/fmicb.2019.00613
Volz, E. M., Le Vu, S., Ratmann, O., Tostevin, A., Dunn, D., Orkin, C., et al. (2018). Molecular epidemiology of HIV-1 Subtype B reveals heterogeneous transmission risk: implications for intervention and control. J. Infect. Dis. 217, 1522–1529. doi: 10.1093/infdis/jiy044
Wertheim, J. O., Kosakovsky Pond, S. L., Forgione, L. A., Mehta, S. R., Murrell, B., Shah, S., et al. (2017). Social and genetic networks of HIV-1 transmission in New York City. PLoS Pathog. 13:e1006000. doi: 10.1371/journal.ppat.1006000
Wertheim, J. O., Leigh Brown, A. J., Hepler, N. L., Mehta, S. R., Richman, D. D., Smith, D. M., et al. (2014). The global transmission network of HIV-1. J. Infect. Dis. 209, 304–313. doi: 10.1093/infdis/jit524
Wertheim, J. O., Murrell, B., Mehta, S. R., Forgione, L. A., Kosakovsky Pond, S. L., Smith, D. M., et al. (2018). Growth of HIV-1 molecular transmission clusters in New York City. J. Infect. Dis. 218, 1943–1953. doi: 10.1093/infdis/jiy431
Wertheim, J. O., Oster, A. M., Hernandez, A. L., Saduvala, N., Banez Ocfemia, M. C., and Hall, H. I. (2016). The international dimension of the U.S. HIV Transmission network and onward transmission of HIV recently imported into the United States. AIDS Res. Hum. Retroviruses 32, 1046–1053. doi: 10.1089/AID.2015.0272
Whiteside, Y. O., Song, R., Wertheim, J. O., and Oster, A. M. (2015). Molecular analysis allows inference into HIV transmission among young men who have sex with men in the United States. AIDS 29, 2517–2522. doi: 10.1097/QAD.0000000000000852
Xu, J., Sonnerborg, A., Gao, L., Wang, P., Bouey, J. Z. H., and Cheng, F. (2020). Delayed treatment for people living with HIV in China, 2004-2016: an analysis of an observational cohort. Int. J. Environ. Res. Public Health 17:1809. doi: 10.3390/ijerph17051809
Xu, J. J., Zhang, M., Brown, K., Reilly, K., Wang, H., Hu, Q., et al. (2010). Syphilis and HIV seroconversion among a 12-month prospective cohort of men who have sex with men in Shenyang. China. Sex Transm. Dis. 37, 432–439. doi: 10.1097/OLQ.0b013e3181d13eed
Zhang, F., Dou, Z., Ma, Y., Zhao, Y., Liu, Z., Bulterys, M., et al. (2009). Five-year outcomes of the China national free antiretroviral treatment program. Ann. Intern. Med. 151, 241–251, W–252. doi: 10.7326/0003-4819-151-4-200908180-00006
Zhang, F., Haberer, J. E., Wang, Y., Zhao, Y., Ma, Y., Zhao, D., et al. (2007). The Chinese free antiretroviral treatment program: challenges and responses. AIDS 21(Suppl. 8), S143–S148. doi: 10.1097/01.aids.0000304710.10036.2b
Zhang, Z., Dai, L., Jiang, Y., Feng, K., Liu, L., Xia, W., et al. (2017). Transmission network characteristics based on env and gag sequences from MSM during acute HIV-1 infection in Beijing. China. Arch. Virol. 162, 3329–3338. doi: 10.1007/s00705-017-3485-z
Zhao, B. (2015). Increase of RT-related transmitted drug resistance in non-CRF01_AE among HIV Type 1–infected men who have sex with men in the 7 cities of China. J. Acquir. Immune Defic. Syndr. 68, 250–255. doi: 10.1097/qai.0000000000000467
Zhao, B., Han, X., Dai, D., Liu, J., Ding, H., Xu, J., et al. (2011). New trends of primary drug resistance among HIV type 1-infected men who have sex with men in Liaoning Province. China. AIDS Res. Hum. Retroviruses 27, 1047–1053. doi: 10.1089/AID.2010.0119
Keywords: HIV-1, transmission cluster, phylodynamics, molecular epidemiology, antiretroviral therapy
Citation: Liu M, Han X, Zhao B, An M, He W, Wang Z, Qiu Y, Ding H and Shang H (2020) Dynamics of HIV-1 Molecular Networks Reveal Effective Control of Large Transmission Clusters in an Area Affected by an Epidemic of Multiple HIV Subtypes. Front. Microbiol. 11:604993. doi: 10.3389/fmicb.2020.604993
Received: 11 September 2020; Accepted: 27 October 2020;
Published: 13 November 2020.
Edited by:Helene Dutartre, UMR 5308 Centre International de Recherche en Infectiologie (CIRI), France
Reviewed by:Bin Su, Capital Medical University, China
Antoinette Van Der Kuyl, University of Amsterdam, Netherlands
Copyright © 2020 Liu, Han, Zhao, An, He, Wang, Qiu, Ding and Shang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hong Shang, email@example.com
†These authors have contributed equally to this work