Ultra-sensitive chemiluminescence imaging DNA hybridization method in the detection of mosquito-borne viruses and parasites

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.


Background
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 [1]. The globally pandemic MBV species encompass dengue virus (DENV), Japanese encephalitis virus (JEV), yellow fever virus (YFV), West Nile virus (WNV) [2] and chikungunya virus (CHIKV) [3]. 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 [4]. Factors like global warming, increased population mobility, and urbanization influence the epidemic scope of MBVs and sometimes even result in sudden outbreak [5]. In 1999, the outbreak of WNV in New York resulted in mass-scale animal mortality [6]. 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 [7]. In 2010, the outbreak of CHIKV in Guangdong province of China led to 173 cases of human infections [8]. From 2005 to 2006, more than 1.3 million cases of CHIKV infections were reported in India [9]. The infected cases of DENV recorded by the World Health Organization (WHO) were approximately 0.9 million [10]. 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 [11]. Each year, more than 20,000 casualties are reported due to severe infections of DENV [12].
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 [17]. 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 [18]. Mosquitoborne 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 [19].
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 [20]. Efficient prevention, control, and treatment of mosquitoborne 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][22][23][24][25], isothermal amplification [26], oligonucleotide microarray [27] and sequencing [28]. 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.

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 aminolabeled 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 (Ther-moFisher), 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) [29], 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-1thio-β-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

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 [30]. 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 streptavidinhorseradish 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 H 2 O 2 (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 selfprogramming 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 probessignal 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 cutoff 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

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   A repeat sequence of (T) 20 with an amino-labeled 3′-end and biotin-labeled 5′-end were used to calibrate the CL signal values during DNA hybridization f A human-original sequence was used as negative control 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. 2 Amplification of serially diluted VLPs by both RT-PCR and PCR. a The layout of the MBVs hybridization capture-chip. The probe P is the QC probe. The probe NC is the negative control probe. b Amplifications of serially diluted VLPs by both RT-PCR and PCR (both based on the GSPs of RT-PCR system A or B). The results showed that although DNase I did not completely digest DNA templates in high concentration VLPs, the extraction of low-concentration VLPs were only detected by RT-PCR and not by PCR, demonstrating that DNA templates had been diluted to sufficiently. c Principle of the CL imaging DNA hybridization method.
Step 1 shows that capture probes are fixed to the aldehyde-chip surface.
Step 2 shows that the denatured RT-PCR products are hybridized on the capture-chip. Steps 3-5 show the principle of CL detection. Biotin is incorporated into reverse strand in RT-PCR amplification. When HRP modified streptavidin is bound, CL signal becomes illuminant by the catalyzed substrates The lengths of in vitro transcribed RNAs were calculated on the length of the inserted gene fragments plus the transcript length of digested PGM-T vector c Copy numbers (copies/μl) = RNA concentration (ng/μl) × Avogadro constant × 10 -9 / (RNA length × 340 g/mol) d GenBank accession numbers (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 10 2 copies/μl except for WNV (10 3 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 10 1 copies/μl; for the others, the values are 10 2 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 10 1.8 PFU/ml and 10 0.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 10 6 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.

Discussion
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 [20]. Many genebased 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 Cut-off value of a probe (grey value) = average signal values of to detect other pathogens (except the particular pathogen detected by this probe) + 3 × standard deviation. A probe signal was considered positive when the value was greater than the cut-off value of a given probe 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 [33]. At present, human vaccines for flavivirus and alphavirus are available; however, these only include prevention of YFV, JEV, and TBEV [34][35][36][37][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-H 2 O 2 , 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][41][42][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 selfpowered 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 [30]. 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 10 1 -10 3 copies/μ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 10 0.8 -10 1.8 PFU/ml, and fits well in a linear regression model, demonstrating the feasibility of quantitative measurement of viral load. The assay is rapid, portable, ultrasensitive, 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 10 8 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