Development of genus-specific universal primers for the detection of flaviviruses
Virology Journal volume 18, Article number: 187 (2021)
Flaviviruses are representative arboviruses carried by arthropods and/or vertebrates; these viruses can pose a public health concern in many countries. By contrast, it is known that a novel virus group called insect-specific flaviviruses (ISFs) also infects arthropods, although no such virus has yet been isolated from vertebrates. The characteristics of ISFs, which affect replication of human-pathogenic flaviviruses within co-infected mosquito cells or mosquitoes without affecting the mosquitoes themselves, mean that we should pay attention to both ISFs and human-pathogenic flaviviruses, despite the fact that ISFs appear not to be directly hazardous to human health. To assess the risk of diseases caused by flaviviruses, and to better understand their ecology, it is necessary to know the extent to which flaviviruses are harbored by arthropods.
We developed a novel universal primer for use in a PCR-based system to detect a broad range of flaviviruses. We then evaluated its performance. The utility of the novel primer pair was evaluated in a PCR assay using artificially synthesized oligonucleotides derived from a template viral genome sequence. The utility of the primer pair was also examined by reverse transcription PCR (RT-PCR) using cDNA templates prepared from virus-infected cells or crude supernatants prepared from virus-containing mosquito homogenates.
The novel primer pair amplified the flavivirus NS5 sequence (artificially synthesized) in all samples tested (six species of flavivirus that can cause infectious diseases in humans, and flaviviruses harbored by insects). In addition, the novel primer pair detected viral genomes in cDNA templates prepared from mosquito cells infected with live flavivirus under different infectious conditions. Finally, the viral genome was detected with high sensitivity in crude supernatants prepared from pooled mosquito homogenates.
This PCR system based on a novel primer pair makes it possible to detect arthropod-borne flaviviruses worldwide (the primer pair even detected viruses belonging to different genetic subgroups). As such, an assay based on this primer pair may help to improve public health and safety, as well as increase our understanding of flavivirus ecology.
Flaviviruses, which are representative arboviruses that infect arthropods and/or vertebrates, can cause public health problems worldwide . Generally, the genus flavivirus is divided into four groups: (i) mosquito-borne flaviviruses, (ii) tick-borne flaviviruses, (iii) insect-specific flaviviruses (ISFs), and (iv) no known vector (NKV) flaviviruses [2,3,4,5]. Of these, mosquito-borne and tick-borne flaviviruses, which include Dengue virus (DENV), West Nile virus (WNV), Yellow fever virus (YFV), Japanese encephalitis virus (JEV), Zika virus (ZIKV), and Tick-borne encephalitis virus (TBEV), are pathogens that play an important role in human diseases such as encephalitis, fever, and hemorrhagic fever [5, 6]. A novel virus group, ISFs was defined recently; this group also infects arthropods, although no ISF has yet been isolated from any vertebrate. ISFs comprise two subgroups: classical ISFs [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23] and dual-host affiliated ISFs [24,25,26,27,28,29,30,31,32]. Because some ISFs enhance or suppress replication of human-pathogenic flaviviruses within co-infected mosquito cells or mosquito laboratory colonies [19, 33,34,35,36], attention should be paid to both ISFs and human-pathogenic flaviviruses, even though ISFs appear not to be directly hazardous to human health. Therefore, to better assess the risk posed by flavivirus-induced diseases in humans, and to increase our understanding of the ecology of flaviviruses, we need to examine the distribution and diversity of flaviviruses using a broad detection system that can detect a wide range of these viruses in arthropods, which are known reservoirs.
Reverse transcription PCR (RT-PCR) is a well-established method used to detect/identify the genomes of micro-organisms such as bacteria, fungi, and viruses [37,38,39,40,41,42]. The RT-PCR method for detecting flaviviruses is well-established and used routinely by many laboratories [43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59]. Viruses belonging to the flavivirus genus have a positive-sense, single-stranded RNA genome (9.2–11 kbp), which encodes three structural (C, prM, and E) and seven non-structural (NS) proteins (NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5) . To date, different research groups have proposed sets of universal primer pairs designed to enable use of RT-PCR to detect NS proteins encoded by a range of flavivirus genomes [54,55,56,57,58,59]. However, although these primer pairs can efficiently detect virus species that are closely related, they are not very good at detecting more distantly related species or completely different species. In particular, these primer pairs may not detect more recently identified flavivirus strains. For example, classical ISFs which comprise a novel virus group, show high levels of genetic deviation from the three major subgroups, i.e., mosquito-borne flaviviruses, tick-borne flaviviruses, and NKV flaviviruses [2,3,4,5]. Thus, it is likely that the above-mentioned universal primers may not detect classical ISFs.
Therefore, the purpose of this study was to develop a detection system that recognizes flavivirus NS5 sequences broadly across the flavivirus genus. First, we identified conserved sites within the NS5 region of the flavivirus genome; such sites are conserved among many different subgroups of viruses, including mosquito -borne flaviviruses, tick-borne flaviviruses, ISFs, and NKV flaviviruses. We then designed a universal primer pair that amplifies this genome segment. To confirm the practical utility of this primer pair, we amplified the NS5 region from eight different flavivirus genomes (mosquito -borne flaviviruses, tick-borne flaviviruses, and ISFs). In addition, we used the primer pair to detect the flavivirus genome specifically in live virus-infected cells (a mosquito-derived cell line). Finally, we used the primer pair to amplify the flavivirus NS5 coding region from crude mosquito supernatants prepared from pooled mosquito homogenates.
Materials and methods
First, a universal primer pair (denoted as pan-flavivirus FW/RV; see the Table 1) was designed to target the NS5 coding region, which is conserved across different species of flavivirus. This was done by performing multiple alignments of NS5 region nucleotide sequences from mosquito-borne flaviviruses, tick-borne flaviviruses, ISFs, and NKV flaviviruses (Fig. 1 and Fig. 2). Sequence alignment was conducted using CLUSTALW and the results were visualized using MEGA  version 10.1.7. Although the predicted amplicon size generated by RT-PCR using this universal primer is similar among species, it remains species-dependent; the predicted amplicon sizes generated from the viral genome template used in this study are as follows: 1018 bp for DENV 2 PDK53 (GenBank accession KU725664); 1021 bp for ZIKV mosquito Haiti/1682/2016 (GenBank accession MF384325); 1027 bp for JEV Nakayama MY 2009 P578662 (GenBank accession HE861351); 1027 bp for WNV isolate BSL26-11 (GenBank accession JQ700442); 1027 bp for YFV strain ES-504 (GenBank accession KY885000); 1018 bp for TBEV strain C11-13 (GenBank accession KP644245); 1009 bp for Aedes flavivirus Narita-21 (GenBank accession AB488408); and 1009 bp for Culex flavivirus Tokyo (GenBank accession AB262759).
PCR amplification of flavivirus NS5 regions
To confirm the utility of the novel universal primer pair (denoted as pan-flavivirus FW/RV), it was used to amplify the NS5 region of representative virus strains belonging to the genus flavivirus (i.e., DENV, ZIKV, JEV, WNV, YFV, TBEV, Aedes flavivirus, and Culex flavivirus) in a PCR; other published universal primer pairs were tested for comparison [54,55,56,57,58]. Individual strains isolated in the 2000s (except JEV) were used to confirm that the universal primers can detect more recent viruses. The following oligonucleotides encoding the NS5 coding region (not the whole sequence) and the 3’-UTR were synthesized artificially (Eurofins Genomics, Tokyo, Japan) for use as templates in the PCR assay: DENV 2 PDK53 (corresponding to bp 8923–10602; accession KU725664); ZIKV mosquito Haiti/1682/2016 (corresponding to bp 8151–10808; accession MF384325); JEV Nakayama/MY/2009/P578662 (corresponding to bp 9036–10770; accession HE861351); WNV BSL26-11 (corresponding to bp 8164–11029; accession JQ700442); YFV ES-504 (corresponding to bp 8120–10699; accession KY885000); TBEV C11-13 (corresponding to bp 8149–10838; accession KP644245); Aedes flavivirus Narita-21 (corresponding to bp 7903–10760; accession AB488408); and Culex flavivirus Tokyo (corresponding to bp 7961–10837; accession AB262759). The PCR assay was performed using DNA polymerase (KOD plus Neo; Toyobo, Osaka, Japan) plus the above primers. Each 25 µl reaction mixture comprised 10 × PCR buffer (2.5 µl), 2 mM dNTPs (2.5 µl), 25 mM MgSO4 (1.5 µl), 10 µM forward primer (0.75 µl), 10 µM reverse primer (0.75 µl), template DNA (1 µl) or distilled water (1 µl; negative control), KOD plus Neo polymerase (0.5 µl), and distilled water (15.5 µl). PCR was conducted using a TaKaRa PCR Thermal Cycler Dice (Takara Bio, Shiga, Japan). The cycling conditions were as follows: initial denaturation at 94 °C for 2 min, followed by 35 amplification cycles (denaturation at 98 °C for 10 s, annealing at 60 °C for 30 s, and extension at 68 °C for 15 or 30 s), and a final extension at 68 °C for 2 min. The extension time was dependent on the primer pair used and on the size of the amplified product, i.e., 15 s for FU1/cFD2 and MA/cFD2, and 30 s for all other primer pairs (please see the Table 1). The amplicons produced by the PCR assay were confirmed by running the products on 1.5% agarose gels, followed by staining with ethidium bromide. The PCR product amplified by the novel universal primers (pan-flavivirus FW/RV) was sequenced.
Amplification of NS5 coding region from infected cells using the RT-PCR method plus the newly designed universal primer pair
Next, the new primer pair was used to detect the flavivirus genome in cells infected with JEV. Here, the NS5 coding region was amplified from cell lysate samples to confirm the specificity of the primer pair for the flavivirus genome. The mosquito cell line (C6/36), which is derived from Aedes albopictus, was cultured (4.0 × 105 cells/well) in 12-well plates at 28 °C in Eagle’s minimum essential medium (Wako Pure Chemical Industries, Osaka, Japan) containing 5% fetal bovine serum and standard antibiotics (penicillin (100 units/mL), streptomycin (100 µg/mL), and amphotericin B (250 ng/mL)). Cells were washed twice with PBS (supplemented with calcium/magnesium; PBS( +)) and then infected with a JEV live attenuated pig vaccine (m strain) (Kyoto Biken Laboratories, Uji, Kyoto) at a multiplicity of infection (m.o.i.) of 0.2, 0.02, and 0.002. Controls were mock-infected with PBS ( +). The m.o.i. was based on the cell number and the titer of the virus used (the virus titer was determined based on the TCID50, which was obtained using C6/36 cells and the Reed-Muench method ). Virus was resuspended in PBS ( +) and incubated with the cells for 2 h at 28 °C. After removing the virus solution, the cells were washed twice with PBS ( +) and cultured at 28 °C in Eagle’s minimum essential medium (Wako Pure Chemical Industries) containing 2% fetal bovine serum and standard antibiotics. Total RNA was extracted from infected cells using a NucleoSpin® RNA Plus (Takara Bio, Shiga, Japan) kit, and single-stranded cDNA was synthesized from the extracted RNA (0.2 µg) using random primers and ReverTra Ace (Toyobo), according to the manufacturer’s instructions. The NS5 coding regions were amplified using a pan-flavivirus FW/RV universal primer pair and DNA polymerase (KOD plus Neo, Toyobo), together with the resultant cDNA template from cells infected with the JEV porcine vaccine strain (m strain), distilled water (negative control), or the genome of Dengue virus 2 PDK53 (accession KU725664) (positive control), as described above, but with 40 amplification cycles. As an additional control for RNA extraction, the actin gene was amplified using the specific primer pair Act-2F (5′-ATGGTCGGYATGGGNCAGAAGGACTC-3′) and Act-8R (5′-GATTCCATACCCAGGAAGGADGG-3′) .
Detection of the viral genome from mosquito homogenate
Next, the primer pair was used to detect the flavivirus genome in crude supernatants prepared from pooled mosquito homogenates containing the JEV vaccine strain (m strain). Mosquitoes (Culex pipiens pallens) were obtained from Sumika Technoservice Corporation (Hyogo, Japan). Mosquito pools (one pool comprised ten mosquitoes) were placed into 50 ml safe-lock microfuge tubes together with zirconia beads and crushed under freezing conditions. Next, 1 ml of PBS (without calcium/magnesium; PBS(-)) was added to the tubes containing the crushed mosquitoes, which were then homogenized on ice. The above-mentioned JEV vaccine strain, which was titrated in C6/36 cells to obtain the TCID50 (as described above), was diluted tenfold and mixed with each homogenized mosquito pool (comprising 10 mosquitoes), followed by centrifugation at 9000 × g (4 °C for 10 min). The final concentration of JEV in the crude supernatant was 105.5–101.5 TCID50/ml. A QIAamp Viral RNA Mini Kit (Qiagen, Hilden, Germany) was used to extract total RNA from the centrifuged homogenate containing JEV (140 µl of whole homogenate was used for RNA extraction; the final volume of extracted RNA was 60 µl). Single-stranded cDNA was synthesized from 10 µl of the extracted RNA and used for cDNA synthesis using random primers and Rever Tra Ace (Toyobo), according to manufacturer’s instructions. The NS5 coding regions were amplified (40 amplification cycles) using the pan-flavivirus FW/RV universal primer pair and DNA polymerase (KOD plus Neo), together with the cDNA template derived from the mosquito cell homogenate.
Measurement of viral copy number by quantitative real-time PCR
In addition, we converted the virus titer in the spiked pooled mosquito samples into genomic copy number by performing quantitative real-time PCR (q-PCR) using cDNA synthesized from RNA extracted from a JEV vaccine strain (serially diluted tenfold from 105.5–100.5 TCID50/ml) and a standard plasmid encoding the JEV genome. Q-PCR was conducted using mix reagent (THUNDERBIRD®Probe qPCR Mix, Toyobo, Osaka) and the above-mentioned cDNA template derived from RNA extracted from the diluted JEV virus, or standard plasmids pEX-K4J2 (Eurofins Genomics) encoding the JEV NS5 region (Nakayama MY 2009 P578662), together with a TaqMan probe and primers (mFU1 and cFD2) targeting the NS5 gene . The TaqMan probe comprised an oligonucleotide with a fluorescent reporter dye attached to the 5’ end, and a non-fluorescent quencher (NFQ) attached to the 3’ end. The sequence of the mFU1 (forward) primer is 5′-TACAACATGATGGGAAAGCGAGAGAAAAA-3′, and that of the cFD2 (reverse) primer is 5′-GTGTCCCAGCCGGCGGTGTCATCAGC-3′. The oligonucleotide sequence of the TaqMan probe is FAM-TCCGTGACATAGCAGGAAAGCAAG-NFQ. Next, quantitative real-time PCR was performed using the CFX96 Real-Time PCR System (Bio-Rad, Hercules, CA) and the following cycling conditions: initial denaturation at 95 °C for 1 min, followed by 40 amplification cycles (denaturation at 95 °C for 15 s and annealing/extension at 58 °C for 45 s). The copy number of the viral genome in each cDNA sample was calculated by reading off the Ct values for cDNAs from a standard curve constructed by plotting the copy number of a plasmid (serially diluted tenfold) encoding the JEV genome versus Ct value (Additional file 1: Fig. S1). The copy number of the viral genome was determined in quadruplicate. The concentration of plasmid encoding the JEV NS5 region (Nakayama MY 2009 P578662) was determined in a NanoDrop ND2000 Spectrophotometer (NanoDrop Technologies, Inc., Wilmington, DE) and converted to molecular copies using the following formula:
Sequence alignment of NS5 regions and primer design
To design primer pairs that detect a broad range of flaviviruses, we performed multiple alignments of NS5 gene sequences from mosquito-borne, tick-borne, ISFs, and NKV flaviviruses. The multiple alignment identified an area of high homology between nucleotide sequences from individual strains of flaviviruses (Fig. 2), which enabled the design of a universal primer pair that detects flaviviruses broadly (pan-flavivirus FW/RV) (Table 1). To overcome some inevitable nucleotide variation (even in this highly homologous region among viral sequences), we gave mix bases such as “R”, “Y”, “H”, “B” or “N” within this primer pair. The nucleotide region bound by the pan-flavivirus RV primer overlaps that recognized by primers EMF1 and CFDJ9977 [54, 55, 58]; however, the lengths and nucleotide sequences of these primers are different from those of our novel primer.
Ability of the universal primers to detect individual flavivirus genomes
To confirm the utility of the primer pair that we designed (pan-flavivirus FW/RV), we used it to amplify the NS5 region of eight different flaviviruses species. We conducted PCR to amplify artificially synthesized oligonucleotides (Type II DENV, ZIKV, JEV, WNV, YFV, TBEV, Aedes flavivirus, and Culex flavivirus) using pan-flavivirus FW/RV and four other primer pairs (FU1/cFD2, MA/cFD2, EMF1/VD8, and FUDJ9166/CFDJ9977) [54,55,56,57,58], which are used commonly to amplify flavivirus genome. As shown in Fig. 3, the primer pairs that were used generated amplicons of the following sizes: FU1/cFD2 and MA/cFD2, < 300 bp; EMF1/VD8, < 700 bp; and FUDJ9166/CFDJ9977 < 900 bp (except for the amplicon of JEV), all of which corresponded to the sizes reported in the literature [54,55,56,57]. The amplicon obtained by PCR using the pan-flavivirus FW/RV primer pair was approximately 1 kbp in size, which corresponds to the expected sizes described in “Primer design” in the “Materials and Methods”. However, except for pan-flavivirus FW/RV, the other primer pairs did not necessarily amplify the NS5 region of all flavivirus species tested. Indeed, FU1/cFD2 amplified Type II DENV, ZIKV, JEV, WNV, YFV, and TBEV, but not ISFs. MA/cFD2 amplified ZIKV, JEV, YFV (weak band) and TBEV (although the product of ZIKV was not of the expected size), but no other species of flavivirus tested. EMF1/VD8 amplified Type II DENV, ZIKV, JEV, and YFV, but no other species tested. FUDJ9166/CFDJ9977 amplified Type II DENV, JEV, WNV, YFV, and ISFs (although the product of JEV was not of the expected size), but not ZIKV or TBEV. In addition, FUDJ9166/CFDJ9977 amplified some additional bands from the genomes of WNV and YFV. The unexpected band observed in the JEV sample is likely an extra-band; this band is larger than the expected size (841 bp) obtained from the template (JEV Nakayama MY 2009 P578662), suggesting that it could be the result of poor annealing to the target sequence. By contrast, pan-flavivirus FW/RV amplified all of the viral genomes and did not generate any non-specific bands. The intensity of the band for each amplicon produced by each primer pair was also different (Fig. 3). The intensity of the PCR products generated by FUDJ9166/CFDJ9977 and pan-flavivirus FW/RV was high; however, the intensity of those generated by the other primer pairs (i.e., FU1/cFD2, MA/cFD2, and EMF1/VD8) were not as high as that generated by pan-flavivirus FW/RV, indicating that these primer pairs are not the preferred choice for further analyses such as gene sequencing. To confirm whether the amplicons generated by the novel universal FW/RV primer are appropriate for gene sequencing, we sequenced the PCR products. The results showed that individual nucleotide sequences, excluding the PCR primers and unread regions at both ends of the PCR products, were 100% identical to the sequences of the templates (Type II DENV (accession KU725664), ZIKV (accession MF384325), JEV (accession HE861351), WNV (accession JQ700442), YFV (accession KY885000), TBEV (accession KP644245), Aedes flavivirus (accession AB488408), and Culex flavivirus (accession AB262759)).
Amplification of the NS5 region of JEV from infected cells
Next, we amplified the viral NS5 region from cells infected with live flavivirus. The mosquito cell line C6/36 was infected with the JEV vaccine strain at an m.o.i. of 0.2, 0.02, and 0.002, or mock-infected with PBS (−). Infected cells were collected, cDNA was synthesized from extracted RNA, and RT-PCR was performed using the universal primer pairs together with synthesized cDNA, distilled water (negative control), or the genome of DENV 2 PDK53 (accession KU725664) (positive control). Specific bands corresponding to the expected amplicon size of JEV were detected in the PCR assay using cDNA templates synthesized from cells at 24 h post-infection, although the intensity of each band was different; the intensity of the bands derived from cells infected at an m.o.i. of 0.2 was high, that from cells infected at an m.o.i. of 0.02 was modest, and that from cells infected at an m.o.i. of 0.002 was low (Fig. 4). These band intensities were higher when the PCR assay was performed using cDNA templates from the cells at 48 h post-infection, although again the intensities differed according to the m.o.i. No bands were observed under mock infection conditions (Fig. 4), suggesting that the primers were specific.
Amplification of the NS5 region of JEV from crude mosquito fluid
Finally, to confirm the utility of the novel primer pairs for amplification of the viral NS5 region from crude insect samples, we tested them using crude homogenates derived from pooled mosquito samples (10 mosquitoes per sample) spiked with different concentrations of the JEV vaccine strain (the virus titer in mosquito homogenates was 105.5–101.5 TCID50/ml). Briefly, cDNA was synthesized from RNA extracted from mosquito homogenates (including JEV), and RT-PCR was conducted using the primer pair; distilled water was used as a negative control (the RT-PCR conditions were the same as those used to amplify virus sequences from the infected mosquito cell line C6/36). Specific bands corresponding to the virus genome were detected in mosquito homogenate samples, with different samples showing different band intensities (Fig. 5). In the case of RNA extracted from the pooled mosquito sample, the detection limit was 102.5 TCID50/ml. We also converted the virus titer into viral genome copy number by performing qPCR using cDNA synthesized from RNA extracted from the JEV vaccine strain (serially diluted tenfold from 105.5–100.5 TCID50/ml) together with a standard plasmid encoding the JEV genome (Additional file 1: Fig. S1). The qPCR assay revealed that the detection limit in crude samples was 102.5 TCID50/ml (equivalent to 104.92 copies/ml of mosquito homogenate), which corresponds to 10–0.43 TCID50 (equivalent to 101.99 copies) for a single reaction in the PCR (Table 2).
Here, we report three main findings with respect to detection of flaviviruses using a novel universal primer pair. We show that: (i) the primer pair amplified the NS5 coding region of all flavivirus species tested; (ii) the primer pair detected the viral genome specifically in cells infected with live flavivirus; and (iii) the primer pair detected the viral genome in crude mosquito homogenates.
Of the primer pairs used in this study, the pair that we designed (pan-flavivirus FW/RV) was the only one able to amplify the NS5 region of all eight species of flavivirus tested; these viruses encompass all major subgroups of mosquito-borne flaviviruses, tick-borne flaviviruses, and ISFs. By contrast, amplification by other primer pairs (FU1/cFD2, MA/cFD2, EMF1/VD8, and FUDJ9166/CFDJ9977) was limited; in particular, they did not amplify the NS5 region of some recently identified flavivirus strains, including ISFs, even though the nucleotide sequences of EMF1 and CFDJ9977 overlap with a region that is also targeted by our pan-flavivirus RV. These results may be because the other primer pairs were designed in the 1990s; therefore, they are unlikely to detect recent strains (the flaviviruses used as oligonucleotide templates for the PCR assays conducted herein were isolated mainly in the 2000s). In addition, classical ISFs identified mainly in the 2000s [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23] (the exception is the ISF prototype “cell fusion agent virus”, which was discovered in 1975)  show large genetic differences from the major subgroups of mosquito-borne or tick-borne flaviviruses. These genetic differences mean that primer pairs FU1/cFD2, MA/cFD2, EMF1/VD8, and FUDJ9166/CFDJ9977, which were designed in the 1990s, are unlikely to detect classical ISFs because they do not recognize gene sequences recently identified as “insect-specific flaviviruses” [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32].
The PCR assay using the pan-flavivirus FW/RV primers yielded bands that were larger in size and higher in intensity than the other well-known primer pairs used in the same assay (Fig. 3 and the Table 1). The fact that the primer pair we designed produces a fragment large enough for nucleotide sequencing means that this primer pair is a potentially powerful tool for genetic analysis of both known and novel flavivirus strains. A recent review by Blitvich and Firth  suggests that some flaviviruses may be novel ISF strains; however, identification of some novel strains was based on ambiguous analyses using inadequate information obtained from genome sequences of limited size (< 300 bp), which are too short for reliable comparison. It is proposed that flaviviruses possessing > 84% nucleotide sequence identity should be classified within the same species . Therefore, in some cases short nucleotide sequences could misidentify “novel species” based on low (< 84%) nucleotide sequence similarity. Thus, to avoid erroneous classification, nucleotide sequences of appropriate length are needed for reliable analyses. The novel universal primer pair will allow reliable genetic analyses because it generates a product > 1000 bp in length, which is suitable for sequence analyses; this property increases the utility of this novel primer pair.
Global detection of flaviviruses in field samples requires a primer pair specific enough to detect the target sequences in the cell lysate. The PCR assay using cDNA templates synthesized from total RNA samples derived from infected mosquito cells revealed that the novel universal primer pair generated specific bands corresponding to the NS5 gene sequences of JEV. In addition, the intensity of the bands produced by pan-flavivirus FW/RV was dependent on the infection conditions (i.e., the m.o.i., the incubation time, and analysis time post-infection), suggesting that the intensity of individual bands reflects the amount of template (copy number) available. This result led us to surmise that our PCR system and specific primer pair can detect not only artificially synthesized oligonucleotides but also cDNA templates prepared from mosquito or tick cell lysates, which harbor variable flavivirus copy numbers. To prove this hypothesis, we used the new primer pair to detect the viral genome in crude mosquito homogenates. The detection limit in samples derived from mosquito homogenate was 102.5 TCID50/ml (equivalent to 104.92 copies/ml of mosquito homogenate or 101.99 copies/PCR reaction); these parameters are consistent with a recent report showing that the viral genome is preserved in wild mosquitoes. Martin et al.  examined the viral load of the cell fusion agent virus, which is a species of ISF, in individual mosquitoes trapped in the field. The viral load in males ranged from 1.25 × 102–5.50 × 106 RNA copies (median > 106 copies) per mosquito, whereas that in females ranged from 5.42 × 103–8.70 × 106 RNA copies (median, > 106 copies) per mosquito (although the viral load in gravid females was lower [3.29 × 102–3.53 × 106 RNA copies per mosquito; median 104 copies]). Based on the parameters reported by Martin et al., and the detection limit of JEV in our experiment, the novel universal primer pair is likely to detect virus in field-collected samples, although the detection limit may differ depending on the type of flavivirus present in the sample. One limitation of this study is that only JEV was used to determine the detection limit of the assay for the flavivirus genome. Further studies are required to determine the detection limit for viruses other than JEV, and to actually detect flaviviruses in field-collected samples.
The present results showing the specificity and detection threshold of the novel universal primer pair suggest that our novel primer pair could efficiently detect flavivirus at levels comparable with those identified in wild mosquitoes, making it a powerful tool for detecting flavivirus genomes in whole arthropod (mosquito and likely tick as well) homogenates.
Globalization and global warming mean that flavivirus-mediated diseases are a public health concern in many parts of the world. In addition, we are unsure about the potential of novel flaviviruses such as ISFs to evolve into major human pathogens. Also, some ISFs can affect the amounts of human-pathogenic virus harbored by arthropods; e.g., superinfection experiments suggest that some ISFs suppress or support replication of mosquito-borne flaviviruses [19, 33,34,35,36]; if this is the case, then ISFs may have an indirect effect on public health. Therefore, it is important to investigate the prevalence of flaviviruses, including novel strains or novel virus groups such as ISFs, in arthropods. This will allow risk assessments regarding possible flavivirus-mediated diseases. It will also increase our understanding of which species have positive or negative effects on replication of other viruses. In addition, unknown flavivirus species are likely to be discovered in hematophagous arthropods in the future. Some of these may even be the cause of diseases of unknown origin. Here, we developed genus-specific universal primers that may make it possible to amplify NS5 sequences from a broad spectrum of flaviviruses belonging to different genetic subgroups. This primer pair will allow us to better investigate the prevalence of flaviviruses in arthropods worldwide, improve public health and safety, and increase our understanding of flavivirus ecology.
Availability of data and materials
All data generated or analyzed during this study are included in the published article and its supplementary figure.
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.
Blitvich BJ, Firth AE. Insect-specific flaviviruses: a systematic review of their discovery, host range, mode of transmission, superinfection exclusion potential and genomic organization. Viruses. 2015;7(4):1927–59.
Gaunt MW, Sall AA, Lamballerie X, Falconar AKI, Dzhivanian TI, Gould EA. Phylogenetic relationships of flaviviruses correlate with their epidemiology, disease association and biogeography. J Gen Virol. 2001;82(Pt 8):1867–76.
Cook S, Holmes EC. A multigene analysis of the phylogenetic relationships among the flaviviruses (Family: Flaviviridae) and the evolution of vector transmission. Arch Virol. 2006;151(2):309–25.
Simmonds P, Becher P, Bukh J, Gould EA, Meyers G, Monath T, Muerhoff S, Pletnev A, Rico-Hesse R, Smith DB, et al. ICTV virus taxonomy profile: flaviviridae. J Gen Virol. 2017;98(1):2–3.
Pierson TC, Diamond MS. The continued threat of emerging flaviviruses. Nat Microbiol. 2020;5(6):796–812.
Roiz D, Vazquez A, Rosso F, Arnoldi D, Girardi M, Cuevas L, Perez-Pastrana E, Sanchez-Seco MP, Tenorio A, Rizzoli A. Detection of a new insect flavivirus and isolation of Aedes flavivirus in Northern Italy. Parasit Vectors. 2012;5:223.
Bolling BG, Eisen L, Moore CG, Blair CD. Insect-specific flaviviruses from Culex mosquitoes in Colorado, with evidence of vertical transmission. Am J Trop Med Hyg. 2011;85(1):169–77.
Bolling BG, Vasilakis N, Guzman H, Widen SG, Wood TG, Popov VL, Thangamani S, Tesh RB. Insect-specific viruses detected in laboratory mosquito colonies and their potential implications for experiments evaluating arbovirus vector competence. Am J Trop Med Hyg. 2015;92(2):422–8.
Calzolari M, Ze-Ze L, Ruzek D, Vazquez A, Jeffries C, Defilippo F, Osorio HC, Kilian P, Ruiz S, Fooks AR, et al. Detection of mosquito-only flaviviruses in Europe. J Gen Virol. 2012;93(Pt 6):1215–25.
Hoshino K, Isawa H, Tsuda Y, Sawabe K, Kobayashi M. Isolation and characterization of a new insect flavivirus from Aedes albopictus and Aedes flavopictus mosquitoes in Japan. Virology. 2009;391(1):119–29.
Hoshino K, Isawa H, Tsuda Y, Yano K, Sasaki T, Yuda M, Takasaki T, Kobayashi M, Sawabe K. Genetic characterization of a new insect flavivirus isolated from Culex pipiens mosquito in Japan. Virology. 2007;359(2):405–14.
Hoshino K, Takahashi-Nakaguchi A, Isawa H, Sasaki T, Higa Y, Kasai S, Tsuda Y, Sawabe K, Kobayashi M. Entomological surveillance for flaviviruses at migratory bird stopover sites in Hokkaido, Japan, and a new insect flavivirus detected in Aedes galloisi (Diptera: Culicidae). J Med Entomol. 2012;49(1):175–82.
Tyler S, Bolling BG, Blair CD, Brault AC, Pabbaraju K, Armijos MV, Clark DC, Calisher CH, Drebot MA. Distribution and phylogenetic comparisons of a novel mosquito flavivirus sequence present in Culex tarsalis Mosquitoes from western Canada with viruses isolated in California and Colorado. Am J Trop Med Hyg. 2011;85(1):162–8.
Cook S, Moureau G, Harbach RE, Mukwaya L, Goodger K, Ssenfuka F, Gould E, Holmes EC, de Lamballerie X. Isolation of a novel species of flavivirus and a new strain of Culex flavivirus (Flaviviridae) from a natural mosquito population in Uganda. J Gen Virol. 2009;90(Pt 11):2669–78.
Goenaga S, Fabbri CM, Garcia JB, Rondan JC, Gardenal N, Calderon GE, Enria DA, Levis SM. New strains of Culex flavivirus isolated in Argentina. J Med Entomol. 2014;51(4):900–6.
Vazquez A, Sanchez-Seco MP, Palacios G, Molero F, Reyes N, Ruiz S, Aranda C, Marques E, Escosa R, Moreno J, et al. Novel flaviviruses detected in different species of mosquitoes in Spain. Vector Borne Zoonotic Dis. 2012;12(3):223–9.
Crabtree MB, Sang RC, Stollar V, Dunster LM, Miller BR. Genetic and phenotypic characterization of the newly described insect flavivirus Kamiti River virus. Arch Virol. 2003;148(6):1095–118.
Hobson-Peters J, Yam AW, Lu JW, Setoh YX, May FJ, Kurucz N, Walsh S, Prow NA, Davis SS, Weir R et al: A new insect-specific flavivirus from northern Australia suppresses replication of West Nile virus and Murray Valley encephalitis virus in co-infected mosquito cells. PLoS One 2013, 8(2):e56534.
Huhtamo E, Moureau G, Cook S, Julkunen O, Putkuri N, Kurkela S, Uzcategui NY, Harbach RE, Gould EA, Vapalahti O, et al. Novel insect-specific flavivirus isolated from northern Europe. Virology. 2012;433(2):471–8.
Sang RC, Gichogo A, Gachoya J, Dunster MD, Ofula V, Hunt AR, Crabtree MB, Miller BR, Dunster LM. Isolation of a new flavivirus related to cell fusing agent virus (CFAV) from field-collected flood-water Aedes mosquitoes sampled from a dambo in central Kenya. Arch Virol. 2003;148(6):1085–93.
Feng Y, Li HB, Zhu J, Zhang YZ, Yang WH, Fan JH, Liang GD, Zhang HL. First isolation of quang binh-like virus from mosquitoes in China. Bing Du Xue Bao. 2014;30(1):57–61.
Zuo S, Zhao Q, Guo X, Zhou H, Cao W, Zhang J. Detection of Quang Binh virus from mosquitoes in China. Virus Res. 2014;180:31–8.
Huhtamo E, Cook S, Moureau G, Uzcategui NY, Sironen T, Kuivanen S, Putkuri N, Kurkela S, Harbach RE, Firth AE, et al. Novel flaviviruses from mosquitoes: mosquito-specific evolutionary lineages within the phylogenetic group of mosquito-borne flaviviruses. Virology. 2014;464–465:320–9.
Huhtamo E, Putkuri N, Kurkela S, Manni T, Vaheri A, Vapalahti O, Uzcategui NY. Characterization of a novel flavivirus from mosquitoes in northern europe that is related to mosquito-borne flaviviruses of the tropics. J Virol. 2009;83(18):9532–40.
Junglen S, Kopp A, Kurth A, Pauli G, Ellerbrok H, Leendertz FH. A new flavivirus and a new vector: characterization of a novel flavivirus isolated from uranotaenia mosquitoes from a tropical rain forest. J Virol. 2009;83(9):4462–8.
Kolodziejek J, Pachler K, Bin H, Mendelson E, Shulman L, Orshan L, Nowotny N. Barkedji virus, a novel mosquito-borne flavivirus identified in Culex perexiguus mosquitoes, Israel, 2011. J Gen Virol. 2013;94(Pt 11):2449–57.
Lee JS, Grubaugh ND, Kondig JP, Turell MJ, Kim HC, Klein TA, O’Guinn ML. Isolation and genomic characterization of Chaoyang virus strain ROK144 from Aedes vexans nipponii from the Republic of Korea. Virology. 2013;435(2):220–4.
Takhampunya R, Kim HC, Tippayachai B, Lee DK, Lee WJ, Chong ST, Kim MS, Lee JS, Klein TA. Distribution and mosquito hosts of Chaoyang virus, a newly reported flavivirus from the Republic of Korea, 2008–2011. J Med Entomol. 2014;51(2):464–74.
Evangelista J, Cruz C, Guevara C, Astete H, Carey C, Kochel TJ, Morrison AC, Williams M, Halsey ES, Forshey BM. Characterization of a novel flavivirus isolated from Culex (Melanoconion) ocossa mosquitoes from Iquitos, Peru. J Gen Virol. 2013;94(Pt 6):1266–72.
Pauvolid-Correa A, Solberg O, Couto-Lima D, Kenney J, Serra-Freire N, Brault A, Nogueira R, Langevin S, Komar N. Nhumirim virus, a novel flavivirus isolated from mosquitoes from the Pantanal, Brazil. Arch Virol. 2015;160(1):21–7.
Rizzo F, Cerutti F, Ballardini M, Mosca A, Vitale N, Radaelli MC, Desiato R, Prearo M, Pautasso A, Casalone C, et al. Molecular characterization of flaviviruses from field-collected mosquitoes in northwestern Italy, 2011–2012. Parasit Vectors. 2014;7:395.
Bolling BG, Olea-Popelka FJ, Eisen L, Moore CG, Blair CD. Transmission dynamics of an insect-specific flavivirus in a naturally infected Culex pipiens laboratory colony and effects of co-infection on vector competence for West Nile virus. Virology. 2012;427(2):90–7.
Kenney JL, Solberg OD, Langevin SA, Brault AC. Characterization of a novel insect-specific flavivirus from Brazil: potential for inhibition of infection of arthropod cells with medically important flaviviruses. J Gen Virol. 2014;95(Pt 12):2796–808.
Newman CM, Cerutti F, Anderson TK, Hamer GL, Walker ED, Kitron UD, Ruiz MO, Brawn JD, Goldberg TL. Culex flavivirus and West Nile virus mosquito coinfection and positive ecological association in Chicago, United States . Vector Borne Zoonotic Dis. 2011;11(8):1099–105.
Kuwata R, Isawa H, Hoshino K, Sasaki T, Kobayashi M, Maeda K, Sawabe K. Analysis of Mosquito-Borne Flavivirus Superinfection in Culex tritaeniorhynchus (Diptera: Culicidae) Cells Persistently Infected with Culex Flavivirus (Flaviviridae). J Med Entomol. 2015;52(2):222–9.
Rodriguez-Lazaro D, Cook N, Ruggeri FM, Sellwood J, Nasser A, Nascimento MS, D’Agostino M, Santos R, Saiz JC, Rzezutka A, et al. Virus hazards from food, water and other contaminated environments. FEMS Microbiol Rev. 2012;36(4):786–814.
Yamamoto Y. PCR in diagnosis of infection: detection of bacteria in cerebrospinal fluids. Clin Diagn Lab Immunol. 2002;9(3):508–14.
Petti CA. Detection and identification of microorganisms by gene amplification and sequencing. Clin Infect Dis. 2007;44(8):1108–14.
Siow WT, Koay ES, Lee CK, Lee HK, Ong V, Ngerng WJ, Lim HF, Tan A, Tang JW, Phua J. The use of polymerase chain reaction amplification for the detection of viruses and bacteria in severe community-acquired pneumonia. Respiration. 2016;92(5):286–94.
Ferrer C, Colom F, Frases S, Mulet E, Abad JL, Alio JL: Detection and identification of fungal pathogens by PCR and by ITS2 and 5.8S ribosomal DNA typing in ocular infections. J Clin Microbiol 2001, 39(8):2873–2879.
Khot PD, Fredricks DN. PCR-based diagnosis of human fungal infections. Expert Rev Anti Infect Ther. 2009;7(10):1201–21.
Lanciotti RS, Kerst AJ, Nasci RS, Godsey MS, Mitchell CJ, Savage HM, Komar N, Panella NA, Allen BC, Volpe KE et al: Rapid detection of west nile virus from human clinical specimens, field-collected mosquitoes, and avian samples by a TaqMan reverse transcriptase-PCR assay. J Clin Microbiol 2000, 38(11):4066–4071.
Pyke AT, Smith IL, van den Hurk AF, Northill JA, Chuan TF, Westacott AJ, Smith GA. Detection of Australasian Flavivirus encephalitic viruses using rapid fluorogenic TaqMan RT-PCR assays. J Virol Methods. 2004;117(2):161–7.
Sanchez-Seco MP, Rosario D, Domingo C, Hernandez L, Valdes K, Guzman MG, Tenorio A. Generic RT-nested-PCR for detection of flaviviruses using degenerated primers and internal control followed by sequencing for specific identification. J Virol Methods. 2005;126(1–2):101–9.
Tang Y, Anne Hapip C, Liu B, Fang CT. Highly sensitive TaqMan RT-PCR assay for detection and quantification of both lineages of West Nile virus RNA. J Clin Virol. 2006;36(3):177–82.
Ayers M, Adachi D, Johnson G, Andonova M, Drebot M, Tellier R. A single tube RT-PCR assay for the detection of mosquito-borne flaviviruses. J Virol Methods. 2006;135(2):235–9.
Kramer LD, Wolfe TM, Green EN, Chiles RE, Fallah H, Fang Y, Reisen WK. Detection of encephalitis viruses in mosquitoes (Diptera: Culicidae) and avian tissues. J Med Entomol. 2002;39(2):312–23.
Lanciotti RS. Molecular amplification assays for the detection of flaviviruses. Adv Virus Res. 2003;61:67–99.
Kauffman EB, Jones SA, Dupuis AP 2nd, Ngo KA, Bernard KA, Kramer LD. Virus detection protocols for west nile virus in vertebrate and mosquito specimens. J Clin Microbiol. 2003;41(8):3661–7.
Lampman RL, Krasavin NM, Szyska M, Novak RJ. A comparison of two West Nile virus detection assays (TaqMan reverse transcriptase polymerase chain reaction and VecTest antigen assay) during three consecutive outbreaks in northern Illinois. J Am Mosq Control Assoc. 2006;22(1):76–86.
Mendez F, Barreto M, Arias JF, Rengifo G, Munoz J, Burbano ME, Parra B. Human and mosquito infections by dengue viruses during and after epidemics in a dengue-endemic region of Colombia. Am J Trop Med Hyg. 2006;74(4):678–83.
Chao DY, Davis BS, Chang GJ. Development of multiplex real-time reverse transcriptase PCR assays for detecting eight medically important flaviviruses in mosquitoes. J Clin Microbiol. 2007;45(2):584–9.
Chang GJ, Trent DW, Vorndam AV, Vergne E, Kinney RM, Mitchell CJ. An integrated target sequence and signal amplification assay, reverse transcriptase-PCR-enzyme-linked immunosorbent assay, to detect and characterize flaviviruses. J Clin Microbiol. 1994;32(2):477–83.
Pierre V, Drouet MT, Deubel V. Identification of mosquito-borne flavivirus sequences using universal primers and reverse transcription/polymerase chain reaction. Res Virol. 1994;145(2):93–104.
Kuno G, Chang GJ, Tsuchiya KR, Karabatsos N, Cropp CB. Phylogeny of the genus Flavivirus. J Virol. 1998;72(1):73–83.
Kuno G. Universal diagnostic RT-PCR protocol for arboviruses. J Virol Methods. 1998;72(1):27–41.
Scaramozzino N, Crance JM, Jouan A, DeBriel DA, Stoll F, Garin D. Comparison of flavivirus universal primer pairs and development of a rapid, highly sensitive heminested reverse transcription-PCR assay for detection of flaviviruses targeted to a conserved region of the NS5 gene sequences. J Clin Microbiol. 2001;39(5):1922–7.
Maher-Sturgess SL, Forrester NL, Wayper PJ, Gould EA, Hall RA, Barnard RT, Gibbs MJ. Universal primers that amplify RNA from all three flavivirus subgroups. Virol J. 2008;5:16.
Kumar S, Nei M, Dudley J, Tamura K. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008;9(4):299–306.
Reed LJ, Muench H. A simple method of estimating fifty percent endpoints. Am J Hyg. 1938;27:493–7.
Staley M, Dorman KS, Bartholomay LC, Fernandez-Salas I, Farfan-Ale JA, Lorono-Pino MA, Garcia-Rejon JE, Ibarra-Juarez L, Blitvich BJ. Universal primers for the amplification and sequence analysis of actin-1 from diverse mosquito species. J Am Mosq Control Assoc. 2010;26(2):214–8.
Stollar V, Thomas VL. An agent in the Aedes aegypti cell line (Peleg) which causes fusion of Aedes albopictus cells. Virology. 1975;64(2):367–77.
Martin E, Tang W, Briggs C, Hopson H, Juarez JG, Garcia-Luna SM, de Valdez MW, Badillo-Vargas IE, Borucki MK, Frank M, et al. Cell fusing agent virus (Flavivirus) infection in Aedes aegypti in Texas: seasonality, comparison by trap type, and individual viral loads. Arch Virol. 2020;165(8):1769–76.
This work was supported by grants from the Oshimo Foundation and, in part, by the Heiwa Nakajima Foundation, the Takahashi Industrial and Economic Research Foundation, the Takeda Science Foundation, the JSPS Core-to-Core Program (B. Asia-Africa Science Platforms), and KAKENHI Grants from the JSPS (grant no. JP 20KK0199).
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Standard curve for the quantitative-PCR assay. A standard curve of Ct values versus copy number. The standard curve was constructed by 10-fold serial dilution (8.16 × 107 to 8.16 copies) of a plasmid encoding the JEV genome. All Ct values are the mean values of quadruplicate measurements, except for the value of 8.16 copies, which is the mean value of triplicate measurements (due to the detection limit).
About this article
Cite this article
Daidoji, T., Morales Vargas, R.E., Hagiwara, K. et al. Development of genus-specific universal primers for the detection of flaviviruses. Virol J 18, 187 (2021). https://doi.org/10.1186/s12985-021-01646-5
- Universal primer