Open Access

Genome analysis revealed novel genotypes and recombination of the human parechoviruses prevalent in children in Eastern China

Gut Pathogens20168:52

https://doi.org/10.1186/s13099-016-0135-z

Received: 12 July 2016

Accepted: 19 October 2016

Published: 8 November 2016

Abstract

Background

Human parechovirus (HPeV) is a genus of virus in the family Picornaviridae, having two species A and B. HPeVs are common infectious agents, usually causing mild diarrhea and respiratory disease in young children.

Results

Here, we collected and sequenced the near complete genome sequences of 17 novel HPeVs from children with diarrhea in eastern China, which showed significant nucleotide sequence divergence. Phylogenetic analysis based on the complete genomes of these HPeV strains revealed that they belonged to seven different genotypes (including three putative novel genotypes), suggesting that HPeVs showed genotype diversity in this area. Recombination analysis indicated that one of the HPeVs, belonging to HPeV-1 with strain name 146Chzj02, was a recombinant generated by inter-genotype recombination among three HPeV strains belonging to three different genotypes, respectively.

Conclusion

Our data revealed the property of the genotype diversity of HPeVs prevalent in children with diarrhea in eastern China, which will be helpful in the future study of the viral evolution of HPeVs and the identification and typing of HPeVs in the clinical laboratory.

Keywords

Human parechovirus Genome sequence Phylogenetic analysis Recombination

Background

Human parechoviruses (HPeVs) belong to Parechovirus genus in the Picornaviridae family, which are non-enveloped, positive-sense RNA viruses with icosahedral capsids [1]. HpeVs have a genome of about 7300nt which flanked by an un-translated region (UTR) at both the 5′ and 3′ ends. The gnome encodes a polyprotein which is post-translationally cleaved by virus proteases to produce the structural (VP0, VP3 and VP1) and non-structural (2A–C and 3A–D) proteins [2].

HPeVs are frequent infectious agents, usually causing mild gastroenteritis and respiratory disease in young children, however, more serious cases, such as flaccid paralysis, encephalitis and myocarditis, have also been reported. HPeV is subdivided into 16 different genotypes [35]. The prototypes of HPeVs, including HPeV-1 and HPeV-2, were originally classified as echo-22 and echo-23, in enterovirus but were regrouped as a separate genus according to their genome organization and biological properties [6, 7]. Among these HPeV genotypes, some were associated with certain specific diseases. HPeV-1 is widely prevalent throughout the world and often found in children with diarrhea and gastroenteritis [3, 8]. HPeV-3 was particularly associated with sepsis syndromes and aseptic meningitis/meningoencephalitis in children [9, 10]. HPeV-4 was isolated from a 5 year old patient with lymphadenitis [11] and might also cause neonatal sepsis [12]. HPeV-6 was discovered from an infant with Reye syndrome [13]. HPeV-7 and HPeV-8 were identified from a patient with non-polio acute flaccid paralysis and a patient with enteritis, respectively [14]. HPeV-10 and -11 were identified in patients with acute gastroenteritis [15, 16]. HPeV-12 was detected in an infant with diarrhea and paralysis [17] while HPeV-14 was discovered from a febrile child [18].

In the present study, we sequenced 17 HPeVs genome sequences from the fecal samples from children <6 years of age with diarrhea, finding three putative new genotypes. Recombination analysis was also performed based on these genomes and one putative recombinant was found.

Methods

Samples and HPeVs detection

From Jan. 2011 to Jun. 2014, a total of 634 stool specimens which were collected from children under 6 years of age with diarrhea who were treated as outpatients were subjected to RT-PCR assay, where the degenerate primers were designed to amplify the VP3/VP1 region of HPeVs, which produced ~360 bp specific bands [19]. The positive PCR products were sequenced and then subjected to BLASTn in GenBank to find whether the HPeV strains belong to putative novel genotypes, where the sequences with <85% similarity to the best BLASTn matches in GenBank were considered to be putative novel genotypes and selected for genome sequencing.

Genome sequencing

RT-PCRs were used to acquire the genome. Primers were designed based on the complete genomes in GenBank which showed high similarities to the 360 bp sequences got in the RT-PCR screening in this study. Briefly, the total RNA was extracted from 200 μl fecal supernatants using QIAamp Viral RNA Mini Kit (Qiagen, Germany) in accordance with the manufacturer’s protocol. RT-PCR was performed by using TaKaRa RNA PCR kit (TaKaRa, Japan). The parameters of PCR amplification included an initial incubation at 94 °C for 5 min, followed by 30 cycles of denaturation at 94 °C for 40 s, annealing for 40 s at a temperature varied depending on different primers, and extension at 72 °C for 1 min, with a final incubation at 72 °C for 10 min. The specific bands were excised from the gel and sequenced by Sanger method (Shanghai Sangon, China).

Phylogenetic analysis

To investigate the relationship between the HPeVs in this study and those with complete genomes in GenBank, all the HPeV genomes were retirved from GenBank and phylogenetic analysis was performed based on the complete genomes. The HPeVs genome sequences from GenBank were pre-analyzed and those sequences showing less than 1% divergence from each other were considered as the same strain and only one of them were included in the phylogenetic analysis. Sequence alignment was performed using CLUSTAL W with the default settings [20]. A phylogenetic tree with 1000 bootstrap resamples of the alignment data sets was generated using the neighbor-joining method in MEGA5.0 [21]. Sequence divergence based on the complete genome of HPeVs was analyzed using AliGROOVE software [22]. The scores were ranging from −1, indicating full random similarity, to +1, non-random similarity.

Virus recombination detection

Detection of potential recombinant sequences, identification of potential parental sequences, and localization of possible recombination break points were determined using the Recombination Detection Program (RDP3.0) [23]. The potentially significant recombination events were further confirmed by SimPlot software [24] and phylogenetic analysis based on the exchanged sequences of the recombinants.

Results

Phylogenetic characterization of the novel HPeVs

RT-PCR screening indicated 83 (13.1%) of the total 634 stool specimens were positive for HPeVs genes. Sequence analysis indicated that 17 sequences showed <85 % similarity to the previous HPeV strains in GenBank, which were then selected for further genome sequencing. In order to characterize these strains, the nearly complete genomes of the 17 putative novel HPeV strains were sequenced and submitted to GenBank with GenBank nos. KT879915–KT879931. The heterogeneity of sequence divergence based on the complete genomes of the 17 HPeVs was indicated in Fig. 1a, which showed that the 17 shared 71–88% identities among themselves, suggesting these HPeVs showed the property genotype diversity. Phylogenetic analysis based on the complete genome of the 17 HPeVs in the present study and those HPeVs with complete genome available in GenBank indicated that the 17 HPeVs identified in this study belonged to 7 different genotypes (Fig. 1b). Seven of them belonged to HPeV-1, being grouped into 4 different clusters within the HPeV-1 clade. One strain clustered with an HPeV-2 strain (HM996978), sharing 82.2% sequence similarity. Two strains clustered together within the HPeV-5 clade, sharing 80.1% sequence identity with each other. One strain belonged to HPeV-6, clustering and sharing 85% sequence identity the other two previous HPeV-6 strains in the same clade. Five of the strains clustered together into a separate group with 74–78% sequence similarities with all the other HPeVs, suggesting they belonged to a new genotype (HPeV-17). One strain clusters separately, sharing 75–83% sequence similarities to all the other HPeVs, representing a new genotype (HPeV-18). The rest one strain formed the deepest branch, sharing 73–77% sequence identities with all the other HPeVs strains, and belonged to a putative new genotype, HPeV-19.
Fig. 1

Sequence divergence and phylogenetic analysis based on the complete genome of HPeVs in the present study. a Nucleotide sequence divergence of the 17 HPeVs in this study was established using AliGROOVE software. The scores were ranging from 0 to +1 indicating the non-random similarity. The darker blue indicated the higher non-randomized accordancy between pairwise sequence comparisons. Strain names in this study were listed on top and the right hand side of the matrix; b neighbor joining phylogenetic tree was constructed with Mega5.0 from multiple alignments of the complete genome of the 17 HPeVs in the present study and other 47 HPeVs with complete genome sequence available in GenBank. Two Ljungan virus strains were included as outgroup. Bootstrap values are indicated at the nodes. The scale bar indicates the number of substitutions per position for a unit branch length. Viruses identified in the present study were displayed with diamonds

Recombination of HPeVs

One potentially significant recombination event was found with a high degree of confidence. Figure 2a, b indicated the bootscan and similarity plots, respectively, of the recombination event which belonged to an intra-genotype recombination occurred between the major parent of an HPeV-1 strain (FM178558) and two minor parents including an HPeV-6 strain (AB252582) and an HPeV-3 strain (JX682576), leading to the recombinant 146Chzj02. Figure 1a, b displayed that this recombination event occurred within the P3 gene of HPeVs. Among the four strains related with this recombination event, FM178558 belongs to HPeV-1 which was isolated from human feces from Netherland [25]. AB252582 was isolated from cerebrospinal fluid from a patient in Japan [13], while JX682576 was isolated from feces of a newborn with gastroenteritis in Germany [26]. To confirm the recombination event, the relevant strains were analyzed by neighbor joining trees using MEGA5.0. Figure 2c, d indicated the trees constructed on the recombinant region and the non-recombinant regions, respectively. The recombinant strain 146Chzj02 clustered closely with its major and minor parental strains in the two phylogenetic trees, respectively. Taken together, our data clearly indicated the existence of the putative recombination event among AB252582, JX682576, and FM178558, which belonged to an inter-genotype recombination.
Fig. 2

Identification of the recombination event based on complete genome of the potential recombinant 146Chzj02. a BOOTSCAN evidence for the recombination origin on the basis of pairwise distance, modeled with a window size 200, step size 20, and 100 Bootstrap replicates using RDP3.0 software; b similarity analysis of complete genomes of potential recombinant 146Chzj02 and its three potential parental strains using Simplot Version 3.5.1 software by a sliding window of 200 nucleotides moving in steps of 20 nucleotides. The shadowy area indicated the recombination event. c Neighbor joining tree constructed using the recombinant region (5694–6798nt); d neighbor joining tree constructed using the non-recombinant region. The nucleotide position is indicated based on the genome sequence of 146Chzj02

To verify the whether recombinant is naturally present in the fecal sample, two sets of specific primers were designed basing on the genome sequence of 146Chzj02, which amplified two fragments including the two recombination breakpoints, respectively. PCR amplified specific bands were sequenced and sequencing results of the two types of PCR amplicons were identical to the originally determined genome, suggesting that 146Chzj02 is natural recombinant.

Discussion

Although many reports indicated that HPeVs can cause infection in young children and lead to various diseases including diarrhea, respiratory diseases, encephalomyelitis, meningitis, myocarditis, lymphadenopathy, hemolytic uremic syndrome, and sudden infant death syndrome, no causal relationship has been established except with HPeV-3 being primarily associated with neonatal sepsis [27, 28]. In the present study, we sequenced the complete genomes of 17 HPeV strains from children <6 years old with diarrhea in eastern China. Up to now, there are totally 16 different types of HPeV, of which type 6 to 16 were identified as genotype by genetic analysis based on sequence divergence [6, 7]. However, we did not include a control group in our study, which provide no sufficient proof of HPeV association to the diarrhea of these children included in this study. Besides, mere detection of HPeV in clinical samples does not confirm or establish its role in diarrhea, as even the asymptomatic individuals may continue to shed HPeVs for longer periods of time [29]. Phylogenetic analysis based on the complete genome sequences of 17 HPeV strains indicated they belonged to 7 different types, which suggested that HPeVs in eastern China showed property of genotype diversity. HPeV1 is the most frequently detected type worldwide, which is also the major genotype in the present study. Considering such a vast diversity of HPeVs prevalent in children population in eastern China, a highly sensitive RT-PCR assay should be designed covering more genotypes of HPeVs in laboratory diagnosis.

Recombination is a rather common phenomenon in picornaviruses [30] and understanding recombination will be helpful in study of the evolution of pathogens and the development of viral vaccine. Like other picornaviruses, recombination in HPeVs has also been found, with breaking points at the capsid encoding and non-structural encoding regions [3134], among these recombination events most belonged to intra-genotype recombination. In the present study, a recombination event was detected in the non-structural region, which was inter-genotype recombination occurred between 3 HPeV strains belong three different genotypes and led to an HPeV-1 recombinant. This finding would support that the most frequent breaking points for recombination in HPeVs flanked the capsid-encoding region [2, 3, 25]. Our data indicated that the recombinant (146Chzj02) and its three parental strains (including AB252582, JX682576, and FM178558) were from four different countries, respectively, which implies that this recombination event might have occurred among the ancestors of the three parental HPeV strains long time ago, when they were simultaneously circulated in a single area in that meantime.

Conclusion

We sequenced 17 complete genome sequences of HPeV from the feces of Children <6 years old with diarrhea in eastern China, which belonged to seven different genotypes including three putative novel genotypes based on phylogenetic analysis. Recombination analysis indicated that one of the HPeVs in the current study belonged to putative recombinant produced by inter-genotype recombination occurred among 3 different genotypes of HPeVs.

Our data provided the genetic diversity of HPeVs from children with diarrhea in eastern China, which will be helpful in the future study of the viral evolution of HPeVs and the identification and typing of HPeVs in the clinical laboratory.

Abbreviations

HPeV: 

human parechovirus

PCR: 

polymerase chain reaction

UTR: 

un-translated region

Declarations

Authors’ contributions

XZ conceived the study. XZ, YS, and YX performed all the experiments. XZ wrote the paper. All authors read and approved the final manuscript.

Acknowledgements

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All the data supporting the results are included in the article.

Ethics approval and consent to participate

This study was approved by the ethics committee of Nanjing Lishui People’s Hospital with Reference no. 20110023. Informed consent was collected from all patients’ parents or guardians.

Funding

This work was supported by the National Natural Sciences Foundation of China No. 31302107. The funding body did not participate in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)
Nanjing Lishui People’s Hospital

References

  1. Joki-Korpela P, Hyypiä T. Parechoviruses, a novel group of human picornaviruses. Ann Med. 2001;33:466–71.View ArticlePubMedGoogle Scholar
  2. Williams CH, Panayiotou M, Girling GD, Peard CI, Oikarinen S, Hyöty H, Stanway G. Evolution and conservation in human parechovirus genomes. The Journal of general virology. 2009;90(Pt 7):1702–12.View ArticlePubMedGoogle Scholar
  3. Benschop KSM, de Vries M, Minnaar RP, Stanway G, van der Hoek L, Wolthers KC, Simmonds P. Comprehensive full-length sequence analyses of human parechoviruses: diversity and recombination. J Gen Virol. 2010;91(Pt 1):145–54.View ArticlePubMedGoogle Scholar
  4. Al-Sunaidi M, Williams CH, Hughes PJ, Schnurr DP, Stanway G. Analysis of a new human parechovirus allows the definition of parechovirus types and the identification of RNA structural domains. J Virol. 2007;81:1013–21.View ArticlePubMedGoogle Scholar
  5. Chieochansin T, Vichiwattana P, Korkong S, Theamboonlers A, Poovorawan Y. Molecular epidemiology, genome characterization, and recombination event of human parechovirus. Virology. 2011;421:159–66.View ArticlePubMedGoogle Scholar
  6. Hyypiä T, Horsnell C, Maaronen M, Khan M, Kalkkinen N, Auvinen P, Kinnunen L, Stanway G. A distinct picornavirus group identified by sequence analysis. Proc Natl Acad Sci USA. 1992;89:8847–51.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Wigand R, Sabin AB. Properties of ECHO types 22, 23 and 24 viruses. Archiv für die gesamte Virusforschung. 1961;11:224–47.View ArticlePubMedGoogle Scholar
  8. Ito M. Isolation and identification of a novel human parechovirus. J Gen Virol. 2004;85:391–8.View ArticlePubMedGoogle Scholar
  9. Renaud C, Harrison CJ. Human parechovirus 3: the most common viral cause of meningoencephalitis in young infants. Infect Dis Clin North Am. 2015;29:415–28.View ArticlePubMedGoogle Scholar
  10. Bissel SJ, Auer RN, Chiang C-H, Kofler J, Murdoch GH, Nix WA, Painter M, Richer M, Sartelet H, Wang G, Wiley CA. Human parechovirus 3 meningitis and fatal leukoencephalopathy. J Neuropathol Exp Neurol. 2015;74:767–77.View ArticlePubMedGoogle Scholar
  11. Benschop KSM, Schinkel J, Luken ME, van den Broek PJM, Beersma MFC, Menelik N, van Eijk HWM, Zaaijer HL, VandenBroucke-Grauls CMJE, Beld MGHM, Wolthers KC. Fourth human parechovirus serotype. Emerg Infect Dis. 2006;12:1572–5.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Jääskeläinen AJ, Kolehmainen P, Kallio-Kokko H, Nieminen T, Koskiniemi M, Tauriainen S, Lappalainen M. First two cases of neonatal human parechovirus 4 infection with manifestation of suspected sepsis, Finland. J Clin Virol. 2013;58:328–30.View ArticlePubMedGoogle Scholar
  13. Watanabe K, Oie M, Higuchi M, Nishikawa M, Fujii M. Isolation and characterization of novel human parechovirus from clinical samples. Emerg Infect Dis. 2007;13:889–95.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Li L, Victoria J, Kapoor A, Naeem A, Shaukat S, Sharif S, Alam MM, Angez M, Zaidi SZ, Delwart E. Genomic characterization of novel human parechovirus type. Emerg Infect Dis. 2009;15:288–91.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Pham NTK, Trinh QD, Takanashi S, Abeysekera C, Abeygunawardene A, Shimizu H, Khamrin P, Okitsu S, Mizuguchi M, Ushijima H. Novel human parechovirus. Sri Lanka. Emerging Infectious Diseases. 2010;16:130–2.View ArticleGoogle Scholar
  16. Pham NTK, Takanashi S, Tran DN, Trinh QD, Abeysekera C, Abeygunawardene A, Khamrin P, Okitsu S, Shimizu H, Mizuguchi M, Ushijima H. Human parechovirus infection in children hospitalized with acute gastroenteritis in Sri Lanka. J Clin Microbiol. 2011;49:364–6.View ArticlePubMedGoogle Scholar
  17. Alam MM, Khurshid A, Shaukat S, Sharif S, Rana MS, Angez M, Naeem M, Zaidi SSZ. Identification of human parechovirus genotype, HPeV-12, in a paralytic child with diarrhea. J Clin Virol. 2012;55:339–42.View ArticlePubMedGoogle Scholar
  18. Benschop K, Thomas X, Serpenti C, Molenkamp R, Wolthers K. High prevalence of human parechovirus (HPeV) genotypes in the Amsterdam region and identification of specific HPeV variants by direct genotyping of stool samples. J Clin Microbiol. 2008;46:3965–70.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Harvala H, Robertson I, McWilliam Leitch EC, Benschop K, Wolthers KC, Templeton K, Simmonds P. Epidemiology and clinical associations of human parechovirus respiratory infections. J Clin Microbiol. 2008;46:3446–53.View ArticlePubMedPubMed CentralGoogle Scholar
  20. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG. Clustal W and Clustal X version 2.0. Bioinform. 2007;23:2947–8.View ArticleGoogle Scholar
  21. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Kück P, Meid SA, Groß C, Wägele JW, Misof B. AliGROOVE–visualization of heterogeneous sequence divergence within multiple sequence alignments and detection of inflated branch support. BMC Bioinformatics. 2014;15:294.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P. RDP3: a flexible and fast computer program for analyzing recombination. Bioinform. 2010;26:2462–3.View ArticleGoogle Scholar
  24. Salminen MO, Carr JK, Burke DS, McCutchan FE. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res Hum Retroviruses. 1995;11:1423–5.View ArticlePubMedGoogle Scholar
  25. Zoll J, Galama JMD, van Kuppeveld FJM. Identification of potential recombination breakpoints in human parechoviruses. J Virol. 2009;83:3379–83.View ArticlePubMedPubMed CentralGoogle Scholar
  26. Eis-Hübinger AM, Eckerle I, Helmer A, Reber U, Dresbach T, Buderus S, Drosten C, Müller A. Two cases of sepsis-like illness in infants caused by human parechovirus traced back to elder siblings with mild gastroenteritis and respiratory symptoms. J Clin Microbiol. 2013;51:715–8.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Boivin G, Abed Y, Boucher FD. Human parechovirus 3 and neonatal infections. Emerg Infect Dis. 2005;11:103–5.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Levorson RE, Jantausch BA, Wiedermann BL, Spiegel HML, Campos JM. Human parechovirus-3 infection: emerging pathogen in neonatal sepsis. Pediatr Infect Dis J. 2009;28:545–7.View ArticlePubMedGoogle Scholar
  29. Kolehmainen P, Oikarinen S, Koskiniemi M, Simell O, Ilonen J, Knip M, Hyöty H, Tauriainen S. Human parechoviruses are frequently detected in stool of healthy Finnish children. J Clin Virol. 2012;54:156–61.View ArticlePubMedGoogle Scholar
  30. Lukashev AN. Recombination among picornaviruses. Rev Med Virol. 2010;20:327–37.View ArticlePubMedGoogle Scholar
  31. Chen H, Zheng X, Chen X, Shi T, Yao Y, Yuan Q, Chen Q, Yu S. Diversity and recombination of human parechovirus in children with acute gastroenteritis in Guangzhou, China. J Med Virol. 2015;87:296–302.View ArticlePubMedGoogle Scholar
  32. Chang JT, Yang CS, Chen YS, Chen BC, Chiang AJ, Chang YH, Tsai WL, Lin YS, Chao D, Chang TH. Genome and infection characteristics of human parechovirus type 1: the interplay between viral infection and type I interferon antiviral system. PLoS ONE. 2015;10:e0116158.View ArticlePubMedPubMed CentralGoogle Scholar
  33. Zhu R, Luo L, Zhao L, Deng J, Wang F, Sun Y, Song Q, Ding Y, Qian Y. Characteristics of the mosaic genome of a human parechovirus type 1 strain isolated from an infant with pneumonia in China. Infect Genet Evol. 2015;29:91–8.View ArticlePubMedGoogle Scholar
  34. Thoi TC, Than VT, Kim W. Whole genomic characterization of a Korean human parechovirus type 1 (HPeV1) identifies recombination events. J Med Virol. 2014;86:2084–91.View ArticlePubMedGoogle Scholar

Copyright

© The Author(s) 2016