- Research article
- Open Access
Rapid detection of human mastadenovirus species B by recombinase polymerase amplification assay
BMC Microbiologyvolume 19, Article number: 8 (2019)
As an important component of the causative agent of respiratory tract infections, enteric and eye infections, Human mastadenoviruses (HAdVs) species B spread easily in the crowd. In this study, we developed a recombinase polymerase amplification (RPA) assay for rapidly detecting HAdVs species B which was comprised of two different formats (real-time and lateral-flow device).
This assay was confirmed to be able to detect 5 different HAdVs species B subtypes (HAdV-B3, HAdV-B7, HAdV-B11, HAdV-B14 and HAdV-B55) without cross-reactions with other subtypes and other respiratory tract pathogens. This RPA assay has not only highly sensitivity with low detection limit of 50 copies per reaction but also short reaction time (< 15 min per detection). Furthermore, the real-time RPA assay has excellent correlation with real-time PCR assay for detection of HAdVs species B presented in clinical samples.
Thus, the RPA assay developed in this study provides an effective and portable approach for the rapid detection of HAdVs species B.
Human mastadenoviruses (HAdVs) are an important cause of human disease, most of which were associated with upper and lower respiratory tract infections, enteric or eye infections. In addition, hepatitis, meningoencephalitis, cystitis and myocarditis, even noninflammatory conditions, such as obesity were also found to be clinical manifestations of HAdVs infections [1,2,3]. Recently, the putative oncogenicity of certain types of the HAdVs in human was proposed based on the results from mammalian animal models and human malignant diseases [4,5,6]. So far, there have already been identified out more than 60 serotypes of AdVs in the world , about one-third of which were confirmed to be associated with human infections . HAdVs infection happened not only in children but also in adults. But more than 80% of infections of HAdVs occured in children, especially those younger than 4 years old . HAdVs are divided into seven species (HAdV-A to G) on the basis of DNA genome homology . The serotypes most frequently associated with respiratory infection are belonged to species B, species C and species E. Some of HAdV-B and HAdV-C, such as HAdV-B3, HAdV-B7, HAdV-B11, HAdV-B14, HAdV-B55 and HAdV-C2 are common causative agents of respiratory tract infections in children [1, 10]. HAdV-F, such as HAdV-F41, is frequently implicated as a cause of infectious diarrhea worldwide .
HAdVs infections are easily spread among the crowd. Timely diagnosis of HAdVs infections will undoubtedly help to prevent theirs transmission. Although viral isolation was the gold standard to identify adenoviral infection, it was time-consuming and needed skilled performers. On the other hand, immunological assay, such as direct immunofluorescence, was rapid to detect viruses. However, this assay usually did not show enough sensitivity and needed special antiserum. With the progress in molecular detection techniques, PCR or real-time PCR assays have rendered rapid and sensitive tools for detection, typing, and monitoring of adenoviral infections [12,13,14,15]. However, field application was impeded in resource-limited setting due to potential contamination and needing expensive equipment. Novel isothermal amplification techniques provided substituted tools for HAdVs detection obviating the need for a thermal cycler. Several isothermal amplification assays for detecting HAdVs have been developed [16,17,18]. All of these assays adopted loop-mediated isothermal amplification (LAMP) methods. In this study, we present a recently developed isothermal amplification method, named recombinase polymerase amplification (RPA), for detecting HAdVs species B. RPA was the proprietary technology of TwistDx Inc.UK. The RPA reaction exploits recombinase, single-stranded DNA binding (SSB) protein and polymerase to initiate the amplification reaction at low temperature (around 37 °C) within minutes. Three TwistAmp amplification formulations were supplied by the provider. TwistAmp basic kit contains the above basic components for the amplification. Amplicons will typically be assessed by an endpoint method, such as gel electrophoresis. TwistAmp exo kit which contains a powerful nuclease (Exonuclease III) can process TwistAmp exo probes during the amplification reaction itself and generate a real-time readout. TwistAmp nfo kit which includes a nfo nuclease (Endonuclease IV) can generate new polymerase extension substrates with suitable nof probes. The antigenic label on the 5′ end of the probe (typically 6-Carboxyfluorescein, FAM or fluorescein isothiocyanate, FITC) becomes conjoined with an antigenic label on the 5′ end of the opposing amplification primer (typically biotin or DIG) and this association can be detected in a ‘sandwich’ assay such as lateral-flow (LF) strips. To resolve the need for resource-limited settings, portable real-time fluorescence scanner device and lateral flow device (LFD) were applied to detect the RPA products in present study.
Viral isolates and clinical specimens
The six virus strains of HAdVs (HAdV-C2, HAdV-B3, HAdV-B7, HAdV-B11, HAdV-B14, and HAdV-B55) isolated from patient’s samples in our institute using Hep-2 cells. The cells were grown in RPMI-1640 (Gabico, USA) supplemented with 10% fetal bovine serum (FBS) (Gabico, USA) and 1% penicillin/streptomycin (Gabico, USA). The cells were cultured at 37 °C in a humidified atmosphere of air containing 5% CO2. Other isolates (positive for the following pathogens including legionella pneumophila, streptococcus pneumoniae, neisseria meningitides and haemophilus influenza) and clinical specimens (positive for the following pathogens including HAdV-F41, mycoplasma pneumoniae, chlamydia pneumoniae, human bocavirus, and herpes simplex virus) identified by pathogen specific real-time PCR assays in our institute were used to evaluate the specificity. The clinical specimens, a total of 186 pharyngeal swabs were obtained from the children who suffered from acute respiratory infection within 3 days after the onset of clinical symptoms in 2013–2016. All these isolates and samples were stored at − 70 °C until used for further detection and analysis.
Primers and probes design and screen
Oligonucleotide primers used in the RPA assays were manually designed based on the conserved region of the hexon gene of HAdVs species B according to the instruction manual provided by TwistDx (Cambridge, United Kingdom), the manufacturer of the RPA reaction. Five sets of primers for candidates were designed and pairs were screened by observing their performance on a 2% agarose gel. After the outer primer pairs were determined, internal probes were screened through a similar process. The set of primers and probe with the brightest band on the agarose gel were chosen for further test.
To generate a detectable product on lateral flow strips, the 5′ end of the forward primer was labeled with a 5′-biotin. The nfo probe was labeled on its 5′ end with a FITC group, a C3 spacer (SpC3) on the 3′ end and a tetrahydrofuran residue (THF) which replaces an internal base. For real-time RPA detection, same primers were used except for without of biotinylation at the 5′ end of the forward primer. Exo probe was conjugated with the FAM and the black hole quencher (BHQ) to the T-bases at internal positions with a THF located in the central part of the two fluorescent groups and a SpC3 labeled on the 3′ end. All oligonucleotides were synthesized by Sangon Biotech (Shanghai, China).
Preparation of plasmid and DNA extraction
PCR product of the hexon gene sequence was cloned into a pMD19-T vector. The plasmid was extracted from culture using the TaKaRa MiniBEST Plasmid Purification Kit (TaKaRa, China) following the manufacturer’s instruction. The concentration was measured by spectrophotometry at 260 nm. The expected copy number of the target gene was calculated according to a previously described formula . The purified plasmid was used as a template to optimize the RPA assays. Genomic DNA from multiple strains and clinical samples was extracted from samples using QIAamp® MinElute® Virus Spin Kit (Qiagen) according to the manufacturer’s instruction and eluted in a final volume of 50 μl. Extracted DNA was stored at − 70 °C until further use.
Real-time PCR assay for HAdVs was adapted from previous published procedure with slight modifications . In brief, real-time PCR reactions were carried out by using Premix Ex Taq Kit (TaKaRa, China). Twenty microliter reaction volume contained 10 μl of 2 × PCR Buffer, 200 nM of primer mix, 100 nM of probe mix, and 5 μl of extracted DNA. Real-time PCR cycling was performed on LightCycler Real-Time PCR system (Roche) as follows: a 30 s heat activation of the Taq polymerase at 95 °C followed by 40 amplification cycles of 5 s at 95 °C and 20 s at 56 °C each (annealing-extension step).
LFD-RPA assay was performed in a 50 μl volume with the TwistAmp nfo kit (TwistDx, Cambridge, UK) containing 420 nM each primer, 120 nM target-specific RPA nfo probe, and 1x rehydration buffer. After mixing these components with dry reagent pellet, 5 μl of DNA was added to the reaction mixture, and 2.5 μl of magnesium acetate (Mg(OAc)2, 280 mM) was pipetted into the tube lids. The lids were closed and the magnesium acetate was spun down into the reaction mixture to trigger the reaction. The reaction mix was placed into the heating block with brief mixing and centrifugation after 3–4 min of incubation. Three temperatures (37, 39, and 42 °C) and six reaction times (5, 10, 15, 20 25 and 30 min) were tested during RPA assay. RPA products were detected using LFD (Ustar Biotech Co. Ltd. Hangzhou, China) according to previous procedure . Briefly, the reaction tube was placed into cartridge and was immobilized by closing the cartridge. Then the cartridge was inserted into detection chamber whose handle was then closed to cut open running buffer reservior and reaction tube. Finally the detection result was read by naked eyes after 5 to 10 min.
Real-time RPA assay
Real-time RPA reaction was set up similar to that of LFD-RPA assay except for using a TwistAmp Exo kit (TwistDX, Cambridge, UK), primers without biotinylation, and nfo probe. The reaction tubes were mixed, centrifuged and then placed into RAA isothermal amplification and detection device (Qitian Co, Wuxi, China). The reaction was performed at 39 °C for 30 min, with brief mixing and centrifugation of reaction tubes after 3–4 min of incubation.
Sensitivity and specificity of the RPA reactions
The sensitivity of the assay was evaluated using 10-fold serial dilutions of plasmid standards (ranging from 1 × 100–1 × 104copies/μl). The specificity of the RPA was determined by analyzing different strains or clinical specimens positive for the pathogens including differential types of HAdVs (HAdV-C2, HAdV-F41), mycoplasma pneumoniae, chlamydia pneumoniae, legionella pneumophila, streptococcus pneumoniae, neisseria meningitides, haemophilus influenza and herpes simplex virus.
Clinical sample analysis and sequencing
To evaluate the performance characteristics of the real-time RPA assay in clinical sample detection, all of the 186 pharyngeal swab specimens were subjected to real-time RPA assay with the parallel analysis by real-time PCR. The positive specimens were further being molecular typing using the protocol described previously .
Optimization of the RPA conditions
According to the standard mentioned methods, one set of primers, one nfo and one exo probe were selected (Table 1) for further optimization. RPA was conducted in 15 min at 37, 39 and 42 °C using 1000 copies of plasmid DNA respectively. There were no significant differences in the LFD-RPA and real-time RPA results for amplification at different temperature. Therefore, the temperature at 39 °C was chosen for the subsequent RPA assays. Six different times comprising of 5, 10, 15, 20, 25 and 30 min were used in the LFD-RPA reaction. A signal on the lateral flow strip can be detected after 5 min of amplification reaction at 39 °C. The intensity of RPA amplicons at 10, 15, 20 and 25 min showed no significant differences. However, nonspecific band occurred in the negative controls when reaction time reaching 30 min. To avoid false positive result, 20 min was selected as an optimal reaction time for RPA assay.
Analytical sensitivity of the RPA assay
The detection threshold of the RPA was determined by using a diluted series of plasmid DNA ranging from 5 to 50,000 copies per reaction. The amplified products were analyzed by LFD and real-time fluorescence monitor. Repeatability of the method was assessed by testing each dilution 5 times. The results were compared with that of real-time PCR. Real-time PCR and LFD-RPA assays detected as low as 5 copies per reaction while real-time RPA assay could detect as low as 50 copies per reaction (Table 2 and Fig. 1).
Analytical specificity of RPA assay
The analytical specificity of the RPA assay was evaluated by testing DNA extracts from various pathogens that could cause similar respiratory symptoms. The results showed that the positive reactions were only observed in the extracts of HAdVs species B while testing results for all other pathogens including mycoplasma pneumoniae, chlamydia pneumoniae, legionella pneumophila, streptococcus pneumoniae, neisseria meningitides, haemophilus influenza, human bocavirus and herpes simplex virus were negative. These results indicate that the analytical specificity of the RPA assay was 100%. To evaluate the applicability of RPA assay to detect HAdVs, seven isolates of HAdVs including five species B (HAdV-B3, HAdV-B7, HAdV-B11, HAdV-B14 and HAdV-B55), one species C (HAdV-C2),one species F (HAdV-F41) were tested. The results revealed HAdVs species B could be detected by both LFD RPA and real-time RPA assays while other species could not be detected (Fig. 2).
Comparison of real-time RPA assay and real-time PCR to detect HAdVs
Using a total of 186 pharyngeal swabs collected from the children who suffered from acute respiratory infection, the clinical performance characteristics of the real-time RPA assay were assessed with using real-time PCR as reference. As showed in Table 3, 35 pharyngeal swabs were detected positive by real-time PCR assay with Ct values ranging from 15.6 to 35.3. Of the 35 samples detected positive in real-time PCR, 33 were detected positive by real-time RPA assay, resulting in a sensitivity of 94.3%. Of the 151 samples tested negative in real-time PCR, all were also tested negative by real-time RPA assay, providing a specificity of 100%. Two samples tested positive by real-time PCR but negative by real-time RPA were identified by sequencing as HAdV-C2. Other 33 samples tested positive by both methods were identified by sequencing contained 22 HAdV-B3, 7 HAdV-B7, 2 HAdV-B11, and 2 HAdV-B55.
HAdV species B is one of the most common pathogens associated with acute respiratory infections, particularly in infants and children. It was difficult for clinicians to distinguish HAdVs infections from others, because clinical symptoms caused by HAdVs infections were similar to those caused by other respiratory pathogens. Thus, rapid, sensitive and specific diagnostic assay for detecting HAdVs species B are urgently needed and would play an important role on initial clinical treatment, avoidance of antibiotic misuse, as well as prevention of virus transmission, especially in developing countries. RPA was a new isothermal amplification method with short turnaround time, low temperature, high sensitivity, and specificity which has been used to detect various viruses [22,23,24]. In this study, we developed a novel molecular method to detect HAdVs species B by combining RPA with real-time fluorescence detection and LFD technique. Since RPA assay can be conducted at a comparatively low constant temperature (from 37 to 42 °C), it could be carried out with a portable real-time fluorescence scanner device without thermal cycler. On the other hand, the read-out can also be observed by naked eyes with LFD which avoids cross-contamination between specimens [20, 25]. These characteristics of RPA assay for HAdVs species B detection developed in this study is suitable for use as a point-of-care molecular diagnostic assay or in low-equipment setting.
Three isothermal amplification methods for the detection of HAdVs were developed based on LAMP technique which can obtain the results within 45 to 60 min [16,17,18]. Compared with LAMP, RPA assay performed in shorter reaction time. The result can be read out within 30 min (20 min RPA plus 5 min on LFD). LAMP requires four to six primers to initialize reaction, while RPA-LFD assay only requires three primers. At the time this article was written, no specialized software is available for design RPA primers. It could be further simplified RPA primer design if specialized software would be developed in the future. Previous studies have shown the sensitivity of RPA is parallel to that of LAMP assay while its specificity is higher than that of LAMP . The developed LFD-RPA assay has shown high sensitivity for detection of HAdVs (5 copies/reaction) which is equal to real-time PCR assay, while higher than LAMP assay (50 to 10,000 copies/reaction) [17, 18]. Although the detection limit of real-time RPA is slightly higher than that of LFD-RPA, this level of sensitivity is sufficient to detect HAdVs species B in clinical samples of pediatric patients which the mean virus loads were 1.7 × 105 copies/ml . No positive results were obtained from a large panel of DNA samples from various viruses and bacteria causing respiratory infection. However, we observed nonspecific amplification bands when reaction time exceeded 30 min in LFD-RPA assay. This may due to the clotting of proteins and primers since all tested samples appeared nonspecific bands when reaction time exceeded 30 min. Similar results were also observed in the detection for yellow fever virus . The false-positive results of LFD-RPA limited its field application. So, we only evaluated characteristics of the real-time RPA assay clinically. The results showed real-time RPA assay was able to identify out 94.3% HAdVs positive samples, thereby demonstrating a high degree of specificity. It should be pointed out that this assay could not distinguish different serotypes although 7 serotype HAdVs were tested in this study and could not detect all serotypes. Since there are more than 60 serotypes ADV, it was difficult to design the primers for all serotypes and degenerate primers and probes should be considered to increase its range of application. So far, only one on-chip RPA assay to detect HAdV-F41 was established with the limit of detection of 35 genomic units/μl . Although we successfully identified out the most of clinical samples which were HAdVs positive as confirmed by real-time PCR assay in this study, more samples of other serotypes need to be tested in the future.
In addition, sample preparation techniques capable of extracting DNA in a manner compatible with the RPA assay need to be developed before this assay being used under field conditions. Alkaline lysis and heat treatment seem to be simple and appropriate methods since they have been successfully used in LAMP assay for detection HAdVs . RPA can tolerate inhibitory substances and temperature variations. Whether those methods are fit for RPA assay needs to be tested in the future .
We have developed a very rapid and sensitive isothermal RPA assay for the detection of HAdVs species B infections. This assay has great potential to be a point-of-care diagnostic tool if suitable sample preparation method has been developed.
Loop-mediated isothermal amplification
Lateral flow device
Polymerase chain reaction
Recombinase polymerase amplification
Lion T. Adenovirus infections in immunocompetent and immunocompromised patients. Clin Microbiol Rev. 2014;27(3):441–62.
Robinson CM, Singh G, Henquell C, Walsh MP, Peigue-Lafeuille H, Seto D, Jones MS, Dyer DW, Chodosh J. Computational analysis and identification of an emergent human adenovirus pathogen implicated in a respiratory fatality. Virology. 2011;409(2):141–7.
Walsh MP, Seto J, Jones MS, Chodosh J, Xu W, Seto D. Computational analysis identifies human adenovirus type 55 as a re-emergent acute respiratory disease pathogen. J Clin Microbiol. 2010;48(3):991–3.
Hohlweg U, Dorn A, Hosel M, Webb D, Buettner R, Doerfler W. Tumorigenesis by adenovirus type 12 in newborn Syrian hamsters. Curr Top Microbiol Immunol. 2004;273:215–44.
Kosulin K, Haberler C, Hainfellner JA, Amann G, Lang S, Lion T. Investigation of adenovirus occurrence in pediatric tumor entities. J Virol. 2007;81(14):7629–35.
Kosulin K, Hoffmann F, Clauditz TS, Wilczak W, Dobner T. Presence of adenovirus species C in infiltrating lymphocytes of human sarcoma. PLoS One. 2013;8(5):e63646.
Robinson CM, Singh G, Lee JY, Dehghan S, Rajaiya J, Liu EB, Yousuf MA, Betensky RA, Jones MS, Dyer DW, et al. Molecular evolution of human adenoviruses. Sci Rep. 2013;3:1812.
Cesario TC. Viruses associated with pneumonia in adults. Clin Infect Dis. 2012;55(1):107–13.
Barrero PR, Valinotto LE, Tittarelli E, Mistchenko AS. Molecular typing of adenoviruses in pediatric respiratory infections in Buenos Aires, Argentina (1999-2010). J Clin Virol. 2012;53(2):145–50.
Matsushima Y, Shimizu H, Kano A, Nakajima E, Ishimaru Y, Dey SK, Watanabe Y, Adachi F, Suzuki K, Mitani K, et al. Novel human adenovirus strain, Bangladesh. Emerg Infect Dis. 2012;18(5):846–8.
Roy S, Vandenberghe LH, Kryazhimskiy S, Grant R, Calcedo R, Yuan X, Keough M, Sandhu A, Wang Q, Medina-Jaszek CA, et al. Isolation and characterization of adenoviruses persistently shed from the gastrointestinal tract of non-human primates. PLoS Pathog. 2009;5(7):e1000503.
van Heerden J, Ehlers MM, Heim A, Grabow WO. Prevalence, quantification and typing of adenoviruses detected in river and treated drinking water in South Africa. J Appl Microbiol. 2005;99(2):234–42.
Heim A, Ebnet C, Harste G, Pring-Akerblom P. Rapid and quantitative detection of human adenovirus DNA by real-time PCR. J Med Virol. 2003;70(2):228–39.
He JW, Jiang S. Quantification of enterococci and human adenoviruses in environmental samples by real-time PCR. Appl Environ Microbiol. 2005;71(5):2250–5.
Huang ML, Nguy L, Ferrenberg J, Boeckh M, Cent A, Corey L. Development of multiplexed real-time quantitative polymerase chain reaction assay for detecting human adenoviruses. Diagn Microbiol Infect Dis. 2008;62(3):263–71.
Li F, Zhao LQ, Deng J, Zhu RN, Sun Y, Liu LY, Li YY, Qian Y. Detecting human adenoviruses in respiratory samples collected from children with acute respiratory infections by loop-mediated isothermal amplification. Zhonghua er ke za zhi = Chinese journal of pediatrics. 2013;51(1):52–7.
Ziros PG, Kokkinos PA, Allard A, Vantarakis A. Development and evaluation of a loop-mediated isothermal amplification assay for the detection of adenovirus 40 and 41. Food Environ Virol. 2015;7(3):276–85.
Wakabayashi T, Yamashita R, Kakita T, Kakita M, Oshika T. Rapid and sensitive diagnosis of adenoviral keratoconjunctivitis by loop-mediated isothermal amplification (LAMP) method. Curr Eye Res. 2004;29(4–5):219–24.
Min BS, Noh YJ, Shin JH, Baek SY, Min KI, Ryu SR, Kim BG, Park MK, Choi SE, Yang EH, et al. Assessment of the quantitative real-time polymerase chain reaction using a cDNA standard for human group a rotavirus. J Virol Methods. 2006;137(2):280–6.
Cui L, Ge Y, Qi X, Xu G, Li H, Zhao K, Wu B, Shi Z, Guo X, Hu L, et al. Detection of severe fever with thrombocytopenia syndrome virus by reverse transcription-cross-priming amplification coupled with vertical flow visualization. J Clin Microbiol. 2012;50(12):3881–5.
Lu X, Erdman DD. Molecular typing of human adenoviruses by PCR and sequencing of a partial region of the hexon gene. Arch Virol. 2006;151(8):1587–602.
Ma S, Li X, Peng B, Wu W, Wang X, Liu H, Yuan L, Fang S, Lu J. Rapid detection of avian influenza a virus (H7N9) by lateral flow dipstick-recombinase polymerase amplification. Biol Pharm Bull. 2018;41(12):1804-8.
Tan KK, Azizan NS, Yaacob CN, Seri NA, Samsudin NI, Teoh BT, Sam SS, AbuBakar S. Operational utility of the reverse-transcription recombinase polymerase amplification for detection of dengue virus. BMC Infect Dis. 2018;18(1):169.
Vasileva Wand NI, Bonney LC, Watson RJ, Graham V, Hewson R. Point-of-care diagnostic assay for the detection of Zika virus using the recombinase polymerase amplification method. J Gen Virol. 2018;99(8):1012–26.
Ge Y, Wu B, Qi X, Zhao K, Guo X, Zhu Y, Qi Y, Shi Z, Zhou M, Wang H, et al. Rapid and sensitive detection of novel avian-origin influenza a (H7N9) virus by reverse transcription loop-mediated isothermal amplification combined with a lateral-flow device. PLoS One. 2013;8(8):e69941.
Yang Y, Qin X, Wang G, Jin J, Shang Y, Zhang Z. Development of an isothermoal amplification-based assay for rapid visual detection of an Orf virus. Virol J. 2016;13:46.
Escadafal C, Faye O, Sall AA, Faye O, Weidmann M, Strohmeier O, von Stetten F, Drexler J, Eberhard M, Niedrig M, et al. Rapid molecular assays for the detection of yellow fever virus in low-resource settings. PLoS Negl Trop Dis. 2014;8(3):e2730.
Kunze A, Dilcher M, Abd El Wahed A, Hufert F, Niessner R, Seidel M. On-Chip isothermal nucleic acid amplification on flow-based Chemiluminescence microarray analysis platform for the detection of viruses and Bacteria. Anal Chem. 2016;88(1):898–905.
Kersting S, Rausch V, Bier FF, von Nickisch-Rosenegk M. Rapid detection of plasmodium falciparum with isothermal recombinase polymerase amplification and lateral flow analysis. Malar J. 2014;13:99.
This study was funded by the National Major Science & Technology Projects for Infectious Disease Control and Prevention (2017ZX10103008), Jiangsu Province Science & Technology Demonstration Project for Emerging Infectious Diseases Control and Prevention (No. BE2015714), Jiangsu Provincial Key Medical Discipline (ZDXKA2016008), the Natural Science Foundation of Jiangsu Province (BK20161583), and the “333” Projects of Jiangsu Province (BRA2017552). The funders contributed to the study design and data collection. The funders had no roles in the study design and data collection.
Availability of data and materials
All data generated or analyzed during current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
All procedures conducted in this study involving human participants and materials were approved by the Jiangsu Provincial Center for Disease Control Ethics Committee, reference number: JSCDCLL2013014. Written consent was obtained from their parents or grandparents (when parents worked outside the hometown, grandparents were legal guardians according to the law of China).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.