Structure-Based Primer Design Minimizes the Risk of PCR Failure Caused by SARS-CoV-2 Mutations

The coronavirus disease 2019 (COVID-19) has caused and is still causing tremendous damage to the global economy and human health. Qualitative reverse transcription-PCR (RT-qPCR) is the golden standard for COVID-19 test. However, the SARS-CoV-2 variants may not only make vaccine less effective but also evade RT-qPCR test. Here we suggest an innovative primer design strategy for the RT-qPCR test of SARS-CoV-2. The principle is that the primers should be designed based on both the nucleic acid sequence and the structure of the protein encoded. The three nucleotides closest to the 3′ end of the primer should be the codon which encodes the tryptophan in the structure core. Based on this principle, we designed a pair of primers targeting the nucleocapsid (N) gene. Since tryptophan is encoded by only one codon, any mutation that occurs at this position would change the amino acid residue, resulting in an unstable N protein. This means that this kind of SARS-CoV-2 variant could not survive. In addition, both our data and previous reports all indicate that the mutations occurring at other places in the primers do not significantly affect the RT-qPCR result. Consequently, no SARS-CoV-2 variant can escape detection by the RT-qPCR kit containing the primers designed based on our strategy.

Introduction

Since its outbreak in December 2019 in Wuhan, China, the coronavirus disease-2019 (COVID-19) soon developed into a global pandemic (Liu R. et al., 2020; Wang et al., 2020a). As of July 14, 2021, more than 187.09 million COVID-19 cases and 4.04 million deaths had been confirmed, giving an overall mortality rate of 2.16%. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified to be the causative agent. Compared to its sister virus SARS-CoV which caused the 2003 SARS epidemic (8,422 cases), SARS-CoV-2 is much more contagious (Liu Y. et al., 2020; Wu et al., 2020). Its R0 value is 2.3, but it could be as high as 5.78 (Bulut and Kato, 2020). Rapid diagnosis and isolation are the first and most effective steps to inhibit the spread of SARS-CoV-2 before the vaccines, and therapeutic drugs are sufficiently effective to control the pandemic (Wang et al., 2020d; Younes et al., 2020). At present, the diagnostic methods for COVID-19 mainly include nucleic acid detection (RT-qPCR), antigen testing, and imaging technology (Jin et al., 2020). Because RT-qPCR has the characteristics of short window period, high sensitivity, and high specificity, it has become the golden standard for COVID-19 test in various institutions around the world (Liu R. et al., 2020; Udugama et al., 2020; Wang et al., 2020b).

It is noteworthy that SARS-CoV-2, as an RNA virus, utilizes RNA-dependent RNA polymerase (RdRp) for the replication and transcription of its genome (Gao et al., 2020). Despite the existence of a “proofreading” mechanism based on nsp14 exonuclease (ExoN) to remove mistakenly incorporated nucleotides, variant strains are still constantly accumulating (Lippi et al., 2020; Wang et al., 2020c). According to the GISAID database, thousands of single-nucleotide variants have been identified across different SARS-CoV-2 strains isolated. Four main types of variants—D614G, Cluster 5, VOC 2020/12/01, and 501Y.V2—have acquired wide attention. The current epidemiological studies have shown that, although these variants do not increase the severity of the disease (assessed by hospitalization time and 28-day mortality) and reinfection, they lead to higher morbidity and hospitalization and more death cases (Abdullahi et al., 2020; Davies et al., 2020; Korber et al., 2020; Tegally et al., 2020; Oude Munnink et al., 2021). It also raises a concern that, as more and more mutations accumulate, some variants may evade the RT-qPCR test, leading to more and more false-negative results (Lippi et al., 2020; Wang et al., 2020b; Rahman et al., 2021). A previous study showed that using degenerate primers in PCR is helpful to detect the variants (Li et al., 2015; Campos and Quesada, 2017; Li et al., 2020), but using degenerate primers brings high cost and low amplification efficiency, which prevents its wide use, especially in cases of detecting positive critical samples (Souvenir et al., 2007; Maher-Sturgess et al., 2008). Therefore, it is an urgent task to find a way to minimize the risk of PCR failure caused by SARS-CoV-2 mutations.

A phylogenetic analysis of the genome sequence of SARS-CoV-2 revealed that ORF1ab, E, and N genes are highly conserved in sarbecoviruses and have been used as target genes by Centers for Disease Control and Prevention (CDCs) in various countries for RT-qPCR test (Guan et al., 2020; Ramirez et al., 2020; Rahman et al., 2021). By far, all these target genes have undergone mutations, with the N gene having the greatest number of mutations on the targets of primers and probes which have been widely used around the world to diagnose COVID-19 (Wang et al., 2020d). Here we systematically analyze the effect of the mismatches between the target sequences and primers or probes on PCR test. Consistent with previous reports, the mismatches occurring at the 3′ end of the primer have the greatest impact on the PCR reaction (Bru et al., 2008; Stadhouders et al., 2010). Therefore, placing the 3′ end of the primers at the site where any mutation is lethal to the virus would greatly reduce the risk of PCR failure caused by SARS-CoV-2 mutations. It has been reported that tryptophan, as the largest amino acid, has a special bulky side chain and plays an important role in protein structure stability (Santiveri and Jiménez, 2010; Bielecki et al., 2014). Among the 20 amino acids that make up proteins, only tryptophan and methionine are encoded by a non-degenerate codon. Combining this information, we suggest an innovative structure-based primer design strategy for the RT-qPCR test of SARS-CoV-2. Experiments showed that the primers designed based on this strategy have the same specificity and sensitivity as the primers designed by US and Chinese CDCs. Mutations affecting the RT-qPCR test always lead to the production of an unstable nucleocapsid, which means that these kinds of SARS-CoV-2 variants cannot survive.

Materials and Methods

Bacterial Strains, Plasmids, and Culture Conditions

The N genes of seven coronaviruses (SARS-CoV, MERS-CoV, HCoV229E, HCoV-OC43, HCoV-NL63, HCoV-HKU1, and SARS-CoV-2) that infect humans were synthesized in Beijing Genomics Institute and then inserted into a pUC57 vector, respectively. The N gene of SARS-CoV-2 was cloned into a PET28a vector. All these plasmids were then transformed into Escherichia coli DH5α cells for amplification. Plasmids containing the mutant N protein gene were constructed using the improved QuikChange method (Xia and Xun, 2017). The partially overlapping primers were used to generate the PCR products for site-directed mutagenesis (SDM). Then, Fast-Digest DpnI was added into the SDM–PCR reaction mix directly without purification. The mixture was incubated at 37°C for 1–3 h to completely remove the residual template plasmid. Furthermore, 5–10 μl DpnI-treated SDM–PCR product was added into the competent cells for DNA transformation. The positive clones were identified by DNA sequencing. Plasmid DNA was extracted by using Plasmid Mini Kit (Omega) and saved at -20°C for future use. For protein expression, plasmids containing the SARS-CoV-2 N gene and mutants were transformed into E. coli BL21 (DE3).

Preparation of the RT-qPCR Template

The plasmids containing the N genes of the other six types of coronaviruses and human genome were prepared. These plasmids and N gene pseudovirus were added to the throat swab from healthy human and placed into the preservative solution (25 mM Tris–HCl, pH 7.6, 1 mM EDTA, 20 mM guanidine thiocyanate) to simulate inactivating clinically positive samples. Then, 200 μl of preservation solution was taken, and the RNA/DNA was extracted using MiniBEST Viral RNA/DNA Extraction Kit (TaKaRa). The DNA template is diluted with enzyme-free water to 10⁷copies/ml, and the RNA is diluted with enzyme-free water to 10⁷, 10⁶, 10⁵, 10⁴, 10³, 700, 500, 400, and 300 copies/ml.

Effect of Mismatch Between Primers and Template on RT-qPCR Reaction

Three pairs of primers designed by the US CDC were chosen as the research object, with the PUC57-N plasmid as the template. A series of mutations was designed on the primers to create mismatches with the template DNA. The mutations change the first, the second, and the third nucleotide of the 3′ end of the primers and one nucleotide in the middle of the primers, respectively. At each position the nucleotide was mutated into another three mismatched nucleotides (for example, A was mutated to T, C, and G). Then, the Ct values for the RT-qPCR reactions containing mutated primers and original primers were compared when the template concentration is 10⁷ copies/ml.

Effect of Mismatch Between Probes and Template on RT-qPCR Reaction

In TaqMan RT-qPCR reaction, fluorescence is produced when Taq DNA Polymerase cuts off the fluorescent group-modified nucleotide at the 5’ end of the probe. Mutations were designed on the template DNA to make mismatches between the first nucleotide at the 5’ end and in the middle of the probe, respectively. At each position, the nucleotide was mutated into another three mismatched nucleotides (for example, A was mutated to T, C, and G). Then, the Ct values for the RT-qPCR reactions containing the mismatched DNA template and probe and the reaction containing the original DNA template and probe were compared, with the template concentration at 10⁷ copies/ml.

Analysis of the Interactions Between a Specific Amino Acid Residue and Surroundings in N Protein

The interactions between a specific amino acid site (W108, W132, W301, M322) and surroundings was analyzed to evaluate their role in stabilizing protein structure. The molecular graphics figures were generated using PyMOL (http://www.pymol.org).

Evaluation of the Sensitivity of Structure-Based Primers

To evaluate the sensitivity of the primers designed based on structure (NAm1 and NAm2), Ct values for the PCR reactions containing our primers and the N gene targeting primers designed by US and Chinese CDCs at four RNA template concentrations (10⁷, 10⁶, 10⁵, and 10⁴ copies/ml) were compared. All the RT-qPCR reactions contain a pair of primers at 400 nM and a probe at 200 nM. The RT-qPCR assays were performed on Analytik jena qTOWER3.

Evaluation of the Specificity of Structure-Based Primers

To date, a total of seven types of coronaviruses have been identified as human pathogens (Zhu et al., 2020). Well-designed primers should be able to distinguish between different coronaviruses, especially between the highly pathogenic and the less pathogenic viruses. For specificity evaluation, the plasmids containing the N genes of the other six types of coronaviruses and human genome were prepared and used as templates at 10⁷ copies/ml in RT-qPCR reactions. The Ct values for RT-qPCR reactions containing different combinations of primers and templates were measured and compared. The data are obtained from triplicate experiments.

Determination of the Limit of Detection

Under the same reaction system (different primer pairs), RT-qPCR was performed with five template concentrations of 10³, 700, 500, 400, and 300 copies/ml, respectively, and each concentration was repeated 20 times.

Result Judgment

When Ct ≤38, it is judged as a positive sample; Ct >38 or no Ct is judged as a negative sample. The template concentration with a detection rate greater than 95% (19 times) was defined as the limit of detection of the primer pair.

Protein Expression and Purification

For the protein expression of N-terminal domain (NTD), C-terminal domain (CTD), and their mutants, E. coli BL21(DE3) cells containing expression plasmids were grown in LB medium supplemented with 50 μg/ml kanamycin at 37°C, 200 rpm. When the OD600 reached 0.8, the temperature was lowered to 16°C, and then a final concentration of 0.5 mM isopropyl β-D-thiogalactopyranoside (IPTG) was added for overnight induction. For the expression of the N protein and its mutants, when the OD600 reached 0.8, IPTG at a final concentration of 0.5 mM was added to induce expression for 3 h.

For protein purification, cells were harvested by centrifugation at 5,000 g for 20 min. Cell pellet was resuspended in lysis buffer (25 mM Tris–HCl, pH 8.0, 500 mM NaCl) and then lysed by sonication on ice. After centrifugation at 28,370 g for 50 min at 4°C, the supernatant was loaded onto a ““nickel chelating sepharose affinity column (GE Healthcare) equilibrated with lysis buffer in advance. The column was then washed with wash buffer (25 mM Tris–HCl, pH 8.0, 500 mM NaCl, 50 mM imidazole) and then eluted with elution buffer (25 mM Tris–HCl, pH 8.0, 500 mM NaCl, 250 mM imidazole). The eluted sample was further purified by size-exclusion chromatography using Superdex 200 (GE Healthcare) in 25 mM Tris–HCl, pH 8.0, and 500 mM NaCl. Finally, SDS-PAGE was used to assess protein purity.

Determination of Protein Thermal Stability

The effects of mutations on specific amino acid residues were assessed by measuring the thermal stability of the wild-type and mutant NTDs and CTDs. For this purpose, the measurement was performed by using Protein Thermal Shift™ Kit on Applied Biosystems QuantStudio 3. Then, 20 μl of protein sample at 0.4 µg/µl was heated from 25 to 99°C at a rate of 0.15°C/s. Data collected from triplicate experiments was used to calculate T_m values using Protein Thermal Shift™ Software 1.4 by fitting the data in the region of analysis to the Boltzmann equation.

Statistics Analysis

All data were collected in at least three independent experiments and were presented as mean ± standard deviation of the results of triplicate experiments. The Student’s t-test was used to perform the statistical analyses. Statistical significance was assessed based on the p-value: *p < 0.05, **p < 0.01, and ***p < 0.001.

Results

The Mismatches That Occur at the 3′ End of the Primer Have the Greatest Impact on PCR Reaction

To find the best strategy for designing primers which make a PCR reaction minimally affected by mutations, the effect of mutations at different positions has been extensively evaluated. The impact of the mutations is indicated by the difference of the cycle threshold (Ct) value with or without mismatch. Our data showed that, for all the six primers tested, when the mismatch occurred at the 3′ end the Ct value difference (ΔCt) reached the maximum (Figure 1). In the worst cases, for example, when the first nucleotide A at the 3′ end of the N3 forward primer was replaced by C or G, the PCR reaction completely failed (Figure 1F). This means that if the virus had a mutation at this position, the RT-qPCR test would always give a false-negative result even if the viral load in the sample is high. The ΔCt decreased as the mismatch moves away from the 3′ end of the primers (Figure 1). When the mismatch occurred in the middle of the primer, ΔCt is within ±1, indicating that the PCR test is basically unaffected.

FIGURE 1

Figure 1 Effect of mismatch between primers and template on PCR reaction. (A, B) Ct values detected when the US N1 forward primer (A) and reverse primer (B) contained mismatch with the template. (C, D) Ct values detected when the US N2 forward primer (C) and reverse primer (D) contained mismatch with the template. (E, F) Ct values detected when the US N3 forward primer (E) and reverse primer (F) contained mismatch with the template. The purple square, blue dot, cyan upper triangle, and magenta lower triangle represent the first, second, third, and 10th (middle position) base of the 3′ end of the primer, respectively. Student’s t-test by SPSS 15 was used for data analysis. A p-value < 0.05 was considered statistically significant. *p < 0.05; **p < 0.01; ***p < 0.001.

By using a similar procedure, we also evaluated the impact of the mismatch occurring at different positions on the probes. The exciting thing is that the PCR reaction is not sensitive to the mismatch on the probes (ΔCt is within ±1) (Figure 2). This means that a strategy to design primers for a robust PCR reaction not sensitive to virus mutation is to place the 3′ end of the primers at the positions where any mutation is fatal to the virus. In this way, the mutations that significantly affect the RT-qPCR test will never happen.

FIGURE 2

Figure 2 Effect of mismatches between probes and template on PCR reaction. (A) Ct values detected when the US N1 probe contained mismatch with the template. (B) Ct values detected when the US N2 probe contained mismatch with the template. (C) Ct values detected when the US N3 probe contained mismatch with the template. Student’s t-test by SPSS 15 was used for data analysis. A p-value < 0.05 was considered statistically significant. *p < 0.05; **p < 0.01.

In order to achieve this goal, the primers should be designed based on both the nucleic acid sequence and the structure of the protein encoded. The three nucleotides closest to the 3′ end of the primer should be the codon which encodes tryptophan in the structure core. Tryptophan is encoded by only one codon, so any mutation at that location alters the amino acid residue. Since tryptophan plays an important role in stabilizing protein structure (Ohmae et al., 2001), the mutation of tryptophan will lead to the destruction of protein structure stability. This means that this kind of SARS-CoV-2 variant could not survive.

Primers Targeting the N Gene Were Designed Based on Both Gene Sequence and Protein Structure

N gene is the major target of the RT-qPCR test for COVID-19 diagnosis. It encodes the N protein, which consists of a NTD and a CTD, and forms a dimer in solution through the dimerization of CTD (Peng et al., 2020). NTD contains two tryptophan residues, W108 and W132, which are located right in the structure core. W108 forms a hydrogen bond with Q58 and L56 and forms a hydrophobic interaction with A55, Y87, R88, R89, R107, G129, and F171, respectively (Figure 3A). W132 forms a hydrogen bond with A125 and forms a hydrophobic interaction with G85, G124, Y123, Y86, L113, L121, and P122, respectively (Figure 3B). A pair of primers were designed, with their 3′ end located on the codon of the two tryptophans, respectively. CTD contains only one tryptophan residue W301. Fortunately, there is M322 which is at the right distance from W301. Methionine also has only one codon, so we choose it as the suboptimal choice and design the second pair of primers. Both W301 and M322 are located on the dimer interface. W301 forms a hydrogen bond with Y298, I304, and A305 and forms a hydrophobic interaction with K299, Q303, Y296, and I291 in the same protomer and with A311 and S312 in the other protomer, respectively (Figure 3C). M322 forms a hydrophobic interaction with W330, L331, T329, and l353 in the same protomer and with L353 in the other protomer, respectively (Figure 3D). Based on the structural information of the N protein, we selected W108, W132, W301, and M322 for primer design because these residue pairs are not only critical for the stability of the protein but also of the right distance on the sequence for a quick PCR reaction. We expect that any mutation that changes one of these amino acid residues would result in an inactive N protein.

FIGURE 3

Figure 3 W108, W132, W301, and M322 are in the structure core of the N protein. (A) The interaction diagram of W108 in the N protein structure. (B) The interaction diagram of W132 in the N protein structure. (C) The interaction diagram of W301 in the N protein structure. (D) The interaction diagram of M322 in the N protein structure.

Effect of Mutations on the Thermal Stability of N Protein

Single-base mutations are much more likely to occur than double- or triple-base mutations during virus replication, for example, the famous P4715L in ORF1ab (nucleotide 14,143, C to T) and D614G in S (nucleotide 23,403, A to G) all resulted from a single-base mutation. The probability of double- or triple-base mutations on the codon of the same amino acid is very small, so we tried to investigate if all the single-nucleotide mutations in the codon of W108, W132, W301, and M322 would produce an inactive N protein. These mutations totally produce 21 mutants of the N protein domain (except for the stop codon). The thermal transition midpoint temperature (T_m) of mutant NTDs and CTDs were measured and compared with the corresponding wild-type domains. Compared with the wild-type NTD, the T_m values of W108C, W108S, W108L, W108R, W108G, W132C, W132S, W132L, W132R, and W132G mutant domains decreased by 15.02 ± 0.09, 16.04 ± 0.06, 12.8 ± 0.18, 17.82 ± 0.02, 22.19 ± 0.19, 14.48 ± 0.08, 14.25 ± 0.08, 18.65 ± 0.02, 19.74 ± 0.28, and 16.05 ± 0.02°C, respectively (Figure 4A). Compared with the wild-type CTD, the T_m values of W301C, W301S, W301L, W301R, W301G, M322V, M322K, M322R, M322T, and M322I mutant domains decreased by 27.76 ± 0.04, 28.02 ± 0.05, 30.13 ± 0.11, 27.74 ± 0.17, 29.42 ± 0.16, 2.955 ± 0.01, 7.095 ± 0.10, 8.82 ± 0.05, 3.79 ± 0.05, and 1.38 ± 0.01°C, respectively. Unexpectedly, the T_m value of M322L mutant domain increased by 1.085°C (Figure 4B).

FIGURE 4

Figure 4 Point mutations affect the thermal stability of the corresponding domains. (A) Thermal transition midpoint temperature difference of the mutant domains at positions W108 and W132. (B) Thermal transition midpoint temperature difference of the mutant domains at positions W301 and M322. Student’s t-test was used to conduct statistical analysis, and differences were considered significant when p <0.05, **p < 0.01, and ***p < 0.001.

All the mutations at W108, W132, and W301 and most mutations at M322 caused a significant decrease of the thermal stability of the NTD or CTD. It was indicated that all these three tryptophan residues are critical for the thermal stability of N protein. Any mutation occurring on any of these tryptophan would decrease the T_m value to below 37°C, which means that these mutants would have no normal functions inside the human body.

The N Protein Mutants Are Poorly Expressed at 37°C

Since the mutations on W108, W132, and W301 cause the T_m values of the corresponding domain to be lower than 37°C, we tested if the full-length protein containing these mutations could be expressed normally. The effect of these mutations on the N protein was evaluated by comparing their expression and solubility with the wild-type protein. The results showed that all mutations on W108, W132, and W301 caused a significant decrease in protein expression level and solubility. The mutant W301 has the lowest solubility (Figures 5A, B). By contrast, the mutations on M322 do not affect the protein expression level (Figure 5B). In addition, all these mutant proteins have a poorer behavior in comparison with the wild-type protein. These results, combined with the data of thermal stability of the mutant domains, indicate that the sites that we selected are critical to the function of the N protein. Therefore, SARS-CoV-2 variants containing these kinds of mutations would be outcompeted by other strains inside the human body and have no chance to influence the RT-qPCR test.

FIGURE 5

Figure 5 The yield of N protein mutants expressed at 37°C. (A) The yield of the mutants at positions 108 and 132. (B) The yield of the mutants at positions 301 and 322. Statistical significance is indicated as compared with wild type using t-test. *p < 0.05; **p < 0.01.

Performance of the Structure-Based Primers

We have shown that if a PCR reaction contains structure-based primers it will be very difficult for a virus to evade the RT-qPCR test by mutation. Now we need to know if the primers designed by this strategy perform as well as the primers designed by both US and Chinese CDCs. Based on the structure of the N protein, we designed two pairs of primers, N Anti-mutation 1 (NAm1) and N Anti-mutation 2 (NAm2), targeting NTD and CTD, respectively. The 3′ end of the forward primer of NAm1 is located on the nucleotide sequence encoding W108, and the 3′ end of the reverse primer is located on the nucleotide sequence encoding W132. For NAm2, the 3′ end of the forward and reverse primers are on the nucleotide sequences encoding W301and M322, respectively. The sequences and T_m values of the primers and probes are shown in Table 1.

TABLE 1

Table 1 The sequences of N anti-mutation (NAm) primers and probes.

The sensitivity of the primers and probes were evaluated by PCR reaction containing a certain amount of the template. The Ct values obtained for NAm1 and NAm2 were compared with the average Ct values obtained for four pairs of primers US N1, US N2, US N3, and CHN N. The results showed that the differences of the Ct values were all within ±1 at four template concentrations (Figure 6A). Both NAm1 and NAm2 meet the limit of detection (500 copies/ml) of the primers designed by the US and Chinese CDCs (Table 2). This means that NAm1 and NAm2 have the same sensitivity as the primers designed by the US and Chinese CDCs.

FIGURE 6

Figure 6 (A) Sensitivity test of N anti-mutation (NAm) primers. The black, green, blue, cyan, orange, and purple curves represent US N1, N2, N3, CHN N, NAm1, and NAm2 primers, respectively. (B) Specificity test of NAm primers. The blue and purple curves represent the amplification curves of NAm1 and NAm2 when the SARS N gene was used as the template, and the flat lines represent the amplification curves when the MERS-CoV, HCoV-229E, HCoV-OC43, HCoV-NL63, and HCoV-HKU1 N gene and human genome were used as templates.

TABLE 2

Table 2 Positive rate for primer pairs at different concentrations.

The specificity of the primers and probes were evaluated by PCR reaction containing a certain amount of the N gene of SARS-CoV-2, SARS-CoV, MERS-CoV, HcoV-229E, HCoV-OC43, HCoV-NL63, and HCoV-HKU1 and human genome. In the PCR assay, when the concentration of the SARS-CoV N gene template was 10⁷ copies/ml, the Ct values for NAm1 and NAm2 were 24.82 and 22.78, respectively. However, when the same amount of MERS-CoV, HCoV-229E, HCoV-OC43, HCoV-NL63, and HCoV-HKU1 N genes and human genome was used as template, no amplification occurs (Figure 6B). Since the SARS virus has not been identified in the world since 2003, we thought that the specificity of NAm1 and NAm2 is both excellent.

Discussion

The pandemic of COVID-19 is still raging in many countries. It not only makes great damages to many developed countries but also makes catastrophe in developing countries like India and other countries in south Asia (Li et al., 2020). However, except for intravenous remdesivir and dexamethasone, which have modest effects in moderate to severe COVID-19, no strong clinical evidence supports the efficacy of any other drugs against SARS-CoV-2 (Asselah et al., 2021). Vaccines are the most effective strategy to prevent infectious diseases, but the limited production means it may take years to vaccinate the majority of the population on earth. Therefore, before effective treatment and vaccines can contain the pandemic successfully, detection and isolation of people infected with SARS-CoV-2 are the first and most effective steps to curb the spread of SARS-CoV-2.

By far, RT-qPCR is the golden standard for COVID-19 test. The success of RT-qPCR depends first and foremost on the performance of the primers and further on the sequence matching between the primers and SARS-CoV-2 genome. The primers currently used are all based on the earliest version of SARS-CoV-2 genome sequence. However, SARS-CoV-2 constantly undergoes mutations. The mutations in S protein may change the antigenicity, leading to immune escape and affecting vaccine design (Naqvi et al., 2020; Silveira et al., 2021). There is no doubt that, once mutations occur in the RT-qPCR-targeted region, the RT-qPCR test may give a false-negative result¹⁹. Since thousands of single-nucleotide variants have been identified by far and the number of mutations is still growing rapidly, the possibility for the appearance of some variants that can evade the RT-qPCR test is also becoming increasingly higher. To avoid this situation, we suggest an innovative primer design strategy such that the primers should be designed based on both the nucleic acid sequence and the structure of the protein encoded (Figure 7). By using this strategy, we designed two pairs of primers targeting the N gene. Experiments showed that these primers have the same sensitivity and high specificity as the primers designed by US and Chinese CDCs. More importantly, any mutation affecting the performance of these primers (especially NAm1) would lead to the production of a protein losing structural stability and thus an inactive virus variant. N protein contains five tryptophans (W52, W108, W132, W301, and W330) and seven methionines (M1, M101, M210, M234, M317, M322, and M411). By analyzing the three-dimensional structure of N protein, we can easily figure out which residues are in the structure core and choose them for primer design. By this means, we need not try every Trp and Met to find the best place, thus greatly reducing the workload. This strategy would be much more convenient when the target protein is much bigger. Therefore, the primer pairs that we designed were the product of integrating the gene sequence and the three-dimensional structure of the protein rather than simply based on the gene sequence and degeneracy of codons.

FIGURE 7

Figure 7 Primer design principle based on protein structure and amino acid codon.

Since SARS-CoV-2 has many genes, it is not difficult to design primers to target other genes based on the same principle and form a multi-channel PCR system. This way, the risk of PCR failure caused by SARS-CoV-2 mutations is minimized. Furthermore, we believe that this strategy can also be used to design primers for the diagnosis of other viruses and even bacterial or fungal pathogens. Since all pathogens have numerous variants and only a small fraction of them are sequenced, the new strategy for primer design will no doubt ease the concerns that the RT-qPCR test may fail to detect certain variants.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

LG conceived and directed this research, is the guarantor of this work, has full access to all of the data in the study, and takes responsibility for the integrity of the data and the accuracy of the data analysis. HD designed the methods. HD, SW, and JZ performed most of the experiments and wrote the manuscript. KZ, FZ, and HW helped with the experiments. SX and WH helped with the project and the writing. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by the Shandong Provincial Key Research and Development Program (2020CXGC011305).

Conflict of Interest

JZ and SX were employed by Shandong Shtars Medical Technology Co., Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors thank Zhifeng Li and Jing Zhu at the Analysis and Testing Center of SKLMT (State Key Laboratory of Microbial Technology, Shandong University) for their assistance in the use of the qTOWER3 instrument for the RT-qPCR assays.

References

Abdullahi, I. N., Emeribe, A. U., Ajayi, O. A., Oderinde, B. S., Amadu, D. O., Osuji, A. I. (2020). Implications of SARS-CoV-2 Genetic Diversity and Mutations on Pathogenicity of the COVID-19 and Biomedical Interventions. J. Taibah Univ. Med. Sci. 15 (4), 258–264. doi: 10.1016/j.jtumed.2020.06.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Asselah, T., Durantel, D., Pasmant, E., Lau, G., Schinazi, R. F. (2021). COVID-19: Discovery, Diagnostics and Drug Development. J. Hepatol. 74 (1), 168–184. doi: 10.1016/j.jhep.2020.09.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Bielecki, M., Wójtowicz, H., Olczak, T. (2014). Differential Roles of Tryptophan Residues in Conformational Stability of Porphyromonas Gingivalis HmuY Hemophore. BMC Biochem. 15 (1), 2. doi: 10.1186/1471-2091-15-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Bru, D., Martin-Laurent, F., Philippot, L. (2008). Quantification of the Detrimental Effect of a Single Primer-Template Mismatch by Real-Time PCR Using the 16S rRNA Gene as an Example. App Environ. Microb. 74 (5), 1660–1663. doi: 10.1128/aem.02403-07

CrossRef Full Text | Google Scholar

Bulut, C., Kato, Y. (2020). Epidemiology of COVID-19. Turk J. Med. Sci. 50 (SI-1), 563–570. doi: 10.3906/sag-2004-172

PubMed Abstract | CrossRef Full Text | Google Scholar

Campos, M. J., Quesada, A. (2017). Strategies to Improve Efficiency and Specificity of Degenerate Primers in PCR. Methods Mol. Biol. 1620, 75–85. doi: 10.1007/978-1-4939-7060-5_4

PubMed Abstract | CrossRef Full Text | Google Scholar

Davies, N. G., Barnard, R. C., Jarvis, C. I., Kucharski, A. J., Munday, J., Pearson, C. A. B., et al. (2021). Estimated Transmissibility and Severity of Novel SARS-CoV-2 Variant of Concern 202012/01 in England. medRxiv 2020, 20248822. doi: 10.1101/2020.12.21.20248640

CrossRef Full Text | Google Scholar

Gao, Y., Yan, L., Huang, Y., Liu, F., Zhao, Y., Cao, L., et al. (2020). Structure of the RNA-Dependent RNA Polymerase From COVID-19 Virus. Science 368 (6492), 779–782. doi: 10.1126/science.abb7498

PubMed Abstract | CrossRef Full Text | Google Scholar

Guan, W. D., Chen, L. P., Ye, F., Ye, D., Wu, S. G., Zhou, H. X., et al. (2020). High-Throughput Sequencing for Confirmation of Suspected 2019-Ncov Infection Identified by Fluorescence Quantitative Polymerase Chain Reaction. Chin. Med. J. (Engl) 133 (11), 1385–1386. doi: 10.1097/CM9.0000000000000792

PubMed Abstract | CrossRef Full Text | Google Scholar

Jin, Y., Yang, H., Ji, W., Wu, W., Chen, S., Zhang, W., et al. (2020). Virology, Epidemiology, Pathogenesis, and Control of COVID-19. Viruses 12 (4), 372. doi: 10.3390/v12040372

CrossRef Full Text | Google Scholar

Korber, B., Fischer, W. M., Gnanakaran, S., Yoon, H., Theiler, J., Abfalterer, W., et al. (2020). Tracking Changes in SARS-CoV-2 Spike: Evidence That D614G Increases Infectivity of the COVID-19 Virus. Cell 182 (4), 812–827.e19. doi: 10.1016/j.cell.2020.06.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Lippi, G., Simundic, A. M., Plebani, M. (2020). Potential Preanalytical and Analytical Vulnerabilities in the Laboratory Diagnosis of Coronavirus Disease 2019 (COVID-19). Clin. Chem. Lab. Med. 58 (7), 1070–1076. doi: 10.1515/cclm-2020-0285

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, K., Shrivastava, S., Stockwell, T. B. (2015). Degenerate Primer Design for Highly Variable Genomes. Methods Mol. Biol. 1275, 103–115. doi: 10.1007/978-1-4939-2365-6_7

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Gayle, A. A., Wilder-Smith, A., Rocklov, J. (2020). The Reproductive Number of COVID-19 is Higher Compared to SARS Coronavirus. J. Travel Med. 27 (2), taaa021. doi: 10.1093/jtm/taaa021

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, R., Han, H., Liu, F., Lv, Z., Wu, K., Liu, Y., et al. (2020). Positive Rate of RT-PCR Detection of SARS-CoV-2 Infection in 4880 Cases From One Hospital in Wuhan, China, From Jan to Feb 2020. Clin. Chim. Acta 505, 172–175. doi: 10.1016/j.cca.2020.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, D., Zhang, J., Li, J. (2020). Primer Design for Quantitative Real-Time PCR for the Emerging Coronavirus SARS-CoV-2. Theranostics 10 (16), 7150–7162. doi: 10.7150/thno.47649

PubMed Abstract | CrossRef Full Text | Google Scholar

Maher-Sturgess, S. L., Forrester, N. L., Wayper, P. J., Gould, E. A., Hall, R. A., Barnard, R. T., et al. (2008). Universal Primers That Amplify RNA From All Three Flavivirus Subgroups. Virol. J. 5, 16. doi: 10.1186/1743-422X-5-16

PubMed Abstract | CrossRef Full Text | Google Scholar

Naqvi, A. A. T., Fatima, K., Mohammad, T., Fatima, U., Singh, I. K., Singh, A., et al. (2020). Insights Into SARS-CoV-2 Genome, Structure, Evolution, Pathogenesis and Therapies: Structural Genomics Approach. Biochim. Biophys. Acta Mol. Basis Dis. 1866 (10), 165878. doi: 10.1016/j.bbadis.2020.165878

PubMed Abstract | CrossRef Full Text | Google Scholar

Ohmae, E., Sasaki, Y., Gekko, K. (2001). Effects of Five-Tryptophan Mutations on Structure, Stability and Function of Escherichia Coli Dihydrofolate Reductase. J. Biochem. 130 (3), 439–447. doi: 10.1093/oxfordjournals.jbchem

PubMed Abstract | CrossRef Full Text | Google Scholar

Oude Munnink, B. B., Sikkema, R. S., Nieuwenhuijse, D. F., Molenaar, R. J., Munger, E., Molenkamp, R., et al. (2021). Transmission of SARS-CoV-2 on Mink Farms Between Humans and Mink and Back to Humans. Science 371 (6525), 172–177. doi: 10.1126/science.abe5901

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Y., Du, N., Lei, Y., Dorje, S., Qi, J., Luo, T., et al. (2020). Structures of the SARS-CoV-2 Nucleocapsid and Their Perspectives for Drug Design. EMBO J. 39 (20), e105938. doi: 10.15252/embj.2020105938

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahman, M. S., Hoque, M. N., Islam, M. R., Islam, I., Mishu, I. D., Rahaman, M. M., et al. (2021). Mutational Insights Into the Envelope Protein of SARS-CoV-2. Gene Rep. 22, 100997. doi: 10.1016/j.genrep.2020.100997

PubMed Abstract | CrossRef Full Text | Google Scholar

Ramirez, J. D., Munoz, M., Hernandez, C., Florez, C., Gomez, S., Rico, A., et al. (2020). Genetic Diversity Among SARS-CoV2 Strains in South America may Impact Performance of Molecular Detection. Pathogens 9 (7), 580. doi: 10.3390/pathogens9070580

CrossRef Full Text | Google Scholar

Santiveri, C. M., Jiménez, M. A. (2010). Tryptophan Residues: Scarce in Proteins But Strong Stabilizers of β-Hairpin Peptides. Biopolymer 94 (6), 779–790. doi: 10.1002/bip.21436

CrossRef Full Text | Google Scholar

Silveira, M. M., Moreira, G., Mendonca, M. (2021). DNA Vaccines Against COVID-19: Perspectives and Challenges. Life Sci. 267, 118919. doi: 10.1016/j.lfs.2020.118919

PubMed Abstract | CrossRef Full Text | Google Scholar

Souvenir, R., Buhler, J., Stormo, G., Zhang, W. (2007). An Iterative Method for Selecting Degenerate Multiplex PCR Primers. Methods Mol. Biol. 402, 245–268. doi: 10.1007/978-1-59745-528-2_12

PubMed Abstract | CrossRef Full Text | Google Scholar

Stadhouders, R., Pas, S. D., Anber, J., Voermans, J., Mes, T. H. M., Schutten, M. (2010). The Effect of Primer-Template Mismatches on the Detection and Quantification of Nucleic Acids Using the 5′ Nuclease Assay. J. Mol. Diagn. 12 (1), 109–117. doi: 10.2353/jmoldx.2010.090035

PubMed Abstract | CrossRef Full Text | Google Scholar

Tegally, H., Wilkinson, E., Giovanetti, M., Iranzadeh, A., Fonseca, V., Giandhari, J., et al. (2020). Emergence and Rapid Spread of a New Severe Acute Respiratory Syndrome-Related Coronavirus 2 (SARS-CoV-2) Lineage With Multiple Spike Mutations in South Africa. medRxiv 2020, 20248640. doi: 10.1101/2020.12.21.20248640

CrossRef Full Text | Google Scholar

Udugama, B., Kadhiresan, P., Kozlowski, H. N., Malekjahani, A., Osborne, M., Li, V. Y. C., et al. (2020). Diagnosing COVID-19: The Disease and Tools for Detection. ACS Nano 14 (4), 3822–3835. doi: 10.1021/acsnano.0c02624

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Horby, P. W., Hayden, F. G., Gao, G. F. (2020a). A Novel Coronavirus Outbreak of Global Health Concern. Lancet 395 (10223), 470–473. doi: 10.1016/S0140-6736(20)30185-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Hozumi, Y., Yin, C., Wei, G. W. (2020d). Mutations on COVID-19 Diagnostic Targets. Genomics 112 (6), 5204–5213. doi: 10.1016/j.ygeno.2020.09.028

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C., Liu, Z., Chen, Z., Huang, X., Xu, M., He, T., et al. (2020b). The Establishment of Reference Sequence for SARS-CoV-2 and Variation Analysis. J. Med. Virol. 92 (6), 667–674. doi: 10.1002/jmv.25762

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, R., Hozumi, Y., Yin, C., Wei, G. W. (2020c). Decoding SARS-CoV-2 Transmission and Evolution and Ramifications for COVID-19 Diagnosis, Vaccine, and Medicine. J. Chem. Inf Model 60 (12), 5853–5865. doi: 10.1021/acs.jcim.0c00501

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, A., Niu, P., Wang, L., Zhou, H., Zhao, X., Wang, W., et al. (2020). Mutations, Recombination and Insertion in the Evolution of 2019-Ncov. bioRxiv 2020, 971101. doi: 10.1101/2020.02.29.971101

CrossRef Full Text | Google Scholar

Xia, Y., Xun, L. (2017). Revised Mechanism and Improved Efficiency of the QuikChange Site-Directed Mutagenesis Method. Methods Mol. Biol. 1498, 367–374. doi: 10.1007/978-1-4939-6472-7_25

PubMed Abstract | CrossRef Full Text | Google Scholar

Younes, N., Al-Sadeq, D. W., Al-Jighefee, H., Younes, S., Al-Jamal, O., Daas, H. I., et al. (2020). Challenges in Laboratory Diagnosis of the Novel Coronavirus SARS-CoV-2. Viruses 12 (6), 582. doi: 10.3390/v12060582

CrossRef Full Text | Google Scholar

Zhu, N., Zhang, D., Wang, W., Li, X., Yang, B., Song, J., et al. (2020). A Novel Coronavirus From Patients With Pneumonia in China 2019. N. Engl. J. Med. 382 (8), 727–733. doi: 10.1056/NEJMoa2001017

PubMed Abstract | CrossRef Full Text | Google Scholar