Self-inhibiting percolation and viral spreading in epithelial tissue

  1. Xiaochan Xu
  2. Bjarke Frost Nielsen
  3. Kim Sneppen  Is a corresponding author
  1. Niels Bohr Institute, University of Copenhagen, Denmark
  2. Novo Nordisk Foundation Center for Stem Cell Medicine, reNEW, University of Copenhagen, Denmark
  3. PandemiX Center, Department of Science and Environment, Roskilde University, Denmark
  4. High Meadows Environmental Institute, Princeton University, United States
5 figures and 1 additional file

Figures

Figure 1 with 2 supplements
Cell states during SARS-CoV-2 infection in human tracheal/bronchial epithelial cells.

(a) 6162 cells (Fiege et al., 2021) covering samples of mock-infected (0 h), 24 hpi (hours post-infection), and 48 hpi visualized with UMAP. (b) Average expression of representative viral genes, …

Figure 1—figure supplement 1
Expression of antiviral genes at different stages of early SARS-CoV-2 infection.

Once the infection has proceeded for 2 days, cells in the A state express all the antiviral genes highly. The typical a state can be identified from the 0 hr panels where some cells express low to …

Figure 1—figure supplement 2
Detected viral genes and expression of IFNs at different stages of early SARS-CoV-2 infection.

Once the infection has proceeded for 2 days, viral genes (cov.*) are detected at a high level in the V cells. The IFN genes are mainly expressed in N cells, amounting to 0.6% at 24hr and to 4.8% …

Figure 2 with 1 supplement
NOVAa model.

(a) The cell state transitions are included in the NOVAa model. The straight black arrows indicate transitions between cell states. The curved yellow arrows indicate the effects of IFNs on activating …

Figure 2—figure supplement 1
Stochastic conversion.

Average final-state cell fractions in a variant of the model with stochastic conversion, i.e. in which each conversion of a neighboring cell takes place with a probability pconv<1. In this figure, pconv=0.5. Note …

Figure 3 with 1 supplement
Cluster size distribution.

(a) The distribution P(s) of cluster sizes of infected cells (s=(N+V)/L2) for different values of pa, simulated by starting with one infected cell in a 2D square lattice of linear extent L=2000. (b) The exponents …

Figure 3—figure supplement 1
Cell fractions within the different states over times, at a range of pa values (below the critical value pc), for two interferon spreading radii, R=1 and R=5.
Figure 4 with 2 supplements
Range of IFN.

(a) Typical cluster for an R=1 simulation at pa~pc=0.278. (b) The dependence of pc with R, approximately reproduced by a fit pc~3R. For comparison, the OVA model as well as percolation has pc~3R/2. In all cases, when …

Figure 4—figure supplement 1
We introduced the simplified OVA model, which does not include an N state, but rather gives each susceptible cell a probability of spontaneously converting to an antiviral state in each time step.

OVA Model: In each time step, L2 updates are made. At each update, one of the L×L sites is randomly selected. If this site is in the O state, nothing happens. If the site is V state then its …

Figure 4—figure supplement 2
Comparison of the Gaussian Model and the main NOVAa model at identical parameters.
Author response image 1

Additional files

Download links