Complete genome sequences of Clostridioides difficile surveillance isolates representing the top 10 ribotypes from the Emerging Infections Program, United States, 2016.

Ten Clostridioides difficile isolates representing the top 10 ribotypes collected in 2016 through the Emerging Infections Program underwent long-read sequencing to obtain high-quality reference genome assemblies. These isolates are publicly available through the CDC & FDA Antibiotic Resistance Isolate Bank.

A novel invasive Streptococcus pyogenes variant sublineage derived through recombinational replacement of the emm12 genomic region.

Group A streptococcal strains potentially acquire new M protein gene types through genetic recombination (emm switching). To detect such variants, we screened 12,596 invasive GAS genomes for strains of differing emm types that shared the same multilocus sequence type (ST). Through this screening we detected a variant consisting of 16 serum opacity factor (SOF)-positive, emm pattern E, emm82 isolates that were ST36, previously only associated with SOF-negative, emm pattern A, emm12. The 16 emm82/ST36 isolates were closely interrelated (pairwise SNP distance of 0-43), and shared the same emm82-containing recombinational fragment. emm82/ST36 isolates carried the sof12 structural gene, however the sof12 indel characteristic of emm12 strains was corrected to confer the SOF-positive phenotype. Five independent emm82/ST36 invasive case isolates comprised two sets of genetically indistinguishable strains. The emm82/ST36 isolates were primarily macrolide resistant (12/16 isolates), displayed at least 4 different core genomic arrangements, and carried 11 different combinations of virulence and resistance determinants. Phylogenetic analysis revealed that emm82/ST36 was within a minor (non-clade 1) portion of ST36 that featured almost all ST36 antibiotic resistance. This work documents emergence of a rapidly diversifying variant that is the first confirmed example of an emm pattern A strain switched to a pattern E strain.

The Impact of SIV-Induced Immunodeficiency on Clinical Manifestation, Immune Response, and Viral Dynamics in SARS-CoV-2 Coinfection.

Persistent and uncontrolled SARS-CoV-2 replication in immunocompromised individuals has been observed and may be a contributing source of novel viral variants that continue to drive the pandemic. Importantly, the effects of immunodeficiency associated with chronic HIV infection on COVID-19 disease and viral persistence have not been directly addressed in a controlled setting. Here we conducted a pilot study wherein two pigtail macaques (PTM) chronically infected with SIVmac239 were exposed to SARS-CoV-2 and monitored for six weeks for clinical disease, viral replication, and viral evolution, and compared to our previously published cohort of SIV-naïve PTM infected with SARS-CoV-2. At the time of SARS-CoV-2 infection, one PTM had minimal to no detectable CD4+ T cells in gut, blood, or bronchoalveolar lavage (BAL), while the other PTM harbored a small population of CD4+ T cells in all compartments. Clinical signs were not observed in either PTM; however, the more immunocompromised PTM exhibited a progressive increase in pulmonary infiltrating monocytes throughout SARS-CoV-2 infection. Single-cell RNA sequencing (scRNAseq) of the infiltrating monocytes revealed a less activated/inert phenotype. Neither SIV-infected PTM mounted detectable anti-SARS-CoV-2 T cell responses in blood or BAL, nor anti-SARS-CoV-2 neutralizing antibodies. Interestingly, despite the diminished cellular and humoral immune responses, SARS-CoV-2 viral kinetics and evolution were indistinguishable from SIV-naïve PTM in all sampled mucosal sites (nasal, oral, and rectal), with clearance of virus by 3-4 weeks post infection. SIV-induced immunodeficiency significantly impacted immune responses to SARS-CoV-2 but did not alter disease progression, viral kinetics or evolution in the PTM model. SIV-induced immunodeficiency alone may not be sufficient to drive the emergence of novel viral variants.

Real-Time Analysis of SARS-CoV-2-Induced Cytolysis Reveals Distinct Variant-Specific Replication Profiles.

The ability of each new SARS-CoV-2 variant to evade host humoral immunity is the focus of intense research. Each variant may also harbor unique replication capabilities relevant for disease and transmission. Here, we demonstrate a new approach to assessing viral replication kinetics using real-time cell analysis (RTCA). Virus-induced cell death is measured in real time as changes in electrical impedance through cell monolayers while images are acquired at defined intervals via an onboard microscope and camera. Using this system, we quantified replication kinetics of five clinically important viral variants: WA1/2020 (ancestral), Delta, and Omicron subvariants BA.1, BA.4, and BA.5. Multiple measures proved useful in variant replication comparisons, including the elapsed time to, and the slope at, the maximum rate of cell death. Important findings include significantly weaker replication kinetics of BA.1 by all measures, while BA.5 harbored replication kinetics at or near ancestral levels, suggesting evolution to regain replicative capacity, and both an altered profile of cell killing and enhanced fusogenicity of the Delta variant. Together, these data show that RTCA is a robust method to assess replicative capacity of any given SARS-CoV-2 variant rapidly and quantitatively, which may be useful in assessment of newly emerging variants.

Reoccurring Escherichia coli O157:H7 Strain Linked to Leafy Greens-Associated Outbreaks, 2016-2019.

Genomic characterization of an Escherichia coli O157:H7 strain linked to leafy greens-associated outbreaks dates its emergence to late 2015. One clade has notable accessory genomic content and a previously described mutation putatively associated with increased arsenic tolerance. This strain is a reoccurring, emerging, or persistent strain causing illness over an extended period.

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.

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%.

Breakthrough Gastrointestinal COVID-19 and Intrahost Evolution Consequent to Combination Monoclonal Antibody Prophylaxis.

Breakthrough gastrointestinal COVID-19 was observed after experimental SARS-CoV-2 upper mucosal infection in a rhesus macaque undergoing low-dose monoclonal antibody prophylaxis. High levels of viral RNA were detected in intestinal sites contrasting with minimal viral replication in upper respiratory mucosa. Sequencing of virus recovered from tissue in 3 gastrointestinal sites and rectal swab revealed loss of furin cleavage site deletions present in the inoculating virus stock and 2 amino acid changes in spike that were detected in 2 colon sites but not elsewhere, suggesting compartmentalized replication and intestinal viral evolution. This suggests suboptimal antiviral therapies promote viral sequestration in these anatomies.

Intra-Host SARS-CoV-2 Evolution in the Gut of Mucosally-Infected Chlorocebus aethiops (African Green Monkeys).

In recent months, several SARS-CoV-2 variants have emerged that enhance transmissibility and escape host humoral immunity. Hence, the tracking of viral evolutionary trajectories is clearly of great importance. Little is known about SARS-CoV-2 evolution in nonhuman primate models used to test vaccines and therapies and to model human disease. Viral RNA was sequenced from rectal swabs from Chlorocebus aethiops (African green monkeys) after experimental respiratory SARS-CoV-2 infection. Two distinct patterns of viral evolution were identified that were shared between all collected samples. First, mutations in the furin cleavage site that were initially present in the virus as a consequence of VeroE6 cell culture adaptation were not detected in viral RNA recovered in rectal swabs, confirming the necessity of this motif for viral infection in vivo. Three amino acid changes were also identified; ORF 1a S2103F, and spike D215G and H655Y, which were detected in rectal swabs from all sampled animals. These findings are demonstrative of intra-host SARS-CoV-2 evolution and may identify a host-adapted variant of SARS-CoV-2 that would be useful in future primate models involving SARS-CoV-2 infection.

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%.

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).

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.

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.

Multispecies Outbreak of Verona Integron-Encoded Metallo-ß-Lactamase-Producing Multidrug Resistant Bacteria Driven by a Promiscuous Incompatibility Group A/C2 Plasmid.

BACKGROUND: Antibiotic resistance is often spread through bacterial populations via conjugative plasmids. However, plasmid transfer is not well recognized in clinical settings because of technical limitations, and health care-associated infections are usually caused by clonal transmission of a single pathogen. In 2015, multiple species of carbapenem-resistant Enterobacteriaceae (CRE), all producing a rare carbapenemase, were identified among patients in an intensive care unit. This observation suggested a large, previously unrecognized plasmid transmission chain and prompted our investigation. METHODS: Electronic medical record reviews, infection control observations, and environmental sampling completed the epidemiologic outbreak investigation. A laboratory analysis, conducted on patient and environmental isolates, included long-read whole-genome sequencing to fully elucidate plasmid DNA structures. Bioinformatics analyses were applied to infer plasmid transmission chains and results were subsequently confirmed using plasmid conjugation experiments. RESULTS: We identified 14 Verona integron-encoded metallo-ß-lactamase (VIM)-producing CRE in 12 patients, and 1 additional isolate was obtained from a patient room sink drain. Whole-genome sequencing identified the horizontal transfer of blaVIM-1, a rare carbapenem resistance mechanism in the United States, via a promiscuous incompatibility group A/C2 plasmid that spread among 5 bacterial species isolated from patients and the environment. CONCLUSIONS: This investigation represents the largest known outbreak of VIM-producing CRE in the United States to date, which comprises numerous bacterial species and strains. We present evidence of in-hospital plasmid transmission, as well as environmental contamination. Our findings demonstrate the potential for 2 types of hospital-acquired infection outbreaks: those due to clonal expansion and those due to the spread of conjugative plasmids encoding antibiotic resistance across species.

Complete and Circularized Genome Assemblies of the Kroppenstedtia eburnea Genus Type Strain and the Kroppenstedtia pulmonis Species Type Strain with MiSeq and MinION Sequence Data.

Kroppenstedtia eburnea DSM 45196T and Kroppenstedtia pulmonis W9323T are aerobic, Gram-positive, filamentous, chemoorganotrophic thermoactinomycetes. Here, we report on the complete and circular genome assemblies generated using Illumina MiSeq and Oxford Nanopore Technologies MinION reads. Putative gene clusters predicted to be involved in the production of secondary metabolites were also identified.

Division of the genus Chryseobacterium: Observation of discontinuities in amino acid identity values, a possible consequence of major extinction events, guides transfer of nine species to the genus Epilithonimonas, eleven species to the genus Kaistella, and three species to the genus Halpernia gen. nov., with description of Kaistella daneshvariae sp. nov. and Epilithonimonas vandammei sp. nov. derived from clinical specimens.

The genus Chryseobacterium in the family Weeksellaceae is known to be polyphyletic. Amino acid identity (AAI) values were calculated from whole-genome sequences of species of the genus Chryseobacterium, and their distribution was found to be multi-modal. These naturally-occurring non-continuities were leveraged to standardise genus assignment of these species. We speculate that this multi-modal distribution is a consequence of loss of biodiversity during major extinction events, leading to the concept that a bacterial genus corresponds to a set of species that diversified since the Permian extinction. Transfer of nine species (Chryseobacterium arachidiradicis, Chryseobacterium bovis, Chryseobacterium caeni, Chryseobacterium hispanicum, Chryseobacterium hominis, Chryseobacterium hungaricum,, Chryseobacterium pallidum and Chryseobacterium zeae) to the genus Epilithonimonas and eleven (Chryseobacterium anthropi, Chryseobacterium antarcticum, Chryseobacterium carnis, Chryseobacterium chaponense, Chryseobacterium haifense, Chryseobacterium jeonii, Chryseobacterium montanum, Chryseobacterium palustre, Chryseobacterium solincola, Chryseobacterium treverense and Chryseobacterium yonginense) to the genus Kaistella is proposed. Two novel species are described: Kaistella daneshvariae sp. nov. and Epilithonimonas vandammei sp. nov. Evidence is presented to support the assignment of Planobacterium taklimakanense to a genus apart from Chryseobacterium, to which Planobacterium salipaludis comb nov. also belongs. The novel genus Halpernia is proposed, to contain the type species Halpernia frigidisoli comb. nov., along with Halpernia humi comb. nov., and Halpernia marina comb. nov.

Conserved Patterns of Symmetric Inversion in the Genome Evolution of Bordetella Respiratory Pathogens.

Whooping cough (pertussis), primarily caused by Bordetella pertussis, has resurged in the United States, and circulating strains exhibit considerable chromosome structural fluidity in the form of rearrangement and deletion. The genus Bordetella includes additional pathogenic species infecting various animals, some even causing pertussis-like respiratory disease in humans; however, investigation of their genome evolution has been limited. We studied chromosome structure in complete genome sequences from 167 Bordetella species isolates, as well as 469 B. pertussis isolates, to gain a generalized understanding of rearrangement patterns among these related pathogens. Observed changes in gene order primarily resulted from large inversions and were only detected in species with genomes harboring multicopy insertion sequence (IS) elements, most notably B. holmesii and B. parapertussis While genomes of B. pertussis contain >240 copies of IS481, IS elements appear less numerous in other species and yield less chromosome structural diversity through rearrangement. These data were further used to predict all possible rearrangements between IS element copies present in Bordetella genomes, revealing that only a subset is observed among circulating strains. Therefore, while it appears that rearrangement occurs less frequently in other species than in B. pertussis, these clinically relevant respiratory pathogens likely experience similar mutation of gene order. The resulting chromosome structural fluidity presents both challenges and opportunity for the study of Bordetella respiratory pathogens.IMPORTANCE Bordetella pertussis is the primary agent of whooping cough (pertussis). The Bordetella genus includes additional pathogens of animals and humans, including some that cause pertussis-like respiratory illness. The chromosome of B. pertussis has previously been shown to exhibit considerable structural rearrangement, but insufficient data have prevented comparable investigation in related species. In this study, we analyze chromosome structure variation in several Bordetella species to gain a generalized understanding of rearrangement patterns in this genus. Just as in B. pertussis, we observed inversions in other species that likely result from common mutational processes. We used these data to further predict additional, unobserved inversions, suggesting that specific genome structures may be preferred in each species.