Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Reliability of Spike Gene Target Failure for ascertaining SARS-CoV-2 lineage B.1.1.7 prevalence in a hospital setting

View ORCID ProfileJosé Afonso Guerra-Assunção, Paul A. Randell, Florencia A. T. Boshier, View ORCID ProfileMichael A. Crone, View ORCID ProfileJuanita Pang, Tabitha Mahungu, Paul S. Freemont, Judith Breuer
doi: https://doi.org/10.1101/2021.04.12.21255084
José Afonso Guerra-Assunção
1Great Ormond Street Institute of Child Health, University College London (UCL), London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for José Afonso Guerra-Assunção
  • For correspondence: a.guerra@ucl.ac.uk
Paul A. Randell
2North West London Pathology, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Florencia A. T. Boshier
1Great Ormond Street Institute of Child Health, University College London (UCL), London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Michael A. Crone
3Department of Infectious Disease, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael A. Crone
Juanita Pang
1Great Ormond Street Institute of Child Health, University College London (UCL), London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Juanita Pang
Tabitha Mahungu
4Department of Infection, Royal Free London NHS Foundation Trust, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Paul S. Freemont
3Department of Infectious Disease, Imperial College London, London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Judith Breuer
1Great Ormond Street Institute of Child Health, University College London (UCL), London, UK
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Data/Code
  • Preview PDF
Loading

Abstract

The appearance of the SARS-CoV-2 lineage B.1.1.7 in the UK in late 2020, associated with faster transmission, sparked the need to find effective ways to monitor its spread. The set of mutations that characterise this lineage include a deletion in position 69 and 70 of the spike protein, which is known to be associated with Spike Gene Target Failure (SGTF) in a commonly used three gene diagnostic qPCR assay. The lower cost and faster turnaround times compared to whole genome sequencing make the use of qPCR for monitoring of the variant spread an attractive proposition. However, there are several potential issues with this approach. Here we use 826 SARS-CoV-2 samples collected in a hospital setting as part of the Hospital Onset COVID Infection (HOCI) study where qPCR was used for viral detection, followed by whole genome sequencing (WGS), to identify the factors to consider when using SGTF to infer lineage B.1.1.7 prevalence in a hospital setting, with potential implications for locations where this variant has recently been introduced.

Introduction

The emergence of the UK Variant of Concern (VOC) 202012/01, also known as lineage B.1.1.7, in South East England during the latter part of 2020 has been associated with rising rates of transmission [1] and potentially with increased disease severity [2].

Increased prevalence of this variant is correlated with an increase in the number of qPCR-based community diagnostic tests that fail to detect spike (S) gene amplicons [3]. So called spike gene target failure (SGTF) has been shown to be due to a deletion of amino acids 69 and 70 in the spike protein leading to failure of probes to bind to the S gene amplicon in one commercial qPCR assay. Detection of other amplicon targets, including the nucleocapsid (N) and ORF1ab genes is unaffected. In community studies SGTF shows good correlation of SGTF with lineage B.1.1.7 after mid-November, but is less accurate before that date [3].

To determine how well SGTF corresponded to VOC in patients hospitalised over the same period we made use of data collected as part of the COVID-19 UK Hospital Onset COVID Infection (COG-UK-HOCI). Only one of 15 hospitals in this trial is using an assay that involves S gene target detection in inpatients while another uses it for healthcare workers only. In this dataset we present data on 826 samples collected between weeks 41 and 53 of 2020. This includes 535 hospital patients collected from week 41, and 291 healthcare workers collected from week 49.

Results

An initial analysis based on unreported S gene results was undertaken. The numbers of samples with and without SGTF from two laboratories are shown in Table 1 and the increasing proportion of SGTF, over time from weeks 45-53 is shown in Figure 1. Among samples with SGTF the proportion due to lineage B.1.1.7, as determined by sequencing, started at 0% in week 43, rising to 20% by week 46 and to nearly 100% by week 53 (Figure S1 and S2). The majority (37/49) of SGTF that were not lineage B.1.1.7 were lineage B.1.258, which is known to carry the same 69/70 two amino acid deletion. SGTF was also observed in 12 samples of other lineages that do not have 69/70 S gene deletion (Table 2).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 1

Number of samples from each laboratory with a breakdown by lineage status as identified by pangolin.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table 2

Breakdown of samples that demonstrate SGTF in this dataset by pangolin lineage and presence of deletion of spike amino acids 69 and 70.

Figure 1
  • Download figure
  • Open in new tab
Figure 1

Number of samples sequenced in each week (A) for each of the hospital trusts in this study. Samples that are lineage B.1.1.7 are shown in blue, while samples from other lineages are shown in green. Samples that demonstrate SGTF are shown in a lighter colour than those that do not demonstrate SGTF. Panel B shows the fraction of samples in each category.

Surprisingly, an apparent signal from the S-gene (based on the autogenerated CT values) was present in 17% (24 out of 143) of samples later confirmed as lineage B.1.1.7s by sequencing, all of which had deletions at amino acids 69/70 (Figure 1 and S3). Closer scrutiny of the amplification curves revealed that in these samples the S-gene curves were not consistent with genuine amplification curves as demonstrated in Figure 2. Although the curves crossed the auto-threshold, these would be reported as S-gene not detected if individual targets were being reported. Without manual curation these samples were falsely identified as not lineage B.1.1.7, resulting in underestimation of the true lineage B.1.1.7 prevalence. To further investigate this phenomenon, we compared the qPCR CT values for each of the three genes for lineage B.1.1.7 with apparent SGTF with those samples demonstrating the anticipated SGTF (Figure 3). Samples with the spurious signal for the S gene had significantly lower N and ORF1ab CT values (compared to samples with SGTF) and within this group the S gene CT values were always significantly higher than the CT values for the other genes (p<0.001).

Figure 2
  • Download figure
  • Open in new tab
Figure 2

Amplification curves obtained for a subset of samples in this study. Each panel corresponds to a different gene: A) ORF1 ab, B) N gene, C) S gene, where some flattened amplification curves that cross the automatic threshold for detection can be seen.

Figure 3
  • Download figure
  • Open in new tab
Figure 3

Comparison between the CT values of lineage B.1.1.7 samples that demonstrate SGTF in comparison to those that do not. Each gene in the assay is represented by a different boxplot (N in red, ORF1ab in blue, S in green). The number of samples in this dataset for each category is shown above each boxplot pair. The RFH health trust has no lineage B.1.1.7 samples that do not show SGTF. In the IHT samples, it is clear that the CT value of the S gene is significantly higher than the other two genes in the assay (p < 0.001).

To investigate whether the viral load was a factor in this phenomenon a subset of samples including lineage B.1.1.7 samples that had previously not shown SGTF, samples that had shown SGTF as expected and samples from other lineages without the relevant deletion were retested neat and when serially diluted. The spurious S gene signal was not reliably reproducible and did not appear to be related to the viral load (data not shown).

Discussion

SGTF has been used to identify lineage B.1.1.7 in the community and forms the basis for many studies looking at its transmission and severity [1,4–6]. Our data suggest that both false positive and false negatives can occur which may skew the positive predictive value of SGTF. In this study, lineage B.1.258, which also carries the spike deletion 69/70, was prevalent locally and we observed false positive lineage B.1.1.7 calls in 9.6% of cases. The incidence of false positives decreased as the more transmissible B.1.1.7 outcompeted B.1.258. By week 51 more than 90% of SGTF was B.1.1.7 (Figure S2).

We also found that there is a risk that spurious S gene results may be observed (as was the case with one of the assays used) and without careful review and consideration of the PCR results the prevalence of SGTF (and consequently lineage B.1.1.7) may be underestimated. This could be due to imperfect colour compensation on qPCR machines other than those made by the manufacturer of the assay, further complicated by the lack of controls for each single target of the multiplex assay to check for colour bleedthrough. The qPCR data from the manufacturer’s machine is also interpreted by separated software developed by the manufacturer and it is unclear how thresholding is performed and how the software handles any noise or curves that do not have an exponential shape.

This phenomenon is evidently dependent on the PCR machine and dyes used. This is reflected in the fact that the false negative results are restricted to just one of the two sites studied (one site used the manufacturer’s qPCR machine, the other a different machine with colour compensation performed with the manufacturer’s calibration plates).

While this phenomenon was only observed in one of the two set of samples analysed here, it has also been independently observed for lineage B.1.1.7 and SGTF in Portugal. [7]. Additionally, a locally developed assay that mimicked SGTF with a different dye did not show the same behaviour (data not shown).

Thus, in summary, SGTF is an important surrogate marker for the VOC 202012/01, and useful for large scale epidemiology studies. However, care needs to be taken where sample numbers are small, for example in care homes or hospitals for example where the intention is to link lineage B.1.1.7 to outcomes or in scenarios where lineage B.1.1.7 variants are in the minority. In these cases, the use of SGTF can be misleading and sequencing is recommended.

Data Availability

Sequencing data was deposited when generated to the appropriate SARS-CoV-2 public repositories (COG-UK and GISAID) and is publicly available.

Materials and Methods

Sequencing

Whole genome sequences for SARS-CoV-2 were generated following a positive qPCR test using either Illumina or ONT nanopore technologies, according with COG-UK [8] processes and deposited in the appropriate repositories. Lineages were determined using Pangolin [9]. The deletion of Spike amino acids 69 and 70 was independently inspected after alignment to reference using minimap2 [10]. The data was aggregated, analysed and visualised with the R statistical framework.

RFH qPCR method

Viral samples are inactivated using RNA Inactivation buffer, and RNA is extracted using an automated in-house method, based on a protocol previously described [11].

The extracted RNA undergoes a RT-qPCR assay and is quantified on a ABI QuantStudio5 system together with a positive control (TaqPath COVID-19 Control) and an extraction control.

Presence of viral targets is indicated by an S shape curve with CT < 37, together with concordant duplicates and the presence of the MS2 control trace at CT < 32

IHT qPCR methods

A sample volume of 200 μl was used for RNA extraction using the Maxwell HT Viral TNA kit (Promega) with a custom extraction protocol [12] on the CyBio FeliX liquid handler (Analytik Jena) with an elution volume of 50 μl. Subsequent RT-qPCR was performed using the TaqPath™ COVID-19 CE- IVD RT-PCR Kit (ThermoFisher Scientific) according to the manufacturer’s instructions and thermocycled on a qTower3 (Analytik Jena). Colour compensation was performed using FAM, VIC, ABY and JUN dyes (ThermoFisher Scientific) according to the qTower3’s calibration instructions. When determining the CT values for the different targets the auto-threshold generated by the analyser was used. Samples were reported as SARS CoV-2 RNA detected when at least 2 targets were detected with typical exponential growth curves. Samples not meeting these criteria were reported as not detected if no targets were present or underwent confirmatory testing if there was only one target. As the overall result was reported (and not individual targets) the autogenerated CT values were not manually curated.

Supplementary Figures

Figure S1
  • Download figure
  • Open in new tab
Figure S1

Total number of samples in this study stratified by week of sample collection and laboratory where they originated from. (A) shows the total number of samples for each week. (B) illustrates the fraction of samples that demonstrate SGTF (blue) or do not demonstrate SGTF (red) across time. (C) illustrates the fraction of samples that were classified by pangolin as lineage B.1.1.7 (blue) or other lineage (red) across time.

Figure S2
  • Download figure
  • Open in new tab
Figure S2

Number of SGTF instances by week (A) for each of the hospital trusts in this study, and fraction of SGTF samples per week that are lineage B.1.1.7 (B). SGTF events are separated between those samples that are not lineage B.1.1.7 (red) and those that are lineage B.1.1.7 (blue).

Figure S3
  • Download figure
  • Open in new tab
Figure S3

Number of lineage B.1.1.7 samples by week (A) for each of the hospital trusts in this study, and fraction of lineage B.1.1.7 samples per week that cause SGTF (B). Lineage B.1.1.7 events are separated between those samples that did not cause SGTF (red) and those that have undetectable S gene (blue).

Acknowledgements

Data for this study was collected as part of the COG-UK HOCI study. COG-UK HOCI is part of COG-UK. COG-UK is supported by funding from the Medical Research Council (MRC) part of UK Research & Innovation (UKRI), the National Institute of Health Research (NIHR) and Genome Research Limited, operating as the Wellcome Sanger Institute.

References

  1. [1].↵
    E. Volz, S. Mishra, M. Chand, J.C. Barrett, R. Johnson, L. Geidelberg, W.R. Hinsley, D.J. Laydon, G. Dabrera, Á. O’Toole, R. Amato, M. Ragonnet-Cronin, I. Harrison, B. Jackson, C.V. Ariani, O. Boyd, N.J. Loman, J.T. McCrone, S. Gonçalves, D. Jorgensen, R. Myers, V. Hill, D.K. Jackson, K. Gaythorpe, N. Groves, J. Sillitoe, D.P. Kwiatkowski, T.C.-19 G.U. ( COG-U. Consortium, S. Flaxman, O. Ratmann, S. Bhatt, S. Hopkins, A. Gandy, A. Rambaut, N.M. Ferguson, Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data, MedRxiv. (2021) 2020.12.30.20249034. https://doi.org/10.1101/2020.12.30.20249034.
  2. [2].↵
    R. Challen, E. Brooks-Pollock, J.M. Read, L. Dyson, K. Tsaneva-Atanasova, L. Danon, Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: matched cohort study, BMJ. 372 (2021) 579. https://doi.org/10.1136/bmj.n579.
  3. [3].↵
    PHE: Investigation of novel SARS-CoV-2 variant of concern 202012/01 - technical briefing 4, 13 January 2021, GOV.UK. (2021). https://www.gov.uk/government/publications/phe-investigation-of-novel-sars-cov-2-variant-of-concern-20201201-technical-briefing-4-13-january-2021 (Accessed February 4, 2021).
  4. [4].↵
    N.G. Davies, S. Abbott, R.C. Barnard, C.I. Jarvis, A.J. Kucharski, J.D. Munday, C.A.B. Pearson, T.W. Russell, D.C. Tully, A.D. Washburne, T. Wenseleers, A. Gimma, W. Waites, K.L.M. Wong, K. van Zandvoort, J.D. Silverman, C.C.-19 W. Group1‡, C.-19 G.U. (COG-U. Consortium‡, K. Diaz-Ordaz, R. Keogh, R.M. Eggo, S. Funk, M. Jit, K.E. Atkins, W.J. Edmunds, Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England, Science. (2021). https://doi.org/10.1126/science.abg3055.
  5. [5].
    M. Patone, K. Thomas, R. Hatch, P.S. Tan, C. Coupland, W. Liao, P. Mouncey, D. Harrison, K. Rowan, P. Horby, P. Watkinson, J. Hippisley-Cox, Analysis of severe outcomes associated with the SARS-CoV-2 Variant of Concern 202012/01 in England using ICNARC Case Mix Programme and QResearch databases, MedRxiv. (2021) 2021.03.11.21253364. https://doi.org/10.1101/2021.03.11.21253364.
  6. [6].↵
    N.G. Davies, C.I. Jarvis, C.C.-19 W. Group, W.J. Edmunds, N.P. Jewell, K. Diaz-Ordaz, R.H. Keogh, Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7, MedRxiv. (2021) 2021.02.01.21250959. https://doi.org/10.1101/2021.02.01.21250959.
  7. [7].↵
    Tracking SARS-CoV-2 VOC 202012/01 (lineage B.1.1.7) dissemination in Portugal: insights from nationwide RT-PCR Spike gene drop out data - SARS-CoV-2 coronavirus, Virological. (2021). https://virological.org/t/tracking-sars-cov-2-voc-202012-01-lineage-b-1-1-7-dissemination-in-portugal-insights-from-nationwide-rt-pcr-spike-gene-drop-out-data/600 (Accessed March 2, 2021).
  8. [8].↵
    An integrated national scale SARS-CoV-2 genomic surveillance network, Lancet Microbe. 1 (2020) e99–e100. https://doi.org/10.1016/S2666-5247(20)30054-9.
    OpenUrl
  9. [9].↵
    A. Rambaut, E.C. Holmes, Á. O’Toole, V. Hill, J.T. McCrone, C. Ruis, L. du Plessis, O.G. Pybus, A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology, Nat. Microbiol. 5 (2020) 1403–1407. https://doi.org/10.1038/s41564-020-0770-5.
    OpenUrl
  10. [10].↵
    H. Li, Minimap2: pairwise alignment for nucleotide sequences, Bioinformatics. 34 (2018) 3094–3100. https://doi.org/10.1093/bioinformatics/bty191.
    OpenUrlCrossRefPubMed
  11. [11].↵
    R. Boom, C.J. Sol, M.M. Salimans, C.L. Jansen, P.M. Wertheim-van Dillen, J. van der Noordaa, Rapid and simple method for purification of nucleic acids, J. Clin. Microbiol. 28 (1990) 495–503. https://doi.org/10.1128/JCM.28.3.495-503.1990.
    OpenUrlAbstract/FREE Full Text
  12. [12].↵
    M.A. Crone, M. Priestman, M. Ciechonska, K. Jensen, D.J. Sharp, A. Anand, P. Randell, M. Storch, P.S. Freemont, A role for Biofoundries in rapid development and validation of automated SARS-CoV-2 clinical diagnostics, Nat. Commun. 11 (2020) 4464. https://doi.org/10.1038/s41467-020-18130-3.
    OpenUrl
Back to top
PreviousNext
Posted April 14, 2021.
Download PDF
Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Reliability of Spike Gene Target Failure for ascertaining SARS-CoV-2 lineage B.1.1.7 prevalence in a hospital setting
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Reliability of Spike Gene Target Failure for ascertaining SARS-CoV-2 lineage B.1.1.7 prevalence in a hospital setting
José Afonso Guerra-Assunção, Paul A. Randell, Florencia A. T. Boshier, Michael A. Crone, Juanita Pang, Tabitha Mahungu, Paul S. Freemont, Judith Breuer
medRxiv 2021.04.12.21255084; doi: https://doi.org/10.1101/2021.04.12.21255084
Digg logo Reddit logo Twitter logo Facebook logo Google logo LinkedIn logo Mendeley logo
Citation Tools
Reliability of Spike Gene Target Failure for ascertaining SARS-CoV-2 lineage B.1.1.7 prevalence in a hospital setting
José Afonso Guerra-Assunção, Paul A. Randell, Florencia A. T. Boshier, Michael A. Crone, Juanita Pang, Tabitha Mahungu, Paul S. Freemont, Judith Breuer
medRxiv 2021.04.12.21255084; doi: https://doi.org/10.1101/2021.04.12.21255084

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Infectious Diseases (except HIV/AIDS)
Subject Areas
All Articles
  • Addiction Medicine (164)
  • Allergy and Immunology (417)
  • Anesthesia (93)
  • Cardiovascular Medicine (867)
  • Dentistry and Oral Medicine (159)
  • Dermatology (98)
  • Emergency Medicine (251)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (398)
  • Epidemiology (8597)
  • Forensic Medicine (4)
  • Gastroenterology (391)
  • Genetic and Genomic Medicine (1775)
  • Geriatric Medicine (170)
  • Health Economics (376)
  • Health Informatics (1252)
  • Health Policy (625)
  • Health Systems and Quality Improvement (472)
  • Hematology (198)
  • HIV/AIDS (380)
  • Infectious Diseases (except HIV/AIDS) (10354)
  • Intensive Care and Critical Care Medicine (554)
  • Medical Education (193)
  • Medical Ethics (51)
  • Nephrology (214)
  • Neurology (1692)
  • Nursing (97)
  • Nutrition (252)
  • Obstetrics and Gynecology (330)
  • Occupational and Environmental Health (451)
  • Oncology (934)
  • Ophthalmology (265)
  • Orthopedics (104)
  • Otolaryngology (172)
  • Pain Medicine (115)
  • Palliative Medicine (40)
  • Pathology (256)
  • Pediatrics (541)
  • Pharmacology and Therapeutics (257)
  • Primary Care Research (210)
  • Psychiatry and Clinical Psychology (1788)
  • Public and Global Health (3877)
  • Radiology and Imaging (629)
  • Rehabilitation Medicine and Physical Therapy (324)
  • Respiratory Medicine (525)
  • Rheumatology (208)
  • Sexual and Reproductive Health (171)
  • Sports Medicine (159)
  • Surgery (191)
  • Toxicology (36)
  • Transplantation (101)
  • Urology (76)