1. Introduction
Arboviral diseases are potent global public health problems due to their significant negative health and socio-economic impacts on many countries, particularly those in resource-limited settings. The major arboviruses belong to the families
Bunyaviridae, Togaviridae, and
Flaviviridae [
1] and are mainly transmitted to humans by hematophagous arthropods (such as mosquitoes, ticks, sandflies), from an animal reservoir or from an infected individual. Amongst the mosquito, the
Aedes and
Culex mosquitoes are considered the most medically important mosquito vectors with most of the arbovirues being transmitted by the Aedes mosquitoes, namely,
Aedes aegypti and Aedes albopictus [
2,
3,
4]. The mosquito-borne arboviruses, such as, yellow fever virus (YFV), Dengue virus (DENV), Zika virus (ZIKV), West Nile virus (WNV), and chikungunya virus (CHIKV), are globally distributed [
5,
6]; growing human populations, urbanization, climatic changes, international traveling, and trades contribute to the extensive spread of these viruses to areas where
Aedes mosquitoes have been infested [
2]. High mutation rates among RNA viruses such as these arboviruses provide conditions for adaptive evolution to new mosquito species and often gain a high degree of receptivity and infectivity which may facilitate arbovirus disease emergence [
7]. DENV and CHIKV are arboviruses of the family
Flaviviridae and
Togaviridae respectively. DENV is the most prevalent arbovirus, present in more than 100 countries within tropical and subtropical regions of the world [
8]. The virus is an enveloped single-stranded, positive-sense RNA virus, comprising four antigenically distinct serotypes, DENV-1 to DENV-4 that exhibit 65% to 70% sequence homology [
2,
3]. CHIKV has now been identified in over 110 countries in Asia, Africa, Europe and the Americas [
9] and is also a positive-sense single-stranded RNA virus with a genome of approximately 11.8 kb and comprises of two open reading frames (ORFs), ORF1 and ORF2, flanked by a 5′ cap and a polyadenylated tail at the 3′UTR [
10].
Moreover, when the multiple serotypes of DENV circulate concurrently with CHIKV, there is a higher risk for more severe forms of the disease such as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS) [
11]. This highlights the importance of accurate and early diagnosis in humans, determination of the circulating viral serotypes, and entomological surveillance at any given location [
6].
Between 2014–2017, 28.4% of dengue cases were serologically confirmed in India [
12]. In the year 2016, followed by 2017, and 2019 maximum number of chikungunya laboratory- confirmed cases were reported in India. Highest confirmed cases were reported in Karnataka, Delhi, and Maharashtra [
13]. In Burkina Faso, DENV is the most widespread arbovirus with the highest prevalence, incidence, and significant morbidities and mortality. In 2016, more than 1061 cases of dengue were reported in Burkina Faso with a case fatality rate of 1.2% [
14]. In 2017, a prevalence of 28.54% of dengue fever was reported among pregnant women in Ouagadougou [
15]. Another sero-epidemiological study of CHIKV was conducted in Ouagadougou, Burkina Faso, with blood samples collected in 2015, and a seroprevalence of 29.1% was reported [
16]. These consecutive reports in different populations (general and pregnant populations) suggest a considerable dengue burden in Burkina Faso. The risk of arbovirus transmission in the country is real because of its tropical climate, which favors the multiplication of
Aedes. Almost 80% of the Burkinabe population engages in farming and/or other agricultural practices. This may constitute a major factor of arbovirus transmission because, during the rainy season, farmers settle in the fields close to the forest and are in closer contact with vectors, which increases the risk of arbovirus transmission [
17,
18].
In the absence of a vaccine and specific antiviral treatment against the DENV and CHIKV,
the most effective means of disease control remain surveillance for early detection of cases in order to intervene with public health measures to contain the cases and thus control the infection.
Serological testing for these infections has well-documented limitations: antibodies immunoglobulin M (IgM) may not be detectable early in the course of infection, a rise inimmunoglobulin G (IgG) between acute and convalescent samples can only provide a retrospective diagnosis, and anti-flavivirus antibodies may cross-react with one another [
19]. IgM and IgG ELISA tests are widely used for rapid serological diagnosis but have the limitation of the inability to identify the circulating viral serotypes [
20]. Virus isolation and amplification in susceptible cell lines is the gold standard for the detection of serotypes, but it is not an appropriate clinical diagnostic assay in early infection since it is laborious and time consuming [
21]. For the prevention and detection of local transmission of DENV and CHIKV via both human and mosquito surveillance, developing a highly specific, sensitive and less expensive detection method with rapid outcomes appear particularly important. Multiplexed PCR-based assays often failed because multiple primers presented in high concentrations interact with each other unless they are exquisitely designed. Non-specific interference of oligonucleotides (DNA and RNA) is also thought to limit further multiplexed PCR [
22]. Many molecular methods have been developed to improve the simultaneous detection of several arboviruses around the world [
11,
23,
24,
25,
26]. The present study was conducted to develop a one-step multiplex
reverse transcription polymerase chain reaction (mRT-PCR) and
Real-time polymerase chain reaction (mRT-qPCR) for the rapid and simultaneous detection and serotyping of DENV and CHIKV to be used in arbovirus surveillance in Burkina Faso.
2. Materials and Methods
2.1. Viruses
Isolates of DENV and CHIKV were used for the development of multiplex PCRs and their analytical performances determination. Reference strains of DENV1-4 (1 isolate per serotype) were obtained from ATCC (VR-1856
™ DENV-1; VR-1584 ™ DENV-2; VR-1256_FD ™ DENV-3; VR-1490 ™ DENV-4) and CHIKV was isolated from a PCR positive human sera sample [
27]. The viruses were propagated in C6/36
Aedes albopictus cells prior to viral RNA extraction. Briefly, 20 µl of isolated sera sample was mixed with DMEM with 2% FBS, filtered and inoculated in confluent 24 well plates of C6/36 cells followed by incubation for 6 days. 90-100% confluent 6-well plates of Vero cells were used for CHIKV and DENV amplifications. 50 µL of virus stocks prepared was added to the flask and homogenized. The flask was incubated at 37°C for 1h, the volume of the flask was made to 25 mL with DMEM 2% FBS supplemented with 1% Pen/Strep and incubated until cytopathic effect was 90%. Infected cell supernatants were harvested in a 15 mL Falcon tubes and centrifuged at 2500 rpm for 5 min and stored at -80 °C. All procedures were carried out using sterile techniques and in a biosafety cabinet.
2.2. Patient Sera Samples
A total of 130 human serum samples were used in the present study after obtaining informed consent from the patients in both India and Burkina Faso. Among the 130 serum samples, 32 were positive for DENV (n=16) and CHIKV (n=16) and 32 from uninfected individuals (negative for both viruses) were part of a study funded by government of India [
28]. Also, 33 of the 130 serum samples were positive only for DENV and 33 serum samples were negative for DENV kindly offered by the National Institute of Public Health, Burkina Faso. The sera samples were well characterized by ELISA and singleplex RT-PCR methods prior to use for the validation (after the proof-of-concept step) of the new multiplex RT-PCR and RT-qPCR methods. The protocol of this study was reviewed and approved by ethical committees of the institutional ethical committees of the Health Science Research, Burkina Faso N° A026-2023/CEIRES/IRSS and ICGEB, New Delhi, ICGEB/IEC/2014/01 version 3.
2.3. Primer Design
Nucleotide sequences for the complete genome of each serotype of the Dengue virus were downloaded from the National Center for Biotechnology Information (NCBI) database and aligned to identify highly conserved regions using MEGA Software version 11. DENV1-4 primers were designed manually, forward conserved primer was designed for all serotypes, and reverse primers of specific serotypes were designed using only specific conserved regions of different serotypes of DENV. A new DENV-3 forward (DENV-3 q F) and the reverse DENV-3 R) have been designed by multi-alignment of another genome of four DENV serotypes for RT-qPCR use. The representative sequences of CHIKV were downloaded from Virus Pathogen Database and Analysis Resource (ViPR) and primers were designed to target the E1 gene using Snap Gene version 5.3.1. All RT-PCR primers were designed using similar parameters so that they would have similar melting temperatures (Tm). The potentials for dimerization and secondary structures were analyzed using Oligo Evaluator
TM online software (data not shown). The length of all amplicons was (200 bp DENV-1, 367 bp DENV-2, 1359 bp DENV-3, 118 bp DENV-4, and 574 bp CHIKV). The specificity of all primer sequences was further confirmed using primer BLAST (NCBI). Detailed information on primers is provided in
Table 1 and genome nucleotide accession numbers are in provided in Supplementary information.
2.4. Viral RNA Extraction
Viral RNA was extracted from 150 µL cell culture supernatant and clinical samples using a NucleoSpin RNA Virus kit (MACHEREY-NAGEL Gmbh Co.KG. Germany), according to the manufacturer’s instructions and eluted in 50 µL of diethyl pyrocarbonate (DEPC) treated water. A purified RNA was quantified using Nanodrop 2000 and cryopreserved at -80 °C until further processing.
2.5. One-Step Multiplex RT-PCR Amplification
Multiplex RT-PCR assay was optimized for simultaneous detection and serotyping of DENV and CHIKV RNAs in cell culture supernatants and human sera. mRT-PCR was performed using PrimeScript One Step RT-PCR Kit (Takara bio-INC). A total volume of 15µL of reaction mixture consisting of 100 ng of extracted RNA, 7.5 µL of 2X RT buffer, 0.6 µL of the enzyme, and 0.45µL (10µM) of forward conserved primer (DENV-1, DENV-2, DENV-3, DENV-4), 0.11µL (10 µM) of forward primer CHIKV, 0.3µL (10 µM) of Reverse primer DENV-1 and DENV-4, 0.11µL (10 µM) of Reverse primer of DENV-2, DENV-3, and CHIKV, the reaction was complete to 15µL with DEPC water. DEPC water was used for negative control. The thermal cycling profile of this assay consists of a 30 min Reverse Transcriptase (RT) step which is performed at 50°C, 95°C for 5 min, and then 35 cycles of 95°C for 30 s, 62°C for 20 s, 72°C for 30s and final extension 72°C for 5 min. mRT-PCR was performed using Applied Biosystems ProFlex™ 3 x 32-well PCR System machine. The PCR products were then analyzed by gel electrophoresis. 15 µL of PCR product were loaded into a 1% (W/V) agarose gel in 1x Tris-Acetate-EDTA buffer with a 1Kbp ladder as molecular weight marker. To detect coinfections between DENV 1-4 and CHIKV, individual serotypes RNA of reference sample were mixed in equal quantities and 3µL was used in mRT-PCR detection.
2.6. Multiplex One-Step Real-Time RT-PCR
A real-time, one-step, multiplex SYBR Green I RT-PCR assay was also developed for the detection of DENV-1 to DENV-4 and CHIKV. This assay was performed in PikoReal 96 Real-Time PCR System machine (Thermo Scientific, United States) using a one-step QuantiTect SYBR Green kit (Qiagen, Hilden, Germany). All mRT-qPCR reactions were performed in 50 µL reactions with 25 µL of SYBR Green, 0.5 µL of the enzyme, 1.4 µL (10µM) of forward conserved primer (DENV-1, DENV-2, DENV-4), 0.46 µl (10 µM) of DENV-3 q forward and reverse primer, 0.47µL (10 µM) of CHIKV forward and reverse primer, 0.5 µL (10 µM) of reverse primer DENV-1 and DENV-4, 0.44 µL (10 µM) of DENV-2 reverse primer, and 100 ng of RNA template. 15µL of reaction was used in triplicate for all viruses. DEPC water was used for negative control. The RT-PCR conditions for the real-time RT-PCR consisted of a 30-minute RT step at 50°C and 10 min of Taq polymerase activation at 95°C, followed by 40 cycles of PCR at 95°C for 30s (denaturation), 62°C for 30s (annealing), 72°C for 30s (extension) and final extension at 60°C for 30s. The melting curve temperature ranged from 60°C to 95°C. The result was positive if cycle threshold (Ct) values is equal to or less than 33 cycles. If Ct is more than 33 cycles, the result was considered negative. The limit of sensitivity of the assay was carried out with known quantitative RNA standards prepared using the one-step QuantiTect SYBR Green method. Briefly, the RNA template of each virus was serially diluted ten-fold with the known concentration of DENV 1-4, CHIKV, and 4µL of the diluted RNA were added to the mRT-qPCR reaction tube and amplified in triplicate. Standard curves have been determined using dilution 102 to 100 Copie/µL, for detection limits determination 10-1 and 10-2 have been added. Melting-curve analysis was performed after PCR amplification to verify that the correct product was amplified by examining its specific Melting temperature (Tm) that was also used to serotype DENV.
After optimization, the positive RNA of DENV and CHIKV and negative RNA from pre-collected sera were assayed in the same condition to evaluate the multiplex RT-PCR and RT-qPCR for their diagnostic potential in a clinical sample and validation. In Burkina, the clinical sample was performed using Applied Biosystems (SimpliAmp Thermo Fisher Scientific) thermal cycler for mRT-PCR and CFX96 Real-Time System (BIO-RAD) for mRT-qPCR.
2.7. Sequence Analysis
After gel electrophoresis, the amplicon of the different viruses was cut and purified using NucleoSpin® Gel and PCR Clean-up Kit (MACHEREY-NAGEL Gmbh Co.KG. Germany) according to the manufacturer’s instructions and eluted in 25 µl NE buffer. 20 𝜇L of purified PCR product was sent to Macrogen (Seoul, South Korea) for Deoxyribonucleic acid (DNA) sequencing by the Sanger method, using the amplification primers. Forward and reverse sequence was assembled by BioEdit Software version 7.2.5.0, to determine consensus sequence and blast using the NCBI BLAST tool to confirm the specificity and accuracy of the amplification of different viruses RNA fragments.
Author Contributions
Conceptualization, L.R.W.B., M.K.G., I.S., and S.S.; methodology, L.R.W.B., S.A.I., D.K.V. and A.K.; validation, L.R.W.B., S.A.I. and S.S.; writing-Original Draft, L.R.W.B., S.A.I., M.K.G., I.S., and S.S; supervision, S.A.I. and S.S. All authors have read and agreed to the published version of the manuscript.
Figure 1.
The DENV genome and primers targeted regions used in this multiplex detection study. Structural proteins are targeted in this study for DENV: DENV conserve forward is located in the capsid (C) protein, Reverse DENV-1 and DENV-4 in the C protein, and Reverse DENV-2 in precursor membrane prM protein. DENV-3 q forward and Reverse DENV-3 in envelope protein E; UTR: untranslated region.
Figure 1.
The DENV genome and primers targeted regions used in this multiplex detection study. Structural proteins are targeted in this study for DENV: DENV conserve forward is located in the capsid (C) protein, Reverse DENV-1 and DENV-4 in the C protein, and Reverse DENV-2 in precursor membrane prM protein. DENV-3 q forward and Reverse DENV-3 in envelope protein E; UTR: untranslated region.
Figure 2.
Comparison sensibility and specificity of singleplex and multiplex RT-PCR. Electrophoresis of 1% agarose gel loaded with 15µL of singleplex RT-PCR and multiplex RT-PCR product showing the specificity of the primers. (A) Singleplex RT-PCR: Lane 1, M (1kb DNA marker); Lane 2, DENV-1; Lane 3, DENV-2; Lane 4, DEN-3; Lane 5, DENV-4; Lane 6, CHIKV; Lane 7, Negative Control (NC). (B) Multiplex RT-PCR: Lane 1, M (1kb DNA marker); Lane 2, mDENV-1; Lane 3, mDENV-2; Lane 4, mDEN-3; Lane 5, mDENV-4; Lane 6, mCHIKV; Lane 7, NC. m: multiplex, bp: basis pair, NC: negative control.
Figure 2.
Comparison sensibility and specificity of singleplex and multiplex RT-PCR. Electrophoresis of 1% agarose gel loaded with 15µL of singleplex RT-PCR and multiplex RT-PCR product showing the specificity of the primers. (A) Singleplex RT-PCR: Lane 1, M (1kb DNA marker); Lane 2, DENV-1; Lane 3, DENV-2; Lane 4, DEN-3; Lane 5, DENV-4; Lane 6, CHIKV; Lane 7, Negative Control (NC). (B) Multiplex RT-PCR: Lane 1, M (1kb DNA marker); Lane 2, mDENV-1; Lane 3, mDENV-2; Lane 4, mDEN-3; Lane 5, mDENV-4; Lane 6, mCHIKV; Lane 7, NC. m: multiplex, bp: basis pair, NC: negative control.
Figure 3.
Coinfection and concurrent infection detection between DENV 1-4 and CHIKV. (A) concurrent infection using RNA mixture of DENV-1 + DENV-2, lane 2; DENV-1 + DENV-3, lane 3; DENV-1 + DENV-4, lane 4; M is 1kb DNA marker. (B) coinfection using RNA mixture of DENV-1 + CHIKV, lane 1; DENV-2 + DENV-3, lane 2; DENV-2 +DENV-4, lane 3; DENV-2 + CHIKV, lane 4; DENV-3 +DENV-4, lane 5; DENV-3 + CHIKV, lane 6; DENV-4 + CHIKV, lane 7; DENV-1 +DENV-3 + CHIKV, lane 8; M is 1kb DNA marker. Electrophoresis of 1% agarose gel loaded with 15µL product showing the specificity of the primers in coinfection detection with mRT-PCR.
Figure 3.
Coinfection and concurrent infection detection between DENV 1-4 and CHIKV. (A) concurrent infection using RNA mixture of DENV-1 + DENV-2, lane 2; DENV-1 + DENV-3, lane 3; DENV-1 + DENV-4, lane 4; M is 1kb DNA marker. (B) coinfection using RNA mixture of DENV-1 + CHIKV, lane 1; DENV-2 + DENV-3, lane 2; DENV-2 +DENV-4, lane 3; DENV-2 + CHIKV, lane 4; DENV-3 +DENV-4, lane 5; DENV-3 + CHIKV, lane 6; DENV-4 + CHIKV, lane 7; DENV-1 +DENV-3 + CHIKV, lane 8; M is 1kb DNA marker. Electrophoresis of 1% agarose gel loaded with 15µL product showing the specificity of the primers in coinfection detection with mRT-PCR.
Figure 4.
Amplification of dengue virus DENV 1-4 and CHIKV in mRT-qPCR. (A) Real-time RT-PCR amplification. (B) Melt curve analysis. Serial 10-fold dilutions were used in amplification by mRT-qPCR in triplicate. The results are positive if the Ct value is equal to or less than 33 cycles, otherwise negative if the Ct value is > 33 cycles. The melting curve temperature ranged from 60°C to 95°C. The color curve represents different virus dilutions in triplicate, The detection limit of this assay for DENV-1, DENV-2, DENV-3, DENV-4, and CHIKV are 10-2 copies/µL. DEPC water was used for negative control.
Figure 4.
Amplification of dengue virus DENV 1-4 and CHIKV in mRT-qPCR. (A) Real-time RT-PCR amplification. (B) Melt curve analysis. Serial 10-fold dilutions were used in amplification by mRT-qPCR in triplicate. The results are positive if the Ct value is equal to or less than 33 cycles, otherwise negative if the Ct value is > 33 cycles. The melting curve temperature ranged from 60°C to 95°C. The color curve represents different virus dilutions in triplicate, The detection limit of this assay for DENV-1, DENV-2, DENV-3, DENV-4, and CHIKV are 10-2 copies/µL. DEPC water was used for negative control.
Table 3.
Diagnostic performance of mRT-qPCR and mRT-PCR in clinical samples.
Table 3.
Diagnostic performance of mRT-qPCR and mRT-PCR in clinical samples.
Country |
Sample |
Total |
mRT-qPCR + N (%) |
mRT-qPCR - N (%) |
mRT-PCR + N (%) |
mRT-PCR - N (%) |
India |
DENV positive |
16 |
16 (100.0) |
0 (0.0) |
14 (87.5) |
2 (12.5) |
DENV negative |
16 |
0 (0.0) |
16 (100.0) |
0 (0.0) |
16 (100.0) |
CHIKV positive |
16 |
16 (100.0) |
0 (0.0) |
15 (93.7) |
1 (6.3) |
CHIKV negative |
16 |
0 (0.0) |
16 (100.0) |
0 (0.0) |
16 (100.0) |
Burkina Faso |
DENV positive |
33 |
33 (100.0) |
0 (0.0) |
30 (90.9) |
3 (9.1) |
DENV negative |
33 |
0 (0.0) |
33 (100.0) |
0 (0.0) |
33 (100.0) |
Sensitivity |
100% |
91 % |
Specificity |
100% |
100% |