Volume 28, Number 7—July 2022
SARS-CoV-2 Delta–Omicron Recombinant Viruses, United States
To detect new and changing SARS-CoV-2 variants, we investigated candidate Delta–Omicron recombinant genomes from Centers for Disease Control and Prevention national genomic surveillance. Laboratory and bioinformatic investigations identified and validated 9 genetically related SARS-CoV-2 viruses with a hybrid Delta–Omicron spike protein.
Emerging variants of SARS-CoV-2 are characterized and monitored closely by national genomic surveillance. In addition to sequencing efforts from US public health, academic, and commercial laboratories, the Centers for Disease Control and Prevention (CDC) collects and sequences SARS-CoV-2 specimens from 64 partners across state, tribal, local, and territorial public health agencies through the National SARS-CoV-2 Strain Surveillance program (https://www.cdc.gov/coronavirus/2019-ncov/variants/cdc-role-surveillance.html) and funds SARS-CoV-2 sequencing through a nationwide network of commercial laboratory testing companies. To date, these efforts have contributed 1.8 million SARS-CoV-2 genomes from the United States to public repositories. The purpose of this genomic surveillance system is to detect and respond dynamically to new and changing SARS-CoV-2 variants (1).
Recombination is an evolutionary mechanism frequently observed in coronaviruses (2,3), and it can lead to rapid accumulation of mutations and heightened transmissibility (4). SARS-CoV-2 recombination events have also been found to arise disproportionately in the spike gene (Y. Turkahia et al., unpub. data, https://www.biorxiv.org/content/10.1101/2021.08.04.455157V1). Recombination between Alpha and Delta SARS-CoV-2 variants has been documented (5–7).
Given the divergence of the Delta and Omicron variant genomes, as well as the known immune-escape properties of Omicron (8,9), a Delta–Omicron recombinant strain could alter the landscape of vaccine and therapeutic effectiveness. In early 2022, viruses resulting from recombination between Delta and Omicron were reported, but further inspection indicated that these claims seemed to have resulted from laboratory artifact or co-infections (10). With this study, we identified candidate Delta–Omicron recombinant genomes from the CDC national genomic surveillance and attempted to rule out laboratory contamination or sequencing error.
We identified 9 candidate recombinant sequences (Table) from CDC national genomic surveillance dataset made publicly available in GenBank and GISAID EpiCoV (https://www.gisaid.org). Using Bolotie, a rapid interclade recombination detection method (3), we identified these sequences as candidate recombinant genomes, having 1 parent in Delta (clade 21J) and 1 in Omicron (clade 21K). Bolotie describes a single breakpoint between nucleotide positions 22035 and 22577 (referenced to GenBank accession no. NC_045512.2); there are no differentiating mutations between clades 21J and 21K within this range. These sequences (EPI_ISL_8720194, EPI_ISL_9147438, EPI_ISL_9147935, EPI_ISL_8981459, EPI_ISL_8981824, EPI_ISL_9088187 [A. Bolze et al., unpub .data, https://www.medrxiv.org/content/medrxiv/early/2022/03/12/2022.03.09.22272113.full.pdf], EPI_ISL_8981712, EPI_ISL_10389339, EPI_ISL_10389336) contain hallmark mutation sets from both Omicron and Delta SARS-CoV-2 lineages, changing from Delta-associated substitutions to Omicron-associated substitutions between spike protein amino acids 158 and 339 (Appendix Figure 1, panel A). This breakpoint is distinct from the 2 clusters of apparent Delta–Omicron recombinants identified in the United Kingdom (https://github.com/cov-lineages/pango-designation/issues/422 and https://github.com/cov-lineages/pango-designation/issues/441), which have a breakpoint upstream of spike in the ORF1ab gene (Appendix Figure 1, panel A), and these samples show a singular breakpoint, unlike concurrently observed Delta–Omicron recombinants in France (P. Colson et al., unpub. data, https://www.medrxiv.org/content/10.1101/2022.03.03.22271812V1).
To rule out Delta and Omicron co-infection, laboratory contamination, and bioinformatic error, we examined the raw read data from the 9 candidate recombinants created from molecular loop and amplicon-based sequencing strategies. Two of these specimens were readily available from the original diagnostic laboratory, and extracted RNA was shipped to CDC for confirmatory sequencing. We used Illumina (https://www.illumina.com) and PacBio (https://www.pacb.com) sequencing of 2 whole-genome amplicon strategies, as well as spike-gene amplification followed by Oxford Nanopore (https://nanoporetech.com) sequencing (Appendix). All sequencing strategies yielded functionally identical consensus sequences compared with the corresponding original sequencing strategies.
Nextclade (11) classified the 9 whole genomes as 21K (Omicron/BA.1). We then split the genomes at position 22150 (within the predicted recombination site range). Nextclade classified the first 22150 base fragment as clade 21J (Delta) and the remainder as clade 21K (Omicron/BA.1). Pangolin version 3.1.20 (pangoLEARN 1.2.123, Scorpio 0.3.16, https://cov-lineages.org) assigned a lineage of none to the full-genome sequences. Pangolin classified the first 22150 base fragment of each recombinant as AY.43 (Delta), although the call was not supported by Scorpio. Inspection of this region revealed closer homology to AY.119.2 (Delta) sequences because of mutations orf1ab:A2855V and orf1ab:A6248S, which are common to AY.119 lineages, and orf1ab:K4176N, which is found in a subset of AY.119.2 (Delta) sequences. The remaining sequence fragment from nt 22151 to the 3′ end was classified by pangolin as BA.1.1 (Omicron). This observation has been documented in the PANGO-designations repository (https://github.com/cov-lineages/pango-designation/issues/439) and is under review for potential lineage assignment.
Detailed sequence analysis confirmed the 2 resequenced specimens as true recombinants and indicated no evidence of co-infection or contamination. Comparison with a representative AY.119.2 (Delta) specimen indicated characteristic Delta mutations (C21618G, C21846T, G21987A, and deletion 22029–22034) at >99% frequency (>600× coverage for Oxford Nanopore, >1,800× coverage for PacBio, >1000× coverage for Illumina) in the 5′ end of the recombinant (Appendix Figure 1, panel B). The 2 BA.1.1 (Omicron) deletions at the beginning of the spike gene (21765–21770 and 21987–21995) and the characteristic Omicron 9-base insertion after nt 22205 were not present in read data, consistent with a Delta origin for the 5′ end of the spike gene. After position 22577, the mutation profiles mirrored that of a representative BA.1.1 (Omicron) specimen (Appendix Figure 1, panel B). Analysis of individual Oxford Nanopore reads showed characteristic Delta mutations co-occurring with Omicron single-nucleotide variants on the same reads (sharing Delta 22029–22034 deletion and Omicron 22673 T>C; Appendix Figure 2). The translated spike protein is a hybrid, containing characteristic amino acids from both Delta and Omicron parents with a breakpoint between the N terminal domain and receptor-binding domain of spike S1 protein (Appendix Figure 1, panel A).
To visualize the parents of the recombinant genomes, we split all candidate recombinant genomes at position 22150, within the predicted breakpoint, and used Nextclade (11) to place each genome fragment (1–22150 and 22151 through the 3′ end) onto a reference tree. We visualized the 2 trees as a tanglegram tree with Auspice (12). Nucleotides 1–22150 clustered with clade 21J (Delta) sequences, and the remaining fragment of the genome clustered with 21K (Omicron/BA.1) (Appendix Figure 3).
Our results provide evidence of a recombinant SARS-CoV-2 genome containing a hybrid spike protein derived from a Delta (AY.119.2)–Omicron (BA.1.1) recombination event. However, the ability to effectively identify and confirm additional recombinant viruses remains challenging because of the range of sequence quality available in the public domain. These limitations are a result of amplification inefficiency and consensus-calling algorithmic error, as well as cases of co-infection or potential sample contamination.
Comparative phenotypic characterization of virus isolates from the recombinant cluster was not possible because all specimens were chemically inactivated. In the spike protein, there are no additional amino acid substitutions within the receptor-binding domain compared with BA.1.1 (Omicron) lineage viruses. Recombinant viruses with this hybrid spike protein were detected over the course of 6 weeks, but the number of cases resulting from those viruses remains low. Most cases were identified within the mid-Atlantic region of the United States. However, epidemiologic linkage cannot be determined because CDC does not collect identifying information for these samples.
Systematic virus surveillance is essential for long-term monitoring of SARS-CoV-2 evolution. Given the potential public health consequences of new variants emerging from recombination, investigations involving laboratory and bioinformatic components, such as the one presented here, are critical for correctly identifying and tracking these viruses.
Ms. Lacek is an informatics health scientist with the Influenza Division, National Center for Immunization and Respiratory Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia. Her primary research interests include viral genomics, next-generation sequencing, and genome assembly.
We thank the public health program and laboratory staff members who contribute to National SARS-CoV-2 Strain Surveillance, including the Association of Public Health Laboratories, commercial laboratory staff members, laboratories submitting tagged baseline SARS-CoV-2 sequences via GISAID EpiCov and GenBank, CDC COVID-19 Pandemic Response Laboratory Task Force, Strain Surveillance and Emerging Variants Bioinformatics Working Group, Natalie Thornburg, and colleagues in the UK Health Security Agency. Furthermore, we gratefully acknowledge the authors from the originating laboratories responsible for obtaining the specimens and the submitting laboratories where genetic sequence data were generated and shared via GISAID. A full acknowledgment table can be found with the EPI_SET ID EPI_SET_20220412tv, in the Data Acknowledgement Locator under GISAID resources (https://www.gisaid.org/epi-search).
All authors have completed and submitted the International Committee of Medical Journal Editors form for disclosure of potential conflicts of interest. No potential conflicts of interest were disclosed.
- National Center for Immunization and Respiratory Diseases (NCIRD) DoVD. Science brief: emerging SARS-CoV-2 variants [cited 2022 Feb 17]. https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/scientific-brief-emerging-variants.html
- Lai MM, Baric RS, Makino S, Keck JG, Egbert J, Leibowitz JL, et al. Recombination between nonsegmented RNA genomes of murine coronaviruses. J Virol. 1985;56:449–56. DOIPubMedGoogle Scholar
- Varabyou A, Pockrandt C, Salzberg SL, Pertea M. Rapid detection of inter-clade recombination in SARS-CoV-2 with Bolotie. Genetics. 2021;218:iyab074.
- Li X, Giorgi EE, Marichannegowda MH, Foley B, Xiao C, Kong X-P, et al. Emergence of SARS-CoV-2 through recombination and strong purifying selection. Sci Adv. 2020;6:
eabb9153. DOIPubMedGoogle Scholar
- Sekizuka T, Itokawa K, Saito M, Shimatani M, Matsuyama S, Hasegawa H, et al. Genome recombination between Delta and Alpha variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Jpn J Infect Dis. 2022;
JJID.2021.844; Epub ahead of print. DOIPubMedGoogle Scholar
- Escalera-Zamudio M, Sánchez HGC, Arredondo LJD, Gutiérrez BRE. Proposal to redesignate B.1.631 as recombinant lineage XB [cited 2022 Feb 18]. https://virological.org/t/re-proposal-to-redesignate-b-1-631-as-recombinant-lineage-xb/746
- Jackson B, Rambaut A, Pybus OG, Robertson DL, Connor T, Loman NJ, et al. Recombinant SARS-CoV-2 genomes involving lineage B.1.1.7 in the UK [cited 2022 Feb 18]. https://virological.org/t/recombinant-sars-cov-2-genomes-involving-lineage-b-1-1-7-in-the-uk/658.
- Liu L, Iketani S, Guo Y, Chan JF, Wang M, Liu L, et al. Striking antibody evasion manifested by the Omicron variant of SARS-CoV-2. Nature. 2022;602:676–81. DOIPubMedGoogle Scholar
- Chakraborty C, Bhattacharya M, Sharma AR. Emerging mutations in the SARS-CoV-2 variants and their role in antibody escape to small molecule-based therapeutic resistance. Curr Opin Pharmacol. 2022;62:64–73. DOIPubMedGoogle Scholar
- Kreier F. Deltacron: the story of the variant that wasn’t. Nature. 2022;602:19. DOIPubMedGoogle Scholar
- Aksamentov I, Roemer C, Hodcroft EB, Neher RA. Nextclade: clade assignment, mutation calling and quality control for viral genomes. J Open Source Softw. 2021;6:3773. DOIGoogle Scholar
- Hadfield J, Megill C, Bell SM, Huddleston J, Potter B, Callender C, et al. Nextstrain: real-time tracking of pathogen evolution. Bioinformatics. 2018;34:4121–3. DOIPubMedGoogle Scholar
TableCite This Article
Original Publication Date: May 12, 2022
Table of Contents – Volume 28, Number 7—July 2022
|EID Search Options|
|Advanced Article Search – Search articles by author and/or keyword.|
|Articles by Country Search – Search articles by the topic country.|
|Article Type Search – Search articles by article type and issue.|
Please use the form below to submit correspondence to the authors or contact them at the following address:
Clinton R. Paden, Centers for Disease Control and Prevention, 1600 Clifton Rd NE, Mailstop H18-6, Atlanta, GA 30329-4027, USA