Skip to main content

New sexually transmitted HIV infections from 2016 to 2050 in Guangdong Province, China: a study based on a dynamic compartmental model

Abstract

Background

In Guangdong Province, China, there is lack of information on the HIV epidemic among high-risk groups and the general population, particularly in relation to sexual transmission, which is a predominant route. The new HIV infections each year is also uncertain owing to HIV transmission from men who have sex with men (MSM) to women, as a substantial proportion of MSM also have female sexual partnerships to comply with social demands in China.

Methods

A deterministic compartmental model was developed to predict new HIV infections in four risk groups, including heterosexual men and women and low- and high-risk MSM, in Guangdong Province from 2016 to 2050, considering HIV transmission from MSM to women. The new HIV infections and its 95% credible interval (CrI) were predicted. An adaptive sequential Monte Carlo method for approximate Bayesian computation (ABC-SMC) was used to estimate the unknown parameter, a mixing index. We calibrated our results based on new HIV diagnoses and proportions of late diagnoses. The Morris and Sobol methods were applied in the sensitivity analysis.

Results

New HIV infections increased during and 2 years after the COVID-19 pandemic, then declined until 2050. New infections rose from 8,828 [95% credible interval (CrI): 6,435–10,451] in 2016 to 9,652 (95% CrI: 7,027–11,434) in 2019, peaking at 11,152 (95% CrI: 8,337–13,062) in 2024 before declining to 7,084 (95% CrI: 5,165–8,385) in 2035 and 4,849 (95% CrI: 3,524–5,747) in 2050. Women accounted for approximately 25.0% of new HIV infections, MSM accounted for 40.0% (approximately 55.0% of men), and high-risk MSM accounted for approximately 25.0% of the total. The ABC-SMC mixing index was 0.504 (95% CrI: 0.239–0.894).

Conclusions

Given that new HIV infections and the proportion of women were relatively high in our calibrated model, to some extent, the HIV epidemic in Guangdong Province remains serious, and services for HIV prevention and control are urgently needed to return to the levels before the COVID-19 epidemic, especially in promoting condom-based safe sex and increasing awareness of HIV prevention to general population.

Peer Review reports

Background

HIV remains a worldwide public health threat, while the efforts for prevention and control efforts have made substantial progress, including reducing of HIV-related morbidity and mortality and increasing life expectancy and access to antiretroviral therapy. The number of people living with HIV was 38.4 million according to the Joint United Nations Programme on HIV/AIDS (UNAIDS) at the end of 2021 [1], and the number of new HIV infections was approximately 1.5 million [2]. The number of people living with HIV was estimated to be 1.25 million in China in 2018, with approximately 80,000 new infections that year [3]. As of 2022, the number of people living with HIV was over 80 thousand in Guangdong Province, which ranked fourth out of 31 provinces [4].

Guangdong Province is located in southern China and has a relatively advanced economy and a higher proportion of sexually active individuals [5]. The number of new HIV diagnoses each year is high and increasing owing to the complex mode of sexual transmission, which has been the main route since 2009 [6], accounting for 52.4%, and increasing continually to 90.3% in 2014. However, the number of new HIV infections predicted in Guangdong Province has not been reported in published peer-reviewed articles. It is important to predict the number of new HIV infections in provinces where incidences of HIV may be different for prevention and control planning, especially since heterosexual transmission is the dominant transmission mode in China.

HIV is increasingly difficult to prevent and control due to intricate modes of sexual transmission in China, including heterosexual transmission from heterosexual men or men who have sex with men (MSM) in heterosexual acts and homosexual transmission. It is observed that HIV began to spread to the general population from key populations in 2007 [4, 7]. Furthermore, impacted by the traditional culture including marital pressures and filial expectations [8, 9], the mode of sexual transmission makes the HIV epidemic complex and diverse due to mixed sexual acts in marital or nonmarital partnerships, commercial or noncommercial partnerships, heterosexual or homosexual partnerships, etc [10].

There are different methods to estimate the number of new HIV infections including mathematical models such as compartmental models and the workbook method; however, a compartmental model is often used to estimate the number of infections [11], which is also recommended by UNAIDS [12,13,14]. Dynamic compartmental models that go beyond the workbook method [15] in long-term prediction are also classical mathematical models of infectious disease transmission, but there are two considerable challenges in developing a compartmental model to predict new HIV infections, including defining risk groups related to MSM depending on the characteristics of a local epidemic and defining a criterion or a reference for model calibration in addition to parameter estimation.

One definition of MSM considers sex of their partners in sexual acts in addition to their number of sexual partnerships in predictive models of new HIV infections; some MSM may have female sexual partnerships in some countries and are also known as men who have sex with men and women (MSMW). The marriage rate among Chinese MSM ranges from 31.2 to 70% [16,17,18,19] and is even higher during the lifetime [20], and the number of wives (tongqi) of MSM is between 1 million and 16 million [21]. However, MSM included only two subgroups depending on their number of sexual partnerships in previous compartmental models [22, 23]. This was inconsistent with the characteristics of an epidemic in these countries as it ignored the risk of HIV transmission [24,25,26], which was higher among MSMW than among MSM only [27], even though MSM are often classified into two subgroups that make interpretability insufficient between the number of predictions and the number of case reports [3]. MSMW may transmit HIV [26, 28]. With the HIV prevalence among MSM increasing, the HIV incidence among their female partners increased by 5.3 times from 2002 to 2010 [29] because MSM failed to disclose their sexual orientation and HIV seropositivity to their wives or female sexual partnerships to satisfy the social and familial expectations [28, 30] in general, thus reducing the awareness of HIV prevention [31] during heterosexual acts with lower condom use than during homosexual acts due to reproductive purposes [18, 32] or other reasons [33,34,35].

Another challenge is how to determine a criterion for model calibration for a compartmental model. Given that the problem of late diagnosis is prevalent worldwide [36,37,38,39] and continues to be prevalent [40], model calibration must also consider the proportion in addition to the number of new HIV diagnoses each year, which was often employed to calibrate owing to a lack of other data. The number predicted must be almost higher than the number of new HIV diagnoses each year in most countries due to late diagnoses [37, 38, 41, 42], whose proportion may increase owing to the impacts of the COVID-19 pandemic [43]; otherwise, the predicted number will be underestimated to some extent. Hence, it is critical to determine the proportion of people diagnosed with HIV infections for model calibration, and the sources should be reasonable and credible and include published studies, AIDS expert interviews or internal estimations from the Centers for Disease Control and Prevention (CDC). This study aimed to predict the number of new HIV infections by developing a compartmental model for the population aged 15 and over considering two sexual transmission routes to further understand the HIV epidemic. Model calibration included the proportion of late diagnoses in addition to the number of new HIV diagnoses each year (almost equal to the number of new cases reported from medical and health institutions at all levels in China because of compulsory guidelines mandated by the National Notifiable Disease Reporting System).

Methods

We developed a deterministic compartmental model with four states to predict the HIV epidemic in the population aged 15 and over in Guangdong Province from 2016 to 2050. The model predicted the number of new HIV infections and its 95% credible interval (CrI). We used the Morris and Sobol methods to analyse the sensitivity of the model parameters. The predicted number was calibrated by comparing it with the number of yearly new HIV diagnoses (almost equal to the number of new cases reported in China) and the potential proportion of late diagnoses reported by AIDS experts in the Guangdong Provincial CDC or published peer-reviewed articles [37]. The model structure, population size of risk groups, model parameters, sensitivity analysis and model calibration are as follows.

Model structure

Considering births and deaths, the basic structure of a compartmental model with four states is shown in Fig. 1. Four ordinary differential equations describing the states are shown in formulas (1)–(4), denoting four states of the HIV epidemic, four risk groups, and two sexual transmission routes among MSM. The four states were S (susceptible), I (infected), D (diagnosed), and T (treated); the infectivity of people in the I and D states and the antiretroviral treatment failure of those in the T state were evaluated. We defined MSM as the high-risk MSM who satisfied at least one of two conditions [22, 23]: (1) more than 10 sexual partners over the past six months [44] and (2) rates of inconsistent condom use over the past six months more than 50.0%, which were approximately 50.0% in Guangdong Province [45]; otherwise, they were classified as low-risk MSM. The four risk groups were heterosexual men, heterosexual women, low-risk MSM, and high-risk MSM. The two sexual transmission routes were heterosexual and homosexual. The probability of HIV acquisition is denoted as \(\lambda \left(t\right)\) in the S population per year. The I(t) is denoted as \(\lambda \left(t\left)S\right(t\right)\). The number of people engaging in potentially high-risk sexual acts per year is denoted as \(N\left(t\right)=S\left(t\right)+I\left(t\right)\), with I (t) including D (t) and T (t). \({\upzeta }\left(t\right)\) is the entry rate per year, which is the sum of the natural birth rate and natural death rate. \(d\left(S,I,D,T\right)\) is the death rate in the states for the risk groups.

Fig. 1
figure 1

A schematic diagram of the compartmental model of HIV transmission

$$\frac{{{X_{j,1,i}}}}{{dt}} = \sum\limits_{j = 1}^4 {{\zeta _j}} {X_{j,1,i}} - \left( {\sum\limits_{k = {\text{2,3}},{4^\prime }} {\lambda _{i,k}^j} \left( t \right)} \right){X_{j,1,i}} - {d_{j,1}}{X_{j,1,i}}$$
(1)
$$\frac{{X}_{j,2,i}}{dt}=\left(\sum _{k=\text{2,3},{4}^{{\prime }}}{\lambda }_{i,k}^{j}\left(t\right)\right){X}_{j,1,i}-{d}_{j,2}{X}_{j,2,i}$$
(2)
$$\frac{{X}_{j,3,i}}{dt}={\delta }_{j}{X}_{j,2,i}-{d}_{j,3}{X}_{j,3,i}$$
(3)
$$\frac{{X}_{j,4,i}}{dt}={\psi }_{j}{X}_{j,3,i}-{d}_{j,4}{X}_{j,4,i}$$
(4)
$$j=\text{1,2},\text{3,4}, i=1, 2$$

j denotes the four risk groups, j\(=\text{1,2},\dots ,4\), 1 indicates heterosexual men, 2 heterosexual women, 3 low-risk MSM, and 4 high-risk MSM. i denotes HIV sexual transmission routes, \(i=\text{1,2}\), 1 indicates heterosexual sexual acts, and 2 indicates homosexual sexual acts. k denotes the four states, \(k=\text{1,2},\dots ,4\)’, where 1 is S, 2 is I, 3 is D, 4 is T, and \(4\)’ is treatment failure.

Population sizes of the risk groups

The number of people in the risk groups, including heterosexual men, heterosexual women, low-risk and high-risk MSM, and low-risk and high-risk MSMW, whose risk groups among MSM or MSMW were classified based on the number of male sexual partners among MSM, were calculated on the basis of the population size and the proportion of people aged 15 and over by sex from then the 2020 Guangdong Statistical Yearbook, which includes data from 2019 and before [46]. Data from 2020 were used with the assumption that the population size has remained relatively steady. Presuming that the population size of low-risk MSM was equal to that of high-risk MSM [45], the proportion of MSM was 5.0% of all men aged 15 and over, and the proportion of heterosexual sexual acts among low- and high-risk MSM was 31.2% [16]. By excluding people already living with HIV, we estimated that the population sizes of heterosexual men, low- and high-risk MSM, and low- and high-risk MSMW were 47.89 million, 1.3 million, 1.3 million, 0.4 million, and 0.4 million, respectively. The population size of heterosexual women was 43.73 million, which excluded people living with HIV and a number of lesbians, whose proportion was presumed to be 5.0%, as the risk of HIV transmission is extremely rare in this group.

Model parameters

Model parameters included three parts depending on whether they were known or unknown (unknown parameters: the probability of HIV acquisitions and a mixing index). The known parameters that we collected (see the Supplementary material) included demographic, behavioural, biological and epidemiological data, coming from peer-reviewed published articles, domestic government reports, AIDS expert interviews, and the viewpoints of key experts on World AIDS Day; these parameters were also employed in other studies and meta-analyses. Demographic parameters included the population sizes of the four groups of people aged 15 and over and the number of people who entered and left this population yearly because of natural population growth. Behavioural parameters included sexual partners, condom use rates and the effectiveness of condoms, the last two of which differed between heterosexual and homosexual sexual acts. Biological parameters included the probability of HIV transmission per sexual act and death rates in different states. Epidemiological parameters included HIV prevalence, rates of new diagnoses, and rates of antiretroviral therapy.

The first unknown parameter was the probability of HIV acquisition, which was a function of the probability of not acquiring HIV during a high-risk sexual act, the number of sexual acts, the number of sexual partners, condom use rates and the effectiveness of condoms. We presumed that HIV transmission in our model occurred only via sextual transmission, including heterosexual and homosexual acts, due to sexual transmission being the dominant route since 2009 in Guangdong Province; the rate of sexual transmission increased continually to 90.3% in 2014. We also presumed that heterosexual men engaged in only heterosexual sexual acts, and women were only infected by men who could also be MSMW. These men were indistinguishable from heterosexual men participating in heterosexual sexual acts due to the lack of a relevant compartment for heterosexual women in the compartmental model. The probability of acquiring HIV in a high-risk sexual act among MSM was calculated based on the total number of homosexual sexual acts, of which insertive and receptive anal sex were indistinguishable. People living with HIV who have experienced antiretroviral treatment failure can spread HIV, or they are not infectious and should achieve a suppressed viral load. The function and its calculation are detailed in the Appendix.

The last unknown parameter was a mixing index that we used approximate Bayesian computation (ABC) [47, 48] to estimate. The index was a randomized mixing level of opting for male or female sexual partners among MSM, also called the assortativity, ranging from 0 to 1 [49], where 0 denoted that MSM would only opt for male sexual partners and 1 denoted that MSM would opt for male or female sexual partners completely at random. The number of HIV infections approached to the number of HIV diagnoses as the randomized mixing level increased [49]. Finally, an adaptive sequential Monte Carlo (SMC) sampler [50] was adopted considering the cost of calculations, given that a value for the mixing index could not be acquired from an epidemiological research or published peer-reviewed studies, which only presented an unconvincing presumed value of 0.5 [23].

Sensitivity analysis

Considering the cost of calculations, we combined a qualitative global method and a quantitative global method, the Morris and Sobol methods, to analyse the sensitivity of model parameters to the number of new HIV infections predicted in our model. The first step qualitatively determined the sensitive parameters via the Morris method [51] and then calculated the quantitative effect of those identified parameters on the number of new HIV infections predicted via the Sobol method [52]. The sensitivity analysis excluded the parameter of the mixing index because it was certainly sensitive. The outputs of the Sobol method included the total-order indices and the first-order indices on qualitatively identified parameters from the absolute values of elementary effects (|EEs|) of the Morris method.

Model calibration

The model calibration data included the following: the yearly new HIV diagnoses based on the system of case reports on government websites [53, 54] on World AIDS Day from 2016 to 2022; and the proportion of HIV diagnoses based on published peer-reviewed studies [41, 55, 56], the viewpoints of a national AIDS expert [3] or provincial AIDS expert interviews [57] from the Guangdong Provincial CDC. Moreover, the number of new HIV infections may have increased during the COVID-19 epidemic due to temporary disruptions in health services and changes in sexual risk acts, including interrupted antiretroviral treatment, reduced HIV testing, and a higher proportion of unprotected sexual acts. Considering the proportion of late diagnoses in Guangdong Province, the difference was measured by determining the relative ratios of the number of new predicted HIV infections divided by the number of new HIV diagnoses each year, which was deemed the “goodness of fit” of the model, and a value larger than 1.2 was applied as a measure of model calibration [37, 41, 58]. A 95% credible interval (CrI) of the posterior parameter from SMC-ABC was adopted to calculate the 95% CrI of the number of new HIV infections. Uncertainty intervals (UI) of it on heterosexual men and women were calculated by realistic scenarios of the minimum and maximum considering sensitivity parameters identified in our sensitivity analysis. To incorporate the potential impact of the COVID-19 pandemic, we made adjustments to the sensitive parameters in our model.

Results

Parameter estimation

The value of the mixing index estimated from ABC-SMC was 0.504 (95% CrI: 0.239, 0.894). The prior presumed distribution of the index was uniformly distributed [0, 1], the sampling sequences were set at 25, the total sample size of each sequence was 10,000, each sampling size was 10, and the effective sample size was also estimated.

Sensitivity analysis

The |EEs| from the Morris method for the number of sexual partners, probabilities of acquiring HIV during high-risk sexual acts, the proportion of condom use and the population sizes of the four states were 3,811.9, 2,441.1, 136.0 and 440.4, respectively. The first three parameters were identified as sensitive parameters, as the population sizes were relatively stable. The Sobol method revealed the following ranking for the three parameters: the number of sexual partners, the probability of acquiring HIV, and the proportion of condom use (refer to Supplementary Fig. 1 for more information). The total-order indices for these parameters were calculated to be 0.73, 0.29, and 0.0009, respectively, while the first-order indices were found to be 0.70, 0.27, and 0.0007, respectively.

Model calibration

The ratios between the predicted number of new HIV infections and the number of new diagnoses per year are also shown as relative ratios, with further adjustments of the sensitivity parameters on the basis of the unknown parameter estimation for the population that experienced sexual transmission in Guangdong Province. The relative ratios were all larger than 1.2 and were calculated according to the proportion of HIV diagnoses or late diagnoses in Guangdong Province (Table 1). The ratios during the COVID-19 pandemic were much larger than those before the pandemic.

Table 1 Model calibration results: Relative ratios larger than 1.2

Number of new infections

The number of new HIV infections in our compartmental model showed that the prediction period was comprehensively divided into five stages on the basis of changing trends (see Fig. 2), including before the COVID-19 epidemic (2016–2019), the COVID-19 epidemic (2020–2022), two or three years after the COVID-19 epidemic (2023–2025), the decade after the COVID-19 epidemic (2026–2035), and the period until 2050. The number of new HIV infections slightly increased from 2016 to 2019, sharply increased from 2020 to 2022, continuously increased from 2023 to 2025, declined yearly from 2026 to 2035, and remained at a stable lower level and slowly decreased until 2050. The number and proportion of new HIV infections predicted for the four risk groups were as follows.

The number of new HIV infections predicted each year in Guangdong Province was found to be 8,828 [95% credible interval (CrI): 6,435–10,451] in 2016, slightly increasing to 9,652 (95% CrI: 7,027–11,434) in 2019, continually jumping to 11,152 (95% CrI: 8,337–13,062) in 2024, beginning to decrease in 2025, then persistently declining to 7084 (95% CrI: 5,165–8,385) in 2035 and continuing to decline to 4,849 (95% CrI: 3,524–5,747) until 2050. The proportion of the predicted number of new HIV infections totally fell by approximately 45.0% between 2016 and 2050, as shown in Table 2. The increase in the predicted number of new HIV infections was much larger during the COVID-19 pandemic and the two years after the end of the pandemic, and then the number continuously declined until 2050. The changing tendency of the number of new HIV infections predicted for the four risk groups was in line with the changing tendency of the total number predicted each year (see Fig. 2).

Table 2 Predicted number of new HIV infections in four susceptible populations in Guangdong Province from 2016 to 2050

The number of new HIV infections and its proportion in the four risk groups from our model varied, as shown in Table 2. The predicted number of the new HIV infections between 2016 and 2050 among heterosexual women was 2,067 (UI: 2,359–3,016) in 2016, which increased to 2,835 (UI: 2,259–2,879) in 2019, jumped to 4,346 (UI: 3,612–4,735) in 2025, and then decreased to 1,228 (UI: 1,008–2,145) in 2050; among heterosexual men, the number was 2,839 (UI: 2,237–3,243) in 2016, which increased to 3,189 (UI: 2,698–3,842) in 2018, and then decreased to 1,751 (UI: 1,427–2,738) in 2050; among high-risk MSM, the number was 2,288 (95% CrI: 660–3,393) in 2016, which slightly increased to 2,568 (95% CrI: 756–3,798) in 2025 and then decreased to 1,046 (95% CrI: 304–1,549) in 2050; and among low-risk MSM, the number was 1,094 (95% CrI: 329–1,612) in 2016, which increased to 1,209 (95% CrI: 357–1,787) in 2019 then to 1,487 (448–2,193) in 2026 and decreased to 824 (95% CrI: 242–1,219) in 2050. The predicted number of new HIV infections can be broken down as follows: women accounted for approximately 25.0% of the total, homosexual transmissions accounted for approximately 40.0% of the total, high-risk MSM accounted for approximately 25.0% of the total, and MSM accounted for 55.0% of men.

Fig. 2
figure 2

Predicted number (95% CrI) of new HIV infections in the population aged 15 and over from 2016 to 2050 in Guangdong Province, China

Discussion

Despite including the period of the COVID-19 pandemic in this prediction, the predicted number of new HIV infections was similar to the number during the HIV epidemic in Guangdong based on the relative ratios from 2016 to 2022. As measures and strategies for HIV prevention and control were carried out regularly before the COVID-19 pandemic in 2016 and 2019, the number of new HIV infections also increased slightly, similar to the number of new HIV diagnoses, which showed a changing trend. Overall, the predicted number increased during the COVID-19 pandemic and for approximately two or three years after the pandemic. The proportion of women out of the total predicted number peaked at 39% during this period. Then, the predicted number decreased continually until 2050.

A range of calibrated criteria, the relative ratios, between the predicted number in this model and HIV diagnoses had values larger than 1.2, fitted with the proportions of HIV diagnoses from Guangdong Provincial AIDS expert interviews on World AIDS Day, which was 71.3% in 2019 [58] and 78.7% in 2021 [57], and the interpretations or viewpoints of national experts regarding the predicted results in China evaluated by Spectrum, which was recommended worldwide [3]. The predicted number of new HIV infections jumped during the COVID-19 pandemic, with the calibrated criterion ranging from 1.4 to 1.7 due to the interruption of or delay in HIV services provided by health care workers who had been diverted for COVID-19 [59] to different extents. These services included HIV testing and counselling, referrals, timely antiretroviral treatment, and the promotion in condom use. Previous studies that used simulation models and cross-sectional data reported that the COVID-19 pandemic should have led to an increase in the number of new HIV infections [43, 60]. HIV diagnoses in another province of China declined by approximately 37.0% during the COVID-19 pandemic compared to the same period before the pandemic, according to a prediction model [61]. Furthermore, there were no convincing data to imply that the variety of sexual behaviours in risk groups differed from that before the COVID-19 pandemic [62].

The predicted number of new HIV infections increased continually for two or three years after the end of the COVID-19 pandemic. Recognizing that high-risk individuals may have been unaware that they had HIV due to reductions in HIV testing, there should have been an increase in HIV transmission because those individuals who worried about contracting COVID-19 did not go to designated institutions for HIV testing. As the COVID-19 epidemic was mitigated, HIV testing failed to promptly reduce HIV transmission, but risk behaviours for HIV may have been unchanged and potentially increased [63] when the COVID-19 epidemic was mitigated. Cross-sectional data showed that the HIV-positive rate of blood donors was much higher than that before the COVID-19 pandemic [64]. The predicted new HIV infections and the trend after the COVID-19 pandemic showed brief increases and then steadily declined after measures and strategies for prevention and control were reinstated, especially scaled up HIV testing [65], timely antiretroviral treatment for newly diagnosed infections, and viral suppression.

Although we used proportions of late diagnoses in addition to new HIV diagnoses from 2016 to 2022 in Guangdong Province to calibrate our model, it was and remains uncertain how many new HIV infections there were and how many were diagnosed in China based on existing research and the viewpoints of national or provincial core experts. The diagnosed proportion of new HIV infections was 68.9% (61.5%, 78.3%) [55] based on predictions made in 2018 by Spectrum regarding the HIV epidemic in China. The diagnosis rate of HIV infections remains low, and the HIV epidemic may be increasing continuously, which is in contrast to predictions from the workbook method and core experts indicating that the HIV epidemic has been stable since 2007 in China [3].

Nevertheless, regarding HIV transmission from MSMW to general women [28], the sex compositions and proportions of heterosexual and homosexual sexual acts associated with the number of new HIV infections from our model were similar to those related to the number of new HIV diagnoses from 2016 to 2019 in Guangdong. The female proportion of new HIV diagnoses and the proportion of heterosexual transmission of new HIV diagnoses in Guangdong during recent years were approximately 30.0% and over 60.0% [66], respectively, which were similar to the national proportions, with the proportions were approximately 30.0% and 70.0% among people living with HIV, respectively [3]. The female proportion of the new HIV infections was higher during the COVID-19 epidemic, and after two years, for the proportion of heterosexual acts among MSMW may increase due to a decline of approximately one-fifth in the number of male sexual partners among MSMW in that period [43], which further could lead to increase for heterosexual transmission among MSMW.

The model parameters and population size in our model were selected based on more reliable sources. The ranges of the parameters were set according to highly cited studies and the viewpoints of core experts, and the prior distribution of the mixing index was reasonable. The model parameters of the HIV epidemic in Guangdong in our model were preferable, and they were replaced with parameters from another province with a similar HIV epidemic were replaced if these parameters for Guangdong were not acquired; otherwise, the values from a meta-analysis were used for China. For some parameters, including the effectiveness of condoms, we opted for values from highly cited studies.

Our modelling study had three strengths. First, the assumptions in the model were in line with the characteristics of the HIV epidemic in Guangdong Province and Chinese social norms and pressures on MSM, including low- and high-risk subgroups classified according to the number of sexual partnerships and the higher proportion of heterosexual sexual acts in China [17, 20]. Second, the calibrated criteria in the model also included the proportion of late diagnoses, which was approximately 30.0% and even higher during the COVID-19 epidemic, in addition to the new HIV diagnoses each year in Guangdong Province. Finally, this is the first compartmental model to publicly predict the number of new HIV infections in members of the population with susceptibility to sexual transmission from 2016 to 2050 in Guangdong Province. However, this study also had four limitations. First, our model failed to identify the risk groups of sexual partnerships, MSMW or heterosexual men, among females with HIV infections. The problem could be further explored by defining of subgroups for females in our model. Second, our model failed to subclassify according to CD4 cell counts, considering that viral suppression may be more reasonably defined as a state of being “Treated”, starting in 2016 for all people with HIV initiating antiretroviral treatment regardless of clinical stage and CD4 cell count. Third, the mixing index failed to randomize between low- and high- risk MSM but did randomize between heterosexual and homosexual sexual acts because HIV transmission from MSM to women may be substantial, impacted by Chinese traditional cultures. Finally, the calibration criterion for our model would also lead to an underestimation of new HIV infections, while we have taken into account the proportion of late diagnoses. And determining an upper value for model calibration is challenging due to the lack of available data.

Conclusion

We developed a deterministic compartmental model to predict the number of new HIV infections from 2016 to 2050. The presumption regarding MSM who engaged in heterosexual acts and the criterion of model calibration may all be in line with the complicated characteristics of the epidemic in Guangdong Province to some degree. The results from the model may simulate a realistic epidemic in Guangdong Province by model calibration. The predicted number slightly increased between 2016 and 2019, was much larger during the COVID-19 pandemic until two years after the end of the pandemic, and then continually declined until 2050. Overall, the HIV epidemic in Guangdong Province remains serious, and it is urgently important to restore services for HIV control and prevention after the COVID-19 pandemic.

Data availability

No datasets were generated or analysed during the current study.

References

  1. In danger: UNAIDS global AIDS update 2022. Geneva: Joint United Nations Programme on HIV/AIDS (UNAIDS). 2022: 354.

  2. UNAIDS core epidemiology slides. Geneva: Joint United Nations Programme on HIV/AIDS (UNAIDS). 2022: 2.

  3. Wu ZY, Mcgoogan JM, Detels R. The enigma of the human immunodeficiency virus (HIV) epidemic in China. Clin Infect Dis. 2021;72(5):876–81.

    Article  PubMed  Google Scholar 

  4. Wu ZY, Chen JF, Scott SR, Mcgoogan JM. History of the HIV epidemic in China. Curr HIV/AIDS Rep. 2019;16(6):458–66.

    Article  PubMed  PubMed Central  Google Scholar 

  5. National Bureau of Statistics of China. 2022 China statistical yearbook: population by province or regions in China in 2022. Beijing: China Statistics Press; 2022.

  6. Li Y, Lin P, Wang Y, Fu X, Ye R, Li J, Feng M, Yu G, Yang F, Yan J, et al. HIV/AIDS surveillance in Guangdong province, 2000–2013. South China J Prev Med. 2015;41(2):101–6.

    Google Scholar 

  7. Lu L, Jia M, Ma Y, Yang L, Chen Z, Ho DD, et al. The changing face of HIV in China. Nature. 2008;455(7213):609–11.

    Article  CAS  PubMed  Google Scholar 

  8. Li X, Zhang B, Wang J, Li Y, Li X, Yu P, Liu M, Liu X. Sexual health status of women who have regular sexual relations with men who have sex with men in mainland China. BMC Public Health. 2017;17(1):168.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Koo FK, Chow EP, Gao L, Fu X, Jing J, Chen L, et al. Socio-cultural influences on the transmission of HIV among gay men in rural China. Cult Health Sex. 2014;16(3):302–15.

  10. Lin Z, Li Y, Fu X, Liu J, Guan Z, Li J, Yang F, Lin P. Study on modes of HIV transmission through heterosexual contact in Guangdong Province, 2015–2016. Chin J AIDS STD. 2017;23(11):1002–5.

    Google Scholar 

  11. Mahy M, Lewden C, Brinkhof MW, Dabis F, Tassie JM, Souteyrand Y, et al. Derivation of parameters used in Spectrum for eligibility for antiretroviral therapy and survival on antiretroviral therapy. Sex Transm Infect. 2010;86(Suppl 2):ii28–34.

    PubMed  Google Scholar 

  12. Khalifa A, Stover J, Mahy M, Idele P, Porth T, Lwamba C. Demographic change and HIV epidemic projections to 2050 for adolescents and young people aged 15–24. Glob Health Action. 2019;12:1662685.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Case KK, Johnson LF, Mahy M, Marsh K, Supervie V, Eaton JW. Summarizing the results and methods of the 2019 Joint United Nations Programme on HIV/AIDS HIV estimates. AIDS. 2019;33(Suppl 3):S197–201.

    Article  PubMed  Google Scholar 

  14. Stover J, Glaubius R, Kassanjee R, Dugdale CM. Updates to the Spectrum/AIM model for the UNAIDS 2020 HIV estimates. J Int AIDS Soc. 2021;24(Suppl 5):e25778.

  15. Lyerla R, Gouws E, García-Calleja JM, Zaniewski E. The 2005 workbook: an improved tool for estimating HIV prevalence in countries with low level and concentrated epidemics. Sex Transm Infect. 2006;82(suppl 3):iii41–4.

    PubMed  PubMed Central  Google Scholar 

  16. Yun K, Xu JJ, Reilly KH, Zhang J, Jiang YJ, Wang N, Shang H. Prevalence of bisexual behaviour among bridge population of men who have sex with men in China: a meta-analysis of observational studies. Sex Transm Infect. 2011;87(7):563–70.

    Article  CAS  PubMed  Google Scholar 

  17. Dai Z, Zhong X, Peng B, Zhang Y, Liang H, Peng H, Zhong XH, Liu X, Huang A. Socio-demographics, sexual behaviours, and use of HIV prevention services among men who have sex with men and women in Western China. Int J STD AIDS. 2016;27(2):133–40.

    Article  PubMed  Google Scholar 

  18. Chow EP, Wilson DP, Zhang L. What is the potential for bisexual men in China to act as a bridge of HIV transmission to the female population? Behavioural evidence from a systematic review and meta-analysis. BMC Infect Dis. 2011;11:242.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Zhang B, Li X, Shi T, Yang L, Zhang J. Initial estimates of population size and a HIV prevalence among MSM or biosexual men in China. Chin J STD/AIDS Pre Cont. 2002;8(4):197–9.

    Google Scholar 

  20. 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  PubMed  PubMed Central  Google Scholar 

  21. Yan F, Shi CX, Zhang Y, Chen J, Williams AB, Li X. I’ll Change his sexual orientation, I don’t think about HIV: a qualitative study to explore attitudes, behaviors, and experiences among wives of men who have sex with men in Mainland China. J Assoc Nurses AIDS Care. 2020;31(4):428–38.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Booton RD, Ong JJ, Lee A, Liu A, Huang W, Wei C, Tang W, Ma W, Vickerman P, Tucker JD, et al. Modelling the impact of an HIV testing intervention on HIV transmission among men who have sex with men in China. HIV Med. 2021;22(6):467–77.

    Article  PubMed  Google Scholar 

  23. 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. 2019;23(2):523–33.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Cummings B, Horth R, McFarland W, Lane T, Young PW, Nalá R, Rutherford G, Raymond HF. Bisexual and bidirectional: assessing the potential for HIV bridging in Mozambique. AIDS Behav. 2018;22(7):2189–98.

    Article  PubMed  Google Scholar 

  25. Hong H, Xu GZ, Zhang DD. Prevalence of HIV and sexually transmitted diseases among female and male partners of men who have sex with men in Ningbo, China. Int J Gynaecol Obstet. 2014;125(1):83.

    Article  PubMed  Google Scholar 

  26. Cai R, Zhao J, Cai W, Chen L, Richardus JH, de Vlas SJ. HIV risk and prevention behaviors in men who have sex with men and women: a respondent-driven sampling study in Shenzhen, China. AIDS Behav. 2014;18(8):1560–8.

    PubMed  Google Scholar 

  27. Wu X, Dou Q, Bao Y, Zhang Y, Tan H. Characteristics of HIV’s sexual behavior and their effect on the secondary transmission rate. J Cent South Univ (Med Sci). 2016;41(7):715–20.

    Google Scholar 

  28. He Q, Wang Y, Lin P, Liu Y, Yang F, Fu X, et al. Potential bridges for HIV infection to men who have sex with men in Guangzhou, China. AIDS Behav. 2006;10(Suppl 4):S17–23.

    Article  PubMed  Google Scholar 

  29. Chow EP, Wilson DP, Zhang L. Estimating HIV incidence among female partners of bisexual men in China. Int J Infect Dis. 2012;16(5):e312–320.

    Article  PubMed  Google Scholar 

  30. Chen F, Ding F, Lin X, Wang X, He H, Huang W, Dong Y, Yang Q, Zhang H. [The status of HIV disclosure to primary sexual partners and sexual behaviors among HIV-infected men who have sex with men of 3 cities in China]. Zhonghua Yu Fang Yi Xue Za Zhi. 2014;48(11):969–73.

    PubMed  Google Scholar 

  31. Yan H, Cao W, Mo P, Huan X, Wang Z, Lin X, Wang X, Gu L, Wang P, Agudile E, et al. Prevalence and associated factors of HIV serostatus disclosure to regular female sex partners among HIV-positive men who have sex with both men and women in China. AIDS Care. 2019;31(8):1026–34.

    Article  PubMed  Google Scholar 

  32. Chen L, Jiang T, Wang H, Hong H, Ge R, Tang H, Wang S, Xu K, Chai C, Ma Q, et al. Development and validation of a risk score for predicting inconsistent condom use with women among men who have sex with men and women. BMC Public Health. 2023;23(1):734.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Shen H, Tang S, Mahapatra T, Tucker JD, Huang S, Yang B, Zhao J, Detels R, Tang W. Condomless vaginal intercourse and its associates among men who have sex with men in China. PLoS ONE. 2016;11(4):e0154132.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Pandya A, Pandya S, Patil B, Merchant S. Invisible others: sexual health vulnerabilities of wives of men who have sex with men. Sex Cult. 2011;16(1):76–89.

    Article  Google Scholar 

  35. Guo H. A study of the sexual behaviour characteristics among married men who have sex with men in Beijing. Beijing: Chinese CDC; 2010.

    Google Scholar 

  36. The Lancet HIV. Time to tackle late diagnosis. Lancet HIV. 2022;9(3):e39.

  37. Jiang H, Liu J, Tan Z, Fu X, Xie Y, Lin K, Yan Y, Li Y, Yang Y. Prevalence of and factors associated with advanced HIV disease among newly diagnosed people living with HIV in Guangdong Province, China. J Int AIDS Soc. 2020;23(11):e25642.

    Article  PubMed  PubMed Central  Google Scholar 

  38. Nishiura H. Estimating the incidence and diagnosed proportion of HIV infections in Japan: a statistical modeling study. PeerJ. 2019;7:e6275.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Wong CS, Wei L, Kim YS. HIV Late Presenters in Asia: management and public health challenges. AIDS Res Treat. 2023;2023:9488051.

  40. Collins S, Namiba A, Sparrowhawk A, Strachan S, Thompson M, Nakamura H. Late diagnosis of HIV in 2022: why so little change? HIV Med. 2022;23(11):1118–26.

    Article  PubMed  Google Scholar 

  41. Lin Z, Li Y, Fu X, Liu J, Guan Z, Li J, Yang F, Lin P. Characteristics and influencing factors for late diagnosed HIV/AIDS cases in Guangdong Province, 2010–2016. China J Prev Med. 2017;43(6):501–6.

    Google Scholar 

  42. Castel C, Sommen C, Strat YL, Alioum A. A multi-state Markov model using notification data to estimate HIV incidence, number of undiagnosed individuals living with HIV, and delay between infection and diagnosis: illustration in France, 2008–2018. Stat Methods Med Res. 2021;30(11):2382–98.

    Article  PubMed  Google Scholar 

  43. Booton RD, Fu G, Macgregor L, Li J, Ong JJ, Tucker JD, Turner KM, Tang W, Vickerman P, Mitchell KM. The impact of disruptions due to COVID-19 on HIV transmission and control among men who have sex with men in China. J Int AIDS Soc. 2021;24(4):e25697.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Zhang L, Chow EP, Wilson DP. Distributions and trends in sexual behaviors and HIV incidence among men who have sex with men in China. BMC Public Health. 2012;12:546.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Fu X, Lin P, Wang Y, Li J, Yu G, Li Y. Prevalence of human immunodeficiency virus and syphilis among MSM in Guangdong from 2009 to 2013. Chin Prev Med. 2014;15(5):169–72.

    Google Scholar 

  46. Statistics Bureau of Guangdong Province, Survey office of Ntional Bureau of Statistics in Guangdong. 2020 Guangdong statistical yearbook: 2019 permanent population of Guangdong. Beijing: China Statistics; 2020.

    Google Scholar 

  47. Sisson SA, Fan Y, Beaumont MA. Handbook of approximate bayesian computation. Boca Raton: CRC; 2019.

    Google Scholar 

  48. Marin J-M, Pudlo P, Robert CP, Ryder RJ. Approximate bayesian computational methods. Stat Comput. 2011;22(6):1167–80.

    Article  Google Scholar 

  49. Garnett GP, Mertz KJ, Finelli L, Levine WC, St Louis ME. The Transmission dynamics of Gonorrhoea: modelling the reported behaviour of infected patients from Newark, New Jersey. Philos Trans R Soc Lond B Biol Sci. 1999;354(1384):787–97.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Del Moral P, Doucet A, Jasra A. An adaptive sequential Monte Carlo method for approximate bayesian computation. Stat Comput. 2012;22(5):1009–20.

    Article  Google Scholar 

  51. Campolongo F, Cariboni J, Saltelli A. An effective screening design for sensitivity analysis of large models. Environ Model Softw. 2007;22(10):1509–18.

    Article  Google Scholar 

  52. Saltelli A, Annoni P, Azzini I, Campolongo F, Ratto M, Tarantola S. Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index. Comput Phys Commun. 2010;181(2):259–70.

    Article  CAS  Google Scholar 

  53. Important News. Guangzhou: Guangdong provincial Center for Disease Control and Prevention. http://cdcp.gd.gov.cn/.

  54. News of basic health service. Guangzhou: People’s Government of Guangdong Province. http://www.gd.gov.cn/zwgk/zdlyxxgkzl/ylws/content/post_4056940.html.

  55. Zhao Y, Han M, Ma Y, Li D. Progress towards the 90-90-90 targets for controlling HIV — China, 2018. China CDC Wkly. 2019;1(1):4–7.

    Article  Google Scholar 

  56. Sun C, Li J, Liu X, Zhang Z, Qiu T, Hu H, Wang Y, Fu G. HIV/AIDS late presentation and its associated factors in China from 2010 to 2020: a systematic review and meta-analysis. AIDS Res Ther. 2021;18(1):96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Zhang H. Sexual transmission as a dominant route and new diagnoses increasing among older people. Guangzhou: YangCheng Evening News. 2021.https://baijiahao.baidu.com/s?id=1717928090025091240&wfr=spider&for=pc.

  58. Guangdong Provincial CDC: Internal report. Guangzhou: Guangdong Provincial Center for Disease Control and Prevention (CDC). 2019.

    Google Scholar 

  59. Mirzaei H, Moradi Y, Abbaszadeh S, Nasiri N, Mehmandoost S, Khezri M, Tavakoli F, Sharifi H. The impact of COVID-19 on disruptions of HIV-related services: a rapid review. Med J Islam Repub Iran. 2022;36:98.

    PubMed  PubMed Central  Google Scholar 

  60. Zang X, Krebs E, Chen S, Piske M, Armstrong WS, Behrends CN, Del Rio C, Feaster DJ, Marshall BDL, Mehta SH, et al. The potential epidemiological impact of coronavirus disease 2019 (COVID-19) on the human immunodeficiency virus (HIV) epidemic and the cost-effectiveness of linked, opt-out HIV testing: a modeling study in 6 US cities. Clin Infect Dis. 2021;72(11):e828–34.

    Article  CAS  PubMed  Google Scholar 

  61. Shi L, Tang W, Hu H, Qiu T, Marley G, Liu X, Chen Y, Chen Y, Fu G. The impact of COVID-19 pandemic on HIV care continuum in Jiangsu, China. BMC Infect Dis. 2021;21(1):768.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Moitra E, Tao J, Olsen J, Shearer RD, Wood BR, Busch AM, LaPlante A, Baker JV, Chan PA. Impact of the COVID-19 pandemic on HIV testing rates across four geographically diverse urban centres in the United States: an observational study. Lancet Reg Health Am. 2022;7:100159.

    PubMed  Google Scholar 

  63. Pampati S, Emrick K, Siegler AJ, Jones J. Changes in sexual behavior, PrEP adherence, and access to sexual health services because of the COVID-19 pandemic among a cohort of PrEP-using MSM in the South. J Acquir Immune Defic Syndr. 2021;87(1):639–43.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Cardona-Arias JA, Higuita-Gutiérrez LF. Impact of COVID-19 pandemic on seroprevalence of HIV, HBV, HCV and HTLV I-II in a blood bank in Medellín, Colombia, 2019–2022. Trop Med Infect Dis. 2023;8(2):118.

  65. Zhang KC, Fang Y, Cao H, Chen H, Hu T, Chen YQ, Zhou X, Wang Z. The impacts of the COVID-19 pandemic on HIV testing utilization among men who have sex with men in China: cross-sectional online survey. JMIR Public Health Surveill. 2022;8(5):e30070.

    Article  PubMed  PubMed Central  Google Scholar 

  66. News. Guangzhou: Health Commission of Guangdong Province. http://wsjkw.gd.gov.cn/zwyw_gzdt/content/post_2171879.html.

Download references

Acknowledgements

We thank Prof. Zhang Jinxin for his valuable advice. We also thank the AIDS experts from Guangdong Provincial CDC for their heuristic viewpoints.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

RY conceived the study, developed the model, collected the parameter data, analyzed the sensitivity of the parameters, calibrated the model, prepared the results, and wrote the first draft of the manuscript. YSL gave valuable advice and revised the first manuscript draft. JG supervised the study and revised the first manuscript draft. All authors contributed to writing and approving the final manuscript.

Corresponding author

Correspondence to Jing Gu.

Ethics declarations

Ethics approval and consent to participate

Not applicable. Our data are mainly from peer-reviewed published articles, domestic government reports, AIDS expert interviews, etc. All methods in the study were carried out in accordance with relevant guidelines and regulations. Since there was almost involved in data of individuals for the study, informed consent and ethical approval were not required.

Consent for publication

Not applicable.

Competing interests

The authors have no actual or potential conflicts of interest related to this study to declare.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ye, R., Lai, Y. & Gu, J. New sexually transmitted HIV infections from 2016 to 2050 in Guangdong Province, China: a study based on a dynamic compartmental model. BMC Public Health 24, 1307 (2024). https://doi.org/10.1186/s12889-024-18735-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-024-18735-z

Keywords