Abstract
High Ct-values falling in the grey zone are frequently encountered in SARS-CoV-2 detection by real-time reverse transcription PCR (rRT-PCR) and have brought urgent challenges in diagnosis of samples with low viral load. Based on the single-stranded DNA reporter trans-cleavage activity by Cas12a upon target DNA recognition, we create a Specific Enhancer for detection of PCR-amplified Nucleic Acids (SENA) to confirm SARS-CoV-2 detection through specifically targeting its rRT-PCR amplicons. SENA is highly sensitive, with its limit of detection being at least 2 copies/reaction lower than that of the corresponding rRT-PCR, and highly specific, which identifies both false-negative and false-positive cases in clinic applications. SENA provides effective confirmation for nucleic acid amplification-based molecular diagnosis, and may immediately eliminate the uncertainty problems of rRT-PCR in SARS-CoV-2 clinic detection.
One Sentence Summary CRISPR-Cas12a-based COVID-19 diagnosis.
Main Text
Since December 2019, the outbreak of COVID-19, caused by the infection of the SARS-CoV-2, has rapidly spread throughout the world (1), and is now a global pandemic. Till June 1, 2020, the outbreak has affected 216 countries, areas and territories, infected 6 million people, and caused more than 370 thousand of death (2). One of the greatest public health concerns in combating against the pandemic is a prompt response to meet the urgent demand for rapid and accurate diagnosis of the virus. Currently, nucleic acid amplification-based molecular diagnostics (MDx) is the most accurate, fast and affordable and thus the preferred method for diagnosis of SARS-CoV-2 infection, and the real-time reverse transcription PCR (rRT-PCR) kits have been successfully developed by quite a few laboratories and commercial companies (3). However, since its clinical application at least five months ago in China, the diagnostic performance of rRT-PCR for SARS-CoV-2 has brought some urgent challenges, particularly with the frequently encountered high Ct-value designated “grey zone” associated uncertain negative or positive readouts (4-8). Besides of “human error” factors such as misconducted sampling, unqualified reagents and uncalibrated diagnostic equipment, inefficient RT reaction and PCR amplification of clinical samples with very low viral loads are likely the major intrinsic causative factors for the fuzzy rRT-PCR readouts and uncertain diagnosis. Although repetitive sampling and assays are implemented for final confirmation of the diagnosis, these “trouble shooting” efforts are time-consuming and may still fail to detect the low viral load samples from some mild or asymptomatic patients, or from the recovering patients, resulting in false-negative diagnosis that may cause serious public concerns in battling against the pandemic (Figure 1).
With the characterization of non-specific trans-cleavage activities against single-stranded nucleic acids in several CRISPR-associated (Cas) proteins, e.g., Cas12, Cas13 and Cas14 (9-15), Clustered Regularly Interspaced Short Palindromic Repeats Diagnostics (CRISPR-Dx) technology (16, 17) was established and has been developing rapidly. The underline mechanism for CRISPR-Dx, as illustrated by the Cas12a-based HOLMES system for example (18), is based on the efficient trans-cleavage activity against a fluorophore quencher (FQ)-labeled single-stranded DNA reporter by Cas12a triggered upon target DNA recognition, which is guided by a specific CRISPR RNA (crRNA), generating exponentially increasing fluorescence signal within several minutes. With this mechanism, here, we design a Specific Enhancer for detection of PCR-amplified Nucleic Acids (SENA) to improve both the detection sensitivity and specificity against the pre-amplified targeted SARS-CoV-2 genomic fragments (Figure 1). Briefly, the COVID-19 clinic samples are firstly analyzed by the well-established rRT-PCR assays and amplicons with uncertain readouts are then verified by SENA in a physically isolated space, avoiding contamination of the PCR laboratory during pipetting.
To prepare appropriate crRNAs for SENA detection, we firstly determined the amplicon sequences from several commercial rRT-PCR kits used in China and then designed specific crRNAs corresponding to each of the distinct amplicons (Table S1). Candidate crRNAs were prepared and analyzed in a SENA system comprised of Cas12a, crRNAs, FQ-reporter and the rRT-PCR products using templates of either the positive or negative controls, and the best crRNAs, i.e., the lowest fluorescence with the negative control and highest with the positive control, were chosen for further SENA assays (data not shown). If a kit detects more than one gene, e.g., Orf1ab (“abbreviated as O”), E and N genes, corresponding crRNAs were mixed in a SENA system to enhance the detection signal.
The performance of SENA was quantitatively characterized via a systematic titration upon rRT-PCR amplicons employing pure SARS-CoV-2 RNA standards comprised of the O and N fragments as the templates. As it is aware that the viral nucleic acids extracted from patients’ samples such as nasopharyngeal swabs usually contain some biological and chemical contaminants that might inhibit the enzyme activities for reverse transcription and PCR reactions (19) and is likely one of the causal effects attributed to the low efficiency of rRT-PCR in clinical analysis. In order to mimic the clinical sampling for the titration experimentation, the RNA standards were serially diluted in buffer prepared by mixing the nucleic acid extracts from 40 COVID-19 negative people, generating RNA templates ranging from 0.025 to 25 copies per reaction (Rx).
Due to the Poisson distribution property of sampling, replica variations become extremely significant when the template copies in individual reaction are designed to be low, i.e., less than 3-4 copies/Rx, near the limit of detection (LoD) for rRT-PCR (20, 21), and extremely low, i.e., equal to and less than 1 copy/Rx. To overcome this sampling ambiguity problem, we performed 9 replicas for groups with 1 and 0.5 RNA template copies/Rx while 6 replicas for each of the rest concentrations. In addition, although the rRT-PCR assay supplier, Shanghai BioGerm (BJ), who follows the Chinese CDC recommended primer sets (Table S1), recommends 40 cycles of PCR amplification, we set 45 cycles as routine aiming at recording maximum exact Ct values if possible. After rRT-PCR reaction, all amplicons were subjected to 3 individual SENA reactions, i.e., N-SENA, O-SENA and mix-SENA with crRNAs targeting O gene, N gene and both, respectively.
Consistent with the theoretical analysis (20) and rigorous experimentation (21), along with the decrease of the RNA templates to less than 3 copies/Rx, the rRT-PCR Ct values in some replicas, primarily that corresponding to the N gene, passed 38 (the cut-off for positive as recommended by the rRT-PCR kit suppliers) but were less than 40, which should be considered as entering the “grey zone”. The Ct values increased steadily when the concentration of the RNA templates further decreased, with more and more replicas showing one or both Ct values entering the “grey zone” and eventually all became “negative”, i.e., greater than 40 or even 45 (Figure 2A, Table S2 and Figure S1). Employing Ct=38 as the cut-off for “positive” detection, we estimated the LoD for O and N genes with 95% confidence interval (CI) of this set of rRT-PCR assay as 3.3≤4.0≤6.1 and 4.0 ≤ 4.1 ≤ 4.4, respectively (Figure S2). Most likely due to the influence of the complex combination of the targeted viral genomic fragments and the clinic sampling background, the LoD determined in this study was clearly higher than the published value of 2.0≤2.5≤3.7, which analyzed single target in a pure system (21).
The rRT-PCR amplicons were further analyzed by SENA detection with the measurement of the fluorescence signals for each corresponding replica. After comparison of the parameters of slope (increase of fluorescence/min) versus FC (the fold of change of fluorescence between that of the sample over that of the negative control at certain time point), we defined a parameter, FCratio, which is the ratio of the FC at 10 min to that at 5 min after the initiation of fluorescence reading (Figure S1 and Table S2). We also found that in the cases with low concentrations of templates, the rRT-PCR efficiency of the two target genes (i.e., O and N) were different so as the SENA detection (Figures S1 and S2). In order to verify the existence of specific amplicons of SARS-CoV-2 nucleic acids in an individual rRT-PCR reaction, all of the amplicons of the replicas with RNA templates ranging from 0.125 to 2 copies/Rx were subjected to next generation sequencing (NGS) analysis. The results were found to be completely consistent with the perspective results of both O-SENA and mix-SENA. In addition, with mix-SENA, not only the signals are generally more significant than that of the O-SENA detection but also may resolve some of the ambiguity readouts found with N-SENA (Figure S1 and Table S2). Based on these results, the mixFCratio was demonstrated as the most sure-proof index for rRT-PCR confirmation, and we empirically estimated that mixFCratio≥1.145 for positive cut-off, and mixFCratio≤1.020 for negative cut-off (Figure 2A). Of course, these two parameters are subject to further verification and adjustment along with the increase of tested samples. Because SENA is rRT-PCR based, the same methodology for determining the rRT-PCR LoD was used to estimate that of SENA by this set of data, corresponding to both individual O and N fragments (Figure S2) and in combination as indicated by the mix-SENA (Figure 2B). As expected, the N-SENA LoD (3.7≤4.3≤4.8 with 95% CI) is very close to that of the N-Ct of rRT-PCR, while the LoD of O-SENA (1.1 ≤ 1.3 ≤ 1.7 with 95% CI) is significantly lower than that of O-Ct (Figure S2). Although the LoD of mix-SENA (1.2≤1.6≤2.1 with 95% CI) is slightly higher than that of O-SENA (Figure 2B and Figure S2), it is apparently caused by its capable of confirming some of the ambiguous amplicons in the extremely low concentration cases (Figure S2) and thus, mixFCratio is chosen for clinic applications.
SENA was further verified in a few hospitals, testing various clinic specimens and samples under different scenarios (Figure 1) and employing few more commercial rRT-PCR diagnosis kits in addition to BJ which was used in the titration experiment (Table S1). Totally 295 clinic samples or specimens (mainly pharyngeal swabs) collected from 282 individuals were tested by rRT-PCR followed by SENA detection (Table S3). Except for asymptomatic carriers, all the cases of uncertain analytic and false positive or negative readouts of rRT-PCR diagnosis were encountered and finally confirmed or corrected by SENA detection.
Specifically, samples from 139 patients of Ruijin Hospital (RJ, Shanghai, China) were assayed by rRT-PCR employing diagnostic kits of ZJ and HD, 137 of which had consistent readouts by both rRT-PCR kits, indicating two positive, 123 negative and 12 suspected that fell in the “grey zone” (Table S3). SENA detection of these samples revealed not only the 12 suspected as negative but also identified one more positive among the original 123 negative individuals, clearly a case of false negative diagnosis (Table S3). Besides, distinct rRT-PCR assay results, positive by HD but negative by ZJ were shown for samples collected from 2 close contacts of COVID-19 patients and apparently asymptomatic (ref to Table S3). However, the amplicons of both ZJ and HD were shown as negative via SENA detection. All these ambiguous rRT-PCR amplicons (17 samples, ref to Table S3) were finally analyzed by NGS, and the results were consistent with the SENA. Noticeably, the rRT-PCR false-negative COVID-19 patient was symptomatically mild at the point of admission with all the clinic laboratory tests negative but turned positive after 24 hours. On the other hand, although those 12 suspected patients had respiratory infection symptoms, they were finally excluded from COVID-19 according to the latest guideline for diagnosis and treatment from China National Health commission (the 6th edition). Similarly, in Shenzhen Second People’s Hospital (SZII, Shenzhen, China), 5 uncertain rRT-PCR readouts for O gene were found among 139 individuals. Three of them had Ct value of 39.47, 39.7 and 40.56, respectively but the following SENA detection gave mixFCratio values less than 1.0 for all of them, indicating all negative. The other two individuals had Ct values of 38.87 and 39.22, while their mixFCratio values were 1.581 and 1,609, respectively, indicating positive for both. In addition, there were another three individuals with Ct values larger than 40 for O gene and 36.09, 35.88 and 37.98 for N gene, respectively; however, the following SENA detection showed mixFCratio values were 1.39, 1.55 and 1.21, respectively, indicating all positive. All these amplicons were further confirmed by NGS analysis (Table S3), obtaining consistent results with those of SENA. Consistently, the three SENA-negative individuals were finally excluded from SARS-CoV-2 infection after being rechecked by rRT-PCR after 24 hours (Table S3). Based on above data, it is clear, SARS-CoV-2 infection suspects with either rRT-PCR Ct values falling in the “grey zone” or with clear patient-contact epidemiological history but negative rRT-PCR tests, are strongly recommended to perform SENA detection to minimize the possibility of misdiagnosis. On the other hand, in case an rRT-PCR-positive suspect does not demonstrate any COVID-19 clinic symptoms and/or signs, SENA detection is also strongly recommended to eliminate either false-positive diagnosis or misdiagnosis of the so-called “asymptomatic carrier” or “asymptomatic patient”.
Besides of preventing false-negative or false-positive diagnosis, the highly sensitive property of SENA may also assist in providing evidence of viral clearance for COVID-19 recovering patients. A female patient in Dongfang Hospital (DF, Shanghai, China) was confirmed as COVID-19 positive by both rRT-PCR and CT scanning and showed ground-glass opacities mixed with consolidation along the subpleural area (Figure 3). Accordingly, the SENA test was positive with the mixFCratio of 1.43. After the hospitalization, the patient was further analyzed by rRT-PCR at two time points, obtaining all negative results with bilateralnasal and pharyngeal swab specimens. However, the mixFCratios of SENA for some of her specimens were 1.64, 1.36 and 1.00, respectively, indicating that the virus was contained and yet to be cleared. On the sixth day, both rRT-PCR and the corresponding SENA detection for all of her specimens were negative and these results were confirmed by NGS and consistent with her normal CT scanning results (Figure 3). Thus, she was discharged from the hospital and safely back to home. Similar cases were found in Jinan of Shandong Province, China, where the fecal samples from two recovering COVID-19 patients were tested negative by rRT-PCR but clearly positive by SENA (Table S3). Considering a certain percentage of the recovered patients discharged from hospitals were reported to be re-detectable positive (RP) (22), the incomplete clearance of the SARS-CoV-2 virus ahead of discharge might be one of the possible causes. Therefore, it could be necessary to consider more sensitive detection approaches such as SENA as a potential index of viral clearance.
To reconfirm and/or improve the cut-off values for SENA mixFCratio, the ambiguous Ct values were re-estimated using the regression functions derived from the rRT-PCR assays with titrated standard RNA templates (Figures S4 and S5), and then the Ct values (both estimated and detected) were plotted against the corresponding mixFCratios (Figure 4). Combining the data from both RNA standards and clinic samples, it is clear that SENA detection is of both high sensitivity, identifying real positive samples with Ct values as high as more than 43 (approaching 50 as estimated), and high specificity, identifying real negative samples with Ct values as low as 39. Therefore, SENA can effectively eliminate uncertain diagnosis of rRT-PCR assays for SARS-CoV-2 infection. In addition, the cut-off value for SENA mixFCratio remains unchanged as 1.145 for positive diagnosis while slightly increased to 1.068 for negative (Figure 4), which is supposed to further increase along with the clinic applications.
Instead of developing a closed CRISPR-Dx system, which, ideally, should be comprised of both target nucleic acids amplification and CRISPR-Cas-based trans-cleavage assays, SENA was created here to match the commercially available and widely applied rRT-PCR kits, and to solve the uncertainty challenge of the rRT-PCR “grey zone” in COVID-19 diagnosis. Although there are dozens of commercial rRT-PCR kits using distinct PCR primers and probes, their corresponding SENA detection kits can be easily designed and prepared following the instructions provided in this study. Considering the fact that qPCR is the most popular MDx system and SENA is simple to operate, SENA thus has the potential to be widely used in various scenarios to solve the uncertainty problems of qPCR as well as other nucleic acid amplification-based MDx. Of course, to minimize the possibility of aerosol contamination during opening of the PCR tubes and pipetting, physical separation of the SENA detection center from the clinical PCR laboratory is an absolute requisition. Abided by this rule, SENA has been demonstrated convenient and effective in several hospitals and centers for disease control and prevention as part of their laboratory routine in combination with rRT-PCR for more sensitive and accurate detection of SARS-CoV-2 infection. Therefore, SENA is a useful technique that meets the urgent needs of combating COVID-19 pandemic.
Data Availability
All data is available in the main text or the supplementary materials.
Competing interests
G.Z. and J.W. are cofounders of Tolo Biotechnology Co., Ltd. W.H., G.Z. and J.W. have filed patent applications relating to the work in this manuscript.
Data and materials availability
“All data is available in the main text or the supplementary materials.”
Supplementary Materials for
This PDF file includes:
Ethics statement
Materials and Methods
Figs. S1 to S5
Tables S1 to S3
Ethics statement
The human pharyngeal swab samples were collected from residual samples that had been obtained for clinical respiratory virus detection. All procedures complied with the Measures for the Ethical Review of Biomedical Research Involving Human Subjects issued by the National Health and Family Planning Commission of The People’s Republic of China. All protocols were approved by the Ethics Committees of Ruijin Hospital, Dongfang Hospital and Shenzhen Second People’s Hospital. All procedures were performed in biosafety level 2 facilities.
Materials and Methods
1. RT-qPCR reactions
Reactions were conducted in a 25-μl reaction mixture following the instructions offered by the commercial suppliers of the reaction kits (ref to Supplementary Table S1). Usually, the reaction cycle parameters were set as reverse transcription at 50 °C for 10 min, denaturation at 95 °C for 5 min, then followed by 45 cycles of amplification, i.e., 95 °C for 10 s and 55 °C for 40 s.
2. SENA reagents and detection
2.1 Preparation of SENA reaction reagents
Candidate crRNA guide sequences with the “TTN” PAM sequences were designed as shown in Supplementary Table S1, and the crRNAs were prepared following the procedures previously described (11). Except the target DNA, the 2× SENA reagent comprises of 2× NEB buffer 3.1, 500 nM LbCas12a (Tolo Biotech.), 1 μM synthesized crRNAs (for each specified target and thus varies according to the RT-qPCR kit supplier, Table S1), 1 μM FQ-reporter, and 1 U/μl RNase inhibitor (TaKaRa).
2.2 SENA detection
To avoid the aerosol contamination of the MDx laboratory, after the RT-qPCR reaction, their products must be transferred to a physically isolated room to perform SENA detection. It is also important to choose proper SENA detection reagents corresponding to the RT-qPCR kits. To prepare a 20-μl SENA reaction system, with corresponding positive and negative controls, 2-μl PCR products and 8-μl RNase-free H2O were mixed with 10-μl 2× SENA reagent, and the mixture was then measured on an appropriate fluorescence reader with FAM fluorescence collected following the programs: 48 °C 30 s per cycle, 20 cycles. Both the slope and the Fluorescence Change (FC) can be calculated at any time points as desired.
3. Next generation sequencing (NGS)
The RT-qPCR products were purified by AMPure XP beads (Beckman Coulter Life Sciences, US), and libraries were then constructed following the procedures of end repair, dA-tailing and adaptor ligation, with the StepWise DNA Lib Prep Kit for Illumina (ABclonal, China). After PCR amplification, samples were sequenced on Illumina Miniseq to produce 2×150 bp paired-end reads. After adaptor trimming and quality trimming, the clean reads were mapped to the reference genome of SARS-CoV-2 (MN908947.3) using Bowtie2 (22).
4. Systematic titration and regression analyses
4.1 Systematic titration experimentation
4.1.1 The standard RNA templates
The SARS-CoV-2 RNA standards were purchased from Genewell (Shenzhen, China). According to the supplier’s information, three plasmids containing the whole sequences of N and E genes, and partial sequence of the Orf1ab, i.e., from 13237 to 13737 of the SARS-CoV-2 complete genome (MN908947.3), were transcribed in vitro individually. The RNA products were mixed with equal molar, aliquoted with addition of 1 μg of human RNA per tube, and subject to lyophilization and subsequent quantification with digital PCR. This SARS-CoV-2 RNA standard dry powder containing 1808 copies of O gene, 179.5 copies of N gene and 1160 copies of E gene was dissolved with 10 μl RNase-free water to obtain the original stock solution (estimated 180.8 copies/μl of O, 179.5 copies/μl of N and 116 copies/μl of E).
4.1.2 The preparation of the serially diluted RNA templates
Pharyngeal swab samples were collected from 40 adult patients in Shenzhen Second Peoples’ Hospital by the Clinic Diagnosis Laboratory and the nucleic acids of each sample were extracted with the pre-packaged nucleic acid extraction kit (Da’An Gene., Ltd., Guangzhou, China), according to the manufacturer’s instructions, ended up with 55-μl extracts per sample. After RT-qPCR assays employing 5 μl of the extracts from each sample, all of the samples were shown to be SARS-CoV-2 negative. The remaining 50-μl extracts of each sample were mixed together and 5 μl of the mixture was once again analyzed by RT-qPCR and confirmed to be SARS-CoV-2 negative. Then, the mixed nucleic acid extract was used as the dilution buffer (totally about 2 ml) for serial dilution of the SARS-CoV-2 RNA standard stocks, generating desired concentrations (i.e., 5, 2, 1, 0.8, 0.6, 0.4, 0.2, 0.1, 0.05, 0.025, 0.01, 0.005 copies/μl), and 5 μl of each of the diluted solutions were used as templates for RT-qPCR analysis, forming gradient template concentrations (i.e., 25, 10, 5, 4, 3, 2, 1, 0.5, 0.25, 0.125, 0.05, 0.025 copies/Rx).
4.1.3 Replica setting
We analyzed 9 replicas for each of the concentrations of 1 and 0.5 RNA template copies/Rx while 6 replicas for each of the rest concentrations (Supplementary Figure S1).
4.1.4 RT-qPCR reactions
RT-qPCR reaction kit was supplied by BJ (Table S1), who follows the primer sets recommended by Chinese CDC. Instead of the recommended 40 cycles of PCR amplification, we set 45 cycles as routine aiming at recording the maximum exact Ct values if possible.
4.1.5 SENA detection
Every amplicon of the RT-qPCR reactions was subjected to SENA detection. Specially, for this experimentation, three sets of SENA reagents were individually used, O-SENA contains the crRNA targeting the O sequence, N-SENA contains the N-targeting crRNA, and the mix-SENA contains crRNAs for both sequences.
4.2 Choice of FCratio as the standard readout for SENA detection and mix-SENA as the standard reagent for clinic application
Three readout parameters were compared as shown in Figure S1 (original data in Table S2). The slope (increase of fluorescence/min) represents the reaction rate of Cas12a trans-cleavage activity, but in this experiment, it represents neither the initial rate of the enzyme under limited substrate condition nor the pure first-order reaction rate varies according to the substrate concentration, particularly, demonstrated in cases of high template concentration, in which, the slope goes down along with the increase of the template concentration. In addition, when the substrate concentration is low, the slope of SENA is hard to be distinguished from that of the negative control, ending with ambiguous cut-offs. The FC (the fold of change of fluorescence between that of the sample over that of the negative control at certain time point, usually 5-30 min) does show clear differences between the positive amplicons from that of negative control, and it also shows certain quantitation character particularly at the low concentration templates cases. Because of these properties, FC has been a parameter used by a few users of CRISPR-Dx (24). However, it seems that the absolute value of the FC usually varies along the reaction time and sometimes it is influenced by the change of the fluorescence signal of the negative control. Besides, it is difficult to determine the “best choice” of the FC recorded at certain time points, which may cause confusing in clinic applications. It is clear, we need a stable readout which reflects the dynamic process and the quantitative correlation of SENA reaction with the low concentration of the templates on one hand and should be robust and accurate for clinical diagnosis on the other hand. We defined FCratio, which is the ratio between FC’s of SENA detection at 10 min versus 5 min after the beginning of the fluorescence reading. It not only measures the fluorescence change of SENA against the negative control background so that the quantity of the amplicons, particularly at the low concentration range may be represented, but also normalizes the slope of the fluorescence change of SENA to eliminate the complex background differences. As shown in Figure S2, FCratio significantly amplified the positive signals and represents the quantitation of the amplicons at low template concentrations to certain extent.
As shown in Figure S1, the capacity of the three SENA reactions are compared. It is obvious that N-SENA is the least sensitive one, while although O-SENA seems much more sensitive than that of N-SENA and largely comparable to that of mix-SENA, its signal at the very low template concentration range seems uncertain in some cases. Therefore, for clinic application, mix-SENA is the choice.
4.3 Quality analyses of the RT-qPCR titration data
The quality of the Ct values vs the concentration of the standard templates for RT-qPCR of the systematic titration experiment were analyzed both empirically and statistically. Firstly, as shown in Table S2 and Figure S1, the amplification efficiencies of the two genes represented by the valuable Ct readouts were different. The apparently lower sensitivity of the N gene amplification is in contrast with the clinic experiences and might due to the difference in the property of the templates used in different experiments (laboratory standard RNA template vs clinical real viral template). Secondly, although linear regression can be readily made between the Ct-values and the log10 (conc) (Figure S3A) as that of the previously published tests, the quality of the regression as judged by the R2s (Figure S3A) and the residues (Figure S3B) are clearly suboptimal likely due to the limited number of replicas in the experimentation. On the other hand, this titration was designed with taking at least the two most fundamental limitation factors about the sensitivity of COVID-19 RT-qPCR diagnosis into consideration, i.e., the sampling ambiguity and the influence of the biological/chemical contaminants from the clinic samples. Therefore, the data will be used for determination of the LoD for RT-qPCR and SENA (Materials and Methods 4.4 and Figure S2), and the regression function will be used to estimate the “apparent Ct” values of the samples with Ct-values greater than 45 (no amplification signal) but their SENA detection is positive (Materials and Methods 4.5 and Figure S4).
4.4 Determination of limit of detection (LoD) for RT-qPCR and SENA
The LoD values for RT-qPCR (N-Ct and O-Ct) and SENA (N-FCratio, O-FCratio and mix-FCratio) were estimated based on the systematic titration employing standard RNA templates (Figure S1 and Table S2). The fractions of positive replicates versus the number of target molecules (copies) per reaction for N and O gene of COVID-19 were plotted and used the sigmoidal function (1) to fit the data via R (software version 3.5.0). The 95% confidence intervals were derived by bootstrapping the model residues and were visualized by R (software version 3.5.0) with built-in ggplot2 library (21).
4.5 Regression of RT-qPCR Ct-values and the SENA FCratio versus the concentration of the templates employing the data from the systematic titration
4.5.1 Regression of RT-qPCR Ct-values with the concentration of the templates (copies/Rx)
Since PCR product increased exponentially with the initial concentration of the sample (x), and the Ct value of RT-qPCR parameter (y) was inversely correlated with the initial concentration, especially in the range of low copy number (low template concentration) samples, the power function equation (a<1) should be suitable for the data fitting. However, some of the experimental groups included very low initial sample concentrations (<1 copy/Rx), those amplification efficiencies should be different (particularly affected by sampling ambiguity) from that of the groups with high initial template concentration. Therefore, the power function formula with four parameters (Y=aXn+bXm) was used to match all the experimental group data to obtain a more accurate data model (Functions 1, Figure S4).
4.5.2 Regression of SENA FCratio with the concentration of the templates (copies/Rx)
The exponential function (first order association kinetics of the interaction between a substrate and an enzyme, Y=a+b(1-ecX)) is used to fit the data of FCratio against the concentration of the templates. At low concentration (especially when the concentration is less than 2 copies/Rx), the FCratio is positively correlated with the template concentrations. However, when the template concentration reaches to 2 copies/Rx and more, the FCratio does not increase accordingly and the curve tend to be flatted out. In addition, as FCratio have been already normalized by the fluorescence signal of the negative background, it is stable and, in this case, we give the parameter Y0 being set as a constant value between 0.9 to 1 (Functions 2, Figure S5).
4.5.3 Regression of RT-qPCR Ct-values with the SENA mix-FCratio values
In practice, quite significant portions of the clinic positive samples detected by mix-SENA with their FCratio readings higher than the positive cut-off, but with a negative PCR Ct value (>40-45, depending on the scenario). Under certain circumstances, people may be interested to learn the copy number of the templates for the corresponding RT-qPCR assays or even the “probable” Ct values of these assays. With the aid of the above-mentioned two regression functions (1 in Figure S4 and 2 in Figure S5), these data could be estimated. One may firstly substitute the Y in Function 2 by the measured FCratio value and the corresponding X can be calculated representing the “estimated concentration of the template). Then this X value can be used to estimate the corresponding “estimated Ct-value” as the Y of Function 1. We estimated all the ambiguous Ct-values of the positive amplicons and plotted them against their corresponding mix-FCratio (Figure 4). All the real and estimated Ct values for both N and O genes are plotted against the corresponding FCratio values of mix-SENA as X axis. An exponential decay function (with X0 = 1; When X≤X0, Ct=∞; otherwise, one phase decay) fits well to all the data (R2=0.9238) and is used for analyzing the clinic data and adjust the cut-off values accordingly (Figure 4).
5. Detection of antibodies
The detection of anti-SARS-CoV2 antibodies was executed by the point-of-care microfluidic platform integrating a home-made fluorescence detection analyzer (Suxin, Shanghai, China). A total of 10-μl plasma was added into the loading chamber of microchip followed by the addition of 70-μl sample dilution buffer. After incubation for 15 min at room temperature, the microchips were loaded onto the fluorescence detection analyzer, and fluorescence signal was detected from the analyzer, following the manufacturer’s instruction.
Supplementary Figures
Acknowledgments
We thank Yucai Wang and Linxian Li for discussions and support. J.W. is supported by the National Natural Science Foundation of China (31922046); W.H. is supported by the National Key R&D Program of China (2019YFA0906000) and NSFC (81772737); D.W. is supported by NSFC (81972830) and Shanghai “Rising Stars of Medical Talent” Youth Development Program-Clinical Laboratory Practitioners Program; F.S. is supported by NSFC (31670757); G.Z. is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB19040200).
References
References
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.