Skip to main content

Phylogenetic analysis highlights the role of older people in the transmission of HIV-1 in Fuyang, Anhui Province, China



The proportion of older HIV-1 infected people in China has increased rapidly in recent years. Elucidation of the transmission characteristics of this high-risk population subgroup is helpful for the development of tailored interventions.


A phylogenetic analysis was performed that uses available HIV-1 pol sequences amplified with nested RT-PCR from plasma samples of all newly diagnosed participants spanning from October 2017 to September 2018 in Fuyang, Anhui Province. Transmission clusters were identified as two or more sequences that shared a corresponding node with an aLRT-SH value ≥90 in the maximum-likelihood phylogenetic tree and had an overall mean genetic distance of ≤1.5%. A local transmission cluster was defined as a cluster that had more than 80% of its sequences from Fuyang. The role of older people in local HIV-1 transmission was determined using an integration of molecular and demographic data.


Of 362 available sequences, 14 subtypes, and 28 local transmission clusters were identified. It was found that the proportion of older people in the local transmission cluster (69/77, 89.61%) was much higher than that of younger people (46/114, 40.35%) (χ2 test, P < 0.001). In the pretreatment drug resistance analysis, the proportion of sequences with PDRMs in the local transmission cluster was not significantly different between the older people group (57.14%, 4/7) and non-old-aged group (11.11%, 1/9) (Fisher’s exact test, P > 0.05).


By combining phylogenetic analyses with demographic data, more detailed information was provided about the local transmission structure in Fuyang. These findings suggested that older people play an important role in local transmission, and more tailored interventions for this population subgroup are urgently needed.

Peer Review reports


It has been commonly assumed that older people are not a high-risk group for human immunodeficiency viruses (HIV) infection [1]. However, when HIV has become a manageable chronic disease, older people living with HIV, defined as HIV-infected adults who are aged 50 and over, [2] hit a record high of approximately 5.8 million accounting for 15.8% of 36.7 million HIV infections in 2015, according to The Joint United Nations Program on HIV and AIDS (UNAIDS) report [3]. The report noted that older people also remained sexually active and exhibited a moderate to high level of interest toward sex well into their 80s [4,5,6]. Moreover, they were even more likely to have multiple partners and use condoms inconsistently [7]. These facts led to the realization that older people were becoming a high-risk group. Also, according to the Center for Disease Control and Prevention (CDC) annual data, the proportion of older people newly diagnosed with HIV-1 infection increased in recent years in some areas of China [8, 9].

Fuyang, neighboring Henan Province and located in the northwest of the Anhui Province, is an underdeveloped city with an area of 9775 km2 that represents this trend. According to a recent local HIV-1 epidemiologic survey, this city was once defined as an HIV-1 epidemic focus where infections have been primarily caused by commercial blood donations in the early phase of HIV-1 epidemics, [10, 11] and the number of newly reported HIV cases reached to its highest of 506 in 2004 due to the large-scale screening in previous blood donors. After implementing the policy of standardized management of blood stations, the reported HIV/AIDS cases decreased to its lowest of 143 in 2011. Meanwhile, the proportion of HIV infected through sexual contact was beginning to be higher than that of blood donations, and it has been determined that currently the HIV-1 epidemic is primarily driven by sexual transmissions, which accounted for more than 95%. By 2014, a total of 3521 HIV-1 infected cases (486 HIV infections, 1917 AIDS cases and 1118 deaths) had been reported, and the HIV prevalence was 0.031% (2403/7,823,000). Moreover, this survey also highlighted that the proportion of older newly HIV-1 infections rapidly increased from 19.0 to 32.1% in 2009–2014, though the proportion of local older people in Fuyang which had a residential population of 7.6–8.4 million has stabilized at 23% on average and didn’t change significantly in last decade. Elucidation of the transmission characteristics of this high-risk population subgroup is helpful for the development of tailored interventions.

Phylogenetic analysis is a useful tool that can be reliably used to define clusters of closely related sequences reflecting actual transmission, [12] while traditional epidemiological approaches based on behavioral information cannot. When combined with epidemiological and clinical data, more detailed transmission information can be identified that could be of public health relevance [13,14,15]. This method had been successfully used to implement tailored preventive interventions for different high-risk population subgroups, such as commercial sex workers, [16] men who have sex with men (MSM), [17] and injecting drug users (IDUs) [18].

The aim of this study is to reveal the local HIV-1 transmission structure of Fuyang by using HIV-1 pol sequences that have been proven to be sufficiently variable to permit the phylogenetic reconstruction of transmissions [19]. This study particularly focuses on older people with the goal of providing information regarding tailored preventive strategies.


Study population and specimens

All newly HIV-1 diagnoses in Fuyang from October 2017 to September 2018 were enrolled in this study. Blood samples were collected by the CDC of Fuyang for the purpose of baseline drug resistance genotyping. All 369 antiretroviral-naive participants signed written informed consents prior to sample donations, and 340 persons completed standardized questionnaires that included demographic data, while 29 persons denied. This study was reviewed and approved by the ethics committees of the Beijing Institute of Microbiology and Epidemiology and the Anhui Center for Disease Prevention and Control.

Sequence data

Total viral RNA was extracted from plasma isolated from blood samples using the Roche High Pure Viral RNA kit (Roche, REF: 11858882001). A nested reverse transcriptase polymerase chain reaction (RT-PCR) was applied to amplify a pol fragment of the 1062 bp length region (from 2253 to 3314 according to HXB2 calibrator) spanning the protease gene and partial reverse transcriptase gene using the One Step RNA PCR Kit (Takara, RR055A) and the ExTaq Kit (Takara, RR902A) with sets of primers and thermal cycling conditions as described previously [20]. The PCR products were purified and subjected to direct DNA sequencing on an Applied Biosystems 3730 Sequencer.

Phylogenetic analyses

Phylogenetic analysis was used to infer the HIV-1 subtypes. The sequences were then submitted to a quality control program on The Los Alamos National Laboratory (LANL) HIV database ( to check for laboratory contamination. Then, all the available sequences were aligned using reference strains of group M and the circulating recombinant forms (CRFs) in China using the ClustalW multiple alignment program included in the BioEdit software. They were then manually edited. A phylogenetic tree was constructed using the neighbor-joining (NJ) method based on the Kimura 2-parameter model with 1000 bootstrap replicates. For sequences of unknown HIV-1 subtype, the jumping profile hidden Markov model (jpHMM) ( was applied to screen for recombination breakpoints.

For the purpose of transmission cluster analysis, we followed a methodology described previously and included a data set of similar sequences from LANL HIV database by identifying the ten closest matches with BLAST for each of the Fuyang sequences in the study [21]. The final data set containing 1645 sequences (n = 362 from Fuyang and n = 1283 from LANL sequences) was then aligned to find the best fitting nucleotide substitution model using the FindModel tool ( A maximum likelihood (ML) phylogenetic tree was constructed in the PhyML 3.0 online program using the GTR nucleotide substitution model and nearest neighbor interchanges (NNI) to estimate the tree topology by approximate likelihood ratio test Shimodaira-Hasegawa (aLRT-SH) like branch support [22]. In addition, the pairwise genetic distance was computed using the Tamura-Nei model in the MEGA6 software. Transmission clusters were picked out in ClusterPicker1.2 when the relationship among two or more sequences fulfilled the following two stringent conditions: the tree nodes with an aLRT-SH value ≥90 and an overall mean genetic distance less than 0.015 nucleotide substitutions per site [21, 23]. Clusters that had more than 80% sequences from Fuyang were defined as local transmission clusters [21].

Pretreatment drug resistance analysis

For the analysis of drug resistance mutations (DRMs), all available sequences were submitted to the Genotypic Resistance Interpretation Algorithm of the Stanford drug resistance database ( Pretreatment drug resistance mutations (PDRMs) were identified according to the WHO 2009 update of the DRMs list [24].

Statistical analyses of demographic data

The demographic data were all based on a self-reported standardized questionnaire administered at the time of diagnosis. The variables of gender, age, education, marriage, transmission category, and HIV-1 subtype were evaluated for the purpose of depicting the local HIV-1 infection status. To further focus on elderly people, persons aged 50 or older were considered to be in the older people group, and persons younger than 50 years old were conversely classified into the non-old-aged group. All the variables that differed between these two groups were examined. Additionally, the contribution of older people to local HIV-1 transmission was also accessed by comparing the number of older people and non-old-aged people enrolled in local transmission cluster and non-local transmission cluster. All statistical analyses were tested using Pearson’s χ2 test or Fisher’s exact test based on a two-sided α of 0.05 using R version 3.5.1.


Description of the demographic information

After excluding seven people due to the poor quality of their pol sequences, the final study population consisted of 333 people with demographic data available and 29 persons without demographic information.

Overall, the study population fraction of males to females was 4.29:1 (270/63). In reference to the transmission category, as many as 60.66% (202/333) of HIV-1 infected cases acquired the virus by heterosexual contact, 38.44% (128/333) of HIV-1 infected cases were MSM, 0.3% (1/333) were infected by injection drug use or heterosexual contact, 0.3% (1/333) were IDU, and 0.3% (1/333) of cases had an unknown transmission category. The median age of all of the participants was 44 years (Interquartile range: 31–59), of which the older people group accounted for 37.84% (126/333), and the proportion of persons in the non-old-aged group was 62.16% (207/333). The demographic comparison of these two groups is shown in Table 1. There was no significant difference in gender composition, but with regard to educated and married people, older people were more likely to be divorced (29.37% vs. 12.56%) and undereducated (the proportion of whom had a degree under high school: 92.06% vs. 68.08%). A more important fact is that the transmission category also exhibited statistically significant difference (P < 0.001). Older people were more likely to be infected through heterosexual contact (113/126, 89.68%), in contrast, the proportion of MSM (117/207, 56.52%) was higher than whom infected through heterosexual contact (89/207, 43.00%) in non-old aged group. Although there is the possibility of misclassification of MSM as heterosexuals and it has been confirmed in investigations conducted in UK, [25, 26] older people influenced by Chinese traditional culture especially in underdeveloped areas of China are more likely to be heterosexual, and it is also found that some of them are even involved in drug-used heterosexual activities [27].

Table 1 Demographic differences of older people group vs. non-old aged group with HIV-1 infection in Fuyang city

HIV-1 genotype distribution

Of the 362 available sequences, a total of 14 subtypes consisting of three pure subtypes, seven CRFs, and four unique recombinant forms (URFs) were identified. The distribution of HIV-1 subtypes was as follows: CRF07_BC (122/362, 33.70%), B (115/362, 31.77%), CRF01_AE (89/362, 24.59%), CRF08_BC (7/362, 1.93%), C (2/362, 0.55%), G (2/362, 0.55%), CRF55_01B (1/362, 0.28%), CRF67_01B (2/362, 0.55%), CRF68_01B (1/362, 0.28%), CRF85_BC (1/362, 0.28%), and URF (20/362, 5.52%) (Fig. 1b). CRF07_BC, B and CRF01_AE primarily contributed to the local HIV-1 epidemic. With respect to the constituent ratio of these three major circulating subtypes, the older people group was significantly different from the non-old-aged group (χ2 test, P < 0.001). Older people were primarily infected with the B subtype virus (62.70%, 79/126), followed by CRF07_BC (21.43%, 27/126) and CRF01_AE (11.11%, 14/126). However, CRF07_BC (42.51%, 88/207) dominated in the non-old-aged group, followed by CRF01_AE (33.33%, 69/207) and the B subtype (12.08%, 25/207) (Fig. 1c).

Fig. 1

a. Map of China highlighting Fuyang, Anhui Province with the population and area labeled below. b. Pie chart in the lower-left corner representing the HIV-1 subtype distribution of the 362 newly diagnosed people, with available sequences from October 2017 to September 2018 in Fuyang. c. Differences in the constituent ratios of three major circulating subtypes (B subtype, CRF01_AE, and CRF07_BC) between the older people group and the non-old-aged group. “***” represents P < 0.001

Older people played an important role in local transmission

A total of 86 transmission clusters consisting 485 sequences (n = 207 from Fuyang and n = 278 reference sequences) was identified. Of them, 28 were local transmission clusters, with 14 clusters in the B subtype, 9 clusters in CRF07_BC, 2 clusters in CRF01_AE and 3 clusters in other subtypes. The distribution of non-local transmission clusters among subtypes was significantly different (χ2 test, P < 0.001), with 35 clusters in CRF01_AE, 16 clusters in CRF07_BC, 5 clusters in other subtypes and only 2 clusters in B subtype.

For the purpose of explicating the importance role of older people in local transmission, the proportion of older people and non-old-aged people in local transmission cluster was compared. With 207 people from Fuyang identified in transmission cluster, 77 were older people, 114 were non-old-aged people and 16 people with unknown age. It was found that the proportion of older people in local transmission cluster (69/77, 89.61%) was much higher than that of non-old-aged people (46/114, 40.35%) (χ2 test, P < 0.001).

Interestingly, more than half of older people in local transmission cluster (43/69, 62.32%) were identified in three large clusters (11, 12, and 22 sequences). Of these large clusters, one was all older people (13/13) and the other two clusters contained 84.62% (11/13) and 86.36% (19/22) of older people (Fig. 2). The median age of these clusters was homogeneous and is as follows: 73 years (13/13, Interquartile range: 69–77), 68 years (11/13, Interquartile range: 59–70), and 70 years (19/22, Interquartile range: 61.75–75). On the contrary, only one large local transmission cluster containing 10 non-old-aged people was identified. Taken together, these results suggested that older people played an important role in local transmission.

Fig. 2

Transmission clusters in the maximum likelihood (ML) phylogenetic trees. A total of 86 transmission clusters consisting 485 sequences (n = 207 from Fuyang and n = 278 reference sequences) were identified by constructing the ML tree. The results are shown by subtype groups: a. B subtype, b. CRF01_AE, c. CRF07_BC and d. Other subtypes. The branch support (aLRT-SH) above 90% is shown at the corresponding nodes, and corresponding branch colors indicate cluster classifications (Red represents the local transmission cluster; Blue represents the non-local transmission cluster). Square colors indicate the sequences from different population subgroups (Pink: older people; Light blue: non-old-aged people; Light green: reference sequences; Grey: unknown-aged people)

Pretreatment drug resistance

The overall rate of PDRMs in this study population was 4.42% (16/362), which was as low as the percentage reported in a previous survey from the Anhui Province [20]. The rate was 4.35% (9/207) in non-old-aged group, while it reached 5.56% (7/126) among older people (Fisher’s exact test, P > 0.05). As shown in Table 2, all seven older people shared the K103 N mutation that confers resistance to nonnucleoside reverse transcriptase inhibitors (NNRTIs). In contrast, the PDRMs were more complicated among the nine non-old-aged people, and they had eight mutations that conferred resistance to NNRTIs, nucleoside reverse transcriptase inhibitors (NRTIs), and protease inhibitors (PIs). Moreover, the proportion of sequences with PDRMs in local transmission cluster was 57.14% (4/7) in older people group, while it was only 11.11% (1/9) in non-old-aged group (Fisher’s exact test, P > 0.05).

Table 2 All participants with pretreatment drug resistance mutations


This analysis investigated the complicated HIV-1 epidemic in Fuyang that was identified to have 14 subtypes within 1 year. In contrast with the result of a nationwide molecular epidemiological survey indicating CRF01_AE was the most prevalent strain in China after 2007, [28] it was found that CRF07_BC predominated in this area, followed by subtype B that accounted for as much as 31.77% of the strains. It is rare to have such a high proportion of subtype B after the early-mid 1990s when commercial plasma collection in central China resulted in an explosion of subtype B infections among former plasma donors (FPDs) [29]. According to a recent local HIV-1 epidemiologic survey, the epidemic in Fuyang was once affected by FPDs. It is now primarily driven by sexual contact transmission [11]. Therefore, these results suggested that subtype B strains had possibly spread out of the FPDs into the general population through sexual contact. This hypothesis is consistent with the result of our previous study conducted in Henan Province [30].

With further phylogenetic analysis, subtype B sequences also showed a tendency to be in local transmission clusters. Moreover, of the 69 subtype B sequences in local transmission clusters, 72.46% (50/69) sequences were obtained from older people infected through heterosexual contact. When combined with clusters in CRF01_AE and CRF07_BC, a more significant result was found. This was that older people were more likely to be in local transmission cluster. In particular, only one reference sequence from LANL HIV database was included in three large clusters in total. On the contrary, non-old-aged people were more likely to be clustered with reference sequences obtained from other cities in China to form national transmission clusters. This suggested that older people predominated in the local HIV-1 transmission. This discrepancy could be attributed to social and economic factors. Fuyang is a relatively economically underdeveloped city compared to first-line cities, such as Beijing and Guangzhou, where transmission clusters among older people have rarely been identified [31, 32]. Older people in Fuyang prefer to settle down, while non-old-aged people tend to work outside the city. Statistically, about 2.8 million Fuyang residents were migrant workers [11]. The possibility cannot be excluded that some non-old-aged people were seasonal migrants who were infected outside Fuyang and were diagnosed locally. As a consequence, non-old-aged infected people might engage in multiregional virus exchange, while older people are more likely to contribute to local transmission. Future studies that track the floating population may help delineate the local epidemic more precisely.

A sexual behavior survey conducted in Guangxi, China, indicated that older people had little knowledge of HIV-1 preventions and did not perceive themselves to be at risk of HIV-1 [7]. Their lower education, which was found in this study, may be one of the major factors that contribute to the increasingly high proportion of older infected people and active transmission among this group. Moreover, another important factor often overlooked is the lack of emotional care, especially in older people who are divorced. A previous study found that older people were more vulnerable to loneliness and social isolation, [33] and some of them may negatively internalize or experience stigma when they are infected with HIV-1 [34]. As a result, they may engage in high-risk behaviors. This suggests that preventive interventions only are not adequate and more emotional and social support is also needed. Actually, we have not only taken actions to publicize the HIV knowledge by media like TV and radio, but also given lectures on HIV knowledge and humanistic care in communities and rural areas. By now, the treatment target is estimated to be achieved to “75–94-92”.

Though 94% of people with diagnosed HIV infection have received sustained antiretroviral therapy (ART), we should also focus on the increasing PDRMs due to extensive use of ART. In this study, the sequences of the older people had a prevalence of PDRMs (5.56%). Furthermore, all of the older people having PDRMs shared K103 N TDR site, and 51.95% (4/7) of them were in local transmission clusters. This implied that prevention from local transmission of K103 N which has high-level resistance to the free drugs (efavirenz and nevirapine) for HIV infections in China was in need, otherwise it will give rise to a reduction in efficacy of first-line ART. Added to this, older people were often associated with age-related chronic diseases, such as cardiovascular disease, cancer, hypertension, and diabetes. Such circumstances might raise the base-line treatment complexity and increase healthcare costs when coupled with PDRMs. Therefore, tailored preventive interventions for older people are urgently needed.

The generalizability of these results is subjected to certain limitations. For instance, although all newly diagnosed participants in Fuyang within 1 year were recruited, our results only show a representation of central cities in China. There are an increasing number of infected older people cases that have been reported for many other regions. Therefore, the contribution of older people to HIV-1 transmission should draw attention. The phylogenetic analysis applied in this study provides a useful tool for investigations of population-based HIV-1 transmission. However, our results are limited by the disadvantages of cross-sectional samples, which are unable to assess transmission dynamics. Therefore, longitudinal observational studies need to be conducted in the future to obtain more information regarding transmission among older people.

Despite these limitations, this is the first investigation to reveal the transmission structure among older people in China using an integration of molecular, clinical, and demographic data. These findings suggest that older people play an important role in local transmission, and as a previous survey suggested, limited knowledge about HIV-1 results in a lack of motivation to seek testing [35]. Therefore, the contribution of older people may be underestimated to some extent. Importantly, it is urgent that more tailored interventions be made and implemented. Otherwise, the treatment of older HIV-1 infected people will pose a great challenge in the future.


This study highlighted that older people which was commonly assumed not to be a high-risk group for HIV-1 infection played an important role in the contribution of local HIV-1 transmission in Fuyang, Anhui Province, China. As the increasing numbers of older HIV-1 infected people were reported according to UNAIDS data, more attentions to older people should also be paid in other regions. By combining with demographic data, phylogenetic analysis was a useful tool can be applied to provide more tailored interventions which were definitely needed to the older people population subgroup.

Availability of data and materials

The gene sequences were deposited in the GenBank with the accession number: MK458938 – MK459299.



Antiretroviral therapy


the Center for Disease Control and Prevention


Circulating recombinant forms


Drug resistance mutations


Former plasma donors


Human immunodeficiency viruses


Injecting drug users


the jumping profile hidden Markov model


Los Alamos National Laboratory


Maximum likelihood


Men who have sex with men




Nearest neighbor interchanges


Nonnucleoside reverse transcriptase inhibitors


Nucleoside reverse transcriptase inhibitors


Pretreatment drug resistance mutations


Protease inhibitors


Reverse transcriptase polymerase chain reaction


the Joint United Nations Program on HIV and AIDS


Unique recombinant forms


  1. 1.

    Altschuler J, Katz AD, Tynan M. Developing and implementing an HIV/AIDS educational curriculum for older adults. Gerontologist. 2004;44(1):121–6.

    Article  Google Scholar 

  2. 2.

    Emlet CA, Brennan DJ, Brennenstuhl S, Rueda S, Hart TA, Rourke SB. The impact of HIV-related stigma on older and younger adults living with HIV disease: does age matter? AIDS Care. 2015;27(4):520–8.

    Article  Google Scholar 

  3. 3.

    UNAIDS. Get on the Fast-Track: The life-cycle approach to HIV. http://wwwunaidsorg/sites/default/files/media_asset/Get-on-the-Fast-Track_enpdf 2016:92-98.

  4. 4.

    Beckman N, Waern M, Gustafson D, Skoog I. Secular trends in self reported sexual activity and satisfaction in Swedish 70 year olds: cross sectional survey of four populations, 1971–2001. BMJ (Clinical research ed). 2008;337:a279.

    Article  Google Scholar 

  5. 5.

    Karraker A, Delamater J, Schwartz CR. Sexual frequency decline from midlife to later life. J Gerontol Ser B Psychol Sci Soc Sci. 2011;66(4):502–12.

    Article  Google Scholar 

  6. 6.

    Galinsky AM. Sexual touching and difficulties with sexual arousal and orgasm among U.S. older adults. Arch Sex Behav. 2012;41(4):875–90.

    Article  Google Scholar 

  7. 7.

    Zhao Q, Mao Y, Li X, Shen Z, Zhou Y. Age differences in sexual risk behaviors and related factors among people living with HIV in Guangxi, China. AIDS Care. 2018;30(4):523–30.

    Article  Google Scholar 

  8. 8.

    Wang W, Zhang X, Hao XG. HIV/AIDS epidemiological characteristics of the elderly in Quzhou city from 2002 to 2016. Chin J AIDS STD (in Chinese. 2017;23(12):1105–6.

    Google Scholar 

  9. 9.

    Ye ZM, Wang DY, Zhao LN, Su DH, Chen XY, Wang RQ, Zhang HM, Jin X. Epidemiological characteristics of HIV /AIDS patients aged 50 and above in Wenzhou, 2012-2016. Mod Prev Med (in Chinese). 2018, 45(5):776–80.

  10. 10.

    Xu JQ, Wang JJ, Han LF, Xu C, Ruan YH, Xu ZH, Chen X, Liu ZD, Wang J, Su B, et al. Epidemiology, clinical and laboratory characteristics of currently alive HIV-1 infected former blood donors naive to antiretroviral therapy in Anhui Province, China. Chin Med J. 2006;119(23):1941–8.

    Article  Google Scholar 

  11. 11.

    Hu B, Ni J, He L, Wu ZY, Wu RT, Xu J. Epidemic situation of HIV /AIDS in Fuyang city, 1995-2014. Mod Prev Med (in Chinese). 2017;44(5):917–21.

    Google Scholar 

  12. 12.

    Grabowski MK, Redd AD. Molecular tools for studying HIV transmission in sexual networks. Curr Opin HIV AIDS. 2014;9(2):126–33.

    Article  Google Scholar 

  13. 13.

    Morgan E, Oster AM, Townsell S, Peace D, Benbow N, Schneider JA. HIV-1 infection and transmission networks of younger people in Chicago, Illinois, 2005-2011. Public Health Rep (Washington, DC : 1974). 2017;132(1):48–55.

    Article  Google Scholar 

  14. 14.

    Chalmet K, Staelens D, Blot S, Dinakis S, Pelgrom J, Plum J, Vogelaers D, Vandekerckhove L, Verhofstede C. Epidemiological study of phylogenetic transmission clusters in a local HIV-1 epidemic reveals distinct differences between subtype B and non-B infections. BMC Infect Dis. 2010;10(1):262.

    Article  Google Scholar 

  15. 15.

    Yerly S, Junier T, Gayet-Ageron A, Amari EBE, von Wyl V, Günthard HF, Hirschel B, Zdobnov E, Kaiser L, Study atSHC. The impact of transmission clusters on primary drug resistance in newly diagnosed HIV-1 infection. AIDS (London, England). 2009;23(11):1415–23 1410.1097/QAD.1410b1013e32832d32840ad.

    Article  Google Scholar 

  16. 16.

    Dennis AM, Murillo W, de Maria Hernandez F, Guardado ME, Nieto AI, Lorenzana de Rivera I, Eron JJ, Paz-Bailey G. Social network-based recruitment successfully reveals HIV-1 transmission networks among high-risk individuals in El Salvador. J Acquir Immune Defic Syndr (1999). 2013;63(1):135–41.

    Article  Google Scholar 

  17. 17.

    Oster AM, Pieniazek D, Zhang X, Switzer WM, Ziebell RA, Mena LA, Wei X, Johnson KL, Singh SK, Thomas PE, et al. Demographic but not geographic insularity in HIV transmission among young black MSM. AIDS (London, England). 2011;25(17):2157–65 2110.1097/QAD.2150b2013e32834bfde32839.

    Article  Google Scholar 

  18. 18.

    Temereanca A, Oprea C, Wertheim JO, Ianache I, Ceausu E, Cernescu C, Mehta SR, Ruta S. HIV transmission clusters among injecting drug users in Romania. Rom Biotechnol Lett. 2017;22(1):12307–15.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Hué S, Clewley JP, Cane PA, Pillay D. HIV-1 pol gene variation is sufficient for reconstruction of transmissions in the era of antiretroviral therapy. AIDS (London, England). 2004;18(5):719–28.

    Article  Google Scholar 

  20. 20.

    Shen Y, Su B, Wu J, Qin Y, Jin L, Miao L, Liu A, Cheng X. The prevalence of transmitted HIV drug resistance among MSM in Anhui province, China. AIDS Res Ther. 2014;11:19.

    Article  Google Scholar 

  21. 21.

    Arimide DA, Abebe A, Kebede Y, Adugna F, Tilahun T, Kassa D, Assefa Y, Balcha TT, Bjorkman P, Medstrand P. HIV-genetic diversity and drug resistance transmission clusters in Gondar, northern Ethiopia, 2003-2013. PLoS One. 2018;13(10):e0205446.

    Article  Google Scholar 

  22. 22.

    Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.

    CAS  Article  Google Scholar 

  23. 23.

    Hassan AS, Pybus OG, Sanders EJ, Albert J, Esbjornsson J. Defining HIV-1 transmission clusters based on sequence data. AIDS (London, England). 2017;31(9):1211–22.

    Article  Google Scholar 

  24. 24.

    Bennett DE, Camacho RJ, Otelea D, Kuritzkes DR, Fleury H, Kiuchi M, Heneine W, Kantor R, Jordan MR, Schapiro JM, et al. Drug resistance mutations for surveillance of transmitted HIV-1 drug-resistance: 2009 update. PLoS One. 2009;4(3):e4724.

    Article  Google Scholar 

  25. 25.

    Ragonnet-Cronin M, Hue S, Hodcroft EB, Tostevin A, Dunn D, Fawcett T, Pozniak A, Brown AE, Delpech V, Brown AJL. Non-disclosed men who have sex with men in UK HIV transmission networks: phylogenetic analysis of surveillance data. The lancet HIV. 2018;5(6):e309–16.

    Article  Google Scholar 

  26. 26.

    Hué S, Brown AE, Ragonnet-Cronin M, Lycett SJ, Dunn DT, Fearnhill E, Dolling DI, Pozniak A, Pillay D, Delpech VC, et al. Phylogenetic analyses reveal HIV-1 infections between men misclassified as heterosexual transmissions. AIDS (London, England). 2014;28(13):1967–75.

    Article  Google Scholar 

  27. 27.

    Tan HH, Wong ML, Chan RK. An epidemiological and knowledge, attitudes, beliefs and practices study of sexually transmitted infections in older men. Singap Med J. 2006;47(10):886–91.

    CAS  Google Scholar 

  28. 28.

    Li X, Li W, Zhong P, Fang K, Zhu K, Musa TH, Song Y, Du G, Gao R, Guo Y, et al. Nationwide trends in molecular epidemiology of HIV-1 in China. AIDS Res Hum Retrovir. 2016;32(9):851–9.

    CAS  Article  Google Scholar 

  29. 29.

    He X, Xing H, Ruan Y, Hong K, Cheng C, Hu Y, Xin R, Wei J, Feng Y, Hsi JH, et al. A comprehensive mapping of HIV-1 genotypes in various risk groups and regions across China based on a nationwide molecular epidemiologic survey. PLoS One. 2012;7(10):e47289.

    CAS  Article  Google Scholar 

  30. 30.

    Li L, Sun G, Zhong P, Han J, Li T, Jia D, Liu Y, Sun C, Wang Z, Li H, et al. HIV-1 Thai B strain has spread out of former plasma donors into general population through sexual contact in Henan, China. J Med Virol. 2016;88(4):614–21.

    CAS  Article  Google Scholar 

  31. 31.

    Wang X, Wu Y, Mao L, Xia W, Zhang W, Dai L, Mehta SR, Wertheim JO, Dong X, Zhang T, et al. Targeting HIV prevention based on molecular epidemiology among deeply sampled subnetworks of men who have sex with men. Clin Infect Dis. 2015;61(9):1462–8.

    CAS  Article  Google Scholar 

  32. 32.

    Chen S, Cai W, He J, Vidal N, Lai C, Guo W, He H, Chen X, Fu L, Peeters M, et al. Molecular epidemiology of human immunodeficiency virus type 1 in Guangdong province of southern China. PLoS One. 2012;7(11):e48747.

    CAS  Article  Google Scholar 

  33. 33.

    Wang JK, Zhao XD. Family functioning, social support, and quality of life for Chinese empty nest older people with depression. Int J Geriatr Psychiatry. 2012;27(11):1204–6.

    Article  Google Scholar 

  34. 34.

    Lin C, Li L, Ji G, Jie W. Emotional social support and access to care among older people living with HIV in rural China. Int J Geriatr Psychiatry. 2015;30(10):1041–7.

    Article  Google Scholar 

  35. 35.

    Brooks JT, Buchacz K, Gebo KA, Mermin J. HIV infection and older Americans: the public health perspective. Am J Public Health. 2012;102(8):1516–26.

    Article  Google Scholar 

Download references


The authors thank all of the participants, peer workers, and the staff of the Fuyang Center for Disease Prevention and Control for the demographic information collection, sample preparation, and testing.


This study was supported by the NSFC (81773493), the State Key Laboratory of Pathogen and Biosecurity (AMMS), Anhui provincial Program of Prevent medicine & Public Health (Grant#2017jk002), National Grand Program on Key Infectious Disease Control (Grant#2017ZX10201101). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information




LY and LL designed the study; ZY and WJ participated in the process of sequences editing and phylogenetic analyses; ZY, WX, SY, LH and HJ performed statistical analysis; YX, DX, HB and SB collected sample and the demographic data; ZY, WJ, XH, LJ, LY and LL participated in the writing process. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Bin Su or Yongjian Liu or Lin Li.

Ethics declarations

Ethics approval and consent to participate

This study was reviewed and approved by the ethics committees of the Beijing Institute of Microbiology and Epidemiology and the Anhui Center for Disease Prevention and Control. All participants signed written informed consents prior to sample donations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wu, J., Zhang, Y., Shen, Y. et al. Phylogenetic analysis highlights the role of older people in the transmission of HIV-1 in Fuyang, Anhui Province, China. BMC Infect Dis 19, 562 (2019).

Download citation


  • HIV-1
  • Older people
  • Transmission clusters
  • Phylogenetic
  • Epidemiology