Abstract
Introduction The initial phase of the coronavirus disease 2019 (COVID-19) pandemic caused significant interruptions to human immunodeficiency virus (HIV) prevention and healthcare services. In Australia, these interruptions coincided with a large reduction in HIV diagnoses but it is unknown if this reflects a reduction in incidence or barriers to testing. We aimed to understand the changes in HIV transmission among men who have sex with men in Australia and the impact of disruptions to HIV prevention, testing and treatment services due to the COVID-19 pandemic and the public health response by incorporating time-sensitive factors in real-world situations.
Methods We developed a mathematical model to estimate monthly HIV incidence between January 2020 and August 2022. We obtained aggregated monthly data for sexual partners, condom use, HIV testing, pre-exposure prophylaxis (PrEP) use, and migration. Three scenarios were simulated: 1) a COVID scenario with all changes in place; 2) a no COVID-19 scenario where input parameters remained at pre-COVID-19 values, and 3) a no COVID-19 scenario with continued PrEP scale-up.
Results In the absence of the COVID-19 pandemic, the estimated number of cumulative infections from January 2020 to August 2022 would have been 1,266 [95% Percentile Interval (PI): 1,100–1,466] compared to 915 [95% PI: 729–1,181] for the COVID-19 scenario (a 27.7% reduction). The largest reduction in infections (44.3%) occurred in 2020 with 273 (95% PI: 221– 344) infections versus 490 (95% PI: 435–552) in the no-COVID scenario. There was a rebound with 387 infections (95% PI: 307–502) in 2021, followed by a reduction to a stable level by August 2022. Our model identified reductions in sexual partners as the leading factor contributing to the change in HIV infections and diagnoses (−24.8% and −10.5%, respectively).
Conclusions A substantial reduction in new HIV infections and diagnoses in Australia occurred during the early stages of the COVID-19 pandemic, which was largely due to reduced HIV transmission. A rebound in infections as sexual partnerships increased between 2021-2022, highlights the imperative to maintain vigorous response efforts and take advantage of the gains made to virtually eliminate HIV transmission in Australia.
Introduction
The coronavirus disease 2019 (COVID-19) not only resulted in millions of confirmed cases and deaths worldwide [1] but also severely interrupted healthcare use [2]. This led to fewer clinic visits, admissions, tests and treatments [2], potentially increasing the burden of infectious diseases including human immunodeficiency virus (HIV), tuberculosis, and malaria [3].
In the early stages of the COVID-19 pandemic, Australian state and federal governments introduced border closure and lockdown measures to control the spread of COVID-19 [4, 5]. Physical restrictions including physical-distancing rules, stay-at-home orders, and closure of non-essential businesses were also introduced [4]. The COVID-19 vaccine rollout in Australia was launched in late February 2021 [6]. Most restrictions were lifted in December 2021 when more than 90% of people over the age of 16 were fully vaccinated [7], and the international borders were reopened to all fully vaccinated visa holders in February 2022 [8]. These measures significantly influenced sexual behaviours and the use of HIV services [4, 5]. Studies reported a substantial reduction in HIV tests during lockdowns and a dramatic decrease in sexual contacts during the first month of the national lockdown [4, 5]. These changes likely contributed to the large reduction (40%) in the number of HIV diagnoses in Australia, decreasing from 895 in 2019 to 541 in 2021 [9]. However, it is unknown if this reduction in HIV diagnoses is mainly related to a fall in testing due to the decreased access to HIV services [4] or a fall in new infections due to reduced sexual contacts [5, 10]. Another factor could have been border closures and the reduced number of migrant arrivals from overseas. Furthermore, the pandemic affected the use of preventive measures such as pre-exposure prophylaxis (PrEP) and condoms. A previous study in Australia showed that a large proportion (42%) of men who have sex with men (MSM) discontinued PrEP use in 2020 [11]. Similar impacts were reported in other high-income settings with sexual minority men in the United States reporting decreased condomless anal sex with casual partners during the early stages of the COVID-19 pandemic [12]. However, the subsequent influence of these changes on the HIV epidemic remains unclear. To better understand the impacts of the COVID-19 pandemic on the HIV epidemic, an evaluation incorporating all factors related to HIV transmission affected by COVID-19 is required.
The impacts of COVID-19 on the HIV epidemic [13–16] have been assessed through modelling. However, many studies were published in the early stage of the pandemic and relied on assumed impacts on factors related to HIV transmission, rather than actual directly observed impacts. Most modelling studies assumed constant impacts and often overlooked the dynamic impacts of specific policies such as lockdowns and border closures on HIV transmission [13–16]. This highlights the need for more nuanced studies that examine the dynamics of the HIV epidemic in response to the COVID-19 response.
MSM are one of the key populations affected by HIV [17]. In Australia, male-to-male sex is the primary HIV risk exposure, accounting for 70% of new HIV infections between 2015 and 2019 [9] and most HIV infections among MSM are attributable to sex with non-romantic partners (casual sexual partners and friends with benefits) [18]. Detailed datasets are available for this population in Australia. This study aimed to use a mathematical model to evaluate the impact of COVID-19 restrictions on the HIV epidemic among MSM in Australia, considering the empirically recorded monthly changes in sexual behaviours, HIV testing, and migration during the acute phase of the COVID-19 pandemic.
Methods
We first developed a simple but general and flexible HIV transmission model which we then applied to the MSM population in Australia. This Simple HIV Model is written in R version 4.3.2. and is publicly available under an open access license [19]. No ethical approval or consent was required as the study only used publicly available data.
Simple HIV Model
Our Simple HIV Model expands previously developed models [20–22]. The model uses difference equations to describe the HIV cascade and the effect of changes in sexual partners, condom use, HIV testing, and PrEP use on HIV transmission. At a pre-specified timestep t, Nt is the number of people living with HIV in the population, It is the number of new infections, and βt is the transmission probability from a person with an unsuppressed HIV viral load to a susceptible partner not on PrEP. Therefore, if no interventions were implemented, HIV-positive individuals would transmit Ntβt infections at timestep t. To represent the steps of the HIV cascade, we defined dt to be the proportion of people living with HIV who are diagnosed, τt to be the proportion of those diagnosed on antiretroviral therapy (ART), σt to be the proportion of those on ART with suppressed viral load (< 200 copies per ml). If ϕ is the reduction in transmission due to viral suppression then HIV-positive individuals transmit Ntβt(1 − dtτtσtϕ) infections at timestep t. To account for PrEP use, we defined ωt to be the proportion of the susceptible population on PrEP, and ∊p to be the efficacy of PrEP in preventing infection. The number of infections at time step t then becomes Ntβt(1 − dtτtσtϕ)(1 − ωt∊p). With these definitions, a simple risk equation was constructed to estimate the number of new infections due to transmission to susceptible sexual partners:
To account for changes in the population due to migration and mortality, we defined Mt to be the number of people who enter Australia living with HIV, Et to be the number of people who leave Australia living with HIV, and μt to be the probability of all-cause mortality in people living with HIV. Then the number of people living with HIV in the population at timestep t+1 would be (Nt + It + Mt − Et)(1 − μt). These definitions allowed us to construct Equation 2:
Before calculating Nt and It, an initial value for the transmission probability (β0) at time t = 0 is required. This is obtained using initial values for the population size, new infections, and each parameter in the following equation:
Effect of changing sexual behaviours
To consider the impacts of changes in sexual behaviours, we incorporated the effects of changes in condom use and sexual partners. If the proportion of condom use changes to ct at timestep t compared to an initial level c0 with an effectiveness of ∊c, the relative change in the transmission probability is given by . If the number of sexual partners changes to pt of the time compared to a pre-covid level of p0, the relative change in the transmission probability is given by
.
By applying these relative changes to the initial transmission probability, we obtain:
Effect of changes in testing rates
The proportion diagnosed over time is affected by changes in testing rates. If the proportion of the population tested at time-step t is Tt, the number of people living with HIV diagnosed at the next time step is given by:
With Equation 2, this can be used to calculate the proportion diagnosed at the next time step:
We only included diagnosed immigrants (Mt) and emigrants (Et) in the model assuming the number undiagnosed is small compared to the overall population and has a minimal impact on the results.
Defining the number of new diagnoses at the initial timestep to be D0, we solved Equation 7 for the proportion of the population tested at the initial timestep (T0) before calculating dt+1:
Application of model
For this study, we focused on the MSM population in Australia and used a monthly timestep starting from December 2019. We used notifications from national surveillance as a proxy for the initial number of diagnosed cases (D0) [9]. The changes in each parameter were captured over time following the implementation of restrictions and border closures in response to the COVID-19 outbreak in March/April 2020 [23]. Up to that time we assumed input parameters remained at their December 2019 values. To estimate the impact of the COVID-19 pandemic on sexual partners, we used the reported number of non-romantic partners (casual sexual partners and friends with benefits) as a proxy for the number of sexual partners at both time 0 (p0) and time t (pt). To estimate the impact of the COVID-19 pandemic on condom use, we used the reported proportion of MSM with casual partners reporting no condomless anal intercourse in the past four weeks as a proxy for the proportion of MSM reporting condom use at both time 0 (c0) and time t (ct). Constant parameters are presented in Table 1. The input variables ωt, Tt, ct, and pt are time varying and were estimated using data from Kirby Institute’s monitoring of PrEP uptake in Australia, the Flux cohort study, and the Gay Community Periodic Survey [5, 24–26] (details provided in the Supplementary Material).
Scenarios
We simulated three scenarios from December 2019 to August 2022. The first scenario was the “COVID-19 scenario”, where all changes in monthly parameter values occur and reflect the impact of COVID-19. The second scenario was the “no COVID-19 scenario”, a counterfactual scenario where the input parameters remained at their values in December 2019 and corresponded to pre-COVID-19 values. The third scenario was the “no COVID-19 plus PrEP scenario”, an alternative counterfactual scenario where PrEP scale-up and the corresponding decrease in condom use prior to 2019 continued during 2020-2022, with the other parameters remaining at their December 2019 value. In this scenario, we assumed that the proportion of MSM on PrEP increased from 36.5% to 58.2% over this period following trends for the number of people dispensed PrEP [26] and the proportion using condoms from 73.1% to 64.9% following the trends reported in the Gay Community Periodic Survey [25] (Table S1). Further details of each scenario are provided in the Supplementary Material.
Further simulations were conducted by separately modelling the impact of COVID-19 on each time series variable to assess their individual effects on the HIV epidemic. Migration was considered as a single variable incorporating both MSM living with HIV who enter (Mt) and leave (Et) Australia.
Model calibration, validation, and uncertainty analysis
We calibrated the transmission probability at the initial timestep (β0) to fit the average monthly number of new HIV infections among MSM in 2019 using Equation 3 [9]. We then calibrated the testing rate at the initial timestep (T0) to fit both the average monthly number of new HIV infections among MSM in 2019 and the average monthly number of MSM diagnosed with HIV in 2019 using Equation 7 [9]. The annual numbers of MSM diagnosed with HIV from 2020 to 2022 were then used to validate the results of the model. To account for uncertainty, all three scenarios were simulated 1000 times by randomly sampling from each parameter range. We assumed that the uncertainty range for parameters including PrEP use (ωt), sexual partners (pt), condom use (ct), and testing rate (Tt) at each timestep were ±5%. The results from each simulation were then used to derive a 95% percentile interval (95% PI) for the model outputs.
Results
After calibration, the estimated number of annual diagnoses in the COVID-19 scenario closely aligned with data from national surveillance between 2020-2022 (Figure 1), providing a validation for our model estimates.
The estimated diagnoses from the model in 2022 were adjusted to reflect the annual number.
Number of new infections, diagnoses, and MSM with HIV in each scenario
As shown in Table 2 and Figure 2, in the absence of the COVID-19 pandemic (no COVID-19 scenario), the estimated number of cumulative infections from 2020 to 2022 would have been the highest, at 1,266 [95% PI: 1,100–1,466]. Figure 2 shows that there was a stable level of monthly new infections with a slightly decreasing trend. In contrast, there was a 27.7% reduction in the number of cumulative infections for the COVID-19 scenario (915, [95% PI: 729–1,181]) (Table 2). The largest reduction in infections (44.3%) occurred in 2020 with 273 (95% PI: 221–344) infections versus 490 (95% PI: 435–552) for the no COVID-19 scenario. There was a rebound with 387 infections (95% PI: 307–502) in 2021, coinciding with a rebound in sexual partners (see Supplementary Material Figure S1), followed by a reduction to a stable level by August 2022 (Figure 2). An increase in PrEP use (no COVID-19 plus PrEP scenario) would have resulted in a slight reduction in infections compared to the no COVID-19 scenario, but not as substantial as observed in the COVID-19 scenario. The estimated number of cumulative infections from 2020 to 2022 in the no COVID-19 plus PrEP scenario was 1,151 (95% PI: 998–1,339), with a 9.1% reduction compared to the no COVID-19 scenario (Table 2). The monthly difference in the number of new infections between the no COVID-19 scenario and the no COVID-19 plus PrEP scenario expanded over time, with the largest difference in infections (21.1%) occurring in August 2022 (Figure 2). The uncertainty ranges for the three scenarios started to overlap in April 2021.
Bands indicate 95% percentile intervals.
Similarly, in the no COVID-19 scenario, the number of cumulative diagnoses over the simulation period would have been 1,613 (95% PI: 1,041–2,221) compared to 1,301 (95% PI: 758–1,870) for the COVID-19 scenario (a 19.3% reduction). The projected number of cumulative diagnoses in no COVID-19 plus PrEP scenario was intermediate to the estimates of the other two scenarios, at 1,579 (95% PI: 1,006–2,189) (Table 2). The number of MSM with HIV would have been slightly higher in the no COVID-19 scenario than that in the COVID-19 scenario, increasing from 21,273 (95% PI: 20,285–22,281) to 21,724 (95% PI: 20,757–22,725) by August 2022 (Table 2). Whereas in the no COVID-19 plus PrEP scenario the number of MSM with HIV was intermediate to the estimates of the other two scenarios, at 21,617 (95% PI: 20,664–22,623) (Table 2). Further details are provided in the Supplementary Material.
Impacts of each time series variable on the HIV epidemic
Our model identified sexual partners as the biggest driver of the changes in new HIV infections (Table S3, Figure S2), with a reduction of 24.8% over the simulation period compared to the no COVID-19 scenario. The other individual variables (HIV testing, condom use, PrEP use, and migration) had smaller effects on infections (2.5% increase, 1.7% reduction, 2.5% increase, and 0.1% reduction, respectively). In terms of HIV diagnoses, the substantial decrease was primarily due to a combination of the reduction in sexual partners (10.5% reduction) and HIV testing (7.5% reduction) over the study period. Further details are provided in the Supplementary Material.
Discussion
Our study used a mathematical model to evaluate the impact of COVID-19 on the HIV epidemic among MSM in Australia. Our simulations estimated that COVID-19-related interruptions to HIV testing, PrEP use, and migration as well as the associated impact on partner numbers and condom use resulted in a substantial decrease in both diagnoses (19.3%) and new HIV infections (27.7%). The largest reduction in new infections (44.3%) occurred in 2020, followed by a rebound in 2021 and then a reduction to a stable level in 2022, at a similar level to what might happened if COVID-19 did not occur, which suggests that the direct impact of COVID-19 on HIV transmission dynamics may have diminished over time. Our study provides new insights into the impacts of COVID-19 on the HIV epidemic by capturing monthly changes in time-sensitive factors in real-world situations including sexual partners, HIV testing, condom use, PrEP use, and migration, compared to previous studies that relied on assumed data or constant impacts [13–16].
Our study identified reductions in sexual partners, proxied by the number of non-romantic partners, and HIV testing as the primary drivers behind the observed decline in HIV diagnoses. This finding highlights the importance of maintaining HIV testing services when public health restrictions are in place in response to a severe pandemic. Telehealth-delivered, web-based, and home-based HIV testing services could be promoted to ensure continuing HIV testing [27–29]. Also, national surveillance reported an increasing trend in the proportion of late diagnoses among MSM between 2019 and 2022, indicating a COVID-19 lag effect on the HIV epidemic [9].
We also found that physical distancing restrictions and the resulting reduction in sexual partners had the biggest impact on new infections, with a 24.8% decrease alone compared to a 27.7% reduction when all changes were included, counteracting the increased risk of infection due to disrupted HIV testing and PrEP use. Our findings are consistent with previous modelling. A study from the United States found that compared to condom use, PrEP use, and HIV testing, the sexual partner numbers were the biggest driver of the change in new HIV infections and a theoretical 50% decrease in the sexual partner numbers could reduce the number of new HIV infections by 24% [14]. Another modelling study showed that the reduced number of sexual partners could offset the impact of HIV service disruptions on new HIV infections [30]. However, these studies only used theoretical scenarios and did not consider empirically observed monthly changes. Compared to these studies, our study incorporated data on the actual dynamics of these time-sensitive factors in real-world situations.
This study projected that continued PrEP scale-up would have slightly reduced the number of new HIV infections with an intensified impact over time, which is consistent with findings from previous studies [31, 32]. This finding highlights the importance of the continuity and expansion of PrEP use. Despite the COVID-19 pandemic disrupting the increasing trend, PrEP use had fully recovered by October 2021 and has continued to increase steadily until the latest available data in June 2023 [33].
There are several limitations in our study. First, results from the current study should be interpreted with caution due to the simplicity of the model. However, our validation results indicate that our model estimates the number of new HIV diagnoses over the entire simulation period well. Also, the robustness of the time series data for variables such as PrEP use, HIV testing, condom use, and sexual partners supported the reliability of the results from our model. Second, we did not have monthly data on ART use and viral suppression to incorporate into the model. Several studies have highlighted the significant impacts that the changes in ART use and viral suppression during the COVID-19 pandemic can have on the HIV epidemic [14–16]. However, previous studies in Australia showed that the proportion of viral suppression and access to ART essentially continued uninterrupted, even during the lockdown restrictions [34, 35]. This finding suggests that the uncertainty introduced by the lack of monthly data on ART use and viral suppression is likely minimal. Third, our results may not fully capture the impacts of key variables on the HIV epidemic during March and April 2020, as the Flux study only started in May 2020. We assumed that the value of these variables in April 2020 was the same as that in May 2020. Given that restrictions and the border closure were introduced in mid-March 2020, we assumed the impact on these variables for March was half of that in April. Our model may underestimate the impact of restrictions on these variables in March and April 2020, as some restrictions were eased in May 2020 [36]. Fourth, we considered the overall impact of lockdown measures on the HIV epidemic. However, these measures varied significantly between states and territories. These regional differences in public health responses could have had varying impacts on factors related to HIV transmission, potentially influencing our results. Future analyses could explore these jurisdictional variations to better understand the HIV epidemic at the state and territory level during the pandemic. Fifth, potential discrepancies arise because we used the proportion of MSM with casual partners reporting no condomless anal intercourse in the past four weeks as a proxy for condom use and the number of sexual partners to represent overall sexual activity rather than consider sexual acts, as a proxy for sexual activity. Finally, in our model, we assumed that changes in condom use and PrEP use are independent. This assumption may not fully reflect real-world behaviours, as those taking up PrEP are also the ones who are reducing condom use. By applying relative risks independently, the model may overestimate HIV transmission risk and underestimate the impact of PrEP use.
Our study is the first study to incorporate the time varying effects of COVID-19 restrictions on the HIV epidemic among MSM using real-world data. Our study demonstrates that COVID-19 restrictions contributed to a genuine reduction in new HIV infections in Australia. The generalisability and flexibility of the Simple HIV model means it is adaptable to other settings or populations.
Conclusions
A substantial reduction in new HIV infections and diagnoses in Australia occurred during the early stages of the COVID-19 pandemic, which was largely due to reduced HIV sexual partnerships. A rebound in infections between 2021-2022, highlights the imperative to maintain vigorous response efforts and take advantage of the gains made to virtually eliminate HIV transmission in Australia. This outcome shows how changes in time-sensitive factors in response to large-scale public health policies or interventions can affect HIV epidemics. By understanding this relationship, public health authorities can better prepare for and mitigate the effects of similar disruptions on the HIV epidemic in the future.
Competing interests
None declared.
Authors’ contributions
RW and RTG conceived the study. RW and RTG developed the model. RW set up and ran the scenarios. RW drafted the manuscript. JAK, MH, SM, and RTG validated the model input. All authors were involved in writing and revising the manuscript.
Funding
RW is supported by an Australian Government Research Training Program Scholarship. The Kirby Institute is funded by the Australian Government Department of Health and is affiliated with the Faculty of Medicine, UNSW Sydney, Australia.
Data Availability Statement
All data produced in the present work are contained in the manuscript and Supplementary Materials. The Simple HIV Model code used to produce results for this study along with model input parameters and summary results are available online at https://github.com/The-Kirby-Institute/Simple_HIV_model/
Supporting Information
Supporting Information file 1: Supplementary Material
docx. Additional methodological details and results are provided in the file.
Acknowledgements
We would like to thank Dr Doug Fraser for the information on the PrEP uptake in Australia.
Footnotes
E-mail addresses of authors: RW: rweng{at}kirby.unsw.edu.au, J.A K: akwon{at}kirby.unsw.edu.au, MH: mhammoud{at}kirby.unsw.edu.au, BC: brent{at}napwha.org.au, NS: nick.scott{at}burnet.edu.au, SM: smcgregor{at}kirby.unsw.edu.au, R.T G: rgray{at}kirby.unsw.edu.au