| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
1 Department of Microbiology, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok 10700, Thailand
2 Department of Immunology, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok 10700, Thailand
3 Department of Pathology, Faculty of Medicine Siriraj Hospital, Mahidol University, Bangkok 10700, Thailand
Correspondence
Prasert Auewarakul
sipaw{at}mahidol.ac.th
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Published online ahead of print on 12 May 2008 as DOI 10.1099/vir.0.2008/002469-0.
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
). Although, most sequence variations are neutral and offer no competitive advantage, the quasispecies provides multiple variants, which can be readily selected if there is a change in selective pressure. Transmission of the virus to a new host species provides new selective pressures and can result in the expansion of the ‘best fit’ minor variants for adaptation to the new environment.
In general, human and avian influenza A viruses differ in their recognition of host cell receptors. Human influenza viruses preferentially recognize 
-2,6-linked sialic acid, while avian influenza viruses recognize -2,3-linked sialic acid (Rogers & Paulson, 1983). However, the highly pathogenic avian influenza viruses subtype H5N1 can transmit directly from avian species to humans (Subbarao et al., 1998; Tran et al., 2004). Even though the H5N1 viruses can infect and cause severe disease in humans, they do not bind the -2,6-linked sialic acid receptor with high affinity (Ha et al., 2001; Stevens et al., 2006). This property is believed to be one of the major factors that prevent the H5N1 virus from transmitting efficiently amongst humans and causing a pandemic. Amino acid substitutions in the haemagglutinin (HA) gene can lead to the altered receptor-binding preference of the virus from -2,3-linked to -2,6-linked sialic acid (Auewarakul et al., 2007; Stevens et al., 2006; Yamada et al., 2006). This would enable avian H5N1 viruses to recognize human-type host cell receptors and could potentially enable the virus to transmit efficiently within the human population and cause a catastrophic pandemic. Therefore, it is extremely important to monitor the viral changes that may lead to the emergence of pandemic viruses.
The earliest possible detection of selected mutants is by looking at sequences within viral quasispecies before they expand and become dominant virus. In order to detect mutants that might have altered phenotypes we studied the viral sequences at the level of quasispecies. We conducted a study where the viral sequence was directly amplified, cloned and sequenced from a nasopharyngeal aspirate or tissue specimens. The specimens were obtained from fatal human cases in Thailand.
Selection at the protein level can be measured by 
(dN/dS), in which dN=non-synonymous substitution rate (non-synonymous changes per non-synonymous site) and dS=synonymous substitution rate (synonymous changes per synonymous site). If amino acid changes provide better fitness, the mutations will be fixed at a higher rate than synonymous mutations. This results in dN>dS and >1. Originally, was calculated as an average for the whole gene, which does not allow sensitive detection of individual amino acid residue under positive selection. Subsequently codon-based models that allow the ratio to vary amongst sites were developed (Nielsen & Yang, 1998; Yang & Nielsen, 2000). These models describe ratio distribution amongst sites: M0 assumes one constant for all sites; M3 classifies sites into discrete classes with different ; M7 allows to vary according to a distribution that represents negative or neutral selection; M8 adds on top of M7 a discrete class for sites with positive selection (>1). Likelihood ratio test is used as a statistical test of goodness-or-fit to compare the two models and test whether the more complex model, e.g. M8 or M3, fits the dataset significantly better than the simpler model, e.g. M7 or M0. An empirical Bayes approach is then used to calculate the posterior probability that each site is from a particular site class, and sites with high posterior probabilities coming from the class with >1 are inferred to be under positive selection (Yang et al., 2005).
|
METHODS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
). Patient A was a 5-year-old boy who had progressive viral pneumonia that led to respiratory failure and death by 12 days after the onset of illness. In December 2005, he developed a fever, stomach ache, nausea and vomiting and was admitted to the hospital 10 days later when he developed dyspnea and a chest radiography showed patchy infiltration at the right middle lobe. The symptom of dyspnea worsened to respiratory failure and the pulmonary infiltration spread to both lungs on the following day. The diagnosis of avian influenza was suspected on day 12. The use of the antiviral drug oseltamivir was started and the patient expired on the same day. The patient was not known to have had direct contact with any sick or dying birds, but he played in the yard where the birds were often present. Nasopharyngeal aspirate was collected on day 12 after the onset of illness.
Autopsy was conducted in the second (patient B) and third patients (patient C). Patient B, previously reported (Uiprasertkul et al., 2007), was a 48-year-old man who had progressive viral pneumonia in October 2005. He had a fever, cough, runny nose, myalgia and chest pain at the onset of illness. Dyspnea developed on day 2 of the illness and a chest radiograph showed interstitial infiltrations at the right upper and left middle lung fields and a mass-like infiltration at the right middle lung field. The diagnosis of avian influenza was suspected on day 4 of the illness after a history of direct contact with dying chickens was revealed. Respiratory secretions were then sent to national laboratories and they were confirmed positive for influenza (H5N1) virus. The patient died on day 6 of the illness.
Patient C, previously reported (Uiprasertkul et al., 2005), was a 6-year-old boy who had progressive viral pneumonia in January 2004. He was initially treated with multiple broad-spectrum anti-microbial agents. Virological diagnosis of H5N1 infection was made on day 7 of the illness. After oseltamivir became available in Thailand, he was treated on day 15 of his illness with this agent until he died. He was also treated with methylprednisolone on day 15 until death and with granulocyte colony-stimulating factor for leukopenia from day 5 to 10 of the illness. The patient died on day 17 of the illness.
The use of the patients' specimens was approved by the Ethics Committee of the Faculty of Medicine Siriraj Hospital.
Viral RNA, cloning and quasispecies analysis.
For patient A, total RNA from the nasopharyngeal specimen was extracted according to the manufacturer's protocol (QIAmp RNA mini kit; Qiagen). For patients B and C, total RNAs were extracted by using Trizol from paraffin-embedded blocks of lung and intestine tissue samples and then purified using Qiagen RNAeasy kit according to the manufacturer's instructions.
A fragment of the HA gene covering the receptor-binding site (nt 413–905) was amplified from RNA extracted from the specimen by using the high fidelity enzyme Pfu (Promega) and the primers HHAf2 (5'-GGTCCAGTCATGAAGCCTCA-3') and HA-H5r12 (5'-TTTATCGCCCCCATTGGAGT-3'). The PCR product was cloned into pGEM T-Easy. One hundred clones of each sample were picked up and sequenced.
The selective pressures acting on the receptor region were estimated by using the CODEML program in the PAML package. We used models M7 and M8, where M7 contains 10 categories to describe amongst sites, all constrained to be <1; M8 differs from M7 only in that it estimates for an extra class of sites (p10) at which can be >1 (Yang, 1997). Models were compared using a likelihood ratio test and the Bayes Empirical Bayes (BEB) method was used for a posteriori estimation of individual codons under positive selection (Yang et al., 2005).
|
RESULTS |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
) within this sequenced region spanning 143 aa residues (128–270). At these positions, the sequence majority in all the samples contains similar amino acids as in the consensus sequence of all human H5N1 viruses available in the GenBank database, except for positions 133 and 138 where mutant amino acids were found in higher frequencies. Total numbers of synonymous and non-synonymous substitutions, as well as lists of non-synonymous substitutions in all the RNA samples are shown in Table 2. Some of these mutations with low frequencies were not picked up as positive-selection site by the BEB analysis. Amongst the positive-selection sites, six sites were repeatedly found when each patient was analysed individually. These sites were 133, 138, 161, 186, 222 and 227 (H3 numbering system). Of these all but the 133 and 161 positions are in the known receptor-binding domain. The receptor-binding site of HA at the tip of HA1 globular domain is composed of three secondary structure elements: the 190 helix (residues 190–198), the 130 loop (residues 135–138) and the 220 loop (residues 221–228), forming the sides of each site; and the base made up of the conserved residues Tyr98, Trp153, His183 and Tyr195 (Skehel & Wiley, 2000). Although, the positions 133 and 230 are not in the receptor-binding domain, they flank the 130 and 220 loops and they showed positive selection in our analyses, suggesting that mutations at these sites might also contribute to the receptor-binding adaptation. A138V, N186K and S227N mutations were previously reported to confer -2,6-linked sialic acid binding to H5N1 virus (Auewarakul et al., 2007; Gambaryan et al., 2006; Shinya et al., 2005; Yamada et al., 2006). The mutant sequences in our analyses are similar to these mutations at positions 138 and 227, whereas our mutation at position 186 is N186D. Amongst the 17 positive-selection sites that were identified in our study, eight sites are not known to be related to the receptor-binding domain. But, these sites all showed only low mutation frequencies. In other words, all strong positive-selection sites with high mutation frequencies are related to the receptor-binding domain.
|
|
). The higher frequency of positive-selection sites in patient C may be related to the fact that this patient died on day 17 of the illness, while the samples from the other two patients were collected earlier in the course of the illness. The longer duration of infection provided a longer period under the selective pressure and might cause the virus to gain more adaptation. In patient C, there were markedly more positive-selection sites that are not related to the receptor-binding domain, including the N-linked glycosylation site at position 158. This suggests that other selective pressures, such as the immune response, might be involved in late phases of the disease.
When each type of tissue from patients B and C was analysed, four and seven positive-selection sites were identified in the lung and intestine, of which two and four sites are related to the receptor-binding domain, respectively. (Table 1). There are some differences in the positive-selection sites from different tissues. In particular, the S227N mutation was found in high frequency (36.7 %) in the intestine of patient C and in low frequency (3 %) in the intestine of patient B, while it was absent in the lung tissue from both cases. This suggested a strong selection for this mutation and a compartmentalization of the viral population within the patients.
Most of the identified positive-selection sites concentrated in the N-terminal 2/3 part of the sequences. Sixteen sites were identified in the region covering the receptor-binding domain from residue 128 to 230, whereas only two sites were identified in the rest of the sequences (residue 231–270).
Mutations that have been reported to alter receptor-binding specificity of H5 are A138V, N186K, Q196R, S227N, Q226L and G228S (Auewarakul et al., 2007; Gambaryan et al., 2006; Shinya et al., 2005; Stevens et al., 2006; Yamada et al., 2006). Our analyses provide new candidate mutations that may affect the receptor-binding specificity. Whether these mutations have functional consequences is under investigation.
|
DISCUSSION |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
, 2007), it is possible that they constitute an unknown epitope.
The observed differences in the viral quasispecies in lung and intestine suggest that there might be a compartmentalization of viral infection and that the selective pressure might be different among different tissues. Human lung has been shown to contain -2,3-linked sialic acid in alveolar epithelial cells, whereas human intestinal epithelium lacks this receptor. In human intestine the -2,3-linked sialic acid receptor was identified only on neurons (Yao et al., 2008). It is not clear whether the difference in the sialic acid receptor distribution between the two tissues contributed to the different selection of viral sequences.
Our data demonstrate adaptation of the receptor-binding domain of H5N1 virus in infected human tissues. If allowed to be transmitted further to other human hosts, the mutants would be likely to be selected further and expand, and eventually cause emergence of a potential pandemic virus. Understanding the adaptation is therefore of upmost importance. Adaptation of avian influenza virus to human host involves multiple mechanisms. However, the receptor usage preference is likely to be a major step in the adaptation process. Although, there have been reports on mutations that altered receptor-binding specificity of H5, those reported mutations only conferred partial switching from -2,3-linkage tropism to dual tropism (Auewarakul et al., 2007; Gambaryan et al., 2006; Stevens et al., 2006; Yamada et al., 2006). It is likely that these mutations are not sufficient and full switching to -2,6-linkage tropism is probably needed for an efficient transmission in the human population. Such mutation that can cause a complete switching of H5 is not known. The fact that our analyses could pick up mutations that are known to change the receptor-binding property of H5N1 viruses indicates that they can be used to screen and search for mutations with significant functional effects. Our analyses offer an approach to find candidate mutations, which should be studied further for determining functionality. Finding mutations with pandemic potential before the actual emergence of such viruses will provide genetic markers for vigilant monitoring, which will hopefully help us to avoid the pandemic.
|
ACKNOWLEDGEMENTS |
|---|
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMETHODSRESULTSDISCUSSIONREFERENCES |
|---|
Eigen, M. (1996). On the nature of virus quasispecies. Trends Microbiol 4, 216–218.[CrossRef][Medline]
Gambaryan, A., Tuzikov, A., Pazynina, G., Bovin, N., Balish, A. & Klimov, A. (2006). Evolution of the receptor binding phenotype of influenza A (H5) viruses. Virology 344, 432–438.[CrossRef][Medline]
Ha, Y., Stevens, D. J., Skehel, J. J. & Wiley, D. C. (2001). X-ray structures of H5 avian and H9 swine influenza virus hemagglutinins bound to avian and human receptor analogs. Proc Natl Acad Sci U S A 98, 11181–11186.
Kaverin, N. V., Rudneva, I. A., Ilyushina, N. A., Varich, N. L., Lipatov, A. S., Smirnov, Y. A., Govorkova, E. A., Gitelman, A. K., Lvov, D. K. & Webster, R. G. (2002). Structure of antigenic sites on the haemagglutinin molecule of H5 avian influenza virus and phenotypic variation of escape mutants. J Gen Virol 83, 2497–2505.
Kaverin, N. V., Rudneva, I. A., Govorkova, E. A., Timofeeva, T. A., Shilov, A. A., Kochergin-Nikitsky, K. S., Krylov, P. S. & Webster, R. G. (2007). Epitope mapping of the hemagglutinin molecule of a highly pathogenic H5N1 influenza virus by using monoclonal antibodies. J Virol 81, 12911–12917.
Nielsen, R. & Yang, Z. (1998). Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. Genetics 148, 929–936.
Rogers, G. N. & Paulson, J. C. (1983). Receptor determinants of human and animal influenza virus isolates: differences in receptor specificity of the H3 hemagglutinin based on species of origin. Virology 127, 361–373.[CrossRef][Medline]
Shinya, K., Hatta, M., Yamada, S., Takada, A., Watanabe, S., Halfmann, P., Horimoto, T., Neumann, G., Kim, J. H. & other authors (2005). Characterization of a human H5N1 influenza A virus isolated in 2003. J Virol 79, 9926–9932.
Skehel, J. J. & Wiley, D. C. (2000). Receptor binding and membrane fusion in virus entry: the influenza hemagglutinin. Annu Rev Biochem 69, 531–569.[CrossRef][Medline]
Stevens, J., Blixt, O., Tumpey, T. M., Taubenberger, J. K., Paulson, J. C. & Wilson, I. A. (2006). Structure and receptor specificity of the hemagglutinin from an H5N1 influenza virus. Science 312, 404–410.
Subbarao, K., Klimov, A., Katz, J., Regnery, H., Lim, W., Hall, H., Perdue, M., Swayne, D., Bender, C. & other authors (1998). Characterization of an avian influenza A (H5N1) virus isolated from a child with a fatal respiratory illness. Science 279, 393–396.
Tran, T. H., Nguyen, T. L., Nguyen, T. D., Luong, T. S., Pham, P. M., Nguyen, V. C., Pham, T. S., Vo, C. D., Le, T. Q. & other authors (2004). Avian influenza A (H5N1) in 10 patients in Vietnam. N Engl J Med 350, 1179–1188.
Uiprasertkul, M., Puthavathana, P., Sangsiriwut, K., Pooruk, P., Srisook, K., Peiris, M., Nicholls, J. M., Chokephaibulkit, K., Vanprapar, N. & Auewarakul, P. (2005). Influenza A H5N1 replication sites in humans. Emerg Infect Dis 11, 1036–1041.[Medline]
Uiprasertkul, M., Kitphati, R., Puthavathana, P., Kriwong, R., Kongchanagul, A., Ungchusak, K., Angkasekwinai, S., Chokephaibulkit, K., Srisook, K. & other authors (2007). Apoptosis and pathogenesis of avian influenza A (H5N1) virus in humans. Emerg Infect Dis 13, 708–712.[Medline]
Yamada, S., Suzuki, Y., Suzuki, T., Le, M. Q., Nidom, C. A., Sakai-Tagawa, Y., Muramoto, Y., Ito, M., Kiso, M. & other authors (2006). Haemagglutinin mutations responsible for the binding of H5N1 influenza A viruses to human-type receptors. Nature 444, 378–382.[CrossRef][Medline]
Yang, Z. (1997). PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 13, 555–556.
Yang, Z. & Nielsen, R. (2000). Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol 17, 32–43.
Yang, Z., Wong, W. S. & Nielsen, R. (2005). Bayes empirical bayes inference of amino acid sites under positive selection. Mol Biol Evol 22, 1107–1118.
Yao, L., Korteweg, C., Hsueh, W. & Gu, J. (2008). Avian influenza receptor expression in H5N1-infected and noninfected human tissues. FASEB J 22, 733–740.
Received 28 March 2008;
accepted 2 May 2008.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |
