key: cord-0820845-akurajpf
authors: Legrand, Loïc; Vabret, Astrid; Dina, Julia; Petitjean‐Lecherbonnier, Joëlle; Stéphanie, Gouarin; Cuvillon, Delphine; Tripey, Valérie; Brouard, Jacques; Freymuth, François
title: Epidemiological and phylogenic study of human metapneumovirus infections during three consecutive outbreaks in Normandy, France
date: 2011-01-24
journal: J Med Virol
DOI: 10.1002/jmv.22002
sha: 6acc9cc8d7764e9e7c9a4f109525e64b971c1d21
doc_id: 820845
cord_uid: akurajpf

Human metapneumovirus (hMPV) is responsible for respiratory tract disease, particularly in the young and elderly population. An epidemiological and phylogenic study was performed on children admitted to hospital with an acute lower respiratory tract infection (LRI). Data were obtained and analyzed over three consecutive winters, from 2002–2003 to 2004–2005. Each year during the winter period, from November to March, 2,415 nasal swabs were tested by a direct immunofluorescence assay (DFA) for influenza viruses A and B, respiratory syncytial virus, parainfluenza viruses, and adenoviruses. Rhinoviruses, enteroviruses, and coronaviruses OC43 and 229E were detected by RT‐PCR. A RT‐PCR designed for the M gene was performed on negative samples for hMPV detection and phylogenic analyses. For the three consecutive winters, hMPV represented 10%, 22.6%, and 8.8% of virus‐negative samples, respectively. In most cases, clinical symptoms indicated a LRI with a final diagnosis of bronchiolitis. During the winter of 2003–2004, all viral clusters (A1, A2, B1, and B2) that circulated in France shifted progressively from the A group to the B group. This study determined the prevalence of hMPV in Normandy, its clinical impact and permitted the analysis of the molecular evolution during the successive outbreaks. J. Med. Virol. 83:517–524, 2011. © 2011 Wiley‐Liss, Inc.

Human metapneumovirus (hMPV) was first isolated by a Dutch group in 2001 [van den Hoogen et al., 2001] . Since then, the virus has been detected in all continents [Nissen et al., 2002; Peret et al., 2002; Stockton et al., 2002; Cuevas et al., 2003; Freymuth et al., 2003; Peiris et al., 2003; Fodha et al., 2004] . Human metapneumovirus is the single human pathogen within the genus Metapneumovirus, subfamily Pneumovirinae, family Paramyxoviridae. Serological studies suggest that the virus has been circulating in humans for at least the past 60 years [van den Hoogen et al., 2001; Ebihara et al., 2003; Mejias et al., 2004] . However, recent evidence from molecular studies suggests that hMPV may have diverged from the avian metapneumovirus type C as early as 200 years ago [de Graaf et al., 2008; Yang et al., 2009] .

Human metapneumovirus causes upper and lower respiratory tract diseases such as bronchiolitis or pneumonia Williams et al., 2004 Williams et al., , 2006 Bouscambert-Duchamp et al., 2005; Gray et al., 2006; Sarasini et al., 2006] . Human infections occur mainly during the winter period and affect, in particular, young children, elderly people, and immunocompromised patients Bastien et al., 2003a; Falsey et al., 2003] . It was shown serologically that children under 5 years old were the most vulnerable and susceptible to contracting the virus [van den Hoogen et al., 2001; Ebihara et al., 2003; Regev et al., 2006] . Furthermore, hMPV can be detected in respiratory samples by reverse transcription-PCR (RT-PCR) [Mackay et al., 2003; Maertzdorf et al., 2004; Bellau-Pujol et al., 2005] or real-time RT-PCR [Hopkins et al., 2008; Kaida et al., 2008] . In addition, several monoclonal antibodies were developed for immunofluorescence assays (IFA) [Ebihara et al., 2005; Landry et al., 2005; Gerna et al., 2006; Aslanzadeh et al., 2008] . Molecular studies have shown that hMPV strains can be divided into two major groups: A and B, and 4 minor groups: A1, A2, B1, and B2 [Bastien et al., 2003b; Boivin et al., 2004; .

This paper describes studies concerning hMPV infection and reports on a genetic analysis of hMPV strains that were detected in infected children during three consecutive winters in Normandy, France.

The study was performed on 2,415 children admitted to the University Hospital of Caen and the hospital of Flers, for an acute lower respiratory tract infection (LRI) from November to March during three consecutive years (November 2002 to March 2004 . All patients underwent a nasal aspirate which was processed immediately by conventional methods (see below). The remaining samples were frozen at À708C and used for subsequent molecular analyses. All samples were collected and tested in conformation with the ethical legislation in place.

The conventional methods included a direct immunofluorescence assay (DFA). This immunofluorescence assay used commercially available monoclonal antibodies (Imagen 1 , Dako, Glostrup, Denmark) and permitted the detection of influenza virus types A and B, respiratory syncytial virus (RSV), parainfluenza virus types 1, 2 and 3, and adenoviruses. When the DFA was negative, samples were inoculated onto HuH7 cell lines grown in 48-well microplates as described previously . Briefly, after four days of incubation, HuH7 cells were harvested and used for DFA with the same antiviral antibodies. When a cytopathic effect was observed, supernatants of the relevant wells were harvested and used in a multiplex PCR for the detection of rhinovirus, enterovirus and human coronaviruses 229E, OC43 [Bellau-Pujol et al., 2005] .

RNA was extracted from nasal aspirates with the QIAamp Viral RNA minikit (Qiagen, Hilden, Germany) according to the manufacturer's instructions. RT-PCR was likewise performed with the one-step RT-PCR kit (Qiagen) using primers designed for the M gene. The forward primer hMPV.M1 (5 0 -CCCTTTGTTTCAGGCC-AAC-3 0 ) and the reverse primer hMPV.M2 (5 0 -GCA-GCTTCAACAGTAGCTG-3 0 ) provided a 416 pb PCR product. The thermocycler conditions were 30 min at 508C, 15 min at 958C following by 40 cycles of 1 min at 958C, 1 min at 608C and 1 min at 728C. Finally, an extension at 728C was performed during 10 min. The result was analyzed on 2% agarose gel.

The PCR sensitivity for each group was established using two fragments of the M gene containing the amplification product (Eurofin-MWG-Operon, Ebersberg, Germany) representing NL00-1 (AF371337) and NL99-1 (AY525843), respectively. NL00-1 was considered like group A prototype sequence and NL99-1 like group B prototype sequence. Synthesized DNA were cloned in a pCR2.1 vector and transcription was performed using T7 polymerase (Applied Biosystems/Ambion, Courtaboeuf, France). After DNase treatment (Applied Biosystems/ Ambion), RNA concentration was determined by measuring the absorbance at 260 nm with the Nanodrop spectrophotometer system. Serial 10-fold dilutions were performed to determine the sensitivity of the assay in copy of the viral RNA per reaction.

An hemi-nested PCR was performed with the reverse primer MPV.M2 and a new forward primer MPV.M3 (5 0 -AGGCCAACACACCACCAG-3 0 ). The amplification was realized through the use of an AmpliTaq (Applied Biosystem, Foster City, CA) with the following program: 5 min at 958C followed by 40 cycles of 30 sec at 958C, 30 sec at 608C, and 45 sec at 728C and a final elongation at 728C for 10 min. The 405 pb product obtained was first purified with the enzyme ExoSap-IT (Amersham Biosciences, Uppsala, Sweden) according to the manufacturer's instructions and then quantified on a spectrometer. Subsequently, the purified PCR product was sequenced on a Beckman Ceq TM 8000 genetic analysis system with the CEQ DTCS Quick start kit (Beckman Coulter, Brea, CA). The sequences were assembled and analyzed with BioEdit and ClustalW software. A phylogenic study was carried out using Mega 3.1 software with the neighbor joining method. A bootstrap was performed with a replicate rate of 1,000.

The nucleotide sequences reported in this study were submitted to the GenBank under accession numbers H384124 to H384159.

The RT-PCR allowed detection of about 30 and 50 copies per reaction for the A and B group, respectively. This method was performed on virus-negative nasal aspirates. The frequency of hMPV infection in hospitalized children could be assessed only during the first and last winters because the study performed during the winter of 2003-2004 was limited to children under 2 years of age; this study was described by Freymuth et al. [2006] . In order to compare the hMPV prevalence to the other common cold viruses, the period from November to March was studied exclusively. In 2002 and 2004, there were 498 and 642 samples, respectively, that contained a respiratory virus out of the 1,040 and 1,112 nasal aspirates tested (Table I) . Human metapneumovirus was detected in 59 (10%) and 52 (8.8%) of the virus-negative specimens analyzed, respectively. As a function of age, there were 34.2%, 22.5%, 16.4%, and 26.9% cases of infected children between 0 and 6 months, 6-12 months, 12-24 months, and older than 24 months, respectively (Table II) .

During the winters of 2002-2004, other viruses were scrutinized for by a DFA test or by using culture models or molecular tools. Human metapneumovirus was the fourth most commonly detected virus (8.8%). Human respiratory syncytial virus (hRSV) was the most predominant and represented 45.6% of virus-positive specimens followed by rhinoviruses (12.7%) and the influenza A virus (11.3%). The other detected viruses were parainfluenza viruses 1, 2, and 3 (5.6%); adenoviruses (5.3%); coronaviruses (4.8%); enteroviruses (3%); and the influenza B virus (2.8%). In children under 2 years old, hMPV (8.5%) became the third most frequently encountered virus after hRSV (50.7%) and rhinoviruses (14.2%). Influenza A virus represented about 7.7% of the cases studied; adenoviruses 5.5%; parainfluenza viruses 1%, 2%, and 3 5%; coronaviruses 4.4%; enteroviruses 2.8%; and influenza B virus 1.2% (Table I) .

Clinical information was available from 46 out of the 52 hMPV-infected children during the winter of 2004. The mean age of the infected children was 21 months (Table II ) and the male/female ratio was equal. Children were hospitalized about 3.5 days after the onset of symptoms. The most frequent clinical signs at admission were rhinitis (84.8%), cough (71.7%), respiratory distress syndrome (30.4%), and otitis (17.4%). Nineteen children (41.3%) presented with an upper respiratory tract infection and 26 children (56.5%) with a LRI. The final clinical diagnosis was bronchiolitis in 17 cases (37%), broncho-pneumonia in six cases (13%), and acute exacerbation of asthma in three cases (6.6%). Signs of gastroenteritis were present in 17 (30.4%) of hMPVinfected children and rotavirus or adenovirus were present in nine cases. Seven hMPV-infected children (15.2%) required oxygen therapy.

During the winter of 2002-2003, the hMPV and hRSV epidemics appeared at the end of November and they peaked almost immediately at 3-4 weeks later. In contradistinction, during the winters of 2003-2004 and 2004-2005 , the hMPV and hRSV epidemics appeared Data in brackets concern children less than 2 years old.

b

These data were obtained only in children less than 2 years old. These data were obtained only in children less than 2 years old.

Phylogenic Study of Human Metapneumovirus later and continued more slowly. Human metapneumovirus seemed to have a constant presence without any obvious peak in activity (Fig. 1) .

A phylogenetic analysis was carried out on the M gene of 107 hMPV strains isolated during the three consecutive outbreaks (Fig. 2) . The hMPV strains were distributed into two groups, A and B, and into four relevant subgroups: A1, A2, B1, and B2. The percentages of identity in the M gene between the hMPV group A and B strains varied from 84.2% to 88.8%. The percentages of genetic similitude between the strains of subgroups A1 and A2 and between the strains of B1 and B2 subgroups were 94.1-96.3% and 93.8-95.7%, respectively.

Human metapneumovirus strains of different groups and subgroups were associated with each outbreak (Fig. 3a) . However, there was often one dominant subgroup in each epidemic. Human metapneumovirus belonging to the A2 subgroup and B1 subgroup were predominant during the winters of 2002-2003 and 2004-2005, respectively (Fig. 3b) . For both these periods, two strains (Cae03-6-31 and Cae05-40-27) were more prominent and represented greater than 50% of the isolated hMPV. During the winter of 2003-2004, the four subgroups were present equally (Fig. 3c) .

The present study reports on an epidemiological, clinical, and phylogenic study of hMPV outbreak during three consecutive winters in Normandy, France. Human metapneumovirus infection was assessed in nasal aspirates received in the laboratory during the months of November to March in order to perform a comparison of hMPV rate with the other common cold viruses. It was reported that, the hMPV infection season is longer than that of hRSV or influenza viruses [van den Hoogen et al., 2003; Esper et al., 2004; Mackay et al., 2006] . As such, the number of positive cases may have been underestimated although it was shown previously that hMPV circulates predominantly during the winterspring period [Bastien et al., 2003a; van den Hoogen et al., 2003; Esper et al., 2004; Williams et al., 2004; Pitoiset et al., 2010] . The number and the percentage of hMPV detected in virus-negative samples during the epidemics of 2002-2003, 2003-2004, and 2004-2005 were 59 (10%), 17 (22.6%), and 52 (8.8%), respectively. The incidence of hMPV infections are reported in many other studies Williams et al., 2004] . It is also suggested that the number of hMPV infections is variable from one winter to another [Falsey et al., 2003; Mackay et al., 2006; Williams et al., 2006] .

Human metapneumovirus was the fourth most prominent virus encountered during the present study and became the third most prominent if the results are focused only on children under 2 years old. Threequarters of the hMPV were detected on children under 2 years old and more than one third on those less than 6 months old. Human respiratory syncytial virus was the most predominant virus observed with an increase in prevalency in children under 2 years old. Influenza A virus was more commonly encountered in the population of over 2 years old. Unlike the hMPV and the other isolated viruses, the influenza A and B viruses show a distinct prevalence between children under 2 years old and those over 2 years old.

In comparison with hRSV and influenza viruses A and B, the outbreak of hMPV occurs over the same period and appears to be closely related to the hRSV epidemic profile. In fact, hMPV and hRSV outbreaks showed the same variation during the winter months [Boivin et al., 2003; Esper et al., 2004] . The influenza epidemic profile, which might be a more independent entity, did not correlate with hMPV or hRSV. During these three winters, the influenza A virus is more prominent than influenza B virus and the percentage was appreciatively the same for that described in Europe [Paget et al., 2003; Meijer et al., 2006] .

Only the negative samples were taken because a low percentage of co-infection was reported previously [Peiris et al., 2003; van den Hoogen et al., 2003; Mullins et al., 2004] and thus, even if this strategy introduced a selection bias, it allows a more reliable evaluation of the clinical impact of hMPV in affected children. Unlike upper respiratory tract infections, human metapneumovirus infection is normally associated with an acute LRI as observed in the present study and in other previous reports Freymuth et al., 2003; Williams et al., 2004] . Indeed, some patients present with only gastroenteritis in the absence of altered respiratory signs. No difference in severity of clinical manifestation was observed between the infection with group A or group B.

Cultivation of hMPV is very difficult to achieve and, at the start of the present study, no DFA test was available commercially. The RT-PCR was designed for the M gene to detect hMPV and to have an equivalent sensitivity for both genotypes. The M gene appeared to be suitable for hMPV detection, since it has already been validated [ Bellau-Pujol et al., 2005; Chano et al., 2005 ] and a genetic discrimination has been determined [van den Hoogen et al., 2001; Bastien et al., 2003b; Carr et al., 2005; Chano et al., 2005] . Moreover, the M gene seems to be a good candidate for the detection of other viruses like influenza viruses or coronaviruses [Donofrio et al., 1992; Vabret et al., 2001] . This PCR appears to be adapted to hMPV detection and characterization since the four hMPV clusters described were amplified successfully Mackay et al., 2004] .

In the present study, four hMPV subtypes termed, ''A1,'' ''A2,'' ''B1,'' and ''B2'' were detected by RT-PCR. The rate of the different subgroups was not the same for each winter and the strain composition varied during each outbreak. The strains circulating in France were described in four previous studies Bouscambert-Duchamp et al., 2005; Foulongne et al., 2006; Pitoiset et al., 2010] . The first two studies reported only A1 and A2 subtypes of hMPV, while Foulongne et al. [2006] detected four subtypes of the virus with an A1 subgroup in predominance. The last study [Pitoiset et al., 2010] reports on a phylogenic analysis carried out on 100 strains isolated over 8 years, from 2002 to 2009 with alternative subgroup predominance. During the first winter 2002-2003, the A group was the most prevalent, with 87.5% of A2 subgroup being the most prominent. This repartition was also found in many other parts of the world [Esper et al., 2004; Mackay et al., 2004; Ludewick et al., 2005; Pitoiset et al., 2010] . Two studies described another pattern of circulation of the virus during the winter of [2002] [2003] in which a prevalent B1 subgroup was noted [Gerna et al., 2005; Gray et al., 2006] . For the winter of [2003] [2004] , no predominant group could be identified. The incidence of groups A and B were approximately the same, although the number of strains belonging to group B and A were nine and eight, respectively. The B2 and A2 subgroup were the most frequently identified with five strains for each cluster. During this same period in Italy (winter of [2003] [2004] , the A2 subtype was also the most prominent strain, even if other strains were present Gerna et al., 2005] . The rate of occurrence of group B was confirmed during the winter of 2004-2005 and represented 96% of all hMPV detected, including one hMPV B2 strain. Over the three winter periods, 46.5% of B1, 37.5% of A2, 8% of A1, and 8% of B2 strains were observed. These data are closely related to those proposed by Mackay et al. [2006] in a study carried out in Australia during the same period. In addition, both predominant subgroups found in 2002 and 2004 were the same (viz., A2 and B1).

The RNA fragment studied is small (405 bp), but some mutations that were found in this gene allowed us to discriminate between, and compare, the different isolates. The sequence analysis for the three epidemic periods showed two distinct genetic profiles. The first profile, found during the winter of 2002-2003 and 2004-2005 , seemed to be composed of two predominant strains showing exactly the same sequence (Cae03-6-31 and Cae05-40-27) which represented at least 52% of hMPV detected, increasing to 79.5% in 2002-2003. During these two winters (2002-2003 and 2004-2005) , the foremost strains belonged to the A2 and B1 subgroups, respectively. This finding corroborates previous observations [Mackay et al., 2006] . For the winter period 2003-2004, the isolated strains belonged to the four subgroups and, at the most, two identical sequences were found. Even if it is impossible to consider these isolates like different strains because the analysis was assessed on too short a RNA fragment, such heterogeneity was not observed during the other winter periods and might be considered as the beginning of a shift from group A to group B as confirmed by the data of 2004. A shift of this nature was reported previously [Agapov et al., 2006; Mackay et al., 2006; Oliveira et al., 2009] Despite the fact that the 2003-2004 study is incomplete and that the period studied was only from November to March, the rate of hMPV can still be evaluated each year and the molecular evolution of the viruses can be analyzed. Human metapneumovirus is circulating in Normandy at a rate equivalent to that described previously Bastien et al., 2003a; Boivin et al., 2003; Esper et al., 2004] . The present clinical study demonstrates that hMPV is responsible for LRI in the majority of cases studied. All the reported hMPV subgroups are present in France but not with the same incidence each year. The results confirm the genetic evolution of hMPV in Normandy with a switch from the A group to the B group being noted during the winter of 2003-2004. Finally, molecular tools developed for the M gene appear to be useful for diagnosis and molecular analysis, since the four subgroups are detected and can be discriminated between in a significant manner.

Genetic variability of human metapneumovirus infection: Evidence of a shift in viral genotype without a change in illness

Prospective evaluation of rapid antigen tests for diagnosis of respiratory syncytial virus and human metapneumovirus infections

Human metapneumovirus infection in the Canadian population

Sequence analysis of the N, P, M and F genes of Canadian human metapneumovirus strains

Development of three multiplex RT-PCR assays for the detection of 12 respiratory RNA viruses

Virological features and clinical manifestations associated with human metapneumovirus: A new paramyxovirus responsible for acute respiratory-tract infections in all age groups

Global genetic diversity of human metapneumovirus fusion gene

Human metapneumovirus infections in hospitalized children

Detection of human metapneumovirus RNA sequences in nasopharyngeal aspirates of young French children with acute bronchiolitis by real-time reverse transcriptase PCR and phylogenetic analysis

Human metapneumovirusassociated respiratory tract infections in the Republic of Ireland during the influenza season of

Epidemiological survey of human metapneumovirus infection in a large pediatric tertiary care center

Human metapneumovirus and respiratory syncytial virus

Evolutionary dynamics of human and avian metapneumoviruses

Detection of influenza A and B in respiratory secretions with the polymerase chain reaction

Seroprevalence of human metapneumovirus in Japan

Detection of human metapneumovirus antigens in nasopharyngeal secretions by an immunofluorescent-antibody test

A 1-year experience with human metapneumovirus in children aged <5 years

Human metapneumovirus infections in young and elderly adults

Detection of human metapneumovirus in two Tunisian children

Human metapneumovirus infection in young children hospitalized with respiratory tract disease

Presence of the new human metapneumovirus in French children with bronchiolitis

Replication of respiratory viruses, particularly influenza virus, rhinovirus, and coronavirus in HuH7 hepatocarcinoma cell line

Comparison of multiplex PCR assays and conventional techniques for the diagnostic of respiratory virus infections in children admitted to hospital with an acute respiratory illness

Simultaneous detection and typing of human metapneumovirus strains in nasopharyngeal secretions and cell cultures by monoclonal antibodies

Changing circulation rate of human metapneumovirus strains and types among hospitalized paediatric patients during three consecutive winter-spring seasons

Multi-year study of human metapneumovirus infection at a large US Midwestern Medical Referral Center

Detection and characterisation of human metapneumovirus from children with acute respiratory symptoms in north-west England

Evaluation of real-time RT-PCR compared with conventional RT-PCR for detecting human metapneumovirus RNA from clinical specimens

Detection of human metapneumovirus in clinical samples by immunofluorescence staining of shell vial centrifugation cultures prepared from three different cell lines

Human metapneumovirus genetic variability

Molecular assays for detection of human metapneumovirus

Use of the P gene to genotype human metapneumovirus identifies 4 viral subtypes

Genetic diversity of human metapneumovirus over 4 consecutive years in Australia

Realtime reverse transcriptase PCR assay for detection of human metapneumoviruses from all known genetic lineages

Epidemiological and virological assessment of influenza activity in Europe, during the 2004-2005 winter

Human metapneumovirus: A not so new virus

Human metapneumovirus infection among children hospitalized with acute respiratory illness

Evidence of human metapneumovirus in Australian children

Epidemiology and genetic variability of human Metapneumovirus during a 4-year-long study in southeastern Brazil

Heterogeneous influenza activity across Europe during the winter of

Children with respiratory disease associated with metapneumovirus in Hong Kong

Respiratory tract reinfections by the new human Metapneumovirus in an immunocompromised child

Characterization of human metapneumoviruses isolated from patients in North America

Human metapneumovirus genotypes and severity of disease in young children (n¼100) during a 7-year study in Dijon hospital, France

Characterization of human metapneumovirus infections in Israel

Detection and pathogenicity of human metapneumovirus respiratory infection in pediatric Italian patients during a winter-spring season

Human metapneumovirus as a cause of community-acquired respiratory illness

Direct diagnosis of human respiratory coronaviruses 229E and OC43 by the polymerase chain reaction

Human metapneumovirus and lower respiratory tract disease in otherwise healthy infants and children

The role of human metapneumovirus in upper respiratory tract infections in children: A 20-year experience

Genetic diversity and evolution of human metapneumovirus fusion protein over twenty years