Sample pooling strategy was intended to use standard molecular diagnostic laboratory equipment and commercially available RNA isolation kits and SARS Cov2 PCR tests. The pooling experiments were performed according to the instructions for particular tests from manufacturers. The underlying assumption was to pool multiple spesimens and process them with only one RNA isolation event and a single PCR test without significant loss of test sensitivity. The presented approach leads to effective cost drop up to six-fold per sample in comparison to individual patient testing. To test the effectiveness and sensitivity of that strategy, we performed two separate sets of experiments. In the first, one SARS CoV-2 positive sample was pooled with eight negative. In the second, one positive sample was merged with five SARS CoV-2 negative. The basic idea was to combine and concentrate several samples to the maximal volume, which can be extracted with the single QIAamp Mini column. Obtained Ct values for pooled samples were compared with Ct values for SARS CoV-2 positive samples processed individually. The difference between these Ct's reflects the possible loss of initial viral RNA templates number in merged samples.
Conducted experiments showed the limits of sensitivity of the proposed strategy of pooling samples and Z-Path-COVID-19-CE Genesig Real-Time RT-PCR test itself. The choice of that particular test was motivated by the extensive documentation prepared by the manufacturer. The kit's instruction contains such information as analytical sensitivity and specificity as well as repeatability, inter-instrument reproducibility, operator reproducibility, daily reproducibility and accuracy. The declared effectiveness of SARS CoV-2 RNA detection with a 100% replicate detection at the level of about five virus copies per reaction was also confirmed in our research. Seven tenfold dilutions of standardised SARS-CoV-2 RNA templates starting from 1,336x106 to 1,34 copies per reaction were amplified. The calculated standard curve showed excellent linearity r2=0,991, Figure 2.
All the swabs used in the project (both COVID 19 positive and negative) were processed (extracted and tested) individually to confirm their status. Also, all supernatants left after pooled samples concentrating were subject to further testing for SARS CoV-2 RNA presence, and no signal was detected.
To calculate the most efficient pool size a web-based application for pooling as described at https://www.chrisbilder.com/shiny was used. The hierarchical testing included an experimental prevalence rate of 5%, an assay lower limit of detection of 1 to 3 RNA copies/µL, an assay sensitivity of 97%, an assay specificity of 100%, a two-stage pooling algorithm, and a range of pool sizes of 3 to 20 samples. These calculations predicted the optimal testing configuration of an initial pool size of 5 followed by individual testing. The expected number of tests for the optimal testing configuration is 2.10. This leads to an expected number of tests per individual of 2.10 / 5 = 0.42. Thus, two-stage hierarchical testing reduces the expected number of tests by 58% when compared to individual testing. For the overall implementation of the algorithm, the sensitivity is 0.9409, the specificity is 1.0000, the positive predictive value is 1.0000, and the negative predictive value was 0.9969.
Six and nine of pooled samples were tested in two different experiments in eight different reaction sets. Results were summarised in Table 1 and Table 2.
Table 1. Six samples pooling tests results.
Sample name
|
Average Ct
|
∆Ct1
|
Calculated starting templates quantity
|
pooled group #1
|
30.34 ± 0.42
|
1.57
|
169.82 ± 42.64
|
CoV-2 positive #1
|
28.78 ± 0.06
|
|
445.73 ±16.31
|
pooled group # 2
|
26.08±0.26
|
1.41
|
2459.16 ± 407.82
|
CoV-2 positive #2
|
24.68±0.06
|
|
5921.69±219.45
|
pooled group # 3
|
26.08±0.65
|
1.25
|
2562.31±898.64
|
CoV-2 positive #3
|
24.83±0.09
|
|
5383.47±295.71
|
pooled group # 4
|
31.36±0.14
|
0.50
|
87.74±7.68
|
CoV-2 positive #4
|
30.85±0.14
|
|
120.44±10.41
|
pooled group # 5
|
35.14±0.31
|
1.38
|
8.14±1.56
|
CoV-2 positive #5
|
33.76±0.32
|
|
19.47±4.06
|
pooled group # 6
|
-
|
-
|
-
|
CoV-2 positive #6
|
36.13±0.30
|
|
4.37±0.79
|
pooled group # 7
|
-
|
-
|
-
|
CoV-2 positive #7
|
36.06±0.16
|
|
4.52±0.47
|
pooled group # 8
|
31.58±0.45
|
1.86
|
78.09±22.88
|
CoV-2 positive #8
|
29.72±0.02
|
|
246.25±2.70
|
Average ∆Ct
|
1.33±0.46
|
|
1 The difference between Ct values of the particular pooled group and corresponding positive sample.
Two of eight positive samples (#6 and #7) showed the Ct value below the Genesig Real-Time RT-PCR test guaranteed limit of detection, which is 4.64 viral copies per reaction, Table 1, Table 2. However, both samples were analysed six times, and the repeatability of the test results was 100%. Positive status of these samples could not be confirmed with two other COVID 19 IVD multiplex tests (data not shown).
Table 2. Nine samples pooling test results.
Sample name
|
Average Ct
|
∆Ct1
|
Calculated starting templates quantity
|
pooled group #1
|
31.46±0.31
|
1.80
|
83.05±15.12
|
CoV-2 positive #1
|
29.66±0.16
|
|
256.22±24.97
|
pooled group # 2
|
27.34±0.24
|
2.62
|
1108.34±164.61
|
CoV-2 positive #2
|
24.72±0.01
|
|
5749.46±43.17
|
pooled group # 3
|
28.25±0.37
|
2.59
|
630.24±154.01
|
CoV-2 positive #3
|
25.67±0.05
|
|
3162.27±92.19
|
pooled group # 4
|
32.89±0.14
|
2.59
|
33.36±2.94
|
CoV-2 positive #4
|
30.29±0.18
|
|
171.58±19.55
|
pooled group # 5
|
37.30±1.72
|
3.01
|
2.69±2.45
|
CoV-2 positive #5
|
34.29±0.37
|
|
13.99±3.21
|
pooled group # 6
|
-
|
-
|
-
|
CoV-2 positive #6
|
36.80±0.06
|
|
2.83±0.11
|
pooled group # 7
|
-
|
-
|
-
|
CoV-2 positive #7
|
35.71±0.85
|
|
6.21±3.60
|
pooled group # 8
|
32.85±0.24
|
2.85
|
34.31±4.91
|
CoV-2 positive #8
|
30.00±0.14
|
|
205.49±17.99
|
|
Average ∆Ct
|
2.58±0.42
|
|
1 The difference between Ct values of the particular pooled group and corresponding positive sample.
Calculated average delta Ct value, showed a statistical difference between pooled and individual positive samples in examined sets. The average 1.33 cycles difference in six samples pooling, confirmed the 2.5 x loss of the initial amount of SARS CoV-2 RNA during the processing merged samples. The initial viral loss is even higher in nine samples pooling sets. The value of average delta Ct equal 2.58 can be translated to sixfold loss (5.97x) of the viral RNA through the pooling and concentrating process.
The analytical sensitivity of the Genesig Real-Time RT-PCR test is equal to 4.64 SARS-CoV-2 RNA copies in reaction. Basing on test's sensitivity and confirmed loss of RNA during the concentrating process, the minimal initial viral RNA templates number which could be detected is 11,6 viral copies in case of pooling of six samples, and 27.7 in nine samples merging. Many COVID 19 IVD multiplex real-time PCR test declare minimal analytical sensitivity on ten RNA copies in reaction. Pooling, concentrating and processing of six samples in single RNA extraction followed by SARS CoV-2 test, seemed to be an interesting alternative for lowering the already enormous COVID 19 molecular diagnostics costs.