Evaluating the trade-off between transmissibility and virulence of SARS-CoV-2 by mathematical modeling

1 Introduction

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), an RNA virus, can be transmitted by droplets that escape lungs through coughing or sneezing and infects humans (direct transmission), or is deposited in surfaces and infects humans when in contact with this contaminated surface (indirect transmission) [1] [2]. Coronavirus disease 2019 (CoViD-19) caused by SARS-CoV-2 infection was declared a pandemic by the World Health Organization on March 11, 2020. In general, the fatality rate in elder patients (60 years or more) is much higher than those with 60 years or less [3].

Since the worldwide outbreak of CoViD-19, like all RNA-based viruses, SARS-CoV-2 mutates faster but tends to change more slowly than HIV or influenza viruses [4]. “A virus with one or several new mutations (small changes during the replication) is referred to as a variant of the original virus” [4]. Currently, three variants challenge the vaccination’s control efforts – B.1.1.7 (United Kingdom), B.1.351 (South Africa), and P.1 (Brazil) (see [5] and references therein). Many countries closed their frontiers, mostly travelers from United Kingdom, South Africa, and Brazil, to control these variants. It is currently accepted that the more virulent variants of SARS-CoV-2 are also more transmissible based on the increased deaths caused by these strains [6].

Based on a mathematical model, we evaluate the trade-off between transmissibility and virulence. We hypothesize that the transmissibility is measured by the transmission rate (ultimately, by the basic reproduction number R₀), and the virulence is assessed by the ratio between asymptomatic and symptomatic individuals. Using the data collected from São Paulo State (Brazil), we estimate the basic reproduction number considering two virulence levels. We take as our framework the CoViD-19 epidemic in São Paulo State to evaluate the trade-off between transmissibility and virulence.

2 Material and methods

The SARS-CoV-2 transmission model is described by the system of ordinary differential equations presented in Appendix A, named the SQEAPMDR model. We drop out the pulses in equation (A.2) in Appendix A and transfer them to the boundary conditions. Hence, equations for susceptible individuals are for isolated individuals Q_j and infected individuals, with j = y, o, and for recovered individuals, where N_j = S_j + Q_j + E_j + A_j + P_j + M_j + D_j, and N = N_y + N_o + I obeys with the initial number of population being N (0) = N₀ = N_0y + N_0o, where N_0y and N_0o are the size of young and elder subpopulations at t = 0. Table 1 summarizes the model compartments (variables).

The initial and boundary conditions supplied to the system of equations (1), (2) and (3) are given by equations (A.3) and (A.5) in Appendix A.

Table 2 summarizes the model parameters. The description of the assigned values to the CoViD-19 natural history parameters can be found in [7]. The values of the parameters not shown in the table are estimated considering two SARS-CoV-2 virulence levels.

To estimate the transmission rates β_1j, β_2j, and β_3j, j = y, o, the proportion in isolation u, and the protection factor ε, we use the accumulated number of severe CoViD-19 cases Ω given by which are the exit from class P, and entering into class D.

To estimate the additional mortality rates α_j, j = y, o, we use the number of deaths due to severe CoViD-19 cases given by Π = Π_y + Π_o, where In the estimation of the additional mortality rates, we must bear in mind that the time at which new cases and deaths were registered does not have direct correspondence; instead, they are delayed by Δ days, that is, Π (t + Δ) = αD(t). Among exits from compartment D (γ_2j, µ, and α_j), we are counting only the deaths caused by severe CoViD-19.

In Appendix B, the steady-state of the system of equations in terms of fractions corresponding to equations (1), (2) and (3) was analyzed to obtain the basic reproduction number R₀, which is given by equation (B.10) in Appendix B. The fractions written as and result in where N_0y and N_0o are the initial numbers of young and elder subpopulations with N₀ = N_0y + N_0o, and Letting z_y = z_o = 0 (R_2y = R_2o = 0), we retrieve the basic reproduction number obtained in [8].

We choose two key parameters to evaluate the trade-off between transmissibility and virulence.

1. Transmissibility – The force of infection λ given by equation (A.1) in Appendix A is the per-capita incidence rate, and λS_y (λS_o) is the total number of the new cases in the young (elder) subpopulation per unit of time. Notwithstanding, the intensity of infection is proportional to the amount of virus released by infectious individuals and the capacity to infect susceptible individuals during contacts between them. The transmission rates carry on this fact besides other characteristics not considered here (social network, movement, demography, genetic, nutritional and health conditions, etc.). For this reason, we assume that the basic reproduction number R₀ calculated by the estimated transmission rates measures the transmissibility of SARS-CoV-2.

2. Virulence – Enhancing the virus’s capacity to infect target cells, an increased number of cells are infected. Consequently, the possibility of manifesting mild or severe disease increases, increasing the symptomatic CoViD-19 cases. For this reason, we translate the virulence with the ratio between asymptomatic and symptomatic individuals.

Notice that another effect of higher virulence is the increased amount of virus released by infectious individuals in the environment, which must increase the transmission. This fact is currently accepted for CoViD-19 epidemic [6]. We consider two broadly separated values for the proportion of asymptomatic individuals and estimate the transmission rates to calculate the basic reproduction number R₀.

3 Results

At the beginning of the epidemic, we have two data sets: Severe CoViD-19 cases (those in hospitals were tested and confirmed) and deaths. We study the transmissibility and virulence of SARS-CoV-2 considering the data collection from São Paulo State, Brazil [9]. Due to the lack of mass testing (PCR and serology) when SARS-CoV-2 transmission begins, the available severe CoViD-19 curve D is considered the epidemic curve.

To evaluate the trade-off between transmissibility and virulence, we assumed that the basic reproduction number R₀ measures the transmissibility, and the proportions of asymptomatic individuals l_y and l_o assess the virulence. Letting l_y = l_o, we consider two values: l₁ = 0.8 (lower virulence, with the ratio between asymptomatic and symptomatic individuals being 4 : 1) and l₂ = 0.2 (higher virulence, with the ratio 1 : 4).

To estimate the transmission rates, the proportion in isolation, and the protective factor (x stands for one of these parameters), we calculate where Ω (x, t_i) is the accumulated severe CoViD-19 cases calculated from equation (5) and Ω^ob (t_i) is the accumulated severe CoViD-19 registered cases at day t_i. We search for the value of x minimizing the Sum. To evaluate Ω, the numerical solutions of equations (1), (2) and (3) are obtained using the initial conditions given by equation (A.4) in Appendix A. To estimate the additional mortality rates, we substitute Ω (x, t_i) by Π (α, t_i) given by equation (6) and Ω^ob (t_i) by Π^ob (t_i).

The observed data from February 26 to November 31 [9] is partitioned into two sets – The first set is used to estimate the model parameters (input set), and the estimated model is then confronted with the second set to assess its prediction ability (test set) [10]. Hence, we estimate the model parameters using severe CoViD-19 (Ω^ob) and deaths (Π^ob) data from February 26 to May 13, 2020 (input data set). Then, we use the estimated parameters to predict the epidemic under interventions (quarantine and relaxation) from May 14 to November 31 (test data set) and compare the outcomes with the observed data.

We estimate the model parameters using the data from February 26 to May 13 following the estimation procedure presented in [7]. For higher virulence (l₂ = 0.2), we obtain (see Figure 1 below):

• Download figure
• Open in new tab

Figure 1:

The estimated accumulated curve Ω and the observed data Ω^ob (a), and the accumulated deaths Π and the observed Π^ob (b) for the higher (continuous line) and lower (dashed line) virulence. Estimations were done considering data from São Paulo State.

1. Transmission rates (from February 26 to April 3) – Estimated values are β_y = 0.69 and β_o = 0.79 (both in days⁻¹), resulting in the basic reproduction number R₀ = 7.02. When the CoViD-19 epidemic was initiated, we did not have any control mechanisms, which is characterized as the natural epidemic. The effect of isolation on March 24 appears around 9 days later.

2. Proportion in partial quarantine (from March 24 to April 13) – Estimated value is k = 0.48. The effect of adopting protective measures on April 4 appears around 9 days later.

3. Protective measures (from April 4 to May 13) – Estimated value is ε = 0.53.

4. Additional mortality rates (from March 16 to May 28) – Estimated values are α_y = 0.0018 and α_o = 0.0071 (both in days⁻¹) fixing Δ = 15 days. The first death due to CoViD-19 was on March 16, and we considered 15 more observed data.

The estimated parameters are summarized in Table 3. For lower virulence (l₂ = 0.8), we transport the estimated parameters from [7]. The period from February 26 to May 13 characterizes the beginning of the epidemic under partial quarantine (isolation) and protective measures adopted by individuals. Based on the two virulence levels’ estimated parameters considering the input data set, we study the epidemic’s prediction ability during the following test data set.

Figure 1 shows the estimated curve of accumulated CoViD-19 cases Ω and the observed data Ω^ob (a), and the accumulated deaths Π and the observed Π^ob (b) for the higher (continuous line) and lower (dashed line) virulence. Estimations were done considering São Paulo State’s data from February 26 to May 13, 2020, for Ω, and from March 16 to May 28 for Π.

The natural epidemic is characterized by the estimated transmission rates from February 26 to April 3. Since then, the observed data correspond to the epidemic with interventions (quarantine and relaxation). Figure 2 shows the extended accumulated Covid-19 cases Ω (a) and the corresponding natural epidemic curves D (b) for the higher (continuous line) and lower (dashed line) virulence.

• Download figure
• Open in new tab

Figure 2:

The extended accumulated Covid-19 cases Ω (a) and the corresponding natural epidemic curves D (b) for the higher (continuous line) and lower (dashed line) virulence.

Figure 2(a) shows that at the end of the epidemic’s first wave, Ω approaches the asymptote (plateaux) 9.46 × 10⁵ for lower virulence and 3.01 × 10⁶ for higher virulence. The curve of deaths Π (not shown) has the sigmoid shape as Ω and approaches the asymptote 5.73 × 10⁴ for lower virulence, and 1.95 × 10⁵ for higher virulence. Figure 2(b) shows that the natural epidemic’s peak is 3.86 × 10⁵, occurring on May 3, 2020, for lower virulence, and 1.17 × 10⁶ occurring on May 12 for higher virulence. In the natural epidemic (Figures 1 until April 3, and 2), the estimated curves of Ω and Π are quite the same in the initial phase for lower and higher virulence strains; however, the latter approaches to asymptotic values increased by around 330%. The peak of the higher virulent SARS-CoV-2 epidemic curve is increased by around 300% from the lower one; hence we take the peaks of curves and not the asymptotic values to compare epidemic outcomes.

Let us evaluate the estimated parameters’ prediction ability considering two virulence levels on the epidemic under interventions in São Paulo State. Besides the input data set, we include the test data set and show the observed data Ω^ob and Π^ob from February 26 to November 31. The model is solved numerically during the input and test data sets using the estimated parameters given in Table 3. The resulting epidemiological scenarios correspond to the quarantine epidemic and are shown in Figures 3 to 5.

• Download figure
• Open in new tab

Figure 3:

The extended curves of the accumulated CoViD-19 cases Ω (a) and deaths Ω (b) for higher (continuous line) and lower (dashed line) virulence. The curves correspond to the epidemic under isolation.

Figure 3 shows the curves of the accumulated CoViD-19 cases Ω (a) and deaths Π (b) for higher (continuous line) and lower (dashed line) virulence. Notice that Figure 3 is the extension of Figure 1. The lower and higher virulence estimations’ prediction ability is analyzed below by comparing the epidemic outcomes with the observed data Ω^ob and Π^ob.

From Figure 3, for lower virulence, the accumulated CoViD-19 cases at the end of the epidemic under isolation’s first wave Ω and deaths Π are 3.85 × 10⁵ (41%) and 2.37 × 10⁴ (41.3%), and for higher virulence, they are 1.26 × 10⁶ (41.9%) and 8.29 × 10⁴ (42.5%). The percentage between parentheses is the ratio to the natural epidemic. The partial quarantine decreased by around 40% of the severe CoViD-19 cases and deaths from the natural epidemic in lower and higher virulence. For the higher virulence, the asymptotic values of Ω and Π increased by around 340% from the lower virulence, close to the natural epidemic’s increase.

Figure 4 shows the curves of asymptomatic (A) and pre-symptomatic (P) individuals (a), and the estimated curves of mild (M) and severe (D) CoViD-19 cases (b) for higher (continuous line) and lower (dashed line) virulence. Figure 4 also shows the curve of carriers capable of transmitting the virus (A + P) and the curve of apparent CoViD-19 (M + D).

• Download figure
• Open in new tab

Figure 4:

The estimated curves of asymptomatic (A) and pre-symptomatic (P) individuals (a), and the estimated curves of mild (M) and severe (D) CoViD-19 cases (b) for higher (continuous line) and lower (dashed line) virulence.

Figure 4 illustrates the impact of control mechanisms on flattening the natural epidemic curve. For lower virulence, the peaks of asymptomatic (A), pre-symptomatic (P), mild (M), and severe (D) CoViD-19 are 2.23 × 10⁶, 1.95 × 10⁵, 5.86 × 10⁵, and 0.66 × 10⁵, and for higher virulence, the peaks are 4.85 × 10⁵, 6.64 × 10⁵, 1.91 × 10⁶, and 1.89 × 10⁵. (For the lower virulence, the peak of A is increased by 460%, and the peaks of P, M, and D are reduced by 340%, 326%, and 286% from the higher virulence.) The asymptomatic and pre-symptomatic individuals are carriers, which difficult the control of infection. Comparing the lower and higher virulence, the ratio between the peaks of carriers (A + P) is 212%, and the ratio between peaks of apparent CoViD-19 (M + D) is 31%.

Figure 5 shows the curves of transmitting SARS-CoV-2 individuals (A + P + z_yM_y + z_oM_o) (a), and the curves of harboring SARS-CoV-2 individuals (E, A, P, M, D) (b) for higher (continuous line) and lower (dashed line) virulence.

• Download figure
• Open in new tab

Figure 5:

The curves of transmitting SARS-CoV-2 individuals (a), and the curves of harboring SARS-CoV-2 individuals (b) for higher (continuous line) and lower (dashed line) virulence.

From Figure 5(a), for lower virulence, the peaks of transmitting SARS-CoV-2 individuals (A +P +z_yM_y +z_oM_o) for lower and higher virulence are 2.68 × 10⁶ and 1.99 × 10⁶, occurring on June 19 and July 21. The ratio between peaks of all infectious individuals with lower and higher virulence is 135%. From Figure 5(b), for lower virulence, the peaks of exposed (E), releasing virus (A + P + M + D), and harboring virus (E + A + P + M + D) individuals are 1.43 × 10⁶, 3.06 × 10⁶, and 4.41 × 10⁶, and for higher virulence, the peaks are 1.21 × 10⁶ (85%), 3.21 × 10⁶ (105%), and 4.36 × 10⁶ (99%). The percentage between parentheses is the ratio to the lower virulence. The individuals releasing the virus are slightly higher in a higher virulence, but all individuals harboring the virus are practically equal in lower and higher virulence transmissions. Notice that exposed individuals (the new infection cases) are higher in lower virulence due to a higher R₀. The peaks occur close to those observed in Figure 5(a).

Figures 3-5 showed that the enhanced virulence increased the deaths due to CoViD-19 by around 300%; however, the transmissibility measured by R₀ was reduced by around 76%. The higher virulence increased by around 300% the number of mild and severe CoViD-19 cases, potentially more transmissible but in isolation; however, the number of asymptomatic, presymptomatic, and a fraction of mild CoViD-19 individuals transmitting the virus is reduced by 74%. The SARS-CoV-2 transmission decreases in the populational point of view, and the number of deaths rises as the virulence increases. Initially, the model variables (compartments) have quite the same values, but the peaks of the more virulent SARS-CoV-2 occur around one month later. In particular, the number of deaths in the more virulent epidemic is three times higher than the less one, although both strains presented the same fatalities during the first three months.

Let us assess the prediction ability by comparing the observed data with the lower and higher virulence epidemiological scenarios. Remember that the lower and higher virulent SARS-CoV-2, considering the input data set comprising the quarantine period, provided quite similar fitted Ω and Π curves (see Figure 1). However, from Figure 3, the lower virulent epidemic fits very well the subset of test data from May 13 to June 30 (Ω) or July 20 (Π), showing that quarantine lasted until that date. (The first observed point detaching from the lower virulence curve Π occurs on July 21, while the severe CoViD-19 curve Ω separates on July 1, 20 days earlier. In the estimation of mortality rates α_y and α_o, we used Δ = 15 days in Π (t + Δ) = αD(t), not Ω.) São Paulo State initiated the relaxation by mid-June 2020, with the effect on the epidemic curve D appearing 9 days later. Hence, the subset of test data since July 1 must be described by a model incorporating intermittent pulses of release [8], which is not dealt with here. Therefore, the less virulent SARS-CoV-2 transmission explains better the CoViD-19 epidemic until the beginning of the relaxation.

The relation Π (t + Δ) = αD(t) does not sustain due to the mass test including mild and eventually asymptomatic and pre-symptomatic individuals in the observed data. This fact is apparent in Figure 3, where the registered SARS-CoV-2 infections increase more quickly than the observed deaths. The supplementary deaths since July 21 can be explained by the beginning of the relaxation 20 days earlier. Moreover, the registered number of deaths in São Paulo State is 58873 on February 25, 2021, which is higher than 57300 predicted by the low virulent natural epidemic (see underneath Figure 2). However, São Paulo State is still in partial quarantine. Hence, based on the model, the supplementary deaths can not be explained only by the relaxation but by the circulation of more virulent variant(s) of the original SARS-CoV-2.

4 Discussion and conclusion

The persistent transmission of SARS-CoV-2 should result in a high number of mutations and variants. In this pool of variants, some acquired enhanced ability to infect cells, increasing the risk of death. Moreover, more viruses will be released in the environment when more cells are infected, increasing the transmission. However, the severe symptomatic individuals are isolated in hospitals, and mild cases are detected by test and advised to isolate. Contrarily, the asymptomatic and pre-symptomatic individuals (carriers) are not seen and circulate freely. Additionally, the carriers are in close contact with other susceptible individuals, contrarily to individuals manifesting any suspicious symptoms (not necessarily CoViD-19). This fact was explained by the model, showing that higher virulence causes increased deaths but is less transmissible – By the fact that R₀ is lower, the epidemic curve increases more slowly, but the curve of deaths increases much more quickly than the less virulent one.

In the first months of the CoViD-19 outbreak, the absence of a mass test against SARS-CoV-2 infection and clinical follow-up of infected individuals did not permit the discrimination between asymptomatic and pre-symptomatic individuals. Considering this initial period, we explained the CoViD-19 data from São Paulo State using two broadly separated values for the ratio between asymptomatic and symptomatic individuals (4 : 1 and 1 : 4). For each proportion l, we estimated the transmission rates and calculated the basic reproduction number R₀.

Figures 3-5 showed two broadly separated virulence to explain CoViD-19 data from São Paulo State. The partial quarantine and protective measures flattened the epidemic curve in São Paulo State, but the transmission was maintained at a relatively higher level (the effective reproduction number reached one around 80 days later [7]). This epidemiological scenario can be explained by spreading a less virulent original SARS-CoV-2 or a pool of less virulent variants dominating during the quarantine period. However, since July 1, the observed number of deaths was systematically above the estimated curve of Π considering quarantine only. The effects of relaxation initiated in São Paulo State appeared on July 1, but it was not fully implemented. We do not have a realistic assessment of the relaxation neither the transmission of SARS-CoV-2 among the isolated individuals (in [7], a decreased transmission in the isolated population is allowed). At the end of February 2021, São Paulo State still maintained part of the population in isolation, but 58873 deaths on February 25, 2021, surpassed 57300 predicted by the natural epidemic under quarantine. This higher number of deaths under isolation may be a consequence of the appearance of more virulent variants (but not so high as l = 0.2 shown in Figures 1-5) during the one-year SARS-CoV-2 transmission. Hence, based on the model’s analysis, we can assume that virulent(s) SARS-CoV-2 can be found in the current pool of mutants, resulting in a continuously increasing number of severe CoViD-2 and deaths. However, more detection of higher virulence severe cases and deaths is synonymous with more transmissibility is a wrong statement according to the model.

The epidemiological scenario in Amazonas State (Brazil) can shed more light on the virulent variants’ role in association with relaxation. In Amazonas State, the occurrence of deaths in January-February, 2021, is higher than that occurred during all 2020 [11]. This increased number can not be explained only by the partial relaxation once a significant number of deaths caused by the P.1 variant was found. Hence, Amazonas State is another example of a pool containing high virulent SARS-CoV-2 transmission.

Finally, during the epidemic, the fast mutations in RNA-virus can result in a pool of variants. This pool can behave on average with lower or higher virulence than the original SARS-CoV-2. In the short-run epidemic, if a pool of variants acts on average as a high (low) virulent strain, the number of severe cases and deaths increases (decreases), but the transmission is decreased (increased). However, this trade-off between transmissibility and virulence shows that, in the long-run epidemic, lower virulence with higher transmissible variant(s) will prevail. Therefore, in the short-run epidemic control, strict isolation of symptomatic and perhaps suspicious individuals must be adopted to avoid high virulent SARS-CoV-2 transmission causing more death risk.

Funding

This research received no specific grant from any funding agency, commercial or not-for-profit sectors.

Conflicts of interest/Competing interests

Not applicable.

Availability of data and material

The data that support the findings of this study are openly available in SP contra o novo coronavírus (Boletim completo) at https://www.seade.gov.br/coronavirus/.

Author contributions

Hyun Mo Yang: Conceptualization, Methodology, Formal analysis, Writing - Original draft preparation, Validation, Supervision. Luis Pedro Lombardi Junior: Software, Data Curation, Visualization, Validation. Ariana Campos Yang: Conceptualization, Validation, Investigation.

A The SQEAPMDR model

One of the main aspects of CoViD-19 is increased fatality in the elder subpopulation. For this reason, a population is divided into two groups, composed of young (60 years old or less, denoted by subscript y) and elder (60 years old or more, denoted by subscript o) subpopulations. This community’s vital dynamic is described by the per-capita rates of birth (ϕ) and mortality (µ), and φ is the aging rate, that is, the flow from young subpopulation y to elder subpopulation o. Another aspect is the presence of the pre-symptomatic individuals, that is, individuals without symptoms transmitting SARS-CoV-2 before the onset of the disease [12].

Hence, for each subpopulation j (j = y, o), individuals are divided into seven compartments: susceptible S_j, isolated (quarantine) Q_j, exposed and incubating E_j, asymptomatic A_j, pre-symptomatic (or pre-diseased) individuals P_j, symptomatic individuals with mild CoViD-19 M_j, and severe CoViD-19 D_j. However, all young and elder individuals in classes A_j, M_j, and D_j enter into the same recovered class R (this is the 8^th class, but common to both subpopulations). Hence, the SQEAPMDR model has 15 compartments.

The natural history of CoViD-19 is the same for young (j = y) and elder (j = o) subpopulations. We assume that individuals in the asymptomatic (A_j), pre-diseased (P_j), and a fraction z_j of mild CoViD-19 (M_j) classes are transmitting the virus. Other infected classes ((1 − z_j) M_j and D_j) are under voluntary or forced isolation. Susceptible individuals are infected at a rate λ_jS_j (known as the mass action law [13]), where λ_j is the per-capita incidence rate (or force of infection) defined by λ_j = λ (δ_jy + ψδ_jo), with λ being where δ_ij is the Kronecker delta, with δ_ij = 1 if i = j, and 0, if i = j; and β_1j, β_2j and β_3j are the transmission rates, that is, the rates at which a virus encounters a susceptible people and infects him/her. In [8], a particular model was analyzed letting z_y = z_o = 0 and χ_y = χ_o = 1.

Susceptible individuals are infected at a rate λ_j and enter into class E_j. After an average period 1/σ_j in class E_j, where σ_j is the incubation rate, exposed individuals enter into the asymptomatic class A_j (with probability l_j) or pre-diseased class P_j (with probability 1 − l_j). After an average period 1/γ_j in class A_j, where γ_j is the recovery rate of asymptomatic individuals, asymptomatic individuals acquire immunity (recovered) and enter into recovered class R. Possibly asymptomatic individuals can manifest symptoms at the end of this period, and a fraction 1 − χ_j enters into mild CoViD-19 class M_j. For symptomatic individuals, after an average period 1/γ_1j in class P_j, where γ_1j is the infection rate of pre-diseased individuals, pre-diseased individuals enter into severe CoViD-19 class D_j (with probability 1 − k_j) or mild CoViD-19 class M_j (with probability k_j). Individuals in class D_j acquire immunity after a period 1/γ_2j, where γ_2j is the recovery rate of severe CoViD-19, and enter into recovered class R or die under the disease-induced (additional) mortality rate α_j. Individuals in mild CoViD-19 class M_j acquire immunity after a period 1/γ_3j, where γ_3j is the recovery rate of mild CoViD-19, and enter into recovered class R.

The dynamic equations for E, A, P, M, D, and R are obtained through the balance between inflow and outflow in each compartment (see equations (2) and (3) in the main text). For the equations corresponding to S and Q, we consider a unique pulse in isolation at time t = τ, described by u_jS_jδ (t − τ), with j = y, o. The fraction of persons in isolation is u_j and δ (x) is the Dirac delta function, that is, δ (x) =∞, if x = 0, otherwise, δ (x) = 0, with . The isolation is translated as These equations are simulated permitting intermittent interventions to the boundary conditions.

The initial conditions (simulation time t = 0) supplied to the system of equations (1), (2), and (3) are, for j = y, o, where is a non-negative number. For instance, means that there is not any exposed person (young and elder) at the beginning of the epidemic. In this paper, we use plus Q_y (0) = Q_o (0) = R(0) = 0, where t = 0 corresponds to the calendar time when the first case was confirmed (February 26 for São Paulo State). For São Paulo State, N_0y = 37.8 million and N_0o = 6.8 million. (See [7] for details in the initial conditions’ setup.)

The isolation implemented at τ = 27 (corresponding to calendar time March 24, 2020) is described by the boundary conditions with τ⁻ = lim_t→τ t (for t < τ), and τ ⁺ = lim_τ←t t (for t > τ).

B The steady-state analysis of the SQEAPMDR model

The system of equations (1), (2) and (3) does not reach steady-state (time varying population). However, this system in term of fractions attains steady-state. Defining the fraction x_j = X_j/N, for j = y, o, with X_j = {S_j, Q_j, E_j, A_j, P_j, M_j, D_j, R}, we have using equation (4), and the system of equations (1), (2) and (3) in terms of fractions become, for susceptible individuals, for isolated and infected individuals, and for recovered individuals, where λ is the force of infection given by equation (A.1) re-written as and This new system of equations has a steady-state, that is, the number of individuals in all classes varies with time. However, their fractions attain a steady-state (the sum of derivatives of all classes is zero).

The trivial (disease-free) equilibrium point P ⁰ of the new system of equations (B.6), (B.7) and (B.8) is given by for j = y and o, where with .

Let us assess the stability of P ⁰ by applying the next generation matrix theory considering the vector of variables x = (e_y, a_y, p_y, m_y, e_o, a_o, p_o, m_o) [14]. We apply method proposed in [15] and proved in [16]. To obtain the basic reproduction number, diagonal matrix V is considered. Hence, the vectors f and v are and where the superscript T stands for the transposition of a matrix, from which we obtain the matrices F and V (see [14]) evaluated at the trivial equilibrium P ⁰, which were omitted. The next generation matrix FV ⁻¹ is and the characteristic equation corresponding to FV ⁻¹ is with the basic reproduction number R₀ being given by where the initial fractions and are given by equation (B.9), and the partial basic reproduction numbers R_1y, R_2y, R_1o, and R_2o are given by equation (8) in the main text. The spectral radius ρ (FV ⁻¹) is the biggest solution of a third-degree polynomial, not easy to evaluate. The procedure proposed in [15] allows us to obtain the threshold R₀ as the sum of coefficients of the characteristic equation. Hence, the trivial equilibrium point P ⁰ is locally asymptotically stable if R₀ < 1.