key: cord-0753778-bftc5sfx
authors: Sarkar, Shovon Lal; Alam, A. S. M. Rubayet Ul; Das, Prosanto Kumar; Pramanik, Md. Hasan Ali; Al-Emran, Hassan M.; Jahid, Iqbal Kabir; Hossain, M. Anwar
title: Development and validation of cost-effective one-step multiplex RT-PCR assay for detecting the SARS-CoV-2 infection using SYBR Green melting curve analysis
date: 2022-04-20
journal: Sci Rep
DOI: 10.1038/s41598-022-10413-7
sha: e48b7583710ebc23003b4ea054ee27fe72ffea7d
doc_id: 753778
cord_uid: bftc5sfx

TaqMan probe-based commercial real-time (RT) PCR kits are expensive but most frequently used in COVID-19 diagnosis. The unprecedented scale of SARS-CoV-2 infections needs to meet the challenge of testing more persons at a reasonable cost. This study developed a simple and cost-effective alternative diagnostic method based on melting curve analysis of SYBR green multiplex assay targeting two virus-specific genes along with a host-specific internal control. A total of 180 randomly selected samples portioning into two subsets based on crude and high-quality RNA extraction were used to compare this assay with a nationwide available commercial kit (Sansure Biotech Inc., (Hunan, China)), so that we could analyze the variation and validity of this in-house developed method. Our customized-designed primers can specifically detect the viral RNA likewise Sansure. We separately optimized SYBR Green RT-PCR reaction of N, E, S, and RdRp genes based on singleplex melting curve analysis at the initial stage. After several rounds of optimization on multiplex assays of different primer combinations, the optimized method finally targeted N and E genes of the SARS-CoV-2 virus, together with the β-actin gene of the host as an internal control. Comparing with the Sansure commercial kit, our proposed assay provided up to 97% specificity and 93% sensitivity. The cost of each sample processing ranged between ~2 and ~6 USD depending on the purification level of extracted RNA template. Overall, this one-step and one-tube method can revolutionize the COVID-19 diagnosis in low-income countries.

In brief, the primer sets were carefully designed against specifically chosen viral genes. We initially standardized the RT-PCR condition for each gene and then targeted multiple genes for conducting multiplex PCR by combining the virus-specific primer sets along with a house-keeping gene of humans as the internal control. Optimization of the annealing temperature and concentrations of those primers and using melting curve analysis mainly sifted out the non-specificity issues of the SYBR green method. We then validated the assay with a total of 180 clinical samples where an equal number of samples were extracted either by crude or column-based RNA extraction system. To assess the effectiveness of our assay, the Sansure kit was considered as the gold standard here. The overall representation of the workflow is delineated in Fig. 1 .

Primer sequence validation. The oligo-analyzer tool showed no stable secondary structure, hairpins, homodimers, and cross-dimer formed within the primer set sequences. The result of the Primer-BLAST tool from the NCBI showed that all the genes' primer set only matched with the expected target size of the SARS-CoV-2 virus genome (reference strain Wuhan-Hu-1), indicating that the performance of the SYBR Green assay would be specific enough for SARS-CoV-2 detection. Searching against the human genome and other pathogens showed no similarity, bolstering that the viral genome segments alone would be amplified (Supplementary  Table S1 ). We found no mutation that engendered the new variants within the primer-binding region of the viral genome. Mutations with low frequency were identified within the last five bases of spike forward and reverse primers, while RdRp, N, and E targeted primers were void of end site mutations. The result of the searches were represented in the Supplementary Table S1 . Development of singleplex SYBR green assay. Using JUST_N1, JUST_E1, JUST_S1, JUST_RdRp1, β-actin, and GAPDH primer set individually, and singleplex assays were at first performed for four randomly selected clinical samples (three positive low-Ct and one negative, previously detected by Sansure RT-PCR kit). Specific desired band position for each gene target was observed in the electrophoresis gel as mentioned in the Table 1 (Table 1 and Supplementary Fig. S1 ). In the melting curve analysis for positive samples, the amplicon of JUST_N1, JUST_E1, JUST_S1, and JUST_RdRp1 primer set produced a specific melting temperature (Tm) peak at 82.32 ± 0.17 °C, 79.40 ± 0.31 °C, 76.52 ± 0.17 °C, and 77.57 ± 0.17 °C, respectively ( Fig. 2 and Supplementary Fig. S2 ). However, a distinct dissociation curve generated for the negative samples suggested a non-specific signal or primer-dimer formation (Supplementary Table S2 ). The amplicon of housekeeping gene β-actin and GAPDH produced a specific Tm peak at 85.78 ± 0.24 °C and 87.59 ± 0.18 °C (Fig. 2) and 82.95 ± 0.036 °C (Supplementary Fig. S2 ), respectively, for both positive and negative samples. Exact details of the information were shown in Supplementary Table S2 .

Optimization of multiplex SYBR green assay. Considering two positive clinical samples, we used JUST_N1, JUST_E1, JUST_S1, JUST_RdRp1, β-actin, and GAPDH primer sets to conduct duplex assay in nine different combinations: JUST_N1 + JUST_E1, JUST_N1 + JUST_S1, JUST_N1 + JUST_RdRp1, JUST_ E1 + JUST_S1, JUST_E1 + JUST_RdRp1, JUST_S1 + JUST_RdRp1, JUST_N1 + GAPDH, JUST_E1 + GAPDH, and JUST_S1 + GAPDH; triplex assay in five different combinations: JUST_N1 + JUST_E1 + JUST_S1, JUST_ N1 + JUST_RdRp1 + JUST_S1, JUST_N1 + JUST_RdRp1 + JUST_E1, JUST_N1 + JUST_E1 + β-actin, and JUST_N1 + JUST_E1 + GAPDH. The quadruplex assay was performed using primer sets combined with JUST_ N1 + JUST_E1 + JUST_S1 + JUST_RdRp1 for one of those positive samples (Supplementary Table S2 ).

In the duplex assays, the amplicon of JUST_N1 + JUST_E1; JUST_N1 + JUST_S1; JUST_N1 + JUST_RdRp1; JUST_E1 + JUST_RdRp1; JUST_N1 + GAPDH; JUST_E1 + GAPDH, and JUST_S1 + GAPDH produced specific www.nature.com/scientificreports/ In the triplex assays, the amplicon of JUST_N1 + JUST_E1 + JUST_S1; JUST_N1 + JUST_RdRp1 + JUST_S1, and JUST_N1 + JUST_RdRp1 + JUST_E1 produced two specific Tm peaks at (82.42 °C and 76.02 °C), (82.27 °C and 77.23 °C), and (82.27 °C and 79.07 °C), respectively. Notably, three specific amplicons were appropriately produced in the triplex assay of JUST_N1 + JUST_E1 + β-actin and generated Tm peak at 81.82 °C, 78.49 °C, and (85.78 °C and 87.73 °C -a signature Tm peak for β-actin), respectively. In contrast, the amplicon of the housekeeping gene GAPDH combined with JUST_N1 and JUST_E1 in the triplex assay was not amplified properly and produced two specific Tm peaks at 82.20 °C and 78.43 °C, respectively. Finally, in the quadruplex assay of JUST_E1_F  ATT CGT TTC GGA AGA GAC  AGG  63  This study  117  02-29  26,254-26,274   JUST_E1_R  CGC ACA CAA TCG AAG CGC  62.3  This study  106-123  26,350-26,367   JUST_N1_F  ACC CAA TAA TAC TGC GTC  TTGG  63.5  This study  138   135-156  28,409-28,430   JUST_N1_R  GGT AGC TCT TCG GTA GTA  GCC  64.4  This study  253-273  28,527-28,547   JUST_RdRp1_F  GTA CTG ATG TCG TAT ACA  GGGC  63  This study  104   79- www.nature.com/scientificreports/ JUST_N1 + JUST_E1 + JUST_S1 + JUST_RdRp1, only two amplicons were generated and produced two specific Tm peaks at 82.12 °C and 77.54 °C.

Since the amplicon size produced by using the primer set JUST_S1 is small (66 bp), it produced a Tm peak at a lower temperature. A similar Tm peak also appeared due to primer-dimer or non-specific amplified or fragmented products. We thus could not distinguish the S-gene-specific amplified product. There had also arisen some cases of mutations within the last 5 bases of the S targeting primer region (Supplementary Table S1 ) during the experiment. Hence, we eliminated the S targeting primers from the further multiplex study. We also omitted JUST-RdRp1 from multiplex assay for its small product size (104 bp) and less accurate result in the melting curve when run with other gene targeting and host-specific primer sets in the respective triplex assays. Housekeeping gene target GAPDH was omitted from the final multiplex assay since it produced a closer Tm peak with JUST_N1 and was not amplified efficiently with the JUST_E1 combination. Therefore, we selected JUST_N1 + JUST_E1 + β-actin primer combination for the final multiplex SYBR Green assay to detect nucleocapsid and envelope gene of SARS-CoV-2 and beta-actin gene of the host as the internal control. The result is summarized in Supplementary Table S2 and Supplementary Fig. S3 . The schematic representation of the optimized protocol is shown in Fig. 3 .

confirming the primer set for multiplex reaction, we intended to optimize the method for primer concentrations and PCR conditions. For initially establishing optimal SYBR green-based protocol, we used three positive samples and one negative sample as previously confirmed by Sansure. The singleplex qRT-PCR was performed using different sets of primer (JUST_N1, JUST_E1, JUST_S1 and JUST_RdRp1) at different concentrations (50-250 nM) and different annealing temperatures (60-70°C). In our study, the optimal condition of annealing temperature was set at 62 °C and 200 nM primer concentration. We also lessened the annealing and extension time (combinedly 25 s instead of 1 min as per standard protocol) to reduce the primer dimer formation and get a higher quantity of target amplicons.

In this study, the amplicon of JUST_N1 primer and JUST_E1 primer produced a specific Tm peak at 81.8 ± 0.40 °C and 78.4 ± 0.33 °C in the melt curve for positive samples. For internal control, β-actin primer amplified product generated specific Tm peaks in two positions at around 85.8 ± 0.3 °C and 87.8 ± 0.3 °C, that is a signature peak (Fig. 2) . However, a clearly distinct dissociation curve was generated for negative samples that varied with the melt curve as identified for the positive samples, suggesting non-specific signals or the formation of primer-dimers in the melt curve peaks for the negative samples. This is the crucial parameter in the analysis of the specificity of curves for the SYBR Green methodology. www.nature.com/scientificreports/ Assessment of SYBR green assay against Probe-based method (Sansure kit): crude RNA extraction method. We assayed 90 samples from COVID-19 cases with Sansure kit and our optimized SYBR Green-based one-step RT-PCR protocol. The RNA was extracted using crude method beforehand. In both methods, 49 and 8 samples were positive and negative, respectively ( Table 2 and Fig. 4) . We observed rest of the 33 samples with contradictory results (27 false-positive and 6 false-negative) as shown in Table 2 and in Appendix 1 in details. To reconfirm the presence of virus, we performed the targeted PCR amplification of spike protein coding sequence for these 33 conflicting samples and did sanger sequencing for three representative samples: one from concordant positive (164.02), false-positive (164.36), and false-negative (165.104) each group. Then, we identified sequence matching with the SARS-CoV-2 spike in BLAST (> 99%) and MEGA7 based alignment (Supplementary Fig. S5 and S6). Among 33 samples, we designated eight of them as undetermined where the clear amplified product could not be confirmed. We found a 47% (9/19) positive and 67% (4/6) negative results for spike amplified products corroborated with our SYBR-green based assay. On the other hand, a 33% (2/6) positive and 53% (10/19) negative results accorded with the Sansure kit method ( Table 2 , Supplementary Fig. S5 , and Appendix 1). Table 2 . Results of samples for crude and column-based RNA extraction methods. a Undertermined results were designated based on spike coding sequence targeted amplification where we could not identify a particular band as positive or negative for the gene. (Table 2 and Appendix 2). After spike gene specific amplification, we found three of them negative and one undetermined ( Supplementary Fig. S5 , Table 2 and Appendix 2). The triplicate runs of the multiplex assay showed consistent results for all the samples except four samples (Appendix 2). . We designated nine samples as undetermined through this process where the gel image could not represent exact results ( Table 2 ). The sanger results showed the similarity of the sequence to the SARS-CoV-2 ( Supplementary Fig. S6 ).

The values of the sensitivity, specificity, positive likelihood ratio, negative likelihood ratio, positive predictive value, negative predictive value, and kappa index of the SYBR Green assay for both crude and column-based RNA extraction approaches were presented in Table 3 . By column-based extraction, the sensitivity and the specificity of the SYBR green assay raised to 7.74% (96.83-89.09)% and 69.73% (92.59-22.86)%, respectively, when we considered Sansure solely as the gold standard. The kappa index of the crude vs column-based extraction also varied (0.134 vs 0.89). Considering the results of spike PCR amplification and further Sanger sequencing of representative samples for the false-positive and false-negative samples, the values increased (Table 3) . To be more precise, the difference in specificity increased 38.31% (92.86-54.55)%. The approximate cost of detection for the crudely extracted sample was < $2, but for the column-based detection, its cost was ≤ $6 (Supplementary Table S5 ).

In crude extraction methods, both N and E gene-specific peaks were found in 32 out of 76 concordant samples, and either N or E gene was detected in rest of the 44 samples (Supplementary Table S4 in 10 0 ,10 -1 , 10 -2 , 10 -3 , 10 -4 , and 10 -5 dilution of template RNA (Fig. 5) . The R 2 value of the N and E were 0.99 and 0.98, respectively, for gene-specific RT-PCR showing the curve fitting the data well. The slope value of the N gene based curve (-3.55) showed a good correlation between viral RNA numbers and Ct values (Fig. 5) . For E gene, we did not get a better slope value (-3.0), given that this gene is present in low copy number in virus-infected cell. Since the N or E gene RNA transcripts within the cell other than a part of viral genome do present within the samples, we cannot correlate our data with viral copy number directly. Table 3 . Operating characteristics of in-house SYBR green assay. PPV positive predictive value, NPV negative predictive value, LR + positive likelihood ratio, LR − negative likelihood ratio. a Considering Sansure for condodant samples and S gene amplification for contradictory false-positive and negative samples as the standard, respectively. The sample number here would be 82 (n = 82) due to undetermined results of eight samples in spike amplification step. b Considering Sansure as standard for 86 concordant samples and S gene amplification as standard for the four contradictory samples. The sample number would be 89 (n-89) here due to one undetermined result in PCR and gel run. 

COVID-19 will continue worldwide in the foreseeable future. Apart from regular diagnostic purposes in the health sector, regular monitoring of the disease is necessary for many other sectors through valid and recognized scientific testing (i.e., real-time PCR). Even after the waning of the pandemic in the future, the COVID-19 testing will continue due to safety issues, research purposes, and the possibility of emerging vaccine-escape mutants. Developed countries can test a vast number of samples every day to control virus spread. In contrast, developing countries are still facing the obstacle of diagnosing a considerable number of patients. The testing cost is also higher in developing countries because of no raw material manufacturing industries and cross-border custom tax issues. Nevertheless, laboratories and diagnostic centers prefer the high-price probe-based method to detect SARS-CoV-2. The main reason is the absence of an alternative, low-cost, and comparably effective reverse transcription real-time PCR (RT-PCR) assay. We thus developed and validated a multiplex SYBR Green method to mimic the nationwide regularly used Taqman-based commercial kit (Sansure Biotech Ltd.) assay, even at a low-priced budget.

Notably, most previous studies were performed in two-step processes (cDNA preparation and RT-PCR reaction in different tubes). We performed those in a one-step reverse transcription RT-PCR using a single tube (cDNA preparation and PCR in the same tube). For instance, Gomez et al. (2020) and Dorlass et al. (2020) converted extracted RNA into cDNA first and then performed the PCR with the cDNA template 7,12 . This one-step approach eases the workflow, expedites overall processing, and reduces the chance of contamination.

Four customized primer sets against N, E, RdRp, and S genes of SARS-CoV-2 were employed for two reasons. The N and S are the most abundantly expressed transcripts during viral replication 15 , and the E and RdRp are the least evolving genes due to minimized effects of the selective pressure 16 (Fig. 1) . We could detect the specific genes of only SARS-CoV-2 with these primer sets, provided that other human coronaviruses might present in the clinical samples. However, we found a partial matching of the JUST_RdRp1, JUST_N1, and JUST_E1 primers with the SARS coronavirus Tor2 genome, but SARS was not detected in humans and other animals after 2005 17 . Moreover, our designed primers targeted the SARS-CoV-2 solely through screening out the mutations present among the GISAID global viral sequences within the last five bases of the 3´-end of the primer. The primers were also screened against the co-variants and clade-featured mutations that independently originated in different parts of the world, and notably, those mutants emerged after most of the previous SYBR green-based methods 18 . The finally selected primers (JUST_N1 and JUST_E1) to detect SARS-CoV-2 in the optimized multiplex reaction were also been found effective against the recently circulating and emerged variants of concern (VoC), such as Delta and Omicron (Supplementary Table S1 ). Thus, this assay can effectively detect the variants. www.nature.com/scientificreports/ Another supremacy of our study is the multiplexing of viral N and E genes with host-specific β-actin gene. As an internal control, β-actin here was used for the same purpose as in TaqMan assays for confirming the human sample and whether the sampling has been collected accurately from the nasopharyngeal swab. This design engendered a unique approach and was not present in any existing SYBR Green-based study for COVID-19 diagnosis. However, the minor drawback of a multiplex SYBR green system is unable to predict the presence or absence of any particular gene or primer-dimer by specific Ct-value. Melting curve-based detection, wherein the dissociation points of three genes are distinct, and any flat curve on that melting point predicted the absence of the specific gene in the sample, herein has solved the issue. Hence, our study is superior to other studies in the imitation of standard TaqMan assay to separate the disparate COVID-19 positive and negative samples. Pereira-Gómez et al. (2021) used primer sets against ORF1b-nsp14 or N gene targets using melting curve analysis to differentiate primer-dimer from the target product amplification 12 . Park et al. (2020) established a multiplex cost-effective SYBR Green-based method using primer sets targeting RdRp, N, E, and S genes while suggesting melting curve analysis for confirming the specific target product amplification. However, their works did not include the internal control particular primers in the multiplex reaction, so those assays will not be appropriate for standard diagnostic purpose 10 .

The peak(s) in the melt curve shows the presence of a particular gene-targeting amplicon(s), and the derivative reporter is a measurement for the approximate quantity of that target gene(s). The melting curve analysis can also help identify mutations occurring within the target region of the amplicon 12 SYBR Green RT-PCR assays after crude and column-based RNA extraction exhibited an entirely different picture in terms of specificity and positive likelihood ratio (LR +). Even after considering spike amplification and Sanger, we found the specificity range up to ~ 75% at 95% CI (Supplementary Table S3 ). The reason is that these statistical parameters heavily rely upon the false-positive value. We found more false-positive results for our SYBR Green assay in crude extraction system-based assay. There could be two-fold reasons for this issue: false-negative results for those samples in Sansure kit and/or S-gene amplification and a non-specific product amplification in the SYBR Green assay. In a crude extraction system, the chance of getting positive results in Sansure kit-based PCR will be lower because the possibility of inhibiting the fluorescent probe in Sansure by extraction systembased inhibitors is higher. The lowest Ct value of those contradictory samples for both crude and column-based (≥ 35) was the proof (Appendices 1 and 2) . It is sometimes difficult to get concrete positive results in conventional PCR where RT-PCR failed. The nine undetermined results in 31 spike gene-specific conventional PCR indicated such an issue. The undetermined results in the gel run also matched with the RT-PCR invalid result where the guaranteed amplification could not be ensured. Because of these reasons, more positive results might be found in SYBR Green RT-PCR compared to Sansure and Spike gene-based amplification.

Although the fluorescence probe-based RT-PCR is advantageous for specific target gene detection, any interfering molecule for the fluorophores and inhibitors to polymerase's 5′-3′ exonuclease activity decreases its efficiency. The Sansure kit has such issues that can lead to underestimating viral RNA quantity or even falsenegative results 22 even with high-performance extracted RNA in some cases 23 The low amount of viral RNA and extraction-related fluorescence-inhibitors made it more challenging to identify the viruses in Sansure assay than SYBR Green or conventional PCR.

The RT-PCR has incredibly advanced molecular diagnostics; conventional PCR (cnPCR) can still show comparable efficiency 24 The cnPCR was found to detect SARS-CoV-2 with an adequate yield as an alternative to RT-PCR 25 . For the guaranteed presence of the virus, the cnPCR is used in combination with the Sanger sequencing. Spike-gene specific amplification following Sanger sequencing of representative cnPCR positive samples thus confirmed the presence of virus among the contradictory samples (false-positive and false-negative). Spike gene amplification was considered since it is the second-highest expressed gene and is also less mutating than the N gene 26 . A small amplicon (228 bp) was targeted to identify even the fragmented RNA genome due to sampling issues and shorter transcripts. Besides, the cnPCR worked well in some cases where the probe-based method failed due to fluorescence-blocking inhibitors 27 . We determined 17 instances for the Sansure negative samples showing positive or undetermined results in the cnPCR ( Table 2) .

Standard curve analysis showed a slightly less correlation and efficiency than desired for E gene (slope = − 3.0, R 2 = 0.98) in 10 0 -10 -5 dilution series. The reason might be the presence of non-specific RNA after extraction and inhibitors with extracted RNA, presence of E gene at low concentration, and primer-dimer formation. The values for the N gene fell within the standard range (slope = − 3.55, R 2 = 0.99) as it is present in a higher number than the E gene.

Our target is not to quantify instead to detect the presence or absence of virus in the clinical samples that we have accomplished here. Nevertheless, this study has some limitations. In melt curve analysis, the dominance of N and E gene presence was observed in column-based and crude RNA extraction, respectively. Other researchers, however, reported a similar trend despite working with a very small percentage of the samples 28 www.nature.com/scientificreports/ lower amount of RNA template, multiplexing of three genes in SYBR Green, and arises of possible mutations during the study in highly dynamic nucleocapsid protein 26,28 might reduce the chance to amplify N gene-specific region. Moreover, the sensitivity and the specificity were not wholly determined because of not validating all the studied samples further in cell culture. However, we endeavored to determine an identical sensitivity and specificity while comparing the deviated samples between SYBR Green and Sansure assay ( Table 2 ). There is a huge chance to recover lower concentrated RNA in column-based extraction after elution. Still, in the quick extract method, the recovery rate of lower concentrated RNA is low, giving false-negative results ( Table 2) 14 .

Another subtle limitation of our technology is that only expert personnel will analyze the results since it is based on melting curve analysis. Finally, PCR/Sanger-based validation was found more accurate when the Ct value was less (usually < 30), and false-positive cases would be higher if the PCR test was used as the sole gold standard. SYBR Green technique was suitable for detecting SARS-CoV-2 in clinical samples; and our user-friendly and affordable protocol detected clinical SARS-CoV-2 as efficiently as the standard costly Taqman protocol. Since this SYBR green-based method of SARS-CoV-2 detection was performed well after optimization and the results can be easily predicted by melting curve analysis, we recommend this method as an easier, cheaper, and reliable alternative to the costly probe-based method to increase the testing capacity for low-and middle-income countries where reagent supply is limited and high testing capacity is desired.

Primer designing for SARS-CoV-2 strains. We designed our primers by targeting N, E, S, and RdRp protein coding nucleotide sequences based on the aligned sequence data of all circulating SARS-CoV-2 genome sequences available in GISAID. We used CoVariants (https:// covar iants. org/ varia nts) and National Genomics Data Center (https:// bigd. big. ac. cn/ ncov/ varia tion/ annot ation) for checking whether the mutation sites fall within the primer binding region or not. We also checked the primers against the genome of Homo sapiens, six common human coronaviruses HCoV-OC43, HCoV-229E, HCoV-NL63, HCoV-HKU1, MERS-CoV, and SARS-CoV, and finally main respiratory and opportunistic viruses and pathogens in PrimerBLAST (Supplementary Table S1 ). As the internal control, we primarily targeted three house-keeping genes of the human genome, GAPDH and beta-actin. We designed the SARS-CoV-2 specific primers in a way to maintain the annealing temperature at ~ 60 °C for efficient qPCR detection and then used oligo-analyzer tool to check the primers for stem-loops and highly energetic dimer formation (< − 9 kcal/mol). The amplicons' size was distinct so that we could identify separate products in the melting curve (Table 1) .

Clinical sample selection and ethical consideration. The sample size was calculated based on the prevalence of positive samples using the following equation Randomly (random number generator using Microsoft Excel inc.) selected 6 positive and 4 negative samples were selected from left-over samples once in a week and tested by our proposed method besides the routine TaqMan based method. In 9 weeks from 10 October 2020 to 30 November 2020, we collected a total sample size of 3836 where 560 samples were positive. Using the above equation with a 10% margin of error and 95% confidence level for the samples and finally a total of 90 nasopharyngeal swab samples (sample set A) were selected. Similarly, another set of 90 samples (sample set B) was collected in between August and September of 2021 from 12,154 samples (2330 positive cases). All patients' samples were selected from the continuous surveillance of COVID-19 at the Genome Center, Jashore University of Science and Technology (JUST). The study was approved by the ethical review committee (ERC) of Jashore University of Science and Technology, Bangladesh (Reference: ERC/FBS/JUST/2020-45, Date: 06/10/2020). We performed all experiments according to the relevant guidelines and regulations. The participants were informed about the study, provided their informed consent, and we used their left-over samples from the routine surveillance.

Crude RNA extraction. According to the manufacturer's instruction, the viral RNA was extracted from the sample set' A (Appendix 1) using QuickExtract™ RNA Extraction Kit (Lucigen, USA) (cat number: QER090150). Briefly, a particular volume of the samples (5-100 µl) VTM was separated from patients swab sample inside Biosafety Cabinet (Class II) and mixed with an equal volume of ice-cold QuickExtract RNA Extraction Solution (5-100 µl) in Eppendorf tube and vortex-mixed for 1 min, followed by immediate transfer into ice.

Column based RNA extraction. To better ensure our method compared to the commercial kit, we extracted RNA from sample set B (Appendix 2) using a more efficient column-based RNA extraction system using AFCPrep™ Viral RNA extraction Kit (cat. AFC-VRNA-048), according to the manufacturer's instructions. Briefly, after vortexing the sample tube, we transferred 140 µl of sample into a nuclease-free Eppendorf tube and then we added 560 µl of AFCL buffer, vortexed, and incubated at room temperature for 10 min. After that, we added 560 µl ethanol (96 ~ 100%) in the tube and mixed well by brief vortexing. Afterward, we combined a RNA extraction (RE) Column with a Collection Tube and transferred up to 700 µl of sample mixture to the RE Column, centrifuged for 1 min at 12,000 RPM, discarded the flow-through, and transferred the remaining sample mixture to the RE Column and centrifuging under the same conditions as before. We then discarded the flow-through, and reused the Collection Tube. Then, we performed similar steps with Wash Buffer-1(500 µl), Wash Buffer-2 (700 µl, washed twice), and blank RE Column to dry. Finally, we combined the RE Column with www.nature.com/scientificreports/ a new Elution Tube, added 55 µl of RNase-free water, and centrifuged at 15,000 RPM for 3 min to elute the viral RNA.

Commercial fluorescence-based RT-PCR. In TaqMan probe-based RT-qPCR method, commercially available SARS-CoV-2 nucleic acid detection kits (Sansure Biotech, China) were compared with the in-house SYBR Green kit. The kit contains SARS-CoV-2 ORF1ab and N genes, and as internal control, human IRC genes (i.e., Rnase P). The reaction conditions and procedures were applied according to the protocol described elsewhere and all reactions were performed in duplicate to confirm reproducibility. of the amplicons to optimize the PCR conditions and primer set concentration. The target SARS-CoV-2 genes included N1, S1, E1, and RdRp1. In addition, either GAPDH or β-actin of humans, a housekeeping gene, was used as an internal positive control. For each reaction, 5 μl of extracted RNA template was used for the SARS-CoV-2 specific target primer sets and human internal positive control primer set. Gradient RT-qPCR was performed with an increasing annealing temperature (Tm) from 60 to 70 °C based on the melting temperature of each primer set. Different primer concentrations (forward, reverse; 50 nM, 75 nM; 100 nM, 150 nM; 200 nM, 250 nM)were used to optimize the engagement of the primer sets. The basic RT-qPCR conditions were set according to the Luna® Universal SYBR green One-Step RT-qPCR Kit (New England Biolabs Inc, MA). Detailed protocol is provided in the manufacturer's recommended procedures (https:// inter natio nal. neb. com/ produ cts/ e'005-luna-unive rsal-one-step-rt-qpcr-kit# Produ ct% 20Inf ormat ion). We mixed 10 μl master mix (2X) with 1 μl enzyme mix (20X), 200 nM forward primer, and 250 nM reverse primer for each gene in single primer set. The PCR conditions were set at 58 °C incubation for 12 min for reverse transcription followed by the initial denaturation at 95 °C for 3 min, then cycle denaturation at 95 °C for 15 s and an extension at 62 °C for 25 s for 45 cycles. For the melt curve analysis, we set 95 °C for 15 s, then 60 °C for 1 min followed by 95 °C for 15 s. The ramp rate of last transformation of 60 °C to 95 °C was set at 0.05 °C/seconds. For passive reference, ROX was used as the passive reference and SYBR dye was used to check the fluorescence in the Real Time PCR machine (QuantStudio 3.0, Applied Biosystem). Low Ct valued samples (< 20 as detected in TaqMan method) and nuclease-free water were included in every qPCR run as a positive control and negative control, respectively.

Optimization of multiplex RT-PCR. For the multiplexing, two sets of SARS-CoV-2 specific primers, JUST_N1 and JUST_E1, as well as, β-actin (human control) primer set were used since the amplified products had distinct melting curve peaks-that eased the identification of different products. The primer concentration was optimized for multiplex RT-qPCR as follows: JUST_N1 forward 300 nM, JUST_N1 reverse 400 nM, JUST_E1 forward 150 nM, JUST_E1 reverse 175 nM, β-actin forward 150 nM, and β-actin reverse 200 nM. 10 μl master mix (2×) with 1 μl enzyme mix (20X) of Luna® Universal One-Step RT-qPCR Kit (New England Bioscience) with 5 μl of template RNA was used under the PCR conditions as performed for the singleplex one. We set the criteria to determine the positive and negative results: presence of any peak for either N or E genes in the melting curve. The results were not interpreted based on the derivative reporter value of the melting curve for each gene since a peak with low value but good shape was considered in this case. The overall protocol of the assay has been presented within a single kit-based manual style (Fig. 3) .

Generation of standard curve and melting curve analysis. For preparing the standard curve, template RNA was diluted tenfold for 5 times (10 0 to 10 -5 ). The diluted template RNA was used to conduct both single-plex RT-qPCR assay for JUST_N1 and JUST_E1. The derivative reporter values of the melting curve as measured for the N and E genes were analyzed to compare between the SYBR positive-Sansure positive and SYBR positive-Sansure negative groups. We targeted to set a threshold based on the derivative reporter value from the melting curve.

Assay precision and validation of RT-PCR method. We performed technical replicates (triplicate) for the samples that were extracted by column purification method. To ensure amplification of the correct RT-PCR Automated RNA extraction protocol for the contradictory samples. For the contradictory samples (false-positive and false-negative), we here used automated RNA extraction for better efficiency. MagMAX™ Viral/Pathogen Nucleic Acid Isolation Kit (Cat. No. A42352) for better extraction of contradictory samples according to the manufacturer's instructions. Briefly, we added 10 µl Proteinase K, 1000 µl Wash buffer, 1000 µl 80% ethanol, and 500 µl 80% ethanol to each well of different rows of the sample plate. Afterwards, we added 200 µl sample and 550 µl Binding Bead Mix (530 µl Binding Solution + 20 µl Total Nucleic Acid Magnetic Beads) into wells of a row A and tip comb in the last row of the sample plate. Besides, we added 40 µl Elution solution to each well of the elution strip. Selecting the program MVP_DUO on the instrument KingFisher™ Duo Prime and starting the run, then loaded the Elution Strip and Sample Plate into position when prompted by the instrument. After the protocol was complete (~ 25 min after start), we immediately removed the Elution Strip from the instrument and recovered RNA elute.

Pneumonia of unknown aetiology in Wuhan, China: potential for international spread via commercial air travel

Analytical sensitivity and efficiency comparisons of SARS-CoV-2 RT-qPCR primer-probe sets

Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR

Types of assays for SARS-COV-2 testing: A review

Comparison of SYBR Green and TaqMan methods in quantitative real-time polymerase chain reaction analysis of four adenosine receptor subtypes

A new SYBR Green real-time PCR to detect SARS-CoV-2

Lower cost alternatives for molecular diagnosis of COVID-19

A one-step, one-tube real-time RT-PCR based assay with an automated analysis for detection of SARS-CoV-2

SYBR green one-step qRT-PCR for the detection of SARS-CoV-2 RNA in saliva

Optimization of primer sets and detection protocols for SARS-CoV-2 of coronavirus disease 2019 (COVID-19) using PCR and real-time PCR

Development of a laboratory-safe and low-cost detection protocol for SARS-CoV-2 of the coronavirus disease 2019 (COVID-19)

Evaluation of SYBR Green real time PCR for detecting SARS-CoV-2 from clinical samples

Evaluation of the diagnostic performance of nine commercial RT-PCR kits for the detection of SARS-CoV-2 in Colombia

Current nucleic acid extraction methods and their implications to point-of-care diagnostics

The architecture of SARS-CoV-2 transcriptome

RdRp mutations are associated with SARS-CoV-2 genome evolution

Tracing the SARS-coronavirus

Centre for Disease Control and Prevention. SARS-CoV-2 Variant Classifications and Definitions

Color multiplexing hybridization probes using the apolipoprotein E locus as a model system for genotyping

KRAS genotyping by digital PCR combined with melting curve analysis

Development of Sybr green I-based melting curve method for HER2I655V polymorphism detection in breast cancer

PCR inhibition in qPCR, dPCR and MPS-mechanisms and solutions

Comparison of the diagnostic efficacy between two PCR test kits for SARS-CoV-2 nucleic acid detection

Quantitative real-time PCR is not more sensitive than "conventional

Validation of conventional PCR-like alternative to SARS-CoV-2 detection with target nucleocapsid protein gene in naso-oropharyngeal samples

Evolutionary dynamics of SARS-CoV-2 nucleocapsid protein and its consequences

PCR inhibitors-occurrence, properties and removal

A novel point mutation in the N gene of SARS-CoV-2 may affect the detection of the virus by reverse transcription-quantitative PCR

A rapid and cost-effective multiplex ARMS-PCR method for the simultaneous genotyping of the circulating SARS-CoV-2 phylogenetic clades

Sanger sequencing of the amplified products (spike gene). For purifying the desired PCR product from the low melting 1% agarose gel, we used Wizard® SV gel and PCR Clean-Up system. The representative amplicons were then subjected to Sanger sequencing with BigDye™ Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific) in Applied Biosystems SeqStudio genetic analyzer as per the optimized protocol of Islam et al. (2021) . The trace or chromatogram (.ab1) files from the Sanger sequencing were analyzed using the Sequencing Analysis Software V6.0 (Thermofisher, USA). NCBI BLAST was performed initially and the alignment to SARS-CoV-2 spike gene was also checked in MEGA7 (https:// www. megas oftwa re. net/). Statistical analysis. Sensitivity, Specificity, kappa index, Positive predictive value (PPV) and Negative predictive value (NPV) were calculated by analyzing data in IBM SPSS Statistics. Positive Likelihood ratio (LR +) and Negative likelihood ratio (LR-) were calculate using MedCalc Software Ltd. Diagnostic test evaluation calculator. https:// www. medca lc. org/ calc/ diagn ostic_ test. php (Version 20.027; accessed February 16, 2022). Simple linear regression analysis and standard curve were generated RStudio using packages: readxl, ggplot2, ggpubr and ggpmisc.

We acknowledge GISAID for sharing the sequence data and IDT for giving the opportunity to use the tools for validating primers in silico. We also acknowledge the Ministry of Health and Family Welfare, Bangladesh, for giving us permission for the SARS-CoV-2 diagnosis.

This study was funded by Jashore University of Science and Technology Research Grant supported by the University Grants Commission of Bangladesh (#FoBST-06).

The authors declare no competing interests.

Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-022-10413-7.Correspondence and requests for materials should be addressed to I.K.J. or M.A.H.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.