key: cord-0820249-jioroc9x authors: Li, Ziyan; Chen, Xue; Huang, Zili; Zhou, Jing; Liu, Rui; Lv, Yi title: Multiplex Nucleic Acid Assay of SARS-CoV-2 via a Lanthanide Nanoparticle-Tagging Strategy date: 2021-09-08 journal: Anal Chem DOI: 10.1021/acs.analchem.1c02657 sha: ba32f8d74acbb324749b42bed0b0ff0372edc452 doc_id: 820249 cord_uid: jioroc9x [Image: see text] Early diagnosis, early isolation, and early treatment are efficient solutions to control the COVID-19 pandemic. To achieve the accurate early diagnosis of SARS-CoV-2, a multiplex detection strategy is required for the cross-validation to solve the problem of “false negative” of the existing gold standard assay. Here, we present a multicomponent nucleic acid assay platform for SARS-CoV-2 detection based on lanthanide nanoparticle (LnNP)-tagging strategy. For targeting SARS-CoV-2’s RNA fragments ORF1ab gene, RdRp gene, and E gene, three LnNP probes can be used simultaneously to identify three sites in one sample through elemental mass spectrometry detection with limits of detection of 1.2, 1.3, and 1.3 fmol, respectively. With the multisite cross-validation, we envision that this multiplex and sensitive detection platform may provide an effective strategy for SARS-CoV-2 fast screening with a high accuracy. T he pandemic of novel coronavirus disease (COVID-19) 1 by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to a surge of a series of serious public health, medical, and society issues. Early diagnosis, early isolation, and early treatment are efficient solutions to control the COVID-19 pandemic. As one of the diagnostic criteria, computerized tomography images can only provide a diagnosis accuracy of 76.4%. 2 Antibody-based serological tests are not suitable for the early diagnosis 3,4 because they are not produced until 3 weeks after infection. 5 The two currently used nucleic acid pathogenic diagnoses are nucleic acid sequencing and nucleic acid fluorescence quantitative polymerase chain reaction (qPCR). 6 Nucleic acid sequencing is relatively time-consuming and expensive. 7 As the current gold standard, qPCR falls in a dilemma of "false negative". 8, 9 Up to 20 to 67% of patients are tested positive only after multiple repeated tests, 10 as also described in National Health Commission of China's new coronavirus pneumonia diagnosis and treatment plan. 11 To achieve the accurate diagnosis of SARS-CoV-2, a multiplexed detection strategy is required for multiple characteristic-related markers for the cross-validation. The U.S. Centers for Disease Control (CDC) gives examples of recommended gene selection and multicomponent testing methods. 12 World Health Organization (WHO) recommends three qPCR reactions to confirm COVID-19 infection, while a total of 3 qPCR tests were performed sequentially for each patient sample. 13 Besides, electrochemical biosensors, 14 colorimetric sensors, 15 a loop-mediated isothermal amplification (LAMP) combined with clustered regularly interspaced short palindromic repeats (CRISPER) system, 16 ,17 microfluidic assay, 6 and point-of-care testing (POCT) platforms 18 have been developed for single or multicomponent detection. 19−23 Despite great success, 24 these optical and electrochemical strategies are practically limited by spectral overlap-caused interference for the simultaneous multicomponent detection. 25 Elemental mass spectrometry, especially inductively coupled plasma mass spectrometry (ICPMS), 26−28 is emerging as a powerful platform for the multiplex bioassay 29−31 due to its high sensitivity, 32 wide dynamic range, 33 and multiplex sensing ability. 34−36 The multiplex bioassay by ICPMS 37 possesses many superiorities, such as a low sample consumption, 38 a short analysis time, a simplified procedure, and a low cost. 38−40 Zhang's group developed a multiplex nucleic acid assay and immunoassay by lanthanide macrocyclic compound tagging. 41, 42 Wang et al. performed a multiplex assay for protease markers. 43 Simultaneous detection of multiple microRNAs by the MNAzyme-catalyzed amplification assay was realized, 44 while single-particle mode ICPMS was applied for a threecomponent homogeneous immunoassay of pancreatic cancer. 45 Multicomponent marker and disease subtype analysis provides the possibility of the cross-validation through a simultaneous monitoring of multiple parameters, 46 which plays a significant role in improving the accuracy of the disease diagnosis. However, the existing assay formats are not compatible with the direct detection of SARS-CoV-2, which is mainly limited by the low sensitivity of macrocyclic compound element labels. Herein, we proposed a multiplexed nucleic acid assay for improving the detection rate of SARS-CoV-2 via lanthanide nanoparticle (NP)-tagging strategy. Three lanthanide NP (LnNPs, Ln = Tb, Ho, and Eu)-based nanoprobes were designed and applied to detect SARS-CoV-2's ORF1ab gene, RdRp gene, and E gene. Three kinds of capture DNA attached to magnetic beads (MBs) could hybridize with corresponding targets, and the half remaining chain could also hybridize with the as-prepared probe DNA-modified NPs. LnNP tags were digested in the enhancer solution and released lanthanide ions from NPs to those of the β-NTA-Tb 3+ /Ho 3+ /Eu 3+ −TOPO complex. 47 The LnNPs are ideal multiplex labels 48 because lanthanides are numerous and similar in nature, 49 and there are 55 to 64 atoms exist per nm −3 NP. 50 Thanks to multiatoms inside each NP, a high sensitivity is obtained without the need of complex nucleic acid amplification procedures. Apparatus and Reagents. A PerkinElmer NexION 350 inductively coupled plasma mass spectrometer was used to detect 153 Eu, 159 Tb, and 165 Ho in standard mode. The operating conditions are listed in Table S1 . Apparatus for characterization and detection is provided in the Supporting Information. Other details about reagents and materials used are also provided in it. DNA and RNA sequences used in this work are also listed in Table S2 , Supporting Information. Synthesis of −COOH-Functionalized NaLnF 4 (Ln = 153 Eu, 159 Tb, and 165 Ho) NPs. Carboxylic acid-functionalized NaLnF 4 (Ln = 153 Eu, 159 Tb, and 165 Ho) NPs were synthesized as bioprobes adopting a one-step solvothermal route by employing polyacrylic acid (PAA) as the surfactant and a capping ligand, both of which provide excess −COOH for further bioconjugation and improving NPs' hydrophilicity ( Figure 1a ). Details are given in the Supporting Information. Preparation of Probe DNA-NaLnF 4 NPs (O/R/E-Tb/Ho/ Eu NPs). To prepare probe DNA-NaLnF 4 NPs as a signal reporter, NH 2 -modified probe DNA was directly attached with the −COOH group coating on the NaLnF 4 surface through the EDC−NHS (N-ethyl-N′-(3-(dimethyl-amino)propyl)carbodiimide-N-hydroxy succinimide)reaction. Specific details of the labeling process are given in the Supporting Information. After ligation, NPs were blocked with 2% bovine serum albumin (BSA) in TBS buffer to prevent protein's nonspecific binding and the residue carboxylic acid groups were blocked with excess tris in TBS buffer. Subsequently, the as-prepared probe DNA-NaLnF 4 NPs were diluted to 1.8 mL per OD DNA and made ready to use. Preparation of Capture DNA-MBs (Capture O/R/E-MB). The process of attaching oligonucleotides on MBs was carried out according to the instructions of Dynabeads M-280. Ahead of all, pretreatment to avoid RNase contamination for RNA manipulation was performed, then biotinylated capture DNA was immobilized, specific details are given in Supporting Information. Following that the MBs with 2% BSA were blocked to prevent nonspecific binding, and the DNA-coated MBs were redispersed in 100 μL Tris-HCl buffer stored at 4°C for further use. Analysis of ORF1ab by Sandwich Assay-ICPMS. Before multicomponent RNA fragment detection, the feasibility and conditions of a single component assay were surveyed by utilizing ORF1ab as a model SARS-CoV-2 positive gene. To fulfill it, a certain amount of capture O-MBs was used for each sample with blocking ahead, 0 and 2 pmol ORF1ab were added under different conditions. Subsequently, probe O-NaTbF 4 NPs participate and react for fully chain hybridization. After the reaction is completed, 200 μL enhancer solution was utilized to release ions from the formed structure for ICPMS quantification. Simultaneous Analysis of ORF1ab/RdRp/E by Sandwich Assay-ICPMS. The RNA sandwich hybridization system was 7.5 μL capture DNA-MBs. Because of the specific recognition between capture O/R/E-MB, targets O/R/E, and probe DNA-NaLnF 4 NPs (O/R/E-Tb/Ho/Eu NPs), a sandwich-like structure was formed with MBs as both solid- Figure S1 ) for the three lanthanides we adopted, we can convert the amount of target according to signal strength. Preparation of Viral RNA Extraction for Standard Recovery Experiment. The real sample matrix of three throat swab samples for the recovery experiment was collected shows lattice fringes with an observed d-spacing of 0.300 nm, which is in good agreement with the lattice spacing in the (101) planes of hexagonal NaTbF 4 (PDF #27-0809). The Xray powder diffraction (XRD) patterns of the NPs exhibit peak position and intensities that can be well indexed in accord with NaTbF 4 crystals ( Figure S2a ). Elements and corresponding proportions of a kind of particles by X-ray photoelectron spectroscopy (XPS) analysis give a peak area of Na, Tb, and F about 1:1:4 ( Figure S3a ). The one-step method also wraps the carboxyl group in the surface of the synthesis of the NPs, and the zeta potential of the naked NPs is positive. The NPs wrapped in the PAA encircled by the method exhibit negative electrical properties (Figure 2d) . Similarly, in the Fourier transform infrared spectroscopy (FTIR), the PAA-modified NPs exhibit a carboxyl group characteristic peak compared to bare NPs (Figure 2f ). Those of the other two NPs were similar, and detailed properties characterized by TEM and high resolution TEM (HRTEM) (Figures S4, S5 Characterization of Probe DNA-NaLnF 4 NPs. Characterization of Probe DNA-NaLnF 4 NPs is essential for our increased understanding of bioconjugation. This modification process can be characterized by multiple means. To understand the contents of the DNA attached to the nanoprobe we added for each sample, we used the high-resolution ICPMS (HR-ICPMS) to quantitatively analyze 31 P on the nanoprobes. 51 We use the same length of DNA as used for a standard curve of the P element, then make a quantity of probe DNA (Figure 2c) , and it is concluded that each sample substantially corresponds to 8 pmol of probe added, with a similar connecting efficiency of each NP. The modification process was also validated by zeta potential with noticeable change (Figure 2d ). In addition, dynamic light scattering (DLS) determined the increase of the particle size after DNA ligation. Characterization of the other NPs using zeta potential ( Figure S6 ) and DLS ( Figure S7 ) is given in Supporting Information. Characterization of Capture DNA-MBs. The capture DNA-MBs were prepared by coupling all three capture DNA with the same molar concentration onto SA-MBs in a one-pot process because it has similar analysis performance compared to separate coupling ( Figure S9 ). In theory, each type of capture DNA strand accounted for one-third of the total attached capture DNA. Changes in zeta-potential ( Figure S10a Figure S11 ). Another more intuitive evidence is that under TEM, it shows naked MBs without immobilization in Figure 3a . After hybridization, three kinds of NPs were attached. Several NaTbF 4 , NaHoF 4 , and NaEuF 4 could be observed clearly (Figure 3b) . A similar phenomenon could also be observed in the SEM image ( Figure S12 ) with its inset graph giving the compositional analysis by energydispersive X-ray spectroscopy (EDS). There emerged elemental peaks of Tb, Ho, and Eu in the EDS diagram, and small particles could be found on MBs (Figure 3c ). The above evidence demonstrated successful hybridization, and a sandwich structure was formed. To further prove this strategy's feasibility, five artificial samples containing different concentrations of ORF1ab, RdRp, and E were prepared and hybridized with capture O/R/E-MBs and O/R/E-Tb/Ho/Eu NPs. With the increase in ORF1ab, RdRp, and E concentrations, the signal intensities of 159 Tb, 165 Ho, and 153 Eu also increased (Figure 3d) . Optimization of Experimental Conditions. Experimentally, one of the RNA fragments, ORF1ab with the corresponding tagging NP NaTbF 4 , was chosen as the model RNA for optimizing experimental conditions. It has a similar length to the other two components of the test object, so the corresponding characteristics such as Tm are similar. Based on this, various experimental parameters were optimized using blocking methods, reaction time and temperature, and the number of components of this system (capture DNA, capture-MB, and probe-DNA). The specific details about the reaction are given in Supporting Information. Effects of Different Blocking Methods. NPs will have nonspecific adsorption of some molecules due to their nanoscale effects. We use proteins to block the remaining sites of the NPs that have been bound to DNA that may be nonspecifically adsorbed. Figure 4a reveals that there has been a steep decline in the signal when the blocking condition is not suitable for this system. The blocking condition we finally selected is 2% BSA, 30 min. Effects of the Amount of Capture-DNA Added to MBs. In a further set of experiments, we explored the reaction ratio of MBs to captured DNA (Figure 4b ). According to the binding capacity instructions of the MBs, we have selected a range of DNA concentrations to achieve the maximum binding efficiency and at the same time to ensure that DNA addition is not too excessive to cause insufficient washing in which free DNA might lead to interference on subsequent sandwich hybridization reactions. Finally, the condition we selected is the group with the highest signal-to-noise ratio, that is, when the concentration of capture DNA is 2 μM, the volume ratio of DNA to MBs is 10:1. Effects of Hybridization Time. Hybridization time is a common condition that has a considerable impact on reaction efficiency. To guarantee enough time and control nonspecific adsorption, we explored the reaction efficiency in the time range of 60−180 min, as shown in Figure 4c . Taking signal growth and blank control as dual considerations, as time grows, the signal value first shows an upward trend and then gradually stabilizes, while the blank value first decreases due to the competition of the target molecule and then increases due to time. Considering the signal-to-noise ratio, we choose 150 min as the best hybridization time. Effects of the Titer for As-prepared Probe-DNA. What can be seen in Figure 4d is the influence of probes of different titers on the experimental results. For 60 μL NPs and 1OD DNA per tube; the corresponding condition is the diluted volume of each probe, ranging from 0.6 to 3.0 mL. From the changing trend of the signal and blank value in the figure, when the probe is more concentrated and the titer is relatively low, the NPs on the excess probe will produce more nonspecific adsorption on the reaction substrate. As the temperature increases, the blank gradually decreases, but if the probe is too dilute, the concentration of probe DNA used to bind the target will decrease, and the reaction speed and efficiency will subside, too, making it impossible to obtain high signal values. Taking these effects into consideration, we assumed that the titer 1 to 1.8 mL with the highest signal-to-noise ratio was selected for applying under subsequent experimental conditions. Effects of the As-prepared MB-Capture Amount. As the unit of the sandwich structure, MB-capture and Probe-DNA have similar influence trends on the reaction. Under the premise of ensuring the amount for reaction, as a reaction matrix, the increase of MBs is positively correlated with the increase in nonspecific adsorption (Figure 4e ). Therefore, we choose a moderate amount, 2.5 μL per sample, and perform subsequent reactions under this condition. Effect of the Hybridization Temperature Program. We studied the appropriate time for the interaction. Under the same reaction time, when the selected reaction temperature is below the melting temperature (T m ) of 15−20°C of the target sequence, the reaction would be more efficient; the lower the temperature, the more stable the double-stranded binding DNA. To explore the checks and balances of these two factors, we have designed a series of hybridization temperature programs as shown in Figure 4f . When the temperature is 43°C , the best signal-to-noise ratio can be obtained. Finally, 43°C was chosen as the optimum reaction temperature. Specificity and Selectivity. When multiple components are the target to be analyzed, it is necessary to ensure that there will be no mutual influence between them. We experimentally prove the mechanism of this strategy and demonstrate its specificity. We have prepared a series of target gene fragment solutions with different components, as shown in Figure 5a . Similarly, there may be a variety of nontarget nucleic acid fragments in real samples. Due to the high sequence homology of RNA, we detected a series of sequences similar to the three targets to prove the selectivity of the established strategy. The results are shown in Figure 5b . Experimental groups were containing single-mismatched RNA fragments and triplemismatched RNA fragments of both ORF1ab/RdRp/E. What is more, random RNA and random DNA were also detected, and the results were similar to those of blank groups, indicating a decent selectivity. Analytical Performance. After optimizing a series of reaction conditions, we explored the analytical performance of the established analytical method as listed in Figure 6a , the three targets can be mixed and quantified and have similar growth trends. In the inset figures of Figure 6b −d, after we perform logarithmic processing and linear fitting of the data, decent linear relationships were obtained at 20−4000 fmol for ORF1ab, RdRp, and E with equations y = 112,372 lg x − 11,479, y = 112,372 lg x − 11,479, and y = 29,967 lg x − 24,277. The limits of detection (3σ/k, LODs) for ORF1ab, RdRp, and E were 1.2, 1.3, and 1.3 fmol, respectively. For investigating the self-magnification of this strategy, the analytical performance was compared with those of other multicomponent simultaneous analysis methods by ICPMS. We use lanthanide NPs as element markers. Compared with the traditional macrocyclic compound-lanthanide element tagging, it has a more comprehensive linear range and three orders of magnitude higher sensitivity; compared with the DSN amplification strategy, the LOD was decreased by 40-fold (Table 3 ). This shows that even working with the simplest sandwich hybrid structure, our strategy can achieve higher sensitivity and high multiplicity without additional amplification reactions. Recoveries and Determined Results. Healthy human throat swab samples were collected, purified, collected, spiked, and analyzed to evaluate the method's applicability. To investigate the influence of the matrix, three target RNAs numbered as 1#, 2#, and 3# were spiked in RNA extraction solution simultaneously at the concentrations of 0, 200, and 2000 fmol. Moreover, the spiked samples were subjected to the proposed method. Quantitative recovery outcomes are all listed in Table 2 . As shown, 181 fmol of ORF1ab,189 fmol of RdRp, and 183 fmol of E were simultaneously found. The recoveries of 90, 94, and 92% were gained when 200 fmol RNA was added and extracted in the 1# throat swab sample. The spiked recovery results of other matrix samples at different concentrations are also within the range of 90−112%. Thus, it has added robust evidence for tolerance of the RNA virus extraction sample matrix and showed prospects of this strategy for the actual application of the experiments. In this work, we have established a multiplexed sandwich RNA analytical method based on NP self-enhanced RNA analysis to detect three specific RNA gene sites of SARS-CoV-2 simultaneously. Three lanthanide NP probes synthesized by the one-step solvent thermal method are modified by a carboxyl group and further coupled with DNA to perform signal report units. Meanwhile, with the chelating ability of β-NTA and the magnetic properties of SA-MBs, the lanthanide tag as the multiplexed probe can be quickly released in the form of ions to form a uniform solution. Thus, high-sensitivity detection can be realized through the self-amplification of NPs. If other nucleic acid amplification strategies are introduced into the system, the sensitivity can be further improved. We expect that this detection platform can increase the detection rate of SARS-CoV-2 through the simultaneous detection of multiple genes, and this strategy can be extended to the sensitive simultaneous detection of other nucleic acids or biomarkers. Proc. 2020 National Health Commission of the People's Republic of China and Traditional Chinese Medicine of the People's Republic of China (2020) Guidelines for the Diagnosis and Treatment of Coronavirus Disease CDC. CDC 2019-Novel Coronavirus (2019-nCoV) Real-Time RT-PCR Diagnostic Panel SOP# CDC-006-00019, 2020. (13) CDC. CDC 2019-nCoV Real-Time RT-PCR Diagnostic Panel Product Information Sheet SOP# CDC-006-00006 The authors declare no competing financial interest. We gratefully acknowledged the National Natural Science Foundation of China (no. 22074098 and 22074096), the Talents Program of Sichuan Province (no. 903), Sichuan Analytical Chemistry pubs.acs.org/ac Article pubs.acs.org/ac Article