Dynamically linking influenza virus infection kinetics, lung injury, inflammation, and disease severity

Abstract
Introduction
Results
Discussion
Materials and methods
Appendix 1
Appendix 2
Appendix 3
Data availability
References
Article and author information
Metrics

Abstract

Influenza viruses cause a significant amount of morbidity and mortality. Understanding host immune control efficacy and how different factors influence lung injury and disease severity are critical. We established and validated dynamical connections between viral loads, infected cells, CD8⁺ T cells, lung injury, inflammation, and disease severity using an integrative mathematical model-experiment exchange. Our results showed that the dynamics of inflammation and virus-inflicted lung injury are distinct and nonlinearly related to disease severity, and that these two pathologic measurements can be independently predicted using the model-derived infected cell dynamics. Our findings further indicated that the relative CD8⁺ T cell dynamics paralleled the percent of the lung that had resolved with the rate of CD8⁺ T cell-mediated clearance rapidly accelerating by over 48,000 times in 2 days. This complimented our analyses showing a negative correlation between the efficacy of innate and adaptive immune-mediated infected cell clearance, and that infection duration was driven by CD8⁺ T cell magnitude rather than efficacy and could be significantly prolonged if the ratio of CD8⁺ T cells to infected cells was sufficiently low. These links between important pathogen kinetics and host pathology enhance our ability to forecast disease progression, potential complications, and therapeutic efficacy.

Introduction

Over 15 million respiratory infections and 200,000 hospitalizations result from influenza A viruses (IAVs) each year (Thompson et al., 2004; Simonsen et al., 2000; Taubenberger and Morens, 2008; Medina and García-Sastre, 2011). The incidence and severity of IAV infections increases when new strains emerge and/or when there is a lack of prior immunity. A robust immune response is crucial for resolving viral infections, but immune-mediated pathology can exacerbate disease (Duan and Thomas, 2016; Moskophidis and Kioussis, 1998; Mauad et al., 2010; Rygiel et al., 2009; La Gruta et al., 2007). High viral loads can also play a role in disease progression (Boon et al., 2011; de Jong et al., 2006), but these do not always correlate with the strength of the host response or with disease severity (Granados et al., 2017; Marathe et al., 2016; Toapanta and Ross, 2009; Smith et al., 2019; Gao et al., 2013). An understanding of how viral loads, host immune responses, and disease progression are related is critical to identify disease-specific markers that may help predict hospitalization or other complications.

During IAV infection in both humans and animals, viral loads increase rapidly for the first 1–2 days of infection before reaching a peak (e.g., as in Smith et al., 2018; Baccam et al., 2006; Miao et al., 2010; Toapanta and Ross, 2009; Carrat et al., 2008; Srivastava et al., 2009; Gao et al., 2013). If the host is infected with a novel strain or has no prior immunity, viral loads in the lung then begin to decline, first slowly (sometimes <1 log₁₀ TCID₅₀/d ) then rapidly (> 4-5 log₁₀ TCID₅₀/d) (Smith et al., 2018; Gao et al., 2013). We previously quantified this biphasic viral decline with a mathematical model, which indicated that the rate of infected cell clearance increases as the density of infected cells decreases (Smith et al., 2018). The timing of the second, rapid viral decay phase coincides with the expansion of CD8⁺ T cells, which are the primary cell responsible for clearing infected cells and resolving the infection (Zhang and Bevan, 2011; Chen et al., 2018; McMichael et al., 1983; La Gruta and Turner, 2014; Kreijtz et al., 2011), and, to a lesser extent, neutralizing antibodies (Chen et al., 2018; Kreijtz et al., 2011; Fang and Sigal, 2005; Wang et al., 2015; McMichael et al., 1983) and cytotoxic CD4⁺ T cells (Wilkinson et al., 2012). For the CD8⁺ T cell response, in particular, it remains unclear whether the efficacy of these cells is dictated by their own density (Graw and Regoes, 2009; Wiedemann et al., 2006), infected cell density (Merrill, 1982; Caramalho et al., 2009; Halle et al., 2016), or both (Gadhamsetty et al., 2014). While quantifying dynamically changing CD8⁺ T cell efficacy is difficult in vitro and in vivo, the question is ripe for in silico investigation. Several modeling studies have described CD8⁺ T cell-mediated infected cell clearance for various viral infections in humans and animals, including IAV, HIV, and LCMV (e.g., as in Ganusov et al., 2011; De Boer and Perelson, 1998; Perelson and Bell, 1982; Miao et al., 2010; Cao et al., 2016; Lee et al., 2009; Price et al., 2015; Graw and Regoes, 2009; Merrill, 1982; Gadhamsetty et al., 2014; Baral et al., 2019; Bonhoeffer et al., 2000; Conway and Perelson, 2015). Some of these studies have also attempted to link this response and other immune responses to inflammation or disease severity (Price et al., 2015; Baral et al., 2019; Manchanda et al., 2014), but have not yet found the appropriate mathematical relation with the available data. In addition, for IAV infections in particular, varied efficiency of the CD8⁺ T cell response throughout the course of infection, their early antigen-specific and lung-resident responses (Yoon et al., 2007; McGill and Legge, 2009; Keating et al., 2018), and the balanced consequences on viral loads and host pathology have not yet been investigated in detail.

A better understanding of infected cell clearance may also yield insight into the damage induced to the lung and the ensuing immunopathology during IAV infection. In general, widespread alveolar disease is observed in patients who succumb to the infection (Bautista et al., 2010; Shieh et al., 2010; Gao et al., 2013), and CT scans show bilateral and multi-focal areas of lung damage (e.g. as in Kohr et al., 2010; Li et al., 2011; Ajlan et al., 2009; Soto-Abraham et al., 2009; Gill et al., 2010; Li et al., 2012). Further, hospitalized patients that died as a result from infection with the 2009 H1N1 influenza virus had large numbers of CD8⁺ T cells present in their lung tissue (Mauad et al., 2010). Large pulmonary populations of CD8⁺ T cells contribute to lung injury by targeting IAV-infected cells (reviewed in Duan and Thomas, 2016) and inducing ‘bystander damage’ to uninfected epithelial cells (van de Sandt et al., 2017). However, their population levels do not necessarily indicate active infected cell removal or immunopathology as their extended presence in the lung also contributes to surveillance and repair (Matheu et al., 2013; Sun et al., 2009). Macrophages and neutrophils also persist in the lung following viral clearance, and their role in inflammation and tissue damage is well documented (reviewed in Watanabe et al., 2019; La Gruta et al., 2007). A favorable outcome requires a balance between immune-mediated protection and pathology, and the viral-immunological landscape that drives disease and the markers that distinguish mild from severe influenza are complex. This is particularly true in humans (Tang et al., 2019) where obtaining high quality data from the lung is challenging, and viral loads in the upper respiratory tract do not always reflect the lower respiratory tract environment (Feikin et al., 2017; Ai et al., 2020; Wong et al., 2020; Li et al., 2012; Ruuskanen et al., 2011).

The accumulation of damage to the epithelium during IAV infection, either from virally-induced cell lysis or immune-mediated effects, is relatively understudied. We and others have modeled the lung damage and inflammation inflicted during pulmonary infections (e.g., as in Smith et al., 2011b; Reynolds et al., 2006; Dunster et al., 2014; Ramirez-Zuniga et al., 2019; Baral et al., 2019) but did so without sufficient data to constrain the model formulations. The difficulty in measuring the dynamics of infected cells and in establishing how damage correlates to specific host responses has been the primary impediment. However, recent technological advances, including the use of reporter viruses (Manicassamy et al., 2010; Tran et al., 2013; Karlsson et al., 2015; Luker and Luker, 2010) and lung histomorphometry (Marathe et al., 2016; Marathe et al., 2017; Sartorius et al., 2007; Zachariadis et al., 2006; Boyd et al., 2020) within animal models, have provided an opportunity to acquire these types of measurements. Whole lung histomorphometry, which is broadly defined as a quantitative microscopic measurement of tissue shape, has recently increased in use due to the ability to directly stain, visualize, and quantify areas of infected cells. Even with these techniques, quantitative data over the course of infection is not currently available. Having data such as these should help reconcile potential nonlinearities in infected cell clearance and provide insight into the accumulated lung inflammation and damage, which are generally thought to be markers of disease severity. In addition, it should bolster the development of robust, predictive computational models, which have historically lacked constraint to these types of data.

In general, disease severity may not be directly correlated to viral loads or specific immunological components. In humans infected with IAV, viral loads typically increase prior to the onset of systemic symptoms, which peak around 2–3 d post-infection (pi) (Carrat et al., 2008; Lee and Lee, 2010; Han et al., 2018). Some symptoms (e.g., fever) are cytokine mediated (Monto et al., 2000), but respiratory symptoms often last longer and can remain following viral clearance (Lee and Lee, 2010; Eccles, 2005). This is particularly true when there is scarring of the lung tissue (Li et al., 2012). In the murine model, weight loss is used as an indicator of disease progression and severity, where greater weight loss corresponds to more severe disease (Trammell and Toth, 2011; Bouvier and Lowen, 2010; Parzych et al., 2013). Animal weights generally drop slowly in the first several days during an IAV infection and more rapidly as the infection begins to resolve (Srivastava et al., 2009; Huang et al., 2017). This is unlike viral load dynamics in these animals, which increase rapidly during the first 0–3 d pi then remain relatively constant prior to resolution (Smith et al., 2018). Because weight loss can occur following resolution of the infection, immune-mediated pathology is thought to play a role (Lauder et al., 2011; Xu et al., 2004; Ostler et al., 2002; Parzych et al., 2013; Duan and Thomas, 2016). Host and pathogen factors, such as age, viral proteins, and inoculum size, have also been shown to influence disease progression (Toapanta and Ross, 2009; Lu et al., 2018; Smith et al., 2019; McAuley et al., 2007). While the precise contribution of different factors to IAV-associated disease and mortality remain elusive, having tools that can link these with disease severity and decipher correlation from causation would improve our ability to effectively predict, understand, and treat the disease.

To gain deeper insight into the dynamics of viral resolution during primary pulmonary influenza infection and investigate the connection between viral loads, damage, inflammation, and disease severity, we developed and validated a mathematical model that describes viral kinetics and CD8⁺ T cell-mediated infected cell clearance. The model accurately predicted the dynamics of effector and memory CD8⁺ phenotypes, agreed with our previous results that infected cells are cleared in a density-dependent manner with CD8 efficiency rapidly accelerating by over 48,000 times in 2 d, and illuminated tradeoffs between the innate and adaptive immune responses. Our model predicted that infection duration is dependent on the magnitude of CD8⁺ T cells rather than their efficacy, which we verified by depleting CD8⁺ T cells. Quantitative whole lung histomorphometry showed that infected areas increase during the first 6 d of infection before clearing rapidly, that the relative number of CD8⁺ T cells corresponded to the cleared area, and corroborated the model-predicted infected cell dynamics. These data also showed that the infection-induced lung injury is distinct from inflammation, and that each correlates to different immune responses and could be predicted using independent mathematical equations. In addition, the data revealed nonlinear connections between these two pathological measurements and disease severity. Together, our data, model, and analyses provide a robust quantification of the density-dependent nature of CD8⁺ T cell-mediated clearance, and the critical connections between these cells and the dynamics of viral loads, infected cells, lung injury, inflammation, and disease severity.

Results

Virus and CD8⁺ T cell kinetics

In animals infected with 75 TCID₅₀ PR8, virus rapidly infects cells or is cleared within 4 hr pi and is undetectable in all animals at this time (Figure 1). Virus then increases exponentially and peaks after ∼2 d pi. Following the peak, virus enters a biphasic decline. In the first phase (2–6 d pi), virus decays slowly and at a relatively constant rate (0.2 log₁₀ TCID₅₀/d) (Smith et al., 2018). In the second phase (7–8 d pi), virus is cleared rapidly with a loss of 4–5 log₁₀ TCID₅₀ in 1–2 d (average of –3.8 log₁₀ TCID₅₀/d) (Smith et al., 2018). CD8⁺ T cells remain at their baseline level from 0–2 d pi, where ∼15% are located in the circulating blood, ∼75% in the lung vasculature, and ∼10% in the lung parenchyma (Figure 1—figure supplement 1) with the majority in the lung as CD43^–CD127⁺. IAV-specific CD8⁺ T cells that are recently primed (CD25⁺CD43⁺) or effector cells (CD25^–CD43⁺) (Keating et al., 2018) begin appearing in the lung and increase slightly from 2–3 d pi but remain at low levels. These cells are thought to proliferate within the lung at least once by 4 d pi (McGill and Legge, 2009), and their populations briefly contract (3–5 d pi) before expanding rapidly (5–8 d pi). During the expansion phase, >95% of recovered cells are in the lung (∼70% in the parenchyma and ∼30% in the vasculature) with CD25⁺CD43⁺ and CD25^–CD43⁺ as the predominant phenotypes. This is in accordance with prior studies that showed these phenotypes are IAV-specific and that their expansion dynamics in the lung were not altered by removing blood-borne CD8⁺ T cells from the analysis (Anderson et al., 2012; Keating et al., 2018).

Figure 1 with 2 supplements see all

Download asset Open asset

Schematic and fit of the CD8⁺ T cell viral kinetic model.

(A) Schematic of the CD8⁺ T cell model in Equation (1)-(6). In the model, target cells ( $T$ ) are infected at rate $β V$ . Infected cells ( $I_{1}$ ) undergo an eclipse phase and transition to become productively-infected cells ( $I_{2}$ ) at rate $k$ . Virus ( $V$ ) is produced by infected cells at rate $p$ and is cleared at rate $c$ . Infected cells are cleared at rate $δ$ by non-specific mechanisms and at rate $δ_{E} (I_{2}, E)$ by effector CD8⁺ T cells ( $E$ ; denoted CD8_E). The dashed lines represent interactions between infected cells and CD8_E. Initial CD8_E influx ( $ξ (E) = ξ / (K_{E} + E)$ ) is proportional to infected cells and is limited by CD8_E with half-saturation constant $K_{E}$ . CD8_E expansion ( $η$ ) occurs proportional to infected cells $τ_{E}$ days ago. Memory CD8⁺ T cell ( $E_{M}$ ; denoted CD8_M) generation occurs at rate $ζ$ and proportional to CD8_E $τ_{M}$ days ago. (B) Fit of the CD8⁺ T cell model (Equation (1)-(6)) to virus and total CD8⁺ T cells from the lungs of mice infected with 75 TCID₅₀ PR8 (10 mice per time point). The total number of CD8⁺ T cells is $\hat{E} = E + E_{M} + {\hat{E}}_{0}$ . (C) Total CD8⁺ T cells in the lung parenchyma (gray circles) and overlay of the model predicted values ( $E + E_{M}$ ). (D) Fit of the model to virus, CD25^–CD43⁺ CD8⁺ T cells (cyan diamonds; $E$ ), and CD43^–CD127⁺CD103⁺CD69⁺ CD8⁺ T cells (blue diamonds; $E_{M}$ ) (five mice per time point). The solid and dashed black lines are the optimal solutions and the gray shading is are the model solutions using parameter sets within the 95% CIs. Parameters are given in Table 1. Data are shown as mean ± standard deviation.

CD8⁺ T cell expansion corresponds to the second viral decay phase with sixty percent of mice clearing the infection by 8 d pi and the other forty percent by 9 d pi (Figure 1). Most CD43⁺CD8⁺ T cell phenotypes decline following viral clearance (8–10 d pi) but do not return to their baseline level by 12 d pi. Long-lived antigen-specific memory phenotypes down regulate CD43 (Jones et al., 1994; Harrington et al., 2000; Onami et al., 2002) and gradually increase substantially beginning at 9 d pi with most as CD25^–, which is qualitatively similar to other studies (Harrington et al., 2000). At 12 d pi, ∼55% of CD8⁺ T cells remain in the lung parenchyma and ∼20% in the circulating blood indicating exit from the lung (Figure 1—figure supplement 1).

Viral kinetic model with density-dependent CD8⁺ T cell-mediated clearance

We previously described the influenza viral load kinetics and biphasic decline during primary pulmonary infection using a density-dependent (DD) model (Smith et al., 2018), which assumes that the rate of infected cell clearance increases as the density of infected cells decreases (i.e., $δ_{d} (I_{2}) = δ_{d} / (K_{δ} + I_{2})$ ). Because the rapid decay of virus is thought to be due to the clearance of infected cells by CD8⁺ T cells, it is unknown if early CD8⁺ T cell presence contributes to infected cell clearance, and no model has captured the entire CD8⁺ T cell time course, we developed a model that describes the dynamics of these cells and their efficiency in resolving the infection (Equation (1)-(6); Figure 1A). The model includes equations for effector ( $E$ , denoted CD8_E) and memory ( $E_{M}$ , denoted CD8_M) CD8⁺ T cells, and two mechanisms of infected cell clearance. The first mechanism is from unspecified, innate mechanisms, which is relatively constant ( $δ$ ) and primarily acts during the first viral decay phase (2–6 d pi). The second is the CD8_E-mediated infected cell clearance, which occurs at a rate that increases as the density of infected cells decreases ( $δ_{E} (I_{2}, E) = δ_{E} E / (K_{δ_{E}} + I_{2})$ ) and primarily acts during the rapid, second viral decay phase (7–8 d pi). Excluding this density dependence entirely resulted in a significant and premature decline in viral loads, which disagreed with the experimental data. We also tested whether the density dependence could be included in the CD8⁺ T cell expansion rate rather than in the infected cell clearance rate (see Appendix 1) as other models have done (e.g., as in Baral et al., 2019; Bonhoeffer et al., 2000; Conway and Perelson, 2015). This modification yielded a close fit to the CD8⁺ T cell data at 6–10 d pi but not at the early time points (Appendix 1—figure 1). In addition, the viral load data was underestimated at 7 d pi causing the solution to miss the rapid decline between 7–8 d pi and result in a statistically poorer fit. Thus, retaining the density-dependence in the rate of infected cell clearance most accurately captured the entire dataset. The model includes terms for the initial CD8_E influx at 2–3 d pi ( $ξ I_{2} / (K_{E} + E)$ ) and for CD8_E expansion ( $η E I_{2} (t - τ_{E})$ ), which accounts for the larger increase between 5–8 d pi. To capture the contraction of CD8⁺ T cells between these times (3–5 d pi), it was necessary to assume that the initial CD8_E influx is regulated by their own population (i.e., $ξ (E) = ξ / (K_{E} + E)$ ). In both terms, the increase is proportional to the number of infected cells. Although memory CD8⁺ T cells were not the primary focus here, it was necessary to include the CD8_M population because CD8⁺ T cells are at a significantly higher level at 10–12 d pi than at 0 d pi (Figure 1B). Fitting the model simultaneously to viral loads and total CD8⁺ T cells from the (non-perfused) lungs of infected animals illustrated the accuracy of the model (Figure 1B). The resulting parameter values, 95% confidence intervals (CIs), ensembles, and histograms are given in Table 1, Figure 2, and Figure 2—figure supplements 1–2.

Figure 2 with 2 supplements see all

Download asset Open asset

Parameter ensembles and histograms.

Parameter ensembles (A) and histograms (B) resulting from fitting the CD8⁺ T cell viral kinetic model (Equation (1)-(6)) to viral titers and total CD8⁺ T cells from mice infected with 75 TCID₅₀ PR8. (A) The rates of infected cell clearance by non-specific mechanisms (δ) and by CD8_E ( $δ_{E}$ ) are slightly negatively correlated. Correlations were also present between the rates of CD8_E clearance ( $d_{E}$ ), CD8_E expansion (η), and CD8_M generation (ζ). The axes limits reflect imposed bounds. Additional ensemble plots are in Figure 2—figure supplements 1–2. (B) The histograms show that the majority of parameters are well-defined with the exception of the eclipse phase transition rate ( $k$ ), one of the half-saturation constants ( $K_{E}$ ), and the CD8_M generation delay ( $τ_{M}$ ).

Table 1

CD8⁺ T cell model parameters.

Parameters and 95% confidence intervals obtained from fitting the CD8⁺ T cell model (Equation (1)-(6)) to viral titers and total CD8⁺ T cells (‘Total CD8’) or viral titers, CD25^–CD43⁺CD8⁺ T cells, and CD43^–CD127⁺CD103⁺CD69⁺CD8⁺ T cells (‘Specific CD8_E,M Phenotypes’) from mice infected with 75 TCID₅₀ PR8. CD8_E and CD8_M denote effector ( $E$ ) and memory ( $E_{M}$ ) CD8⁺ T cells, respectively. The total number of CD8⁺ T cells is $\hat{E} = E + E_{M} + \hat{E_{0}}$ and is denoted by CD8.

Parameter	Description	Units	Total CD8		Specific CD8_E,Mphenotypes
Parameter	Description	Units	Value	95% CI	Value	95% CI
$β$	Virus infectivity	${TCID}_{50}^{- 1} d^{- 1}$	6.2 × 10^–5	[5.3 × 10^–6, 1.0 × 10^–4]	3.7 × 10^–5	[1.1 × 10^–5, 9.4 × 10^–5]
$k$	Eclipse phase transition	$d^{- 1}$	4.0	[4.0, 6.0]	5.1	[4.0, 6.0]
$p$	Virus production	${TCID}_{50} {cell}^{- 1} d^{- 1}$	1.0	[5.8 × 10^–1 × 1.1 × 10²]	1.5	[0.73, 13.6]
$c$	Virus clearance	$d^{- 1}$	9.4	[5.6, 9.5 × 10²]	12.1	[5.8, 17.5]
$δ$	Infected cell clearance	$d^{- 1}$	2.4 × 10^–1	[1.0 × 10^–1, 6.6 × 10^–1]	3.0 × 10^–1	[1.9 × 10^–1, 5.9 × 10^–1]
$δ_{E}$	Infected cell clearance by CD8_E	$cells {CD8}_{E}^{- 1} d^{- 1}$	1.9	[3.3. × 10^–1, 2.0]	5.7	[1.7, 8.5]
$K_{δ_{E}}$	Half-saturation constant	$cells$	4.3 × 10²	[1.0 × 10², 2.9 × 10⁵]	1.3 × 10²	[1.0 × 10¹, 8.6 × 10²]
$ξ$	CD8_E infiltration	$C D 8_{E}^{2} {c e l l}^{- 1} d^{- 1}$	2.6 × 10⁴	[1.3 × 10², 8.7 × 10⁴]	2.9 × 10³	[8.0 × 10², 1.7 × 10⁴]
$K_{E}$	Half-saturation constant	${CD8}_{E}$	8.1 × 10⁵	[1.0 × 10³, 1.0 × 10⁶]	2.2 × 10⁶	[1.1 × 10⁶, 8.3 × 10⁶]
$η$	CD8_E expansion	${cell}^{- 1} d^{- 1}$	2.5 × 10^–7	[1.6 × 10^–8, 6.7 × 10^–7]	3.5 × 10^–7	[2.3 × 10^–7, 5.2 × 10^–7]
$τ_{E}$	Delay in CD8_E expansion	$d$	3.6	[2.1, 5.9]	3.3	[2.6, 3.8]
$d_{E}$	CD8_E clearance	$d^{- 1}$	1.0	[5.1 × 10^–2, 2.0]	1.1	[3.3 × 10^–1, 2.5]
$ζ$	CD8_M generation	$C D 8_{M} C D 8_{E}^{- 1} d^{- 1}$	2.2 × 10^–1	[1.0 × 10^–2, 9.4 × 10^–1]	3.2 × 10^–2	[1.0 × 10^–2, 2.2 × 10^–1]
$τ_{M}$	Delay in CD8_M generation	$d$	3.5	[3.0, 4.0]	3.3	[2.1, 5.3]
${\hat{E}}_{0}$	Baseline CD8 or CD8_E	$CD8$ or ${CD8}_{E}$	4.2 × 10⁵	[3.3. × 10⁵, 5.3 × 10⁵]	6.0 × 10²	[3.1 × 10², 1.8 × 10²]
${\hat{E}}_{M 0}$	Baseline CD8_M	${CD8}_{M}$	-	-	6.5 × 10²	[3.2 × 10², 1.5 × 10³]
$T (0)$	Initial uninfected cells	$cells$	1 × 10⁷	-	1 × 10⁷	-
$I_{1} (0)$	Initial infected cells	$cells$	75	-	75	-
$I_{2} (0)$	Initial infected cells	$cells$	0	-	0	-
$V (0)$	Initial virus	${TCID}_{50}$	0	-	0	-
$E (0)$	Initial CD8_E	${CD8}_{E}$	0	-	0	-
$E_{M} (0)$	Initial CD8_M	${CD8}_{M}$	0	-	0	-

Plotting the model predicted dynamics for the lung-specific CD8⁺ T cells ( ${CD8}_{L} = E + E_{M}$ ) illustrated the accuracy of the model in predicting their dynamics within the lung parenchyma without fitting to these data (Figure 1C). One benefit of using the total CD8⁺ T cells is that the model automatically deduces the dynamics of effector-mediated killing and memory generation without needing to specify which phenotypes might be involved in these processes as they may be dynamically changing. However, the rates of expansion, contraction, and infected cell clearance may be different if only certain phenotypes are engaged. Thus, we examined whether the model could fit the dynamics of the predominant effector (CD25^–CD43⁺) and memory (CD43^–CD127⁺CD103⁺CD69⁺) phenotypes. Re-fitting the model to these data suggested that no changes to model formulation were needed and there were only small alterations to select CD8-specific parameter values (Figure 1D; Table 1).

Plotting the model ensembles revealed a correlation between the two infected cell clearance parameters ( $δ$ and $δ_{E}$ ; Figure 2B), which represent the efficacy of the non-specific immune response and the CD8⁺ T cell response, respectively. Performing a sensitivity analysis showed that the viral load dynamics do not change substantially when these parameters are increased or decreased (Appendix 2—figure 1) . However, decreasing the rate of non-specific infected cell clearance (i.e., lower $δ$ ) resulted in a significant increase in the number of CD8_E due to the small increase in the number of infected cells (Figure 1). Even with a larger CD8_E population, recovery was delayed by only ~0.1 d. Given the correlation between $δ$ and $δ_{E}$ (Figure 2B), a more efficient CD8_E response (i.e., higher $δ_{E}$ ) may be able to overcome this short delay in resolution. The lack of sensitivity to changes in the infected cell clearance parameters is in contrast to the DD model, where the viral dynamics were most sensitive to perturbations in $δ_{d}$ (Appendix 2—figure 1; Smith et al., 2018), which encompasses multiple processes. With CD8⁺ T cells explicitly included in the model, the infection duration was most sensitive to changes in the rate of CD8_E expansion ( $η$ ) (Figure 1, Appendix 1—figure 1; discussed in more detail below).

Examining the parameter ensembles and sensitivity analysis also yielded insight into how other model parameters affect the CD8⁺ T cell response. The rates of CD8_E expansion ( $η$ ) and clearance ( $d_{E}$ ) were slightly correlated, indicating a balance between these two processes (Figure 2A). This correlation and the sensitivity of $η$ produced model dynamics that were also sensitive to changes in the CD8_E clearance rate ( $d_{E}$ ) (Appendix 3—figure 2). As expected, the rates of CD8_M generation ( $ζ$ ) and CD8_E clearance ( $d_{E}$ ) were correlated (Figure 2A). It has been estimated that approximately 5–10% of effector CD8⁺ T cells survive to become a long-lasting memory population (Kaech et al., 2002). Despite the inability to distinguish between CD8_E and CD8_M in the total CD8⁺ T cell data, the model predicts that 17% of CD8_E transitioned to a memory class by 15 d pi. When considering only the CD25^–CD43^+/– effector and memory phenotypes, the model estimates this value to be ~7%. Additional insight into the regulation of the CD8⁺ T cell response, results from the model fitting, and a comparison of the DD model and the CD8⁺ T cell model are included in Appendix 2.

Density-dependent infected cell clearance

Given the accuracy of the model, we next sought to gain further insight into the nonlinear dynamics of CD8⁺ T cell-mediated infected cell clearance. Plotting the clearance rate ( $δ_{E} (I_{2}, E) = δ_{E} E / (K_{δ_{E}} + I_{2})$ ) as a function of infected cells ( $I_{2}$ ) and CD8_E ( $E$ ) (Figure 3) confirmed that there is minimal contribution from CD8_E-mediated clearance to viral load kinetics or infected cell kinetics prior to 7 d pi (Figure 3A–C, markers a–b). At the initiation of the second decay phase (7 d pi), the clearance rate is $δ_{E} (I_{2}, E) = 3.5 / d$ (Figure 3A–B, marker c). As the infected cell density declines toward the half-saturation constant ( $K_{δ_{E}} = 4.3 \times 10^{2} cells$ ), the clearance rate increases rapidly to a maximum of 4830/d (Figure 3A–C, markers d–g). The model predicts that there are 6 × 10⁵ infected cells remaining at 7 d pi, which can be eliminated by CD8_E in 6.7 h.

Figure 3

Download asset Open asset

Density-dependent infected cell clearance by CD8⁺ T cells and their impact on recovery time.

(A) The rate of CD8_E-mediated infected cell clearance ( $δ_{E} (I_{2}, E) = δ_{E} E / (K_{δ_{E}} + I_{2})$ ) plotted as a function of infected cells ( $I_{2}$ ) and effector CD8⁺ T cells ( $E$ ; CD8_E). The colored markers (denoted a–g) indicate the infected cell clearance rate that corresponds to different time points during the infection for the best-fit solution. (B) Values of $δ_{E} (I_{2}, E)$ for the indicated time points associated with the markers a–g. (C) Corresponding locations of the various $δ_{E} (I_{2}, E)$ values (markers a–g) on the best-fit solution of the CD8⁺ T cell model for virus ( $V$ ), infected cells ( $I_{2}$ ), and CD8_E ( $E$ ). (D) Solutions of the CD8⁺ T cell model (Equation (1)-(6)) for virus ( $V$ ) and total CD8⁺ T cells ( $\hat{E}$ ) using the best-fit parameters (black line) and when varying the CD8_E expansion rate ( $η$ ; magenta lines) to illustrate how different total CD8⁺ T cell magnitudes alter infection duration. The magenta lines are solutions from when the percent ${\hat{E}}_{\max}$ relative to ${\hat{E}}_{\max}$ from the best-fit solution was 42% (solid line), 39.2% (dash-dotted line), 39.1% (dashed line), or 37% (dotted line). (E) The time at which infected cells reach the half-saturation constant ( $I_{2} = K_{δ_{E}}$ ; gray circles) and the infection duration (time where $\log_{10} V = 0$ ; black diamonds) are shown for the various CD8⁺ T cell magnitudes. The gray line between these points is the time required to eliminate $K_{δ_{E}}$ infected cells and achieve complete resolution of the infection ( $\log_{10} V = 0$ ). (F) Fit of the CD8⁺ T cell model (Equation (1)-(6)) to viral loads and CD8⁺ T cells (magenta diamonds) following depletion at −2, 0, 3, and 7 d pi (magenta arrows). The best model (Supplementary file 2) resulted in fewer target cells ( $T_{0}$ ), a lower CD8_E influx ( $ξ$ ), and a higher CD8_E expansion rate ( $η$ ). All other parameters were fixed to the best-fit value in Table 1. The solid lines are $\hat{E} = E + E_{M} + {\hat{E}}_{0}$ and the dashed lines are $E$ for the cases where CD8⁺ T cells were depleted (magenta) and where they were not depleted (black). (G) Comparison of the log₁₀ ratio of virus to CD69⁺CD8⁺ T cells with and without CD8⁺ T cell depletion (magenta and black, respectively). All data are shown as mean ± standard deviation.

To explore how recovery time is altered by varying CD8_E levels, we examined the resulting dynamics from increasing or decreasing the rate of CD8_E expansion ( $η$ ). When $η$ was increased by 50%, the CD8_E population increased by a substantial 670% (Appendix 3—figure 1). However, this was insufficient to significantly shorten the infection (8.4 d versus 7.8 d). The infection duration could be reduced if CD8_E expansion began earlier (i.e., smaller $τ_{E}$ ; Appendix 3—figure 2). Although recovery is not significantly expedited by a larger CD8_E population, our model predicted that the infection would be dramatically prolonged if these cells are sufficiently depleted (Figure 3D–E and Figure 1). This in silico analysis revealed a bifurcation in recovery time, such that the infection is either resolved within ∼15 d pi or may last up to ∼45 d if CD8_E are below a critical magnitude required to resolve the infection (Figure 3D–E). The critical number of total CD8⁺ T cells needed for successful viral clearance was ${\hat{E}}_{m a x}^{c r i t} = 7.4 \times 10^{5}$ CD8, which was 39.2% of the maximum number of CD8⁺ T cells obtained from the best-fit solution (i.e., ${\hat{E}}_{m a x} = 1.9 \times 10^{6}$ CD8). This corresponds to 17% of CD8_E (i.e., $E_{m a x}^{c r i t} = 2.3 \times 10^{5}$ CD8_E; Figure 3D–E). The model analysis indicated that decreasing the total number of CD8⁺ T cells by as little as 0.1% from this critical level (i.e., 39.2% to 39.1%) lengthened the infection from 15 d pi to 25 d pi (Figure 3D).

Dynamical changes during CD8⁺ T cell depletion

To further test our model formulation and identify how viral load kinetics are altered with dynamically changing CD8⁺ T cells, we infected groups of mice with 75 TCID₅₀ and depleted these cells at −2 d, 0 d, 3 d, and 7 d pi with an anti-CD8α antibody (clone 2.43; Figure 3F). CD8⁺ T cells were reduced with over 99% efficiency and only 1.3% remained 2 d after depletion in the absence of infection (i.e., at 0 d pi; Figure 3F). CD8⁺ T cells remained >1 log₁₀ lower throughout the infection. Unexpectedly, the corresponding viral loads were significantly lower at 2 d pi (p=0.02) and consistently lower at other time points early in the infection (4 d pi (p=0.1) and 6 d pi (p=0.2)). As expected and predicted by our mathematical model, viral loads were significantly higher at 8 d pi (4.68 log₁₀ TCID₅₀ compared to 1.42 log₁₀ TCID₅₀; p=0.01). By 10 d pi, a sufficient number of CD8⁺ T cells were present and all animals had cleared the infection (Figure 3F). Interestingly, the number of CD8⁺ T cells at 10 d pi was only slightly higher than at 8 d pi (5.25 log₁₀ versus 5.02 log₁₀; p=0.064) further highlighting the density-dependent dynamics described above.

Given that viral loads were lower at early time points and that the anti-CD8α antibody is known to cause concentration-dependent changes in CD8_E differentiation, activation, and efficacy (Cross et al., 2019) in addition to resulting in death of the cells that would initiate activation of and removal by other immune cells, we did not expect our model to match the data without modulation of parameter values. However, we did expect that no changes to the model formulation would be required. In total, we tested >30 ‘models’ where 1–4 parameters were altered and found one model that was significantly better according to the AIC (Supplementary file 1). In that model, the initial number of target cells ( $T_{0}$ ) was 2.5x lower (~4 × 10⁶ versus 1 × 10⁷ cells), the rate of initial CD8_E influx ( $ξ$ ) was 2x lower (1.3 × 10⁴ versus $2.6 \times 10^{4} {C D 8}_{E}^{2} {c e l l}^{- 1} d^{- 1}$ ), and the rate of CD8_E expansion ( $η$ ) was 4x higher (1 × 10^–6 versus $2.5 \times 10^{- 7} {c e l l}^{- 1} d^{- 1}$ ). The second best model had a lower cost but was penalized by an additional parameter. That model suggested similar results but replaced the effect on $T_{0}$ with a combination of a lower virus production rate ( $p$ ) and higher non-specific infected cell clearance rate ( $δ$ ). Both models resulted in fewer infected cells and, thus, have approximately equivalent interpretations.

To further examine these findings, we plotted the ratio of virus to activated (CD69⁺) CD8⁺ T cells (Figure 3G). At 2 d pi, the ratio was similar for the depleted and mock-treated groups (p=0.41) suggesting that the depletion-induced reduction in virus was proportional to the reduction in activated CD8⁺ T cells. However, at 4 d pi, the ratio in the depleted groups was significantly higher (p=0.0028) for otherwise similar levels of virus. This signifies that there were disproportionately low numbers of CD8⁺ T cells and, thus, reduced influx (i.e., lower $ξ$ ). By 6 d pi, the ratios were again similar (p=0.25) implying a higher expansion of these cells (i.e., higher $η$ ). Investigation into the predicted target cell reduction is included below.

Modeling lung injury dynamics

To investigate the dynamics of infected cells and their clearance by CD8⁺ T cells, we quantified these cells and the progression and resolution of lung injury using serial whole lung histomorphometry (Figure 4A). Antigen-positive areas of the lung (‘active’ lesions) were first detectable at 2 d pi (Figure 4A–B), which coincides with the peak in viral loads. The infected area continued to increase in a nonlinear manner until 6 d pi, whereas viral loads remained high until 7 d pi (Figure 1B). At this time, resolution of the infection began and the infected area declined at a rate of ∼28.7%/d between 6–7 d pi (Figure 4—figure supplement 1). Few to no infected cells were present at 8 d pi (Figure 4A). Correspondingly, virus was undetectable in most animals by 8 d pi (Figure 1B). Because the percent active lesion is a reflection of the influenza-positive cells, we examined whether the CD8⁺ T cell model accurately predicted these dynamics. In the model, the accumulated infection is defined by the cumulative area under the curve (CAUC) of the productively infected cells ( $I_{2}$ ). Plotting the percent active lesion against the CAUC of $I_{2}$ shows that the model accurately reflects the cumulative infected cell dynamics and, thus, the infection progression within the lung (Figure 4B). Plotting the CAUC of $I_{2}$ for all parameter sets in the 95% CIs and from fitting the model to specific phenotypes further illustrates the accuracy by showing that the heterogeneity in the histomorphometry data is captured (Figure 4B), which is larger than the heterogeneity in viral loads (Figure 1B). Plotting the predicted active lesion dynamics for the model parameterized to the data where CD8⁺ T cells were depleted suggested that there was a ∼22% reduction in the active lesion (Figure 4B).

Figure 4 with 2 supplements see all

Download asset Open asset

Histomorphometry and inflammation of the IAV-infected lung.

(A) Whole lung sections with histomorphometry showing the areas of influenza NP-positive ‘active’ lesions (red), previously infected ‘inactive’ lesions with minimal antigen-positive debris (green), or mixed active and inactive regions (orange) throughout the infection. Representative images from each group are shown. (B) Percent active lesion (red squares) plotted together with the cumulative area under the curve (CAUC) of the predicted infected cell dynamics ( $I_{2}$ ) obtained from fitting the CD8⁺ T cell model. The linear decline in the active lesion (−28.7%/d; see Figure 4—figure supplement 1) was used to estimate the decline after 6 d pi. (C) Percent inactive lesion (green squares) plotted together with the percent maximum CD8_E ( $E / E_{\max}$ ) obtained from fitting the CD8⁺ T cell model. (D) The total lesion (blue squares) is the addition of the active and inactive lesions. To include all measurements on the same scale, the CAUC of $I_{2}$ was multiplied by a scaling factor of 14.2% per 1 × 10⁷ cells, and the percent maximum CD8_E was multiplied by a scaling factor of 0.46%. (E) Fit of Equation (7) to the alveolar (white triangles) and interstitial (yellow triangles) inflammation scores. The solid black, dashed black, and solid magenta lines are the curves generated using the best-fit parameters obtained from fitting the model to the total CD8⁺ T cells, the CD25^–CD43⁺ and CD43^–CD127⁺CD103⁺CD69⁺ CD8⁺ T cells, and the total CD8⁺ T cells during CD8 depletion, respectively. The gray shading are the curves generated using the 95% CI parameters from fitting the model to the total CD8⁺ T cells. All data are shown as mean ± standard deviation.

Antigen-negative, previously-infected or damaged areas of the lung (‘inactive’ lesions) are evident beginning at 5 d pi (Figure 4A,C). This resolution of the infection continued from 5-8 d pi, causing a 15.1%/d increase in the area of inactive lesions (Figure 4—figure supplement 1). Following this, healing of the injured lung is apparent as the inactive lesioned area declines (−14.5%/d from 9–10 d pi; Figure 4C and Figure 4—figure supplement 1). These dynamics generally parallel the CD8⁺ T cell dynamics but are nonlinearly correlated (Figure 4—figure supplement 1). Fitting a line to the CD8⁺ T cell data from 5-8 d pi indicated that the influx rate of all phenotypes is 4.94 × 10⁵ cells/d, of lung-specific phenotypes is 3.97 × 10⁵ cells/d, and of the CD25^–CD43⁺ effector phenotype is 1.98 × 10⁵ cells/d (Figure 4—figure supplement 1). Thus, the model estimates that, on average, every 100,000 total, lung-specific, or CD25^–CD43⁺ CD8⁺ T cells clear ~3.1%, ~3.8%, or ~7.6% of the infected areas within the lung, respectively. During the CD8⁺ T cell contraction phase, a similar linear regression analysis suggested that these cells decline at rates of ~ 4.13 × 10⁵ CD8/d (total), ~ 2.35 × 10⁴ CD8/d (lung-specific), and ~ 3.83 × 10⁴ CD8/d (CD25^–CD43⁺) (Figure 4—figure supplement 1). Similar to the relation discussed above, the dynamics of the damaged areas of the lung corresponded to the dynamics of the percent maximum CD8_E (i.e., $E / E_{m a x}$ ) in the model (Figure 4C). Our model suggested a ∼37% reduction in the inactive lesion during CD8⁺ T cell depletion (Figure 4C). Adding the predicted dynamics for the active and inactive lesions agrees with the dynamics of the total lesion, but is slightly underestimated when using the model fit to CD25^–CD43⁺ CD8⁺ T cells (Figure 4D). In addition, the predicted total lesion was reduced for the case where CD8⁺ T cells were depleted (Figure 4D).

Modeling lung inflammation dynamics

In addition to measuring virus-induced lung damage, lung inflammation was scored (Figure 4E). Both alveolar and interstitial inflammation begin to increase at 2 d pi with the sharpest increase between 1–3 d pi. Inflammation continues to increase until 6 d pi with a maximum score of 3.5–4.0 out of 5. Resolution of inflammation was not evident during the time course of our data, which concluded at 10 d pi. This is in contrast to the lung damage inflicted by the virus, which begins declining at 8 d pi and shows that ∼15% of the damage was repaired by 10 d pi. Inflammation was linearly correlated to inflammatory macrophages and to neutrophils (Figure 4—figure supplement 2), which were excluded from the model. However, the knowledge that the model’s predicted infected cell dynamics are accurate and that these cells can produce cytokines and chemokines that attract cells like macrophages and neutrophils suggested that our model could be used to estimate inflammation. To do this, we fit Equation (7) to the inflammation scores while keeping all other parameters fixed to their best-fit values (Table 1). The results suggested that the equation captured the inflammation dynamics, and that the contribution from the initial infected cell class ( $I_{1}$ ) was $α_{1} = 4.27 per 10^{7} cells / d$ and the contribution from the productively infected cells ( $I_{2}$ ) was $α_{2} = 0.87 per 10^{7} cells / d$ . For the case where CD8⁺ T cells were depleted, Equation (7) estimated the inflammation scores to be reduced to ∼1.5 (Figure 4E).

Weight loss relates nonlinearly to lung injury and inflammation

To monitor disease progression, weight loss was measured daily throughout the course of infection (Figure 5). During the first 5 d pi, animals gradually lost ∼4% of their initial weight. This was followed by a sharper drop (8%) at 6 d pi. Animal weights increased slightly at 7 d pi (∼6%) before reaching peak weight loss (10–14%) at 8 d pi. Following virus resolution, the animals’ weights began to restore as the inactive lesions resolved (9–10 d pi; Figure 5A). Weight loss was reduced during CD8⁺ T cell depletion (Figure 5A). Plotting weight loss against the percent total (active and inactive) lesioned area of the lung shows the similarity in their dynamics (Figure 5A) and revealed a nonlinear relation (Figure 5B). To further quantify their relationship, we fit the saturating function $L (W) = l_{m a x} W^{n} / (W^{n} + K_{w}^{n})$ to these data, where $L$ is the percent total lesioned area, $W$ is percent weight loss, l_max is the maximum rate of the interaction, $K_{w}$ is the half-saturation constant, and $n$ is the Hill coefficient. This function provided a close fit to the data ( $R^{2} = 0.92$ ; black line in Figure 5B) with best-fit parameters $l_{m a x} = 39.7 %$ total lesioned area, $K_{w} = 2.58 %$ weight loss, and $n = 5.24$ . Plotting the estimated percent total lesion during CD8⁺ T cell depletion (Figure 4D) together with the measured weight loss (Figure 5A) also supported the nonlinear relation. Using the same best-fit parameters ( $K_{w} = 2.58 %$ weight loss and $n = 5.24$ ) together with the predicted maximum percent total lesion ( $l_{m a x} = 17.73 %$ total lesioned area) showed the accuracy of this function (magenta line in Figure 5D). Weight loss was also nonlinearly correlated to the alveolar and interstitial inflammation scores (Figure 5C), although the relation was slightly different compared to that for the lung injury data. Fitting the same function to these data independently for alveolar and interstitial inflammation also provided a close fit ( $R^{2} = 0.98$ and $R^{2} = 0.97$ ; black lines in Figure 5D). Best-fit parameters for alveolar inflammation were $l_{a_{\max}} = 3.63$ score, $K_{a_{w}} = 1.95 %$ weight loss, and $n_{a} = 3.65$ . Best-fit parameters for interstitial inflammation were $l_{i_{\max}} = 3.40$ score, $K_{i_{w}} = 1.96 %$ weight loss, and $n_{i} = 3.15$ . For the CD8-depleted case, plotting the estimated inflammation (Figure 4E) together with the measured weight loss (Figure 5A) again supported the nonlinear relation. In addition, using the same best-fit parameters ( $K_{a_{w}} = 1.95 %$ weight loss, $n_{a} = 3.65$ ; $K_{i_{w}} = 1.96 %$ weight loss, $n_{i} = 3.15$ ) together with the predicted maximum inflammation score ( $l_{i_{\max}} = 1.41$ score) showed the accuracy (magenta line in Figure 5D).

Figure 5

Download asset Open asset

Weight loss dynamics and its relation to lung injury and inflammation.

(A) The percent total (active and inactive) lesion (blue squares) plotted together with the percent weight loss (black diamonds) to illustrate their similar dynamics. (B) Fit of a saturating function ( $L (W) = l_{\max} W^{n} / (K_{w}^{n} + W^{n})$ ; black line) to the mean percent total lesioned area ( $L$ ) and mean weight loss ( $W$ ) for all time points (black circles). The best-fit parameters were $l_{\max} = 39.7 %$ total lesioned area, $K_{w} = 2.6 %$ weight loss, and $n = 5.2$ . The magenta circles are the predicted percent total lesioned area (Figure 4D) with the corresponding weight loss during CD8 depletion, and the magenta line is the prediction using the best-fit parameters ( $K_{w} = 2.6 %$ weight loss and $n = 5.2$ ) together with the maximum predicted percent total lesion ( $l_{\max} = 17.73 %$ total lesioned area). (C) The alveolar (white triangles) and interstitial (yellow triangles) inflammation plotted together with the percent weight loss (black diamonds). (D) Fit of a saturating function (black line) to the mean alveolar (white circles) and interstitial (yellow circles) inflammation scores and mean weight loss for all time points. The best-fit parameters for alveolar inflammation (white circles; black line) were $l_{a_{m a x}} = 3.63$ score, $K_{a_{w}} = 1.95 %$ weight loss, and $n_{a} = 3.65$ , and for interstitial inflammation were $l_{i_{m a x}} = 3.40$ score, $K_{i_{w}} = 1.96 %$ weight loss, and $n_{i} = 3.15$ . The magenta circles are the predicted percent inflammation score (Figure 4E) with the corresponding weight loss during CD8 depletion, and the magenta line is the prediction using the best-fit parameters ( $K_{a_{w}} = 1.95$ % weight loss, $n_{a} = 3.65$ ; $K_{i_{w}} = 1.96 %$ weight loss, $n_{i} = 3.15$ ) together with the maximum predicted inflammation ( $l_{m a x} = 1.41$ score).

Discussion

Influenza A virus infections pose a significant threat to human health, and it is crucial to understand and have tools that can predict how the virus spreads within the lower respiratory tract, how specific immune responses contribute to infection control, and how these relate to disease progression. Although it has been difficult to directly relate these features and obtain high quality data from the lower respiratory tract in humans, we circumvented the challenge by pairing comprehensive experimental data with robust mathematical models and analyses. Our iterative model-driven experimental approach (Smith, 2018a; Smith, 2018b) revealed important dynamic relations between virus, infected cells, immune cells, lung damage, inflammation, and disease severity (summarized in Figure 6). Identifying these nonlinear connections allows for more accurate interpretations of viral infection data and significant improvement in our ability to predict disease severity, the likelihood of complications, and therapeutic efficacy.

Figure 6

Download asset Open asset

Summary of the connections between the kinetics of virus, infected cells, CD8⁺ T cells, lung injury, lung inflammation, and disease severity.

Summary of the relations between the dynamics of virus, infected cells, CD8⁺ T cells, lung lesions, lung inflammation, and weight loss established by our analysis. Given that viral loads and weight loss are the most easily measured variables, our mathematical model (Equation (1)-(6)) could be used to estimate the kinetics of infected cells, CD8⁺ T cells, and inflammation. The cumulative area under the curve (CAUC) of the productively infected cell dynamics ( $I_{2}$ ) yields an estimate of the percent lung infected (active lesion) while the predicted relative CD8_E dynamics ( $E / E_{m a x}$ ) yields an estimate of the percent lung resolved (inactive lesion). The total amount of lung involved (% lung infected and % lung resolved) or the inflammation scores can then be used to estimate weight loss through the functions $L_{L} (W)$ and $L_{I} (W)$ . These connections could be reversed and weight loss used to estimate viral load kinetics.

Our histomorphometric data and analyses provided a robust quantification of the extent of influenza virus infection within the lung and the dynamics of lung injury and inflammation. The images look strikingly similar to CT scans and postmortem histopathology from infected patients, which also show extensive bilateral damage and cytopathic effects in epithelial cells (Kohr et al., 2010; Li et al., 2011; Ajlan et al., 2009; Soto-Abraham et al., 2009; Gill et al., 2010; Li et al., 2012; Shieh et al., 2010). In addition, our results agree with weekly CT imaging of hospitalized patients infected with the 2009 H1N1 influenza virus showing that the damage worsens during the first week of illness and that cell proliferation begins in the second week (Li et al., 2012). Further investigating the utility of CT imaging could enhance the translatability of our findings. The model’s 95% CI predictions of the lesioned area nicely captured the heterogeneity in the data (Figure 4), which is remarkable given the minimal variability in the virus and CD8⁺ T cell data and model fit (Figure 1). The differential heterogeneity suggests that small changes to lung viral loads can induce significant changes in disease severity. Indeed, this has been observed in both humans and animals when antivirals are administered (Toniolo Neto, 2014; Baloxavir Marboxil Investigators Group et al., 2018; Treanor et al., 2000; Moscona, 2005; Gubareva et al., 2000; Koshimichi et al., 2018) or when host responses are experimentally modulated (Channappanavar et al., 2016; Szretter et al., 2009; Iwasaki and Nozima, 1977). It was also observed in our CD8 depletion data, where our model predicted a significant reduction in the lesioned area (Figure 4) while viral loads were only slightly lower (Figure 3). This may indicate that viral loads are generally not a reliable measure of infection severity, and the knowledge should aid future studies seeking to estimate therapeutic efficacy or the effects of immune modulation. There was a slight over-estimation of the inactive lesion early in the infection (1–5 d pi, dark gray shading; Figure 4) when using the fit to the total CD8⁺ T cells, which was abrogated when using the CD25^-CD43⁺CD8⁺ T cells. However, these cells underestimated the total lesion at 6–7 d pi, suggesting contributions from other CD8⁺ phenotypes. In general, modeling different CD8⁺ T cell phenotypes may help clarify their in vivo functions and better define the rates of CD8-mediated infected cell clearance. Although our results suggest that this is not necessarily required, they may be more critical to consider when predicting the response of preexisting and potentially cross-reactive immunity particularly given that memory resident T cells are only beginning to be understood (Zens et al., 2016; Slütter et al., 2017; Van Braeckel-Budimir et al., 2018; Wakim et al., 2015; Uddbäck et al., 2021).

The new knowledge about the extent of infection within the lung also uncovered the nonlinear connection between disease severity and various measures of lung pathology (Figure 5). This discovery is significant because it suggests there is a unique link to the extent of both lung injury and inflammation, which are often used interchangeably yet have distinct dynamics (Figure 4) and correlate to different immune responses (Figure 4—figure supplements 1–2). In the same vein, both measurements can be estimated using the dynamics of infected cells (Figure 4) but with independent mathematical relations (CAUC of $I_{2}$ versus Equation (7)). Our CD8 depletion data (Figure 3) and other experimental studies using histomorphometry (Marathe et al., 2016) corroborate a relationship between weight loss and lung pathology during IAV infection. For example, animals treated with antivirals in various conditions (single or combination therapy and in immunocompetent or immunosuppressed hosts) demonstrated that, although viral loads are not always significantly reduced, the antiviral-induced reductions in weight loss were paired with decreased infected areas of the lung (Marathe et al., 2016). Further examining lung injury and inflammation kinetics and their connection to weight loss in various experimental infection settings (e.g., different ages or sex, under therapy, other strains or viruses, etc.) and deciphering how innate and humoral responses contribute to these measurements should improve our predictive capabilities. In addition, the nonlinearity of the connections warrants further investigation as it could indicate cooperativity of multiple mechanisms. Studies such as these may be particularly helpful in understanding the exacerbated morbidity and mortality in elderly individuals, where there are lower viral loads, increased weight loss, symptoms, and/or mortality within animal models (Toapanta and Ross, 2009; Smith et al., 2019). Because murine models are imperfect systems with important immunologic and physiologic distinctions from humans (Mestas and Hughes, 2004), uncovering how our results translate to human infection with novel or circulating influenza viruses and identifying the corresponding, tractable human symptom that could be a robust predictor of disease would be ideal. However, this could be arduous due to inherent human and viral heterogeneity and contributions from co-morbidities and complex immunologic histories.

Despite their limitations, animal models are used to make viral and immunological discoveries and to develop and test therapeutics and vaccines (reviewed in Thangavel and Bouvier, 2014; Margine and Krammer, 2014; Smith and McCullers, 2014), and this work demonstrates a significant potential for the easily obtained weight loss data to be used, analyzed, and interpreted within both modeling and experimental studies. In our data and others’ animal data, there is a spike in weight loss, symptom score, and/or inflammation at ∼6 d pi (Music et al., 2014; Kugel et al., 2009; Lu et al., 2018; Huang et al., 2017; Figure 5). Our data suggests that this may be due to CD8_E activity while the infection continues to spread. The increase in animal weights starting at 7 d pi is concurrent with the decline in infected lesions (Figure 5). Interestingly, there was no corresponding decrease in inflammation, and the spike was not apparent in the CD8 depletion data where severity was reduced. The subsequent increase in weight loss following infection resolution (8–9 d pi; Figure 4C) could be attributed to CD8⁺ T cell-mediated pathology and/or to ongoing inflammation. The tighter correlation to inflammation might suggest the latter; however, some human and animal studies indicate that large numbers of CD8⁺ T cells pose a risk of acute lung tissue injury (Duan and Thomas, 2016; Moskophidis and Kioussis, 1998; Mauad et al., 2010; Rygiel et al., 2009; La Gruta et al., 2007). According to our model predictions, excessive CD8⁺ T cell numbers may augment disease progression yet do not improve recovery time (Appendix 3—figure 2). Instead, an earlier onset of CD8_E proliferation (i.e., smaller $τ_{E}$ ) would be required to significantly shorten the infection (Appendix 3—figure 2). This aligns with evidence that hosts with adaptive immune responses primed by vaccine or prior infection recover more rapidly (Duan et al., 2015; Weinfurter et al., 2011). While higher CD8⁺ T cell numbers have little impact on viral kinetics, our model and data agree with clinical and experimental studies from a wide range of host species that impaired CD8⁺ T cell responses can prolong an IAV infection (Wang et al., 2015; Miao et al., 2010; McMichael et al., 1983; Yap and Ada, 1978; Wells et al., 1981; Bender et al., 1992). In some scenarios, such as in immunocompromised hosts, virus can persist for up to several weeks if CD8⁺ T cell-mediated clearance is unsuccessful (Figure 3; Wang et al., 2015; Eichelberger et al., 1991; Hou et al., 1992; Xue et al., 2017; Memoli et al., 2014). The bifurcation in recovery time revealed by the model suggests that this may occur when the number of lung-specific CD8⁺ T cells are less than 2.31 × 10⁵ cells (Figure 3D–E). Our CD8 depletion data show that the precise number will be dependent on other infection variables. Minimally, we would expect this number to vary depending on parameters like the dose, rate of virus replication, degree of prior immunity, and/or the infected cell number and lifespan, which has been noted in another modeling study that detailed similar bifurcating behavior (Yates et al., 2011). Although some previously published models also suggested delayed resolution with depleted CD8⁺ T cell responses (Dobrovolny et al., 2013; Miao et al., 2010; Cao et al., 2016; Lee et al., 2009), this bifurcation has not been observed and their estimated delays in recovery do not amount to the long-lasting infections in immunodeficient patients (Wang et al., 2015; Eichelberger et al., 1991; Hou et al., 1992). Our model’s ability to capture the dynamics when CD8⁺ T cells are depleted is encouraging, and the data are an important reminder that experimental modulation of cell populations (e.g., through depletion or genetic knockouts) is complicated and that the data from such systems should be interpreted cautiously. It also brings into question studies that have used these types of data without validated mathematical models or quantification of other immunological variables.

Our data and analyses provided strong support that the infected cell density impacts the rate at which they are cleared by effector CD8⁺ T cells (CD8_E) (Figures 3–4, Figure 4—figure supplement 1). Although CD8⁺ T cell dynamics can be somewhat replicated when assuming the density dependence lies within their expansion, that assumption could not recover the viral load dynamics (Appendix 1—figure 1). Discriminating between these mechanisms is difficult in vivo, but the ability of the model in Equation (1)-(6) to capture the entire time course of CD8_E dynamics in multiple scenarios is compelling. Regardless of the mechanism, we first detected this density-dependence in a model that excluded specific immune responses (see Appendix 2) (Smith et al., 2018). Simple models like that one are useful to capture nonlinearities, but they cannot distinguish between different mechanisms.

Several factors may contribute to the density-dependent change in the rate of CD8_E-mediated clearance. One possibility is that the slowed rate of clearance at high infected cell densities is due to a ‘handling time’ effect, which represents the time required for an immune cell to remove an infected cell (e.g., as in Smith et al., 2018; Li and Handel, 2014; Le et al., 2014; Graw and Regoes, 2009; Gadhamsetty et al., 2014; Pilyugin and Antia, 2000; Smith et al., 2011b). When CD8_E interact with infected cells, a complex is formed for ∼20–40 min (Mempel et al., 2006; Deguine et al., 2010; Zagury et al., 1975; Yannelli et al., 1986; Perelson and Bell, 1982). Because CD8_E could not interact with other infected cells during this time, the global rate of infected cell clearance would be lowest when infected cells outnumber CD8_E. In addition, contact by a single CD8_E may be insufficient to remove an infected cell (Halle et al., 2016). Infected cell clearance is more frequently observed after interactions with multiple CD8_E, with an average of 3.9 contact events either serially or simultaneously (Halle et al., 2016). Thus, the high density of infected cells early in the infection reduces the probability that a single virus-infected cell would be targeted a sufficient number of times to induce cell death. However, as CD8_E accumulate and the density of infected cells decreases (Figure 3A), the probability of simultaneous interactions will increase. These interactions may also be influenced by CD8_E movement, where their velocity slows at the height of infected cell clearance (6–8 d pi) (Lambert Emo et al., 2016; Matheu et al., 2013). This should reduce the time required to remove an infected cell and, thus, result in a higher efficiency. Moreover, it is possible that spatial limitations also contribute, such that a high infected cell density may hinder CD8⁺ T cells from reaching infected cells on the interior of the infected tissue. Crowding of immune cells at the periphery of infected lesions has been observed in other infections (e.g., in tuberculosis granulomas [Ulrichs et al., 2004; Russell et al., 2009]) and has been suggested in agent-based models of influenza virus infection (Levin et al., 2016).

In addition to illuminating the connections between viral kinetics and pathology, the histomorphometric data validated the model’s infected cell dynamics (Figure 4B). The dynamics of susceptible and infected cells throughout the infection and the accuracy of the target cell limited approximation used within influenza viral kinetic models have been questioned for several years (Smith and Perelson, 2011; Smith, 2018b; Beauchemin and Handel, 2011; Ahmed et al., 2017; Gallagher et al., 2018; Boianelli et al., 2015; Handel et al., 2018). The ability of the model to accurately predict the histomorphometry and CD8⁺ T cell depletion data corroborates the use of this approximation, which assumes a limited number of available target cells and describes their decline by infection only. Although the slowing of viral loads beginning at 2 d pi could be due to a variety of innate immune-mediated mechanisms (e.g., macrophages, neutrophils, type I interferons, etc.), adding more complex dynamics to our model was unnecessary to describe the data. Further, the model and data agree that there are few infected cells during the time when viral loads are growing most rapidly (0–2 d pi; Figures 1B and 4B). We previously used this information to derive approximations for the model and gain deeper insight into how each parameter influences the kinetics (Smith et al., 2010), which has helped numerous studies interpret their results (Miao et al., 2010; Smith et al., 2011a; Smith, 2018b; Holder and Beauchemin, 2011). The data also supports the model’s hypothesis that that there is minimal clearance of infected cells prior to CD8_E expansion (Figure 4). The knowledge of the model’s accuracy and of the spatiality in the lung should aid investigation into the mechanisms that limit virus growth during the early stages of the infection.

Employing targeted model-driven experimental designs to examine and validate theoretical predictions like the ones presented here is pivotal to elucidating the mechanisms of infection spread and clearance (Smith, 2018b). Examining other infections (e.g., coronaviruses), modifications to the dynamics (e.g., lethal doses), and the connection between lung measurements and those more easily acquired from upper respiratory tract will help refine the dynamical links between virus, host immune responses, and disease severity and identify their generalizability. Determining the factors that influence disease severity is vital to understanding the disproportionate mortality in at-risk populations (e.g., elderly) and to improving therapeutic design. This is particularly important because current antivirals alleviate symptoms but do not always effectively lower viral loads (Toniolo Neto, 2014; Baloxavir Marboxil Investigators Group et al., 2018; Treanor et al., 2000; Moscona, 2005; Gubareva et al., 2000; Koshimichi et al., 2018). The predictive capabilities of validated models like the one here should prove useful in forecasting infection dynamics for a variety of scenarios. These tools and analyses provide a more meaningful interpretation of infection data, new ways to interpret weight loss data in animal models, and a deeper understanding of the progression and resolution of the disease, which will undoubtedly aid our ability to effectively combat influenza.

Materials and methods

Ethics statement

Request a detailed protocol

All experimental procedures were performed under protocols O2A-020 or 17–096 approved by the Animal Care and Use Committees at St. Jude Children’s Research Hospital (SJCRH) or the University of Tennessee Health Science Center (UTHSC), respectively, under relevant institutional and American Veterinary Medical Association (AVMA) guidelines. All experimental procedures were performed in a biosafety level two facility that is accredited by the American Association for Laboratory Animal Science (AALAS).

Mice

Adult (6-week-old) female BALB/cJ mice were obtained from Jackson Laboratories (Bar Harbor, ME) or Charles River Laboratories (Wilmington, Massachusetts). Mice were housed in groups of five mice in high-temperature 31.2 cm × 23.5 cm × 15.2 cm polycarbonate cages with isolator lids (SJCRH) or in 38.2 cm × 19.4 cm × 13.0 cm solid-bottom polysulfone individually ventilated cages (UTHSC). Rooms used for housing mice were maintained on a 12:12 hr light:dark cycle at 22 ± 2°C with 50% humidity in the biosafety level two facility at SJCRH (Memphis, TN) or UTHSC (Memphis, TN). Prior to inclusion in the experiments, mice were allowed at least 7 days to acclimate to the animal facility such that they were 7 weeks old at the time of infection. Laboratory Autoclavable Rodent Diet (PMI Nutrition International, St. Louis, MO; SJCRH) or Teklad LM-485 Mouse/Rat Sterilizable Diet (Envigo, Indianapolis, IN; UTHSC) and autoclaved water were available ad libitum. All experiments were performed under an approved protocol and in accordance with the guidelines set forth by the Animal Care and Use Committee at SJCRH or UTHSC.

Infectious agents

Request a detailed protocol

All experiments were done using the mouse adapted influenza A/Puerto Rico/8/34 (H1N1) (PR8).

Infection experiments

Request a detailed protocol

The viral infectious dose (TCID₅₀) was determined by interpolation using the method of Reed and Muench, 1938 using serial dilutions of virus on Madin-Darby canine kidney (MDCK) cells. Mice were intranasally inoculated with 75 TCID₅₀ PR8 diluted in 100 μl of sterile PBS. In a subset of animals, CD8⁺ T cells were depleted by IP injection at days −2, 0, 3, and 7 pi with 100 μg of the rat anti-CD8α antibody (clone 2.43) that was purified from ATCC hybridoma (per manufacturer instructions) in 250 μl of PBS. Depletion efficiency was confirmed in the lung and spleen as described below. Mice were weighed at the onset of infection and each subsequent day to monitor illness and mortality. Mice were euthanized if they became moribund or lost 30% of their starting body weight. For viral load and total CD8⁺ T cell quantification, experiments were repeated three times and in each facility to ensure reproducibility and two complete experiments (ten animals per time point) were used for these studies. For all other experiments, the data was compared to prior results in addition to being repeated at select time points to ensure reproducibility, and one complete experiment (five animals per time point) was used. For pathology scoring and histomorphometry quantitation, five animals per time point were used. Power was calculated using G*Power.

Lung harvesting for viral and cellular dynamics

Request a detailed protocol

For total CD8⁺ T cell quantification, mice were euthanized by CO₂ asphyxiation. To distinguish blood-borne CD8⁺ T cells from those in the lung parenchyma, deeply anesthetized mice (5% isoflurane) were retro-orbitally injected with 3 μg of anti-CD45 antibody (PerCP, clone 30-F11, Biolegend) 3 min prior to euthanasia (Anderson et al., 2012; Keating et al., 2018). Mice were then euthanized by 33% isoflurane inhalation and their lungs perfused with 10 ml PBS prior to removal. For all experiments, lungs were then aseptically harvested, washed three times in PBS, and placed in 500 μl sterile PBS. Whole lungs were digested with collagenase (1 mg/ml, Sigma C0130), and physically homogenized by syringe plunger against a 40 μm cell strainer. Cell suspensions were centrifuged at 4°C, 500xg for 7 min. The supernatants were used to determine the viral titers (TCID₅₀) by serial dilutions on MDCK monolayers. Following red blood cell lysis, cells were washed in MACS buffer (PBS, 0.1M EDTA, 0.01M HEPES, 5 mM EDTA and 5% heat-inactivated FBS). Cells were then counted with trypan blue exclusion using a Cell Countess System (Invitrogen, Grand Island, NY) and prepared for flow cytometric analysis as indicated below.

Lung titers

Request a detailed protocol

For each animal, viral titers were obtained using serial dilutions on MDCK monolayers and normalized to the total volume of lung homogenate supernatant. The resulting viral loads are shown in Figure 1B, provided in Source data 1, and were previously published and utilized for calibration of the density-dependent model (see Appendix 2; Smith et al., 2018).

Flow cytometric analysis

Request a detailed protocol

Flow cytometry (LSRII, BD Biosciences, San Jose, CA (SJCRH) or ZE5 Cell Analyzer, Bio-Rad, Hercules, CA (UTHSC)) was performed on the cell pellets after incubation with 200 μl of a 1:2 dilution of Fc block (human-γ globulin) on ice for 30 min, followed by viability (Biolegend, Zombie Violet Fixable Viability Kit) and surface marker staining with anti-mouse antibodies. For total CD8⁺ T cell, macrophage, and neutrophil quantification, we used antibodies CD3e (BV786, clone 145–2 C11, Biolegend), CD4 (PE-Cy5, clone RM4-5, BD Biosciences), CD8α (BV605, clone 53–6.7, BD Biosciences), Ly6C (APC, clone HK1.4, eBioscience), F4/80 (PE, clone BM8, eBioscience), CD11c (eFluor450, clone N418, eBioscience), CD11b (Alexa700, clone M1/70, BD Biosciences), MHC-II (FITC, clone M5/114.15.2, Invitrogen), CD49b (APCe780, clone DX5, eBioscience), and Ly6G (PerCp-Cy5.5, clone 1A8, Biolegend). The data were analyzed using FlowJo 10.4.2 (Tree Star, Ashland, OR) where viable cells were gated from a forward scatter/side scatter plot and singlet inclusion. Following neutrophil (Ly6G^hi) and subsequent macrophage (CD11c^hiF4/80^hi) exclusion, CD8⁺ T cells were gated as CD3e⁺DX5^–CD4^–CD8⁺.

To distinguish blood borne CD8⁺ T cells from those in the lung parenchyma, we used antibodies CD3e (BV786, clone 145–2 C11, Biolegend), CD4 (FITC, clone RM4-5, Biolegend), CD8α (BV605, clone 53–6.7, BD Biosciences), B220 (APCe780, clone RA3-6B2, eBioscience), CD49b (APCe780, clone DX5, eBioscience), CD62L (PE-Cy7, clone MEL-14, Biolegend), CD69 (PE, clone H1.2F3, Biolegend), CD44 (PE-Dazzle594, clone IM7, Biolegend), CD25 (Alexa700, clone PC61, Biolegend), CD43 (APC, clone 1B11, Biolegend), CD127 (PE-Cy5, clone A7R34, Biolegend), and CD103 (BV711, clone 2E7, Biolegend) or CD314 (NKG2D) (BV711, clone CX5, BD Bioscience). The data were analyzed using FlowJo 10.6.2 (Tree Star, Ashland, OR) where viable cells were gated from a forward scatter/side scatter plot, singlet inclusion, and viability dye exclusion. Blood borne CD8⁺ T cells were gated as CD45⁺ and those in the lung parenchyma as CD45^-. The total in each sub-population was gated as CD3e⁺B220^–DX5^–CD4^–CD8⁺. IAV-specific CD8⁺ T cells were then gated as CD25⁺CD43⁺ (recently activated), CD25^–CD43⁺ (effector) (Keating et al., 2018), and CD43^–CD127⁺CD103⁺CD69⁺ (long-lived memory). Expression of CD44, CD69, CD62L, and NKG2D were also assessed to ensure appropriate classification.

For all experiments, the absolute numbers of CD8⁺ T cells were calculated based on viable events analyzed by flow cytometry as related to the total number of viable cells per sample. The data are provided in Source data 1. We use the kinetics of the total number of CD8⁺ T cells (non-perfused lung; ‘total’), the total number of CD8⁺ T cells in the lung parenchyma (perfused lung; ‘lung-specific’), and virus-specific CD8⁺ T cell phenotypes in the lung parenchyma as defined by surface staining (described above). We chose this approach because the use of tetramers yields epitope-specific cell dynamics that vary in time and magnitude (e.g., as in Toapanta and Ross, 2009; Keating et al., 2018) and would complicate the model and potentially skew the results. In addition, tetramer staining is not available for all epitopes in BALB/cJ mice and they may underestimate the dynamics of IAV-specific cells (Keating et al., 2018). The gating strategies are shown in Figure 1—figure supplement 2.

Lung immunohistopathologic and immunohistochemical (IHC) evaluation

Request a detailed protocol

The lungs from IAV-infected mice were fixed via intratracheal infusion and then immersion in 10% buffered formalin solution. Tissues were paraffin embedded, sectioned, and stained for influenza virus using a primary goat polyclonal antibody (US Biological, Swampscott, MA) against influenza A, USSR (H1N1) at 1:1000 and a secondary biotinylated donkey anti-goat antibody (sc-2042; Santa Cruz Biotechnology, Santa Cruz, CA) at 1:200 on tissue sections subjected to antigen retrieval for 30 min at 98°C. The extent of virus spread was quantified by capturing digital images of whole-lung sections (two dimensional) stained for viral antigen using an Aperio ScanScope XT Slide Scanner (Aperio Technologies, Vista, CA) then manually outlining defined fields. Alveolar areas containing virus antigen-positive pneumocytes were highlighted in red (defined as ‘active’ infection), whereas lesioned areas containing minimal or no virus antigen-positive debris were highlighted in green (defined as ‘inactive’ infection). Lesions containing a mix of virus antigen-positive and antigen-negative pneumocytes were highlighted in orange (defined as ‘mixed’ infection). The percentage of each defined lung field was calculated using the Aperio ImageScope software. Pulmonary lesions in HE-stained histologic sections were assigned scores based on their severity and extent. Representative images from each group and quantitative analyses (five animals/group) of viral spread and lung pathology during IAV infection are shown in Figure 4, and provided in Source data 1.

CD8⁺ T cell model

Request a detailed protocol

To examine the contribution of CD8⁺ T cells to the biphasic viral load decay, we expanded the density-dependent (DD) model (see Appendix 2) to include two mechanisms of infected cell clearance (non-specific clearance ( $δ$ ) and CD8⁺ T cell-mediated clearance ( $δ_{E} (I_{2}, E)$ )) and two CD8⁺ T cell populations: effector ( $E$ , denoted CD8_E) and memory ( $E_{M}$ , denoted CD8_M) CD8⁺ T cells. The model is given by Equation (1)-(6).

\frac{d T}{d t} = - β T V

\frac{d I_{1}}{d t} = β T V - k I_{1}

\frac{d I_{2}}{d t} = k I_{1} - δ I_{2} - δ_{E} (I_{2}, E) I_{2}

\frac{d V}{d t} = p I_{2} - c V

\frac{d E}{d t} = ξ (E) I_{2} + η E I_{2} (t - τ_{E}) - d_{E} E

\frac{d E_{M}}{d t} = ζ E (t - τ_{M})

In this model, target cells become infected with virus at rate $β V$ per day. Once infected, cells enter an eclipse phase ( $I_{1}$ ) before transitioning at rate $k$ per day to a productively-infected state ( $I_{2}$ ). These infected cells produce virus at rate $p$ TCID₅₀ per infected cell per day, and virus is cleared at rate $c$ per day. Virus-producing infected cells ( $I_{2}$ ) are cleared by non-specific mechanisms (e.g., apoptosis and/or innate immune responses) at a constant rate $δ$ per day. Innate immune responses were excluded from the model because the viral load data is linear during the time where they act (2–7 d pi) and, thus, additional equations cannot improve the fit. Productively infected cells are cleared by CD8_E at rate $δ_{E} (I_{2}, E) = δ_{E} E / (K_{δ_{E}} + I_{2})$ per day, where the rate of infected cell clearance is $δ_{E} / K_{δ_{E}}$ per CD8_E per day and $K_{δ_{E}}$ is the half-saturation constant. The CD8_E-mediated clearance rate ( $δ_{E} (I_{2}, E)$ ) is dependent on the density of infected cells and is similar to the infected cell clearance term in the DD model (see Equation (A11); Smith et al., 2018). Similar density-dependent forms have also been used in models that describe the CD8⁺ T cell response to other virus infections (De Boer and Perelson, 1998; Li and Handel, 2014; Gadhamsetty et al., 2014). Models that exclude this density-dependence were examined, but these models resulted in a statistically poor fit to the data as defined by the Akaike Information Criteria (AIC) (Supplementary file 2). This is due in part to the increase of CD8⁺ T cells from 3-5 d pi (Figure 1). We did examine a model that excluded these dynamics and included CD8⁺ T cell expansion as a density-dependent function (e.g., $η E I_{2} (t - τ_{E}) / (K_{I} + I_{2})$ ) while keeping a linear rate of CD8-mediated infected cell clearance ( $δ_{E} E I_{2}$ ) (see Appendix 1). This model was adapted from other published models (Conway and Perelson, 2015; Bonhoeffer et al., 2000; Baral et al., 2019) and can produce similar dynamics from 6-10 d pi, but it was not statistically supported by our data (Supplementary file 1). We further tested the model in Baral et al., 2019, but this model could not fit our data and lacked statistical support (Supplementary file 1).

The model assumes that the initial CD8_E influx in the lung is proportional to infected cells at rate $ξ (E) = ξ / (K_{E} + E)$ CD8_E per cell per day, which is down-regulated by the CD8_E already present in the lung. The associated half-saturation constant is $K_{E}$ . Similar terms for CD8_E regulation have been used in modeling HIV infections (De Boer and Perelson, 1998; Müller et al., 2001) and in models that examined CD8⁺ T cell proliferation mechanisms (De Boer and Perelson, 1995). We also examined whether their influx is proportional to infected cells at a delay (i.e., $ξ (E) I (t - τ_{I})$ ). While this modification better captured the initial increase in CD8⁺ T cells, the additional parameter was not supported. In our model, CD8_E expansion occurs at rate $η$ per infected cell per day with time delay $τ_{E}$ . This term accounts for local CD8_E proliferation in the lung (McGill and Legge, 2009; Lawrence et al., 2005) and migration of CD8_E from secondary lymphoid organs (Zhang and Bevan, 2011; Bedoui and Gebhardt, 2011; Wu et al., 2011; Kang et al., 2011). The delay may signify the time it takes CD8_E to become activated by antigen presenting cells, differentiate, proliferate, and/or migrate to the infection site. The lung CD8_E population declines due to cell death and/or emigration at rate $d_{E}$ per day. These cells transition to CD8_M ( $E_{M}$ ) at rate ζ CD8_M per CD8_E per day after $τ_{M}$ days. The model schematic and fit to the viral load, total, effector, and memory CD8⁺ T cell data are in Figure 1. All code is provided in Source code 1.

Inflammation model

Request a detailed protocol

To estimate the alveolar and interstitial inflammation without modeling other cell classes (e.g., macrophages and neutrophils), we assumed that inflammation in the lung ( $L_{I}$ ) was proportional to the infected cells ( $I_{1}$ and $I_{2}$ ) according to the equation,

\frac{d L_{I}}{d t} = α_{1} I_{1} + α_{2} I_{2},

where $α_{1, 2}$ has units of score/cell/d and defines the inflammation score contribution from each infected cell class. A decay term was excluded because inflammation does not resolve on the timescale of our data.

Parameter estimation

Request a detailed protocol

Given a parameter set $θ$ , the cost $C (θ)$ was minimized across parameter ranges using an Adaptive Simulated Annealing (ASA) global optimization algorithm (Smith et al., 2018) to compare experimental and predicted values of virus ( $V$ ; log₁₀ TCID₅₀/lung) and log₁₀ total CD8⁺ T cells/lung ( $\hat{E} = E + E_{M} + {\hat{E}}_{0}$ , where ${\hat{E}}_{0}$ is the initial number of CD8⁺ T cells at 0 d pi), or log₁₀ TCID₅₀/lung virus ( $V$ ), log₁₀ effector CD8⁺ T cells/lung ( $\hat{E} = E + {\hat{E}}_{0}$ ), and log₁₀ memory CD8⁺ T cells/lung ( ${\hat{E}}_{M} = E_{M} + {\hat{E}}_{M_{0}}$ ). The cost function is defined by,

\begin{array}{ll} C (θ) = \sum_{i, j} (V (θ, t_{i}) - v_{i, j})^{2} + \sum_{i, j} ({\hat{E}}_{M} (θ, t_{i}) - m_{i, j})^{2} + s_{E} [\sum_{i, j} (\hat{E} (θ, t_{i}) - e_{i, j})^{2} + \\ \sum_{i} \sqrt{γ_{i}} {(\frac{\hat{E} (θ, t_{i + 1}) - \hat{E} (θ, t_{i - 1})}{t_{i + 1} - t_{i - 1}} - \frac{1}{γ_{i}} \sum_{j} \frac{e_{i + 1, j} - e_{i - 1, j}}{t_{i + 1} - t_{i - 1}})}^{2}], \end{array}

where $(t_{i}, v_{i, j})$ is the viral load data, $(t_{i}, e_{i, j})$ is the total or effector CD8⁺ T cell data, $(t_{i}, m_{i, j})$ is the memory CD8⁺ T cell data, and $V (θ, t_{i})$ , $\hat{E} (θ, t_{i})$ , and $\hat{E_{M}} (θ, t_{i})$ are the corresponding model predictions. Here, $s_{E} = (v_{max} - v_{min}) / (e_{max} - e_{min})$ is a scaling factor, and $γ_{i} = J_{i + 1} J_{i - 1}$ where $J_{i}$ is the number of observations at time $t_{i}$ . Errors of the log₁₀ data were assumed to be normally distributed. To explore and visualize the regions of parameters consistent with the model, we fit Equation (1)-(6) to 2000 bootstrap replicates of the data. If the fit was within $χ^{2} = 0.05$ of the best-fit and the CD8 derivative was not a statistical outlier as determined by the function isoutlier, then the bootstrap was considered successful (Smith et al., 2011a; Smith et al., 2013; Smith et al., 2018). For each best-fit estimate, we provide 95% confidence intervals (CI) obtained from the bootstrap replicates (Table 1). Calculations were performed either in MATLAB using a custom built ASA algorithm (Smith et al., 2018) or in Python using the simanneal package (Perry, 2018) followed by a L-BFGS-B (Byrd et al., 1995; Zhu et al., 1997) deterministic minimization through SciPy’s minimize function. MATLAB ode15s and dde23 or SciPy integrate.ode using lsoda and PyDDE (Cairns, 2008) were used as the ODE and DDE solvers.

Estimated parameters in the CD8⁺ T cell model included the rates of virus infection ( $β$ ), virus production ( $p$ ), virus clearance ( $c$ ), eclipse phase transition ( $k$ ), non-specific infected cell clearance ( $δ$ ), CD8_E-mediated infected cell clearance ( $δ_{E}$ ), half-saturation constants ( $K_{δ_{E}}$ and $K_{E}$ ), CD8_E infiltration ( $ξ$ ), CD8_E expansion ( $η$ ), delay in CD8_E expansion ( $τ_{E}$ ), CD8_Eclearance ( $d_{E}$ ), CD8_M generation ( $ζ$ ), delay in CD8_M generation ( $τ_{M}$ ), and the baseline number of CD8⁺ T cells ( ${\hat{E}}_{0}$ , ${\hat{M}}_{0}$ ). Bounds were placed on the parameters to constrain them to physically realistic values. Because biological estimates are not available for all parameters, ranges were set reasonably large based on preliminary results and previous estimates (Smith et al., 2018). The rate of infection ( $β$ ) was allowed to vary between $10^{- 6} - 10^{- 1} {TCID}_{50}^{- 1} d^{- 1}$ , and the rate of virus production ( $p$ ) between $10^{- 1} - 10^{3} {TCID}_{50} {cell}^{- 1} d^{- 1}$ . Bounds for the virus clearance rate ( $c$ ) were $1 d^{- 1}$ ( $t_{1 / 2} = 16.7 h$ ) and $10^{3} d^{- 1}$ ( $t_{1 / 2} = 1 \min$ ). To ensure biological feasibility, the lower and upper bounds for the eclipse phase transition rate ( $k$ ) were 4-6 d^-1 as done previously (Smith et al., 2018).

The rate of non-specific infected cell clearance (δ) was given limits of 0.05–1 d^-1. The CD8_E-mediated infected cell clearance rate ( $δ_{E}$ ) varied between $0.01 - 2 cells {CD8}_{E}^{- 1} d^{- 1}$ , and the associated half-saturation constant ( $K_{δ_{E}}$ ) was bounded between 10¹–10⁶ cells. The upper bound of $δ_{E}$ was chosen to maintain the convergence of $δ$ to nonzero values and consistency with prior results where $δ$ approximates the slope of the first decay phase (Smith et al., 2010). Bounds for the rate of CD8_E infiltration (ξ) were $10^{2} - 10^{6} {C D 8}_{E}^{2} {c e l l}^{- 1} d^{- 1}$ , and bounds for the half-saturation constant ( $K_{E}$ ) were 10³-10⁷ CD8_E. The CD8_E expansion rate ( $η$ ) varied between 10^–8–10^–6 cell^–1 d^–1, and the delay in CD8_E expansion ( $τ_{E}$ ) between 2–6 d. The rate of CD8_E clearance ( $d_{E}$ ) had limits of 0.05–2 d^–1. The rate of CD8_M generation ( $ζ$ ) varied between $0.01 - 1 C D 8_{M} C D 8_{E}^{- 1} d^{- 1}$ , and the delay in CD8_M generation ( $τ_{M}$ ) varied between 3–4 d (total CD8⁺ T cell fit) or 2–6 d pi (effector and memory CD8⁺ T cell fit). Larger bounds were examined for this parameter; however, the parameter was non-identifiable in the total CD8⁺ T cell fit and a small range was required for convergence. Bounds for the baseline number of CD8⁺ T cells ( ${\hat{E}}_{0}$ , ${\hat{E}}_{M_{0}}$ ) were set to the upper and lower values of the data at 0 d pi ( $3.0 \times 10^{5} - 5.3 \times 10^{5} CD8$ (total), $2.5 \times 10^{2} - 1.4 \times 10^{3} {CD8}_{E}$ (effector), $3.4 \times 10^{2} - 8.0 \times 10^{2} {CD8}_{M}$ (memory)).

The initial number of target cells ( $T (0)$ ) was set to 10⁷ cells (Smith et al., 2011a). The initial number of infected cells $I_{1} (0)$ was set to 75 cells to reflect an initial dose of 75 TCID₅₀ (Smith et al., 2018). We previously found that estimating $I_{1} (0)$ , fixing $V (0) = 75 {TCID}_{50}$ , or estimating $V (0)$ did not improve the fit and could not be statistically justified (Smith et al., 2018). The initial number of productively infected cells ( $I_{2} (0)$ ), the initial free virus ( $V (0)$ ), and the initial number of CD8_E ( $E (0)$ ) and CD8_M ( $E_{M} (0)$ ) were set to 0. All other parameter estimations were done as described in the text.

Linear regression

Request a detailed protocol

The function polyfit in MATLAB was used to perform linear regression on the percent active lesioned area, the percent inactive lesioned area, and the CD8⁺ T cells during the expansion phase (5–8 d pi) and the contraction phase (9–10 d pi). Linear fits are shown in Figure 4—figure supplement 1.

Cumulative area under the curve

Request a detailed protocol

The function cumtrapz in MATLAB was used to estimate the cumulative area under the curve (CAUC) of the infected cells ( $I_{2}$ ) for the best-fit model solution.

Appendix 1

Alternate CD8⁺ T cell Models

We examined alternate formulations of the CD8⁺ T cell model to further investigate the density-dependence in the CD8⁺ T cell response. Rather than assuming that CD8_E-mediated clearance of infected cells is dependent on their density, the model in Equation (A1)-(A6) assumes that the rate of CD8_E expansion is dependent on the density of infected cells. Similar models have been used to study CD8⁺ T cell responses during HIV infection (Conway and Perelson, 2015; Bonhoeffer et al., 2000) and other viral infections (Baral et al., 2019). The differences between this alternate model and the CD8⁺ T cell model (Equation (1)-(6)) are in bold.

\frac{d T}{d t} = - β T V

\frac{d I_{1}}{d t} = β T V - k I_{1}

\frac{d I_{2}}{d t} = k I_{1} - δ I_{2} - δ_{𝐄𝐚} {𝐄𝐈}_{𝟐}

\frac{d V}{d t} = p I_{2} - c V

\frac{d E}{d t} = \frac{η_{𝐄𝐚}}{𝐊_{𝐄𝐚} + 𝐈_{𝟐}} I_{2} (t - τ_{E}) - d_{E} E

\frac{d E_{M}}{d t} = ζ E (t - τ_{M})

When a linear CD8_E-mediated infected cell clearance rate is included, the CD8_E dynamics between 3–5 d pi cannot be replicated and returned a statistically worse fit via AIC (Supplementary file 1). However, because these cells may not have effector functions and contribute to infected cell clearance at this time (see Discussion), we excluded these data and the term $ξ I_{2} / (K_{E} + E)$ when fitting Equation (A1)-(A6) to the viral load and total CD8⁺ T cell data (Appendix 1—figure 1). The CD8 dynamics are similar to those generated by CD8⁺ T cell model. However, the alternate model underestimates the data at day seven and the sharp decline between 7–8 d. While it is inappropriate to directly compare these models due to the varying number of data points, we assessed their goodness of fit with and without contributions from the data at 3–5 d pi (Supplementary file 1). Regardless of the data inclusion or exclusion, the alternate model had a higher AIC in all contexts and, thus, is not statistically justifiable. Similarly, we examined the model in Baral et al., 2019 (Equation (A7)-(A8)), which also could not replicate our data (Appendix 1—figure 1).

\frac{d I}{d t} = κ_{b} I (1 - \frac{I}{I_{m a x}}) - δ_{E b} I E

\frac{d E}{d t} = η_{b} \frac{I E}{K_{E b} + I} - d_{E b} \frac{I E}{K_{D b} + I}

Appendix 1—figure 1

Download asset Open asset

Fit of alternate CD8⁺ T cell models.

Fit of the CD8⁺ T cell model (Equation (1)-(6); solid black line) compared to the fit of two alternate CD8⁺ T cell models (Equation (A1)-(A6) (solid blue line) and Equation (A7)-(A8) (dashed blue line)) to virus and CD8⁺ T cells (excluding 3–5 d pi; white squares) from the lungs of mice infected with 75 TCID₅₀ PR8 (10 mice per time point). Resulting parameter values were $δ_{E a} = 4.02 \times 10^{- 6} {C D 8}_{E}^{- 1} d^{- 1}$ , $η_{E a} = 3.12 \times 10^{5} {CD8}_{E} / d$ , and $K_{E a} = 9.53 \times 10^{5} i n f e c t e d c e l l s$ , and $κ_{b} = 3.33 d^{- 1}$ , $I_{m a x} = 1.0 \times 10^{6} i n f e c t e d c e l l s$ , $δ_{E b} = 1.53 \times 10^{- 5} {CD8}_{E}^{- 1} d^{- 1}$ , $η_{b} = 1.49 d^{- 1}$ , $K_{E b} = 9.5 i n f e c t e d c e l l s$ , $d_{E b} = 0.95 d^{- 1}$ , and $K_{D b} = 1.90 \times 10^{6} i n f e c t e d c e l l s$ . All other parameters are in Table 1, and the AICs are in Supplementary file 1. Data are shown as mean ± standard deviation.

Appendix 2

Comparison of the Density-Dependent and CD8⁺ T cell Models

We previously developed a density-dependent (DD) viral kinetic model, which describes the biphasic decline of viral loads without inclusion of specific host responses (Smith et al., 2018). This model tracks four populations: susceptible epithelial (‘target’) cells ( $T$ ), two classes of infected cells ( $I_{1}$ and $I_{2}$ ), and virus ( $V$ ) (Smith et al., 2018).

\frac{d T}{d t} = - β T V

\frac{d I_{1}}{d t} = β T V - k I_{1}

\frac{d I_{2}}{d t} = k I_{1} - δ_{d} (I_{2}) I_{2}

\frac{d V}{d t} = p I_{2} - c V

Briefly, in the DD model, virus-producing infected cells ( $I_{2}$ ) are cleared according to the function $δ_{d} (I_{2}) = δ_{d} / (K_{δ} + I_{2})$ , where $δ_{d} / K_{δ}$ is the maximum per day rate of infected cell clearance and $K_{δ}$ is the half-saturation constant (Appendix 2—figure 1). All other terms are common to the CD8 T⁺ cell model (Equation (1)-(6)). The DD model provides a close fit to the viral load data in Figure 1B and replicates the biphasic viral load decline while excluding the dynamics of specific immune responses (Smith et al., 2018). Unsurprisingly, the CD8⁺ T cell model is also capable of reproducing the biphasic viral load decay (Figure 1B and Appendix 2—figure 1). In that model, infected cell clearance is split into terms for non-specific clearance ( $δ$ ) and CD8_E-mediated clearance ( $δ_{E} (I_{2}, E) = δ_{E} E / (K_{δ} + I_{2})$ ) (Appendix 2—figure 1).

Because the CD8⁺ T cell model is more mechanistic than the DD model, most of the correlations between the parameters common to both models (i.e., the rates of virus infectivity ( $β$ ), virus production ( $p$ ), and virus clearance ( $c$ )) were reduced (Appendix 2—figure 1A). In addition, the correlations between the infected cell clearance parameters ( $δ_{d}$ and $K_{δ}$ or $δ_{E}$ and $K_{δ_{E}}$ ) and between the rate of virus infectivity ( $β$ ) and their ratios ( $δ_{d} / K_{δ}$ or $δ_{E} / K_{δ_{E}}$ ) were abolished (Figure 2—figure supplement 1–2). There was a negative correlation between the infected cell clearance parameters ( $δ$ and $δ_{E}$ ; Figure 2B), which may reflect the connection between the efficacy of early immune mechanisms and the CD8⁺ T cell response. This result is in line with experimental evidence that the innate immune responses modulate the activation of adaptive immunity (Iwasaki and Medzhitov, 2015; Luster, 2002; Iwasaki and Medzhitov, 2010; Le Bon and Tough, 2002; Jain and Pasare, 2017).

The differences in model structure between the two models yielded changes in parameter sensitivity and model behavior during the rapid viral clearance phase (Appendix 2—figure 1). In the DD model, the most sensitive parameter is the infected cell clearance, $δ_{d}$ (Appendix 2—figure 1). A 50% decrease in this parameter resulted in a ∼7 d delay in viral resolution (Appendix 2—figure 1; Smith et al., 2018). In the CD8⁺ T cell model, however, viral resolution is delayed by <1 d if the CD8_E-mediated infected cell clearance parameter ( $δ_{E}$ ) is reduced by 50% (Appendix 2—figure 1, Appendix 3—figure 1). The rates of CD8_E expansion ( $η$ ) and decay ( $d_{E}$ ) are sensitive and, thus, significantly influence the viral resolution kinetics (Figures 1–2). A 50% decrease in $η$ results in a ∼6 d delay in recovery (Appendix 2—figure 1, Appendix 3—figure 1) whereas a 48% decrease in $η$ prolongs the infection by ∼30 d (Figure 3D–E). This bifurcation in recovery time is a unique feature of the CD8⁺ T cell model.

Appendix 2—figure 1

Download asset Open asset

Parameter behavior of the density-dependent model and the CD8⁺ T cell model.

(A) Comparison of parameters that were common between the density-dependent model (gray, Equation (A9)-(A12)) and the CD8⁺ T cell model (purple, Equation (1)-(6)). Correlations were evident between parameters relating to the rates of virus infectivity ( $β$ ), virus production ( $p$ ), and virus clearance ( $c$ ). However, the strength of the correlation was significantly reduced in the CD8⁺ T cell model. The eclipse phase parameter ( $k$ ) was not well-defined in either model. (B) In the density-dependent model (gray), the viral kinetics and the infection duration were sensitive to small changes in the infected cell clearance parameter ( $δ_{d}$ ). This parameter was well-defined with a narrow 95% CI. In the CD8⁺ T cell model (purple), changing the CD8_E-mediated infected cell clearance parameter ( $δ_{E}$ ) had little impact on viral kinetics or CD8⁺ T cell kinetics. However, these kinetics were most sensitive to changes in the rate of CD8_E expansion ( $η$ ), which was well-defined with a narrow 95% CI.

Appendix 3

Regulation of the CD8⁺ T cell Response

To further understand the regulation of the CD8⁺ T cell response, we examined the 2-D parameter ensembles (Figure 2A–C, Figure 2—figure supplements 1–2) and the results from the sensitivity analysis (Appendix 3—figure 1–2). Overall, few parameters were correlated. There was an expected, although small, positive correlation between the rate of CD8_E infiltration ( $ξ$ ) and the associated half-saturation constant ( $K_{E}$ ) (Figure 2—figure supplement 1–2), which represents the coordination between CD8_E recruitment and the processes that prevent an overabundance of these cells. Likewise, a negative correlation was detected between the rate of initial CD8_E influx ( $ξ$ ) and the initial number of CD8⁺ T cells ( ${\hat{E}}_{0}$ ) (Figure 2—figure supplements 1–2). The influx rate ( $ξ$ ) was also positively correlated with the delay in CD8_E expansion ( $τ_{E}$ ) (Figure 2—figure supplements 1–2). The rates of CD8_E expansion ( $η$ ) and decay ( $d_{E}$ ) are correlated (Figure 2C), indicating a balance between these two processes. This correlation was expected and reflects the coordination of mechanisms that regulate CD8⁺ T cell numbers, which may be necessary to limit excessive immunopathology while still resolving the infection (Moskophidis and Kioussis, 1998; Duan and Thomas, 2016; La Gruta et al., 2007). Further, because of this correlation and the sensitivity of $η$ (Appendix 3—figure 1), the CD8⁺ T cell kinetics are sensitive to changes in $d_{E}$ (Appendix 3—figure 1). However, increasing the decay rate had less impact on the viral load kinetics, comparatively. Because $d_{E}$ is correlated with both $η$ and the rate of CD8_M generation ( $ζ$ ) (Figure 2C), it naturally follows that $η$ and $ζ$ are correlated (Figure 2—figure supplements 1–2). Changing the rates of virus infectivity ( $β$ ), production ( $p$ ), or clearance ( $c$ ) had little effect on viral load or CD8⁺ T cell kinetics (Appendix 3—figure 1).

Appendix 3—figure 1

Download asset Open asset

Sensitivity of the CD8⁺ T cell model.

Solutions of the CD8⁺ T cell model (Equation (1)-(6)) with the indicated parameter ( $β$ , $k$ , $p$ , $c$ , $δ$ , $δ_{E}$ , $K_{δ_{E}}$ , $ξ$ , $K_{E}$ , or $η$ ) increased (red) or decreased (blue) 50% from the best-fit value (Table 1). CD8_E denotes effector CD8⁺ T cells.

Appendix 3—figure 2

Download asset Open asset

Sensitivity of the CD8⁺ T cell model.

Solutions of the CD8⁺ T cell model (Equation (1)-(6)) with the indicated parameter ( $τ_{E}$ , $d_{E}$ , $ζ$ , $τ_{M}$ , or ${\hat{E}}_{0}$ ) increased (red) or decreased (blue) 50% from the best-fit value (Table 1). CD8_E and CD8_M denote effector and memory CD8⁺ T cells, respectively.

Data availability

All data generated or analyzed during this study are included in the manuscript and supporting files. Source data files have been provided.

References

Preprint
(2017) Mathematical analysis of a mouse experiment suggests little role for resource depletion in controlling influenza infection within host
arXiv.

https://arxiv.org/abs/1705.02565
- Google Scholar
1. Ai T
2. Yang Z
3. Hou H
4. Zhan C
5. Chen C
6. Lv W
7. Tao Q
8. Sun Z
9. Xia L
(2020) Correlation of chest CT and RT-PCR testing for coronavirus disease 2019 (COVID-19) in China: a report of 1014 cases
Radiology 296:E32–E40.
https://doi.org/10.1148/radiol.2020200642
- PubMed
- Google Scholar
(2009) Swine-origin influenza A (H1N1) viral infection: radiographic and CT findings
American Journal of Roentgenology 193:1494–1499.
https://doi.org/10.2214/AJR.09.3625
- PubMed
- Google Scholar
1. Anderson KG
2. Sung H
3. Skon CN
4. Lefrancois L
5. Deisinger A
6. Vezys V
7. Masopust D
(2012) Cutting edge: intravascular staining redefines lung CD8 T cell responses
The Journal of Immunology 189:2702–2706.
https://doi.org/10.4049/jimmunol.1201682
- PubMed
- Google Scholar
(2006) Kinetics of influenza A virus infection in humans
Journal of Virology 80:7590–7599.
https://doi.org/10.1128/JVI.01623-05
- PubMed
- Google Scholar
1. Baloxavir Marboxil Investigators Group
2. Hayden FG
3. Sugaya N
4. Hirotsu N
5. Lee N
6. de Jong MD
7. Hurt AC
8. Ishida T
9. Sekino H
10. Yamada K
11. Portsmouth S
12. Kawaguchi K
13. Shishido T
14. Arai M
15. Tsuchiya K
16. Uehara T
17. Watanabe A
(2018) Baloxavir marboxil for uncomplicated influenza in adults and adolescents
New England Journal of Medicine 379:913–923.
https://doi.org/10.1056/NEJMoa1716197
- PubMed
- Google Scholar
1. Baral S
2. Antia R
3. Dixit NM
(2019) A dynamical motif comprising the interactions between antigens and CD8 T cells may underlie the outcomes of viral infections
PNAS 116:17393–17398.
https://doi.org/10.1073/pnas.1902178116
- PubMed
- Google Scholar
(2010) Clinical aspects of pandemic 2009 influenza A (H1N1) virus infection
The New England Journal of Medicine 362:1708–1719.
https://doi.org/10.1056/NEJMra1000449
- PubMed
- Google Scholar
1. Beauchemin CA
2. Handel A
(2011) A review of mathematical models of influenza A infections within a host or cell culture: lessons learned and challenges ahead
BMC Public Health 11 Suppl 1:S7.
https://doi.org/10.1186/1471-2458-11-S1-S7
- PubMed
- Google Scholar
1. Bedoui S
2. Gebhardt T
(2011) Interaction between dendritic cells and T cells during peripheral virus infections: a role for antigen presentation beyond lymphoid organs?
Current Opinion in Immunology 23:124–130.
https://doi.org/10.1016/j.coi.2010.11.001
- PubMed
- Google Scholar
1. Bender BS
2. Croghan T
3. Zhang L
4. Small PA
(1992) Transgenic mice lacking class I major histocompatibility complex-restricted T cells have delayed viral clearance and increased mortality after influenza virus challenge
Journal of Experimental Medicine 175:1143–1145.
https://doi.org/10.1084/jem.175.4.1143
- PubMed
- Google Scholar
(2015) Modeling influenza virus infection: a roadmap for influenza research
Viruses 7:5274–5304.
https://doi.org/10.3390/v7102875
- PubMed
- Google Scholar
(2000) Risks and benefits of structured antiretroviral drug therapy interruptions in HIV-1 infection
Aids 14:2313–2322.
https://doi.org/10.1097/00002030-200010200-00012
- PubMed
- Google Scholar
1. Boon AC
2. Finkelstein D
3. Zheng M
4. Liao G
5. Allard J
6. Klumpp K
7. Webster R
8. Peltz G
9. Webby RJ
(2011) H5N1 influenza virus pathogenesis in genetically diverse mice is mediated at the level of viral load
mBio 2:e00171.
https://doi.org/10.1128/mBio.00171-11
- PubMed
- Google Scholar
1. Bouvier NM
2. Lowen AC
(2010) Animal models for influenza virus pathogenesis and transmission
Viruses 2:1530–1563.
https://doi.org/10.3390/v20801530
- PubMed
- Google Scholar
1. Boyd DF
2. Allen EK
3. Randolph AG
4. Guo XJ
5. Weng Y
6. Sanders CJ
7. Bajracharya R
8. Lee NK
9. Guy CS
10. Vogel P
11. Guan W
12. Li Y
13. Liu X
14. Novak T
15. Newhams MM
16. Fabrizio TP
17. Wohlgemuth N
18. Mourani PM
19. Wight TN
20. Schultz-Cherry S
21. Cormier SA
22. Shaw-Saliba K
23. Pekosz A
24. Rothman RE
25. Chen KF
26. Yang Z
27. Webby RJ
28. Zhong N
29. Crawford JC
30. Thomas PG
31. PALISI Pediatric Intensive Care Influenza (PICFLU) Investigators
(2020) Exuberant fibroblast activity compromises lung function via ADAMTS4
Nature 587:466–471.
https://doi.org/10.1038/s41586-020-2877-5
- PubMed
- Google Scholar
1. Byrd RH
2. Lu P
3. Nocedal J
4. Zhu C
(1995) A limited memory algorithm for bound constrained optimization
SIAM Journal on Scientific Computing 16:1190–1208.
https://doi.org/10.1137/0916069
- Google Scholar
Software
1. Cairns BJ
(2008)
Pydde, version 0.2.2

Pydde.
1. Cao P
2. Wang Z
3. Yan AW
4. McVernon J
5. Xu J
6. Heffernan JM
7. Kedzierska K
8. McCaw JM
(2016) On the role of CD8⁺ T cells in determining recovery time from influenza virus infection
Frontiers in Immunology 7:611.
https://doi.org/10.3389/fimmu.2016.00611
- PubMed
- Google Scholar
(2009) Visualizing CTL/melanoma cell interactions: multiple hits must be delivered for tumour cell annihilation
Journal of Cellular and Molecular Medicine 13:3834–3846.
https://doi.org/10.1111/j.1582-4934.2008.00586.x
- PubMed
- Google Scholar
1. Carrat F
2. Vergu E
3. Ferguson NM
4. Lemaitre M
5. Cauchemez S
6. Leach S
7. Valleron AJ
(2008) Time lines of infection and disease in human influenza: a review of volunteer challenge studies
American Journal of Epidemiology 167:775–785.
https://doi.org/10.1093/aje/kwm375
- PubMed
- Google Scholar
1. Channappanavar R
2. Fehr AR
3. Vijay R
4. Mack M
5. Zhao J
6. Meyerholz DK
7. Perlman S
(2016) Dysregulated type I interferon and inflammatory Monocyte-Macrophage responses cause lethal pneumonia in SARS-CoV-Infected mice
Cell Host & Microbe 19:181–193.
https://doi.org/10.1016/j.chom.2016.01.007
- PubMed
- Google Scholar
1. Chen X
2. Liu S
3. Goraya MU
4. Maarouf M
5. Huang S
6. Chen JL
(2018) Host immune response to influenza A virus infection
Frontiers in Immunology 9:320.
https://doi.org/10.3389/fimmu.2018.00320
- PubMed
- Google Scholar
1. Conway JM
2. Perelson AS
(2015) Post-treatment control of HIV infection
PNAS 112:5467–5472.
https://doi.org/10.1073/pnas.1419162112
- PubMed
- Google Scholar
1. Cross EW
2. Blain TJ
3. Mathew D
4. Kedl RM
(2019) Anti-CD8 monoclonal antibody-mediated depletion alters the phenotype and behavior of surviving CD8+ T cells
PLOS ONE 14:e0211446.
https://doi.org/10.1371/journal.pone.0211446
- PubMed
- Google Scholar
1. De Boer RJ
2. Perelson AS
(1995) Towards a general function describing T cell proliferation
Journal of Theoretical Biology 175:567–576.
https://doi.org/10.1006/jtbi.1995.0165
- PubMed
- Google Scholar
1. De Boer RJ
2. Perelson AS
(1998) Target cell limited and immune control models of HIV infection: a comparison
Journal of Theoretical Biology 190:201–214.
https://doi.org/10.1006/jtbi.1997.0548
- PubMed
- Google Scholar
1. de Jong MD
2. Simmons CP
3. Thanh TT
4. Hien VM
5. Smith GJ
6. Chau TN
7. Hoang DM
8. Chau NV
9. Khanh TH
10. Dong VC
11. Qui PT
12. Cam BV
13. Ha doQ
14. Guan Y
15. Peiris JS
16. Chinh NT
17. Hien TT
18. Farrar J
(2006) Fatal outcome of human influenza A (H5N1) is associated with high viral load and hypercytokinemia
Nature Medicine 12:1203–1207.
https://doi.org/10.1038/nm1477
- PubMed
- Google Scholar
(2010) Intravital imaging reveals distinct dynamics for natural killer and CD8(+) T cells during tumor regression
Immunity 33:632–644.
https://doi.org/10.1016/j.immuni.2010.09.016
- PubMed
- Google Scholar
(2013) Assessing mathematical models of influenza infections using features of the immune response
PLOS ONE 8:e57088.
https://doi.org/10.1371/journal.pone.0057088
- PubMed
- Google Scholar
(2015) Diverse heterologous primary infections radically alter immunodominance hierarchies and clinical outcomes following H7N9 influenza challenge in mice
PLOS Pathogens 11:e1004642.
https://doi.org/10.1371/journal.ppat.1004642
- PubMed
- Google Scholar
1. Duan S
2. Thomas PG
(2016) Balancing immune protection and immune pathology by CD8(+) T-Cell responses to influenza infection
Frontiers in Immunology 7:25.
https://doi.org/10.3389/fimmu.2016.00025
- PubMed
- Google Scholar
(2014) The resolution of inflammation: a mathematical model of neutrophil and macrophage interactions
Bulletin of Mathematical Biology 76:1953–1980.
https://doi.org/10.1007/s11538-014-9987-x
- PubMed
- Google Scholar
1. Eccles R
(2005) Understanding the symptoms of the common cold and influenza
The Lancet Infectious Diseases 5:718–725.
https://doi.org/10.1016/S1473-3099(05)70270-X
- PubMed
- Google Scholar
(1991) Clearance of influenza virus respiratory infection in mice lacking class I major histocompatibility complex-restricted CD8+ T cells
Journal of Experimental Medicine 174:875–880.
https://doi.org/10.1084/jem.174.4.875
- PubMed
- Google Scholar
1. Fang M
2. Sigal LJ
(2005) Antibodies and CD8+ T cells are complementary and essential for natural resistance to a highly lethal cytopathic virus
The Journal of Immunology 175:6829–6836.
https://doi.org/10.4049/jimmunol.175.10.6829
- PubMed
- Google Scholar
1. Feikin DR
2. Fu W
3. Park DE
4. Shi Q
5. Higdon MM
6. Baggett HC
7. Brooks WA
8. Deloria Knoll M
9. Hammitt LL
10. Howie SRC
11. Kotloff KL
12. Levine OS
13. Madhi SA
14. Scott JAG
15. Thea DM
16. Adrian PV
17. Antonio M
18. Awori JO
19. Baillie VL
20. DeLuca AN
21. Driscoll AJ
22. Ebruke BE
23. Goswami D
24. Karron RA
25. Li M
26. Morpeth SC
27. Mwaba J
28. Mwansa J
29. Prosperi C
30. Sawatwong P
31. Sow SO
32. Tapia MD
33. Whistler T
34. Zaman K
35. Zeger SL
36. O' Brien KL
37. Murdoch DR
38. PERCH Study Group
(2017) Is higher viral load in the upper respiratory tract associated with severe pneumonia? findings from the PERCH study
Clinical Infectious Diseases 64:S337–S346.
https://doi.org/10.1093/cid/cix148
- PubMed
- Google Scholar
(2014) A general functional response of cytotoxic T lymphocyte-mediated killing of target cells
Biophysical Journal 106:1780–1791.
https://doi.org/10.1016/j.bpj.2014.01.048
- PubMed
- Google Scholar
1. Gallagher M
2. Brooke C
3. Ke R
4. Koelle K
(2018) Causes and consequences of spatial Within-Host viral spread
Viruses 10:627.
https://doi.org/10.3390/v10110627
- Google Scholar
(2011) Killing of targets by CD8 T cells in the mouse spleen follows the law of mass action
PLOS ONE 6:e15959.
https://doi.org/10.1371/journal.pone.0015959
- PubMed
- Google Scholar
1. Gao R
2. Bhatnagar J
3. Blau DM
4. Greer P
5. Rollin DC
6. Denison AM
7. Deleon-Carnes M
8. Shieh WJ
9. Sambhara S
10. Tumpey TM
11. Patel M
12. Liu L
13. Paddock C
14. Drew C
15. Shu Y
16. Katz JM
17. Zaki SR
(2013) Cytokine and chemokine profiles in lung tissues from fatal cases of 2009 pandemic influenza A (H1N1): role of the host immune response in pathogenesis
The American Journal of Pathology 183:1258–1268.
https://doi.org/10.1016/j.ajpath.2013.06.023
- PubMed
- Google Scholar
1. Gill JR
2. Sheng ZM
3. Ely SF
4. Guinee DG
5. Beasley MB
6. Suh J
7. Deshpande C
8. Mollura DJ
9. Morens DM
10. Bray M
11. Travis WD
12. Taubenberger JK
(2010) Pulmonary pathologic findings of fatal 2009 pandemic influenza A/H1N1 viral infections
Archives of Pathology & Laboratory Medicine 134:235–243.
https://doi.org/10.5858/134.2.235
- PubMed
- Google Scholar
1. Granados A
2. Peci A
3. McGeer A
4. Gubbay JB
(2017) Influenza and rhinovirus viral load and disease severity in upper respiratory tract infections
Journal of Clinical Virology 86:14–19.
https://doi.org/10.1016/j.jcv.2016.11.008
- PubMed
- Google Scholar
1. Graw F
2. Regoes RR
(2009) Investigating CTL mediated killing with a 3D cellular automaton
PLOS Computational Biology 5:e1000466.
https://doi.org/10.1371/journal.pcbi.1000466
- PubMed
- Google Scholar
(2000) Influenza virus neuraminidase inhibitors
The Lancet 355:827–835.
https://doi.org/10.1016/S0140-6736(99)11433-8
- PubMed
- Google Scholar
1. Halle S
2. Keyser KA
3. Stahl FR
4. Busche A
5. Marquardt A
6. Zheng X
7. Galla M
8. Heissmeyer V
9. Heller K
10. Boelter J
11. Wagner K
12. Bischoff Y
13. Martens R
14. Braun A
15. Werth K
16. Uvarovskii A
17. Kempf H
18. Meyer-Hermann M
19. Arens R
20. Kremer M
21. Sutter G
22. Messerle M
23. Förster R
(2016) In Vivo Killing Capacity of Cytotoxic T Cells Is Limited and Involves Dynamic Interactions and T Cell Cooperativity
Immunity 44:233–245.
https://doi.org/10.1016/j.immuni.2016.01.010
- PubMed
- Google Scholar
1. Han A
2. Poon JL
3. Powers JH
4. Leidy NK
5. Yu R
6. Memoli MJ
(2018) Using the influenza Patient-reported outcome (FLU-PRO) diary to evaluate symptoms of influenza viral infection in a healthy human challenge model
BMC Infectious Diseases 18:353.
https://doi.org/10.1186/s12879-018-3220-8
- PubMed
- Google Scholar
(2018) Progress and trends in mathematical modelling of influenza A virus infections
Current Opinion in Systems Biology 12:30–36.
https://doi.org/10.1016/j.coisb.2018.08.009
- Google Scholar
1. Harrington LE
2. Galvan M
3. Baum LG
4. Altman JD
5. Ahmed R
(2000) Differentiating between memory and effector CD8 T cells by altered expression of cell surface O-glycans
Journal of Experimental Medicine 191:1241–1246.
https://doi.org/10.1084/jem.191.7.1241
- PubMed
- Google Scholar
1. Holder BP
2. Beauchemin CA
(2011) Exploring the effect of biological delays in kinetic models of influenza within a host or cell culture
BMC Public Health 11 Suppl 1:S10.
https://doi.org/10.1186/1471-2458-11-S1-S10
- PubMed
- Google Scholar
1. Hou S
2. Doherty PC
3. Zijlstra M
4. Jaenisch R
5. Katz JM
(1992)
Delayed clearance of Sendai virus in mice lacking class I MHC-restricted CD8 T cells

Journal of Immunology 149:1319–1325.
- PubMed
- Google Scholar
1. Huang CT
2. Hung CY
3. Chen TC
4. Lin CY
5. Lin YC
6. Chang CS
7. He YC
8. Huang YL
9. Dutta A
(2017) Rapamycin adjuvant and exacerbation of severe influenza in an experimental mouse model
Scientific Reports 7:4136.
https://doi.org/10.1038/s41598-017-04365-6
- PubMed
- Google Scholar
1. Iwasaki A
2. Medzhitov R
(2010) Regulation of adaptive immunity by the innate immune system
Science 327:291–295.
https://doi.org/10.1126/science.1183021
- Google Scholar
1. Iwasaki A
2. Medzhitov R
(2015) Control of adaptive immunity by the innate immune system
Nature Immunology 16:343–353.
https://doi.org/10.1038/ni.3123
- PubMed
- Google Scholar
1. Iwasaki T
2. Nozima T
(1977)
Defense mechanisms against primary influenza virus infection in mice: I. the roles of interferon and neutralizing antibodies and Thymus dependence of interferon and antibody production

Journal of Immunology 118:256–263.
- PubMed
- Google Scholar
1. Jain A
2. Pasare C
(2017) Innate control of adaptive immunity: beyond the Three-Signal paradigm
The Journal of Immunology 198:3791–3800.
https://doi.org/10.4049/jimmunol.1602000
- PubMed
- Google Scholar
1. Jones AT
2. Federsppiel B
3. Ellies LG
4. Williams MJ
5. Burgener R
6. Duronio V
7. Smith CA
8. Takei F
9. Ziltener HJ
(1994)
Characterization of the activation-associated isoform of CD43 on murine T lymphocytes

Journal of Immunology 153:3426–3439.
- PubMed
- Google Scholar
1. Kaech SM
2. Hemby S
3. Kersh E
4. Ahmed R
(2002) Molecular and functional profiling of memory CD8 T cell differentiation
Cell 111:837–851.
https://doi.org/10.1016/S0092-8674(02)01139-X
- PubMed
- Google Scholar
1. Kang SS
2. Herz J
3. Kim JV
4. Nayak D
5. Stewart-Hutchinson P
6. Dustin ML
7. McGavern DB
(2011) Migration of cytotoxic lymphocytes in cell cycle permits local MHC I-dependent control of division at sites of viral infection
Journal of Experimental Medicine 208:747–759.
https://doi.org/10.1084/jem.20101295
- PubMed
- Google Scholar
(2015) Visualizing real-time influenza virus infection, transmission and protection in ferrets
Nature Communications 6:6378.
https://doi.org/10.1038/ncomms7378
- PubMed
- Google Scholar
1. Keating R
2. Morris MY
3. Yue W
4. Reynolds CE
5. Harris TL
6. Brown SA
7. Doherty PC
8. Thomas PG
9. McGargill MA
(2018) Potential killers exposed: tracking endogenous influenza-specific CD8⁺ T cells
Immunology and Cell Biology 96:1104–1119.
https://doi.org/10.1111/imcb.12189
- PubMed
- Google Scholar
(2010) Imaging appearance of swine-origin influenza A (novel 2009 H1N1) pneumonia in an immunocompromised patient
Radiology Case Reports 5:361.
https://doi.org/10.2484/rcr.v5i1.361
- PubMed
- Google Scholar
(2018) Safety, tolerability, and pharmacokinetics of the novel Anti-influenza agent baloxavir marboxil in healthy adults: phase I study findings
Clinical Drug Investigation 38:1189–1196.
https://doi.org/10.1007/s40261-018-0710-9
- PubMed
- Google Scholar
(2011) Immune responses to influenza virus infection
Virus Research 162:19–30.
https://doi.org/10.1016/j.virusres.2011.09.022
- PubMed
- Google Scholar
1. Kugel D
2. Kochs G
3. Obojes K
4. Roth J
5. Kobinger GP
6. Kobasa D
7. Haller O
8. Staeheli P
9. von Messling V
(2009) Intranasal administration of alpha interferon reduces seasonal influenza A virus morbidity in ferrets
Journal of Virology 83:3843–3851.
https://doi.org/10.1128/JVI.02453-08
- PubMed
- Google Scholar
(2007) A question of self-preservation: immunopathology in influenza virus infection
Immunology & Cell Biology 85:85–92.
https://doi.org/10.1038/sj.icb.7100026
- PubMed
- Google Scholar
1. La Gruta NL
2. Turner SJ
(2014) T cell mediated immunity to influenza: mechanisms of viral control
Trends in Immunology 35:396–402.
https://doi.org/10.1016/j.it.2014.06.004
- PubMed
- Google Scholar
1. Lambert Emo K
2. Hyun YM
3. Reilly E
4. Barilla C
5. Gerber S
6. Fowell D
7. Kim M
8. Topham DJ
(2016) Live imaging of influenza infection of the Trachea reveals dynamic regulation of CD8+ T cell motility by antigen
PLOS Pathogens 12:e1005881.
https://doi.org/10.1371/journal.ppat.1005881
- PubMed
- Google Scholar
1. Lauder SN
2. Taylor PR
3. Clark SR
4. Evans RL
5. Hindley JP
6. Smart K
7. Leach H
8. Kidd EJ
9. Broadley KJ
10. Jones SA
11. Wise MP
12. Godkin AJ
13. O'Donnell V
14. Gallimore AM
(2011) Paracetamol reduces influenza-induced immunopathology in a mouse model of infection without compromising virus clearance or the generation of protective immunity
Thorax 66:368–374.
https://doi.org/10.1136/thx.2010.150318
- PubMed
- Google Scholar
(2005) Frequency, specificity, and sites of expansion of CD8+ T cells during primary pulmonary influenza virus infection
The Journal of Immunology 174:5332–5340.
https://doi.org/10.4049/jimmunol.174.9.5332
- PubMed
- Google Scholar
(2014) Mathematical modeling provides kinetic details of the human immune response to vaccination
Frontiers in Cellular and Infection Microbiology 4:177.
https://doi.org/10.3389/fcimb.2014.00177
- PubMed
- Google Scholar
1. Le Bon A
2. Tough DF
(2002) Links between innate and adaptive immunity via type I interferon
Current Opinion in Immunology 14:432–436.
https://doi.org/10.1016/S0952-7915(02)00354-0
- PubMed
- Google Scholar
1. Lee HY
2. Topham DJ
3. Park SY
4. Hollenbaugh J
5. Treanor J
6. Mosmann TR
7. Jin X
8. Ward BM
9. Miao H
10. Holden-Wiltse J
11. Perelson AS
12. Zand M
13. Wu H
(2009) Simulation and prediction of the adaptive immune response to influenza A virus infection
Journal of Virology 83:7151–7165.
https://doi.org/10.1128/JVI.00098-09
- PubMed
- Google Scholar
1. Lee CS
2. Lee JH
(2010) Dynamics of clinical symptoms in patients with pandemic influenza A (H1N1)
Clinical Microbiology and Infection 16:389–390.
https://doi.org/10.1111/j.1469-0691.2010.03117.x
- PubMed
- Google Scholar
1. Levin D
2. Forrest S
3. Banerjee S
4. Clay C
5. Cannon J
6. Moses M
7. Koster F
(2016) A spatial model of the efficiency of T cell search in the influenza-infected lung
Journal of Theoretical Biology 398:52–63.
https://doi.org/10.1016/j.jtbi.2016.02.022
- PubMed
- Google Scholar
1. Li P
2. Su DJ
3. Zhang JF
4. Xia XD
5. Sui H
6. Zhao DH
(2011) Pneumonia in novel swine-origin influenza A (H1N1) virus infection: high-resolution CT findings
European Journal of Radiology 80:e146–e152.
https://doi.org/10.1016/j.ejrad.2010.05.029
- PubMed
- Google Scholar
1. Li P
2. Zhang JF
3. Xia XD
4. Su DJ
5. Liu BL
6. Zhao DL
7. Liu Y
8. Zhao DH
(2012) Serial evaluation of high-resolution CT findings in patients with pneumonia in novel swine-origin influenza A (H1N1) virus infection
The British Journal of Radiology 85:729–735.
https://doi.org/10.1259/bjr/85580974
- PubMed
- Google Scholar
1. Li Y
2. Handel A
(2014) Modeling inoculum dose dependent patterns of acute virus infections
Journal of Theoretical Biology 347:63–73.
https://doi.org/10.1016/j.jtbi.2014.01.008
- PubMed
- Google Scholar
1. Lu J
2. Duan X
3. Zhao W
4. Wang J
5. Wang H
6. Zhou K
7. Fang M
(2018) Aged mice are more resistant to influenza virus infection due to reduced inflammation and lung pathology
Aging and Disease 9:358.
https://doi.org/10.14336/AD.2017.0701
- PubMed
- Google Scholar
1. Luker KE
2. Luker GD
(2010) Bioluminescence imaging of reporter mice for studies of infection and inflammation
Antiviral Research 86:93–100.
https://doi.org/10.1016/j.antiviral.2010.02.002
- PubMed
- Google Scholar
1. Luster AD
(2002) The role of chemokines in linking innate and adaptive immunity
Current Opinion in Immunology 14:129–135.
https://doi.org/10.1016/S0952-7915(01)00308-9
- Google Scholar
(2014) Within-host influenza dynamics: a small-scale mathematical modeling approach
Biosystems 118:51–59.
https://doi.org/10.1016/j.biosystems.2014.02.004
- PubMed
- Google Scholar
(2010) Analysis of in vivo dynamics of influenza virus infection in mice using a GFP reporter virus
PNAS 107:11531–11536.
https://doi.org/10.1073/pnas.0914994107
- PubMed
- Google Scholar
1. Marathe BM
2. Wong SS
3. Vogel P
4. Garcia-Alcalde F
5. Webster RG
6. Webby RJ
7. Najera I
8. Govorkova EA
(2016) Combinations of oseltamivir and T-705 extend the treatment window for highly pathogenic influenza A(H5N1) Virus infection in mice
Scientific Reports 6:26742.
https://doi.org/10.1038/srep26742
- PubMed
- Google Scholar
1. Marathe BM
2. Mostafa HH
3. Vogel P
4. Pascua PNQ
5. Jones JC
6. Russell CJ
7. Webby RJ
8. Govorkova EA
(2017) A pharmacologically immunosuppressed mouse model for assessing influenza B virus pathogenicity and oseltamivir treatment
Antiviral Research 148:20–31.
https://doi.org/10.1016/j.antiviral.2017.10.021
- PubMed
- Google Scholar
1. Margine I
2. Krammer F
(2014) Animal models for influenza viruses: implications for universal vaccine development
Pathogens 3:845–874.
https://doi.org/10.3390/pathogens3040845
- PubMed
- Google Scholar
1. Matheu MP
2. Teijaro JR
3. Walsh KB
4. Greenberg ML
5. Marsolais D
6. Parker I
7. Rosen H
8. Oldstone MB
9. Cahalan MD
(2013) Three phases of CD8 T cell response in the lung following H1N1 influenza infection and sphingosine 1 phosphate agonist therapy
PLOS ONE 8:e58033.
https://doi.org/10.1371/journal.pone.0058033
- PubMed
- Google Scholar
1. Mauad T
2. Hajjar LA
3. Callegari GD
4. da Silva LF
5. Schout D
6. Galas FR
7. Alves VA
8. Malheiros DM
9. Auler JO
10. Ferreira AF
11. Borsato MR
12. Bezerra SM
13. Gutierrez PS
14. Caldini ET
15. Pasqualucci CA
16. Dolhnikoff M
17. Saldiva PH
(2010) Lung pathology in fatal novel human influenza A (H1N1) infection
American Journal of Respiratory and Critical Care Medicine 181:72–79.
https://doi.org/10.1164/rccm.200909-1420OC
- PubMed
- Google Scholar
1. McAuley JL
2. Hornung F
3. Boyd KL
4. Smith AM
5. McKeon R
6. Bennink J
7. Yewdell JW
8. McCullers JA
(2007) Expression of the 1918 influenza A virus PB1-F2 enhances the pathogenesis of viral and secondary bacterial pneumonia
Cell Host & Microbe 2:240–249.
https://doi.org/10.1016/j.chom.2007.09.001
- PubMed
- Google Scholar
1. McGill J
2. Legge KL
(2009) Cutting edge: contribution of lung-resident T cell proliferation to the overall magnitude of the antigen-specific CD8 T cell response in the lungs following murine influenza virus infection
The Journal of Immunology 183:4177–4181.
https://doi.org/10.4049/jimmunol.0901109
- PubMed
- Google Scholar
(1983) Cytotoxic T-cell immunity to influenza
New England Journal of Medicine 309:13–17.
https://doi.org/10.1056/NEJM198307073090103
- PubMed
- Google Scholar
1. Medina RA
2. García-Sastre A
(2011) Influenza A viruses: new research developments
Nature Reviews Microbiology 9:590–603.
https://doi.org/10.1038/nrmicro2613
- PubMed
- Google Scholar
1. Memoli MJ
2. Athota R
3. Reed S
4. Czajkowski L
5. Bristol T
6. Proudfoot K
7. Hagey R
8. Voell J
9. Fiorentino C
10. Ademposi A
11. Shoham S
12. Taubenberger JK
(2014) The natural history of influenza infection in the severely immunocompromised vs nonimmunocompromised hosts
Clinical Infectious Diseases 58:214–224.
https://doi.org/10.1093/cid/cit725
- PubMed
- Google Scholar
(2006) Regulatory T cells reversibly suppress cytotoxic T cell function independent of effector differentiation
Immunity 25:129–141.
https://doi.org/10.1016/j.immuni.2006.04.015
- PubMed
- Google Scholar
1. Merrill SJ
(1982) Foundations of the use of an enzyme-kinetic analogy in cell-mediated cytotoxicity
Mathematical Biosciences 62:219–235.
https://doi.org/10.1016/0025-5564(82)90084-0
- Google Scholar
1. Mestas J
2. Hughes CC
(2004) Of mice and not men: differences between mouse and human immunology
The Journal of Immunology 172:2731–2738.
https://doi.org/10.4049/jimmunol.172.5.2731
- PubMed
- Google Scholar
1. Miao H
2. Hollenbaugh JA
3. Zand MS
4. Holden-Wiltse J
5. Mosmann TR
6. Perelson AS
7. Wu H
8. Topham DJ
(2010) Quantifying the early immune response and adaptive immune response kinetics in mice infected with influenza A virus
Journal of Virology 84:6687–6698.
https://doi.org/10.1128/JVI.00266-10
- PubMed
- Google Scholar
(2000) Clinical signs and symptoms predicting influenza infection
Archives of Internal Medicine 160:3243–3247.
https://doi.org/10.1001/archinte.160.21.3243
- PubMed
- Google Scholar
1. Moscona A
(2005) Neuraminidase inhibitors for influenza
New England Journal of Medicine 353:1363–1373.
https://doi.org/10.1056/NEJMra050740
- PubMed
- Google Scholar
1. Moskophidis D
2. Kioussis D
(1998) Contribution of virus-specific CD8+ cytotoxic T cells to virus clearance or pathologic manifestations of influenza virus infection in a T cell receptor transgenic mouse model
Journal of Experimental Medicine 188:223–232.
https://doi.org/10.1084/jem.188.2.223
- PubMed
- Google Scholar
(2001) Small variations in multiple parameters account for wide variations in HIV–1 set–points: a novel modelling approach
Proceedings of the Royal Society of London. Series B: Biological Sciences 268:235–242.
https://doi.org/10.1098/rspb.2000.1358
- Google Scholar
1. Music N
2. Reber AJ
3. Lipatov AS
4. Kamal RP
5. Blanchfield K
6. Wilson JR
7. Donis RO
8. Katz JM
9. York IA
(2014) Influenza vaccination accelerates recovery of ferrets from lymphopenia
PLOS ONE 9:e100926.
https://doi.org/10.1371/journal.pone.0100926
- PubMed
- Google Scholar
1. Onami TM
2. Harrington LE
3. Williams MA
4. Galvan M
5. Larsen CP
6. Pearson TC
7. Manjunath N
8. Baum LG
9. Pearce BD
10. Ahmed R
(2002) Dynamic regulation of T cell immunity by CD43
Journal of Immunology 168:6022–6031.
https://doi.org/10.4049/jimmunol.168.12.6022
- PubMed
- Google Scholar
(2002) Virus clearance and immunopathology by CD8(+) T cells during infection with respiratory syncytial virus are mediated by IFN-gamma
European Journal of Immunology 32:2117–2123.
https://doi.org/10.1002/1521-4141(200208)32:8<2117::AID-IMMU2117>3.0.CO;2-C
- PubMed
- Google Scholar
1. Parzych EM
2. DiMenna LJ
3. Latimer BP
4. Small JC
5. Kannan S
6. Manson B
7. Lasaro MO
8. Wherry EJ
9. Ertl HC
(2013) Influenza virus specific CD8⁺ T cells exacerbate infection following high dose influenza challenge of aged mice
BioMed Research International 2013:876314.
https://doi.org/10.1155/2013/876314
- PubMed
- Google Scholar
1. Perelson AS
2. Bell GI
(1982)
Delivery of lethal hits by cytotoxic T lymphocytes in multicellular conjugates occurs sequentially but at random times

Journal of Immunology 129:2796–2801.
- PubMed
- Google Scholar
Software
1. Perry M
(2018)
Simanneal, version 0.4.1

Simanneal.
1. Pilyugin SS
2. Antia R
(2000) Modeling immune responses with handling time
Bulletin of Mathematical Biology 62:869–890.
https://doi.org/10.1006/bulm.2000.0181
- PubMed
- Google Scholar
1. Price I
2. Mochan-Keef ED
3. Swigon D
4. Ermentrout GB
5. Lukens S
6. Toapanta FR
7. Ross TM
8. Clermont G
(2015) The inflammatory response to influenza A virus (H1N1): An experimental and mathematical study
Journal of Theoretical Biology 374:83–93.
https://doi.org/10.1016/j.jtbi.2015.03.017
- PubMed
- Google Scholar
(2019) Mathematical modeling of energy consumption in the acute inflammatory response
Journal of Theoretical Biology 460:101–114.
https://doi.org/10.1016/j.jtbi.2018.08.033
- PubMed
- Google Scholar
1. Reed LJ
2. Muench H
(1938) A simple method of estimating fifty per cent endpoints12
American Journal of Epidemiology 27:493–497.
https://doi.org/10.1093/oxfordjournals.aje.a118408
- Google Scholar
(2006) A reduced mathematical model of the acute inflammatory response: I. Derivation of model and analysis of anti-inflammation
Journal of Theoretical Biology 242:220–236.
https://doi.org/10.1016/j.jtbi.2006.02.016
- PubMed
- Google Scholar
1. Russell DG
2. Cardona PJ
3. Kim MJ
4. Allain S
5. Altare F
(2009) Foamy macrophages and the progression of the human tuberculosis granuloma
Nature Immunology 10:943.
https://doi.org/10.1038/ni.1781
- PubMed
- Google Scholar
(2011) Viral pneumonia
Lancet 377:1264–1275.
https://doi.org/10.1016/S0140-6736(10)61459-6
- PubMed
- Google Scholar
1. Rygiel TP
2. Rijkers ES
3. de Ruiter T
4. Stolte EH
5. van der Valk M
6. Rimmelzwaan GF
7. Boon L
8. van Loon AM
9. Coenjaerts FE
10. Hoek RM
11. Tesselaar K
12. Meyaard L
(2009) Lack of CD200 enhances pathological T cell responses during influenza infection
The Journal of Immunology 183:1990–1996.
https://doi.org/10.4049/jimmunol.0900252
- PubMed
- Google Scholar
1. Sartorius A
2. Lu Q
3. Vieira S
4. Tonnellier M
5. Lenaour G
6. Goldstein I
7. Rouby JJ
(2007) Mechanical ventilation and lung infection in the genesis of air-space enlargement
Critical Care 11:R14.
https://doi.org/10.1186/cc5680
- PubMed
- Google Scholar
1. Shieh WJ
2. Blau DM
3. Denison AM
4. Deleon-Carnes M
5. Adem P
6. Bhatnagar J
7. Sumner J
8. Liu L
9. Patel M
10. Batten B
11. Greer P
12. Jones T
13. Smith C
14. Bartlett J
15. Montague J
16. White E
17. Rollin D
18. Gao R
19. Seales C
20. Jost H
21. Metcalfe M
22. Goldsmith CS
23. Humphrey C
24. Schmitz A
25. Drew C
26. Paddock C
27. Uyeki TM
28. Zaki SR
(2010) 2009 pandemic influenza A (H1N1): pathology and pathogenesis of 100 fatal cases in the United States
The American Journal of Pathology 177:166–175.
https://doi.org/10.2353/ajpath.2010.100115
- PubMed
- Google Scholar
(2000) The impact of influenza epidemics on hospitalizations
The Journal of Infectious Diseases 181:831–837.
https://doi.org/10.1086/315320
- PubMed
- Google Scholar
(2017) Dynamics of influenza-induced lung-resident memory T cells underlie waning heterosubtypic immunity
Science Immunology 2:eaag2031.
https://doi.org/10.1126/sciimmunol.aag2031
- PubMed
- Google Scholar
(2010) An accurate two-phase approximate solution to an acute viral infection model
Journal of Mathematical Biology 60:711–726.
https://doi.org/10.1007/s00285-009-0281-8
- PubMed
- Google Scholar
(2011a) Effect of 1918 PB1-F2 expression on influenza A virus infection kinetics
PLOS Computational Biology 7:e1001081.
https://doi.org/10.1371/journal.pcbi.1001081
- PubMed
- Google Scholar
(2011b) Mathematical model of a three-stage innate immune response to a pneumococcal lung infection
Journal of Theoretical Biology 276:106–116.
https://doi.org/10.1016/j.jtbi.2011.01.052
- PubMed
- Google Scholar
(2013) Kinetics of coinfection with influenza A virus and Streptococcus pneumoniae
PLOS Pathogens 9:e1003238.
https://doi.org/10.1371/journal.ppat.1003238
- PubMed
- Google Scholar
(2018) Influenza Virus Infection Model With Density Dependence Supports Biphasic Viral Decay
Frontiers in Microbiology 9:1554.
https://doi.org/10.3389/fmicb.2018.01554
- PubMed
- Google Scholar
1. Smith AM
(2018a) Validated models of immune response to virus infection
Current Opinion in Systems Biology 12:46–52.
https://doi.org/10.1016/j.coisb.2018.10.005
- PubMed
- Google Scholar
1. Smith AM
(2018b) Host-pathogen kinetics during influenza infection and coinfection: insights from predictive modeling
Immunological Reviews 285:97–112.
https://doi.org/10.1111/imr.12692
- PubMed
- Google Scholar
(2019) Influenza virus inoculum volume is critical to elucidate age-dependent mortality in mice
Aging Cell 18:e12893.
https://doi.org/10.1111/acel.12893
- PubMed
- Google Scholar
1. Smith AM
2. McCullers JA
(2014) Secondary bacterial infections in influenza virus infection pathogenesis
Current Topics in Microbiology and Immunology 385:327–356.
https://doi.org/10.1007/82_2014_394
- PubMed
- Google Scholar
1. Smith AM
2. Perelson AS
(2011) Influenza A virus infection kinetics: quantitative data and models
Wiley Interdisciplinary Reviews: Systems Biology and Medicine 3:429–445.
https://doi.org/10.1002/wsbm.129
- PubMed
- Google Scholar
(2009) Pathological changes associated with the 2009 H1N1 virus
New England Journal of Medicine 361:2001–2003.
https://doi.org/10.1056/NEJMc0907171
- Google Scholar
(2009) Host genetic background strongly influences the response to influenza a virus infections
PLOS ONE 4:e4857.
https://doi.org/10.1371/journal.pone.0004857
- PubMed
- Google Scholar
1. Sun J
2. Madan R
3. Karp CL
4. Braciale TJ
(2009) Effector T cells control lung inflammation during acute influenza virus infection by producing IL-10
Nature Medicine 15:277–284.
https://doi.org/10.1038/nm.1929
- PubMed
- Google Scholar
1. Szretter KJ
2. Gangappa S
3. Belser JA
4. Zeng H
5. Chen H
6. Matsuoka Y
7. Sambhara S
8. Swayne DE
9. Tumpey TM
10. Katz JM
(2009) Early control of H5N1 influenza virus replication by the type I interferon response in mice
Journal of Virology 83:5825–5834.
https://doi.org/10.1128/JVI.02144-08
- PubMed
- Google Scholar
1. Tang BM
2. Shojaei M
3. Teoh S
4. Meyers A
5. Ho J
6. Ball TB
7. Keynan Y
8. Pisipati A
9. Kumar A
10. Eisen DP
11. Lai K
12. Gillett M
13. Santram R
14. Geffers R
15. Schreiber J
16. Mozhui K
17. Huang S
18. Parnell GP
19. Nalos M
20. Holubova M
21. Chew T
22. Booth D
23. Kumar A
24. McLean A
25. Schughart K
(2019) Neutrophils-related host factors associated with severe disease and fatality in patients with influenza infection
Nature Communications 10:3422.
https://doi.org/10.1038/s41467-019-11249-y
- PubMed
- Google Scholar
1. Taubenberger JK
2. Morens DM
(2008) The pathology of influenza virus infections
Annual Review of Pathology: Mechanisms of Disease 3:499–522.
https://doi.org/10.1146/annurev.pathmechdis.3.121806.154316
- PubMed
- Google Scholar
1. Thangavel RR
2. Bouvier NM
(2014) Animal models for influenza virus pathogenesis, transmission, and immunology
Journal of Immunological Methods 410:60–79.
https://doi.org/10.1016/j.jim.2014.03.023
- PubMed
- Google Scholar
1. Thompson WW
2. Shay DK
3. Weintraub E
4. Brammer L
5. Bridges CB
6. Cox NJ
7. Fukuda K
(2004) Influenza-associated hospitalizations in the united states
Jama 292:1333–1340.
https://doi.org/10.1001/jama.292.11.1333
- PubMed
- Google Scholar
1. Toapanta FR
2. Ross TM
(2009) Impaired immune responses in the lungs of aged mice following influenza infection
Respiratory Research 10:112.
https://doi.org/10.1186/1465-9921-10-112
- PubMed
- Google Scholar
1. Toniolo Neto J
(2014) Neuraminidase inhibitors for preventing and treating influenza in healthy adults and children
Sao Paulo Medical Journal 132:256–257.
https://doi.org/10.1590/1516-3180.20141324T2
- PubMed
- Google Scholar
1. Trammell RA
2. Toth LA
(2011)
Markers for predicting death as an outcome for mice used in infectious disease research

Comparative Medicine 61:492–498.
- PubMed
- Google Scholar
1. Tran V
2. Moser LA
3. Poole DS
4. Mehle A
(2013) Highly sensitive real-time in vivo imaging of an influenza reporter virus reveals dynamics of replication and spread
Journal of Virology 87:13321–13329.
https://doi.org/10.1128/JVI.02381-13
- PubMed
- Google Scholar
1. Treanor JJ
2. Hayden FG
3. Vrooman PS
4. Barbarash R
5. Bettis R
6. Riff D
7. Singh S
8. Kinnersley N
9. Ward P
10. Mills RG
(2000) Efficacy and safety of the oral neuraminidase inhibitor oseltamivir in treating acute influenza: a randomized controlled trial. US Oral Neuraminidase Study Group
JAMA 283:1016–1024.
https://doi.org/10.1001/jama.283.8.1016
- PubMed
- Google Scholar
(2021) Long-term maintenance of lung resident memory T cells is mediated by persistent antigen
Mucosal Immunology 14:92–99.
https://doi.org/10.1038/s41385-020-0309-3
- PubMed
- Google Scholar
1. Ulrichs T
2. Kosmiadi GA
3. Trusov V
4. Jörg S
5. Pradl L
6. Titukhina M
7. Mishenko V
8. Gushina N
9. Kaufmann SH
(2004) Human tuberculous granulomas induce peripheral lymphoid follicle-like structures to orchestrate local host defence in the lung
The Journal of Pathology 204:217–228.
https://doi.org/10.1002/path.1628
- PubMed
- Google Scholar
(2018) Repeated Antigen Exposure Extends the Durability of Influenza-Specific Lung-Resident Memory CD8⁺ T Cells and Heterosubtypic Immunity
Cell Reports 24:3374–3382.
https://doi.org/10.1016/j.celrep.2018.08.073
- PubMed
- Google Scholar
1. van de Sandt CE
2. Bárcena M
3. Koster AJ
4. Kasper J
5. Kirkpatrick CJ
6. Scott DP
7. de Vries RD
8. Herold S
9. Rimmelzwaan GF
10. Kuiken T
11. Short KR
(2017) Human CD8⁺ T Cells Damage Noninfected Epithelial Cells during Influenza Virus Infection In Vitro
American Journal of Respiratory Cell and Molecular Biology 57:536–546.
https://doi.org/10.1165/rcmb.2016-0377OC
- PubMed
- Google Scholar
(2015) Antibody-targeted vaccination to lung dendritic cells generates tissue-resident memory CD8 T cells that are highly protective against influenza virus infection
Mucosal Immunology 8:1060–1071.
https://doi.org/10.1038/mi.2014.133
- PubMed
- Google Scholar
1. Wang Z
2. Wan Y
3. Qiu C
4. Quiñones-Parra S
5. Zhu Z
6. Loh L
7. Tian D
8. Ren Y
9. Hu Y
10. Zhang X
11. Thomas PG
12. Inouye M
13. Doherty PC
14. Kedzierska K
15. Xu J
(2015) Recovery from severe H7N9 disease is associated with diverse response mechanisms dominated by CD8⁺ T cells
Nature Communications 6:6833.
https://doi.org/10.1038/ncomms7833
- PubMed
- Google Scholar
(2019) The role of macrophages in the resolution of inflammation
The Journal of Clinical Investigation 129:2619–2628.
https://doi.org/10.1172/JCI124615
- PubMed
- Google Scholar
1. Weinfurter JT
2. Brunner K
3. Capuano SV
4. Li C
5. Broman KW
6. Kawaoka Y
7. Friedrich TC
(2011) Cross-reactive T cells are involved in rapid clearance of 2009 pandemic H1N1 influenza virus in nonhuman primates
PLOS Pathogens 7:e1002381.
https://doi.org/10.1371/journal.ppat.1002381
- PubMed
- Google Scholar
(1981)
Recovery from a viral respiratory infection. I. Influenza pneumonia in normal and T-deficient mice

Journal of Immunology 126:1036–1041.
- PubMed
- Google Scholar
(2006) Cytotoxic T lymphocytes kill multiple targets simultaneously via spatiotemporal uncoupling of lytic and stimulatory synapses
PNAS 103:10985–10990.
https://doi.org/10.1073/pnas.0600651103
- PubMed
- Google Scholar
1. Wilkinson TM
2. Li CK
3. Chui CS
4. Huang AK
5. Perkins M
6. Liebner JC
7. Lambkin-Williams R
8. Gilbert A
9. Oxford J
10. Nicholas B
11. Staples KJ
12. Dong T
13. Douek DC
14. McMichael AJ
15. Xu XN
(2012) Preexisting influenza-specific CD4+ T cells correlate with disease protection against influenza challenge in humans
Nature Medicine 18:274–280.
https://doi.org/10.1038/nm.2612
- PubMed
- Google Scholar
1. Wong HYF
2. Lam HYS
3. Fong AH
4. Leung ST
5. Chin TW
6. Lo CSY
7. Lui MM
8. Lee JCY
9. Chiu KW
10. Chung TW
11. Lee EYP
12. Wan EYF
13. Hung IFN
14. Lam TPW
15. Kuo MD
16. Ng MY
(2020) Frequency and Distribution of Chest Radiographic Findings in Patients Positive for COVID-19
Radiology 296:E72–E78.
https://doi.org/10.1148/radiol.2020201160
- PubMed
- Google Scholar
1. Wu H
2. Kumar A
3. Miao H
4. Holden-Wiltse J
5. Mosmann TR
6. Livingstone AM
7. Belz GT
8. Perelson AS
9. Zand MS
10. Topham DJ
(2011) Modeling of influenza-specific CD8+ T cells during the primary response indicates that the spleen is a major source of effectors
Journal of Immunology 187:4474–4482.
https://doi.org/10.4049/jimmunol.1101443
- PubMed
- Google Scholar
1. Xu L
2. Yoon H
3. Zhao MQ
4. Liu J
5. Ramana CV
6. Enelow RI
(2004) Cutting edge: pulmonary immunopathology mediated by antigen-specific expression of TNF-alpha by antiviral CD8+ T cells
Journal of Immunology 173:721–725.
https://doi.org/10.4049/jimmunol.173.2.721
- PubMed
- Google Scholar
1. Xue KS
2. Stevens-Ayers T
3. Campbell AP
4. Englund JA
5. Pergam SA
6. Boeckh M
7. Bloom JD
(2017) Parallel evolution of influenza across multiple spatiotemporal scales
eLife 6:e26875.
https://doi.org/10.7554/eLife.26875
- PubMed
- Google Scholar
(1986)
Reorientation and fusion of cytotoxic T lymphocyte granules after interaction with target cells as determined by high resolution cinemicrography

Journal of Immunology 136:377–382.
- PubMed
- Google Scholar
1. Yap KL
2. Ada GL
(1978) Cytotoxic T cells in the lungs of mice infected with an influenza A virus
Scandinavian Journal of Immunology 7:73–80.
https://doi.org/10.1111/j.1365-3083.1978.tb00428.x
- PubMed
- Google Scholar
(2011) Virus replication strategies and the critical CTL numbers required for the control of infection
PLOS Computational Biology 7:e1002274.
https://doi.org/10.1371/journal.pcbi.1002274
- PubMed
- Google Scholar
1. Yoon H
2. Legge KL
3. Sung SS
4. Braciale TJ
(2007) Sequential activation of CD8+ T cells in the draining lymph nodes in response to pulmonary virus infection
Journal of Immunology 179:391–399.
https://doi.org/10.4049/jimmunol.179.1.391
- PubMed
- Google Scholar
(2006) gammadelta T cells regulate the early inflammatory response to bordetella pertussis infection in the murine respiratory tract
Infection and Immunity 74:1837–1845.
https://doi.org/10.1128/IAI.74.3.1837-1845.2006
- PubMed
- Google Scholar
1. Zagury D
2. Bernard J
3. Thierness N
4. Feldman M
5. Berke G
(1975) Isolation and characterization of individual functionally reactive cytotoxic T lymphocytes: conjugation, killing and recycling at the single cell level
European Journal of Immunology 5:818–822.
https://doi.org/10.1002/eji.1830051205
- Google Scholar
(2016) Vaccine-generated lung tissue–resident memory T cells provide heterosubtypic protection to influenza infection
JCI Insight 1:e85832.
https://doi.org/10.1172/jci.insight.85832
- Google Scholar
1. Zhang N
2. Bevan MJ
(2011) CD8(+) T cells: foot soldiers of the immune system
Immunity 35:161–168.
https://doi.org/10.1016/j.immuni.2011.07.010
- PubMed
- Google Scholar
1. Zhu C
2. Byrd RH
3. Lu P
4. Nocedal J
(1997) Algorithm 778: L-BFGS-B: Fortran subroutines for large-scale bound-constrained optimization
ACM T Math Software 23:550–560.
https://doi.org/10.1145/279232.279236
- Google Scholar

Article and author information

Author details

Margaret A Myers

Department of Pediatrics, University of Tennessee Health Science Center, Memphis, United States

Contribution
Software, Formal analysis, Investigation, Visualization, Writing - original draft

Contributed equally with
Amanda P Smith

Competing interests
No competing interests declared
Amanda P Smith

Department of Pediatrics, University of Tennessee Health Science Center, Memphis, United States

Contribution
Data curation, Formal analysis, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

Contributed equally with
Margaret A Myers

Competing interests
No competing interests declared
Lindey C Lane

Department of Pediatrics, University of Tennessee Health Science Center, Memphis, United States

Contribution
Data curation

Competing interests
No competing interests declared
David J Moquin

Department of Anesthesiology, Washington University School of Medicine, St. Louis, United States

Contribution
Data curation

Competing interests
No competing interests declared
Rosemary Aogo

Department of Pediatrics, University of Tennessee Health Science Center, Memphis, United States

Contribution
Formal analysis, Writing - review and editing

Competing interests
No competing interests declared
Stacie Woolard

Flow Cytometry Core, St. Jude Children's Research Hospital, Memphis, United States

Contribution
Formal analysis

Competing interests
No competing interests declared
Paul Thomas

Department of Immunology, St. Jude Children’s Research Hospital, Memphis, United States

Contribution
Methodology, Writing - review and editing

Competing interests
No competing interests declared
Peter Vogel

Department of Pathology, St. Jude Children's Research Hospital, Memphis, United States

Contribution
Data curation, Formal analysis, Visualization, Methodology, Writing - original draft, Writing - review and editing

Competing interests
No competing interests declared
Amber M Smith

Department of Pediatrics, University of Tennessee Health Science Center, Memphis, United States

Contribution
Conceptualization, Resources, Data curation, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing - original draft, Writing - review and editing

For correspondence
amber.smith@uthsc.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-7092-6904

Funding

National Institute of Allergy and Infectious Diseases (AI139088)

Margaret A Myers
Amanda P Smith
Lindey C Lane
Rosemary Aogo

National Institute of Allergy and Infectious Diseases (AI125324)

Margaret A Myers
Amanda P Smith
Lindey C Lane
David J Moquin
Amber M Smith

National Institute of Allergy and Infectious Diseases (AI100946)

Amber M Smith

American Lebanese Syrian Associated Charities (Internal Funding)

Margaret A Myers
Amanda P Smith
Lindey C Lane
David J Moquin
Stacie Woolard
Paul Thomas
Peter Vogel
Amber M Smith

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This work was supported by NIH grants AI100946, AI125324, and AI139088 and ALSAC. We thank Alan Perelson for his helpful comments, and Robert Michael for technical assistance.

Ethics

Animal experimentation: All experimental procedures were performed under protocols O2A-020 or 17-096 approved by the Animal Care and Use Committees at St. Jude Children's Research Hospital (SJCRH) or the University of Tennessee Health Science Center (UTHSC), respectively, under relevant institutional and American Veterinary Medical Association (AVMA) guidelines. All experimental procedures were performed in a biosafety level 2 facility that is accredited by the American Association for Laboratory Animal Science (AALAS).

Version history

Preprint posted: February 19, 2019 (view preprint)
Received: March 28, 2021
Accepted: July 14, 2021
Accepted Manuscript published: July 20, 2021 (version 1)
Version of Record published: August 17, 2021 (version 2)

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.