Human bocavirus (HBoV) 1 is considered an important respiratory pathogen, while the role of HBoV2-4 in clinical disease remains somewhat controversial. Since, they are characterized by a rapid evolution, worldwide surveillance of HBoVs’ genetics is necessary. This study explored the prevalence of HBoV genotypes in pediatric patients with respiratory tract infection in Croatia and studied their phylogeny. Using multiplex PCR for 15 respiratory viruses, we investigated 957 respiratory samples of children up to 18 years of age with respiratory tract infection obtained from May 2017 to March 2021 at two different hospitals in Croatia. Amplification of HBoV near-complete genome or three overlapping fragments was performed, sequenced, and their phylogenetic inferences constructed. HBoV was detected in 7.6% children with a median age of 1.36 years. Co-infection was observed in 82.2% samples. Sequencing was successfully performed on 29 HBoV positive samples, and all belonged to HBoV1. Croatian HBoV1 sequences are closely related to strains isolated worldwide, and no phylogenetic grouping based on mono- or co-infection cases or year of isolation was observed. Calculated rates of evolution for HBoV1 were 10⁻⁴ and 10⁻⁵ substitutions per site and year. Recombination was not detected among sequences from this study.

Purpose. This study investigated the HPIV3 circulating strains in Croatia and whether the other parts of HPIV3 genome (F gene and HN 582 nucleotides fragment) could be equally suitable for genetic and phylogenetic analysis. Methodology. Clinical materials were collected in period 2011–2015 from children suffering from respiratory illnesses. In positive HPIV3 samples viral genome was partially amplified and sequenced for HN and F genes. Obtained sequences were analysed by phylogenetic analysis and genetic characterization was performed. Results. All samples from this study belonged to subcluster C and over a short period of time, genetic lineage C3a gained prevalence over the other C genetic lineages, from 39 % in 2011 to more than 90 % in 2013 and 2014. Phylogenetic classifications of HPIV3 based on the entire HN gene, HN 582 nt fragment and entire fusion (F) gene showed identical classification results for Croatian strains and the reference strains. Molecular analysis of the F and HN glycoproteins, showed their similar nucleotide diversity (Fcds P=0.0244 and HNcds P=0.0231) and similar Ka/Ks ratios (F Ka/Ks=0.0553 and HN Ka/Ks=0.0428). Potential N-glycosylation sites, cysteine residues and antigenic sites are generally strongly conserved in HPIV3 glycoproteins from both our and the reference samples. Conclusion. The HPIV3 subclaster C3 (genetic lineage C3a) became the most detected circulating HPIV3 strain in Croatia. The results indicated that the HN 582 nt and the entire F gene sequences were as good for phylogenetic analysis as the entire HN gene sequence.

Small hydrophobic (SH) gene is one of the mostly diverse genomic regions of human respiratory syncytial virus (HRSV). Its coding region constitutes less than 50% of the complete gene length, enabling SH gene to be highly variable and the SH protein highly conserved. In standard HRSV molecular epidemiology studies, solely sequences of the second hypervariable region of the glycoprotein gene (HVR2) are analyzed. To what extent do the strains identical in HVR2 differ elsewhere in genomes is rarely investigated. Our goal was to investigate whether diversity and inter-genotype differences observed for HVR2 are also present in the SH gene. We sequenced 198 clinical samples collected within a limited area and time frame. In this HRSV collection, rapid and significant changes in HVR2 occurred. Over 20% of strains from this pool (containing HRSV genotypes NA1, ON1, GA5, BA9 and BA10) would be incorrectly assumed to be identical to another strain if only the HVR2 region was analysed. The majority of differences found in SH gene were located in the 5′ untranslated region (UTR). Seven indels were detected, one was genotype GA5 specific. An in-frame deletion of 9 nucleotides (coding for amino acids 49–51) was observed in one of group A strains. Fifteen different SH protein sequences were detected; 68% of strains possessed the consensus sequence and most of others differed from the consensus in only one amino acid (only 4 strains differed in 2 amino acids). The majority of differing amino acids in group A viruses had the same identity as the corresponding amino acids in group B strains. When analysis was restricted to strains with identical HVR2 nucleotide sequences and differing SH protein sequences, 75% of differences observed in the SH ectodomain were located within region coding for amino acids 49–51. Basing HRSV molecular epidemiology studies solely on HVR2 largely underestimates the complexity of circulating virus populations. In strain identification, broadening of the genomic target sequence to SH gene would provide a more comprehensive insight into viral pool versatility and its evolutionary processes.

Background: The families Paramyxoviridae and Pneumoviridae comprise a broad spectrum of viral pathogens that affect human health. The matrix (M) protein of these viruses has a central role in their life cycle. In line with this, molecular characteristics of the M proteins from variable viruses that circulated in Croatia were investigated. Methods: Sequences of the M proteins of human parainfluenza virus (HPIV) 1–3 within the family Paramyxoviridae, human metapneumovirus (HMPV), and human respiratory syncytial virus from the family Pneumoviridae were obtained and analyzed. Results: M proteins were very diverse among HPIVs, but highly conserved within each virus. More variability was seen in nucleotide sequences of M proteins from the Pneumoviridae family. An insertion of 8 nucleotides in the 3′ untranslated region in 1 HMPV M gene sequence was discovered (HR347-12). As there are no samples with such an insertion in the database, this insertion is of interest and requires further research. Conclusion: While we have confirmed that M proteins were conserved among individual viruses, any changes that are observed should be given attention and further researched. Of special interest is inclusion of HPIV2 M proteins in this analysis, as these proteins have not been studied to the same extent as other paramyxoviruses.