You will be shown a paragraph from a biomedical research paper, followed by a question about the paper. Please extract text from the paragraph to answer the question. Do not include any text in your repsonse other than the answer.

Paragraph: We conducted plaque assays on six-well plate monolayers of three immortalized mammalian kidney cell lines: [1] Vero (African green monkey) cells, which are IFN-defective and thus limited in antiviral capacity (Desmyter et al., 1968) ; [2] RoNi/7.1 (Rousettus aegyptiacus) cells which demonstrate idiosyncratic induced interferon responses upon viral challenge  (Kuzmin et al., 2017; Arnold et al., 2018; Biesold et al., 2011; Pavlovich et al., 2018) ; and [3] PaKiT01 (Pteropus alecto) cells which constitutively express IFN-a (Zhou et al., 2016; Crameri et al., 2009) . To intensify cell line-specific differences in constitutive immunity, we carried out infectivity assays with GFP-tagged, replication-competent vesicular stomatitis Indiana viruses: rVSV-G, rVSV-EBOV, and rVSV-MARV, which have been previously described (Miller et al., 2012; Wong et al., 2010) . Two of these viruses, rVSV-EBOV and rVSV-MARV, are recombinants for which cell entry is mediated by the glycoprotein of the bat-evolved filoviruses, Ebola (EBOV) and Marburg (MARV), thus allowing us to modulate the extent of structural, as well as immunological, antiviral defense at play in each infection. Previous work in this lab has demonstrated incompatibilities in the NPC1 filovirus receptor which render PaKiT01 cells refractory to infection with rVSV-MARV (Ng and Chandrab, 2018, Unpublished results) , making them structurally antiviral, over and above their constitutive expression of IFN-a. All three cell lines were challenged with all three viruses at two multiplicities of infection (MOI): 0.001 and 0.0001. Between 18 and 39 trials were run at each cell-virus-MOI combination, excepting rVSV-MARV infections on PaKiT01 cells at MOI = 0.001, for which only eight trials were run (see Materials and methods; Figure 1 -figure supplements 1-3, Supplementary file 1). Because plaque assays restrict viral transmission neighbor-to-neighbor in two-dimensional cellular space (Howat et al., 2006) , we were able to track the spread of GFP-expressing virus-infected cells across tissue monolayers via inverted fluorescence microscopy. For each infection trial, we monitored and re-imaged plates for up to 200 hr of observations or until total monolayer destruction, processed resulting images, and generated a time series of the proportion of infectious-cell occupied plate space across the duration of each trial (see Materials and methods). We used generalized additive models to infer the time course of all cell culture replicates and construct the multi-trial dataset to which we eventually fit our mechanistic transmission model for each cell line-virus-specific combination ( Figure 1; Figure 1 -figure supplements 1-5). All three recombinant vesicular stomatitis viruses (rVSV-G, rVSV-EBOV, and rVSV-MARV) infected Vero, RoNi/7.1, and PaKiT01 tissue cultures at both focal MOIs. Post-invasion, virus spread rapidly across most cell monolayers, resulting in virus-induced epidemic extinction. Epidemics were less severe in bat cell cultures, especially when infected with the recombinant filoviruses, rVSV-EBOV and rVSV-MARV. Monolayer destruction was avoided in the case of rVSV-EBOV and rVSV-MARV infections on PaKiT01 cells: in the former, persistent viral infection was maintained throughout the 200 hr duration of each experiment, while, in the latter, infection was eliminated early in the time series, preserving a large proportion of live, uninfectious cells across the duration of the experiment. We assumed this pattern to be the result of immune-mediated epidemic extinction (Figure 1) . Patterns from MOI = 0.001 were largely recapitulated at MOI = 0.0001, though at somewhat reduced total proportions (Figure 1-figure supplement 5 ). A theoretical model fit to in vitro data recapitulates expected immune phenotypes for bat cells We next developed a within-host model to fit to these data to elucidate the effects of induced and constitutive immunity on the dynamics of viral spread in host tissue ( Figure 1 ). The compartmental within-host system mimicked our two-dimensional cell culture monolayer, with cells occupying five distinct infection states: susceptible (S), antiviral (A), exposed (E), infectious (I), and dead (D). We modeled exposed cells as infected but not yet infectious, capturing the 'eclipse phase' of viral integration into a host cell which precedes viral replication. Antiviral cells were immune to viral infection, in accordance with the 'antiviral state' induced from interferon stimulation of ISGs in tissues adjacent to infection (Stetson and Medzhitov, 2006) . Because we aimed to translate available data into modeled processes, we did not explicitly model interferon dynamics but instead scaled the rate of cell progression from susceptible to antiviral (r) by the proportion of exposed cells (globally) in the system. In systems permitting constitutive immunity, a second rate of cellular acquisition of antiviral status (") additionally scaled with the global proportion of susceptible cells in the model. Compared with virus, IFN particles are small and highly diffusive, justifying this global signaling assumption at the limited spatial extent of a six-well plate and maintaining consistency with previous modeling approximations of IFN signaling in plaque assay (Howat et al., 2006) . To best represent our empirical monolayer system, we expressed our state variables as proportions (P S , P A , P E , P I , and P D ), under assumptions of frequency-dependent transmission in a wellmixed population (Keeling and Rohani, 2008) , though note that the inclusion of P D (representing the proportion of dead space in the modeled tissue) had the functional effect of varying transmission with infectious cell density. This resulted in the following system of ordinary differential equations:

We defined 'induced immunity' as complete, modeling all cells as susceptible to viral invasion at disease-free equilibrium, with defenses induced subsequent to viral exposure through the term r. By contrast, we allowed the extent of constitutive immunity to vary across the parameter range of " > 0, defining a 'constitutive' system as one containing any antiviral cells at disease-free equilibrium. In fitting this model to tissue culture data, we independently estimated both r and "; as well as the cell-to-cell transmission rate, b, for each cell-virus combination. Since the extent to which constitutively-expressed IFN-a is constitutively translated into functional protein is not yet known for bat hosts (Zhou et al., 2016) , this approach permitted our tissue culture data to drive modeling inference: even in PaKiT01 cell lines known to constitutively express IFN-a, the true constitutive extent of the system (i.e. the quantity of antiviral cells present at disease-free equilibrium) was allowed to vary through estimation of ": For the purposes of model-fitting, we fixed the value of c, the return rate of antiviral cells to susceptible status, at 0. The small spatial scale and short time course (max 200 hours) of our experiments likely prohibited any return of antiviral cells to susceptible status in our empirical system; nonetheless, we retained the term c in analytical evaluations of our model because regression from antiviral to susceptible status is possible over long time periods in vitro and at the scale of a complete organism (Radke et al., 1974; Rasmussen and Farley, 1975; Samuel and Knutson, 1982) . Before fitting to empirical time series, we undertook bifurcation analysis of our theoretical model and generated testable hypotheses on the basis of model outcomes. From our within-host model system (Equation 1-5), we derived the following expression for R 0 , the pathogen basic reproduction number (Supplementary file 2):

Pathogens can invade a host tissue culture when R 0 >1. Rapid rates of constitutive antiviral acquisition (") will drive R 0 <1: tissue cultures with highly constitutive antiviral immunity will be therefore resistant to virus invasion from the outset. Since, by definition, induced immunity is stimulated following initial virus invasion, the rate of induced antiviral acquisition (r) is not incorporated into the equation for R 0 ; while induced immune processes can control virus after initial invasion, they cannot prevent it from occurring to begin with. In cases of fully induced or absent immunity (" ¼ 0), the R 0 equation thus reduces to a form typical of the classic SEIR model:

At equilibrium, the theoretical, mean field model demonstrates one of three infection states: endemic equilibrium, stable limit cycles, or no infection ( Figure 2) . Respectively, these states approximate the persistent infection, virus-induced epidemic extinction, and immune-mediated epidemic extinction phenotypes previously witnessed in tissue culture experiments ( Figure 1 ). Theoretically, endemic equilibrium is maintained when new infections are generated at the same rate at which infections are lost, while limit cycles represent parameter space under which infectious and susceptible populations are locked in predictable oscillations. Endemic equilibria resulting from cellular regeneration (i.e. births) have been described in vivo for HIV (Coffin, 1995) and in vitro for herpesvirus plaque assays (Howat et al., 2006) , but, because they so closely approach zero, true limit cycles likely only occur theoretically, instead yielding stochastic extinctions in empirical time series. Bifurcation analysis of our mean field model revealed that regions of no infection (pathogen extinction) were bounded at lower threshold (Branch point) values for b, below which the pathogen was unable to invade. We found no upper threshold to invasion for b under any circumstances (i.e. b high enough to drive pathogen-induced extinction), but high b values resulted in Hopf bifurcations, which delineate regions of parameter space characterized by limit cycles. Since limit cycles so closely approach zero, high bs recovered in this range would likely produce virus-induced epidemic extinctions under experimental conditions. Under more robust representations of immunity, with higher values for either or both induced (r) and constitutive (") rates of antiviral acquisition, Hopf bifurcations occurred at increasingly higher values for b, meaning that persistent infections could establish at higher viral transmission rates ( Figure 2 ). Consistent with our derivation for R 0 , we found that the Branch point threshold for viral invasion was independent of changes to the induced immune parameter (r) but saturated at high values of " that characterize highly constitutive immunity ( Figure 3) . We next fit our theoretical model by least squares to each cell line-virus combination, under absent, induced, and constitutive assumptions of immunity. In general, best fit models recapitulated expected outcomes based on the immune phenotype of the cell line in question, as described in the general literature (Table 1 Ironically, the induced immune model offered a slightly better fit than the constitutive to rVSV-MARV infections on the PaKiT01 cell line (the one cell line-virus combination for which we know a constitutively antiviral cell-receptor incompatibility to be at play). Because constitutive immune assumptions can prohibit pathogen invasion (R 0 <1), model fits to this time series under constitutive assumptions were handicapped by overestimations of ", which prohibited pathogen invasion. Only by incorporating an exceedingly rapid rate of induced antiviral acquisition could the model guarantee that initial infection would be permitted and then rapidly controlled. In all panel (A) plots, the rate of induced immune antiviral acquisition (r) was fixed at 0.01. Panel (B) depicts dynamics under variably induced immunity, ranging from absent (left: r=0) to high (right: r=1). In all panel (B) plots, the rate of constitutive antiviral acquisition (") was fixed at 0.0001 Branch point curves are represented as solid lines and Hopf curves as dashed lines. White space indicates endemic equilibrium (persistence), gray space indicates limit cycles, and black space indicates no infection (extinction). Other parameter values for equilibrium analysis were fixed at: b = .025, m = .001, s = 1/6, c = 0. Special points from bifurcations analyses are listed in Supplementary file 3. In fitting our theoretical model to in vitro data, we estimated the within-host virus transmission rate (b) and the rate(s) of cellular acquisition to antiviral status (r or r + ") ( Table 1 ; Supplementary file 4). Under absent immune assumptions, r and " were fixed at 0 while b was estimated; under induced immune assumptions, " was fixed at 0 while r and b were estimated; and under constitutive immune assumptions, all three parameters (r, ", and b) were simultaneously estimated for each cell-virus combination. Best fit parameter estimates for MOI=0.001 data are visualized in conjunction with br and b -" bifurcations in (r) and (B) the constitutive immunity rate of antiviral acquisition ("). Panels show variation in the extent of immunity, from absent (left) to high (right). Branch point curves are represented as solid lines and Hopf curves as dashed lines. White space indicates endemic equilibrium (persistence), gray space indicates limit cycling, and black space indicates no infection (extinction). Other parameter values for equilibrium analysis were fixed at: b = .025, m = .001, s = 1/6, a = 1/6, c = 0. Special points from bifurcations analyses are listed in Supplementary file 3. space corresponding to theoretical limit cycles, consistent with observed virus-induced epidemic extinctions in stochastic tissue cultures. In contrast to Vero cells, the induced immunity model offered the best fit to all RoNi/7.1 data, consistent with reported patterns in the literature and our own validation by qPCR ( Table 1; Arnold et al., 2018; Kuzmin et al., 2017; Biesold et al., 2011; Pavlovich et al., 2018) . As in Vero cell trials, we estimated highest b values for rVSV-G infections on RoNi/7.1 cell lines but here recovered higher b estimates for rVSV-MARV than for rVSV-EBOV. This reversal was balanced by a higher estimated rate of acquisition to antiviral status (r) for rVSV-EBOV versus rVSV-MARV. In general, we observed that more rapid rates of antiviral acquisition (either induced, r, constitutive, ", or both) correlated with higher transmission rates (b). When offset by r, b values estimated for RoNi/7.1 infections maintained the same amplitude as those estimated for immune-absent Vero cell lines but caused gentler epidemics and reduced cellular mortality (Figure 1) . RoNi/7.1 parameter estimates localized in the region corresponding to endemic equilibrium for the deterministic, theoretical model (Figure 4) , yielding less acute epidemics which nonetheless went extinct in stochastic experiments.

Question: Which cells are IFN-defective and therefore limited in  antiviral capacity?
demonstrate idiosyncratic induced interferon responses upon viral challenge