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.

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.

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.

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.

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.

The Hidden Genomics of Chlamydia trachomatis.

The application of whole-genome sequencing has moved us on from sequencing single genomes to defining unravelling population structures in different niches, and at the -species, -serotype or even -genus level, and in local, national and global settings. This has been instrumental in cataloguing and revealing a huge a range of diversity in this bacterium, when at first we thought there was little. Genomics has challenged assumptions, added insight, as well as confusion and glimpses of truths. What is clear is that at a time when we start to realise the extent and nature of the diversity contained within a genus or a species like this, the huge depth of knowledge communities have developed, through cell biology, as well as the new found molecular approaches will be more precious than ever to link genotype to phenotype. Here we detail the technological developments and insights we have seen during the relatively short time since we began to see the hidden genome of Chlamydia trachomatis.

Massive expansion of Ubiquitination-related gene families within the Chlamydiae.

Gene loss, gain, and transfer play an important role in shaping the genomes of all organisms; however, the interplay of these processes in isolated populations, such as in obligate intracellular bacteria, is less understood. Despite a general trend towards genome reduction in these microbes, our phylogenomic analysis of the phylum Chlamydiae revealed that within the family Parachlamydiaceae, gene family expansions have had pronounced effects on gene content. We discovered that the largest gene families within the phylum are the result of rapid gene birth-and-death evolution. These large gene families are comprised of members harboring eukaryotic-like ubiquitination-related domains, such as F-box and BTB-box domains, marking the largest reservoir of these proteins found among bacteria. A heterologous type III secretion system assay suggests that these proteins function as effectors manipulating the host cell. The large disparity in copy number of members in these families between closely related organisms suggests that nonadaptive processes might contribute to the evolution of these gene families. Gene birth-and-death evolution in concert with genomic drift might represent a previously undescribed mechanism by which isolated bacterial populations diversify.

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.

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.

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.

Random transposon mutagenesis of Verrucomicrobium spinosum DSM 4136(T).

The Verrucomicrobia are a bacterial group of growing interest due to their environmental ubiquity as free-living and host-associated microbes. They also exhibit an unusual compartmentalized cell plan, shared with members of neighboring phyla that include the Planctomycete bacteria. However, Verrucomicrobia are currently difficult to study, due to a lack of available genetic tools that would permit robust testing of hypotheses formulated from ecological and genomic data. To our knowledge, there are no published studies describing the transformation of exogenous DNA into any members of the Verrucomicrobia (or the neighboring phylum containing Planctomycetes). Here, we present a procedure for the transformation of DNA into Verrucomicrobium spinosum DSM 4136(T) via electroporation and the first description of a random transposon mutant library in this organism. We anticipate that this approach could be applied successfully to other Verrucomicrobia, providing opportunities to test the role of predicted gene function in ecological interactions and identify genes associated with the distinctive Planctomycete-Verrucomicrobial cell plan.

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.

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.

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.

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.