# Genetic variability of human metapneumo- and bocaviruses in children with respiratory tract infections

Vasiliki Pogka, Afroditi Moutousi, Athanasios Kossyvakis, Antonios Kalliaropoulos, Dionyssios N Sgouras, Maria Giannaki, Andreas Mentis

► To cite this version:

Vasiliki Pogka, Afroditi Moutousi, Athanasios Kossyvakis, Antonios Kalliaropoulos, Dionyssios N Sgouras, et al.. Genetic variability of human metapneumo- and bocaviruses in children with respiratory tract infections. Influenza and Other Respiratory Viruses, 2014, 8 (1), pp.107-115. 10.1111/irv.12185 . pasteur-02195905

HAL Id: pasteur-02195905

https://riip.hal.science/pasteur-02195905

Submitted on 26 Jul 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Genetic variability of human metapneumo- and bocaviruses in children with respiratory tract infections

Vasiliki Pogka,ª Afroditi Moutousi,ª Athanasios Kossyvakis,ª Antonios Kalliaropoulos,ª
Dionyssios N. Sgouras,ᵇ Maria Giannaki,ᶜ Andreas F. Mentisª,ᵇ

*Correspondence*: Andreas F. Mentis, Hellenic Pasteur Institute, Vas. Sofias 127, 11521 Athens, Greece.

E-mail: mentis@pasteur.gr

Accepted 10 August 2013. Published Online 7 November 2013.

**Objectives** The genotypic analysis of human metapneumo-
(HMPV) and boca-(HBoV) viruses circulating in Greece and their
comparison to reference and other clinical strains.

**Design** Genetic analysis of representative strains over three consecutive winter seasons of the years 2005–2008.

**Setting** Representative positive specimens for HMPV and HBoV from paediatric patients of healthcare units and hospitals in Southern Greece with influenza-like illness or other respiratory tract infections.

**Sample** Seven to ten positive specimens for either HMPV or HBoV from each winter period. In total, 24 specimens positive for HMPV and 26 for HBoV, respectively.

**Main outcome measures** Sequence diversity of HMPV and HBoV strains by sequencing the complete G and VP1/VP2 genes, respectively.

**Results** In total, 24 HMPV strains were found to have a 92–100% nucleotide and a 85.9–100% amino acid identity. Phylogenetic

analysis based on the number of amino acid differences, revealed
circulation of 4 different subclusters belonging to genetic lineage B2.
Similarly, analysis of 26 HBoV strains indicated that 22 clustered
within genotype St2, 2 into genotype St1 and the remaining 2
formed a third cluster derived from potential recombination
between different St1 genotype strains. St2 HBoV genotype was
observed throughout the whole observation period whereas St1 only
during the second and the third winter period. Higher levels of
heterogeneity were observed between HMPV compared to HBoV
strains.

**Conclusions** Phylogenetic analysis revealed circulation of one single lineage (B2) for HMPV viruses and predominance of St2 genotype for HBoV viruses. A possible recombination between St1 genotype strains of HBoV was observed.

**Keywords** Bocavirus, children, metapneumovirus, phylogenetic analysis, respiratory viruses.

*Please cite this paper as:* Pogka *et al.* (2014) Genetic variability of human metapneumo- and bocaviruses in children with respiratory tract infections. Influenza and Other Respiratory Viruses 8(1), 107–115.

Introduction

Newly discovered viruses, such as human metapneumovirus (HMPV) and human bocavirus (HBoV) are important pathogens causing respiratory tract infection (RTI) in susceptible populations, particularly in children and the elderly.<sup>1–4</sup>

HMPV is a negative single stranded RNA virus belonging to the family of *Paramyxoviridae*, which was first isolated in 2001<sup>5</sup> from nasopharyngeal aspirates. It is responsible for about 5–15% of the worldwide respiratory tract infections, affecting both young children and adults, causing symptoms ranging from mild disease of the upper respiratory tract system to severe bronchiolitis and pneumonia.<sup>2</sup> In infants, HMPV incidence is even higher, reaching 25%.<sup>6</sup> Genetic

analysis of HMPV isolates revealed two ‘major’ lineages (A
and B) and four ‘minor’ lineages (A1, A2, B1 and B2), based
on the sequence variability of the attachment (G) and fusion
(F) surface glycoproteins. The existence of two further
sublineages, A2a and A2b, has further been suggested.<sup>7</sup>

HBoV, a single stranded DNA virus belonging to the
family of *Parvoviridae*, was discovered in 2005⁸ and has been
detected in nasopharyngeal aspirates, sera and blood samples
of patients with respiratory tract infection and in faecal
specimens of patients with acute respiratory illness and/or
gastroenteritis. HBoV affects, mostly, children younger than
2 years old, with its incidence ranging from 2·7 to 19% in
respiratory and from 0·8 to 9·1% in faecal samples, obtained
from patients with respiratory infection and gastroenteritis,
respectively.⁹ Four species of HBoV have been recently

<footnote>ªNational Influenza Reference Laboratory of Southern Greece, Hellenic Pasteur Institute, Athens, Greece. ᵇLaboratory of Medical Microbiology, Hellenic Pasteur Institute, Athens, Greece. ᶜDepartment of Microbiology, Aghia Sophia Children's Hospital, Athens, Greece.</footnote>
discovered with HBoV1 being predominantly a respiratory pathogen, whereas HBoV2, HBoV3 and HBoV4 being found mainly in stool.¹⁰ Nucleotide sequence analysis of the two capsid proteins (VP1 and VP2), which show high variability compared with coding sequences of the two viral nonstructural proteins (NS1 and NP1), has further divided HBoV1 into genotypes St1 and St2. These two major HBoV genotypes correspond to the original St1 (Stockholm 1) and St2 (Stockholm 2) isolates.⁸

Epidemiological and phylogenetic analysis of HMPV and HBoV has been carried out extensively in other countries but not in Greece, where only a few studies concerning newly recognized respiratory viruses have been published. Using molecular methods for identification, we recently reported the incidence of 13 respiratory pathogens, including HMPV and HBoV, in children presenting with influenza-like illness (ILI).¹¹ In this study, we aimed to determine the genetic diversity of the HMPV and HBoV clinical strains identified within the Greek population by comparison with reference, as well as, wild strains described in other countries.

# Patients and methods

## Clinical specimens

Clinical specimens including nasopharyngeal aspirates and/or throat swabs were collected over the winter seasons (November to May) of 2005/2006, 2006/2007 and 2007/2008 by paediatricians in healthcare units of Southern Greece and in collaborating paediatric hospitals of Athens. Specimens were sent to the National Influenza Reference Laboratory of Southern Greece in viral transport medium (Mediaproducts BV, Groningen, The Netherlands). Our sample group consisted in total of 3306 specimens from paediatric patients aged 0–18 years-old, 1272 out of whom suffered from ILI (sudden onset of symptoms, with high fever and cough in the

absence of other diagnosis) and the remaining 2034 from other upper and lower respiratory tract infections which did not fall within the definition of ILI. The presence of HMPV RNA and HBoV DNA concerning patients with ILI has been described in our previous study.¹¹ In the present study, the same PCR protocols were used to reveal the presence of HMPV and HBoV in patients with other upper and lower respiratory tract infections (data not shown). Moreover, all specimens found negative for HBoV1 in our previous study¹¹ were further screened for the three recently identified species, HBoV2, HBoV3 and HBoV4, as previously described.¹²

## Nucleotide sequence analysis of HMPV and HBoV strains

In total, presence of HMPV RNA was detected in 188 (5·7%) and HBoV DNA was detected in 193 (5·8%) cases of the 3306 patients screened. Seven to ten positive specimens for either HMPV or HBoV belonging to each winter period (in total 50 samples) were further studied using specific primers mapping highly variable regions of the viral genome. The selected specimens amounting to 13% of total positive strains corresponded to the beginning (November-December), middle (January-February) and end (March-May) of each winter season. More specifically, 24 samples positive for HMPV were further amplified to obtain the genomic region corresponding to the 969 bp DNA fragment, encompassing the full-length of G protein coding region (240 amino acids).¹³ Likewise, 26 samples positive for HBoV were amplified utilizing 4 overlapping sets of primers, spanning the whole length (672 amino acids) of the coding region corresponding to the VP1 and VP2 proteins (Table 1).¹⁴

PCR products were purified using the QIAquick PCR purification kit (Qiagen, Hilden, Germany) and the Min-EluteTM gel extraction kit (Qiagen) and subjected to direct sequencing in both directions utilizing the GenomeLab

<table><caption>**Table 1.** Primers sequences for sequencing analysis used in this study</caption><thead><tr><td>Primer</td><td>Sequence (5'-3')</td><td>Target</td><td>Position (5'-3')</td><td>Reference citation</td></tr></thead><tbody><tr><td>HMPV F1</td><td>TACAAAACAAGAACATGGACAAG</td><td>G</td><td>6183-6206*</td><td>13</td></tr><tr><td>HMPV R1</td><td>GAGATAGACATTAACAGTGGATT</td><td>G</td><td>7127-7149*</td><td>13</td></tr><tr><td>HBoV F1</td><td>GATAACTGACGAGGAAATGCT</td><td>VP1/VP2</td><td>3009-3029**</td><td>14</td></tr><tr><td>HBoV R1</td><td>AGTATGTCCATGGAGTTGTGA</td><td>VP1/VP2</td><td>3711-3731**</td><td>14</td></tr><tr><td>HBoV F2</td><td>TTCAGAATGGTCACCTCTACA</td><td>VP1/VP2</td><td>3639-3659**</td><td>14</td></tr><tr><td>HBoV R2</td><td>CTGTGCTTCCGTTTTGTCTTA</td><td>VP1/VP2</td><td>4266-4286**</td><td>14</td></tr><tr><td>HBoV F3</td><td>AACTTTGACTGTGAATGGGTTA</td><td>VP1/VP2</td><td>4172-4193**</td><td>14</td></tr><tr><td>HBoV R3</td><td>AAATAGTGCCCTGGAGGATGAT</td><td>VP1/VP2</td><td>4767-4787**</td><td>14</td></tr><tr><td>HBoV F4</td><td>ACCAAGGGCTGACAAACACA</td><td>VP1/VP2</td><td>4711-4730**</td><td>14</td></tr><tr><td>HBoV R4</td><td>TGTACAACAACAACACATTAAAAG</td><td>VP1/VP2</td><td>5276-5299**</td><td>14</td></tr></tbody></table>

*Numbering according to reference FJ168778.

**Numbering according to reference DQ000496.
DTCS-quick start sequencing kit (Beckman Coulter, Brea, CA, USA) on a CEQTM 8000 genetic analyser (Beckman Coulter).

## Phylogenetic analysis

Sequences obtained were examined in terms of closest homology sequence using BLAST software (http://www.ncbi.nlm.nih.gov/BLAST/). Multiple sequence alignments of representative strains identified in Europe, in addition to GenBank sequence data from original HMPV and HBoV prototype strains were made by BioEdit sequence alignment editor. Phylogenetic analysis of HMPV and HBoV-positive samples was performed using the molecular evolutionary genetics analysis (MEGA) software, version 5.15 Phylogenetic trees were constructed by the neighbour-joining method (NJ, tree algorithm inferred with the Kimura 2-parameter substitution model of sequence evolution)16,17 while, a bootstrap resampling analysis was performed (1.000 replicates) to test tree robustness.18 Moreover, comparison of the mean sequence diversity in pairwise alignments was applied on both our strains and the representative prototypes. More specifically, the aligned sequences were compared using uncorrected p-distances (proportion of observed differences), handling gaps with complete deletion.

## Nucleotide sequences accession numbers

All 50 sequences determined in the present report have been deposited in the GenBank/EMBL/DDBJ database under accession numbers JQ513461-JQ513510. The reference strains used in this study have been previously assigned with the following accession numbers: a) for the Netherlands HMPV prototypes AY296040 (NL/94/01), AY296034 (NL/1/99), AF371337 (NL/1/00), AY296021 (NL/00/17) b) for Stockholm HBoV prototypes DQ000495 (St1), DQ000496 (St2).

# Results

## Phylogenetic analysis of HMPV sequences

Twenty-four HMPV-positive samples, including 9 isolates detected in 2006, 8 detected in 2007 and 7 isolates detected in 2008, were characterized at the molecular level and a phylogenetic tree of the complete G gene was generated. Comparison with the prototype strains from the Netherlands13 revealed that all Greek-HMPV isolates clustered with NL/94/01 in genetic lineage B2 (Figure 1) and presented a nucleotide (nt) identity of 92.0–100.0% and an amino acid (aa) identity of 85.9–100.0% to each other.

Phylogenetic analysis of the studied strains classified them within 4 separate groups based on the number of aa differences they shared. More specifically, group 1 included HMPV strains that circulated during the first winter period and shared an aa substitution at position 238 (P238S), group

2 included strains which circulated during the second winter season and shared an aa substitution at position 200 (N200K) and group 3 included strains that circulated during the whole period under investigation and revealed 5 aa substitutions (D84Y, L105P, T166P, E190K and S198P). Concerning group 4, sequence alignment revealed the insertion of 2 additional amino acids in specimen 1106/08 (one Lys and one Glu) and 4 amino acids in specimen 181/07 (two Lys and two Glu). These two HMPV strains also shared a number of amino acid substitutions in relation to all the other Greek specimens analysed in this study (Figure 2).

Intergroup analysis exposed 4 aa substitutions between groups 2 and 3 (I22L, E63D, D208N and A213T), 7 aa substitutions between groups 1, 2 and 3 (T39S, L56S, Q93R, S109L, T149A, G155R and P220S) and 1 aa substitution between groups 1, 3 and 4 (D111G/N). Comparing the Greek strains with the HMPV B2 reference sequence, they all shared 8 aa substitutions throughout the whole G protein (M55T/K, A98T, A99T, Y133H/R, Y187N/K, V188G/R, S201N/G, H221Q/Y).

To summarize, sequencing analysis suggested that among HMPV strains, only those belonging to group 3 circulated throughout the whole period under study. Intragroup analysis within group 3 also suggested specific variations within the strains that circulated during each winter season. More specifically, all 2007 strains revealed a proline to threonine aa substitution at position 73 (P73T) and a proline to serine aa substitution at position 121 (P121S). The P121S substitution was also observed in all strains from 2006 and in 3 strains from 2008. Within each winter period, no sequence variations were observed among strains in each group.

## Phylogenetic analysis of HBoV sequences

Twenty-six representative HBoV-positive samples were analysed using the complete VP1/VP2-gene region. Eight of these samples were detected in 2006, 8 in 2007 and the remaining 10 samples in 2008. Phylogenetic analysis of the Greek-HBoV strains indicated that 2 (7.7%) of them clustered into genotype St1 and 22 (84.6%) into genotype St2 (Figure 3). The remaining 2 (7.7%) strains studied, HBoVs/Athens.GRC/35.07 and HBoVs/Athens.GRC/108.07, which clustered neither with St1 nor with St2 genotype, formed a third group which shared 98.2–98.5% nt identity with St1 and 98.4–98.7% nt identity with St2. All of the Greek strains showed an nt identity of 98.2–100.0% and a deduced aa identity of 95.6–100.0% with each other and appeared to be most closely related to the Tunisian (TUN2207) and the Chinese (CQ201102) strains.

More specifically, for specimens belonging to St1 and St2 groups, analysis of intra- and intergenotype nt and aa sequences revealed different levels of homology. When comparing strains belonging to the above genotypes (19 different strains), diversity at the nt level varied between 0.05
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Evolutionary relationships of Greek HMPV strains with reference genotypes based on the complete coding sequence of G gene. The evolutionary history was inferred using the Neighbor-Joining method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. Only values above 60% are shown. Strains isolated in 2006 (▲), 2007 (●) and in 2008 (■) are depicted. Reference strains are indicated in italics.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 3.** Evolutionary relationships of the Greek HBoV strains with reference genotypes using the complete coding sequence of VP1/VP2 gene. Evolutionary history was inferred using the Neighbour-Joining method. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. Only values above 60% are shown. Strains isolated in 2006 (▲), 2007 (●) and in 2008 (■) are depicted. Reference strains are indicated in italics.</figcaption></figure>
and 0·99% and at the aa level between 0·15 and 0·45%. When comparing only the two St1 strains, diversity at the nt level was 0·55% and at the aa level was 0·15%. Likewise, for St2 strains (17 different strains), the highest divergence was up to 0·55% at the nt and 0·3% at the aa level.

Interestingly, 5 aa substitutions were observed in both strains belonging to the third HBoV group and included L40S, A340P, G415S, T521N and F540Y while three additional mutations (P430L, L433F and P541L) were revealed in the strain HBoVs/Athens.GRC/108.07. The boot scanning analysis suggested that the strains HBoVs/Athens.GRC/35.07 and HBoVs/Athens.GRC/108.07 resulted from recombination between a HBoVs/Athens.GRC/27.07-like and a HBoVs/Athens.GRC/24.08-like strain (belonging to St1 group), with a breakpoint located at position 1260 bp of the complete VP1/VP2 region (Figure 4). Bidirectional sequencing, as well as, the fact that the recombination breakpoint was located within the 3rd amplicon of VP1/VP2 region suggested that this result was not due to mixed infection.

In conclusion, sequencing data suggested that at least 3 groups of HBoV cocirculated in Southern Greece during the study period. Similar observations were also made by Chieochansin et al.,¹⁴ Lin et al¹⁹ and Xu et al.²⁰ In each of the above study populations, including ours, the representative strain of group 1 and group 2 was the prototype strains St1 and St2, respectively. The Greek strains, which clustered within group 1, were found to share a common aa substitution at position 590 (S590T, with respect to VP1 numbering) with all the St1 reference strains used in the phylogenetic analysis, compared with group 2 (Figure 3), an observation also made by others.¹⁴,¹⁹ Group 3 included 4 strains which came from Taiwan,¹⁹ Thailand¹⁴ and Greece. Additional substitutions at positions 415 and 540 (G415S

and F540Y) found in the VP1 protein of Greek strains HBoVs/Athens.GRC/35.07 and HBoVs/Athens.GRC/108.07 had also been reported in two Italian strains.²¹

Discussion

This is the first study on the genetic diversity of HMPV and HBoV strains in children with respiratory tract infections in Greece, including paediatric patients who suffered from ILI¹¹ as well as other respiratory infections.

Since its initial documentation in the Netherlands in 2001, HMPV has been reported to circulate worldwide.²²⁻²⁴ We only found the B2 HMPV lineage circulating among the Greek population during the winter season of the three years examined, although cocirculation of all four viral lineages during the same period is not unusual.²⁵,²⁶ Reports from Croatia, Ireland and France have also indicated predominant circulation of lineage B2, during the winter periods of 2005–2006, 2006–2007 and the years 2006–2008, respectively.²⁷⁻²⁹ Gaunt et al³⁰ suggested a potential back-shift to lineage B2 in Scotland during the winter season of 2007-2008. The observation of only one circulating lineage among the study population was also mentioned in other studies. In Spain,³¹ Italy³² and China³³,³⁴ circulation of HMPV A2 lineage was reported. However, a number of studies have indicated that circulating predominant strains may vary on an annual basis and that their replacement may occur, on average, every 1–3 years.³⁵,³⁶

As year-round surveillance was not performed and summer infections from HMPV and HBoV cannot be excluded, recirculation of the same HMPV lineage over three consecutive years within the Greek population could be the result of the maintenance of the virus throughout the year. Alternatively, HMPV lineage B2 could have been

<figure><img src="image_4.png" /><figcaption>**Figure 4.** Applied bootscan in Simplot software to detect the recombination point between hBoV strains. The most likely breakpoint was located at position 1260 bp of the complete VP1/VP2 gene between a HBoVs/ Athens.GRC/27.07-like and a HBoVs/ Athens.GRC/24.08-like strain.</figcaption></figure>
reintroduced each year in the Greek population with minor intralineage alterations causing infections to the naïve population.

Phylogenetic analysis of the HBoV full-length VP1 and VP2 proteins was performed on 26 HBoV-positive samples revealing limited sequence variations among them as well as in relation to available sequences of previously described HBoV strains. Genotype St2 was the only HBoV genotype, which circulated during the first winter period of 2005–2006, while both St1 and St2 cocirculated during the second and the third study period, with predominance of genotype St2. Interestingly, the majority of HBoV strains, detected in Italian children during winter periods of the years 2004–2008, also belonged to genotype St2 (62.5%).²¹ An even higher percentage (75.4%) of HBoV strains belonging to genotype St2 during the years 2002–2005 was reported from Germany.³⁷

Our study exposed the circulation of a third group of HBoV strains during the second surveillance period. Lin et al¹⁹ and Xu et al²⁰ also reported the circulation of HBoV strains belonging to a third group different from St1 and St2, which resulted from recombination events between different circulating viral strains. Such genetic recombination events have been reported in parvoviruses³⁸,³⁹ and in different HBoV species.⁴⁰⁻⁴² Our data suggest that the strains HBoVs/ Athens.GRC/35.07 and HBoVs/Athens.GRC/108.07 resulted from recombination events between a HBoVs/Athens.GRC/ 27.07-like and a HBoVs/Athens.GRC/24.08-like strains, with a breakpoint site in the VP1/VP2 region. To date, similar intraspecies recombination events involving HBoV1 have been reported before, with a breakpoint located in the NS1,¹⁹ VP1/VP2²¹ and NS1 and VP1/VP2²⁰ regions. Taken together, these data suggest that intraspecies recombination could occur in HBoV1 at different positions of its genome. Further studies, including whole genome sequencing, are necessary to identify breakpoint patterns in these viruses, thus explor- ing their potential association with viral spread among populations.

Higher similarity observed between HBoV sequences to that observed between HMPV strains may be explained by the widespread assumption that mutation rates correspond to the fidelity of replication. DNA viruses, such as HBoV, mutate slower compared with RNA viruses, such as HMPV, because DNA polymerases possess proof-reading activity, which further reduces mutation rates during DNA replica- tion.⁴³ Moreover, the rate of generation and correction of genomic mutations may depend on a number of factors, including generation time, genomic architecture, replication speed, transmission and environmental effects.

Our study highlighted the heterogeneity of viral strains circulating within the Greek population by determination of the complete sequences of HMPV G gene and HBoV VP1/ VP2 gene. The phylogenetic analysis performed aimed to

monitor the genetic evolution of HMPV and HBoV during three consecutive years and revealed an individual predom- inant lineage for the former and a combination of genotypes for the latter.

## References

1 Falsey AR, Erdman D, Anderson LJ, Walsh EE. Human metapneumo- virus infections in young and elderly adults. J Infect Dis 2003; 187:785-790.

2 Kahn JS. Epidemiology of human metapneumovirus. Clin Microbiol Rev 2006; 19:546-557.

3 Weissbrich B, Neske F, Schubert J et al. Frequent detection of bocavirus DNA in German children with respiratory tract infections. BMC Infect Dis 2006; 6:109.

4 Maggi F, Andreoli E, Pifferi M, Meschi S, Rocchi J, Bendinelli M. Human bocavirus in Italian patients with respiratory diseases. J Clin Virol 2007; 38:321-325.

5 van den Hoogen BG, de Jong JC, Groen J et al. A newly discovered human pneumovirus isolated from young children with respiratory tract disease. Nat Med 2001; 7:719-724.

6 Maggi F, Pifferi M, Vatteroni M et al. Human metapneumovirus associated with respiratory tract infections in a 3-year study of nasal swabs from infants in Italy. J Clin Microbiol 2003; 41:2987-2991.

7 Feuillet F, Lina B, Rosa-Calatrava M, Boivin G. Ten years of human metapneumovirus research. J Clin Virol 2012; 53:97-105.

8 Allander T, Tammi MT, Eriksson M, Bjerkner A, Tiveljung-Lindell A, Andersson B. Cloning of a human parvovirus by molecular screening of respiratory tract samples. Proc Natl Acad Sci U S A 2005; 102:12891-12896.

9 Lindner J, Modrow S. Human bocavirus-a novel parvovirus to infect humans. Intervirology 2008; 51:116-122.

10 Jartti T, Hedman K, Jartti L, Ruuskanen O, Allander T, Soderlund- Venermo M. Human bocavirus-the first 5 years. Rev Med Virol 2011; 22:46-64.

11 Pogka V, Kossivakis A, Kalliaropoulos A et al. Respiratory viruses involved in influenza-like illness in a Greek pediatric population during the winter period of the years 2005-2008. J Med Virol 2011; 83:1841-1848.

12 Kantola K, Sadeghi M, Antikainen J et al. Real-time quantitative PCR detection of four human bocaviruses. J Clin Microbiol 2010; 48:4044-4050.

13 van den Hoogen BG, Herfst S, Sprong L et al. Antigenic and genetic variability of human metapneumoviruses. Emerg Infect Dis 2004; 10:658-666.

14 Chieochansin T, Chutinimitkul S, Payungporn S et al. Complete coding sequences and phylogenetic analysis of Human Bocavirus (HBoV). Virus Res 2007; 129:54-57.

15 Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGAS: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 2011; 28:2731-2739.

16 Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 1987; 4:406-425.

17 Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980; 16:111-120.

18 Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution 1985; 39:783-791.

19 Lin JH, Chiu SC, Lin YC et al. Clinical and genetic analysis of Human Bocavirus in children with lower respiratory tract infection in Taiwan. J Clin Virol 2009; 44:219-224.
20 Xu L, He X, Zhang DM et al. Surveillance and genome analysis of human bocavirus in patients with respiratory infection in guangzhou, china. PLoS One 2012; 7:e44876.

21 Zappa A, Canuti M, Frati E et al. Co-circulation of genetically distinct human metapneumovirus and human bocavirus strains in young children with respiratory tract infections in Italy. J Med Virol 2011; 83:156-164.

22 Arnott A, Vong S, Sek M et al. Genetic variability of human metapneumovirus amongst all ages population in Cambodia between 2007 and 2009. Infect Genet Evol 2011; 15:43-52.

23 Escobar C, Luchsinger V, de Oliveira DB, Durigon E, Chnaiderman J, Avendano LF. Genetic variability of human metapneumovirus isolated from Chilean children, 2003-2004. J Med Virol 2009; 81:340-344.

24 Lamson DM, Griesemer S, Fuschino M, St George K. Phylogenetic analysis of human metapneumovirus from New York State patients during February through April 2010. J Clin Virol 2012; 53:256-258.

25 Gerna G, Sarasini A, Percivalle E et al. Prospective study of human metapneumovirus infection: diagnosis, typing and virus quantification in nasopharyngeal secretions from pediatric patients. J Clin Virol 2007; 40:236-240.

26 Mackay IM, Bialasiewicz S, Jacob KC et al. Genetic diversity of human metapneumovirus over 4 consecutive years in Australia. J Infect Dis 2006; 193:1630-1633.

27 Ljubin-Sternak S, Santak M, Cepin-Bogovic J et al. Detection of genetic lineages of human metapneumovirus in Croatia during the winter season 2005/2006. J Med Virol 2008; 80:1282-1287.

28 Pitoiset C, Darniot M, Huet F, Aho SL, Pothier P, Manoha C. Human metapneumovirus genotypes and severity of disease in young children (n = 100) during a 7-year study in Dijon hospital, France. J Med Virol 2010; 82:1782-1789.

29 Carr MJ, Waters A, Fenwick F, Toms GL, Hall WW, O'Kelly E. Molecular epidemiology of human metapneumovirus in Ireland. J Med Virol 2008; 80:510-516.

30 Gaunt E, McWilliam-Leitch EC, Templeton K, Simmonds P. Incidence, molecular epidemiology and clinical presentations of human metapneumovirus; assessment of its importance as a diagnostic screening target. J Clin Virol 2009; 46:318-324.

31 Camps M, Ricart S, Dimova V et al. Prevalence of human metapneumovirus among hospitalized children younger than 1 year in Catalonia, Spain. J Med Virol 2008; 80:1452-1460.

32 Guido M, Quattrocchi M, Campa A et al. Human metapneumovirus and human bocavirus associated with respiratory infection in Apulian population. Virology 2011; 417:64-70.

33 Mao HW, Yang XQ, Zhao XD. Characterization of human metapneumoviruses isolated in Chongqing, China. Chin Med J (Engl) 2008; 121:2254-2257.

34 Zhang C, Du LN, Zhang ZY et al. Detection and genetic diversity of human metapneumovirus in hospitalized children with acute respiratory infections in Southwest China. J Clin Microbiol 2012; 50:2714-2719.

35 Schildgen V, van den Hoogen B, Fouchier R et al. Human Metapneumovirus: lessons learned over the first decade. Clin Microbiol Rev 2011; 24:734-754.

36 Aberle JH, Aberle SW, Redlberger-Fritz M, Sandhofer MJ, Popow-Kraupp T. Human metapneumovirus subgroup changes and seasonality during epidemics. Pediatr Infect Dis J 2010; 29:1016-1018.

37 Neske F, Blessing K, Tollmann F et al. Real-time PCR for diagnosis of human bocavirus infections and phylogenetic analysis. J Clin Microbiol 2007; 45:2116-2122.

38 Shackelton LA, Hoelzer K, Parrish CR, Holmes EC. Comparative analysis reveals frequent recombination in the parvoviruses. J Gen Virol 2007; 88:3294-3301.

39 Ohshima T, Mochizuki M. Evidence for recombination between feline panleukopenia virus and canine parvovirus type 2. J Vet Med Sci 2009; 71:403-408.

40 Arthur JL, Higgins GD, Davidson GP, Givney RC, Ratcliff RM. A novel bocavirus associated with acute gastroenteritis in Australian children. PLoS Pathog 2009; 5:e1000391.

41 Kapoor A, Slikas E, Simmonds P et al. A newly identified bocavirus species in human stool. J Infect Dis 2009; 199:196-200.

42 Kapoor A, Simmonds P, Slikas E et al. Human bocaviruses are highly diverse, dispersed, recombination prone, and prevalent in enteric infections. J Infect Dis 2010; 201:1633-1643.

43 Duffy S, Shackelton LA, Holmes EC. Rates of evolutionary change in viruses: patterns and determinants. Nat Rev Genet 2008; 9:267-276.