Next Article in Journal
Students in Dormitories Were Not Major Drivers of the Pandemic during Winter Term 2020/2021: A Cohort Study with RT-PCR and Antibody Surveillance in a German University City
Previous Article in Journal
COVID. A New Journal to Affirm the Role of Science between Communication and Social Responsibility
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A Simple and Fast Method to Sequence the Full-Length Spike Gene for SARS-CoV-2 Variant Identification from Patient Samples

1
Microbiologie Fondamentale et Pathogénicité, CNRS-UMR 5234, Université de Bordeaux, F-33000 Bordeaux, France
2
SBM Department, UB’L3 Platform, Université de Bordeaux, F-33000 Bordeaux, France
3
CHU Bordeaux, Virology Laboratory, Université de Bordeaux, F-33000 Bordeaux, France
*
Author to whom correspondence should be addressed.
Submission received: 8 July 2021 / Revised: 26 August 2021 / Accepted: 1 September 2021 / Published: 3 September 2021
(This article belongs to the Topic Burden of COVID-19 in Different Countries)

Abstract

:
Since the beginning of the pandemic, a race has been underway to detect SARS-CoV-2 virus infection (PCR screening, serological diagnostic kits), treat patients (drug repurposing, standard care) and develop a vaccine. After almost a year of active circulation worldwide, SARS-CoV-2 variants have appeared in different countries. Those variants include mutations in multiple regions of the genome, particularly in the spike gene. Because this surface protein is a key player in both the spread of the virus and the efficacy of vaccine strategies, the challenge is to efficiently monitor the appearance of spike mutations in the population. The present work describes a procedure based on the widely available Sanger technology to produce a full-length sequence of the spike gene from patient-derived samples.

1. Introduction

In late 2019, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in Wuhan, Hubei province, China. This local outbreak rapidly developed into the COVID-19 pandemic [1,2]. At the end of March 2021, the total number of cases was estimated to be 128 million (https://coronavirus.jhu.edu, 29 March 2021). Since the middle of 2020, there have been hundreds of thousands of new infections each day, with peaks approaching one million. Unsurprisingly, several new variants have emerged during this expansion, some of which have become variants of interest, or even variants of concern (VOC) [3]. Their prevalence in the population quickly increased and the majority of new infections are now due to these variants rather than the parental strain. Mutations have been described in the 5′ untranslated region (UTR), ORF1ab, ORF3a, Spike, M, ORF8 and N. Therefore, monitoring the genetic evolution in the SARS-CoV-2 virus is becoming a top priority in the management of the COVID-19 pandemic.
The entry of coronaviruses is mediated by the surface glycoprotein spike [4,5,6]. Because of its role in receptor binding, the spike protein was quickly found to be the best candidate for producing neutralizing antibodies [7]. Current vaccines and novel therapeutic strategies under development are based on this spike protein. Accordingly, it is crucial to monitor the appearance of mutations in circulating viruses that might affect collective as well as acquired immunity.
The first described amino acid mutation was a D to G transition in the spike coding region (D614G) [8]. Observed in Europe as early as February 2020, this variation increased viral replication/infectivity [9]. It is now present in almost all circulating viruses and is shared among all of the current variants of interest/concern. Of note, this mutation has also been associated with an increased susceptibility to neutralizing antibodies [10]. Conversely, other mutations, such as the N501Y and E484K mutations located in the receptor-binding domain of spike (RBD), can affect the immunity conferred by polyclonal antibodies. The United Kingdom (UK) variant, also known as VOC 202012/01 (Nextstrain 20I/501Y.V1, PANGO lineage B.1.1.7), emerged in the south-east of the country in December 2020 [11]. In addition to the D614G, it harbored multiple mutations including N501Y, as well as a deletion at positions 69/70, affecting the diagnostic capacity of PCR tests involving the S region (Table 1, Alpha).
Soon thereafter, another VOC was identified with a similar N501Y transition but in South Africa [12]. Known as 20H/Beta or B.1.351, this variant (Table 1, Beta) appeared independently of the alpha virus and presented additional transitions including the spike E484K. Finally, a third VOC was identified in Brazil. Known as P.1 or 20J/501Y.V3, this variant shared some spike mutations seen in the beta virus. In all cases, numerous other mutations can be associated with these major variants, and new mutations are appearing constantly. That is the case, for example, with the BRE-IPP03921 variant, which recently emerged in Brittany, France, and which was identified after a recrudescence of COVID-19 symptoms in patients who tested negative for SARS-CoV-2 by routine PCR screening. More recently, the VOC delta (B.1.617.2), first described in India, is now spreading in the United Kingdom and all around Western Europe [13].
Overall, all variants harbour mutations in the viral spike, and these mutations spread throughout the protein sequence. However, present sequencing protocols of the spike cover two regions (genome positions 21,701 to 22,344 and 22,763 to 23,310), while sequencing of the whole spike is performed by NGS when a variant is suspected. Next generation sequencing would be the best solution to obtain a full-length spike sequence or even a whole genome sequence in order to survey any kind of mutation. Still, the required technology is absent in most southern countries and is not even readily available in every hospital in Europe. In contrast, Sanger is more widely available and makes it possible to obtain 800–1000 base reads per run. However, the SARS-CoV-2 spike is a 1273-amino acid protein encoded by the corresponding 3822-base open reading frame. Accordingly, we developed a procedure based on Sanger technology to sequence the full-length spike in three steps from nasal or saliva screening swab extracts. This procedure is compatible with other patient-derived samples such as bronchoalveolar lavages (BAL). The entire process can be performed in less than 14 h and provides sufficient coverage to identify existing and future novel variants.

2. Materials and Methods

Clinical Sample Collection. Patient samples were obtained from the Bordeaux university Hospital virology laboratory. To validate the protocol, various types of clinical samples were used. Twelve among 15 samples were RNA extracts from patients tested at the hospital COVID-19 centre as positive by qualitative real-time RT-PCR (nasal swabs). Three samples were RNA extracted from broncho-alveolar lavages (BAL). The last one (culture strain) consisted of RNA extracts from culture supernatant of Vero E6 cells infected with a reference strain (WIV04/2019 obtained through the EVAg and Pasteur Institute). Infection was performed in the BSL3 facility of Bordeaux University (UB’L3) as described in our recent work [14]. Ct values of RT-qPCR were determined during patient follow-up in the hospital using SARS-CoV-2 R-gene (BioMérieux, Craponne, France) and GeneFinderTM COVID-19 PLUS RealAmp (Osang Healthcare, Gyeonggi-do, Korea) kits.
RNA extraction. A pure viral RNA kit from Roche was used to extract RNA from samples of either BAL or cell culture supernatant (200 µL). After binding and washes, elution was performed with 50 µL of water. The quality of the extraction was checked by monitoring the absorbance profile using a NanoDrop 2000 (ThermoFisher scientific, Waltham, MA, USA).
Amplification of the spike gene. Primers used for amplification and sequencing were purchased from Eurofins Genomics (Germany), and are reported in Table 2. Amplification of a single fragment of 4200 bp covering the entire spike gene was performed from 10 µL of extracted viral RNA, as previously described [15]. RT-PCR was carried out in 25 μL reaction using forward F1-1 and reverse F1-1 primers (0.4 μM of each primer) in a buffer containing 0.4 mM of each dNTP, 3.2 mM MgSO4 and 0.5 unit of SuperScript™ III RT/Platinum™ Taq Mix (Invitrogen, ThermoFisher Scientific, Waltham, MA, USA). After an initial RT step performed at 50 °C for 30 min, the PCR cycling program consisted of a denaturation at 95 °C for 2 min followed by 30 cycles of 95 °C for 30 s, 60 °C for 15 s and 68 °C for 5 min. The reaction concluded with a final extension at 68 °C for 5 min. A nested PCR was performed with 2 µL of this RT-PCR product at a final volume of 25 μL using 0.8 μM of each primer (forward F1-2 and reverse F1) in a reaction buffer containing 15 mM of MgCl2, 5 mM of dNTPs and 0.25 unit of ThermoSTAR 2 HS Taq (Eurobio, Les Ulis, France). The quality of the PCR products was verified by electrophoresis using a QIAxcel Advanced automate (Qiagen, Germantown, MD, USA) before cleaning using ExoSAP-IT™ PCR Product Cleanup Reagent kit (Applied Biosystems, ThermoFisher Scientific, Waltham, MA, USA).
Sanger Sequencing and sequences analysis. Purified PCR products were sequenced using the BigDye® Terminator v3.1 Cycle Sequencing Kit on an Applied Biosystems 3500XL DX automat (Applied Biosystems, ThermoFisher Scientific, Waltham, MA, USA). For that, ten sequencing reactions were performed, each using 2 µL of nested PCR products, and using primers F1-2 and F1-2 plus eight additional primers (Table 2). Abi files were either aligned to the reference SARS-CoV-2 genome in SnapGene 4.3.11 or uploaded to the SmartGene IDNS software Coronavirus Module (IDNS® version v3_11_0r2(r32344) © SmartGene 2021).

3. Results

The aim of this work was to establish a robust but simple, fast procedure to generate sequencing data on the full-length spike gene using Sanger technology. Since the beginning of 2020, several hundreds of thousands of people have been screened every day for SARS-CoV-2 infection. Accordingly, there are a tremendous number of samples that merit further characterization through spike sequencing. The majority of the collected samples are nasal swabs that are further processed (RNA extraction) for RT-qPCR detection of two to three viral genes. We intended to make use of these PCR positive samples that were already extracted to enable subsequent full-length spike gene sequencing. Sequencing was performed for 16 samples from patients, comprising 8 males and 7 females. Ct values ranged from 15.14 to 35.76 for the N gene, and 16.6 to undetectable (Ct > 40) for the RdRp gene (Table 3). For some samples, E gene screening was also performed and Ct values ranged from 13.53 to 25.64 (Table 3).
First, RNA extracts from nasal swabs were subjected to reverse transcription using primers Forward F1-1 and Reverse F1-1 located in adjacent genes, namely, orf1ab and orf3a, respectively (Table 2). This led to the generation of a 4587-base-pair DNA fragment that was further amplified by nested PCR using primers Forward F1-2 and Reverse F1-2 (Table 2). At this step, the quality of amplification could be assessed if needed through conventional agarose separation or any other electrophoretic method. As seen in Figure 1B, this process from nasal swab extracts allowed us to produce a single band of 4199 bp. Before further processing, samples were enzymatically cleaned with Exo I and rSAP treatments. Sequencing was performed using 2 µL of the PCR mix per primer pair, allowing up to 13 sequencing reactions per sample. We found that six was the minimal number of Sanger reactions necessary to fully cover the 3822 bp spike gene. Deletion events in the 5′ side were unambiguously detected, even with a single read. However, mutations at primer binding sites may affect sequencing efficacy. In addition, the presence of mix bases due to multiple viral populations may interfere with proper base assignment, and higher coverage might be required in hotspot regions. Accordingly, primers were designed to produce a higher overlap (4x coverage) at spike amino acids positions 400 to 800 (22,760–23,962 nt, Figure 1A), because these correspond to the receptor binding domain and represent the most relevant region for variant determination. Thus, we proposed a set of 10 primers to unambiguously sequence the full-length spike gene. Using this protocol, complete sequencing was achieved for all tested samples (Table 3).
Sequencing could be achieved using various types of samples. In addition to nasal swabs, BAL or cell culture viral supernatants were used as alternative starting materials. In the latter case, RNA extraction was performed using 200 µL of sample concentrated to 50 µL of total RNA extract. Then, these alternative extracts were used exactly as previously described for the nasal swab extracts. Analyses of sequences were then performed using the SmartGene IDNS software Coronavirus Module (Figure 2). This allowed us to detect mutations in comparison with the reference strain. Combinations of mutations unambiguously allowed us to classify the viruses as WT or in relation to the circulating alpha and beta variants, with alpha representing the most frequent circulating form in our region at this time. In addition to previously reported mutations using current protocols, several additional mutations could be identified in the spike gene, especially the L452R, a mutation found in delta and kappa variants, as well as in epsilon, known to impact antibody binding and increase transmissibility. Interestingly, sequencing data could be obtained similarly from samples arising from screening facility samples, as well as from intensive care unit derived samples (BAL). Moreover, our data showed that if necessary, hard to sequence samples could alternatively be amplified in the lab (cell culture viral amplification) without introducing additional mutations (at least on early passages).

4. Discussion

Since the original breakout of SARS-CoV-2 at the end of 2019 in China, tremendous effort has been made by the scientific community to better understand the physiopathology of the virus and find effective antivirals and vaccines. Even though several vaccines have been successfully developed in record time, several mutations have arisen which threaten their future efficacy. Altogether, we developed a simple and reliable method to sequence whole spike gene in less than 14 h, directly from RNA extract left over from screening facilities. This protocol is able to unambiguously identify an exhaustive list of known and yet unseen mutations, which will help monitor circulating variants around the world if widely adopted by screening facilities. The surveillance of new variants will certainly be the most significant challenge in future months, and since next generation sequencing is not available in many regions, this simple method, based on Sanger sequencing, appears to be a viable option for global surveillance.

Author Contributions

Sample collection, P.B., M.-E.L. and C.T.; sample preparation, P.R.-P. and M.-L.B. methodology, P.R.-P.; formal analysis, P.R.-P. and M.M.; writing—original draft preparation, P.R.-P. and M.M.; review and editing, M.M. and M.-L.A.; funding acquisition, M.-L.A. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the CNRS and the Bordeaux University for support. This work was supported by the Region Nouvelle-Aquitaine (AMI FLASH Recherche et Innovations COVID). ML Blondot was supported by the Region Nouvelle Aquitaine and UBreACT Mission COVID-19 of Bordeaux University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw sequence data can be obtained directly from the authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhou, P.; Yang, X.L.; Wang, X.G.; Hu, B.; Zhang, L.; Zhang, W.; Si, H.R.; Zhu, Y.; Li, B.; Huang, C.L.; et al. A pneumonia outbreak associated with a new coronavirus of probable bat origin. Nature 2020, 579, 270–273. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zhu, N.; Zhang, D.; Wang, W.; Li, X.; Yang, B.; Song, J.; Zhao, X.; 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. Centers for Disease Control and Prevention. Available online: https://www.cdc.gov/coronavirus/2019-ncov/cases-updates/variant-surveillance/variant-info.html (accessed on 29 March 2021).
  4. Bosch, B.J.; van der Zee, R.; de Haan, C.A.; Rottier, P.J. The coronavirus spike protein is a class I virus fusion protein: Structural and functional characterization of the fusion core complex. J. Virol. 2003, 77, 8801–8811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Wrapp, D.; Wang, N.; Corbett, K.S.; Goldsmith, J.A.; Hsieh, C.L.; Abiona, O.; Graham, B.S.; McLellan, J.S. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 2020, 367, 1260–1263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ou, X.; Liu, Y.; Lei, X.; Li, P.; Mi, D.; Ren, L.; Guo, L.; Guo, R.; Chen, T.; Hu, J.; et al. Characterization of spike glycoprotein of SARS-CoV-2 on virus entry and its immune cross-reactivity with SARS-CoV. Nat. Commun. 2020, 11, 1620. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Piccoli, L.; Park, Y.J.; Tortorici, M.A.; Czudnochowski, N.; Walls, A.C.; Beltramello, M.; Silacci-Fregni, C.; Pinto, D.; Rosen, L.E.; Bowen, J.E.; et al. Mapping Neutralizing and Immunodominant Sites on the SARS-CoV-2 Spike Receptor-Binding Domain by Structure-Guided High-Resolution Serology. Cell 2020, 183, 1024–1042.e1021. [Google Scholar] [CrossRef] [PubMed]
  8. Korber, B.; Fischer, W.M.; Gnanakaran, S.; Yoon, H.; Theiler, J.; Abfalterer, W.; Hengartner, N.; Giorgi, E.E.; Bhattacharya, T.; Foley, B.; et al. Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus. Cell 2020, 182, 812–827.e819. [Google Scholar] [CrossRef] [PubMed]
  9. Zhou, B.; Thao, T.T.N.; Hoffmann, D.; Taddeo, A.; Ebert, N.; Labroussaa, F.; Pohlmann, A.; King, J.; Steiner, S.; Kelly, J.N.; et al. SARS-CoV-2 spike D614G change enhances replication and transmission. Nature 2021, 592, 122–127. [Google Scholar] [CrossRef] [PubMed]
  10. Weissman, D.; Alameh, M.G.; de Silva, T.; Collini, P.; Hornsby, H.; Brown, R.; LaBranche, C.C.; Edwards, R.J.; Sutherland, L.; Santra, S.; et al. D614G Spike Mutation Increases SARS CoV-2 Susceptibility to Neutralization. Cell Host Microbe 2021, 29, 23–31 e24. [Google Scholar] [CrossRef] [PubMed]
  11. Tang, J.W.; Tambyah, P.A.; Hui, D.S. Emergence of a new SARS-CoV-2 variant in the UK. J. Infect. 2021, 82, e27–e28. [Google Scholar] [CrossRef] [PubMed]
  12. Tegally, H.; Wilkinson, E.; Giovanetti, M.; Iranzadeh, A.; Fonseca, V.; Giandhari, J.; Doolabh, D.; Pillay, S.; San, E.J.; Msomi, N.; et al. Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa. medRxiv 2020. [Google Scholar] [CrossRef]
  13. Campbell, F.; Archer, B.; Laurenson-Schafer, H.; Jinnai, Y.; Konings, F.; Batra, N.; Pavlin, B.; Vandemaele, K.; Van Kerkhove, M.D.; Jombart, T.; et al. Increased transmissibility and global spread of SARS-CoV-2 variants of concern as at June 2021. Euro Surveill 2021, 26, 2100509. [Google Scholar] [CrossRef]
  14. Beucher, G.; Blondot, M.-L.; Celle, A.; Pied, N.; Recordon-Pinson, P.; Esteves, P.; Faure, M.; Métifiot, M.; Lacomme, S.; Dacheaux, D.; et al. SARS-CoV-2 transmission via apical syncytia release from primary bronchial epithelia and infectivity restriction in children epithelia. bioRxiv 2021. [Google Scholar] [CrossRef]
  15. Jary, A.; Leducq, V.; Malet, I.; Marot, S.; Klement-Frutos, E.; Teyssou, E.; Soulie, C.; Abdi, B.; Wirden, M.; Pourcher, V.; et al. Evolution of viral quasispecies during SARS-CoV-2 infection. Clin. Microbiol Infect. 2020, 26, 1560.e1–1560.e4. [Google Scholar] [CrossRef]
Figure 1. (A). Nucleotide diversity around the spike coding region. Screenshot of the Nextstrain.org website taking into account global GISAID data from 19 December 2020 to 8 April 2021 (top). The bottom panel represents the region amplified following RT-PCR (grey box encompassed by black arrows). (B). Typical quality control gel obtained after nested PCR amplification, M1 and M2 are DNA ladders, 1 to 4 are amplicons from four patients samples and t1-t2 are negative controls of the two steps amplification.
Figure 1. (A). Nucleotide diversity around the spike coding region. Screenshot of the Nextstrain.org website taking into account global GISAID data from 19 December 2020 to 8 April 2021 (top). The bottom panel represents the region amplified following RT-PCR (grey box encompassed by black arrows). (B). Typical quality control gel obtained after nested PCR amplification, M1 and M2 are DNA ladders, 1 to 4 are amplicons from four patients samples and t1-t2 are negative controls of the two steps amplification.
Covid 01 00028 g001
Figure 2. Smartgene IDNS software Coronavirus Module. This is an example of the spike region covered by the 10 amplifications.
Figure 2. Smartgene IDNS software Coronavirus Module. This is an example of the spike region covered by the 10 amplifications.
Covid 01 00028 g002
Table 1. List of reported mutations for the three VOCs and the 20A.EU2 (corresponding to SARS CoV2 clade which dominated the European second outbreak in 2020) compared to the reference Wuhan strain (WIV04/2019).
Table 1. List of reported mutations for the three VOCs and the 20A.EU2 (corresponding to SARS CoV2 clade which dominated the European second outbreak in 2020) compared to the reference Wuhan strain (WIV04/2019).
StrainWIV04/201920A.EU2AlphaBetaGammaDelta
AA position on spike\#N/A210121212
18L FF
19T R
20T N
26P S
69H -
70V -
80D A
138D Y
142G D
144Y -
154E K
157F -
158R -
190R S
215D G
242L -
243A -
244L -
246R I
417K NT
452L R
477SN
478T K
484E KKQ
501N YYY
570A D
614DGGGGG
655H Y
681P H R
701A V
716T I
950D N
982S A
1027T I
1071Q H
1118D H
1176V F
# = total number of mutations compared to the reference sequence. Deletions are indicated with a minus (−) sign. Signature mutations are highlighted in green, while mutations that could impact viral infectivity and/or immune response are in red.
Table 2. Primers used in this study.
Table 2. Primers used in this study.
Primer NameSequencePositions (Ref. NC_045512.2)Reference
RT-PCRForward F1-1GCATGGTGGACAGCCTTTGT21,220 to 21,239[7]
Reverse F1-1TGGGTTTTTGGAACGGCATT25,787 to 25,806
Nested PCR
+ sequencing
Forward F1-2GGGGTACTGCTGTTATGTCTTT21,422 to 21,443
Reverse F1-2CACCCTTGGAGAGTGCTAGT25,602 to 25,621
SequencingnCoV-2019_76_LEFTAGGGCAAACTGGAAAGATTGCT22,798 to 22,819ANRS AC43
nCoV-2019_76_RIGHTACCTGTGCCTGTTAAACCATTGA23,190 to 23,212
CV_14_FCAGGCTGCGTTATAGCTTGGAA22,851 to 22,872CDC
CV_14_RAACCAGTGTGTGCCATTTGAA24,870 to 24,850
CV_15_FGCAGAAATCAGAGCTTCTGCTAATC24,608 to 24,632
SARS2_B6FCCAGCAACTGTTTGTGGACCTA23,123 to 23,144
SARS SMFTCAATCCATCATTGCCTACAC23,629 to 23,649This work
SARS S3RGAGCAAAGGTGGCAAAACAG24,138 to 24,157
Table 3. Clinical samples characteristics and spike sequence.
Table 3. Clinical samples characteristics and spike sequence.
Sample TypeSampling DateSexCt RTqPCRKnown Spike MutationsOthers Spike MutationsVariant
NRdRpE
Nasal swabFeb 2021M28.7928.79NDdel69-70, del144, N501Y, A570D, D614G, P681H, T716I, S982A, D1118HV6A, Q14Q/PAlpha
Nasal swabFeb 2021F31.4131.44NDdel69-70, del144, N501Y, A570D, D614G, P681R, T716I, S982A, D1118H Alpha
Nasal swabFeb 2021F17.6819.39NDdel69-70, del144, N501Y, A570D, D614G, P681H, T716I, S982A, D1118HV47V/I, A831V, T998T/P, T1006T/PAlpha
Nasal swabFeb 2021M20.2621.62NDdel69-70, del144, N501Y, A570D, D614G, P681H, T716I, S982A, D1118HV6AAlpha
Nasal swabFeb 2021F24.9125.50NDdel69-70, del144, N501Y, A570D, D614G, P681R, T716I, S982A, D1118H Alpha
Nasal swabFeb 2021F21.3720.90NDL18F, D80A, D215G, del242-244, K417N, E484K, N501Y, D614G, A701V Beta
Nasal swabFeb 2021F35.16>40NDdel69-70, del144, N501Y, A570D, D614G, P681H, T716I, S982A, D1118H Alpha
Nasal swabFeb 2021M26.9028.65NDL18F, L452R, N501Y, A653V, H655Y, D796Y, G1219VR346I, Q667H, K1191Nnew
Nasal swabFeb 2021F23.7231.50NDL18F, L452R, N501Y, A653V, H655Y, D796Y, G1219VR346I, Q667H, K1191Nnew
Nasal swabFeb 2021M35.7631.50NDdel69-70, del144, N501Y, A570D, D614G, P681H, T716I, S982A, D1118HI197I/V, A831V, D839D/G, T1238T/IAlpha
Nasal swabFeb 2021M15.1416.5813.53D614GA222VWT
Nasal swabFeb 2021M27.1929.5725.64S477N, D614G, G1219C 20A.EU2
BALApril 2021M17.9019.4316.89del69-70, del144, N501Y, A570D, D614G, P681R, T716I, S982A, D1118H Alpha
BALFeb 2021F15.2616.613.56S477N, D614GF2L, S459F, K1045E20A.EU2
BALFeb 2021M21.9522.1220.29D614GT29I, Q675HWT
Cell cultureMay 2020-NDNDNDno mutationno mutationWT
BAL: branchoalveolar lavage, ND: not determined. Mutations that could impact viral infectivity and/or immune response are in red.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Recordon-Pinson, P.; Blondot, M.-L.; Bellecave, P.; Lafon, M.-E.; Tumiotto, C.; Métifiot, M.; Andreola, M.-L. A Simple and Fast Method to Sequence the Full-Length Spike Gene for SARS-CoV-2 Variant Identification from Patient Samples. COVID 2021, 1, 337-344. https://0-doi-org.brum.beds.ac.uk/10.3390/covid1010028

AMA Style

Recordon-Pinson P, Blondot M-L, Bellecave P, Lafon M-E, Tumiotto C, Métifiot M, Andreola M-L. A Simple and Fast Method to Sequence the Full-Length Spike Gene for SARS-CoV-2 Variant Identification from Patient Samples. COVID. 2021; 1(1):337-344. https://0-doi-org.brum.beds.ac.uk/10.3390/covid1010028

Chicago/Turabian Style

Recordon-Pinson, Patricia, Marie-Lise Blondot, Pantxika Bellecave, Marie-Edith Lafon, Camille Tumiotto, Mathieu Métifiot, and Marie-Line Andreola. 2021. "A Simple and Fast Method to Sequence the Full-Length Spike Gene for SARS-CoV-2 Variant Identification from Patient Samples" COVID 1, no. 1: 337-344. https://0-doi-org.brum.beds.ac.uk/10.3390/covid1010028

Article Metrics

Back to TopTop