1. National Center for Clinical Laboratories, Beijing Hospital, National Center of Gerontology; Institute of Geriatric Medicine, Chinese Academy of Medical Sciences, Beijing, People's Republic of China.
2. Graduate School, Peking Union Medical College, Chinese Academy of Medical Sciences, Beijing 100730, People's Republic of China.
3. Beijing Engineering Research Center of Laboratory Medicine, Beijing, People's Republic of China.
#These authors contributed equally to this manuscript.
In December 2019, a new coronavirus disease (COVID-19) outbreak occurred in Wuhan, China. Severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2), which is the seventh coronavirus known to infect humans, is highly contagious and has rapidly expanded worldwide since its discovery. Quantitative nucleic acid testing has become the gold standard for diagnosis and guiding clinical decisions regarding the use of antiviral therapy. However, the RT-qPCR assays targeting SARS-CoV-2 have a number of challenges, especially in terms of primer design. Primers are the pivotal components of a RT-qPCR assay. Once virus mutation and recombination occur, it is difficult to effectively diagnose viral infection by existing RT-qPCR primers. Some primers and probes have also been made available on the WHO website for reference. However, no previous review has systematically compared the previously reported primers and probes and described how to design new primers in the event of a new coronavirus infection. This review focuses on how primers and probes can be designed methodically and rationally, and how the sensitivity and specificity of the detection process can be improved. This brief review will be useful for the accurate diagnosis and timely treatment of the new coronavirus pneumonia.
Keywords: coronavirus, SARS-CoV-2, quantitative nucleic acid testing, primer design, sensitivity
In December 2019, outbreak of a new coronavirus disease (COVID-19), defined by the International Committee on Taxonomy of Viruses, occurred in Wuhan, China. This virus has since been named severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) . SARS-CoV-2 is highly contagious and has rapidly expanded worldwide since its discovery. As of 12 May, 2020, a total of 4,088,848 cases of SARS-CoV-2 infection have been confirmed worldwide with 283,153 deaths , although this may be underestimated due to inadequate testing in many countries .
SARS-CoV-2 is an enveloped, positive-sense, single-stranded RNA virus. It possesses a comparatively large genome approaching 30 kb. The genome is arranged in the order of a 5′ untranslated region (UTR)-replicase complex (open reading frame [ORF] 1ab)-structural protein [Spike(S)-Envelope(E)-Membrane (M)-Nucleocapsid (N)]-3′UTR and non-structural ORFs . The S-protein has a strong binding affinity to human ACE2 , and this, along with a long incubation time before manifesting symptoms, means that SARS-CoV-2 has high transmissibility with a mortality rate of approximately 3% [6,7]. This also poses considerable challenges regarding the timely treatment and effective prevention and control of COVID-19. To better understand the origin of SARS-CoV-2 and its genetic relationship with other coronaviruses, phylogenetic analyses of coronavirus sequences from various sources can be performed using the neighbour-joining method in MEGA X . ORF 1ab, E, and N genes are highly conserved among sarbecoviruses, which are a subgenus of betacoronavirus. The RNA-dependent RNA polymerase (RdRp; nsp12) positioned on ORF 1ab plays a crucial role in RNA synthesis [9,10] with features of rapid mutation and recombination . As a result, the conserved regions (ORF 1ab, E, and N genes) are usually selected as the standard target genes for primer and probe design .
SARS-CoV-2 has been identified as the seventh coronavirus known to infect humans . The previous six coronaviruses were: Human coronavirus (HCoV)-229E, HCoV-NL63, HCoV-HKU1, HCoV-OC43, Middle East respiratory syndrome coronavirus (MERS-CoV), and severe acute respiratory syndrome coronavirus (SARS-CoV). HCoV-229E and HCoV-NL63 are human alphacoronaviruses, whereas HCoV-HKU1, HCoV-OC43, MERS-CoV, and SARS-CoV are human betacoronaviruses. SARS-CoV-2 is a novel betacoronavirus, and is genetically similar to the SARS-like bat coronaviruses in the subgenus Sarbecovirus . A complete genomic comparison showed that the new coronavirus has the highest homology with Bat SL-CoVZC45 and Bat SL-CoVZXC21 coronaviruses, but differs from the currently known SARS with a similarity of just 79.7% . Although substantial genetic differences exist between the novel coronavirus and other betacoronaviruses, cross-reactions in RT-qPCR or antibody assays for SARS or other betacoronaviruses are possible if the primers and antigenic epitopes are not carefully selected . Therefore, it is necessary to design specific primers and probes in the regions with the lowest feasible similarity to other viruses to avoid false detection of SARS-CoV .
The timely and accurate laboratory testing of samples from cases under investigation are essential parts of the management of emerging infections. Molecular techniques have been used successfully to identify infectious diseases for many years. As can be observed from the identification path of SARS-CoV-2 [18,19], high-throughput sequencing is a powerful tool for the discovery of pathogens, changing the way we respond to infectious disease outbreaks, improving our understanding of disease occurrence, accelerating the identification of pathogens, and promoting data sharing [12,13]. However, the analysis method is labour-intensive and time-consuming and, thus, cannot be widely used as a detection method or screening tool.
As sequence information about SARS-CoV-2 has recently become available, quantitative nucleic acid testing has become the gold standard for clinical decisions regarding the use of antiviral therapy [20-22]. The most recently developed nucleic acid testing method is quantitative real-time PCR (RT-qPCR) [23,24], although RT-qPCR assay design optimization can be a complicated process. In addition to being rapid, the RT-qPCR assay has the advantage of a standardized protocol that can be easily adapted for the detection of other respiratory viruses. RT-qPCR can be performed under uniform amplification conditions; thus, target-specific primer and probe sets can be utilized . Moreover, testing can also be performed using a high-throughput system .
However, there are a number of problems related to the RT-qPCR assays targeting SARS-CoV-2. First, the sensitivity of this method is not high in clinical application, and studies have reported cases of positive CT scan results and negative RT-qPCR results at initial presentation . Some researchers and clinicians have argued that CT imaging should be used to identify SARS-CoV-2 infection, evaluate treatment response and prognosis of COVID-19 [28, 29], but a number of cases have shown progressive multiple peripheral ground-glass opacities in both lungs despite negative RT-qPCR results . Furthermore, an accurate detection of SARS-CoV-2 virus can be made from oropharyngeal, nasopharyngeal swabs and/or lower respiratory tract samples (sputum or tracheal aspirates) ; thus, sampling could also be one of the reasons for low sensitivity. The primers, probes, and reagents used for detection could be another critical factor. Finally, RT-qPCR primers are designed to amplify the target regions of the SARS-CoV-2 genome. However, the novel coronaviruses are RNA viruses, and once mutations and recombination occur, RT-qPCR primers will not be able to effectively amplify the viral sequences. Thus, mutations and recombination can reduce the sensitivity of the RT-qPCR.
Primers are pivotal components of a qPCR assay. Some primers and probes have also been made available on the WHO website for reference. However, no previous review has systematically described how to design primers and probes methodically and rationally in the face of a new coronavirus. Therefore, the focus of the present review is to discuss effective primer designing and improving the sensitivity and specificity of the detection process. This review also discusses the need for accurate diagnosis and timely treatment of the new coronavirus pneumonia.
At present, gene sequences of the six known coronaviruses (HCoV-229E, HCoV-NL63, HCoV-HKU1, HCoV-OC43, MERS-CoV, and SARS-CoV) have been published, and primers and probes can be designed based on the available genomic information. When an unknown coronavirus suddenly emerges and causes a new type of infectious pneumonia, how should we respond? In such a scenario, designing effective primers with high specificity and sensitivity can have immense value for diagnosis, tracking, and tracing. The design of SARS-CoV-2 primers and probes and the step-by-step analysis process to accurately detect the virus are shown in Figure 1.
Zhang et al. were the first to determine the full-length genomic sequence of the SARS-CoV-2 virus . On 26 Dec, 2019, six days after the onset of disease, a 41-year-old man with fever, unproductive cough, pain, and weakness was admitted to the Central Hospital of Wuhan. Bronchoalveolar lavage fluid (BALF) was collected to perform deep meta-transcriptomic sequencing. Primers used for RT-qPCR were designed based on the whole genome of WH-Human 1 coronavirus (MN908947). The detailed reference sequences of the viral RNA (GenBank MN908947) are shown in Table 1. The genome sequence of the causative novel coronavirus was shared through the Global Initiative on Sharing All Influenza Data (GISAID) platform (http://www.gisaid.org/) from 12 Jan 2020 . Subsequently, a number of other SARS-CoV-2 sequences were released; these sequences were found to be almost identical .
Detailed information of the reference sequence MN908947
|ORF||Function||Location (nt)||Length (bp)|
|1a||Encoded nonstructural proteins (nsp1 to nsp11), essential for viral replication, viral assembly, immune response modulation, etc.||--||--|
|1b||Encoded nonstructural proteins (nsp12 to nsp16), essential for viral replication||--||--|
|S||Binding to cell receptor and mediating virus-cell fusion||21,563-25,384||3,822|
|E||Envelope protein, virus assembly, and morphogenesis||26,245-26,472||228|
|M||Membrane protein, virus assembly||26,523-27,191||669|
|N||Nucleocapsid protein, forms complexes with genomic RNA, interaction with M protein for viral assembly||28,274-29,533||1,260|
The analysis process of primer and probe design for SARS-CoV-2 in chronological order. BALF: bronchoalveolar lavage fluid; GISAID: Global Initiative on Sharing All Influenza Data.
On 1 Jan 2020, Corman et al. used the close correlation between the SARS-CoV-2 and SARS coronavirus genomes to artificially synthesize the 2019-nCoV sequence and establish an experimental detection process in the absence of virus isolates . On 10 Jan 2020, Chan et al reported that two patients with fever, respiratory symptoms, and pulmonary infiltrates on chest radiographs were enrolled at the University of Hong Kong-Shenzhen Hospital. Subsequently, between 11 to 15 Jan 2020, five other members of the patients' family also presented to the hospital for clinical assessment . The first set of primers targeted 344 bp of the RdRp gene of all SARS-related coronaviruses. Next, these positive clinical samples were analysed by Nanopore sequencing, and a second set of primers were designed from the SARS-CoV-2 genome sequence and targeted a 158 bp region of the Spike (S) gene of this novel coronavirus . Moreover, Won et al. designed primers based on the sequence of the first patient infected with SARS-CoV-2 in Korea [34,35]. On 22 Jan 2020, the downloaded sequence and those for SARS coronavirus, bat SARS-like coronaviruses, and other representative coronaviruses among the first publicly available sequences in GenBank (Accession number: MN908947) were edited and aligned. Two sequence regions (ORF1b and N) that are highly conserved among sarbecoviruses were selected for primer and probe design in the study by Chu et al. . On 14 Feb 2020, 95 full-length genomic sequences of SARS-CoV-2 strains from the NCBI and GISAID databases were retrieved, and a reference sequence was established. Researchers found that the overall variation in the ORF regions was low; 13 variation sites in the 1a, 1b, S, 3a, M, 8, and N regions were identified, among which positions nt28144 in ORF 8 and nt8782 in ORF 1a showed mutation rates of 30.53% (29/95) and 29.47% (28/95), respectively . There may be selective mutations in SARS-CoV-2, and it is necessary to avoid certain regions when designing primers and probes.
Generally, in the face of unknown pathogenic and highly infectious coronaviruses, the first step involves extracting RNA from patient samples and sequencing the whole genome. However, the whole genome sequencing operation is difficult and costly, and it is impossible to use this method for all patients with suspected infection. Genome sequencing can be used as a detection method in a few samples, and the results can be used as a reference sequence to design primers and probes targeting the target gene, so that RT-qPCR detection can be achieved. However, RNA viruses are susceptible to mutations. Shen et al. recently identified a remarkable level of viral diversity in some infected patients, accounting for a median number of four intra-individual viral variants, which is suggestive of the rapid evolution of SARS-CoV-2 . With the increasing number of cases, results from genome sequencing will also increase. Genome sequences of multiple coronaviruses can be compared to determine the reference sequence. For primer and probe design, it is more reliable and scientifically sound to select conserved regions instead of mutation sites.
Assays are useful only if the correct target has been identified and used in the assay design. Therefore, it is necessary to know as much information as possible related to the RNA genome sequences. The first step involves obtaining the maximum amount of information from sequence databases; it is crucial to refer to the accession or individual transcript number of related sequences for assay design.
To design a PCR primer set, a reference sequence is needed to identify the exact sequence being targeted and from where to select the primer pair candidates. Ideally, the reference sequence consists of not only the regions targeted for amplification, but also all the nucleotide sequences that will be involved in the amplification process. This is essential to ensure that the designed primers do not amplify non-specific products .
When a new coronavirus is identified, the genome sequence of the virus can first be obtained through standard high-throughput sequencing of patient samples, which can be used as a reference sequence. When the coronavirus spreads further, the sequencing results of different coronavirus strains can be uploaded to the GISAID (https://www.gisaid.org) and NCBI (https://www.ncbi.nlm.nih.gov/) databases. Primer 7.0 can be used to compare the reference nucleotide sequence with those of related human isolates and to analyse the variation at different locations. In addition, the ClustalW program of MEGA (7.0.14) and Mafft (7.450) software can be used to conduct homology analyses and multiple sequence alignments. These sequences can be retrieved from experimental data or online databases. In conclusion, the reference sequence is determined by selecting the highest frequency nucleotide at each position.
Usually, degenerate primers are designed in two different instances. The first is when the target sequence is unknown or presents variability. In such a case, the target sequence can be inferred through the amino acid sequence. Degenerate primers can overcome the uncertainty of the target sequence due to the redundancy of the genetic code . The second is when degenerate primers are used to amplify targeted sequences across multiple genotypes by allowing variant primers to simultaneously amplify multiple sequences. Previous studies have detected five different SARS-CoV-2 haplotypes, and the results indicated active genetic recombination . The mutation regions may reduce the accuracy of RT-qPCR detection . Thus, designing degenerate primers and/or probes is an effective method for addressing the challenge of virus mutations.
Primers can be designed as regular or degenerate. Degenerate primer design requires a set of reference sequences. When regular primers are synthesised, only two populations of primers are present in the PCR reaction: the forward and reverse primer populations. When degenerate primers are synthesized, there may be multiple forward and reverse primers. A good principle for designing highly efficient degenerate primers is to establish a limit of degeneracy of up to 64 (the number of different sequences within the mixture primer). It is recommended that the degenerate primers contain R(A/G), Y(C/T), S(C/G), W(A/T), K(G/T), and M(A/C), which have two variants at each position, instead of B(C/G/T), D(A/G/T), H(A/C/T), V(A/C/G), and N(G/A/T/C), which have more than two variants at each position .
Many institutions have studied primer-probe sets for the confirmation of SARS-CoV-2. The molecular diagnosis of COVID-19 is currently carried out by one-step RT-qPCR targeting SARS-CoV-2. Primers and probes have been developed by researchers at China CDC (China), Charité-Universitätsmedizin Berlin Institute of Virology (Germany), The University of Hong Kong (Hong Kong), the National Institute of Health (Thailand), and the Centers for Disease Control and Prevention (CDC, USA) (https://www.who.int/emergencies/diseases/novel-coronavirus-2019/technical-guidance/laboratory-guidance). The primers and probes reported for the SARS-CoV-2 RT-qPCR assays are shown in Table 2.
A conserved region in the structural protein envelope E gene is usually selected for the E gene assay, which normally functions as the first line screening assay. The E genes of the reference (MN908947) and SARS (NC_004718) sequences were aligned using nucleotide BLAST (https://blast.ncbi.nlm.nih.gov/). The results are shown in Figure 2. The E gene primer/probe set described by Corman et al. was found to be more sensitive than those of the N gene and RdRp . The position of the primer targeting the E gene is shown in Figure 2 and indicated with pink text. From the alignment of the E gene between the SARS-CoV-2 reference sequence (MN908947) and SARS-CoV (NC_004718), there was 94% homology and the mismatch positions are marked in red in Figure 2. Corman et al. selected the forward and reverse primers in the complete match position so that both SARS-CoV-2 and SARS-CoV could be detected. It is believed that SARS has been eliminated in humans and the last reported human SARS case was detected in 2004 [48,49]. Samples with positive SARS-CoV results were considered to be infected by SARS-CoV-2 or its related animal coronaviruses. Selecting the common sequence of bat-associated SARS-related viruses is essential to ensure broad sensitivity even in the case of multiple independent acquisitions of variant viruses. Therefore, this pancoronavirus assay is a useful tool for screening sample collections for the presence of all known and potentially unknown coronaviruses. When designing the specific primer/probe targeting the E gene, the mismatch position should be considered (Figure 2).
Compared with the E gene length in SARS-CoV-2, the gene lengths of ORF 1ab and N are longer, and more variation points are present. Therefore, regions in the ORF 1ab and N genes are usually selected for the target-specific forward and reverse primers designed for confirmatory testing. Research has shown that SARS-CoV-2 is distinct from SARS-CoV in the phylogeny of the complete RdRp gene . As shown in Table 2, tests for the RdRp gene in Germany and ORF1b-nsp14 gene in Hong Kong used degenerate probes and primers. There are two main reasons for this detection approach. First, the genetic diversity of 2019-nCoV in humans and animals is yet to be fully determined. Second, many laboratories lack positive controls for SARS-CoV-2. The serially diluted RNA samples extracted from SARS-CoV-infected cells usually work as the positive control (obtained via the European virus archive global (EVAg), https://www.european-virus-archive.com) . SARS-CoV-2 isolation is not recommended as a routine diagnostic procedure. In vitro transcribed RNA of the target gene of SARS-CoV-2 can be cloned into a plasmid as the RNA standard. On one hand, the degenerate probes and primers can detect the positive control (SARS-CoV), and on the other hand, they can detect SARS-CoV-2. For example, Corman et al. formulated a broad-range probe reacting with SARS-CoV and SARS-CoV-2, an additional specific probe that reacts only with SARS-CoV-2. In their study, RdRp and E genes are more sensitive than N genes, and the RdRp gene assay is recommended as a confirmatory assay. However, Chu et al. found that the N gene assay was approximately 10 times more sensitive than the ORF-1b gene assay in detecting positive clinical specimens from two patients via RT-qPCR ; based on their detection performance, the N gene RT-qPCR is recommended as a screening assay, and the ORF1b-nsp14 assay is recommended as a confirmatory assay. Further experimental studies are required to evaluate the sensitivity of the primer/probe set shown in Table 2. For one target gene, there may exist many regions for which the primer/probe is designed. We believe that it is more reasonable to compare the primer/probe sensitivity for the same target gene among different test regions. However, if the regions of multiple target genes are different in each study, the comparability among different studies is lost. Alternatively, the primer/probe sets from different studies can be evaluated under the same experimental conditions with the same test samples.
The USA CDC designed N1, N2, and N3 genes (shown in Table 2) as the target of the SARS-CoV-2 Real-Time RT-PCR Diagnostic Panel. In the SARS-CoV-2 RT-PCR assay, N1 and N2 were designed for the specific detection of SARS-CoV-2, whereas N3 was designed for the universal detection of SARS-like coronaviruses . Experiments found that the CDC N1 and N2 primer/probe sets performed better than the N3 set, and the N3 primer and probe set was removed from the Diagnostic Panel of revision 2 for the CDC 2019-nCoV Real-Time RT-PCR Diagnostic Panel.
Researchers in Korea performed experiments to compare the previously reported primer-probe sets from the WHO for molecular diagnosis. Their results revealed that in the case of targeting the RdRp/ORF1 region, the ORF 1ab (China) set might be more sensitive than others. 2019-nCoV_N2 (USA) and NIID_2019-nCOV_N (Japan) sets may be recommended for the sensitive RT-qPCR assay of the N region. Therefore, an appropriate combination of the ORF 1ab (China), 2019-nCoV_N2 (USA), and NIID_2019-nCOV_N (Japan) sets should be selected for sensitive and reliable laboratory confirmation of SARS-CoV-2 . Other researchers found that the E gene primer/probe set described by Corman et al. and the N2 set developed by the CDC were the most sensitive assays after evaluating seven different primer/probe sets (RdRp, E, and N genes in the study by Corman et al. and N1, N2, and N3 by CDC and RdRp developed by the University of Washington Clinical Virology Lab) .
From the perspective of primers, there are three main ways to improve the sensitivity of RT-qPCR: primer concentrations, degenerate primers, and multi-target detection. First, primer concentrations for probe-based assays were 300-900 nM . The primer concentrations of all tests (Table 2) were within this range (data not shown), and increasing the concentration of primers appropriately may improve the sensitivity of RT-qPCR. For example, Vijgen et al. found that primer concentrations increased up to 400 nm per reaction can improve the sensitivity of the assay . Second, degenerate primers are used to cope with the genetic diversity of SARS-CoV-2. Finally, regarding multi-target detection, screening assays with a single target region are more vulnerable to sequence variations than dual- or triple-target assays [55,56]. Overall, the E gene assay was used for pan-sarbecovirus detection, and the RdRp assay and the N gene assay may eventually be utilised as confirmatory assays.
In addition to sensitivity, specificity is an important indicator when considering primers and probes. Whether there is cross-reactivity is an important aspect of evaluating specificity. The proportion of co-infections has been underestimated due to either a lack of sensitivity or the limitations of virus detection technologies . As co-infections can occur, all patients that meet the suspected case definition should be tested for COVID-19 regardless of whether another respiratory pathogen is found. However, all tests shown in Table 2 indicate no cross-reactivity with HCoV-229E, HCoV-NL63, HCoV-HKU1, HCoV-OC43, MERS-CoV, and other respiratory viruses, indicating that the specificity of the primers and probes listed in Table 2 is high.
A successful PCR reaction depends on a number of factors, most of which are related to the primer design quality. The primer length should be between 16 and 28 nucleotides, as shown in Table 2. Primer melting temperatures (Tm) should be between 50°C and 62°C, depending on the GC content. The GC content of the primer is preferably 45-55%. Moreover, the GC content and Tm value of the upstream and downstream primers should be kept close. The Tm difference between the two primers should be lower than 5°C , and the annealing temperature should be 5°C lower than that for the primer with the lowest Tm. Because the target sequence of the primers is variable for degenerate primers, a good approach is to set annealing temperatures between 50°C and 55°C, although specific primers could be used to estimate the PCR conditions more accurately. For the polymorphisms in the SARS-CoV-2 genomes, there may be mismatches in primers or probes that cause false-negative conditions . The significance of the regular examination of primers and probes has been emphasized . False negative results can be reduced by designing primers and probes with melting temperatures >60°C and run conditions of the assay with annealing temperatures at 55°C to tolerate one to two mismatches when amplifying a coronavirus gene with higher variability.
Information regarding the primers and probes reported for SARS-CoV-2 real-time reverse-transcription PCR assays
|Target genes||Country||Name||Sequence (5' → 3')||Reference sequence||Nucleotide position||Reference|
|RdRp gene||Germany||Forward primer;|
(W=A/T; R=A/G; M=A/C; S=G/C)
|Alignment of sequenced virus genomes||--|||
|ORF 1ab||China||Forward primer:|
|ORF 1ab||China||Forward primer:|
|ORF1b-nsp14||Hong Kong||Forward primer:|
(W=A/T; Y=C/T; R=A/G)
|N gene||China||Forward primer:|
|N gene||Hong Kong||Forward primer,|
|N1 gene||USA||Forward primer:|
|GAC CCC AAA ATCAGCGAA AT|
|N2 gene||USA||Forward primer:|
|TTACAA ACATTGGCCGCA AA|
|N3 gene||USA||Forward primer:|
|GGGAGCCTTGAA TAC ACC AAA A|
|N gene||Thailand||Forward primer:|
|E gene||Germany||Forward primer:|
|RNase P gene (RP)||USA||Forward primer:|
Abbreviations: E, envelope protein gene; M, membrane protein gene; N, nucleocapsid protein gene; ORF, open reading frame; RdRp, RNA-dependent RNA polymerase gene; FAM, 6-carboxyfluorescein; BBQ, blackberry quencher; BHQ-1, black hole quencher 1; TAMRA, carboxytetramethylrhodamine; IBFQ, Iowa Black FQ.
Sequence alignment of the reference (MN908947) and SARS (NC_004718) sequences. Red colour marks the mismatch positions and pink colour marks the forward primer and reverse complementary sequence of the reverse primer in the study by Corman et al. .
The 5′- and 3′-ends of the primers play different roles; although the 5′-end contributes marginally, matching the 3′ end is critical for PCR efficiency amplification. It has been found that a single mismatch in the last three nucleotides counting from the 3′-end is acceptable, but the second impaired position drops efficiency below the detectable level on agarose gel electrophoresis .The base at the terminal 3' bases of the primer generally does not require A nucleotide. More attention should be paid to the last five bases; in particular, no more than three G/C residues should be included in the last five bases and the primer should never end with three consecutive G/C residues, such as GGG or CCC, otherwise the complementary, dimer, or hairpin structure at the 3′ end of the primer may cause the PCR reaction to fail. Multiple pairs of primers in the same amplification system can also increase the probability of primer dimer formation, which can be predicted by related software such as Oligo7 . The bases of primers should be randomly distributed, and it is suggested that there should be no more than four complementary bases between the forward and reverse primers. As reported in previous studies, primer dimers can be prevented by modifying with 2′-O-methyl bases in the penultimate base . The closer the probe is to the upstream primer, the higher the hydrolysis efficiency. However, it should be kept above 3 bp. The 5' end of the probe cannot start with a G base because it can quench the fluorophore . As shown in Table 2, most of the probes do not start with a G base except for that for the N gene developed in Hong Kong.
After designing the primer/probe following the general principle strategy, in silico analysis is important for specificity or exclusivity testing. For example, an alignment was performed with the oligonucleotide primer and probe sequences of the CDC 2019 nCoV RT-qPCR diagnostic panel with all publicly available nucleic acid sequences for 2019-nCoV in GenBank. NCBI Primer-BLAST can be used to confirm primer specificity. It is best to confirm specificity again after primer probe design is completed, especially since the primers and probes of multiple RT-qPCR need to be mixed together. It has also been suggested that primer/probe sets are selected from conserved sequences on stem structures to avoid loop structures , which may be different from the general principle that recommends designing primers to avoid secondary structure in templates. As a general principle, this is correct because the template can be paired to form the stem by itself, which may not be conducive to accurate PCR quantification. However, what we want to emphasize in this review is that, firstly, SARS-CoV-2 is an RNA virus, which is unstable and easily degraded by RNase A . The stem structure is more stable than the loops that are easily cleaved [65,66] and by avoiding loops and targeting stems, the original RNA value can be nearly reached, even if the SARS-CoV-2 RNA is degraded under improper transport and storage conditions before analysis, resulting in broader applicability . In addition, during reverse transcription in the RT-qPCR assay, the temperature is mostly set to 55°C, and this temperature is conducive to the effective opening of RNA secondary structures . During the amplification process, through denaturation, the secondary structures will also be opened. It is important to check that all secondary structures have a Tm (°C) less than the qPCR annealing temperature (normally 55-60°C) to ensure effective amplification . Therefore, we suggest that in the primer design for SARS-CoV-2, it is more effective to design the primer targeting position on the stem, at least to make sure that the 3 ' terminal bases are on the stem (as shown in Figure 3). However, further experiments are needed to explore and verify this.
The secondary structures of the ORF 1ab and N gene fragments amplified using the forward and reverse primers from Table 2 are shown in Figure 3. The secondary structures of the target genes were predicted by Mfold (http://mfold.rna.albany.edu/?qZmfold) using the minimum free energy method. The forward and reverse primers designed by each country avoided the loop positions as much as possible. Researchers also revealed that if the loop size, determined by the distance between paired palindrome elements, is less than 11 base pairs and the stem length, determined by the alignment of the palindrome elements, exceeds 14 hydrogen bonds, then there may be stem loop interference in the primer binding site . According to this, there are no stem loop interferences observed in the binding sites of primers developed in each country (Figure 3).
Secondary structures of the ORF 1ab and N gene fragments amplified using the forward and reverse primers. (A) The forward primer and reverse complementary sequence of the reverse primer for ORF 1ab from China are indicated in red. (B) The forward primer and reverse complementary sequence of the reverse primer for ORF1b from Hong Kong are shown in light blue. (C) The forward primer and reverse complementary sequence of the reverse primer for RdRp from Germany are shown in pink. (D) The forward primer and reverse complementary sequence of the reverse primer for the N1 gene from USA are shown in yellow and the N gene in Thailand are shown in dark blue. (E) The forward primer and reverse complementary sequence of the reverse primer for the N2 gene from USA are shown in yellow, the N gene from Hong Kong are shown in light blue, and the N gene from Japan are shown in green. (G) The forward primer and reverse complementary sequence of the reverse primer for the N3 gene from USA are shown in yellow and the N gene from China in red. Orange represents the overlap. F represents the forward primer and R represents the reverse complementary sequence of the reverse primer.
List of practical primer design tools
|Aim||Name||Software or Website|
|Secondary structure prediction||Mfold||http://mfold.rna.albany.edu/?qZmfold|
|Primer design||Primer Premier||Software|
|JCVI Primer Designer||https://sourceforge.net/projects/primerdesigner/|
|Test the primer specificity||NCBI Primer-BLAST||https://www.ncbi.nlm.nih.gov/tools/primer-blast/|
The primer and probe sequences can also be designed using software tools [70,71]. The commonly used primer design tools are listed in Table 3. Even if the same tools are used, the same suggested primers and probes may not always be the same, but this does not indicate that the quality of primer design is poor. For example, the primers and probes designed for the ORF 1ab of SARS-CoV-2 have different sequences distributed by the Chinese CDC and the USA CDC. Indeed, there are many primer and probe combinations that can meet the requirements of RT-qPCR experiments, so the highest quality combination does not always need to be used. Regardless of how advanced the algorithm of the design tool is, the performance of the probe primer ultimately depends on the experimental data and whether there are non-specific products to be confirmed by gel electrophoresis after PCR or a dissociation curve (which is recommended). The amplification conditions also need to be continuously optimized through experiments. For example, the optimal Tm and optimized oligo concentrations were eventually confirmed by temperature gradient experiments. The parameters given by the design tool are for reference only and should never replace the subsequent confirmation process. From an economic point of view, it is generally recommended to design a probe and three corresponding pairs of primers for a target and determine the best primer combination after experimental screening. The entire flow of the primer design is shown in Figure 4.
It is clear that no diagnostic modality is perfect. In order to obtain RT-qPCR results with high accuracy, each step must be performed satisfactorily in the preanalytical, analytical, and post-analytical steps. It is important to avoid sample contamination and ensure adequate procedures for specimen collection, handling, transport, and storage. For the analytical procedures, many problems, such as active viral recombination, a lack of harmonization of primers and probes, and non-specific PCR annealing, should be addressed. Finally, the healthcare personnel must be educated on the interpretation of RT-qPCR results. Despite high sensitivity, negative RT-qPCR results observed at one or two time points are insufficient to exclude SARS-CoV-2 infection in patients, especially for those with typical clinical presentations or clear epidemic indications. Therefore, negative results do not preclude SARS-CoV-2 infection and should not be used as the sole basis for treatment or other patient management decisions. Negative results must be combined with clinical observations, patient history, and epidemiological information.
Each RT-qPCR run should include both external and internal controls, and laboratories are encouraged to participate in external quality assessment schemes that should be established as soon as possible to monitor analytical quality and harmonize the assays. In many studies [34,44], the internal control primer set targeted the human ribonuclease P (RNase P) protein. This protein is an ubiquitous enzyme in all cells and cellular compartments . The GAPDH gene can also be used as the target gene for the internal control [17,57]. In order to identify regions of the viral genome, further analysis of the molecular targets is needed to achieve the highest possible diagnostic accuracy. The primer/probe set should be evaluated with further experiments, so that more sensitive and specific primers/probes can be developed. Ordering primers and probes to perform validation testing on functionality and potential contaminants is also recommended .
The entire flow of primer design. (A) Establishment of the reference sequence. Different types of samples are collected, RNA is extracted, and then high-throughput sequencing is performed to determine the reference sequence. Alternatively, after multiple sequencing results are aligned, the reference sequence can be determined. (B) Identification of the target genes. By comparing the SARS-CoV-2 genome with other bat-associated SARS-related viral genomes, conserved and non-conserved regions were identified. Conserved regions such as ORF 1ab, E gene and N gene can be used as target genes for designing primers. Primers are designed at the same target region (E gene) of SARS-CoV-2 genome and bat-associated SARS-related viral genes for pan-sarbecovirus detection. Specific primer probes can be designed for specific target gene regions of SARS-CoV-2 (ORF 1ab and N gene) for confirmation experiments. (C) Improving the sensitivity of the RT-qPCR assay. Primers are designed according to the general principles of primers or by using software. Degenerate primers can increase the sensitivity of the reaction. When using primers to detect the target gene, if the target gene is mutated, the universal primer may amplify, but the degenerate primer can make a false negative result into a positive result (Grey heads indicate that false negative results become positive results). After considering the secondary structure of the target gene, the target gene on the stem can prevent RNA from being cleaved by RNase A, which can also increase the detection rate. NCBI Primer-BLAST verifies the specificity of the primers. In the experimental process of the designed primer probe, dual-target detection and multi-target detection can also improve the sensitivity of the reaction. The probe is labelled with different fluorescent dyes, or each target gene is detected in a separate assay with the probes labelled with the same fluorescent dye.
As RT-qPCR is the current gold standard for the etiological diagnosis of SARS-CoV-2 infection, the diagnostic accuracy of this technique should be a foremost prerequisite. This review provides an overview of the most crucial components of a RT-qPCR assay-primer design (summarised in Figure 4). First, it is important to identify the reference sequence using a step-by-step process, and eventually align the sequence, if required. The gene targets include the N, E, ORF 1ab (or RdRP), and S genes. The more conserved E gene is the target for the pan-coronavirus assay, while N and RdRP genes mainly function as targets for confirmatory assays. After considering the basic rules of primer/probe design, and degenerate primers/probes against SARS-CoV can be used as a positive control; moreover, primers that account for the variation in coronavirus sequences can be developed. The single-target method is prone to lead to false negative results once a site in the target sequence is mutated. Thus, dual-target or multi-target detection can be additionally used to improve detection sensitivity. After multiple pairs of primers are designed, regardless of whether they are designed manually or by software, further verification is necessary. In silico analysis of primer and probe sequences is needed to analyse the specificity. However, the experimental conformation is more important. Only by evaluating the diagnostic ability of the designed primers and probes on a number of samples, we can identify primers/probes with good sensitivity and specificity. Furthermore, the reaction conditions of the assay, such as the concentration of the primers and probes, Tm, and annealing temperature, must be optimized for maximum accuracy of the results.
SARS-CoV-2: severe acute respiratory syndrome coronavirus 2; MERS-CoV: Middle East respiratory syndrome coronavirus; SARS-CoV: severe acute respiratory syndrome coronavirus; ORF: Open reading frame; S: Spike; E: Envelope; M: Membrane; N: Nucleocapsid; RT-qPCR: Quantitative real-time PCR; RdRp: RNA-dependent RNA polymerase; BALF: Bronchoalveolar lavage fluid; GISAID: Global Initiative on Sharing All Influenza Data; FAM: 6-carboxyfluorescein; BBQ: Blackberry quencher; BHQ-1: Black hole quencher 1; TAMRA: Carboxytetramethylrhodamine; IBFQ: Iowa Black FQ; CDC: Centers for Disease Control and Prevention; EVAg: European virus archive global.
This work was supported by the “AIDS and Hepatitis, and Other Major Infectious Disease Control and Prevention” Program of China under Grant [No. 2018ZX10102001]. The funder had no role in the study design, data collection, analysis, interpretation, or writing of the paper.
The authors have declared that no competing interest exists.
1. Coronaviridae Study Group of the International Committee on Taxonomy of V. The species Severe acute respiratory syndrome-related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat Microbiol. 2020 10.1038/s41564-020-0695-z
2. WHO. Coronavirus disease 2019 (COVID-19) Situation Report - 113 2020
3. Bustin SA, Nolan T. RT-qPCR Testing of SARS-CoV-2: A Primer. International journal of molecular sciences. 2020 21
4. Kandeel M, Ibrahim AA, Fayez M, Al-Nazawi M. From SARS and MERS CoVs to SARS-CoV-2: moving toward more biased codon usage in viral structural and non-structural genes. J Med Virol. 2020 10.1002/jmv.25754
5. Xu X, Chen P, Wang J, Feng J, Zhou H, Li X. et al. Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission. Science China Life Sciences. 2020;63:457-60
6. Huang C, Wang Y, Li X, Ren L, Zhao J, Hu Y. et al. Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China. Lancet. 2020;395:497-506
7. Li Q, Guan X, Wu P, Wang X, Zhou L, Tong Y. et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia. N Engl J Med. 2020;382:1199-207
8. S K, G S, M L, C K. biology TKJM, evolution. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. 2018;35:1547-9
9. Shannon A, Le NT-T, Selisko B, Eydoux C, Alvarez K, Guillemot J-C. et al. Remdesivir and SARS-CoV-2: Structural requirements at both nsp12 RdRp and nsp14 Exonuclease active-sites. Antiviral Res. 2020 178
10. sciences EAJL. Ribavirin, Remdesivir, Sofosbuvir, Galidesivir, and Tenofovir against SARS-CoV-2 RNA dependent RNA polymerase (RdRp): A molecular docking study. 2020; 117592.
11. Su S, Wong G, Shi W, Liu J, Lai ACK, Zhou J. et al. Epidemiology, Genetic Recombination, and Pathogenesis of Coronaviruses. Trends Microbiol. 2016;24:490-502
12. Guan W-D, Chen L-P, Ye F, Ye D, Wu S-G, Zhou H-X. et al. High-throughput sequencing for confirmation of suspected 2019-nCoV infection identified by fluorescence quantitative polymerase chain reaction. Chin Med J (Engl). 2020 10.1097/CM9.0000000000000792
13. Zhu N, Zhang D, Wang W, Li X, Yang B, Song J. et al. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N Engl J Med. 2020;382:727-33
14. Chu DKW, Pan Y, Cheng SMS, Hui KPY, Krishnan P, Liu Y. et al. Molecular Diagnosis of a Novel Coronavirus (2019-nCoV) Causing an Outbreak of Pneumonia. Clin Chem. 2020 hvaa029
15. Chen L, Liu W, Zhang Q, Xu K, Ye G, Wu W. et al. RNA based mNGS approach identifies a novel human coronavirus from two individual pneumonia cases in 2019 Wuhan outbreak. Emerg Microbes Infect. 2020;9:313-9
16. Che XY, Qiu LW, Liao ZY, Wang YD, Wen K, Pan YX. et al. Antigenic cross-reactivity between severe acute respiratory syndrome-associated coronavirus and human coronaviruses 229E and OC43. The Journal of infectious diseases. 2005;191:2033-7
17. Lu R, Zhao X, Li J, Niu P, Yang B, Wu H. et al. Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet. 2020;395:565-74
18. Wu F, Zhao S, Yu B, Chen YM, Wang W, Song ZG. et al. A new coronavirus associated with human respiratory disease in China. Nature. 2020;579:265-9
19. Chan JF, Yuan S, Kok KH, To KK, Chu H, Yang J. et al. A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster. Lancet. 2020;395:514-23
20. Drexler JF, Kupfer B, Petersen N, Grotto RM, Rodrigues SM, Grywna K. et al. A novel diagnostic target in the hepatitis C virus genome. PLoS medicine. 2009;6:e31
21. Liu R, Han H, Liu F, Lv Z, Wu K, Liu Y. et al. 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. 2020 S0009-8981(20)30112-1
22. Xu J, Zhao S, Teng T, Abdalla AE, Zhu W, Xie L. et al. Systematic Comparison of Two Animal-to-Human Transmitted Human Coronaviruses: SARS-CoV-2 and SARS-CoV. Viruses. 2020;12:E244
23. Cobb BR, Vaks JE, Do T, Vilchez RA. Evolution in the sensitivity of quantitative HIV-1 viral load tests. J Clin Virol. 2011;52(Suppl 1):S77-82
24. Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M. et al. The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments. Clin Chem. 2009;55:611-22
25. van Elden LJ, van Loon AM, van Alphen F, Hendriksen KA, Hoepelman AI, van Kraaij MG. et al. Frequent detection of human coronaviruses in clinical specimens from patients with respiratory tract infection by use of a novel real-time reverse-transcriptase polymerase chain reaction. The Journal of infectious diseases. 2004;189:652-7
26. Pfefferle S, Reucher S, Nörz D, Lütgehetmann M. Evaluation of a quantitative RT-PCR assay for the detection of the emerging coronavirus SARS-CoV-2 using a high throughput system. Euro Surveill. 2020;25:10.2807 /1560-7917.ES.2020.25.9.2000152
27. Xie X, Zhong Z, Zhao W, Zheng C, Wang F, Liu J. Chest CT for Typical 2019-nCoV Pneumonia: Relationship to Negative RT-PCR Testing. Radiology. 2020 200343-
28. Zhao W, Zhong Z, Xie X, Yu Q, Liu J. CT Scans of Patients with 2019 Novel Coronavirus (COVID-19) Pneumonia. Theranostics. 2020;10:4606-13
29. Liu F, Zhang Q, Huang C, Shi C, Wang L, Shi N. et al. CT quantification of pneumonia lesions in early days predicts progression to severe illness in a cohort of COVID-19 patients. Theranostics. 2020;10:5613-22
30. Fang Y, Zhang H, Xie J, Lin M, Ying L, Pang P. et al. Sensitivity of Chest CT for COVID-19: Comparison to RT-PCR. Radiology. 2020 200432
31. Zou L, Ruan F, Huang M, Liang L, Huang H, Hong Z. et al. SARS-CoV-2 Viral Load in Upper Respiratory Specimens of Infected Patients. N Engl J Med. 2020 10.1056/NEJMc2001737
32. Y S, bulletin MJJEsbEslmtEcd. GISAID: Global initiative on sharing all influenza data - from vision to reality. 2017; 22.
33. Corman VM, Landt O, Kaiser M, Molenkamp R, Meijer A, Chu DKW. et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Euro Surveill. 2020;25:2000045
34. Won J, Lee S, Park M, Kim TY, Park MG, Choi BY. et al. Development of a Laboratory-safe and Low-cost Detection Protocol for SARS-CoV-2 of the Coronavirus Disease 2019 (COVID-19). Exp Neurobiol. 2020 10.5607/en20009
35. Park WB, Kwon NJ, Choi SJ, Kang CK, Choe PG, Kim JY. et al. Virus Isolation from the First Patient with SARS-CoV-2 in Korea. J Korean Med Sci. 2020;35:e84-e
36. Wang C, Liu Z, Chen Z, Huang X, Xu M, He T. et al. The establishment of reference sequence for SARS-CoV-2 and variation analysis. J Med Virol. 2020 10.1002/jmv.25762
37. Shen Z, Xiao Y, Kang L, Ma W, Shi L, Zhang L. et al. Genomic diversity of SARS-CoV-2 in Coronavirus Disease 2019 patients. Clin Infect Dis. 2020: ciaa203.
38. Li K, Brownley A. Primer design for RT-PCR. Methods in molecular biology (Clifton, NJ). 2010;630:271-99
39. Campos MJ, Quesada A. Strategies to Improve Efficiency and Specificity of Degenerate Primers in PCR. Methods in molecular biology (Clifton, NJ). 2017;1620:75-85
40. Yi H. 2019 novel coronavirus is undergoing active recombination. Clin Infect Dis. 2020
41. Lippi G, Simundic AM, Plebani M. Potential preanalytical and analytical vulnerabilities in the laboratory diagnosis of coronavirus disease 2019 (COVID-19). Clinical chemistry and laboratory medicine. 2020
42. Institut Pasteur P. Real-time RT-PCR assays for the detection of SARS-CoV-2. 2020.
43. China_CDC. Specific primers and probes for detection 2019 novel coronavirus. 2020.
44. US_CDC. 2019-Novel Coronavirus (2019-nCoV) Real-time rRT-PCR Panel: Primers and Probes. 2020.
45. Nao; N, Shirato; K, Katano; H, Matsuyama; S, Takeda M. Detection of second case of 2019-nCoV infection in Japan. 2020.
46. Health TMoP. Diagnostic detection of Novel coronavirus 2019 by Real time RT- PCR. 2020.
47. Nalla AK, Casto AM, Huang MW, Perchetti GA, Sampoleo R, Shrestha L. et al. Comparative Performance of SARS-CoV-2 Detection Assays using Seven Different Primer/Probe Sets and One Assay Kit. J Clin Microbiol. 2020
48. M W, M Y, H X, W L, B K, B Z. et al. SARS-CoV infection in a restaurant from palm civet. 2005; 11: 1860-5.
49. Dennis Lo YM, Chiu RWK. Racing towards the development of diagnostics for a novel coronavirus (2019-nCoV). Clin Chem. 2020 hvaa038
50. JL R, CM P, EA G, X dL, R C, B C. et al. The European Virus Archive goes global: A growing resource for research. 2018; 158: 127-34.
51. Wang P, Anderson N, Pan Y, Poon L, Charlton C, Zelyas N. et al. The SARS-CoV-2 Outbreak: Diagnosis, Infection Prevention, and Public Perception. Clin Chem. 2020 hvaa080
52. Jung YJ, Park G-S, Moon JH, Ku K, Beak S-H, Kim S. et al. Comparative analysis of primer-probe sets for the laboratory confirmation of SARS-CoV-2. bioRxiv 20200225964775. 2020
53. Bustin S, Huggett J. qPCR primer design revisited. Biomolecular detection and quantification. 2017;14:19-28
54. Vijgen L, Moës E, Keyaerts E, Li S, Van Ranst M. A pancoronavirus RT-PCR assay for detection of all known coronaviruses. Methods in molecular biology (Clifton, NJ). 2008;454:3-12
55. Chudy M, Weber-Schehl M, Pichl L, Jork C, Kress J, Heiden M. et al. Blood screening nucleic acid amplification tests for human immunodeficiency virus Type 1 may require two different amplification targets. Transfusion. 2012;52:431-9
56. Lee HK, Lee CK, Loh TP, Tang JW, Chiu L, Tambyah PA. et al. Diagnostic testing for pandemic influenza in Singapore: a novel dual-gene quantitative real-time RT-PCR for the detection of influenza A/H1N1/2009. The Journal of molecular diagnostics: JMD. 2010;12:636-43
57. Bellau-Pujol S, Vabret A, Legrand L, Dina J, Gouarin S, Petitjean-Lecherbonnier J. et al. Development of three multiplex RT-PCR assays for the detection of 12 respiratory RNA viruses. J Virol Methods. 2005;126:53-63
58. Chuang LY, Cheng YH, Yang CH. Specific primer design for the polymerase chain reaction. Biotechnology letters. 2013;35:1541-9
59. Gunson RN, Collins TC, Carman WF. Practical experience of high throughput real time PCR in the routine diagnostic virology setting. J Clin Virol. 2006;35:355-67
60. Kamau E, Agoti CN, Lewa CS, Oketch J, Owor BE, Otieno GP. et al. Recent sequence variation in probe binding site affected detection of respiratory syncytial virus group B by real-time RT-PCR. J Clin Virol. 2017;88:21-5
61. Cha RS, Zarbl H, Keohavong P, Thilly WG. Mismatch amplification mutation assay (MAMA): application to the c-H-ras gene. PCR methods and applications. 1992;2:14-20
62. Crockett AO, Wittwer CT. Fluorescein-labeled oligonucleotides for real-time pcr: using the inherent quenching of deoxyguanosine nucleotides. Analytical biochemistry. 2001;290:89-97
63. Chen L, Li W, Zhang K, Zhang R, Lu T, Hao M. et al. Hepatitis C Virus RNA Real-Time Quantitative RT-PCR Method Based on a New Primer Design Strategy. The Journal of molecular diagnostics: JMD. 2016;18:84-91
64. Prats-Ejarque G, Blanco JA, Salazar VA, Nogués VM, Moussaoui M, Boix E. Characterization of an RNase with two catalytic centers. Human RNase6 catalytic and phosphate-binding site arrangement favors the endonuclease cleavage of polymeric substrates. Biochimica et biophysica acta General subjects. 2019;1863:105-17
65. Villegas J, Burzio V, Villota C, Landerer E, Martinez R, Santander M. et al. Expression of a novel non-coding mitochondrial RNA in human proliferating cells. Nucleic acids research. 2007;35:7336-47
66. Smith RM, Walton CM, Wu CH, Wu GY. Secondary structure and hybridization accessibility of hepatitis C virus 3'-terminal sequences. J Virol. 2002;76:9563-74
67. Baranauskas A, Paliksa S, Alzbutas G, Vaitkevicius M, Lubiene J, Letukiene V. et al. Generation and characterization of new highly thermostable and processive M-MuLV reverse transcriptase variants. Protein engineering, design & selection: PEDS. 2012;25:657-68
68. Thornton B, Basu C. Rapid and simple method of qPCR primer design. Methods in molecular biology (Clifton, NJ). 2015;1275:173-9
69. Li K, Brownley A, Stockwell TB, Beeson K, McIntosh TC, Busam D. et al. Novel computational methods for increasing PCR primer design effectiveness in directed sequencing. BMC bioinformatics. 2008;9:191
70. Rozen S, Skaletsky H. Primer3 on the WWW for general users and for biologist programmers. Methods in molecular biology (Clifton, NJ). 2000;132:365-86
71. Chen SH, Lin CY, Cho CS, Lo CZ, Hsiung CA. Primer Design Assistant (PDA): A web-based primer design tool. Nucleic acids research. 2003;31:3751-4
72. V G, A V, chemistry ASJTJob. RNase P: variations and uses. 2002; 277: 6759-62.
73. Laboratory testing for coronavirus disease 2019 (COVID-19) in suspected human cases. 2020
Corresponding author: Jinming Li, National Center for Clinical Laboratories, No.1 Dahua Road, Dongdan, Beijing 100730, People's Republic of China. Phone: 86 10 58115053; Fax: 86 10 65212064; E-mail: jmliorg.cn.