Genomic Surveillance for SARS-CoV-2 Variants: Circulation of Omicron Lineages - United States, January 2022-May 2023.

CDC has used national genomic surveillance since December 2020 to monitor SARS-CoV-2 variants that have emerged throughout the COVID-19 pandemic, including the Omicron variant. This report summarizes U.S. trends in variant proportions from national genomic surveillance during January 2022-May 2023. During this period, the Omicron variant remained predominant, with various descendant lineages reaching national predominance (>50% prevalence). During the first half of 2022, BA.1.1 reached predominance by the week ending January 8, 2022, followed by BA.2 (March 26), BA.2.12.1 (May 14), and BA.5 (July 2); the predominance of each variant coincided with surges in COVID-19 cases. The latter half of 2022 was characterized by the circulation of sublineages of BA.2, BA.4, and BA.5 (e.g., BQ.1 and BQ.1.1), some of which independently acquired similar spike protein substitutions associated with immune evasion. By the end of January 2023, XBB.1.5 became predominant. As of May 13, 2023, the most common circulating lineages were XBB.1.5 (61.5%), XBB.1.9.1 (10.0%), and XBB.1.16 (9.4%); XBB.1.16 and XBB.1.16.1 (2.4%), containing the K478R substitution, and XBB.2.3 (3.2%), containing the P521S substitution, had the fastest doubling times at that point. Analytic methods for estimating variant proportions have been updated as the availability of sequencing specimens has declined. The continued evolution of Omicron lineages highlights the importance of genomic surveillance to monitor emerging variants and help guide vaccine development and use of therapeutics.

Travel history among persons infected with SARS-CoV-2 variants of concern in the United States, December 2020-February 2021.

The first three SARS-CoV-2 phylogenetic lineages classified as variants of concern (VOCs) in the United States (U.S.) from December 15, 2020 to February 28, 2021, Alpha (B.1.1.7), Beta (B.1.351), and Gamma (P.1) lineages, were initially detected internationally. This investigation examined available travel history of coronavirus disease 2019 (COVID-19) cases reported in the U.S. in whom laboratory testing showed one of these initial VOCs. Travel history, demographics, and health outcomes for a convenience sample of persons infected with a SARS-CoV-2 VOC from December 15, 2020 through February 28, 2021 were provided by 35 state and city health departments, and proportion reporting travel was calculated. Of 1,761 confirmed VOC cases analyzed, 1,368 had available data on travel history. Of those with data on travel history, 1,168 (85%) reported no travel preceding laboratory confirmation of SARS-CoV-2 and only 105 (8%) reported international travel during the 30 days preceding a positive SARS-CoV-2 test or symptom onset. International travel was reported by 92/1,304 (7%) of persons infected with the Alpha variant, 7/55 (22%) with Beta, and 5/9 (56%) with Gamma. Of the first three SARS-CoV-2 lineages designated as VOCs in the U.S., international travel was common only among the few Gamma cases. Most persons infected with Alpha and Beta variant reported no travel history, therefore, community transmission of these VOCs was likely common in the U.S. by March 2021. These findings underscore the importance of global surveillance using whole genome sequencing to detect and inform mitigation strategies for emerging SARS-CoV-2 VOCs.

Genome Sequences of Hemolytic and Nonhemolytic Listeria innocua Strains from Human, Food, and Environmental Sources.

This report describes genome sequences for nine Listeria innocua strains that varied in hemolytic phenotypes on sheep blood agar. All strains were sequenced using Pacific Biosciences (PacBio) single-molecule real-time (SMRT) chemistry; overall, the average read length of these sequences was 2,869,880 bp, with an average GC content of 37%.

The first complete genome of the simian malaria parasite Plasmodium brasilianum.

Naturally occurring human infections by zoonotic Plasmodium species have been documented for P. knowlesi, P. cynomolgi, P. simium, P. simiovale, P. inui, P. inui-like, P. coatneyi, and P. brasilianum. Accurate detection of each species is complicated by their morphological similarities with other Plasmodium species. PCR-based assays offer a solution but require prior knowledge of adequate genomic targets that can distinguish the species. While whole genomes have been published for P. knowlesi, P. cynomolgi, P. simium, and P. inui, no complete genome for P. brasilianum has been available. Previously, we reported a draft genome for P. brasilianum, and here we report the completed genome for P. brasilianum. The genome is 31.4 Mb in size and comprises 14 chromosomes, the mitochondrial genome, the apicoplast genome, and 29 unplaced contigs. The chromosomes consist of 98.4% nucleotide sites that are identical to the P. malariae genome, the closest evolutionarily related species hypothesized to be the same species as P. brasilianum, with 41,125 non-synonymous SNPs (0.0722% of genome) identified between the two genomes. Furthermore, P. brasilianum had 4864 (82.1%) genes that share 80% or higher sequence similarity with 4970 (75.5%) P. malariae genes. This was demonstrated by the nearly identical genomic organization and multiple sequence alignments for the merozoite surface proteins msp3 and msp7. We observed a distinction in the repeat lengths of the circumsporozoite protein (CSP) gene sequences between P. brasilianum and P. malariae. Our results demonstrate a 97.3% pairwise identity between the P. brasilianum and the P. malariae genomes. These findings highlight the phylogenetic proximity of these two species, suggesting that P. malariae and P. brasilianum are strains of the same species, but this could not be fully evaluated with only a single genomic sequence for each species.

Genomics and metagenomics of Madurella mycetomatis, a causative agent of black grain mycetoma in Sudan.

Madurella mycetomatis is one of the main causative agents of mycetoma, a debilitating neglected tropical disease. Improved understanding of the genomic diversity of the fungal and bacterial causes of mycetoma is essential to advances in diagnosis and treatment. Here, we describe a high-quality genome assembly of M. mycetomatis and results of the whole genome sequence analysis of 26 isolates from Sudan. We demonstrate evidence of at least seven genetically diverse lineages and extreme clonality among isolates within these lineages. We also performed shotgun metagenomic analysis of DNA extracted from mycetoma grains and showed that M. mycetomatis reads were detected in all sequenced samples with the average of 11,317 reads (s.d. +/- 21,269) per sample. In addition, 10 (12%) of the 81 tested grain samples contained bacterial reads including Streptococcus sp., Staphylococcus sp. and others.

Multiple lineages of monkeypox virus detected in the United States, 2021-2022.

Monkeypox is a viral zoonotic disease endemic in Central and West Africa. In May 2022, dozens of non-endemic countries reported hundreds of monkeypox cases, most with no epidemiological link to Africa. We identified two lineages of monkeypox virus (MPXV) among two 2021 and seven 2022 US monkeypox cases: the major 2022 outbreak variant called B.1 and a minor contemporaneously sampled variant called A.2. Analyses of mutations among these two variants revealed an extreme preference for GA-to-AA mutations indicative of human APOBEC3 cytosine deaminase activity among Clade IIb MPXV (previously West African, Nigeria) sampled since 2017. Such mutations were not enriched within other MPXV clades. These findings suggest that APOBEC3 editing may be a recurrent and a dominant driver of MPXV evolution within the current outbreak.

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.

Genomic Surveillance for SARS-CoV-2 Variants: Predominance of the Delta (B.1.617.2) and Omicron (B.1.1.529) Variants - United States, June 2021-January 2022.

Genomic surveillance is a critical tool for tracking emerging variants of SARS-CoV-2 (the virus that causes COVID-19), which can exhibit characteristics that potentially affect public health and clinical interventions, including increased transmissibility, illness severity, and capacity for immune escape. During June 2021-January 2022, CDC expanded genomic surveillance data sources to incorporate sequence data from public repositories to produce weighted estimates of variant proportions at the jurisdiction level and refined analytic methods to enhance the timeliness and accuracy of national and regional variant proportion estimates. These changes also allowed for more comprehensive variant proportion estimation at the jurisdictional level (i.e., U.S. state, district, territory, and freely associated state). The data in this report are a summary of findings of recent proportions of circulating variants that are updated weekly on CDC's COVID Data Tracker website to enable timely public health action.† The SARS-CoV-2 Delta (B.1.617.2 and AY sublineages) variant rose from 1% to >50% of viral lineages circulating nationally during 8 weeks, from May 1-June 26, 2021. Delta-associated infections remained predominant until being rapidly overtaken by infections associated with the Omicron (B.1.1.529 and BA sublineages) variant in December 2021, when Omicron increased from 1% to >50% of circulating viral lineages during a 2-week period. As of the week ending January 22, 2022, Omicron was estimated to account for 99.2% (95% CI = 99.0%-99.5%) of SARS-CoV-2 infections nationwide, and Delta for 0.7% (95% CI = 0.5%-1.0%). The dynamic landscape of SARS-CoV-2 variants in 2021, including Delta- and Omicron-driven resurgences of SARS-CoV-2 transmission across the United States, underscores the importance of robust genomic surveillance efforts to inform public health planning and practice.

Analysis of the initial lot of the CDC 2019-Novel Coronavirus (2019-nCoV) real-time RT-PCR diagnostic panel.

At the start of the COVID-19 pandemic, the Centers for Disease Control and Prevention (CDC) designed, manufactured, and distributed the CDC 2019-Novel Coronavirus (2019-nCoV) Real-Time RT-PCR Diagnostic Panel for SARS-CoV-2 detection. The diagnostic panel targeted three viral nucleocapsid gene loci (N1, N2, and N3 primers and probes) to maximize sensitivity and to provide redundancy for virus detection if mutations occurred. After the first distribution of the diagnostic panel, state public health laboratories reported fluorescent signal in the absence of viral template (false-positive reactivity) for the N3 component and to a lesser extent for N1. This report describes the findings of an internal investigation conducted by the CDC to identify the cause(s) of the N1 and N3 false-positive reactivity. For N1, results demonstrate that contamination with a synthetic template, that occurred while the "bulk" manufactured materials were located in a research lab for quality assessment, was the cause of false reactivity in the first lot. Base pairing between the 3' end of the N3 probe and the 3' end of the N3 reverse primer led to amplification of duplex and larger molecules resulting in false reactivity in the N3 assay component. We conclude that flaws in both assay design and handling of the "bulk" material, caused the problems with the first lot of the 2019-nCoV Real-Time RT-PCR Diagnostic Panel. In addition, within this study, we found that the age of the examined diagnostic panel reagents increases the frequency of false positive results for N3. We discuss these findings in the context of improvements to quality control, quality assurance, and assay validation practices that have since been improved at the CDC.

Complete Genome Sequence of Rickettsia parkeri Strain Black Gap.

A unique genotype of Rickettsia parkeri, designated R. parkeri strain Black Gap, has thus far been associated exclusively with the North American tick, Dermacentor parumapertus. The compete genome consists of a single circular chromosome with 1,329,522 bp and a G+C content of 32.5%.

Complete Genome Sequence of Rhodococcus sp. Strain W8901, a Human Clinical Specimen, Assembled Using MiSeq and MinION Sequence Data.

Rhodococcus sp. strain W8901 is a Gram-positive, aerobic, mycolic acid-containing coccobacillus obtained from a patient with acute lymphocytic leukemia. Here, we report on the complete, circular genome sequence obtained using Illumina MiSeq and Oxford Nanopore Technologies MinION reads in order to better resolve the phylogeny of a rare pathogen.

Characterization of a novel transitional group Rickettsia species (Rickettsia tillamookensis sp. nov.) from the western black-legged tick, Ixodes pacificus.

A previously unrecognized Rickettsia species was isolated in 1976 from a pool of Ixodes pacificus ticks collected in 1967 from Tillamook County, Oregon, USA. The isolate produced low fever and mild scrotal oedema following intraperitoneal injection into male guinea pigs (Cavia porcellus). Subsequent serotyping characterized this isolate as distinct from recognized typhus and spotted fever group Rickettsia species; nonetheless, the isolate remained unevaluated by molecular techniques and was not identified to species level for the subsequent 30 years. Ixodes pacificus is the most frequently identified human-biting tick in the western United States, and as such, formal identification and characterization of this potentially pathogenic Rickettsia species is warranted. Whole-genome sequencing of the Tillamook isolate revealed a genome 1.43 Mbp in size with 32.4 mol% G+C content. Maximum-likelihood phylogeny of core proteins places it in the transitional group of Rickettsia basal to both Rickettsia felis and Rickettsia asembonensis. It is distinct from existing named species, with maximum average nucleotide identity of 95.1% to R. asembonensis and maximum digital DNA-DNA hybridization score similarity to R. felis at 80.1%. The closest similarity at the 16S rRNA gene (97.9%) and sca4 (97.5%/97.6% respectively) is to Candidatus 'Rickettsia senegalensis' and Rickettsia sp. cf9, both isolated from cat fleas (Ctenocephalides felis). We characterized growth at various temperatures and in multiple cell lines. The Tillamook isolate grows aerobically in Vero E6, RF/6A and DH82 cells, and growth is rapid at 28 °C and 32 °C. Using accepted genomic criteria, we propose the name Rickettsia tillamookensis sp. nov., with the type strain Tillamook 23. Strain Tillamook 23 is available from the Centers for Disease Control and Prevention Rickettsial Isolate Reference Collection (WDCM 1093), Atlanta, GA, USA (CRIRC accession number RTI001T) and the Collection de Souches de l'Unité des Rickettsies (WDCM 875), Marseille, France (CSUR accession number R5043). Using accepted genomic criteria, we propose the name Rickettsia tillamookensis sp. nov., with the type strain Tillamook 23 (=CRIRC RTI001=R5043).

Genomic Surveillance for SARS-CoV-2 Variants Circulating in the United States, December 2020-May 2021.

SARS-CoV-2, the virus that causes COVID-19, is constantly mutating, leading to new variants (1). Variants have the potential to affect transmission, disease severity, diagnostics, therapeutics, and natural and vaccine-induced immunity. In November 2020, CDC established national surveillance for SARS-CoV-2 variants using genomic sequencing. As of May 6, 2021, sequences from 177,044 SARS-CoV-2-positive specimens collected during December 20, 2020-May 6, 2021, from 55 U.S. jurisdictions had been generated by or reported to CDC. These included 3,275 sequences for the 2-week period ending January 2, 2021, compared with 25,000 sequences for the 2-week period ending April 24, 2021 (0.1% and 3.1% of reported positive SARS-CoV-2 tests, respectively). Because sequences might be generated by multiple laboratories and sequence availability varies both geographically and over time, CDC developed statistical weighting and variance estimation methods to generate population-based estimates of the proportions of identified variants among SARS-CoV-2 infections circulating nationwide and in each of the 10 U.S. Department of Health and Human Services (HHS) geographic regions.* During the 2-week period ending April 24, 2021, the B.1.1.7 and P.1 variants represented an estimated 66.0% and 5.0% of U.S. SARS-CoV-2 infections, respectively, demonstrating the rise to predominance of the B.1.1.7 variant of concern† (VOC) and emergence of the P.1 VOC in the United States. Using SARS-CoV-2 genomic surveillance methods to analyze surveillance data produces timely population-based estimates of the proportions of variants circulating nationally and regionally. Surveillance findings demonstrate the potential for new variants to emerge and become predominant, and the importance of robust genomic surveillance. Along with efforts to characterize the clinical and public health impact of SARS-CoV-2 variants, surveillance can help guide interventions to control the COVID-19 pandemic in the United States.

Chromosome-Level Genome Sequence of Leishmania (Leishmania) tropica Strain CDC216-162, Isolated from an Afghanistan Clinical Case.

PacBio and Illumina MiSeq platforms were used for genomic sequencing of a Leishmania (Leishmania) tropica strain isolated from a patient infected in Pakistan. PacBio assemblies were generated using Flye v2.4 and polished with MiSeq data. The results represent a considerable improvement of the currently available genome sequences in the GenBank database.

Prevalent, protective, and convergent IgG recognition of SARS-CoV-2 non-RBD spike epitopes.

The molecular composition and binding epitopes of the immunoglobulin G (IgG) antibodies that circulate in blood plasma after severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection are unknown. Proteomic deconvolution of the IgG repertoire to the spike glycoprotein in convalescent subjects revealed that the response is directed predominantly (>80%) against epitopes residing outside the receptor binding domain (RBD). In one subject, just four IgG lineages accounted for 93.5% of the response, including an amino (N)-terminal domain (NTD)-directed antibody that was protective against lethal viral challenge. Genetic, structural, and functional characterization of a multidonor class of "public" antibodies revealed an NTD epitope that is recurrently mutated among emerging SARS-CoV-2 variants of concern. These data show that "public" NTD-directed and other non-RBD plasma antibodies are prevalent and have implications for SARS-CoV-2 protection and antibody escape.

Cutaneous Leishmaniasis Caused by an Unknown Leishmania Strain, Arizona, USA.

We investigated an autochthonous case of cutaneous leishmaniasis caused by a genetically different Leishmania sp. in a patient in Arizona, USA. This parasite was classified into the subgenus Leishmania on the basis of multilocus DNA sequence and phylogenetic analyses of the rRNA locus and 11 reference genes.

Draft Chromosome Sequences of a Clinical Isolate of the Free-Living Ameba Naegleria fowleri.

We present the chromosome sequences of a Naegleria fowleri isolate from a human primary amebic meningoencephalitis (PAM) case. The genome sequences were assembled from Illumina HiSeq and PacBio sequencing data and verified with the optical mapping data. This led to the identification of 37 contigs representing 37 chromosomes in N. fowleri.

Clade-specific chromosomal rearrangements and loss of subtelomeric adhesins in Candida auris.

Candida auris is an emerging fungal pathogen of rising concern due to global spread, the ability to cause healthcare-associated outbreaks, and antifungal resistance. Genomic analyses revealed that early contemporaneously detected cases of C. auris were geographically stratified into four major clades. While Clades I, III, and IV are responsible for ongoing outbreaks of invasive and multidrug-resistant infections, Clade II, also termed the East Asian clade, consists primarily of cases of ear infection, is often susceptible to all antifungal drugs, and has not been associated with outbreaks. Here, we generate chromosome-level assemblies of twelve isolates representing the phylogenetic breadth of these four clades and the only isolate described to date from Clade V. This Clade V genome is highly syntenic with those of Clades I, III, and IV, although the sequence is highly divergent from the other clades. Clade II genomes appear highly rearranged, with translocations occurring near GC-poor regions, and large subtelomeric deletions in most chromosomes, resulting in a substantially different karyotype. Rearrangements and deletion lengths vary across Clade II isolates, including two from a single patient, supporting ongoing genome instability. Deleted subtelomeric regions are enriched in Hyr/Iff-like cell-surface proteins, novel candidate cell wall proteins, and an ALS-like adhesin. Cell wall proteins from these families and other drug-related genes show clade-specific signatures of selection in Clades I, III, and IV. Subtelomeric dynamics and the conservation of cell surface proteins in the clades responsible for global outbreaks causing invasive infections suggest an explanation for the different phenotypes observed between clades.

Complete and Circularized Bacterial Genome Sequence of Gordonia sp. Strain X0973.

Gordonia sp. strain X0973 is a Gram-positive, weakly acid-fast, aerobic actinomycete obtained from a human abscess with Gordonia araii NBRC 100433T as its closest phylogenetic neighbor. Here, we report using Illumina MiSeq and PacBio reads to assemble the complete and circular genome sequence of 3.75 Mbp with 3,601 predicted coding sequences.