Next Article in Journal
Bayesian Analysis of Spatial Model for Frequency of Tornadoes
Next Article in Special Issue
Brochosomes and Other Bioaerosols in the Surface Layer of the Atmosphere of Moscow Metropolis
Previous Article in Journal
Aerosols as Vectors for Contaminants: A Perspective Based on Outdoor Aerosol Data from Kuwait
Previous Article in Special Issue
Composition and Concentration of the Biogenic Components of the Aerosols Collected over Vasyugan Marshes and Karakan Pine Forest at Altitudes from 500 to 7000 m
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Assessment of Airborne Transmission of Human and Animal Influenza Viruses in the Ferret Model

1
Federal Budgetary Research Institution State Research Center of Virology and Biotechnology “Vector” Rospotrebnadzor, Koltsovo 630559, Russia
2
Joint Russian-Vietnamese Tropical Research and Technology Centre, Nguyen Van Huyen, Nghia Do, Cau Giay, Hanoi 650000, Vietnam
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Atmosphere 2023, 14(3), 471; https://doi.org/10.3390/atmos14030471
Submission received: 27 January 2023 / Revised: 18 February 2023 / Accepted: 25 February 2023 / Published: 27 February 2023
(This article belongs to the Special Issue Bioaerosol Composition and Measurement)

Abstract

:
The WHO has developed a tool to assess the risk of newly emerging influenza viruses with pandemic potential (TIPRA). According to TIPRA, the main parameters for assessing the risk of human-to-human transmission of a novel influenza virus are its ability to bind to human cell receptors of the upper respiratory tract (URT) and transmit in model animals. The aim of this study was to quantify airborne transmission of human and animal influenza viruses in the ferret model. The transmission of influenza viruses was studied in the ferret model in an aerobiology chamber. Airborne particles concentration and fractional composition in the aerobiology chamber were measured using an aerosol particle counter and analytical aerosol filters. Viral load in ferret nasal washings and aerosol filters was determined by titration in MDCK cells and quantitative RT-PCR. Genetic analysis of influenza viruses was performed using virus genome sequences obtained by NGS. After intranasal infection, human and animal influenza viruses replicated in the cells of nasal mucosa in ferrets. The level of virus airborne particles contamination provided by infected animals depends on the infectious dose and differs significantly between influenza virus strains. The studied avian influenza viruses show insufficient transmission in the ferret model, while human and swine influenza viruses are highly transmitted in ferrets. We propose a quantitative model of airborne transmission of influenza virus from donor to recipient ferrets. Level of influenza virus transmission in the ferret model correlates with genetic markers of virus receptor specificity and the level of virus airborne particle contamination induced by donor ferrets.

1. Introduction

Influenza is a severe infectious disease and a serious problem for public health. The virus infects up to 15% of the world population annually and causes thousands of deaths worldwide. Influenza pandemics arise every 10–50 years; each pandemic is caused by a new virus strain that can transmit in the human population and has no immunity in humans [1].
Avian influenza viruses and mammalian influenza viruses such as swine influenza are a serious threat, since they can cause a new pandemic. Some avian influenza viruses have already crossed the interspecies barrier and can transmit from birds to humans. Despite the absence of registered cases of mass human-to-human transmission of these strains to date, there is a risk of emergence of new avian influenza viruses that can spread in the human population. Quantification of virus transmission in an animal model is an important component in assessing the pandemic risk of newly emerging avian influenza viruses [2].
Various kinetic [3,4] and stochastic [5,6] models of the epidemiological process are used to analyze the pandemic potential of new viruses. Dose–response relationship and duration of the pathogenic effect are considered input data for quantification of viral and microbial risk [7,8,9].
In 2016, the WHO launched the Tool for Influenza Pandemic Risk Assessment (TIPRA) and reported that virus transmission in animal models is one of the most important elements in evaluating the general risk of emergence of an influenza virus with pandemic potential [10]. Influenza A viruses are found in a wide range of host animals. The virus is often transmitted between one and different animal species, which is accompanied by genetic mutations and rearrangements [11]. These genetic changes can lead to the emergence of a virus with high human-to-human transmissibility. The emergence of a pandemic influenza virus in 2009 in addition to sporadic transmission of some influenza viruses from animals to humans highlight the necessity to monitor and assess the potential risk of occurrence of new influenza viruses that can cause pandemics. Risk assessment is aimed at determining the probability of emergence of an influenza virus with pandemic potential [12]. TIPRA recommends evaluating four viral properties for risk assessment: virus ability to bind to human cell receptors, virus genome features, ability to transmit in animal models, and sensitivity to antiviral treatment.
The basis for influenza virus transmissibility is its interaction with the host cell surface. This interaction is provided by two virus envelope proteins: hemagglutinin (HA) and neuraminidase (NA). HA binds to this cell surface receptor sialic acid and acts as one of the major factors providing species specificity of influenza virus. Avian viruses usually bind to alpha-2′3-linked sialic acid, while human viruses predominantly bind to alpha-2′6-linked sialosides [13]. NA is an enzyme cleaving sialic acids.
This cleavage facilitates virus release from the cell after infection by preventing aggregation of virions [14,15] and virus release from mucins rich in sialic acid [16,17]. The upper respiratory tract (URT) is the primary infection site in humans. The virus infects cells predominantly expressing alpha-2,6-linked sialic acids [18], while URT mucin in humans is rich in sialosides of avian type (alpha-2,3). The balance is observed between receptor-binding ability of HA and receptor cleavage activity of NA. In order to avoid inhibition by mucins, human influenza virus must have HA with low avidity to alpha-2′6-linked sialic acid and NA with high enzymatic activity against alpha-2′3-linked sialic acid. Human influenza virus must possess high ability to bind alpha-2′6-linked saccharides for effective attachment and penetration into the target cell, while effective cleavage of 2′6-linked saccharides by NA is required for virus release and prevention of virus aggregation after release from the cell. Thus, it is reasonable to assume that compatible levels of HA and NA activity are required for effective virus replication and transmission. Such studies are of strong interest for assessment of pandemic potential of influenza viruses. The role of balance of the activities of these two proteins in receptor binding was first assessed by using a new biophysical approach based on bio-layer interferometry in the real-time mode [19].
Among animal models for studying influenza virus, mice are usually unable to transmit infection from animal to animal [20,21,22], which makes the system unsuitable for studying virus transmissibility. Ferrets present a well-established model to study both influenza virus transmission and pathogenesis. They are susceptible to non-adapted swine and human influenza viruses and provide contact-dependent transmission from infected ferrets to an uninfected animal [23]. At the same time, highly pathogenic avian influenza viruses are non-transmissible between ferrets [24]. However, the ferret model has several practical disadvantages: they are expensive, relatively large, and thus require complex housing conditions. In addition, it is difficult to find ferrets that have not been previously infected with influenza. In this regard, guinea pig models of virus infection were developed, which are easy to handle and less expensive [25]. A more affordable and less demanding guinea pig model allows one to conduct more experiments and, thus, increases the statistical power of transmission experiments. Guinea pigs are easily infected by human and avian influenza viruses without adaptation. However, they either lack or exhibit extremely mild clinical symptoms of influenza infection compared to ferrets and humans. Although influenza viruses can replicate at high titers in the respiratory tract of guinea pigs, they usually do not exhibit severe disease symptoms even when infected with strains pathogenic to humans and ferrets [26].
The presence of pronounced clinical disease signs, the specificity of human influenza virus infection alongside the possibility of non-contact influenza virus transmission make ferrets the most suitable model to study virus transmissibility. Natural infection has an obvious advantage, since it provides a realistic infectious dose and results in virus replication kinetics that better imitate the kinetics of natural influenza virus infection in humans [27,28,29,30].
A person infected with the influenza virus contaminates the environment. This infectious material may be transmitted between people in many different ways, and the importance of different routes of infection is not yet known. The virus is shed by the infected host during events such as coughing, sneezing and talking. An “expiratory spray” is formed from particles of different sizes in which the virus is present. The virus is introduced to a new host by inhalation or contact. Contact transmission for the influenza virus is not considered dominant. Droplet transmission is often considered significant because the droplets have a high infectious potential and carry a high infectious dose; however, the droplets cannot reach the target cells in URT by inhalation and the role of the droplets is limited [31] Inhaled airborne infectious particles (≤5 µm) can be shed by patients during coughing. The pandemic and epidemic influenza viruses that have circulated in humans throughout the past century were transmitted mostly via the airborne route. Of all the pathways, perhaps the most interesting is the airborne transmission of influenza; evidence to support the importance of this overall contribution is increasing but is still inconclusive [28,32]. Droplet size distribution in human coughing has shown droplet size to be coughed up between 0.6–15.9 µm and an average modal size of 8.35 µm, although the cough droplet size distribution is multimodal and has three peaks at about 1 µm, 2 µm and 8 µm [33].
In infected ferrets, seasonal A/H1N1 and A/H3N2 viruses and mammal-adapted avian A/H5N1 virus are transmitted via the air in the URT, more specifically, from the nasal respiratory epithelium, and not from the trachea, bronchus or lungs [34]. Viruses are expelled from URT via breathing, sneezing, coughing and vocal activity. In a ferret model, influenza virus transmission between animals was mediated by airborne particles larger than 1.5 μm, consistent with the quantity and size of virus-laden particles released by the donors [35]. Onward transmission by donors was most efficient before fever onset and may continue for 5 days after inoculation. Interestingly, that limited amount of virus-laden particles at submicron size were exhaled by influenza-inoculated donor ferrets, nevertheless, the researchers demonstrated that droplets smaller than 1.5 μm, artificially generated by nebulizer, efficiently infect the recipient ferrets with influenza virus.
Despite numerous studies of the influenza virus on the ferret model [27,28,29,30,36,37,38,39] there remains a knowledge gap in understanding the mechanism of transmissibility of influenza viruses and difficulties in assessing the epidemiological danger of the influenza virus based on data obtained in the ferret model. For an in-depth study of the transmissibility of the influenza virus, we have developed an aerobiological chamber. In this study on an aerobiological chamber, we have made a comparative assessment of the airborne transmission of six influenza viruses of human and animal origin in a ferret model.

2. Materials and Methods

2.1. Viruses

The pandemic influenza virus A/California/07/2009 (H1N1)pdm09 was provided by the Centers for Disease Control and Prevention (CDC, Atlanta, Georgia, USA). Avian influenza viruses A/chicken/Primorsky Krai/1771/2018 (H9N2), A/chicken/Kostroma/1718/2017 (H5N2), A/chicken/Nghe An/08VTC/2018 (H5N1), and A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2), as well as swine influenza virus A/swine/Irkutsk/155/2017 (H3N2), were obtained from the collection of microorganisms of the State Research Center of Virology and Biotechnology “Vector”. Influenza viruses were propagated in 9-day-old chicken embryos.

2.2. Animals, Animal Procedures, and Intranasal ID50

Female ferrets aged 6–8 months were obtained from the breeding facilities of laboratory animals of the State Research Center of Virology and Biotechnology “Vector”. Animals were kept in individually ventilated cages on a standard diet. Animal housing and all procedures and manipulations with animals were performed in accordance with the request of the State Research Center of Virology and Biotechnology “Vector” on 3 June 2022, protocol No. 2 (28 June 2022) of the Bioethics Commission of the State Research Center of Virology and Biotechnology “Vector”. Serological testing was conducted to confirm the absence of virus-specific antibodies in animals before the start of the study and verify virus transmission from infected to healthy animals during the experiment. Donor ferrets were infected intranasally under anesthesia with a combination of Zoletil 100 (Delpharm Tours, France) and Xyla (Interchemie, Estonia); 0.25 mL of viral suspension was administered into each nostril. Animals were divided into four groups of 3 animals each to measure intranasal ID50. Animals of each group received a serial dilution of an influenza virus strain intranasally. The disease was diagnosed in infected animals based on the presence of viral RNA in the nasal washings 4 and 6 days after infection and seroconversion status 21 days after infection. Probit analysis was used to calculate intranasal ID50.

2.3. Aerobiology Chamber, Analytical Filters, and Particle Counter

To study transmissibility of influenza A viruses, an original ELC 04-60 aerobiology chamber (TIEGEL GmbH, Radeberg, Germany) was used (Figure 1). The chamber is designed for experiments in small and medium laboratory animals. The system consists of a technical air preparation module and four animal cabinet-chambers for animal exposure. The chamber is climate-controlled with the temperature range of 4–20 °C and relative humidity of 20–60%. A relative humidity of 30% and air temperature of 20 °C with decreased air pressure (−200 Pa) were maintained in the cabinet-chamber during experiments on influenza virus transmissibility in the ferret model. Aerosol samples from each cabinet-chamber were obtained through two sampling outlets to AFA-BA-3 analytical aerosol filters. Samples from the third air duct entered the aerosol particle counter (TSI AeroTrak APC 9303-01, Shoreview, Minnesota, USA), which operated in the following mode: 5 min/5 min during a 5 h animal exposure.
In order to assess the transmissibility of influenza A viruses in ferrets, animals were divided into two groups: donor group of 12 animals and recipient group of 12 animals. Donor ferrets were divided into four subgroups of 3 animals each a day before the experiment, and each subgroup was infected intranasally with an influenza virus strain at a dose of 1–106 FFU. Three animals of each subgroup of donor ferrets were placed in one of the four ELC 04-60 TIEGEL cabinet-chambers a day after infection. After this, a corresponding group of three naïve recipient ferrets was placed in a different cage of each cabinet-chamber. The distance between the cages was 10 cm, which excluded the direct contact between infected animals and recipients. The airflow was directed from infected donor ferrets to recipient ferrets; the flow speed was 3 cm/s. Simultaneous exposure of animals lasted 5 h daily for six days. During the exposure, air samples were captured onto particle counters and AFA-BA-3 analytical aerosol filters. Virus disease in animals was detected based on the presence of viral RNA in the nasal washings within 14 days after intranasal infection or airborne contact and seroconversion status 21 days after intranasal infection or airborne contact. The experiments were conducted one time.
The airborne infectious dose (AID) inhaled by the recipient ferret was calculated from viral load on two analytical aerosol filters of the chamber, expressed in FFU-equivalent, and minute respiratory volume in ferrets as:
A I D = v F A × M V × t ,
where v F —average viral load on analytical aerosol filter, FFU-equivalent; A —volume of air passing through the analytical aerosol filter during exposure, L.; M V —ferret minute volume, L/min; t —time of exposure, min. Calculation ferret minute volume M V was conducted in accordance with [40].

2.4. Determination of Viral Load in Biological Samples in MDCK Cells Using FFU Assay

Tenfold dilutions of nasal washingss from laboratory animals were prepared and added to MDCK cell monolayer in 96-well plates. Growth medium was removed from flat-bottom plate wells with cell monolayer, cells were washed twice with supportive growth medium (DMEM, 100 μ/mL penicillin, 25 mM HEPES, 0.2% solution of bovine serum albumin, 2 μg/mL of TPCK). Next, dilutions of test samples were added, and cells were incubated for 1 h at 37 °C and 5% CO2. The medium was then removed from wells, and cells were washed with supportive medium once. A total of 150 μL of supportive medium was added to all wells. The medium was removed after 18–20-h incubation at 37 °C and 5% CO2, and 200 μL of 80% cold acetone (−20 °C) was added, followed by incubation for 20–30 min. Acetone was removed, and the wells were dried. Next, 50 µL of diluted 1:1000 mouse monoclonal antibodies to influenza virus nucleoprotein was added to each well. The plate was incubated at 37 °C for 1 h, wells were washed three times with PBS, and secondary rabbit anti-mouse IgG antibodies conjugated to horseradish peroxidase (Abcam) were added at a ratio of 1:1000. After 1 h incubation, wells were washed three times with PBS, and AEC substrate (3-amino-9-ethylcarbazole, Sigma) solution was added. After 30 min incubation, the solution was removed, and the plate was washed with PBS once. Infected red-brown cells were counted, and the virus titer was calculated as the number of focus-forming units (FFU) per 1 mL of nasal washings (FFU/mL). The analysis was performed in duplicate wells.

2.5. Determination of Virus RNA Load in Biological Samples by RT-PCR

Influenza virus RNA was quantified in nasal washingss and analytical filters by quantitative RT-PCR with real-time data acquisition. The RIBO-sorb kit (InterLabService Ltd., Moscow, Russia) was used for RNA isolation. Reverse transcription was carried out using the Reverta-L kit (InterLabService Ltd., Moscow, Russia) in an incubator. The resulting influenza A cDNA fragments were amplified using the AmpliSens Influenza virus A/B-FL kit (InterLabService Ltd., Moscow, Russia). PCR and data registration were conducted on a RotorGene 6000 real-time cycler. The calibration curve of viral RNA Ct values on the concentration of influenza virus in FFU/mL was plotted to calculate the concentration of influenza virus strains in FFU-equivalents/mL. The PCR assay was performed without repeats and was considered correct when all control samples of a given test system were released.

2.6. Statistical Analysis

Statistical data processing was performed using Origin 8.1 data analysis package and Microsoft Excel. The analysis was carried out with log-transformed data due to the insignificance in this study of linear scale errors in exponential processes which, in our case, is virus reproduction.
Linear regression was used to test the relationships of measured viral particle concentrations in nasal washingss to the dose of inoculation due to the ability to test for the trend. Similarly, regression analysis was used to evaluate the calculated intranasal ID50 and the calculated dose–response function. Trends with R2 values above 0.5 were considered satisfactory.

3. Results

3.1. Virus Yield in the URT Cells of Donor Ferrets

Figure 2 shows virus concentration in the nasal washings of ferrets after intranasal infection with different doses of six influenza virus strains. Pandemic influenza virus A/California/07/2009 (H1N1)pdm09 is first detected in the nasal washings of ferrets upon infection with a minimal dose of 5 FFU. Virus concentration in the nasal washings reaches the maximum value of 2.3 × 106 FFU/mL at the infectious dose of 5 × 102 FFU, while further increase in infection dose does not result in the growth in virus concentration on the mucous membrane of the URT in ferrets (Figure 2A). Avian influenza viruses A(H5N2) and A(H5N1) at low doses of intranasal infection induce virus accumulation in the nasal washings at minimal concentrations (<4.5 × 103 FFU/mL), while virus concentration in the nasal washings does not exceed 2.9 × 104 FFU/mL for both strains at infectious dose of >1000 FFU (Figure 2B,D).
The opposite dependence of the virus concentration on the mucus membrane of the URT in ferrets on infectious dose is noted for avian influenza viruses A(H9N2) and swine influenza virus A(H3N2). Virus concentration in the nasal washings decreases with an increase in the dose of intranasal infection. Virus with a maximal concentration of 3.3 × 104 – 3.9 × 105 FFU/mL is detected in the nasal washings upon infection with avian influenza A(H9N2) viruses at infectious doses of <102 FFU. Virus concentration in the washings decreases to 1.4 × 104 FFU/mL upon an increase in intranasal infection dose to 106 FFU (Figure 2C,E). For swine influenza A(H3N2) virus, the maximum virus concentration in the nasal washings is 1.6 × 104 FFU/mL for infectious dose of 102 FFU; it further decreases to 1.1 × 103 FFU/mL with an increase in infectious dose to 105 FFU (Figure 2F).
The obtained results indicate that a decrease in virus concentration on the mucus membrane of the URT is observed with an increase in intranasal infectious dose (IID) for several strains. Such a decrease in virus yield can be explained by induction of non-specific antiviral defense mechanism, for instance, via the interferon pathway. Two avian influenza viruses A/chicken/Nghe An/08VTC/2018 (H5N1) and A/chicken/Kostroma/1718/2017 (H5N2), despite the pronounced specificity to avian alpha-2′3 receptors, demonstrated the ability to infect ferrets at high infectious doses. This is apparently due to the presence of low amounts of receptors to avian influenza on the cell surface of the mucous membrane of the URT in ferrets. On the contrary, absolute predominance of human receptors type alpha-2′6 on the mucous membrane of the URT in ferrets provides high virus yields on the mucus membrane upon infection with low doses of A/California/07/2009 (H1N1)pdm09, A/chicken/Primorsky Krai/1771/2018 (H9N2), and A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2) in ferrets. At the same time, the pandemic virus A/California/07/2009 (H1N1)pdm09 is accumulated in the URT at concentrations 10 times higher than those for A(H9N2) viruses and 100 times higher than those for A(H3N2), A(H5N2), and A(H5N1).

3.2. Genomic Characterization of Influenza Virus Strains Used in the Study

Complete genomes of 6 influenza viruses A/California/07/2009 (H1N1)pdm09, A/chicken/Primorsky Krai/1771/2018 (H9N2) (Y280 genetic line), A/chicken/Kostroma/1718/2017 (H5N2) (clade 2.3.4.4b) [41], A/chicken/Nghe An/08VTC/2018 (H5N1) (clade 2.3.2.1c), A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2) (genetic line Y280) and A/swine/Irkutsk/155/2017 (H3N2) were obtained using NGS for genetic analysis. Genetic analysis of the obtained genome sequences revealed the presence of some molecular markers of virulence and transmissibility in the HA and NS genes of the H5N1, H5N2, and H9N2 viruses. The 158D (H3 numbering) mutation in the HA gene results in the loss of the N-glycosylation site, which is often associated with a change in receptor specificity. The 145G mutation increases binding to alpha 2,6 receptors. The polybasic site in the HA gene greatly increases pathogenicity. Mutations in the NS1 205N and NS2 47A genes lead to a decrease in the antiviral response in mammals. Deletion 80-84 AA in the same NS1 gene increases virulence, and the PDZ binding domain of the NS1 gene ESEV suppresses the antiviral response (H5N1 genetic changes inventory) (Table 1).

3.3. Airborne Infection of Recipient Ferrets

AFA-BA-3 analytical aerosol filters were used to assess biological contamination of air in the aerobiological chamber. Air samples were collected onto analytical filters during the entire period of cohabitation exposure of donor and recipient ferrets. During this time, the virus trapped on the filter was inactivated, so the average virus concentration in the air of the aerobiology chamber was calculated from the viral RNA amount in the analytical filter, which was estimated by quantitative RT-PCR, and expressed in FFU-equivalents in accordance with the calibration curve. We also measured the concentration of aerosol particles in the chamber for three fractions using the particle counter: 0.3–1 µm, 1–5 µm, and >5 µm.
Ferrets create airborne pollution in the chamber. The maximum number of 16–63×103 aerosol particles of 0.3–1 µm per air liter was registered (Table 2). Influenza virus transmissibility in ferrets is known to be provided by aerosol particles of >1.5 µm [35]. The number of particles of the size 1.0–5.0 µm in the aerobiology chamber is an order of magnitude lower, their concentration is in the range of 1.7–6.4 × 103 particles per air liter (Table 2). Apparently, virus particles with a concentration in the air reaching the maximum value on day 4 after infection (0.012 FFU-equivalent/L) and then decreasing to 0.003 FFU/mL by day 6 after infection are distributed in this fraction. On day 4 after infection, there is an average of one viral particle (FFU-equivalent) per approximately 3.47 × 105 aerosol particles of >1 µm, the size corresponding to influenza virus transmissibility in ferrets. On day 6 after infection, the viral load in airborne particles decreases almost seven-fold, and one FFU-equivalent per 2.16 × 106 influenza transmissible airborne particles is observed.
During the period of cohabitation exposure in aerobiology chamber, recipient ferrets receive an airborne infectious dose (AID) of influenza virus from the air contaminated by donor ferrets. As shown in Figure 2, different virus concentrations were detected in the URT washings of donor ferrets infected with different IID that induced the different concentrations of airborne virus in the chamber.
Figure 3 presents the dependence of calculated dose inhaled by recipient ferrets on IID in donor ferrets.
Figure 3 shows that recipient ferrets received AID of 3.4 × 102–5.5 × 103 FFU-equivalents of A/swine/Irkutsk/155/2017 (H3N2) and A/California/07/2009 (H1N1)pdm09 influenza viruses, which is significantly higher than AID50 for this strains (25–52 FFU-equ). Recipient ferrets inhaled about 10 FFU-equivalents of A/chicken/Primorsky Krai/1771/2018 (H9N2) influenza virus, which is close to lower value of AID50 confidence interval at 95% for this strain (13 FFU).
Other influenza virus strains induced chamber contamination to a much lesser extent: recipient ferrets inhaled no more than 8 FFU-equivalents during cohabitation with donor ferrets.
Figure 4 presents the dependence of the probability of infection in recipient ferrets kept in the downward airflow to infectious donor ferrets in the aerobiology chamber on the influenza virus dose inhaled by recipient ferrets from the air contaminated by donor ferrets. The inhaled dose is presented as IID50, which makes it possible to combine the results for all studied strains in one graph.
The dose–response curve was reconstructed using Origin 8.1 software package by approximating the experimental data using the Levenberg–Marquardt algorithm and the function:
P A I D = 1 n = 1 t 1 1 2 1 + e r f l o g A I D n log A I D 50 2 σ 2 ,
where «t» is the time from the start of donors’ impact on recipients, with AID50 and σ as parameters. The values of inhaled AID were normalized by the value of IID50, which were measured for each of six strains in intranasally infected ferrets.
The obtained correlation is characterized by the following parameters: AID50 = 16.75 IID50 (3.14–89.33, p = 0.95) and σ = 1.18 ± 0.43. Thus, we managed to evaluate the dose providing infection in 50% recipient ferrets by airborne transmission from donor ferrets. This dose means that more than 50% of recipient ferrets are infected after inhaling >16.75 IID50.
Thus, transmissibility of influenza virus in the ferret model is highly probable in the case when donor ferrets can maintain a significant virus level in the surrounding air, which provides an inhaled dose of influenza virus sufficient enough to infect recipient ferrets.

4. Discussion

We proposed a ferret model in the aerobiology chamber study and quantified transmission of human and animal influenza viruses using this model to study influenza virus transmissibility. Ferrets are a well-established model for both transmission and pathogenesis of influenza viruses because their clinical symptoms are similar to those of humans [23,24,27,28,42,43,44,45,46,47].
Despite numerous studies on influenza virus transmissibility in the ferret model, there are no clearly established conditions determining the level of viral load in the environment that is required for donor ferrets to provide virus transmission to recipient ferrets. Animal experiments provide reasonable estimates for human susceptibility in dose–response studies [48]. A dose–response function connects a number of pathogens with a probability of infection. Various mathematical dose–response models are used to describe the probability of infection of the body [49]. So, from the hypothesis of independent action for a homogeneous population, a one-parameter exponential dose-effect relationship P(ID) = 1−exp(−pID), is valid [50]. Here, P(ID) is the probability of infecting an organism when a certain number of virions enters it, having a Poisson distribution with an average value equal to ID, the parameter p is the probability of infecting the host with one virion. The probability of infecting the host with a single virion, p, is related to the 50% infectious dose ID50 of the virus to the host by a simple ratio p = ln2/ID50, where ID50 is expressed as the number of virions. The exponential one-parameter dose–response model is useful for microbiology, but researchers widely use the popular two-parameter Poisson beta function model P(ID) = 1 − [1 + (ID/ID50)(21/α − 1)]−α with parameters ID50 and α [51]. The two-parameter Poisson beta function is the best known of the algebraic dose–response models with sigmoid curves. It fits well with another two-parameter lognormal distribution model P(ID) = 1/2[1 + erf(lg(ID/ID50)/21/2σlg)] with parameters ID50 and σlg, where ID50 is an infecting dose corresponding to a 50% response, σlg is a standard deviation for a normally distributed logarithm of ID, “erf” is an error function. The lognormal distribution for representing dose–response data is indistinguishable from the Poisson beta function for the most practical purposes. Mathematically, it can be shown that if the probability of infection (PN) is the product of the set of underlying probabilities of several events (Pis), i.e., P(N) = ∏ P(i), and if each of these underlying P(i) is uniformly distributed within its defined range, then the resulting probability of infection (PN) has a lognormal distribution [52,53]. The parameters of the lognormal distribution model are intuitive clear for researchers and can be measured experimentally. For this reason, we have chosen a lognormal distribution model to describe our experimental results in ferrets.
Our study demonstrates the relationship between influenza airborne and intranasal ID50 in the ferret model which, in turn, substantially determines airborne transmissibility.
Influenza virus transmission between ferrets is highly probable when donor ferrets can maintain the virus level in the environment that provides the inhaled dose of influenza virus of >16.75 IID50 in recipient ferrets for the studied influenza virus strains. Lower level of the virus in the environment, even a single virus, can also induce airborne transmission, but at very small probability.
Airborne particle contamination is provided by donor ferrets due to airborne excretion of influenza virus from the infected animals. The significant difference between IID50 and AID50 value may be explained by the different route of the virus entering into the nasal cavity. The IID50 value was measured for intranasally infected ferrets, whereas AID50 is obtained in ferrets infected via airborne route. In this case, only part of the inhaled AID may be sedimented in nasal mucosa, another part of the inhaled virus is eliminated during exhalation, and the sedimented virus dose in the nose is smaller than the inhaled virus dose [54].
The pattern of influenza virus accumulation on the mucous membrane of the URT in ferrets is unique for each specific strain. A comparative study of intranasal infection in ferrets with six influenza A virus strains showed that three strains demonstrate a decrease in viral concentration on the mucous membrane of the URT with an increase in the dose of intranasal infection. The latter three strains increase the virus yield on the mucosa of the URT with an increasing infectious dose. These three viruses include 2009 pandemic influenza strain A/California/07/2009 (H1N1)pdm09 and avian influenza A(H5) strains A/chicken/Kostroma/1718/2017 (H5N2) and A/chicken/Nghe An/08VTC/2018 (H5N1). The remaining three strains, which caused avian influenza virus outbreaks in poultry farms, namely A/chicken/Primorsky Krai/1771/2018 (H9N2) and A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2), as well as swine influenza A/swine/Irkutsk/155/2017 (H3N2) virus, circulating in a pig farm in the Irkutsk region, demonstrated either weak or strong negative correlation between virus concentration on the mucous membrane of the URT in ferrets and an intranasal infectious dose with an increase in intranasal infectious dose. Non-specific antiviral response was hypothesized to explain a decrease in viral replication upon an increase in infectious dose for the three studied strains of subtypes (H9N2) and A(H3N2). The opposite pattern for other strains, namely an increase in viral yield with an increase in infectious dose, primarily for A(H5) strains, is apparently due to inhibition of the host antiviral response to infection.
This study revealed significant differences in the reproduction of influenza viruses on the mucosa of the URT of ferrets for different strains. Genetic analysis revealed amino acid substitutions, which were correlated with the characteristics of a nonspecific antiviral reaction of the organism. The two studied strains A/chicken/Kostroma/1718/2017 (H5N2) and A/chicken/Nghe An/08VTC/2018 (H5N1) with a positive correlation between intranasal infectious dose and virus yield on the mucous membrane of the URT differ in the presence of the C-terminal amino acid region ESEV in the PDZ-binding motif of NS1, which prevents early apoptosis and inhibits expression of interferon-stimulated genes thus providing a high level of virus replication [55,56]. Furthermore, hemagglutinin in these two viruses differs in the presence of the polybasic site R-X-R/K-R-R-K-R-G at 323–330. This region is susceptible to proteolytic cleavage, found in many A(H5N1) viruses, and associated with a pronounced increase in pathogenicity [57,58].
Two influenza A(H9N2) virus strains do not possess the above genetic markers associated with a high level of viral replication. However, proteins NS1 and NS2 contain two related substitutions (N200S in NS1 and T47A in NS2), which provide decreased antiviral response in the host [36]. This is manifested in a weak correlation between virus accumulation in the URT in ferrets and low negative dependence on the infectious dose.
Swine influenza virus A/swine/Irkutsk/155/2017 (H3N2) is the only strain lacking genetic markers of host’s antiviral response inhibition and showing the maximum negative correlation between viral yield and infectious dose.
In this study, the ability to transmit H1N1pdm09, H3N2 viruses isolated from pigs, H5N1, H5N2 and 2 H9N2 viruses isolated from birds was studied in ferrets. Airborne transmissibility in ferrets has been shown for seasonal influenza viruses H1N1pdm09 A/California/07/2009 and H3N2 A/swine/Irkutsk/155/2017. Ferrets are animal models for human influenza infection and have been shown to be transmissible with seasonal influenza viruses [37]. The H3N2 A/swine/Irkutsk/155/2017 virus is reassortant with the HA, NA genes from the H3N2 seasonal influenza virus and internal genes from the seasonal influenza virus of another H1N1pdm09 subtype. Reassortant H9N2 viruses with internal genes from seasonal influenza viruses have previously been shown to be able to replicate and be transmitted in ferrets [59]. Since the introduction of A/goose/Guangdong/1/96 (H5N1) in China in 1996, influenza viruses of this subtype have spread throughout the world. As a result of reassortment with low pathogenic influenza viruses, there has been a further increase in genetic diversity and the emergence of viruses of different subtypes, which can be designated as H5Nx [60]. It has been established that most studied H5Nx viruses of clade 2.3.2.1. and 2.3.4.4. are not transmitted between mammals [61]. However, some of them have a limited ability to be transmissible by contact in ferrets and other animal models [62]. As a result of the work, it was found that the studied viruses A/chicken/Nghe An/08VTC/2018 (H5N1) and A/chicken/Kostroma/1718/2017 (H5N2) are not transmitted between ferrets by airborne particles. Limited available information about genetic markers of transmissibility often does not allow predicting or evaluating transmissibility and, in some cases, human and avian H5Nx viruses with the same known markers may differ in transmissibility [61]. The main markers of transmissibility are markers that determine the receptor specificity of the virus, the stability of the HA protein when the virus exits the endosome into the cell, as well as markers associated with virus replication. Basically, these markers are located in different positions of the receptor-binding site of the HA protein and the PB2 protein [28,63]. The influence of markers in HA at positions 158 (position associated with the glycosylation site), 224, 226, 318 (H3 numbering) on transmissibility was established in a ferret model study of an artificial H5N2 reassortant containing HA from H5N2 and NA and internal genes from the H3N2 seasonal influenza virus [63]. Obtained reassortant virus with HA from avian influenza and internal genes from seasonal influenza virus became transmissible via airborne particles. Another study showed that the G228S mutation in HA of H5Nx viruses, in addition to those listed, affects the transmissibility of the virus in ferrets by direct contact [64]. The profile of the studied viruses A/chicken/Nghe An/08VTC/2018 (H5N1) and A/chicken/Kostroma/1718/2017 (H5N2) differed between them by one amino acid at position 158 of the HA gene and was as follows: HA 110H, 158D/N, 160A, 224N, 226Q, 228G, 318T, PB2 627E, which corresponded to the “avian” type. At the same time, the remaining genes of the studied H5Nx viruses may contain markers that have not yet been identified that affect transmissibility.
The published literature also shows limited direct contact transmission in ferrets for H9N2 subtype viruses [65]. For the A/chicken/Primorsky Krai/1771/2018 virus studied in this work, airborne transmission was shown, which was established by the presence of specific antibodies in the blood on day 21 after exposure. For the A/chicken/Thanh Hoa/V1S5VTC/2020 virus, no ferret-to-ferret transmission has been shown. The lack of isolation of the virus by the recipient ferret may indicate a low adaptability of this virus to mammals. The main markers of transmissibility of H9N2 viruses are markers that determine receptor specificity. Significance of transmissibility markers is strongly dependent on their position in the genome on their genetic environment. Thus, the decisive role of the Q226L [H3 numbering] mutation was previously demonstrated [65], but modern viruses of the H9N2 subtype almost all contain the Q226L mutation. Then, many mutations were found on another viruses, including the key ones at positions 190, 226, and 227, which determine receptor specificity and avidity of binding to receptors [66]. Further monitoring of H9N2 viruses made it possible to establish that mutations T205A, D208E, V216L, V245I in viruses circulating in 2002-2005 and mutations S119R, D145G, Q156R, A160D, T212I, Q227M, R246K in viruses, circulating in 2011-2017, increased affinity for either human-type receptors or both human- and avian-type receptors [67]. According to GISAID EpiFlu, almost all H9N2 viruses isolated from human infections in 2018-2020 have these mutations. At the same time, some of these mutations were absent in avian H9N2 virus genomes. Genetic analysis of the studied influenza viruses A/chicken/Primorsky Krai/1771/2018 and A/chicken/Thanh Hoa/V1S5VTC/2020 showed that almost all of the mutations associated with increased affinity to human type receptors were present in the HA gene, except for markers in positions 119, 145 and 160, which had different amino acids (119Q, 145N, 160N). It is possible that absence of these markers of host specificity and potential absence of other markers not yet identified make the viruses less adapted for transmission in mammals, in particular, ferrets, compared to bird-to-human transmission. The studied viruses differ from each other by the substitution at position 145N/G in HA, and this difference may be the reason for the observed difference in transmissibility.
Virus replication on the mucous membrane of the URT in ferrets leads to viral contamination of the chamber air. We showed that donor ferrets infected intranasally with the pandemic influenza A/California/07/2009 (H1N1)pdm09 virus create airborne contamination in the aerobiology chamber at the level providing infection in recipient ferrets with the respiratory dose of 4.2 × 103–5.5 × 103 FFU-equivalents, which is 100 times higher than AID50 for this strain. Similarly, donor ferrets intranasally infected with swine influenza virus A/swine/Irkutsk/155/2017 (H3N2) create viral contamination inside the aerobiological chamber at such a level that recipient ferrets receive a respiratory dose of 1.3 × 102–3.4 × 102 FFU-equivalents, which is more than 10 times the value of the AID50 for this strain. Part of the donor ferrets intranasally infected with avian influenza virus A/chicken/Primorsky Krai/1771/2018 (H9N2) create viral contamination inside the aerobiological chamber at such a level that recipient ferrets receive a respiratory dose about 10 FFU-equivalents, which is close to lower value of AID50 95% confidence interval.
The other studied influenza viruses induced chamber contamination to a significantly lesser extent, and recipient ferrets inhaled <8 FFU-equivalents during cohabitation with donor ferrets.
Despite the fact that virus was detected in the nasal washings of ferrets after intranasal infection with each of the six strains, only three strains demonstrated the ability of donor ferrets to contaminate the chamber environment at the level, when the airborne virus concentration was sufficient to infect recipient ferrets.

5. Conclusions

Our results obtained in the ferret model cannot be directly extrapolated to humans. However, they allow for ranking viruses based on their transmissibility between mammals. The obtained results make it possible to develop an evidence base of phenotypic and genetic markers of transmissibility of new animal influenza viruses to carry out timely anti-epidemic measures by constantly monitoring viruses in the environment in order to preserve the life and health of the human population.

Author Contributions

Conceptualization, A.R. and E.G.; data curation, A.G., M.K. and A.D.; formal analysis, A.G., V.M., I.S. and S.S.; investigation, A.G., G.O. and T.T.N.; methodology, A.G., N.D. and A.M.; project administration, A.G., A.K., E.G., R.M. and A.R.; supervision, I.S., E.G., R.M. and A.R.; visualization, A.G. and M.K.; writing—original draft preparation, A.R., A.G., M.K., A.D. and N.K.; writing—review and editing, A.R., A.G., M.K., A.D. and N.K. All authors have read and agreed to the published version of the manuscript.

Funding

State Assignment No. 13/21 (FBRI SRC VB ‘Vector’ Rospotrebnadzor) supported this research.

Institutional Review Board Statement

The animal study protocol was approved by the Ethics Committee of Federal Budgetary Research Institution State Research Center of Virology and Biotechnology “Vector” Rospotrebnadzor (No. 2 of 28.06.22).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AEC3-amino-9-ethylcarbazole
AIDAirborne Infectious Dose
AID5050% Airborne Infectious Dose
cDNAComplementary Deoxyribonucleic Acid
DMEMDulbecco’s Modified Eagle Medium
ELCEuroLabClime
FFUFocus-Forming Units
GISAIDGlobal Initiative on Sharing All Influenza Data
HEPES4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid
ID5050% Infectious Dose
IgGImmunoglobulin G
IIDIntranasal Infectious Dose
IID5050% Intranasal Infectious Dose
MDCKMadin-Darby Canine Kidney
NGSNext-Generation Sequencing
NSNonstructural Protein
PBSPhosphate-Buffered Saline
RNARiboNucleic Acid
RT-PCRReverse Transcription, Polymerase Chain Reaction
TIPRATool for Influenza Pandemic Risk Assessment
TPCKTosyl Phenylalanyl Chloromethyl Ketone
URTThe Upper Respiratory Tract
WHOWorld Health Organization

References

  1. Krammer, F.; Smith, G.J.D.; Fouchier, R.A.M.; Peiris, M.; Kedzierska, K.; Doherty, P.C.; Palese, P.; Shaw, M.L.; Treanor, J.; Webster, R.G.; et al. Influenza. Nat. Rev. Dis. Primers 2018, 4, 1–21. [Google Scholar] [CrossRef] [PubMed]
  2. Wille, M.; Holmes, E.C. The Ecology and Evolution of Influenza Viruses. Cold Spring Harb. Perspect. Med. 2020, 10, a038489. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Anderson, R.M.; May, R.M. Infectious Diseases of Humans: Dynamics and Control; Oxford University Press: Oxford, UK, 1991. [Google Scholar]
  4. Zhang, S.-H. Estimating Transmissibility of Seasonal Influenza Virus by Surveillance Data. J. Data Sci. 2011, 9, 55–64. [Google Scholar] [CrossRef]
  5. Riley, S.; Wu, J.T.; Leung, G.M. Optimizing the Dose of Pre-Pandemic Influenza Vaccines to Reduce the Infection Attack Rate. PLoS Med. 2007, 4, 1032–1040. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Payet, C.; Voirin, N.; Vanhems, P.; Ecochard, R. A statistical model to assess the risk of communicable diseases associated with multiple exposures in healthcare settings. BMC Med. Res. Methodol. 2013, 13, 26. [Google Scholar] [CrossRef]
  7. Haas, C.N.; Rose, J.B.; Gerba, C.P. Quantitative Microbial Risk Assessment, 2nd ed.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2014. [Google Scholar]
  8. Haas, C.N.; Rose, J.B.; Gerba, C.P.; Regli, S. Risk Assessment of Virus in Drinking Water. Risk Anal. 1993, 13, 545–552. [Google Scholar] [CrossRef]
  9. Codeço, C.T. Endemic and epidemic dynamics of cholera: The role of the aquatic Reservoir. BMC Infect. Dis. 2001, 1, 1. [Google Scholar] [CrossRef] [Green Version]
  10. WHO. Tool for Influenza Pandemic Risk Assessment (TIPRA); WHO: Geneva, Switzerland, 2016; pp. 1–60. [Google Scholar]
  11. Webster, R.G.; Monto, A.S.; Braciale, T.J.; Lamb, R.A. Textbook of Influenza, 2nd ed.; Wiley-Blackwell: West Sussex, UK, 2015; pp. 175–189. [Google Scholar]
  12. WHO. Pandemic Influenza Risk Management: A WHO Guide to Inform and Harmonize National and International Pandemic Preparedness and Response; WHO: Geneva, Switzerland, 2017; pp. 1–62. [Google Scholar]
  13. Skehel, J.J.; Wiley, D.C. Receptor binding and membrane fusion in virus entry: The influenza hemagglutinin. Annu. Rev. Biochem. 2000, 69, 531–569. [Google Scholar] [CrossRef]
  14. Palese, P.; Compans, R.W. Inhibition of influenza virus replication in tissue culture by 2-deoxy-2,3-dehydro-N-trifluoroacetylneuraminic acid (FANA): Mechanism of action. J. Gen. Virol. 1976, 33, 159–163. [Google Scholar] [CrossRef]
  15. Baigent, S.J.; Bethell, R.C.; McCauley, J.W. Genetic analysis reveals that both haemagglutinin and neuraminidase determine the sensitivity of naturally occurring avian influenza viruses to zanamivir in vitro. Virology 1999, 263, 323–338. [Google Scholar] [CrossRef] [Green Version]
  16. Cohen, M.; Zhang, X.-Q.; Senaati, H.P.; Chen, H.-W.; Varki, N.M.; Schooley, R.T.; Gagneux, P. Influenza A penetrates host mucus by cleaving sialic acids with neuraminidase. Virol. J. 2013, 10, 321. [Google Scholar] [CrossRef] [Green Version]
  17. Matrosovich, M.N.; Matrosovich, T.Y.; Gray, T.; Roberts, N.A.; Klenk, H.-D. Neuraminidase is important for the initiation of influenza virus infection in human airway epithelium. J. Virol. 2004, 78, 12665–12667. [Google Scholar] [CrossRef] [Green Version]
  18. Shinya, K.; Ebina, M.; Yamada, S.; Ono, M.; Kasai, N.; Kawaoka, Y. Avian flu: Influenza virus receptors in the human airway. Nature 2006, 440, 435–436. [Google Scholar] [CrossRef]
  19. Benton, D.J.; Martin, S.R.; Wharton, S.A.; McCauley, J.W. Biophysical Measurement of the Balance of Influenza A Hemagglutinin and Neuraminidase Activities. J. Biol. Chem. 2015, 290, 6516–6521. [Google Scholar] [CrossRef] [Green Version]
  20. Schulman, J.L.; Kilbourne, E.D. Experimental transmission of influenza virus infection in mice. I. The period of transmissibility. J. Exp. Med. 1963, 118, 257–266. [Google Scholar] [CrossRef]
  21. Matsuoka, Y.; Lamirande, E.W.; Subbarao, K. The Mouse Model for Influenza. Curr. Protoc. Microbiol. 2009, 15, 15G.3.1–15G.3.30. [Google Scholar] [CrossRef]
  22. Edenborough, K.M.; Gilbertson, B.P.; Brown, L.E. A Mouse Model for the Study of Contact-Dependent Transmission of Influenza A Virus and the Factors That Govern Transmissibility. J. Virol. 2012, 86, 12544–12551. [Google Scholar] [CrossRef] [Green Version]
  23. Linster, M.; van Boheemen, S.; de Graaf, M.; Schrauwen, E.J.A.; Lexmond, P.; Mänz, B.; Bestebroer, T.M.; Baumann, J.; van Riel, D.; Rimmelzwaan, G.F.; et al. Identification, characterization, and natural selection of mutations driving airborne transmission of A/H5N1 virus. Cell 2014, 157, 329–339. [Google Scholar] [CrossRef] [Green Version]
  24. Yen, H.-L.; Lipatov, A.S.; Ilyushina, N.A.; Govorkova, E.A.; Franks, J.; Yilmaz, N.; Douglas, A.; Hay, A.; Krauss, S.; Rehg, J.E.; et al. Inefficient Transmission of H5N1 Influenza Viruses in a Ferret Contact Model. J. Virol. 2007, 81, 6890–6898. [Google Scholar] [CrossRef] [Green Version]
  25. Lowen, A.C.; Mubareka, S.; Tumpey, T.M.; García-Sastre, A.; Palese, P. The guinea pig as a transmission model for human influenza viruses. Proc. Natl. Acad. Sci. USA 2006, 103, 9988–9992. [Google Scholar] [CrossRef] [Green Version]
  26. Kwon, Y.K.; Lipatov, A.S.; Swayne, D.E. Bronchointerstitial pneumonia in guinea pigs following inoculation with H5N1 high pathogenicity avian influenza virus. Vet. Pathol. 2009, 46, 138–141. [Google Scholar] [CrossRef] [PubMed]
  27. Matsuoka, Y.; Lamirande, E.W.; Subbarao, K. The Ferret Model for Influenza. Curr. Protoc. Microbiol. 2009, 15, 15G.2.1–15G.2.29. [Google Scholar] [CrossRef] [PubMed]
  28. Herfst, S.; Schrauwen, E.J.; Linster, M.; Chutinimitkul, S.; de Wit, E.; Munster, V.J.; Sorrell, E.M.; Bestebroer, T.M.; Burke, D.F.; Smith, D.J.; et al. Airborne transmission of influenza A/H5N1 virus between ferrets. Science 2012, 336, 1534–1541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Van der Vries, E.; Veldhuis Kroeze, E.J.; Stittelaar, K.J.; Linster, M.; Van der Linden, A.; Schrauwen, E.J.; Leijten, L.M.; van Amerongen, G.; Schutten, M.; Kuiken, T.; et al. Multidrug resistant 2009 A/H1N1 influenza clinical isolate with a neuraminidase I223R mutation retains its virulence and transmissibility in ferrets. PLoS Pathog. 2011, 7, e1002276. [Google Scholar] [CrossRef] [Green Version]
  30. Hamelin, M.-E.; Baz, M.; Bouhy, X.; Beaulieu, E.; Dube, K.; Mallett, C.; Boivin, G. Reduced airborne transmission of oseltamivir-resistant pandemic A/H1N1 virus in ferrets. Antivir. Ther. 2011, 16, 775–779. [Google Scholar] [CrossRef] [Green Version]
  31. Wells, W.F. Airborne Contagion and Air Hygiene. An Ecological Study of Droplet Infections; Harvard University Press: Cambridge, MA, USA, 1955. [Google Scholar]
  32. Killingley, B.; Nguyen-Van-Tam, J. Routes of influenza transmission. Influenza Other Respir. Viruses 2013, 7, 42–51. [Google Scholar] [CrossRef] [Green Version]
  33. Yang, S.; Lee, G.W.; Chen, C.M.; Wu, C.C.; Yu, K.P. The size and concentration of droplets generated by coughing in human subjects. J. Aerosol Med. 2007, 20, 484–494. [Google Scholar] [CrossRef]
  34. Richard, M.; van den Brand, J.M.A.; Bestebroer, T.M.; Lexmond, P.; de Meulder, D.; Fouchier, R.A.M.; Lowen, A.C.; Herfst, S. Influenza A viruses are transmitted via the air from the nasal respiratory epithelium of ferrets. Nat. Commun. 2020, 11, 766. [Google Scholar] [CrossRef] [Green Version]
  35. Zhou, J.; Wei, J.; Choy, K.T.; Sia, S.F.; Rowlands, D.K.; Yu, D.; Wu, C.Y.; Lindsley, W.G.; Cowling, B.J.; McDevitt, J.; et al. Defining the sizes of airborne particles that mediate influenza transmission in ferrets. Proc. Natl. Acad. Sci. USA 2018, 115, 2386–2392. [Google Scholar] [CrossRef] [Green Version]
  36. Imai, H.; Shinya, K.; Takano, R.; Kiso, M.; Muramoto, Y.; Sakabe, S.; Murakami, S.; Ito, M.; Yamada, S.; Le, M.T.; et al. The HA and NS genes of human H5N1 influenza A virus contribute to high virulence in ferrets. PLoS Pathog. 2010, 6, e1001106. [Google Scholar] [CrossRef] [Green Version]
  37. Belser, J.A.; Pulit-Penaloza, J.A.; Maines, T.R. Ferreting Out Influenza Virus Pathogenicity and Transmissibility: Past and Future Risk Assessments in the Ferret Model. Cold Spring Harb. Perspect. Med. 2020, 10, a038323. [Google Scholar] [CrossRef] [Green Version]
  38. Belser, J.A.; Lau, E.H.Y.; Barclay, W.; Barr, I.G.; Chen, H.; Fouchier, R.A.M.; Hatta, M.; Herfst, S.; Kawaoka, Y.; Lakdawala, S.S.; et al. Robustness of the Ferret Model for Influenza Risk Assessment Studies: A Cross-Laboratory Exercise. mBio 2022, 13, e01174-22. [Google Scholar] [CrossRef]
  39. Belser, J.A.; Eckert, A.M.; Tumpey, T.M.; Maines, T.R. Complexities in Ferret Influenza Virus Pathogenesis and Transmission Models. Microbiol. Mol. Biol. Rev. 2016, 80, 733–744. [Google Scholar] [CrossRef] [Green Version]
  40. MacInnes, H.; Zhou, Y.; Gouveia, K.; Cromwell, J.; Lowery, K.; Layton, R.C.; Zubelewicz, M.; Sampath, R.; Hofstadler, S.; Liu, Y.; et al. Transmission of aerosolized seasonal H1N1 influenza A to ferrets. PLoS ONE 2011, 6, e24448. [Google Scholar] [CrossRef] [Green Version]
  41. Marchenko, V.; Goncharova, N.; Susloparov, I.; Kolosova, N.; Gudymo, A.; Svyatchenko, S.; Danilenko, A.; Durymanov, A.; Gavrilova, E.; Maksyutov, R.; et al. Isolation and characterization of H5Nx highly pathogenic avian influenza viruses of clade 2.3.4.4 in Russia. Virology 2018, 525, 216–223. [Google Scholar] [CrossRef]
  42. Fabian, P.; McDevitt, J.J.; DeHaan, W.H.; Fung, R.O.; Cowling, B.J.; Chan, K.H.; Leung, G.M.; Milton, D.K. Influenza virus in human exhaled breath: An observational study. PLoS ONE 2008, 3, e2691. [Google Scholar] [CrossRef] [Green Version]
  43. Koster, F.; Gouveia, K.; Zhou, Y.; Lowery, K.; Russell, R.; MacInnes, H.; Pollock, Z.; Layton, R.C.; Cromwell, J.; Toleno, D.; et al. Exhaled aerosol transmission of pandemic and seasonal H1N1 influenza viruses in the ferret. PLoS ONE 2012, 7, e33118. [Google Scholar] [CrossRef] [Green Version]
  44. Gustin, K.M.; Belser, J.A.; Wadford, D.A.; Pearce, M.B.; Katz, J.M.; Tumpey, T.M.; Maines, T.R. Influenza virus aerosol exposure and analytical system for ferrets. Proc. Natl. Acad. Sci. USA 2011, 108, 8432–8437. [Google Scholar] [CrossRef] [Green Version]
  45. Everett, H.E.; Nash, B.; Londt, B.Z.; Kelly, M.D.; Coward, V.; Nunez, A.; van Diemen, P.M.; Brown, I.H.; Brookes, S.M. Interspecies Transmission of Reassortant Swine Influenza A Virus Containing Genes from Swine Influenza A(H1N1)pdm09 and A(H1N2) Viruses. Emerg. Infect. Dis. 2020, 26, 273–281. [Google Scholar] [CrossRef]
  46. Zitzow, L.A.; Rowe, T.; Morken, T.; Shieh, W.J.; Zaki, S.; Katz, J.M. Pathogenesis of avian influenza A (H5N1) viruses in ferrets. J. Virol. 2002, 76, 4420–4429. [Google Scholar] [CrossRef] [Green Version]
  47. Belser, J.A.; Szretter, K.J.; Katz, J.M.; Tumpey, T.M. Use of animal models to understand the pandemic potential of highly pathogenic avian influenza viruses. Adv. Virus Res. 2009, 73, 55–97. [Google Scholar] [PubMed]
  48. Armstrong, T.W.; Haas, C.N. A quantitative microbial risk assessment model for Legionnaires’ disease: Animal model selection and dose-response modeling. Risk Anal. 2007, 27, 1581–1596. [Google Scholar] [CrossRef] [PubMed]
  49. Brouwer, A.F.; Weir, M.H.; Eisenberg, M.C.; Meza, R.; Eisenberg, J.N.S. Dose-response relationships for environmentally mediated infectious disease transmission models. PLoS Comput. Biol. 2017, 13, e1005481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Peto, S. A Dose-Response Equation for the Invasion of Micro-Organisms. Biometrics 1953, 9, 320–335. [Google Scholar] [CrossRef]
  51. Kitajima, M.; Huang, Y.; Watanabe, T.; Katayama, H.; Haas, C.N. Dose-response time modeling for highly pathogenic avian influenza A (H5N1) virus infection. Lett. Appl. Microbiol. 2011, 53, 438–444. [Google Scholar]
  52. Peleg, M. Microbial Dose-Response Curves and Disinfection Efficacy Models Revisited. Food Eng. Rev. 2021, 13, 305–321. [Google Scholar] [CrossRef]
  53. Peleg, M.; Normand, M.D.; Corradini, M.G. Construction of Food and Water Borne Pathogens’ Dose–Response Curves Using the Expanded Fermi Solution. J. Food Sci. 2011, 76, 82–89. [Google Scholar] [CrossRef]
  54. Rowe, B.R.; Canosa, A.; Meslem, A.; Rowe, F. Increased airborne transmission of COVID-19 with new variants, Implications for health policies. Build. Environ. 2022, 219, 109132. [Google Scholar] [CrossRef]
  55. Liu, H.; Golebiewski, L.; Dow, E.C.; Krug, R.M.; Javier, R.T.; Rice, A.P. The ESEV PDZ-binding motif of the avian influenza A virus NS1 protein protects infected cells from apoptosis by directly targeting scribble. J. Virol. 2010, 84, 11164–11174. [Google Scholar] [CrossRef] [Green Version]
  56. Shepardson, K.; Larson, K.; Cho, H.; Johns, L.L.; Malkoc, Z.; Stanek, K.; Wellhman, J.; Zaiser, S.; Daggs-Olson, J.; Moodie, T.; et al. A Novel Role for PDZ-Binding Motif of Influenza A Virus Nonstructural Protein 1 in Regulation of Host Susceptibility to Postinfluenza Bacterial Superinfections. Viral. Immunol. 2019, 32, 131–143. [Google Scholar] [CrossRef]
  57. Webster, R.G.; Rott, R. Influenza A virus pathogenicity: The pivotal role of hemagglutinin. Cell 1987, 50, 665–666. [Google Scholar] [CrossRef]
  58. Horimoto, T.; Kawaoka, Y. The hemagglutinin cleavability of a virulent avian influenza virus by subtilisin-like endoproteases is influenced by the amino acid immediately downstream of the cleavage site. Virology 1995, 210, 466–470. [Google Scholar] [CrossRef] [Green Version]
  59. Cáceres, C.J.; Rajao, D.S.; Perez, D.R. Airborne Transmission of Avian Origin H9N2 Influenza A Viruses in Mammals. Viruses 2021, 13, 1919. [Google Scholar] [CrossRef]
  60. Claes, F.; Morzaria, S.P.; Donis, R.O. Emergence and dissemination of clade 2.3.4.4 H5Nx influenza viruses-how is the Asian HPAI H5 lineage maintained. Curr. Opin. Virol. 2016, 16, 158–163. [Google Scholar] [CrossRef]
  61. Yamaji, R.; Saad, M.D.; Davis, C.T.; Swayne, D.E.; Wang, D.; Wong, F.Y.K.; McCauley, J.W.; Peiris, J.S.M.; Webby, R.J.; Fouchier, R.A.M.; et al. Pandemic potential of highly pathogenic avian influenza clade 2.3.4.4 A(H5) viruses. Rev. Med. Virol. 2020, 30, e2099. [Google Scholar] [CrossRef] [Green Version]
  62. Kaplan, B.S.; Russier, M.; Jeevan, T.; Marathe, B.; Govorkova, E.A.; Russell, C.J.; Kim-Torchetti, M.; Choi, Y.K.; Brown, I.; Saito, T.; et al. Novel Highly Pathogenic Avian A(H5N2) and A(H5N8) Influenza Viruses of Clade 2.3.4.4 from North America Have Limited Capacity for Replication and Transmission in Mammals. mSphere 2016, 1, e00003-16. [Google Scholar] [CrossRef] [Green Version]
  63. Imai, M.; Watanabe, T.; Hatta, M.; Das, S.C.; Ozawa, M.; Shinya, K.; Zhong, G.; Hanson, A.; Katsura, H.; Watanabe, S.; et al. Experimental adaptation of an influenza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature 2012, 486, 420–428. [Google Scholar] [CrossRef] [Green Version]
  64. Pulit-Penaloza, J.A.; Brock, N.; Pappas, C.; Sun, X.; Belser, J.A.; Zeng, H.; Tumpey, T.M.; Maines, T.R. Characterization of highly pathogenic avian influenza H5Nx viruses in the ferret model. Sci. Rep. 2020, 10, 12700. [Google Scholar] [CrossRef]
  65. Wan, H.; Sorrell, E.M.; Song, H.; Hossain, M.J.; Ramirez-Nieto, G.; Monne, I.; Stevens, J.; Cattoli, G.; Capua, I.; Chen, L.M.; et al. Replication and transmission of H9N2 influenza viruses in ferrets: Evaluation of pandemic potential. PLoS ONE 2008, 3, e2923. [Google Scholar] [CrossRef] [Green Version]
  66. Peacock, T.P.; Sealy, J.E.; Harvey, W.T.; Benton, D.J.; Reeve, R.; Iqbal, M. Genetic determinants of receptor-binding preference and zoonotic potential of H9N2 avian influenza viruses. J. Virol. 2020, 95, e01651-20. [Google Scholar] [CrossRef]
  67. Liu, Y.; Li, S.; Sun, H.; Pan, L.; Cui, X.; Zhu, X.; Feng, Y.; Li, M.; Yu, Y.; Wu, M.; et al. Variation and Molecular Basis for Enhancement of Receptor Binding of H9N2 Avian Influenza Viruses in China Isolates. Front. Microbiol. 2020, 11, 602124. [Google Scholar] [CrossRef] [PubMed]
Figure 1. TIEGEL ELC 04-60 dynamic aerobiology chamber.
Figure 1. TIEGEL ELC 04-60 dynamic aerobiology chamber.
Atmosphere 14 00471 g001
Figure 2. Effect of intranasal infectious dose (IID) on virus load in the nasal washings of donor ferrets, infected with different doses of influenza virus strains: (A) A/California/07/2009 (H1N1)pdm09, (B) A/chicken/Kostroma/1718/2017 (H5N2), (C) A/chicken/Primorsky Krai/1771/2018 (H9N2), (D) A/chicken/Nghe An/08VTC/2018 (H5N1), (E) A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2), (F) A/swine/Irkutsk/155/2017 (H3N2). The X axis—IID in donor ferrets. The Y-axis—virus titer in nasal washings of donor ferrets.
Figure 2. Effect of intranasal infectious dose (IID) on virus load in the nasal washings of donor ferrets, infected with different doses of influenza virus strains: (A) A/California/07/2009 (H1N1)pdm09, (B) A/chicken/Kostroma/1718/2017 (H5N2), (C) A/chicken/Primorsky Krai/1771/2018 (H9N2), (D) A/chicken/Nghe An/08VTC/2018 (H5N1), (E) A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2), (F) A/swine/Irkutsk/155/2017 (H3N2). The X axis—IID in donor ferrets. The Y-axis—virus titer in nasal washings of donor ferrets.
Atmosphere 14 00471 g002
Figure 3. Effect of IID in donor ferrets on estimated AID inhaled by recipients, for donor ferrets infected with influenza virus strains: (A) A/California/07/2009 (H1N1)pdm09, (B) A/chicken/Kostroma/1718/2017 (H5N2), (C) A/chicken/Primorsky Krai/1771/2018 (H9N2), (D) A/chicken/Nghe An/08VTC/2018 (H5N1), (E) A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2), (F) A/swine/Irkutsk/155/2017 (H3N2). The X axis—IID in donor ferrets expressed in lg FFU. The Y axis—AID inhaled by recipient ferrets calculated as lg FFU-equivalents based on RT-PCR of the nasal washings from the analytical aerosol filter. Red line indicates AID50 for each strain. Dotted red lines indicate 95% confidence interval. Red triangle indicates AID when recipient ferrets have become infected. Blue solid dots indicate AID when recipient ferrets have not become infected.
Figure 3. Effect of IID in donor ferrets on estimated AID inhaled by recipients, for donor ferrets infected with influenza virus strains: (A) A/California/07/2009 (H1N1)pdm09, (B) A/chicken/Kostroma/1718/2017 (H5N2), (C) A/chicken/Primorsky Krai/1771/2018 (H9N2), (D) A/chicken/Nghe An/08VTC/2018 (H5N1), (E) A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2), (F) A/swine/Irkutsk/155/2017 (H3N2). The X axis—IID in donor ferrets expressed in lg FFU. The Y axis—AID inhaled by recipient ferrets calculated as lg FFU-equivalents based on RT-PCR of the nasal washings from the analytical aerosol filter. Red line indicates AID50 for each strain. Dotted red lines indicate 95% confidence interval. Red triangle indicates AID when recipient ferrets have become infected. Blue solid dots indicate AID when recipient ferrets have not become infected.
Atmosphere 14 00471 g003
Figure 4. Effect of AID received from the air contaminated by donor ferrets on the probability of infection P(AID) in recipient ferrets. The X axis—AID inhaled by recipient ferrets expressed in IID50 units. The Y axis—proportion of infected recipient ferrets expressed as %. Experimental points for different influenza virus strains were approximated using a sigmoid function with the following parameters: AID50 = 16.75 IID50 and σ = 1.18. Curve parameters were calculated based on the experimental points using Origin 8.1 data analysis package.
Figure 4. Effect of AID received from the air contaminated by donor ferrets on the probability of infection P(AID) in recipient ferrets. The X axis—AID inhaled by recipient ferrets expressed in IID50 units. The Y axis—proportion of infected recipient ferrets expressed as %. Experimental points for different influenza virus strains were approximated using a sigmoid function with the following parameters: AID50 = 16.75 IID50 and σ = 1.18. Curve parameters were calculated based on the experimental points using Origin 8.1 data analysis package.
Atmosphere 14 00471 g004
Table 1. Molecular markers of virulence and transmissibility in the HA and NS genes of the studied H5N1, H5N2, and H9N2 viruses.
Table 1. Molecular markers of virulence and transmissibility in the HA and NS genes of the studied H5N1, H5N2, and H9N2 viruses.
StrainExperiment Data:
Increasing Virus Concentration with Increasing Ferret Infection Dose
HANS1 PDZbm
Suppresses the Immune Response
Prevents Early Apoptosis
NS1 Del 80-84 AA
Increases Virulence
Polybasic Cleavage Site in HA, Greatly Increases PathogenicityAdaptation to Humans and Mammals (According to Epidemiological Data)NS1 N200S (N205S) Together with T47A (T48A) in NS2 Are Associated with Decreased Antiviral Response in Mammals
A/chicken/Nghe An/08VTC/2018 (H5N1)Yes (increase)158DYes (ESEV)YesYesNoNS1 S205N, NS2 A47T
Do not reduce antiviral response
A/chicken/Kostroma/1718/2017 (H5N2)Yes (increase)158NYes (ESEV)NoYesNoNS1 N205I, NS2 T47S
Do not reduce antiviral response
A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2)
No (no change)145GNoNoNoNoNS1 N205S, NS2 T47A,
Decreased antiviral response
A/chicken/Primorsky Krai/1771/2018 (H9N2)No (significant dependency is not defined)145NNoNoNoNoNS1 N205S, NS2 T47A
Decreased antiviral response
A/swine/Irkutsk/155/2017(H3N2)No (decrease) NoNoNoAdaptation to pigsNS1 205N, NS2 A47T
Do not reduce antiviral response
A/California/07/2009 (H1N1)pdm09Yes (increase) NoNoNoStrong adaptation to humanNS1 205N, NS2 A47T
Do not reduce antiviral response
Table 2. Fractional composition of aerosol particles and airborne viral load (FFU-equivalent) in the air of aerobiology chamber at various time periods after intranasal infection of donor ferrets with 10 and 106 FFU of avian influenza A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2) virus.
Table 2. Fractional composition of aerosol particles and airborne viral load (FFU-equivalent) in the air of aerobiology chamber at various time periods after intranasal infection of donor ferrets with 10 and 106 FFU of avian influenza A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2) virus.
StrainDose, FFUAerosol Particle Concentration UnitsConcentration of Aerosol ParticlesConcentration of Virus Particles, FFU-Equivalent/L
0.3–1.0 µm1.0–5.0 µm>5.0 µm
Day after InfectionDay after InfectionDay after InfectionDay after Infection
246246246246
A/chicken/Thanh Hoa/V1S5VTC/2020 (H9N2)10Pcs/L53,003
± 13,512
38,260
± 16,211
63,787
± 32,449
5237
± 1454
4058
± 1842
6437
± 3639
218
± 84
111
± 39
62
± 25
ND*0.01200.0030
%90.7
± 23.1
90.2
± 38.2
90.8
± 46.2
9
± 2.5
9.6
± 4.3
9.2
± 5.2
0.37
± 0.14
0.26
± 0.09
0.09
± 0.03
106Pcs/L23,315
± 8139
20,598
± 10,231
16,396
± 9955
2172
± 834
1711
± 903
1306
± 900
64
± 31
27
± 12
26
± 12
0.00090.00110.0010
%91.2
± 31.9
92.2
± 45.8
92.5
± 56.1
8.5
± 3.3
7.7
± 4
7.4
± 5.1
0.25
± 0.12
0.12
± 0.05
0.15
± 0.07
Unchallenged ferrets (control)0Pcs/L33,761 ± 51373288 ± 580644 ± 210ND*ND*ND*
%89.6 ± 13.68.7 ± 1.51.7 ± 0.6
ND*—not detected, below detection limit (<0.00017 FFU-equivalent/L).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gudymo, A.; Onkhonova, G.; Danilenko, A.; Susloparov, I.; Danilchenko, N.; Kosenko, M.; Moiseeva, A.; Kolosova, N.; Svyatchenko, S.; Marchenko, V.; et al. Quantitative Assessment of Airborne Transmission of Human and Animal Influenza Viruses in the Ferret Model. Atmosphere 2023, 14, 471. https://doi.org/10.3390/atmos14030471

AMA Style

Gudymo A, Onkhonova G, Danilenko A, Susloparov I, Danilchenko N, Kosenko M, Moiseeva A, Kolosova N, Svyatchenko S, Marchenko V, et al. Quantitative Assessment of Airborne Transmission of Human and Animal Influenza Viruses in the Ferret Model. Atmosphere. 2023; 14(3):471. https://doi.org/10.3390/atmos14030471

Chicago/Turabian Style

Gudymo, Andrey, Galina Onkhonova, Alexey Danilenko, Ivan Susloparov, Natalia Danilchenko, Maxim Kosenko, Anastasia Moiseeva, Natalia Kolosova, Svetlana Svyatchenko, Vasily Marchenko, and et al. 2023. "Quantitative Assessment of Airborne Transmission of Human and Animal Influenza Viruses in the Ferret Model" Atmosphere 14, no. 3: 471. https://doi.org/10.3390/atmos14030471

APA Style

Gudymo, A., Onkhonova, G., Danilenko, A., Susloparov, I., Danilchenko, N., Kosenko, M., Moiseeva, A., Kolosova, N., Svyatchenko, S., Marchenko, V., Nhai, T. T., Kuznetsov, A., Gavrilova, E., Maksyutov, R., & Ryzhikov, A. (2023). Quantitative Assessment of Airborne Transmission of Human and Animal Influenza Viruses in the Ferret Model. Atmosphere, 14(3), 471. https://doi.org/10.3390/atmos14030471

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

Article Metrics

Back to TopTop