Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

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

Hyun Mo Yang, Luis Pedro Lombardi Junior, Ariana Campos Yang
doi: https://doi.org/10.1101/2021.02.27.21252592
Hyun Mo Yang
a,bDepartment of Applied Mathematics, State University of Campinas Praça Sérgio Buarque de Holanda, 651; CEP: 13083-859, Campinas, SP, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: hyunyang@ime.unicamp.br
Luis Pedro Lombardi Junior
a,bDepartment of Applied Mathematics, State University of Campinas Praça Sérgio Buarque de Holanda, 651; CEP: 13083-859, Campinas, SP, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ariana Campos Yang
cDivision of Allergy and Immunology, General Hospital of the Medicine School of University of São Paulo Av. Dr. Eneas Carvalho de Aguiar, 255; CEP: 05403-000, São Paulo, SP, Brazil
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

Background At the beginning of 2020, SARS-CoV-2 spread to all continents, and since then, mutations have appeared in different regions of the world. The appearance of more virulent mutations leads to asseverate that they are also more transmissible. We analyzed the lower and higher virulent SARS-CoV-2 epidemics to establish a relationship between transmissibility and virulence based on a mathematical model.

Methods A compartmental mathematical model based on the CoViD-19 natural history encompassing the age-dependent fatality was applied to evaluate the SARS-CoV-2 transmissibility and virulence. The transmissibility was measured by the basic reproduction number R0 and the virulence by the proportion of asymptomatic individuals. The model parameters were fitted considering the observed data from São Paulo State.

Results The numbers of severe CoViD-19 and deaths are three times higher, but R0 is 25% lower in more virulent SARS-CoV-2 transmission than in a less virulent one. However, the number of more virulent SARS-CoV-2 transmitting individuals is 25% lower, mainly due to symptomatic individuals’ isolation, explaining the increased transmission in lower virulence.

Conclusions The quarantine study in São Paulo State showed that the more virulent SARS-CoV-2 resulted in a higher number of fatalities but less transmissible than the less virulent one. One possible explanation for the number of deaths surpassing that predicted by the low virulent SARS-CoV-2 infection could be the transmission of more virulent variant(s).

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 R0), 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 Embedded Image for isolated individuals Qj and infected individuals, with j = y, o, Embedded Image and for recovered individuals, Embedded Image where Nj = Sj + Qj + Ej + Aj + Pj + Mj + Dj, and N = Ny + No + I obeys Embedded Image with the initial number of population being N (0) = N0 = N0y + N0o, where N0y and N0o are the size of young and elder subpopulations at t = 0. Table 1 summarizes the model compartments (variables).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1:

Summary of the model variables, for j = y (young) and j = o (elder).

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.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2:

Summary of the model parameters (j = y, o) and values (rates in days−1 and proportions are dimensionless). (*) Values fitted considering two levels of the SARS-CoV-2 virulence.

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 Embedded Image 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 Embedded Image 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 R0, which is given by equation (B.10) in Appendix B. The fractions written as Embedded Image and Embedded Image result in Embedded Image where N0y and N0o are the initial numbers of young and elder subpopulations with N0 = N0y + N0o, and Embedded Image Letting zy = zo = 0 (R2y = R2o = 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 λSy (λSo) 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 R0 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 R0.

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 R0 measures the transmissibility, and the proportions of asymptomatic individuals ly and lo assess the virulence. Letting ly = lo, we consider two values: l1 = 0.8 (lower virulence, with the ratio between asymptomatic and symptomatic individuals being 4 : 1) and l2 = 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 Embedded Image where Ω (x, ti) is the accumulated severe CoViD-19 cases calculated from equation (5) and Ωob (ti) is the accumulated severe CoViD-19 registered cases at day ti. 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, ti) by Π (α, ti) given by equation (6) and Ωob (ti) by Πob (ti).

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 (l2 = 0.2), we obtain (see Figure 1 below):

Figure 1:
  • 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−1), resulting in the basic reproduction number R0 = 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−1) 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 (l2 = 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.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 3:

The estimated parameters (j = y, o) considering two SARS-CoV-2 virulence levels (rates in days−1 and proportions are dimensionless).

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.

Figure 2:
  • 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 × 105 for lower virulence and 3.01 × 106 for higher virulence. The curve of deaths Π (not shown) has the sigmoid shape as Ω and approaches the asymptote 5.73 × 104 for lower virulence, and 1.95 × 105 for higher virulence. Figure 2(b) shows that the natural epidemic’s peak is 3.86 × 105, occurring on May 3, 2020, for lower virulence, and 1.17 × 106 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.

Figure 3:
  • 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 × 105 (41%) and 2.37 × 104 (41.3%), and for higher virulence, they are 1.26 × 106 (41.9%) and 8.29 × 104 (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).

Figure 4:
  • 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 × 106, 1.95 × 105, 5.86 × 105, and 0.66 × 105, and for higher virulence, the peaks are 4.85 × 105, 6.64 × 105, 1.91 × 106, and 1.89 × 105. (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 + zyMy + zoMo) (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.

Figure 5:
  • 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 +zyMy +zoMo) for lower and higher virulence are 2.68 × 106 and 1.99 × 106, 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 × 106, 3.06 × 106, and 4.41 × 106, and for higher virulence, the peaks are 1.21 × 106 (85%), 3.21 × 106 (105%), and 4.36 × 106 (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 R0. 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 R0 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 R0 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 R0.

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.

Data Availability

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

https://www.seade.gov.br/coronavirus/

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 Sj, isolated (quarantine) Qj, exposed and incubating Ej, asymptomatic Aj, pre-symptomatic (or pre-diseased) individuals Pj, symptomatic individuals with mild CoViD-19 Mj, and severe CoViD-19 Dj. However, all young and elder individuals in classes Aj, Mj, and Dj enter into the same recovered class R (this is the 8th 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 (Aj), pre-diseased (Pj), and a fraction zj of mild CoViD-19 (Mj) classes are transmitting the virus. Other infected classes ((1 − zj) Mj and Dj) are under voluntary or forced isolation. Susceptible individuals are infected at a rate λjSj (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 Embedded Image 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 zy = zo = 0 and χy = χo = 1.

Susceptible individuals are infected at a rate λj and enter into class Ej. After an average period 1/σj in class Ej, where σj is the incubation rate, exposed individuals enter into the asymptomatic class Aj (with probability lj) or pre-diseased class Pj (with probability 1 − lj). After an average period 1/γj in class Aj, 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 Mj. For symptomatic individuals, after an average period 1/γ1j in class Pj, where γ1j is the infection rate of pre-diseased individuals, pre-diseased individuals enter into severe CoViD-19 class Dj (with probability 1 − kj) or mild CoViD-19 class Mj (with probability kj). Individuals in class Dj 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 Mj 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 ujSjδ (t − τ), with j = y, o. The fraction of persons in isolation is uj and δ (x) is the Dirac delta function, that is, δ (x) =∞, if x = 0, otherwise, δ (x) = 0, with Embedded Image. The isolation is translated as Embedded Image 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, Embedded Image where Embedded Image is a non-negative number. For instance, Embedded Image means that there is not any exposed person (young and elder) at the beginning of the epidemic. In this paper, we use Embedded Image plus Qy (0) = Qo (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, N0y = 37.8 million and N0o = 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 Embedded Image with τ− = limt→τ 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 xj = Xj/N, for j = y, o, with Xj = {Sj, Qj, Ej, Aj, Pj, Mj, Dj, R}, we have Embedded Image using equation (4), and the system of equations (1), (2) and (3) in terms of fractions become, for susceptible individuals, Embedded Image for isolated and infected individuals, Embedded Image and for recovered individuals, Embedded Image where λ is the force of infection given by equation (A.1) re-written as Embedded Image and Embedded Image 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 0 of the new system of equations (B.6), (B.7) and (B.8) is given by Embedded Image for j = y and o, where Embedded Image with Embedded Image.

Let us assess the stability of P 0 by applying the next generation matrix theory considering the vector of variables x = (ey, ay, py, my, eo, ao, po, mo) [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 Embedded Image and Embedded Image 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 0, which were omitted. The next generation matrix FV −1 is Embedded Image and the characteristic equation corresponding to FV −1 is Embedded Image with the basic reproduction number R0 being given by Embedded Image where the initial fractions Embedded Image and Embedded Image are given by equation (B.9), and the partial basic reproduction numbers R1y, R2y, R1o, and R2o are given by equation (8) in the main text. The spectral radius ρ (FV −1) is the biggest solution of a third-degree polynomial, not easy to evaluate. The procedure proposed in [15] allows us to obtain the threshold R0 as the sum of coefficients of the characteristic equation. Hence, the trivial equilibrium point P 0 is locally asymptotically stable if R0 < 1.

Footnotes

  • a emails: hyunyang{at}ime.unicamp.br

  • b luispedro jr{at}hotmail.com

  • c arianacy{at}gmail.com

References

  1. [1].↵
    Yuan L, et al. Aerodynamic analysis of SARS-CoV-2 in two Wuhan hospitals. Published April 27, 2020. Available at: https://doi.org/10.1038/s41586-020-2271-3 (accessed May 10, 2020).
  2. [2].↵
    Zhu N, et al. A novel Coronavirus from patients with pneumonia in China, 2019. N. Eng. J. Med. 2020; 381: 1–7. doi: 10.1056/NEJMoa2001017.
    OpenUrlCrossRefPubMed
  3. [3].↵
    WHO, Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19), 16-24 February 2020 (2020).
  4. [4].↵
    WHO, Coronavirus Disease (COVID-19): Virus Evolution. 30 December 2020. Available at: https://www.who.int/news-room/q-a-detail/sars-cov-2-evolution (accessed February 23, 2021).
  5. [5].↵
    CDC, SARS-CoV-2 Variants. Updated Jan. 31, 2021. Available at: https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/variant-surveillance/variant-info.html (accessed February 23, 2021).
  6. [6].↵
    Korber B, et al. Tracking changes in SARS-CoV-2 spike: Evidence that D614G increases infectivity of the COVID-19 virus. Cell 2021; 182: 812–817. doi: https://doi.org/10.1016/j.cell.2020.06.043.
    OpenUrl
  7. [7].↵
    Yang HM, Lombardi Junior LP, Castro FFM, Campos AC. Mathematical modeling of the transmission of SARS-CoV-2 – Evaluating the impact of isolation in São Paulo State (Brazil) and lockdown in Spain associated with protective measures on the epidemic of covid-19. PlosOne 2020; in revision.
  8. [8].↵
    Yang HM, Lombardi Junior LP, Castro FFM, Yang AC. Mathematical model describing Covid-19 in São Paulo State, Brazil – Evaluating isolation as control mechanism and forecasting epidemiological scenarios of release. Epidemiology and Infection 2020 148: e155. doi: 10.1017/S0950268820001600.
    OpenUrlCrossRef
  9. [9].↵
    SEADE, SP contra o novo coronavírus – Boletim completo. Available from: https://www.seade.gov.br/coronavirus/?utmsource=portal&utmmedium=banner&utmcampaign=boletim-completo (accessed August 14, 2020).
  10. [10].↵
    Alvarenga MY, Sameshima K, Baccal á LA, Yang HM. Non-linear analysis of the rhythmic activity in rodent brain. Mathematical Biosciences 1999; 157 (1-2): 287–302.
    OpenUrlPubMed
  11. [11].↵
    Informações COVID-19. http://www.transparencia.am.gov.br/covid-19/monitoramento-covid-19/ (accessed February 27, 2021).
  12. [12].↵
    Arons MM, et al. Presymptomatic SARS-CoV-2 infections and transmission in a skilled nursing facility. The New Engl. Jour. Medicine 2020; April 24, 2020. doi: 10.1056/NEJ-Moa2008457.
    OpenUrlCrossRef
  13. [13].↵
    Anderson RM, May, RM. Infectious Diseases of Human. Dynamics and Control. Oxford, New York, Tokyo: Oxford University Press; 1991:757 p.
  14. [14].↵
    Diekmann O, Heesterbeek JAP, Roberts MG. The construction of next-generation matrices for compartmental epidemic models. J. R. Soc. Interface 2010; 7:873–885.
    OpenUrlCrossRefPubMedWeb of Science
  15. [15].↵
    Yang HM. The basic reproduction number obtained from Jacobian and next generation matrices – A case study of dengue transmission modelling. BioSystems 2014; 126:52–75.
    OpenUrl
  16. [16].↵
    Yang HM, Greenhalgh D. Proof of conjecture in: The basic reproduction number ob-tained from Jacobian and next generation matrices – A case study of dengue transmission modelling. Appl. Math. Comput. 2015; 265:103–107.
    OpenUrl
Back to top
PreviousNext
Posted March 01, 2021.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Evaluating the trade-off between transmissibility and virulence of SARS-CoV-2 by mathematical modeling
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Evaluating the trade-off between transmissibility and virulence of SARS-CoV-2 by mathematical modeling
Hyun Mo Yang, Luis Pedro Lombardi Junior, Ariana Campos Yang
medRxiv 2021.02.27.21252592; doi: https://doi.org/10.1101/2021.02.27.21252592
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Evaluating the trade-off between transmissibility and virulence of SARS-CoV-2 by mathematical modeling
Hyun Mo Yang, Luis Pedro Lombardi Junior, Ariana Campos Yang
medRxiv 2021.02.27.21252592; doi: https://doi.org/10.1101/2021.02.27.21252592

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Epidemiology
Subject Areas
All Articles
  • Addiction Medicine (162)
  • Allergy and Immunology (416)
  • Anesthesia (91)
  • Cardiovascular Medicine (863)
  • Dentistry and Oral Medicine (159)
  • Dermatology (98)
  • Emergency Medicine (251)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (394)
  • Epidemiology (8573)
  • Forensic Medicine (4)
  • Gastroenterology (388)
  • Genetic and Genomic Medicine (1759)
  • Geriatric Medicine (167)
  • Health Economics (373)
  • Health Informatics (1250)
  • Health Policy (622)
  • Health Systems and Quality Improvement (469)
  • Hematology (196)
  • HIV/AIDS (378)
  • Infectious Diseases (except HIV/AIDS) (10318)
  • Intensive Care and Critical Care Medicine (553)
  • Medical Education (192)
  • Medical Ethics (51)
  • Nephrology (213)
  • Neurology (1682)
  • Nursing (97)
  • Nutrition (252)
  • Obstetrics and Gynecology (328)
  • Occupational and Environmental Health (451)
  • Oncology (930)
  • Ophthalmology (264)
  • Orthopedics (102)
  • Otolaryngology (172)
  • Pain Medicine (114)
  • Palliative Medicine (40)
  • Pathology (253)
  • Pediatrics (538)
  • Pharmacology and Therapeutics (254)
  • Primary Care Research (209)
  • Psychiatry and Clinical Psychology (1775)
  • Public and Global Health (3854)
  • Radiology and Imaging (626)
  • Rehabilitation Medicine and Physical Therapy (320)
  • Respiratory Medicine (521)
  • Rheumatology (208)
  • Sexual and Reproductive Health (168)
  • Sports Medicine (158)
  • Surgery (191)
  • Toxicology (36)
  • Transplantation (101)
  • Urology (76)