key: cord-0854962-lhe51ws7 authors: Wong, Sallene; Pabbaraju, Kanti; Wong, Anita; Fonseca, Kevin; Drews, Steven J. title: Development of a real-time RT-PCR assay for detection of resistance to oseltamivir in influenza A pandemic (H1N1) 2009 virus using single nucleotide polymorphism probes date: 2011-05-31 journal: Journal of Virological Methods DOI: 10.1016/j.jviromet.2011.02.014 sha: aa8592ed325e1e4cfa48e85e7557c3e9dea70a39 doc_id: 854962 cord_uid: lhe51ws7 Abstract Resistance to oseltamivir in pandemic (H1N1) 2009 influenza A virus is linked to an amino acid change from histidine (H) to tyrosine (Y) at position 275 in the neuraminidase protein (NA). A real-time one step RT-PCR assay using single nucleotide polymorphism (SNP) probes was developed to detect this mutation in respiratory specimens. The limit of detection was 47.6 copies/reaction for wild-type H275 RNA and 52.9 copies/reaction for the mutant H275Y RNA. The assay did not cross-react with other respiratory pathogens. The clinical sensitivity and specificity of the assay was compared to the gold standard Sanger sequencing method using 25 sensitive, 15 resistant and 20 negative samples. The sensitivity and specificity was 88.0% and 100% respectively with the SOIV_Osel_SEN probe designed to detect the H275 allele and 100% for the SOIV_Osel_RES probe detecting the 275Y allele. The sensitivity of the assay using nine admixtures of sensitive and resistant alleles was 88.9% and 77.8% with the SOIV_Osel_SEN probe and SOIV_Osel_RES probe respectively. The presence of mixed sensitive and resistant alleles in patient samples and mixtures of in vitro RNA were detected reproducibly. This assay can be used for screening of original samples for oseltamivir resistance without the need for culture and phenotypic testing. The influenza pandemic of 2009 saw the emergence of a new influenza A (H1N1) 2009 strain that showed initial resistance to adamantanes and susceptibility to oseltamivir (Dawood et al., 2009) . However, by the summer of 2009, there were reports of isolates with an amino acid change from histidine (H) to tyrosine (Y) at position 275 in the neuraminidase (NA) glycoprotein that was linked to oseltamivir resistance (Chen et al., 2009; Harvala et al., 2010; Janies et al., 2010; Thabet et al., 2010; WHO, 2009a) . Some of these cases were identified in Alberta, Canada starting in July 2009. Although the prevalence of oseltamivir resistance in the 2009-2010 respiratory season has been low, there have been instances in the recent past of dramatic increases in oseltamivir and adamantane resistance in seasonal influenza A (IFVA) over a relatively short time frame in the absence of selective pressure by anti-viral usage (Deyde et al., 2007; Dharan et al., 2009; Hauge et al., 2009; Hurt et al., 2009; Janies et al., 2010; Nelson et al., 2009; Pabbaraju et al., 2010; Sheu et al., 2008) . Resistance to oseltamivir has been reported to result from point mutations in several regions of the NA protein of the virus (e.g., D79G, S247G or S247N, and H275Y) (Sheu et al., 2008) . Phylogenetic studies have revealed that the oseltamivir-resistant seasonal IFVA (H1N1) strain evolved from a susceptible strain (clade 2B) which circulated in the 2007-2008 season by acquiring the H275Y change in the NA gene . Currently, the gold standard method for detection of the C823T nucleotide polymorphism that leads to the H275Y amino acid change is dideoxy Sanger sequencing of the NA gene. Neuraminidase inhibition assays including commercial kits such as the NA-Star ® Influenza Neuraminidase Inhibitor Resistance assay (Applied Biosystems, Foster city, CA, USA) have been described for phenotypic resistance testing (Hurt et al., 2009 ). These phenotypic assays require culturing of the virus, which would not be possible in the case of novel sub-type variants such as the pandemic (H1N1) 2009 prior to characterization of the virus. Both sequencing and phenotypic testing are labor intensive, time-consuming, and expensive. The use of pyrosequencing for monitoring nucleotide changes has also be described, however this requires specialized instrumentation (Deyde et al., 2009 (Deyde et al., , 2010 . A reverse-transcriptase PCR coupled with a restriction fragment length polymorphism assay (RTPCR-RFLP) has been described for seasonal H1N1 viruses but it is a lengthy process and not practical for high volume testing (Guo et al., 2009 determination of oseltamivir resistance for seasonal and pandemic H1N1 influenza A viruses has been reported recently using liquid microarrays on the Luminex platform (Mahony et al., 2010) . Realtime reverse-transcriptase (RT)-PCR assays using single nucleotide polymorphism (SNP) probes for allelic discrimination have been described for the detection of H275Y mutation in seasonal H1N1 viruses (Bolotin et al., 2009; Carr et al., 2008; Operario et al., 2010) and pandemic (H1N1) 2009 viruses (van der Vries et al., 2010) . Rapid methods for the detection of mutations leading to oseltamivir resistance allow diagnostic laboratories to monitor the development and transmission of resistance especially in high risk populations such as immuno-compromised patients and patients requiring antiviral treatment. This study describes the validation of a high throughput, sensitive and specific RT-PCR for detection of the C823T nucleotide polymorphism leading to the H275Y amino acid change in the NA protein of pandemic (H1N1) 2009 virus. This assay could be used in surveillance of tamiflu resistance, detection of a resistant virus and monitoring response to therapy. Respiratory specimens submitted to the Provincial Laboratory for Public Health (ProvLab), Alberta, Canada, between July 6th 2009 and January 4th 2010 for respiratory virus testing, were included in this analysis. The current diagnostic algorithm for testing of influenza A in our laboratory includes primary screening using a real-time RT-PCR assay (InfA) targeting the matrix gene of all human influenza A viruses (Dawood et al., 2009) . Samples that test positive for IFVA are tested subsequently for H1, H3 and pandemic (H1N1) 2009 sub-types (Dawood et al., 2009; Pabbaraju et al., 2009) . Specimens from patients with a significant travel history, epidemiological links to pandemic (H1N1) 2009 virus, severe respiratory infection, immuno-compromised/bone marrow transplant cases and intensive care unit patients were identified based on the requisitions submitted by the physician. Additionally, patients from the community with none of the above predisposing factors for development of resistance were also included in the analysis. Selected pandemic (H1N1) 2009 virus confirmed samples were tested for oseltamivir resistance using the gold standard Sanger sequencing method (n = 306). Sixty-nine specimens were selected for the assay validation, of which 45 were upper respiratory specimens including nasopharyngeal and throat swabs, 24 were lower respiratory specimens including endotracheal aspirate, bronchoalveolar lavage and auger suction. Extraction of total nucleic acid from all clinical specimens was performed on the easyMAG ® automated extraction system (bioMérieux Canada Inc., Québec, Canada) following the manufacturer's instructions and as described previously (Pabbaraju et al., 2009 ). Primers and two minor-groove-binding (MGB) hydrolysis probes were designed based on the available NA gene sequences from GenBank (as of September 21, 2009) using the Primer Express TM software [version 3.0.0, Applied Biosystems (ABI), Foster City, CA, USA]. The two MGB probes were designed for discrimination of sensitive and resistant alleles. The probe specific for detection of the H275 allele that has been reported to lead to oseltamivir sensitivity (SOIV Osel SEN) was labelled with 6carboxyfluorescein (FAM); and a probe specific for the H275Y allele that has been reported to lead to oseltamivir resistance (SOIV Osel RES) was labelled with VIC at the 5 end; both probes had a non-fluorescent black hole quencher (BHQ) at the 3 end. Nucleotide sequences of the primers and probes are indicated in Table 1 . Optimal primer concentrations were determined by checkerboard titrations using concentrations ranging from 400 nM to 800 nM. The optimal probe concentration was determined by testing concentrations ranging from 100 nM to 200 nM. Primers were synthesized at the Core DNA services, University of Calgary (Alberta, Canada) and the probes were purchased from ABI. Full length amplification of the NA gene with the H275 and H275Y alleles from pandemic (H1N1) 2009 virus was performed using primers described previously (Mak et al., 2006) . These amplified fragments were cloned into the plasmid vector pCRII-TOPO ® and transformed into One Shot ® TOP-10 chemically competent Escherichia coli cells (Invitrogen, Carlsbad, CA, USA). The recombinant plasmids with cloned NA fragments were purified using the QIAprep ® Kit (Qiagen), linearized using the restriction enzyme Hind III and used as templates for transcription. RNA was transcribed using the T7 RiboMAX TM Express kit (Promega, Madison, WI, USA) and treated with TURBO ® DNase (Applied Biosystems Canada, Streetsville, Ontario, Canada) for degradation of template DNA. The purified RNA was quantified using the NanoDrop TM spectrophotometer (Thermo Fisher Scientific, Wilmington, Delaware, USA) and the absorbance value was used to calculate the copy number of transcribed RNA. Amplification was carried out using SuperScript TM III Platinum ® One-step Quantitative Kit (Invitrogen) in a final volume of 25 L with 5 L of extracted RNA as template, 800 nM each of the forward and reverse primers, and 200 nM of each probe. Amplification and detection was performed on the 7500 Sequence Detection System (7500SDS, Applied Biosystem TM (ABI), Foster City, CA, USA). The Table 1 Primer and probe sequences designed for detection of the C823T mutation. The SNP position is highlighted in bold and underlined. Oligonucleotide name Oligonucleotide sequence ( amplification conditions were: reverse transcription at 48 • C for 30 min, Taq activation for 10 min at 95 • C (there was no difference in CT values with activation for 2 min), followed by 45 cycles of amplification comprising of denaturation for 15 s at 95 • C, annealing and primer extension for 1 min at 60 • C. A threshold of 0.3 was used for analysis of amplification curves and the baseline was calculated using default settings on the software (SDS version 1.4). 2.6. Sensitivity, dynamic range, specificity and reproducibility of the RT-PCR Analytical sensitivity of the real-time RT-PCR was determined by testing serial dilutions of quantified in vitro transcribed RNA for both the H275 and H275Y alleles in eight replicates. Serial dilutions of RNA spanning eight logs of template concentration were tested in triplicate to determine the dynamic range of the assay. The primers and probes were designed to be specific for the pandemic (H1N1) 2009 virus sequence, however, they were evaluated experimentally for cross-reactivity with the following targets using extracted DNA/RNA from high titre viral and bacterial cultures or positive respiratory samples as template: IFVA (seasonal H1N1, n = 13; seasonal H3N2, n = 11; H1N2, n = 1; recombinant H3N2, n = 1), influenza B (IFVB) (n = 2), parainfluenza (PIV) (types 1, 2, 3, 4B), respiratory syncytial virus (RSV) (A and B), human coronaviruses (hCoV) (229E, OC43, NL63 and HKU1), human metapneumovirus (hMPV) (lineages 1 and 2), human rhinovirus 1B, coxsackievirus A16, echovirus 2, human adenovirus type 3, human bocavirus, Chlamydophila pneumoniae, Legionella pneumophila and Mycoplasma pneumoniae. The inter-assay and intra-assay variability was assessed using three clinical samples each with the H275 and H275Y alleles at three levels of viral load as determined by the cycle threshold (Ct) values. The samples were tested in triplicate on three independent runs. Sixty-nine samples that had been confirmed previously as positive for pandemic (H1N1) 2009 virus by real-time RT-PCR and sequenced for the NA gene were used for testing assay accuracy. These included 25 sensitive (C823) and 15 resistant (823T) samples. Nine samples containing admixtures of C823 and 823T alleles as determined by close observation of the electropherogram were also included. Twenty samples which had tested negative for IFVA viruses were included as negative controls. In vitro transcribed RNA with the C823 and 823T markers were selected at 2 different concentrations to represent high and low viral titre samples. The copy numbers for high titre RNA were 4.76 × 10 4 and 5.29 × 10 4 per reaction for C823 and 823T, respectively. The copy numbers for low titre RNA were 4.76 × 10 2 and 5.29 × 10 2 per reaction for C823 and 823T, respectively. At each concentration, a mixture of RNA was prepared ranging from 5% to 95% of each allele. The RNA mixtures were tested in triplicate for the detection of mixed clones and estimation of sensitivity for detection when both alleles are present. Conventional Sanger sequencing of the region of interest in the NA gene was used as the gold-standard method for this study. Three hundred and six clinical samples that tested positive for pandemic (H1N1) 2009 virus from immuno-compromised, intensive care unit cases, patients who had received oseltamivir treatment, had a travel history and from the community were selected for sequencing of the partial NA gene including the H275Y (C823T) region. A 415 bp fragment of the NA gene was amplified using primers (NA675F nested sw and NA1089R nested sw) (Table 1) , and the amplified product was purified using the EZ-10 PCR purification kit (Bio Basic Inc., Markham, Ontario, Canada). Sequencing was performed on both strands using the BigDye ® Terminator v3.1 Cycle Sequencing Kit (ABI) in the ABI PRISM ® 3130-Avant Genetic Analyzer. Careful manual examination of the electropherogram at base pair 823 was also performed for the presence of mixed bases. 3.1. Dynamic range, sensitivity, specificity and reproducibility of the assay This RT-PCR assay showed a broad dynamic range using in vitro transcribed RNA. Linear amplification of template was observed between 4.76 × 10 1 to 4.76 × 10 8 and 5.29 × 10 2 to 5.29 × 10 9 copies per reaction for H275 and H275Y alleles, respectively ( Fig. 1a and b) . The assay efficiency was 98.8% and 98.0% for the SOIV Osel SEN probe and SOIV Osel RES probe respectively, and the R 2 values comparing template concentrations and Ct values were 1.0 for both probes. Limit-of-detection limit for the assay was determined by testing 10-fold serial dilutions of quantified RNA in triplicate, resulting in 47.6 copies and 52.9 copies for the H275 and the H275Y alleles per reaction, respectively. The primers and probes did not cross-react with high titre respiratory viruses and bacteria listed in the materials and methods section resulting in 100% specificity. The intra-assay coefficient of variation in Ct values for three clinical samples representing a range of viral loads was tested in triplicate on three independent runs. This parameter ranged from 0.08% to 6.06% with the H275 allele and 0.02-2.95% for H275Y samples. The inter-assay variability ranged from 0.56% to 2.13% and 0.27-1.01% with the H275 and H275Y alleles, respectively. Accuracy of the RT-PCR assay was determined using 25 samples with C823, 15 with 823T and nine with a mixture of C and T nucleotides at position 823 by the gold standard Sanger sequencing method. Additionally, 20 samples that were negative for IFVA were also tested. The 24 specimens containing the H275Y allele and H275/H275Y admixtures populations were from 10 patients belong to immuno-compromised and intensive care unit cases identified in the described screening process previously. The 20 IFVA negative samples were not detected by the assay indicating good primer and probe specificity. Of the 25 samples with the H275 allele, the SOIV Osel SEN probe resulted 22 as positive and three samples tested negative (Table 2a) . These samples had been tested previously by a screening assay for influenza targeting the matrix gene (InfA) (Dawood et al., 2009 ). The Ct values by this screening assay for the three discordant samples were 26.44, 35.48, and 36.79. A negative result suggests either a low viral load or sequence variation in the primer or probe binding region. The sequence of the assay detection region was determined for these samples and analyzed for the presence of mutations. The sample with a Ct value of 35.48 had no change in the probe binding region indicating that the sample possibly tested negative as a result of low viral load. The samples with Ct values of 26.44 and 36.79 were found to have a base pair mutation at different positions in the probe binding region (Fig. 2) . Researchers using this region for the design of probes to detect oseltamivir resistance should be aware of the presence of these mutations in some isolates. These 25 samples were not detected by the SOIV Osel RES probe as expected by Fig. 1. (a) Dynamic range of the developed RT-PCR assay. Amplification curves (Delta Rn vs. Cycle number) and standard curve (Ct vs. log10 concentration) for serial dilutions of the plasmid with C823 nucleotide. Assay dynamic range was linear between template concentrations ranging from 4.76 × 10 1 copies/reaction to 4.76 × 10 8 copies/reactions. (b) Dynamic range of the developed RT-PCR assay. Amplification curves (Delta Rn vs. Cycle number) and standard curve (Ct vs. log10 concentration) for serial dilutions of the plasmid with C823T nucleotide. Assay dynamic range was linear between template concentrations ranging from 5.29 × 10 2 copies/reaction to 5.29 × 10 9 copies/reactions. Performance of RT-PCR assay using SOIV Osel SEN probe in comparison to the gold standard Sanger sequencing method. Performance of RT-PCR assay using SOIV Osel RES probe in comparison to the gold standard Sanger sequencing method. the probe design. Fig. 2 shows an alignment of the two SNP probes and mutations detected in this region. Samples with the H275Y allele were present either as 823T (n = 15) or admixtures of C823 and 823T (n = 9). All of the 15 samples with only the 823T were detected by the SOIV Osel RES probe and not by the SOIV Osel SEN probe. Of the 9 admixtures, the SOIV Osel SEN probe detected eight as positive and the SOIV Osel RES probe detected seven as positive. One admixture that was not detected by the SOIV Osel SEN probe (but detected by the SOIV Osel RES probe) showed a major 823T population and a minor C823 population based on peak height in the electropherogram and also had a Ct value of 35.36 for detection of IFVA, indicating a low viral load (Table 2a) . Two samples with admixtures that had minor populations of the 823T nucleotide present as determined by sequencing were not detected by the SOIV Osel RES probe (Table 2b) . Assay sensitivity for the detection of samples with C823 by the SOIV Osel SEN probe was 88% and detection of 823T by the SOIV Osel RES probes was 100%. Sensitivity for the detection of samples with mixtures of C823 and 823T by SOIV Osel SEN probe and SOIV Osel RES probe was 88.9% and 77.8%, respectively. There was no non-specific detection of the H275 or H275Y alleles, thus both probes were 100% specific (Tables 2a and 2b) . Representative amplification curves for the detection of C823, 823T and a mixture of both nucleotides are shown in Fig. 3 . The sensitive and resistant samples showed a single amplification curve; whereas the admixtures had two detectable amplification curves as indicated in Fig. 3. To test the ability of the assay for detection of both sequence variants when they are present in various ratios and at different viral loads, a mixture of known template concentration of resistant and sensitive in vitro RNA was tested. Two mixtures were prepared, one at a low copy number with 4.76 × 10 2 and 5.29 × 10 2 copies/reaction and the other mixture with 4.76 × 10 5 and 5.29 × 10 5 copies/reaction for the sensitive and resistant RNA, respectively. In vitro RNA with the two sequence variants was mixed in ratios from 5% to 95% at both viral loads. In the mixture prepared using higher copies of in vitro RNA, the SOIV Osel RES probe detected the presence of 823T when its proportion ranged from 20% to 100% of the total template concentration; this probe did not give a positive result when 823T was present at concentrations less than 20%. In the same mixtures, the presence of C823 was detected by SOIV Osel SEN probe when present at a concentrations ranging from 5% to 100% of total template as indicated in Table 3a . For mixtures with low copies of in vitro RNA, the SOIV Osel RES probe detected the presence of 823T when its proportion ranged from 20% to 100% of the total template concentration; this probe did not give a positive result when 823T was present at concentrations less than 20%. In the same mixtures, the presence of C823 was detected by SOIV Osel SEN probe when present at a concentrations ranging from 20% to 100% of total template and presence of C823 at lower than 20% was not detected (Table 3b) . Resistance to oseltamivir in seasonal and pandemic influenza A (H1N1) has been linked primarily to an H275Y polymorphism in NA glycoprotein Sheu et al., 2008) . The H275Y amino acid change was first reported in pandemic (H1N1) virus in June 2009 in an immuno-compromised patient treated with oseltamivir (WHO 2009 briefly note 18), and more cases have since been reported sporadically worldwide (Baz et al., 2009; Chen et al., 2009; Le et al., 2010) . As there are still some questions as to whether or not this H275Y polymorphism impairs viral fitness and transmission (Seibert et al., 2010) , we must be prepared for the possibility that oseltamivir-resistant isolates of the pandemic (H1N1) 2009 virus could expand in a manner similar to that seen in the previous seasonal H1N1 viruses Besselaar et al., 2008; Hurt et al., 2009; Lackenby et al., 2008; Meijer et al., 2009; Sheu et al., 2008) . Screening tools should take into account the low prevalence of oseltamivir-resistance in the pandemic (H1N1) 2009 strain and ideally detection should be approached with cost effective and simple screening methods (Dawood et al., 2009; Harvala et al., 2010; Pabbaraju et al., 2009; Thabet et al., 2010) . The current gold standard methods for detection of resistance include phenotypic detection using the NA-Star ® assay from Applied Biosystems or genetic detection of relevant mutations using Sanger sequencing. Both of these methodologies are labor intensive, time-consuming and expensive. Additionally, phenotypic methods require viral culture which may not be feasible for novel sub-type influenza viruses that are restricted to high-level biosafety environments. Other methods such as reverse-transcriptase PCR coupled with restriction fragment length polymorphism assay (RTPCR-RFLP), pyrosequencing, microarray and multiplex RT-PCR with Luminex detection technologies have also been reported for detection of resistance but require specialized equipment and training (Bolotin et al., 2009; Carr et al., 2008; Deyde et al., 2009 Deyde et al., , 2010 Guo et al., 2009; Hindiyeh et al., 2010; Mahony et al., 2010; Operario et al., 2010; van der Vries et al., 2010) . Real-time RT-PCR using probe technology is now commonplace in molecular diagnostic laboratories; these assays are highthroughput and amenable to large scale testing of patient samples with a quick turn-around time for results. Many laboratories are currently equipped to utilize SNP hydrolysis probes for the detec-tion of oseltamivir resistance in pandemic (H1N1) 2009 virus. This approach has been reported previously for the detection of the C823T mutation in seasonal H1N1 (Bolotin et al., 2009; Carr et al., 2008; Operario et al., 2010) and pandemic (H1N1) 2009 viruses (Hindiyeh et al., 2010; van der Vries et al., 2010) . The possible caveats of using SNP probes to detect H275 and H275Y include the presence of mutations in the probe region adjacent to the SNP positions. This may result in false negative results; the presence of such mutations was detected in this study and has been published previously (Hindiyeh et al., 2010) . To address this issue degenerate probes have been utilized by other investigators, however, the issue of cross reactivity and sensitivity was not investigated thoroughly (Hindiyeh et al., 2010) . SNP assays may also be unable to detect minor populations in heterogenous mixtures of H275 and H275Y alleles within some specimens. Given the concerns about how H275Y SNP assays perform in the context of mixed viral populations and mutations in the probe-binding region, the authors feel that the results of a SNP assay should be interpreted with caution and ambiguous results should be confirmed by sequencing. In conclusion, this manuscript describes the development and validation of a SNP assay for the detection of H275Y mutations associated with oseltamivir resistance in pandemic (H1N1) 2009 virus. The strength of this assay is its ability to identify H275Y mutations rapidly without specialized sequencing equipment. This tool is still relevant in the post-pandemic period as the pandemic strain is currently the circulating seasonal H1N1 influenza A strain of the 2010-2011 respiratory season. Emergence of H274Y oseltamivirresistant A(H1N1) influenza viruses in Japan during the 2008-2009 season Emergence of oseltamivir-resistant pandemic H1N1 virus during prophylaxis Widespread oseltamivir resistance in influenza A viruses (H1N1) Development of a novel real-time reversetranscriptase PCR method for the detection of H275Y positive influenza A H1N1 isolates Rapid molecular detection of the H275Y oseltamivir resistance gene mutation in circulating influenza A (H1N1) viruses Oseltamivir-resistant influenza A pandemic (H1N1) 2009 virus Emergence of a novel swine-origin influenza A (H1N1) virus in humans Surveillance of resistance to adamantanes among influenza A(H3N2) and A(H1N1) viruses isolated worldwide Pyrosequencing as a tool to detect molecular markers of resistance to neuraminidase inhibitors in seasonal influenza A viruses Detection of molecular markers of drug resistance in 2009 pandemic influenza A (H1N1) viruses by pyrosequencing Infections with oseltamivir-resistant influenza A(H1N1) virus in the United States Rapid identification of oseltamivir-resistant influenza A(H1N1) viruses with H274Y mutation by RT-PCR/restriction fragment length polymorphism assay The emergence of oseltamivir-resistant pandemic influenza A(H1N1) 2009 virus amongst hospitalised immunocompromised patients in Scotland Oseltamivirresistant influenza viruses A (H1N1), Norway Rapid detection of influenza A pandemic (H1N1) 2009 virus neuraminidase resistance mutation H275Y by real-time reverse transcriptase PCR Emergence and spread of oseltamivir-resistant A(H1N1) influenza viruses in Oceania Selection for resistance to oseltamivir in seasonal and pandemic H1N1 influenza and widespread co-circulation of the lineages Emergence of resistance to oseltamivir among influenza A(H1N1) viruses in Europe A community cluster of oseltamivir-resistant cases of 2009 H1N1 influenza Development of a novel bead-based multiplex PCR assay for combined subtyping and oseltamivir resistance genotyping (H275Y) of seasonal and pandemic H1N1 influenza A viruses Longitudinal analysis of genotype distribution of influenza A virus from Oseltamivir-resistant influenza virus A (H1N1) The origin and global emergence of adamantane resistant A/H3N2 influenza viruses Neuraminidase inhibitor susceptibility testing in human influenza viruses: a laboratory surveillance perspective Highly sensitive and quantitative detection of the H274Y oseltamivir resistance mutation in seasonal A/H1N1 influenza virus Design and validation of real-time reverse transcription-PCR assays for detection of pandemic (H1N1) 2009 virus Adamantane resistance in seasonal human influenza A viruses from Calgary Oseltamivir-resistant variants of the 2009 pandemic H1N1 influenza A virus are not attenuated in the guinea pig and ferret transmission models Surveillance for neuraminidase inhibitor resistance among human influenza A and B viruses circulating worldwide from Evaluation of a rapid molecular algorithm for detection of pandemic influenza A (H1N1) 2009 virus and screening for a key oseltamivir resistance (H275Y) substitution in neuraminidase First oseltamivir-resistant pandemic (H1N1) 2009 influenza from EMR Oseltamivir resistance in immunocompromised hospital patients pandemic (H1N1) 2009 breifing note 18 Genetic makeup of amantadine-resistant and oseltamivir-resistant human influenza A/H1N1 viruses We would like to thank the staff in Molecular Diagnostics and Virology departments at Provincial Laboratory for Public Health (ProvLab) for their technical assistance in undertaking routine respiratory virus testing. In particular, we thank Ms. Kara Gill and Ms Danielle Zarra for their excellent technical support for sequencing of samples.