Ultra-sensitive chemiluminescence imaging DNA hybridization method in the detection of mosquito-borne viruses and parasites
Parasites & Vectors volume 10, Article number: 44 (2017)
Mosquito-borne viruses (MBVs) and parasites (MBPs) are transmitted through hematophagous arthropods-mosquitoes to homoiothermous vertebrates. This study aims at developing a detection method to monitor the spread of mosquito-borne diseases to new areas and diagnose the infections caused by MBVs and MBPs.
In this assay, an ultra-sensitive chemiluminescence (CL) detection method was developed and used to simultaneously detect 19 common MBVs and MBPs. In vitro transcript RNA, virus-like particles (VLPs), and plasmids were established as positive or limit of detection (LOD) reference materials.
MBVs and MBPs could be genotyped with high sensitivity and specificity. The cut-off values of probes were calculated. The absolute LODs of this strategy to detect serially diluted in vitro transcribed RNAs of MBVs and serially diluted plasmids of MBPs were 102–103copies/μl and 101–102copies/μl, respectively. Further, the LOD of detecting a strain of pre-quantified JEV was 101.8–100.8PFU/ml, fitted well in a linear regression model (coefficient of determination = 0.9678).
Ultra-sensitive CL imaging DNA hybridization was developed and could simultaneously detect various MBVs and MBPs. The method described here has the potential to provide considerable labor savings due to its ability to screen for 19 mosquito-borne pathogens simultaneously.
Mosquito-borne viruses (MBVs) are transmitted from hematophagous arthropod/mosquitoes to homoiothermous vertebrates. As competent hosts, arthropods can support large virus inocula, and transmit them from an infected donor to a recipient during blood-feeding. Three major families in MBVs include the Togaviridae, Flaviviridae, and Bunyaviridae . The globally pandemic MBV species encompass dengue virus (DENV), Japanese encephalitis virus (JEV), yellow fever virus (YFV), West Nile virus (WNV)  and chikungunya virus (CHIKV) . In China, MBVs severely affect people living in tropical regions and rural areas. There are more than 70 identified flaviviruses, of which 35 can affect humans . Factors like global warming, increased population mobility, and urbanization influence the epidemic scope of MBVs and sometimes even result in sudden outbreak . In 1999, the outbreak of WNV in New York resulted in mass-scale animal mortality . In 2012, the outbreak of WNV in Texas resulted in 1,868 human infections (89 deaths) and an economic cost of more than $47 million . In 2010, the outbreak of CHIKV in Guangdong province of China led to 173 cases of human infections . From 2005 to 2006, more than 1.3 million cases of CHIKV infections were reported in India . The infected cases of DENV recorded by the World Health Organization (WHO) were approximately 0.9 million . However, actual infected cases worldwide are likely to be more than 390 million; 96 million infected patients show the clinical or sub-clinical symptoms of DENV infection . Each year, more than 20,000 casualties are reported due to severe infections of DENV .
In addition, malaria parasites and lymphatic filariae, borne by mosquitoes, are currently the most common parasites in the world [13, 14]. Malaria and lymphatic filariasis have been considered as the world’s highest priority control and eradication parasitic diseases due to their wide prevalence [15, 16]. In tropical regions, these two parasites share the same propagation host and may also co-infect humans . Mosquito-borne parasites (MBPs) infecting humans include Plasmodium vivax, Plasmodium malariae, Plasmodium falciparum and Plasmodium ovale. Malaria is rampant in Asia, Africa, and Latin America, with half of the world’s population living in malaria-risk areas. In 2006, 250 million people were infected by malaria and approximately 1 million deaths were reported, among which children (<5 years) and pregnant women were most affected . Mosquito-borne lymphatic filariae include Wuchereria bancrofti, Brugia malayi and Brugia timori. Approximately 1.3 billion people are living in the lymphatic filariasis-risk areas. Following malaria, lymphatic filariasis is the second most prevalent arthropod-borne disease worldwide, which affects 128 million people each year in Asia, Africa, Western Pacific, and parts of America .
The spread of the mosquito-borne pathogens from the epidemic areas to new ones poses a serious threat to the monitoring, prevention, and control of disease . Efficient prevention, control, and treatment of mosquito-borne diseases rely on the rapid detection of the mosquito-borne pathogens and accurate diagnosis. Currently, methods for detecting MBVs and MBPs are mainly based on genetic assays, which include real-time PCR [21–25], isothermal amplification , oligonucleotide microarray  and sequencing . However, these methods were designed for detecting only certain species of MBVs or MBPs. A method that can simultaneously detect common MBVs and MBPs has not yet been reported. In this assay, a reliable and portable, ultra-sensitive chemiluminescence (CL) imaging DNA hybridization was developed and used to simultaneously detect common MBVs and MBPs. In vitro transcript RNA, virus-like particles, and plasmids were established as positive or limit of detection (LOD) reference materials. The specificity and sensitivity of the method was validated in cultivated MBVs, simulated mosquito samples, reference materials, and non-MBVs.
Specimen collection and processing
Cultivated MBVs strains were obtained from the Institute of Microbiology Epidemiology of the Academy of Military Medical Sciences and National Institute for Viral Disease Control and Prevention, China CDC. Lysates (TIANamp Virus RNA Kit, Tiangen Biotech Beijing Co., Ltd.) were added to all the viral samples prior to experiments. Total RNAs were extracted using TIANamp Virus RNA Kit as described in the manufacturer’s protocol (Tiangen Biotech Beijing, China) and stored at -70 °C until use. Uninfected Aedes vigilax were obtained from the Institute of Microbiology Epidemiology of Academy of Military Medical Sciences and stored at -70 °C until use.
Primers and probes design and evaluation
In this assay, 12 MBVs and 7 MBPs were chosen for detection: DENV1-4, JEV, YFV, WNV, St. Louis encephalitis virus (SLEV), CHIKV, eastern equine encephalitis virus (EEEV), western equine encephalitis virus (WEEV), Venezuelan equine encephalitis virus (VEEV), P. vivax, P. malariae, P. falciparum, P. ovale, W. bancrofti, B. malayi and B. timori. Complete genomic sequences of these MBVs, 18S rRNA gene sequences of malaria parasites, SspI repeat DNA sequence of W. bancrofti, and Hhal repeat region sequences of B. malayi and B. timori were downloaded from the GenBank database. The sequences were aligned using AlignX (a component of the Vector NTI Advance 10.3.0) to compare the homology between potential targets belong to the same genus. Gene-specific primers (GSPs) were derived from published literature or designed from conserved regions using Primer Premier 5 (PREMIER Biosoft International, Palo Alto, USA). Reverse primers used in Multiplex RT-PCR or PCR were biotin-labelled at 5′ end to form a CL reaction. The capture probes, which were used to capture the target DNA amplification fragments, were designed by Oligo 7 (Molecular Biology Insights, Inc., Colorado Springs, USA). Continuous homologous sequences of congeners were excluded. A human-original sequence was used as negative control probe. A repeat sequence of (T)12 with an amino-labeled 3′-end was connected to the 3′-end of all the probes to fix the repeat sequence with the aldehyde-chip surface. A repeat sequence of (T)20 with a biotin-labeled 5′-end and an amino-labeled 3′-end served as a quality control (QC) probe. All the primers and probes were confirmed by NCBI BLAST and synthesized by Sangon Biotech Co., Ltd. (Shanghai, China).
To evaluate the efficiency of primers, serially diluted viral RNAs or plasmid DNAs were amplified by RT-PCR or PCR and resolved by 2% agarose gel electrophoresis. Furthermore, the most efficient pairs of primers were biotin labeled in reverse primers, and RT-PCR or PCR was used to evaluate the DNA hybridization efficiency of the capture probes.
Preparation of in vitro transcript RNA
To verify the LOD of this strategy and calibrate the concentration of virus-like particles (VLPs), in vitro transcribed RNAs of these MBVs were prepared. First, plasmids were constructed by cloning target gene fragments to PGM-T vector as described in the manufacturer’s protocol (Tiangen Biotech Beijing Co., Ltd.). Following digestion with FastDigest Sal I restriction endonuclease (ThermoFisher, Waltham, USA), the plasmids were purified by agarose gel and retrieved. Subsequently, in vitro transcribed RNAs were prepared with TranscriptAid T7 High Yield Transcription Kit (ThermoFisher), digested by deoxyribonuclease (DNase) I, and purified by phenol-chloroform extraction method as described in the manufacturer’s protocol (ThermoFisher). Further, in vitro transcribed RNAs were purified and digested by DNase I using RNeasy Mini Kit (QIAGEN, Mississauga, Canada) as described in the manufacturer’s protocol. The purified RNAs were mixed with 2 × RNA Loading Dye Solution (ThermoFisher) and resolved by 2% agarose gel electrophoresis. Finally, in vitro transcribed RNAs were quantified by UV spectrophotometry and subject to serial dilutions at a magnitude of 10 and stored at -70 °C until use.
Preparation of virus-like particles
The maturation protein gene and coat protein gene of MS2 bacteriophage were amplified by MS2-BamH I and MS2-Hind III primers (modified from reference literature) , and the sequence was as follows: MS2-BamH I, 5′-AAC GGA TCC CAT GGC TAT CGC TGT AGG TAG CC-3′, MS2-Hind III, 5′-CAT AAG CTT CTT CGA CAT GGG TAA TCC TCA TGT T-3′. The bases in italics represent endonuclease sites of BamH I and Hind III. The PET-MS2 vector was constructed using the MS2 amplified fragment and PET-28a (+) vector (Novagen, Merck KGaA, Darmstadt, Germany) and ligated using T4 DNA ligase (New England BioLabs® Inc., Ipswich, USA) after digestion with restriction enzymes (BamH I and Hind III, ThermoFisher), according to the manufacturer’s protocol. The Hind III restriction enzyme site and protection bases were added to the 5′ end of MBVs’ GSPs (both forward and reverse), and target genes of these mosquito-borne viruses were amplified and digested by Hind III restriction enzyme. PET-MS2 vectors were also digested by Hind III and dephosphorylated by alkaline phosphatase (New England BioLabs® Inc.). Target gene fragments of MBVs and PET-MS2 vectors were ligated using T4 DNA ligase and cloned to BL21 (DE3) competent cells (Tiangen Biotech Beijing Co., Ltd.), according to the manufacturer’s instructions. These bacterial plasmids were induced by isopropyl-1-thio-β-galactopyranoside (IPTG) to express VLPs and lysed by sonication. The VLPs were precipitated by PEG8000 and purified by 15–45% (w/v) sucrose density gradient ultracentrifugation at 32,000× rpm for 6 h/4 °C (Beckman, MAL-80 ROTOR, Brea, USA). The aliquot containing VLPs was dialyzed against PBS at 4 °C for 24 h to remove sucrose. The purified VLPs were used as positive RNA references materials.
Multiplex asymmetric RT-PCR and PCR amplification
The GSPs of MBVs and MBPs were used in two multiplex asymmetric RT-PCR systems and one multiplex asymmetric PCR system, respectively. For optimal amplification efficiency, orthogonal experiment was used to optimize the concentration of the GSPs. Briefly, the concentration ratio of the forward and reverse primers were set to 1:5. The high, medium, and low concentrations of primers were set based on the results of preliminary experiments. The selected primer concentrations were determined by evaluation of the amplification efficiencies of the target genes. Each RT-PCR was performed in a 25 μl reaction volume, containing 12.5 μl of 2× One Step Buffer, 1.0 μl of PrimeScript One-step Enzyme Mix (DRR055A, Takara Biotechnology (Dalian) Co., Ltd., Dalian, China), 5 μl of total RNA template, and specific primer mix. RT-PCR system A amplified the target genes of mosquito-borne flavivirus included DENV, JEV, YFV, SLEV, and WNV. System A contained GSPs of DC10418 (0.2 μmol), JEV-3UTR-F2 (0.28 μmol), YFV-3UTR-F1 (0.4 μmol), WNV-3UTR-F2 (0.2 μmol), SLE-3UTR-F1 (0.1 μmol), and CDC10564 (2 μmol). RT-PCR system B amplified the target genes of mosquito-borne alphavirus included EEEV, WEEV, VEEV, and CHIKV. System B contained GSPs of EEE-E1-F (0.4 μmol), EEE-E1-R (2 μmol), CHIKV-F1 (0.2 μmol), CHIKV-R2 (1 μmol), WEE-F (0.08 μmol), WEE-R (0.4 μmol), VEEV-cap-F3 (0.12 μmol), and VEEV-cap-R3 (0.6 μmol). PCR was performed in a 25 μl reaction volume containing 2.5 μl of 10× EX Taq Buffer, 0.5 μl of TaKaRa Ex Taq HS (RR006, Takara Biotechnology (Dalian) Co., Ltd.), 200 μmol of each dNTP, 5 μl of total DNA template, and specific primer mix. PCR system amplified target genes of plasmodium and lymphatic filariae. PCR system contained GSPs of MAL-F1 (0.16 μmol), MAL-R1 (0.8 μmol), BM/WBF (0.08 μmol), WBR (0.4 μM), MGB-HhaI-For (0.16 μmol), and MGB-HhaI-Rev (0.8 μmol). Amplifications were performed on a Veriti 96-Well Thermal Cycler PCR system (Applied Biosystems, Foster City, USA). RT-PCR was performed under the following conditions: 30 min at 50 °C; 2 min at 94 °C; 45 cycles of 20 s at 94 °C, 20 s at 55 °C, and 20 s at 72 °C; and a final extension of 5 min at 72 °C. PCR was performed under the following conditions: 2 min at 95 °C; 45 cycles of 20 s at 94 °C, 20 s at 55 °C, and 20 s at 72 °C; and a final extension of 5 min at 72 °C.
DNA Hybridization and signal detection
Captured probes of MBVs and MBPs (50 μmol) were used to establish the MBVs and MBPs capture-chip, and the spotting was repeated thrice in the vertical direction on a aldehyde modified glass slide (Baio Technology Shanghai Co., Ltd., Shanghai, China) surface with uniform proportional printing buffer (5% glycerol, 0.1% sodium dodecyl sulfate (SDS), 6× saline-sodium citrate buffer (SSC), and 2% (wt/vol) Ficoll 400) as described in previous studies . The QC probe, used at 12.5 μmol final concentration, was spotted and repeated six to seven times in the horizontal direction to calibrate the CL signal values. Each aldehyde slide was divided into 10 blocks (11 × 11 mm) by a waterproof film to detect 10 different samples.
After amplification, the two multiplex RT-PCR reaction products, amplified using the same template, were blended. After 5 min of denaturation at 95 °C, the 5 μl of RT-PCR mixtures or PCR products were immediately placed on ice for 5 min and mixed with 5 μl of DNA hybridization buffer (8× SSC, 0.6% SDS, 10% formylamine, and 10× Denhardt). DNA hybridization mixtures (10 μl) were hybridized on the MBVs or MBPs capture-chip for 1 h at 45 °C. The capture-chip was washed for 30 s each with 1× SSC and 0.2% SDS, 0.2× SSC, and 0.1× SSC at room temperature. The capture-chip was incubated with 15 μl of streptavidin-horseradish peroxidase (Str-HRP, Sigma-Aldrich, St. Louis, USA) for 25 min at 37 °C. Subsequently, the capture-chip was washed with PBST (phosphate buffer, 0.05% Tween 20) 10 s at room temperature. Finally, 10 μl of pre-mixed CL HRP substrate luminal solution and H2O2 (Millipore Corporation, Boston, USA) were added to the capture-chip and detected immediately with a micro-light level imaging system (developed in-house, patent number CN201330010609) (The working principle of this CL imaging DNA hybridization method were showed in Fig. 2c).
Determination of cut-off values of capture probes
In this assay, low background CL imaging and self-programming software were employed as the imaging method and the interpretation software was used to read the grey values of the probes. The grey values recorded by the software were set to be in the range of 0–256. The CL signal values of probes = mean of signal values of probes - signal values of background. The CL signal value was calibrated as follows: the calibrated CL signal value of a probe = mean of CL signal value of the probe / mean of CL signal value of QC probe (the same detection block) × 100. The cut-off values of probes were calculated as follows: cut-off value of a probe = average signal values of to detect other pathogens (except the pathogen detected by this probe) + 3 × standard deviation. Probe signal was considered positive when the value was greater than the cut-off value of the probe.
Specificity and sensitivity
Due to the lack of MBVs and MBPs samples, the specificity of the strategy was evaluated only by 6 cultivated MBVs strains and an inactivated JEV Vaccine (Table 1). The other MBVs and MBPs were evaluated by artificial VLPs, in vitro transcribed RNAs, or plasmids. Further, a panel of non-MBVs, which include tick-borne encephalitis virus (TBEV), avian influenza A (H5N1), 2009 influenza A (H1N1) (PH1N1), seasonal influenza A (H3N2), seasonal influenza A (H1N1), Hepatitis B virus (HBV), Hepatitis C virus (HCV), adenovirus AD2, AD3, AD40, rubella virus, mumps virus, and respiratory syncytial viruses (RSV) HK6 and B, were also used to determine the specificity of the strategy.
To evaluate the LOD of this strategy, in vitro transcribed RNAs or plasmids were quantified using UV spectrophotometer, and the copy numbers were calculated, serially diluted to 10-fold, and were detected. Purified VLPs were also diluted and detected to determine the relative concentrations as compared to in vitro transcribed RNA.
Simulated detection of mosquito samples
To test simulated mosquito samples, uninfected mosquito pools of 50 Aedes vigilax were spiked with 50 μl 106 copies/μl part of single or mixed VLPs or plasmids. Then the spiked mosquito pools were homogenized in 500 μl diluent containing 20% fetal bovine serum in phosphate-buffered saline using T10 basic ULTRA-TURRAX® small dispersing instrument (IKA, Staufen, Germany). The homogenized samples were centrifuged for 5 min at 12,000× rpm. Supernatants were removed to clean microfuge tubes. RNA or DNA was extracted from the supernatant using TIANamp Virus RNA Kit (Tiangen Biotech Beijing Co., Ltd.) or TIANamp Genomic DNA Kit (Tiangen Biotech Beijing Co., Ltd.), respectively, as described in the manufacturer’s protocol. Subsequently, the extracted RNAs or DNAs were amplified by RT-PCR or PCR systems and detected by the CL imaging DNA hybridization assay.
Primers and probes design and evaluation
The amplification efficiency of 81 GSPs, obtained from published literatures, and 82 GSPs designed for this assay were precisely evaluated by serially diluted viral RNAs or plasmid DNAs. Biotin-labeled reverse primers were used for multiplex asymmetric RT-PCR or PCR to evaluate DNA hybridization efficiency of 223 capture probes. Finally, 20 GSPs and 21 probes were chosen to amplify the MBVs and MBPs targeted genes and capture the target DNA amplification fragments (Tables 2 and 3).
Construction of in vitro transcribed RNAs and VLPs
Nine in vitro transcribed RNAs of MBVs were prepared to determine the LOD of the MBVs detection assay and calibrate the concentration of VLPs. After purification, in vitro transcribed RNAs were verified by agarose gel electrophoresis (Fig. 1a) and quantified by UV spectrophotometry. The copy numbers of these in vitro transcribed RNAs were calculated as: copy numbers (copies/μl) = RNA concentration (ng/μl) × Avogadro constant × 10-9/ (RNA length × 340 g/mol) (Table 4). The in vitro transcribed RNAs were serial diluted to 10-fold and stored at -70 °C until use.
Twelve VLPs containing the amplified target RNAs were also prepared as the positive RNA references of the MBVs detection strategy. RNase A and DNase I digestion on the prepared VLPs showed that the nucleic acid bands between 1,000 and 2,000 bp tolerated both enzymes due to their packaging in the internal section of VLPs (Fig. 1b). To determine whether the target RNAs had been packaged in the internal regions, RNase A and DNase I digested VLPs were serially diluted (10-fold), total RNA was extracted, amplified by both multiplex RT-PCR and PCR (both based on the GSPs of RT-PCR system A or B), and detected by MBVs hybridization capture-chip. The results showed that DNase I did not completely digested the DNA templates from high concentration VLPs; the extractions from low-concentration VLPs could only be detected by RT-PCR but not PCR (Fig. 2b), thereby demonstrating that the DNA templates had been diluted to sufficiently low concentrations that cannot be detected.
Determination of cut-off value of capture probes
In this assay, the cut-off values of a probe were calculated as follows, cut-off = average signal value to detect other pathogens (except the particular pathogen detected by this probe) + 3 × standard deviation. The cut-off values of all the probes are shown in Table 5.
Evaluation of specificity
The specificity of the strategy was evaluated using cultivated MBVs strains, inactivated JEV vaccine strain, artificial VLPs, in vitro transcribed RNAs, and plasmids. The results showed that this assay genotyped MBVs and MBPs with no cross-reactions between the samples and reference materials (Fig. 3a, b). The CL signal value (i.e. that had not been calibrated) of each probe was shown in 3-D bar chart (Fig. 3c, d). Furthermore, a panel of non-MBVs did not yield any positive results with this strategy, suggesting the method is indeed highly specific (Additional file 1: Figure S1).
Determination of LOD
To determine the absolute LOD of this strategy, serially diluted in vitro transcribed RNAs or plasmids, which had been quantified by the copy numbers, served as the reference materials. The absolute LODs for the MBVs were 102 copies/μl except for WNV (103 copies/μl), (Fig. 4a). The calibrated CL signal values of the in vitro transcribed RNAs are described in a line chart (Fig. 4c). Further, the absolute LODs for P. falciparum, P. ovale, and P. malariae are 101 copies/μl; for the others, the values are 102 copies/μl (Fig. 5a). The calibrated CL values of these plasmid DNAs are also described in a line chart (Fig. 5b)
A strain of JEV virus derived from National Institute for Viral Disease Control and Prevention was also employed as a reference material to determinate the LOD of this strategy. The results showed that the LOD for JEV was between 101.8 PFU/ml and 100.8 PFU/ml (Fig. 4b). A linear regression model based the calibrated CL signal values against the JEV concentrations shows satisfactory fitness, with the coefficient of determination being 0.9678 (Fig. 4d).
Simulated detection of mosquito samples
To test simulated mosquito samples, uninfected mosquito pools of 50 Aedes vigilax were spiked with 106 copies/μl of six single VLPs and plasmids or mixed VLPs and plasmids, followed by the CL imaging DNA hybridization assay. The results showed that this assay could detect all mimic samples of single MBVs, single MBPs, mixed mosquito-borne flavivirus, mixed mosquito-borne alphavirus, and mixed MBPs (Fig. 6). Therefore, the proposed assay can potentially test targeted MBVs and MBPs in actual mosquito samples.
As a competent vector, arthropods can support the sufficiently large inoculum of the virus and transmit the MBVs from an infected donor to a recipient host during blood-feeding. MBVs are also transmitted between male and female mosquitoes or from a female to the offspring. Furthermore, malaria and lymphatic filariasis, transmitted by hematophagous arthropods, are still the most prevalent mosquito-borne parasitic diseases in humans [13, 14]. A continuous effort to monitor the MBVs is critical to prevent sudden outbreaks . Many gene-based diagnostic methods have been used to detect MBVs and MBPs. However, to the best of our knowledge, a portable, ultra-sensitive CL DNA hybridization method to simultaneously detect common MBVs and MBPs has not been reported. Because of the lack of antiviral drugs [10, 31, 32], the use of vaccines has emerged as an efficient strategy to control the spread of MBV diseases . At present, human vaccines for flavivirus and alphavirus are available; however, these only include prevention of YFV, JEV, and TBEV [34–38], while the human vaccines for the prevention of DENV, WNV, CHIKV, SLEV, WEEV, EEEV and VEEV are still not approved for clinical trials [33, 39]. As a result, it has become imperative to develop reliable methods for rapid diagnosis, timely prevention, and efficient control of the mosquito-borne diseases. Rapid diagnosis is not only the indispensable method for controlling the infectious diseases, but also the prerequisite for timely prevention and efficient treatment. Infectious diseases are only able to be effectively prevented and treated after rapid and accurate diagnosis. Therefore, the spread of these infectious diseases, for which the drugs are not yet established, can be effectively reduced by controlling their vectors and accurate monitoring. In this assay, we developed a portable, ultra-sensitive CL imaging DNA hybridization method to detect various MBVs and MBPs simultaneously. This method is based on horseradish peroxidase catalyzing luminol-H2O2, and its performance proved to be highly sensitive and specific for the 12 MBVs and 7 MBPs that were included in the current study. The traditional epidemic of these pathogens involves many tropical and sub-tropical areas of Asia, Africa, the Americans, Europe and the Pacific archipelago [40–43]. Currently, there may not be a region in the world where all 19 pathogens are present. So the inclusion of the 19 pathogens in the current study appears exceeding the needs of clinical detections. However, we believe our study delivers a proof-of-concept study where such a diagnostic assay may test for more pathogens and save significant labour and time when a patient presents with symptoms of unknown etiology. Furthermore, given the fact that mosquito-borne diseases may occur in non-endemic areas due to the uncertainties of global human migration, for example, the first imported cases of Zika virus in China and Europe [44, 45], a rapid and ready-to-be-applied method will have significant clinical implications. For point of care testing, a self-powered portable CL CCD imager was developed. A self-programming software was employed to interpret the software readings and to study the grey values of the probes. However, other commercial CL imagers based on CCD imaging technology can also be used for CL imaging. Despite needing a CL imager, this CL imaging method has lower cost and faster detection speed than our previous visual method which was based on quantum dot-catalyzed silver deposition .
In the current study, due to the lack of pathogen nucleic acid samples, we validated the specificity of our method using the positive and negative references. The absolute LOD, which was evaluated using serially diluted in vitro transcribed RNAs and plasmids, was in the range of 101–103copies/μl. A strain of JEV was also employed as a reference material to determine the absolute LOD. Results from the two repeat tests showed that the absolute LOD for JEV was in the range of 100.8–101.8PFU/ml, and fits well in a linear regression model, demonstrating the feasibility of quantitative measurement of viral load. The assay is rapid, portable, ultra-sensitive, and high-throughput. The entire detection time, spanning from sample extraction to result reading, is between 6–8 h.
Since most of the targets in this assay are potent pathogens, we did not collect actual samples or all the nucleic acids. Instead, 14 artificial target pathogen gene sequences were constructed by overlap extension method and cloned into plasmid vectors as well as those target pathogens having nucleic acids samples (datas not shown). Subsequently, RNAs of these target viruses were in vitro transcribed from these plasmids. We observed that DNA templates from the plasmids were not completely removed by DNase I digestion. Consequently, RNeasy Mini Kit (QIAGEN), a purification method based on spin column, was also employed to purify and digest these in vitro transcribed RNAs on column. Although residual DNA templates were not completely removed from in vitro transcribed RNAs, they were negligible after the concentrations of in vitro transcribed RNAs were diluted to lower than 108 copies/μl due to the twice DNase I digestion (dates not shown). Similar issue was observed during the preparation of VLPs of these target MBVs. However, to maintain the integrity, the VLPs were not digested by DNase I on the spin column. Consequently, the VLPs, as RNA positive reference materials, were also diluted to a sufficiently low concentration to overcome the influence of residual DNA templates as the concentration of residual DNAs were only 103–104-fold lower than that of the RNAs within VLPs.
In the present assay, a uniform reverse primer located in the 3′UTR region of flavivirus genome was employed in the multiplex asymmetric RT-PCR. Consequently, relatively shorter probes were used for flavivirus to eliminate cross signals among flavivirus species due to the homology of the genomic sequences of the 3′UTRs. On the other hand, our method did not distinguish between B. malayi and B. timori because the target sequences (i.e. HhaIrepeat sequence) are almost identical. Furthermore, methods based on DNA hybridization inevitably appeared cross signal with off target capture probes, especially for the samples containing extremely high concentrations of the target pathogens. Consequently, the cut-off values of capture probes showed in Table 5 need to be adjusted and refined to test the different concentrations of the target pathogen for future application. At present, there are many commercial extraction kits available to extract genomic RNA and DNA from diverse origins, so the proposed strategy can potentially detect these MBVs and MBPs regardless of the origins of these pathogens. In the proof-of-concept experiment, we spiked uninfected mosquito pools with VLPs or plasmids, demonstrated that our CL imaging DNA hybridization method could detect target MBVs and MBPs in mosquito samples. Therefore, further research should be focused on extending the method into the clinical settings for the epidemiological investigations of MBVs and MBPs.
Ultra-sensitive CL imaging DNA hybridization was developed and could simultaneously detect MBVs and MBPs. This method is based on horseradish peroxidase catalyzing luminol-H2O2, and its performance proved to be rapid, portable, high-throughput, highly sensitive and specific. The entire operation time, spanning from sample extraction to result reading, is between 6–8 h. The method described here has the potential to provide considerable labor savings due to its ability to screen for 19 mosquito-borne pathogens simultaneously.
Eastern equine encephalitis virus
Hepatitis B virus
Hepatitis C virus
Japanese encephalitis virus
Limit of detection
Respiratory syncytial viruses
Sodium dodecyl sulfate
St. Louis encephalitis virus
Saline-sodium citrate buffer
Tick-borne encephalitis virus
Venezuelan equine encephalitis virus
Western equine encephalitis virus
World Health Organization
West Nile virus
Yellow fever virus
Hubalek Z. Mosquito-borne viruses in Europe. Parasitol Res. 2008;103 Suppl 1:S29–43.
Weissenbock H, Hubalek Z, Bakonyi T, Nowotny N. Zoonotic mosquito-borne flaviviruses: worldwide presence of agents with proven pathogenicity and potential candidates of future emerging diseases. Vet Microbiol. 2010;140(3–4):271–80.
Ho PS, Ng MM, Chu JJ. Establishment of one-step SYBR green-based real time-PCR assay for rapid detection and quantification of chikungunya virus infection. Virol J. 2010;7:13.
Palani G, Padmanabhan PP, Ramesh K, Asadullah KS, Sambasivam M, Arunagiri K, Krishnasamy K. Development and evaluation of flavi-immunoglobulin M capture enzyme-linked immunosorbent assay. Indian J Pathol Micr. 2013;56(3):269–71.
Gould EA, Higgs S. Impact of climate change and other factors on emerging arbovirus diseases. T Roy Soc Trop Med H. 2009;103(2):109–21.
Lanciotti RS, Roehrig JT, Deubel V, Smith J, Parker M, Steele K, et al. Origin of the West Nile virus responsible for an outbreak of encephalitis in the northeastern United States. Science. 1999;286(5448):2333–7.
Murray KO, Ruktanonchai D, Hesalroad D, Fonken E, Nolan MS. West Nile virus, Texas, USA, 2012. Emerg Infect Dis. 2013;19(11):1836–8.
Wu D, Wu J, Zhang Q, Zhong H, Ke C, Deng X, et al. Chikungunya outbreak in Guangdong Province, China, 2010. Emerg Infect Dis. 2012;18(3):493–5.
Arankalle VA, Shrivastava S, Cherian S, Gunjikar RS, Walimbe AM, Jadhav SM, et al. Genetic divergence of Chikungunya viruses in India (1963–2006) with special reference to the 2005–2006 explosive epidemic. J Gen Virol. 2007;88(Pt 7):1967–76.
World Health Organization and the Special Program for Research and Training in Tropical Diseases (TDR). Dengue Guidelines for Diagnosis, Treatment, Prevention and Control: New Edition. 2009.
Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496(7446):504–7.
Gubler DJ. The economic burden of dengue. Am J Trop Med Hyg. 2012;86(5):743–4.
Birkholtz LM, Bornman R, Focke W, Mutero C, de Jager C. Sustainable malaria control: transdisciplinary approaches for translational applications. Malaria J. 2012;11:431.
Manguin S, Bangs MJ, Pothikasikorn J, Chareonviriyaphap T. Review on global co-transmission of human Plasmodium species and Wuchereria bancrofti by Anopheles mosquitoes. Infect Genet Evol. 2010;10(2):159–77.
Kyelem D, Biswas G, Bockarie MJ, Bradley MH, El-Setouhy M, Fischer PU, et al. Determinants of success in national programs to eliminate lymphatic filariasis: a perspective identifying essential elements and research needs. Am J Trop Med Hyg. 2008;79(4):480–4.
Molyneux DH, Zagaria N. Lymphatic filariasis elimination: progress in global programme development. Ann Trop Med Parasit. 2002;96 Suppl 2:S15–40.
Buck AA, Anderson RI, MacRae AA. Epidemiology of poly-parasitism. I. Occurrence, frequency and distribution of multiple infections in rural communities in Chad, Peru, Afghanistan, and Zaire. Tropenmed Parasitol. 1978;29(1):61–70.
World Malaria Report 2008. http://www.who.int/malaria/publications/atoz/9789241563697/en/. Accessed 1 Jan 2008.
Lymphatic filariasis. http://www.who.int/mediacentre/factsheets/fs102/en/. Accessed 1 Feb 2016.
Gould EA, Higgs S, Buckley A, Gritsun TS. Potential arbovirus emergence and implications for the United Kingdom. Emerg Infect Dis. 2006;12(4):549–55.
Chien LJ, Liao TL, Shu PY, Huang JH, Gubler DJ, Chang GJ. Development of real-time reverse transcriptase PCR assays to detect and serotype dengue viruses. J Clin Microbiol. 2006;44(4):1295–304.
Hull R, Nattanmai S, Kramer LD, Bernard KA, Tavakoli NP. A duplex real-time reverse transcriptase polymerase chain reaction assay for the detection of St. Louis encephalitis and eastern equine encephalitis viruses. Diagn Micr Infec Dis. 2008;62(3):272–9.
Kang X, Li Y, Liu H, Lin F, Cai X, Sun T, et al. A duplex real-time reverse transcriptase polymerase chain reaction assay for detecting western equine and eastern equine encephalitis viruses. Virol J. 2010;7:284.
Mishra K, Raj DK, Hazra RK, Dash AP, Supakar PC. The development and evaluation of a single step multiplex PCR method for simultaneous detection of Brugia malayi and Wuchereria bancrofti. Mol Cell Probe. 2007;21(5–6):355–62.
Rao RU, Weil GJ, Fischer K, Supali T, Fischer P. Detection of Brugia parasite DNA in human blood by real-time PCR. J Clin Microbiol. 2006;44(11):3887–93.
Lanciotti RS, Kerst AJ. Nucleic acid sequence-based amplification assays for rapid detection of West Nile and St. Louis encephalitis viruses. J Clin Microbiol. 2001;39(12):4506–13.
Xiao-Ping K, Yong-Qiang L, Qing-Ge S, Hong L, Qing-Yu Z, Yin-Hui Y. Development of a consensus microarray method for identification of some highly pathogenic viruses. J Med Virol. 2009;81(11):1945–50.
Hu Z, Nordstrom H, Nowotny N, Falk KI, Sandstrom G. Anchored pan dengue RT-PCR and fast sanger sequencing for detection of dengue RNA in human serum. J Med Virol. 2010;82(10):1701–10.
Huang J, Yang CM, Wang LN, Meng S, Deng W, Li JM. A novel real-time multiplex reverse transcriptase-polymerase chain reaction for the detection of HIV-1 RNA by using dual-specific armored RNA as internal control. Intervirology. 2008;51(1):42–9.
Zhang Y, Liu Q, Wang D, Chen S, Wang S. Simultaneous detection of oseltamivir- and amantadine-resistant influenza by oligonucleotide microarray visualization. PLoS One. 2013;8(2):e57154.
Plummer EM, Shresta S. Mouse models for dengue vaccines and antivirals. J Immunol Methods. 2014;410:34–8.
Reichert E, Clase A, Bacetty A, Larsen J. Alphavirus antiviral drug development: scientific gap analysis and prospective research areas. Biosecur Bioterror. 2009;7(4):413–27.
Metz SW, Pijlman GP. Arbovirus vaccines; opportunities for the baculovirus-insect cell expression system. J Invertebr Pathol. 2011;107(Suppl):S16–30.
Frierson JG. The yellow fever vaccine: a history. Yale J Biol Med. 2010;83(2):77–85.
Camacho LA, Freire Mda S, Leal Mda L, Aguiar SG, Nascimento JP, Iguchi T, et al. Immunogenicity of WHO-17D and Brazilian 17DD yellow fever vaccines: a randomized trial. Rev Saude Publ. 2004;38(5):671–8.
Hombach J, Solomon T, Kurane I, Jacobson J, Wood D. Report on a WHO consultation on immunological endpoints for evaluation of new Japanese encephalitis vaccines, WHO, Geneva, 2–3 September, 2004. Vaccine. 2005;23(45):5205–11.
Hoke CH, Nisalak A, Sangawhipa N, Jatanasen S, Laorakapongse T, Innis BL, et al. Protection against Japanese encephalitis by inactivated vaccines. New Engl J Med. 1988;319(10):608–14.
Amicizia D, Domnich A, Panatto D, Lai PL, Cristina ML, Avio U, Gasparini R. Epidemiology of tick-borne encephalitis (TBE) in Europe and its prevention by available vaccines. Hum Vacc Immunother. 2013;9(5):1163–71.
Ishikawa T, Yamanaka A, Konishi E. A review of successful flavivirus vaccines and the problems with those flaviviruses for which vaccines are not yet available. Vaccine. 2014;32(12):1326–37.
Petersen LR, Powers AM. Chikungunya: epidemiology. F1000Res. 2016;5:82.
Wang JL, Li TT, Huang SY, Cong W, Zhu XQ. Major parasitic diseases of poverty in mainland China: perspectives for better control. Infect Dis Poverty. 2016;5(1):67.
Caraballo H, King K. Emergency department management of mosquito-borne illness: malaria, dengue, and West Nile virus. Emerg Med Pract. 2014;16(5):1–23. quiz 23–24.
Mackenzie JS, Williams DT. The zoonotic flaviviruses of southern, south-eastern and eastern Asia, and Australasia: the potential for emergent viruses. Zoonoses Public Hlth. 2009;56(6–7):338–56.
Zhong YB, Liu XQ, Deng YC, Xu PH, Zhong GR, Zhang W. First case of laboratory-confirmed Zika virus infection imported into China. Chinese Med J-peking. 2016;129(16):2013–4.
Tappe D, Rissland J, Gabriel M, Emmerich P, Gunther S, Held G, et al. First case of laboratory-confirmed Zika virus infection imported into Europe, November 2013. Euro Surveill. 2014;19(4).
We thank the Institute of Microbiology Epidemiology of Academy of Military Medical Sciences, National Institute for Viral Disease Control and Prevention and the 302nd Hospital of the Chinese People’s Liberation Army for providing cultivated samples of MBVs, non-MBVs and Uninfected mosquito pools.
This work was supported by the National Science and Technology Major Project (No. 2013ZX10004-101 and No. 2015ZX09J15105-001-001), the Science and Technology Major of Guangdong Province (No. 2012A080203005) and the Science and Technology Program of Tianjin (No. 15PTGCSY00150). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Availability of data and materials
The datasets supporting the conclusions of this article are included within Additional file 1.
YZ and QL wrote the main manuscript text; BZ participated in the experiment; XW participated in the experimental results discussed; SC and SW guided and developed the experimental program. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Figure S1. Determination of the specificity of the MBVs detection by detection of non-MBVs. A panel of non-MBVs were detected by this MBVs detection strategy to determine the specificity. There were no positive results (i.e. detection of these non-MBVs). YFV and EEEV were used as positive references. (XML 7 kb)
About this article
Cite this article
Zhang, Y., Liu, Q., Zhou, B. et al. Ultra-sensitive chemiluminescence imaging DNA hybridization method in the detection of mosquito-borne viruses and parasites. Parasites Vectors 10, 44 (2017). https://doi.org/10.1186/s13071-017-1975-1