Xander F.
van Kooten‡
,
Yana
Rozevsky‡
,
Yulia
Marom
,
Efrat
Ben Sadeh
and
Amit
Meller§
*
Department of Biomedical Engineering, Technion-IIT, Haifa 32000, Israel. E-mail: ameller@technion.ac.il
First published on 7th March 2022
The majority of RNA based COVID-19 diagnostics employ enzymatic amplification to achieve high sensitivity, but this relies on arbitrary thresholding, which complicates the comparison of test results and may lead to false outcomes. Here we introduce solid-state nanopore sensing for label-free quantification of SARS-CoV-2 RNA in clinical nasal swab samples. This PCR-free method involves reverse transcribing a target gene on the viral RNA before enzymatically digesting all but the resulting dsDNA. Ratiometric quantification of RNA abundance is achieved by single-molecule counting and length-based nanopore identification of dsDNA from a SARS-CoV-2 gene and a human reference gene. We graded nasal swab samples from >15 subjects and find that the SARS-CoV-2 ratiometric nanopore index correlates well with the reported RT-qPCR threshold cycle for positive classified samples. Remarkably, nanopore analysis also reports quantitative positive outcomes for clinical samples classified as negative by RT-qPCR, suggesting that the method may be used to diagnose COVID-19 in samples that may evade detection. We show that the sample preparation workflow can be implemented using a compact microfluidic device with integrated thermal control for semi-automated processing of extremely small sample volumes, offering a viable route towards automated, fast and affordable RNA quantification in a small and portable device.
Single-molecule sensing can overcome some of the limitations of amplification-based detection, provided that it offers sufficient specificity and sensitivity, and preferably avoids non-linear and error-prone amplification steps.16–18 Nanopores are an emerging class of single-molecule biosensors, developed primarily for single-molecule DNA sequencing. Owing to their high sensitivity, in recent years nanopore biosensors have been adapted for sensing clinical biomarkers in biofluids, such as circulating DNAs and even proteins.19–22 Specifically, solid-state nanopores (ssNPs) can be used to quantify specific RNA and mRNA molecules with high accuracy and sensitivity by converting them to double-stranded DNA using a lossless biochemical assay.23 Reverse transcription quantitative nanopore sensing (RT-qNP) allows direct quantification of gene expression, bypassing the need for PCR amplification.23 To date, however, the use of ssNP biosensors for quantitative molecular detection of SARS-CoV-2 in clinical samples has not been reported. Here we present a proof-of-principle study to assess the ability of ssNPs to quantify multiple RNA types directly extracted from clinical patient samples. We show that the nanopore can distinguish between viral RNA and a human reference gene, thus allowing direct counting of the two genes in a rapid and label-free molecular test. Importantly, the ratio of the two RNA events can serve as a quantitative measure for the viral load in the sample, showing good correlation with the Ct value obtained by RT-qPCR run in a certified diagnostic lab. Remarkably, our method reports quantitative RNA loads in clinical samples classified as negative SARS-CoV-2 by RT-qPCR, in some cases indicating a clear positive result.
The purely additive library preparation of the sensing method described here lends itself for liquid handling in a microfluidic device, as it does not involve any centrifugation or purification steps. This may open the possibility for fully automated and fast sample processing and subsequent nanopore analysis. To demonstrate this, we introduce a design for a device capable of preparing clinical sample volumes for purely electrical analysis by ssNPs, hence permitting future development of a fully portable compact sensing technology. The fluidic device performs all the necessary functions for reverse transcription of RNA and digestion of the undesired analytes, using pressure-driven flow with capillary synchronization valves in a disposable lab-on-chip format. We use this device to process both analytical samples containing SARS-CoV-2 RNA and clinical samples from confirmed COVID-19 patients, and show that nanopore quantification of these samples yields results in line with those processed in a standard well-plate format.
In this study we count the ratio of dsDNA molecules reverse transcribed from two genes: SARS-CoV-2 RdRP, and the human housekeeping gene RPP30. The two dsDNA fragments are produced in the same reaction by including specific reverse primers for the second strand synthesis, resulting in 107 bp and 758 bp long molecules for the two genes, respectively. This design allows straightforward discrimination between the two strands using ssNPs ranging from 3.5 nm to roughly 6 nm in diameter. Fig. 1e shows typical dsDNA translocations corresponding to the two target genes reverse transcribed from synthetic RNAs sources (see ESI Fig. S1 and S2†). The event diagram heat map shows the density of single-molecule events, each represented by its IB and tD values. Notably, in this method validation step, results were obtained in two consecutive experiments using the same nanopore device (G = 9.33 nS, corresponding to a ∼4 nm pore). The device was washed thoroughly with clean buffer in between measurement until no events were observed for at least 10 minutes. The two dsDNA lengths produced two distinct clusters of translocation events, where the longer RPP30 produced longer and deeper amplitude events, as expected. Double-Gaussian fitting of the histograms yields a clear separation of peaks for IB (Fig. 1f), but less so for tD, due to the broader distribution of events (Fig. 1g). A good separation of populations could nevertheless be obtained by using a Gaussian Mixture Model, which uses two-dimensional information to cluster events. Importantly, the ability to serially perform independent analyses of the two synthetic and pure dsDNA components using the same ssNP, establishes the ability of the nanopore to quantify multiple components for a ratiometric approach.
The events capture rate (R), which is proportional the dsDNA concentration in the cis chamber, depends on the size and shape of the nanopore used,25 as well as on the efficiency of the sample collection and preparation. As such, the events rate of RdRP alone is an unsuitable metric for comparing between samples analyzed using different nanopore devices. Instead, the reference gene, which acts as an internal calibrator for the event rate, should be used to calculate a normalized metric. Accordingly, we define the SARS-CoV-2 ratiometric nanopore index as:
![]() | (1) |
Fig. 2b displays the analysis of 4 clinical samples, using 4 different nanopore devices. We start by analyzing a pre-COVID-19 nasal swab taken before the outbreak of the pandemic, serving as a true negative clinical sample (S0). As expected, despite the presence of both primer sets in the reaction, this sample displays a single population of events with relatively low IB value of 0.4, which we attribute to the human RPP30 gene. We note that the other reaction byproducts, such as digested ssDNA or digested proteins, do not produce additional signals. The measurement was repeated using a different nanopore device, yielding the same result (see ESI Fig. S10†). In contrast, samples S1, S2 and S3, with decreasing SARS-CoV-2 RNA abundance as indicated by their RT-qPCR cycle threshold (Ct) values, exhibit two distinct populations of events. We ascribe long, low-amplitude events to RdRP (labelled in red) and longer, deeper events to RPP30. The upper part of each panel displays typical sets of events, the middle part displays the 2D GMM clustering of the events and the lower part shows the arrival time histogram for both populations. We find that while the event distributions may vary among different nanopore devices, the GMM classification is able to differentiate between the two populations in each experiment. This highlights the importance of co-analyzing the reference gene as an internal control for sample integrity and for nanopore benchmarking.
About a dozen additional clinical samples were analyzed in a similar fashion. We find that in nine of the samples analyzed, as the Ct value increased from 16 to nearly 34, the corresponding ratiometric nanopore index (eqn (1)) decreased from 0.61 to 0.42 (Fig. 2d, red circles). Histograms of IB and tD for all samples are provided in ESI Fig. S3–S6.† We find that the nanopore results and the relative RNA expression levels measured by single-molecule counting can be related by a non-linear Hill function: with the following parameters: b = 0.43 ± 0.01; m = 0.611 ± 0.005; r = −8.37 ± 1.57; x1/2 = 25.6 ± 0.7. Here x1/2 represents the RT-qPCR cycle number at which the nanopore identifies equal abundance of the viral gene (RdRP) and the reference human gene (RPP30). The strong non-linear factor r is likely related to the fact that we compare an exponential PCR amplification process to unamplified molecule-by-molecule counting, which is strictly linear. Notably, at both extremes, namely at very high viral load (Ct < 15 cycles) and very low viral load (Ct > 35), the results appear to be less dependent on the number of cycles. This is a consequence of the RT-qPCR process and is likely not due to the nanopore sensitivity limits, as these limit values correspond to RNP values of 0.75 and 1.57 for the base (b) and maximum (m) values, respectively. In a typical translocation batch, one would then have to misclassify roughly 25% to 34% events, which is highly unlikely in these experiments.24
Next, we analyzed six additional samples that were classified as negative, as they produced no detectable signal in RT-qPCR after 40 cycles. The analysis of each sample was repeated at least twice using independent ssNP devices (ESI Fig. S3, S5 and S6†). Four out of these six samples yielded a single population of events, and were therefore assigned a value of RNP = 0 (Fig. 2d, green circles). However, two other samples (ESI Fig. S3, S5,† and Fig. 2c) returned a positive outcome in three repeats despite being classified as negative by RT-qPCR. These samples are marked as blue circles in Fig. 2d. We estimated the expected Ct values for each false negative case by the Hill function, corresponding to moderate SARS-CoV-2 abundance (Ct = 23–25). The data used for the quantitative analysis is summarized in Table 1 and ESI Fig. S3–S5.†
Sample | C t value | R RPP30 (s−1) | R RdRP (s−1) | R NP |
---|---|---|---|---|
S0 | — | 1.27 ± 0.02 | 0 | 0 |
S1 | 16.03 | 0.36 ± 0.02 | 0.57 ± 0.02 | 0.61 ± 0.02 |
S2 | 27.22 | 0.26 ± 0.01 | 0.25 ± 0.01 | 0.49 ± 0.01 |
S3 | 33.62 | 0.87 ± 0.02 | 0.63 ± 0.03 | 0.42 ± 0.04 |
S4 | N.D. | 1.60 ± 0.04 | 1.75 ± 0.05 | 0.52 ± 0.06 |
S5 | N.D. | 0.12 ± 0.01 | 0.16 ± 0.01 | 0.57 ± 0.01 |
S6 | 24.86 | 0.15 ± 0.01 | 0.18 ± 0.01 | 0.54 ± 0.01 |
S7 | N.D. | 5.25 ± 0.06 | 0 | 0 |
S8 | N.D. | 0.35 ± 0.03 | 0 | 0 |
S9 | N.D. | 0.74 ± 0.03 | 0 | 0 |
S10 | N.D. | 0.48 ± 0.01 | 0 | 0 |
S11 | 20.28 | 0.81 ± 0.04 | 1.13 ± 0.04 | 0.58 ± 0.04 |
S12 | 31.77 | 0.37 ± 0.01 | 0.30 ± 0.02 | 0.45 ± 0.02 |
S13 | 14.3 | 0.07 ± 0.01 | 0.11 ± 0.01 | 0.61 ± 0.01 |
S14 | 13.89 | 0.29 ± 0.03 | 0.46 ± 0.03 | 0.61 ± 0.03 |
S15 | 22.07 | 0.14 ± 0.01 | 0.19± 0.01 | 0.57 ± 0.01 |
To ensure controlled mixing with minimal dispersion of incubated sample, the liquids are sequentially introduced into an initially empty (air-filled) hydrophobic channel. However, this approach is likely to trap air bubbles,28 which expand dramatically at elevated temperatures and can push liquids back to the reservoirs. To overcome this, we implemented capillary valves at each junction, which synchronize interfaces arriving at different times by retaining them until the applied pressure exceeds the Laplace (‘burst’) pressure, which is defined by the geometry.28–30 The operation of a capillary synchronization junction is shown in Fig. 3d, where a filling pressure Pf is initially applied to the reservoir containing the RNA sample and RT mix. Once both interfaces have reached the junction (ii), a pressure Pb, which exceeds the Laplace pressure, is applied, causing the valve to burst and the flows to merge (iii). Apart from preventing bubbles and allowing temperature-sensitive reagents to be stored away from the ‘hot zone’ of the chip, capillary valves offer a way towards full automation, by presenting incoming flows with a ‘soft’ block that can withstand an applied pressure. A similar constriction geometry is added as a retention valve at the end of every incubation chamber, to pin the fluid front and prevent capillary creeping during incubation.
To demonstrate the performance of the microfluidic device for nanopore quantification of viral samples, we compared ssNP sensing of an additional clinical sample (S6, Ct = 25) after sample preparation in the microfluidic device and in a vial. Fig. 3 and ESI Fig. S7–S9† show the nanopore results for samples processed using the conventional and microfluidic workflows (Fig. 3e and f, respectively). Remarkably, despite being run on pores with a different conductance, the microfluidic and conventional workflow yield nearly identical relative CoV-2 rates (RNP = 0.55 ± 0.01 and 0.54 ± 0.01, respectively). This demonstrates the feasibility of microfluidic sample processing and highlights the sensitivity and robustness of the nanopore-based workflow.
A fully autonomous device would further require integration of an in-line nanopore sensor to eliminate the manual sample transfer,33 and should support pre-loading of reagents in e.g. thermoelectrically cooled reservoirs. Furthermore, an investigation of reaction efficiency in each stage (currently totaling 60 min) is expected to show that incubation times can be substantially reduced, making the microfluidic workflow suitable for fast-turnaround testing scenarios.
Our studies open an avenue towards a lab-on-chip device for biomedical and clinical needs, demonstrated in the context of COVID-19, but potentially also relevant to other pathogens as well as cancer diagnostics based on RNA biomarkers. Here, we illustrated the strength of the technique by diagnosing a range of SARS-CoV-2 positive and negative clinical samples. In some cases, our results point to false negative RT-qPCR classification (Fig. 2 and ESI Fig. S3, S5†). While the disparity between nanopore and RT-qPCR results in these cases remains to be investigated, false negative RT-qPCR outcomes are known to occur, perhaps most notably in samples taken during disease progression.10,11
We believe that the accuracy of the nanopore-based test can be largely attributed to the fact that both the viral and human gene are co-processed and co-analyzed, and the results are always presented as a ratio of the two. This cancels out sample-to-sample variations due to e.g. sample source concentration, enzymatic process efficiency and sampling errors. Furthermore, the lack of amplification eliminates any potential amplification bias or error. In the nanopore test, true negative samples (pre-COVID-19 nasal swabs) clearly display a single event population corresponding to the human reference gene RPP30 (see Fig. 2 and ESI Fig. S6†). The positive signal of RPP30 provides a powerful confirmation that the clinical sample was acquired effectively. In contrast, the appearance of two event populations in samples with significant relative counts is strong evidence of the presence of the viral gene. This translates to high nanopore accuracy, as shown in Fig. 2 and ESI Fig. S3–S5.†
The workflow shown here is expected to be compatible with various RNA sampling approaches, including direct detection from saliva without any extraction31 or even from sewage.32 The combination of a microfluidic processing device with nanopore sensing has an important impact towards establishing rapid lab-on-chip devices for medical diagnostics in general.
Extracted full NA clinical SARS-CoV-2 samples were provided by Rambam Medical Center (Haifa, Israel). All samples were approved by Rambam Health Corporation Ethics Committee (Haifa, Israel) and were diagnosed as positive or negative, based on RT-qPCR results (Ct values). Nasal swabs taken both from nostrils and throat were collected by healthcare providers. The samples were mixed with lysis buffer and then whole nucleic acids (NA) were extracted using automated NA instruments, STARMag 96 × 4 Universal Cartridge kit’ (Seegene Inc., South Korea) or magLEAD (Precision System Science). RT-qPCR was performed on the E, RdRP, and N genes of SARS-CoV-2 (Sarbeco probes).3,34 Reactions with 8 μl of samples were heated to 50 °C for 30 minutes for reverse transcription, denatured in 95 °C for 10 minutes, and then 40 cycles of amplification were carried out at 95 °C for 15 seconds and 55 °C for 32 seconds. Fluorescence was measured using the FAM parameters. As a result, the cycle threshold (Ct) was determined except for the undetermined cases in which no signal was observed for 40 PCR cycles or more (classified as “negative”). Information on Ct values is presented in Table 1. The clinical samples were divided into 5 μl aliquots and stored at −80 °C.
For the Nanopore-based analysis 5 μl aliquots of the same NA extracted samples used for the RT-qPCR, were diluted to final volume of 10 μl, with the following buffer: 20 mM Tris–HCl 10 mM (NH4)2SO4, 10 mM KCl, 0.1% Triton® X-100 (pH 8.8 @ 25 °C), 6 mM MgSO4, 1.4 mM dNTPs, 0.2 mM of each primer, 3 U WarmStart® RTx (NEB) and 6 U Bst 2.0 WarmStart™ DNA polymerase (NEB). No loop primers were used in the reaction to avoid amplification. The reaction was carried out at 62 °C for 10 min. A positive control sample composed of 0.5 × 106 copies synthetic SARS-CoV-2 and 50 ng of extracted total RNA from HCT116 cells, was run in parallel. Afterwards, cDNA samples were treated with 20 U of Exonuclease I (NEB) at 37 °C for 15 min, followed by 2.5 U RNase I (Thermo Fisher Scientific) at 50 °C for 20 min, and finally 0.2 U of ProK (NEB). The sample was diluted 40-fold in nanopore buffer.
Clinical sample processing for the microfluidic workflow: the clinical samples were prepared using the same concentrations as detailed above, but only 1 μl of each sample was introduced into the device. After processing the samples with the fluidic device, the sample was diluted 20-fold in nanopore buffer.
Chips were cleaved from the wafer and cleaned using a 2:
1 solution of H2SO4/H2O2. They were then glued to a Teflon insert using Ecoflex 5 (Smooth-ON) and immersed in buffer (1 M KCl, 40 mM Tris–HCl, 1 mM EDTA). An Ag/AgCl electrode was placed in both the cis and trans reservoir and nanopores were drilled using dielectric breakdown (DB).36 The locally thinned SiNx ensured that the pore was formed at the center of the membrane. DB voltage pulses with an amplitude of 8–9 V and a duration of 225 ms were applied at 1 s intervals, and the ionic current was probed at 300 mV between pulses. Once the current exceeded 0.4 nA, the pulse amplitude was reduced to 1–3 V until the desired nanopore conductance was reached.
Afterward, the nanopores were kept under a low probing voltage (0.15 to 0.3 V) in a buffer solution (1 M KCl, 40 mM Tris–HCl, 1 mM EDTA, pH 7.5) to obtain a stable open-pore current. During the experiment, translocation events were monitored using an Axon 200B amplifier, filtered at 100 kHz, and acquired using a custom LabVIEW software (National Instruments). After collecting the data, we performed offline analysis using a custom LabVIEW program to extract the dwell time (tD), current blockage (IB), and arrival time (ta) of each translocation event according to an electrical threshold. Collisional events (short current spikes due to unsuccessful translocations) could be filtered by setting a current threshold in the analysis of IB < 0.8.
Optical monitoring of the nanopores was performed as follows: 500 mM CaCl2 was added to the Cis chamber, 500 nM of RHOD-2 (Molecular Probes®) and 10 mM EGTA (Ca2+ chelator) were added to the trans chamber. When negative potential bias is applied, Ca2+ ions driven to the pore and binds to RHOD-2, releasing strong fluorescence signal in the pore location. The signal is radially declines as Ca2+ ions bind to EGTA. A custom wide field microscope with a collimated 560 nm laser (iFlex-Viper, PointSource) coupled through a single-mode optical fiber is focused on a high NA objective (Olympus Plan Apochromat 60×/1.45) and expanded to illuminate the membrane position at the image plane. The emitted light is collected by the same objective and focused onto an EMCCD camera (Andor iXon 887).36
The PDMS/glass device was placed on the copper heat diffuser of two thermoelectric heaters (Melcor CPI-4-31-045L) wired in series and controlled in a closed loop (3040, Newport). The cold side was bonded with thermal grease to a heat sink which was cooled using a chiller (PolyScience). After pipetting the reagents into the reservoirs, we connected the pressure lines. The pressure in the reservoirs was set by a controller (MFCS-EZ, Fluigent). Each reservoir was connected to two parallel on/off valves (075-T2NC12-32 M, Bio-Chem) acting as a two-way selector valve. The valves were controlled by a 6-channel valve controller (VC-6, Warner Instrument Corp.). For visualization purposes, the device was imaged with an upright stereomicroscope (AZ100, Nikon).
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1nr07787b |
‡ Equally contributing authors. |
§ Current address: The WYSS Institute, Harvard University, Boston, MA 02115 USA. |
This journal is © The Royal Society of Chemistry 2022 |