# Epidemiological Consequences of Viral Interference: A Mathematical Modeling Study of Two Interacting Viruses

^{1}School of Health Professions, Eastern Virginia Medical School, Norfolk, VA, United States^{2}Department of Physics & Astronomy, Texas Christian University, Fort Worth, TX, United States

Some viruses have the ability to block or suppress growth of other viruses when simultaneously present in the same host. This type of viral interference or viral block has been suggested as a potential interaction between some respiratory viruses including SARS-CoV-2 and other co-circulating respiratory viruses. We explore how one virus' ability to block infection with another within a single host affects spread of the viruses within a susceptible population using a compartmental epidemiological model. We find that population-level effect of viral block is a decrease in the number of people infected with the suppressed virus. This effect is most pronounced when the viruses have similar epidemiological parameters. We use the model to simulate co-circulating epidemics of SARS-CoV-2 and influenza, respiratory syncytial virus (RSV), and rhinovirus, finding that co-circulation of SARS-CoV-2 and RSV causes the most suppression of SARS-CoV-2. Paradoxically, co-circulation of SARS-CoV-2 and influenza or rhinovirus results in almost no change in the SARS-CoV-2 epidemic, but causes a shift in the timing of the influenza and rhinovirus epidemics.

## 1. Introduction

A novel coronavirus first detected in Wuhan, China in late 2019 spread rapidly around the world (Chen et al., 2020; Wu et al., 2020). While many people infected with the virus remain asymptomatic, a significant number of people develop COVID-19 (Dong et al., 2020; Verity et al., 2020)—a disease characterized by severe respiratory distress that leads to hospitalization and potentially death (Goyal et al., 2020; Jiang et al., 2020; Sun et al., 2020). The first waves of SARS-CoV-2 infection occurred in the late spring and summer months when there are few co-circulating respiratory viruses. However, subsequent waves have occurred at times when we normally expect seasonal circulation of viruses such as influenza and respiratory syncytial virus (RSV) (Choe et al., 2019; DeGroote et al., 2020). There was a fear that the combined effect of SARS-CoV-2 and illness due to these other respiratory viruses could lead to an excessive burden on the healthcare system (Kissler et al., 2020; Xu and Li, 2020), but in many parts of the world, the usual seasonal respiratory virus epidemics did not materialize (Casalegno et al., 2021; Eisen et al., 2021; Le Hingrat et al., 2021; Liu et al., 2021; Tempia et al., 2021; Williams et al., 2021). It's not clear if the suppression of other respiratory viruses is due to an interaction between SARS-CoV-2 and these other viruses or due to non-pharmaceutical interventions imposed to control the spread of SARS-CoV-2 (Oh et al., 2021; Redlberger-Fritz et al., 2021; Wagatsuma et al., 2021).

There are observational studies that have noted a lack of SARS-CoV-2 coinfections with other respiratory viruses within a single host (Blasco et al., 2020; Kim et al., 2020; Nowak et al., 2020; Xing et al., 2020; Alhumaid et al., 2021). While this could be attributed to the limited circulation of other respiratory viruses, other studies suggest that the limited number of coinfections might be due to virus-virus interactions within the host (Dee et al., 2021; Stowe et al., 2021). A recent modeling study has suggested that many respiratory viruses actually suppress SARS-CoV-2 infections (Pinky and Dobrovolny, 2020), potentially preventing patients from contracting SARS-CoV-2 while they are infected with another virus

It is unclear what this type of viral interference means for the eventual co-circulation of SARS-CoV-2 with other respiratory viruses. Studies have noted that outbreaks of other respiratory viruses influence each other's timing. For example, several studies found that circulation of the 2009 pandemic influenza virus was delayed in some regions, coinciding with an outbreak of rhinovirus, suggesting that rhinovirus infections delayed the influenza outbreak (Anestad and Nardbo, 2009, 2011; Casalegno et al., 2010a). Studies also noted that the 2009 pandemic influenza outbreak delayed outbreaks of several other respiratory viruses (Casalegno et al., 2010b; Pascalis et al., 2012; Yang et al., 2012; Sun et al., 2014), but it is not just outbreaks of novel viruses that alter outbreaks of circulating viruses. Influenza and RSV have similar seasonality yet they are never seen to peak at the same time, suggesting some kind of interference (Anestad, 1982; Velasco-Hernandez et al., 2015). Interaction between annual seasonal influenza outbreaks and rhinovirus often seems to have an effect on rhinovirus prevalence (Nickbakhsh et al., 2019) and early annual influenza A outbreaks have been observed to shift the timing of RSV, influenza B, and coronavirus outbreaks in the same year (van Asten et al., 2016). More generally, a recent study investigated covariance of a variety of respiratory virus outbreaks finding both negative and positive covariance between a number of viral pairs (Mair et al., 2019).

Epidemiological models have also been developed in an attempt to understand how viruses interact at the population level. Many of these models assume interaction between the viruses based on some form of cross-immunity (Ackleh and Allen, 2005; Allen and Kirupaharan, 2005; Nuno et al., 2005; Saunders, 1981; Andreason, 2018; Garba et al., 2013; Ackleh and Salceanu, 2014; Bhattacharyya et al., 2015; Baguelin and Eggo, 2018; Almaraz and Gomez-Corral, 2019; Amador et al., 2019; Gutierrez-Jara et al., 2019), finding competitive exclusion or at least a diminished presence of one virus depending on the strength and mathematical formulation of the cross-immunity. Other interactions between viruses have also been considered. One epidemiological model examined the interaction of two viruses (influenza and RSV) under the assumption of one virus being competitively stronger than the other, finding that interaction of the two viruses caused emergence of a biennial fluctuation in the size of epidemics (Velasco-Hernandez et al., 2015). A study of different mechanisms of interaction between two influenza strains (H1N1 and H3N2) determined that the interactions only caused observable interference when the attack rates of the strains were high (Ackerman et al., 1990).

In all these models, the cross-protective immunity or other interaction is assumed to be equivalent for both viruses and typically lasts for the duration of the epidemic, whereas interaction through within-host viral interference leads to a one-sided effect that provides only temporary protection against one of the viruses. In this article, we investigate the case of an asymmetric interaction between viruses based on the viral suppression of SARS-CoV-2 by other respiratory viruses. In our model, infection with one virus (influenza, for example) and exposure to the second (SARS-CoV-2) will temporarily prevent infection with the second virus (SARS-CoV-2), but infection with the second virus (SARS-CoV-2) and exposure to the first (influenza) will lead to coinfection. We find that this type of interaction leads to a smaller epidemic of the suppressed virus, but can also alter the timing of the epidemic of the dominant virus.

## 2. Materials and Methods

### 2.1. Mathematical Model

Spread of an infectious disease through a population is often modeled using ordinary differential equation (ODE) compartment models. The common basis for many of these epidemiological models is the Susceptible-Exposed-Infectious-Recovered (SEIR) model where each class of the population is represented by a compartment and movement of people between compartments is governed by ODEs (Hethcote, 2000). We extend the SEIR compartmental model to include two viruses,

The model diagram is given in Figure 1. Initially, the population is susceptible to both viruses (*S*). People can be infected by either virus with force of infection β_{i} (*i* = 1, 2) and move into the exposed compartments (*E*_{i}). Our assumption is that virus 1 (non-SARS-CoV-2 virus such as influenza) blocks infection with virus 2 (SARS-CoV-2) within the host, so people exposed to virus 1 transition to the infectious compartment (*I*_{1}). Upon recovery from virus 1, they are immune to virus 1, but are now susceptible to virus 2 (${S}_{2}^{(1)}$). Once infected with the second virus, they move through the exposed (${E}_{2}^{(1)}$) and infectious (${I}_{2}^{(1)}$) stages, after which they are fully recovered (*R*) and immune to both viruses. People who are infected first with virus 2 can develop a coinfection. Thus, people in both the exposed and infected compartments for the second virus (*E*_{2} and *I*_{2}) can become coinfected. We assume that coinfected exposed (*E*_{3}) people have moved to the infected phase for virus 2, so can infect others with this virus, but are still only exposed to virus 1. After some transition time, they become fully coinfected (*I*_{3}) and capable of transmitting both viruses. After some time, coinfected patients also recover and are immune to both viruses. Those patients exposed to virus 2 who do not become coinfected move to the infected compartment (*I*_{2}) and recover. At this point, they are immune to virus 2, but are still susceptible to virus 1 (${S}_{1}^{(2)}$). Once infected with virus 1, they move through the corresponding exposed (${E}_{1}^{(2)}$) and infected (${I}_{1}^{(2)}$) compartments to fully recover.

**Figure 1**. An extended SEIR model that includes two interacting viruses. Infection with virus 1 temporarily blocks infection with virus 2, but infection with virus 2 can lead to coinfection upon exposure to virus 1. The transmission rate of each virus is β_{i}. The duration of the exposed period is 1/*k*_{i} and the duration of the infectious period is 1/δ_{i}.

We make a number of simplifying assumptions to keep the model tractable. We assume that the population remains fixed, with total population size *N*, and that births and deaths are negligible over the course of the epidemic. We also do not track deaths due to the infectious diseases. For an epidemic model with a single virus, this would mean that the “recovered” compartment represents all people who are no longer able to participate in the infection, whether it is because they have recovered and are immune or because they have died. In this model, we have compartments with partial recovery where people are immune to one virus, but still susceptible to the other. If someone dies from the first infection, then they cannot be infected with the second virus, so we are over-counting the number of people with sequential infections. This might be substantial, particularly for people infected first with SARS-CoV-2, which has a fairly high death rate (Salzberger et al., 2020). (This assumption is relaxed in the Supplementary Material). We assume that order of infection does not change any of the infection parameters for a particular virus. In this model, these are the force of infection β_{i}, the incubation period 1/*k*_{i}, and the recovery rate δ_{i} (*i* = 1, 2). We do, however, explore the possibility that coinfection can lead to a longer recovery time.

### 2.2. Simulations

For simulations, we use influenza epidemic parameters (Spencer et al., 2021) as a baseline to investigate interacting viruses (Table 1). We then consider what might happen during a potential next wave of SARS-CoV-2 infection that could coincide with seasonal circulation of other respiratory viruses. We use parameters for influenza, respiratory syncytial virus (RSV), and rhinovirus (RV) derived from fits to several years of seasonal outbreaks in the US (Spencer et al., 2021). Note that, particularly for influenza, this means that we are not modeling any specific strain of influenza. For SARS-CoV-2, we use parameters derived from fits to data from the initial outbreak in China (Anderson et al., 2020; Bentout et al., 2020). Parameters are given in Table 1. This leaves us with one unknown parameter: the coinfection recovery rate δ_{3}. Since the recovery rates for SARS-CoV-2 coinfections are unknown, we consider both a faster and slower coinfection recovery rate: δ_{3}= 0.1 /d or δ_{3}= 0.01 /d. Simulations are run using python's odeint differential equation solver (Bell et al., 2022).

**Table 1**. SARS-CoV-2 (Anderson et al., 2020; Bentout et al., 2020), influenza virus (H1N1) (Spencer et al., 2021), respiratory syncytial virus (RSV) (Spencer et al., 2021), and rhinovirus (RV) (Spencer et al., 2021) parameters for simulation.

### 2.3. Sensitivity Analysis

We use the partial rank correlation coefficient (PRCC) to assess which model parameters contribute most to changes in epidemic outcome. We allow parameter values to range between ±10% of their base value and use 1,000 different randomly selected parameter combinations to calculate the PRCC. The partial rank correlation coefficient is close to ±1 if there is a high degree of correlation between the independent and dependent variables, with positive values indicating a positive correlation (both increase or decrease together) and negative values indicating a negative correlation. The output variables examined for PRCC are total number of virus 1 infections, total number of virus 2 infections, and the number of coinfections.

## 3. Results

### 3.1. Epidemics With Identical Viruses

We use a standard compartmental SEIR-type ODE model described by Equation (1). In the model, people infected with virus 1 cannot be infected by virus 2, whereas people infected with virus 2 can be coinfected. The model also allows for sequential infections without any asymmetry, i.e., a person can get the other virus after recovering from their initial infection. Infection is initiated by 100 infected individuals for each of the viruses, *I*_{i}(0) = 100. We assume an initial susceptible population, *S*(0), of 331,002,651 (the population of the US). All other classes are initially set to zero. We start by examining the interaction of two identical viruses using parameters for influenza infection (Table 1).

Figure 2A shows the projected epidemics caused by two identical viruses circulating in a population, but with virus 1 having the ability to block infection with virus 2. For comparison, the dashed black line shows the epidemic if only one of these viruses is circulating on its own. For both viruses in the co-circulating epidemic, the number of people infected with a single virus (either 1 or 2) is lower than if the virus was circulating without any interaction. The number of people infected with virus 2, however, is more suppressed than the number of people infected with virus 1. Even if we include the coinfected population, there is an overall lowering of the number of people infected with virus 2 in the population, while the total number of people infected with virus 1 stays the same (Figure 2B).

**Figure 2**. Interacting viral epidemics. Infection is initiated by 100 virus 1 and virus 2 infected individuals assuming the initial number of susceptible population is the current population of the US, i.e., 331,002,651. **(A)** The dashed line shows the trajectory of a single epidemic, while the remaining lines show the number of people infected during a co-circulating epidemic, with red giving the number of people infected with virus 1, blue giving the number of people infected with virus 2, and cyan giving the coinfected population. **(B)** The fraction of the population infected at the end of the epidemic for a single virus (black) or co-circulation of virus 1 (red) and virus 2 (blue). The left bars do not include co-infected people, while the right bars include co-infected people.

We use partial rank correlation coefficient (PRCC) to assess how different model parameters affect the outcome of the epidemic. We consider three outcomes: the total number of people infected with virus 1, the total number of people infected with virus 2, and the total number of coinfections. Results of the PRCC are shown in Figure 3. The results for the numbers of people infected with either virus 1 or virus 2 are fairly symmetric, with the total number of people infected with virus 1 depending most strongly on virus 1 infection and recovery rates (β_{1} and δ_{1}), and the total number of people infected with virus 2 depending most strongly on virus 2 infection and recovery rates (β_{2} and δ_{2}). The asymmetry in the model is most apparent when looking at the parameter dependence of the total number of coinfections. The total number of coinfections depends more strongly on the infection rate and recovery rate of the second virus than on the infection and recovery rates of the first virus. The only way to become coinfected is to first be infected with the second virus, so the more virus 2 infections and the longer they last, the more likely we are to see coinfections.

**Figure 3**. Partial rank correlation coefficients for model parameters. Bars indicate the coefficients, which measure the sensitivity of total number of virus 1 infections **(A)**, virus 2 infections **(B)**, and number of coinfections **(C)** on model parameters.

### 3.2. Epidemics With Different Viruses

In reality the viruses that circulate at the same time are described by different epidemiological parameters. To better understand how differences in viral spread enhance or suppress the population-level interaction of the two viruses, we examine the epidemic when the infection rate (β_{1} or β_{2}) of either virus changes. We use three measures of the epidemic curve to assess how the epidemic changes: the peak number of infected, the time of peak, and the cumulative number of infected.

Results of this analysis are shown in Figure 4. The left and center columns show changes in peak number of infected (top row), time of peak (center row), and total infected (bottom row) measured relative to their value in a single virus epidemic. The left column shows results for virus 1 (the dominant virus; Figures 4A,D,G) while the center column shows results for virus 2 (the suppressed virus; Figures 4B,E,H). In each figure we vary either the infection rate of virus one (β_{1}) or the infection rate of virus 2 (β_{2}). Thus, varying β_{1} increases the peak number of infected for virus 1 (red line; Figure 4A), but has little visible effect on the peak number of infected for virus 2 (red line; Figure 4B). Varying β_{2} has the opposite effect (blue lines; Figures 4A,B). Altogether, we see that varying β_{i} causes changes primarily in the epidemic of virus *i*, where there is a minimum β_{i} needed for the presence of an epidemic. Beyond this threshold, as β_{i} increases, there is a higher peak number of infections (Figures 4A,B), an earlier time of peak (Figures 4D,E) and higher total number of infections (Figures 4G,H). There is little change in the epidemic time course of one virus as the infection rate of the other virus is varied.

**Figure 4**. Co-circulating viruses with different infection rates. We change the infection rate of virus 1 (red lines) or virus 2 (blue lines), leaving other parameters the same, and plot the peak number of infected **(A–C)**, time of epidemic peak **(D–F)**, and total number of infected **(G–I)** for virus 1 **(A,D,G)** and virus 2 **(B,E,H)**. All values are presented as relative to their values for a single virus epidemic. **(C,F,I)** The graphs highlight asymmetries in response of the two viruses by subtracting blue lines from corresponding red lines in the left and center columns.

There is also very little obvious asymmetry in the graphs; changing the infection rate of the dominant and suppressed virus appears to have more or less the same outcome. However, if we look carefully, there is a small dip in virus 2 response to changes in β_{1} (near relative β = 1) that is not apparent in the response of virus 1 to changes in β_{2} of the left column. To highlight this asymmetry, we subtract the direct responses (virus *i* response to change in β_{i}) of virus 1 and virus 2, as well as subtract the indirect responses (virus *i* response to change in β_{j}) of virus 1 and virus 2, with the results given in Figures 4C,F,I. Any asymmetry in the response to changes in infection rate is confined to a small range near equal infection rates when the viruses have similar epidemiological characteristics. The dominant virus has a slightly higher peak number of infections, higher total infections and earlier time of peak than the suppressed virus, no matter which infection rate is varied (either β_{1} or β_{2}). Thus, the within host suppression of one virus by another is reflected in decreased population-level spread of the suppressed virus.

### 3.3. SARS-CoV-2 and Other Respiratory Viruses

Based on the results of a recent modeling study that suggests that SARS-CoV-2 is suppressed by several common respiratory viruses (Pinky and Dobrovolny, 2020), we investigated interacting epidemics of SARS-CoV-2 with influenza, respiratory syncytial virus (RSV), and rhinovirus (RV). We start both epidemics at the same time with the same number of infected people (100), which might not be entirely realistic. Model predictions of the epidemics are shown in Figure 5. Note that these simulations use parameters describing spread of each of these viruses in a typical epidemic season without any sort of non-pharmaceutical interventions (NPI) such as mask-wearing, social distancing, or closures/lockdowns. Thus, the predictions are more representative of what could happen in coming years rather than what occurred in 2020 or even in 2021 when some of these NPIs were at least partially in place.

**Figure 5**. Common respiratory viruses co-circulating with SARS-CoV-2. Graphs show the model predictions for interacting epidemics of SARS-CoV-2 and influenza **(A,B)**, RSV **(C,D)**, and rhinovirus **(E,F)** assuming that SARS-CoV-2 is suppressed during coinfections with each of these viruses. **(A,C,E)** assume a coinfection recovery rate of 0.1 /d while **(B,D,F)** assume a coinfection recovery rate of 0.01 /d.

In the case of influenza and rhinovirus, the influenza and RV epidemics are predicted to start largely after SARS-CoV-2 has run its course, so there is little effect of the interaction on the time course of the SARS-CoV-2 epidemic. There is, however, a noticeable effect on the time course of the influenza and RV epidemics, which are shifted, and occur earlier than if either virus was circulating alone. This shift is seen for both values of coinfection recovery rate in the case of influenza, but only for the lower recovery rate in the case of RV. The shift occurs when coinfected individuals transmit the infection for a longer period of time than singly infected individuals. While there are few coinfected individuals, the longer recovery time allows them to infect a few more people, speeding up the epidemic.

In the case of SARS-CoV-2 and RSV interacting epidemics, the epidemics are predicted to peak at nearly the same time, so there is a more visible effect of the within-host viral block. The peak of the SARS-CoV-2 epidemic is reduced and delayed as compared to the single epidemic. Additionally, there is a clear change in the slope of the SARS-CoV-2 epidemic curve during the peak of the RSV epidemic, suggesting that the presence of RSV is preventing some SARS-CoV-2 infections.

### 3.4. Extensions to the Base Model

Two possible extensions of the model are examined in the Supplementary Material, but are found to have a small effect on model predictions. The first is including death of patients, a factor that could be particularly important when considering SARS-CoV-2 and influenza, which have higher case fatality rates than many other respiratory viruses (Zhang and Zhao, 2012; Liu et al., 2020; Fan et al., 2021). Inclusion of death decreases the number of people available for a sequential secondary infection, so there is a reduction in the total number of people infected with either virus (Supplementary Figure 1). We also allowed for virus 1-infected patients to become coinfected with virus 2, assuming that the virus 2 infection would be suppressed enough to be asymptomatic, but would stimulate enough immune response to prevent further virus 2 infections for the remainder of the epidemic (Supplementary Figure 3). This assumption leads to a larger number of coinfected patients and even larger suppression of virus 2 since a subset of virus 2 infected people do not transmit the infection (Supplementary Figure 4).

## 4. Discussion

We explored an epidemic involving two viruses that interact via suppression of one virus if the two viruses occur in the same host. Such viral interference is seen in many viral coinfections (Whitakerdowling and Youngner, 1987; Kumar et al., 2018), and was hypothesized as a possible explanation for an apparent lack of SARS-CoV-2 coinfections with other respiratory viruses (Pinky and Dobrovolny, 2020). Our model predicts that the within host suppression of one virus by another manifests in a decrease in prevalence of the suppressed virus at the population level.

One of the big concerns of public health officials is the potential for a surge of sick patients that might overwhelm hospital capacity (Li et al., 2020; Moghadas et al., 2020). There has been particular concern about co-circulation of SARS-CoV-2 and other respiratory viruses, like influenza, that can also cause surges in hospitalization rates (Bertolani et al., 2018). Should both epidemics peak at the same time, hospitals might not have sufficient resources to treat all those in need. Our simulations indicate that there will be a decrease in the prevalence of the suppressed virus, which helps with the hospital capacity issue, although the effect is small and might not be enough to prevent overwhelming of healthcare systems.

Hospitalization rates will also be affected by the clinical severity of coinfections. If coinfections have a protective effect, there could be less strain on medical resources; if coinfections result in more severe clinical outcomes, hospitals could be even more strained. Some studies have found that coinfections decrease disease severity (Martin et al., 2011, 2013) or are at least no more severe than mono-infections (Brand et al., 2012; Martin et al., 2013; Asner et al., 2014; Rotzen-Ostlund et al., 2014; Mexico Emerging Infect Dis, 2019; Xiang et al., 2021), although others have found that coinfections can be more severe than mono-infections (Waner, 1994; Goka et al., 2015; Alosaimi et al., 2021; Musuuza et al., 2021). In the case of SARS-CoV-2 coinfections, studies are equally mixed about the severity of clinical disease, with some studies indicating a protective effect (Chekuri et al., 2021; Goldberg et al., 2021), others showing a worsening of clinical outcomes (Alosaimi et al., 2021; Stowe et al., 2021), and still others showing no significant difference between SARS-CoV-2 coinfections and mono-infections (Cheng et al., 2021; Guan et al., 2021). A meta-analysis of SARS-CoV-2 coinfections with influenza indicated no overall increase in mortality associated with coinfections, but found that SARS-CoV-2/influenza coinfections had decreased mortality in China and increased mortality in other regions (Guan et al., 2021), suggesting that other factors besides the characteristics of the two pathogens might be involved in determining severity.

Using parameters for influenza, RSV, and rhinovirus, we found that RSV had the most impact on the prevalence of SARS-CoV-2 when the two viruses co-circulate. Of the viruses investigated here, SARS-CoV-2 and RSV have the most similar infection rates (β). According to the analysis of section 3.2, the largest changes occur when the two viruses have similar growth rates, so this is why our model predicts the largest effect for the RSV/SARS-CoV-2 combination. In reality, all of these viruses, and more, are co-circulating together, making it difficult to isolate the interaction of one pair (de Celles et al., 2022). In our simulations, we assume that influenza, RSV, and rhinovirus all suppress SARS-CoV-2, as was suggested in a mathematical modeling study (Pinky and Dobrovolny, 2020). These modeling results have been recently supported by experiments showing suppression of SARS-CoV-2 by rhinovirus in human respiratory epithelium (Dee et al., 2021) and suppression of SARS-CoV-2 viral loads in mice first infected with influenza (Achdout et al., 2021), as well as in ferrets first infected with influenza (Bao et al., 2021). Our model could also be used to simulate the converse, SARS-CoV-2 suppression of other respiratory viruses, although for influenza and rhinovirus, we expect to see little change in the epidemic dynamics since their infection rates are quite different from the infection rate of SARS-CoV-2.

Two extensions of our interaction model, inclusion of patient death and incomplete viral block, are considered in the Supplementary Material, but other factors might also affect the predicted dynamics. For example, asymptomatic infections contribute to spread of influenza (Furuya-Kanamori et al., 2016; Cui et al., 2019), RSV (Munywoki et al., 2015), rhinovirus (Martin et al., 2018), and SARS-CoV-2 (Al-Sadeq and Nasrallah, 2020; Dobrovolny, 2020), but are not explicitly accounted for in our model. Additionally, people infected with SARS-CoV-2 are asked to isolate, keeping them from further spreading the infection (Hellewell et al., 2020). While symptomatic patients with other respiratory viruses will often avoid typical daily activities, isolation measures are not as strict, so this would add an additional asymmetry between the dynamics of SARS-CoV-2 and other respiratory viruses. There is the possibility of other, less direct, interactions between respiratory virus infections. For example, a number of studies have observed a correlation between influenza vaccination and decreased COVID-19 severity in patients (Amato et al., 2020; Marín-Hernández et al., 2021; Zanettini et al., 2020), suggesting a possible interaction between the influenza immune response and SARS-CoV-2 virus. Perhaps the biggest factor not considered in the model is mitigation strategies for reducing spread of SARS-CoV-2, such as wearing of masks and social distancing. These NPIs not only change the time course of the SARS-CoV-2 epidemic, but also help stem the spread of other respiratory viruses (Ngonghala et al., 2020) and are largely thought to be responsible for the lack of typical seasonal respiratory viruses in 2020 (Oh et al., 2021; Redlberger-Fritz et al., 2021; Wagatsuma et al., 2021).

Despite these limitations, our model provides insight into how co-circulating viruses will spread through the population if they participate in viral interference at the within-host level. The model suggests a reduction in the total number of cases of the suppressed virus, although the effect is small, but also suggests a shift in the timing of the peak of the dominant virus is possible.

## Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

## Author Contributions

LP and HD: conceptualization, methodology, validation, and writing—review and editing. LP: software and formal analysis. HD: writing—original draft preparation, supervision, and project administration. Both authors have read and agreed to the published version of the manuscript.

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.830423/full#supplementary-material

## References

Achdout, H., Vitner, E. B., Politi, B., Melamed, S., Yahalom-Ronen, Y., Tamir, H., et al. (2021). Increased lethality in influenza and SARS-CoV-2 coinfection is prevented by influenza immunity but not SARS-CoV-2 immunity. *Nat. Commun.* 12, 5819. doi: 10.1038/s41467-021-26113-1

Ackerman, E., Longini, I., Seaholm, S., and Hedin, A. (1990). Simulation of mechanisms of viral interference in influenza. *Int. J. Epidemiol.* 19, 444–454. doi: 10.1093/ije/19.2.444

Ackleh, A., and Allen, L. (2005). Competitive exclusion in SIS and SIR epidemic models with total cross immunity and density-dependent host mortality. *Discr. Cont. Dyn. Syst. B* 5, 175–188. doi: 10.3934/dcdsb.2005.5.175

Ackleh, A. S., and Salceanu, P. L. (2014). Robust uniform persistence and competitive exclusion in a nonautonomous multi-strain SIR epidemic model with disease-induced mortality. *J. Math. Biol.* 68, 453–475. doi: 10.1007/s00285-012-0636-4

Alhumaid, S., Al Mutair, A., Al Alawi, Z., Alshawi, A. M., Alomran, S. A., Almuhanna, M. S., et al. (2021). Coinfections with bacteria, fungi, and respiratory viruses in patients with SARS-CoV-2: a systematic review and meta-analysis. *Pathogens* 10, 809. doi: 10.3390/pathogens10070809

Allen, L. J., and Kirupaharan, N. (2005). Asymptotic dynamics of deterministic and stochastic epidemic models with multiple pathogens. *Int. J. Numer. Anal. Model.* 2, 329–344.

Almaraz, E., and Gomez-Corral, A. (2019). Number of infections suffered by a focal individual in a two-strain SIS model with partial cross-immunity. *Math. Methods Appl. Sci.* 42, 4318–4330. doi: 10.1002/mma.5652

Alosaimi, B., Naeem, A., Hamed, M. E., Alkadi, H. S., Alanazi, T., Al Rehily, S. S., et al. (2021). Influenza co-infection associated with severity and mortality in COVID-19 patients. *Virol. J.* 18, 127. doi: 10.1186/s12985-021-01594-0

Al-Sadeq, D. W., and Nasrallah, G. K. (2020). The incidence of the novel coronavirus SARS-CoV-2 among asymptomatic patients: a systematic review. *Int. J. Infect. Dis.* 98, 372–380. doi: 10.1016/j.ijid.2020.06.098

Amador, J., Armesto, D., and Gomez-Corral, A. (2019). Extreme values in SIR epidemic models with two strains and cross-immunity. *Math. Biosci. Eng.* 16, 1992–2022. doi: 10.3934/mbe.2019098

Amato, M., Werba, J. P., Frigerio, B., Coggi, D., Sansaro, D., Ravani, A., et al. (2020). Relationship between influenza vaccination coverage rate and COVID-19 outbreak: an italian ecological study. *Vaccines* 8, 535. doi: 10.3390/vaccines8030535

Anderson, R. M., Heesterbeek, H., Klinkenberg, D., and Hollingsworth, T. D. (2020). How will country-based mitigation measures influence the course of the covid-19 epidemic? *Lancet* 395, 931–934. doi: 10.1016/S0140-6736(20)30567-5

Andreason, V. (2018). Epidemics in competition: partial cross-immunity. *Bull. Math. Biol.* 80, 2957–2977. doi: 10.1007/s11538-018-0495-2

Anestad, G. (1982). Interference between outbreaks of respiratory syncytial virus and influenza virus infection. *Lancet* 1:502. doi: 10.1016/s0140-6736(82)91466-0

Anestad, G., and Nardbo, S. A. (2009). Interference between outbreaks of respiratory viruses. 1*Euro Surveil*. 14, 19359. doi: 10.2807/ese.14.41.19359-en

Anestad, G., and Nardbo, S. A. (2011). Virus interference. Did rhinoviruses activity hamper the progress of the 2009 influenza A (H1N1) pandemic in Norway? *Med. Hypoth.* 77, 1132–1134. doi: 10.1016/j.mehy.2011.09.021

Asner, S. A., Science, M. E., Tran, D., Smieja, M., Merglen, A., and Mertz, D. (2014). Clinical disease severity of respiratory viral co-infection versus single viral infection: a systematic review and meta-analysis. *PLoS ONE* 9, e99392. doi: 10.1371/journal.pone.0099392

Baguelin, L. O. M., and Eggo, R. M. (2018). Influenza interaction with cocirculating pathogens and its impact on surveillance, pathogenesis, and epidemic profile: a key role for mathematical modelling. *PLoS Pathog.* 14, e1006770. doi: 10.1371/journal.ppat.1006770

Bao, L., Deng, W., Qi, F., Lv, Q., Song, Z., Liu, J., et al. (2021). Sequential infection with H1N1 and SARS-CoV-2 aggravated COVID-19 pathogenesis in a mammalian model, and co-vaccination as an effective method of prevention of COVID-19 and influenza. *Signal Transduct. Target. Ther.* 6:200. doi: 10.1038/s41392-021-00618-z

Bell, P., Burovski, E., Charlong, J., Gommers, R., Picus, M., Reddy, T., et al. (2022). *Scipy Version 1.7.1*.

Bentout, S., Chekroun, A., and Kuniya, T. (2020). Parameter estimation and prediction for coronavirus disease outbreak 2019 (COVID-19) in Algeria. *AIMS Publ. Health* 7, 306–318. doi: 10.3934/publichealth.2020026

Bertolani, A., Fattore, G., and Pregliasco, F. (2018). The hospitalization burden of influenza: just the tip of the iceberg? *Glob. Region. Health Tech. Assess.* doi: 10.1177/2284240318777148

Bhattacharyya, S., Gesteland, P. H., Korgenski, K., Bjornstad, O. N., and Adler, F. R. (2015). Cross-immunity between strains explains the dynamical pattern of paramyxoviruses. *Proc. Natl. Acad. Sci. U.S.A.* 112, 13396–13400. doi: 10.1073/pnas.1516698112

Blasco, M. L., Buesa, J., Colomina, J., Forner, M. J., Galindo, M. J., Navarro, J., et al. (2020). Co-detection of respiratory pathogens in patients hospitalized with coronavirus viral disease—2019 pneumonia. *J. Med. Virol.* 92, 1799–1801. doi: 10.1002/jmv.25922

Brand, H. K., de Groot, R., Galama, J. M., Brouwer, M. L., Teuwen, K., Hermans, P. W., et al. (2012). Infection with multiple viruses is not associated with increased disease severity in children with bronchiolitis. *Pediatr. Pulmonol.* 47, 393–400. doi: 10.1002/ppul.21552

Casalegno, J., Ottmann, M., Duchamp, M. B., Escuret, V., Billaud, G., Frobert, E., et al. (2010a). Rhinoviruses delayed the circulation of the pandemic influenza A (H1N1) 2009 virus in France. *Clin. Microbiol. Infect.* 16, 326–329. doi: 10.1111/j.1469-0691.2010.03167.x

Casalegno, J., Ottmann, M., Duchamp, M. B., Valette, M., Morfin, F., and Lina, B. (2010b). Impact of the 2009 influenza A(H1N1) pandemic wave on the pattern of hibernal respiratory virus epidemics, France, 2009. *Euro Surveil.* 15:19485. doi: 10.2807/ese.15.06.19485-en

Casalegno, J.-S., Ploin, D., Cantais, A., Masson, E., Bard, E., Valette, M., et al. (2021). Characteristics of the delayed respiratory syncytial virus epidemic, 2020/2021, Rhone Loire, France. *Euro Surveil.* 26, 2100630. doi: 10.2807/1560-7917.ES.2021.26.29.2100630

Chekuri, S., Szymczak, W. A., Goldstein, D. Y., Nori, P., Rolon, R. M., Spund, B., et al. (2021). SARS-CoV-2 coinfection with additional respiratory virus does not predict severe disease: a retrospective cohort study. *J. Antimicrob. Chemother.* 76, 12–19. doi: 10.1093/jac/dkab244

Chen, N., Zhou, M., Dong, X., Qu, J., Gong, F., Han, Y., et al. (2020). Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. *Lancet* 395, 507–513. doi: 10.1016/S0140-6736(20)30211-7

Cheng, Y., Ma, J., Wang, H., Wang, X., Hu, Z., Li, H., et al. (2021). Co-infection of influenza A virus and SARS-CoV-2: a retrospective cohort study. *J. Med. Virol.* 93, 2947–2954. doi: 10.1002/jmv.26817

Choe, Y. J., Smit, M. A., and Mermel, L. A. (2019). Seasonality of respiratory viruses and bacterial pathogens. *Antimicrob. Resist. Infect. Control* 8, 125. doi: 10.1186/s13756-019-0574-7

Cui, J., Zhang, Y., Feng, Z., Guo, S., and Zhang, Y. (2019). Influence of asymptomatic infections for the effectiveness of facemasks during pandemic influenza. *Math. Biosci. Eng.* 16, 3936–3946. doi: 10.3934/mbe.2019194

de Celles, M. D., Goult, E., Casalegno, J.-S., and Kramer, S. C. (2022). The pitfalls of inferring virus-virus interactions from co-detection prevalence data: application to influenza and SARS-CoV-2. *Proc. R. Soc. B* 289, 20212358. doi: 10.1098/rspb.2021.2358

Dee, K., Goldfarb, D. M., Haney, J., Amat, J. A. R., Herder, V., Stewart, M., et al. (2021). Human rhinovirus infection blocks severe acute respiratory syndrome coronavirus 2 replication within the respiratory epithelium: implications for COVID-19 epidemiology. *J. Infect. Dis.* 224, 31–38. doi: 10.1093/infdis/jiab147

DeGroote, N. P., Haynes, A. K., Taylor, C., Killerby, M. E., Dahl, R. M., Mustaquim, D., et al. (2020). Human parainfluenza virus circulation, United States, 2011-2019. *J. Clin. Virol.* 124, 104261. doi: 10.1016/j.jcv.2020.104261

Dobrovolny, H. M. (2020). Modeling the role of asymptomatics in infection spread with application to SARS-CoV-2. *PLoS ONE* 15, e0236976. doi: 10.1371/journal.pone.0236976

Dong, Y., Zhou, H., Li, M., Zhang, Z., Guo, W., Yu, T., et al. (2020). A novel simple scoring model for predicting severity of patients with SARS-CoV-2 infection. *Transbound. Emerg. Dis.* 67, 2823–2829. doi: 10.1111/tbed.13651

Eisen, A. K. A., Gularte, J. S., Demoliner, M., de Abreu Goes Pereira, V. M., Heldt, F. H., Filippi, M., et al. (2021). Low circulation of influenza A and coinfection with SARS-CoV-2 among other respiratory viruses during the COVID-19 pandemic in a region of southern Brazil. *J. Med. Virol.* 93, 4392–4398. doi: 10.1002/jmv.26975

Fan, G., Yang, Z., Lin, Q., Zhao, S., Yang, L., and He, D. (2021). Decreased case fatality rate of COVID-19 in the second wave: a study in 53 countries or regions. *Transbound. Emerg. Dis.* 68, 213–215. doi: 10.1111/tbed.13819

Furuya-Kanamori, L., Cox, M., Milinovich, G. J., Magalhaes, R. J. S., Mackay, I. M., and Yakob, L. (2016). Heterogeneous and dynamic prevalence of asymptomatic influenza virus infections. *Emerg. Infect. Dis.* 22, 1052–1056. doi: 10.3201/eid2206.151080

Garba, S., Safi, M., and Gumel, A. (2013). Cross-immunity-induced backward bifurcation for a model of transmission dynamics of two strains of influenza. *Nonlin. Anal. Real World Appl.* 14, 1384–1403. doi: 10.1016/j.nonrwa.2012.10.003

Goka, E. A., Vallely, P. J., Mutton, K. J., and Klapper, P. E. (2015). Single, dual and multiple respiratory virus infections and risk of hospitalization and mortality. *Epidemiol. Infect.* 143, 37–47. doi: 10.1017/S0950268814000302

Goldberg, E., MHasegawa, K., Lawrence, A., Kline, J. A., and Camargo, C. A. (2021). Viral coinfection is associated with improved outcomes in emergency department patients with SARS-CoV-2. *Western J. Emerg. Med.* 22, 1262–1269. doi: 10.5811/westjem.2021.8.53590

Goyal, P., Choi, J. J., Pinheiro, L. C., Schenck, E. J., Chen, R., Jabri, A., et al. (2020). Clinical characteristics of Covid-19 in New York City. *N. Engl. J. Med.* 382, 2372–2374. doi: 10.1056/NEJMc2010419

Guan, Z., Chen, C., Li, Y., Yan, D., Zhang, X., Jiang, D., et al. (2021). Impact of coinfection with SARS-CoV-2 and influenza on disease severity: a systematic review and meta-analysis. *Front. Publ. Health* 9, 773130. doi: 10.3389/fpubh.2021.773130

Gutierrez-Jara, J., Cordova-Lepe, F., and Munoz-Quezada, M. (2019). Dynamics between infectious diseases with two susceptibility conditions: a mathematical model. *Math. Biosci.* 309, 66–77. doi: 10.1016/j.mbs.2019.01.005

Hellewell, J., Abbott, S., Gimma, A., Bosse, N., Jarvis, C., Russell, T. W., et al. (2020). Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. *Lancet Glob. Health* 8, E488–E496. doi: 10.1016/S2214-109X(20)30074-7

Hethcote, H. W. (2000). The mathematics of infectious diseases. *SIAM Rev.* 42, 599–653. doi: 10.1137/S0036144500371907

Jiang, F., Deng, L., Zhang, L., Cai, Y., Cheung, C. W., and Xia, Z. (2020). Review of the clinical characteristics of coronavirus disease 2019 (COVID-19). *J. Gen. Intern. Med.* 35, 1545–1549. doi: 10.1007/s11606-020-05762-w

Kim, D., Quinn, J., Pinsky, B., Shah, N. H., and Brown, I. (2020). Rates of co-infection between SARS-CoV-2 and other respiratory pathogens. *JAMA* 323, 2085–2086. doi: 10.1001/jama.2020.6266

Kissler, S. M., Tedijanto, C., Goldstein, E., Grad, Y. H., and Lipsitch, M. (2020). Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period. *Science* 368, 860–868. doi: 10.1126/science.abb5793

Kumar, N., Sharma, S., Barua, S., Tripathi, B. N., and Rouse, B. T. (2018). Virological and immunological outcomes of coinfections. *Clin. Microbiol. Rev.* 31, e00111-17. doi: 10.1128/CMR.00111-17

Le Hingrat, Q., Bouzid, D., Choquet, C., Laurent, O., Lescure, F.-X., Timsit, J.-F., et al. (2021). Viral epidemiology and SARS-CoV-2 co-infections with other respiratory viruses during the first COVID-19 wave in Paris, France. *Influenza Other Resp. Viruses* 15, 425–428. doi: 10.1111/irv.12853

Li, R., Pei, S., Chen, B., Song, Y., Zhang, T., Yang, W., et al. (2020). Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV2). *Science* 368, 489–493. doi: 10.1126/science.abb3221

Liu, J.-Q., Xu, J.-W., Sun, C.-Y., Wang, J.-N., Wang, X.-T., Chen, X., et al. (2020). Age-stratified analysis of SARS-CoV-2 infection and case fatality rate in China, Italy, and South Korea. *Eur. Rev. Med. Pharmacol. Sci.* 24, 12575–12578. doi: 10.26355/eurrev_202012_24054

Liu, P., Xu, M., Cao, L., Su, L., Lu, L., Dong, N., et al. (2021). Impact of COVID-19 pandemic on the prevalence of respiratory viruses in children with lower respiratory tract infections in China. *Virol. J.* 18, 159. doi: 10.1186/s12985-021-01627-8

Mair, C., Nickbakhsh, S., Reeve, R., McMenamin, J., Reynolds, A., Gunson, R. N., et al. (2019). Estimation of temporal covariances in pathogen dynamics using Bayesian multivariate autoregressive models. *PLoS Comput. Biol.* 15, e1007492. doi: 10.1371/journal.pcbi.1007492

Marín-Hernández, D., Schwartz, R. E., and Nixon, D. F. (2021). Epidemiological evidence for association between higher influenza vaccine uptake in the elderly and lower COVID-19 deaths in Italy. *J. Med. Virol.* 93, 64–65. doi: 10.1002/jmv.26120

Martin, E. T., Fairchok, M. P., Stednick, Z. J., Kuypers, J., and Englund, J. A. (2013). Epidemiology of multiple respiratory viruses in childcare attendees. *J. Infect. Dis.* 207, 982–989. doi: 10.1093/infdis/jis934

Martin, E. T., Kuypers, J., Chu, H. Y., Foote, S., Hashikawa, A., Fairchok, M. P., et al. (2018). Heterotypic infection and spread of rhinovirus A, B, and C among childcare attendees. *J. Infect. Dis.* 218, 848–855. doi: 10.1093/infdis/jiy232

Martin, E. T., Kuypers, J., Wald, A., and Englund, J. A. (2011). Multiple versus single virus respiratory infections: viral load and clinical disease severity in hospitalized children. *Influenza Other Respir. Viruses* 6, 71–77. doi: 10.1111/j.1750-2659.2011.00265.x

Moghadas, S. M., Shoukat, A., Fitzpatrick, M. C., Wells, C. R., Sah, P., Pandey, A., et al. (2020). Projecting hospital utilization during the COVID-19 outbreaks in the United States. *Proc. Natl. Acad. Sci. U.S.A.* 117, 9122–9126. doi: 10.1073/pnas.2004064117

Munywoki, P. K., Koech, D. C., Agoti, C. N., Bett, A., Cane, P. A., Medley, G. F., et al. (2015). Frequent asymptomatic respiratory syncytial virus infections during an epidemic in a rural Kenyan household cohort. *J. Infect. Dis.* 212, 1711–1718. doi: 10.1093/infdis/jiv263

Musuuza, J. S., Watson, L., Parmasad, V., Putman-Buehler, N., Christensen, L., and Safdar, N. (2021). Prevalence and outcomes of co-infection and superinfection with SARS-CoV-2 and other pathogens: a systematic review and meta-analysis. *PLoS ONE* 16:e0251170. doi: 10.1371/journal.pone.0251170

Ngonghala, C. N., Iboi, E., Eikenberry, S., Scotch, M., MacIntyre, C. R., Bonds, M. H., et al. (2020). Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel coronavirus. *Math. Biosci.* 325, 108364. doi: 10.1016/j.mbs.2020.108364

Nickbakhsh, S., Mair, C., Matthews, L., Reeve, R., Johnson, P. C., Thorburn, F., et al. (2019). Virus-virus interactions impact the population dynamics of influenza and the common cold. *Proc. Natl. Acad. Sci. U.S.A.* 116, 27142–27150. doi: 10.1073/pnas.1911083116

Nowak, M. D., Sordillo, E. M., Gitman, M. R., and Paniz, A. E. (2020). Co-infection in SARS-CoV-2 infected patients: where are influenza virus and rhinovirus/enterovirus? *J. Med. Virol.* 92, 1699–1700. doi: 10.1002/jmv.25953

Noyola, D. E., Hunsberger, S., Valdés Salgado, R., Powers, J. H. III., Galindo-Fraga, A., Ortiz-Hernández, A. A., et al. (2019). Comparison of rates of hospitalization between single and dual virus detection in a mexican cohort of children and adults with influenza-like illness. *Open Forum Infect. Dis.* 6, ofz424. doi: 10.1093/ofid/ofz424

Nuno, M., Feng, Z., Martcheva, M., and Castillo-Chavez, C. (2005). Dynamics of two-strain influenza with isolation and partial cross-immunity. *SIAM J. Appl. Math.* 65, 964–982. doi: 10.1137/S003613990343882X

Oh, D.-Y., Buda, S., Biere, B., Reiche, J., Schlosser, F., Duwe, S., et al. (2021). Trends in respiratory virus circulation following COVID-19-targeted nonpharmaceutical interventions in Germany, January — September 2020: analysis of national surveillance data. *Lancet Reg. Health Eur.* 6, 100112. doi: 10.1016/j.lanepe.2021.100112

Pascalis, H., Temmam, S., Turpin, M., Rollot, O., Flahault, A., Carrat, F., et al. (2012). Intense co-circulation of non-influenza respiratory viruses during the first wave of pandemic influenza pH1N1/2009: a cohort study in Reunion Island. *PLoS ONE* 7, e44755. doi: 10.1371/journal.pone.0044755

Pinky, L., and Dobrovolny, H. M. (2020). SARS-CoV-2 coinfections: could influenza and the common cold be beneficial? *J. Med. Virol.* 92, 2623–2630. doi: 10.1016/j.jtbi.2019.01.011

Redlberger-Fritz, M., Kundi, M., Aberle, S. W., and Puchhammer-Stockl, E. (2021). Significant impact of nationwide SARS-CoV-2 lockdown measures on the circulation of other respiratory virus infections in Austria. *J. Clin. Virol.* 137, 104795. doi: 10.1016/j.jcv.2021.104795

Rotzen-Ostlund, M., Eriksson, M., Lindell, A. T., Allander, T., Wirgart, B. Z., and Grillner, L. (2014). Children with multiple viral respiratory infections are older than those with single viruses. *Acta Paediatr.* 103, 100–104. doi: 10.1111/apa.12440

Salzberger, B., Buder, F., Lampl, B., Ehrenstein, B., Hitzenbichler, F., and Hanses, F. (2020). Epidemiology of SARS-CoV-2 infection and COVID-19. *Internist*. 61, 782–788. doi: 10.1007/s00108-020-00834-9

Spencer, J. A., Shutt, D. P., Moser, S. K., Clegg, H., Wearing, H. J., Mukundan, H., et al. (2021). Distinguishing viruses responsible for influenza-like illness. *medRxiv*. doi: 10.1101/2020.02.04.20020404

Stowe, J., Tessier, E., Zhao, H., Guy, R., Muller-Pebody, B., Zambon, M., et al. (2021). Interactions between SARS-CoV-2 and influenza, and the impact of coinfection on disease severity: a test-negative design. *Int. J. Epidemiol.* 50, 1124–1133. doi: 10.1093/ije/dyab081

Sun, L., Shen, L., Fan, J., Gu, F., Hu, M., An, Y., et al. (2020). Clinical features of patients with coronavirus disease 2019 from a designated hospital in Beijing, China. *J. Med. Virol.* 92, 2055–2066 doi: 10.1002/jmv.25966

Sun, Y., Zhu, R., Zhao, L., Deng, J., Wang, F., Ding, Y., et al. (2014). Effect of human rhinovirus infection in pediatric patients with influenza-like illness on the 2009 pandemic influenza A(H1N1) virus. *Chin. Med. J.* 127, 1656–1660. doi: 10.3760/cma.j.issn.0366-6999.20132386

Tempia, S., Walaza, S., Bhiman, J. N., McMorrow, M. L., Moyes, J., Mkhencele, T., et al. (2021). Decline of influenza and respiratory syncytial virus detection in facility-based surveillance during the COVID-19 pandemic, South Africa, January to October 2020. *Euro Surveil.* 26, 200160. doi: 10.2807/1560-7917.ES.2021.26.29.200160

van Asten, L., Bijkerk, P., Fanoy, E., van Ginkel, A., Suijkerbuijk, A., van der Hoek, W., et al. (2016). Early occurrence of influenza A epidemics coincided with changes in occurrence of other respiratory virus infections. *Influenza Other Resp. Dis.* 10, 14–26. doi: 10.1111/irv.12348

Velasco-Hernandez, J. X., Nunez-Lopez, M., Comas-Garcia, A., Noyola Cherpitel, D. E., and Capistran Ocampo, M. (2015). Superinfection between influenza and RSV alternating patterns in San Luis Potosi State, Mexico. *PLoS ONE* 10, e0115674. doi: 10.1371/journal.pone.0115674

Verity, R., Okell, L. C., Dorigatti, I., Winskill, P., Whittaker, C., Imai, N., et al. (2020). Estimates of the severity of coronavirus disease 2019: a model-based analysis. *Lancet Inf. Dis.* 20, 669–677. doi: 10.1016/S1473-3099(20)30243-7

Wagatsuma, K., Koolhof, I. S., Shobugawa, Y., and Saito, R. (2021). Decreased human respiratory syncytial virus activity during the COVID-19 pandemic in Japan: an ecological time-series analysis. *BMC Infect. Dis.* 21, 734. doi: 10.1186/s12879-021-06461-5

Waner, J. L. (1994). Mixed viral infections: detection and management. *Clin. Microbiol. Rev.* 7, 143–151.

Whitakerdowling, P., and Youngner, J. (1987). Viral interference-dominance of mutant viruses over wild-type virus in mixed infections. *Microbiol. Rev.* 51, 179–191. doi: 10.1128/MMBR.51.2.179-191.1987

Williams, T. C., Sinha, I., Barr, I. G., and Zambon, M. (2021). Transmission of paediatric respiratory syncytial virus and influenza in the wake of the COVID-19 pandemic. *Euro Surveil.* 26, 2100186. doi: 10.2807/1560-7917.ES.2021.26.29.2100186

Wu, F., Zhao, S., Yu, B., Chen, Y.-M., Wang, W., Song, Z.-G., et al. (2020). A new coronavirus associated with human respiratory disease in China. *Nature* 579, 265–271. doi: 10.1038/s41586-020-2008-3

Xiang, X., Wang, Z., Ye, L., He, X., Wei, X., Ma, Y., et al. (2021). Co-infection of SARS-COV-2 and influenza a virus: a case series and fast review. *Curr. Med. Sci.* 41, 51–57. doi: 10.1007/s11596-021-2317-2

Xing, Q., Li, G-j., Xing, Y-h., Chen, T., Li, W-j., Ni, W., et al. (2020). Precautions are needed for COVID-19 patients with coinfection of common respiratory pathogens. *Lancet*. doi: 10.2139/ssrn.3550013

Xu, S., and Li, Y. (2020). Beware of the second wave of COVID-19. *Lancet* 395, 1321–1322. doi: 10.1016/S0140-6736(20)30845-X

Yang, Y., Wang, Z., Ren, L., Wang, W., Vernet, G., Paranhos-Baccala, G., et al. (2012). Influenza a/H1N1 2009 pandemic and respiratory virus infections, Beijing, 2009-2010. *PLoS ONE* 7, e45807. doi: 10.1371/journal.pone.0045807

Zanettini, C., Omar, M., Dinalankara, W., Imada, E. L., Colantuoni, E., Parmigiani, G., et al. (2020). Influenza vaccination and COVID-19 mortality in the USA. *MedRxiv*. 9:427. doi: 10.1101/2020.06.24.20129817

Keywords: SARS-CoV-2, influenza, respiratory syncytial virus, rhinovirus, viral interference, mathematical model

Citation: Pinky L and Dobrovolny HM (2022) Epidemiological Consequences of Viral Interference: A Mathematical Modeling Study of Two Interacting Viruses. *Front. Microbiol.* 13:830423. doi: 10.3389/fmicb.2022.830423

Received: 07 December 2021; Accepted: 14 February 2022;

Published: 11 March 2022.

Edited by:

Jianwei Wang, Chinese Academy of Medical Sciences and Peking Union Medical College, ChinaReviewed by:

Jacques Izopet, Université Toulouse III Paul Sabatier, FranceXavier López-Labrador, Fundación Para el Fomento de la Investigación Sanitaria y Biomédica de la Comunitat Valenciana (FISABIO), Spain

Copyright © 2022 Pinky and Dobrovolny. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hana M. Dobrovolny, h.dobrovolny@tcu.edu