Next Article in Journal
Cholesteryl Ester Transfer Protein Inhibition Reduces Major Adverse Cardiovascular Events by Lowering Apolipoprotein B Levels
Next Article in Special Issue
SARS-CoV-2 Spatiotemporal Genomic and Molecular Analysis of the First Wave of the COVID-19 Pandemic in Macaé, the Brazilian Capital of Oil
Previous Article in Journal
Multi-Omics Studies Unveil Extraciliary Functions of BBS10 and Show Metabolic Aberrations Underlying Renal Disease in Bardet–Biedl Syndrome
Previous Article in Special Issue
Role of the Pangolin in Origin of SARS-CoV-2: An Evolutionary Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Differentiation of SARS-CoV-2 Variants Using RT-qPCRs by Targeting Recurrent Mutation Sites: A Diagnostic Laboratory Experience from Multi-Center Regional Study, August 2020–December 2021, Poland

Department of Biomedical Research, National Medicines Institute, 00-725 Warsaw, Poland
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(16), 9416; https://doi.org/10.3390/ijms23169416
Submission received: 28 July 2022 / Revised: 16 August 2022 / Accepted: 18 August 2022 / Published: 20 August 2022
(This article belongs to the Special Issue Genomic Variation and Epidemiology of SARS-CoV-2)

Abstract

:
Rapid identification of SARS-CoV-2 variants is essential for epidemiological surveillance. RT-qPCR-based variant differentiation tests can be used to quickly screen large sets of samples for relevant variants of concern/interest; this study was conducted on specimens collected at 11 centers located in Poland during routine SARS-CoV-2 diagnostics between August 2020 and December 2021. A total of 1096 samples (with CT < 30) were screened for Alpha, Beta, Delta, Kappa and Omicron variants using commercial assays targeting repeat mutation sites. Variants were assigned to 434 (39.6%) specimens; the remaining 662 (60.4%) samples were not classified (no tested mutations detected). Alpha (n = 289; 66.59%), Delta (n = 115; 26.5%), Kappa (n = 30; 6.91%) and Omicron (n = 2; 0.46%) variants were identified and their distribution changed over time. The first Alpha variant appeared in October 2020, and it began to gradually increase its proportion of the virus population by June 2021. In July 2021, it was replaced by the Delta variant, which already dominated by the end of the year. The first Kappa was detected in October 2021, while Omicron was found in December 2021. The screening of samples allowed the determination of epidemiological trends over a time interval reflecting the national COVID-19 waves.

1. Introduction

The ongoing more than two years pandemic of coronavirus disease 2019 (COVID-19) caused by coronavirus 2 of the severe acute respiratory syndrome (SARS-CoV-2) constitutes a major global health concern [1]. Since the first identification of SARS-CoV-2 in December 2019 [2], the World Health Organization has recorded more than 560 million confirmed cases of COVID-19, including over 6 million deaths [3]; this pandemic is not only a serious public health problem but also an economic and social one [4]. In addition, COVID-19 remains a major epidemiological challenge due to the worldwide spread of SARS-CoV-2 and because the continued evolution of this virus indicates that it is likely to be around for many years to come [5,6].
Epidemiological investigations are increasingly using genome sequence data of disease-causing pathogens. For the first time, such a large amount of genome sequence data on an etiological agent of illness was generated and made publicly available just after the actual emergence of SARS-CoV-2 and its variants [5,6,7]. Under conditions of neutral genetic drift, the rate of SARS-CoV-2 evolution is estimated to be approximately 1 × 10−3 mutations/nucleotide yearly [8]. Although most mutations are insignificant, there are also some that confer selective benefits, e.g., related to changes in binding to host cell receptors, and this results in increased infectivity and transmissibility, as well as weak antigen-antibody interaction, reduced neutralization and finally, immune escape [9]; they are so advantageous that, through adaptive evolution within the host, they can increase their frequency in the overall viral population, leading to new lineages. The emergence of SARS-CoV-2 mutants with features dangerous to public health has attracted attention since late 2020; they have been classified as a variant of interest (VOI) and a variant of concern (VOC) in order to establish global surveillance and research efforts. These variants have specific multiple changes to the Wuhan-Hu strain isolated in China during the first outbreak in late 2019. The mutations may make them insensitive to vaccines and can lead to increased infectivity, risk of hospitalization and reduced treatment efficacy, as well as causing diagnostic failures [6,10]. These effects are mainly related to changes occurring in the SARS-CoV-2 spike (S) protein [11].
Among the most relevant variants, the highlighted ones are Alpha (lineage B.1.1.7), Beta (lineage B.1.351), Gamma (lineage P.1), Delta (lineage B.1.627.2) and Omicron (lineage B.1.1.529). All these variants belonged to the VOCs and spread worldwide [6,10,12]. The Alpha variant appeared for the first time in late September 2020 in the United Kingdom (UK) [13], and by the end of the year, it had become the dominant lineage in this country [14]; this variant presents three mutations of interest (N501Y, D614G and P681H) in the S protein [10]. The N501Y alteration increases the affinity of the virus for the host cell ACE2 receptor, and the D614G and the P681H mutations are responsible for enhanced viral transmission [15]. The Beta variant was first identified in October 2020 in South Africa and soon became the dominant strain, constituting the second wave in the region [16]; this variant contains five substitutions of interest in the S protein (K417N, E484K, N501Y, D614G and A701V), three of which are the same as in the Alpha [10]. The K417N and the E484K mutations increase the binding affinity to human ACE2, and the A701V affects antibody neutralization [15,17]. The Gamma variant was initially detected in Brazil [18] and Japan [19]; this variant is characterized by five mutations of interest in the S protein (K417T, E484K, N501Y, D614G H655Y), which, in addition to the last one, are also present in the Beta [10]. The H655Y change enhances the transmissibility of the virus [20]. The Delta variant was first found in India in December 2020 [21], and it includes four substitutions of interest in the S protein, three new ones, L452R and T478K and P681R, and the D614G earlier observed in Alpha, Beta and Gamma [10]. These mutations reduce the ability to neutralize antibodies to recognize the virus and block infection [22]. The Omicron variant was identified for the first time in November 2021 in South Africa [23]. There are several dozen mutations of interest in the S protein in the Omicron sub-lineages. Some of them have been previously identified in the other VOCs [10,23]. The Omicron has been documented in epidemiological studies to be associated with an increased risk of re-infection [24].
The identification of SARS-CoV-2 variants emphasized the need to develop widely available methods to detect and track mutations. Next-generation sequencing (NGS) has been recognized as the gold standard for screening SARS-CoV-2 variants [25]. Although NGS can characterize a variant with the highest accuracy, its implementation is time and labor-consuming, relatively expensive and difficult to access in most medical laboratories. In contrast, commercial PCR-based assays to identify mutations for variant differentiation are rapid, simple to perform, not too costly and therefore, can be used in routine diagnostics [26].
The aim of the study was a retrospective analysis of SARS-CoV-2 positive samples identified at the National Institute of Medicine (NMI) to differentiate virus variants using RT-qPCR kits targeting repeat mutation sites. The study focused on the determination of (i) the time of first appearance of the variant, (ii) the prevalence of variants, and finally, (iii) the changes in SARS-CoV-2 epidemiology over almost 1.5 years.

2. Results

2.1. Study Specimen Characteristics

2.1.1. Background Analysis

A total of 33,780 clinical samples from hospitals (H), mobile collection sites (MCS) and epidemiological surveillance of individuals in home isolation/quarantine (I/Q) were obtained between 1 August 2020 and 31 December 2021; this corresponded to 17,025 (50.4%) and 16,755 (49.6%) specimens during the study period of the year, respectively. The characteristics of the sample collection are shown in Table 1.
A total of 3217 (9.52%) SARS-CoV-2 positive cases were recognized; this included 2444 (14.36%) and 773 (4.61%) specimens recovered in 2020 and 2021, respectively. The most of them noticed in November 2020 (n = 1029; 19.17%), October 2020 (n = 758; 14.57%) and March 2021 (n = 246; 13.28%). The prevalence of SARS-CoV-2 positive cases from the particular centers varied with n = 99 (3.17%) for H1, n = 98 (21.97%) for H2, n = 698 (12.31%) for H3, n = 22 (7.12%) for H14, n = 164 (22.34%) for H5, n = 8 (9.2%) for H6, n = 62 (31.96%) for H7, n = 746 (52.57%) for MCS1, n = 223 (6.81%) for MCS2, n = 303 (17.63%) for I/Q1 and n = 21 (44.68%) for I/Q2 during the test period of 2020. Next year, it was n = 90 (1.18%), n = 54 (4.1%), n = 81 (9.57%), n = 9 (1.76%), n = 86 (10.94%), n = 0 (0%), n = 91 (20.04%), n = 216 (32.1%), n = 121 (2.95%), n = 8 (1.96%) and n = 17 (34%) from the center/collection, respectively.

2.1.2. Final Collection

The final collection (specimens with all tested viral gene cycle thresholds below thirty, CT < 30 as recommended by the test manufacturer) contained 1096 samples in total (34% of all positive cases in the range of CT = 40) (Table 1); it included 589 and 507 SARS-CoV-2 positive cases in 2020 and 2021, which corresponded to 24.1% and 65.59% of all positive specimens during the study periods, respectively. In 2020, this accounted for n = 27 (27.27%), n = 27 (27.55%), n = 147 (21.06%), n = 9 (40.91%), n = 40 (24.39%), n = 3 (37.5%), n = 16 (25.81%), n = 169 (22.65%), n = 30 (13.45%), n = 111 (36.63%) and n = 10 (47.62%) from the center/collection, respectively. The following year, it was n = 65 (72.22%), n = 12 (22.22%), n = 62 (76.54%), n = 8 (88.89%), n = 46 (53.49%), n = 0 (0%), n = 47 (51.65%), n = 147 (68.06%), n = 104 (85.95%), n = 2 (25%) and n = 14 (82.35%) in each facility, respectively.

2.1.3. Demographics of the Individuals

The characteristics of the final study cohort are shown in Table 2. The patient group consisted of 554 females (50.55%) and 542 males (49.45%) ranging in age from 0 to 94 years, with a median age of 48 years.
The most common group were individuals aged 36–64 with 502 (45.8%) (n = 270 and n = 232 in 2020 and 2021, respectively), followed by those aged 19–35 with 256 (23.36%) (n = 153 and n = 103 in 2020 and 2021, respectively), ≥65 with 254 (23.18%) (n = 124 and n = 130 in 2020 and 2021, respectively), 6–18 with 52 (4.75%) (n = 28 and n = 24 in 2020 and 2021, respectively) and 0–5 with 32 (2.92%) cases (n = 14 and n = 18 in 2020 and 2021, respectively).

2.2. Identification of SARS-CoV-2 Variants

A set of 1096 samples were tested for variant differentiation. Results of this analysis are shown in Table 3.
The study was performed for n = 70, n = 46, n = 130, n = 155, n = 188, n = 87, n = 55, n = 157, n = 40, n = 11, n = 1, n = 8, n = 2, n = 4, n = 25, n = 62 and n = 54 samples from each month over almost 1.5 years. The monthly distribution of SARS-CoV-2 variants is shown in Figure 1 and Table 4.
In total, 436 (39.78%) specimens were assigned a variant. The remaining 660 (60.22%) samples were not classified (no tested mutations detected). Alpha (n = 289; 66.59%), Delta (n = 115; 26.5%), Kappa (n = 30; 6.91%) and Omicron (n = 2; 0.46%) variants were detected in the tested specimens; the Beta variant was not identified. There was no statistical significant relationship between gender of patients in whom variants were detected (female, n = 554/1096 vs. male, n = 542/1096; p = 0.767).
Only mutations associated with the Alpha variant were detected in 69 (11.71%) specimens from 2020. The remaining 520 (88.26%) samples were not categorized. The first case of the Alpha variant was identified in a sample collected in 21 October 2020 from a patient hospitalized in H3.
Most cases of the Alpha variant were recorded in MCS1 (n = 30; 17.75%), but the highest percentage was in H2 (n = 11; 40.74%). In both hospital and outpatient samples from 2020, the amount of this variant was similar (p = 0.569) with n = 34/269 (12.64%) and n = 35/320 (10.94%) specimens, respectively.
A variety of variants were detected in specimens collected in 2021. Specific mutations were detected for 365 (72%) samples. The remaining 142 (28%) specimens were not classified. The Alpha variant was dominant and represented 220 (43.4%) probes. Between March and June 2021, the variant accounted for 92.36% (n = 145), 77.5% (n = 31), 75.% (n = 9) and 100% (n = 1) of the tested specimens in this time interval, respectively. Most cases of the Alpha variant were reported in MCS1 with n = 100 and MCS2 with n = 35 (68.03% and 33.65% of all samples from these centers, respectively) and the majority of them were collected in March (n = 88; 88% and n = 24; 72.73%, respectively); however, in H2 hospital the variant was present in 83.33% (n = 10) of the studied probes from this center.
The Alpha variant was replaced by the Delta variant in July 2021. A total of 115 (22.7%) samples represented the Delta variant. The first case of this variant was detected in a specimen recovered in 9 July 2021 from patient hospitalized in H3; this variant appeared in samples from several hospitals: H1(n = 34; 52.31%), H2 (n = 2; 16.67%), H3 (n = 22; 35.48%) and H5 (n = 9; 19.57%), mobile collection site MCS2 (n = 44; 42.31%), and epidemiological surveillance center I/Q2 (n = 4; 28.57%). The percentage of Delta variant in relation to other variants present in particular center was the highest in H1 (52.31%), H3 (35.48%) and MCS2 (42.31%). 27.92% (n = 67/240) of all hospital cases included in the study from 2021 were identified as Delta variant, while the non-hospital samples were only 17.98% (n = 48/267) (p = 0.034).
The first case of the Kappa variant was identified in sample taken 7 October 2021 from patient hospitalized in H1. A total of 30 (5.92%) samples represented this variant. The variant occurred at a frequency of 12%, 35.48% and 9.26% from October to December 2021, respectively; it appeared in few hospitals: H1 (n = 5; 7.69%), H3 (n = 10; 16.13%) and H5 (n = 1; 2.17%), and outpatients from MCS2 (n = 11; 10.58%) and I/Q2 (n = 3; 21.43%). Mutation specific to the Kappa variant accounted for 5.92% of all variant-differentiated samples from 2021.
Only two cases of the Omicron variant were reported in the study, both in December 2021 in MCS2; they consisted 1.92% of the mutation-tested samples in this center and 0.4% of all from 2021. The first one was identified in a specimen collected in 13 December 2021.

3. Discussion

The vaccination is one of the most effective ways to protect from infectious diseases. The development of vaccines against COVID-19 proved successful in preventing symptomatic infection and illness; however, despite the availability of vaccines, many countries experienced multiple COVID-19 waves, largely caused and driven by emerging new VOCs [27]; they prolong the continuation of the disease, leading to patient deaths, economic losses and a further burden on public health systems [4]. Globally, five COVID-19 waves were observed during the study period, from 1 August 2020 to 31 December 2021 [3]. The approximate peak disease periods are July–August 2020 (first wave), October 2020–January 2021 (second wave), March–April 2021 (third wave), July–August 2021 (fourth wave) and finally November–December 2021 (fifth wave with a peak in January 2022). While the first wave was due to the spread of a wild type (WT) of SARS-CoV-2, the following ones were caused by the emergence and dissemination of further different VOCs.
The situation was slightly different in Poland, where three waves of COVID-19 were observed at the time [28]. The first illness peak appeared in October–November 2020 (first wave), followed by another in February–March 2021 (second wave), and then October–December 2021 (third wave). The number of positive SARS-CoV-2 cases identified during the routine diagnostics in the NMI reflected a nationwide trend, with peaks in the first and second waves, as we already noted in our previously report [29]. A renewed increase in the number of positive samples occurred in the autumn/winter period (full data not shown, partial data in Figure 1), coinciding with the time of the third wave of COVID-19 in the country. A total of 3217 SARS-CoV-2 positive cases were detected, of which 1096 (34%) samples were included in the retrospective study. The final collection was represented by 589 (24.1%) and 507 (65.59%) specimens from 2020 and 2021, respectively. These data indicate that fewer patients in the first year of the study were diagnosed in the active phase of infection with high viral load (CT < 30) than in the second one; this result seems to be related to the health care burden, delays in diagnosis, and the availability of tests for SARS-CoV-2; however, it may also be due to the fact that SARS-CoV-2 was a new virus in the human population in 2020, so people attributed the initial symptoms of infection to other diseases, and therefore got examined later. A statistically significant majority of the tested samples were from adults (p < 0.0001) in roughly the same number of each gender, which is consistent with COVID-19 population studies [30,31]. Some papers report that males present an increased risk of severe acute respiratory syndrome infection and higher mortality [32,33]. In our study, we could not confirm this information because the laboratory did not possess clinical data to follow the diseases and their outcomes.
There are only few publications tracing SARS-CoV-2 variants in Poland [29,34,35,36,37]. Most of them based on analysis of viral genomic data deposited at the Global Initiative on Sharing All Influenza Data (GISAID) [5]; it currently contains 83,764 sequences from Poland [5]. Although the NGS is the gold standard for differentiating SARS-CoV-2 variants, due to its long turnaround time and high cost per sample, its use is less feasible by diagnostic laboratories [26]. Therefore, RT-qPCR assays detecting VOCs-specific mutations appear to be a better option for timely clinical applications and population surveillance [26,38,39,40,41]. The present study applied just this method to test a set of more than 1000 samples for SARS-CoV-2 variants. Most of the specimens from 2020 were not classified and likely represented the WT virus. These results are consistent with national data, where the original virus is assigned to the first wave of COVID-19 in the country [42]. The first VOC, Alpha variant, was identified in a sample collected in October 2020. Interestingly, previous data reported that this variant appeared in Poland in December 2020 [37,42]. Since its emergence, the Alpha variant has gradually increased its proportion in the following months, as we already observed in our previous survey [29]. Between March and June 2021, the variant accounted for more than 92% of all recognized cases; this finding also reflects national data, in which the Alpha variant was responsible for the second wave of the epidemic [42]. In July, it was replaced by the Delta variant (the first case identified at the beginning of the month) and was the most commonly identified variant until the end of 2021 (although, it should be noted that the May–September period was represented by a small number of samples). Country data indicate that its initial appearance was in April 2021, but similar to our study, it began to dominate from July and caused the third wave of COVID-19 [42]. The Delta variant was slightly more frequent in hospitalized patients than in outpatients, which was also observed by other authors [43,44]. Two cases of the Omicron variant were reported in December 2021. Based on national data, this variant appeared similar to our research in the last month of 2021 and triggered the fourth wave of disease the following year [42]. Of the investigated VOCs in this study, only the Beta variant was not detected, although other reports indicate few cases of it in the country in 2021 [35,42]; however, one assay identified the Kappa variant, which belongs to the VOIs. The first case of this variant was found in October 2021 and it was represented in the sample collection until the end of the year.
Our study shows that screening by RT-qPCR-based variant differentiation tests allows the identification of epidemiological trends (both temporal and variant dominance) reflecting national COVID-19 waves. Although some of the specimens were not classified by this method, it enabled the testing of a large set of samples. Based on the global data of SARS-CoV-2 spread, the non-categorized specimens from 2020 probably belonged to the WT variant; however, for those from 2021, it seems appropriate to call them “other” considering that they may possess some mutations to the Wuhan-Hu strain that were not identified by the applied assays. Potentially, we should also use the term “likely” for the detected variants (e.g., Alpha-like variant), since only specific mutations were found, and their presence and/or absence of other changes was not verified by genome sequencing. Although, as the NGS-based national data indicated, the variant matches were complete.
A limitation of the study was that not all samples were examined, but only those for which the CT values of the analyzed genes in RT-qPCRs were less than 30. Therefore, in the remaining cases, a specific variants were not defined, potentially underestimating their numbers; this is evidenced by our earlier study, that included all positive samples from the mobile collection site from this study (MCS1) [29]; however, it should be noted that the prevalence of the Alpha variant in the overall virus population remained similar in both investigations.

4. Materials and Methods

4.1. Specimen Collection

SARS-CoV-2-positive samples collected from patients between 1 August 2020 and 31 December 2021 were analyzed. These included inpatients from seven medical centers (4 third- and 3 s-level hospitals; H1–H7), and outpatients from two drive/walk-thru mobile collection sites (MCS1—drive-thru mobile collection site, MCS2—walk-thru mobile collection site) and epidemiological surveillance of individuals in home isolation/quarantine (I/Q1—samples from towns in Mazovia, I/Q2—samples from Warsaw) in Mazovia, Poland. The specimens were nasopharyngeal or combined pharyngeal and nasal swabs obtained by the NMI as part of routine examinations for SARS-CoV-2 infection. The samples were transported in virus-dedicated media, along with the patient’s basic demographics and consent to use the material for scientific research obtained through a laboratory form. All specimens were coded to prevent identification of personal information during the testing, and after the diagnostic process they were stored at −80 °C, pending a retrospective study.

4.2. Nucleic Acids Extraction Methods and SARS-CoV-2 Detection

Virus RNA acids were isolated using two procedures, column-based extraction by Viral DNA/RNA kit (A&A Biotechnology, Gdansk, Poland) and magnetic-based method via NucleoMag Pathogen kit (Machery-Nagel, Duren, Germany) with modifications to the manufacturers’ instructions, as described previously [45].
The MutaPLEX® Coronavirus (SARS-CoV-2) kit (Immundiagnostik AG, Bensheim, Germany) detecting S/RdRP and E genes of SARS-CoV-2 was used for in vitro diagnostics. RT-qPCRs were carried out with Applied Biosystems QuantStudio 6 Pro Real-Time PCR system (Life Technologies Holdings Pte Ltd., Singapore). Both the test parameters and reading rules (based on CT values) were performed according to the test manufacturer’s recommendations.

4.3. Differentiation of SARS-CoV-2 Variants

The samples with CT < 30 in all SARS-CoV-2 genes detected during routine diagnostics were tested for some VOCs identification using three commercial variant-differentiation assays. The tests were carried out according to the manufacturers’ guidelines, using an Applied Biosystems QuantStudio 6 Pro Real-Time PCR System instrument. To determine the time of first emergence of the variant in the study collection, all samples found at least one month before the initial worldwide appearance of the VOC were examined.

4.3.1. Alpha Variant Identification

Samples previously identified as SARS-CoV-2 positive isolated from patients between August 2020 and December 2021 were tested for Alpha using the Bosphore® SARS-CoV-2 Variant Detection Kit v1 assay (Anatolia Geneworks, Istanbul, Turkey); this assay uses appropriately designed primers and detects the A570D, P681H and Y144del mutations that have been observed in the UK variants.

4.3.2. Beta, Delta and Kappa Variants Identification

The specimens recovered from patients between September 2020 and December 2021 were tested for detection of the Beta, Delta and Kappa variants. The ID™ SARS-CoV-2/VOC evolution Pentaplex kit (ID Solutions, Grabels, France) was used to find these variants. The primers used in the assay serve to amplify specific target sequences. In addition to the variant-specific L452R, E484K and E484Q mutations, the kit also identifies unmodified N2, RdRp1 and RdRp2 SARS-CoV-2 genes, which allowed verification of the used samples. The manufacturer’s data show that the limit of detection is 10 copies per PCR reaction for SARS-CoV-2 and K417N mutation and 40 copies/PCR for L452R and E484K mutations.

4.3.3. Omicron Variant Identification

The samples collected from October 2021 to December 2021 were checked for Omicron variant. Detection of this variant was performed using PKamp VariantDetect SARS-CoV-2 RT-PCR KIT Combination 6 (PerkinElmer, Waltham, MA, USA); this kit detects L452R, P681H and P681R mutations using primers and TaqMan probes, which allows the identification of the Alpha and Delta variants in addition to Omicron. Furthermore, it also provides verification of SARS-CoV-2 positive samples by detecting the N/ORF1ab virus genes. The detection limit of this test amounts to 200 copies per PCR reaction for tested mutation.

4.4. Statistical Analysis

The statistical analysis of data and figure was performed with GraphPad Prism 7.0 software (GraphPad Software, Inc., San Diego, CA, USA). Two-sided chi-square or Fisher’s exact test was used in p-value calculation at alpha < 0.05.

5. Conclusions

The present study showed that differentiation of SARS-CoV-2 variants by RT-qPCR is highly effective in retrospective studies of a large set of samples and reflects national epidemiological trends (prevalence and types of virus variants). The obtained results entitle this method to application in real-time VOCs monitoring. Changes in the functional properties of new SARS-CoV-2 variants can lead to differences in the clinical course of the infection. Therefore, rapid identification of such variants in population allows timely introduction of appropriately adjusted prevention and adaptation of medical intervention in the form of an optimal procedure and treatment method.

Author Contributions

Conceptualization, M.K. and A.B.; methodology, K.W., M.K., J.W., A.K., I.W. and A.B.; statistical analysis, J.W.; tables and figure preparation, K.W., J.W. and A.B.; writing—original draft preparation, K.W. and A.B.; writing—review and editing, K.W. and A.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financed by internal funding from the National Medicines Institute.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and was a part of routine diagnostics in a medical laboratory. The samples used for scientific research were analyzed anonymously. The genetic material of the viruses was the object of interest and human material was not analyzed.

Informed Consent Statement

Informed consent including General Data Protection Regulation (GDPR) was obtained from all subjects involved in the study.

Data Availability Statement

All relevant data are within the manuscript.

Acknowledgments

We would like to thank the laboratory workers who performed routine diagnostics on the samples.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Coronavirus Disease 2019. Available online: https://www.who.int/emergencies/diseases/novelcoronavirus-2019 (accessed on 13 June 2022).
  2. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Xiang, Z.; Huang, B.; Shi, W.; Lu, R.; et al. A novel coronavirus from patients with pneumonia in China, 2019. N. Engl. J. Med. 2020, 382, 727–733. [Google Scholar] [CrossRef] [PubMed]
  3. World Health Organization. WHO Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int (accessed on 22 July 2022).
  4. Jackson, J.K. Global Economic Effects of COVID-19: Overview—CRS Reports. Available online: https://crsreports.congress.gov/product/pdf/R/R46270#:~:text=Estimates%20indicate%20the%20COVID%2D19,slower%20rate%20projected%20for%202022 (accessed on 13 June 2022).
  5. Global Initiative on Sharing All Influenza Data (GISAID). Munich: Global Initiative on Sharing All Influenza Data. c2008–2022. Available online: https://www.gisaid.org (accessed on 13 June 2022).
  6. World Health Organization. Tracking SARS-CoV-2 Variants as of 7 June 2022. Available online: https://www.who.int/activities/tracking-SARS-CoV-2-variants (accessed on 13 June 2022).
  7. Wu, F.; Zhao, S.; Yu, B.; Chen, Y.-M.; Wang, W.; Song, Z.-G.; Hu, Y.; Ta, Z.-W.; Tian, J.-H.; Pei, Y.-Y.; et al. A new coronavirus associated with human respiratory disease in China. Nature 2020, 579, 265–269. [Google Scholar] [CrossRef] [Green Version]
  8. Duchene, S.; Featherstone, L.; Haritopoulou-Sinanidou, M.; Rambaut, A.; Lemey, P.; Baele, G. Temporal signal and the phylodynamic threshold of SARS-CoV-2. Virus Evol. 2020, 6, veaa061. [Google Scholar] [CrossRef] [PubMed]
  9. Harvey, W.T.; Carabelli, A.M.; Jackson, B.; Gupta, R.K.; Thomson, E.C.; Harrison, E.M.; Ludden, C.; Reeve, R.; Rambaut, A.; COVID-19 Genomics UK (COG-UK) Consortium; et al. SARS-CoV-2 variants, spike mutations and immune escape. Nat. Rev. Microbiol. 2021, 19, 409–424. [Google Scholar] [CrossRef] [PubMed]
  10. European Centre for Disease Prevention and Control. SARS-CoV-2 Variants of Concern as of 9 June 2022. Available online: https://www.ecdc.europa.eu/en/covid-19/variants-concern (accessed on 13 June 2022).
  11. Di Caro, A.; Cunha, F.; Petrosillo, N.; Beeching, N.; Ergonul, O.; Petersen, E.; Koopmans, M. SARS-CoV-2 escape mutants and protective immunity from natural infections or immunizations. Clin. Microbiol. Infect. 2021, 27, 823–826. [Google Scholar] [CrossRef]
  12. Lazarevic, I.; Pravica, V.; Miljanovic, D.; Cupic, M. Immune Evasion of SARS-CoV-2 Emerging Variants: What Have We Learnt So Far? Viruses 2021, 13, 1192. [Google Scholar] [CrossRef]
  13. Rambaut, A.; Loman, N.; Pybus, O.; Barclay, W.; Barrett, J.; Carabelli, A.; Connor, T.; Peacock, T.; Robertson, D.L.; Volz, E.; et al. Preliminary Genomic Characterisation of an Emergent SARS-CoV-2 Lineage in the UK Defined by a Novel Set of Spike Mutations. Available online: https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-ukdefined-by-a-novel-set-of-spike-mutations/563 (accessed on 13 June 2022).
  14. Volz, E.; Mishra, S.; Chand, M.; Barrett, J.C.; Johnson, R.; Geidelberg, L.; Hinsley, W.R.; Laydon, D.J.; Dabrera, G.; O’Toole, A.; et al. Assessing Transmissibility of SARS-CoV-2 Lineage B.1.1.7 in England. Nature 2021, 593, 266–269. [Google Scholar] [CrossRef]
  15. Choi, J.Y.; Smith, D.M. SARS-CoV-2 Variants of Concern. Yonsei Med. J. 2021, 62, 961–968. [Google Scholar] [CrossRef]
  16. Tegally, H.; Wilkinson, E.; Giovanetti, M.; Iranzadeh, A.; Fonseca, V.; Giandhari, J.; Doolabh, D.; Pillay, S.; San, E.J.; Msomi, N.; et al. Detection of a SARS-CoV-2 variant of concern in South Africa. Nature 2021, 592, 438–443. [Google Scholar] [CrossRef]
  17. Mohammadi, M.; Shayestehpour, M.; Mirzaei, H. The impact of spike mutated variants of SARS-CoV2 [Alpha, Beta, Gamma, Delta, and Lambda] on the efficacy of subunit recombinant vaccines. Braz. J. Infect. Dis. 2021, 25, 101606. [Google Scholar] [CrossRef]
  18. Faria, N.R.; Morales Claro, I.; Candid, D.; Moyses Franco, L.A.; Andrade, P.S.; Coletti, T.M.; Silva, C.A.M.; Sales, F.C.; Manuli, E.R.; Aguiar, R.S.; et al. Genomic Characterisation of an Emergent SARS-CoV-2 Lineage in Manaus: Preliminary Findings. Available online: https://virological.org/t/genomiccharacterisation-of-an-emergent-sars-cov-2-lineagein-manaus-prelimina (accessed on 17 June 2022).
  19. Fujino, T.; Nomoto, H.; Kutsuna, S.; Ujiie, M.; Suzuki, T.; Sato, R.; Fujimoto, T.; Kuroda, K.; Wakita, T.; Ohmagari, N. Novel SARS-CoV-2 variant in travelers from Brazil to Japan. Emerg. Infect. Dis. 2021, 27, 1243–1245. [Google Scholar] [CrossRef] [PubMed]
  20. Escalera, A.; Gonzalez-Reiche, A.S.; Aslam, S.; Mena, I.; Laporte, M.; Pearl, R.L.; Fossati, A.; Rathnasinghe, R.; Alshammary, H.; van de Guchte, A.; et al. Mutations in SARS-CoV-2 variants of concern link to increased spike cleavage and virus transmission. Cell Host Microbe 2022, 30, 373–387.e7. [Google Scholar] [CrossRef] [PubMed]
  21. Joshi, M.; Kumar, M.; Srivastava, V.; Kumar, D.; Rathore, D.; Pandit, R.; Chaitanya, G. First Detection of SARS-CoV-2 Delta Variant (B.1.617.2) in the Wastewater of (Ahmedabad) India. Available online: https://www.medrxiv.org/content/10.1101/2021.07.07.21260142v1 (accessed on 17 June 2022).
  22. Liu, Z.; VanBlargan, L.A.; Bloyet, L.M.; Rothlauf, P.W.; Chen, R.E.; Stumpf, S.; Zhao, H.; Errico, J.M.; Theel, E.S.; Liebeskind, M.J.; et al. Identification of SARS-CoV-2 spike mutations that attenuate monoclonal and serum antibody neutralization. Cell Host Microbe 2021, 29, 477–488.e4. [Google Scholar] [CrossRef] [PubMed]
  23. Wang, L.; Cheng, G. Sequence analysis of the emerging SARS-CoV-2 variant Omicron in South Africa. J. Med. Virol. 2022, 94, 1728–1733. [Google Scholar] [CrossRef]
  24. Pulliam, J.R.C.; van Schalkwyk, C.; Govender, N.; von Gottberg, A.; Cohen, C.; Groome, M.J.; Dushoff, J.; Mlisana, K.; Moultrie, H. Increased risk of SARS-CoV-2 reinfection associated with emergence of Omicron in South Africa. Science 2022, 376, eabn4947. [Google Scholar] [CrossRef]
  25. World Health Organization. Genomic Sequencing of SARS-CoV-2: A Guide to Implementation for Maximum Impact on Public Health. Available online: https://www.who.int/publications/i/item/9789240018440 (accessed on 13 June 2022).
  26. European Centre for Disease Prevention and Control (ECDC). Methods for the Detection and Identification of SARS-CoV-2 Variants. Available online: https://www.ecdc.europa.eu/en/publications-data/methods-detection-and-identification-SARSCoV-2-variant (accessed on 13 June 2022).
  27. Dutta, A. COVID-19 waves: Variant dynamics and control. Sci. Rep. 2022, 12, 9332. [Google Scholar] [CrossRef]
  28. World Health Organization. WHO Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int/region/euro/country/pl (accessed on 22 July 2022).
  29. Kozinska, A.; Wegrzynska, K.; Komiazyk, M.; Walory, J.; Wasko, I.; Baraniak, A. Viral Etiological Agent(s) of Respiratory Tract Infections in Symptomatic Individuals during the Second Wave of COVID-19 Pandemic: A Single Drive-Thru Mobile Collection Site Study. Pathogens 2022, 11, 475. [Google Scholar] [CrossRef]
  30. Louca, S. SARS-CoV-2 infections in 165 countries over time. Int. J. Infect. Dis. 2021, 111, 336–346. [Google Scholar] [CrossRef]
  31. Martín, U.; Bacigalupe, A.; Jiménez Carrillo, M. COVID-19 and gender: Certainties and uncertainties in monitoring the pandemic. Rev. Esp. Salud Publica 2021, 95, 33927179. [Google Scholar]
  32. Haitao, T.; Vermunt, J.V.; Abeykoon, J.; Ghamrawi, R.; Gunaratne, M.; Jayachandran, M.; Narang, K.; Parashuram, S.; Suvakov, S.; Garovic, V.D. COVID-19 and Sex Differences: Mechanisms and Biomarkers. Mayo Clin. Proc. 2020, 95, 2189–2203. [Google Scholar] [CrossRef]
  33. Mukherjee, S.; Pahan, K. Is COVID-19 Gender-sensitive? J. Neuroimmune Pharmacol. 2021, 16, 38–47. [Google Scholar] [CrossRef] [PubMed]
  34. Komiazyk, M.; Walory, J.; Gawor, J.; Ksiazek, I.; Gromadka, R.; Baraniak, A. Case Report of COVID-19 after full vaccination: Viral loads and anti-SARS-CoV-2 antibodies. Diagnostics 2021, 11, 1815. [Google Scholar] [CrossRef] [PubMed]
  35. Charkiewicz, R.; Nikliński, J.; Biecek, P.; Kiśluk, J.; Pancewicz, S.; Moniuszko-Malinowska, A.M.; Flisiak, R.; Krętowski, A.J.; Dzięcioł, J.; Moniuszko, M.; et al. The first SARS-CoV-2 genetic variants of concern (VOC) in Poland: The concept of a comprehensive approach to monitoring and surveillance of emerging variants. Adv. Med. Sci. 2021, 66, 237–245. [Google Scholar] [CrossRef] [PubMed]
  36. Serwin, K.; Aksak-Wąs, B.; Parczewski, M. Phylodynamic Dispersal of SARS-CoV-2 Lineages Circulating across Polish-German Border Provinces. Viruses 2022, 14, 884. [Google Scholar] [CrossRef]
  37. Gogolewski, K.; Miasojedow, B.; Sadkowska-Todys, M.; Stepień, M.; Demkow, U.; Lech, A.; Szczurek, E.; Rabczenko, D.; Rosińsk, M.; Gambin, A. Data-driven case fatality rate estimation for the primary lineage of SARS-CoV-2 in Poland. Methods 2022, 203, 584–593. [Google Scholar] [CrossRef] [PubMed]
  38. Ham, R.E.; Smothers, A.R.; Che, R.; Sell, K.J.; Peng, C.A.; Dea, D. Identifying SARS-CoV-2 Variants of Concern through Saliva-Based RT-qPCR by Targeting Recurrent Mutation Sites. Microbiol. Spectr. 2022, 10, e0079722. [Google Scholar] [CrossRef]
  39. Umunnakwe, C.N.; Makatini, Z.N.; Maphanga, M.; Mdunyelwa, A.; Mlambo, K.M.; Manyaka, P.; Nijhuis, M.; Wensing, A.; Tempelman, H.A. Evaluation of a commercial SARS-CoV-2 multiplex PCR genotyping assay for variant identification in resource-scarce settings. PLoS ONE 2022, 17, e0269071. [Google Scholar] [CrossRef]
  40. Castro, G.M.; Sicilia, P.; Bolzon, M.L.; Lopez, L.; Barbás, M.G.; Pisano, M.B.; Ré, V.E. Tracking SARS-CoV-2 Variants Using a Rapid Typification Strategy: A Key Tool for Early Detection and Spread Investigation of Omicron in Argentina. Front. Med. 2022, 9, 851861. [Google Scholar] [CrossRef]
  41. Borillo, G.A.; Kagan, R.M.; Marlowe, E.M. Rapid and Accurate Identification of SARS-CoV-2 Variants Using Real Time PCR Assays. Front. Cell. Infect. Microbiol. 2022, 12, 894613. [Google Scholar] [CrossRef]
  42. Visualization of COVID-19 Spread. Available online: https://covidhub.psnc.pl/eng/mapa/#! (accessed on 23 July 2022).
  43. Sheikh, A.; McMenamin, J.; Taylor, B.; Robertson, C. SARS-CoV-2 Delta VOC in Scotland: Demographics, risk of hospital admission, and vaccine effectiveness. Lancet 2021, 397, 2461–2462. [Google Scholar] [CrossRef]
  44. Shiehzadegan, S.; Alaghemand, N.; Fox, M.; Venketaraman, V. Analysis of the Delta Variant B.1.617.2 COVID-19. Clin. Pract. 2021, 11, 778–784. [Google Scholar] [CrossRef] [PubMed]
  45. Komiazyk, M.; Walory, J.; Kozinska, A.; Wasko, I.; Baraniak, A. Impact of the Nucleic Acid Extraction Method and the RT-qPCR Assay on SARS-CoV-2 Detection in Low-Viral Samples. Diagnostics 2021, 11, 2247. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The timeline of distribution of SARS-CoV-2 variants. The light blue color represents specimens for which the tested mutations were not identified, red corresponds to Alpha, dark blue—Delta, green—Kappa, yellow—Omicron variant. The values in the columns indicate the number of samples.
Figure 1. The timeline of distribution of SARS-CoV-2 variants. The light blue color represents specimens for which the tested mutations were not identified, red corresponds to Alpha, dark blue—Delta, green—Kappa, yellow—Omicron variant. The values in the columns indicate the number of samples.
Ijms 23 09416 g001
Table 1. Centers included in the study with the numbers of total samples and negative/positive cases recovered at each site in 2020 and 2021.
Table 1. Centers included in the study with the numbers of total samples and negative/positive cases recovered at each site in 2020 and 2021.
Center 1Type of
Center/
Collection
20202021
No. of
Samples
No. of
Negative
Cases
No. of
Positive
Cases
No. of
Positive
Cases
with
CT < 30
No. of
Samples
No. of
Negative
Cases
No. of
Positive
Cases
No. of
Positive
Cases
with
CT < 30
H1third level
hospital
312530269927760374299065
H24463489827131712245412
H3567149736981478467478162
H430928722951049898
H5second level hospital734570164407866588646
H68779830000
H719413262164543449147
MCS1drive-thru
walk-thru
1419673746169673435216147
MCS2327430512233041073971121104
I/Q1Epidemio-logical Survei-llance1719141630311140939482
I/Q24726211050291714
1 H—hospital, MCS—mobile collection site, I/Q—home isolation/quarantine.
Table 2. Characteristics of the patients with samples subjected to SARS-CoV-2 variant differentiation by center, gender and age group.
Table 2. Characteristics of the patients with samples subjected to SARS-CoV-2 variant differentiation by center, gender and age group.
No. of Tested Samples 2
Center 1H1H2H3H4H5H6H7MCS1MCS2I/Q1I/Q2
Variables
Female13/3613/277/319/722/263/04/2380/7114/6055/22/4
Male14/2914/1070/310/118/200/012/2489/7616/4456/08/10
Age in years:
0–51/100/05/30/20/00/01/13/10/14/00/0
6–181/101/08/10/02/00/00/07/20/119/00/0
19–354/151/027/154/411/31/00/157/2710/3637/01/2
36–6418/2820/356/245/218/222/05/1980/7919/4842/25/5
≥653/25/951/190/09/210/010/2622/381/819/04/7
1 H—hospital, MCS—mobile collection site, I/Q—home isolation/quarantine. 2 the numbers of tested samples are separated for 2020 and 2021.
Table 3. SARS-CoV-2 variants in each center with the numbers of their samples in 2020 and 2021.
Table 3. SARS-CoV-2 variants in each center with the numbers of their samples in 2020 and 2021.
Center 120202021
No. of
Tested
Samlpes
SARS-CoV-2 VariantNo. of
Tested
Samples
SARS-CoV-2 Variant
AlphaOtherAlpha DeltaKappaOmicron Other
H127324651834508
H227111612102000
H3147111366222221008
H4909860002
H540634461391023
H6303000000
H716313471100036
MCS11693013914710000047
MCS230129104354411212
I/Q11114107200002
I/Q2100101454302
1 H—hospital, MCS—mobile collection site, I/Q—home isolation/quarantine.
Table 4. The monthly prevalence of SARS-CoV-2 variants in particular centers.
Table 4. The monthly prevalence of SARS-CoV-2 variants in particular centers.
Month
of Collection
VariantNo. of Variant Samples in Centers 1
H1H2H3H4H5H6H7MCS1MCS2I/Q1I/Q2
20 AugustOther0015000001540
20 SeptemberOther0017000005240
20 OctoberAlpha10501000040
Other923801030188292
20 NovemberAlpha113020018100
Other911375130644202
20 DecemberAlpha1103030312000
Other632941107771306
21 JanuaryAlpha33210025000
Other1050701434721
21 FebruaryAlpha00610027200
Other0021102111100
21 MarchAlpha8691502882402
Other10007012100
21 AprilAlpha40337030803
Other10107000000
21 MayAlpha31201010100
Other20001000000
21 JuneAlpha00000010000
21 JulyDelta10200000500
21 AugustDelta10100000000
21 SeptemberDelta10101000000
Other00010000000
21 OctoberDelta900010001100
Kappa10000000200
Other10000000000
21 NovemberDelta9112020001103
Kappa20900000803
Other00000000200
21 DecemberDelta1316050001701
Kappa20101000100
Omicron00000000200
Other20000000101
1 H—hospital, MCS—mobile collection site, I/Q—home isolation/quarantine.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wegrzynska, K.; Komiazyk, M.; Walory, J.; Kozinska, A.; Wasko, I.; Baraniak, A. Differentiation of SARS-CoV-2 Variants Using RT-qPCRs by Targeting Recurrent Mutation Sites: A Diagnostic Laboratory Experience from Multi-Center Regional Study, August 2020–December 2021, Poland. Int. J. Mol. Sci. 2022, 23, 9416. https://doi.org/10.3390/ijms23169416

AMA Style

Wegrzynska K, Komiazyk M, Walory J, Kozinska A, Wasko I, Baraniak A. Differentiation of SARS-CoV-2 Variants Using RT-qPCRs by Targeting Recurrent Mutation Sites: A Diagnostic Laboratory Experience from Multi-Center Regional Study, August 2020–December 2021, Poland. International Journal of Molecular Sciences. 2022; 23(16):9416. https://doi.org/10.3390/ijms23169416

Chicago/Turabian Style

Wegrzynska, Karolina, Magdalena Komiazyk, Jaroslaw Walory, Aleksandra Kozinska, Izabela Wasko, and Anna Baraniak. 2022. "Differentiation of SARS-CoV-2 Variants Using RT-qPCRs by Targeting Recurrent Mutation Sites: A Diagnostic Laboratory Experience from Multi-Center Regional Study, August 2020–December 2021, Poland" International Journal of Molecular Sciences 23, no. 16: 9416. https://doi.org/10.3390/ijms23169416

APA Style

Wegrzynska, K., Komiazyk, M., Walory, J., Kozinska, A., Wasko, I., & Baraniak, A. (2022). Differentiation of SARS-CoV-2 Variants Using RT-qPCRs by Targeting Recurrent Mutation Sites: A Diagnostic Laboratory Experience from Multi-Center Regional Study, August 2020–December 2021, Poland. International Journal of Molecular Sciences, 23(16), 9416. https://doi.org/10.3390/ijms23169416

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop