Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2009, p. 1372-1378, Vol. 47, No. 5
0095-1137/09/$08.00+0 doi:10.1128/JCM.01993-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.
,
Wadsworth Center, New York State Department of Health, Albany, New York,1 Wisconsin State Laboratory of Hygiene, 465 Henry Mall, Madison, Wisconsin 537062
Received 15 October 2008/ Returned for modification 12 December 2008/ Accepted 17 March 2009
|
|
|---|
|
|
|---|
Adamantanes are effective only against influenza type A virus; they function by blocking ion flow through the 97-amino-acid viral proton channel (the M2 protein of the matrix gene), which is necessary for replication (1). Adamantanes have been shown to be effective agents for the treatment of influenza symptoms as well as for prophylaxis (13, 23). However, recent publications have documented a rapid increase in the incidence of adamantane-resistant influenza virus strains during the 2005-2006 influenza season (4, 12). Adamantane resistance is associated with mutations in a short region of the M2 protein sequence encompassed in nucleotides 789 to 815 of the matrix gene. Thus far, five resistance mutations, in the codons for amino acids 26, 27, 30, 31, and 34, have been reported. The predominant circulating resistance mutation reported affects amino acid 31, resulting in a serine-to-asparagine (S31N) change (24, 26). During the 2005-2006 influenza season, the Centers for Disease Control and Prevention (CDC) found that more than 90% of A/H3N2 viruses circulating in the United States had the S31N adamantane resistance mutation and recommended that adamantanes not be used for treatment or prophylaxis of influenza until susceptibility can be reestablished (5). The results also indicated that only 4% of A/H1N1 viruses tested were adamantane resistant (12). In Asia during the same period, the levels of resistance were slightly higher: nearly all A/H3N2 and 15% of A/H1N1 viruses tested had an adamantane resistance mutation (12).
The use of NAIs has increased as a result of the decreased efficacy of adamantanes against such a large proportion of circulating influenza viruses. NAI antivirals are analogues of sialic acid, the biological neuraminidase (NA) substrate. The NAIs disrupt viral replication by occupying the NA active site, thus limiting the binding of sialic acid and the subsequent NA cleavage of host cell receptors (21, 28). Although currently not as common as adamantane resistance, NAI resistance can develop as a result of amino acid-altering mutations in the NA gene of influenza virus. Characterized NAI resistance mutations, distributed throughout the NA gene, result in catalytic or structural changes that affect the sialic acid binding site. Several studies have developed models predicting that the increased use of NAI antivirals, such as during an influenza pandemic, may result in an increase in circulating NAI-resistant viruses, which could greatly decrease the overall efficacy of antivirals used in a pandemic situation (10, 20). In Japan, oseltamivir is often prescribed for prophylactic therapy, as well as for treatment of influenza cases. The widespread use of the drug has been cited as a reason for higher rates of resistance in Japan, previously found to be as high as 18% in pediatric cases (18). In the early part of the 2007-2008 influenza season, an oseltamivir-resistant A/H1N1 virus strain was reported first in Norway and subsequently in other countries worldwide (19, 29). In many European nations where oseltamivir use is very limited, a high percentage of the 2007-2008 influenza season's A/H1N1 virus strains had the oseltamivir resistance mutation, H274Y. This finding may indicate that oseltamivir resistance can develop and spread, even in the absence of the selective pressure exerted by excessive use of the drug. Despite variable reports on the transmissibility of drug-resistant influenza viruses (3, 16), the H274Y mutant spread rapidly throughout Europe during the 2007-2008 season. Thus far, the oseltamivir-resistant virus has remained susceptible to zanamivir (7, 22).
Clearly, timely and reliable methods are needed to monitor the development of resistance to influenza antivirals if these drugs are to be used effectively in both seasonal and pandemic situations. Here, we describe such methods for monitoring the mutations that are known to increase resistance to either adamantanes or NAIs. Pyrosequencing allows the rapid screening of the critical region of the influenza A virus matrix gene, for detection of specific adamantane resistance mutations, or of the NA gene, for detection of a specific NA resistance mutation (i.e., H274Y). For more-extensive NAI resistance genotypic analysis, a method that uses traditional Sanger sequencing and that is based on newly compiled databases of influenza virus sequences is described. We also report results, acquired with these methods, from surveillance programs conducted by two state public health laboratories during the 2006-2007 and 2007-2008 influenza seasons in New York (NY) and Wisconsin (WI).
|
|
|---|
(ii) RT-PCR. A Qiagen one-step RT-PCR kit was used throughout. Results were analyzed using Bio-Rad iQ5 software 2.0 (Bio-Rad, Hercules, CA). The crossing threshold (Ct) for each run was set by the user to the middle of the exponential portion of the amplification curve.
For conventional sequencing, RT-PCR was performed and samples were analyzed on 1% agarose gels. Bands of appropriate size were purified with Ultrafree-DA spin columns (Millipore, Bedford, MA) or with a Zymoclean gel DNA recovery kit (Zymo Research, Orange, CA) before sequencing.
(iii) Data set construction. Four different data sets were created for the optimization and selection of PCR primers, as well as for selection of additional sequences used to create positive control oligonucleotides. For targets within the M2 region of the influenza A virus matrix gene, full-length nucleotide gene sequences were selected from the Influenza Virus Resource of NCBI (2); this resource includes data from the Influenza Genome Sequencing Project of the National Institute of Allergy and Infectious Diseases and from GenBank. Where groups of identical sequences were present, representative sequences were selected. The final matrix gene data set included nucleotide sequences from 168 viruses collected worldwide over 72 years. Sequences were aligned, analyzed for conserved regions (>20 nucleotides homologous throughout the data set), and checked for errors using MEGA (Molecular Evolutionary Genetics Analysis) software version 3.1 (Center for Evolutionary Functional Genomics, Tempe, AZ) and Sequencher 4.5 (Gene Codes, Ann Arbor, MI) software. Pyrosequencing primer M2seqR and conventional sequencing primer DDseqM2_F (Table 1) were selected. This data set was also used to verify primer sequences that we selected for use from those described by Bright et al. (4) using the default settings of the CLUSTALW program within the BioEdit sequence alignment editor (Ibis Biosciences, Carlsbad, CA) and Sequencher 4.5.
|
View this table: [in a new window] |
TABLE 1. Influenza virus matrix gene M2 primers
|
|
View this table: [in a new window] |
TABLE 2. Influenza virus NA gene RT-PCR and dideoxy sequencing primers
|
(iv) Pyrosequencing. Pyrosequencing of samples from NY was done by following primarily the procedure of Bright et al. (4), with modifications to allow sequencing in both directions. For pyrosequencing of the forward strand, a second set of reactions, using a biotinylated forward primer (M2_F Biotin), a nonbiotinylated reverse primer (M2_R), and a new pyrosequencing primer, M2seqR (Table 1), were used to produce sequence complementary to that in the original Bright et al. reaction.
All NY M2 pyrosequencing was performed by the Wadsworth Center Molecular Genetics Core Facility, using a PyroMark ID pyrosequencer (Biotage AB, Uppsala, Sweden) according to the manufacturer's protocol in a cyclic dispensation mode (GATC) for 20 cycles. Sequences were analyzed with a 36-bp quality window. The software automatically designated each sample as "pass," "check," or "fail." "Pass" samples were unambiguous, with clearly defined peaks, while "fail" runs were clearly not usable, with no clear peaks above background or only nonsense sequences. "Check" runs required examination by the operator and were usually reliable sequences with some complicating factors, such as spurious peaks, added or deleted bases in homopolymer regions, and wide or double peaks, that could often be corrected by the operator.
Nucleic acid extracted from a cultured sample from a patient infected with H3N2 virus was used as an adamantane-sensitive control for pyrosequencing. Synthetic adamantane-resistant positive controls were designed by alteration of the A/Hong Kong/1/68 virus sequence so that it contained the following five amino acid codons known to confer adamantane resistance: L26F, V27A, A30T, S31N, and G34E (1, 25). Control sequences designated MUTA/HKcontF and MUTA/HKcontR (see Fig. S1 in the supplemental material) were designed to hybridize with pyrosequencing primers M2seqF and M2seqR, respectively, thus allowing sequencing in both directions. 5'-Biotinylated oligonucleotides of these control sequences were synthesized (Integrated DNA Technologies) with an additional 18 nucleotides of marker sequence. Control oligonucleotides (in binding buffer) were used at 0.625 µM in the pyrosequencing protocol.
(v) Dideoxy sequencing. Dideoxy sequencing was performed by the Wadsworth Center Molecular Genetics Core Facility using an ABI BigDye Terminator version 3.1 cycle sequencing kit with an ABI 3130 or 3730 DNA analyzer (Applied Biosystems, Foster City, CA).
For sequencing of the M2 region of the influenza A virus matrix gene, forward primer DDseqM2_F was selected from the matrix gene data set for use with primer M2_R (Table 1). RT-PCR with this primer pair produced a 456-bp amplicon encompassing nucleotides 572 to 1027 of the gene; this amplicon was subsequently sequenced in both directions by using the same primers as for the amplification reaction.
For dideoxy NA gene sequencing, RT-PCR primer pairs for A/H1N1, A/H3N2, and influenza B viruses were chosen (from their respective data sets, as detailed above) to amplify approximately 1.3 kb of the NA gene (Table 2). The amplified segments encompassed approximately 92% of the gene and contained all reported sites of mutation previously shown to decrease susceptibility to NAIs, as well as regions considered to be potential sites of additional, as-yet-uncharacterized resistance mutations. Internal primers are shown in Table 2, along with the amplification primers. Sequences obtained from these reactions were then aligned with Sequencher software, and the complete consensus gene sequence was determined.
WI. (i) Specimens and initial screening. Influenza virus was obtained from patient respiratory samples and/or isolates from cultured patient samples. Total nucleic acid was extracted using a MagNA Pure LC apparatus (Roche Diagnostics, Basel, Switzerland) according to the manufacturer's instructions. All primary samples were typed and HA subtyped using real-time RT-PCR. The NA subtype was inferred from the HA subtype through confirmation of a subset of samples, via sequencing at the CDC, to determine currently circulating influenza A subtypes. The NA subtype was further confirmed by successful human A/H1N1 virus-specific pyrosequencing of the NA gene.
(ii) RT-PCR. A Qiagen one-step RT-PCR kit (Qiagen, Inc., Germantown, MD) was used throughout.
(iii) Pyrosequencing. Pyrosequencing of the WI samples was done by using the same methods as for the NY samples (described above) for unidirectional pyrosequencing of both the M2 and NA genes of the A/H1N1 virus. For each reaction, the biotinylated RT-PCR product size was confirmed by using precast 2% agarose E gels (Invitrogen, Carlsbad, CA), by following the manufacturer's protocol. The remaining sample was used for pyrosequencing reactions. Adamantane- and NAI-sensitive and -resistant controls were selected and validated from nucleic acid extracted from cultured surveillance clinical specimens referred from participating clinical virology laboratories. Primer sequences for pyrosequencing detection of the H274Y mutation in the NA gene were obtained from Deyde et al. (11).
|
|
|---|
Further studies were conducted to determine the real-time PCR Ct value, as an indicator of the viral RNA quantity sufficient to yield quality pyrosequencing results directly from primary specimens. In preliminary studies in WI, a Ct value of 26.0 was chosen as the cutoff for exclusion from pyrosequencing testing for the S31N adamantane resistance mutation (see Fig. S2A in the supplemental material). Four of 33 test samples had Ct values greater than 26 and would not have been pyrosequenced had the cutoff been used. Of the four, three had "fail" pyrosequencing results and one had a "check" result. A subsequent analysis of all 2006-2007 samples revealed that the excluded specimens (those with Ct values greater than 26) accounted for approximately 15% of all surveillance samples collected. Similar pyrosequencing validation studies in WI were conducted for detection of the H274Y NAI resistance mutation (see Fig. S2B in the supplemental material), but a cutoff value was not used, and all samples were pyrosequenced. Of 59 samples, 5 were rated "fail," 8 were rated "check," and 46 were rated "pass." Only the pyrosequencing results that were rated "pass" gave usable results.
Analysis of the bidirectional adamantane resistance mutation pyrosequencing quality ratings for the 2006-2007 NY samples revealed that nearly all primary samples tested generated usable sequence, when tested directly. For NY samples, a cutoff Ct value was not used; all 91 primary specimens were tested. In the forward direction (see Fig. S2C in the supplemental material), all "pass" and "check" samples (82 and 8, respectively) gave usable sequences. Only two primary samples were not pyrosequenced successfully; however, for one of these "fail" samples, sequence was obtained using the reverse primer. Thus, only one primary specimen from the entire season did not give satisfactory results in at least one direction.
The NY validation also included a comparison of pyrosequencing and dideoxy sequencing results from 24 A/H1N1 and 35 A/H3N2 virus samples obtained during 2001 to 2007. The 24 A/H1N1 virus samples consisted of nucleic acid extracts from 16 primary specimens (7 with matched culture isolates) and 8 culture isolates only, all from different patients. No known adamantane resistance mutations were found in H1N1 virus samples obtained from primary extracts (Table 3). One H1N1 virus from a culture extract obtained in September 2006 (no primary specimen extract was available) had a S31N mutation; however, the patient had recently traveled to China. The 35 A/H3N2 virus samples consisted of 22 primary specimens (two with matched culture isolates) and 13 culture isolates only, all from different patients. In 30 (86%) H3N2 virus samples, a S31N resistance mutation was detected (Table 3); of these, 18 were primary specimens and 12 were culture extracts for which no primary sample was available. No other known resistance mutations were detected in any of the samples. Dideoxy sequencing results were identical to pyrosequencing results in all cases (results not shown).
|
View this table: [in a new window] |
TABLE 3. NY influenza A virus adamantane resistance mutations, 2001 to 2007
|
Adamantane resistance surveillance. A total of 236 surveillance samples (primary specimens or culture isolates) identified as influenza A virus were successfully evaluated during the 2006-2007 influenza season for the presence of amino acid changes conferring resistance to the adamantanes, at the NY (90 samples) and WI (146 samples) reference laboratories (Table 4). Of these, 144 (61.0%) were subtyped as influenza A/H1 virus, 81 (34.3%) as influenza A/H3 virus, and 11 (4.7%) could not be subtyped. None of the A/H1 virus samples possessed resistance mutations, while 82.7% of the A/H3 virus samples and 54.5% of the nontyped samples had the S31N mutation.
|
View this table: [in a new window] |
TABLE 4. Adamantane S31N resistance mutations detected in NY and WI specimens
|
NAI resistance surveillance. A total of 54 influenza virus surveillance samples were successfully evaluated in NY during the 2006-2007 influenza season for the presence of mutations conferring resistance to the NAIs (Table 5). Of these, 18 (33.3%) were subtyped as influenza A/H1N1 virus, 17 (31.5%) as influenza A/H3N2 virus, and 19 (35.2%) as influenza B virus. None of these influenza virus samples possessed a NAI resistance mutation.
|
View this table: [in a new window] |
TABLE 5. Neuraminidase H274Y resistance mutations detected in NY and WI specimens
|
|
|
|---|
Pyrosequencing was used to screen specimens for adamantane-resistant influenza A viruses. Pyrosequencing is a relatively new method suitable for rapid acquisition of short reads (100 to 200 bp) of genomic sequence, so as to identify the presence of known mutations. Because the adamantane resistance mutations are clustered in a short span of the influenza M2 gene sequence (see Fig. S1 in the supplemental material), pyrosequencing is an appropriate tool for rapid analysis. Using newly selected data sets, the NY laboratory used a modification of a previously published pyrosequencing method for this genomic region, enabling sequencing of the influenza virus in both directions, as a further confirmatory measure. Additionally, with the design of the new M2 gene dideoxy sequencing assay, we were able to directly compare the two methods. The consistency in results between the pyrosequencing and dideoxy sequencing assays served to validate the pyrosequencing methodology. Occasionally, poor pyrosequencing results were identified by the PyroMark ID software or during postrun analysis. These included spurious peaks, added or deleted bases in homopolymer regions, and wide or double peaks; such errors tend to arise from an overloading of starting material in the pyrosequencing reaction. Reduction of the biotinylated PCR product starting amount to 15 µl usually corrected the problem.
Results also confirmed that pyrosequencing of extracts from fresh (unfrozen) primary specimens during routine influenza surveillance can be a reliable method for the rapid detection of adamantane and NAI resistance in influenza A viruses. When a primary specimen test results in a pyrosequencing quality rating of "fail," the test should be performed on the corresponding virus culture isolate if available. Overall, we could be reasonably confident that samples with Ct values in the high 20s, or below, in the initial influenza A real-time RT-PCR screening assay would be accurately sequenced with the pyrosequencing method. The cutoff value of 26 was chosen to ensure that a high percentage of the samples tested with pyrosequencing would yield usable sequence. In a public health setting, testing the majority of the samples will often provide sufficient information for assessment of an outbreak, or evaluation of resistance in a state or region. The cutoff value was included as a possible guideline for obtaining information from the majority of the samples in a short time with minimal unusable results. While bidirectional pyrosequencing is not a necessity in a surveillance system, it does give increased confidence in results, especially if they are to be reported clinically. Also, for samples for which no sequence can be obtained in one direction, the other direction may give acceptable data. In fact, the reverse pyrosequencing method was slightly more reliable in this regard than the originally published forward method (data not shown).
Pyrosequencing was also used to screen specimens for oseltamivir-resistant A/H1N1 viruses. Because the H274Y mutation had previously been detected and was well characterized (9, 17), primers for pyrosequencing of this region (11) were used for rapid screening to detect this mutation. To monitor for other NAI resistance mutations that could be present in currently circulating influenza A/H1N1, A/H3N2, and B virus strains, primers were developed in NY to allow dideoxy sequencing of the NA gene. This method can be more informative than pyrosequencing, since it generates bidirectional sequence from nearly the entire gene, using two sequencing reactions per direction. Unfortunately, sample extracts in some cases had been frozen and thawed up to three times, promoting RNA degradation; this is a factor that can particularly affect influenza B virus, which is more sensitive to repeated freezing and thawing (J.M.L., unpublished observations). Six point mutations conferring amino acid changes, in addition to the H274Y mutation, were detected with this method. Initial testing found that none inhibited NA (data not shown), and further work on these mutations is ongoing. However, given the current specific interest in the circulating influenza A virus oseltamivir-resistant H274Y viruses, any combination of the A/N1 forward and reverse primers in Table 2 can be used in a single assay to detect the H274Y mutation, which is situated in an area common to all primer set amplicons.
In the early part of the 2007-2008 influenza season, routine monitoring of antiviral resistance resulted in the identification of a number of oseltamivir-resistant A/H1N1 viruses. Continued screening during the season confirmed that all oseltamivir-resistant viruses had an NA gene mutation that results in the H274Y amino acid change, and all were subtype A/H1N1. For the 2007-2008 season, the CDC reported that 10.2% of A/H1N1 virus surveillance strains contained the oseltamivir resistance mutation. For NY and WI, the respective rates were 11.1 and 17.4% of A/H1N1 virus-positive specimens. The rapid emergence of oseltamivir resistance in this subtype, from less than 1% to more than 10%, suggests that an increasingly large proportion of A/H1N1 viruses will not be treatable with oseltamivir and that rapid and reliable methods for determining the antiviral resistance status of influenza viruses are crucial. Initial data for the 2008-2009 influenza season indicate that nearly all A/H1 viruses have the H274Y mutation and are oseltamivir resistant (8).
In addition to the sequencing methods presented here for the detection of resistance mutations, real-time PCR methods sensitive enough to allow detection of the separate, single-base-pair differences that confer resistance to influenza antivirals are needed; several new methods have recently been published (6, 14). Real-time PCR methods capable of this type of detection will allow increasing numbers of laboratories to monitor resistance, since the instrumentation necessary is already available in a large number of clinical laboratories. We plan to expand our two-state surveillance network for influenza antiviral resistance detection to include other state public health laboratories that perform resistance testing on influenza-positive specimens and isolates. Such a network is needed to provide region-specific antiviral test results in support of CDC testing efforts and will be important in any potential pandemic situation.
This work was supported in part by an appointment to the Emerging Infectious Disease Fellowship Program administered by the Association of Public Health Laboratories and was funded in part by the CDC.
Published ahead of print on 25 March 2009. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»