- Volume 7, Issue 6, 2021

Herpes simplex virus serotype 2 (HSV-2) is a ubiquitous human pathogen that causes recurrent genital infections and ulcerations. Many HSV-2 strains with different biological properties have been identified, but only the genomes of HSV-2 strains HG52, SD90e and 333 have been reported as complete and fully characterized sequences. We de novo assembled, annotated and manually curated the complete genome sequence of HSV-2 strain MS, a highly neurovirulent strain, originally isolated from a multiple sclerosis patient. We resolved both DNA ends, as well as the complex inverted repeats regions present in HSV genomes, usually undisclosed in previous published partial herpesvirus genomes, using long reads from Pacific Biosciences (PacBio) technology. Additionally, we identified isomeric genomes by determining the alternative relative orientation of unique fragments in the genome of the sequenced viral population. Illumina short-read sequencing was crucial to examine genetic variability, such as nucleotide polymorphisms, insertion/deletions and sequence determinants of strain-specific virulence factors. We used Illumina data to fix two disrupted open reading frames found in coding homopolymers after PacBio assembly. These results support the combination of long- and short-read sequencing technologies as a precise and effective approach for the accurate de novo assembly and curation of complex microbial genomes.

The occurrence of multiple strains of a bacterial pathogen such as   M. tuberculosis   or   C. difficile   within a single human host, referred to as a mixed infection, has important implications for both healthcare and public health. However, methods for detecting it, and especially determining the proportion and identities of the underlying strains, from WGS (whole-genome sequencing) data, have been limited. In this paper we introduce SplitStrains, a novel method for addressing these challenges. Grounded in a rigorous statistical model, SplitStrains not only demonstrates superior performance in proportion estimation to other existing methods on both simulated as well as real   M. tuberculosis   data, but also successfully determines the identity of the underlying strains. We conclude that SplitStrains is a powerful addition to the existing toolkit of analytical methods for data coming from bacterial pathogens and holds the promise of enabling previously inaccessible conclusions to be drawn in the realm of public health microbiology.

Enterotoxigenic Escherichia coli (ETEC) expressing the colonization pili CFA/I are common causes of diarrhoeal infections in humans. Here, we use a combination of transposon mutagenesis and transcriptomic analysis to identify genes and pathways that contribute to ETEC persistence in water environments and colonization of a mammalian host. ETEC persisting in water exhibit a distinct RNA expression profile from those growing in richer media. Multiple pathways were identified that contribute to water survival, including lipopolysaccharide biosynthesis and stress response regulons. The analysis also indicated that ETEC growing in vivo in mice encounter a bottleneck driving down the diversity of colonizing ETEC populations.

The vaginal microbiota, normally characterized by lactobacilli presence, is crucial for vaginal health. Members belonging to L. crispatus and L. gasseri species exert crucial protective functions against pathogens, although a total comprehension of factors that influence their dominance in healthy women is still lacking. Here we investigated the complete genome sequence and comprehensive phenotypic profile of L. crispatus strain BC5 and L. gasseri strain BC12, two vaginal strains featured by anti-bacterial and anti-viral activities. Phenotype microarray (PM) results revealed an improved capacity of BC5 to utilize different carbon sources as compared to BC12, although some specific carbon sources that can be associated to the human diet were only metabolized by BC12, i.e. uridine, amygdalin, tagatose. Additionally, the two strains were mostly distinct in the capacity to utilize the nitrogen sources under analysis. On the other hand, BC12 showed tolerance/resistance towards twice the number of stressors (i.e. antibiotics, toxic metals etc.) with respect to BC5. The divergent phenotypes observed in PM were supported by the identification in either BC5 or BC12 of specific genetic determinants that were found to be part of the core genome of each species. The PM results in combination with comparative genome data provide insights into the possible environmental factors and genetic traits supporting the predominance of either L. crispatus BC5 or L. gasseri BC12 in the vaginal niche, giving also indications for metabolic predictions at the species level.

Although the beneficial effects of probiotics are likely to be associated with their ability to colonize the gut, little is known about the characteristics of good colonizers. In a systematic analysis of the comparative genomics, we tried to elucidate the genomic contents that account for the distinct host adaptability patterns of Lactobacillus and Bifidobacterium species. The Bifidobacterium species, with species-level phylogenetic structures affected by recombination among strains, broad mucin-foraging activity, and dietary-fibre-degrading ability, represented niche conservatism and tended to be host-adapted. The Lactobacillus species stretched across three lifestyles, namely free-living, nomadic and host-adapted, as characterized by the variations of bacterial occurrence time, guanine–cytosine (GC) content and genome size, evolution event frequency, and the presence of human-adapted bacterial genes. The numbers and activity of host-adapted factors, such as bile salt hydrolase and intestinal tissue-anchored elements, were distinctly distributed among the three lifestyles. The strains of the three lifestyles could be separated with such a collection of colonization-related genomic content (genes, genome size and GC content). Thus, our work provided valuable information for rational selection and gut engraftment prediction of probiotics. Here, we have found many interesting predictive results for bacterial gut fitness, which will be validated in vitro and in vivo.

Pseudomonas corrugata constitute one of the phylogenomic subgroups within the Pseudomonas fluorescens species complex and include both plant growth-promoting rhizobacteria (PGPR) and plant pathogenic bacteria. Previous studies suggest that the species diversity of this group remains largely unexplored together with frequent misclassification of strains. Using more than 1800 sequenced Pseudomonas genomes we identified 121 genomes belonging to the P. corrugata subgroup. Intergenomic distances obtained using the genome-to-genome blast distance (GBDP) algorithm and the determination of digital DNA–DNA hybridization values were further used for phylogenomic and clustering analyses, which revealed 29 putative species clusters, of which only five correspond to currently named species within the subgroup. Comparative and functional genome-scale analyses also support the species status of these clusters. The search for PGPR and plant pathogenic determinants showed that approximately half of the genomes analysed could have a pathogenic behaviour based on the presence of a pathogenicity genetic island, while all analysed genomes possess PGPR traits. Finally, this information together with the characterization of phenotypic traits, allows the reclassification proposal of Pseudomonas fluorescens F113 as Pseudomonas ogarae sp. nov., nom rev., type strain F113T (=DSM 112162T=CECT 30235T), which is substantiated by genomic, functional genomics and phenotypic differences with their closest type strains.

Species of the floating, freshwater fern Azolla form a well-characterized symbiotic association with the non-culturable cyanobacterium Nostoc azollae, which fixes nitrogen for the plant. However, several cyanobacterial strains have over the years been isolated and cultured from Azolla from all over the world. The genomes of 10 of these strains were sequenced and compared with each other, with other symbiotic cyanobacterial strains, and with similar strains that were not isolated from a symbiotic association. The 10 strains fell into three distinct groups: six strains were nearly identical to the non-symbiotic strain, Nostoc ( Anabaena ) variabilis ATCC 29413; three were similar to the symbiotic strain, Nostoc punctiforme , and one, Nostoc sp. 2RC, was most similar to non-symbiotic strains of Nostoc linckia. However, Nostoc sp. 2RC was unusual because it has three sets of nitrogenase genes; it has complete gene clusters for two distinct Mo-nitrogenases and an alternative V-nitrogenase. Genes for Mo-nitrogenase, sugar transport, chemotaxis and pili characterized all the symbiotic strains. Several of the strains infected the liverwort Blasia, including N. variabilis ATCC 29413, which did not originate from Azolla but rather from a sewage pond. However, only Nostoc sp. 2RC, which produced highly motile hormogonia, was capable of high-frequency infection of Blasia. Thus, some of these strains, which grow readily in the laboratory, may be useful in establishing novel symbiotic associations with other plants.

Located at the tip of cell surface glycoconjugates, sialic acids are at the forefront of host–microbe interactions and, being easily liberated by sialidase enzymes, are used as metabolites by numerous bacteria, particularly by pathogens and commensals living on or near diverse mucosal surfaces. These bacteria rely on specific transporters for the acquisition of host-derived sialic acids. Here, we present the first comprehensive genomic and phylogenetic analysis of bacterial sialic acid transporters, leading to the identification of multiple new families and subfamilies. Our phylogenetic analysis suggests that sialic acid-specific transport has evolved independently at least eight times during the evolution of bacteria, from within four of the major families/superfamilies of bacterial transporters, and we propose a robust classification scheme to bring together a myriad of different nomenclatures that exist to date. The new transporters discovered occur in diverse bacteria, including Spirochaetes , Bacteroidetes , Planctomycetes and Verrucomicrobia , many of which are species that have not been previously recognized to have sialometabolic capacities. Two subfamilies of transporters stand out in being fused to the sialic acid mutarotase enzyme, NanM, and these transporter fusions are enriched in bacteria present in gut microbial communities. Our analysis supports the increasing experimental evidence that competition for host-derived sialic acid is a key phenotype for successful colonization of complex mucosal microbiomes, such that a strong evolutionary selection has occurred for the emergence of sialic acid specificity within existing transporter architectures.

Respiratory syncytial viruses (RSVs) are an important cause of mortality worldwide and a major cause of respiratory tract infections in children, driving development of vaccine candidates. However, there are large gaps in our knowledge of the local evolutionary and transmission dynamics of RSVs, particularly in understudied regions such as the Middle East. To address this gap, we sequenced the complete genomes of 58 RSVA and 27 RSVB samples collected in a paediatric cohort in Amman, Jordan, between 2010 and 2013. RSVA and RSVB co-circulated during each winter epidemic of RSV in Amman, and each epidemic comprised multiple independent viral introductions of RSVA and RSVB. However, RSVA and RSVB alternated in dominance across years, potential evidence of immunological interactions. Children infected with RSVA tended to be older than RSVB-infected children [30 months versus 22.4 months, respectively (P value = 0.02)], and tended to developed bronchopneumonia less frequently than those with RSVB, although the difference was not statistically significant (P value = 0.06). Differences in spatial patterns were investigated, and RSVA lineages were often identified in multiple regions in Amman, whereas RSVB introductions did not spread beyond a single region of the city, although these findings were based on small sample sizes. Multiple RSVA genotypes were identified in Amman, including GA2 viruses as well as three viruses from the ON1 sub-genotype that emerged in 2009 and are now the dominant genotype circulating worldwide. As vaccine development advances, further sequencing of RSV is needed to understand viral ecology and transmission, particularly in under-studied locations.

Salmonella enterica serovar Typhimurium is the leading cause of salmonellosis in Australia, and the ability to identify outbreaks and their sources is vital to public health. Here, we examined the utility of whole-genome sequencing (WGS), including complete genome sequencing with Oxford Nanopore technologies, in examining 105 isolates from an endemic multi-locus variable number tandem repeat analysis (MLVA) type over 5 years. The MLVA type was very homogeneous, with 90 % of the isolates falling into groups with a five SNP cut-off. We developed a new two-step approach for outbreak detection using WGS. The first clustering at a zero single nucleotide polymorphism (SNP) cut-off was used to detect outbreak clusters that each occurred within a 4 week window and then a second clustering with dynamically increased SNP cut-offs were used to generate outbreak investigation clusters capable of identifying all outbreak cases. This approach offered optimal specificity and sensitivity for outbreak detection and investigation, in particular of those caused by endemic MLVA types or clones with low genetic diversity. We further showed that inclusion of complete genome sequences detected no additional mutational events for genomic outbreak surveillance. Phylogenetic analysis found that the MLVA type was likely to have been derived recently from a single source that persisted over 5 years, and seeded numerous sporadic infections and outbreaks. Our findings suggest that SNP cut-offs for outbreak cluster detection and public-health surveillance should be based on the local diversity of the relevant strains over time. These findings have general applicability to outbreak detection of bacterial pathogens.

Klebsiella pneumoniae is a frequent cause of nosocomial and severe community-acquired infections. Multidrug-resistant (MDR) and hypervirulent (hv) strains represent major threats, and tracking their emergence, evolution and the emerging convergence of MDR and hv traits is of major importance. We employed whole-genome sequencing (WGS) to study the evolution and epidemiology of a large longitudinal collection of clinical K. pneumoniae isolates from the H301 hospital in Beijing, China. Overall, the population was highly diverse, although some clones were predominant. Strains belonging to clonal group (CG) 258 were dominant, and represented the majority of carbapenemase-producers. While CG258 strains showed high diversity, one clone, ST11-KL47, represented the majority of isolates, and was highly associated with the KPC-2 carbapenemase and several virulence factors, including a virulence plasmid. The second dominant clone was CG23, which is the major hv clone globally. While it is usually susceptible to multiple antibiotics, we found some isolates harbouring MDR plasmids encoding for ESBLs and carbapenemases. We also reported the local emergence of a recently described high-risk clone, ST383. Conversely to strains belonging to CG258, which are usually associated to KPC-2, ST383 strains seem to readily acquire carbapenemases of different types. Moreover, we found several ST383 strains carrying the hypervirulence plasmid. Overall, we detected about 5 % of simultaneous carriage of AMR genes (ESBLs or carbapenemases) and hypervirulence genes. Tracking the emergence and evolution of such strains, causing severe infections with limited treatment options, is fundamental in order to understand their origin and evolution and to limit their spread. This article contains data hosted by Microreact.

The increasing use of PCR for the detection of gastrointestinal pathogens in hospital laboratories in England has improved the detection of Shiga toxin-producing Escherichia coli (STEC), and the diagnosis of haemolytic uraemic syndrome (HUS). We aimed to analyse the microbiological characteristics and phylogenetic relationships of STEC O26:H11, clonal complex (CC) 29, in England to inform surveillance, and to assess the threat to public health. There were 502 STEC belonging to CC29 isolated between 2014 and 2019, of which 416 were from individual cases. The majority of isolates belonged to one of three major sequence types (STs), ST16 (n=37), ST21 (n=350) and ST29 (n=24). ST16 and ST29 were mainly isolated from cases reporting recent travel abroad. Within ST21, there were three main clades associated with domestic acquisition. All three domestic clades had Shiga toxin subtype gene (stx) profiles associated with causing severe clinical outcomes including STEC-HUS, specifically either stx1a, stx2a or stx1a/stx2a. Isolates from the same patient, same household or same outbreak with an established source for the most part fell within 5-SNP single linkage clusters. There were 19 5-SNP community clusters, of which six were travel-associated and one was an outbreak of 16 cases caused by the consumption of contaminated salad leaves. Of the remaining 12 clusters, 9/12 were either temporally or geographically related or both. Exposure to foodborne STEC O26:H11 ST21 capable of causing severe clinical outcomes, including STEC-HUS, is an emerging risk to public health in England. The lack of comprehensive surveillance of this STEC serotype is a concern, and there is a need to expand the implementation of methods capable of detecting STEC in local hospital settings.

The bacterial genotoxin colibactin interferes with the eukaryotic cell cycle by causing dsDNA breaks. It has been linked to bacterially induced colorectal cancer in humans. Colibactin is encoded by a 54 kb genomic region in Enterobacteriaceae . The colibactin genes commonly co-occur with the yersiniabactin biosynthetic determinant. Investigating the prevalence and sequence diversity of the colibactin determinant and its linkage to the yersiniabactin operon in prokaryotic genomes, we discovered mainly species-specific lineages of the colibactin determinant and classified three main structural settings of the colibactin–yersiniabactin genomic region in Enterobacteriaceae . The colibactin gene cluster has a similar but not identical evolutionary track to that of the yersiniabactin operon. Both determinants could have been acquired on several occasions and/or exchanged independently between enterobacteria by horizontal gene transfer. Integrative and conjugative elements play(ed) a central role in the evolution and structural diversity of the colibactin–yersiniabactin genomic region. Addition of an activating and regulating module (clbAR) to the biosynthesis and transport module (clbB-S) represents the most recent step in the evolution of the colibactin determinant. In a first attempt to correlate colibactin expression with individual lineages of colibactin determinants and different bacterial genetic backgrounds, we compared colibactin expression of selected enterobacterial isolates in vitro. Colibactin production in the tested Klebsiella species and Citrobacter koseri strains was more homogeneous and generally higher than that in most of the Escherichia coli isolates studied. Our results improve the understanding of the diversity of colibactin determinants and its expression level, and may contribute to risk assessment of colibactin-producing enterobacteria.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the novel coronavirus responsible for the COVID-19 pandemic, continues to cause a significant public-health burden and disruption globally. Genomic epidemiology approaches point to most countries in the world having experienced many independent introductions of SARS-CoV-2 during the early stages of the pandemic. However, this situation may change with local lockdown policies and restrictions on travel, leading to the emergence of more geographically structured viral populations and lineages transmitting locally. Here, we report the first SARS-CoV-2 genomes from Palestine sampled from early March 2020, when the first cases were observed, through to August of 2020. SARS-CoV-2 genomes from Palestine fall across the diversity of the global phylogeny, consistent with at least nine independent introductions into the region. We identify one locally predominant lineage in circulation represented by 50 Palestinian SARS-CoV-2, grouping with genomes generated from Israel and the UK. We estimate the age of introduction of this lineage to 05/02/2020 (16/01/2020–19/02/2020), suggesting SARS-CoV-2 was already in circulation in Palestine predating its first detection in Bethlehem in early March. Our work highlights the value of ongoing genomic surveillance and monitoring to reconstruct the epidemiology of COVID-19 at both local and global scales.

The COVID-19 pandemic has spread rapidly throughout the world. In the UK, the initial peak was in April 2020; in the county of Norfolk (UK) and surrounding areas, which has a stable, low-density population, over 3200 cases were reported between March and August 2020. As part of the activities of the national COVID-19 Genomics Consortium (COG-UK) we undertook whole genome sequencing of the SARS-CoV-2 genomes present in positive clinical samples from the Norfolk region. These samples were collected by four major hospitals, multiple minor hospitals, care facilities and community organizations within Norfolk and surrounding areas. We combined clinical metadata with the sequencing data from regional SARS-CoV-2 genomes to understand the origins, genetic variation, transmission and expansion (spread) of the virus within the region and provide context nationally. Data were fed back into the national effort for pandemic management, whilst simultaneously being used to assist local outbreak analyses. Overall, 1565 positive samples (172 per 100 000 population) from 1376 cases were evaluated; for 140 cases between two and six samples were available providing longitudinal data. This represented 42.6 % of all positive samples identified by hospital testing in the region and encompassed those with clinical need, and health and care workers and their families. In total, 1035 cases had genome sequences of sufficient quality to provide phylogenetic lineages. These genomes belonged to 26 distinct global lineages, indicating that there were multiple separate introductions into the region. Furthermore, 100 genetically distinct UK lineages were detected demonstrating local evolution, at a rate of ~2 SNPs per month, and multiple co-occurring lineages as the pandemic progressed. Our analysis: identified a discrete sublineage associated with six care facilities; found no evidence of reinfection in longitudinal samples; ruled out a nosocomial outbreak; identified 16 lineages in key workers which were not in patients, indicating infection control measures were effective; and found the D614G spike protein mutation which is linked to increased transmissibility dominates the samples and rapidly confirmed relatedness of cases in an outbreak at a food processing facility. The large-scale genome sequencing of SARS-CoV-2-positive samples has provided valuable additional data for public health epidemiology in the Norfolk region, and will continue to help identify and untangle hidden transmission chains as the pandemic evolves.

Clostridioides difficile is the leading cause of healthcare-associated infectious diarrhoea. However, it is increasingly appreciated that healthcare-associated infections derive from both community and healthcare environments, and that the primary sites of C. difficile transmission may be strain-dependent. We conducted a multisite genomic epidemiology study to assess differential genomic evidence of healthcare vs community spread for two of the most common C. difficile strains in the USA: sequence type (ST) 1 (associated with ribotype 027) and ST2 (associated with ribotype 014/020). We performed whole-genome sequencing and phylogenetic analyses on 382 ST1 and ST2 C. difficile isolates recovered from stool specimens collected during standard clinical care at 3 geographically distinct US medical centres between 2010 and 2017. ST1 and ST2 isolates both displayed some evidence of phylogenetic clustering by study site, but clustering was stronger and more apparent in ST1, consistent with our healthcare-based study more comprehensively sampling local transmission of ST1 compared to ST2 strains. Analyses of pairwise single-nucleotide variant (SNV) distance distributions were also consistent with more evidence of healthcare transmission of ST1 compared to ST2, with 44 % of ST1 isolates being within two SNVs of another isolate from the same geographical collection site compared to 5.5 % of ST2 isolates (P-value=<0.001). Conversely, ST2 isolates were more likely to have close genetic neighbours across disparate geographical sites compared to ST1 isolates, further supporting non-healthcare routes of spread for ST2 and highlighting the potential for misattributing genomic similarity among ST2 isolates to recent healthcare transmission. Finally, we estimated a lower evolutionary rate for the ST2 lineage compared to the ST1 lineage using Bayesian timed phylogenomic analyses, and hypothesize that this may contribute to observed differences in geographical concordance among closely related isolates. Together, these findings suggest that ST1 and ST2, while both common causes of C. difficile infection in hospitals, show differential reliance on community and hospital spread. This conclusion supports the need for strain-specific criteria for interpreting genomic linkages and emphasizes the importance of considering differences in the epidemiology of circulating strains when devising interventions to reduce the burden of C. difficile infections.

As part of the ongoing studies with clinically relevant   Klebsiella   spp., we characterized the genomes of three clinical GES-5-positive ST138 strains originally identified as Klebsiella oxytoca. bla OXY gene, average nucleotide identity and phylogenetic analyses showed the strains to be   Klebsiella michiganensis  . Affiliation of the strains to ST138 led us to demonstrate that the current multi-locus sequence typing scheme for   K. oxytoca   can be used to distinguish members of this genetically diverse complex of bacteria. The strains encoded the kleboxymycin biosynthetic gene cluster (BGC), previously only found in   K. oxytoca   strains and one strain of   Klebsiella grimontii  . The finding of this BGC, associated with antibiotic-associated haemorrhagic colitis, in   K. michiganensis   led us to carry out a wide-ranging study to determine the prevalence of this BGC in   Klebsiella   spp. Of 7170 publicly available   Klebsiella   genome sequences screened, 88 encoded the kleboxymycin BGC. All BGC-positive strains belonged to the   K. oxytoca   complex, with strains of four (  K. oxytoca  ,   K. pasteurii  ,   K. grimontii  ,   K. michiganensis  ) of the six species of complex found to encode the complete BGC. In addition to being found in   K. grimontii   strains isolated from preterm infants, the BGC was found in   K. oxytoca   and   K. michiganensis   metagenome-assembled genomes recovered from neonates. Detection of the kleboxymycin BGC across the   K. oxytoca   complex may be of clinical relevance and this cluster should be included in databases characterizing virulence factors, in addition to those characterizing BGCs.

Members of the bacterial genus Vibrio utilize chitin both as a metabolic substrate and a signal to activate natural competence. Vibrio cholerae is a bacterial enteric pathogen, sub-lineages of which can cause pandemic cholera. However, the chitin metabolic pathway in V. cholerae has been dissected using only a limited number of laboratory strains of this species. Here, we survey the complement of key chitin metabolism genes amongst 195 diverse V. cholerae . We show that the gene encoding GbpA, known to be an important colonization and virulence factor in pandemic isolates, is not ubiquitous amongst V. cholerae . We also identify a putatively novel chitinase, and present experimental evidence in support of its functionality. Our data indicate that the chitin metabolic pathway within V. cholerae is more complex than previously thought, and emphasize the importance of considering genes and functions in the context of a species in its entirety, rather than simply relying on traditional reference strains.