Early estimates of the incidence trend and the reproductive number of the monkeypox epidemic in Brazil

INTRODUCTION

In July of 2022, the World Health Organization declared that the monkeypox outbreak was a public health emergency of international concern (1). Until August 13, 2022, there were 35.414 confirmed cases and seven deaths worldwide. In Brazil, the first case was confirmed on June 8, 2022, and 2.849 have been confirmed cases until August 13, 2022 (2). Incidence trends and transmissibility measures are essential to recognize the potential extension of the epidemic and to guide public health interventions (3). Consequently, we estimated the relative increase and the Rt (effective reproductive number) of the monkeypox incidence in Brazilian states with the highest number of cases.

METHODS

This was an ecological study based on surveillance data of the epidemiological reports posted by the Brazilian National Health Ministry (4). We created a database with the new and accumulated cases (5). Our analysis focused on the four states with more than 100 cases up to August 13, 2022. These states were São Paulo (SP), Rio de Janeiro (RJ), Minas Gerais (MG), and Goiás, which geographically includes the Federal District of Brasilia (GO-FD).

We only analyzed confirmed cases defined as suspected cases with positive PCR results for the Monkeypox virus (6). We estimated the weekly growth of the incidence of Monkeypox using two regression approaches starting one day after the first case was reported. First, we used the count models of Poisson and negative binomial regression and reported the last one when the alpha term differed significantly from zero. We also applied the Prais-Wisten regression to the natural logarithm of the weekly case count. For each state, we reported the relative incidence increase calculated as 1-e^β, where β is the regression coefficient.

To estimate Rt, we used the “estimate_R” package contained in the “EpiEstim” library of the R programming language (7). We used the “parametric_si” configuration of the package. The script is available as supplementary material. One of the parameters needed to estimate the Rt is the serial interval (SI), i.e., the time between the illness onset in the primary case and that in a secondary. Since the microdata of Monkeypox cases in Brazil are unavailable, we used a predefined SI recently calculated in the United Kingdom (8). That SI was obtained from 17 pairs of cases (primary - secondary), exhibiting a mean of 9.8 days and a standard deviation of 9.1 days (8). For each state, we calculated the mean Rt after the 20th day since the outbreak, when we observed a more stable daily Rt.

This study used public data, and the analysis was performed using R (version 4.2.1) and Stata (version 17.0, Stata Corp LP, College Station, TX, USA).

RESULTS

The state with the most cases was SP, followed by RJ and GO-FD (figure 1). However, the highest relative increase of incidence was observed in GO-FD, being 98.1% per week according to the Poisson model and 95.1% using the Prais-Winsten regression. In all states, the increasing trend was statistically significant with both regression approaches (Table).

• Download figure
• Open in new tab

Figure 1:

Brazil maps with Monkeypox incidence by region: North (A), Northeast (B), Center West (C), Southeast (D) and South (E)

The daily Rt oscillated between 1.9 and 3.1 in GO-FD, 1.5 and 2.2 in SP, 1.6 and 2.1 in RJ, and 1.6 and 2.0 in MG (figure 2). The means of Rt in those states were respectively: 2.28, 1.94, 1.77, and 1.71 (Table).

• Download figure
• Open in new tab

Figure 2:

Effective reproductive number (Rt) of monkeypox epidemic in four Brazilian states. GO-FD (A), SP (B), RJ (C) and MG (D).

DISCUSSION

The confirmed Monkeypox cases in Brazil up to August 13, 2022, was more than double of the 1.369 cases that were notified up to July 31, 2022. We observed that the Rt of the four analyzed states may be comparable with some calculated for covid-19 (9). Monkeypox transmission appears to mainly occur by direct contact in the current outbreak (10) has a high viral load on skin lesions (11) and there is evidence of a lesion staying positive on the PCR test for up to 21 days (12). We hypothesize that a prolonged infectious period could explain those Rt values.

Our study had some methodological limitations. The sample size was limited and Rt can be unstable during the early stage of epidemic. Therefore, we restricted our analysis to four states and Rt was estimated after the 20th day of the first case. It is also important to note that we estimated the Rt with a fixed SI obtained from another population (8). Even though the virus is the same involved in the current outbreak, we recommend replicating the Rt estimations using locally obtained SI (when available). Nonetheless, by using this preset SI we could compare the profile of the states with a higher number of Monkeypox cases. In that sense, our Rt estimates correlated well with the observed upward trends across the states, which give a reference of the potential transmissibility of Monkeypox.

This study illustrated the critical and significant growth of Monkeypox cases in Brazil, highlighting the need for effective public health measures such as testing, contact tracing, and vaccination of the risk groups. It also standed out the importance of surveillance to ensure that all cases are timely studied and their potential determinants well documented.

Conflicts of Interest

The authors declare no conflict of interest.

Funding Statement

This work had no specific funding. However, FADQ is beneficiary of a fellowship for research productivity from the National Council for Scientific and Technological Development - CNPq, process/contract identification: 312656/2019-0.