key: cord-0899094-usl8z6ep
authors: Zheng, Wen-zhi; Wei, Tian-li; Ma, Fen-lian; Yuan, Wu-mei; Zhang, Qian; Zhang, Ya-xin; Cui, Hong; Zheng, Li-shu
title: Human polyomavirus type six in respiratory samples from hospitalized children with respiratory tract infections in Beijing, China
date: 2015-10-13
journal: Virol J
DOI: 10.1186/s12985-015-0390-5
sha: 814ec1dc13c5f446a41d1b3f1becb6495897455a
doc_id: 899094
cord_uid: usl8z6ep

BACKGROUND: HPyV6 is a novel human polyomavirus (HPyV), and neither its natural history nor its prevalence in human disease is well known. Therefore, the epidemiology and phylogenetic status of HPyV6 must be systematically characterized. METHODS: The VP1 gene of HPyV6 was detected with an established TaqMan real-time PCR from nasopharyngeal aspirate specimens collected from hospitalized children with respiratory tract infections. The HPyV6-positive specimens were screened for other common respiratory viruses with real-time PCR assays. RESULTS: The prevalence of HPyV6 was 1.7 % (15/887), and children ≤ 5 years of age accounted for 80 % (12/15) of cases. All 15 HPyV6-positive patients were coinfected with other respiratory viruses, of which influenza virus A (IFVA) (8/15, 53.3 %) and respiratory syncytial virus (7/15, 46.7 %) were most common. All 15 HPyV6-positive patients were diagnosed with lower respiratory tract infections, and their viral loads ranged from 1.38 to 182.42 copies/μl nasopharyngeal aspirate specimen. The most common symptoms were cough (100 %) and fever (86.7 %). The complete 4926-bp genome (BJ376 strain, GenBank accession number KM387421) was amplified and showed 100 % identity to HPyV6 strain 607a. CONCLUSIONS: The prevalence of HPyV6 was 1.7 % in nasopharyngeal aspirate specimens from hospitalized children with respiratory tract infections, as analyzed by real-time PCR. Because the coinfection rate was high and the viral load low, it was not possible to establish a correlation between HPyV6 and respiratory diseases. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s12985-015-0390-5) contains supplementary material, which is available to authorized users.

The family Polyomaviridae contains viruses that are among the smallest known to infect humans, in terms of both their particle sizes and genome lengths. The 40-45 nm nonenveloped icosahedral particles carry a circular doublestranded DNA genome of approximately 5000 bp, which is divided into three functional regions: the noncoding control region, the early gene region encoding the large T antigen and the small T antigen, and the late coding region, encoding capsid proteins VP1, VP2, and VP3 [1] .

Human polyomaviruses (HPyVs) have not been associated with any severe acute disease in healthy humans. However, the viral proteins expressed by polyomaviruses (PyVs) can initiate the transformation and immortalization of cultured cells and cause cancer in experimental animals [2] . The HPyV family currently consists of 13 members, including JCPyV [3] , BKPyV [4] , WUPyV [5] , KIPyV [6] , MCPyV [7] , HPyV6 [8] , HPyV7 [8] , TSPyV [9] , HPyV9 [10] , HPyV10 [11] [12] [13] , STLPyV [14] , HPyV12 [15] and NJPyV-2013 [16] .

Previous studies have indicated that a number of HPyVs are associated with human diseases, such as progressive multifocal leukoencephalopathy (JCPyV), hemorrhagic cystitis (BKPyV), Merkel cell carcinoma (MCPyV), and trichodysplasia spinulosa (TSPyV) [3, 7, 9, [17] [18] [19] . A serological study demonstrated that the seroprevalence of these viruses in 1501 plasma samples from healthy adult blood donors was BKPyV (82 %), JCPyV (39 %), KIPyV (55 %), WUPyV (69 %), MCPyV isolate 350 (25 %); and MCPyV isolate 339 (42 %). The seroprevalence of all polyomaviruses in children under 21 years of age (n = 721) was similar to that in the adult population, suggesting that primary exposure to these viruses occurs in early childhood and seems to result in lifelong persistence [20] . Nevertheless, the natural histories of most HPyVs and their prevalence in human diseases are not yet well known.

In 2008, MCPyV was discovered by Feng et al. [7] . HPyV6 and HPyV7 were discovered with MCPyV in skin swabs from the foreheads of healthy volunteers [8] . However, later research could not demonstrate a relationship between HPyV6 and Merkel cell carcinoma or other skin diseases. A phylogenetic analysis of the complete HPyV6 genome indicated that HPyV6 shared a branch with KIPyV and WUPyV. Previous reports have shown that genomic fragments of KIPyV and WUPyV have been regularly detected in nasopharyngeal aspirates of children with respiratory tract infections (RTIs) and are suspected of a causal relationship with respiratory disease. However, the link between these PyVs and respiratory diseases remains speculative [5, 6, 21, 22] .

A real-time PCR assay was established here to determine the prevalence of HPyV6 throughout a time period of 12 months (from October 2011 to September 2012). The prevalence of HPyV6 was 1.7 % (15/887). A complete HPyV6 genome was amplified, sequenced and found to be identical with a HPyV6 isolate from USA.

With the standard curve derived from serial DNA dilutions, the dynamic range of the real-time PCR assay was 10 0 -10 10 copies/μl and the limit of detection was one copy. The coefficient of determination (R 2 = 0.99667) showed a good linear correlation. The TaqMan-based real-time PCR assay to detect HPyV6 did not amplify any other viral pathogen, showing the excellent specificity of this assay. Four different DNA concentrations (10 5 -10 8 copies per reaction) were repeated five times in each run. The maximum coefficient of variation was 0.66 %, which indicates good precision (data not shown).

A total of 887 NPA samples were obtained from 887 children with RTI. The sex ratio (male:female) was 524:363 (1.44:1) and the median age was 24 months (ranging from 3 days to 14 years). The prevalence of HPyV6 was 1.7 % (15/887) in the 887 NPA specimens tested. Of the 15 HPyV6-positive patients, eight were male (8/524, 1.53 %) and seven were female (7/363, 1.93 %), so the prevalence was similar in both sexes (p > 0.05). The ages of the infected patients ranged from 4 days to 13 years, and children ≤ 5 years of age accounted for 80 % (12/15) of the total HPyV6-positive children. The age distribution of the HPyV6-infected children indicated that those aged 37-48 months had the highest infection rate Table 1) .

All 15 HPyV6-positive patients were diagnosed with lower RTI, including bronchopneumonia in nine (60 %), acute bronchitis in three (20 %), bronchitis in two (13.3 %), and pneumonia in one (6.7 %). The most common symptom was cough, which occurred in all 15 patients (100 %). Other clinical presentations included fever (n = 13, 86.7 %), gasping (n = 3, 20 %), vomiting (n = 3, 20 %), and diarrhea (n = 1, 6.7 %). Among the HPyV6-positive specimens, the HPyV6 genome copies ranged from 1.38 to 182.42 copies/μl NPA on a realtime PCR assay ( Table 2 ). The HPyV6 genome copy number was 52.07 copies/μl NPA in children suffering lower RTIs only, and was slightly higher than 29.75 copies/μl NPA in those infected with lower RTIs and other diseases, but these did not differ significantly (p > 0.05, Mann-Whitney U test).

To determine the complete HPyV6 genomic sequence, overlapping genomic fragments were amplified with nested PCR using 15 pairs of virus-specific primers (listed in Additional file 1: Table S2 ), and the complete 4926-bp genome was compiled with BioEdit 9.0 software. The genome sequence of HPyV6 was deposited in GenBank under accession number KM387421. BLAST analysis with the complete HPyV6 genome showed a high level of nucleic acid identity to the six full HPyV6 genomes available in GenBank. Relationsship to the HPyV6 strain 607a was 100 %.

HPyV6, thought to be a skin-tropic polyomavirus, was initially described in 2010. Since then, subgenomic fragments of HPyV6 DNA have been detected in a variety of specimen types, including skin, respiratory secretion samples, and various tumor samples (for instance, The answers to these questions will require more information on the biology and epidemiology of the HPyVs. The genomes and proteins of HPyV6, one of the novel HPyVs, show little sequence homology with previously reported HPyVs (BKPyV and JCPyV). Although HPyV6 encodes a conserved, potentially carcinogenic LTAg, previous studies have shown no association between HPyV6 and tumors [1] . Furthermore, a phylogenetic analysis indicated that LTAg of HPyV6 (KM387421) is only distantly related to its homologues in other cancerassociated HPyVs. The HPyV6, WUPyV, and KIPyV strains formed a clade in the complete genome and VP1 amino acid phylogenies, whether HPyV6 also associate with respiratory infection which need more clinical and experimental evidences to support. HPyV6 has been detected in specimens from the human respiratory tract, but there are as yet insufficient epidemiological data to demonstrate a correlation between HPyV6 and respiratory disease. Because initial infections with most HPyVs occur in infancy, the prevalence of HPyV6 in NPAs from children was detected with real-time PCR.

HPyV6 displayed an overall prevalence of 1.7 % in NPA samples collected from children in a hospital in China, which is similar to its prevalence reported previously (0.5-2 %) [34, 35] . It has not been confirmed that HPyV6 infects humans via the respiratory tract, but the respiratory tract may be a possible route of transmission. In this study, HPyV6 was mainly detected in children less than 5 years of age, and the peak incidence occurred in spring. All 15 HPyV6-positive patients were coinfected with other respiratory viruses, of which IFVA and RSV were the most common. The HPyV6-positive patients were diagnosed with lower RTIs, 60 % had bronchopneumonia, and the most common symptoms were cough and fever.

Although known HPyVs cause disease in patients with immune-system imbalances, they do not seem to cause obvious illnesses in the great majority of infected individuals. However, BKPyV can induce nephropathy in kidney transplantation recipients and JCPyV causes progressive multifocal leukoencephalopathy. Progressively increasing or high viral loads are also associated with high-level viral replication and disease. For instance, progressive multifocal leukoencephalopathy and hemorrhagic cystitis are related to high viral loads of JCPyV and BKPyV, respectively. The present study cannot confirm that HPyV6 is the cause of RTIs in hospitalized children, because the viral loads of HPyV6 were low (1.38-182.42 copies/μl) and the coinfection rate with other respiratory viruses was high. In previous reports, HPyVs were detected in the respiratory tract and skin but, to date, there has been insufficient evidence that HPyV6 is associated with any respiratory tract disease or skin disease [36] . Whether HPyV6 induce any disease requires the analysis of further data.

The detection rate for HPyV6 by real-time PCR assay was 1.7 % in 887 NPA samples collected from hospitalized children with RTI. An association between HPyV6 and respiratory diseases could not been revealed due to the high coinfection rate and the low HPyV6 viral load.

From 1 October, 2011, to 30 September, 2012, 887 nasopharyngeal aspirate (NPA) specimens were collected continuously from 887 hospitalized children with RTI, who ranged in age from 3 days to 14 years, at Beijing Friendship Hospital, Beijing, China. The patients' parents or guardians gave their written informed consent for specimen collection and testing, and the project was approved by the Ethical Committee of the Beijing Friendship Hospital. The NPA specimens were collected and transported immediately to the National Institute for Viral Disease Control and Prevention, China CDC, and stored at −80°C until further processing. All demographic data and clinical findings were recorded on a standard form. Total nucleic acids were extracted from each NPA specimen using the QIAamp MinElute Virus Spin Kit (Qiagen, Beijing, China).

A TaqMan-based real-time PCR assay for the detection of VP1 gene of HPyV6 was designed with Primer Express 3.0 software. The primer sequences used were HPyV6-F 5'-TTAACACCCTTCTTTGTGCTGCTA-3' and HPyV6-R 5'-GCCCAATTATTCAAAGCAGCTAA-3' , and the probe sequence was HPyV6-P FAM-CTGTCACAG GCCTGCTGAGCAATAGATTTC-TAMRA. The specificities of the primers and probe were evaluated in GenBank with BLAST. The primers and probe were synthesized by Invitrogen (Beijing, China). A common reaction mix was prepared for the real-time PCR assays. Briefly, the final 20 μl reaction mix contained 10 μl of TaqMan Gene Expression Master Mix, 1.8 μl of each primer (10 pmol/μl), 0.2 μl of probe (10 pmol/μl), 2 μl of pMD18-T/VP1 plasmid template (plasmid pMD18-T linked to the VP1 gene of HPyV6), and 4.2 μl of H 2 O. The amplification conditions included an initial incubation at 50°C for 2 min and 95°C for 15 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min, using the Mx3005P qPCR System (Agilent Stratagene). Ten-fold serial dilutions of the pMD18-T/VP1 plasmid (from 10 0 to 10 10 copies/μl) were added to the real-time PCR reactions in duplicate. The results were used to generate a standard curve for HPyV6. Specificity was assessed by testing mixed samples of other common HPyVs, including WUPyV, KIPyV, JCPyV, BKPyV, and HPyV7. To test the reproducibility of the assay, we added 10 5 -10 8 copies/μl of pMD18-T/VP1 plasmid to each reaction and each concentration of DNA was repeated five times. In addition, housekeeping gene glyceraldehyd-3-phosphate dehydrogenase (GAPDH) was used as internal control. 2 μl of nucleic acid of each NPA specimen was added to each reaction. Only samples that were positive according to both PCR and DNA sequencing were considered "positive".

Fifteen overlapping fragments of the complete genome of HPyV6 strain BJ376 were PCR amplified (Table 2) with the Takara Ex Taq kit. The 15 overlapping fragments were then cloned into pMD18-T and sequenced (Invitrogen). The nucleotide sequence of the full-length HPyV6 genome was then compiled using BioEdit 9.0 software. The full-length HPyV6 sequence was aligned with the sequences of other HPyVs and other HPyV6 strains available in GenBank with DNAStar software. A neighbor-joining tree was constructed with MEGA 6.0.

The significance of differences between the prevalence rates and viral loads of various groups was tested with Fisher's exact test and the Mann-Whitney U test. All analyses were performed with SPSS 19.0 software. 

Genome analysis of the new human polyomaviruses

Human polyomaviruses in disease and cancer

Cultivation of papova-like virus from human brain with progressive multifocal leucoencephalopathy

New human papovavirus (B.K.) isolated from urine after renal transplantation

Identification of a novel polyomavirus from patients with acute respiratory tract infections

Identification of a third human polyomavirus

Clonal integration of a polyomavirus in human Merkel cell carcinoma

Merkel cell polyomavirus and two previously unknown polyomaviruses are chronically shed from human skin

Discovery of a new human polyomavirus associated with trichodysplasia spinulosa in an immunocompromized patient

A novel human polyomavirus closely related to the african green monkeyderived lymphotropic polyomavirus

Complete genome sequence of a tenth human polyomavirus

Discovery of a novel polyomavirus in acute diarrheal samples from children

Identification of MW polyomavirus, a novel polyomavirus in human stool

Discovery of STL polyomavirus, a polyomavirus of ancestral recombinant origin that encodes a unique T antigen by alternative splicing

Identification of a novel human polyomavirus in organs of the gastrointestinal tract

Identification of a novel polyomavirus in a pancreatic transplant recipient with retinal blindness and vasculitic myopathy

Association between a high BK virus load in urine samples of patients with graftversus-host disease and development of hemorrhagic cystitis after hematopoietic stem cell transplantation

Progressive multifocal leukoencephalopathy and promyelocytic leukemia nuclear bodies: a review of clinical, neuropathological, and virological aspects of JC virus-induced demyelinating disease

SV40 and human cancer: a review of recent data

Seroepidemiology of human polyomaviruses

Presence of the newly discovered human polyomaviruses KI and WU in Australian patients with acute respiratory tract infection

Molecular epidemiology of KI and WU polyomaviruses in infants with acute respiratory disease and in adult hematopoietic stem cell transplant recipients

Rapid and sensitive method using multiplex real-time PCR for diagnosis of infections by influenza A and influenza B viruses, respiratory syncytial virus, and parainfluenza viruses 1, 2, 3, and 4

Development and Assay of RNA Transcripts of Enterovirus Species A to D, Rhinovirus Species A to C, and Human Parechovirus:Assessment of Assay Sensitivity and Specificity of Real-Time Screening and Typing Methods

WU and KI polyomaviruses in respiratory samples from allogeneic hematopoietic cell transplant recipients

Development and evaluation of real-time PCR assays for the detection of the newly identified KI and WU polyomaviruses

Human bocavirus in patients with respiratory tract infection

Pring-Akerblom P. Rapid and quantitative detection of human adenovirus DNA by real-time PCR

Real-time reverse transcriptase PCR assay for detection of human metapneumoviruses from all known genetic lineages

Design and performance testing of quantitative real time PCR assays for influenza A and B viral load measurement

Impact of human coronavirus infections in otherwise healthy children who attended an emergency department

Exploring the prevalence of ten polyomaviruses and two herpes viruses in breast cancer

Prevalence of human polyomaviruses in common and rare types of non-Merkel cell carcinoma skin cancer

Detection of novel polyomaviruses, TSPyV, HPyV6, HPyV7, HPyV9 and MWPyV in feces, urine, blood, respiratory swabs and cerebrospinal fluid

Human polyomaviruses in children undergoing transplantation

No evidence for association of HPyV6 or HPyV7 with different skin cancers

Submit your next manuscript to BioMed Central and take full advantage of: • Convenient online submission • Thorough peer review • No space constraints or color figure charges • Immediate publication on acceptance • Inclusion in PubMed, CAS, Scopus and Google Scholar • Research which is freely available for redistribution Submit your manuscript at www

We thank the Beijing Friendship Hospital, Capital Medical University, for providing the samples, and the laboratory staff who contributed to this study. No compensation was given to subjects for participation in this study. 

Additional file 1: Table S1 . Primers and probes used to detect respiratory viruses. Table S2 . Primers used to amplify HPyV6 with nested PCR (DOCX 28 kb)

The authors declare that they have no competing interests.