Subsection
The 2019 SARS-CoV-2 database (RCoV19) from the National Genomics Data Center (NGDC) contains data for comparison of SARS-CoV-2 subtypes. We obtained the mutation site of each subtype and its mutation frequency. Through screening, we selected mutation sites with mutation frequencies above 0.95 and in close proximity as our screening targets. All mutation sites we used are shown in Fig. 1. The yellow blocks indicate the mutation sites present on the ARMS primers, which were used as our initial screening targets to specifically select subtypes from the majority of the SARS-CoV-2 viruses, while the other blocks indicate the mutations present on the probes of different subtypes, which are able to distinguish between the subtypes that share mutations on the primers.
Plasmid Construction
The sequences of all isoforms were obtained from the open-source database of the New Crown Pneumonia isoform sequences provided by NCBI (www.ncbi.nlm.nih.gov) and GISAID (www.gisaid.org). The sequence numbers are shown in Table 1. After comparison of the published sequences through SnapGene, all sequences contained loci with high mutation rates in agreement with those in RCoV19. We cloned the S gene sequences of all isoforms into the vector pcDNA3.1. Plasmids that served as templates for subsequent experiments were purchased from General Biol (Anhui) Co., Ltd .
The plasmids were centrifuged at 10,000 g for 2 min, and then TE Buffer (Cat. No: BL531A, Biosharp, Beijing, China) was added. The plasmid DNA concentration was measured with a NanoDrop One system (Thermo Scientific) and then diluted with TE buffer. The concentration was confirmed with NanoDrop. These steps were repeated until the concentration of all plasmids was 10 ng/µL. The calibrated plasmids were subsequently diluted 1000× (10 pg/µL) and 10,000× (1 pg/µL) for experiments.
Primer And Probe Design
The sequences of all primers and probes designed using Primer Premier 6, Beacon Designer 8, and SnapGene are shown in Table 1. We aimed to preserve the mutation sites of each isoform while following the basic principles of primer design, so that the mutation sites on the primers are as close to the 3′ end as possible and the mutation sites on the probes are located in the middle of the sequence.
All primers were 18 to 22 nucleotides in length with Tm values of 45–55°C, upstream and downstream consistency as high as possible, a GC content of 40–60%, G/C within the last five bases, no four or more replicates of any base, and an amplicon size of 100–200 bp. The specificity of SARS-CoV‐2 and other organisms was examined using BLAST in NCBI. Primers were purchased from Zhejiang Sunya Biotechnology Co., Ltd.
To enhance the specificity of stem sequence recognition, probe size is 30–40 bp, with a loop of usually 20–30 bp in length which is complementary to the target sequence and a stem of usually 5–7 bp in length with high GC content. They are paired with each other to form the stem structure. The fluorescent groups are labeled at one end of the probe, while quenchers are labeled at the other end. A 5′ FAM fluorescent group and a 3′ Dabcyl quenching group were used to adapt to our molecular hybridization method. Other fluorescent groups can be used according to the research needs. Probes were purchased from General Biol (Anhui) Co., Ltd .
TABLE 1 Template, primer and probe sequences for COVID-19 subtypes
The letters in bold represent hairpins of molecular beacons and the red letters represent point mutated nucleotides A wide center space indicates a missing position.
Pcr And Gel Electrophoresis
To test the amplification using ARMS-PCR primers, we used a Bio-Rad T 100 thermocycler for PCR amplification. Each reaction system contained 4 µL 5× Taq Pro Buffer (#PN101, Vazyme Biotech Co. Ltd., Beijing, China), 0.5 µL forward primer (10 µM), 0.5 µL reverse primer (10 µM), 1 µL template (10 pg/µl), 0.2 µL Taq Pro HS DNA polymerase (Vazyme Biotech Co. Ltd., Beijing, China), and 12.8 µL ddH2O.
The PCR program was as follows: 95°C for 3 min, followed by 30 cycles of 95°C for 10 s and 50°C for 15 s. The PCR products were analyzed using 1.5% agarose electrophoresis and visualized using an Azure Imaging Systems 280 (USA).
Primer Optimization
To optimize primer specificity, we mismatched the base pair at the penultimate or third position at the 3′ end and then detected the specificity of individual primer pairs using a real-time fluorescence qPCR instrument (Bio-Rad CFX Connect) and Vazyme ChamQ SYBR qPCR Master Mix (without Rox). Each reaction system contained 10 µL 2× ChamQ SYBR qPCR Master Mix (Vazyme Biotech Co. Ltd., Beijing, China), 0.4 µL forward primer F (10 µM), 0.4 µL reverse primer (10 µM), 1 µL template (1 pg/µL), and 8.2 µL ddH2O. The PCR program was as follows: 95°C for 3 min, followed by 40 cycles of 95°C for 10 s and 50°C for 15 s. The fluorescence signal was collected during the reaction. Finally, the primers with non-specific samples and no amplification or large Ct value differences were selected for subsequent experiments, and for the primers that did not meet the requirements, the specificity was enhanced by adjusting the sequence or introducing a new mismatch.
Probe Specificity Detection
The probe assay was used to detect samples of different templates by matching the screened primers with the corresponding probes to confirm whether the corresponding probes can specifically distinguish between different target mutants that have the same primers.
qPCR amplifications were carried out in a 20 µL reaction volume containing 4 µL 5× Taq Pro Buffer (Vazyme Biotech Co. Ltd., Beijing, China), 1 µL (1 pg/µL) DNA template, 0.4 µL Taq DNA polymerase (5 U/µL, Vazyme Biotech Co. Ltd., Beijing, China), 0.4 µL (10 µM) each of the forward and reverse primers, 0.2 µL Beacon (10 µM), and 13.6 µL ddH2O.
For real-time fluorescence PCR detection using Vazyme Taq Pro HS DNA polymerase reagent in a Bio-Rad CFX Connect system, the PCR program was as follows: 95°C for 3 min, followed by 40 cycles of 95°C for 10 s and 55°C for 15 s. The fluorescence signal was collected during the reaction. The specificity of the probe was demonstrated by the fact that the amplification curve was only detected when using the corresponding template.
Molecular Hybridization
Molecular hybridization is a method for detecting specific products at time nodes after amplification. After amplification, a fluorescent probe was added to detect the background signal value. Due to the presence of the molecular beacon probe stem–loop structure, the probe did not produce fluorescence in the normal state. After denaturation at 95°C, the probe specifically combined with the amplification product at 55°C, and the emitted fluorescence signal was captured using a fluorescence receiver. The signal value at this time was multiplied by the difference from the background.
First, a preliminary amplification was carried out on that fragment. In this case, we formed multiplex PCR by combining multiple primer pairs, such as A570D-Primer + G339D-Primer + P681R-Primer or T547K-primer + LPPA24-27S-primer. The feasibility of this method was tested using a multiplex PCR protocol. Each reaction system contained 10 µL 5× Taq Pro Buffer (Vazyme Biotech Co. Ltd., Beijing, China), 2.5 µL forward primer (10 µM), 2.5 µL reverse primer (10 µM), 10 µL template (1 pg/µL), 1 µL Taq Pro HS DNA polymerase (Vazyme Biotech Co. Ltd., Beijing, China), and (29 − 5 × N) µL ddH2O, where N = 1–3. The PCR program was as follows: 95°C for 3 min, followed by 35 cycles of 95°C for 10 s and 50°C for 15 s.
The amplification product was added to a black 384-well plate (#w/lid Cell Culture Sterile PS, Thermo Fisher Scientific) with 1 µL of beacon probe (10 µM). The probe with the FAM fluorophore was then detected using a SpectraMax iD5 multi-mode microplate reader at 55°C (excitation at 484 nm, emission at 525 nm). The average value from three replicates was taken as the background, i.e., the fluorescence signal of the reagent itself when the probe was not attached to the template.
After denaturation at 95°C for 3 min, molecular hybridization of the probe to the amplified product was performed at 55°C for 30 s. At this point the products were read again on a SpectraMax iD5 multi-mode microplate reader. The average value from the three repeats was taken. The signal change at this point reflects the detection of the template by the probe. Specific templates will produce a more pronounced increase in signal.
Sensitivity Test
To verify the detection limit of this assay, both RT-PCR and molecular hybridization were used. The real-time fluorescence quantification system and the molecular hybridization system performed as expected based on previous tests. First, the plasmid copy number of 9100 bp at 1 ng/µL = (6.02 × 1023 × 10− 6)/(9100 × 660) = 1011 copies/mL. Real-time fluorescence qPCR was performed using concentrations of 1012 copies/mL, 1010 copies/mL, 108 copies/mL, 106 copies/mL, 104 copies/mL, and 102 copies/mL. Plasmid positive controls were used, and deionized water was used as the negative control as the detection limit reference. Molecular hybridization was performed using concentrations of 1010 copies/mL, 108 copies/mL, 106 copies/mL, 104 copies/mL, and 102 copies/mL, with deionized water as a negative control as the assay-limiting reference.