Emerg Infect DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention22264682331009711-051510.3201/eid1801.110515DispatchDispatchMolecular Evolution of Respiratory Syncytial Virus Fusion Gene, Canada, 2006–2010Molecular Evolution of RSV Fusion GenePapenburgJesseCarbonneauJulieHamelinMarie-ÈveIsabelSandraBouhyXavierOhoumanneNajwaDéryPierrePaesBosco A.CorbeilJacquesBergeronMichel G.De SerresGastonBoivinGuyMcGill University Health Centre, Montréal, Québec, Canada (J. Papenburg);Centre de Recherche du Centre Hospitalier Universitaire de Québec, Québec City, Québec, Canada (J. Papenburg, J. Carbonneau, M.-È. Hamelin, S. Isabel, X. Bouhy, P. Déry, J. Corbeil, M.G. Bergeron, G. Boivin);Institut National de Santé Publique du Québec, Québec (N. Ohoumanne, G. De Serres);McMaster Children’s Hospital, Hamilton, Ontario, Canada (B.A. Paes)Address for correspondence: Guy Boivin, CHUL, Rm RC-709, 2705 Blvd Laurier, Québec City, Québec G1V 4G2, Canada; email: guy.boivin@crchul.ulaval.ca12012181120124

To assess molecular evolution of the respiratory syncytial virus (RSV) fusion gene, we analyzed RSV-positive specimens from 123 children in Canada who did or did not receive RSV immunoprophylaxis (palivizumab) during 2006–2010. Resistance-conferring mutations within the palivizumab binding site occurred in 8.7% of palivizumab recipients and none of the nonrecipients.

Keywords: Human respiratory syncytial viruschildrenfusion (F) proteinphylogenypalivizumabvirusesCanada

Human respiratory syncytial virus (RSV) is the most common cause of acute respiratory tract infections (RTIs) and a major cause of hospital admission and death among children <5 years of age worldwide (1). Risk for severe RSV-associated illness is highest among children born prematurely or with chronic medical disorders (2). Palivizumab immunoprophylaxis is the only available measure to prevent severe RSV disease.

The RSV fusion (RSV-F) surface glycoprotein mediates virus fusion to host cells. It is a major antigenic determinant that elicits neutralizing antibodies and cytotoxic T-lymphocyte immunity (3). Palivizumab (MedImmune, Gaithersburg, MD, USA) is a humanized mouse monoclonal antibody that inhibits RSV-F by binding to a defined epitope (residues 262–276) (4,5). Palivizumab immunoprophylaxis is recommended for the prevention of serious lower RTIs caused by RSV in children at high risk (6). RSV strains with mutations in key amino acid residues within the palivizumab binding site are resistant to this antibody (79); however, little is known about the prevalence of such mutations in clinical samples. Furthermore, despite its role in RSV pathogenesis, immunity, and prevention strategies, few data on RSV-F molecular evolution are available (10,11) because previous phylogenetic studies have focused on the RSV-G glycoprotein (12,13). Therefore, we monitored evolutionary changes in RSV-F, particularly potential resistance mutations in the palivizumab binding site, among strains from children who did and did not receive palivizumab.

The Study

This cohort study was approved by the Centre Hospitalier Universitaire de Québec Research Ethics Board. Participants were <3 years of age and either received medical attention at an outpatient pediatric clinic or were hospitalized at a pediatric center for acute RTI during 4 winter seasons (2006–2010), in Québec City, Québec, Canada.

Clinical data were prospectively collected at study entry and after 1-month follow-up. For all patients, at the first visit a nasopharyngeal aspirate was collected. The aspirate was frozen at −80°C until subsequent testing by a multiplex PCR/DNA hybridization assay that detects RSV genotype-A (RSV-A), RSV-B, and 22 other respiratory viruses (Infiniti RVP assay; Autogenomics, Carlsbad, CA, USA) (14).

RSV infection was identified in aspirates from 467 (63.6%) of 734 hospitalized children (257 RSV-A, 210 RSV-B) and from 147 (48.2%) of 305 outpatient children (85 RSV-A, 62 RSV-B). During 2006–2010, a total of 724 children received palivizumab in the Québec City region (L. Cliche, pers. comm.). RSV-positive samples from all 12 study participants receiving palivizumab and from 100 not receiving palivizumab underwent RSV-F sequencing.

Additionally, F-gene analysis was performed on 11 RSV-positive clinical samples from palivizumab recipients retrospectively identified by using neonatal clinic registries at McMaster Children’s Hospital (Hamilton, Ontario, Canada) and Montréal Children’s Hospital (Montréal, Québec, Canada) during 2009–2010. Clinical data were collected by chart review.

RNA was extracted directly from nasopharyngeal samples by using a QIAmp Viral RNA Mini Kit (QIAGEN, Mississauga, Ontario, Canada). Random primers (Amersham, Piscataway, NJ, USA) and Superscript II RT Kit (Invitrogen, Carlsbad, CA, USA) were used for reverse transcription. PCR amplification was performed with QuantiFast Probe PCR+ROX Vial Kit (QIAGEN); primers and thermocycling conditions are available from G.B. upon request. RSV-F amplicons were sequenced by using an automated sequencer (Applied Biosystems, Foster City, CA, USA).

Along with the newly generated nucleotide sequences from the 123 children (23 palivizumab recipients; 100 nonrecipients), we analyzed 92 clinical RSV-F sequences provided by other investigators (10,11) or available from GenBank. Palivizumab exposure was unknown for these samples, and all originated outside Canada. We also included 10 unpublished RSV-F sequences from a 2004–2005 study of palivizumab recipients in Canada (15). Multiple-sequence alignment was performed with ClustalW in MEGA5 (www.megasoftware.net). A 1,524-bp region (positions 79–1602 in the prototype A2 and B1 RSV-F genes) was translated into amino acid sequences, and the palivizumab binding site (residues 262–276) was assessed for variations. We removed 55 redundant (identical) nucleotide sequences. The final dataset comprised 170 unique sequences, including 89 that were newly generated in this study (GenBank accession nos. JF776691–JF776779).

Phylogenetic reconstructions were computed in MEGA5 (Figure). Strains from Canada tended to group together. However, this segregation likely reflects temporal evolution rather than geographic influence because all RSV-F sequences originating from other countries were collected >10 years ago. Previous RSV-G studies have demonstrated concurrent circulation of related lineages in distant areas (13). Overall, RSV-F was highly conserved (Table 1). RSV-A exhibited more genetic diversity than RSV-B. This difference may reflect sampling bias because more RSV-A sequences originating from diverse locations and years were available; however, the greater genetic diversity of RSV-A compared with RSV-B has also been observed in a study conducted in South Korea (10).

Phylogenetic analysis of 170 near–full-length unique respiratory syncitial virus fusion (RSV-F) gene sequences (nt 79–1602). Panels A and B are detailed phylograms of the RSV-A and RSV-B taxa analyzed, respectively. One bovine RSV-F sequence was added to the dataset (GenBank accession no. AF295543.1) as the outgroup (not shown) and used for rooting the phylograms. Topology was inferred by using the neighbor-joining method, and evolutionary distances were computed by using the maximum-composite likelihood method in MEGA5 software (www.megasoftware.net). The topologic accuracy of the tree was evaluated by using 1,000 bootstrap replicates. Only bootstrap values >50% are shown. Blue text and triangles represent RSV strains isolated in Canada and sequenced at the Centre de Recherche du Centre Hospitalier Universitaire de Québec; red branches indicates a sublineage of RSV-A with a N276S mutation; underlining indicates prototypical RSV A2 and RSV B1 strains. The clinical origin of strains from Canada (prospective study hospitalized patient [H] or clinic patient [C]; 2004–2005 palivizumab study patient [S]; retrospectively identified patient from the Montréal Children’s Hospital [MCH] or McMaster Children’s Hospital [MAC]) is indicated, followed by the year collected, the specimen identifier, and the result of RSV genogroup testing (RSV-A [A], RSV-B [B]) by multiplex PCR/DNA hybridization assay (14). Specimens with a nonsilent mutation at codon 272 have the amino acid substitution identified in brackets. When a taxon represents >1 identical sequences, the number of patients that it represents and the number of palivizumab recipients (PZB) among these patients are noted in parentheses. SA, South Africa; CHN, People’s Republic of China; UK, United Kingdom, JPN, Japan; AUS, Australia; USA, United States; CDC, US Centers for Disease Control and Prevention; NY, New York; SL, St. Louis; WV, West Virginia; MN, Minnesota; MD, Maryland; WA, Washington; SEL, Seoul, South Korea. Scale bars represent substitutions per basepair per the indicated horizontal distance.

Estimated mean nucleotide and amino acid identities of 170 unique respiratory syncytial virus fusion gene sequences*
SequenceMean nucleotide identity, % ± SDMean amino acid identity, % ± SD
Overall, n = 17089.59 ± 8.3895.99 ± 3.25
RSV-A (within group), n = 10597.06 ± 1.3398.85 ± 0.76
RSV-B (within group), n = 6598.65 ± 0.6999.48 ± 0.28
RSV-A vs. RSV-B (between groups), n = 17080.85 ± 0.6692.64 ± 0.56

*Mean nucleotide and amino acid identities were estimated by calculating pairwise distances in MEGA5 software (www.megasoftware.net) and using the maximum-composite likelihood and Poisson correction models, respectively. RSV, respiratory syncytial virus.

Of the 23 RSV-F amino acid sequences from patients who received palivizumab (Table 2), 2 exhibited a mutation at residue 272 (K272Q, K272M), whereas no such mutations were identified in the other 100 strains obtained from those who did not receive palivizumab (8.7% versus 0%; p = 0.03, by Fisher exact test). A sublineage of RSV-A, noteworthy for an N276S substitution in the palivizumab binding site, emerged during 2008–2009 (8 [44%] of 18 RSV-A strains) and became the predominant RSV-A clade during 2009–2010 (25 [100%] of 25 strains) in palivizumab recipients and nonrecipients. This N276S lineage was detected in all 3 Canadian communities studied. No sequences from GenBank or other studies (10,11) harbored any palivizumab binding site mutation.

Characteristics of 23 children in whom clinically significant respiratory syncytial virus respiratory tract infection developed while receiving palivizumab immunoprophylaxis, Canada, 2006–2010*
Location and patient IDAge, mo/ sex†GA at birth, wk + dUnderlying
comorbiditiesNo. doses PZB‡Delay, d§Clinical diagnosesHMultiplex PCR/DNA resultsMutation
Québec City, Québec (2006–2010)
C0607-10239/F32 + 4Prematurity, LBW321BronchiolitisNoRSV-A;
enterovirus type AK272Q
H0607-06424/M38 + 3Congenital myopathy315Pneumonia;
bronchospasmYesRSV-BNF
H0607-13212/M38 + 5Pulmonary artery stenosis57BronchiolitisYesRSV-ANF
H0708-1994/M30 + 4Prematurity, VLBW414BronchiolitisYesRSV-BNF
H0809-03711/F27 + 5Prematurity, ELBW314BronchiolitisYesRSV-ANF
C0809-10556/F29 + 0Prematurity, ELBW, triplet427BronchiolitisNoRSV-AN276S
C0809-10566/M29 + 0Prematurity, ELBW, triplet427BronchiolitisNoRSV-AN276S
C0809-10576/M29 + 0Prematurity, VLBW, triplet427BronchiolitisNoRSV-AN276S
H0910-0044/F39 + 5Choanal hypoplasia116Apnea;
upper RTIYesRSV-AN276S
H0910-1406/M25 + 5Prematurity, ELBW429BronchiolitisYesRSV-BNF
H0910-14413/F26 + 5Prematurity, VLBW47PneumoniaYesRSV-AK272M, N276S
H0910-1509/M28 + 3Prematurity, VLBW412Upper RTI;
acute otitis mediaYesRSV-A;
adenovirus type CN276S
Montréal, Québec (2009–2010)
MCH0910-00115/M40 + 4Total anomalous pulmonary venous return326PneumoniaYesRSV¶N276S
MCH0910-0026/F39 + 0Pulmonary valve stenosis, right aortic arch27BronchiolitisYesRSV¶N276S
MCH0910-0035/M39 + 6Cystic fibrosis324BronchiolitisNoRSV¶N276S
MCH0910-0047/M36 + 2Prematurity, BPD hypotonia46BronchiolitisYesRSV¶N276S
MCH0910-00515/M40 + 4Neuromuscular disorder, recurrent aspirations413Upper RTI;
acute otitis mediaNoRSV¶N276S
MCH0910-0062/M34 + 6Prematurity, LBW114BronchiolitisYesRSV¶N276S
MCH0910-00719/F25 + 0Prematurity, ELBW, BPD319Bronchiolitis, bronchospasmYesRSV¶N276S
MCH0910-0082/F38 + 1Neuromuscular disorder, ventricular septal defect212BronchiolitisYesRSV¶N276S
Hamilton, Ontario (2009–2010)
MAC0910-0011/F34 + 5Prematurity, LBW23BronchiolitisYesRSV¶N276S
MAC0910-0026/F34 + 3Prematurity, LBW, twin125BronchiolitisYesRSV¶N276S
MAC0910-0036/F34 + 3Prematurity,
VLBW, IUGR, twin127BronchiolitisYesRSV¶N276S

*Patient identification (ID) nomenclature: hospitalized (H) or clinic (C) prospective study participant, Montréal Children’s Hospital (MCH) or McMaster Children’s Hospital (MAC) patient. GA, gestational age; PZB, palivizumab; multiplex PCR/DNA, hybridization assay; mutation, mutation in respiratory syncytial virus fusion protein PZB binding site (residues 262–276); RSV, respiratory syncytial virus; LBW, low birthweight (1,500–2,500 g); NF, no mutation found in PZB binding site; VLBW, very low birthweight (1,000–1,499 g); ELBW, extremely low birthweight (<1,000 g); RTI, respiratory tract infection; BPD, bronchopulmonary dysplasia; IUGR, intrauterine growth restriction.
†Median patient age 6.0 mo (range 1–24 mo).
‡Mean ± SD no. palivizumab doses received that winter: 3.0 ± 1.2 doses.
§Median interval between last palivizumab dose and symptom onset: 15.0 d (range 3–27 d).
¶Retrospectively identified participants from Montréal Children’s Hospital or McMaster Children’s Hospital were RSV-positive by direct immunofluorescence assay (Chemicon International, Temecula, CA, USA) and were not tested by the multiplex PCR/DNA hybridization assay (14).

Microneutralization assays were performed as described elsewhere, with minor modifications (9). RSV was incubated (for 2 h at 37°C) with serially diluted palivizumab, then cultured in Vero cells (for 5 d). RSV replication and 50% inhibitory concentrations (IC50) were subsequently determined by F protein quantification by using ELISA. The mean ± SD IC50 of C0910–1006A, a N276S strain, was 0.33 ± 0.04 µg/mL, similar to that of RSV-A2 wild type (IC50 0.46 ± 0.04 µg/mL) and therefore was considered susceptible. Position-272 variants did not grow in culture and were not tested.

Conclusions

We report the prevalence of resistance-conferring mutations in RSV-F among children receiving or not receiving palivizumab. Although infrequent (8.7% of infections in palivizumab recipients), residue 272 mutations were significantly associated with palivizumab exposure and not observed at all in nonexposed patients.

We identified 2 new clinical specimens with position 272 mutations (K272Q and K272M). We cannot exclude the possibility that additional specimens contained mixed viral populations with minor proportions of position-272 mutants not detectable by conventional sequencing methods. Changes at this position (from lysine to asparagine, glutamine, glutamic acid, methionine, or threonine) have produced palivizumab resistance in vitro (9) and in a cotton rat model (7). As previously reported, 1 (10.0%) of 10 sequences from our 2004–2005 study of palivizumab patients carried a 272 mutation (K272E) (15). The K272E substitution is the only substitution also demonstrated to confer resistance to motavizumab, an enhanced-potency monoclonal antibody developed by affinity maturation of palivizumab (9). We could not perform neutralization assays on position-272 variants because they did not grow in culture. This finding suggests that such changes adversely affect viral replicative capacity (7,9).

Phylogenetic analysis demonstrated that mutations in the palivizumab binding site occur in diverse genetic backgrounds; all 3 strains with substitutions at residue 272 grouped to different clades (Figure). Furthermore, these mutant strains caused mild disease treatable in an outpatient clinic and severe illness requiring hospitalization.

From 3 Canadian communities we detected a lineage harboring an N276S mutation in 44.4% of RSV-A sequences from 2008–2009 and 100% from 2009–2010, unrelated to palivizumab exposure. Adams et al. have proposed that N276S led to palivizumab resistance in a clinical specimen (8). However, that sample also comprised a K272E subpopulation. Our microneutralization assay results and unpublished neutralization data using recombinant viruses and clinical isolates (Q. Zhu, pers. comm.) suggest that N276S does not confer resistance.

Although serious RSV RTIs during palivizumab prophylaxis remain uncommon, we observed an 8.7% prevalence of known resistance mutations among 23 medically attended patients receiving palivizumab. These findings underscore the need for continued monitoring of RSV-F evolution.

Suggested citation for this article: Papenburg J, Carbonneau J, Hamelin M-E, Isabel S, Bouhy X, Ohoumanne N, et al. Molecular evolution of respiratory syncytial virus fusion gene, Canada, 2006–2010. Emerg Infect Dis [serial on the Internet]. 2012 Jan [date cited]. http://dx.doi.org/10.3201/eid1801.110515

Acknowledgments

We thank Patricia Fontela for her constructive review of this manuscript; Jane McDonald, Lorraine Piché, and Diane Newby for providing the clinical samples from patients at the Montréal Children’s Hospital; James Mahony and Sylvia Chung for providing the clinical samples from patients at McMaster Children’s Hospital; Susan Steele for performing medical chart reviews at McMaster Children’s Hospital; Yun Kyung Kim and Hoan Jong Lee for sharing previously unpublished RSV-F sequences from the Seoul National University; Hao Qiang Zheng for sharing previously unpublished RSV-F sequences from the US Centers for Disease Control and Prevention; and Marie-Hélène Cavanagh and Laetitia Aerts for technical assistance.

This study was supported by the Canadian Institutes of Health Research and a research grant from Abbott Laboratories Ltd.

Dr Papenburg is a pediatric infectious disease specialist and medical microbiologist who performed this work during a fellowship at the Centre de Recherche en Infectiologie de l’Université Laval. His clinical practice is now at the Montréal Children’s Hospital, McGill University. His research interests include the clinical epidemiology and molecular virology of childhood respiratory tract infections.

ReferencesNair H, Nokes DJ, Gessner BD, Dherani M, Madhi SA, Singleton RJ, Global burden of acute lower respiratory infections due to respiratory syncytial virus in young children: a systematic review and meta-analysis. Lancet. 2010;375:154555 10.1016/S0140-6736(10)60206-120399493Welliver RC Review of epidemiology and clinical risk factors for severe respiratory syncytial virus (RSV) infection. J Pediatr. 2003;143:S1127 10.1067/S0022-3476(03)00508-014615709Collins PL, Crowe JE Respiratory syncytial virus and metapneumovirus. In: Knipe DM, Howley PM, editors. Field’s virology.,5th ed. New York: Lippincott Williams and Wilkins; 2007 p. 1601–46.Crowe JE, Firestone CY, Crim R, Beeler JA, Coelingh KL, Barbas CF, Monoclonal antibody-resistant mutants selected with a respiratory syncytial virus–neutralizing human antibody fab fragment (Fab 19) define a unique epitope on the fusion (F) glycoprotein. Virology. 1998;252:3735 10.1006/viro.1998.94629878616Huang K, Incognito L, Cheng X, Ulbrandt ND, Wu H Respiratory syncytial virus–neutralizing monoclonal antibodies motavizumab and palivizumab inhibit fusion. J Virol. 2010;84:813240 10.1128/JVI.02699-0920519399Committee on Infectious Diseases From the American Academy of Pediatrics: policy statements–modified recommendations for use of palivizumab for prevention of respiratory syncytial virus infections. Pediatrics. 2009;124:1694701 10.1542/peds.2009-234519736258Zhao X, Chen FP, Megaw AG, Sullender WM Variable resistance to palivizumab in cotton rats by respiratory syncytial virus mutants. J Infect Dis. 2004;190:19416 10.1086/42551515529258Adams O, Bonzel L, Kovacevic A, Mayatepek E, Hoehn T, Vogel M Palivizumab-resistant human respiratory syncytial virus infection in infancy. Clin Infect Dis. 2010;51:1858 10.1086/65353420528472Zhu Q, McAuliffe JM, Patel NK, Palmer-Hill FJ, Yang CF, Liang B, Analysis of respiratory syncytial virus preclinical and clinical variants resistant to neutralization by monoclonal antibodies palivizumab and/or motavizumab. J Infect Dis. 2011;203:67482 10.1093/infdis/jiq10021208913Kim YK, Choi EH, Lee HJ Genetic variability of the fusion protein and circulation patterns of genotypes of the respiratory syncytial virus. J Med Virol. 2007;79:8208 10.1002/jmv.2089117457915Zheng H, Storch GA, Zang C, Peret TC, Park CS, Anderson LJ Genetic variability in envelope-associated protein genes of closely related group A strains of respiratory syncytial virus. Virus Res. 1999;59:8999 10.1016/S0168-1702(98)00132-410854168Gilca R, De Serres G, Tremblay M, Vachon ML, Leblanc E, Bergeron MG, Distribution and clinical impact of human respiratory syncytial virus genotypes in hospitalized children over 2 winter seasons. J Infect Dis. 2006;193:548 10.1086/49852616323132Peret TC, Hall CB, Hammond GW, Piedra PA, Storch GA, Sullender WM, Circulation patterns of group A and B human respiratory syncytial virus genotypes in 5 communities in North America. J Infect Dis. 2000;181:18916 10.1086/31550810837167Raymond F, Carbonneau J, Boucher N, Robitaille L, Boisvert S, Wu WK, Comparison of automated microarray detection with real-time PCR assays for detection of respiratory viruses in specimens obtained from children. J Clin Microbiol. 2009;47:74350 10.1128/JCM.01297-0819158263Boivin G, Caouette G, Frenette L, Carbonneau J, Ouakki M, De Serres G Human respiratory syncytial virus and other viral infections in infants receiving palivizumab. J Clin Virol. 2008;42:527 10.1016/j.jcv.2007.11.01218164233