key: cord-0000112-wm4c7rk6
authors: Evseenko, Vasily A; Bukin, Eugeny K; Zaykovskaya, Anna V; Sharshov, Kirill A; Ternovoi, Vladimir A; Ignatyev, George M; Shestopalov, Alexander M
title: Experimental infection of H5N1 HPAI in BALB/c mice
date: 2007-07-27
journal: Virol J
DOI: 10.1186/1743-422x-4-77
sha: 33e17e2a998f1757142ba1df9058e67f360688c4
doc_id: 112
cord_uid: wm4c7rk6

BACKGROUND: In 2005 huge epizooty of H5N1 HPAI occurred in Russia. It had been clear that territory of Russia becoming endemic for H5N1 HPAI. In 2006 several outbreaks have occurred. To develop new vaccines and antiviral therapies, animal models had to be investigated. We choose highly pathogenic strain for these studies. RESULTS: A/duck/Tuva/01/06 belongs to Quinghai-like group viruses. Molecular markers – cleavage site, K627 in PB2 characterize this virus as highly pathogenic. This data was confirmed by direct pathogenic tests: IVPI = 3.0, MLD(50 )= 1,4Log10EID(50). Also molecular analysis showed sensivity of the virus to adamantanes and neuraminidase inhibitors. Serological analysis showed wide cross-reactivity of this virus with sera produced to H5N1 HPAI viruses isolated earlier in South-East Asia. Mean time to death of infected animals was 8,19+/-0,18 days. First time acute delayed hemorrhagic syndrome was observed in mice lethal model. Hypercytokinemia was determined by elevated sera levels of IFN-gamma, IL-6, IL-10. CONCLUSION: Assuming all obtained data we can conclude that basic model parameters were characterized and virus A/duck/Tuva/01/06 can be used to evaluate anti-influenza vaccines and therapeutics.

Influenza A (H5N1) virus now becomes a real threat for humans. Since 1997, when first human case of H5N1 HPAI had been reported, more than 317 people were infected and 191 died [1] . Before 2005 attention was attracted to Thailand, Vietnamese and Indonesian viruses. In the beginning of 2005 outbreak on Quinghai lake occurred [2] . Later "Quinghai-like" viruses spreaded to most part of Russia, European countries and Africa and caused numerous outbreaks. Only in Russia more than 1 million of different species and sorts of poultry died and been slaughtered [3] . Confirmed cases in Azerbaijan, Egypt, Iraq, and Turkey was caused by Quinghai-like viruses. Earlier HPAI viruses were investigated in mice [4, 5] and murine models were successively used for reverse genetics made influenza vaccines [6] . It was shown that H5N1 HPAI viruses could have different pathogenicity for mice [7] . Several molecular markers were choused to explain differences. Multibasic cleavage site with 627K in PB2 designate to highly pathogenic phenotype for mice. Also important role of pulmonary cytokines elevation was highlighted [8] . Combination of adaptation for wild waterfowl and high virulence for mammals makes Quinghai-like viruses presumably pandemic. Also, in future, because of ability for rapid spreading for long distances, this group of viruses can appear in North and South America and cause outbreaks.

Human disease caused by HPAI viruses can be characterized as acute viral pneumonia aggravated by ARDS, toxic shock and multiple organ failure. System dysfunction mediated by hypercytokinemia and high viral load [9] . To be ready for new influenza pandemy it is necessary to use animal models, in vaccine and antivirals studies, which most closely reflect human disease. Isolates from FRSI SRC VB "VECTOR" repository which were characterized previously were examined for MLD 50 , molecular markers of pathogenicity, sensitivity to amantadines and neuraminidase inhibitors, to be candidates for murine model. Among the investigated isolates A/duck/Tuva/01/06 has best features to be used.

Genes of A/duck/Tuva/01/06 were sequenced and analyzed for molecular markers of pathogenicity. Also phylogenetic analysis was performed. Results are presented in figure 1 . A/duck/Tuva/01/06 belongs to group of Qinghailike viruses. HA contains 5 polybasic aminoacids (PQGRRKKKR↓GL) in cleavege site of HA [15] . The receptor binding domen can be characterized as "avian" [16] . High pathogenicity to mammals in general correlates with presence of 627K in PB2 [17] .

The analysis of non-structural protein 1 (NS1) which also could be contributed for high virulence of H5N1 viruses revealed deletion of 5 amino acids similar to those in H5N1 viruses of genotype Z which could be contributed to increased expression of TNF-α and IP-10 protein in primary human macrophages [18] . A/duck/Tuva/01/06 con-tained Glu 92 in the NS1 and contained "avian-like" PDZdomain ligand ESEV [19] . It was shown that the most recent H5N1 strains isolated in Southeast Asia were resistant to amantadine and rimantadine; antiviral drugs targeted the M2 ion channels of influenza A viruses [20, 21] . It was also reported about Oseltamivir resistant H5N1 viruses isolation from humans [22, 23] . To determine the potential sensitivity of studied H5N1 viruses to these antivirals, amino acid sequences of the M2 and NA proteins were analyzed.

Variants of influenza A viruses resistant to amantadine possessed amino acid substitutions at one of 5 residues (26, 27, 30, 31, and 34) in the M2 protein [24, 25] . Sequence analysis did not reveal any mutations associated with resistance to amantadine. Thus all A/duck/Tuva/01/ 06 is potentially sensitive to this class of antiviral agents. Amino acid residues 119, 274, 292 and 294 in the NA protein (numbering according to the HA of H2 subtype) are crucial for the sensitivity of influenza A viruses to neuraminidase inhibitors [26] ; substitution H 274 →Y in the NA conferred resistance to Oseltamivir was observed in clinical H5N1 isolates [25, 26] . Sequence comparison of the NA protein of A/duck/Tuva/01/06 aligned with the NA of N2 subtype of A/Wuhan/359/95 (H3N2) influenza virus showed phenotype potentially sensitive to neuraminidase inhibitors.

A/duck/Tuva/01/06 showed wide cross-reactivity with sera against H5N1 HPAI viruses isolated earlier in South-Eastern Asia. HI results can be found in table 1. These features persuade to use this virus in studies of vaccines made from various H5N1 influenza viruses.

First MID 50 and MLD 50 for A/duck/Tuva/01/06 were determined (table 2) . To determine mean time to death (m.t.d) and 2C. In several cases (9 animals totally) the disease was complicated by severe intestine atony, which can independently lead to death or by pressuring on diaphragm can intensify respiratory failure.

We also determined virus titers in several organ tissues. As it was expected the highest titers was observed in lungs -5,3 log EID 50 . Brain titers were also high -3,4 log EID 50 . In spleen, liver and kidney tissues virus titers were lower then 1 logEID 50 and considered not significant.

We investigated the involvement of several cytokines in immunopathogenesis of experimental H5N1 HPAI infection in mice. Results of ELISA technique revealed alteration of expression both pro-inflammatory and antiinflammatory cytokines after the challenge (figure 3). In general, the most marked changes of cytokine levels were observed before the death of mice.

The minimal concentration of IFN-γ was detected on day 5 (14.3 ± 10.8 pg/ml), however, its levels enlarged about 8-fold (256 ± 27 pg/ml) during the course of the infection when compared with uninfected animals. On days 3 and 5 systemic production of TNF-α was below the detection limit of the assay. A peak was reached on day 7 by the cytokine (24 ± 3.2 pg/ml) and its levels remained elevated on day 8. Interestingly, concentrations of IL-1β in mice after the challenge were significantly lower in comparison with the constitutive expression of the mediator in intact animals. An abrupt decrease of IL-1β was detected on day 3 post infection, but was followed by step increase from day 5. After the 2.5-fold enlargement on day 3 the levels of IL-6 decreased dramatically on day 5, and the highest levels of the cytokine were determined at the end of observation period (133 ± 12 pg/ml). The constitutive production of IL-10 was undetectable. The dynamics of IL-10 showed a gradual growth with the maximum level (92.1 ± 6.0 pg/ml) reached before the death of mice. We observed statistically significant increase of IL-12 after the challenge. Concentrations of the cytokine retained constant in infected mice, except the unexpected decline occurred on day 7. The expression of IL-18 could not be detected throughout the entire period of observation.

Until 2005 avian influenza was regional problem of sev- [4] and question "why had only some wild waterfowl died?" is still unclear. Most of the outbreaks in Russia associated with wild birds. The same time viruses adapted to wild birds are extremely pathogenic for poultry and mice. This "competitive advantage" makes Quinghailike viruses most probable candidate to be precursor for new pandemic influenza virus. At the same time pathogenesis of different (phylogenetical clades) HPAI reveal common causes. The principal causes of rapid mice death after infecting with HPAI are primary viral pneumonia, ARDS, lesions of central nervous system and multiple organ failure. Our data suggest that A/duck/Tuva/01/06 strain of HPAI caused lethal pneumonia and spread systemically to the brain in BALB/c mice. Lesion of respiratory epithelium and following an activation of monocytes/macrophages results in a release of proinflammatory cytokines (TNF-α, IL-6) which are a hallmark of ARDS in murine model [27] . Despite powerful anti-influenza virus effects of TNF-α in lung tissue, as it was described previously [28] , we consider that elevated production of the cytokines seems to be crucial in the pathogenesis of HPAI infection. Moreover, it was shown that lethal H5N1 viruses are resistant to antiviral effects of interferons and TNF-α [29] . Virus-induced overexpression of TNF-α as well as high IFN-γ lead to activation of endothelium and imbalance in blood coagulation system [30] . This may explain the hemorrhagic syndrome as observed in some of animals. To pay attention that IL-12 is a potent inducer of IFN-γ synthesis by blood mononuclear cells [31] , we concluded the same cytokines hyperproduction reflects macrophage overactivation and subsequent hypercytokinemia. This cascade of events Phylogenetic tree based on full length sequencesof HA Figure 2 Phylogenetic tree based on full length sequencesof HA. Nucleotide sequences were analyzed by using the neighborjoining method with 500 bootstraps. The phylogenetic tree was rooted to the HA gene of A/goose/Guangdong/1/96 (H5N1) virus.

including inflammatory mediator production, changes in blood coagulation system and microvascular permeability was denoted as systemic inflammatory response syndrome (SIRS) [32] . On the other hand, we proposed that the prominent production of IL-10 from the early stages of the experimental HPAI infection was the compensatory response to overproduction of proinflammatory cytokines such as TNF-α, IL-6 and IL-12. However, the role of IL-10, which principle function seems to be containment and eventual termination of inflammation [33] , in HPAI pathogenesis is unclear. Also there is an uncertain discrepancy between undetectable expression of IL-18 and high levels of other Th1-cytokines (IFN-γ and IL-12).

Summing up, in our study BALB/c mice infected with HPAI, strain A/duck/Tuva/01/06, appeared to be able to produce the innate immune response, which culminated to the development of shock and subsequent multiple organ failure. The main characteristics of our model are comparable to the previously described fatal cases of H5N1 influenza in humans [10, 11] . Proposed model reflects lesions not only same organs but also mediating levels of some (IFN-γ, IL-6, IL-10) cytokines in terminal conditions.

The implication of different cytokines in immunopathogenesis of experimental HPAI is beyond question. But to understand exact mechanisms, which determine the disease outcome, further experiments remain to be done. 

All experiments were performed in BSL 3+ facilities of FSRI SRC VB "Vector" of Rospotrebnadzor licensed for working with highly pathogenic avian influenza viruses.

Stock of A/duck/Tuva/01/06 was produced in 9 days-old chicken embryos. Allantoic fluid was aliquoted and stored at -80°C. The infectivity of stock viruses was determined in 10 days-old embryonated chicken eggs; titers were calculated by the method of Reed and Muench [10] and were expressed as log 10 of 50% egg infective dose (EID 50 ) in 1 ml of allantoic fluid.

Viral RNA was isolated from virus-containing allantoic fluid with the RNeasy Mini kit (Qiagen, Valencia, CA) as specified by the manufacturer. Uni-12 primer was used for reverse transcription. PCR was performed with a set of primers specific for each gene segment of Influenza A virus [11] . PCR products were purified with the QIAquick PCR purification (Qiagen).

Sequencing was done with Beckman Coulter Genom-eLab™ Methods development kit Dye terminator Cycle Sequencing according instructions of manufacturer. Primers for sequence were obtained from E. Hoffman (SJCRH, Memphis, TN). Sequence products were analyzed on automatic sequence analyzer Beckman Coulter CEQ2000.

Phylogenetical analysis was done on HA full gene sequence DQ861291 using MEGA 2.1 software. Phylogenetical tree was built by Neighbor-Joining method; matrix of distances was counted with p-distance algorithm. Reliability of clades was checked with bootstrap analysis with 500 replications. Other genes in GenBank DQ861291-DQ861295.

Cross-reaction of A/duck/Tuva/01/06 was defined by hemagglutination inhibition test (HI) with 0.5% CRBC [12] with a panel of antisera against H5N1 HPAI. 

Six-week-old inbred male BALB/c mice (vivarium of FRSI SRC VB "Vector"). Animals were placed to individual cages with food and water available ad libitum. To determine the MLD 50 and MID 50 , mice were anaesthetized by diethyl ether inhalation and infected intranasally with 50 µl 10-fold serial dilutions of allantoic fluidin PBS (pH 7,2). Each group contained 10 animals. Animals were observed daily for 15 days for mortality (MLD 50 ) or sacrificed on day 5 after the challenge with following virus detection in the lungs by inoculation of 10 days-old embryonated chicken eggs (MID 50 ). MLD 50 and MID 50 were calculated by the method of Reed and Muench. Animals from group where 1MLD 50 had been observed were taken to determine virus titers in lung, spleen, kidneys, and liver and brain tissues. Mind time to death (m.t.d) was calculated as previously described [13] . Pathogenicity to chickens was determined by IVPI test [14] . All animal studies were performed according protocols approved by Animal Care & Use committee of FSRI SRC VB "Vector".

To determine IFN-γ, TNF-α, IL-6, IL-10, IL1-β, IL-12 we use ELISA R& D Systems kits (Minneapolis, MN, USA). Serum levels of IL-18 were measured using commercial Mouse IL-18 ELISA test kit (MBL, Nagoya, Japan). Detection limits were as follows: TNF-α, less then 5,1 pg/ml; IL1-β, 3,0 pg/ml; IL6, 3,1 pg/ml; IL10, 4,0 pg/ml; IL-18, 25 pg/ml. Sera was taken on 0,3,5,7,8 days and aliquots and stored -80°C upon usage. Day 8 was chosen because m.t.d defined earlier in the work was 8,19 ± 0,18 days. Statistics was performed with Student t-test. Values p < 0,05 considered to be reliable.

World Health Organization: Epidemic and Pandemic Alert and Response (EPR)

HighlyPathogenic H5N1 Influenza Virus Infection in Migratory Birds

Avian influenza in Siberia -2005: Laboratory and epidemiological studies, antiepidemic measures during the epizooty of avian influenza in poultry in Siberian and Ural Federal Regions of the Russian Federation

Neurovirulence in mice of H5N1 influenza virus genotypes isolated from Hong Kong poultry in 2001

Avian influenza (H5N1) viruses isolated from humans in Asia in 2004 exhibit increased virulence in mammals

Efficacy of H5 influenza vaccines produced by reverse genetics in a lethal mouse model

Avian influenza (H5N1) viruses isolated from humans in Asia in 2004 exhibit increased virulence in mammals

Pathogenesis of Hong Kong H5N1 influenza virus NS gene reassortants in mice: the role of cytokines and B-and T-cell responses

Fatal outcome of human influenza A (H5N1) is associated with high viral load and hypercytokinemia

A simple method for estimating fifty percent endpoints

Universal primer set for the full-length amplification of all influenza A viruses

US Department of Health, Education and Welfare Immunology Series 6

Statistical methods in microbiological studies

A color atlas and text on avian influenza. Papi Editore

Influenza: emergence and control

Structure and receptor specificity of the hemagglutinin from an H5N1 influenza virus

The polymerase complex genes contribute to the high virulence of the human H5N1 influenza virus isolate A/ Vietnam/1203/04

H5N1 influenza: a protean pandemic threat

Large-scale sequence analysis of avian influenza isolates

Genesis of a highly pathogenic and potentially pandemic H5N1 influenza virus in eastern Asia

Detection of amantadine-resistant variants among avian influenza viruses isolated in North America and Asia

Avian flu: isolation of drug-resistant H5N1 virus

Oseltamivir resistance during treatment of influenza A (H5N1) infection

Themolecular basis of the specific anti-influenza action of amantadine

Influenza virus M2protein has ion channel activity

Molecular mechanisms of influenza virus resistance to neuraminidase inhibitors

Acute Respiratory Distress Syndrome Induced by Avian Influenza A (H5N1) Virus in Mice

Tumor necrosis factor alpha exerts powerful anti-influenza virus effects in lung epithelial cells

Lethal H5N1 influenza viruses escape host anti-viral cytokine responses

The interactions between inflammation and coagulation

Regulation of human cytolytic lymphocyte responses by interleukin-12

Toward an epidemiology and natural history of SIRS (systemic inflammatory response syndrome)

Interleukin-10 and the interleukin-10 receptor

This work was supported by Bio Industry Initiative (BII) of the US department of State grant ISTC #3436.

The author(s) declare that they have no competing interests.

VE carried out molecular genetic analysis, performed animal studies, design of experiments and drafted manuscript. EB performed immunoassays and obtained data analysis. AZ participated in animal studies. KS assisted in animal studies. VT was responsible for sequence. GI participated in study design and coordination. AS carried out coordination. All authors read and approved the final manuscript.