Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.
Advertisement
Scientific Reports volume 13, Article number: 15913 (2023)
153
2
Metrics details
To understand the epidemiological and genetic characteristics of B19V, a multiple-province surveillance of patients with febrile rash illnesses (FRIs) were conducted in China during 2009 ~ 2021. The clinical specimens of 3,820 FRI patients were collected and tested for B19V DNA. A total of 99 (2.59%) patients were positive for B19V, and 49 (49.49%) were children under 5 years old. B19V infections occurred throughout the year without obvious seasonal pattern. Ten NS1-VP1u sequences and seven genome sequences were obtained in this study, identified as subgenotype 1a. Combined with the globally representative genome sequences, no temporal and geographic clustering trends of B19V were observed, and there was no significant correlation between B19V sequences and clinical manifestations. The evolutionary rate of the B19V genome was 2.30 × 10–4 substitutions/site/year. The number of negative selection sites was higher than that of positive selection sites. It was the first to comprehensively describe the prevalence patterns and evolutionary characteristics of B19V in FRI patients in China. B19V played the role in FRI patients. Children under 5 years old were the main population of B19V infection. Subgenotype 1a was prevalent in FRI patients in China. B19V showed a high mutation rate, while negative selection acted on the genome.
As a member of the family Parvoviridae, human parvovirus B19 (B19V) was a global and common infectious pathogen in humans1. B19V was first identified by Cossart et al.2 during the evaluation of tests for hepatitis B virus surface antigen in 1975. Subsequently, B19V infection has been reported worldwide. The transmission of B19V infection usually occurs through the respiratory route, and infections can also be transmitted vertically from mother to fetus and through the transfusion of blood products and bone marrow transplants3,4,5. B19V infection is common in childhood, and it can also occur throughout adulthood albeit at a lower rate6. The findings showed a significant negative correlation between viremia rates and age, and the positive rate of B19V DNA decreased with age, from 2.24% in 19–30 years to 0.87% in 41–50 years7,8. However, the prevalence of IgG antibodies directed against B19V ranges from 2 to 15% in children 1 to 5 years old, 15 to 60% in children 6 to 19 years old, 30 to 60% in adults, and more than 85% in the geriatric population9,10,11,12. The manifestations of B19V infection depend on the age, immunity and hematologic status of the host1,13. Most B19V infections are asymptomatic. The common clinical manifestation of B19V infection is erythema infectiosum in children, while arthropathy is a more common manifestation of infection in adults, particularly in women14. In immunocompromised hosts, persistent B19V infection presents with pure red cell aplasia and chronic anemia13,15. B19V is now recognized as the only etiologic agent of erythema infectiosum16,17,18. Due to similar eruption symptoms, erythema infectiosum can be confused with rubella19,20,21,22.
B19V is a small nonenveloped single-stranded DNA virus with an approximately 5.5-kb long genome. This DNA encodes six viral proteins, among which the three major proteins are nonstructural protein 1 (NS1), viral protein 1 (VP1) and viral protein 2 (VP2). In addition, there are three small nonstructural proteins of 7.5 kDa, X, and 11 kDa. Based on the phylogenetic analysis of the NS1-VP1u region, B19V was classified into three genotypes: genotype 1, genotype 2, and genotype 323. With the increasing number of genome sequences obtained, B19V was further subdivided into subgenotypes, which were useful for molecular epidemiological studies. Genotype 1 was further divided into two subgenotypes, 1a and 1b, and genotype 3 was divided into two subgenotypes, 3a and 3b24,25.
To date, investigations of B19V have been carried out by many researchers in recent years, including in China. However, most of them focused on viral infection in blood products from healthy donors, organ transplantation and the status of B19V infection in pregnant women26,27,28,29,30. There have been few studies on the prevalence of B19V in patients with febrile rash illnesses (FRIs). From January 2009 to June 2021, active surveillance of FRIs was conducted in ten provinces of China, including Anhui, Beijing, Hebei, Henan, Hunan, Shandong, Shanxi, Shaanxi, Shanghai and Xinjiang Uygur Autonomous Region. Multiple viral pathogens were tested in the clinical specimens obtained from FRI patients, including measles virus, rubella virus, human enterovirus, varicella-zoster virus, dengue virus, B19V, Epstein‒Barr virus, and human herpes virus 6, which covered the most common viral pathogens of FRIs. In this study, the epidemiological and genetic characteristics of B19V were comprehensively analyzed based on the surveillance data of FRI patients in ten provinces in China from 2009 to 2021.
In total, the clinical specimens of 3,820 FRI patients were tested for B19V DNA during 2009 ~ 2021, of which 99 were positive for B19V (Table 1). The positive rate of B19V was 2.59% in FRI patients. The positive rate of B19V was significant higher in females than in males (p = 0.040). The median age of B19V-positive patients was 5 years old (IQR: 1–20). Among the 99 B19V-positive patients, 49 cases (49.49%) were children under 5 years old, 31 cases (31.31%) were adolescents aged 5–17 years old, 19 cases (19.19%) were adults aged 18–59 years old, and no B19V-positive patient was detected in the elderly aged ≥ 60 years old. The positive rate was highest in the age group of 5 ~ 17 years old (3.59%), followed by < 5 years old (2.57%), 18–59 years old (2.08%), and ≥ 60 years old (0.00%). There was a significant difference in the positive rates of B19V in different age groups (p = 0.046). B19V was detected in FRI patients throughout the year, with no obvious seasonal pattern (p = 0.155). The positive rate of B19V in autumn was highest (3.51%), followed by winter (3.23%), summer (2.70%) and spring (2.02%). The B19V-positive rates in FRI patients were slightly higher in autumn and winter than in spring and summer, but the difference was not significant. In terms of regional distribution, the positive rate of B19V in FRI patients in Northern China (3.50%) was higher than that in Southern China (0.44%), and the difference was statistically significant (p = 0.001).
A total of 10 NS1-VP1u sequences were obtained from the clinical specimens of 99 FRI patients with B19V infection in this study, including 4 from Shanxi, 3 from Shanghai, 2 from Shaanxi, and 1 from Henan. These sequences were preliminarily identified by a BLAST search in the NCBI Nucleotide GenBank database (https://blastncbi.nlm.nih.gov/) as subgenotype 1a of genotype 1, and the nucleotide similarity was 98.7% ~ 100.0%.
Based on the phylogenetic tree constructed with the NS1-VP1u dataset, B19V sequences could be divided into three genotypes (Fig. 1). The average inter-genotype P-distances of B19V ranged from 0.069 to 0.103, and the average intra-genotype P-distances were 0.015 for genotype 1, 0.021 for genotype 2, and 0.025 for genotype 3. The average inter-subgenotype P-distances in the same genotype were 0.031 ~ 0.032, and the average intra-subgenotype P-distances were 0.011 ~ 0.022. Combined with the globally representative B19V sequences downloaded from GenBank, no trend in temporal or geographical clustering of B19V was observed in the phylogenetic tree. In addition, the sequences obtained from healthy donor blood and the clinical specimens of the B19V patients with different clinical manifestations were also interleaved in the phylogenetic tree.
Phylogenetic tree of the B19V sequence based on the NS1-VP1u region. A phylogenetic tree was constructed with the maximum likelihood method. Bootstrap values above 80% are shown (1000 replicates). The blue diamond represents the reference sequences of genotype or subgenotype; the green solid circle represents the Chinese representative sequences; the red solid circle represents the sequence obtained in this study.
A total of 29 Chinese representative sequences of B19V were analyzed in this study, including 10 sequences from FRI patients and 19 sequences downloaded from GenBank. The 29 Chinese sequences belonged to subgenotypes 1a, 1b and 3b, of which 79.31% (23/29) were subgenotype 1a (Supplementary Table S1). All sequences obtained from FRI patients belonged to subgenotype 1a. In addition, the sequences of subgenotype 1a were also obtained from healthy donor blood and the clinical specimens of patients with kidney transplants, HIV, and pneumoniae, while the sequences of subgenotypes 1b and 3b were all sourced from healthy donor blood.
A total of 7 B19V genome sequences were obtained in this study, including 2 from Shanxi, 2 from Shanghai, 2 from Shaanxi, and 1 from Henan. The 7 genome sequences all belonged to subgenotype 1a, and the nucleotide similarity was 98.5% ~ 99.7%. The phylogenetic tree was conducted based on the genome dataset of B19V. The topologies of these sequences were almost identical in the ML trees of the genome and NS1-VP1u region (Supplementary Figure S1). The average intergenotype P-distances ranged from 0.090 to 0.130, and the average intragenotype P-distances of B19V were 0.017 for genotype 1, 0.026 for genotype 2, and 0.041 for genotype 3. The average inter-subgenotype P-distances in the same genotype were 0.050 ~ 0.055, and the average intra-subgenotype P-distances were 0.016 ~ 0.031. In terms of genetic distance, genotype 1 was far from genotypes 2 and 3, while genotype 3 showed higher genetic diversity than genotype 1 and genotype 2. The genetic diversity of genotype 1 was the lowest, although it was widely prevalent in the world.
There were some genotype-specific amino acid sites in the six proteins of B19V, such as 26 in NS1, 7 in 7.5 kDa, 10 in VP1, 2 in X, 4 in VP2, and 13 in 11 kDa (Supplementary Table S2). The amino acid mutation rates of the six proteins from high to low were 11 kDa (39.36%, 37/94), X (23.46%, 19/81), 7.5 kDa (18.92%, 14/74), NS1 (12.67%, 85/671), VP1 (11.78%, 92/781) and VP2 (8.30%, 46/554). The amino acid mutation rates of the three major proteins were lower than those of the three small proteins.
Because 6 sequences missed the date information of case onset and sample collection, 52 out of 58 sequences in the genome dataset were used for analysis of B19V. By using BEAST package (version 1.10.4), no time trend for B19V was observed in the BEAST tree (Supplementary Figure S2). And, the evolutionary rate of the B19V genome was estimated to be 2.30 × 10–4 substitutions/site/year (95% HPD: 1.60 × 10–4 ~ 3.04 × 10–4 substitutions/site/year). Among them, the evolutionary rate of the subgenotype 1a B19V genome was 2.64 × 10–4 substitutions/site/year (95% HPD: 2.05 × 10–4 ~ 3.26 × 10–4 substitutions/site/year). In the B19V genome, the evolutionary rates and nucleotide similarity of genes encoding different proteins were variable. The evolutionary rates of NS1, VP1 and VP2 were similar, slightly higher than that of the B19V genome. 11 kDa was the fastest evolving small nonstructural protein in the genome, with an evolutionary rate of up to 3.58 × 10–4 substitutions/site/year. The 7.5 kDa protein was the slowest evolving protein in the genome, with an evolutionary rate of 9.03 × 10–5 substitutions/site/year (Table 2).
To better understand the evolutionary dynamics of B19V, selection pressures on six proteins in the genome were also analyzed in this study. The ω values for six proteins were all less than 1, indicating that these proteins were all under negative selection pressure. Through the analysis with the FEL, SLAC and FUBAR methods, negative selection sites were found in the coding regions of six proteins, and their number was far higher than that of positive selection sites (Table 3). The number of negative selection sites in the coding regions of the three major proteins varied with different methods. In addition, through the analysis with the FEL, FUBAR and MEME methods, a small number of positive selection sites were also found in the coding genes of NS1, VP1 and VP2, with a P value of < 0.05 (MEME and FEL methods) or a posterior probability of > 0.95 (FUBAR method) (Table 3). A total of 4 positive selection sites were found in the NS1 protein, among which 3 sites (V151 → G, I181 → M/L, A206 → N) were located at the N-terminus of the NS1 protein (MEME), and one site (F554 → S) was located at the C-terminus of the NS1 protein (FUBAR). Notably, a nucleotide substitution (T1661C) in the NS1 protein (F554 → S) was found to be a nonsynonymous mutation with positive selection pressure, while the same substitution was a synonymous mutation at the first nucleotide of codon 65 (TTG → CTG) in the 7.5 kDa protein, where it was found to be subject to negative selection pressure. Four positive selection sites were also found in the VP1 protein, of which two sites (E4 → T/K/Q/N/A, P96 → S) were located in VP1u, and two sites (VP1:S548, T716 → S; VP2:S321, T489 → S) were located in the C-terminus of the VP1/VP2 protein (MEME). In particular, a highly active mutation was found at codon 4 in VP1u (FUBAR).
B19V was one of the viral pathogens in FRI patients, especially for children under 5 years of age. It was frequently detected in measles- and rubella-negative patients in countries with measles/rubella elimination or near elimination20. For example, in Campinas, Brazil, a setting of low measles and rubella virus transmission, 2.4% of FRI patients under 40 years of age were laboratory confirmed to be positive for B19V31. In Bulgaria, 56.18% of measles- and rubella-negative patients were positive for B19V viral DNA, and 48.97% were positive for B19V-IgM32. In a survey of children under 15 years of age with measles- and rubella-like illness in Iran, 10.81% of children were positive for B19V viral DNA, while 19.21% and 18.87% were positive for B19V-IgM and B19V-IgG, respectively20. The positive rate of B19V among FRI patients varies widely and depends on geographical region, detection method, measles and rubella control status and age groups tested. In China, the positive rates of B19V in FRI patients ranged from 0.43 to 9.52%33,34,35,36,37. In this study, the positive rate of B19V in FRI patients was 2.59%, similar to the previous reports. Surveillance data in this study showed that nearly half of B19V patients were children under 5 years old. B19V infection in FRI patients occurred throughout the year without obvious seasonal patterns, but the B19V-positive rates in FRI patients were slightly higher in autumn and winter than in spring and summer in China. In contrast, a study in Iran found slightly different results from ours, as it found that the majority of B19V infections were concentrated between March and June, with May being the most dominant month of the year20. The positive rate of B19V in FRI patients in Northern China was significantly higher than that in Southern China. Since there are very few epidemiological studies on B19V in FRI patients in China and even in the world, further studies are needed to fully characterize the prevalence of B19V in the future.
To date, three genotypes of B19V have been recognized as distinct circulating lineages worldwide. Compared with genotypes 2 and 3, genotype 1 was the predominant genotype of B19V, circulating in all continents, especially subgenotype 1a. Similar to the predominant genotype of B19V circulating globally, most of the B19V sequences detected in China were genotype 1, including the sequences obtained in FRI patients in this study. So far, the Chinese B19V sequences of genotype 1 have been obtained from healthy donor blood and patients with FRIs, kidney transplants, HIV, etc. In addition, a few B19V sequences of subgenotypes 1b and 3b were detected from healthy donor blood27.
No B19V sequence of genotype 2 has been detected in China. Genotype 2 was sporadically detected in a few countries of Europe, including Finland, Germany, and Italy38,39,40. Because it was very rarely detected as circulating virus and frequently detected as persistent in tissues of elderly individuals, genotype 2 was considered an ancient genotype41. The data suggested that a global replacement of genotype 2 with genotype 1 may have occurred in the 1960s41. Genotype 3 was highly dominant in West Africa, which has been regarded as an endemic region for genotype 3 B19V infection25. However, it now tends to spread outside of Africa, especially subgenotype 3b, which is increasingly being found in some countries of Asia, Europe, and America41.
Despite being a DNA virus, B19V showed a high mutation rate. At the nucleotide level of the genome, the intergenotype genetic distances among the three genotypes of B19V ranged from 0.090 to 0.130. The phylogenetic analysis of B19V showed that genotype 1 was far from genotypes 2 and 3 in terms of genetic distance. Compared with genotypes 2 and 3, genotype 1 showed low diversity although a large number of B19V genotype 1 strains are widespread around the world. The genetic variation within genotype 1 is usually less than 2%, and that of genotype 2 and genotype 3 is usually between 3 and 10%, similar to the previous report41. Moreover, no temporal or geographic clustering trend among three genotypes was observed in phylogenetic tree of B19V.
B19V infection is associated with a variety of clinical manifestations. In addition to erythema infectiosum, B19V can also cause aplastic anemia, transient aplastic crisis, hepatic inflammation, myocarditis, arthritis, and neurological disorders. The clinical manifestations of B19V infection depend on the interplay among the pathogenetic potential of the virus, its adaptation to different cellular environments, and the physiological and immune status of the infected individuals41,42. The current evidence suggests that all genotypes have been found in patients with B19V-related symptoms as well as asymptomatic individuals43. In this study, all the B19V sequences obtained from the patients with FRIs belonged to genotype 1. Except for erythema infectiosum, genotype 1 could also be detected in patients with different B19V clinical manifestations, while genotype 3 was mostly detected in patients who had undergone liver/kidney transplantation or patients with leukemia44,45. Genotype 2 has only been detected in the bone marrow of older patients with cytopenias of unknown origin or chronic lymphocytic leukemia44,46. After acute infection, residual viral DNA of B19V can remain in tissue for decades or even throughout life47,48. However, no clear correlations were observed between sequences and clinical manifestations in the phylogenetic tree of B19V.
In the evolutionary analysis, the average evolutionary rate of the B19V genome in this study (2.30 × 10–4 substitutions/site/year) was similar to that of the subgenotype 1a B19V genome (2.64 × 10–4 substitutions/site/year). These two evolutionary rates were in the range of 10–4 substitutions/site/year, which was consistent with the previous studies25,49,50,51. The evolutionary rates of NS1, VP1 and VP2 were similar, slightly higher than that of the B19V genome. The two small nonstructural proteins, 11 kDa and 7.5 kDa, were the fastest and slowest evolving proteins in the genome, respectively. The amino acid mutation rates of the three major proteins were lower than those of the three small proteins, suggesting that the high mutation rate may enable progressive diversification of viral isolates50,51. The selection pressure analysis in this study revealed that ω values in different coding regions were less than 1. Negative selection acted on the B19V genome, including the six proteins, indicating that the evolution of B19V was under strong negative selection, similar to previously published results49,50. Only a few positive selection sites were found in NS1, VP1 and VP2 in this study. While no positive selection site was found in the functional domain of NS1 protein, including the RI motif-replication-initiator motif (codons 79–147), NTP-binding and helicase domains (codons 320–416), and putative transactivation domains: TAD1 (codons 416–424), TAD2 (codons 523–531), and TAD3 (codons 566–576)50. Four positive selection sites in VP1 codons (4, 96, 548 and 716) were found in the dominant VP1/VP2 immune epitopes (VP1-F1: aa 2–100, VP1-F3 = VP2: aa 227–781), which could elicit a long-lasting immune response52. These positive selection sites in the antigenic domain of VP1/VP2 may be driven by immune response pressures.
B19V played the role in FRI patients in China. Children under 5 years old are the main population of B19V infection. Genotype 1 was prevalent in FRI patients with B19V infection in China, especially subgenotype 1a. B19V showed a high mutation rate in viral evolution, while negative selection acted on most of the B19V genome. These findings would provide a better understanding of the burden of B19V infections and the development of potential vaccines against B19V.
Through our multiple-province surveillance of FRI patients for more than 10 years, the data obtained in this study was the first to at least partially reflect the prevalence of B19V in FRI patients in China. However, due to the low viral yield in clinical samples and failure to grow in cell cultures of B19V, only a few NS1-VP1u sequences and genome sequences were obtained in FRI patients. The limited epidemiological data and sequences of B19V may not comprehensively reflect the prevalence of B19V in China. Therefore, a systematic and continuous surveillance of B19V is crucial in future.
Active FRI surveillance was conducted from January 2009 to June 2021 in multiple provinces in China. The surveillance enrolled FRI patients of all ages. The patients with the clinical manifestation of fever (> 37.5 °C, lasting more than one day) and systemic or localized skin or mucosal rash were enrolled in sentinel hospitals after they or their parents/guardians provided informed consent. Detailed clinical and epidemiological information was collected through a standardized case reporting form by the staff of sentinel hospitals. This study was approved by the Ethics Committee of the National Institute for Viral Disease Control and Prevention, Chinese Centers for Disease Control and Prevention and all methods were performed in accordance with the relevant guidelines. In this study, only clinical specimens including throat swabs or sera were collected and detected from FRI patients and no human experimentation was involved.
The epidemiological dataset of FRI patients with B19V infection was established and analyzed in this study. Because age data did not conform to a normal distribution, the age of FRI patients with B19V infection is presented as the median and interquartile range (IQR). The age groups were divided as follows: < 5, 5 ~ 17, 18 ~ 59, and ≥ 60 years old. In order to analyze the regional differences in the prevalence of B19V, ten provinces were divided into Northern China (Xinjiang Uygur Autonomous Region, Beijing, Hebei, Henan, Shaanxi, Shandong, Shanxi) and Southern China (Anhui, Hunan, Shanghai) (Supplementary Table S3). Descriptive statistical analyses were performed to analyze the positive rates of B19V by gender, age group, season and region during 2009 ~ 2021. Chi-square tests were used to analyze the significant differences in frequency data in SPSS 26.0, and p < 0.05 was considered to be statistically significant.
Clinical specimens, including throat swab or serum specimens, were collected from FRI patients. According to the manufacturer’s instructions, viral DNA was extracted from clinical specimens of FRI patients by using a viral DNA/RNA nucleic acid extraction and purification kit (Xi’an Tianlong Technology Co., Ltd, China) and stored at −80 °C until testing. The preliminary screening of the eight FRI-associated viruses mentioned above was performed by real-time PCR53. The 994 nt NS1-VP1u region (positions 2117 to 3110 in NC_000883.2) and the near complete coding B19V genome (positions 616 to 5174 in NC_000883.2) were amplified by nested PCR from the B19V-positive specimens by using Platinum PCR SuperMix (Invitrogen, USA). The pairs of primers for NS1-VP1u and the genome are listed in Supplementary Table S4. The first and second rounds of PCR amplification were performed under the following conditions: denaturation at 94 °C for 4 min; amplification for 40 cycles of 94 °C for 1 min, 45 °C for 1 min, and 72 °C for 1 min 40 s; and extension for 10 min at 72 °C at the end of the reaction. The PCR-positive products of B19V were identified by 1.5% agarose gel electrophoresis and sequenced by Sangon Bioengineering (Shanghai) Co., Ltd. The raw B19V sequence was edited and assembled by Sequencher software (Version 5.4.5).
A total of 85 representative NS1-VP1u sequences were downloaded from GenBank, which consisted of 6 reference sequences of the B19V genotype/subgenotype, 19 sequences from China, and 60 sequences from other countries. Fifty-one representative B19V genome sequences were downloaded from GenBank from 13 countries worldwide during 1973 ~ 2021. These representative sequences were derived from plasma from healthy blood donors and the clinical specimens of B19V patients with different clinical manifestations, including erythema infectiosum, arthritis, aplastic crisis, hepatic inflammation, myocarditis, etc.
Multiple nucleotide sequences, including the sequences obtained in this study and the representative sequences downloaded from GenBank, were aligned and edited using MAFFT software (Version 7.4.5) and MEGA software (Version 7.0). Subsequently, the NS1-VP1u dataset and the genome dataset were generated and used for genotype identification and viral evolutionary analysis of B19V, respectively (Supplementary Table S1).
Based on NS1-VP1u and the genome dataset of B19V, phylogenetic trees were constructed by the maximum likelihood (ML) approach in MEGA software. The phylogenetic tree was tested with 1000 bootstrap replications, and bootstrap values greater than 80% are indicated on the trees. The P-distances within and between genotypes/subgenotypes at the nucleotide level and the amino acid mutation rates in six proteins were calculated by using MEGA software, and nucleotide similarities were calculated by using BioEdit software (Version 7.0.5.1).
The genome dataset was analyzed with the BEAST package (version 1.10.4) to estimate the evolutionary rate of B19V and the BEAST tree of B19V was also conducted in this study. Initially, the nucleotide substitution model was determined by jModelTest2 (version 2.1.6) online software. With Bayesian skyline as coalescent tree priors, three different molecular clock models were implemented in the BEAST analysis, including strict clock, uncorrelated exponential relaxed clock, and uncorrelated log-normal relaxed clock. Log marginal likelihoods were determined by stepping-stone sampling (SS). The best fit models for the genome and six proteins of B19V were determined by Bayesian factor analysis. The Markov Chain Monte Carlo method was performed for 50 million generations and sampled so that 1,000 trees were generated. Finally, the convergence of the chains and the effective sample size (> 200) were determined by Tracer software (version 1.7.1), and the uncertainty of the parameter estimates was assessed by a 95% HPD interval.
Selection pressure was also analyzed in six proteins of B19V based on the genome dataset. DnaSP6 software (Version 6.0) was used to calculate the ω values (ω = dN/dS) of B19V, where dN represents the nonsynonymous substitution rate and dS represents the synonymous substitution rate. Gene-specific and site-specific selection pressures were measured as the value dN/dS at each codon site and estimated using four different codon-based maximum-likelihood methods (FEL, SLAC, FUBAR and MEME), and Tamura-Nei model (TrN) or Hasegawa-Kishino-Yano model (HKY85) were used as nucleotide substitution models50. All methods were obtained available at the DataMonkey online version of the Hyphy package (www.datamonkey.org) with significance levels set at p < 0.05, and the posterior probability of the FUBAR algorithm was > 0.95.
Seven genome sequences and three NS1-VP1u sequences of B19V in this study were submitted to GenBank with accession numbers OR533486 ~ OR533495 (Supplementary Table S1).
This study was approved by the ethical review committee of the National Institute for Viral Disease Control and Prevention, Chinese Centers for Disease Control and Prevention. The informed consents were signed by patients or their legal guardians.
The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
Heegaard, E. D. & Brown, K. E. Human parvovirus B19. Clin. Microbiol. Rev. 15, 485–505. https://doi.org/10.1128/CMR.15.3.485-505.2002 (2002).
Article PubMed PubMed Central Google Scholar
Cossart, Y. E., Field, A. M., Cant, B. & Widdows, D. Parvovirus-like particles in human sera. Lancet 1, 72–73. https://doi.org/10.1016/s0140-6736(75)91074-0 (1975).
Article CAS PubMed Google Scholar
Azzi, A., Morfini, M. & Mannucci, P. M. The transfusion-associated transmission of parvovirus B19. Transfus. Med. Rev. 13, 194–204. https://doi.org/10.1016/s0887-7963(99)80033-9 (1999).
Article CAS PubMed Google Scholar
Broliden, K. Parvovirus B19 infection in pediatric solid-organ and bone marrow transplantation. Pediatr. Transplant 5, 320–330. https://doi.org/10.1034/j.1399-3046.2001.00035.x (2001).
Article CAS PubMed Google Scholar
Rogo, L. D., Mokhtari-Azad, T., Kabir, M. H. & Rezaei, F. Human parvovirus B19: a review. Acta Virol. 58, 199–213. https://doi.org/10.4149/av_2014_03_199 (2014).
Article CAS PubMed Google Scholar
Young, N. S. & Brown, K. E. Parvovirus B19. N Engl. J. Med. 350, 586–597. https://doi.org/10.1056/NEJMra030840 (2004).
Article CAS PubMed Google Scholar
Abdelrahman, D., Al-Sadeq, D. W., Smatti, M. K., Taleb, S. A., AbuOdeh, R. O., Al-Absi E. S., et al. Prevalence and Phylogenetic Analysis of Parvovirus (B19V) among Blood Donors with Different Nationalities Residing in Qatar. Viruses, 13: (2021) doi:https://doi.org/10.3390/v13040540.
Kooistra, K., Mesman, H. J., de Waal, M., Koppelman, M. H. G. M. & Zaaijer, H. L. Epidemiology of high-level parvovirus B19 viraemia among Dutch blood donors, 2003–2009. Vox Sang. 100, 261–266. https://doi.org/10.1111/j.1423-0410.2010.01423.x (2011).
Article CAS PubMed Google Scholar
Anderson, L. J. et al. Detection of antibodies and antigens of human parvovirus B19 by enzyme-linked immunosorbent assay. J. Clin. Microbiol. 24, 522–526. https://doi.org/10.1128/jcm.24.4.522-526.1986 (1986).
Article CAS PubMed PubMed Central Google Scholar
Cohen, B. J. & Buckley, M. M. The prevalence of antibody to human parvovirus B19 in England and Wales. J. Med. Microbiol. 25, 151–153. https://doi.org/10.1099/00222615-25-2-151 (1988).
Article CAS PubMed Google Scholar
Kelly, H. A. et al. The age-specific prevalence of human parvovirus immunity in Victoria Australia compared with other parts of the world. Epidemiol Infect 124, 449–457. https://doi.org/10.1017/s0950268899003817 (2000).
Article CAS PubMed PubMed Central Google Scholar
Tsujimura, M. et al. Human parvovirus B19 infection in blood donors. Vox. Sang. 69, 206–212. https://doi.org/10.1111/j.1423-0410.1995.tb02595.x (1995).
Article CAS PubMed Google Scholar
Qiu, J., Soderlund-Venermo, M. & Young, N. S. Human Parvoviruses. Clin. Microbiol. Rev. 30, 43–113. https://doi.org/10.1128/CMR.00040-16 (2017).
Article PubMed Google Scholar
White, D. G. et al. Human parvovirus arthropathy. Lancet 1, 419–421. https://doi.org/10.1016/s0140-6736(85)91145-6 (1985).
Article CAS PubMed Google Scholar
El-Mahallawy, H. A., Mansour, T., El-Din, S. E., Hafez, M. & Abd-el-Latif, S. Parvovirus B19 infection as a cause of anemia in pediatric acute lymphoblastic leukemia patients during maintenance chemotherapy. J. Pediatr. Hematol. Oncol. 26, 403–406. https://doi.org/10.1097/00043426-200407000-00001 (2004).
Article PubMed Google Scholar
Anderson, M. J., Lewis, E., Kidd, I. M., Hall, S. M. & Cohen, B. J. An outbreak of erythema infectiosum associated with human parvovirus infection. J. Hyg. (Lond) 93, 85–93. https://doi.org/10.1017/s0022172400060964 (1984).
Article CAS PubMed Google Scholar
Naides, S. J. Erythema infectiosum (fifth disease) occurrence in Iowa. Am. J. Public Health 78, 1230–1231. https://doi.org/10.2105/ajph.78.9.1230 (1988).
Article CAS PubMed PubMed Central Google Scholar
Plummer, F. A. et al. An erythema infectiosum-like illness caused by human parvovirus infection. N Engl. J. Med. 313, 74–79. https://doi.org/10.1056/NEJM198507113130203 (1985).
Article CAS PubMed Google Scholar
Yermalovich, M. A., Hubschen, J. M., Semeiko, G. V., Samoilovich, E. O. & Muller, C. P. Human parvovirus B19 surveillance in patients with rash and fever from Belarus. J. Med. Virol. 84, 973–978. https://doi.org/10.1002/jmv.23294 (2012).
Article PubMed Google Scholar
Rezaei, F. et al. Prevalence and genotypic characterization of Human Parvovirus B19 in children with measles- and rubella-like illness in Iran. J. Med. Virol. 88, 947–953. https://doi.org/10.1002/jmv.24425 (2016).
Article CAS PubMed Google Scholar
Los de Angeles Ribas, M.,Tejero, Y.,Cordero, Y.,Perez, D.,Sausy, A.,Muller, C. P., et al. Identification of human parvovirus B19 among measles and rubella suspected patients from Cuba. J. Med. Virol., 91:1351-1354. (2019) doi:https://doi.org/10.1002/jmv.25444.
Grolhier, C. et al. When a viral eruption hides another one: intrafamilial outbreak of parvovirus B19 and measles virus co-infections: case report. BMC Infect Dis. 20, 496. https://doi.org/10.1186/s12879-020-05183-4 (2020).
Article PubMed PubMed Central Google Scholar
Servant, A. et al. Genetic diversity within human erythroviruses: identification of three genotypes. J. Virol. 76, 9124–9134. https://doi.org/10.1128/jvi.76.18.9124-9134.2002 (2002).
Article CAS PubMed PubMed Central Google Scholar
Toan, N. L. et al. Phylogenetic analysis of human parvovirus B19, indicating two subgroups of genotype 1 in Vietnamese patients. J. Gen. Virol. 87, 2941–2949. https://doi.org/10.1099/vir.0.82037-0 (2006).
Article CAS PubMed Google Scholar
Parsyan, A., Szmaragd, C., Allain, J. P. & Candotti, D. Identification and genetic diversity of two human parvovirus B19 genotype 3 subtypes. J. Gen. Virol. 88, 428–431. https://doi.org/10.1099/vir.0.82496-0 (2007).
Article CAS PubMed Google Scholar
Li X, Lin Z, Liu J, Tang Y, Yuan X, Li N, et al. Overall prevalence of human parvovirus B19 among blood donors in mainland China: A PRISMA-compliant meta-analysis. Medicine (Baltimore), 99:e19832. (2020) doi:https://doi.org/10.1097/MD.0000000000019832.
Jia, J. et al. Existence of various human parvovirus B19 genotypes in Chinese plasma pools: genotype 1, genotype 3, putative intergenotypic recombinant variants and new genotypes. Virol. J. 13, 155. https://doi.org/10.1186/s12985-016-0611-6 (2016).
Article CAS PubMed PubMed Central Google Scholar
Cheng, P., Jian, Q., Fu, Z. & Ma, Y. Parvovirus B19-associated severe anemia in adult liver transplant recipients: a case series and review of the literature. Surg. Infect. (Larchmt) 23, 848–856. https://doi.org/10.1089/sur.2022.186 (2022).
Article PubMed Google Scholar
Huang, Q. et al. Parvovirus B19 infection in kidney transplant recipients: a prospective study in a teaching hospital in Shanghai, China. Transpl. Immunol. 74, 101667. https://doi.org/10.1016/j.trim.2022.101667 (2022).
Article CAS PubMed Google Scholar
Zhou, Y. et al. Detection of cytomegalovirus, human parvovirus B19, and herpes simplex virus-1/2 in women with first-trimester spontaneous abortions. J. Med. Virol. 87, 1749–1753. https://doi.org/10.1002/jmv.24218 (2015).
Article CAS PubMed Google Scholar
de Moraes, J. C. et al. Etiologies of rash and fever illnesses in Campinas Brazil. J. Infect. Dis. 204(Suppl 2), S627-636. https://doi.org/10.1093/infdis/jir490 (2011).
Article PubMed Google Scholar
Toshev, A., Ivanova, S., Kovaleva, V., Andonova, L. & Mihneva, Z. Detection of Human parvovirus B19 (HPVB19) in serum samples from fever-rash ill individuals during the rubella outbreak (2005) in Bulgaria. Biotechnol. Biotechnol. Equip. 28, 1103–1107. https://doi.org/10.1080/13102818.2014.967746 (2014).
Article CAS PubMed PubMed Central Google Scholar
Qiu, Q. et al. Viral pathogenic spectrum and epidemiological features of rash and fever syndrome in Shanxi province, China, during 2009–2015. Chinese J. Virol. 33, 169–175. https://doi.org/10.13242/j.cnki.bingduxuebao.003117 (2017).
Article Google Scholar
Wang, L. et al. The epidemiological characteristics of rash and fever syndrome in Gansu Province from 2009 to 2019. Chinese J. Dis. Control Prevent. 25, 1102–1106. https://doi.org/10.16462/j.cnki.zhjbkz.2021.09.019 (2021).
Article CAS Google Scholar
Ju, X., Xu, A., Fang, Q. & Huang, J. Etiology study on rash and fever illness in Guangdong Province, 2010–2012. Chinese J. Dis. Control and Prevent. 17, 670–673 (2013).
Google Scholar
Fei, Y. et al. Study on the etiology of rash and fever illness in Pudong New Area of Shanghai from 2010 to 2017. Chinese J. Dis. Control and Prevent. 23, 550–554. https://doi.org/10.16462/j.cnki.zhjbkz.2019.05.012 (2019).
Article Google Scholar
Chen M, Wang Y, Zhang Z, Zou L, Dou X, Huang F. Detection of rash and fever illness associated viruses in suspected measles cases in Beijing from 2017 to 2018. Int. J. Virol. (2019) 26:
Nguyen, Q. T., Wong, S., Heegaard, E. D. & Brown, K. E. Identification and characterization of a second novel human erythrovirus variant, A6. Virology 301, 374–380. https://doi.org/10.1006/viro.2002.1585 (2002).
Article CAS PubMed Google Scholar
Schneider, B. et al. Simultaneous persistence of multiple genome variants of human parvovirus B19. J. Gen. Virol. 89, 164–176. https://doi.org/10.1099/vir.0.83053-0 (2008).
Article CAS PubMed Google Scholar
Toppinen, M. et al. The landscape of persistent human DNA viruses in femoral bone. Forensic. Sci. Int. Genet. 48, 102353. https://doi.org/10.1016/j.fsigen.2020.102353 (2020).
Article CAS PubMed Google Scholar
Gallinella, G. Parvovirus B19 Achievements and Challenges. ISRN Virol. 2013, 1–33. https://doi.org/10.5402/2013/898730 (2013).
Article CAS Google Scholar
Servant-Delmas, A., Lefrere, J. J., Morinet, F. & Pillet, S. Advances in human B19 erythrovirus biology. J. Virol. 84, 9658–9665. https://doi.org/10.1128/JVI.00684-10 (2010).
Article CAS PubMed PubMed Central Google Scholar
Jain, A. & Kant, R. Genotypes of erythrovirus B19, their geographical distribution & circulation in cases with various clinical manifestations. Indian J. Med. Res. 147, 239–247. https://doi.org/10.4103/ijmr.IJMR_1816_16 (2018).
Article CAS PubMed PubMed Central Google Scholar
Sanabani, S., Neto, W. K., Pereira, J. & Sabino, E. C. Sequence variability of human erythroviruses present in bone marrow of Brazilian patients with various parvovirus B19-related hematological symptoms. J. Clin. Microbiol. 44, 604–606. https://doi.org/10.1128/JCM.44.2.604-606.2006 (2006).
Article CAS PubMed PubMed Central Google Scholar
Alves, A. D. et al. Persistence of Parvovirus B19 in liver from transplanted patients with acute liver failure. Future Microbiol. 15, 307–317. https://doi.org/10.2217/fmb-2019-0224 (2020).
Article CAS PubMed Google Scholar
Garcia Sde, O. et al. Parvovirus among patients with cytopenia of unknown origin in Brazil: a case-control study. J. Clin. Microbiol. 49, 1578–1580. https://doi.org/10.1128/JCM.00077-11 (2011).
Article PubMed Google Scholar
S. W, N.S. Y, K.E. B. Prevalence of parvovirus B19 in liver tissue: no association with fulminant hepatitis or hepatitis-associated aplastic anemia. J. Infect Dis.,187(10):1581-1586. (2003) doi:https://doi.org/10.1086/374781
Norja, P. et al. Bioportfolio: lifelong persistence of variant and prototypic erythrovirus DNA genomes in human tissue. Proc. Natl. Acad. Sci. U S A 103, 7450–7453. https://doi.org/10.1073/pnas.0602259103 (2006).
Article ADS CAS PubMed PubMed Central Google Scholar
Shackelton, L. A. & Holmes, E. C. Phylogenetic evidence for the rapid evolution of human B19 erythrovirus. J. Virol. 80, 3666–3669. https://doi.org/10.1128/JVI.80.7.3666-3669.2006 (2006).
Article CAS PubMed PubMed Central Google Scholar
Stamenkovic, G. G. et al. Substitution rate and natural selection in parvovirus B19. Sci. Rep. 6, 35759. https://doi.org/10.1038/srep35759 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Seetha, D. et al. Molecular-genetic characterization of human parvovirus B19 prevalent in Kerala State India. Virol J. 18, 96. https://doi.org/10.1186/s12985-021-01569-1 (2021).
Article CAS PubMed PubMed Central Google Scholar
Musiani, M. et al. Immunoreactivity against linear epitopes of parvovirus B19 structural proteins. Immunodominance of the amino-terminal half of the unique region of VP1. J. Med. Virol. 60, 347–352 (2000).
Article CAS PubMed Google Scholar
Cui, A. et al. Development of a multiplex one-step real-time RT-PCR assay for the simultaneous detection of eight viruses associated with febrile rash illnesses. Biosafety and Health 2, 89–94. https://doi.org/10.1016/j.bsheal.2020.04.003 (2020).
Article Google Scholar
Download references
We thank staff members of the FRI surveillance network laboratories and sentinel hospitals for assistance with field investigation, administration, data and specimen collection, and laboratory testing.
This work was supported by grants from the Key Technologies Research and Development Program of the National Ministry of Science (2009ZX10004202, 2013ZX10004202, 2018ZX10713002 and 2018ZX10713001-003).
National Health Commission (NHC) Key Laboratory of Medical Virology and Viral Diseases, WHO WPRO Regional Reference Measles/Rubella Laboratory, National Institute for Viral Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, 102206, People’s Republic of China
Haoran Jiang, Qi Qiu, Yangzi Zhou, Yan Zhang, Wenbo Xu & Aili Cui
School of Public Health and Management, Shandong First Medical University & Shandong Academy of Medical Sciences, Jinan, 250117, People’s Republic of China
Haoran Jiang & Xiaomei Li
Shanghai Municipal Center for Disease Control and Prevention, Institute of Infectious Disease Prevention and Control, 1380 Zhongshan West Road, Xuhui District, Shanghai, 200336, People’s Republic of China
Qi Qiu
NHC Key Laboratory of Systems Biology of Pathogens, Institute of Pathogen Biology, Chinese Academy of Medical Sciences & Peking Union Medical College, Beijing, 100730, People’s Republic of China
Yangzi Zhou
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
H.J. and A.C. prepared the manuscript; H.J., Q.Q., and Y.Z. performed the identification and sequencing; H.J., A.C., Y.Z. and X.L. analyzed the data; A.C., Q.Q. and H.J. designed the study; A.C. and W.X. acquired funding, and supervised the study. A.C. reviewed the manuscript. All authors read and approved the final manuscript.
Correspondence to Aili Cui or Xiaomei Li.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Jiang, H., Qiu, Q., Zhou, Y. et al. The epidemiological and genetic characteristics of human parvovirus B19 in patients with febrile rash illnesses in China. Sci Rep 13, 15913 (2023). https://doi.org/10.1038/s41598-023-43158-y
Download citation
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-023-43158-y
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.
Advertisement
© 2023 Springer Nature Limited
Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.