Estimating the burden of mpox among MSM in South Africa
=======================================================

* Ruth McCabe
* Leigh F. Johnson
* Lilith Whittles

## Abstract

Despite seeing few cases during the 2022-23 mpox global outbreak, recent reports of 22 cases among men-who-have-sex-with-men (MSM) in South Africa, including three deaths, have raised concerns about underreported community transmission. We used a Monte Carlo simulation model to estimate the true epidemic size, considering increased severity of mpox among MSM living with advanced HIV (MSMLAHIV), documented overrepresentation of people living with HIV (PLHIV) among mpox cases, and HIV prevalence in South Africa. We estimate there have been between 220-450 cases among MSM LHIV in South Africa, implying a total of 290-560 cases among all MSM. We provide an upper bound of 750-1,600 cases as a sensitivity analysis where the prevalence of HIV among mpox patients is the same as population prevalence among MSM in South Africa. Estimates in both scenarios suggest a substantial number of undetected cases, with case ascertainment rates estimated between 1% and 8%. Our findings underscore the need for enhanced surveillance, targeted public health interventions, and awareness campaigns to mitigate the outbreak’s impact at a population-level.

## Introduction

Since 2022, the global outbreak of mpox has caused over 97,000 cases across 116 countries (1). In contrast to outbreaks in historically endemic countries, transmission mostly occurred via sexual networks, with cases heavily concentrated among men who have sex with men (MSM). Prior to 2024, South Africa was largely unaffected by the global outbreak of mpox, recording only five cases, all of which were linked to international travel (2). However, between May and July 2024 the country has reported 22 cases and three deaths (2). All 22 cases occurred in men aged between 17 and 43 (3), most of whom self-identify as MSM (4). Genomic sequencing was performed on a subset of five cases, all of which were confirmed to be MPXV subclade IIb – linking them to the 2022-23 global outbreak, rather than the ongoing epidemic of MPXV clade I in the Democratic Republic of Congo and bordering countries (5). The global outbreak of MPXV clade IIb has been associated with a low case fatality ratio (CFR) 0.21% (n=203/97,745; 95% CI: 0.15% - 0.24%), so the three recent deaths in South Africa have prompted concern about the true burden of mpox. Furthermore, the recent cases had no history of international travel or attendance at high-risk gatherings, and are largely without epidemiological links, fuelling concern of sustained, undetected community transmission (4,6).

Globally, PLHIV have been disproportionately affected by mpox, and are over-represented among both cases and deaths (7). Increased mpox severity and fatality rates among people living with advanced HIV (PLAHIV) have been well-documented (8), which is of particular concern due to typically higher prevalence of HIV among MSM compared to the general population. South Africa is among the countries most affected by the HIV worldwide, with an estimated 18.3% (95% CI: 15.6%-20.5%) of the adult population living with HIV, which rises to 29.4% (25.5% - 33.5%) in MSM (9,10). At least two-thirds of the recent South African mpox cases (n=15/22) occurred in people living with advanced HIV – including the individuals with fatal outcomes (6).

Here, by leveraging the relationship between deaths and cases via the CFR, we estimate the potential total number of mpox cases underlying the observed number of cases and deaths thus far (July 2024), accounting for the high prevalence of HIV in South Africa and its effect of mpox severity.

## Methods

To estimate the total number of mpox cases in MSM, *C*, given observed deaths, all of which were in MSMLAHIV (*D* = *D**AHIV*), we first calculated the expected number of cases among MSMLAHIV, *C**AHIV*, using the corresponding observed mpox CFR among PLAHIV, *P* (*D**AHIV* |*C**AHIV*): ![Formula][1]</img>  

We then used this estimate to calculate the expected number of cases among all MSMLHIV, *C**HIV*, by considering the proportion of mpox cases with HIV whose infection is in an advanced stage: ![Formula][2]</img>  

Finally, we estimated the total number of mpox cases among MSM, using the prevalence of HIV among mpox cases: ![Formula][3]</img>  

### Model parametrisation

#### Case fatality ratio in MSMLAHIV

A global case series of mpox among PLHIV found clear evidence of higher CFRs in people with lower CD4+ T-cell counts (8). Of the 27 mpox deaths recorded among a sample of 382 cases with HIV coinfection during the 2022-23 global outbreak, 23 occurred in the 85 individuals CD4+ T-cell count <100/mm3, corresponding to a CFR of 27.1% (95% exact binomial confidence interval (CI): 18.0% -37.8%). The remaining four deaths were recorded in the 94 individuals with CD4+ T-cell counts between 100-200/mm3, equating to a CFR of 15.1% (95% CI: 10.2% - 21.2%) in those with CD4+ T-cell count <200/mm3. No deaths occurred in the 204/382 mpox cases with CD4+ T-cell count ≥200/mm3. The definition of ‘advanced HIV’ used in the South African mpox case reports is unclear but is typically defined as CD4+ T-cell count <200/mm3, which we take as our central value. However, we perform sensitivity analysis using an alternative definition based on CD4+ T-cell count <100/mm3.

#### Proportion of mpox cases LHIV with advanced HIV

We assume that the proportion of MSMLHIV whose HIV infection is at an advanced stage is the same among mpox cases as in the MSM population, specifically, that the probability of having advanced HIV can be considered independent of the probability of having mpox, conditional on HIV infection. To quantify the proportion of MSM LHIV whose infection is in an advanced stage, we use estimates from the Thembisa 4.7, a transmission-dynamic deterministic compartmental model of HIV in South Africa (9), which has been used to produce HIV estimates for national policy development including those reported by UNAIDS (11). The model predicts that 6.6% of adults living with HIV in South Africa in 2024 will have a CD4+ T-cell count below 200/mm3.

For our sensitivity analysis of the definition of ‘advanced HIV’, we assume that 3.3% of the population will have CD4+ T-cell count <100/mm3, (i.e. that CD4+ T-cell counts are distributed roughly uniformly below 200/mm3). This assumption aligns with the global case series of mpox in PLAHIV, in which mpox cases with CD4+ T-cell counts <100/mm3 comprised 47.4% (n=84/179; 95% CI 40.0% - 55.1%) of mpox cases with CD4+ T-cell counts<200/mm3 (8).

#### Proportion of mpox cases LHIV

Globally, PLHIV have observed a higher incidence of mpox cases, likely due to shared sexual network risk (7). A conservative lower bound for the prevalence of HIV among mpox cases would therefore be the prevalence of HIV among MSM, *P* (*N**HIV*), which the Thembisa model estimates at 29.5% (95% CI: 25.5% - 33.5%) in South Africa. However, this method may overstate the overall number of cases and so we consider this as a sensitivity analysis. Instead, our central scenario considers that the odds of LHIV given mpox infection are proportional to the odds of LHIV. ![Formula][4]</img>  

We estimated the odds ratio, *α*, using mpox case data from the USA in which HIV prevalence was 56.4% (n=5,350/9,492; 95% CI: 55.4% - 57.4%) (12), compared to prevalence of 12.4% (95% CI: 10.9% - 13.8%) among all MSM (13). Serological evidence has shown a low prevalence of undiagnosed mpox in the US, which supports the use of the data to form an unbiased assessment of risk-patterns (14). These data suggest the odds of having HIV were *α* = 9.2 (95% CI: 7.9 – 10.6) times higher among mpox cases than in the general MSM population.

We multiplied the odds of HIV infection among all MSM by *α* to obtain a central estimate of the prevalence of HIV among mpox cases in South Africa: ![Formula][5]</img>  

#### Modelled scenarios

We created a Monte Carlo simulation model to quantify the potential number of mpox cases in South Africa i) in MSMLAHIV; ii) in all MSMLHIV; and iii) in MSM in total. Conditional probabilities were modelled using beta distributions, parameterised from the literature, either directly from reported sample numbers where possible, or estimated via method-of-moments from mean and 95% confidence intervals (Table 1). We generated 10,000 samples from parametric distributions to propagate the statistical uncertainty inherent in the parameter estimates throughout the calculations, allowing us to present results in terms of mean and 95% confidence intervals. Analyses were performed in R v4.2.2; model code is available at [https://github.com/mrc-ide/mpox_sa](https://github.com/mrc-ide/mpox_sa).

View this table:
[Table 1:](http://medrxiv.org/content/early/2024/08/14/2024.08.13.24311919/T1)

Table 1: Model parameters, descriptions and sources

## Results

Table 2 summarises the results of the central scenario and associated sensitivity analyses.

View this table:
[Table 2:](http://medrxiv.org/content/early/2024/08/14/2024.08.13.24311919/T2)

Table 2: 
Estimated potential number of mpox cases in South Africa among MSMLAHIV, MSMLHIV, and in total under four scenarios with corresponding parameter assumptions. Results are presented as mean (95% confidence interval) where applicable

### Estimated mpox cases among MSM LHIV

The three mpox deaths observed among PLAHIV in South Africa are consistent with an underlying burden of 21 cases (95% CI: 15 – 29) (Table 2). The 15 cases that have been observed in PLAHIV up to June 2024 are therefore broadly consistent with modelled expectations, subject to case under-ascertainment. Taking the observed cases as a lower bound of the true total implies a case-ascertainment rate of 74.9% (95% CI: 51.2% - 100%) among MSMLAHIV. Under the assumption that 6.6% of HIV-positive MSM with mpox have CD4+ T-cell counts < 200/mm3, we estimate total cases in PLHIV to be 313 (95% CI: 227 – 444).

### Estimated total mpox cases in MSM

Allowing for overrepresentation of PLHIV among mpox cases, in line with relative odds observed in the US, suggests an HIV prevalence of 79.1% (95% CI: 74.9% - 82.9%) among MSM with mpox in South Africa, compared to 29.4% (95% CI: 25.5% - 33.5%) in all MSM. Taking this population prevalence of HIV into account, we estimate the total burden of mpox at 395 (95% CI: 287 - 562) cases. These estimates correspond to CFR at the high end of that observed in the global outbreak: 0.5%-1.0%, reflecting the high burden of HIV in South Africa, and an mpox case ascertainment of around 3.9%-7.7% to date.

### Sensitivity analyses

Using the alternative definition of ‘advanced HIV’ as CD4+ T-cell count <100/mm3, the observed three deaths are compatible with 11 cases (95% CI: 8 - 17) among MSMLAHIV, which would imply that case ascertainment among MSMLAHIV has been extremely high, with 99.4% (95% CI: 90.4% - 100%) of infections in this group already observed. Taking the 15 observed cases in PLAHIV as a lower bound for the true number of cases, the assumption that 3.3% of MSMLHIV have CD4+ T-cell counts <100/mm3 is compatible with a higher number of cases than the CD4+ T-cell count <200/mm3 definition, suggesting an underlying case burden of 458 (95% CI: 455 – 503) among MSMLHIV and 579 (95% CI 549 - 640) in total.

In the conservative scenario, where we assume that prevalence of HIV among mpox cases is in line with background HIV prevalence among MSM in South Africa at 29.4% (95% CI: 25.5% - 33.5%), we estimate the total number of cases in the outbreak to be 1,068 (95% CI: 750 - 1,543) if the reported deaths had CD4+ T-cell counts <200/mm3, or up to 1,565 (95% CI: 1,361 - 1,832) under the stricter definition of CD4+ T-cell counts <100/mm3. These figures are 2-3-fold higher than in the central scenario and are consistent with an overall infection fatality ratio of around 0.2-0.3%, similar to that observed globally, and suggest that only 1%-3% of recent cases in South Africa have been detected.

## Discussion

Using recent reports of 15 mpox cases and three deaths among MSMLAHIV in South Africa, we used a simple probabilistic model to estimate the underlying burden of unobserved cases, both in total and among MSMLHIV.

A key advantage of our pragmatic approach is that we extrapolate from published literature on mpox epidemics in similar populations, as well as incorporating context-specific factors, such as the high prevalence of HIV among MSM in South Africa. A further strength is that we incorporate parametric uncertainty via the use of Monte Carlo simulations, allowing us to present a comprehensive view of the range of outcomes compatible with the observed number of cases and deaths. Finally, we present plausible upper and lower bounds for the true burden of mpox in South Africa by conducting a range of sensitivity analyses; we use different definitions of advanced HIV based on CD4 count thresholds, and present both cautious and neutral assumptions about the prevalence of HIV among mpox cases relative to the general population.

To enable estimation of the mpox case burden, we have made simplifying assumptions that should be borne in mind when interpreting the results. The uncertainty associated with each of these assumptions has been propagated through to our final estimates, which explains the relatively large confidence intervals that we report here. Firstly, the accuracy of our results is dependent on any biases in published parameter estimates; for example, if many cases among PLAHIV have been undetected, then the CFR could be lower than estimated in the global case series (8), which would serve to increase estimates of case numbers presented here. Concerns around biased estimates of severity in PLAHIV are mitigated by the fact that mpox symptoms are also more severe in this population, suggesting ascertainment is likely high (8). However, it could also be the case that the true number of deaths from mpox is under-reported due to ambiguity in the documentation of causes of death, which may suggest the observed CFR is an under-estimate. Secondly, given the case patterns observed globally, as well as recently in South Africa, our study limits the population at risk to MSM, and does not consider transmission to other populations, such as female partners. Consequently, our estimates may be an underestimate of the full burden of mpox in South Africa. Thirdly, we assume the severity of clade IIb mpox is unchanged since global outbreak, and that national differences in the observed CFR are attributable to the underlying prevalence of AHIV among MSM. Fourth, we assume that prevalence of AHIV is the same in mpox patients as in among all MSM LWHIV – however individuals with advanced HIV may be less likely to engage in sex, and so be less likely to contract mpox, compared to other HIV-infected individuals, which would also serve to increase estimates of the underlying number of cases. Moreover, recent studies have suggested that the percentage of men LAHIV in Sub-Saharan Africa could be as high as 13% (15). Adopting this higher estimate in our model would approximately halve the expected number of mpox cases. Conversely, LAHIV may make contracting mpox more likely if sexual contact with an infected individual does occur, due to inhibited mucosal immunity. Finally, we note that while information on the HIV status on the initial 20 cases was reported, this was not mentioned for the most recent two cases. We assume that these individuals are not LAHIV but note that if this were true then the estimated number of mpox cases among MSM would increase.

To place our findings in the broader context of global mpox outbreaks among MSM, in countries where epidemics have risen and peaked, likely due to a combination of network saturation effects and behavioural adaptation, the final epidemic size has been estimated to range from 0.1% to 1.5% of the MSM population (16), noting that determining accurate population denominators for MSM can be challenging. In South Africa, the MSM population is estimated at approximately 310,000 individuals (10), suggesting the potential final size of the mpox outbreak could reach up to 4,650 cases, or 3-10x the current estimated size. This underscores the urgency for targeted public health interventions, surveillance and risk communication and communication engagement (RCCE) to manage the outbreak effectively and mitigate its impact.

## Conclusions

A combination of evidence, including 22 mpox cases with no reported travel history, no epidemiological links, and no attendance at high-risk events/venues, coupled with three observed deaths, despite the low case fatality reported globally, suggests undetected community transmission of mpox in South African MSM. We estimate a current epidemic size of between 280-570 cases in our central scenario, or 750-1,600 cases in the conservative scenario assuming PLHIV are not overrepresented among mpox cases – corresponding to a detection rate of between 1 in 14 and 1 in 100 cases. The true size of the epidemic could be confirmed with serological studies, but in the meantime, awareness campaigns to encourage testing, symptomatic abstention from sex, and vaccination where available are recommended to mitigate the risks of future transmission.

## Data Availability

Analyses were performed in R v4.2.2; model code is available at https://github.com/mrc-ide/mpox_sa. All data are publicly available and cited in the manuscript.

## Contributions

LKW conceived the study, designed the methodology, and performed the analysis. LKW and RM drafted the report. LFJ advised on parameter values and revised the report.

## Declaration of interests

All authors declare no competing interests.

## Funding

LKW acknowledges funding from the Wellcome Trust (grant number 218669/Z/19/Z). RM and LKW acknowledge funding from the Medical Research Council (MRC) Centre for Global Infectious Disease Analysis (MR/X020258/1) funded by the UK MRC and carried out in the frame of the Global Health EDCTP3 Joint Undertaking supported by the EU. The views expressed are those of the authors and not necessarily those of MRC or Wellcome.

## License

For the purpose of open access, the author has applied a ‘Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising.

## Acknowledgments

LKW would like to acknowledge Dr Alyssa Bilinski’s ‘Napkin Math’ approach, which greatly guided her thinking.

*   Received August 13, 2024.
*   Revision received August 13, 2024.
*   Accepted August 14, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.World Health Organization. 2022-23 Mpox (Monkeypox) Outbreak: Global Trends [Internet]. [cited 2023 Sep 14]. Available from: [https://worldhealthorg.shinyapps.io/mpx\_global/#2\_Global\_situation\_update](https://worldhealthorg.shinyapps.io/mpx\_global/#2_Global_situation_update)
    
    

2.  2.Mahase E. Mpox: South Africa reports 20 new cases and three deaths. BMJ. 2024 Jul 10;386:q1529.
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjE3OiIzODYvanVsMTBfNS9xMTUyOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzE0LzIwMjQuMDguMTMuMjQzMTE5MTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

3.  3.Africa Centres for Disease Control and Prevention (Africa CDC). Mpox situation in Africa [Internet]. 2024 Jul. Available from: [https://www.africacdc.org](https://www.africacdc.org)
    
    

4.  4.Mpox Outbreak: Current Status and Ongoing Response Efforts(11 July 2024) [Internet]. NICD. 2024 [cited 2024 Aug 5]. Available from: [https://www.nicd.ac.za/mpox-outbreak-current-status-and-ongoing-response-efforts/](https://www.nicd.ac.za/mpox-outbreak-current-status-and-ongoing-response-efforts/)
    
    

5.  5.Deadlier strain of mpox spreads to multiple African countries [Internet]. [cited 2024 Aug 5]. Available from: [https://www.science.org/content/article/deadlier-strain-mpox-spreads-more-african-countries](https://www.science.org/content/article/deadlier-strain-mpox-spreads-more-african-countries)
    
    

6.  6.Disease Outbreak News: Mpox - South Africa (09 July 2024) - South Africa | ReliefWeb [Internet]. 2024 [cited 2024 Jul 22]. Available from: [https://reliefweb.int/report/south-africa/disease-outbreak-news-mpox-south-africa-09-july-2024](https://reliefweb.int/report/south-africa/disease-outbreak-news-mpox-south-africa-09-july-2024)
    
    

7.  7.Mitjà O, Ogoina D, Titanji BK, Galvan C, Muyembe JJ, Marks M, et al. Monkeypox. Lancet Lond Engl. 2023 Jan 7;401(10370):60–74.
    
    

8.  8.Mitjà O, Alemany A, Marks M, Mora JIL, Rodríguez-Aldama JC, Silva MST, et al. Mpox in people with advanced HIV infection: a global case series. The Lancet. 2023 Mar 18;401(10380):939–49.
    
    

9.  9.Johnson L, Dorrington R. Thembisa version 4.7: A model for evaluating the impact of HIV/AIDS in South Africa.
    
    

10. 10.Joint United Nations Programme on HIV/AIDS. UNAIDS DATA 2022. 2022.
    
    

11. 11.Johnson LF, May MT, Dorrington RE, Cornell M, Boulle A, Egger M, et al. Estimating the impact of antiretroviral treatment on adult mortality trends in South Africa: A mathematical modelling study. PLOS Med. 2017 Dec 12;14(12):e1002468.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1002468&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F14%2F2024.08.13.24311919.atom) 

12. 12.Eustaquio PC, Salmon-Trejo LAT, McGuire LC, Ellington SR. Epidemiologic and Clinical Features of Mpox in Adults Aged >50 Years — United States, May 2022–May 2023. Morb Mortal Wkly Rep. 2023 Aug 18;72(33):893–6.
    
    

13. 13.Crepaz N, Hess KL, Purcell DW, Hall HI. Estimating national rates of HIV infection among men who have sex with men, persons who inject drugs and heterosexuals in the United States. AIDS Lond Engl. 2019 Mar 15;33(4):701–8.
    
    

14. 14.Minhaj FS, Singh V, Cohen SE, Townsend MB, Scott H, Szumowski J, et al. Prevalence of Undiagnosed Monkeypox Virus Infections during Global Mpox Outbreak, United States, June– September 2022 - Volume 29, Number 11—November 2023 - Emerging Infectious Diseases journal - CDC. [cited 2024 Jul 25]; Available from: [https://wwwnc.cdc.gov/eid/article/29/11/23-0940\_article](https://www.nc.cdc.gov/eid/article/29/11/23-0940_article)
    
    

15. 15.aidsmap.com [Internet]. 2024 [cited 2024 Aug 13]. Nearly two million people in sub-Saharan Africa may have advanced HIV disease. Available from: [https://www.aidsmap.com/news/mar-2024/nearly-two-million-people-sub-saharan-africa-may-have-advanced-hiv-disease](https://www.aidsmap.com/news/mar-2024/nearly-two-million-people-sub-saharan-africa-may-have-advanced-hiv-disease)
    
    

16. 16.Endo A,  Jung S mok, Miura F. Mpox emergence in Japan: ongoing risk of establishment in Asia. The Lancet. 2023 Jun 10;401(10392):1923–4.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(23)00766-3&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=37269865&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F14%2F2024.08.13.24311919.atom)

 [1]: /embed/graphic-1.gif
 [2]: /embed/graphic-2.gif
 [3]: /embed/graphic-3.gif
 [4]: /embed/graphic-4.gif
 [5]: /embed/graphic-5.gif