Skip to main content
  • Research article
  • Open access
  • Published:

A mathematical model of biomedical interventions for HIV prevention among men who have sex with men in China



The new HIV treatment guidelines in China recommend antiretroviral therapy (ART) for all people living with HIV, but significant gaps in implementation still exist. Pre-exposure prophylaxis (PrEP) can effectively reduce the risk of HIV transmission among men who have sex with men (MSM). This study assessed the epidemiological impact and cost effectiveness of PrEP, enhanced biomedical interventions and their combination among MSM in China.


A deterministic mathematical model was developed and projected over 20 years to assess the impact of the PrEP, biomedical interventions and their combinations. Incidence and prevalence of HIV were measured, and cost-effectiveness was assessed using incremental cost (international dollars, Int.$) per quality-adjusted life year (QALY) gained.


A total of 0.78 million new HIV infections were estimated to occur over the next 20 years if no additional interventions are implemented among MSM. The PrEP-only strategy covering 25–75% of HIV-negative high-risk MSM can prevent 0.09–0.20 million (12.1–25.7%) new infections, at a cost of 17,277–18,452 Int.$/QALY. The optimal cost-effectiveness path is from test-and-treat to the combination strategy of test-and-treat and PrEP. Some strategies could almost eliminate new HIV infections over the next 20 years.


PrEP, test-and-treat, and their combinations among MSM are effective and cost-effective relative to current policy. PrEP is an important and cost-effective addition to current policy in China.

Peer Review reports


HIV/AIDS continues to be a serious public health problem in China, especially among key risk populations. It is estimated that 780,000 people were living with HIV (PLWH) in China in 2011, of whom 32.5% were men who have sex with men (MSM) [1]. HIV prevalence among MSM in China increased sharply from 5.5% in 2009 to 8.0% in 2015 [2], and rapid upward trends in HIV incidence were also found among this population, increasing from 0.39 (0.21–0.56) per 100 person-years in 2000 to 0.98 (0.70–1.25) in 2010 [3]. The percentage of reported cases attributed to MSM also increased from 2.5% in 2006 to 25.8% in 2014 [4] and MSM are now the fastest growing group of PLWH in China.

The Chinese government has increased policy attention and scaled up responses to the HIV/AIDS epidemic [5], including the recent revision of national treatment criteria, which now recommend initiation of ART for all PLWH, regardless of World Health Organization (WHO) clinical stage and at any CD4 cell count [6]. However, in real-world settings, significant gaps still exist [7] to achieve the Joint United Nations Programme on HIV/AIDS (UNAIDS) 90–90-90 targets by 2020. By the end of 2015, China had not reached the 90–90-90 targets, with only 68% of PLWH diagnosed, 67% of diagnosed PLWH in treatment and 91% of people in treatment attaining viral load targets [7]. For example, the lifetime- and annual- HIV testing rates among MSM in China are only 47 and 38% respectively [8]. It is estimated that 61.6–87% of all HIV cases among MSM remain undiagnosed [9, 10]. Around 56% [11] of diagnosed HIV-positive MSM received ART. Because of these gaps in intervention effectiveness, new strategies need to be considered, and recently attention has turned to pre-exposure prophylaxis (PrEP) – the use of tenofovir and emtricitabine among HIV-negative individuals. Large scale randomized controlled trials have consistently shown that both daily PrEP [12] and on-demand PrEP [13] can effectively reduce the risk of HIV transmission among high-risk MSM by 44–92% depending on adherence. New WHO guidelines [14] recommend offering PrEP to people at substantial risk of HIV infection, but this option has not yet been made available in China.

It is important to assess the potential impact of PrEP strategies and their combination with other biomedical interventions in the Chinese context. In this study, we used mathematical modeling to investigate the epidemiological impact and cost-effectiveness of various biomedical interventions and their combinations on HIV transmission over the next 20 years in China. This study assesses the benefits both of full implementation of current policies and the timely introduction of novel policies, and makes recommendations for future HIV policy responses in China.


We developed a deterministic compartmental mathematical model that divides the population at risk into five groups based on HIV serostatus, disease stage and CD4 count, and further divides the population into three groups based on serostatus and treatment state. This model builds on previous models in China [15], the USA [16] and sub-Saharan Africa [17], structures disease stages to reflect WHO treatment recommendations [18], and incorporates a stage for acute infection. A separate “treatment” compartment for HIV-negative subjects allows the incorporation of PrEP into the model. The compartmental structure of the model is shown in Fig. 1, with mathematical symbols removed for simplicity. Arrows between boxes indicate population flows, with arrows not terminating in a compartment indicating mortality. Flow into and out of each compartment is described by an ordinary differential equation (ODE), with each arrow representing a single flow process contributing to a single ODE, and described by a single rate parameter. The rate of flow from compartments 1–3 (uninfected) to compartments 4 or 5 (acute infection) is determined dynamically by a force of infection term composed from the infectiousness of any HIV-positive sexual contact, rates of sexual partnership, and probability that any single partnership is infectious. The full system of ODEs and detailed formulae for the forces of infection are given in Additional file 1 Appendix Section 2.2.

Fig. 1
figure 1

HIV transmission compartmental model structure

The compartment structure shown in Fig. 1 was replicated for high-risk and low-risk MSM, and the total population of MSM was divided between high-risk and low-risk at a ratio of 1:4 (20% of the population was high-risk) in a total population of 3,625,000. This division was based on annual sexual partnerships, with high-risk men having 15 per year and low-risk men having 2.6, for a balanced average of 5 partnerships per year [19]. All population estimates and parameters were drawn from published literature and are summarized in supplementary Table S2. In this model without ART, the asymptomatic disease stage (CD4 > 200) will last for about 8 years [20]. International randomized controlled trials (RCTs) [21, 22] demonstrated a 93–100% reduction in HIV transmission risk for ART. ART effectiveness was assumed to be 90% to reflect the rate of viral suppression among those on ART in China [7]. PrEP effectiveness was assumed to be 60% to reflect the adherence rates from demonstration PrEP studies in China [23, 24]. Some assumptions were also made about reductions in sexual partnerships at advanced disease stages and upon identification of HIV-positive serostatus (Table S2).

The system of ODEs developed to describe Fig. 1 was solved numerically in Matlab 2016 using a difference equation method with a time step of one month. The model was projected forward 20 years to 2037 and data on HIV incidence, HIV prevalence and HIV mortality extracted. The model was calibrated against HIV prevalence data from 2005 to 2015, using a simple search procedure to identify 400 parameter sets based on the deviance criterion [25]. The 400 parameter sets that closest calibrated to observed data were used as uncertainty ranges, and epidemiological outcomes were estimated as the mean of the 400 uncertainty runs, weighted by the inverse deviance of each run. Formulae for epidemiological outcomes and heuristics for model simulation and calibration are described in detail in Additional file 1 Appendix sections 2.4–2.6.

Modeled policy scenarios

We projected a base-case model that assumed current Chinese HIV treatment guidelines were followed for the next 20 years with no change in testing uptake and treatment entry rates. Against this base case we modeled several alternative scenarios and all of their combinations:

  1. 1.

    Scenario 1: Test-and-treat strategy that is fully compliant with the WHO 90–90-90 recommendations, effectively representing the current HIV treatment guidelines implemented at much greater effectiveness than the current standard. This scenario assumed annual testing rates of 90% for all MSM, with an ART utilization rate of 90% for all diagnosed PLWH, and 90% of all PLWH receiving ART achieving viral suppression (i.e., 90% ART effectiveness).

  2. 2.

    Scenario 2–4: PrEP for high-risk MSM with coverage of 25, 50% or 75% of the entire population of identified HIV-negative high-risk MSM.

  3. 3.

    Scenario 5–7: PrEP for high-risk MSM (with the same three levels of PrEP coverage) combined with expanded annual VCT of 90% MSM.

  4. 4.

    Scenario 8–10: PrEP for high-risk MSM (with the same three levels of PrEP coverage) combined with the test-and-treat strategy.

PrEP was modeled separately with and without widespread VCT because current levels of HIV testing amongst MSM are low, and PrEP implemented against this low backdrop of VCT may have limited effectiveness. Test-and-treat to the 90–90-90 standard was modeled as a separate policy scenario to the base case to reflect low current testing and treatment entry rates, especially in rural areas with limited ability to finance the required interventions [7]. We defined elimination of HIV as HIV incidence of less than 1 per 1000 person-years [17].

Health economic analysis

We estimated the cost of all policy scenarios and conducted a full health economic analysis. Policy scenario costs were estimated on the basis of known costs of testing, ART treatment and average cost of treating AIDS and associated opportunistic infections. PrEP costs were assumed on the basis of the domestic price of Truvada. Costs were only estimated for that proportion of the population assumed to be covered by a given intervention. We compared the estimated cost of each policy with the HIV/AIDS-related healthcare costs averted by reduced infections or slowed disease progression, relative to the base case only. Cost-effectiveness analysis was conducted following standard protocols [26], calculating the average cost-effectiveness ratio (CER) for each policy scenario relative to the base case, and then calculating the incremental cost-effectiveness ratio (ICER) of each scenario relative to the next most cost-effective scenario. All costs were reported in 2017 International Dollars (Int.$), calculated using the International Monetary Fund gross domestic product (GDP) deflator and implied purchasing power parity conversion rates. All cost and health outcomes (e.g., quality-adjusted life year, QALY) were calculated with an annual discount rate of 3%. All cost data and quality of life multipliers were drawn from published literature and are summarized in supplementary Table S2.

Sensitivity analysis

We conducted two scenario-based sensitivity analyses. We tested the impact of the PrEP and combined scenarios when PrEP adherence was 30, 60% or 90%, to explore the possible effectiveness with various adherence rates. Secondly, we tested the impact of the PrEP and combined scenarios after adjusting for 100% risk compensation, in which all high-risk MSM in the PrEP program completely ceased all forms of safer sex behavior, to test the potential impact on a PrEP program of risk compensation. Full details and results of these sensitivity analyses are given in Additional file 1: Appendix section 3.


Baseline HIV epidemic projections

The model estimates of the HIV prevalence and incidence rate among Chinese MSM in the past eleven years showed good fit with the national estimates and surveillance data [2] (Fig. 2), and all selected models were not statistically significantly different from the calibration data on a goodness-of-fit test based on their calibration deviance statistic. Without any additional intervention, HIV prevalence will increase to 11.2% in 2037 and there will be 0.78 million cumulative new infections among MSM over the next 20 years (Table 1). In this scenario the HIV incidence rate will be 0.73 per 100 person-years in 2037. Among all scenarios, only combination scenarios (scenarios 9–10) will achieve HIV elimination: the combination strategy of test-and-treat and PrEP (50% of high-risk MSM), by 2037; and the combination strategy of test-and-treat and PrEP (75% of high-risk MSM), by 2033.

Fig. 2
figure 2

Estimated HIV prevalence and incidence between 2005 and 2037 in the status quo and the ten intervention scenarios. Left) HIV prevalence; Right) HIV incidence

Table 1 Benefits, costs and savings in healthcare costs of the individual and combined interventions, 2018~2037

Impact of biomedical interventions

The costs and health benefits (QALYs) for all intervention scenarios are shown in Table 1 and Fig. 3. Among all scenarios, the most cost-effective strategy was test-and-treat, which could prevent 57.9% (0.45 million) of the total new infections estimated under the status quo scenario, with an average CER of 1754 Int.$/QALY.

Fig. 3
figure 3

Cost-effectiveness frontier for all strategies. PrEP, pre-exposure prophylaxis; QALY, quality-adjusted life years. Zone 1 (areas above the red dashed line): Interventions located in zone 1 cost more than three times per capita gross domestic product per QALY gained, and are defined as not cost effective. Zone 2 (areas between the red dashed line and the blue dashed line): Interventions located in zone 2 cost between one and three times per capita gross domestic product per QALY gained, and are defined as cost effective [40, 41]. Zone 3 (areas below the blue dashed line): Interventions located in zone 3 cost less than per capita gross domestic product per QALY gained, and are defined as very cost effective [40, 41]. The black dashed line represents the optimal cost-effectiveness path.

With PrEP effectiveness of 60%, the PrEP-only strategies for high-risk MSM (scenario 2–4) can prevent 0.09–0.20 million infections, 12.1–25.7% of the total new infections, yielding an average CER of 17,277–18,452 Int.$/QALY.

Compared to the base-case scenario, the combination strategies of PrEP and expanded VCT (scenario 5–7) could prevent 0.24–0.36 million new infections, 30.8–47.0% of the total new infections, at a cost of 13,835–18,110 Int.$ per QALY gained. With a combination of test-and-treat and PrEP for 75% of high-risk MSM (scenario 10), 0.58 million new infections (i.e., 75.2% of the total new infections) could be prevented, at a cost of 12,218 Int.$ per QALY gained. When resources are available, the optimal cost-effectiveness path is from test-and-treat to the combination strategy of test-and-treat and PrEP (25% of high-risk MSM), with an ICER of 72,545 Int.$/QALY; followed by the same combination strategy of test-and-treat and PrEP, but with higher PrEP coverage, costing 65,428–73,500 Int.$ per QALY gained depending on coverage level.

Table 1 also shows the saving in healthcare cost for all intervention scenarios compared to the base case. The test-and-treat strategy was the most cost-saving strategy, with a much greater saving in healthcare cost than the incremental costs. Among all combination strategies, the savings in healthcare costs can cover 63.7–100% of the incremental cost, and they can avert 0.08–0.13 million more infections compared to test-and-treat only strategy.

Sensitivity analysis

The sensitivity ranges of the epidemiological impact and cost-effectiveness are also shown in Table 1. The test and treat strategy would avert 46.4–66.5% of the total new HIV infections, whereas PrEP covering 25%/50%/75% of high-risk MSM would avert 8.7–16.7%, 15.1–26.9%, and 19.6–33.3% of the total new infections, respectively. With a combination strategy of test-and-treat and PrEP (75% of high-risk MSM), 69.7–80.5% of the total new infections over 20 years would be prevented. All ten strategies (scenario 1–10) remain cost-effective compared to the base case under all ranges of sensitivity values (Table 1), costing less than three times per capita gross domestic product per QALY gained.

We also tested the robustness of our findings to changes in the assumption about PrEP effectiveness, which depends on PrEP adherence (Additional file 1 Figure S3 and Table S5 in the supplement). With a decrease in PrEP effectiveness from 60 to 30%, under the PrEP 25%/50%/75% high-risk MSM strategies (scenarios 2–4), the total number of new infections prevented will drop from 0.09–0.20 million to 0.05–0.12 million (12.1–25.7% vs. 7.0–14.9%, additional file 1 Figure S3 and Table S5); the cost per QALY gained will increase from 17,277–18,452 Int.$ to 29,902–30,414 Int.$ (Table S5 in the supplement). The results with an increase in PrEP effectiveness from 60 to 90% are described in the supplement.

In the analysis of risk compensation, it was assumed that all PrEP users in scenario 2–10 completely stopped using condoms with their sex partners. With an assumed PrEP effectiveness of 60%, the number of infections prevented will reduce by around 22.2–31.2% in scenarios 2–4 with risk compensation compared to the scenarios without risk compensation (Table S6 (c)), suggesting that risk compensation could be a significant hindrance to the impact of this strategy in situations where PrEP effectiveness is not ideal.

Finally, we modeled the impact of all interventions over a 30-year time horizon, shown in section 3.2 of the supplementary materials.


This study used a compartmental mathematical model to estimate the future trends of the HIV epidemic among MSM in China and to evaluate the cost-effectiveness of individual and combined biomedical interventions. The base case scenario estimated that without scale-up of current treatment strategies and implementation of new interventions the HIV epidemic will not be well controlled in MSM. The provision of PrEP to 25% of high-risk MSM in combination with a fully implemented test-and-treat strategy could prevent more than 0.5 million new HIV infections by 2037. With the combined strategy of test-and-treat and PrEP among high-risk MSM it is possible to achieve HIV elimination by 2033. All individual and combined biomedical interventions were shown to be cost-effective, and two strategies were also shown to be cost-saving.

UNAIDS proposed the ambitious 90–90-90 targets to end the HIV epidemic by 2030. However, in this study, we found that the 90–90-90 strategy alone may not be sufficient to eliminate HIV transmission. Under this strategy, HIV incidence at 2037 is estimated to be 2.4 per 1000 person-years, which is much higher than the incidence threshold for HIV elimination. This result confirms findings from other studies [15, 27, 28] and shows that novel HIV prevention strategies are greatly warranted among MSM in China.

This study also confirms other findings about the effectiveness and cost-effectiveness of PrEP strategies among MSM [27, 29,30,31], extends these findings to a middle-income setting and compares them against test-and-treat strategies. This study also confirms a recent finding about the epidemiological impact and effectiveness of PrEP strategies among Chinese MSM [32]. In this study, we further investigated the conditions for HIV elimination and tested the impact of PrEP adherence and risk compensation on the main findings. Both of these studies show the importance of scaling up existing strategies and implementing novel intervention strategies in China. Our effectiveness and cost-effectiveness estimates may be underestimated because: 1) the PrEP strategy is a comprehensive programme, rather than solely biomedical. Although not modeled in this study, the PrEP programme also includes components of regular HIV testing, behavioral counselling, motivational interviewing and risk reduction, which could lead to further gains in health benefits; 2) in this study, we only analyzed a daily PrEP strategy, but an on-demand strategy would reduce the average number of days that the participants need to take PrEP, lowering the cost of PrEP interventions and significantly improving their acceptability; and 3) the patents for PrEP drugs expired in December 2017 [30], leading to potential future price reductions that could not be incorporated into the cost estimates for this study.

Previously published clinical trials [12, 33] have shown that the protective effects of PrEP relied heavily on drug adherence, which will likely be lower in actual implementation than in trials. In this study we assumed a PrEP effectiveness of 60%, which reflected the adherence rates in PrEP demonstration studies in China [23]. In the scenario-based sensitivity analysis, we tested possible effectiveness of both increased (90%) and decreased (30%) adherence rates. Even under the worst adherence scenario of 30% PrEP effectiveness, the PrEP only strategies can still prevent 7.0–14.9% of HIV infections over a 20 year time horizon. This shows that, while efforts to improve adherence will be important in practical implementation of PrEP in China, failure to achieve trial-level adherence will not fundamentally undermine the effectiveness of the policy itself.

Another concern about PrEP is risk compensation. A recent meta-analysis found evidence of an increase in condomless sex and STI prevalence among PrEP users [34]. One study found that 35% of MSM in China were more likely to have unprotected sex while taking PrEP [35]. Our risk-compensation analyses also found that when PrEP adherence is low, its effectiveness is sensitive to risk compensation. However, with good PrEP adherence (90%), even a 100% risk compensation in condom use would not affect the effectiveness of the intervention for HIV prevention. Therefore, if significant levels of risk compensation are expected amongst the target population, intensive efforts at improving adherence will be needed to realize the full benefits of the PrEP policy. Furthermore, we did not model multiple risk-compensation situations, in which MSM receiving PrEP increase both unsafe sex and number of partners. An appropriately designed and well-implemented demonstration study is needed in China to fully answer these questions, but the possibility of risk compensation should not be taken as a reason not to consider implementing PrEP in China, as our modeling shows that only the most extreme forms of risk compensation will undermine the strategy, and even then it will likely still realize some benefits. However, it is always important to monitor risk compensation among users in the context of PrEP scale-up to assess the impact of PrEP on the sexual health (including HIV and other STI) of users and to inform preventive strategies.

Our study has several limitations. First, heterosexual transmission was not discussed in our model, although 17–35% of Chinese MSM are married [36]. Second, we only targeted up to 75% of high-risk MSM for PrEP. Targeting a larger percentage or expanding to low-risk MSM would lead to a higher impact on the HIV epidemic but a reduced cost-effectiveness. We chose the 75% coverage strategy because awareness of PrEP is low among Chinese MSM [37], and a targeted program aimed at high-risk MSM could be more achievable and feasible. Third, we also did not include false-positive results, resistance to ART, and adverse events of PrEP in our model, but unless extreme negative events happen during the intervention period, their effects on the impact of interventions are not likely to be substantial. Fourth, three of our PrEP scenarios (scenarios 5–7) include the assumption of high rates of voluntary counselling and testing, that might be infeasible in the current environment in China. However, to account for this we also considered three basic scenarios (scenarios 2–4) where PrEP was offered to high-risk MSM, but these MSM did not increase their testing rates to high levels. Even in these more realistic scenarios, we found PrEP was effective. Finally, our study used commonly-accepted cost-effectiveness thresholds of 1-3x per capita GDP for identifying which strategies would be cost-effective, but recent research has shown that these thresholds may not be appropriate for low- and middle-income countries, where use of these thresholds may lead to adoption of lower priority interventions [38]. Using a proposed alternative threshold for China of approximately 50% of GDP [39] would potentially change our conclusion regarding the cost-effectiveness of many of the interventions we considered, but even under this much stricter experimental alternative cost-effectiveness threshold, the test-and-treat strategy, and test-and-treat with 25% PrEP, remain cost-effective. This shows the enormous public health potential of improved HIV interventions in China.


Our findings suggest that both test-and-treat and PrEP can be effective tools for HIV prevention among MSM, and are predicted to be cost-effective even with poor PrEP effectiveness. HIV elimination is achievable in the next 20 years with a combination of these biomedical interventions. PrEP for HIV prevention should be considered for implementation among high-risk MSM in China.



Antiretroviral therapy


Cost-effectiveness ratio


Gross domestic product


Human Immunodeficiency Virus


Incremental cost-effectiveness ratio


Men who have sex with men


Ordinary differential equation


People living with HIV


Pre-exposure prophylaxis


Quality-adjusted life year


Randomized controlled trials


Joint United Nations Programme on HIV/AIDS


United States of America


World Health Organization


  1. Ministry of Health of China. Joint Unitied nations Programme on HIV/AIDS. In: World health organization: 2011 estimates for the HIV/AIDS epidemic in China; 2011.

    Google Scholar 

  2. Tang S, Tang W, Meyers K, Chan P, Chen Z, Tucker JD. HIV epidemiology and responses among men who have sex with men and transgender individuals in China: a scoping review. BMC Infect Dis. 2016;16(1):588.

    Article  Google Scholar 

  3. Zhang L, Chow EPF, Jing J, Zhuang X, Li X, He M, Sun H, Li X, Gorgens M, Wilson D, et al. HIV prevalence in China: integration of surveillance data and a systematic review. Lancet Infect Dis. 2013;13(11):955–63.

    Article  Google Scholar 

  4. National Health and Family Planning Commission of China: 2015 China AIDS response progress report. 2015.

    Google Scholar 

  5. Wu Z, Sullivan SG, Wang Y, Rotheram-Borus MJ, Detels R. Evolution of China's response to HIV/AIDS. Lancet (London, England). 2007;369(9562):679–90.

    Article  Google Scholar 

  6. Pintye J, Baeten JM, Celum C, Mugo N, Ngure K, Were E, Bukusi E, John-Stewart G, Heffron RA. Maternal tenofovir disoproxil fumarate use during pregnancy is not associated with adverse perinatal outcomes among HIV-infected east African women: a prospective study. J Infect Dis. 2017;216(12):1561-568.

    Article  Google Scholar 

  7. Wu Z. The progress and challenges of promoting HIV/AIDS 90–90-90 strategies in China. The Chinese Journal of Disease Control and Prevention. 2016;20(12):1187–9.

    Google Scholar 

  8. Zou H, Hu N, Xin Q, Beck J. HIV testing among men who have sex with men in China: a systematic review and meta-analysis. AIDS Behav. 2012;16(7):1717–28.

    Article  Google Scholar 

  9. Chow EP, Wilson DP, Zhang L. The next era of HIV in China: rapidly spreading epidemics among men who have sex with men. J Acquir Immune Defic Syndr. 2010;55(4):e32–3 author reply e33–34.

    Article  Google Scholar 

  10. Hsieh YH, Ruan Y, Chen CW, Shi W, Li D, Luo F, Shao Y. HIV prevalence and underreporting of men who have sex with men in Beijing. Int J STD AIDS. 2012;23(8):606–7.

    Article  Google Scholar 

  11. Hu X, Chen F, Ding F, Lin X, Wang X, He H, Huang W, Yang Q, Zhang H. Coverage of HIV related follow-up intervention and antiretroviral treatment and its correlation among HIV-positive men who have sex with men of 3 cities in China. Zhonghua yu fang yi xue za zhi [Chinese journal of preventive medicine]. 2015;49(11):945–9.

    Google Scholar 

  12. Grant RM, Lama JR, Anderson PL, McMahan V, Liu AY, Vargas L, Goicochea P, Casapia M, Guanira-Carranza JV, Ramirez-Cardich ME, et al. Preexposure chemoprophylaxis for HIV prevention in men who have sex with men. N Engl J Med. 2010;363(27):2587–99.

    Article  CAS  Google Scholar 

  13. Molina JM, Capitant C, Spire B, Pialoux G, Cotte L, Charreau I, Tremblay C, Le Gall JM, Cua E, Pasquet A, et al. On-demand Preexposure prophylaxis in men at high risk for HIV-1 infection. N Engl J Med. 2015;373(23):2237–46.

    Article  CAS  Google Scholar 

  14. World Health Organization. Guideline on when to start antiretroviral therapy and on pre-exposure prophylaxis for HIV. Geneva, Switzerland; 2015.

  15. Li J, Gilmour S, Zhang H, Koyanagi A, Shibuya K. The epidemiological impact and cost-effectiveness of HIV testing, antiretroviral treatment and harm reduction programs. AIDS. 2012;26(16):2069–78.

    Article  Google Scholar 

  16. Long EF, Brandeau ML, Owens DK. The cost-effectiveness and population outcomes of expanded HIV screening and antiretroviral treatment in the United States. Ann Intern Med. 2010;153(12):778–89.

    Article  Google Scholar 

  17. Granich RM, Gilks CF, Dye C, De Cock KM, Williams BG. Universal voluntary HIV testing with immediate antiretroviral therapy as a strategy for elimination of HIV transmission: a mathematical model. Lancet (London, England). 2009;373(9657):48–57.

    Article  Google Scholar 

  18. Organization WH. Consolidated Guidelines on the Use of Antiretroviral Drugs for Treating and Preventing HIV Infection. Geneva: Recommendations for a Public Health Approach; 2013.

    Google Scholar 

  19. Zhang L, Fung Chow EP, Wilson DP. Men who have sex with men in China have relatively low numbers of sexual partners. Infectious disease reports. 2011;3(1):e10.

    Article  Google Scholar 

  20. Lodi S, Phillips A, Touloumi G, Geskus R, Meyer L, Thiebaut R, Pantazis N, Amo JD, Johnson AM, Babiker A, et al. Time from human immunodeficiency virus seroconversion to reaching CD4+ cell count thresholds <200, <350, and <500 cells/mm(3): assessment of need following changes in treatment guidelines. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America. 2011;53(8):817–25.

    Article  CAS  Google Scholar 

  21. Cohen MS, Chen YQ, McCauley M, Gamble T, Hosseinipour MC, Kumarasamy N, Hakim JG, Kumwenda J, Grinsztejn B, Pilotto JH, et al. Antiretroviral therapy for the prevention of HIV-1 transmission. N Engl J Med. 2016;375(9):830–9.

    Article  CAS  Google Scholar 

  22. Rodger AJ, Cambiano V, Bruun T, Vernazza P, Collins S, van Lunzen J, Corbelli GM, Estrada V, Geretti AM, Beloukas A, et al. Sexual activity without condoms and risk of HIV transmission in Serodifferent couples when the HIV-positive partner is using suppressive antiretroviral therapy. Jama. 2016;316(2):171–81.

    Article  Google Scholar 

  23. Gao F, Zhong X, Peng B, Zhang Y, Liang H, Zou Y, Dai J, Liu Q, Song B, Huang W, et al. Pre-exposure prophylaxis for HIV high-risk groups: an analysis of clinical compliance. Acad J Second Mil Univ. 2012;33(11):1186–93.

    Google Scholar 

  24. Zhang Y, Zhong X, Peng B, Huang W, Huang A. Cohort retention and predictors in HIV pre-exposure prophylaxis clinical trial among men who have sex with men. Journal of Third Military Medical University. 2013;35(22):2469–72.

    Google Scholar 

  25. Cameron AC, Trivedi PK. Regression analysis of count data, second edition. Edn. Cambridge. New York, NY: Cambridge University Press; 2013.

    Book  Google Scholar 

  26. Gold MR. Cost-effectiveness in health and medicine. New York: Oxford University Press; 1996.

    Google Scholar 

  27. Punyacharoensin N, Edmunds WJ, De Angelis D, Delpech V, Hart G, Elford J, Brown A, Gill ON, White RG. Effect of pre-exposure prophylaxis and combination HIV prevention for men who have sex with men in the UK: a mathematical modelling study. Lancet HIV. 2016;3(2):e94–e104.

    Article  Google Scholar 

  28. Phillips AN, Cambiano V, Miners A, Lampe FC, Rodger A, Nakagawa F, Brown A, Gill ON, De Angelis D, Elford J, et al. Potential impact on HIV incidence of higher HIV testing rates and earlier antiretroviral therapy initiation in MSM. Aids. 2015;29(14):1855–62.

    Article  CAS  Google Scholar 

  29. Juusola JL, Brandeau ML, Owens DK, Bendavid E. The cost-effectiveness of preexposure prophylaxis for HIV prevention in the United States in men who have sex with men. Ann Intern Med. 2012;156(8):541–50.

    Article  Google Scholar 

  30. Nichols BE, Boucher CA, van der Valk M, Rijnders BJ, van de Vijver DA. Cost-effectiveness analysis of pre-exposure prophylaxis for HIV-1 prevention in the Netherlands: a mathematical modelling study. Lancet Infect Dis. 2016;16(12):1423–9.

    Article  Google Scholar 

  31. Cambiano V, Miners A, Dunn D, McCormack S, Ong KJ, Gill ON, Nardone A, Desai M, Field N, Hart G, et al. Cost-effectiveness of pre-exposure prophylaxis for HIV prevention in men who have sex with men in the UK: a modelling study and health economic evaluation. Lancet Infect Dis. 2018;18(1):85–94.

    Article  Google Scholar 

  32. Zhang L, Peng P, Wu Y, Ma X, Soe NN, Huang X, Wu H, Markowitz M, Meyers K. Modelling the epidemiological impact and cost-effectiveness of PrEP for HIV transmission in MSM in China. AIDS Behav. 2018.

  33. Kashuba AD, Patterson KB, Dumond JB, Cohen MS. Pre-exposure prophylaxis for HIV prevention: how to predict success. Lancet (London, England). 2012;379(9835):2409–11.

    Article  Google Scholar 

  34. Traeger MW, Schroeder SE, Wright EJ, Hellard ME, Cornelisse VJ, Doyle JS, Stoove MA. Effects of pre-exposure prophylaxis for the prevention of HIV infection on sexual risk behavior in men who have sex with men: a systematic review and meta-analysis. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America. 2018.

  35. Zhou F, Gao L, Li S, Li D, Zhang L, Fan W, Yang X, Yu M, Xiao D, Yan L, et al. Willingness to accept HIV pre-exposure prophylaxis among Chinese men who have sex with men. PLoS One. 2012;7(3):e32329.

    Article  CAS  Google Scholar 

  36. Lin H, Ding Y, Liu X, Wu Q, Shen W, He N. High prevalence of HIV infection and bisexual networks among a sample of men who have sex with men in eastern China. PLoS One. 2015;10(6):e0129300.

    Article  Google Scholar 

  37. Zhang Y, Peng B, She Y, Liang H, Peng HB, Qian HZ, Vermund SH, Zhong XN, Huang A. Attitudes toward HIV pre-exposure prophylaxis among men who have sex with men in western China. AIDS Patient Care STDs. 2013;27(3):137–41.

    Article  CAS  Google Scholar 

  38. Ochalek J, Lomas J, Claxton K. Cost per DALY averted thresholds for low- and middle-income countries: evidence from cross country data. Working Papers; 2015.

    Google Scholar 

  39. Woods B, Revill P, Sculpher M, Claxton K. Country-level cost-effectiveness thresholds: initial estimates and the need for further research. Value in health : the journal of the International Society for Pharmacoeconomics and Outcomes Research. 2016;19(8):929–35.

    Article  Google Scholar 

  40. World Health Organization. Making choices in health. Geneva: WHO guide to cost-effective- ness analysis; 2003.

    Google Scholar 

  41. WHO Commission on Macroeconomics and Health. Macro- economics and health. Geneva, Switzerland: investing in health for economic development; 2001.

    Google Scholar 

Download references


Not applicable.


The study was supported by the National Natural Science Foundation of China under a young scientists’ grant (81803334), a Major Infectious Disease Prevention and Control of the National Science and Technique Major Project (2018ZX10715004), and a National Centre for Global Health and Medicine Research Development grant (30A-1007).

Availability of data and materials

All data generated or analysed during this study are included in this published article and its supplementary information files.

Author information

Authors and Affiliations



JL, SG, and JTFL were responsible for the conception and design of this study. SG and JL worked on the development, coding, and analysis of the mathematical model. JL, SG and LP wrote the first draft. JL, LP, JG, HZ, YR, CH and YH collected and collated data for the inputs to the model. JTFL, YR, JG, HZ and YH provided scientific discussion and comments. JL, SG, JG, YR, LP, HZ, CH, YH, and JTFL read, revised and approved the final manuscript.

Corresponding author

Correspondence to Stuart Gilmour.

Ethics declarations

Ethics approval and consent to participate

Not applicable, as this was a mathematical model.

Consent for publication

Not applicable.

Competing interests

The author declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Supplementary Appendix. A complete description of the mathematical methods used in development of the model. (DOCX 5387 kb)

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, J., Peng, L., Gilmour, S. et al. A mathematical model of biomedical interventions for HIV prevention among men who have sex with men in China. BMC Infect Dis 18, 600 (2018).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: