Evaluation of the Interactions between Mumps Virus and Guinea Pig.

Mumps is a highly contagious viral disease that can be prevented by vaccination. In the last decade, we have encountered repeated outbreaks of mumps in highly vaccinated populations, which call into question the effectiveness of available vaccines. Animal models are crucial for understanding virus-host interactions, and viruses such as mumps virus (MuV), whose only natural host is the human, pose a particular challenge. In our study, we examined the interaction between MuV and the guinea pig. Our results present the first evidence that guinea pigs of the Hartley strain can be infected in vivo after intranasal and intratesticular inoculation. We observed a significant viral replication in infected tissues up to 5 days following infection and induction of cellular and humoral immune responses as well as histopathological changes in infected lungs and testicles, without clinical signs of disease. Transmission of the infection through direct contact between animals was not possible. Our results demonstrate that guinea pigs and guinea pig primary cell cultures represent a promising model for immunological and pathogenetic studies of the complex MuV infection. IMPORTANCE Understanding of mumps virus (MuV) pathogenesis and the immune responses against MuV infection is limited. One of the reasons is the lack of relevant animal models. This study explores the interaction between MuV and the guinea pig. We demonstrated that all tested guinea pig tissue homogenates and primary cell cultures are highly susceptible to MuV infection and that α2,3-sialylated glycans (MuV cellular receptors) are being abundantly expressed at their surface. The virus remains in the guinea pig lungs and trachea for up to 4 days following intranasal infection. Although asymptomatic, MuV infection strongly activates both humoral and cellular immune response in infected animals and provides protection against virus challenge. Infection of the lungs and testicles after intranasal and intratesticular inoculation, respectively, is also supported by histopathological changes in these organs. Our findings give perspective for application of guinea pigs in research on MuV pathogenesis, antiviral response, and vaccine development and testing.

Mass spectrometry-based investigation of measles and mumps virus proteome.

BACKGROUND: Measles (MEV) and mumps virus (MUV) are enveloped, non-segmented, negative single stranded RNA viruses of the family Paramyxoviridae, and are the cause of measles and mumps, respectively, both preventable by vaccination. Aside from proteins coded by the viral genome, viruses are considered to contain host cell proteins (HCPs). The presence of extracellular vesicles (ECVs), which are often co-purified with viruses due to their similarity in size, density and composition, also contributes to HCPs detected in virus preparations, and this has often been neglected. The aim was to identify which virus-coded proteins are present in MEV and MUV virions, and to try to detect which HCPs, if any, are incorporated inside the virions or adsorbed on their outer surface, and which are more likely to be a contamination from co-purified ECVs. METHODS: MUV, MEV and ECVs were purified by ultracentrifugation, hydrophobic interaction chromatography and immunoaffinity chromatography, proteins in the samples were resolved by SDS-PAGE and subjected to identification by MALDI-TOF/TOF-MS. A comparative analysis of HCPs present in all samples was carried out. RESULTS: By proteomics approach, it was verified that almost all virus-coded proteins are present in MEV and MUV particles. Protein C in MEV which was until now considered to be non-structural viral protein, was found to be present inside the MeV virions. Results on the presence of HCPs in differently purified virus preparations imply that actin, annexins, cyclophilin A, moesin and integrin ß1 are part of the virions. CONCLUSIONS: All HCPs detected in the viruses are present in ECVs as well, indicating their possible function in vesicle formation, or that most of them are only present in ECVs. Only five HCPs were constantly present in purified virus preparations, regardless of the purification method used, implying they are likely the integral part of the virions. The approach described here is helpful for further investigation of HCPs in other virus preparations.

Variability analysis and inter-genotype comparison of human respiratory syncytial virus small hydrophobic gene.

BACKGROUND: 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. METHODS: We sequenced 198 clinical samples collected within a limited area and time frame. In this HRSV collection, rapid and significant changes in HVR2 occurred. RESULTS: 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. CONCLUSIONS: 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.

Genetic diversity of human metapneumovirus in hospitalized children with acute respiratory infections in Croatia.

Human metapneumovirus (HMPV) is recognized as a global and frequent cause of acute respiratory tract infections among people of all ages. The objectives of this study were molecular epidemiology and evolutionary analysis of HMPV strains which produced moderate and severe acute respiratory tract infections in children in Croatia during four consecutive seasons (2011-2014). A total of 117 HMPV-positive samples collected from hospitalized pediatric patients presenting with acute respiratory tract infections and tested by direct immunofluorescence assay were first analyzed by amplifying a part of the F gene. Sixteen samples were further analyzed based on complete F, G, and SH gene sequences. HMPV genome was identified in 92 of 117 samples (78%) and the circulation of multiple lineages of HMPV was confirmed. In 2011, 2012, and 2014, subgroups A2 and B2 co-circulated, while B1 gained prevalence in 2013 and 2014. The study established the presence of a novel subcluster A2c in Croatia. This subcluster has only recently been detected in East and Southeast Asia. This study provides new insights into epidemiology and genetic diversity of HMPV in this part of Europe.

Genetic Variability and Sequence Relatedness of Matrix Protein in Viruses of the Families Paramyxoviridae and Pneumoviridae.

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.

First recorded case of paramyxovirus infection introduced into a healthy snake collection in Croatia.

BACKGROUND: In the present study, we describe the first paramyxovirus infection in a snake collection in Croatia caused by an introduction of new snakes that were not previously tested and didn't show any signs of disease. CASE PRESENTATION: In less than a month after introduction into a healthy colony, new snakes began to show respiratory symptoms (i.e. mouth opening, wheezing, etc.) and died within a month and a half after antibiotic therapy was applied. The same symptoms and a high mortality rate were then observed in in-contact snakes from other collections belonging to different snake families. CONCLUSIONS: Two entries of new snakes in different time periods were recorded and recognized as possible sources of infection. We stress the need for veterinary health control and monitoring of snakes prior to transportation as well as implementing obligatory quarantine measures to minimize the risk of infection among newly established snake groups.

Genetic diversity among human parainfluenza virus type 2 isolated in Croatia between 2011 and 2014.

The dynamics and evolution of the human parainfluenza virus type 2 (HPIV2) in Croatia, and also globally, are largely unknown. Most HPIV2 infections are treated symptomatically outside the hospital setting. Thus, the diagnosis is missing making it difficult to follow the genetic variation and evolution of the HPIV2. This study explores hospitalized HPIV2 cases in Croatia during 4-year period (2011-2014). Most cases in this period were reported in October or November (68.75%) and most of patients were under 2 years of age (81.25%). For molecular analyses, we used the F and HN gene sequences and showed that although both regions are equally suitable for phylogenetic analyses it would be advantageous to use regions longer than 2 kb for HPIV2 analyses of isolates which are spatially and temporally closely related. We show here that the dominant cluster in this area was cluster G3 while only one strain isolated in this period was positioned in the distant cluster G1a. Further monitoring of the HPIV2 will determine whether cluster G3 will remain dominant or it will be overruled by cluster G1a. This will be important for the surveillance of virus circulation in population and significance of the viral infection. J. Med. Virol. 88:1733-1741, 2016. © 2016 Wiley Periodicals, Inc.

Common position of indels that cause deviations from canonical genome organization in different measles virus strains.

BACKGROUND: The canonical genome organization of measles virus (MV) is characterized by total size of 15 894 nucleotides (nts) and defined length of every genomic region, both coding and non-coding. Only rarely have reports of strains possessing non-canonical genomic properties (possessing indels, with or without the change of total genome length) been published. The observed mutations are mutually compensatory in a sense that the total genome length remains polyhexameric. Although programmed and highly precise pseudo-templated nucleotide additions during transcription are inherent to polymerases of all viruses belonging to family Paramyxoviridae, a similar mechanism that would serve to non-randomly correct genome length, if an indel has occurred during replication, has so far not been described in the context of a complete virus genome. METHODS: We compiled all complete MV genomic sequences (64 in total) available in open access sequence databases. Multiple sequence comparisons and phylogenetic analyses were performed with the aim of exploring whether non-recombinant and non-evolutionary linked measles strains that show deviations from canonical genome organization possess a common genetic characteristic. RESULTS: In 11 MV sequences we detected deviations from canonical genome organization due to short indels located within homopolymeric stretches or next to them. In nine out of 11 identified non-canonical MV sequences, a common feature was observed: one mutation, either an insertion or a deletion, was located in a 28 nts long region in F gene 5' untranslated region (positions 5051-5078 in genomic cDNA of canonical strains). This segment is composed of five tandemly linked homopolymeric stretches, its consensus sequence is G6-7C7-8A6-7G1-3C5-6. Although none of the mononucleotide repeats within this segment has fixed length, the total number of nts in canonical strains is always 28. These nine non-canonical strains, as well as the tenth (not mutated in 5051-5078 segment), can be grouped in three clusters, based on their passage histories/epidemiological data/genetic similarities. There are no indications that the 3 clusters are evolutionary linked, other than the fact that they all belong to clade D. CONCLUSIONS: A common narrow genomic region was found to be mutated in different, non-related, wild type strains suggesting that this region might have a function in non-random genome length corrections occurring during MV replication.

Identification of mumps virus protein and lipid composition by mass spectrometry.

BACKGROUND: Mumps virus is a negative-sense, single stranded RNA virus consisting of a ribonucleocapsid core enveloped by a lipid membrane derived from host cell, which causes mumps disease preventable by vaccination. Since virus lipid envelope and glycosylation pattern are not encoded by the virus but dependent on the host cell at least to some extent, the aim of this work was to analyse L-Zagreb (L-Zg) mumps virus lipids and proteins derived from two cell types; Vero and chicken embryo fibroblasts (CEF). Jeryl Lynn 5 (JL5) mumps strain lipids were also analysed. METHODS: Virus lipids were isolated by organic phase extraction and subjected to 2D-high performance thin layer chromatography followed by lipid extraction and identification by matrix-assisted laser desorption/ionization mass spectrometry (MALDI MS). Virus samples were also subjected to gel electrophoresis under denaturating conditions and protein bands were excised, in-gel trypsinized and identified by MS as well as tandem MS. RESULTS: Results showed that lipids of both mumps virus strains derived from Vero cells contained complex glycolipids with up to five monosaccharide units whereas the lipid pattern of mumps virus derived from CEF was less complex. Mumps virus was found to contain expected structural proteins with exception of fusion (F) protein which was not detected but on the other hand, V protein was detected. Most interesting finding related to the mumps proteins is the detection of several forms of nucleoprotein (NP), some of which appear to be C-terminally truncated. CONCLUSIONS: Differences found in lipid and protein content of mumps virus demonstrated the importance of detailed biochemical characterization of mumps virus and the methodology described here could provide a means for a more comprehensive quality control in vaccine production.

Stability, biophysical properties and effect of ultracentrifugation and diafiltration on measles virus and mumps virus.

Measles virus and mumps virus (MeV and MuV) are enveloped RNA viruses used for production of live attenuated vaccines for prophylaxis of measles and mumps disease, respectively. For biotechnological production of and basic research on these viruses, the preparation of highly purified and infectious viruses is a prerequisite, and to meet that aim, knowledge of their stability and biophysical properties is crucial. Our goal was to carry out a detailed investigation of the stability of MeV and MuV under various pH, temperature, shear stress, filtration and storage conditions, as well as to evaluate two commonly used purification techniques, ultracentrifugation and diafiltration, with regard to their efficiency and effect on virus properties. Virus titers were estimated by CCID50 assay, particle size and concentration were measured by Nanoparticle tracking analysis (NTA) measurements, and the host cell protein content was determined by ELISA. The results demonstrated the stability of MuV and MeV at pH <9 and above pH 4 and 5, respectively, and aggregation was observed at pH >9. Storage without stabilizer did not result in structural changes, but the reduction in infectivity after 24 hours was significant at +37 °C. Vortexing of the viruses resulted in significant particle degradation, leading to lower virus titers, whereas pipetting had much less impact on virus viability. Diafiltration resulted in higher recovery of both total and infectious virus particles than ultracentrifugation. These results provide important data for research on all upstream and downstream processes on these two viruses regarding biotechnological production and basic research.

Genomic diversity of mumps virus and global distribution of the 12 genotypes.

The WHO recently proposed an updated nomenclature for mumps virus (MuV). WHO currently recognizes 12 genotypes of MuV, assigned letters from A to N (excluding E and M), which are based on the nucleotide sequences of small hydrophobic (SH) and haemagglutinin-neuraminidase (HN) genes. A total of 66 MuV genomes are available in GenBank, representing eight of the 12 genotypes. To complete this dataset, whole genomes of seven isolates representing six genotypes (D, H, I, J, K and L) and one unclassified strain were sequenced. SH and HN genes of other representative strains were also sequenced. The degree of genetic divergence, predicted amino acid substitutions in the HN and fusion (F) proteins and geographic distributions of MuV strains were analysed based on the updated dataset. Nucleotide heterogeneity between genotypes reached 20% within the SH gene, with a maximum of 9% within the HN gene. The geographic and chronologic distributions of the 12 genotypes were summarised. This review contributes to our understanding of strain diversity for wild type MuV, and the results support the current WHO nomenclature.

Early Evolution of Human Respiratory Syncytial Virus ON1 Strains: Analysis of the Diversity in the C-Terminal Hypervariable Region of Glycoprotein Gene within the First 3.5 Years since Their Detection.

OBJECTIVE: Characterization of the phylogeny and diversity of human respiratory syncytial virus (HRSV) genotype ON1 that occurred during its early evolution (within the first 3.5 years since the detection of the first ON1 strains). ON1 strains have a 72-nucleotide-long in-frame duplication within the second hypervariable domain of the glycoprotein gene (HVR2). METHODS: All available HVR2 sequences of strains belonging to the ON1 genotype published prior to June 20, 2014 were collected. Multiple sequence alignments, phylogeny, phylogeography, sequence clustering and putative protein analyses were performed. RESULTS: The worldwide spread and diversification of ON1 strains are presented. Only in a minority of ON1 strains do the two replicas remain identical, and various ON1 strains possess common differences between the first and the second copy (segments A and B). Mutations of the progenitor sequence were more frequent in segment B, a higher overall diversity on the protein level and more putative glycosylation sites exist in segment B, and, unlike in segment A, positive selection acts on that protein region. CONCLUSIONS: The fast spread of the novel HRSV genotype ON1 has been accompanied by its rapid concurrent diversification. Differences in variability of the two replicas within HVR2 were detected, with C-terminal replica being more variable.

Critical factors for the replication of mumps virus in primary chicken embryo fibroblasts defined by the use of design of experiments (DoE).

Live attenuated vaccines against mumps virus (MuV) have been traditionally produced by passaging the virus in the embryonated chicken eggs or primary chicken embryo fibroblasts (CEFs). Virus propagation on these cell substrates enables successful virus attenuation and retains it sufficiently antigenic to induce lasting protective immunity in humans. The aim of this study was to identify critical factors for MuV replication in primary CEFs grown on a small-scale level in order to explore possibilities for improvements in the virus replication and yield. The effect of differently prepared cells, culturing conditions, and infection conditions on virus yield was estimated by employing statistical design of experiments (DoE) methodology. Our results show that the preparation of primary CEFs and the way of their infection substantially impact virus yield and are critical for efficient MuV replication. These process parameters should be considered in further process optimization. We also demonstrate the applicability of DoE in optimization of virus replication as a crucial step in obtaining high virus yields.

A study of the genetic variability of human respiratory syncytial virus in Croatia, 2006-2008.

Human respiratory syncytial virus (HRSV) is a common etiological agent of acute lower respiratory tract disease in infants. The molecular epidemiology of HRSV in Croatia over four consecutive seasons (from 2006 to 2008) was investigated. A total of 72 HRSV samples were chosen from 696 screened cases in a pediatric clinic in Zagreb. Molecular characterization of HRSV revealed the predominance of HRSV group B viruses in the first two epidemic seasons and HRSV group A viruses in the next two seasons. According to the phylogenetic analysis, NA1 and BA9 were the predominant circulating HRSV genotypes detected during the study. Overall, 82.9% of all HRSV A strains belonged to the NA1 genotype. The HRSV B genotype BA9, detected in two consecutive seasons (2006 and 2007), was the predominant circulating HRSV B genotype, accounting for 80.6% of all HRSV B strains. This study provides data on the circulation pattern of HRSV genotypes in Croatia and their molecular characterization.

Is Better Standardization of Therapeutic Antibody Quality in Emerging Diseases Epidemics Possible?

During the ongoing COVID-19 epidemic many efforts have gone into the investigation of the SARS-CoV-2-specific antibodies as possible therapeutics. Currently, conclusions cannot be drawn due to the lack of standardization in antibody assessments. Here we describe an approach of establishing antibody characterisation in emergent times which would, if followed, enable comparison of results from different studies. The key component is a reliable and reproducible assay of wild-type SARS-CoV-2 neutralisation based on a banking system of its biological components - a challenge virus, cells and an anti-SARS-CoV-2 antibody in-house standard, calibrated to the First WHO International Standard immediately upon its availability. Consequently, all collected serological data were retrospectively expressed in an internationally comparable way. The neutralising antibodies (NAbs) among convalescents ranged from 4 to 2869 IU mL-1 in a significant positive correlation to the disease severity. Their decline in convalescents was on average 1.4-fold in a one-month period. Heat-inactivation resulted in 2.3-fold decrease of NAb titres in comparison to the native sera, implying significant complement activating properties of SARS-CoV-2 specific antibodies. The monitoring of NAb titres in the sera of immunocompromised COVID-19 patients that lacked their own antibodies evidenced the successful transfusion of antibodies by the COVID-19 convalescent plasma units with NAb titres of 35 IU mL-1 or higher.

Clinical and molecular characterization of a parechovirus type 1 outbreak in neonates in Croatia.

During July 2009 an outbreak in neonates represented with gastrointestinal and respiratory symptoms was observed at the Neonatal Postintensive Care Unit, Clinical Hospital Center, Zagreb. Human parechovirus type 1 (HPeV1) was isolated from seven patients, one of whom was asymptomatic. All but one were premature neonates with serious underlying conditions, and all recovered fully after several days. In order to characterize the HPeV1s, sequencing of the VP1/2A region was conducted on six isolates from the outbreak and four isolates detected in Croatia in 2008 and 2007. The analysis of sequence similarity showed that the nucleotide identity between the prototype strain (Harris) and HPeV1 isolated in Croatia was 76.5-77.5%. Croatian strains from 2007 and 2009 clustered together with strains from the Netherlands and Germany detected in 2003 and 2006, respectively, while strains from 2008 clustered with the strain from Finland detected in 2000. Change of the dominant strains each year may suggest antigenic variation as a result of viral response to specific immunity of the target population.

Comparative analysis of CE-SSCP to standard RFLP-CE-FLA method in quantification of known viral variants within an RNA virus quasispecies.

RNA viruses display the highest replication error rate in our biosphere, leading to highly diverse viral populations termed quasispecies. The gold standard method for detection and quantification of variants in a quasispecies is cloning and sequencing, but it is expensive, laborious and time consuming. Therefore, other mutation detection approaches, including SSCP, are often used. In this study, we demonstrate development and the usage of a CE-SSCP method for quantification of two nearly identical viral variants in heterogenic population of a mumps virus strain and its comparison to RFLP-CE-fragment length analysis (RFLP-CE-FLA). Analyzed PCR fragments were of the same size (245 bp) with one difference in their nucleotide sequence. The limit of detection of both methods was at 5% of the minor variant. When PCR amplicons of the two variants were pooled, methods' results were very similar. On the contrary, the quantification results of samples in which variants were mixed prior to PCR showed substantial difference between the two methods. Our results indicate that although both methods can be used for detection and monitoring of a specific mutation within a viral population, caution should be taken when quantitative analysis of complex samples is based solely on results of one method.

Influence of Ribavirin on Mumps Virus Population Diversity.

Frequent mumps outbreaks in vaccinated populations and the occurrence of neurological complications (e.g., aseptic meningitis or encephalitis) in patients with mumps indicate the need for the development of more efficient vaccines as well as specific antiviral therapies. RNA viruses are genetically highly heterogeneous populations that exist on the edge of an error threshold, such that additional increases in mutational burden can lead to extinction of the virus population. Deliberate modulation of their natural mutation rate is being exploited as an antiviral strategy and a possibility for rational vaccine design. The aim of this study was to examine the ability of ribavirin, a broad-spectrum antiviral agent, to introduce mutations in the mumps virus (MuV) genome and to investigate if resistance develops during long-term in vitro exposure to ribavirin. An increase in MuV population heterogeneity in the presence of ribavirin has been observed after one passage in cell culture, as well as a bias toward C-to-U and G-to-A transitions, which have previously been defined as ribavirin-related. At higher ribavirin concentration, MuV loses its infectivity during serial passaging and does not recover. At low ribavirin concentration, serial passaging leads to a more significant increase in population diversity and a stronger bias towards ribavirin-related transitions, independently of viral strain or cell culture. In these conditions, the virus retains its initial growth capacity, without development of resistance at a whole-virus population level.

Population Variability Generated during Rescue Process and Passaging of Recombinant Mumps Viruses.

Recombinant mumps viruses (MuVs) based on established vaccine strains represent attractive vector candidates as they have known track records for high efficacy and the viral genome does not integrate in the host cells. We developed a rescue system based on the consensus sequence of the L-Zagreb vaccine and generated seven different recombinant MuVs by (a) insertion of one or two additional transcription units (ATUs), (b) lengthening of a noncoding region to the extent that the longest noncoding region in MuV genome is created, or (c) replacement of original L-Zagreb sequences with sequences rich in CG and AT dinucleotides. All viruses were successfully rescued and faithfully matched sequences of input plasmids. In primary rescued stocks, low percentages of heterogeneous positions were found (maximum 0.12%) and substitutions were predominantly obtained in minor variants, with maximally four substitutions seen in consensus. ATUs did not accumulate more mutations than the natural MuV genes. Six substitutions characteristic for recombinant viruses generated in our system were defined, as they repetitively occurred during rescue processes. In subsequent passaging of primary rescue stocks in Vero cells, different inconsistencies within quasispecies structures were observed. In order to assure that unwanted mutations did not emerge and accumulate, sub-consensus variability should be closely monitored. As we show for Pro408Leu mutation in L gene and a stop codon in one of ATUs, positively selected variants can rise to frequencies over 85% in only few passages.

Prevalence and Molecular Characterization of Human Bocavirus Detected in Croatian Children with Respiratory Infection.

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-4 and 10-5 substitutions per site and year. Recombination was not detected among sequences from this study.