Congruity of genomic and epidemiological data in modelling of local cholera outbreaks.

Cholera continues to be a global health threat. Understanding how cholera spreads between locations is fundamental to the rational, evidence-based design of intervention and control efforts. Traditionally, cholera transmission models have used cholera case-count data. More recently, whole-genome sequence data have qualitatively described cholera transmission. Integrating these data streams may provide much more accurate models of cholera spread; however, no systematic analyses have been performed so far to compare traditional case-count models to the phylodynamic models from genomic data for cholera transmission. Here, we use high-fidelity case-count and whole-genome sequencing data from the 1991 to 1998 cholera epidemic in Argentina to directly compare the epidemiological model parameters estimated from these two data sources. We find that phylodynamic methods applied to cholera genomics data provide comparable estimates that are in line with established methods. Our methodology represents a critical step in building a framework for integrating case-count and genomic data sources for cholera epidemiology and other bacterial pathogens.

Choclo virus (CHOV) recovered from deep metatranscriptomics of archived frozen tissues in natural history biorepositories.

BACKGROUND: Hantaviruses are negative-stranded RNA viruses that can sometimes cause severe disease in humans; however, they are maintained in mammalian host populations without causing harm. In Panama, sigmodontine rodents serve as hosts to transmissible hantaviruses. Due to natural and anthropogenic forces, these rodent populations are having increased contact with humans. METHODS: We extracted RNA and performed Illumina deep metatranscriptomic sequencing on Orthohantavirus seropositive museum tissues from rodents. We acquired sequence reads mapping to Choclo virus (CHOV, Orthohantavirus chocloense) from heart and kidney tissue of a two-decade old frozen museum sample from a Costa Rican pygmy rice rat (Oligoryzomys costaricensis) collected in Panama. Reads mapped to the CHOV reference were assembled and then validated by visualization of the mapped reads against the assembly. RESULTS: We recovered a 91% complete consensus sequence from a reference-guided assembly to CHOV with an average of 16X coverage. The S and M segments used in our phylogenetic analyses were nearly complete (98% and 99%, respectively). There were 1,199 ambiguous base calls of which 93% were present in the L segment. Our assembled genome varied 1.1% from the CHOV reference sequence resulting in eight nonsynonymous mutations. Further analysis of all publicly available partial S segment sequences support a clear relationship between CHOV clinical cases and O. costaricensis acquired strains. CONCLUSIONS: Viruses occurring at extremely low abundances can be recovered from deep metatranscriptomics of archival tissues housed in research natural history museum biorepositories. Our efforts resulted in the second CHOV genome publicly available. This genomic data is important for future surveillance and diagnostic tools as well as understanding the evolution and pathogenicity of CHOV.

Severe Acute Respiratory Syndrome Coronavirus-2 Lambda Variant Collected from a Child from Arkansas and Sequenced.

An eleven-year-old tested positive for SARS-CoV-2 Lambda variant. Sequencing was performed on the Oxford Nanopore and the Illumina NextSeq 500. Both platforms identified all 7 of the synonymous mutations in the sample, while all 28 nonsynonymous mutations were identified from Oxford Nanopore and 20 nonsynonymous mutations were identified from Illumina.

Genomic characterization of Staphylococcus aureus isolates causing osteoarticular infections in otherwise healthy children.

BACKGROUND: Pediatric osteoarticular infections are commonly caused by Staphylococcus aureus. The contribution of S. aureus genomic variability to pathogenesis of these infections is poorly described. METHODS: We prospectively enrolled 47 children over 3 1/2 years from whom S. aureus was isolated on culture-12 uninfected with skin colonization, 16 with skin abscesses, 19 with osteoarticular infections (four with septic arthritis, three with acute osteomyelitis, six with acute osteomyelitis and septic arthritis and six with chronic osteomyelitis). Isolates underwent whole genome sequencing, with assessment for 254 virulence genes and any mutations as well as creation of a phylogenetic tree. Finally, isolates were compared for their ability to form static biofilms and compared to the genetic analysis. RESULTS: No sequence types predominated amongst osteoarticular infections. Only genes involved in evasion of host immune defenses were more frequently carried by isolates from osteoarticular infections than from skin colonization (p = .02). Virulence gene mutations were only noted in 14 genes (three regulating biofilm formation) when comparing isolates from subjects with osteoarticular infections and those with skin colonization. Biofilm results demonstrated large heterogeneity in the isolates' capacity to form static biofilms, with healthy control isolates producing more robust biofilm formation. CONCLUSIONS: S. aureus causing osteoarticular infections are genetically heterogeneous, and more frequently harbor genes involved in immune evasion than less invasive isolates. However, virulence gene carriage overall is similar with infrequent mutations, suggesting that pathogenesis of S. aureus osteoarticular infections may be primarily regulated at transcriptional and/or translational levels.

Global landscape of SARS-CoV-2 genomic surveillance and data sharing.

Genomic surveillance has shaped our understanding of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants. We performed a global landscape analysis on SARS-CoV-2 genomic surveillance and genomic data using a collection of country-specific data. Here, we characterize increasing circulation of the Alpha variant in early 2021, subsequently replaced by the Delta variant around May 2021. SARS-CoV-2 genomic surveillance and sequencing availability varied markedly across countries, with 45 countries performing a high level of routine genomic surveillance and 96 countries with a high availability of SARS-CoV-2 sequencing. We also observed a marked heterogeneity of sequencing percentage, sequencing technologies, turnaround time and completeness of released metadata across regions and income groups. A total of 37% of countries with explicit reporting on variants shared less than half of their sequences of variants of concern (VOCs) in public repositories. Our findings indicate an urgent need to increase timely and full sharing of sequences, the standardization of metadata files and support for countries with limited sequencing and bioinformatics capacity.

Global landscape of SARS-CoV-2 genomic surveillance, public availability extent of genomic data, and epidemic shaped by variants.

Genomic surveillance has shaped our understanding of SARS-CoV-2 variants, which have proliferated globally in 2021.We collected country-specific data on SARS-CoV-2 genomic surveillance, sequencing capabilities, public genomic data from multiple public repositories, and aggregated publicly available variant data. Then, different proxies were used to estimate the sequencing coverage and public availability extent of genomic data, in addition to describing the global dissemination of variants. We found that the COVID-19 global epidemic clearly featured increasing circulation of Alpha since the start of 2021, which was rapidly replaced by the Delta variant starting around May 2021. SARS-CoV-2 genomic surveillance and sequencing availability varied markedly across countries, with 63 countries performing routine genomic surveillance and 79 countries with high availability of SARS-CoV-2 sequencing. We also observed a marked heterogeneity of sequenced coverage across regions and countries. Across different variants, 21-46% of countries with explicit reporting on variants shared less than half of their variant sequences in public repositories. Our findings indicated an urgent need to expand sequencing capacity of virus isolates, enhance the sharing of sequences, the standardization of metadata files, and supportive networks for countries with no sequencing capability.

Tracing Transmission of Sin Nombre Virus and Discovery of Infection in Multiple Rodent Species.

Sin Nombre orthohantavirus (SNV), a negative-sense, single-stranded RNA virus that is carried and transmitted by the North American deer mouse Peromyscus maniculatus, can cause infection in humans through inhalation of aerosolized excreta from infected rodents. This infection can lead to hantavirus cardiopulmonary syndrome (HCPS), which has an ∼36% case-fatality rate. We used reverse transcriptase quantitative PCR (RT-qPCR) to confirm SNV infection in a patient and identified SNV in lung tissues in wild-caught rodents from potential sites of exposure. Using viral whole-genome sequencing (WGS), we identified the likely site of transmission and discovered SNV in multiple rodent species not previously known to carry the virus. Here, we report, for the first time, the use of SNV WGS to pinpoint a likely site of human infection and identify SNV simultaneously in multiple rodent species in an area of known host-to-human transmission. These results will impact epidemiology and infection control for hantaviruses by tracing zoonotic transmission and investigating possible novel host reservoirs. IMPORTANCE Orthohantaviruses cause severe disease in humans and can be lethal in up to 40% of cases. Sin Nombre orthohantavirus (SNV) is the main cause of hantavirus disease in North America. In this study, we sequenced SNV from an infected patient and wild-caught rodents to trace the location of infection. We also discovered SNV in rodent species not previously known to carry SNV. These studies demonstrate for the first time the use of virus sequencing to trace the transmission of SNV and describe infection in novel rodent species.

Genotypic and phenotypic profiles of antibiotic-resistant bacteria isolated from hospitalised patients in Bangladesh.

OBJECTIVES: Characterisation of resistance phenotype and genotype is crucial to understanding the burden and transmission of antimicrobial resistance (AMR). This study aims to determine the spectrum of AMR and associated genes encoding aminoglycoside, macrolide and ß-lactam classes of antimicrobials in bacteria isolated from hospitalised patients in Bangladesh. METHODS: 430 bacterial isolates from patients with respiratory, intestinal, wound infections and typhoid fever, presenting to clinical care from 2015 to 2019, were examined. They included Escherichia coli (n = 85); Staphylococcus aureus (n = 84); Salmonella typhi (n = 82); Klebsiella pneumoniae (n = 42); Streptococcus pneumoniae (n = 36); coagulase-negative staphylococci (n = 28); Enterococcus faecalis (n = 27); Pseudomonas aeruginosa (n = 26); and Acinetobacter baumannii (n = 20). Reconfirmation of these clinical isolates and antimicrobial susceptibility tests was performed. PCR amplification using resistance gene-specific primers was done, and the amplified products were confirmed by Sanger sequencing. RESULTS: 53% of isolates were multidrug-resistant (MDR), including 97% of Escherichia coli. There was a year-wise gradual increase in MDR isolates from 2015 to 2018, and there was an almost twofold increase in the number of MDR strains isolated in 2019 (P = 0.00058). Among the 5 extended-spectrum ß-lactamases investigated, CTX-M-1 was the most prevalent (63%) followed by NDM-1 (22%); Escherichia coli was the major reservoir of these genes. The ermB (55%) and aac(6')-Ib (35%) genes were the most frequently detected macrolide and aminoglycoside resistance genes, respectively. CONCLUSION: MDR pathogens are highly prevalent in hospital settings of Bangladesh.

Emergence in late 2020 of multiple lineages of SARS-CoV-2 Spike protein variants affecting amino acid position 677.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein (S) plays critical roles in host cell entry. Non-synonymous substitutions affecting S are not uncommon and have become fixed in a number of SARS-CoV-2 lineages. A subset of such mutations enable escape from neutralizing antibodies or are thought to enhance transmission through mechanisms such as increased affinity for the cell entry receptor, angiotensin-converting enzyme 2 (ACE2). Independent genomic surveillance programs based in New Mexico and Louisiana contemporaneously detected the rapid rise of numerous clade 20G (lineage B.1.2) infections carrying a Q677P substitution in S. The variant was first detected in the US on October 23, yet between 01 Dec 2020 and 19 Jan 2021 it rose to represent 27.8% and 11.3% of all SARS-CoV-2 genomes sequenced from Louisiana and New Mexico, respectively. Q677P cases have been detected predominantly in the south central and southwest United States; as of 03 Feb 2021, GISAID data show 499 viral sequences of this variant from the USA. Phylogenetic analyses revealed the independent evolution and spread of at least six distinct Q677H sub-lineages, with first collection dates ranging from mid-August to late November 2020. Four 677H clades from clade 20G (B.1.2), 20A (B.1.234), and 20B (B.1.1.220, and B.1.1.222) each contain roughly 100 or fewer sequenced cases, while a distinct pair of clade 20G clusters are represented by 754 and 298 cases, respectively. Although sampling bias and founder effects may have contributed to the rise of S:677 polymorphic variants, the proximity of this position to the polybasic cleavage site at the S1/S2 boundary are consistent with its potential functional relevance during cell entry, suggesting parallel evolution of a trait that may confer an advantage in spread or transmission. Taken together, our findings demonstrate simultaneous convergent evolution, thus providing an impetus to further evaluate S:677 polymorphisms for effects on proteolytic processing, cell tropism, and transmissibility.

Coevolving Plasmids Drive Gene Flow and Genome Plasticity in Host-Associated Intracellular Bacteria.

Plasmids are important in microbial evolution and adaptation to new environments. Yet, carrying a plasmid can be costly, and long-term association of plasmids with their hosts is poorly understood. Here, we provide evidence that the Chlamydiae, a phylum of strictly host-associated intracellular bacteria, have coevolved with their plasmids since their last common ancestor. Current chlamydial plasmids are amalgamations of at least one ancestral plasmid and a bacteriophage. We show that the majority of plasmid genes are also found on chromosomes of extant chlamydiae. The most conserved plasmid gene families are predominantly vertically inherited, while accessory plasmid gene families show significantly increased mobility. We reconstructed the evolutionary history of plasmid gene content of an entire bacterial phylum over a period of around one billion years. Frequent horizontal gene transfer and chromosomal integration events illustrate the pronounced impact of coevolution with these extrachromosomal elements on bacterial genome dynamics in host-dependent microbes.

Genomics of the Argentinian cholera epidemic elucidate the contrasting dynamics of epidemic and endemic Vibrio cholerae.

In order to control and eradicate epidemic cholera, we need to understand how epidemics begin, how they spread, and how they decline and eventually end. This requires extensive sampling of epidemic disease over time, alongside the background of endemic disease that may exist concurrently with the epidemic. The unique circumstances surrounding the Argentinian cholera epidemic of 1992-1998 presented an opportunity to do this. Here, we use 490 Argentinian V. cholerae genome sequences to characterise the variation within, and between, epidemic and endemic V. cholerae. We show that, during the 1992-1998 cholera epidemic, the invariant epidemic clone co-existed alongside highly diverse members of the Vibrio cholerae species in Argentina, and we contrast the clonality of epidemic V. cholerae with the background diversity of local endemic bacteria. Our findings refine and add nuance to our genomic definitions of epidemic and endemic cholera, and are of direct relevance to controlling current and future cholera epidemics.

Draft Genome Sequences of Chlamydiales Bacterium STE3 and Neochlamydia sp. Strain AcF84, Endosymbionts of Acanthamoeba spp.

Chlamydiales bacterium STE3 and Neochlamydia sp. strain AcF84 are obligate intracellular symbionts of Acanthamoeba spp. isolated from the biofilm of a littoral cave wall and gills from striped tiger leaf fish, respectively. We report the draft genome sequences of these two environmental chlamydiae affiliated with the family Parachlamydiaceae.

Characterization of a thaumarchaeal symbiont that drives incomplete nitrification in the tropical sponge Ianthella basta.

Marine sponges represent one of the few eukaryotic groups that frequently harbour symbiotic members of the Thaumarchaeota, which are important chemoautotrophic ammonia-oxidizers in many environments. However, in most studies, direct demonstration of ammonia-oxidation by these archaea within sponges is lacking, and little is known about sponge-specific adaptations of ammonia-oxidizing archaea (AOA). Here, we characterized the thaumarchaeal symbiont of the marine sponge Ianthella basta using metaproteogenomics, fluorescence in situ hybridization, qPCR and isotope-based functional assays. 'Candidatus Nitrosospongia ianthellae' is only distantly related to cultured AOA. It is an abundant symbiont that is solely responsible for nitrite formation from ammonia in I. basta that surprisingly does not harbour nitrite-oxidizing microbes. Furthermore, this AOA is equipped with an expanded set of extracellular subtilisin-like proteases, a metalloprotease unique among archaea, as well as a putative branched-chain amino acid ABC transporter. This repertoire is strongly indicative of a mixotrophic lifestyle and is (with slight variations) also found in other sponge-associated, but not in free-living AOA. We predict that this feature as well as an expanded and unique set of secreted serpins (protease inhibitors), a unique array of eukaryotic-like proteins, and a DNA-phosporothioation system, represent important adaptations of AOA to life within these ancient filter-feeding animals.

Mitigating Anticipated Effects of Systematic Errors Supports Sister-Group Relationship between Xenacoelomorpha and Ambulacraria.

Xenoturbella and the acoelomorph worms (Xenacoelomorpha) are simple marine animals with controversial affinities. They have been placed as the sister group of all other bilaterian animals (Nephrozoa hypothesis), implying their simplicity is an ancient characteristic [1, 2]; alternatively, they have been linked to the complex Ambulacraria (echinoderms and hemichordates) in a clade called the Xenambulacraria [3-5], suggesting their simplicity evolved by reduction from a complex ancestor. The difficulty resolving this problem implies the phylogenetic signal supporting the correct solution is weak and affected by inadequate modeling, creating a misleading non-phylogenetic signal. The idea that the Nephrozoa hypothesis might be an artifact is prompted by the faster molecular evolutionary rate observed within the Acoelomorpha. Unequal rates of evolution are known to result in the systematic artifact of long branch attraction, which would be predicted to result in an attraction between long-branch acoelomorphs and the outgroup, pulling them toward the root [6]. Other biases inadequately accommodated by the models used can also have strong effects, exacerbated in the context of short internal branches and long terminal branches [7]. We have assembled a large and informative dataset to address this problem. Analyses designed to reduce or to emphasize misleading signals show the Nephrozoa hypothesis is supported under conditions expected to exacerbate errors, and the Xenambulacraria hypothesis is preferred in conditions designed to reduce errors. Our reanalyses of two other recently published datasets [1, 2] produce the same result. We conclude that the Xenacoelomorpha are simplified relatives of the Ambulacraria.

High quality reference genomes for toxigenic and non-toxigenic Vibrio cholerae serogroup O139.

Toxigenic Vibrio cholerae of the O139 serogroup have been responsible for several large cholera epidemics in South Asia, and continue to be of clinical and historical significance today. This serogroup was initially feared to represent a new, emerging V. cholerae clone that would lead to an eighth cholera pandemic. However, these concerns were ultimately unfounded. The majority of clinically relevant V. cholerae O139 isolates are closely related to serogroup O1, biotype El Tor V. cholerae, and comprise a single sublineage of the seventh pandemic El Tor lineage. Although related, these V. cholerae serogroups differ in several fundamental ways, in terms of their O-antigen, capsulation phenotype, and the genomic islands found on their chromosomes. Here, we present four complete, high-quality genomes for V. cholerae O139, obtained using long-read sequencing. Three of these sequences are from toxigenic V. cholerae, and one is from a bacterium which, although classified serologically as V. cholerae O139, lacks the CTXφ bacteriophage and the ability to produce cholera toxin. We highlight fundamental genomic differences between these isolates, the V. cholerae O1 reference strain N16961, and the prototypical O139 strain MO10. These sequences are an important resource for the scientific community, and will improve greatly our ability to perform genomic analyses of non-O1 V. cholerae in the future. These genomes also offer new insights into the biology of a V. cholerae serogroup that, from a genomic perspective, is poorly understood.

The history, genome and biology of NCTC 30: a non-pandemic Vibrio cholerae isolate from World War One.

The sixth global cholera pandemic lasted from 1899 to 1923. However, despite widespread fear of the disease and of its negative effects on troop morale, very few soldiers in the British Expeditionary Forces contracted cholera between 1914 and 1918. Here, we have revived and sequenced the genome of NCTC 30, a 102-year-old Vibrio cholerae isolate, which we believe is the oldest publicly available live V. cholerae strain in existence. NCTC 30 was isolated in 1916 from a British soldier convalescent in Egypt. We found that this strain does not encode cholera toxin, thought to be necessary to cause cholera, and is not part of V. cholerae lineages responsible for the pandemic disease. We also show that NCTC 30, which predates the introduction of penicillin-based antibiotics, harbours a functional ß-lactamase antibiotic resistance gene. Our data corroborate and provide molecular explanations for previous phenotypic studies of NCTC 30 and provide a new high-quality genome sequence for historical, non-pandemic V. cholerae.

Publisher Correction: Genomic insights into the 2016-2017 cholera epidemic in Yemen.

In the HTML version of this Letter, the affiliations for authors Andrew S. Azman, Dhirendra Kumar and Thandavarayan Ramamurthy were inverted (the PDF and print versions of the Letter were correct); the affiliations have been corrected online.

Genomic insights into the 2016-2017 cholera epidemic in Yemen.

Yemen is currently experiencing, to our knowledge, the largest cholera epidemic in recent history. The first cases were declared in September 2016, and over 1.1 million cases and 2,300 deaths have since been reported1. Here we investigate the phylogenetic relationships, pathogenesis and determinants of antimicrobial resistance by sequencing the genomes of Vibrio cholerae isolates from the epidemic in Yemen and recent isolates from neighbouring regions. These 116 genomic sequences were placed within the phylogenetic context of a global collection of 1,087 isolates of the seventh pandemic V. cholerae serogroups O1 and O139 biotype El Tor2-4. We show that the isolates from Yemen that were collected during the two epidemiological waves of the epidemic1-the first between 28 September 2016 and 23 April 2017 (25,839 suspected cases) and the second beginning on 24 April 2017 (more than 1 million suspected cases)-are V. cholerae serotype Ogawa isolates from a single sublineage of the seventh pandemic V. cholerae O1 El Tor (7PET) lineage. Using genomic approaches, we link the epidemic in Yemen to global radiations of pandemic V. cholerae and show that this sublineage originated from South Asia and that it caused outbreaks in East Africa before appearing in Yemen. Furthermore, we show that the isolates from Yemen are susceptible to several antibiotics that are commonly used to treat cholera and to polymyxin B, resistance to which is used as a marker of the El Tor biotype.