Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Adamantane-Resistant Influenza A Viruses in the World (1902–2013): Frequency and Distribution of M2 Gene Mutations

  • Guoying Dong ,

    dongguoying88@aliyun.com (GD); gj28@ibp.ac.cn (GJ); hehx@ioz.ac.cn (HH)

    Affiliations College of Global Change and Earth System Science, Beijing Normal University, Beijing, China, National Research Center For Wildlife Born Diseases, Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China, National Laboratory of Biomacromolecules, Institute of Biophysics of Chinese Academy of Sciences, Beijing, China

  • Chao Peng,

    Affiliation College of Global Change and Earth System Science, Beijing Normal University, Beijing, China

  • Jing Luo,

    Affiliation National Research Center For Wildlife Born Diseases, Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China

  • Chengmin Wang,

    Affiliation National Research Center For Wildlife Born Diseases, Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China

  • Le Han,

    Affiliation College of Global Change and Earth System Science, Beijing Normal University, Beijing, China

  • Bin Wu,

    Affiliation National Research Center For Wildlife Born Diseases, Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China

  • Guangju Ji ,

    dongguoying88@aliyun.com (GD); gj28@ibp.ac.cn (GJ); hehx@ioz.ac.cn (HH)

    Affiliation National Laboratory of Biomacromolecules, Institute of Biophysics of Chinese Academy of Sciences, Beijing, China

  • Hongxuan He

    dongguoying88@aliyun.com (GD); gj28@ibp.ac.cn (GJ); hehx@ioz.ac.cn (HH)

    Affiliation National Research Center For Wildlife Born Diseases, Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China

Abstract

Adamantanes (amantadine and rimantadine) have been used to prevent and treat influenza A virus infections for many years; however, resistance to these drugs has been widely reported in the world. To investigate the frequency and distribution of M2 gene mutations in adamantane-resistant influenza variants circulated in the world between 1902 and 2013, 31251 available M2 protein sequences from different HA-subtype influenza A viruses (H1–H17) were analyzed and adamantane resistance-associated mutations were compared (L26F, V27A, A30T, A30V, S31N, G34E, and L38F). We find that 45.2% (n = 14132) of influenza A (H1–H17) viruses circulating globally were resistant to adamantanes, and the vast majority of resistant viruses (95%) bear S31N mutations. Whereas, only about 1% have V27A mutations and other mutations (L26F, A30T, G34E, and L38F) were extremely rare (their prevalence appeared to be < 0.2%). Our results confirm that H1, H3, H5, H7, H9, and H17 subtype influenza A viruses exhibit high-level resistance to adamantanes. In contrast, the appearance of adamantane-resistant mutants in H2, H4, H6, H10, and H11 subtypes was rare. However, no adamantane resistance viruses were identified among other HA subtypes (H8, H12–H16). Our findings indicate that the frequency and distribution of adamantane-resistant influenza variants varied among different HA subtypes, host species, years of isolation, and geographical areas. This comprehensive study raises concerns about the increasing prevalence of adamantane-resistant influenza A viruses and highlights the importance of monitoring the emergence and worldwide spread of adamantane-resistant variants.

Introduction

Influenza A virus is a common cause of respiratory infections, infecting humans, birds, pigs, horses, as well as other species [1]. Eighteen Hemagglutinin (HA, H1 to H18) and eleven Neuraminidase (NA, N1 to N11) subtypes of influenza viruses can cause epidemics and pandemics through antigenic drift and antigenic shift, respectively [1]. Considering the significant threat that influenza A viruses pose to global economy and human health, preparedness for a potential influenza pandemic is a global priority [2].

The primary means to protect humans against influenza virus infection is by vaccination. However, in the absence of an effective and widely available vaccine, protection from a pandemic and/or long-term endemicity would rely largely on the prophylactic and therapeutic properties of antiviral treatment. Therefore, effective anti-influenza drugs offer a valuable addition to the use of vaccines, particularly for those not caused by drug-resistant viruses [3]. Currently, two classes of anti-influenza drugs, adamantanes (amantadine and its methyl derivative rimantadine) and neuraminidase inhibitors (oseltamivir and zanamivir), are used for prophylaxis and treatment of influenza A infections [4]. These FDA-approved drugs are effective for influenza A viruses, including highly pathogenic H5N1 viruses [5]. Their therapeutic efficacy, however, has come into question with the emergence of drug-resistant viruses. Adamantane drugs, also known as M2 channel blockers, could block the ion channel formed by the M2 protein encoded by the M gene of influenza A viruses and inhibit early stages of virus replication [6]. As adamantane derivatives, amantadine (SymmetrelTM) and rimantadine (FlumadineTM) have been used successfully for influenza A virus infection for more than 30 years because of their wide availability and low cost [7, 8]. The effectiveness of prophylaxis with these drugs varied between 80% and 90%, and their use can reduce the duration of illness by 1.5 days if treatment is begun within 48 hours of infection [9]. However, the main drawbacks of ultilizing adamantane derivatives are that drug-resistant variants develop rapidly [10]. Presently, animal experiments and clinical studies have indicated the frequent occurrence of amantadine-resistant influenza variants after exposure to the drugs. During clinical use in humans, adamantane resistance emerged in about 30% of patients after 2–3 days of the treatment [11]. Moreover, drug-resistant variants can replicate as efficiently as sensitive ones, and they can also transmit efficiently from one individual to other individuals without apparent loss of pathogenicity [11]. Additionally, complete cross-resistance exists between amantadine and rimantadine because of the common mechanism of action [12]. To reduce the emergence and spread of drug-resistant influenza viruses, amantadine-oseltamivir combination chemotherapy has been evaluated and their synergistic antiviral effects have been confirmed [13]. Therefore, understanding the mechanism responsible for the emergence and spread of antiviral resistance is important for controlling seasonal and pandemic influenza.

After four decades of effective use of adamantanes in the prophylaxis and treatment of influenza, global resistance to these drugs has increased dramatically. In the period from 1991 to 1995, the global frequency of adamantane resistance among A (H3N2) viruses was as low as 0.8% [14]. However, in recent years, a significant worldwide increase in adamantane resistance has been reported for human H1N1 and H3N2 seasonal influenza viruses and H5N1 avian influenza viruses [1517]. For example in China, the frequency of amantadine-resistant H1N1 variants has greatly increased from 28% during 2004–2005 to 72% in 2005–2006 [18]. Similarly, in Taiwan and Canada, the frequencies of adamantane-resistance raised from 0 in 2004–2005 to 24% and 33% in 2005–2006, respectively [18]. In the United States, resistance to adamantine among A (H3N2) viruses circulating had reached 92% in 2005 [19]. Specially, adamantane resistance had reached 100% during 2005–2006 in influenza A (H3N2) isolates from some Southeast Asian countries [18]. It is worth mentioned that the most recent H5N1 strains isolated from humans and poultry in Vietnam and Thailand show resistance to the adamantanes [20]. The rapid emergence and spread of adamantane-resistant variants of influenza H1N1, H3N2, recent H7N9 viruses from China, and H5N1 viruses have led to changes in CDC recommendations for the use of adamantanes in the control of influenza A virus infections [21].

The genetic basis for resistance to adamantanes has been well characterised and amino acid substitutions at residues L26, V27, A30, S31, and G34 within the transmembrane domain of the M2 protein are associated with adamantane resistance [12, 22]. Abed et al. constructed recombinant H1N1 influenza A viruses with the commonly observed L26F, V27A, A30T, S31N, G34E, and V27A/S31N mutations in the M2 gene by using reverse genetics, and showed all of these mutations are able to cause amantadine resistance [23]. So far, the known mutations that confer resistance against adamantanes are L26F, V27A, A30T, A30V, S31N, G34E and L38F in the M2 transmembrane region [6, 22]. Bright et al. found that out of 7000 influenza A isolates collected from countries in Africa, Asia, Europe, the Americas and Oceania between 1994 and 2005, the frequency of adamantane resistance mutations rose from 0.4% in 1994–1995 to 12.3% in 2003–2004 [16]. Most drug-resistant influenza viruses contain one of these amino acid changes, but variants with dual mutations have also been described [12]. In particular, the S31N substitution renders the virus highly resistant to adamantanes, and the incidence of clinical isolates bearing this mutation has jumped from a small percentage to approximately 97% in recent years [18, 19, 24]. Given the uncertainty of the therapeutic value about adamantanes, particular to a pandemic and/or long-term endemicity, it is necessary to study deeply whether resistance to adamantanes is prevalent in various influenza viruses in the world. In this present study, we analyzed the incidence rate, host distribution, temporal distribution, and geographic distribution of adamantane-resistance among H1–H17 subtype influenza A viruses isolated worldwide from 1902 to 2013, and also tried to clarify the possible reasons of high adamantane-resistance incidence.

Methods

Sequence data

All amino acid sequence data of the M2 proteins used in this study were obtained from the influenza virus sequence database available in GenBank (Influenza Virus Resource [http://www.ncbi.nlm.nih.gov/genomes/FLU/FLU.html], accessed on 8 June and 31 October 2013). Firstly, the shorter and identical sequences were collapsed using the web servers of Influenza Virus Resource in the National Center for Biotechnology Information (NCBI), and the remaining sequences were then utilized for further analysis. Secondly, a total of 31251 M2 protein sequences (nucleotide positions 1–27 and 716–979) were selected and downloaded representing the spectrum of influenza A diversity from 1902 to 2013, together with information about the subtype, host, location and year of isolation. Subtypes H1–H17 of influenza A virus from different host species, including human, avian, swine, equine, and canine viruses, were included. As a result, M2 protein sequences of 12579 H1, 539 H2, 9414 H3, 1001 H4, 2856 H5, 1208 H6, 1147 H7, 124 H8, 1151 H9, 548 H10, 421 H11, 133 H12, 74 H13, 11 H14, 7 H15, 35 H16, and 3 H17 influenza A viruses were used for the analysis. All GenBank accession numbers of adamantane resistant viruses were listed in the S1S6 Tables through primary classification in host, subtype, and isolation location.

Sequence Analysis

All of these sequence data were aligned by Clustal W (BioEdit version 7.0.5) (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). Substitutions of residues Leu26, Val27, Ala30, Ser31, Gly34, His37, Leu38, and Trp41 in the transmembrane region of the M2 ion channel protein were used to screen amantadine-resistant mutants [12, 25]. The frequencies of the appearance of recognized adamantine resistance-associated amino acids (L26F, V27A, A30T, A30V, S31N, G34E, and L38F) were identified and tabulated. Among 31251 sequences, the number of strains with mutations was shown, and their detail information was listed in the S1S6 Tables.

Results

The frequency of adamantane resistant influenza variants

A total of 31251 influenza A viruses were screened for the specific mutations known to correlate with resistance to adamantanes. Our analysis based on the M2 protein sequence data available in the GenBank showed a diversity in the frequency of adamantane resistant influenza variants among different HA subtypes, host species, years of isolation, and geographical areas (S1S6 Tables). In the present study, a total of 14132 (45.22%) of the 31251 influenza A viruses contained one or two amino acid substitutions in the M2 protein known to cause resistance to adamantanes [12, 22]. The prevalence of adamantane resistance over the time period of the study is shown in Table 1. These results suggested that the frequency of emergence of drug-resistant strains varied among different HA subtypes. Among these isolates, we identified adamantane-resistant influenza variants from H1–H7, H9–H11, and H17 subtypes. However, analysis of viruses isolated from 1902 to 2013 did not reveal H8 and H12–H16 subtype influenza isolates with amino acid substitutions in the transmembrane region of M2 protein corresponding to resistance to adamantanes (Table 1). In other words, no adamantane resistant viruses were identified among the influenza H8, H12, H13, H14, H15, and H16 subtype. In contrast, adamantane-resistant strains frequently occurred among viruses of six HA subtypes with the frequencies of H1 (8777/12579, 69.78%), H3 (4107/9414, 43.63%), H5 (804/2856, 28.15%), H7 (144/1147, 12.55%), H9 (265/1151, 23.02%), and H17 (3/3, 100%). The remaining HA subtype influenza viruses were identified with <0.5% adamantane resistant variants except H6 subtype influenza viruses. Analysis of sequence data from H6 subtype strains indicated 2.07% (25/1208) resistant variants with M2 protein-associated mutations that would confer resistance to adamantanes [22]. Notably, the numbers of resistant variants differed among HA/NA subtype (Table 1). For example, a total of 8193 influenza A (H1N1) viruses and 583 influenza A (H1N2) isolates were found to be adamantane-resistant. Similarly, 4095 influenza A (H3N2), 7 influenza A (H3N1) and 5 influenza A (H3N8) viruses were found to contain mutations conferring resistance to adamantanes (Table 1). Our study further suggested that the frequency of emergence of adamantane-resistant strains varied among different HA or HA/NA subtypes.

thumbnail
Table 1. The frequencies of adamantane-resistant variants among all influenza A (H1-H17) viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.t001

The frequencies of adamantane resistance mutations of viruses from influenza H1–H17 subtypes are also analyzed in this study. The critical residues conferring resistance to adamantanes were shown in Table 2. Results from our analyses indicated that the viruses with L26F, V27A, A30T, G34E, and L38F mutations were found sporadically. However, the majority of drug-resistant influenza variants contained point mutations resulting in a serine to-asparagine change at amino acid 31 (S31N) of the M2 protein that conferred adamantane resistance [22]. Among these resistant viruses, 96.7% (8490/8777) of the H1, 93.6% (3844/4107) of the H3, 66.7% (2/3) of the H4, 83.3% (670/804) of the H5, 92% (23/25) of the H6, 86.1% (124/144) of the H7, 87.5% (232/265) of the H9, 100% (1/1) of the H10, 50% (1/2) of the H11, and 100% (3/3) of the H17 strains demonstrated the S31N substitution (Table 2). Furthermore, influenza H1, H3 and H5 variants contained the L26F, V27A, A30T, L26F/S31N, and V27A/S31N adamantane resistant mutations (Table 2). These resistant strains possessed amino acid substitutions in the M2 protein at three positions–L26F, V27A, and S31N [22]. The G34E and G34E/S31N mutations emerged only in influenza A H1 viruses, whereas the L38F mutation only occurred in the influenza A H2 virus. In addition, the L38F/S31N and A30V/S31N mutations were only found in influenza A H3 viruses. Our study further revealed adamantane resistance mutations appeared across different HA subtypes.

thumbnail
Table 2. The frequency of adamantane resistance mutations of viruses from different influenza HA subtypes.

https://doi.org/10.1371/journal.pone.0119115.t002

Host distribution of adamantane resistant influenza variants

Adamantane resistance is caused by a single point mutation at the following amino acid positions in the transmembrane region of the M2 protein: 26, 27, 30, 31, 34, and 38. The known mutations that confer adamantane resistance are L26F, V27A, A30T, A30V, S31N, G34E and L38F [12, 22]. To better understand host distribution of adamantane resistance, we conducted the M2 protein sequence analyses focusing on positions 26, 27, 30, 31, 34, and 38 for amino acid substitutions associated with naturally occurring adamantane resistance. The numbers of adamantane resistant influenza variants with drug resistance mutation(s) among different hosts are shown in Table 3 and Table 4. Of those, adamantane resistant influenza A H1 and H3 variants mainly emerged in human and swine. In contrast, influenza A H5, H7 and H9 subtype viruses with adamantane resistance primarily occurred in avian. Particularly, influenza A H4, H6, H10, and H11 subtype adamantane resistant strains were only detected in avian species, whereas influenza H2 and H17 subtype adamantane resistant variants were only found in human and bat, respectively (Table 3). Adamantane-resistant mutations known as L26F, V27A, A30T, S31N, G34E, and L38F were detected at all six sites (positions 26, 27, 30, 31, 34, and 38) of M2 protein in human influenza viruses (Table 4). We also found influenza strains with double adamantane resistance mutations (L26F/S31N, V27A/S31N, A30T/S31N, A30V/S31N, G34E/S31N, and L38F/S31N) at positions 26 and 31, 27 and 31, 30 and 31, 31 and 34 or 31 and 38 in the human population. Adamantane-resistant mutations (L26F, V27A, A30T, S31N) were found at positions 26, 27, 30, and 31 in avian and swine influenza viruses. Double resistance mutations at positions 26 and 31 (L26F/S31N), 27 and 31 (V27A/S31N) were also found in these hosts (Table 4). However, No G34E resistant mutation at position 34 was detected in avian influenza viruses. Similarly, Avian and swine influenza viruses also had no L38F adamantane resistant mutation at position 38. Notably, adamantane resistance mutations in all hosts were detected most frequently at position 31 (S31N), followed by position 27 (V27A). The most common drug resistance mutation was S31N in human, avian, and swine influenza viruses with the frequencies of 98.7% (10813/10954), 88.1% (909/1032), and 77.4% (1595/2062), respectively. In addition, a small number of S31N adamantane resistance mutations were also detected in canine, cheetah, ferret, mink, cat, equine, bat, and environment influenza viruses (Table 4). In S1S6 Tables we show more adamantane resistant variants carrying the S31N resistance mutation at position 31 in H3N2 and H1N1 human influenza viruses after 2003. Adamantane resistant strains of avian influenza H5N1 viruses have also been found more frequently since 2000. Our study further suggested that the frequency of emergence of adamantane-resistant strains varied among different hosts and adamantane resistance mutations appeared across different hosts.

thumbnail
Table 3. Host distribution of adamantane resistant influenza variants analyzed in this study among the different HA subtypes.

https://doi.org/10.1371/journal.pone.0119115.t003

thumbnail
Table 4. Host distribution of adamantane-resistance mutations analyzed in this study.

https://doi.org/10.1371/journal.pone.0119115.t004

Geographical distribution of adamantane-resistant influenza variants

To further understand the geographical distribution of adamantane resistance, the M2 protein sequences related to drug-resistant phenotypes were analyzed focusing solely on geographic location of influenza variants with adamantane-resistant mutations. The geographic analysis summarized in Table 5 indicated a substantially rising percentage of adamantane-resistant influenza viruses isolated from the world including specific countries in Africa, Asia, Europe, America, and Oceania. Adamantane resistant influenza A H1 variants were distributed widely in the world. Of these variants, 52.1% (4572/8777) distributed in several countries in America (United States, Canada, Mexico, and Nicaragua) and 26.4% (2315/8777) were detected in Asia (China, Singapore, Thailand, Hong Kong, South Korea, Cambodia, Viet Nam, Japan, and Taiwan). In addition, about one-fifth of the influenza H1 variants distributed in Europe (United Kingdom, Spain, Greece, Germany, Finland, Russia, and Netherlands) (Table 5). Adamantane resistant influenza A H3 variants were obtained mainly from 26 countries, showing a wide distribution across the world. Of those, 34.8% (1427/4107) isolates from Asia and 56.4% (2316/4107) variants from America were adamantane resistant. All influenza A H6 and H9 variants with adamantane resistance only distributed in Asia (China, Taiwan, and Hong Kong). For influenza A H5 adamantane resistance variants, a great number of them were distributed in Asia including China, Thailand, Cambodia, Viet Nam, and Indonesia, while a small number of them were found in Africa (Egypt). Importantly, the largest proportion of Asian drug-resistant avian influenza viruses of H5 and H9 subtypes occurred in China (Table 5). Influenza A H7 adamantane resistance variants distributed primarily in China of Asia and United States of America, with the highest frequency of resistance found in China (64.6%). Similarly, geographical analyses of adamantane-resistant mutations were also conducted in this study. The numbers of adamantane-resistant mutations from all of the HA subtypes are shown in Table 6. The L26F mutation distributed primarily in Hong Kong and New Zealand, the V27A mutation mainly distributed in China and Indonesia, and the A30T mutation distributed primarily in China. The S31N mutation distributed across the world (Table 6). This is the most common mutation known to confer resistance to adamantanes. In contrast, the G34E and L38F mutations only distributed in United States and South Korea, respectively. For double adamantane resistance mutations, the L26F/S31N mutations distributed mainly in Thailand, United Kingdom, and United States. The majority of the V27A/S31N mutations emerged in United States, Spain, South Korea, and Indonesia (Table 6). The A30T/S31N mutations were found in Nicaragua and Australia, and the G34E/S31N mutations occured in United Kingdom and Canada. The L38F/S31N and A30V/S31N mutations distributed only in Taiwan and United States, respectively. Our study further suggested adamantane resistance mutations appeared across different geographic regions and the frequency of emergence of adamantane-resistant strains varied among different geographical areas.

thumbnail
Table 5. Geographical distribution of adamantane-resistant influenza variants.

https://doi.org/10.1371/journal.pone.0119115.t005

thumbnail
Table 6. Geographical distribution of adamantane-resistant mutations analyzed.

https://doi.org/10.1371/journal.pone.0119115.t006

Temporal distribution of adamantane resistant influenza variants

To better understand time distribution of adamantane resistant influenza variants, the M2 protein sequence analyses were conducted focusing attention on years of isolation of influenza viruses with adamantane resistance. The numbers and frequencies of adamantane resistant variants from H1N1, H3N2, H5N1, H6, H7, and H9N2 subtype influenza viruses among different years of isolation are given in Table 712. The trend for overall data showed a considerable increase. And the frequency of emergence of adamantane-resistant variants from different HA subtypes varied among different years of isolation during 2001–2013. For H1 subtype, temporal distribution of adamantane-resistant influenza H1N1 variants during 1933–2013 revealed substantial increases in rates for USA, China, Singapore, United Kingdom, Spain, Canada, Thailand, Hong Kong, Mexico, Australia, Nicaragua, Finland, Greece, Russia and Taiwan (Table 7). Trends of rate changes from these countries were not considerable before the spikes occurred in 2009, which indicate the existence of such spikes and their considerable effects on the overall trends. As shown in Table 7, the highest level of adamantane resistance was continually seen during 2009–2013 at 100% in many countries (e.g. Singapore, Spain, Portugal, Nicaragua, Chile). For H3 subtype, temporal distribution of adamantane-resistant variants from H3N2 influenza viruses during 1968–2013 showed substantial increases in rates from many countries and regions across the five continents: USA, China, Hong Kong, Singapore, Taiwan, Nicaragua, Vietnam, Thailand, Uganda, Japan, Russia, Germany, Kenya, and Australia (Table 8). The frequency of adamantane-resistant H3N2 influenza viruses has increased sharply since 2003, particularly in China. The proportion of isolates resistant to adamantanes continued to increase during 2005–2006. Specifically, isolates collected during this period in Vietnam, Spain, Mexico, Cuba, and New Zealand exhibited 100% drug resistance. During 2003–2009, resistants frequency increased considerably in mainland China, Hong Kong, Taiwan and Thailand. During 2005–2013, there were substantial increases in the rates in USA and Singapore. In contrast, no adamantane resistance was observed in United Kingdom during this period. During 2008–2010, the overall rate change in Kenya, Uganda, Russia, Netherlands, Mexico, and Japan was considerable.

thumbnail
Table 7. Temporal distribution of adamantane-resistant influenza variants from H1N1 subtypes viruses during 1933–2013.

https://doi.org/10.1371/journal.pone.0119115.t007

thumbnail
Table 8. Temporal distribution of adamantane-resistant influenza variants from H3N2 subtypes viruses during 1968–2013.

https://doi.org/10.1371/journal.pone.0119115.t008

thumbnail
Table 9. Temporal distribution of adamantane-resistant influenza variants from H5N1 subtype viruses during 1959–2013.

https://doi.org/10.1371/journal.pone.0119115.t009

thumbnail
Table 10. Temporal distribution of adamantane-resistant influenza variants from H6 subtypes viruses during 1999–2013.

https://doi.org/10.1371/journal.pone.0119115.t010

thumbnail
Table 11. Temporal distribution of adamantane-resistant influenza variants from H7 subtypes viruses during 1902–2013.

https://doi.org/10.1371/journal.pone.0119115.t011

thumbnail
Table 12. Temporal distribution of adamantane-resistant influenza variants from H9N2 subtypes viruses during 1998–2013.

https://doi.org/10.1371/journal.pone.0119115.t012

For the H5 subtype, temporal distribution of adamantane-resistant influenza H5N1 variants during 1959–2013 indicated considerable rate changes from six distinctly important countries including China, Thailand, Cambodia, Vietnam, and Indonesia of Asia and Egypt of Africa. Furthermore, rates among these countries were also considerably different (Table 9). For instance, the frequency of H5N1 adamantane-resistant viruses obtained from Indonesia began increasing in 2003, with a considerable spike (100%) between 2007 and 2008. This trend of having spikes was followed by an increasing frequency in Thailand and Vietnam in 2003. Notably, the percentage of adamantane resistant H5N1 viruses gathered from Cambodia is particularly high and consistently reached 100% during 2004–2012. For H6 subtype, temporal distribution of adamantane-resistant H6 influenza variants during 1999–2013 suggested that H6N1 and H6N2 mutants with considerable rate changes emerged mainly in Taiwan and Vietnam, respectively (as shown in Table 10). No H6N2 resistant viruses were detected in samples collected from China. For H7 subtype, temporal distribution of adamantane-resistant H7 influenza variants during 1902–2013 suggested that H7N9, H7N7, and H7N6 mutants with considerable rate changes emerged mainly in China in 2013, whereas H7N3, H7N2, and H7N1 mutants were respectively detected in Pakistan, USA, and Italy by 2005 (as shown in Table 11). For H9 subtype, temporal distribution of H9N2 adamantane resistant viruses during 1998–2013 displayed considerable changes in rates in Asia for specific countries or areas such as China, Hong Kong, Bangladesh, and United Arab (Table 12). In 2009, one or more drug-resistant virus was identified in five countries: China (20/23, 87.0%), South Korea (1/24, 4.2%), Vietnam (3/4, 75.0%), Bangladesh (4/7, 57.1%), and Iran (2/2, 100%). However, no resistance was found in Hong Kong, Japan, Israel, and United Arab in the same period. By comparison, the frequency of resistant H9N2 variants isolated from China was consistently near 10% during 2002–2005.

In conclusion, our study further suggested the frequency of emergence of adamantane-resistant influenza strains varied among different years of isolation.

Discussion

Anti-influenza drugs play an important role in a comprehensive approach on controlling influenza A virus infections. Adamantanes had been excellent anti-influenza medicines until the recent emergence of resistant viruses. In this study, 31251 different subtype influenza A viruses (H1–H17) isolated in the world from 1902 to 2013 were assessed for resistance to adamantanes. On the basis of M2 protein sequence analysis, our study reveals a continuing worldwide increase in adamantane resistance, and suggests that the frequency of emergence of drug-resistant influenza variants varied among different HA subtypes, host species, years of isolation, and geographical areas. This study, which evaluates the frequency and distribution of 14132 adamantane-resistant influenza viruses obtained worldwide, is the largest and most comprehensive report on adamantane resistance to date.

The frequency of resistance to adamantanes among circulating influenza A viruses has dramatically increased over the past few years [16]. In this study we confirm that H1, H3, H5, H7, H9, and H17 subtype influenza A viruses with the specific resistance-associated mutations in their M2 genes exhibited high-level resistance to adamantanes. In contrast, the appearance of H2, H4, H6, H10, and H11 adamantane-resistant mutants was a rare event. However, no adamantane resistance viruses were identified among other HA subtypes (H8, H12–H16) (Table 1). Our findings extend those reported in 2005 [16] and 2006 [19] and document a continuous and proportional increase in influenza A (H1, H3, H5, H7, H9) variants showing resistance to adamantanes. The global incidence of adamantane resistance among A (H3N2) has increased dramatically. This escalating trend in circulating H3N2 drug-resistant variants was first observed among viruses isolated from Asia in 2000 and then in other regions of the world in 2005 (Table 8). A significant increase in resistance was also detected in many countries throughout the world among A (H1N1) viruses collected during 2009–2012 (Table 7). Although adamantane resistant H5N1 variants are present in Asia, their distribution appeared to be largely limited to Thailand, Cambodia, and Vietnam, It is worth noting that most H5N1 viruses from Indonesia and China are sensitive to adamantine [5] (Table 9). The apparent geographical disparity in the susceptibility of H5N1 isolates to adamantane is unexplained. Similarly, H9N2 viruses isolated from Asian countries including China, Hong Kong, Bangladesh, and United Arab show significant increases in drug-resistance frequencies (Table 12). Particularlly, H7N9 resistant mutants only distributed in mainland China and Taiwan of Asia also show a significant increase in the incidence (100%) of adamantane resistance in 2003 (Table 11). Our findings that epidemic and pandemic strains of influenza were identified first in Asia accord with the results from other previous studies [26].

Adamantane resistance in the influenza A virus is associated with six amino acid substitutions in the M2 protein according to previous reports [16]. The known mutations that confer adamantane resistance are L26F, V27A, A30T (A30V), S31N, G34E, and L38F [12, 22]. It was reported that mutated viruses may either lose the ability to bind M2 ion channel blockers, as with the S31N or A30T amino acid substitutions [27], or bind the blockers but retain M2 function, as with amino acid replacements L26F or V27A at residue 26 or 27 [27]. In our study described here, the most common adamantane resistance mutation was S31N in influenza A viruses from human, avian, and swine (Table 4). Most adamantane-resistant influenza variants (95%) bear S→N amino acid mutations at position 31 of the M2 protein, whereas only about 1% have V→A mutations at aa 27. Drug-resistance mutations at other amino acid positions (L26F, A30T, G34E, and L38F) are extremely rare (their prevalence appeared to be < 0.2%) (Table 2). Among these mutations, the S31N is the most frequently reported [14, 19, 28], suggesting that variants containing the S31N substitution might possess a significant advantage on viral replication or transmission, leading to more efficient circulation. Additionally, the high levels of the dual L26I and S31N resistance mutations mainly detected in H5N1 variants obtained from Thailand, Vietnam, and Cambodia indicate that viruses carrying this dual motif are stably selected (S3 Table). Given that the dual L26I and S31N motif are identified firstly in China in 2002 in our study, it appears that China might have been the location of introduction or generation of viruses. Our results further suggest that M2 drug-resistant mutations could have occurred spontaneously before these drugs were developed, implying that greater caution is needed in the use of adamantanes.

It is known that influenza viruses resistant to amantadine and rimantadine can emerge quickly when these drugs are being used to control influenza outbreaks [11, 29]. Adamantanes have been considered first-line drugs for the prophylaxis and treatment of influenza A virus infections. The widespread use of adamantanes has been associated with the rapid emergence of resistant viruses which are as genetically stable, virulent and transmissible as the wild-type virus [30]. Here we report the results of a comprehensive study focusing on adamantane resistance in influenza A viruses circulating worldwide. Our results reveal a significant increase in adamantane-resistance frequencies in influenza A H1N1, H3N2, H5N1, H6, H7, and H9N2 viruses circulating during 2001–2013 (Table 712). The high levels of adamantane-resistance found in our study are comparable with other studies monitoring resistance worldwide [17, 18, 19, 31]. For example, adamantane resistance in A (H3N2) strains continued to increase in many countries from 2005 to 2007 (e.g. New Zealand, USA, Malaysia, and China) [31] (Table 8). However, resistance in A (H1N1) strains had been quite variable during that time period in some countries such as Japan, Australia, USA and Canada [15, 32, 33] (Table 7). The high frequency of adamantane-resistant variants indicates that continuous global surveillance and rapid identification of mutants are essential to monitor the emergence and spread of drug resistance, and to help with making informed decisions about antiviral usage in control of influenza virus infections.

To our knowledge, this study is the first attempt to evaluate the frequency and distribution of adamantane resistant influenza variants throughout the world over a period of 111 years. Results from our study provide further evidence that amino acid residues 26, 27, 30, 31, 34, and 38 in the M2 protein play a major role in determining the resistance phenotype of adamantanes of influenza A viruses and highlight the necessity of monitoring the susceptibility of influenza A viruses to antiviral drugs. Our findings raise concerns about the increasing prevalence of adamantane-resistant influenza variants and draw attention to the importance of tracking the emergence and worldwide spread of drug-resistant variants.

Supporting Information

S1 Table. Accession numbers and adamantane resistance mutations of M2 protein sequences from H1 subtype influenza A viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.s001

(XLSX)

S2 Table. Accession numbers and adamantane resistance mutations of M2 protein sequences from H3 subtype influenza A viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.s002

(XLSX)

S3 Table. Accession numbers and adamantane resistance mutations of M2 protein sequences from H5 subtype influenza A viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.s003

(XLSX)

S4 Table. Accession numbers and adamantane resistance mutations of M2 protein sequences from H7 subtype influenza A viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.s004

(XLSX)

S5 Table. Accession numbers and adamantane resistance mutations of M2 protein sequences from H9 subtype influenza A viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.s005

(XLSX)

S6 Table. Accession numbers and adamantane resistance mutations of M2 protein sequences from H2, H4, H6, H10, H11, and H17 subtype influenza A viruses used in this study.

https://doi.org/10.1371/journal.pone.0119115.s006

(XLSX)

Acknowledgments

We appreciate Dr. Hong Zhang from Z-Biomed for helping with the language expression of the manuscript. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author Contributions

Conceived and designed the experiments: HH GJ GD. Performed the experiments: GD CP. Analyzed the data: JL CW LH BW. Wrote the paper: GD.

References

  1. 1. Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y. Evolution and ecology of influenza A viruses. Microbiol Rev. 1992; 56: 152–179. pmid:1579108
  2. 2. Webster RG, Hulse D. Controlling avian flu at the source. Nature 2005; 35: 415–416.
  3. 3. Gerdil C. The annual production cycle for influenza vaccine. Vaccine. 2003; 21: 1776–1779. pmid:12686093
  4. 4. Englund JA. Antiviral therapy of influenza. Semin Pediatr Infect Dis. 2002; 13: 120–128. pmid:12122950
  5. 5. Hurt AC, Selleck P, Komadina N, Shaw R, Brown L. Susceptibility of highly pathogenic A (H5N1) avian influenza viruses to the neuraminidase inhibitors and adamantanes. Antiviral Res. 2007; 73: 228–231. pmid:17112602
  6. 6. Cheung CL, Rayner JM, Smith GJ, Wang P, Naipospos TS, et al. Distribution of mantadine-resistant H5N1 avian influenza variants in Asia. J Infect Dis. 2006; 193: 1626–1629. pmid:16703504
  7. 7. Dolin R, Reichman RC, Madore HP, Maynard R, Linton PN, et al. Controlled trial of amantadine and rimantadine in the prophylaxis of influenza A infection. N Engl J Med. 1982; 307: 580–584. pmid:7050702
  8. 8. Belshe RB, Burk B, Newman F, Cerruti RL, Sim IS. Resistance of influenza A virus to amantadine and rimantadine: results of one decade of surveillance. J Infect Dis. 1989; 159: 430–435. pmid:2915166
  9. 9. Reuman PD, Bernstein DI, Keefer MC, Young EC, Sherwood JR, et al. Efficacy and safety of low dosage amantadine hydrochloride as prophylaxis for influenza A. Antiviral Res. 1989; 11: 27–40. pmid:2712549
  10. 10. Hayden FG. Amantadine and rimantadine—clinical aspects, In Richman D. D. (ed.), Antiviral drug resistance: John Wiley and Sons; 1996. pp. 59–77.
  11. 11. Hayden FG, Hay AJ. Emergence and transmission of influenza A viruses resistant to amantadine and rimantadine. Curr Top Microbiol Immunol. 1992; 176: 119–130. pmid:1600749
  12. 12. Belshe RB, Smith MH, Hall CB, Betts R, Hay AJ. Genetic basis of resistance to rimantadine emerging during treatment of influenza virus infection. J Virol. 1988; 62: 1508–1512. pmid:3282079
  13. 13. Ilyushina NA, Hoffmann E, Salomon R, Webster RG, Govorkova EA. Amantadine oseltamivir combination therapy for H5N1 influenza virus infection in mice. Antivir Ther. 2007; 12: 363–370. pmid:17591026
  14. 14. Ziegler T, Hemphill ML, Ziegler ML, Perez-Oronoz G, Klimov AL, et al. Low incidence of rimantadine resistance in field isolates of influenza A viruses. J Infect Dis. 1999; 180: 935–939. pmid:10479115
  15. 15. Saito R, Suzuki Y, Li D, Zaraket H, Sato I, et al. Increased incidence of adamantane-resistant influenza A(H1N1) and A(H3N2) viruses during the 2006–2007 influenza season in Japan. J Infect Dis. 2008; 197: 630–632. (Author’s reply, 197:632–633.) pmid:18275281
  16. 16. Bright R A, Medina MJ, Xu X, Perez-Oronoz G, Wallis TR, et al. Incidence of adamantine resistance among influenza A (H3N2) viruses isolated worldwide from 1994 to 2005: a cause for concern. Lancet. 2005; 366: 1175–1181. pmid:16198766
  17. 17. He G, Qiao J, Dong C, He C, Zhao L, et al. Amantadine resistance among H5N1 avian influenza viruses isolated in Northern China. Antivir Res. 2008; 77: 72–76. pmid:17897729
  18. 18. Deyde VM, Xu X, Bright RA, Shaw M, Smith CB, et al. Surveillance of resistance to adamantanes among influenza A(H3N2) and A(H1N1) viruses isolated worldwide, J Infect Dis. 2007; 196: 249–257. pmid:17570112
  19. 19. Bright RA, Shay DK, Shu B, Cox NJ, Klimov AI. Adamantane resistance among influenza A viruses isolated early during the 2005–2006 influenza season in the United States. JAMA. 2006; 295:891–894. pmid:16456087
  20. 20. World Health Organization Global Influenza Program Surveillance Network. Evolution of H5N1 avian influenza viruses in Asia. Emerg Infect Dis. 2005; 11: 1515–1521. pmid:16318689
  21. 21. Bright R A, Shay D, Bresee J, Klimov A, Cox N, et al. High levels of adamantane resistance among influenza A (H3N2) viruses and interim guidelines for use of antiviral agents—United States, 2005–2006 influenza season. MMWR Morb Mortal Wkly Rep. 2006; 55: 44–46. pmid:16424859
  22. 22. Hay AJ, Zambon MC, Wolstenholme AJ, Skehel JJ, Smith MH. Molecular basis of resistance of influenza A viruses to amantadine. J Antimicrob Chemother. 1986; 8 (suppl B): 19–29.
  23. 23. Abed Y, Goyette N, Boivin G. Generation and characterization of recombinant influenza A (H1N1) viruses harboring amantadine resistance mutations. Antimicrob Agents Chemother. 2005; 49: 556–559. pmid:15673732
  24. 24. Hayden FG. Antiviral resistance in influenza viruses-implications for management and pandemic response. N Engl J Med. 2006; 354: 785–788. pmid:16495389
  25. 25. Wu B, Wang CM, Dong GY, Luo J, Zhao BH, et al. Molecular characterization of H1N1 influenza A viruses from human cases in North America. Chinese Sci Bull. 2009; 54: 2179–2192.
  26. 26. Cox NJ, Brammer TL, Regnery HL. Influenza: global surveillance for epidemic and pandemic variants. Eur J Epidemiol. 1994; 10: 467–470. pmid:7843358
  27. 27. Astrahan P, Kass I, Cooper MA, Arkin I. A novel method of resistance for influenza against a channel-blocking antiviral drug. Proteins. 2004; 55: 251–257. pmid:15048819
  28. 28. Tang JW, Ngai KL, Wong JC, Lam WY, Chan PK. Emergence of adamantane-resistant influenza A(H3N2) viruses in Hong Kong between 1997 and 2006, J Med Virol. 2008; 80: 895–901. pmid:18360910
  29. 29. Shiraishi K, Mitamura K, Sakai-Tagawa Y, Goto H, Sugaya N, et al. High frequency of resistant viruses harboring different mutations in amantadine-treated children with influenza. J Infect Dis. 2003; 188: 57–61. pmid:12825171
  30. 30. Bean WJ, Threlkeld SC, Webster RG. Biologic potential of amantadine-resistant influenza A virus in an avian model. J Infect Dis. 1989; 159: 1050–1056. pmid:2723453
  31. 31. Barr IG, Hurt AC, Iannello P, Tomasov C, Deed N, et al. Increased adamantane resistance in influenza A (H3) viruses in Australia and neighbouring countries in 2005. Antiviral Res. 2007; 73: 112–117. pmid:16963130
  32. 32. Barr IG, Hurt AC, Deed N, Iannello P, Tomasov C, et al. The emergence of adamantane resistance in influenza A(H1) viruses in Australia and regionally in 2006. Antiviral Res. 2007; 75: 173–176. pmid:17350111
  33. 33. Deyde VM, Garten R, Sheu TG, Gubareva LV, Klimov AI. Reply to Saito et al. J Infect Dis. 2008; 197: 632–633.