- Volume 9, Issue 9, 2023

Clostridioides difficile , the leading cause of antibiotic-associated diarrhoea worldwide, is a genetically diverse species which can metabolise a number of nutrient sources upon colonising a dysbiotic gut environment. Trehalose, a disaccharide sugar consisting of two glucose molecules bonded by an α 1,1-glycosidic bond, has been hypothesised to be involved in the emergence of C. difficile hypervirulence due to its increased utilisation by the RT027 and RT078 strains. Here, growth in trehalose as the sole carbon source was shown to be non-uniform across representative C. difficile strains, even though the genes for its metabolism were induced. Growth in trehalose reduced the expression of genes associated with toxin production and sporulation in the C. difficile R20291 (RT027) and M120 (RT078) strains in vitro, suggesting an inhibitory effect on virulence factors. Interestingly, the R20291 TreR transcriptional regulatory protein appeared to possess an activator function as its DNA-binding ability was increased in the presence of its effector, trehalose-6-phosphate. Using RNA-sequencing analysis, we report the identification of a putative trehalose metabolism pathway which is induced during growth in trehalose: this has not been previously described within the C. difficile species. These data demonstrate the metabolic diversity exhibited by C. difficile which warrants further investigation to elucidate the molecular basis of trehalose metabolism within this important gut pathogen.

Respiratory diphtheria is a serious infection caused by toxigenic Corynebacterium diphtheriae , and disease transmission mainly occurs through respiratory droplets. Between 2017 and 2019, a large diphtheria outbreak among forcibly displaced Myanmar nationals densely settled in Bangladesh was investigated. Here we utilized whole-genome sequencing (WGS) to characterize recovered isolates of C. diphtheriae and two co-circulating non-diphtheritic Corynebacterium (NDC) species – C. pseudodiphtheriticum and C. propinquum. C. diphtheriae isolates recovered from all 53 positive cases in this study were identified as toxigenic biovar mitis, exhibiting intermediate resistance to penicillin, and formed four phylogenetic clusters circulating among multiple refugee camps. Additional sequenced isolates collected from two patients showed co-colonization with non-toxigenic C. diphtheriae biovar gravis, one of which exhibited decreased susceptibility to the first-line antibiotics and harboured a novel 23-kb multidrug resistance plasmid. Results of phylogenetic reconstruction and virulence-related gene contents of the recovered NDC isolates indicated they were likely commensal organisms, though 80.4 %(45/56) were not susceptible to erythromycin, and most showed high minimum inhibition concentrations against azithromycin. These results demonstrate the high resolution with which WGS can aid molecular investigation of diphtheria outbreaks, through the quantification of bacterial genetic relatedness, as well as the detection of virulence factors and antibiotic resistance markers among case isolates.

Non-typhoidal Salmonella are extremely diverse and different serovars can exhibit varied phenotypes, including host adaptation and the ability to cause clinical illness in animals and humans. In the USA, Salmonella enterica serovar Kentucky is infrequently found to cause human illness, despite being the top serovar isolated from broiler chickens. Conversely, in Europe, this serovar falls in the top 10 serovars linked to human salmonellosis. Serovar Kentucky is polyphyletic and has two lineages, Kentucky-I and Kentucky-II; isolates belonging to Kentucky-I are frequently isolated from poultry in the USA, while Kentucky-II isolates tend to be associated with human illness. In this study, we analysed whole-genome sequences and associated metadata deposited in public databases between 2017 and 2021 by federal agencies to determine serovar Kentucky incidence across different animal and human sources. Of 5151 genomes, 90.3 % were from isolates that came from broilers, while 5.9 % were from humans and 3.0 % were from cattle. Kentucky-I isolates were associated with broilers, while isolates belonging to Kentucky-II and a new lineage, Kentucky-III, were more commonly associated with cattle and humans. Very few serovar Kentucky isolates were associated with turkey and swine sources. Phylogenetic analyses showed that Kentucky-III genomes were more closely related to Kentucky-I, and this was confirmed by CRISPR-typing and multilocus sequence typing (MLST). In a macrophage assay, serovar Kentucky-II isolates were able to replicate over eight times better than Kentucky-I isolates. Analysis of virulence factors showed unique patterns across these three groups, and these differences may be linked to their association with different hosts.

For the last two decades, the human infection frequency of Escherichia coli O157 (O157) in Scotland has been 2.5-fold higher than in England and Wales. Results from national cattle surveys conducted in Scotland and England and Wales in 2014/2015 were combined with data on reported human clinical cases from the same time frame to determine if strain differences in national populations of O157 in cattle could be associated with higher human infection rates in Scotland. Shiga toxin subtype (Stx) and phage type (PT) were examined within and between host (cattle vs human) and nation (Scotland vs England and Wales). For a subset of the strains, whole genome sequencing (WGS) provided further insights into geographical and host association. All three major O157 lineages (I, II, I/II) and most sub-lineages (Ia, Ib, Ic, IIa, IIb, IIc) were represented in cattle and humans in both nations. While the relative contribution of different reservoir hosts to human infection is unknown, WGS analysis indicated that the majority of O157 diversity in human cases was captured by isolates from cattle. Despite comparable cattle O157 prevalence between nations, strain types were localized. PT21/28 (sub-lineage Ic, Stx2a+) was significantly more prevalent in Scottish cattle [odds ratio (OR) 8.7 (2.3–33.7; P<0.001] and humans [OR 2.2 (1.5–3.2); P<0.001]. In England and Wales, cattle had a significantly higher association with sub-lineage IIa strains [PT54, Stx2c; OR 5.6 (1.27–33.3); P=0.011] while humans were significantly more closely associated with sub-lineage IIb [PT8, Stx1 and Stx2c; OR 29 (4.9–1161); P<0.001]. Therefore, cattle farms in Scotland were more likely to harbour Stx2a+O157 strains compared to farms in E and W (P<0.001). There was evidence of limited cattle strain migration between nations and clinical isolates from one nation were more similar to cattle isolates from the same nation, with sub-lineage Ic (mainly PT21/28) exhibiting clear national association and evidence of local transmission in Scotland. While we propose the higher rate of O157 clinical cases in Scotland, compared to England and Wales, is a consequence of the nationally higher level of Stx2a+O157 strains in Scottish cattle, we discuss the multiple additional factors that may also contribute to the different infection rates between these nations.

Respiratory syncytial virus (RSV), or human orthopneumovirus, is a major cause of acute lower respiratory infection (ALRI), particularly in young children, causing significant morbidity and mortality. We used pathogen genomics to characterize the population structure and genetic signatures of RSV isolates circulating in children in New South Wales between 2016 and 2018 and to understand the evolutionary dynamics of these strains in the context of publicly available RSV genomes from the region and globally. Whole-genome phylogenetic analysis demonstrated the co-circulation of a few major RSV clades in the paediatric population from Sydney. The whole-genome-based genotypes A23 (RSV-A ON1-like genotype) and B6 (RSV-B BA9-like genotype) were the predominant RSV-A and RSV-B genotypes circulating during the study period, respectively. These genotypes were characterized with high levels of diversity of predicted N- and O-linked glycosylation patterns in both the G and F glycoproteins. Interestingly, a novel 72-nucleotide triplication in the sequence that corresponds to the C-terminal region of the G gene was identified in four of the A23 genotype sequenced in this study. Consistently, the population dynamics analysis demonstrated a continuous increase in the effective population size of A23 and B6 genotypes globally. Further investigations including functional mapping of mutations and identifying the impact of sequence changes on virus fitness are highly required. This study highlights the potential impact of an integrated approach that uses WG-based phylogeny and studying selective pressure events in understanding the emergence and dissemination of RSV genotypes.

Pantoea stewartii , a plant pathogen, is primarily transmitted through contaminated seeds and insect vectors, with the corn flea beetle (Chaetocnema pulicaria) being the primary carrier. P. stewartii is a bacterium belonging to the order Enterobacterales and can lead to crop diseases that have a significant economic impact worldwide. Due to its high potential for spread, P. stewartii is classified as a quarantine organism in numerous countries. Despite its impact on agriculture, the limited genome sequences of P. stewartii hamper understanding of its pathogenicity and host specificity, and the development of effective control strategies. In this study, a P. stewartii strain (C10109_Jinnung) was discovered in the faecal matter of the Critically Endangered western ground parrot/kyloring (Pezoporus flaviventris) in Australia, which to our knowledge is the first reported P. stewartii genome from a bird source. Whole-genome sequencing and phylogenomic analysis of strain C10109_Jinnung, obtained from a captive psittacine, provides new insights into the genetic diversity and potential transmission route for the spread of P. stewartii beyond insects and plants, where P. stewartii is typically studied. Our findings provide new insights into the potential transmission route for spread of P. stewartii and expand the known transmission agents beyond insects and plants. Expanding the catalogue of P. stewartii genomes is fundamental to improving understanding of the pathogenicity, evolution and dissemination, and to develop effective control strategies to reduce the substantial economic losses associated with P. stewartii in various crops and the potential impact of endangered animal species.

A multidrug-resistant strain of Klebsiella pneumoniae (Kp) sequence type (ST) 1788, an otherwise uncommon ST worldwide, was isolated from 65 patients at 11 hospitals and 11 general practices across South and West Wales, UK, between February 2019 and November 2021. A collection of 97 Kp ST1788 isolates (including 94 from Wales) was analysed to investigate the diversity and spread across Wales and to identify molecular marker(s) to aid development of a strain-specific real-time PCR. Whole genome sequencing (WGS) was performed with Illumina technology and the data were used to perform phylogenetic analyses. Pan-genome analysis of further Kp genome collections was used to identify an ST1788-specific gene target; a real-time PCR was then validated against a panel of 314 strains and 218 broth-enriched screening samples. Low genomic diversity was demonstrated amongst the 94 isolates from Wales. Evidence of spread within and across healthcare facilities was found. A yersiniabactin locus and the KL2 capsular locus were identified in 85/94 (90.4 %) and 94/94 (100 %) genomes respectively; bla SHV-232, bla TEM-1, bla CTX-M-15 and bla OXA-1 were simultaneously carried by 86/94 (91.5 %) isolates; 4/94 (4.3 %) isolates also carried bla OXA-48 carbapenemase. Aminoglycoside and fluoroquinolone resistance markers were found in 94/94 (100 %) and 86/94 (91.5 %) isolates respectively. The ST1788-specific real-time PCR was 100 % sensitive and specific. Our analyses demonstrated recent clonal expansion and spread of Kp ST1788 in the community and across healthcare facilities in South and West Wales with isolates carrying well-defined antimicrobial resistance and virulence markers. An ST1788-specific marker was also identified, enabling rapid and reliable preliminary characterization of isolates by real-time PCR. This study confirms the utility of WGS in investigating novel strains and in aiding proactive implementation of molecular tools to assist infection control specialists.

Aeromonas veronii is an opportunistic pathogen that affects both fish and mammals, including humans, leading to bacteraemia, sepsis, meningitis and even death. The increasing virulence and drug resistance of A. veronii are of significant concern and pose a severe risk to public safety. The Type I restriction-modification (RM) system, which functions as a bacterial defence mechanism, can influence gene expression through DNA methylation. However, little research has been conducted to explore its origin, evolutionary path, and relationship to virulence and drug resistance in A. veronii . In this study, we analysed the pan-genome of 233 A . veronii strains, and the results indicated that it was 'open', meaning that A. veronii has acquired additional genes from other species. This suggested that A. veronii had the potential to adapt and evolve rapidly, which might have contributed to its drug resistance. One Type I methyltransferase (MTase) and two complete Type I RM systems were identified, namely AveC4I, AveC4II and AveC4III in A. veronii strain C4, respectively. Notably, AveC4I was exclusive to A. veronii C4. Phylogenetic analysis revealed that AveC4I was derived from horizontal gene transfer from Thiocystis violascens and exchanged genes with the human pathogen Comamonas kerstersii . Single molecule real-time sequencing was applied to identify the motif methylated by AveC4I, which was unique and not recognized by any reported MTases in the REBASE database. We also annotated the functions and pathways of the genes containing the motif, revealing that AveC4I may control drug resistance in A. veronii C4. Our findings provide new insight on the mechanisms underlying drug resistance in pathogenic bacteria. By identifying the specific genes and pathways affected by AveC4I, this study may aid in the development of new therapeutic approaches to combat A. veronii infections.

The spread of multidrug-resistant tuberculosis (MDR-TB) is a growing problem in many countries worldwide. Resistance to one of the primary first-line drugs, rifampicin, is caused by mutations in the Mycobacterium tuberculosis rpoB gene. So-called borderline rpoB mutations confer low-level resistance, in contrast to more common rpoB mutations which confer high-level resistance. While some borderline mutations show lower fitness in vitro than common mutations, their in vivo fitness is currently unknown. We used a dataset of 394 whole genome sequenced MDR-TB isolates from Bangladesh, representing around 44 % of notified MDR-TB cases over 6 years, to look at differences in transmission clustering between isolates with borderline rpoB mutations and those with common rpoB mutations. We found a relatively low percentage of transmission clustering in the dataset (34.8 %) but no difference in clustering between different types of rpoB mutations. Compensatory mutations in rpoA, rpoB, and rpoC were associated with higher levels of transmission clustering as were lineages two, three, and four relative to lineage one. Young people as well as patients with high sputum smear positive TB were more likely to be in a transmission cluster. Our findings show that although borderline rpoB mutations have lower in vitro growth potential this does not translate into lower transmission potential or in vivo fitness. Proper detection of these mutations is crucial to ensure they do not go unnoticed and spread MDR-TB within communities.

Streptococcus pneumoniae (pneumococcus) is a leading vaccine-preventable cause of childhood invasive disease. Nigeria has the second highest pneumococcal disease burden globally, with an estimated ~49 000 child deaths caused by pneumococcal infections each year. Ten-valent pneumococcal conjugate vaccine (GSK; PCV10) was introduced in December 2014 in a phased approach. However, few studies have characterized the disease-causing pneumococci from Nigeria. This study assessed the prevalence of serotypes, antibiotic susceptibility and genomic lineages using whole genome sequencing and identified lineages that could potentially escape PCV10 (GSK). We also investigated the potential differences in pneumococcal lineage features between children with and without sickle cell disease. A collection of 192 disease-causing pneumococcal isolates was obtained from Kano (n=189) and Abuja (n=3) states, Nigeria, between 1 January 2014 and 31 May 2018. The majority (99 %, 190/192) of specimens were recovered from children aged 5 years or under. Among them, 37 children had confirmed or traits of sickle cell disease. Our findings identified 25 serotypes expressed by 43 Global Pneumococcal Sequence Clusters (GPSCs) and 85 sequence types (STs). The most common serotypes were 14 (18 %, n=35), 6B (16 %, n=31), 1 (9 %, n=17), 5 (9 %, n=17) and 6A (9 %, n=17); all except serotype 6A are included in PCV10 (GSK). PCV10 (SII; PNEUMOSIL) and PCV13 formulations include serotypes 6A and 19A which would increase the overall coverage from 67 % by PCV10 (GSK) to 78 and 82 %, respectively. The pneumococcal lineages were a mix of globally spreading and unique local lineages. Following the use of PCV10 (GSK), GPSC5 expressing serotype 6A, GPSC10 (19A), GPSC26 (12F and 46) and GPSC627 (9L) are non-vaccine type lineages that could persist and potentially expand under vaccine-selective pressure. Approximately half (52 %, 99/192) of the pneumococcal isolates were resistant to the first-line antibiotic penicillin and 44 % (85/192) were multidrug-resistant. Erythromycin resistance was very low (2 %, 3/192). There was no significant difference in clinical manifestation, serotype prevalence or antibiotic resistance between children with and without traits of or confirmed sickle cell disease. In summary, our findings show that a high percentage of the pneumococcal disease were caused by the serotypes that are covered by currently available vaccines. Given the low prevalence of resistance, macrolide antibiotics, such as erythromycin, should be considered as an option to treat pneumococcal disease in Nigeria. However, appropriate use of macrolide antibiotics should be vigilantly monitored to prevent the potential increase in macrolide resistance.

  Bradyrhizobium   comprises a diverse group of bacteria with various lifestyles. Although best known for their nodule-based nitrogen-fixation in symbiosis with legumes, a select group of bradyrhizobia are also capable of photosynthesis. This ability seems to be rare among rhizobia, and its origin and evolution in these bacteria remain a subject of substantial debate. Therefore, our aim here was to investigate the distribution and evolution of photosynthesis in   Bradyrhizobium   using comparative genomics and representative genomes from closely related taxa in the families Nitrobacteraceae, Methylobacteriaceae, Boseaceae and Paracoccaceae. We identified photosynthesis gene clusters (PGCs) in 25 genomes belonging to three different   Bradyrhizobium   lineages, notably the so-called Photosynthetic,   B. japonicum   and   B. elkanii   supergroups. Also, two different PGC architectures were observed. One of these, PGC1, was present in genomes from the Photosynthetic supergroup and in three genomes from a species in the   B. japonicum   supergroup. The second cluster, PGC2, was also present in some strains from the   B. japonicum   supergroup, as well as in those from the   B. elkanii   supergroup. PGC2 was largely syntenic to the cluster found in   Rhodopseudomonas palustris   and   Tardiphaga  . Bayesian ancestral state reconstruction unambiguously showed that the ancestor of   Bradyrhizobium   lacked a PGC and that it was acquired horizontally by various lineages. Maximum-likelihood phylogenetic analyses of individual photosynthesis genes also suggested multiple acquisitions through horizontal gene transfer, followed by vertical inheritance and gene losses within the different lineages. Overall, our findings add to the existing body of knowledge on   Bradyrhizobium  ’s evolution and provide a meaningful basis from which to explore how these PGCs and the photosynthesis itself impact the physiology and ecology of these bacteria.

Nitrofurantoin is a broad-spectrum first-line antimicrobial used for managing uncomplicated urinary tract infection (UTI). Loss-of-function mutations in chromosomal genes nfsA, nfsB and ribE of Escherichia coli are known to reduce nitrofurantoin susceptibility. Here, we report the discovery of nitrofurantoin heteroresistance in E. coli clinical isolates and a novel genetic mechanism associated with this phenomenon. Subpopulations with lower nitrofurantoin susceptibility than major populations (hereafter, nitrofurantoin-resistant subpopulations) in two E. coli blood isolates (previously whole-genome sequenced) were identified using population analysis profiling. Each isolate was known to have a loss-of-function mutation in nfsA. From each isolate, four nitrofurantoin-resistant isolates were derived at a nitrofurantoin concentration of 32 mg l−1, and a comparator isolate was obtained without any nitrofurantoin exposure. Genomes of derived isolates were sequenced on Illumina and Nanopore MinION systems. Genetic variation between isolates was determined based on genome assemblies and read mapping. Nitrofurantoin minimum inhibitory concentrations (MICs) of both blood isolates were 64 mg l−1, with MICs of major nitrofurantoin-susceptible populations varying from 4 to 8 mg l−1. Two to 99 c.f.u. per million demonstrated growth at the nitrofurantoin concentration of 32 mg l−1, which is distinct from that of a homogeneously susceptible or resistant isolate. Derived nitrofurantoin-resistant isolates had 11–66 kb deletions in chromosomal regions harbouring nfsB, and all deletions were immediately adjacent to IS1-family insertion sequences. Our findings demonstrate that the IS1-associated large-scale genetic deletion is a hitherto unrecognized mechanism of nitrofurantoin heteroresistance and could compromise UTI management. Further, frequencies of resistant subpopulations from nitrofurantoin-heteroresistant isolates may challenge conventional nitrofurantoin susceptibility testing in clinical settings.