Heat inactivation of the Murine Hepatitis Virus (MHV), a surrogate for human SARS-CoV-2
Given the hazard classification of the SARS-CoV-2 as a Hazard Group 3 pathogen by the Advisory Committee on Dangerous Pathogens (ACDP), a surrogate model mouse hepatitis coronavirus (MHV-A59) [16], was chosen for viral inactivation studies. MHV is a prototype of the group II coronavirus (beta-coronavirus) in common with severe acute respiratory syndrome coronaviruses such as SARS-CoV-2 [17]. Commonly found in laboratory mice, MHV-A59 is a Hazard Group 1 pathogen, can be replicated to high titers in cell culture and has been used as an experimental model to understand the biology of coronaviruses [16-20]. Furthermore, MHV-A549 is a widely accepted coronavirus surrogate model and has been used to assess the survival of coronaviruses in the environment as well as sensitivity to biocidal agents and heat [20].
To examine the effect of heat on virus viability, we prepared high titer cell-free MHV-A549 samples in the same VTM as used in clinical samples. Samples were subjected to heat inactivation in a thermal cycler for 5, 10 or 15 minutes at 60°C, 65°C and 75°C, followed by a tissue culture infectious dose (TCID50) assay to measure sample viability. In the absence of heat treatment, high levels of infectious virus were observed in the TCID50 assay for samples containing cell-free virus of 3.5x106 TCID50/ml (average of n=3 experiments). Following heat treatment, no live virus was detected under any of the conditions examined while the control conditions confirmed virus infectivity. This demonstrates that heat alone was sufficient to render high titers of the surrogate viral samples non-infectious.
Large-scale inactivation of SARS-CoV-2 in OP/NP clinical samples from COVID-19 patients and the general public
A wet heating system from the catering industry (Electrolux Skyline Combi Oven) was identified for the bulk HI of 500 samples per hour while held within their UN3373 secondary containment. Efficient heat penetration to all samples throughout the oven was achieved by the use of a racking system comprising 25 sample bags per tray, with 10 trays per oven, loaded on a roll-in/roll-out trolley system for easy handling (Fig. 1a and Supplementary Fig. S1). Experiments were performed to optimise a heating cycle to achieve a sustained temperature of at least 65°C for 10 minutes for all of the samples with minimal overheating in a time efficient manner. Twenty temperature probes were packed inside mock OP/NP sample vials along with a swab, and packed in a manner mimicking clinical samples (Fig. 1b). These thermocouple swab vials along with and an additional 230 unused OP/NP swab test kits, assembled as if used, were placed in the oven, with the thermocouple swab vials distributed 2 per tray throughout the oven. The most optimal heating cycle in terms of achieving a minimum heat exposure of 65°C or above for at least 10 minutes in all thermocouple swab vials was identified and confirmed in three independent experiments. Using the optimised parameters (Supplementary Table 1), it took approximately 15 minutes for all samples to reach 65°C, and no sample exceeded 80°C (Fig. 1c and 1d and Supplementary Fig. S2). The overall cycle time was 30 minutes, meaning 500 samples could be subject to the bulk HI protocol per hour per oven.
Next, a range of different OP/NP vial and swab combinations were tested using the thermocouple swabs created to reflect the diversity of test kits used across the COVID-19 testing network (Supplementary Fig. S3). This showed that all of the OP/NP vial and swab types tested achieved the same required heat inactivation conditions, and did not exceed 80°C (Supplementary Fig. S3). Additionally, as the CCTC also analysed pooled OP/NP tests supplied by the University of Cambridge containing up to 10 swabs in a single tube [21], we also confirmed that vials containing multiple OP/NP swabs performed similarly to vials containing single swabs. We observed a 4 minute time difference before reaching the required 65°C, but no additional over-heating (Supplementary Fig. S4).
Having established a suitable temperature profile for HI of OP/NP swab vials in bulk using the industrial steam oven, we then confirmed viral inactivation in this precise setting using the surrogate virus MHV-A549. Clinical samples typically contain a mix a cell free virus and virus infected epithelial cells. We therefore prepared mock OP/NP swab viral samples containing either cell-free MHV-A59 viral particles or MHV-A59 infected 17Cl-1 cells in VTM. Swab vials were packaged in UN3373 secondary containment before treatment according to the bulk HI protocol in the industrial steam oven. As before, TCID50 assays confirmed the absence of infectious virus in all of the MHV-A549 samples following heat treatment. Two independent experiments with multiple biological repeats were performed and identical results were obtained (Fig. 1e). Thus, the bulk HI protocol comprising of a 30 minute cycle achieving a 10 minute hold at 65°C or greater, reduces virus infectivity by up to 6 log10 TCID50/ml.
Detection of SARS-CoV-2 in clinical OP/NP samples following heat exposure
As shown in previous experiments, efficient bulk HI of samples for a minimum of 10 minute heat treatment at 65°C required a temperature ramping programme whereby samples typically reached temperatures of 75°C to 80°C for around 20 minutes. The nature of the bulk HI process means some sample to sample variability in heat exposure is unavoidable. We therefore examined the effect of heat and duration of the heat step on assay sensitivity in clinical OP/NP swab samples. Randomly selected clinical OP/NP swab samples received at the CCTC were retained following standard testing [6] and reserved at 4°C until their RT-qPCR data was available. SARS-CoV-2 positive samples were then identified and used for experiments to examine a range of temperatures and incubation times using a thermal cycler or a small dry air oven. Samples were re-tested in the RT-qPCR assay and Cq values compared to control conditions of heat inactivation at 65°C for 10 minutes, reflecting the HI step of the standard CCTC assay protocol. The control conditions were included instead of comparison to the clinical test result to avoid time-induced changes.
An initial study was performed to understand the intrinsic variability in our assay when retesting known SARS-CoV-2 positive samples, by repeating replicate tests of a panel of samples on three consecutive days (Supplementary Fig. S5). Using an analysis of variance to estimate the pooled within-sample variability we calculated a variance of 0.33 Cq. This value was used as a reference point to determine whether treatment-based differences of subsequent experiments were within or outside the expected assay variability.
We then compared incubation temperatures of 75°C and 80°C for 10 minutes using SARS-CoV-2 positive clinical samples chosen in the Cq range of 17 to 33. Incubation at 75°C for 10 minutes had little effect on SARS-CoV-2 detection, with a 0.39 +/- 0.4 Cq variation which is in line with the known assay variability of 0.33 Cq (Fig. 2a and Supplementary Table 2). Incubation of the same samples at 80°C for 10 minutes however resulted in a small loss in sensitivity (i.e. a higher Cq) for SARS-CoV-2 detection by RT-qPCR with a mean increase in Cq value of 0.69 +/- 0.55 (Fig. 2a and Supplementary Table 2). Within the OP/NP swab samples tested at different temperatures, some were more affected than others. This variation was independent of the level of SARS-CoV-2 detection (i.e. the Cq value). For example in Fig. 2b, sample O28 and O16 (initial Cq values of 18 and 32 respectively) show no effect of heat, while samples O27 and R8 (initial Cq values of 24 and 33 respectively) show a 2 Cq difference in SARS Cov-2 detection.
Next, we evaluated the effect of incubation time at a constant temperature of 65°C using, SARS-CoV-2 positive clinical samples chosen in the Cq range of 21 to 36 Cq. Incubation of samples for 20 minutes showed no effect on SARS-CoV-2 detection, seen as a mean change in Cq value of (0.01 +/- 0.5), however 60 and 90 minute incubation showed an increase in Cq value of up to 2 Cq (Fig 2c and Supplementary Fig. S6). This increase in Cq is above the intra assay variability of 0.33Cq and therefore indicates some loss of sensitivity for these longer incubation times (Supplementary Table 3). Similar to the observations on temperature range, the impact appears to be sample-dependent as opposed to Cq-dependent. For example, samples R1 and R2 having Cq values of 22 and 34 respectively, show only a 1Cq variation with time, while samples O1 and R4 having Cq values of 23 and 28 respectively, show over a 2 Cq variation with time (Fig. 2d and Supplementary Table 2).
In summary, higher temperatures and longer incubation times compared to the existing HI protocol of 10 min at 65°C result in higher Cq values in the RT-qPCR assay. The degree to which the Cq value is shifted appears to be sample-specific rather than Cq-specific. The most likely explanation being variation in OP/NP sample quality such as the presence of RNAses which may be activated during HI, leading to degradation of SARS-CoV-2 RNA [12]. Importantly, the detection of SARS-CoV-2 in all positive samples is maintained.
Validation of a bulk HI protocol for use in an approved clinical SARS-CoV-2 diagnostic assay.
In establishing a bulk HI protocol for viral inactivation, we found that some sample-dependent and unpredictable small loss in assay sensitivity was inevitable. In order to understand if this could be detrimental to the clinical test, we performed a large-scale concordance study in which a representative set of clinical samples were tested in both the standard procedure and the potential new bulk HI workflow in parallel. Three independent experiments were carried out in which clinical OP/NP swab samples were retained following the first step of clinical testing (sample transfer from swab vial), and vials were then re-capped and re-packaged in fresh UN3733 containment corresponding exactly to the original packing (specifically comprising a UN3733 leak proof bag, an inner sealable bag and an absorbent pad). Re-packed samples were subject to the bulk HI protocol (Supplementary Table 1) followed by RNA extraction and RT-qPCR following the CCTC standard procedure with the exception of omitting the laboratory-based HI step following sample lysis. Three independent experiments were run, each using 200-250 randomly chosen clinical OP/ON swab samples. Results were assessed in terms of test result concordance, by comparing the approved diagnostic assay result to the experimental bulk HI protocol test result for each sample (Fig. 3a, Supplementary Tables 4 and 5). The total sample set over the three experiments comprised 630 samples, 105 of which were given a SARS-CoV-2 positive result in the approved diagnostic assay, with RT-qPCR Cq values ranging from 17 to 38. The remaining test results for this sample set were negative (523 samples) or void (i.e. result not valid, 2 samples). Overall 94% of samples achieved the same test result using the bulk HI protocol, which included all but one positive sample with a Cq below 36, and 514 of the Negative samples (Fig. 3a). Positive samples with Cq values of 36 and above were variably detected between the two processes at approximately equal proportions, with 6 samples having concordance results as ‘Positive to Negative’ and 9 samples having results of ‘Negative to Positive’ (Supplementary Table 5). This variation reflects the limit of RT-qPCR assay sensitivity itself, where detection of such very weak positive samples is inherently inconsistent. When weak positive data from Cq 36 and above are excluded from the analysis, an overall concordance of 99% is seen. For samples that were identified as positive in both protocols, the bulk HI protocol Cq level was typically within 2 Cq of the approved assay result (Fig. 3b) and the average change of 0.81 +/- 1.08 Cq was very comparable to the value of 0.8 +/- 0.7 Cq seen for samples incubated for 60 minutes at 65°C (Supplementary Table 4 and Table 3 respectively), as indicated by earlier studies, some samples were seen in which the impact on Cq was greater, with 14 samples out of the 105 positive samples showing a deviation above 2Cq (Fig. 3b). As seen previously the Cq difference is typically independent of the actual Cq value (Fig. 3c and Supplementary Figure 7), however since reporting of SARS-CoV-2 test results is limited to positive, negative or void, the actual Cq value obtained in the assay has no impact on the reporting of positive samples as is the purpose of the assay. Overall the concordance data using clinical samples demonstrates that the bulk HI protocol is highly comparable to the existing approved assay, detects positive results with the same degree of reproducibility as the existing assay, and is therefore suitable for use as clinical diagnostic test.
In the concordance study, OP/NP swab samples were tested as soon as possible after the HI step of the original clinical test. In reality there was an unavoidable lag time of approximately 1 to 4 hours between the bulk HI step and the OP/NP sample being added to RNA extraction lysis buffer (at which point RNA integrity is preserved). This was due to the logistics of the experiment, i.e. the time it took to repack the sample swab vials and follow the bulk HI protocol. In the real-world environment of a high-throughput COVID-19 Testing Centre, the normal scenario is that samples are processed as quickly as possible and immediately after heat inactivation. However, unforeseen circumstances around general logistics or other incidents may arise and extend this time. Therefore, it was critical to determine the length of time following HI that samples can be processed and still preserve the quality of the diagnostic test. Similarly, as a contingency plan against oven malfunction or operator error that could result in uncertainty on whether bulk HI was successful, a fall-back process was needed. In our case, this would be to revert to the standard assay in which samples are heated in the lab after addition of RNA extraction lysis buffer.
To address the risk of time delay, the concordance study already described above had an additional test group in which the samples being retested were retained and then tested for a second time but after a 10 hour delay after the HI cycle. Concordance for this data set was comparable to the samples processed without a delay, i.e. in which the samples were tested within 1 to 4 hours following HI (Fig. 3b, Supplementary Table 5 and Supplementary Fig. S8). In addition, a separate study of 184 clinical OP/ON swab samples was performed in which sample processing occurred 24 hours after HI. For these samples we saw the same result in that concordance is preserved for SARS-CoV-2 positive samples below a Cq of 36 (Supplementary Fig. S9, Supplementary Table 5 and Fig. 3b).
To address the risk of a failed HI cycle and the resultant requirement to revert to the standard assay protocol, we carried out further experiments. Here we subjected clinical OP/NP swab samples to the bulk HI protocol, but then followed the standard assay as if no HI cycle had been used. This meant that the samples were heated twice during the workflow (Supplementary Fig. S10). The time taken to do this work meant that samples were retested using the experimental conditions one day after the clinical test. Nonetheless, we saw 88.8% concordance in the test results of the 171 samples processed, and as before, the majority of discordant test results were above a Cq value of 36 (Supplementary Fig. S10).
Finally, at the time this work was carried out, the CCTC was contributing to the genomic surveillance effort to monitor SARS-CoV-2 variants of concern by providing extracted RNA from positive tests to the Wellcome Sanger Institute for sequencing. Sequencing quality was not changed by the bulk HI workflow, as illustrated by an equivalent base coverage prior to and following implementation of the bulk HI protocol at the CCTC (Supplementary Fig. S11).