Next Article in Journal
Adeno-Associated Vector-Delivered CRISPR/SaCas9 System Reduces Feline Leukemia Virus Production In Vitro
Next Article in Special Issue
Intracellular Life Cycle Kinetics of SARS-CoV-2 Predicted Using Mathematical Modelling
Previous Article in Journal
Genetic Characterization of Highly Pathogenic Avian Influenza A(H5N8) Virus in Pakistani Live Bird Markets Reveals Rapid Diversification of Clade 2.3.4.4b Viruses
Previous Article in Special Issue
Comparative Computational Modeling of the Bat and Human Immune Response to Viral Infection with the Comparative Biology Immune Agent Based Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets

by
Naveen K. Vaidya
1,2,3,*,
Angelica Bloomquist
1,2,3 and
Alan S. Perelson
4
1
Department of Mathematics and Statistics, San Diego State University, San Diego, CA 92182, USA
2
Computational Science Research Center, San Diego State University, San Diego, CA 92182, USA
3
Viral Information Institute, San Diego State University, San Diego, CA 92182, USA
4
Los Alamos National Laboratory, Theoretical Biology and Biophysics Group, Los Alamos, NM 87545, USA
*
Author to whom correspondence should be addressed.
Submission received: 22 June 2021 / Revised: 27 July 2021 / Accepted: 9 August 2021 / Published: 18 August 2021
(This article belongs to the Special Issue Mathematical Modeling of Viral Infection)

Abstract

:
The pre-clinical development of antiviral agents involves experimental trials in animals and ferrets as an animal model for the study of SARS-CoV-2. Here, we used mathematical models and experimental data to characterize the within-host infection dynamics of SARS-CoV-2 in ferrets. We also performed a global sensitivity analysis of model parameters impacting the characteristics of the viral infection. We provide estimates of the viral dynamic parameters in ferrets, such as the infection rate, the virus production rate, the infectious virus proportion, the infected cell death rate, the virus clearance rate, as well as other related characteristics, including the basic reproduction number, pre-peak infectious viral growth rate, post-peak infectious viral decay rate, pre-peak infectious viral doubling time, post-peak infectious virus half-life, and the target cell loss in the respiratory tract. These parameters and indices are not significantly different between animals infected with viral strains isolated from the environment and isolated from human hosts, indicating a potential for transmission from fomites. While the infection period in ferrets is relatively short, the similarity observed between our results and previous results in humans supports that ferrets can be an appropriate animal model for SARS-CoV-2 dynamics-related studies, and our estimates provide helpful information for such studies.

1. Introduction

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a novel coronavirus that causes the infectious disease COVID-19 [1,2,3,4]. The ongoing global pandemic of COVID-19 has spread to almost all countries globally, with more than 178 million confirmed cases and more than 3.8 million deaths as of 21 June 2021 [5]. While tremendous efforts have been put into the control of COVID-19 outbreaks, these have mainly been non-pharmaceutical. A number of different vaccines have been approved. However, the disease is still spreading because of vaccine hesitancy, the lack of availability of vaccines in most parts of the world, and the difficulties of getting vaccines into arms when it is available. The emergence of SARS-CoV-2 variants amplifies already existing substantial threats to global public health [6,7].
The viral kinetics of SARS-CoV-2 during infection within an individual is poorly understood. The viral load data from human infection is rarely available during the early infection period [8,9,10,11,12]. In the past, animal models have successfully substituted for human hosts to study the within-host dynamics of many viruses, such as HIV, HCV, and HBV [13,14,15,16,17]. In particular, respiratory viruses that infect humans have been widely studied using ferrets [18,19,20,21,22,23,24,25,26], and recent studies have shown that SARS-CoV-2 can successfully transmit and replicate in ferrets [1,27,28,29,30,31]. As in previous studies [32,33,34,35,36,37,38,39], mathematical models fit to data from animals can be a powerful tool to characterize within-host viral kinetics of SARS-CoV-2. Such information is helpful to design, evaluate, and identify antiviral agents that can control SARS-CoV-2 infection within COVID-19 patients [1].
There are limited studies on within-host SARS-CoV-2 viral dynamics [8,9,40,41,42,43,44,45,46,47]. While previous studies have provided important insights into the viral dynamics [40,42,43,44], immune responses [8,9,45], and potential antiviral therapy [41,48] within human hosts [8,9,40,41,42,44,45] and macaques [40,43], many aspects of within-host dynamics still remain uncertain. Furthermore, most of the data used in these studies were obtained post-symptom onset, and none of these models dealt with the dynamics within the ferret, which has been considered one of the important animal models for experimental studies of SARS-CoV-2 [1,27,28,29,30,31].
In this study, we used mathematical models to characterize SARS-CoV-2 infection dynamics in two groups of ferrets [1]: one infected with SARS-CoV-2/F13/environment/2020/Wuhan isolated from an environmental sample collected in the Huanan Seafood Market in Wuhan (F13-E) and another isolated from a person infected with SARS-CoV-2/CTan/human/2020/Wuhan (CTan-H). Specifically, we used the data from SARS-CoV-2-infected ferrets [1] in our models to estimate viral kinetic parameters, such as the infection rate, the viral clearance rate, the infectious virus proportion, the infected cell death rate, and the viral production rate. In addition, we further used our model to calculate the basic reproduction number, pre-peak infectious viral growth rate, post-peak infectious viral decay rate, pre-peak infectious viral doubling time, post-peak infectious virus half-life, and the total cell loss in the respiratory tract. Our results provide evidence-based quantitative insights into within-host viral dynamics in animal models that can significantly benefit the control of SARS-CoV-2 infection in humans.

2. Materials and Methods

2.1. Experimental Data

The data used in this study were obtained from digitizing a published experimental infection of ferrets by SARS-CoV-2 [1]. As mentioned above, two different SARS-CoV-2 viruses were considered: F13-E (isolated from an environmental sample) and CTan-H (isolated from a human patient). Six ferrets (three animals in each group) were inoculated intranasally with 10 5 plaque-forming unit (PFU) of F13-E virus (animal: F13-E-1, F13-E-2, and F13-E-3) and CTan-H virus (animal: CTan-H-1, CTan-H-2, and CTan-H-3) in a volume of 1 mL. Viral RNA copies and viral titer (infectious virus) were recorded in nasal washes collected on days 2, 4, 6, 8, and 10 post infection (p.i.) from each animal. For our modeling, the data lying below the lower limit of detection (i.e., 3 l o g 10 viral RNA copies per mL and 1 l o g 10 PFU per mL) were taken as half of the limit of detection.

2.2. Mathematical Model

To model the data containing viral RNA copies and viral titer (infectious virus), we used three viral dynamic models (Model 1: basic viral dynamics; Model 2: viral dynamics with eclipse phase; and Model 3: viral dynamics with immune response), similar to the ones for influenza [36,49], but we considered two types of virus population: infectious ( V i ) and non-infectious ( V n ). As SARS-CoV-2 was found to primarily replicate in the upper respiratory tract (nasal turbinate, soft palate, tonsils) of ferrets [1], we considered cells in the upper respiratory tract as targets of SARS-CoV-2 infection. Infectious SARS-CoV-2, V i , infects target cells, T , at rate β T V i , where β is a rate   constant . Infected cells, I , are assumed to die at per capita rate δ and produce new free virus particles at rate p per cell. A portion, α , of the newly produced free viruses are assumed to be infectious, and the remaining 1 α free viruses are non-infectious. The parameter c represents the clearance rate of both infectious and non-infectious viruses.
The system of equations in Model 1 is
d T d t = β T V i ,             T 0   =   T 0 , d I d t = β T V i δ I ,             I 0   =   I 0 ,   d V i d t = α p I c V i ,                 V i 0   =   V i 0 , d V n d t =   1 α p I c V n ,               V n 0   =   V n 0 .
In Model 2, we consider two classes of infected cells: one in the eclipse phase, I e , which do not produce virus yet, and another in the productive phase, I p , which are actively producing virus. The newly infected cells remain in the eclipse phase for an average duration of 1 / k . With these mechanisms, the system of equations in Model 2 becomes
d T d t = β T V i ,             T 0   =   T 0 , d I e d t = β T V i k I e ,             I e 0   =   I e 0 , d I p d t = k I e δ I p ,             I p 0   =   I p 0 ,   d V i d t = α p I p c V i ,                 V i 0   =   V i 0 , d V n d t =   1 α p I p c V n ,               V n 0   =   V n 0 .
To develop Model 3, we extend Model 1 by incorporating an immune response. There are many potential immune responses, and incorporating all of the effects due to immune responses requires highly complex models. However, as a representative of immune response-models, we consider a simple approach because of limited data. SARS-CoV-2-specific adaptive immune response has not been properly documented in ferrets during the infection period (i.e., 10 days). Antibodies against SARS-CoV-2 in these ferrets have been measured, but only at day 13 and day 20 post infection [1]. Therefore, we consider the nonspecific innate immune response, which provides the first-line defense [50]. Specifically, we include type I Interferon (mainly IFN- α / β ), which induces the expression of many IFN-stimulated antiviral proteins in the neighboring cells, making them refractory to infection [50]. We follow the modeling idea used by Pawelek et al. [36] to model such immune effect. We represent the IFN level by F and assume IFN makes uninfected cells refractory to infection with rate σ . The IFN level grows in proportion to the infected cells at rate g and decays at the per capita rate of ω . For simplicity, we assume the cells in the refractory state do not revert to the susceptible state in this relatively short infection period.
The system of equations representing Model 3 is
d T d t = β T V i σ T F ,             T 0   =   T 0 , d I d t = β T V i δ I ,             I 0   =   I 0 ,   d V i d t = α p I c V i ,                 V i 0   =   V i 0 , d V n d t =   1 α p I c V n ,               V n 0   =   V n 0 ,   d F d t =   g I ω F ,               F 0   =   F 0 .

2.3. Parameters and Data Fitting

Cells in the upper respiratory tract are the main target for SARS-CoV-2 infection [1]. Based on the surface area per epithelial cell and the total area of epithelial cells lining the nasal turbinate of the human upper respiratory tract, Baccam et al. [49] calculated the total cells in the upper respiratory tract of humans to be 4 × 10 8 cells. For simplicity, we scaled this value based on the ratio of the average weight of a human (75 kg) to the average weight of a ferret (2 kg) and thus assume there are 10 7 epithelial cells in ferret upper respiratory tract. Only a fraction of these cells expresses the ACE2 receptor and enzyme TMPRSS2 needed to cleave the spike protein so that the virus can fuse to the cell membrane and successfully enter and infect a cell. Estimates of this fraction range from 1% [51] to 20% [52]. We will assume 10% and take T 0 = 10 6 cells. However, we note that the number of initial target cells chosen, T0, does not affect the results of this study because redefining T T / T 0 ,   I I / T 0 ,   I e I e / T 0 ,   I p I p / T 0 , and p p T 0 will keep the systems unaltered, implying that the only parameter affected by a different choice of T0 is p. Based on a previous study [41], we set g = 1 per day and ω = 0.4 per day [41].
As infection was established with a virus inoculum, we set I 0 = I e 0 = I p 0 = 0 . To estimate the initial viral load inoculated, we obtained the regression using viral load and corresponding virus titer from all animals at all time points. Then, using initial virus titer of 10 5 PFU in the regression formula, we obtained 9 log10 viral RNA copies, which we set to be the initial virus inoculated. The remaining parameters ( V i 0 , β ,   δ , α , p , c ), ( V i 0 , β ,   δ , α , p , c , k ), and ( V i 0 , β ,   δ , α , p , c , σ ) in Model 1, 2, and 3, respectively, were estimated by fitting the models to experimental data from each ferret individually. We solved the model equations numerically using the built-in solver ode45 in MATLAB 2020a (The MathWorks, Inc.). The ode45 is based on an explicit Runge–Kutta (4,5) formula, the Dormand–Prince pair. The predicted log10 values of the total viral load ( V i + V n ) and the infectious virus concentration ( V i ) were fitted to the corresponding log-transformed data via a nonlinear least square regression method in which the sum of the squared residuals (SSR), i.e., the difference between the model predictions and the corresponding experimental values, was minimized. We used the following formula to calculate SSR:
S S R = t d = 1 N d [ { log 10 V 2 t d     log 10 V ¯ 2 t d } 2   +   { log 10 V i 2 t d     log 10 V ¯ i 2 t d } 2 ] ,
where V t d   =   V i t d   +   V n t d represents the viral load at time t d predicted by the model, V i t d represents the infectious virus concentration at time t d predicted by the model, and V ¯ t d and V ¯ i t d are the values of viral RNA copies and the viral titer, respectively, in the experimental measurement. N d = 5 is the total number of time points at which measurements were taken, and the sum is taken over the set of measurement time points. In this study, a total of 10 data points, including two different types (viral load and infectious virus) was available for each fitting.
We also performed the global minimum search by providing various initial guesses distributed uniformly across the reasonable parameter limits. To overcome any possible stiffness of the model equations, we performed data fitting using ode15s, which is an ODE solver designed for stiff equations, and we did not find any difference in our fitting results. For each best-fit parameter, we also provide 95% confidence intervals (CI), which were computed from 250 replicates, by bootstrapping the residuals [53,54]. Since parameter values were constrained to be positive, potentially causing their distributions to be skewed to the right [49], we computed the geometric mean to report the average parameter values.
We also performed an analysis to determine whether the model parameters we expect to estimate are uniquely identifiable for the available viral load and infectious virus data. We note that while parameters may be correlated for the basic TIV (Target cell, Infected cell, and Virus particles) model [55], our model is a generalization that includes infectious virus, and we fitted to both the total viral load as in the TIV model as well as to the measured infectious titer. Thus, we used more data than in the TIV model. To address the identifiability issue, we generated the new data at the time of sampling from the model simulations with estimated known parameters. Then, we fitted the generated data with our fitting routine using different initial guesses as in our global minimization calculations. We found that the parameters, including both p and β , used to generate simulated data were recovered by our fitting procedure. We also repeated the process by adding 5%, 10%, 15%, and 20% random errors to the generated data. In this case, we obtained estimated values of the parameters approximately close to what we used for data generations but with slight variations, as expected. In addition, to quantify the correlation of parameters fitted using our data sets with two different natures (viral load and infectious virus), we also computed Pearson correlation coefficients for each model and each animal using the parameters obtained from replicates by bootstrapping the residuals. Our results show that the parameters estimated in our case are not correlated (Supplementary Table S1). Therefore, we do not expect identifiability issues in our fitting, while we acknowledge some uncertainty due to limitations on the models and data.

2.4. Sensitivity Analysis

To study how the model dynamics are affected by uncertainty present in model parameters, we considered a widely varying parameter space for each model. With these parameter spaces, we performed a global sensitivity analysis using Latin Hypercube sampling for each of the three models.

3. Results

3.1. Data Fitting to Three Viral Dynamics Models and Comparison

For each ferret, we fitted each of the three models to the experimental data. In the fitting of Model 1, we estimated the initial infectious virus concentration ( V i 0 ), the infection rate constant ( β ), the death rate of infected cell ( δ ), the proportion of infectious virus among newly produced virus ( α ), the virus production rate ( p ), and the virus clearance rate ( c ). In Models 2 and 3, we estimated the additional parameters k (the rate of transfer from eclipse phase to infectious phase) and σ (the IFN-induced antiviral efficacy), respectively, along with those parameters estimated in Model 1. The best estimates of these parameters, along with 95% confidence intervals and geometric mean across ferrets, are given in Table 1 (Model 1), Table 2 (Model 2), and Table 3 (Model 3). With these parameters, the model predictions agree well with the experimental data on both viral RNA copies and viral titer (Figure 1—Model 1, Figure 2—Model 2, Figure 3—Model 3).
To examine which model can better fit, we computed the Sum of the Squared Residuals (SSR) and the Akaike Information Criterion (AIC) values. As discussed in the method section, SSR indicates how close the model prediction is to the data, and AIC allows one to compare models given the goodness of fit for a given model, number of data points, and number of parameters. Our model comparison shows that the basic model is approximately as good as the other two extended models (Model 1: 4.24 (SSR), 61.36 (AIC); Model 2: 4.14 (SSR), 151.16 (AIC); Model 3: 3.88 (SSR), 150.52 (AIC)), and thus the additional parameters k and σ , introduced in Models 2 and 3, respectively, did not improve the data fitting. Since the models provide similar dynamics, and the basic model has the best (lowest) AIC value, we present most of the results to follow based on Model 1. However, we observed that Models 2 and 3 also provide a similar conclusion about the viral dynamic properties.
The model fits revealed a trend in SARS-CoV-2 dynamics within a ferret in which the viral load initially decayed for about 17 h and then increased until it reached a peak around 5 days post infection. However, the concentration of infectious virus increased monotonically until it reached a peak. After the peak, both viral load and infectious virus decreased to a level below the limit of detection at around 8 days post infection. At the peak, the viral load reached about 8 log10 RNA copies per mL, among which ~4 log10 (~0.01%) are infectious viruses. The general pattern observed in SARS-CoV-2 dynamics is similar to that for influenza virus [36,49].

3.2. Important Viral Dynamics Parameters

The average estimate for the infection rate, β , was 6.93 × 10 5 per infectious virus per day. As indicated by our estimate of α = 1.75 × 10 4 , only about 0.02% of the total newly produced virus were infectious, and most of the viral particles were non-infectious. This low estimate of the newly produced infectious virus proportion is consistent with other data [56]. We also obtained that the infection in ferret was initiated with about 5 infectious viruses per mL on average (i.e., V i 0   ~   5 /mL), while the experiment was initiated with the intranasal inoculation of 105 PFU of stock viruses in 1 mL. Therefore, only a very few infectious viruses may be needed at the infection site (upper respiratory tract) to initiate and establish SARS-CoV-2 infection.
The infected cell death rate, δ , and virus clearance rate, c , estimated by our model were 5.20 per day and 6.10 per day, respectively. These estimates indicate that infected cells turned over in about 5 h on average, and the free virus gets cleared from the body in about 4 h on average. Using Model 2, we estimated k (the rate of cell transfer from the eclipse phase to the productively infected phase) to be 7.88 per day, which implies that the infected cells started producing virus about 3 h after they were infected. Similarly, we used Model 3 to estimate the rate constant for cells becoming refractory due to IFN as σ = 1.97 × 10 6 per day.

3.3. Basic Reproduction Number

The basic reproduction number, R 0 , represents the average number of infectious viruses generated by a single infectious virus when target cells are not limiting [57,58,59]. For Model 1, R 0 can be calculated using the formula [49]
R 0 = α p β T 0 c δ .
The basic reproduction number, R 0 , for SARS-CoV-2 estimated by our modeling varies among individual ferrets from 1.52 to 2.69, with an average value of 2.05 (Table 4). For within-host virus dynamics, this information can be applied to determine whether a virus can establish infection [57,58,59]. In general, if R 0 < 1 , the infection will die out, and if R 0 > 1 , the infection will spread [57,58,59]. Consistent with this theoretical result, our estimate also showed each R 0 to be greater than 1 as the infection was established in each ferret. In addition, R 0 estimate can also be applied to calculate the effectiveness of antiviral agents, such as a vaccine, required to stop the viral replication and avoid infection. For example, the average value of R 0 = 2.05 indicates that at least 51%, i.e., ( 1 1 / R 0 ) , efficacy is required for antiviral agents blocking viral entry or production to stop the infection.

3.4. Viral Kinetic Properties

We also used our model to estimate other kinetic properties, including pre-peak growth rate and doubling time as well as post-peak decay rate and half-life of infectious virus (Table 4), during the SARS-CoV-2 infection in ferrets. During the pre-peak phase, the concentration of infectious virus grew exponentially, with a per capita growth rate of 2.29 per day (Table 4), resulting in a doubling time of 7.27 h. Similarly, during the post-peak phase, the infectious virus decayed with a per capita decay rate of 1.95 per day (Table 4) and a half-life of 8.54 h. In addition to these rates, we also estimated the net target cell loss in the upper respiratory tract of ferrets due to SARS-CoV-2 infection (Table 4). On average, approximately 78% of target cells were lost in the upper respiratory tract during the entire infectious period of SARS-CoV-2 in ferrets.

3.5. Comparison between F13-E and CTan-H Viruses

We compared the viral dynamics generated due to infection by the two types of viruses, F13-E and CTan-H. As discussed above, these virus types represent viruses in the environment (seafood market of Wuhan) and viruses within an infected person, respectively. Successful establishment of infection by both viruses ( R 0 > 1 in both cases) implies that environmental transmission may also be an important route in addition to direct transmission. We observed a similar trend of viral dynamics in both cases (Figure 1), and the estimates of the parameters did not show any statistically significant difference between the two virus infections (t-test, p-value 0.05 ). Furthermore, a comparison of the kinetic parameters in all three models (Table 1, Table 2 and Table 3) revealed that none of the parameters were statistically different between two viruses (t-test, p-value 0.05 ). These results combined indicate that the environmental virus (F13-E) was as infectious as the within-human virus (CTan-H), causing similar cell loss and other properties of infection.

3.6. Sensitivity Analysis

Even though we estimated model parameters using experimental data, there is uncertainty among the estimated parameters, and the limited number of animals used in the experiment may not be representative of a larger population. Therefore, we performed a global sensitivity analysis of each model parameter using the Latin Hypercube sampling technique. Specifically, we computed partial rank correlation coefficients (PRCC) for each parameter of all three models corresponding to total viral-peak, infectious viral-peak, time to viral-peak, and infected target cell loss (Figure 4). In general, the infection rate ( β ), the virus production rate ( p ), the proportion of infectious virus ( α ), and the rate of transfer from eclipse phase to infectious phase ( k ) positively impacted total viral-peak, infectious viral-peak, and cell loss and negatively impacted time to viral-peak. The viral clearance rate ( c ) , the infected cell death rate ( δ ), and the IFN-induced antiviral efficacy ( σ ) have both positive and negative effects depending on the property being examined (Figure 4). As observed from all the models, compared to the basic viral dynamic parameters ( c ,   p ,   α ,   δ ,   β ), the eclipse phase related parameter ( k ) and the immune-response-related parameter ( σ ) have significantly lower impact on the viral dynamics (Figure 4), supporting the results above that additional parameters in eclipse phase model (Model 2) and immune response model (Model 3) did not improve the data-fit compared to the basic viral dynamics model (Model 1).

3.7. Initial Infectious Virus

As mentioned earlier, animals were inoculated intranasally with ~ 10 5 PFU in a volume of 1 mL, but our estimates showed a relatively low level of initial infectious virus at the infection site (8.18 PFU/mL for F-13E group and 2.46 PFU/mL for CTan-H group). To understand this discrepancy, we performed further data fitting by fixing the initial infectious virus population at the infection site (upper respiratory tract) at V i 0 = 10 2 ,   10 3 , 10 4 PFU/mL (Figure 5). We found that a low level of initial infectious virus is needed to obtain reasonable fits to the data (F-13E-2: SSR = 3.21, 3.99, 7.65, 10.65, AIC = 28.64, 30.80, 37.32, 40.63, for V i 0 = 9.25 ,   10 2 ,   10 3 , 10 4 PFU/mL, respectively; CTan-H-2: SSR = 2.99, 3.64, 5.79, 8.21, AIC = 27.92, 29.89, 34.53, 38.03, for V i 0 = 9.12 ,   10 2 ,   10 3 , 10 4 PFU/mL, respectively). There may be a number of factors that have contributed to this discrepancy. The viral stocks used for inoculation were prepared in Vero E6 cells and the measurement of their infectiousness was done by plaque assay using Vero E6 cells [1]. Vero E6 cells are a stable cell line derived from African green monkey kidney cells. The cells have lost the ability to secrete IFN due to spontaneous gene deletions [60,61], indicating Vero cells should be more susceptible to infection than the target cells in the ferret. Thus, PFUs measured in Vero cells may be considerably greater than the corresponding infectivity in ferrets. Additionally, the infectiousness at the infection site within ferrets may be different from the one estimated in in-vitro experiments, as the virus may adapt to the in-vivo cell environment [62]. Lastly, we note that the infection of SARS-CoV-2 has been successfully established in ferrets with widely varying initial infectious doses and the infection has been reported to depend on the route of transmission [63,64].

4. Discussion

Despite the development of vaccines against SARS-CoV-2, the threat of the COVID-19 pandemic remains. Thus, there is a continuing and urgent need to administer effective antiviral agents for SARS-CoV-2 treatment. Animal models are often used for experimental trials, while developing antiviral agents and ferrets are considered suitable hosts for respiratory diseases, including SARS-CoV-2. Here, we used mathematical models to characterize SARS-CoV-2 infection within ferrets infected with F13-E (SARS-CoV-2 isolated from an environmental sample) or infected with CTan-H (SARS-CoV-2 isolated from an infected individual). Our modeling results, including a global sensitivity analysis, have provided important insights into the viral dynamics and viral kinetics in ferrets and allowed us to compare them with infection in humans. Such information is useful for performing experiments in animals to design antiviral therapies for humans.
The temporal pattern of SARS-CoV-2 dynamics in ferrets is similar to the infection in humans [10,40,41,43,44] as well as to that of some other acute viral infections, such as influenza [36,49], though the infection in ferrets lasts for a significantly shorter duration than in humans (~10 days in ferrets [1] vs >20 days in humans [8,41,44]). Despite the short duration of infection, many of the viral dynamics parameters we estimated in ferrets, including infection rate, virus production rate, viral clearance rate, infected cell death rate, duration in eclipse phase, and IFN-induced antiviral efficacy, are similar to those in humans [8,40,44]. This similarity indicates that ferrets can be an appropriate animal model for experimental trials related to viral dynamics properties. It is worth noting that while we estimated the eclipse phase length and IFN-induced antiviral activities using extended models, our model comparison based on data fitting and the global sensitivity analysis showed that these additional parameters (impact of eclipse phase and IFN-induced antiviral activities) were not significantly impactful for infection dynamics in ferrets as compared to other basic viral dynamics parameters. Minimal effects of SARS-CoV-2-specific antibodies were also found in a previous study [9].
Because of the similarity in the basic viral dynamics parameters between SARS-CoV-2 in ferrets and humans and in other acute viral infections, such as influenza, the basic reproduction number and other viral kinetic properties (pre-peak growth rate, doubling time, post-peak decay rate, and half-life of infectious virus) are in a similar range [41,44]. However, the kinetic properties are quite different from other viral infections, such as HBV and EIAV, which result in chronic infection in the host [38,39]. For example, our estimate of the basic reproduction number in SARS-CoV-2 infection in ferrets ( R 0 = 2.05 ) is lower than HBV infection in chimpanzee ( R 0 > 15 ) [39] and EIAV infection in horses ( R 0 > 15 ) [38]. Additionally, in EIAV [38], the infection cell death rate is negligible (~0.06 per day), resulting in only 0.7% cell loss. In contrast, SARS-CoV-2 infection causes significantly higher cell death (~5.20 per day) in ferrets, comparable to human influenza (~4 per day [49]), resulting in approximately 78% target cell loss in the upper respiratory tract of ferrets.
We also note that the initial decline in total virus observed in our simulations using the best-fit parameters is significantly larger compared to declines predicted by model fits to data with other viruses such as influenza. In the experiment modeled here, ferrets were infected with 105 PFU in 1 mL of fluid [1]. Comparing the early measurements of the total virus with the measurements of infectious virus, we deduced that 1 PFU corresponds to approximately 104 virions, which is in line with published estimates [56]; thus, we fixed V 0 at 109 per mL. The first measured values were at day 2, and by then, viral loads were of order 106 per mL, suggesting that there must be a dramatic decline in the viral load over the first two days, which our model replicates. The viral load decline very early in infection can be better estimated if more frequent sampling is done, especially within the first two days post-infection.
Our comparative analysis of viral dynamics initiated with a virus isolated from the environment and from a human host shows that both viruses can establish infection within the ferret, and the viral dynamics parameters in both cases are similar. This suggests the possibility of SARS-CoV-2 transmission from fomites. In both cases, we estimated that only about 0.02% of the newly produced virus was infectious, and the remaining portion was non-infectious. In addition, our estimate of V i 0 and our sensitivity analysis suggest that a small number of infectious SARS-CoV-2 at the infection site (upper respiratory tract) initiated infection despite a higher virus titer in the inoculum. These results indicate that ferrets need to be exposed to a large amount of virus for successful infections, consistent with previous studies, in which a small viral dose was not sufficient to infect all the inoculated animals in their experiment [63].
We acknowledge some limitations of our study. Our parameter estimates are based on a limited number of measurements from a small number of animals. Therefore, our estimates of parameters may have uncertainty that may arise from variation across animals. However, we performed a global sensitivity analysis using a wide parameter space, which suggests robustness of the viral dynamics and viral kinetic properties. While our fitting did not support the need to include an immune response in the model, the anti-SARS-CoV-2 antibody only measured for the viral level was undetectable in these ferrets, i.e., after 13 days or 20 days. More data with frequent measurements of immune responses can help improve our modeling and discern more accurately if an immune response is needed to clear this virus in ferrets. We also note that the parameter estimates we obtained may not reflect in viral kinetics in humans quantitatively because of extremely short-time dynamics in ferrets compared to humans. More focus on the qualitative conclusion of our study may be appropriate to predict human infection until these models are extensively tested in human infection with frequent pre-peak data. The extension of our modeling study by including potential antiviral therapy in ferrets can be a future avenue in this research direction.
We recognize that the viral strains used in the study we modeled are not the current variants of concern and that it would be valuable to model the viral kinetics of these strains. However, new variant viruses, such as Alpha (B.1.1.7) and Beta (B.1.351), are being studied [30,31], and when appropriate data is available, our models fit to the data may provide insights into their within-host kinetics.
In summary, we implemented three different viral dynamic models to analyze the data from SARS-CoV-2 infection in ferrets. Estimated parameters, predicted viral kinetic properties, and the global sensitivity analysis over a wide parameter space included in this study may provide useful information for the development of therapies utilizing ferrets as animal models for experimental trials.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/v13081635/s1, Table S1: Pearson Correlation Coefficient of Estimated Parameters from Replicates by Bootstrapping the Residuals.

Author Contributions

N.K.V.: conceptualization, methodology, formal analysis, investigation, writing—original draft preparation, writing—review and editing; A.B.: methodology, formal analysis, investigation; A.S.P.: investigation, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NSF grants DMS-1951793, DMS-1616299, DMS-1836647, and DEB-2030479 (NKV) and NIH grants R01-AI028433 and R01-OD011095 (ASP) and the UGP award from San Diego State University (NKV). The funders had no role in study design, analysis, decision to publish, or preparation of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The study did not report any data. Only digitized published data were used.

Acknowledgments

This work was funded by NSF grants DMS-1951793, DMS-1616299, DMS-1836647, and DEB-2030479 (NKV) and NIH grants R01-AI028433 and R01-OD011095 (ASP) and the UGP award from San Diego State University (NKV). Portions of this work were under the auspices of the U.S. Department of Energy under contract 89233218CNA000001. Authors would like to thank Forest Rohwer for useful discussion.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Shi, J.; Wen, Z.; Zhong, G.; Yang, H.; Wang, C.; Huang, B.; Liu, R.; He, X.; Shuai, L.; Sun, Z.; et al. Susceptibility of ferrets, cats, dogs, and other domesticated animals to SARS-coronavirus 2. Science 2020, 368, 1016–1020. [Google Scholar] [CrossRef] [Green Version]
  2. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.M.; Wang, W.; Song, Z.G.; Hu, Y.; Tao, Z.W.; Tian, J.H.; Pei, Y.Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
  3. Chan, J.F.; Kok, K.H.; Zhu, Z.; Chu, H.; To, K.K.; Yuan, S.; Yuen, K.Y. Genomic characterization of the 2019 novel human-pathogenic coronavirus isolated from a patient with atypical pneumonia after visiting Wuhan. Emerg. Microbes Infect. 2020, 9, 221–236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Zou, L.; Ruan, F.; Huang, M.; Liang, L.; Huang, H.; Hong, Z.; Yu, J.; Kang, M.; Song, Y.; Xia, J.; et al. SARS-CoV-2 Viral Load in Upper Respiratory Specimens of Infected Patients. N. Engl. J. Med. 2020, 382, 1177–1179. [Google Scholar] [CrossRef]
  5. WHO. WHO Coronavirus (COVID-19) Dashboard; WHO: Geneva, Switzerland, 2021. [Google Scholar]
  6. Baric, R.S. Emergence of a Highly Fit SARS-CoV-2 Variant. N. Engl. J. Med. 2020, 383, 2684–2686. [Google Scholar] [CrossRef] [PubMed]
  7. Lauring, A.S.; Hodcroft, E.B. Genetic Variants of SARS-CoV-2-What Do They Mean? JAMA 2021, 325, 529–531. [Google Scholar] [CrossRef] [PubMed]
  8. Goyal, A.; Cardozo-Ojeda, E.F.; Schiffer, J.T. Potency and timing of antiviral therapy as determinants of duration of SARS-CoV-2 shedding and intensity of inflammatory response. Sci. Adv. 2020, 6, eabc7112. [Google Scholar] [CrossRef] [PubMed]
  9. Yang, S.; Jerome, K.R.; Greninger, A.L.; Schiffer, J.T.; Goyal, A. Endogenously Produced SARS-CoV-2 Specific IgG Antibodies May Have a Limited Impact on Clearing Nasal Shedding of Virus during Primary Infection in Humans. Viruses 2021, 13, 516. [Google Scholar] [CrossRef] [PubMed]
  10. Kim, J.Y.; Ko, J.H.; Kim, Y.; Kim, Y.J.; Kim, J.M.; Chung, Y.S.; Kim, H.M.; Han, M.G.; Kim, S.Y.; Chin, B.S. Viral Load Kinetics of SARS-CoV-2 Infection in First Two Patients in Korea. J. Korean Med. Sci. 2020, 35, e86. [Google Scholar] [CrossRef]
  11. Pan, Y.; Zhang, D.; Yang, P.; Poon, L.L.M.; Wang, Q. Viral load of SARS-CoV-2 in clinical samples. Lancet Infect. Dis. 2020, 20, 411–412. [Google Scholar] [CrossRef]
  12. Wölfel, R.; Corman, V.M.; Guggemos, W.; Seilmaier, M.; Zange, S.; Müller, M.A.; Niemeyer, D.; Jones, T.C.; Vollmar, P.; Rothe, C.; et al. Virological assessment of hospitalized patients with COVID-2019. Nature 2020, 581, 465–469. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Berggren, K.A.; Suzuki, S.; Ploss, A. Animal Models Used in Hepatitis C Virus Research. Int. J. Mol. Sci. 2020, 21, 3869. [Google Scholar] [CrossRef] [PubMed]
  14. Ma, Z.M.; Stone, M.; Piatak, M.; Schweighardt, B.; Haigwood, N.L.; Montefiori, D.; Lifson, J.D.; Busch, M.P.; Miller, C.J. High specific infectivity of plasma virus from the pre-ramp-up and ramp-up stages of acute simian immunodeficiency virus infection. J. Virol. 2009, 83, 3288–3297. [Google Scholar] [CrossRef] [Green Version]
  15. Asabe, S.; Wieland, S.F.; Chattopadhyay, P.K.; Roederer, M.; Engle, R.E.; Purcell, R.H.; Chisari, F.V. The size of the viral inoculum contributes to the outcome of hepatitis B virus infection. J. Virol. 2009, 83, 9652–9662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Kumar, R.; Orsoni, S.; Norman, L.; Verma, A.S.; Tirado, G.; Giavedoni, L.D.; Staprans, S.; Miller, G.M.; Buch, S.J.; Kumar, A. Chronic morphine exposure causes pronounced virus replication in cerebral compartment and accelerated onset of AIDS in SIV/SHIV-infected Indian rhesus macaques. Virology 2006, 354, 192–206. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Kumar, R.; Torres, C.; Yamamura, Y.; Rodriguez, I.; Martinez, M.; Staprans, S.; Donahoe, R.M.; Kraiselburd, E.; Stephens, E.B.; Kumar, A. Modulation by morphine of viral set point in rhesus macaques infected with simian immunodeficiency virus and simian-human immunodeficiency virus. J. Virol. 2004, 78, 11425–11428. [Google Scholar] [CrossRef] [Green Version]
  18. Maher, J.A.; DeStefano, J. The ferret: An animal model to study influenza virus. Lab. Anim. 2004, 33, 50–53. [Google Scholar] [CrossRef] [Green Version]
  19. Shi, J.; Deng, G.; Kong, H.; Gu, C.; Ma, S.; Yin, X.; Zeng, X.; Cui, P.; Chen, Y.; Yang, H.; et al. H7N9 virulent mutants detected in chickens in China pose an increased threat to humans. Cell Res. 2017, 27, 1409–1421. [Google Scholar] [CrossRef] [Green Version]
  20. Stittelaar, K.J.; de Waal, L.; van Amerongen, G.; Veldhuis Kroeze, E.J.; Fraaij, P.L.; van Baalen, C.A.; van Kampen, J.J.; van der Vries, E.; Osterhaus, A.D.; de Swart, R.L. Ferrets as a Novel Animal Model for Studying Human Respiratory Syncytial Virus Infections in Immunocompetent and Immunocompromised Hosts. Viruses 2016, 8, 168. [Google Scholar] [CrossRef]
  21. Zhang, Q.; Shi, J.; Deng, G.; Guo, J.; Zeng, X.; He, X.; Kong, H.; Gu, C.; Li, X.; Liu, J.; et al. H7N9 influenza viruses are transmissible in ferrets by respiratory droplet. Science 2013, 341, 410–414. [Google Scholar] [CrossRef]
  22. Imai, M.; Watanabe, T.; Hatta, M.; Das, S.C.; Ozawa, M.; Shinya, K.; Zhong, G.; Hanson, A.; Katsura, H.; Watanabe, S.; et al. Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature 2012, 486, 420–428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Herfst, S.; Schrauwen, E.J.; Linster, M.; Chutinimitkul, S.; de Wit, E.; Munster, V.J.; Sorrell, E.M.; Bestebroer, T.M.; Burke, D.F.; Smith, D.J.; et al. Airborne transmission of influenza A/H5N1 virus between ferrets. Science 2012, 336, 1534–1541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Herfst, S.; Mok, C.K.P.; van den Brand, J.M.A.; van der Vliet, S.; Rosu, M.E.; Spronken, M.I.; Yang, Z.; de Meulder, D.; Lexmond, P.; Bestebroer, T.M.; et al. Human Clade 2.3.4.4 A/H5N6 Influenza Virus Lacks Mammalian Adaptation Markers and Does Not Transmit via the Airborne Route between Ferrets. mSphere 2018, 3, e00405–e00417. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Van den Brand, J.M.; Haagmans, B.L.; Leijten, L.; van Riel, D.; Martina, B.E.; Osterhaus, A.D.; Kuiken, T. Pathology of experimental SARS coronavirus infection in cats and ferrets. Vet. Pathol. 2008, 45, 551–562. [Google Scholar] [CrossRef] [Green Version]
  26. Martina, B.E.; Haagmans, B.L.; Kuiken, T.; Fouchier, R.A.; Rimmelzwaan, G.F.; Van Amerongen, G.; Peiris, J.S.; Lim, W.; Osterhaus, A.D. Virology: SARS virus infection of cats and ferrets. Nature 2003, 425, 915. [Google Scholar] [CrossRef]
  27. Francis, M.E.; Richardson, B.; McNeil, M.; Rioux, M.; Foley, M.K.; Ge, A.; Pechous, R.D.; Kindrachuk, J.; Cameron, C.M.; Richardson, C.; et al. Male sex and age biases viral burden, viral shedding, and type 1 and 2 interferon responses during SARS-CoV-2 infection in ferrets. bioRxiv 2021. [Google Scholar] [CrossRef]
  28. Proud, P.C.; Tsitoura, D.; Watson, R.J.; Chua, B.Y.; Aram, M.J.; Bewley, K.R.; Cavell, B.E.; Cobb, R.; Dowall, S.; Fotheringham, S.A.; et al. Prophylactic intranasal administration of a TLR2/6 agonist reduces upper respiratory tract viral shedding in a SARS-CoV-2 challenge ferret model. EBioMedicine 2021, 63, 103153. [Google Scholar] [CrossRef]
  29. Ryan, K.A.; Bewley, K.R.; Fotheringham, S.A.; Slack, G.S.; Brown, P.; Hall, Y.; Wand, N.I.; Marriott, A.C.; Cavell, B.E.; Tree, J.A.; et al. Dose-dependent response to infection with SARS-CoV-2 in the ferret model and evidence of protective immunity. Nat. Commun. 2021, 12, 81. [Google Scholar] [CrossRef]
  30. Ryu, D.K.; Song, R.; Kim, M.; Kim, Y.I.; Kim, C.; Kim, J.I.; Kwon, K.S.; Tijsma, A.S.; Nuijten, P.M.; van Baalen, C.A.; et al. Therapeutic effect of CT-P59 against SARS-CoV-2 South African variant. Biochem. Biophys. Res. Commun. 2021, 566, 135–140. [Google Scholar] [CrossRef]
  31. Sourimant, J.; Lieber, C.M.; Aggarwal, M.; Cox, R.M.; Wolf, J.D.; Yoon, J.J.; Toots, M.; Ye, C.; Sticher, Z.; Kolykhalov, A.A.; et al. 4’-Fluorouridine is a broad-spectrum orally efficacious antiviral blocking respiratory syncytial virus and SARS-CoV-2 replication. bioRxiv 2021. [Google Scholar] [CrossRef]
  32. Vaidya, N.K.; Ribeiro, R.M.; Miller, C.J.; Perelson, A.S. Viral dynamics during primary simian immunodeficiency virus infection: Effect of time-dependent virus infectivity. J. Virol. 2010, 84, 4302–4310. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Vaidya, N.K.; Ribeiro, R.M.; Perelson, A.S.; Kumar, A. Modeling the Effects of Morphine on Simian Immunodeficiency Virus Dynamics. PLoS Comput. Biol. 2016, 12, e1005127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Vaidya, N.K.; Ribeiro, R.M.; Liu, P.; Haynes, B.F.; Tomaras, G.D.; Perelson, A.S. Correlation Between Anti-gp41 Antibodies and Virus Infectivity Decay During Primary HIV-1 Infection. Front. Microbiol. 2018, 9, 1326. [Google Scholar] [CrossRef] [PubMed]
  35. Mutua, J.M.; Perelson, A.S.; Kumar, A.; Vaidya, N.K. Modeling the effects of morphine-altered virus specific antibody responses on HIV/SIV dynamics. Sci. Rep. 2019, 9, 5423. [Google Scholar] [CrossRef]
  36. Pawelek, K.A.; Huynh, G.T.; Quinlivan, M.; Cullinane, A.; Rong, L.; Perelson, A.S. Modeling within-host dynamics of influenza virus infection including immune responses. PLoS Comput. Biol. 2012, 8, e1002588. [Google Scholar] [CrossRef] [Green Version]
  37. Ciupe, S.M.; Ribeiro, R.M.; Nelson, P.W.; Dusheiko, G.; Perelson, A.S. The role of cells refractory to productive infection in acute hepatitis B viral dynamics. Proc. Natl. Acad. Sci. USA 2007, 104, 5050–5055. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Schwartz, E.J.; Vaidya, N.K.; Dorman, K.S.; Carpenter, S.; Mealey, R.H. Dynamics of lentiviral infection in vivo in the absence of adaptive immune responses. Virology 2018, 513, 108–113. [Google Scholar] [CrossRef]
  39. Ciupe, S.M.; Vaidya, N.K.; Forde, J.E. Early events in hepatitis B infection: The role of inoculum dose. Proc. R. Soc. B 2021, 288, 20202715. [Google Scholar] [CrossRef]
  40. Wang, S.; Pan, Y.; Wang, Q.; Miao, H.; Brown, A.N.; Rong, L. Modeling the viral dynamics of SARS-CoV-2 infection. Math. Biosci. 2020, 328, 108438. [Google Scholar] [CrossRef]
  41. Gonçalves, A.; Bertrand, J.; Ke, R.; Comets, E.; de Lamballerie, X.; Malvy, D.; Pizzorno, A.; Terrier, O.; Rosa Calatrava, M.; Mentré, F.; et al. Timing of Antiviral Treatment Initiation is Critical to Reduce SARS-CoV-2 Viral Load. CPT Pharmacomet. Syst. Pharmacol. 2020, 9, 509–514. [Google Scholar] [CrossRef]
  42. Néant, N.; Lingas, G.; Le Hingrat, Q.; Ghosn, J.; Engelmann, I.; Lepiller, Q.; Gaymard, A.; Ferré, V.; Hartard, C.; Plantier, J.C.; et al. Modeling SARS-CoV-2 viral kinetics and association with mortality in hospitalized patients from the French COVID cohort. Proc. Natl. Acad. Sci. USA 2021, 118. [Google Scholar] [CrossRef] [PubMed]
  43. Gonçalves, A.; Maisonnasse, P.; Donati, F.; Albert, M.; Behillil, S.; Contreras, V.; Naninck, T.; Marlin, R.; Solas, C.; Pizzorno, A.; et al. SARS-CoV-2 viral dynamics in non-human primates. PLoS Comput. Biol. 2021, 17, e1008785. [Google Scholar] [CrossRef] [PubMed]
  44. Kim, K.S.; Ejima, K.; Iwanami, S.; Fujita, Y.; Ohashi, H.; Koizumi, Y.; Asai, Y.; Nakaoka, S.; Watashi, K.; Aihara, K.; et al. A quantitative model used to compare within-host SARS-CoV-2, MERS-CoV, and SARS-CoV dynamics provides insights into the pathogenesis and treatment of SARS-CoV-2. PLoS Biol. 2021, 19, e3001128. [Google Scholar] [CrossRef] [PubMed]
  45. Dogra, P.; Ruiz-Ramírez, J.; Sinha, K.; Butner, J.D.; Peláez, M.J.; Rawat, M.; Yellepeddi, V.K.; Pasqualini, R.; Arap, W.; Sostman, H.D.; et al. Innate Immunity Plays a Key Role in Controlling Viral Load in COVID-19: Mechanistic Insights from a Whole-Body Infection Dynamics Model. ACS Pharmacol. Transl. Sci. 2021, 4, 248–265. [Google Scholar] [CrossRef]
  46. Perelson, A.S.; Ke, R. Mechanistic Modeling of SARS-CoV-2 and Other Infectious Diseases and the Effects of Therapeutics. Clin. Pharmacol. Ther. 2021, 109, 829–840. [Google Scholar] [CrossRef]
  47. Iwanami, S.; Ejima, K.; Kim, K.S.; Noshita, K.; Fujita, Y.; Miyazaki, Y.; Kohno, S.; Miyazaki, K.; Morimoto, S.; Nakaoka, S.; et al. Detection of significant antiviral drug effects on COVID-19 with reasonable sample sizes in randomized controlled trials: A modeling study combined with clinical data. PLoS Med. 2021, 18, e1003660. [Google Scholar] [CrossRef] [PubMed]
  48. Czuppon, P.; Débarre, F.; Gonçalves, A.; Tenaillon, O.; Perelson, A.S.; Guedj, J.; Blanquart, F. Success of prophylactic antiviral therapy for SARS-CoV-2: Predicted critical efficacies and impact of different drug-specific mechanisms of action. PLoS Comput. Biol. 2021, 17, e1008752. [Google Scholar] [CrossRef]
  49. Baccam, P.; Beauchemin, C.; Macken, C.A.; Hayden, F.G.; Perelson, A.S. Kinetics of influenza A virus infection in humans. J. Virol. 2006, 80, 7590–7599. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Maines, T.R.; Szretter, K.J.; Perrone, L.; Belser, J.A.; Bright, R.A.; Zeng, H.; Tumpey, T.M.; Katz, J.M. Pathogenesis of emerging avian influenza viruses in mammals and the host innate immune response. Immunol. Rev. 2008, 225, 68–84. [Google Scholar] [CrossRef]
  51. Sungnak, W.; Huang, N.; Bécavin, C.; Berg, M.; Queen, R.; Litvinukova, M.; Talavera-López, C.; Maatz, H.; Reichart, D.; Sampaziotis, F.; et al. SARS-CoV-2 entry factors are highly expressed in nasal epithelial cells together with innate immune genes. Nat. Med. 2020, 26, 681–687. [Google Scholar] [CrossRef] [Green Version]
  52. Hou, Y.J.; Okuda, K.; Edwards, C.E.; Martinez, D.R.; Asakura, T.; Dinnon, K.H.; Kato, T.; Lee, R.E.; Yount, B.L.; Mascenik, T.M.; et al. SARS-CoV-2 Reverse Genetics Reveals a Variable Infection Gradient in the Respiratory Tract. Cell 2020, 182, 429.e414–446.e414. [Google Scholar] [CrossRef]
  53. Efron, B.; Tibshirani, R.J. An Introduction to the Bootstrap; Chapman & Hall/CRC: New York, NY, USA, 1998. [Google Scholar]
  54. Bates, D.M.; Watts, D.G. Nonlinear Regression Analysis: Its Applications; Wiley: New York, NY, USA, 1988. [Google Scholar]
  55. Miao, H.; Xia, X.; Perelson, A.S.; Wu, H. On Identifiability of Nonlinear Ode Models and Applications in Viral Dynamics. SIAM Rev. Soc. Ind. Appl. Math. 2011, 53, 3–39. [Google Scholar] [CrossRef]
  56. Sender, R.; Bar-On, Y.M.; Gleizer, S.; Bernshtein, B.; Flamholz, A.; Phillips, R.; Milo, R. The total number and mass of SARS-CoV-2 virions. Proc. Natl. Acad. Sci. USA 2021, 118. [Google Scholar] [CrossRef]
  57. Ho, D.D.; Huang, Y. The HIV-1 vaccine race. Cell 2002, 110, 135–138. [Google Scholar] [CrossRef] [Green Version]
  58. Bonhoeffer, S.; May, R.M.; Shaw, G.M.; Nowak, M.A. Virus dynamics and drug therapy. Proc. Natl. Acad. Sci. USA 1997, 94, 6971–6976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Ribeiro, R.M.; Dixit, N.M.; Perelson, A.S. Modelling the in vivo growth rate of HIV: Implications for vaccination. In Multidisciplinary Approaches to Theory in Medicine; Paton, R., McNamara, L.A., Eds.; Elsevier: Amsterdam, The Netherland, 2006. [Google Scholar]
  60. Desmyter, J.; Melnick, J.L.; Rawls, W.E. Defectiveness of interferon production and of rubella virus interference in a line of African green monkey kidney cells (Vero). J. Virol. 1968, 2, 955–961. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Mosca, J.D.; Pitha, P.M. Transcriptional and posttranscriptional regulation of exogenous human beta interferon gene in simian cells defective in interferon synthesis. Mol. Cell. Biol. 1986, 6, 2279–2283. [Google Scholar] [CrossRef]
  62. Klimstra, W.B.; Tilston-Lunel, N.L.; Nambulli, S.; Boslett, J.; McMillen, C.M.; Gilliland, T.; Dunn, M.D.; Sun, C.; Wheeler, S.E.; Wells, A.; et al. SARS-CoV-2 growth, furin-cleavage-site adaptation and neutralization using serum from acutely infected hospitalized COVID-19 patients. J. Gen. Virol. 2020, 101, 1156–1169. [Google Scholar] [CrossRef] [PubMed]
  63. Karimzadeh, S.; Bhopal, R.; Nguyen Tien, H. Review of infective dose, routes of transmission and outcome of COVID-19 caused by the SARS-COV-2: Comparison with other respiratory viruses. Epidemiol. Infect. 2021, 149, e96. [Google Scholar] [CrossRef]
  64. Monchatre-Leroy, E.; Lesellier, S.; Wasniewski, M.; Picard-Meyer, E.; Richomme, C.; Boué, F.; Lacôte, S.; Murri, S.; Pulido, C.; Vulin, J.; et al. Hamster and ferret experimental infection with intranasal low dose of a single strain of SARS-CoV-2. J. Gen. Virol. 2021, 102. [Google Scholar] [CrossRef]
Figure 1. Model 1 fitting to the data. Prediction of the viral dynamics model 1 (solid line) along with the data (filled circle) containing the total viral load (RNA copies per mL) and the infectious virus titer (PFU per mL), collected from ferrets infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2.
Figure 1. Model 1 fitting to the data. Prediction of the viral dynamics model 1 (solid line) along with the data (filled circle) containing the total viral load (RNA copies per mL) and the infectious virus titer (PFU per mL), collected from ferrets infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2.
Viruses 13 01635 g001
Figure 2. Model 2 fitting to the data. Prediction of the viral dynamics model 2 (solid line) along with the data (filled circle) containing the total viral load (RNA copies per mL) and the infectious virus titer (PFU per mL), collected from ferrets infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2.
Figure 2. Model 2 fitting to the data. Prediction of the viral dynamics model 2 (solid line) along with the data (filled circle) containing the total viral load (RNA copies per mL) and the infectious virus titer (PFU per mL), collected from ferrets infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2.
Viruses 13 01635 g002
Figure 3. Model 3 fitting to the data. Prediction of the viral dynamics model 3 (solid line) along with the data (filled circle) containing the total viral load (RNA copies per mL) and the infectious virus titer (PFU per mL), collected from ferrets infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2.
Figure 3. Model 3 fitting to the data. Prediction of the viral dynamics model 3 (solid line) along with the data (filled circle) containing the total viral load (RNA copies per mL) and the infectious virus titer (PFU per mL), collected from ferrets infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2.
Viruses 13 01635 g003aViruses 13 01635 g003b
Figure 4. Sensitivity analysis. Partial rank correlation coefficients (PRCC) from the Latin hypercube sampling method of the total viral-peak size (1st column), infectious viral peak size (2nd column), time to viral-peak (3rd column), and cell loss (4th column) for Model 1 (1st row), Model 2 (2nd Row), and Model 3 (3rd row).
Figure 4. Sensitivity analysis. Partial rank correlation coefficients (PRCC) from the Latin hypercube sampling method of the total viral-peak size (1st column), infectious viral peak size (2nd column), time to viral-peak (3rd column), and cell loss (4th column) for Model 1 (1st row), Model 2 (2nd Row), and Model 3 (3rd row).
Viruses 13 01635 g004aViruses 13 01635 g004b
Figure 5. Role of Initial infectious viruses. Fitting data with initial virus at the infection site ( V i 0 ) estimated vs fixed. Shown here are two representative animals from each of the two groups infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2. The best-fit curves corresponding to V i 0 = 9.25 PFU/mL (F13-E-2) and V i 0 = 9.12 PFU/mL (CTan-H-2) were obtained by allowing V i 0 to be estimated by data fitting process, while other best-fit curves were obtained by fixing the values of V i 0 at 10 2 , 10 3 , and 10 4 PFU/mL.
Figure 5. Role of Initial infectious viruses. Fitting data with initial virus at the infection site ( V i 0 ) estimated vs fixed. Shown here are two representative animals from each of the two groups infected with F13-E SARS-CoV-2 and CTan-H SARS-CoV-2. The best-fit curves corresponding to V i 0 = 9.25 PFU/mL (F13-E-2) and V i 0 = 9.12 PFU/mL (CTan-H-2) were obtained by allowing V i 0 to be estimated by data fitting process, while other best-fit curves were obtained by fixing the values of V i 0 at 10 2 , 10 3 , and 10 4 PFU/mL.
Viruses 13 01635 g005aViruses 13 01635 g005b
Table 1. Parameters of Model 1 estimated from the data fitting. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Table 1. Parameters of Model 1 estimated from the data fitting. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Ferret V i 0
(PFU/mL)
β
(PFU/mL)
δ
(/Day)
α p
(/Cell/Day)
c
(/Day)
F13-E SARS-CoV-2 virus
F13-E-1 2.64 3.55 × 10 5 4.05 2.33 × 10 4 7934 7.13
0.0002 ,   51.75 1.16 × 10 5 ,   9.99 × 10 5 1.05 ,   12.79 4.42 × 10 5 ,   0.92 753 ,   75,412 1.93 ,   49.91
F13-E-2 9.25 4.48 × 10 5 5.05 1.63 × 10 4 8761 6.30
0.12 ,   116.55 1.66 × 10 5 ,   9.98 × 10 5 1.02 ,   12.66 5.48 × 10 5 ,   0.93 653 ,   61,520 2.52 ,   49.98
F13-E-3 22.39 7.05 × 10 5 5.99 1.90 × 10 4 5102 6.64
0.0016 ,   274.25 2.34 × 10 5 ,   1.00 × 10 4 1.11 ,   11.56 4.84 × 10 5 ,   0.91 765 ,   54,913 3.28 ,   49.98
G. Ave. (F13-E) 8.18 4.82 × 10 5 4.97 1.93 × 10 4 7078 6.68
CTan-H SARS-CoV-2 virus
CTan-H-1 0.11 6.50 × 10 5 4.89 1.79 × 10 4 4152 3.67
0.002 ,   80.79 1.38 × 10 5 ,   1.00 × 10 4 0.91 ,   19.99 2.38 × 10 5 ,   0.89 632 ,   94,492 2.22 ,   49.96
CTan-H-2 9.12 9.90 × 10 5 5.52 3.06 × 10 4 1461 3.50
0.05 ,   231.2 2.26 × 10 5 ,   1.00 × 10 4 1.12 ,   14.59 8.57 × 10 5 ,   0.93 618 ,   32,596 3.04 ,   49.89
CTan-H-3 14.13 9.99 × 10 5 5.99 7.32 × 10 5 16,641 13.41
0.48 ,   736.69 3.37 × 10 5 ,   1.00 × 10 4 1.04 ,   19.98 1.63 × 10 5 ,   0.92 723 ,   125,883 2.75 ,   49.99
G. Ave. (CTan-H) 2.46 8.63 × 10 5 5.45 1.59 × 10 4 4656 5.57
Overall G. Ave. 4.48 6.46 × 10 5 5.20 1.75 × 10 4 5741 6.10
p-value≥0.05≥0.05≥0.05≥0.05≥0.05≥0.05
Table 2. Parameters of Model 2 estimated from the data fitting. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Table 2. Parameters of Model 2 estimated from the data fitting. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Ferret V i 0
(PFU/mL)
β
(PFU/mL)
δ
(/Day)
α p
(/Cell/Day)
c
(/Day)
k
(/Day)
F13-E SARS-CoV-2 virus
F13-E-1 3.94 5.21 × 10 5 3.25 2.07 × 10 4 7018 7.99 11.96
1.36 × 10 5 ,   32.58 1.83 × 10 5 ,   2.67 × 10 4 1.78 ,   8.87 4.53 × 10 5 ,   1.05 × 10 3 4980 ,   22,125 2.62 ,   18.21 2.83 ,   19.20
F13-E-2 11.09 8.28 × 10 5 4.03 1.65 × 10 4 3512 3.74 8.71
2.54 × 10 6   ,     71.45 3.07 × 10 5 ,   2.79 × 10 4 1.77 ,   9.92 5.44 × 10 5 ,   6.19 × 10 4 1353 ,   15,579 2.12 ,   13.94 2.68 ,   19.98
F13-E-3 26.80 1.40 × 10 4 3.67 1.89 × 10 4 1545 4.28 10.99
1.09 × 10 3 ,   99.58 5.34 × 10 5 ,   5.00 × 10 4 1.47 ,   9.95 4.55 × 10 5 ,   7.79 × 10 4 390 ,   8088 1.83 ,   21.05 2.26 ,   19.99
G. Ave. (F13-E) 10.54 8.45 × 10 5 3.64 1.86 × 10 4 3365 5.04 10.46
CTan-H SARS-CoV-2 virus
CTan-H-1 0.24 1.06 × 10 4 4.98 1.44 × 10 4 5091 3.65 3.93
3.26 × 10 4 ,   10.26 1.98 × 10 5 ,   4.98 × 10 4 2.13 ,   9.98 3.17 × 10 5 ,   1.38 × 10 3 1946 ,   39,011 2.38 ,   14.84 1.72 ,   19.52
CTan-H-2 21.76 8.02 × 10 5 6.42 3.50 × 10 4 2143 3.98 12.52
3.48 × 10 3 ,   99.97 1.63 × 10 6 ,   4.87 × 10 4 1.77 ,   9.99 1.01 × 10 4 ,   7.49 × 10 3 273 ,   6474 1.28 ,   16.24 2.58 ,   19.99
CTan-H-3 21.03 3.19 × 10 4 4.41 5.07 × 10 5 4971 9.38 4.25
0.03 ,   95.20 1.26 × 10 4 ,   5.00 × 10 4 1.54 ,   9.89 1.82 × 10 5 ,   1.97 × 10 4 4962 ,   5376 2.51 ,   30.18 1.61 ,   19.99
G. Ave. (CTan-H) 4.80 1.49 × 10 4 5.20 1.37 × 10 4 3785 5.15 5.93
Overall G. Ave. 7.11 1.12 × 10 4 4.35 1.60 × 10 4 3569 5.09 7.88
p-value 0.05 0.05 0.05 0.05 0.05 0.05 0.05
Table 3. Parameters of Model 3 estimated from the data fitting. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Table 3. Parameters of Model 3 estimated from the data fitting. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Ferret V i 0
(PFU/mL)
β
(PFU/mL)
δ
(/Day)
α p
(/Cell/Day)
c
(/Day)
σ
(/IFN/Day)
F13-E SARS-CoV-2 virus
F13-E-1 2.53 4.74 × 10 5 2.42 2.20 × 10 4 3980 5.43 2.99 × 10 7
6.13 × 10 4 ,   35.91 1.07 × 10 5 ,   1.88 × 10 4 1.66 ,   8.09 5.63 × 10 5 ,   1.39 × 10 3 3945 ,   12,712 2.51 ,   20.02 4.73 × 10 12 ,   2.54 × 10 6
F13-E-2 14.35 4.32 × 10 5 2.22 1.62 × 10 4 3984 4.12 1.05 × 10 6
8.38 × 10 3 ,     92.33 1.51 × 10 5 ,   2.11 × 10 4 1.70 ,   4.86 4.77 × 10 5 ,   7.36 × 10 4 2987 ,   13,169 2.12 ,   8.61 3.21 × 10 12 ,   5.72 × 10 6
F13-E-3 95.79 3.79 × 10 5 2.26 1.97 × 10 4 3951 6.25 1.12 × 10 6
14.37 ,   96.68 9.60 × 10 6 ,   1.36 × 10 4 1.18 ,   4.72 6.58 × 10 5 ,   9.20 × 10 4 3942 ,   4063 2.23 ,   29.06 8.35 × 10 8 ,   8.29 × 10 6
G. Ave. (F13-E) 15.16 4.27 × 10 5 2.30 1.92 × 10 4 3972 5.19 7.06 × 10 7
CTan-H SARS-CoV-2 virus
CTan-H-1 1.72 1.55 × 10 5 4.62 1.76 × 10 4 14,140 4.07 3.38 × 10 6
0.01 ,   11.12 5.55 × 10 6 ,   9.04 × 10 5 2.06 ,   9.96 3.04 × 10 5 ,   6.92 × 10 4 13,996 ,   16,093 2.68 ,   11.78 1.18 × 10 7 ,   9.98 × 10 6
CTan-H-2 16.22 1.90 × 10 5 1.86 3.18 × 10 4 3983 3.72 6.72 × 10 6
4.83 × 10 3 ,   97.54 6.85 × 10 6 ,   1.89 × 10 4   1.48 ,   5.96 8.02 × 10 5 ,   1.22 × 10 3 2955 ,   4070 2.15 ,   23.20 6.40 × 10 11 ,   1.00 × 10 5
CTan-H-3 35.21 2.37 × 10 5 2.96 6.03 × 10 5 14,100 3.48 7.30 × 10 6
0.23 ,   99.99 9.95 × 10 6 ,   2.81 × 10 4 1.41 ,   9.92 1.52 × 10 5 ,     2.95 × 10 4 14,066 ,   14,121 2.00 ,   30.67 2.07 × 10 7 ,   1.00 × 10 5
G. Ave. (CTan-H) 9.95 1.91 × 10 5 2.94 1.50 × 10 4 9270 3.75 5.50 × 10 6
Overall G. Ave. 12.28 2.85 × 10 5 2.60 1.70 × 10 4 6068 4.41 1.97 × 10 6
p-value 0.05 0.05 0.05 0.05 0.05 0.05 0.05
Table 4. Basic reproduction number ( R 0 ), pre-peak infectious viral growth rate ( r g ), post-peak infectious viral decay rate ( r d ), pre-peak infectious viral doubling time ( t 2 ), and post-peak infectious virus half-life ( t 1 / 2 ) for each ferret. The growth rate and decay rates are computed using the eigenvalues of the Jacobian: r   =   c + δ   +   c + δ 2   +   4 c δ R t 1 / 2 ,   where R t is approximation to the pre-peak and post-peak reproduction numbers. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Table 4. Basic reproduction number ( R 0 ), pre-peak infectious viral growth rate ( r g ), post-peak infectious viral decay rate ( r d ), pre-peak infectious viral doubling time ( t 2 ), and post-peak infectious virus half-life ( t 1 / 2 ) for each ferret. The growth rate and decay rates are computed using the eigenvalues of the Jacobian: r   =   c + δ   +   c + δ 2   +   4 c δ R t 1 / 2 ,   where R t is approximation to the pre-peak and post-peak reproduction numbers. G. Ave. stands for geometric average. p-values indicate the result of comparison between F13-E and CTan-H groups.
Ferret R 0 r g r d t 2 t 1 / 2 Cell Loss (%)
F13-E SARS-CoV-2 virus
F13-E-12.272.662.176.267.6885.76
F13-E-22.012.342.047.108.1479.98
F13-E-31.721.951.778.529.3869.95
G. Ave. (F13)1.992.301.997.248.3778.28
CTan-H SARS-CoV-2 virus
CTan-H-12.692.702.166.167.6991.48
CTan-H-2 2.322.261.857.378.9986.48
CTan-H-31.521.941.748.559.5459.35
G. Ave. (CTan-H)2.122.281.917.298.7177.72
Overall G. Ave.2.052.291.957.278.5478.00
p-value 0.05 0.05 0.05 0.05 0.05 0.05
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vaidya, N.K.; Bloomquist, A.; Perelson, A.S. Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets. Viruses 2021, 13, 1635. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081635

AMA Style

Vaidya NK, Bloomquist A, Perelson AS. Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets. Viruses. 2021; 13(8):1635. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081635

Chicago/Turabian Style

Vaidya, Naveen K., Angelica Bloomquist, and Alan S. Perelson. 2021. "Modeling Within-Host Dynamics of SARS-CoV-2 Infection: A Case Study in Ferrets" Viruses 13, no. 8: 1635. https://0-doi-org.brum.beds.ac.uk/10.3390/v13081635

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop