Development and comparison of a novel multiple cross displacement amplification (MCDA) assay with other nucleic acid amplification methods for SARS-CoV-2 detection

MCDA target gene selection

To identify target genes with highly conserved regions and a suitable GC-content for MCDA, 1,216 SARS-CoV-2 genomes deposited in GISAID (all available complete, high coverage genomes (>29,000 bp) with low coverage flags excluded up until March 26, 2020)⁶ were aligned against the SARS-CoV-2 reference genome: NC_045512.2 using Snippy (v4.3.6) with the -ctgs flag and default settings (https://github.com/tseemann/snippy). A sliding window approach was then applied to identify conserved 300bp windows with GC content >43%, low SNP number, and low total SNP number (total SNPs was defined as the number of strains with a SNP in a given window). Three 300 bp conserved regions were identified and selected for MCDA primer design with two regions in ORF1ab (NC_045512.2: 515-831 and 12968-13288) and one in the N gene (NC_045512.2: 28345-28647).

MCDA primer design

For each region, 4 sets of MCDA primers were designed as previously described⁵. Each primer set consisted of 2 cross-primers (CP1/CP2), 2 displacement primers (F1/F2) and 6 amplification primers (C1/C2, D1/D2, R1/R2) (Supplementary Table 1). Non-specific primer binding was assessed using BLASTN against 14 non-SARS-CoV-2 coronaviruses used in Lamb et al.⁷, human genome (hg19) and 11 other common bacterial and viral respiratory pathogen/microbiome species.

Preparation of DNA/RNA standards

For each region, ∼500 bp gene fragments for ORF1ab (NC_045512.2: 416-931 and 12869-13388) and N gene (NC_045512.2: 28246-28747) were synthesised with an additional 100 bp up and downstream of the target region (ThermoFisher) (Supplementary Table 2). Each fragment contained a T7 promoter for transcription and M13 adapters for amplification. Seven DNA standards from synthesised gene fragments were prepared (100,000, 10,000, 5000, 1000, 500, 250 and 100 copies/µl).

For RNA, 1 pg of synthesised gene fragments were transcribed overnight at 37°C using T7 RNA polymerase (Sigma). Overnight DNA digestion was performed using the turbo DNA free kit (ThermoFisher) and further treated with DNase I (NEB) until all traces of DNA were removed. Complete DNA removal was confirmed after each round of DNase treatment using RT-PCR with the SensiFAST SYBR kit (Bioline) and F1/R1 MCDA primers (Supplementary Table 1). The transcribed RNA was serially diluted and used as input. Since copy number cannot be determined, lowest detectable dilutions were used for sensitivity comparison.

Initial evaluation of MCDA primer sets and optimisation of isothermal amplification temperature

MCDA reactions were performed using the WarmStart LAMP (DNA and RNA) kit (NEB) which contains a warmstart RTx reverse transcriptase and Bst2.0 polymerase for simultaneous reverse transcription and isothermal amplification. Antarctic thermolabile UDG was also added in each reaction to prevent carryover contamination.

For each primer set, a primer mix containing: 0.4 μM of F1 and F2, 0.8 μM of C1 and C2, 1.2 μM of R1, R2, D1 and D2 and 2.4 μM of CP1 and CP2 was used. Standard desalting purified primers were used for the initial evaluation and optimisation tests while HPLC purified primers were used for sensitivity and speed comparison against LAMP and RT-PCR tests.

For the initial evaluation of each MCDA primer set, a 10 µl reaction was used and contained: 5 µl of 2x WarmStart master mix, 0.2 µl of fluorescent dye, 1.2 µl of MCDA primer mix, 0.2 µl of 1U/µl Antarctic thermolabile UDG (NEB), 0.7 µl of 10 mM dUTP, 1.7 µl of H₂O and 1 µl of DNA template (1000 copies/µl). Specificity of MCDA primer sets were also evaluated using purified human genomic DNA (Sigma) and a microbial community DNA standard (Zymo Research). MCDA reactions were performed in triplicates in the Rotor-Gene Q (Qiagen) with isothermal amplification at either 60°C, 63°C or 65°C for 1 h and real time fluorescence detection every 60 seconds, followed by enzyme inactivation at 95°C for 5 min and a final melt curve from 50°C – 99°C to ensure correct MCDA product.

Comparison of MCDA, LAMP and RT-PCR

To compare the speed and sensitivity (limit of detection) of MCDA, LAMP and RT-PCR, published primers targeting the same SARS-CoV-2 MCDA ORF1ab (NC_045512.2: 416-931) and N (NC_045512.2: 28246-28747) regions were used (Supplementary Table 1). For LAMP, two published primer sets from Zhang et al.³ which targeted the same region as our MCDA were compared. For RT-PCR, there were no suitable published primers pairs which targeted the same ORF1ab region, therefore only primers submitted by the National Institute of Health, Thailand against the N gene was compared⁸. All primers were HPLC-purified grade.

MCDA, LAMP and RT-PCR were tested in three independent runs (biological replicates) using the same aliquot of DNA/RNA. Each run contained 3 technical replicates. The limit of detection was defined as the highest dilution where all 9 replicates (3 biological replicates x 3 technical replicates) were detected.

To reduce between run variations, 10 µl MCDA and LAMP reactions were set up and performed simultaneously in the same run. MCDA reactions were prepared as described above. For LAMP, each 10 µl reaction contained: 5 µl 2x WarmStart master mix (NEB), 0.2 µl fluorescent dye (NEB), 0.2 µl 1U/µl Antarctic thermolabile UDG (NEB), 0.7 µl 10 mM dUTP, 1 µl LAMP primer mix, 1.9 µl of H₂O and 1 µl of DNA/RNA template. Each LAMP primer mix contained 16 µM FIP and BIP, 2 µM F3 and B3 and 4 µM LF and LB. MCDA and LAMP isothermal amplification was performed at 65°C as described above.

For RT-PCR using DNA templates, 10 µl reactions containing 5 µl SensiFAST probe No-ROX mix (Bioline), 0.5 µl primer mix (40 µM F and R, 10 µM probe), 3.5 µl of H₂O and 1 µl DNA template were used. The cycling conditions were 95°C for 2 min, followed by 45 cycles of 95°C for 15 secs and 55°C for 30 secs.

For RT-PCR with RNA templates, 10 µl reactions were set up containing 5 µl SensiFAST probe No-ROX One-Step mix (Bioline), 0.5 µl primer mix (40 µM F and R, 10 µM probe), 0.1 µl reverse transcriptase (Bioline), 3.4 µl of H₂O and 1 µl RNA template. Reverse transcription was performed at 45°C for 20 minutes followed by RT-PCR amplification as described above for DNA.

To compare the speed of RT-PCR, cycle threshold (Ct) was converted to time using the following equation: Time = (Ct x 50 sec) + 120 sec. The detection time required for RT-PCR was calculated based on the cycling conditions (45 sec per cycle plus an initial 120 sec hold) and the ramp rate for the Rotor-gene Q (5 sec per cycle). The ramp rate for the Rotor gene Q is 15°C/s for heating and 20°C/s for cooling according to the manufacture’s technical information (https://www.qiagen.com/us/resources/download.aspx?id=2120af5e-8daf-4184-b277-aeb6ef5bbc05&lang=it-IT).

Development of MCDA assays for SARS-CoV-2 detection

Three 300 bp conserved regions suitable for MCDA primer design were identified from the genome alignment of 1,216 SARS-CoV-2 strains. Two regions, designated as region 1 and 2, belonged to the ORF1ab gene at NC_045512.2 position 515-831 and 12968-13288 respectively. One region, designated as region 3, corresponded to the N-gene at NC_045512.2 position 28345-28647. Four MCDA primer sets for each region was designed and evaluated (Supplementary Table 1).

Each MCDA primer set was initially tested at 3 isothermal amplification temperatures (60°C, 63°C and 65°C) using 1000 DNA copies/reaction as the starting template. As seen in Figure 1, regardless of the primer sets used, the slowest amplification time was observed at 60°C. Amplification at 63°C and 65°C were similar and 65°C was chosen as the isothermal amplification temperature used.

• Download figure
• Open in new tab

Figure 1:

Initial evaluation of MCDA primer sets at 3 different isothermal amplification temperature (60°C, 63°C and 65°C). Four MCDA primer sets were designed for each target region chosen. A-C: Region 1 ORF1ab: 515-831 D-E: Region 2 ORF1ab: 12968-13288 F-H: Region 3 N gene:28345-28647.

To maintain MCDA assay robustness against SNPs which may affect MCDA primer binding and amplification efficiency, primer sets from two different regions were chosen for further development as a duplex assay. Amplification of region 2 was the slowest for all primer sets (Figure 1 D-F) compared with region 1 and 3, taking between 15-25 minutes at 65°C. Primer sets in region 2 also had very high variation between technical replicates (data not shown). Therefore region 2 was removed from further evaluation.

Region 3 amplification of the N gene was the fastest with primer set 2 followed closely by primer set 3 (Figure 1G-I). Primer set 1 and 4 were the slowest for region 3 and were therefore eliminated from further testing. We also observed that primer set 2 had tighter technical replicates compared to primer set 3 (data not shown), thus region 3 primer set 2 was chosen as our final MCDA primer set for further sensitivity and specificity testing.

Within region 1, primer set 3 was the slowest with fluorescence appearing at ∼35 min (Figure 1C). This primer set was removed from further consideration. Primer set 1 was the fastest primers to amplify region 1 and was chosen for inclusion in our MCDA assay.

Therefore, the final primer sets chosen for MCDA SARS-CoV-2 detection was region 1 (ORF1ab) primer set 1 and region 3 (N gene) primer set 2 (Figure 2). Both primer sets showed no non-specific amplification when tested against human and microbial community genomic DNA.

• Download figure
• Open in new tab

Figure 2:

The nucleotide sequences and position of the final (A) ORF1ab and (B) N gene COVID-19 MCDA primer sets chosen in this study. Right and left arrows show sense and complementary sequences, respectively while coloured text indicate the position of primers: F1/F2 displacement primers in purple, P1/P2 primers in green, C1/C2 amplification primers in blue, D1/D2 amplification primers in yellow and R1/R2 amplification primers in red.

Sensitivity and time to detection comparison of MCDA, LAMP and RT-PCR

The sensitivity and speed for MCDA, LAMP and RT-PCR were then compared for two SARS-CoV-2 genes. For the N gene (region 3), detection by MCDA was consistently faster than LAMP, by ∼10-13 minutes, for most DNA dilutions tested (Table 1). The average fastest detection time for MCDA was 5.2 minutes at 100,000 copies/reaction while for LAMP it was 15 minutes. MCDA also had higher sensitivity with the limit of detection at 1000 copies/reaction while for LAMP it was 5000 copies/reaction. This limit of detection was equivalent to a RT-PCR Ct value of 32.4 and 30.3. respectively. A greater number of sporadic detections for higher dilutions were also observed for MCDA compared to LAMP. MCDA was also significantly faster than RT-PCR by ∼20 minutes for lower dilutions and 10 minutes for higher dilutions. At 100,000 copies/reaction, the detection time for RT-PCR was 23 minutes. However, RT-PCR had the highest sensitivity with the limit of detection at 100 copies/reaction. Similar results were also observed using RNA template with RT-PCR being the most sensitive technique, detecting RNA at 10⁻⁶ dilution. For MCDA, the lowest RNA dilution detected was 10⁻⁴, with only sporadic detection for LAMP at this dilution.

For ORF1ab, LAMP was more sensitive than MCDA with the limit of detection at 500 copies/reaction and had more sporadic detection at higher dilutions (Table 2). The limit of detection for MCDA was 5000 copies/reaction. LAMP was also slightly faster than MCDA by ∼1-3 minutes for concentrations above the limit of detection. The fastest time to detection for ORF1ab LAMP was 8.4 minutes while for MCDA it was 9.7 minutes.