key: cord-0744840-7ehvcubd authors: Pitoiset, Cécile; Darniot, Magali; Huet, Frédéric; Aho, Serge‐Ludwig; Pothier, Pierre; Manoha, Catherine title: Human metapneumovirus genotypes and severity of disease in young children (n = 100) during a 7‐year study in Dijon hospital, France date: 2010-09-01 journal: J Med Virol DOI: 10.1002/jmv.21884 sha: e1895ec2cdc422394f8d20d542f598e7de461d8c doc_id: 744840 cord_uid: 7ehvcubd Human metapneumovirus is a cause of respiratory tract infections at all ages. Our objectives were to analyze the distribution of the A and B genotypes over 7 years in Dijon and to investigate a possible association between hMPV genotypes and disease severity. During 2002–2009, we genotyped the 100 isolates from children <3 years old with hMPV. Phylogenetic analysis indicated a change in the distribution of hMPV genotype over the years. Severity was then measured by detailed clinical evaluation. The hospitalization rate was greater when genotype B was involved 72.5% versus 53.3% (P = 0.054). Those infected with genotype B tended to have a higher clinical score, as measured by Vicente et al. 2006 (P = 0.07). We showed that, although clinical severity is not clearly associated with hMPV genotype in this study, pathological signs on chest X‐ray were observed more often in B subgroup (P < 0.01). J. Med. Virol. 82:1782–1789, 2010. © 2010 Wiley‐Liss, Inc. After being discovered in the Netherlands in 2001 [van den Hoogen et al., 2001] , human metapneumovirus (hMPV) was identified in 3-10% of hospitalized young children with acute respiratory tract infection in many countries worldwide as reviewed by Kahn [2006] . hMPV is a major pathogen that causes acute respiratory illness (ARI) in individuals of all ages. Higher morbidity is observed in young children, the elderly [Boivin et al., 2002; Falsey et al., 2003; and immunocompromised adults [Sumino et al., 2005; Williams et al., 2005; O'Gorman et al., 2006] . The clinical characteristics of hMPV infection in children range from mild upper respiratory tract to severe lower respiratory tract disease. The main clinical symptoms include rhinorrhea, cough, wheezing, and exacerbation of asthma. hMPV isolates are separated by phylogenetic analysis into two major genotypes (A and B), subdivided in four subgroups (A1, A2, B1, and B2) [Biacchesi et al., 2003; . Subgroup A2 has then been subdivided in two lineages A2a and A2b [Huck et al., 2006b] . Each season, hMPV circulation may vary, and hMPV genotypes circulate simultaneously in variable proportions [Williams et al., 2004; Gerna et al., 2005; Mackay et al., 2006] . In this study, we analyzed the distribution of the two hMPV genotypes and determined whether the clinical characteristics and severity of ARI episodes in children under 3 years of age correlated with a particular genotype. The Hospital of Dijon is of medium size with a capacity of 1,700 beds, 20,000 children are admitted to the emergency department every year. During [Manoha et al., 2008] . This assay does not include pre-screening culture methods. Respiratory syncytial virus (RSV), influenza (FLU), and parainfluenza (PI) were detected by IFA (Argène, Varilhes, France) throughout the study. RNA was extracted from the samples using TRIzol (Invitrogen, Cergy-Pontoise, France). To detect hMPV, two previously described sets of primers were used. From December 2002 to December 2004, the first set of primers amplified a part of the L gene (171 bp) [van den Hoogen et al., 2001] . From 2005 to 2009, PCR was based on a second primer set which amplified a part of the F gene (750 bp) . The analysis focused on all subjects under 3 years old who had positive results for hMPV from December 2002 to March 2009. However, those who were co-infected with other respiratory viruses (RSV, FLU, or PI) were excluded from the analysis. Records were also excluded when medical physicians suspected respiratory bacterial infections and positive results were found. The clinical characteristics of the 100 included patients were collected mainly retrospectively through reviews of patients' medical file by a reviewer blinded to the genotyping result. The severity of the illness was assessed by two different clinical scores, CSS1 and 2. As described by Vicente et al. [2006] , CSS1 ranged from 0 to 3 on the basis of the need for hospitalization, oxygen saturation <90%, and stay in an intensive care unit stay. As the clinical characteristics of hMPV infection resemble those of RSV, the second CSS was similar to that described by Martinello et al. [2002] for RSV and ranged from 0 to 6. Two points were given if the patient required mechanical ventilation during the illness, and one point was given for each of the following: hospital admission, hospitalization for >5 days, oxygen saturation <90%, and use of supplemental oxygen. PCR-amplified fragments were purified and directly sequenced in both directions using Big Dye Terminator (Applera Corporation, Foster City, CA) on an automatic DNA sequencer (Applera Corporation, Foster City, CA). Sequence comparisons were performed on a 102-nt segment of the L gene [van den Hoogen et al., 2001] and/or a 449-nt segment of the F gene . Bionumerics software (Applied Maths, Sint-Martens-Latem, Belgium) was used to determine the hMPV genotypes. Multiple-sequence alignments including some hMPV gene sequences available from GenBank were performed. Phylogenetic trees were generated by means of the UPGMA method. Bootstrap support was determined by 100 resamplings of the sequences. Results were analyzed using STATA version 8 (StataCorp, College Station, TX). Statistical differences in the distribution of hMPV genotypes over the years were determined by the w 2 -test. The clinical characteristics and severity of illness associated with genotype A and B were compared using one-way analysis' test or Kruskal-Wallis' test for continuous variables and the w 2 -test for categorical data, as appropriate. Multivariate logistic regression analysis was conducted to estimate independence of potential confounding factors on chest radiographs, that is, signs of respiratory distress, chronic respiratory disease, gastroesophageal reflux disease (GERD), and cardiopathy. 1A ). Among the 100 samples, 60% belonged to genotype A and 40% belonged to genotype B: of the 60A, seven were classified as A1 and 53 as A2; of the 40B, 12 were classified as B1 and 28 as B2 (Figs. 2 and 3). For the tree based on the L-gene fragment (Fig. 2) , which combined one third of the sequences, nucleotide identity between groups A and B was 81.4-82.4%, while it was 91.2-96.1% within group A. No B1 sequence was included in this tree, the nucleotide identity within subgroup B2 was 97.1-99.5%. For the tree based on the F-gene fragment (Fig. 3) , nucleotide identity between groups A and B was 83.2-85.1%. Within genotype A, nucleotide identity was 92.7-100%, within subgroup A1 (97.3-99.6%) and within subgroup A2 (96.4-100%); within genotype B, nucleotide identity was (95.0-95.8%), within subgroup B1 (98.9-100%), and within subgroup B2 (97.7-100%). Subgroup A was thus more divergent and the existence of two lineages within subtype A2, labelled A2a and A2b, as recently described by Huck et al. [2006b] , was confirmed. As the PCR on the L gene used at the beginning of the study (2002) (2003) (2004) did not allow us to distinguish lineages, only 27 A2 isolates were classified as either A2a or A2b. Five isolates (24%) belonged to A2a (97.9-99.4%), and 22 isolates (76%) belonged to A2b (98.0-99.6%). One strain isolated at the end of February 2008 (Dij-199FEB-08) did not fit into this dichotomy. The A genotype was predominant during the first two seasons of the study, then both genotypes A and B circulated together with alternative predominance of one or the other genotype (Fig. 1B) . There was not only a change in the predominant genotype but also a change within subgroups (Fig. 1C) . A2 was detected throughout the study. The A1 subgroup was detected only during the first two seasons of the study, whereas the B1 subgroup emerged in 2003-2004, peaked the next winter then remained at a low level from 2005-2006 to 2008-2009 . The B2 subgroup was detected mainly during the last three seasons of the study; its predominance was opposite to that of B1. Strains of genotypes A and B circulated simultaneously throughout the studied seasons although the relative distribution of each genotype changed over the years (w 2 -test, P ¼ 0.001). As lower sensitivity has been reported for the primer set (L6/L7) on the reference strains NL/17//00(A2) and NL/1/99 (B1) [Maertzdorf et al., 2004] , we may have missed some A2 or B1 isolates during the first two seasons. Even when we excluded this period from the analysis, the relative distribution of the genotypes (w 2 -test, P < 0.05) and of the subgroups (w 2 -test, P ¼ 0.001) changed significantly in different years. According to Genotype The sex and age were similar in both groups, with a median age of 5 months for children infected with hMPV genotype A and 6.5 months for children infected with hMPV genotype B (Table I) . When we compared the clinical characteristics associated with each genotype, we found a similar pattern. Major clinical symptoms for both groups included cough (A: 78%; B 70%), rhinorrhea (A: 53%; B 40%), fever (A: 61%; B: 67%), and wheezing The proportion of hospitalizations was higher for genotype B but without reaching statistical significance (A: 53%, B: 72.5%, w 2 -test, P ¼ 0.054) and the median duration of hospitalization was similar (A: 1.5 days; B: 2.5 days; Table I ). Pathological findings on chest X-ray were more frequent in genotype B (A: 15%; B 44%, P < 0.01), with most of the children showing bronchial wall thickening, alveolar condensation, perihilar haziness, or thoracic distension. There was no effect of confounding factors on chest radiographs as calculated by logistic regression analysis. We did not find any association between genotype and any critical factors of severity such as oxygen saturation <90%, intensive care unit stay or use of supplemental oxygen. Both clinical severity scores were slightly higher for genotype B, but the differences were not statistically significant (CSS1: P ¼ 0.07, CSS2: P ¼ 0.24). Of note, the one and only child who did not present any risk factors (cardiopathy, chronic respiratory disease, J. Med. Virol. DOI 10.1002/jmv etc.) but needed to be hospitalized in an intensive care unit had a B2 genotype. When patients with underlying diseases (cardiopathy, pulmonary disease and GERD) were removed from the analysis, the rate of hospitalization still tended to be higher for genotype B compared to genotype A (A: 48%; B 68%; P ¼ 0.089). When the four subgroups were compared, subgroup A1 was associated with the lowest rate and duration of hospitalization, percentage of LRTD as well as lowest values of CSS (Table II ). The differences, however, did not reach statistical significance. A1 appeared to induce a milder disease than did the other subgroups. When we restricted the analysis to patients admitted to hospital (n ¼ 61), clinical severity among genotypes was similar except for pharyngitis, which was more common in children with genotype B (P ¼ 0.041) than in those with genotype A. hMPV are an important cause of ARI in young children, and in immunocompromised or elderly adults. To date, studies on hMPV have mainly focused on clinical characterization of hMPV infection, but only few reports have studied the influence of genotype on the severity of infection. In this report, we characterized a large cohort of hMPV isolates from children over 7 consecutive years and detailed the clinical characteristics of hMPV infection in relation to genotype. Percentages are given relative to cases of genotype A (n ¼ 60) or B (n ¼ 40). For the following item, radiological signs in chest, only children who had a chest X-ray have been taken into account. Chronic respiratory disease was measured as binary variable and included bronchopulmonary dysplasia and one case of lung agenesis; GERD, gastroesophageal reflux disease. Asthma in the family was defined as a history of medical diagnosis of asthma in first-degree-relatives. Significant differences are in bold. Differences close to significance are also indicated. Detection over a long period showed that most of the infections occur between December and April; these results are similar to those showing a higher frequency in late winter and spring [van den Hoogen et al., 2001; Cane et al., 2003; Esper et al., 2004; Robinson et al., 2005; Agapov et al., 2006; Rafiefard et al., 2008] . We also report on the domination of the A genotype over the B genotype (60/40). These results are consistent with those of Rafiefard et al. [2008] who also identified hMPV genotypes over a long period (5 years) in Sweden. More precisely, we showed that genotype A predominated from 2002 to 2004, then both genotypes A and B cocirculated. In 2002-2003 a large number of specimens were positive for genotype A in our study as was the case in several countries. This result may be moderated because of the use of a primer set (originally published by Van den Hoogen 2001) located on L gene which has been shown to amplify hMPV strains of the different sublineages with variable sensitivities [Maertzdorf et al., 2004; Sarasini et al., 2006] . The assay on the F gene presented a higher sensitivity (100Â) than the L assay for the four subgroups. The gain in sensitivity obtained with the second pair over the first PCR was probably sufficient to allow the detection of metapneumovirus in children presenting to the hospital with respiratory symptoms and therefore relatively high viral loads, but some B strains may have been missed during the first two seasons. A change in the predominant genotype was observed over the years as well as a change in subgroups. Within genotypes, subgroup A2 and B2 accounted for most cases. A1 cases were detected in the early years, mainly 2001-2002 after which only a few A1 cases were reported worldwide in [Ludewick et al., 2005 Mackay et al., 2006; Huck et al., 2006b; Rafiefard et al., 2008] . A1 cases were no longer reported from 2004 to 2006 Matsuzaki et al., 2008] . Concurrently, subgroup B1 emerged to become the predominant subgroup in 2004-2005 in South Africa [Ludewick et al., 2005] , Australia , Japan [Matsuzaki et al., 2008] , and France as shown in this study. This underlines the fact that the prevalence of B1 was similar worldwide. Similar observations on changes in predominant circulating strains have led to suggestions that there is evasion of pre-existing immunity [Ludewick et al., 2005; Agapov et al., 2006] , which can be attributed to changes in immunity of the population in response to antigenic differences. Despite the stability of lineages over the years, a region of considerable variation (amino acids 260-300) on the fusion protein, analogous to an major antigenic region of the RSV fusion protein, has been identified [Yang et al., 2009] . This relative instability could be a way for the virus to escape the immune response. However, we cannot exclude the possibility that the dynamics of the virus population may vary from one genotype to another and from one year to another depending on environmental factors. Within lineage A2, two clusters, designated A2a and A2b, have been identified [Huck et al., 2006a,b] , with approximately 70% of the strains detected from 2002 to 2004 classified as A2b. Distinction between A2b and A2a strains was reported during one season in Croatia 2005 /2006 [Ljubin-Sternak et al., 2008 . In our study, RT-PCR analysis of the F fragment gene used from 2004 to 2009 allowed us to show that A2a and A2b were prevalent during this period and co-circulated. A majority of isolates belonged to subgroup A2b. At the end of the study, one isolate was on a separate branch from A2a and A2b. It may represent a new variant within genotype A and may signify an attempt of the virus to escape the current immune response of the population. Percentages are given relative to cases of A1, A2, B1, or B2 subgroups. For the following item, radiological signs in chest, only children who had a chest X-ray have been taken into account. Significative differences are in bold. Differences close to significance are also indicated. In the present study, we sought to determine whether the hMPV genotype was associated with disease severity. We did not find any major differences in the clinical pattern of hMPV infection in children <3 for each of the genotypes of hMPV suggesting that variations in viral genotype do not appear to be a major causal factor of differences in severity of the illness. Abnormalities on chest radiographs have been shown to occur frequently in children infected with hMPV [Esper et al., 2003; Williams et al., 2004; Wilkesmann et al., 2006] ; in the present study, we showed a higher frequency for genotype B. Three teams have screened a high number of isolates to look for correlations between genotype and disease severity in small children. None reported any significant differences between clinical data according to genotypes [Agapov et al., 2006; Sloots et al., 2006; Vicente et al., 2006 ]. However, a difference in the severity index as well as in the diagnosis of pneumonia was reported [Vicente et al., 2006] , suggesting that hMPV genotype A may be more pathogenic than hMPV genotype B. Although we studied the same population, that is children <3 years old, and applied the same severity index as described by Vicente, we did not find a correlation between genotype and severity score. The conclusion of these studies, which included a large number of strains was that there is no clear link between viral genotype and susceptibility to more severe hMPV infection. Our study was based on metapneumovirus, and most of the main causes of coinfections, namely RSV, Flu, and PI, were excluded. We did not test for human bocavirus, adenoviruses, coronavirus, and rhinoviruses. We probably missed some co-infections, but the clinical impact of viral coinfection involving hMPV and these viruses is unclear. Sloots et al. [2006] found no significant difference in disease severity between hMPV infected patients with a coinfection and those without. Only hMPV-RSV dual infections are thought to be more severe [Semple et al., 2005 ]. An association between hMPV infection and asthma has been reported by some authors who suggested that hMPV is a causative agent of acute wheezing in young children [Jartti et al., 2002; von Linstow et al., 2004; Manoha et al., 2007] . We did not find any difference in asthma between children infected with hMPV genotype A and children infected with genotype B. However, when we compared the four subgroups, asthma was only observed in patients with subgroup A2 and B2. Rhinoviruses are frequently associated with LRTD and asthmatic exacerbations. We found no significant difference regarding asthma according to genotype, but we cannot exclude the possibility that the presence of another virus, in particular rhinovirus, contributed to this clinical presentation. Recently, Matsuzaki et al. [2008] compared the clinical characteristics of A2, B1, and B2 hMPV infections for 85 outpatients under 15 years old over three seasons and showed that wheezing was statistically less frequent in A2 than in B1 and B2. They also reported a higher frequency of laryngitis in children with B1. Our study included samples collected over several seasons, which comprehensively characterize the clinical presentation of hMPV infection in children who present to hospital emergency departments. However, no clear relationship between one hMPV genotype and disease severity has been established. One possibility could be that the virulence is maintained by evolution in subgroup predominance, leading to progression from rather mild infections often due to the A1 subgroup to lower RTI mainly due to A2 and B2. However, other factors such as the viral load [Bosis et al., 2008] may be involved in disease severity. Investigation of a large panel of isolates as well as predicting factors of illness severity are required to gain insights into the molecular epidemiology of hMPV and viral pathogenesis. Genetic variability of human metapneumovirus infection: Evidence of a shift in viral genotype without a change in illness Genetic diversity between human metapneumovirus subgroups Virological features and clinical manifestations associated with human metapneumovirus: A new paramyxovirus responsible for acute respiratory-tract infections in all age groups Association between high nasopharyngeal viral load and disease severity in children with human metapneumovirus infection Human metapneumovirus in a haematopoietic stem cell transplant recipient with fatal lower respiratory tract disease Human metapneumovirus infection in the United States: Clinical manifestations associated with a newly emerging respiratory infection in children A 1-year experience with human metapneumovirus in children aged <5 years Human metapneumovirus infections in young and elderly adults Changing circulation rate of human metapneumovirus strains and types among hospitalized pediatric patients during three consecutive winter-spring seasons Human metapneumovirus infection in a hematopoietic stem cell transplant recipient with relapsed multiple myeloma and rapidly progressing lung cancer Novel human metapneumovirus sublineage Metapneumovirus and acute wheezing in children Epidemiology of human metapneumovirus Detection of genetic lineages of human metapneumovirus in Croatia during the winter season Human metapneumovirus genetic variability Genetic diversity of human metapneumovirus over 4 consecutive years in Australia Real-time reverse transcriptase PCR assay for detection of human metapneumoviruses from all known genetic lineages Epidemiological and clinical features of hMPV, RSV and RVs infections in young children Rapid and sensitive detection of metapneumovirus in clinical specimens by indirect fluorescence assay using a monoclonal antibody Correlation between respiratory syncytial virus genotype and severity of illness Clinical impact of human metapneumovirus genotypes and genotype-specific seroprevalence in Yamagata Human metapneumovirus in adults: A short case series Epidemiologic characteristics and seasonal distribution of human metapneumovirus infections in five epidemic seasons in Seasonality and clinical features of human metapneumovirus infection in children in Northern Alberta Detection and pathogenicity of human metapneumovirus respiratory infection in pediatric Italian patients during a winterspring season Dual infection of infants by human metapneumovirus and human respiratory syncytial virus is strongly associated with severe bronchiolitis Detection of severe human metapneumovirus infection by real-time polymerase chain reaction and histopathological assessment Differences in clinical severity between genotype A and genotype B human metapneumovirus infection in children Human metapneumovirus and respiratory syncytial virus in hospitalized Danish children with acute respiratory tract infection Human metapneumovirus infections cause similar symptoms and clinical severity as respiratory syncytial virus infections Human metapneumovirus and lower respiratory tract disease in otherwise healthy infants and children A prospective study comparing human metapneumovirus with other respiratory viruses in adults with hematologic malignancies and respiratory tract infections Genetic diversity and evolution of human metapneumovirus fusion protein over twenty years This work was supported by the Centre Hospitalier Universitaire of Dijon. We thank Jean-Baptiste Bour for helpful comments, Emmanuel Florentin for statistics on virology laboratory data and Philip Bastable for editorial assistance.