Next Article in Journal
A Nectin1 Mutant Mouse Model Is Resistant to Pseudorabies Virus Infection
Next Article in Special Issue
Detection of SARS-CoV-2 Nucleocapsid, Spike, and Neutralizing Antibodies in Vaccinated Japanese
Previous Article in Journal
Mass Spectrometry Approaches for SARS-CoV-2 Detection: Harnessing for Application in Food and Environmental Samples
Previous Article in Special Issue
Early and Longitudinal Humoral Response to the SARS-CoV-2 mRNA BNT162b2 Vaccine in Healthcare Workers: Significance of BMI, Adipose Tissue and Muscle Mass on Long-Lasting Post-Vaccinal Immunity
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 14
Issue 5
10.3390/v14050873
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

SARS-CoV-2 Variants in Paraguay: Detection and Surveillance with an Economical and Scalable Molecular Protocol

by
Magaly Martinez
1,†,
Phuong-Vi Nguyen
2,†,
Maxwell Su
2,
Fátima Cardozo
3,
Adriana Valenzuela
3,
Laura Franco
1,
María Eugenia Galeano
1,
Leticia Elizabeth Rojas
1,
Chyntia Carolina Díaz Acosta
1,
Jonás Fernández
1,
Joel Ortiz
1,
Florencia del Puerto
1,
Laura Mendoza
3,
Eva Nara
1,
Alejandra Rojas
4,*,‡ and
Jesse J. Waggoner
2,5,*,‡
1
Departamento de Biología Molecular y Biotecnología, Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción, San Lorenzo 111241, Paraguay
2
Department of Medicine, Division of Infectious Diseases, Emory University, 1760 Haygood Drive NE, Room E-169, Bay E-1, Atlanta, GA 30322, USA
3
Departamento de Salud Pública, Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción, San Lorenzo 111241, Paraguay
4
Departamento de Producción, Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción, Dr. Cecilio Báez y Dr. Villamayor, Campus Universitario, San Lorenzo 111241, Paraguay
5
Department of Global Health, Rollins School of Public Health, Atlanta, GA 30322, USA
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
‡
These authors contributed equally to this work.
Viruses 2022, 14(5), 873; https://doi.org/10.3390/v14050873
Submission received: 2 March 2022 / Revised: 18 April 2022 / Accepted: 19 April 2022 / Published: 22 April 2022
(This article belongs to the Collection SARS-CoV-2 and COVID-19)
Browse Figures
Versions Notes

Abstract

:
SARS-CoV-2 variant detection relies on resource-intensive whole-genome sequencing methods. We sought to develop a scalable protocol for variant detection and surveillance in Paraguay, pairing rRT-PCR for spike mutations with Nanopore sequencing. A total of 201 acute-phase nasopharyngeal samples were included. Samples were positive for the SARS-CoV-2 N2 target and tested with the Spike SNP assay to detect mutations associated with the following variants: alpha (501Y), beta/gamma (417variant/484K/501Y), delta (452R/478K), and lambda (452Q/490S). Spike SNP calls were confirmed using amplicon (Sanger) sequencing and whole-genome (Nanopore) sequencing on a subset of samples with confirmed variant lineages. Samples had a mean N2 Ct of 20.8 (SD 5.6); 198/201 samples (98.5%) tested positive in the Spike SNP assay. The most common genotype was 417variant/484K/501Y, detected in 102/198 samples (51.5%), which was consistent with the P.1 lineage (gamma variant) in Paraguay. No mutations (K417 only) were found in 64/198 (32.3%), and K417/484K was identified in 22/198 (11.1%), consistent with P.2 (zeta). Seven samples (3.5%) tested positive for 452R without 478K, and one sample with genotype K417/501Y was confirmed as B.1.1.7 (alpha). The results were confirmed using Sanger sequencing in 181/181 samples, and variant calls were consistent with Nanopore sequencing in 29/29 samples. The Spike SNP assay could improve population-level surveillance for mutations associated with SARS-CoV-2 variants and inform the judicious use of sequencing resources.

1. Introduction

Amidst the global pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), numerous variants have emerged due to mutations in the positive-sense RNA genome. Variants of concern (VOCs) and variants of interest (VOIs) bear mutations that impact detection, treatment, clinical severity, transmission, and/or immune protection from prior infection or vaccination [1]. Whole-genome sequencing is the reference method for detecting and tracking variants, but this is a time- and resource-intensive process with the capacity varying markedly between regions [2,3,4,5]. It was estimated that 5% of all samples need to be sequenced to reliably detect variants at a prevalence of 0.1–1% [3], yet for many regions, this goal remains unattainable. Limited available sequencing capacity and overwhelming case numbers have motivated the development of real-time reverse transcription PCRs (rRT-PCRs) for the detection and surveillance of specific mutations found among VOCs/VOIs [6,7,8,9,10,11,12,13,14,15,16,17,18,19,20]. For these reasons, our group designed the Spike SNP assay, which is a multiplex rRT-PCR designed to detect specific mutations in the receptor binding motif of the Spike receptor binding domain (RBD; Figure 1) [20].
In Paraguay, the first confirmed SARS-CoV-2 infection was reported on 7 March 2020 [21], but 18 months into the pandemic, <0.1% of samples from the country could be sequenced and entered in the GISAID database (www.gisaid.org, accessed 1 March 2022) [22]. Published reports on case series in Paraguay have been limited to descriptive analyses of symptoms and disease outcomes [23,24]. As such, little is known about variants that circulated in the country during the beginning of the pandemic in 2020. However, data from 2021 show that variants originally detected in Brazil, such as the P.1 (gamma variant) and P.2 (zeta) lineages, circulated widely until the middle of that year (GISAID database). Nanopore sequencing on the portable MinION sequencer from ONT (Oxford Nanopore Technologies, Oxford, UK) became available in Asunción, the capital of Paraguay, to detect and monitor SARS-CoV-2 variants for the first time, where this technology provides lower-cost and more-rapid whole-genome sequencing compared to other methods [2,25]. The objective of the current study was to develop an economical and scalable protocol based on the Spike SNP assay and Nanopore sequencing to identify and monitor VOCs/VOIs in acute-phase samples in Paraguay.

2. Methods

2.1. Clinical Samples

Nasopharyngeal (NP) samples were obtained from residual baseline surveillance samples at the Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción (IICS-UNA) Laboratory. Samples were de-identified, aliquoted, and stored at −80 °C until nucleic acid extraction. A convenience set of 201 SARS-CoV-2 positive samples from November 2020 to April 2021 was chosen for this study. Samples tested positive for the envelope gene target of the Charité protocol [26] with a cycle threshold (Ct) < 34 and were collected from different cities in Paraguay from symptomatic (n = 187) and asymptomatic (n = 14) individuals. This study was reviewed and approved by the IICS Scientific and Ethics Committee (P38/2020) and the Emory Institutional Review Board (study 00110736). Patient consent was waived under approval of IICS Scientific and Ethics Committee (P38/2020) due to all samples being de-identified remaining of diagnostics samples.

2.2. Spike SNP Design and Performance

Hydrolysis probes were designed to detect mutations associated with emerging variants of concern: L452Q (CAG > CTG) and F490S (UUU > UCU, C.37, lambda), E484Q (GAA > CAA, B.1.617.1, kappa), T478K (ACA > AAA, B.1.617.2, delta), and a modified T478K containing one mixed base in codon 477 for the detection of omicron. All new probes were designed as described [20], without locked nucleic acids (LNAs). Unmodified probes for 484K and 501Y were also designed and evaluated. New probes were individually evaluated using SARS-CoV-2 genomic RNA in singleplex assays. Probes that provided the most sensitive signal based on the cycle threshold (Ct) value, with preserved specificity, were selected. Subsequently, probes were tested in combination and compared side-by-side with the corresponding singleplex assays. Probe evaluations were performed by utilizing variant control samples that were kindly provided by the National Institutes of Health Variant Task Force and contained the following SARS-CoV-2 variants: alpha, beta, gamma, kappa, lambda, delta, and omicron. Variant control samples were heat-inactivated clinical samples that were prepared in the same manner as all clinical samples and then tested across a range of N2 Ct values from 20–36. Variants and lineages had been confirmed via sequencing (Helix, San Mateo, CA, USA).
Nucleic acids were extracted from all samples on an EMAG instrument (bioMérieux, Durham, NC) from 200 µL of NP swab and eluted in 50 µL. All rRT-PCRs were performed on a Rotor-Gene Q instrument (Qiagen, Germantown, MD, USA) using 20 µL reactions of the Luna Probe One-Step RT-qPCR Kit (New England Biolabs, Ipswich, MA, USA) and 5 µL of nucleic acid eluate. Each sample was tested in a duplex reaction for the SARS-CoV-2 nucleocapsid 2 target and RNase P (N2RP), as well as two different reactions of the Spike SNP assay, including probes for (1) K417, 452R, 484K, and 501Y, as described [20]; and (2) 452Q, 490S, and 484Q. During the study, two probes for 478K were designed, and all samples with 452R were tested in singleplex reactions with the original probe to search for the delta variant. Analytical performance of the second Spike SNP reaction was also evaluated following the inclusion of this probe. As sample collection pre-dated the emergence of omicron, samples were not re-tested with the revised 478K probe that was designed to also detect the mutation in this variant. The N2RP assay was performed as described [20,27]. For full Spike SNP methods, please see the Supplementary Materials. This study was performed in accordance with the STARD 2015 reporting guidelines (accessed 16 January 2022).

2.3. Sequencing

All rRT-PCR products from the first Spike SNP reaction (348-base-pair amplicon) were stored at 4 °C and later shipped to GeneWiz (South Plainfield, NJ, USA) for Sanger sequencing. This was performed to confirm the Spike SNP calls in the current study. The amplicon sequence was considered to be high quality and included for analysis if it achieved a quality score ≥ 40 and had base calls for all probe targets in at least one direction. A quality score ≥ 40, as defined by the company, signifies a mean probability of error ≤ 0.01% for each base call across the amplicon.
A representative subset of 29 samples underwent complete genome sequencing to confirm the variant lineage and identify mutations outside of the Spike SNP amplicon. Samples were selected to include the ancestral strain; alpha, gamma, and zeta variants; and one sample with the L452R mutation. Complete genome sequencing was performed following the Artic Network protocol using the V3 primer scheme [28]. The genomic libraries were prepared and then loaded on an R9.4.1 flow cell (ONT FLO-MIN106) for sequencing with a MinION device (MIN-101B, ONT, Oxford, UK). Basecalling, demultiplexing, and trimming were carried out using the Guppy toolkit integrated into the MinKNOW v.4.3.4 Nanopore software ( ONT, Oxford, UK). Genome assembly and consensus sequences were obtained with Nanopore EPI2ME v.3.4.1(ONT, Oxford, UK) (using Flye and Medaka respectively). The Nextclade tool was used for clade definition and mutation calling [29] and lineages were assigned using the Pangolin web server [30]. Sequences were deposited in the GISAID database.

2.4. Statistics

Basic statistical analyses were performed using GraphPad Prism version 9.0.2 (GraphPad Software, San Diego, CA, USA).

3. Results

3.1. Spike SNP Optimization

For the initial testing, a Spike SNP assay was developed that contained probes for 452Q, 484Q, and 490S. With the global emergence of B.1.617.2 and the lack of 484Q detection, the 478K probe was incorporated into the assay in place of 484Q. The addition of the original or revised 478K probe had no impact on 452Q or 490S detection (Figure 2), and despite overlapping sequences, the 478K and 484Q probes demonstrated no interference when combined in the same reaction (Figure S1). New probes for 484K and 501Y without LNA bases provided sensitive detection at both sites, though Ct values with the unmodified 484K probe were approximately two cycles later than those generated with the LNA probe (Figure S2). Testing of clinical samples in the current study was performed with the original 484K and 501Y probes to maintain consistency with existing protocols. All primers and probes used in or developed for this study are shown in Table 1, and their respective locations across the amplicon are displayed in Figure 1.

3.2. SNP Detection in Clinical Samples

Clinical samples were collected between November 2020 and April 2021. A total of 113 participants were female (56.8% of the 199 with gender data available), and the age and day of symptoms at sample collection are shown in Table 2, along with the N2RP and Spike SNP assay results. The distribution of N2 target Ct values is shown in Figure 3 (Ct values provided in Table S1). The most common genotype identified was 417variant/484K/501Y (102/198 samples, 51.5%), which was consistent with B.1.351 (beta) and P.1 (gamma). Large numbers of cases were also detected with a signal for K417 only (64/198, 32.3%), which was consistent with a non-variant/ancestral lineage, and K417/484K (22/198, 11.1%), which was observed in P.2 (zeta), among others. Seven samples (3.5%) had the genotype K417/452R but tested negative with the 478K probe and were thus not consistent with B.1.617.2.
Samples with the 417variant/484K/501Y genotype had significantly lower N2 Ct values (mean 19.3, standard deviation 5.0) than all other samples (22.0, 5.4; p < 0.001; data not shown) and samples with K417 only (22.1, 5.4; p = 0.002; Figure 3B). No differences in patient age or day of symptoms were observed based on the genotype (p > 0.05 for all comparisons).

3.3. Amplicon and Whole-Genome Sequencing

High-quality amplicon sequences were obtained from 181 samples (181/201, 90.0%) using Sanger sequencing, and the results confirmed Spike SNP genotype calls for all samples (Table 2). Characteristic findings for samples with 484K, with and without 501Y, are shown (Figure 4). All samples with variant calls at position 417 had the mutation AAG > ACG conferring K417T were observed in P.1. One sample had 484K detected in the Spike SNP assay (genotype K417/484K) but E484 was detected based on the consensus Sanger sequence. On review of the tracing, there was a mixture of G and A bases at the first position of codon 484, which is consistent with the evidence of a mixed viral population in this sample. Notably, this sample had a late Ct value for 484K (38.4) relative to K417 (34.3) when compared to other samples with this genotype (mean Ct difference 0.4, standard deviation 0.2). In addition to mutations targeted in the Spike SNP assay, four samples (2.2%) had a mutation conferring T478R (ACA > AGA) and one sample (0.5%) had a mutation resulting in N501T (AAU > ACU).
Whole-genome sequences for 29 samples were obtained using the Artic protocol. Ten different lineages were identified (Table 3) [30]. Nanopore sequencing confirmed Spike SNP genotype calls for all samples. Samples bearing genotype 417variant/484K/501Y belonged to the P.1 lineage and P.1.2 sub-lineage. The genotype of only K417 was found in samples of lineages B.1.1.28, B.1.1.33, B.1.1.277, and N3, all within the 20B Nextstrain clade and lineage B.1.499 (clade 20C) [29]. The P.2 lineage that also belongs to the 20B clade showed the K417/484K genotype. One sample with genotype K417/501Y was classified as B.1.1.7 (alpha). A lineage from 2020 that was still circulating during 2021, namely, A.2.5.2, was found in one sample with genotype K417/452R.

4. Discussion

Among individuals with a SARS-CoV-2 infection in Paraguay, the Spike SNP assay provided accurate detection of targeted mutations in the RBD and, in conjunction with the employed sequencing protocol, confirmed the predominance of B.1.1.28, P.2, and P.1 in this population between November 2020 and April 2021, as well as the relative absence of B.1.1.7. Given the frequent travel between Paraguay and Brazil, the B.1.1.28 lineage was expected to account for a high proportion of cases in late 2020 and early 2021, to be later replaced by P.1 and P.2. Samples with P.1 had lower Ct values in our population. A similar finding was initially documented in Brazil and likely resulted in increased transmissibility of the gamma variant with the replacement of other lineages [31,32]. The P.1 and P.2 lineages were readily differentiated from one another in the Spike SNP assay, and in a real-world scenario, the N2RP and Spike SNP assays could be performed sequentially to diagnose SARS-CoV-2 infection, identify a set of common mutations, and select samples for downstream sequencing [4,5,20].
For the current study, new probes were generated for mutations that occurred within the target region and are associated with emerging VOCs/VOIs, including delta, lambda, and omicron variants (Figure 1). As new designs do not involve separate primer–probe sets, the assay can be readily modified and reconfigured to include probes that match circulating variants in a particular region. All new probes were designed without LNA bases yet demonstrated similar performances to previous designs (Figure S1). Although LNA bases were shown to be important for the detection of K417 mutations, avoiding their use in other probes improves access to and reduces costs for these oligos. An additional benefit to the current design is that up to four mutations are targeted in a single reaction, which has been found to improve variant classification [9]. Other published rRT-PCR protocols for the detection of a similar number of mutations generate two or more amplicons in single-reaction [12,15,17] or multiple-reaction designs [6,8,14,16], and such protocols may be more difficult to modify in response to the emergence of new variants.
As shown in the current study, the Spike SNP assay can be paired with whole-genome sequencing on a MinION device to provide an economical and tractable molecular protocol for SARS-CoV-2 variant detection. Nanopore sequencing has become a common method worldwide and particularly in locations where the use of short-read sequencing platforms is less feasible. MinION sequencing was well described for SARS-CoV-2 and provides formal confirmation of detected variants and lineages, which cannot be achieved using rRT-PCR alone [4,5,25,31]. However, the Spike SNP assay may provide more sensitive detection of specific mutations than can be achieved with current sequencing protocols. Amplicon sequencing was performed in the current study to confirm Spike SNP calls for assay validation. While this provides less information than whole-genome sequencing, it serves to differentiate B.1.351 from P.1 and identify additional mutations in the receptor-binding motif. In our study, four individuals had the T478R mutation; the whole genome of two of them was sequenced and assigned to the B.1.499 lineage. This lineage circulated widely in Argentina during 2020, but the mutation T478R was not reported. While this mutation is distinct from T478K found in the B.1.617.2 (delta) and B.1.1.529 (omicron), mutations at position 478 have developed in vitro and confer decreased antibody neutralization, indicating that such mutations may benefit the virus through immunity escape [33,34]. Finally, one individual was identified with a possible mixed infection or a minor variant bearing 484K [35,36]. Such infections do not appear to cause more severe clinical disease [35], but this case highlights one potential benefit of rRT-PCR surveillance if key mutations can be detected that are not observed in the consensus sequence.
The degree of multiplexing possible in the Spike SNP assay is limited by the number of channels available in common real-time PCR instruments (4–5), as well as the defined 348-base-pair target region. Based on the results of testing for 478K and 484Q in a single reaction, overlapping probes did not impact the performance (Figure S1). However, multiplexed detection of key mutations outside of the Spike SNP target region would require additional primer–probe sets and assay optimization. A limitation of the current study was the utilization of retrospective convenience samples to evaluate the testing protocol. Although new data on SARS-CoV-2 variants in Paraguay were generated, we cannot calculate variant prevalence from these data.
In conclusion, the Spike SNP assay provided accurate detection of mutations that were associated with VOCs/VOIs in Paraguay. This can be implemented in SARS-CoV-2 testing protocols to triage samples for whole-genome sequencing in the MinION platform or direct amplicon sequencing, which provides additional information. The assay can be implemented in any laboratory performing rRT-PCR to improve the surveillance for these mutations and inform the judicious use of scarce sequencing resources.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v14050873/s1, File S1: Spike SNP Performance; Table S1: Spike SNP and N2RP line data for 201 clinical samples tested from Paraguay; Figure S1: 478K and 484Q probes do not interfere despite overlapping sequences; Figure S2: Performance of unmodified hydrolysis probes compared to probes that include LNA bases.

Author Contributions

Conceptualization, M.M., F.C. and A.R.; Data curation, M.M., A.R. and J.J.W.; Formal analysis, M.M., L.M., A.R. and J.J.W.; Funding acquisition, M.M., A.R. and J.J.W.; Investigation, M.M., P.-V.N., M.S., F.C., A.V., L.F., M.E.G., L.E.R., C.C.D.A., J.F., J.O., F.d.P., L.M., E.N., A.R. and J.J.W.; Methodology, M.M., P.-V.N., M.S., F.C., A.V., L.F., M.E.G., L.E.R., C.C.D.A., J.F., J.O., F.d.P., L.M., E.N., A.R. and J.J.W.; Project administration, M.M. and A.R.; Resources, A.R.; Supervision, M.M.; Writing—original draft, M.M., A.R. and J.J.W.; Writing—review and editing, M.M., P.-V.N., M.S., F.C., A.V., L.F., M.E.G., L.E.R., C.C.D.A., J.F., J.O., F.d.P., L.M., E.N., A.R. and J.J.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by funding from Consejo Nacional de Ciencia y Tecnología with FEEI support, CONACYT, Paraguay, PINV20-239 (M.M.); National Institutes of Health grants U54 EB027690 02S1, U54 EB027690 03S1, U54EB027690 03S2, and UL1TR002378; and the Doris Duke Charitable Foundation, Clinical Scientist Development Award 2019089 (J.J.W.).

Institutional Review Board Statement

This study was reviewed and approved by the IICS Scientific and Ethics Committee (P38/2020) and the Emory Institutional Review Board (study 00110736).

Informed Consent Statement

Patient consent was waived under approval of IICS Scientific and Ethics Committee (P38/2020) due to all samples being de-identified remaining of diagnostics samples.

Data Availability Statement

All supporting data necessary for interpretation are provided herein. Complete genome sequences were deposited in GISAID and accession numbers are shown in Table 3. Amplicon sequence tracings are available upon request.

Acknowledgments

The authors thank the staff at the Instituto de Investigaciones en Ciencias de la Salud, Universidad Nacional de Asunción, and collaborating laboratories: Letizia Carpinelli (CYRLAB), Fátima Ovando (Hospital de Clínicas, UNA), Graciela Riera (MeyerLab), and Stefan Goertzen (Hospital Loma Plata) for their support during the study period and the COVID-19 pandemic in Paraguay. We also thank members of the Variant Task Force at Emory University, namely, Aunradha Rao, Leda Bassit, and Morgan Greenleaf, for the provision of samples with characterized variant strains for assay evaluation.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Centers for Disease Control and Prevention SARS-CoV-2 Variant Classifications and Definitions. Available online: https://www.cdc.gov/coronavirus/2019-ncov/variants/variant-info.html#Concern (accessed on 2 September 2021).
  2. Gardy, J.L.; Loman, N.J. Towards a genomics-informed, real-time, global pathogen surveillance system. Nat. Rev. Genet. 2018, 19, 9–20. [Google Scholar] [CrossRef] [PubMed]
  3. Vavrek, D.; Speroni, L.; Curnow, K.J.; Oberholzer, M.; Moeder, V.; Febbo, P.G. Genomic surveillance at scale is required to detect newly emerging strains at an early timepoint. medRxiv 2021. [Google Scholar] [CrossRef]
  4. Pan American Health Organization. Recommendations for Reporting and Notification of Sars-Cov-2 Variants of Concern and Variants of Interest; Pan American Health Organization: Washington, WA, USA, 2021. [Google Scholar]
  5. World Health Organization. Guidance for Surveillance of Sars-Cov-2 Variants; World Health Organization: Geneva, Switzerland, 2021. [Google Scholar]
  6. Banada, P.; Green, R.; Banik, S.; Chopoorian, A.; Streck, D.; Jones, R.; Chakravorty, S.; Alland, D. A simple RT-PCR melting temperature assay to rapidly screen for widely circulating SARS-CoV-2 variants. J. Clin. Microbiol. 2021, 59, e00845-21. [Google Scholar] [CrossRef] [PubMed]
  7. Bedotto, M.; Fournier, P.E.; Houhamdi, L.; Levasseur, A.; Delerce, J.; Pinault, L.; Padane, A.; Chamieh, A.; Tissot-Dupont, H.; Brouqui, P.; et al. Implementation of an in-house real-time reverse transcription-PCR assay for the rapid detection of the SARS-CoV-2 Marseille-4 variant. J. Clin. Virol. 2021, 139, 104814. [Google Scholar] [CrossRef]
  8. Blairon, L.; Cupaiolo, R.; Piteus, S.; Beukinga, I.; Tre-Hardy, M. The challenge of screening SARS-CoV-2 variants of concern with RT-qPCR: One variant can hide another. J. Virol. Methods 2021, 297, 114248. [Google Scholar] [CrossRef]
  9. Camp, J.V.; Buchta, C.; Jovanovic, J.; Puchhammer-Stockl, E.; Benka, B.; Griesmacher, A.; Aberle, S.W.; Goerzer, I. RT-PCR based SARS-CoV-2 variant screening assays require careful quality control. J. Clin. Virol. 2021, 141, 104905. [Google Scholar] [CrossRef]
  10. Heijnen, L.; Elsinga, G.; de Graaf, M.; Molenkamp, R.; Koopmans, M.P.G.; Medema, G. Droplet digital RT-PCR to detect SARS-CoV-2 signature mutations of variants of concern in wastewater. Sci. Total Environ. 2021, 799, 149456. [Google Scholar] [CrossRef]
  11. La Rosa, G.; Mancini, P.; Bonanno Ferraro, G.; Veneri, C.; Iaconelli, M.; Lucentini, L.; Bonadonna, L.; Brusaferro, S.; Brandtner, D.; Fasanella, A.; et al. Rapid screening for SARS-CoV-2 variants of concern in clinical and environmental samples using nested RT-PCR assays targeting key mutations of the spike protein. Water Res. 2021, 197, 117104. [Google Scholar] [CrossRef]
  12. Norz, D.; Grunwald, M.; Olearo, F.; Fischer, N.; Aepfelbacher, M.; Pfefferle, S.; Lutgehetmann, M. Evaluation of a fully automated high-throughput SARS-CoV-2 multiplex qPCR assay with built-in screening functionality for del-HV69/70- and N501Y variants such as B.1.1.7. J. Clin. Virol. 2021, 141, 104894. [Google Scholar] [CrossRef]
  13. Tierling, S.; Kattler, K.; Vogelgesang, M.; Pfuhl, T.; Lohse, S.; Lo Porto, C.; Schmitt, B.; Nastasja, S.; Salhab, A.; Smola, S.; et al. Rapid base-specific calling of SARS-CoV-2 variants of concern using combined RT-PCR melting curve screening and SIRPH technology. Open Forum Infect. Dis. 2021, 8, ofab364. [Google Scholar] [CrossRef]
  14. Vega-Magana, N.; Sanchez-Sanchez, R.; Hernandez-Bello, J.; Venancio-Landeros, A.A.; Pena-Rodriguez, M.; Vega-Zepeda, R.A.; Galindo-Ornelas, B.; Diaz-Sanchez, M.; Garcia-Chagollan, M.; Macedo-Ojeda, G.; et al. RT-qPCR assays for rapid detection of the N501Y, 69–70del, K417N, and E484K SARS-CoV-2 mutations: A screening strategy to identify variants with clinical impact. Front. Cell. Infect. Microbiol. 2021, 11, 672562. [Google Scholar] [CrossRef]
  15. Vogels, C.B.F.; Breban, M.I.; Ott, I.M.; Alpert, T.; Petrone, M.E.; Watkins, A.E.; Kalinich, C.C.; Earnest, R.; Rothman, J.E.; Goes de Jesus, J.; et al. Multiplex qPCR discriminates variants of concern to enhance global surveillance of SARS-CoV-2. PLoS Biol. 2021, 19, e3001236. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, H.; Jean, S.; Eltringham, R.; Madison, J.; Snyder, P.; Tu, H.; Jones, D.M.; Leber, A.L. Mutation-specific SARS-CoV-2 PCR screen: Rapid and accurate detection of variants of concern and the identification of a newly emerging variant with spike L452R mutation. J. Clin. Microbiol. 2021, 59, e0092621. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, H.; Miller, J.A.; Verghese, M.; Sibai, M.; Solis, D.; Mfuh, K.O.; Jiang, B.; Iwai, N.; Mar, M.; Huang, C.; et al. Multiplex SARS-CoV-2 genotyping reverse transcriptase PCR for population-level variant screening and epidemiologic surveillance. J. Clin. Microbiol. 2021, 59, e0085921. [Google Scholar] [CrossRef] [PubMed]
  18. Yaniv, K.; Ozer, E.; Shagan, M.; Lakkakula, S.; Plotkin, N.; Bhandarkar, N.S.; Kushmaro, A. Direct RT-qPCR assay for SARS-CoV-2 variants of concern (Alpha, B.1.1.7 and Beta, B.1.351) detection and quantification in wastewater. Environ. Res. 2021, 201, 111653. [Google Scholar] [CrossRef]
  19. Zelyas, N.; Pabbaraju, K.; Croxen, M.A.; Lynch, T.; Buss, E.; Murphy, S.A.; Shokoples, S.; Wong, A.; Kanji, J.N.; Tipples, G. Precision response to the rise of the SARS-CoV-2 B.1.1.7 variant of concern by combining novel PCR assays and genome sequencing for rapid variant detection and surveillance. Microbiol. Spectr. 2021, 9, e0031521. [Google Scholar] [CrossRef]
  20. Babiker, A.; Immergluck, K.; Stampfer, S.D.; Rao, A.; Bassit, L.; Su, M.; Nguyen, V.; Stittleburg, V.; Ingersoll, J.M.; Bradley, H.L.; et al. Single-amplicon, multiplex real-time RT-PCR with tiled probes to detect SARS-CoV-2 spike mutations associated with variants of concern. J. Clin. Microbiol. 2021, 59, e014462. [Google Scholar] [CrossRef]
  21. Ministerio de Salud Pública y Bienestar Social. Primer Caso del Nuevo Coronavirus en El Paraguay. Available online: https://www.mspbs.gov.py/portal/20535/primer-caso-del-nuevo-coronavirus-en-el-paraguay.html (accessed on 1 September 2021).
  22. Ministerio de Salud Pública y Bienestar Social. Available online: https://www.mspbs.gov.py/index.php (accessed on 25 August 2021).
  23. Sobarzo, P.; Rolón López, J.C.; Narváez Serra, P.F.; López Cañete, S.A. Clinical characteristics of the first 60 patients with SARS CoV-2 admitted to the National Hospital for the period June–August 2020. Rev. Virtual Soc. Parag. Med. Int. 2021, 8, 69–77. [Google Scholar] [CrossRef]
  24. Vargas-Correa, A.; Mereles, E.F.; Segovia, C.N.; Gimenez, A.A.; Santacruz, L.; Ojeda, M.L.; Kunzle, E.H.; Samudio, M. Clinical-epidemiological characteristics of patients confirmed with COVID-19 from the Department of Alto Paraná, Paraguay. Rev. Salud Pública Parag. 2021, 11, 54–61. [Google Scholar] [CrossRef]
  25. Brinkmann, A.; Ulm, S.L.; Uddin, S.; Forster, S.; Seifert, D.; Oehme, R.; Corty, M.; Schaade, L.; Michel, J.; Nitsche, A. AmpliCoV: Rapid whole-genome sequencing using multiplex PCR amplification and real-time Oxford nanopore MinION sequencing enables rapid variant identification of SARS-CoV-2. Front. Microbiol. 2021, 12, 651151. [Google Scholar] [CrossRef]
  26. Corman, V.M.; Landt, O.; Kaiser, M.; Molenkamp, R.; Meijer, A.; Chu, D.K.W.; Bleicker, T.; Brunink, S.; Schneider, J.; Schmidt, M.L.; et al. Detection of 2019 novel coronavirus (2019-nCoV) by real-time RT-PCR. Eur. Surveill. Bull. Eur. Mal. Transm. Eur. Commun. Dis. Bull. 2020, 25, 2000045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Waggoner, J.J.; Stittleburg, V.; Pond, R.; Saklawi, Y.; Sahoo, M.K.; Babiker, A.; Hussaini, L.; Kraft, C.S.; Pinsky, B.A.; Anderson, E.J.; et al. Triplex real-time RT-PCR for severe acute respiratory syndrome coronavirus 2. Emerg. Infect. Dis. 2020, 26, 1633–1635. [Google Scholar] [CrossRef]
  28. Quick, J. nCoV-2019 Sequencing Protocol v3 (LoCost) V.3. Available online: https://www.protocols.io/view/ncov-2019-sequencing-protocol-v3-locost-bh42j8ye (accessed on 1 February 2021).
  29. Hadfield, J.; Megill, C.; Bell, S.M.; Huddleston, J.; Potter, B.; Callender, C.; Sagulenko, P.; Bedford, T.; Neher, R.A. Nextstrain: Real-time tracking of pathogen evolution. Bioinformatics 2018, 34, 4121–4123. [Google Scholar] [CrossRef] [PubMed]
  30. Rambaut, A.; Holmes, E.C.; O’Toole, A.; Hill, V.; McCrone, J.T.; Ruis, C.; du Plessis, L.; Pybus, O.G. A dynamic nomenclature proposal for SARS-CoV-2 lineages to assist genomic epidemiology. Nat. Microbiol. 2020, 5, 1403–1407. [Google Scholar] [CrossRef] [PubMed]
  31. Faria, N.R.; Mellan, T.A.; Whittaker, C.; Claro, I.M.; Candido, D.D.S.; Mishra, S.; Crispim, M.A.E.; Sales, F.C.S.; Hawryluk, I.; McCrone, J.T.; et al. Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil. Science 2021, 372, 815–821. [Google Scholar] [CrossRef] [PubMed]
  32. Naveca, F.G.; Nascimento, V.; de Souza, V.C.; Corado, A.L.; Nascimento, F.; Silva, G.; Costa, A.; Duarte, D.; Pessoa, K.; Mejia, M.; et al. COVID-19 in Amazonas, Brazil, was driven by the persistence of endemic lineages and P.1 emergence. Nat. Med. 2021, 27, 1230–1238. [Google Scholar] [CrossRef]
  33. Singh, J.; Rahman, S.A.; Ehtesham, N.Z.; Hira, S.; Hasnain, S.E. SARS-CoV-2 variants of concern are emerging in India. Nat. Med. 2021, 27, 1131–1133. [Google Scholar] [CrossRef]
  34. 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. [Google Scholar] [CrossRef]
  35. Francisco, R.D.S., Jr.; Benites, L.F.; Lamarca, A.P.; de Almeida, L.G.P.; Hansen, A.W.; Gularte, J.S.; Demoliner, M.; Gerber, A.L.; de C Guimarães, A.P.; Antunes, A.K.E.; et al. Pervasive transmission of E484K and emergence of VUI-NP13L with evidence of SARS-CoV-2 co-infection events by two different lineages in Rio Grande do Sul, Brazil. Virus Res. 2021, 296, 198345. [Google Scholar] [CrossRef]
  36. Lythgoe, K.A.; Hall, M.; Ferretti, L.; de Cesare, M.; MacIntyre-Cockett, G.; Trebes, A.; Andersson, M.; Otecko, N.; Wise, E.L.; Moore, N.; et al. SARS-CoV-2 within-host diversity and transmission. Science 2021, 372, eabg082. [Google Scholar] [CrossRef]
Viruses 14 00873 g001 550
Figure 1. Graphical display of the 348-base-pair Spike SNP target, including the nucleic acid sequence and amino acid translation of the Wuhan-Hu-1 strain (NC_045512.2). The locations of the forward and reverse primer sequences and all hydrolysis probes are shown. Probe colors correspond to the fluorophores on each probe and the detection channel in the Rotor-Gene Q instrument. Detected nucleotide changes are displayed within the probe, and amino acid mutations are shown below the sequence. The probe for K417 is marked (*), as this is the only probe for which variant sequences result in signal loss and this has only been successfully designed with locked nucleic acid (LNA) bases.
Figure 1. Graphical display of the 348-base-pair Spike SNP target, including the nucleic acid sequence and amino acid translation of the Wuhan-Hu-1 strain (NC_045512.2). The locations of the forward and reverse primer sequences and all hydrolysis probes are shown. Probe colors correspond to the fluorophores on each probe and the detection channel in the Rotor-Gene Q instrument. Detected nucleotide changes are displayed within the probe, and amino acid mutations are shown below the sequence. The probe for K417 is marked (*), as this is the only probe for which variant sequences result in signal loss and this has only been successfully designed with locked nucleic acid (LNA) bases.
Viruses 14 00873 g001
Viruses 14 00873 g002 550
Figure 2. Addition of the 478K probe did not impact the signal from 452Q or 490S probes. Amplification curves are shown for 100-fold dilutions of C.37 (lambda) and B.1.617.2 (delta) samples tested together in a single run. The results were similar for (A) 452Q and (C) 490S mutations in C.37 with the probe for 478K (solid curves) or without (dotted curves). Results for B.1.617.2 were similar with the 478K probe (B) used in singleplex (dotted curves) or triplex reactions (solid curves). Nonspecific signals were not detected in any channel.
Figure 2. Addition of the 478K probe did not impact the signal from 452Q or 490S probes. Amplification curves are shown for 100-fold dilutions of C.37 (lambda) and B.1.617.2 (delta) samples tested together in a single run. The results were similar for (A) 452Q and (C) 490S mutations in C.37 with the probe for 478K (solid curves) or without (dotted curves). Results for B.1.617.2 were similar with the 478K probe (B) used in singleplex (dotted curves) or triplex reactions (solid curves). Nonspecific signals were not detected in any channel.
Viruses 14 00873 g002
Viruses 14 00873 g003 550
Figure 3. Distribution of N2 Ct values for (A) all samples and (B) the four most common genotypes detected. (A) N2 Ct values for all tested samples are displayed. Three samples that tested negative in the Spike SNP assay are highlighted (dark triangles) and had N2 Ct values of 32.6, 36.6, and 35.1. (B) N2 Ct values for the four most common genotypes detected. Bars represent mean Ct and standard deviation, and p-values for tested comparisons are displayed above the graph.
Figure 3. Distribution of N2 Ct values for (A) all samples and (B) the four most common genotypes detected. (A) N2 Ct values for all tested samples are displayed. Three samples that tested negative in the Spike SNP assay are highlighted (dark triangles) and had N2 Ct values of 32.6, 36.6, and 35.1. (B) N2 Ct values for the four most common genotypes detected. Bars represent mean Ct and standard deviation, and p-values for tested comparisons are displayed above the graph.
Viruses 14 00873 g003
Viruses 14 00873 g004 550
Figure 4. Spike SNP and amplicon sequencing results were concordant. Amplification curves and amplicon sequence tracings are shown for samples with mutations encoding (A) 484K and 501Y (purple curves) or (B) 484K only with N501 (cranberry). A beta variant control sample is shown for reference (blue curves). Bold squares highlight the SNP mutations in the amplicon sequence.
Figure 4. Spike SNP and amplicon sequencing results were concordant. Amplification curves and amplicon sequence tracings are shown for samples with mutations encoding (A) 484K and 501Y (purple curves) or (B) 484K only with N501 (cranberry). A beta variant control sample is shown for reference (blue curves). Bold squares highlight the SNP mutations in the amplicon sequence.
Viruses 14 00873 g004
Table
Table 1. Spike SNP assay primers and probes used in the current study.
Table 1. Spike SNP assay primers and probes used in the current study.
NameSequence aConcentration (nM) bLocation (5′–3′) c
Primers
SpikeSNP Ext dTGAAGTCAGACAAATCGCTCC40022,777–22,797
SpikeSNP R1 dTGGTGCATGTAGAAGTTCAAAAG40023,103–23,125
LNA Probes
K417 dFAM-CAAACTGGA+A+A+G+ATTGCTG-3IABkFQ20022,802–22,820
484K dCy5-ACCTTGTAATGGTGT+T+A+AAGGTTTT-3IAbRQSp20022,996–23,020
501Y dHEX-CCCAC+T+T+ATGGTGTTGG-3IABkFQ20023,057–23,073
Unmodified Probes
452R dCFR610-ATAATTACCGGTATAGATTGTTTAGGAAGT-BHQ-220022,908–22,937
484KQ670-CACCTTGTAATGGTGTTAAAGGTTTTAA-BHQ-220022,995–23,022
484QQ670-CACCTTGTAATGGTGTTCAAGGTTTTAA-BHQ-220022,995–23,022
501YCFO560-CAACCCACTTATGGTGTTGGTTAC-BHQ-120023,054–23,077
452QCFR610-ATAATTACCAGTATAGATTGTTTAGGAAGT-BHQ-220022,908–22,937
490SFAM-ACTCTCCTTTACAATCATATGGTTTCC-BHQ-120023,028–23,054
478K (original)CFO560-CGGTAGCAAACCTTGTAATGGTG-BHQ-120022,987–23,009
478K (revised)CFO560-CGGTARCAAACCTTGTAATGGTG-BHQ-120022,987–23,009
Abbreviations: BHQ, black hole quencher; CFO560, CAL Fluor Orange 560; CFR610, CAL Fluor Red 610; Cy5, Cyanine 5; FAM, Fluorescein; IABkFQ/IAbRQSp, Iowa Black quenchers; Q670, Quasar 670; a probe sequences listed 5′-fluorophore-sequence-quencher-3′; “+” before base indicates a locked nucleic acid, bold and underlined base indicates the targeted single nucleotide change; b concentration in the final reaction mixture; c location in Wuhan-Hu-1 complete genome sequence (NCBI Reference Sequence: NC_045512.2); d previously published, presented for clarity.
Table
Table 2. Genotypes detected in the Spike SNP assay, patient demographics, and the percent confirmed using amplicon sequencing.
Table 2. Genotypes detected in the Spike SNP assay, patient demographics, and the percent confirmed using amplicon sequencing.
Categoryn (%)Age, Years, Mean (SD) aDay of Symptoms, Mean (SD) bCt, Mean (SD) cSequencing Confirmed, n/N (%) d
N2 detected201 (100)42.6 (16.4)4.8 (3.0)20.8 (5.6)─
Spike SNP detected198 (98.5)42.5 (16.4)4.8 (3.0)20.6 (5.3)181/181 (100)
Genotype e
417var/484K/501Y102 (51.5)42.8 (14.9)4.5 (2.4)19.3 (5.0)96/96 (100)
K41764 (32.3)42.6 (18.8)5.1 (3.7)22.1 (5.4)56/56 (100)
K417/484K22 (11.1)37.4 (14.2)4.8 (2.7)22.2 (5.4)20/20 (100) f
K417/452R7 (3.5)54.0 (15.8)5.0 (4.6)19.4 (4.0)7/7 (100)
K417/490S1 (0.5)41415.21/1 (100)
K417/501Y1 (0.5)68519.71/1 (100)
K417/484K/501Y1 (0.5)191032.1NA
Abbreviations: var, variant; SD, standard deviation; a age available for 194/201 (96.5%) samples; b day of symptoms available for 136/201 (67.7%) samples; c Ct value in the N2 assay; d number with sequence confirmed (n) over the number of amplicons with high-quality sequence and base calls at each SNP position (N); e displayed as % of Spike SNP positive samples; f one sample had a mixed base (A/G) at the position targeted with the 484K probe.
Table
Table 3. Spike mutations and lineages identified using genome sequencing.
Table 3. Spike mutations and lineages identified using genome sequencing.
Sample CodeSpike SNP
Genotype		Nanopore Sequencing
GISAID
Accession No.		Nextstrain CladeWHO LabelPango LineageSpike Deletions/Amino Acid Substitutions a
PY21-100417var/484K/501YEPI_ISL_407189220J (V3)
gamma		P.1L18F/T20N/P26S/D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-9417var/484K/501YEPI_ISL_244477820J (V3)
gamma		P.1L18F//T20N/P26S/D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-57417var/484K/501YEPI_ISL_413747620J (V3)
gamma		P.1L18F/T20N/P26S/D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-58417var/484K/501YEPI_ISL_407189720J (V3)
gamma		P.1L18F/T20N/P26S7D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-84417var/484K/501YEPI_ISL_407189820J (V3)
gamma		P.1L18F/T20N/P26S/D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-16417var/484K/501YEPI_ISL_407189320J (V3)
gamma		P.1L18F/T20N/P26S/D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-11417var/484K/501YEPI_ISL_407190020J (V3)
gamma		P.1L18F/T20N/P26S/D138Y/L189F/R190S/K417T/E484K/N501Y/D614G/H655Y/T1027I/V1176F
PY21-90417var/484K/501YEPI_ISL_407189920J (V3)
gamma		P.1.2L18F/T20N/P26S/T95I/D138Y/R190S/K417T/E484K/N501Y/D614G/H655Y/Q675H/T1027I/V1176F/V1228L
PY21-150K417EPI_ISL_223489920BB.1.1.33D614G
PY21-153K417EPI_ISL_223489120BB.1.1.33D614G
PY21-140K417EPI_ISL_223488020BB.1.1.28D614G/V1176F
PY21-145K417EPI_ISL_408460520BB.1.1.28D614G/V1176F
PY21-146K417EPI_ISL_223488120BB.1.1.28D614G/V1176F
PY21-98K417EPI_ISL_244478820BB.1.1.28D614G/V1176F
PY21-152K417EPI_ISL_223488520BB.1.1.28T22I/D614G/V1176F
PY21-144K417EPI_ISL_407190120BB.1.1.28D614G/N1135Y/V1176F
PY21-43K417EPI_ISL_407189520BB.1.1.28A575S/D614G/V1176F
PY21-44K417EPI_ISL_408460420BB.1.1.28T572I/D614G/V1176F
PY21-47K417EPI_ISL_407189420BB.1.1.28A575S/D614G/V1176F
PY21-148K417EPI_ISL_223489420BB.1.1.28V6F/L18F/D614G/V1176F
PY21-147K417EPI_ISL_223488720BB.1.1.28Q14H/D614G/D1153A/V1176F
PY21-154K417EPI_ISL_223489320BB.1.1.277D614G/V1176F
PY21-52K417/484KEPI_ISL_244478020B
(zeta)		P.2E484K/D614G/V1176F
PY21-157K417/484KEPI_ISL_223487920B
(zeta)		P.2L5F/E484K/D614G/A771S/V1176F
PY21-143K417EPI_ISL_223488920BN.3T76I/D614G
PY21-141K417EPI_ISL_223488320CB.1.499T478R/D614G
PY21-149K417EPI_ISL_223489720CB.1.499T478R/D614G
PY21-102K417/452REPI_ISL_407189619BA.2.5.2L452R/D614G
PY21-156K417/501YEPI_ISL_408460620I (V1)
alpha		B.1.1.7DEL69-70/DEL144N501Y/A570D/D614G/P681H/T716I/S982A/D1118H
Abbreviations: var, variant; WHO, World Health Organization; a bold text indicates amino acid substitutions detected using the Spike SNP assay.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Martinez, M.; Nguyen, P.-V.; Su, M.; Cardozo, F.; Valenzuela, A.; Franco, L.; Galeano, M.E.; Rojas, L.E.; Díaz Acosta, C.C.; Fernández, J.; et al. SARS-CoV-2 Variants in Paraguay: Detection and Surveillance with an Economical and Scalable Molecular Protocol. Viruses 2022, 14, 873. https://doi.org/10.3390/v14050873

AMA Style

Martinez M, Nguyen P-V, Su M, Cardozo F, Valenzuela A, Franco L, Galeano ME, Rojas LE, Díaz Acosta CC, Fernández J, et al. SARS-CoV-2 Variants in Paraguay: Detection and Surveillance with an Economical and Scalable Molecular Protocol. Viruses. 2022; 14(5):873. https://doi.org/10.3390/v14050873

Chicago/Turabian Style

Martinez, Magaly, Phuong-Vi Nguyen, Maxwell Su, Fátima Cardozo, Adriana Valenzuela, Laura Franco, María Eugenia Galeano, Leticia Elizabeth Rojas, Chyntia Carolina Díaz Acosta, Jonás Fernández, and et al. 2022. "SARS-CoV-2 Variants in Paraguay: Detection and Surveillance with an Economical and Scalable Molecular Protocol" Viruses 14, no. 5: 873. https://doi.org/10.3390/v14050873

APA Style

Martinez, M., Nguyen, P. -V., Su, M., Cardozo, F., Valenzuela, A., Franco, L., Galeano, M. E., Rojas, L. E., Díaz Acosta, C. C., Fernández, J., Ortiz, J., del Puerto, F., Mendoza, L., Nara, E., Rojas, A., & Waggoner, J. J. (2022). SARS-CoV-2 Variants in Paraguay: Detection and Surveillance with an Economical and Scalable Molecular Protocol. Viruses, 14(5), 873. https://doi.org/10.3390/v14050873

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