Implication of SARS-CoV-2 evolution in the sensitivity of RT-qPCR diagnostic assays.

Since the initial outbreak of COVID-19 and identification of a novel enveloped RNA betacoronavirus, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), reports suggest the virus might be evolving, albeit at a lower rate than influenza viruses.1, 2, 3 Nevertheless, in light of the rapid and pandemic-scale spread of COVID-19, mutations in SARS-CoV-2 raise new diagnostic challenges, including the redesign of the oligonucleotide sequences in use in RT-qPCR assays to circumvent potential primer–sample mismatches.

We report an analysis of all high-coverage SARS-CoV-2 genome sequences (1825 in total) deposited in the Global Initiative on Sharing All Influenza Data (GISAID) database as of March 30, 2020. We aligned the sequences against the reference sequence obtained from the Wuhan seafood market pneumonia virus isolate, Wuhan-Hu-1 (NC_045512). Subsequently, we annotated in the alignments the binding sites of 33 oligonucleotides developed by different centres and shared by WHO for use in the RT-qPCR detection of SARS-CoV-2 from human samples. We then calculated the nucleotide diversity (π) in the binding region of each oligonucleotide (figure ).

Figure

Genomic position and nucleotide diversity in the binding sites of oligonucleotides in use in RT-qPCR diagnostic assays across 1825 SARS-CoV-2 genomes isolated from humans in different countries

(A) Graphical representation of a SARS-CoV-2 genome (29 903 nucleotides, NC_045512), showing genes (blue) and oligonucleotide primer binding sites (forward primers in dark green, reverse primers in light green, and double-dye oligonucleotide probes in red). (1) Pasteur_RdRp_IP2, (2) ChinaCDC_Orf1ab, (3) Pasteur_RdRp_IP4, (4) Charité_RdRp, (5) HKU_Orf1ab, (6) Charité_E, (7) USACDC_N1, (8) USACDC_N3, (9) ChinaCDC_N, (10) HKU_N, and (11) USACDC_N2. (B) Nucleotide diversity (π) in the binding sites of the oligonucleotides in use in RT-qPCR diagnostic assays across 1825 SARS-CoV-2 genomes isolated from humans in different countries. SARS-CoV-2=severe acute respiratory syndrome coronavirus 2.

The analysis revealed that 79% (26 of 33) of the primer binding sites used in the RT-qPCR assays were mutated in at least one genome. The substitution of three nucleotides (GGG substituted to AAC) at the beginning of the binding site of the forward primer designed by the Chinese National Institute for Viral Disease Control and Prevention in the gene encoding for the nucleocapsid phosphoprotein was of relevance. The AAC variant was in 14% (258 of 1825) of the genomes isolated and sequenced in 24 different countries (Australia, Belgium, Brazil, Chile, China, Czech Republic, Denmark, England, Finland, France, Germany, Greece, Hungary, Iceland, Italy, Mexico, the Netherlands, Peru, Portugal, Spain, Switzerland, the USA, Vietnam, and Wales).

Notwithstanding the possibility of sequencing errors, the fact that some variants were consistently found in different sequencing experiments at independent test sites suggests that these genetic variants are true variants. Moreover, the observation that at least one of the previously designed primers is now likely to be ineffective at detecting up to 14% of the virus variants in circulation strengthens the need to continue optimising the oligonucleotides in use in assays being developed. Oligonucleotide optimisation will be facilitated by global sharing of SARS-CoV-2 genomes and the frequently updated reports on sequence analysis that are available on the GISAID website.