- Open Access
Prevalence of adenovirus in children with acute respiratory tract infection in Lanzhou, China
Virology Journalvolume 10, Article number: 271 (2013)
Human adenovirus (HAdV) is an important agent causing respiratory tract infection in children. Information on the epidemiological and clinical features of HAdV is limited in children with acute respiratory tract infections (ARTIs) in China, especially those of a novel genotype, Ad55.
In total, 1169 nasopharyngeal aspirates were collected from children younger than 14 years with ARTIs between November 2006 and November 2009. The polymerase chain reaction (PCR) was used to screen HAdVs. All PCR-positive products were sequenced.
74 of 1169 (6.33%) specimens were positive for HAdVs. Among positive cases, AdV3 (58/74) was detected most frequently, followed by AdV11 (10/74), AdV2 (2/74), AdV7 (2/69), AdV6 (1/74), and AdV1 (1/74). AdV55 was found in one case. The incidence of HAdV infection peaked in children aged 3–7 years. The most common clinical diagnosis was upper respiratory infection, and the most common syndrome was fever and cough.The comparison of HAdV and RSV group revealed that Children infected with group AdV were significant older than children infected with group RSV, had more fever but less frequently wheezing, and cough, crackles, and cyanosis, The duration of hospitalization between the AdV group and RSV group was not significant, but a greater frequency of LRTIs was observed in RSV group.
HAdV is an important viral agent in children with ARTIs in Lanzhou City, China. Multiple HAdV serotypes co-circulated with Ad3, which was predominant in this 3-year study. The novel AdV55 genotype was found in one case. No fixed seasonal rhythm could be identified.
Acute respiratory tract infections (ARTIs) are a major health problem worldwide, with high morbidity and mortality rates. Human adenovirus (HAdV) is not only an important cause of mild upper respiratory tract illness, but is also associated with more serious diseases, such as severe pneumonia. HAdV is responsible for 5–10% of lower respiratory tract infection (LRTI) in children [1, 2]. Additionally, the occurrence of fatal outcomes and chronic pulmonary sequelae associated with HAdV infection has been reported frequently [3–5].
HAdV infections can occur endemically throughout the year or in epidemics. In most studies, HAdV have been isolated endemically throughout the year [1, 6, 7], but numerous outbreaks of ARTI caused by HAdV have been reported during the last decade in many countries. Particular species and serotypes are more commonly associated with disease, syndromes, epidemiological settings, and demographic risk groups ; however, the prevalence of HAdV in Lanzhou remains unknow.
HAdV belongs to the Mastadenovirus genus and consists of seven species (A–G). To date, 67 types have been identified[8–13]. Whether these new types represent truly new strains or simply recombination events without clinical significance, needs to be further investigated. Respiratory diseases are generally associated with species B (serotypes 3 and 7), C (serotypes 1, 2, and 5), and E (serotype 4) viruses, whereas viruses in species F cause gastroenteritis, and viruses in subspecies D are often associated with keratoconjunctivitis. A study from Cuba recently reported that species D can also be isolated in respiratory specimens as a dominant genotype . In China, AdV3 and AdV7 have been reported as causative agents in epidemic outbreaks of respiratory disease. In 2009, a study from the Chinese Centers for Disease Control (CDC) reported an HAdV acute respiratory disease outbreak in Qishan, Shanxi Province, and AdV55 (QS-DLL), a new adenovirus not previously reported worldwide, was identified as the cause of these infections. The novel HAdV was shown to arise from a hexon recombination between AdV11 and AdV14 . This outbreak of AdV55 in China prompted us to review HAdV-associated ARTIs in Lanzhou City, China.
The objective of this study was to describe the epidemiological and clinical features of HAdV ARTIs over a period of 3 years.
Materials and methods
Study participants and samples
Between November 2006 and November 2009, 1169 nasopharyngeal aspirates (NPAs), consisting of 764 inpatient specimens and 405 outpatient specimens, were collected from 1169 children with ARTIs in the First Hospital of Lanzhou University, Gansu Province, China.
Informed consent was obtained from the parents of all children who provided specimens. The study protocol was approved by the hospital’s ethics committee (The Ethics Committee of the First Hospital of Lanzhou University).
All NPA specimens were collected and transported immediately to the laboratory at the National Institute for Viral Disease Control and Prevention, China CDC, and stored at -80°C until further analysis. Demographic data and clinical findings were recorded.
Viral DNA and RNA were extracted from 140 μL of each NPA specimen using the QIAamp viral DNA and RNA mini kits (Qiagen, Shanghai, China) according to the manufacturer’s protocol. cDNA was synthesized using random hexamer primers with Superscript II RH- reverse transcriptase (Invitrogen, Carlsbad, CA, USA).
In screening for HAdV, we used forward (TTCCCCATGGCICAYAACAC) and reverse (CCCTGGTAKCCRATRTTGTA) primers that target the partial region of the hexon gene to amplify a 482-bp fragment, as described elsewhere . The polymerase chain reaction (PCR) was performed under the following conditions: 94°C for 8 min, followed by 35 cycles at 94°C for 30 s, 50°C for 30 s, and 72°C for 45 s, and a final extension at 72°C for 10 min.
Detection of other respiratory viruses
All specimens were also screened for human respiratory syncytial virus (HRSV), human metapneumovirus (HMPV), influenza viruses A and B (IFVA, IFVB), parainfluenza virus types 1 to 3 (PIV1–3), human rhinoviruses (HRVs), and human coronaviruses (NL63 , HKU1) using a standard reverse-transcription PCR technique[16–19], and human Bocavirus (Bcov) using traditional PCR methods.
Nucleotide sequence analysis
The whole hexon gene was amplified using forward primer hexon-s and reverse primer hexon-as [11, 21] in our AdV11-positive specimens; the predicted product was 3449 bp. PCR was performed under the following conditions: 94°C for 1 min, followed by 35 cycles at 94°C for 30 s, 55°C for 30 s, and 68°C for 5 min, and a final extension at 68°C for 10 min.
All positive PCR products were purified using the QIAquick PCR purification kit (Qiagen), cloned into the pGEM-T Easy vector (Promega, Madison, WI, USA), and sequenced by HUADA Gene Company(Beijing, China). The nucleotide and deduced amino acid sequences of the hexon gene were compared with reference strains downloaded from GenBank. Phylogenetic analyses were conducted using MEGA software (ver. 5.05). Sequences were determined and analyzed using the DNAMAN software package .
The statistical significance of differences between the various groups was tested using the chi-squared test and Fisher’s exact test. A P value of <0.05 was considered to indicate statistical significance. All analyses were performed using the SPSS software (ver. 13.0.)
Between November 2006 and November 2009, 1169 patients were enrolled. The ages of the children with ARTIs ranged from 1 day to 168 months. The vast majority of patients (90.95%) were ≤ 5 years old. The ratio of boys to girls was 1.7:1 and that of outpatients to inpatients was 1:1.9.
Detection of HAdVs and other respiratory viruses
In total, 74 of 1169 (6.33%) specimens were positive for HAdV. The 74 HAdV-positive cases consisted of 16 obtained in 2006–2007, 26 in 2007–2008, and 32 obtained in 2008–2009. The rate of adenovirus coinfection with other respiratory viruses was 62% (48/74); the most common coinfected virus was HRSV (15), followed by HPIV (11), Bcov (7), IFVB (5), IFVA (3), HRV (4), HMPV (4), HKU1 (1), and 2NL63 (1). Adenoviruses represented 11.68% of all Viral agents identified.
Of the 74 HAdV-positive cases (Figure 1), AdV3 (58/74) was detected most frequently, followed by AdV11 (10/74), AdV2 (2/74), AdV7 (2/69), AdV6 (1/74), and AdV1 (1/74). AdV3 accounted for 78.38% of HAdV positive specimens, suggesting an epidemic of ARTIs due to AdV3 during the period from November 2006 to September 2009. AdV11 was isolated sporadically throughout the study period.
Identification of the novel HAdVs - AdV55
In total, 10 AdV11 viruses were detected in 1169 specimens, and the whole hexon gene was amplified successfully in one specimen (LZY153) collected in February 2009. The identity between LZY153 and QS-DLL (FJ643676) was 99.96%. Whole hexon gene sequence analysis showed intraspecies recombination between HAdV-11 and HAdV-14(data not show), similar to that seen for QS-DLL (FJ643676).
The monthly distribution of respiratory adenovirus infection is shown in Figure 1. No HAdV outbreak was identified in Lanzhou City, China, during 2006–2009. HAdV were detected in every month except May to October 2007, February and August 2008, and August, October, and November 2009 (Figure 2).
The mean age of the 74 AdV-positive patients was 21 (range, 0.5–144) months; 79.7% (59/74) of patients were younger than 5 years. Furthermore, 52.17% of AdV infections occurred in children <1 year of age; 11 patients were younger than 6 months, including two newborns. The peak incidence occurred at 3–4 years of age (Table 1). The ratio of males to females was 2:1.
Clinical characteristics of HAdVs infection in children
Information on clinical characteristics was available for 74 patients (Table 2). Upper respiratory infection (URTI) was the most common presentation (32/74, 43.24%), followed by bronchopneumonia (21/74, 28.38%), bronchiolitis (14/74, 18.91%), and bronchitis (7/74, 9.46%). Fever was recorded in 70.27% (52/74) of the patients; prolonged (duration > 10 days) and high (temperature > 39°C) fever was recorded in 32.43% (24/74) of the patients. Cough was observed in 66.22% (49/74) of patients, wheezing in 18(58.1%), and cyanosis in two patients. Gastrointestinal (GI) discomfort, such as diarrhea, was observed in 14 (18.92%).
We compared clinical characteristics between the AdV group and the RSV group (Table 3). Subjects coinfected with RSV and AdV were removed from further analysis, leaving 58 AdV-positive patients and 88 RSV-positive patients. Children infected with AdV were significantly older than children infected with RSV (chi-squared test, p < 0.001) and had fever more frequently (p = 0.030), but wheezing, cough, crackles, and cyanosis less frequently (Table 3). No significant difference in GI discomfort or the duration of hospitalization was observed between groups. But significant differences were observed in the frequency of URTI and LRTI between the two groups, and a greater frequency of LRTIs in RSV group(p < 0.001).
Chest radiographs were obtained from 18 AdV-positive patients; all showed abnormal findings. Peribronchial and perihilar infiltration were observed in 11 (11/18) of the children, coarse lung markings were observed in five (5/18), and one patient who was infected with HAdVs type 7 had pleural effusions and homogenous consolidation. Complete blood cell counts were performed for 40 of the 74 HAdVs-positive children. Most white blood cell (WBC) counts were within the normal range; for 25% (10/40) of the patients, WBC counts were >10 × 109 leukocytes/L. The mean WBC count was 8.41(range, 3.3–19.9 × 109 ) leukocytes/L. The overall hospitalization rate was 51.35% (38/74), with a median stay of 12.75 (range, 3–36) days.
In the present study, 1169 children with ARTIs were enrolled; more than half of the children (73.55%) were infected with at least one viral agent, and HAdV were detected in 74 children (6.33%), accounting for 11.68% of total virus isolates. Similar incidences of HAdV have been reported in Cuba, Taiwan, Brazil, and South Korea [1, 14, 22, 23]. However, higher frequencies of HAdV were observed in Argentina, where rates reached 14.3% in 168 hospitalized children with LRTI during the period from 1991 through 1995. HAdV infections can occur endemically throughout the year or in epidemics. Epidemics of adenoviral respiratory disease are common during the winter and spring. In the present study, although the HAdV detection rate was relatively high in winter and spring, in contrast to HRSV infections [25, 26], no fixed seasonal rhythm of HAdV was identified (Figure 2).
HAdV are responsible for 4–10% of LRTIs in children . Serotypes 1–3, 5, and 7 are frequently recovered from children with ARTIs, but outbreaks of severe infections in healthy children are most frequently reported with AdV7, followed by AdV3 and AdV21 [3, 27, 28]. In hospitalized Korean children with LRTIs due to HAdV, the serotypes recovered most frequently between 1990 and 1998 were AdV7 (41%), AdV3 (15%), and AdV2 (15%) . A two-decade study from Taiwan showed that AdV3 and AdV7 were the major detected types; AdV3 was consistently isolated throughout the survey period . HAdV types 3 and 7 are also two of the most important etiological agents of pneumonia in China; a previous study from China showed that AdV3 was the most prevalent genome type in children with pneumonia over the 27-year period from 1962 to 1988. AdV3 and AdV7 accounted for 69–100% of HAdV strains isolated from patients with pneumonia in Beijing, China . In the present study, six serotypes (types 1, 2, 3, 6, 7, 11) were observed between November 2006 and November 2009 in Lanzhou City, China. AdV3 was predominant during this period, accounting for 78.38% of HAdV. Genome recombination plays an important role in the molecular evolution of HAdV, leading to newly emerging strains that have tropism changes or become virulent [8–13, 30]. Interestingly, recombination between AdV11 and AdV14 was found in our AdV11 isolates. AdV11, which has rarely been reported worldwide in ARTIs, accounted for 13.51% of HAdV isolates in our study. AdV14, which belongs to the species B2 , has been reported to cause outbreaks of severe respiratory disease in many places recently . AdV11, which also belongs to subgenus B2, was isolated more frequently in China between 1965 and 1985 . Since then, AdV11 infections had not been reported for more than 20 years until an acute respiratory outbreak caused by AdV55, showing a hexon recombination between AdV11 and AdV14, occurred in Shanxi Province, China, in the spring of 2008. Our study confirmed that AdV55 was present and circulating in Lanzhou, China.
A previous study indicated that the incidence of HAdV infection was inversely related to age . In our study, HAdV was detected most commonly in children aged 3–7 years (10.00%), and the majority of cases (79.73%) were < 5 years old. These results are consistent with those of studies from Greater Manchester, UK, and Korea, which also revealed that most HAdV infections occurred among children less than 5 years of age . However, in the present study, HAdVs infection was also found in 11 infants younger than 6 months, including two newborns. This result is in contrast to that of a previous study from Taiwan, which found that all HAdVs-positive patients except one infant were older than 6 months ; this phenomenon was explained by the reasoning that infants < 6 months of age are normally protected by maternal antibodies, although the exact mechanism requires further study. Our study also revealed that patients with AdV3 were significantly older than patients infected with other virus types (median age, 30 vs. 8 months), which is also consistent with a previous study from Taiwan. Males were more frequently infected than females with all HAdV serotypes in our study, which is also consistent with the results of a previous study . All of these issues require further study.
The most common clinical diagnosis was URTI (43.24%), and 56.76% of our HAdV-positive patients had lower respiratory tract conditions, including pneumonia, bronchitis, and bronchiolitis. HAdVs were less frequent causes of ARTIs in children than were HRSV, HPIV, and HRV in our study. Furthermore, significant differences were observed in the frequency of URTI and LRTI between the RSV group and HAdV group, and a greater frequency of URTIs in HAdV group. Most LRTIs caused by HAdV are mild and indistinguishable from LRTIs caused by other respiratory viruses; however, HAdV can result in fatal outcomes or residual sequelae. In a previous study from Buenos Aires, Argentina, the overall fatality rate during hospitalization of patients with AdV7 reached 28.6%. Another study from China conducted over a 33-year period (1958–1990) reported a 15.5% fatality rate of HAdV-associated infantile pneumonia. The symptoms manifested during infection with AdV3 and AdV7 have become less severe since the mid-1980s, and no fatal case has been reported since 1986 . In our 3-year study, the RTIs associated with AdV3 resulted in no mortality. This finding could be because of a stabilization of the circulating HAdV type 3 in China, with subsequent increases in herd immunity and reductions in HAdV infections. AdV7 is most frequently reported elsewhere and may be related to severe infections in healthy children . In the present study, only two patients were AdV7-positive; pleural effusions and homogenous consolidation were identified on chest radiographs in one of these patients, with negative results by bacterial culturing. AdV11, which has rarely been reported in respiratory specimens, was second in frequency only to AdV3 in this study and was detected in 14.86% of AdV-positive samples. Sample LZY153, which was confirmed to contain the new AdV55 genotype with hexon recombination between AdV11 and AdV14, was collected in February 2009 from a 12-month-old boy who was admitted to the hospital because of bronchitis and eventually recovered. Given that our sample contained only one confirmed AdV55 case, the clinical significance of this genotype is difficult to analyze. Furthermore, interpretation of the general impact of HAdV infection is also restricted by the difficulty of patient follow-up for the detection of lung sequelae. Long-term surveillance of infection with different HAdV genotypes, especially newly identified genotypes, is needed to explore their clinical effects.
Of the 10 AdV-11 positive samples, only one (LZY153) was successfully amplified for the whole 3449 bp long hexon gene and was identified as AdV 55, which was assumed to be caused by the low amplification efficiency by our PCR for the long fragment and/or the degradation of the DNA in the specimens (the NPAs collected at 2006–2009 ). AdV 55 has been confirmed as a novel recombinant between AdV-11 and AdV-14 at the hexon gene . And it is difficult to identify the other nine AdV11-positive samples as AdV55 by the short fragment of the hexon gene. However, this novel recombinant AdV55 needs further study to explore its clinical significance after more cases with AdV55 were reported and confirmed.
In conclusion, the present study confirmed that multiple HAdV types co-circulated in Lanzhou City, China. AdV3 predominated during November 2006 to November 2009. We also found that AdV55 existed and circulated in China. No fixed seasonal rhythm of HAdV was identified. The incidence of HAdV infection peaked in children aged 3–7 years, but a relatively high infection rate (3.29%) was observed in infants younger than 6 months, which should be given more attention. Further molecular surveillance of HAdV should be undertaken continuously to explore their epidemiological and clinical characteristics.
The present study showed that HAdV is an important viral agent in children with ARTIs in Lanzhou City, China. Multiple HAdV types co-circulated with Ad3, which was predominant in this 3-year study. The novel AdV55 genotype was found in one case. No fixed seasonal rhythm could be identified.
Chen HL, et al.: Respiratory adenoviral infections in children: a study of hospitalized cases in southern Taiwan in 2001–2002. J Trop Pediatr 2004,50(5):279-284. 10.1093/tropej/50.5.279
Yun BY, et al.: Viral etiology and epidemiology of acute lower respiratory tract infections in Korean children. Pediatr Infect Dis J 1995,14(12):1054-1059. 10.1097/00006454-199512000-00005
Hong JY, et al.: Lower respiratory tract infections due to adenovirus in hospitalized Korean children: epidemiology, clinical features, and prognosis. Clin Infect Dis 2001,32(10):1423-1429. 10.1086/320146
Simila S, et al.: Chronic lung damage caused by adenovirus type 7: a ten-year follow-up study. Chest 1981,80(2):127-131. 10.1378/chest.80.2.127
Simila S, Ylikorkala O, Wasz-Hockert O: Type 7 adenovirus pneumonia. J Pediatr 1971,79(4):605-611. 10.1016/S0022-3476(71)80307-4
Kajon AE, et al.: Molecular epidemiology of adenovirus acute lower respiratory infections of children in the south cone of South America (1991–1994). J Med Virol 1996,48(2):151-156. 10.1002/(SICI)1096-9071(199602)48:2<151::AID-JMV6>3.0.CO;2-8
Choi EH, et al.: The association of newly identified respiratory viruses with lower respiratory tract infections in Korean children, 2000–2005. Clin Infect Dis 2006,43(5):585-592. 10.1086/506350
Liu EB, et al.: Genetic analysis of a novel human adenovirus with a serologically unique hexon and a recombinant fiber gene. PLoS One 2011,6(9):e24491. 10.1371/journal.pone.0024491
Walsh MP, et al.: Computational analysis of two species C human adenoviruses provides evidence of a novel virus. J Clin Microbiol 2011,49(10):3482-3490. 10.1128/JCM.00156-11
Robinson CM, et al.: Computational analysis and identification of an emergent human adenovirus pathogen implicated in a respiratory fatality. Virology 2011,409(2):141-147. 10.1016/j.virol.2010.10.020
Yang Z, et al.: Genomic analyses of recombinant adenovirus type 11a in China. J Clin Microbiol 2009,47(10):3082-3090. 10.1128/JCM.00282-09
Matsushima Y, et al.: Genome sequence of a novel virus of the species human adenovirus d associated with acute gastroenteritis. Genome Announc 2013,1(1):e00068-12.
Matsushima Y, et al.: Novel human adenovirus strain. Bangladesh. Emerg Infect Dis 2012,18(5):846-848. 10.3201/eid1805.111584
Belsy A, et al.: Molecular characterization of adenoviral infections in Cuba: report of an unusual association of species D adenoviruses with different clinical syndromes. Arch Virol 2009,154(4):619-627. 10.1007/s00705-009-0338-4
Yan H, et al.: Development of RT-multiplex PCR assay for detection of adenovirus and group A and C rotaviruses in diarrheal fecal specimens from children in China. Kansenshogaku Zasshi 2004,78(8):699-709.
Vabret A, et al.: Detection of the new human coronavirus HKU1: a report of 6 cases. Clin Infect Dis 2006,42(5):634-639. 10.1086/500136
Bastien N, et al.: Human coronavirus NL63 infection in Canada. J Infect Dis 2005,191(4):503-506. 10.1086/426869
Bellau-Pujol S, et al.: Development of three multiplex RT-PCR assays for the detection of 12 respiratory RNA viruses. J Virol Methods 2005,126(1–2):53-63.
Vabret A, et al.: Direct diagnosis of human respiratory coronaviruses 229E and OC43 by the polymerase chain reaction. J Virol Methods 2001,97(1–2):59-66.
Hierholzer JC, et al.: Detection of adenovirus in clinical specimens by polymerase chain reaction and liquid-phase hybridization quantitated by time-resolved fluorometry. J Clin Microbiol 1993,31(7):1886-1891.
Zhu Z, et al.: Outbreak of acute respiratory disease in China caused by B2 species of adenovirus type 11. J Clin Microbiol 2009,47(3):697-703. 10.1128/JCM.01769-08
Moura PO, et al.: Molecular epidemiology of human adenovirus isolated from children hospitalized with acute respiratory infection in Sao Paulo. Brazil. J Med Virol 2007,79(2):174-181. 10.1002/jmv.20778
Kim YJ, et al.: Genome type analysis of adenovirus types 3 and 7 isolated during successive outbreaks of lower respiratory tract infections in children. J Clin Microbiol 2003,41(10):4594-4599. 10.1128/JCM.41.10.4594-4599.2003
Carballal G, et al.: Adenovirus type 7 associated with severe and fatal acute lower respiratory infections in Argentine children. BMC Pediatr 2002, 2: 6. 10.1186/1471-2431-2-6
Jin Y, et al.: Newly identified respiratory viruses associated with acute lower respiratory tract infections in children in Lanzou, China, from 2006 to 2009. Clin Microbiol Infect 2012,18(1):74-80. 10.1111/j.1469-0691.2011.03541.x
Zhang RF, et al.: Human respiratory syncytial virus in children with acute respiratory tract infections in China. J Clin Microbiol 2010,48(11):4193-4199. 10.1128/JCM.00179-10
Kajon AE, et al.: A new genome type of adenovirus 3 associated with severe lower acute respiratory infection in children. J Med Virol 1990,30(1):73-76. 10.1002/jmv.1890300116
Li QG, et al.: Molecular epidemiology of adenovirus types 3 and 7 isolated from children with pneumonia in Beijing. J Med Virol 1996,49(3):170-177. 10.1002/(SICI)1096-9071(199607)49:3<170::AID-JMV3>3.0.CO;2-1
Lin KH, et al.: A two decade survey of respiratory adenovirus in Taiwan: the reemergence of adenovirus types 7 and 4. J Med Virol 2004,73(2):274-279. 10.1002/jmv.20087
Walsh MP, et al.: Computational analysis identifies human adenovirus type 55 as a re-emergent acute respiratory disease pathogen. J Clin Microbiol 2010,48(3):991-993. 10.1128/JCM.01694-09
Lewis PF, et al.: A community-based outbreak of severe respiratory illness caused by human adenovirus serotype 14. J Infect Dis 2009,199(10):1427-1434. 10.1086/598521
Li QG, Hambraeus J, Wadell G: Genetic relationship between thirteen genome types of adenovirus 11, 34, and 35 with different tropisms. Intervirology 1991,32(6):338-350.
Cooper RJ, et al.: The epidemiology of adenovirus infections in Greater Manchester, UK 1982–96. Epidemiol Infect 2000,125(2):333-345. 10.1017/S0950268899004550
Cheng CC, et al.: Molecular and clinical characteristics of adenoviral infections in Taiwanese children in 2004–2005. Eur J Pediatr 2008,167(6):633-640. 10.1007/s00431-007-0562-4
This work was supported in part by the China Mega-Project for Infectious Disease (2011ZX10004-001), and National Natural Science Foundation of China(2012, 81260001).
The authors declare that they have no competing interests.
Rong-fang Zhang carried out RT-PCR, virus detection and sequencing studies. Zhi-ping Xie, Kun-long Yan, Han-chun Gao, Jing-rong Song, Xin-hui Yuan was involved in clinical diagnosis of patients and collecting samples. Yu Jin and Rong-fang Zhang drafted the manuscript.Yu Jin, Yun-de Hou, Zhao-jun Duan conceived the study, made the arrangements for obtaining samples from the hospitals , and finalized the manuscript. All authors read and approved the final manuscript.
Yu Jin, Rong-fang Zhang contributed equally to this work.