Extraction-free protocol combining proteinase K and heat inactivation for detection of SARS-CoV-2 by RT-qPCR

About the Authors:

Valeria Genoud

Contributed equally to this work with: Valeria Genoud, Martin Stortz

Roles Conceptualization, Investigation, Methodology, Project administration, Supervision, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (VG); vlevi12@gmail (VL)

Affiliation: Facultad de Ciencias Exactas y Naturales, Departamento de Química Biológica, Universidad de Buenos Aires, Buenos Aires, Argentina

Martin Stortz

Contributed equally to this work with: Valeria Genoud, Martin Stortz

Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliations Facultad de Ciencias Exactas y Naturales, Departamento de Fisiología, Biología Molecular y Celular, Universidad de Buenos Aires, Buenos Aires, Argentina, Facultad de Ciencias Exactas y Naturales, Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

Ariel Waisman

Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: LIAN-CONICET, FLENI, Belén de Escobar, Argentina

Bruno G. Berardino

Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliations Facultad de Ciencias Exactas y Naturales, Departamento de Química Biológica, Universidad de Buenos Aires, Buenos Aires, Argentina, Facultad de Ciencias Exactas y Naturales, Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

ORCID logo https://orcid.org/0000-0002-8954-5546

Paula Verneri

Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliation: Facultad de Ciencias Exactas y Naturales, Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

Virginia Dansey

Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

Affiliations Facultad de Ciencias Exactas y Naturales, Departamento de Química Biológica, Universidad de Buenos Aires, Buenos Aires, Argentina, Facultad de Ciencias Exactas y Naturales, Unidad de Microanálisis y Métodos Físicos en Química Orgánica (UMYMFOR), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

ORCID logo https://orcid.org/0000-0003-3611-6577

Melina Salvatori

Roles Investigation, Methodology

Affiliation: Instituto de Investigaciones Biomédicas en Retrovirus y SIDA (INBIRS), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

Federico Remes Lenicov

Roles Conceptualization, Investigation, Methodology, Supervision, Validation, Writing – original draft, Writing – review & editing

Affiliation: Instituto de Investigaciones Biomédicas en Retrovirus y SIDA (INBIRS), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

Valeria Levi

Roles Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (VG); vlevi12@gmail (VL)

Affiliations Facultad de Ciencias Exactas y Naturales, Departamento de Química Biológica, Universidad de Buenos Aires, Buenos Aires, Argentina, Facultad de Ciencias Exactas y Naturales, Instituto de Química Biológica de la Facultad de Ciencias Exactas y Naturales (IQUIBICEN), CONICET-Universidad de Buenos Aires, Buenos Aires, Argentina

ORCID logo https://orcid.org/0000-0002-1666-4675

Introduction

The outbreak of the coronavirus disease 2019 (COVID-19), caused by SARS-CoV-2, was declared a pandemic by World Health Organization on March 11, 2020. Confronted by the lack of a vaccine or specific treatment, early detection, isolation and contact tracing are fundamental strategies for reducing the spread of this disease [1].

RT-qPCR is the gold standard method for SARS-CoV-2 determinations [2]. These assays typically involve collecting the clinical specimen usually with nasopharyngeal swabs, extracting the RNA from the sample and analyzing the presence of the viral RNA through RT-qPCR.

Nowadays, the RNA extraction step is considered a bottleneck of SARS-CoV-2 RT-qPCR determinations [3]. This procedure can be automatized but commercial extraction robots are very expensive and may not be affordable for most clinical labs in low-income countries. In contrast, manual RNA-extraction is a quite cumbersome procedure that includes several washing and centrifuging steps at a Biosafety level 2 laboratory and relatively expensive extraction kits. In Argentina, a trained lab technician processes approximately 24 samples every 1.5 hours. Additionally, the skyrocketing increment on the clinical use of RT-qPCR due to the COVID-19 pandemic led to a worldwide shortage of RNA purification reagents [4,5] thus, many research groups around the world have concentrated their efforts on either simplifying or eliminating the extraction step [6–16].

Björn Reinius at the Karolinska Institute in Sweden [17] optimized an extraction-free sample preparation protocol that only includes a very simple heat-inactivation step at 95°C during 5 min (HID method). Through HID RT-qPCR, they were able to detect with high sensitivity the relatively short N1 and N2 regions of the viral RNA which are preserved after the heating step. This group also performed an exhaustive analysis of media formulation compatible with the HID method [17]. More recently, the Centers for Disease Control and Prevention (CDC) has approved a similar thermal-inactivation protocol (1 min at 95°C) that can be used with Quantabio UltraPlex 1-Step ToughMix (4X) and is sufficient for viral inactivation [18]. This protocol also detects N1 and N2 sequences of the viral RNA.

Many protocols traditionally used for extraction-free preparation of clinical samples included treatments with proteinase-K (PK), a commonly used protease that degrade RNases and hence, it helps to preserve RNA integrity [19,20]. A recent work suggested that PK may facilitate extraction-free determinations of SARS-CoV-2 although a major loss in sensitivity with a ~6-units shift in CT values on the determination of the envelope (E) gene [21] may prevent its uses in clinical analysis. Other recent preprint [22] tested a protocol based on PK and heat inactivation using a reduced number of samples (N = 17) in UTM medium with relative low CT values, and claim a good sensitivity on the determination of the E gene. Although these results might seem contradictory with recent findings showing suboptimal detection of the E amplicon on HID samples in UTM medium [17], the efficiency of the PCR is very dependent on the detection kit, transport medium brand, thermocycler, etc.

Here, we combined these strategies and optimized the PK+HID method, a very simple protocol for extraction-free RT-qPCR determinations of SARS-CoV-2. Our results position PK+HID as an inexpensive, versatile and reliable alternative to bypass the traditional RNA extraction procedure and suggest that it could substantially contribute to improve the testing capacity in countries with limited resources and facilities.

Materials and methods

Sample collection

Nasopharyngeal swabs were collected during screening of non-hospitalized symptomatic cases and close contacts to these cases at diverse hospital and clinical centers at the Buenos Aires area (AMBA, Área Metropolitana de Buenos Aires), Argentina. Swabs were deposited in 2–5 ml of saline solution. These clinical samples were processed at the Instituto de Investigaciones Biomédicas en Retrovirus y SIDA (INBIRS, Buenos Aires, Argentina), a specialized center currently dedicated to SARS-CoV-2 diagnosis by RT-qPCR. For routine diagnosis, RNA extraction was performed from 300 μl of swab samples using a Chemagic 360-D (Perkin-Elmer), an automated extraction equipment, according to the manufacturer’s instructions.

This study used the remaining volume of anonymized samples that had been collected for clinical diagnosis. Under these circumstances, and for studies involving development of COVID-19 diagnostic tools, our IRB (the Research Ethics Committee from Fundación Huésped in Buenos Aires, Argentina) deemed unnecessary to obtain informed consent from the patients.

PK+HID sample preparation

Nasopharyngeal swabs deposited in 2–5 ml of saline solution were maintained up to 24 h at 4°C. PK+HID samples were prepared 10–24 h after RNA extraction unless indicated otherwise. 10 μl of a solution 10 mg/ml Proteinase K (Promega, Madison, WI, USA or FlashPrep® ARN SARS-CoV-2 Highway, Inbio Highway, Tandil, Argentina) were placed in 0.2 ml PCR tubes. Then, 90 μl of the nasopharyngeal swabs samples previously vortexed were transferred to each of these PCR tubes. Please, notice that the PK solution should be added to the tubes before transferring the samples and not the other way around. This procedure minimizes the work-hazard since tubes with samples are not reopen until inactivated. It also provides flexibility in the workload of clinical labs because the preparation of PCR tubes with PK solution could be done outside the Biosafety Level 2 facility in a clean, nucleic acid-free area.

For experiments with variable PK concentrations, 10 μl of PK stock solutions with the adequate concentration were added. Samples were incubated at 55°C for 15 min and then at 98°C for 5 min in a thermal cycler with heated lid. Finally, inactivated samples were cooled to 4°C and kept at this temperature until RT-qPCR analysis.

RT-qPCR

One-step reverse transcription and qPCR of N1 and N2 viral regions were performed using the GoTaq® Probe 1-Step RT-qPCR System (Promega, Madison, WI, USA) and the 2019-nCoV CDC Diagnostic Panel (Integrated DNA Technologies, Coralville, IA, USA) according to manufacturer’s instructions. The N1 and N2 viral regions and the human RP (internal control) were detected in singleplex 20 μl-volume reactions using 5 μl of purified RNA or 5 μl of PK+HID samples unless specifically stated in the text. Thermal cycling steps were: 45°C for 15 min, 95°C for 2 min, and 45 cycles of 95°C for 3 s and 55°C for 30 s. RT-qPCR was performed on a StepOnePlus Real-Time PCR System (Applied Biosystems).

One-step reverse transcription and qPCR of N and ORF1ab viral genes were performed following routine diagnostics procedures at INBIRS using the DisCoVery SARS-CoV-2 RT-PCR Detection Kit (Safecare Biotech Hangzhou, China) according to manufacturer’s instructions. The N and ORF1ab viral regions and the human RP internal control were detected in multiplex reactions. Thermal cycling steps were: 50°C for 5 min, 95°C for 30 s, and 45 cycles of 95°C for 5 s and 60°C for 34 s. RT-qPCR was performed on a CFX96 Real-Time System.

One-step reverse transcription and qPCR of N, E and RdRp viral genes were performed following routine diagnostics procedures at INBIRS using the GeneFinderTM COVID-19 Plus RealAmp Kit (Gene Finder, Korea) according to manufacturer’s instructions. The N, E, RdRp viral regions and the human RP internal control were detected in multiplex reactions. Thermal cycling steps were: 50°C for 20 min, 95°C for 5 min, and 45 cycles of 95°C for 15 s and 58°C for 1 min. RT-qPCR was performed on a CFX96 Real-Time System (BioRad).

Infectivity assays

Vero E6 cells were grown in Dulbecco’s modified Eagle’s medium (DMEM) supplemented with 10% fetal bovine serum, 1% penicillin/streptomycin, 0.5 mg/ml Amphotericin B and 1% L-glutamine (all from Gibco). Cells were seeded on a 24-well plate at 3.5 × 105 cells/ml one day prior to infection. Nasopharyngeal swab samples were diluted 1:2 in DMEM, filtered through a 0.22 mm-pore-filter and added to cell cultures for 1 h. Then, cultures were washed and incubated for 3 days. When the cultures showed cytopathic effects, the supernatant was collected for detection of SARS-CoV-2 by RT-qPCR, and the culture was terminated. In the absence of cytopathic effects, the supernatant was passed to a new culture. When no cytopathic effect was observed after the 3rd passage, the supernatant was collected for detection of SARS-CoV-2 by RT-qPCR. All infection experiments were conducted in a biosafety level 3 laboratory.

Loop-mediated isothermal amplification (LAMP)

Reverse transcription, isothermal amplification of the viral genes ORF1a, ORF1b, N and E, and colorimetric detection were performed with COVID-19 Neokit Tecnoami (Neokit S.A.S., Buenos Aires, Argentina), according to manufacturer’s instructions. Briefly, reactions were performed in 0.2 ml PCR tubes mixing 35 μl of reagent and 8 μl of the nasopharyngeal samples. The reactions were run in a thermal cycler with heated lid at 64°C for 1 h. Diagnostics was made according to final color of the solution (violet, negative; blue, positive, S1 Fig). The limit of detection reported by the provider is 12.5 copies whereas the clinical sensitivity, specificity and accuracy are 94, 100 and 96%, respectively.

Data analysis

qPCR amplification curves obtained in those experiments using DisCoVery and GeneFinder detection kits were analyzed with the CFX Manager software (BioRad) to obtain CT values. Data obtained with 2019-nCoV CDC were analyzed with the free software LinRegPCR [23] to obtain CT values. This software was also used to calculate PCR efficiencies and the initial amount of template copies (in arbitrary units).

Results

The treatment of nasopharyngeal swab samples with proteinase-K (PK) improves the performance of the heat inactivation protocol

We used a reduced number of SARS-CoV-2 positive (N = 6) and negative (N = 3) nasopharyngeal swab samples collected in saline solution to test if their preincubation with PK before the thermal inactivation step improves the performance of the HID extraction-free method. Although previous results [17] suggested that saline solution is suboptimal for direct RT-qPCR detection, this is one of the transport medium recommended by CDC and it is also commonly used to collect clinical samples in Argentina.

In these experiments, we mixed up 90 μl of the samples with 10 μl of PK 10 mg/ml (final concentration 1 mg/ml) in 0.2 ml PCR tubes and incubated them in a thermocycler during 15 min at 55°C followed by 5 min at 98°C (PK+HID samples). As controls, we used heat-inactivated aliquots omitting the preincubation with PK at 55°C (HID samples) and RNA samples obtained through a standard extraction protocol (purified RNA samples) from the same nasopharyngeal swabs. These samples were analyzed by RT-qPCR using the 2019-nCoV CDC Diagnostic Panel which detects the N1 and N2 viral regions and the human RNase P (RP) gene as internal control.

Fig 1a shows that all positive and negative samples were identified using either HID or PK+HID treatments. As examples, Fig 1b shows representative amplification curves obtained for a positive PK+HID sample illustrating that the complex matrix of PK+HID samples did not affect the expected sigmoidal-like shape of the curves. Additionally, the pretreatment with PK consistently resulted in lower CT values than those obtained with HID samples and closer to those determined with purified RNA samples (Fig 1a and S1 Table).

[Figure omitted. See PDF.]

Fig 1. Proteinase K improves the performance of the heat inactivation method in RT-qPCR determinations of SARS-CoV-2.

(a-b) Positive (#8, #11, #12, #16, #19, #20) and negative (#1, #2, #3) nasopharyngeal swab samples were processed by heat inactivation (HID, 98°C for 5 min); proteinase K treatment followed by heat inactivation (PK+HID, 55°C for 15 min and 98°C for 5 min) or subjected to RNA extraction (purified RNA). The viral N1 and N2 genes and the human RNase P gene (RP) were amplified and detected by RT-qPCR. (a) CT values obtained from RT-qPCR analysis of the same samples prepared by the three different methods. (b) Representative amplification curves for each gene obtained for one of the positive samples. (c) Amplification efficiencies (EPCR) and initial amount of amplicon copies (n) in PK+HID samples relative to the corresponding HID samples. The median of each measurement is represented with a line in the bars and the lengths of these bars represent the standard error. (d) Positive nasopharyngeal swab samples were subjected to treatment with different concentrations of proteinase K (PK) followed by heat inactivation (55°C for 15 min and 98°C for 5 min). The viral N1 and N2 genes and the human RP gene were amplified and detected by RT-qPCR. ΔCT represents the mean difference between CT values obtained in each analyzed condition and the corresponding to HID samples. Mean ± SEM values are represented (N = 5).

https://doi.org/10.1371/journal.pone.0247792.g001

The analysis of the amplification curves shows that the PCR efficiency (EPCR) was the lowest for the HID treatment followed by PK+HID and then by the RNA extraction treatment (Fig 1c). This result suggests that the PK treatment may contribute to degrade unknown components of the sample matrix that interfere with RT-qPCR. Additionally, we observed a higher PCR amplification efficiency for the N1 viral gene compared to N2, in line with previous results [17,24]. We also analyzed the amount of amplicon copies in PK+HID relative to that determined in HID samples and observed that this ratio was higher than 1 for N2 and RP suggesting that PK treatment also protected these RNA regions from RNase´s action (Fig 1d).

PK+HID involves the addition of the PK buffer and thus a 10% dilution of the saline medium used to collect nasopharyngeal swabs. Therefore, we next evaluate if the improvement of PK+HID treatment compared to heat-inactivation is only due to the dilution of the saline medium. We analyzed five positive samples processed by the standard PK+HID protocol or by HID but adding to the samples the same amount of PK buffer without PK (HID´ samples). S2 Fig shows that CT values for PK+HID samples were lower compared to HID´ samples; the analyses of amplification curves indicate that nPK+HD was higher than nHID´ supporting that PK contributes with the preservation of the integrity of the viral RNA. Therefore, the dilution of saline solution could not completely explain the better performance of PK+HID. Nevertheless, the efficiency of PCR amplification for HID´ samples was slightly lower only for N2 amplicon suggesting that the dilution of saline solution in PK+HID samples also contributes to increase the efficiency of PCR.

We next tested if using a higher volume of PK+HID samples could improve SARS-CoV-2 detection in the RT-qPCR determinations. Although the sensitivity of the assay is expected to be higher with more copies of the viral RNA template, a previous report also suggested that increasing too much this volume may have the opposite effect due to RT-qPCR inhibitors in the sample matrix [17]. We observed that the amplification curves shifted to lower cycles when using 10 μl of sample instead of 5 μl for the RT-qPCR determination resulting in lower CT values for the N1 region (0.3 to 1 units, S2 Table). This result suggests that a higher sample volume could contribute to the detection of SARS-CoV-2 in samples with low viral titers.

Finally, we tested if the performance of the PK+HID method depends on PK concentration in the range of 0.1 to 1 mg/ml. Fig 1d shows that CT values obtained for the viral and human amplicons did not change within this PK concentration range. These values were also lower than those observed in the absence of PK, confirming the improving role of PK in this extraction-free method.

A comparative RT-qPCR analysis with purified RNA samples shows an efficient detection of N1 SARS-CoV-2 region in PK+HID samples

We next ran RT-qPCR determinations of SARS-CoV-2 in PK+HID and purified RNA samples prepared from positive and negative nasopharyngeal swab samples. Based on the results showed in the previous section, we used 1 mg/ml PK for the PK-treatment and 8 μl of PK+HID samples (i.e. between the volumes assayed in the previous section) in the following RT-qPCR assays. For this particular experiment, we only analyzed the N1 viral gene and the RP human internal control to maximize the number of samples loaded in the PCR multiwell plate.

In these experiments, we correctly identified all positive (27/27) and negative (6/6) samples. Fig 2 also shows that the CT values determined for the viral sequence N1 in PK+HID positive samples strongly correlated with those values measured in purified RNA samples supporting that the PK+HID treatment robustly preserves the information on the viral loads in the swab specimens.

[Figure omitted. See PDF.]

Fig 2. PK+HID method exhibits a similar performance than RNA extraction in RT-qPCR determinations of the SARS-CoV-2 N1 gene.

Positive nasopharyngeal swab samples were processed for RNA extraction (purified RNA) or subjected to treatment with PK 1 mg/ml (55°C for 15 min) followed by heat inactivation at 98°C for 5 min (PK+HID). The viral N1 gene and the human RP gene were amplified and detected by RT-qPCR. CT values obtained for the same samples prepared by the two different methods (N = 27), the line obtained by regression of the data (continuous lines) and 95% confidence bands (pink) are represented. The parameter values obtained from the fitting were: slope = 1.17 ± 0.05 and intercept = -3.4 ± 1.2 (N1 amplicon); slope = 1.0 ± 0.1 and intercept = 1.7 ± 3.0 (RP amplicon).

https://doi.org/10.1371/journal.pone.0247792.g002

The PK+HID method can be used in combination with different multiplexing RT-qPCR detection kits allowing the detection of SARS-CoV-2 with high sensitivity and specificity

In the previous sections, we showed that the PK+HID method resulted in an efficient alternative to bypass RNA extraction when detecting the N1 and N2 regions of the viral RNA. However, most laboratories use PCR kits that target other viral genes often combined in multiplex setups. Therefore, we next evaluated if the improved PK+HID method could be used with multiplexing detections kits.

We first used DisCoVery (Safecare Biotech Hangzhou), a SARS-CoV-2 multiplex RT-qPCR detection kit that targets the SARS-CoV-2 genes N and ORF1ab and human RNase P gene. Relevantly, many worldwide providers offer detection kits targeting one or both of these genes.

Fig 3 shows that the PK+HID treatment of the samples allows the detection of viral and internal control genes. In the case of viral amplicons, the CT values obtained with PK+HID samples linearly correlated with those obtained with purified RNA samples. The slopes for both N and ORF1ab were ~1 and the y-intercepts were ~0 suggesting that these sequences are detected equally well in samples prepared with the PK+HID method or through RNA extraction. This observation is also valid for samples with relatively high CT values (>35) indicating that weak positive samples are still detected with the PK+HID method.

[Figure omitted. See PDF.]

Fig 3. PK+HID method exhibits a similar performance than RNA extraction in RT-qPCR determinations of the SARS-CoV-2 N and ORF1ab genes.

Positive nasopharyngeal swab samples were processed for RNA extraction (purified RNA) or subjected to treatment with PK 1 mg/ml (55°C for 15 min) followed by heat-inactivation at 98°C for 5 min (PK+HID). The viral N and ORF1ab genes and the human RP gene were amplified and detected by RT-qPCR. (a) CT values obtained for the same samples prepared by the two different methods, the line obtained by regression of the data (continuous lines) and 95% confidence bands (pink) are represented. The parameter values obtained from the fitting were: slope = 1.02 ± 0.04 and intercept = 0.1 ± 0.9 (N amplicon); slope = 1.02 ± 0.04 and intercept = 0.8 ± 0.9 (ORF1ab amplicon) and slope = 0.4 ± 0.2 and intercept = 15 ± 6 (RP amplicon). Notice that CT values determined for RP spans a smaller range. (b) Amplification efficiencies obtained in PK+HID samples (EPK+HID) relative to the corresponding purified RNA samples (Epurified RNA). The mean values obtained for Epurified RNA were: 1.98 ± 0.03 (N), 2.02 ± 0.06 (ORF1ab) and 1.94 ± 0.03 (RP). The median of each measurement is represented with a line in the bars and the lengths of these bars represent the standard error (N = 16).

https://doi.org/10.1371/journal.pone.0247792.g003

We correctly identified 44/46 positive and 3/3 negative samples. S2 Table compiles the CT values of the false-negative samples (see Discussion). Additionally, the internal control (IC) of two other PK+HID samples was not detected making mandatory repeating the assay; these invalid samples were not included in S3 Table. We determined the amplification efficiency of the viral regions N and ORF1ab and the human gene RP and found that these values were similar for PK+HID and purified RNA samples (Fig 3b). We also compared the amplification curves obtained using this detection kit and swab samples that were either treated by HID or PK+HID or subjected to the RNA extraction procedure (S3 Fig). PK treatment lowered the CT values in comparison to those obtained with HID only, with a mean + SEM difference in CT values of 2.7 + 0.7 (N), 3.1 + 1.0 (ORF1ab) and 1.4 + 0.2 (RP). These results, obtained with a different detection kit, agree with those showed in Fig 1 also supporting that the PK treatment improves the performance of the heat-inactivation method.

The work of Reinius and colleagues [17] confirmed the virus inactivation after HID treatment. As a control, we also run infectivity assays as described in Materials and Methods using 4 nasopharyingeal swab samples with relatively high viral loads (CT values 22.11; 21.32; 22.74 and 19.18 for N amplicon). After the 3rd passage of cell culture, we did not verify any cytopathic effect and the viral RNA was not detected by RT-qPCR showing that the PK+HID treatment also inactivates the virus.

We also analyzed if PK+HID samples can be preserved at 4°C since this could be a common practice in clinical labs. We processed nasopharyngeal samples using the PK+HID treatment and run RT-qPCR determinations after different incubation times at 4°C (0, 20 and 40 h). S4 Fig shows that the CT values for N and ORF1ab genes did not change during the first 20 h and increased at the longest incubation time suggesting that PK+HID samples can be conserved at 4°C for up to 20 h.

A common practice in many clinical labs involves freezing/thawing reagents several times. Therefore, we tested if repeated cycles of freezing/thawing of the PK stock solution may affect the performance of PK+HID. To address this issue, we processed 10 positive samples by PK+HID using a PK solution that was previously subjected to 10 cycles of freezing/thawing or a PK solution thawed only once. The mean differences in CT values between these samples were 0.17 ± 0.15 (N), -0.05 ± 0.20 (ORF1ab) and 0.01 ± 0.14 (RP) (Ndata = 10) showing that the PK stock solution is stable for up to 10 cycles of freezing/thawing.

To analyze if the PK+HID method can be combined with other multiplexing kits, we prepared PK+HD and purified RNA samples from 94 randomly-selected clinical specimens and ran RT-qPCR determinations using the GeneFinder detection kit which amplifies N, E and RdRp regions of the viral RNA and the human RP gene as internal control.

The analysis of CT values measured in PK+HID and RNA samples showed stronger correlations for the N amplicon in comparison to E and RdRp amplicons (Fig 4a). Additionally, some of the N-positive RNA samples (also considered SARS-CoV-2 positive in clinical diagnosis [24]) showed negative results for the E and RdRp amplicons (3/38 and 4/38, respectively); this difference increased in PK+HID samples with 8 more samples showing negative results for both amplicons while being N-positive. As mentioned above, other authors also reported a poor detection and large right-shifts in CT values for E amplicon in extraction-free samples [17].

[Figure omitted. See PDF.]

Fig 4. Performance of PK+HID method in RT-qPCR determinations of the SARS-CoV-2 N, RdRp and E genes in randomly-selected clinical specimens.

Nasopharyngeal swab samples were processed for RNA extraction (purified RNA) or subjected to treatment with PK 1 mg/ml (55°C for 15 min) followed by heat-inactivation at 98°C for 5 min (PK+HID). The viral N, E and, RdRp genes and the human RP gene were amplified and detected by RT-qPCR. (a) CT values obtained for the same samples prepared by the two different methods, the line obtained by regression of the data (continuous lines) and 95% confidence bands (pink) are represented. RP plot also includes the data obtained in SARS-CoV-2 negative samples and spans a smaller CT range. The parameter values obtained from the fitting were: slope = 1.2 ± 0.1 and intercept = -5 ± 3 (N amplicon); slope = 0.7 ± 0.2 and intercept = 9 ±4 (RdRp amplicon); slope = 0.7 ± 0.1 and intercept = 7 ± 3 (E amplicon) and slope = 0.4 ± 0.1 and intercept = 16 ± 3 (RP amplicon). (b) Amplification efficiencies obtained in PK+HID samples (EPK+HID) relative to the corresponding purified RNA samples (Epurified RNA). The mean values obtained for Epurified RNA were: 1.94 ± 0.02 (N), 1.97 ± 0.02 (RdRp), 1.93 ± 0.03 (E) and 1.90 ± 0.03 (RP). The median of each measurement is represented with a line in the bars and the lengths of these bars represent the standard error (N = 8).

https://doi.org/10.1371/journal.pone.0247792.g004

We also observed that the amplification efficiency of viral amplicons N, E and RdRp was similar for PK+HID and purified RNA samples whereas EPCR of the internal control RP was slightly lower for PK+HID samples (Fig 4b). These results support that the PK+HID treatment contributes to inactivate inhibitors of the PCR reaction.

To quantitatively analyze the performance of the PK+HID treatment combined with the use of GeneFinder kit, we estimated the accuracy, sensitivity and specificity of detection as previously described [25] obtaining values of 0.95, 0.92 ± 0.05 and 0.96 ± 0.04, respectively.

In the validation assay described above, swab samples were first analyzed by the usual RT-qPCR procedure including the RNA extraction step; 10 to 24 h later, we processed the same nasopharyngeal swabs by PK+HID. Therefore, we asked if the extra-time that swab samples remain at 4 C could lead to a partial degradation of viral RNA and consequently, to an apparently lower sensitivity of the PK+HID method. To test this hypothesis, we performed the validation experiment running simultaneously the RNA extraction and PK-HID procedures using 106 positive and 106 negative randomly-selected samples and analyzed the RNA and PK+HID samples by RT-qPCR with the detection kit DisCoVery (S5 Fig); unfortunately, we could not run these experiments with the GeneFinder detection kit used above because the provider run out of stock during the pandemic outbreak. The analyses of the data showed that the accuracy was 0,990 and the sensitivity and specificity were both 0.99 ± 0.01. These results strongly support that PK+HID is a reliable protocol for extraction-free RT-qPCR SARS-CoV-2 determinations.

In this last validation assay, we also analyzed the amplification curves of positive samples with relatively low viral loads (CT values >30) to estimate the relative amount of viral amplicon copies in the PK+HID and RNA samples volumes used in the RT-qPCR determination. S5 Fig shows that the relative amount of viral N and ORF1ab amplicons were in the same order in PK+HID vs RNA samples; particularly, the mean ratio nPK+HID/nRNA was 0.7 ± 0.2 (N) and 1.7 ± 0.6 (ORF1ab). This result suggests a similar performance of PK+HID or RNA extraction combined with this detection kit.

Finally, we asked whether PK+HID could be combined with a SARS-CoV-2 LAMP detection kit since they are now widely used in many countries because their simplicity and associated lower costs. S1 Fig includes some examples of positive and negative results obtained combining PK+HID with LAMP. We could detect 29/29 and 5/6 positive and negative samples (previously processed by RT-qPCR with DisCoVery, 15 ≤ CTs ≤ 34); these proof-of-principle experiments suggest that PK+HID could be combined with this type of detection kits.

We also evaluated 25 positive and 53 negative samples treated either by PK+HID or RNA extraction followed by detection with LAMP (S4 Table) and estimated that the accuracy, sensitivity and specificity measured when LAMP is combined with PK+HID in comparison to RNA extraction were 0.95, (0.92 + 0.06) and (0.96 + 0.03), respectively. These results suggest that processing the swab samples with PK+HID instead of RNA extraction introduces minimal losses of sensitivity and specificity when using this LAMP detection kit.

Discussion

The gold standard method for detection of SARS-CoV-2 in nasopharyngeal swab samples is RT-qPCR, a technique with exquisite sensitivity and high specificity commonly applied to analyze purified RNA samples.

The RNA extraction step has nowadays become the bottleneck in COVID-19 testing, especially during pandemic outbreaks when personnel, resources and biosafety compatible facilities to perform the testing are scarce.

The fast evolution of COVID-19 pandemic and the necessity of performing massive and rapid RT-qPCR determinations to test, trace and isolate SARS-CoV-2 positive patients encouraged the scientific community to evaluate alternative methods to bypass the laborious RNA extraction step. We mentioned above heat inactivation methods that make use of a fast thermal-treatment of the samples and allow the identification of SARS-CoV-2 positive samples targeting N1 and N2 viral genes [17]. These results are consistent with other works that also tested the performance of heat-inactivated respiratory samples in RT-qPCR assays [8,11,26–29].

Here, we combined heat inactivation with a pretreatment with PK, protease routinely used in nucleic acid preparations and thus widely available in most clinical laboratories. We mentioned previous works using PK to treat nasopharyngeal swab samples; the protease is also starting to be used in the analysis of saliva samples for the detection of SARS-CoV-2 in extraction-free RT-qPCR determinations [21,22,30,31].

In line with these reports, our data shows that preincubation with PK improves the performance of the heat-inactivation protocol in RT-qPCR assays targeting N1 and N2 viral genes with the 2019-nCoV CDC Diagnostic Panel. However, a validation with a larger number of samples should be done to compare the sensitivities of both methods. We verified that the performance of the PK+HID does not depend on PK concentration in the range of 0.1–1 mg/ml; this robustness of the PK+HID method is appealing since small variations in the PK concentration will not significantly affect the results of SARS-CoV-2 determinations and also allows reducing even more the costs of sample preparation for RT-qPCR assays.

Using the 2019-nCoV CDC Diagnostic Panel, we could correctly identify the SARS-CoV-2 positive and negative analyzed samples and, for the positive samples, we obtained CT values in PK+HID samples very close to those of purified RNAs (2.7 units for N1) whereas heat inactivation without PK treatment produces higher shifts in CTs (5 units for N1).

Remarkably, we demonstrated that the PK+HID method could be used in combination with other RT-qPCR kits that target different viral genes in multiplex RT-qPCR reactions. Specifically, we used DisCoVery (targeting N and ORF1ab genes) and Genefinder (targeting N, E and RdRp genes). The detection of these genes was unexpected since the HID method did not provide satisfactory results in the detection of larger amplicons [17]. We found that CT values of both N and ORF1ab viral amplicons detected with DisCoVery were similar to those measured for the purified RNA samples, in line with our observations with the CDC kit.

Although the detection of the E and RdRp genes with the Genefinder kit was weaker than the detection of the N gene, the relatively high detection sensitivity of the PK+HID method estimated from the analysis of 94 randomly selected samples suggests that the PK+HID method can be used as an alternative, extraction-free RT-qPCR analysis of nasopharyngeal swab samples in clinical diagnostics.

S3 Table compiles the CT values of all false-negative samples collected during the course of our work. Most of these samples show high CT values suggesting that they correspond to swab samples with relatively low viral loads.

We should emphasize that the success of an extraction-free method strongly depends on the correct combination of transport medium and detection kit. Merindol et al [32] found that saline solution swab medium impaired the performance of the AllplexTM 2019-nCoV Assay. Although, in our hands, saline solution swab medium could be used with 2019-nCoV CDC, GeneFinder™ COVID-19 Plus RealAmp and DisCoVery SARS-CoV-2RT-PCR detection kits. Since the extraction-free protocol may reduce the efficiency of the PCR, it is critical to select a PCR kit that preserves sensitivity. While most bibliography, as well as our own findings point to N as the most sensitive viral region, ORF1ab gene also showed a good performance with our extraction-free protocol.

Conclusions

Our results indicate that the combination of proteinase K preincubation with heat inactivation is a very robust and reliable procedure for extraction-free RT-qPCR determinations of SARS-CoV-2. We should mention that preparation of the 94 PK+HID samples required to complete the 96-well PCR plate (2 wells are needed for the positive and negative controls) takes approximately 60 min by one operator whereas manual extraction of the same number of samples using commercial kits would take 5–6 hours of work in a Biosafety Level 2 area with continuous use of a benchtop centrifuge. Additionally, PK treatment involves an estimated cost of 0.37 USD per sample whereas the cost of standard RNA extraction procedures is clearly superior (6–10 USD per sample in Argentina). We should also highlight that the HID method combined with the 2019-nCoV CDC Diagnostic Panel is an inexpensive option since it does not require extra reagents as PK+HID.

The simplicity and relatively short period of time required for the protocol and the full-availability and low price of the required reagents place PK+HID as a cheap, simple and fast alternative to traditional RNA extraction protocols.

Indeed, the use of this method was authorized by the National Administration of Drugs, Foods and Medical Devices of Argentina (ANMAT) on December 14, 2020 and it is currently being used in some public and private diagnostic labs.

Supporting information

S1 Fig. PK+HID can be combined with a detection kit based on LAMP reactions.

Representative colorimetric reactions run using COVID-19 Neokit Tecnoami, based on LAMP. Diagnostics is made according to final color of the solution: blue, positive; violet, negative.

https://doi.org/10.1371/journal.pone.0247792.s001

(PDF)

S2 Fig. Effect of saline solution dilution in PK+HID method.

Five positive nasopharyngeal swab samples (#1 to #5) were processed by adding 10 μl of proteinase K 10mg/ml (PK+HID samples) or 10 μl in proteinase K buffer (HID’ samples) and subjected to thermal incubations (55°C for 15 min and 98°C for 5 min). The viral N1 and N2 genes and the human RNase P gene (RP) were amplified and detected by RT-qPCR. (a) CT values obtained from RT-qPCR analysis of the same samples prepared by both different methods. (b) Ratio between relative amplicon amounts (n) of PK+HID and HID’ samples. The median of each measurement is represented with a line in the bars and the lengths of these bars represent the standard error. (c) Amplification efficiencies (EPCR). The median of each measurement is represented with a line in the bars and the lengths of these bars represent the standard error.

https://doi.org/10.1371/journal.pone.0247792.s002

(PDF)

S3 Fig. Performances of HID and PK+HID methods targeting N and ORF1ab amplicons.

Twelve positive nasopharyngeal swab samples (#1 to #12) were processed by heat inactivation (HID samples), proteinase K treatment followed by heat inactivation (PK+HID samples) or RNA extraction (purified RNA samples). The viral N and ORF1ab genes and the human RNase P gene (RP) were amplified and detected by RT-qPCR. CT values obtained from RT-qPCR analysis of the same samples prepared by the three different methods.

https://doi.org/10.1371/journal.pone.0247792.s003

(PDF)

S4 Fig. Stability of PK+HID samples at 4°C.

Ten positive nasopharyngeal swab samples were processed by PK treatment followed by heat inactivation. These PK+HID samples were kept at 4°C. The viral N and ORF1ab genes and the human RNase P gene (RP) were amplified and detected by RT-qPCR at different times of incubation at 4°C (0, 20 and 40 h). Differences between CT values obtained at different times and time = 0 h. Thin, light lines represent individual sample measurements. Strong lines represent mean + SEM values.

https://doi.org/10.1371/journal.pone.0247792.s004

(PDF)

S5 Fig. Performance of PK+HID method targeting N and ORF1ab amplicons.

Randomly-selected nasopharyngeal swab samples were simultaneously processed for RNA extraction (purified RNA) or subjected to treatment with PK 1 mg/ml (55°C for 15 min) followed by heat-inactivation at 98°C for 5 min (PK+HID). The viral N and ORF1ab genes and the human RP gene were amplified and detected by RT-qPCR. (a) CT values obtained for the same samples prepared by the two different methods; the line obtained by regression of the data (continuous lines). RP plot also includes the data obtained in SARS-CoV-2 negative samples whereas N and ORF1ab plots include those data of negative samples that provided CT values above the positive-negative cut-off. (b) Ratio between the relative amplicon copy amounts detected using PK+HID or RNA purification for those nasopharyngeal swab samples that presented CT values >30 for any of the viral genes in purified RNA samples (Ndata = 10 and 12 for N and ORF1ab genes, respectively).

https://doi.org/10.1371/journal.pone.0247792.s005

(PDF)

S6 Fig.

https://doi.org/10.1371/journal.pone.0247792.s006

(PNG)

S1 Table. Mean variations in CT values in HID and PK+HID samples compared to purified RNA samples.

Values are expressed as mean ± standard error.

https://doi.org/10.1371/journal.pone.0247792.s007

(PDF)

S2 Table. CT values obtained in RT-qPCR analysis of the same HID+PK samples introducing 5 or 10 μl of sample volume as input for RT-qPCR.

https://doi.org/10.1371/journal.pone.0247792.s008

(PDF)

S3 Table. CT values of false-negative samples detected with the PK+HID method, using the results obtained with purified RNA samples as reference.

https://doi.org/10.1371/journal.pone.0247792.s009

(PDF)

S4 Table. CT values of the positive samples processed by RNA extraction and analyzed by RT-qPCR targeting N and ORF1ab viral genes and human RP gene, and diagnostic results obtained with the same samples processed by RNA extraction or PK+HID and analyzed using a LAMP detection kit.

https://doi.org/10.1371/journal.pone.0247792.s010

(PDF)

Acknowledgments

V. L. thanks J. C. Reboreda and A. Quaglino for their support.

Citation: Genoud V, Stortz M, Waisman A, Berardino BG, Verneri P, Dansey V, et al. (2021) Extraction-free protocol combining proteinase K and heat inactivation for detection of SARS-CoV-2 by RT-qPCR. PLoS ONE 16(2): e0247792. https://doi.org/10.1371/journal.pone.0247792

References

1. Kucharski AJ, Klepac P, Conlan AJK, Kissler SM, Tang ML, Fry H, et al. Effectiveness of isolation, testing, contact tracing, and physical distancing on reducing transmission of SARS-CoV-2 in different settings: a mathematical modelling study. Lancet Infect Dis. 2020. Epub 2020/06/20. pmid:32559451.

2. Guglielmi G. The explosion of new coronavirus tests that could help to end the pandemic. Nature. 2020;583(7817):506–9. Epub 2020/07/19. pmid:32681157.

3. Esbin MN, Whitney ON, Chong S, Maurer A, Darzacq X, Tjian R. Overcoming the bottleneck to widespread testing: a rapid review of nucleic acid testing approaches for COVID-19 detection. Rna. 2020;26(7):771–83. Epub 2020/05/03. pmid:32358057

4. Akst J. RNA Extraction Kits for COVID-19 Tests Are in Short Supply in US. The Scientist Magazine. 2020 Mar 11, 2020.

5. Slabodkin G. FDA chief warns of supply ’pressure’ on reagents for coronavirus tests. MedTech Dive Magazine. 2020 March 12, 2020.

6. Brown J RA L., Shah D H K. Validation of an extraction-free RT-PCR protocol for detection of SARS-CoV2 RNA. medRxiv (Posted as a preprint). 2020. Epub May 01, 2020.

7. Sahajpal NS M A. K., A , Ananth S, Kothandaraman A, Hedge M, Chaubey A, et al. Clinical validation of innovative, low cost, kit-free, RNA processing protocol for RT-PCR based COVID-19 testing. medRxiv (Posted as a preprint). 2020. Epub Jul 30, 2020.

8. Lubke N, Senff T, Scherger S, Hauka S, Andree M, Adams O, et al. Extraction-free SARS-CoV-2 detection by rapid RT-qPCR universal for all primary respiratory materials. J Clin Virol. 2020;130:104579. Epub 2020/08/17. pmid:32795959

9. Barza R, Patel P, Sabatini L, Singh K. Use of a simplified sample processing step without RNA extraction for direct SARS-CoV-2 RT-PCR detection. J Clin Virol. 2020;132:104587. Epub 2020/09/09. pmid:32898817.

10. Ragan KB, Bhadra S, Choi JH, Towers D, Sullivan CS, Ellington AD. Comparison of media and standards for SARS-CoV-2 RT-qPCR without prior RNA preparation. medRxiv. 2020. Epub 2020/08/15. pmid:32793925

11. Kriegova E, Fillerova R, Kvapil P. Direct-RT-qPCR Detection of SARS-CoV-2 without RNA Extraction as Part of a COVID-19 Testing Strategy: From Sample to Result in One Hour. Diagnostics (Basel). 2020;10(8). Epub 2020/08/23. pmid:32824767

12. Israeli O, Beth-Din A, Paran N, Stein D, Lazar S, Weiss S, et al. Evaluating the efficacy of RT-qPCR SARS-CoV-2 direct approaches in comparison to RNA extraction. Int J Infect Dis. 2020;99:352–4. Epub 2020/08/14. pmid:32791207

13. Ulloa S, Bravo C, Parra B, Ramirez E, Acevedo A, Fasce R, et al. A simple method for SARS-CoV-2 detection by rRT-PCR without the use of a commercial RNA extraction kit. J Virol Methods. 2020;285:113960. Epub 2020/08/25. pmid:32835738

14. Adams NM, Leelawong M, Benton A, Quinn C, Haselton FR, Schmitz JE. COVID-19 diagnostics for resource-limited settings: Evaluation of "unextracted" qRT-PCR. J Med Virol. 2020. Epub 2020/08/12. pmid:32779772

15. Mancini F, Barbanti F, Scaturro M, Errico G, Iacobino A, Bella A, et al. Laboratory management for SARS-CoV-2 detection: a user-friendly combination of the heat treatment approach and rt-Real-time PCR testing. Emerg Microbes Infect. 2020;9(1):1393–6. Epub 2020/06/20. pmid:32552549.

16. Hasan MR, Mirza F, Al-Hail H, Sundararaju S, Xaba T, Iqbal M, et al. Detection of SARS-CoV-2 RNA by direct RT-qPCR on nasopharyngeal specimens without extraction of viral RNA. PLoS One. 2020;15(7):e0236564. Epub 2020/07/25. pmid:32706827

17. Smyrlaki I, Ekman M, Lentini A, Rufino de Sousa N, Papanicolaou N, Vondracek M, et al. Massive and rapid COVID-19 testing is feasible by extraction-free SARS-CoV-2 RT-PCR. Nat Commun. 2020;11(1):4812. pmid:32968075

18. Prevention CfDCa. CDC 2019-Novel Coronavirus (2019-nCoV) Real-Time RT-PCR Diagnostic Panel In: Diseases DoV, editor. Atlanta, GA2020.

19. Sung H, Yong D, Ki CS, Kim JS, Seong MW, Lee H, et al. Comparative Evaluation of Three Homogenization Methods for Isolating Middle East Respiratory Syndrome Coronavirus Nucleic Acids From Sputum Samples for Real-Time Reverse Transcription PCR. Ann Lab Med. 2016;36(5):457–62. Epub 2016/07/05. pmid:27374711

20. Espy MJ, Patel R, Paya CV, Smith TF. Comparison of three methods for extraction of viral nucleic acids from blood cultures. J Clin Microbiol. 1995;33(1):41–4. Epub 1995/01/01. pmid:7699063

21. Mallmann L, Schallenberger K, Demolliner M, Eisen AKA, Hermann BS, Heldt FH, et al. Pre-treatment of the clinical sample with Proteinase K allows detection of SARS-CoV-2 in the absence of RNA extraction. bioRxiv (Posted as a preprint). 2020. Epub May 09, 2020.

22. Marzinotto S, Mio C, Cifu A, Verardo R, Pipan C, Schneider C, et al. A streamlined approach to rapidly detect SARS-CoV-2 infection, avoiding RNA extraction. medRxiv (Posted as a preprint). 2020. Epub Apr 11, 2020. pmid:33343767

23. Ruijter JM, Ramakers C, Hoogaars WM, Karlen Y, Bakker O, van den Hoff MJ, et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 2009;37(6):e45. Epub 2009/02/25. pmid:19237396

24. Vogels CBF, Brito AF, Wyllie AL, Fauver JR, Ott IM, Kalinich CC, et al. Analytical sensitivity and efficiency comparisons of SARS-CoV-2 RT-qPCR primer-probe sets. Nat Microbiol. 2020. Epub 2020/07/12. pmid:32651556.

25. Akobeng AK. Understanding diagnostic tests 1: sensitivity, specificity and predictive values. Acta Paediatr. 2007;96(3):338–41. Epub 2007/04/05. pmid:17407452.

26. Alcoba-Florez J, Gonzalez-Montelongo R, Inigo-Campos A, de Artola DG, Gil-Campesino H, The Microbiology Technical Support T, et al. Fast SARS-CoV-2 detection by RT-qPCR in preheated nasopharyngeal swab samples. Int J Infect Dis. 2020;97:66–8. Epub 2020/06/04. pmid:32492531.

27. Bruce EA, Huang ML, Perchetti GA, Tighe S, Laaguiby P, Hoffman JJ, et al. Direct RT-qPCR detection of SARS-CoV-2 RNA from patient nasopharyngeal swabs without an RNA extraction step. PLoS Biol. 2020;18(10):e3000896. pmid:33006983 competing interests. While DJS was employed at IXIS LLC at the time of this study, his employment there did not create a competing interest. Further, IXIS LLC had no involvement in this study.

28. Fomsgaard AS, Rosenstierne MW. An alternative workflow for molecular detection of SARS-CoV-2—escape from the NA extraction kit-shortage, Copenhagen, Denmark, March 2020. Euro Surveill. 2020;25(14). Epub 2020/04/16. pmid:32290902

29. Ladha A J J., Abudayyeh O, Gootenberg J, Zhang F. A 5-min RNA preparation method for COVID-19 detection with RT-qPCR. medRxiv (Posted as a preprint). 2020. Epub May 08, 2020. pmid:20055947.

30. Chu AW, Chan WM, Ip JD, Yip CC, Chan JF, Yuen KY, et al. Evaluation of simple nucleic acid extraction methods for the detection of SARS-CoV-2 in nasopharyngeal and saliva specimens during global shortage of extraction kits. J Clin Virol. 2020;129:104519. Epub 2020/07/07. pmid:32629187

31. Michel D, Danzer KM, Gross R, Conzelmann C, Muller JA, Freischmidt A, et al. Rapid, convenient and efficient kit-independent detection of SARS-CoV-2 RNA. J Virol Methods. 2020:113965. Epub 2020/09/07. pmid:32891677

32. Merindol N, Pepin G, Marchand C, Rheault M, Peterson C, Poirier A, et al. SARS-CoV-2 detection by direct rRT-PCR without RNA extraction. J Clin Virol. 2020;128:104423. Epub 2020/05/18. pmid:32416598