Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

The genomic epidemiology of Neisseria meningitidis carriage from a randomised controlled trial of 4CMenB vaccination in an asymptomatic adolescent population

Background

Oropharyngeal carriage of Neisseria meningitidis is frequent during adolescence, representing a major source of invasive meningococcal disease. This study examined the impact of a serogroup B vaccination (Bexsero, GSK 4CMenB) programme on adolescent N. meningitidis carriage using genomic data.

Methods

A total 34,489 oropharyngeal samples were collected as part of a state-wide cluster randomised-controlled trial in South Australia during 2017 and 2018 (NCT03089086). Samples were screened for the presence of N. meningitidis DNA by porA PCR prior to culture. Whole genome sequencing was performed on all 1772 N. meningitidis culture isolates and their genomes were analysed.

Findings

Unencapsulated meningococci were predominant at baseline (36.3% of isolates), followed by MenB (31.0%), and MenY (20.5%). Most MenB were ST-6058 from hyperinvasive cc41/44, or ST-32 and ST-2870 from cc32. For MenY, ST-23 and ST-1655 from cc23 were prevalent. Meningococcal carriage was mostly unchanged due to the vaccination programme; however, a significant reduction in ST-53 capsule-null meningococci prevalence was observed in 2018 compared to 2017 (OR = 0.52; 95% CI: 0.30-0.87, p = 0.0106). This effect was larger in the vaccinated compared to the control group (OR = 0.37; 95% CI: 0.12-0.98, p = 0.0368).

Interpretation

While deployment of the 4CMenB vaccination did not alter the carriage of hyperinvasive MenB in the vaccinated population, it altered the carriage of other N. meningitidis sequence types following the vaccination program. Our findings suggest 4CMenB vaccination is unlikely to reduce transmission of hyperinvasive N. meningitidis strains and therefore ongoing targeted vaccination is likely a more effective public health intervention.

Funding

This work was funded by GlaxoSmithKline Biologicals SA.

Genomic-wide analysis approach revealed genomic similarity for environmental Mexican S. Oranienburg genomes

As the human population grows, an increase in food trade is needed. This elevates the risk of epidemiological outbreaks. One of the prevalent pathogens associated with food production in Mexico has been Salmonella Oranienburg. Effective surveillance systems require microbial genetic knowledge. The objective of this work is to describe the genetic composition of Mexican S. Oranienburg genomes. For that, 53 strains from different environmental sources were isolated and sequenced. Additionally, 109 S. Oranienburg genomes were downloaded. Bioinformatic analyses were used to explore the clonal complex and genomic relatedness. A major clonal group formed by ST23 was identified comprising four STs. 202 SNPs were found the maximum difference among isolates. Virulence genes for host invasion and colonization as rpoS, fimbria type 1, and, T3SS were found common for all isolates. This study suggests that Mexican S. Oranienburg strains are potential pathogens circulating continuously in the region between host and non-host environments.

Exploring the genetic basis of natural resistance to microcins

Enterobacteriaceae produce an arsenal of antimicrobial compounds including microcins, ribosomally produced antimicrobial peptides showing diverse structures and mechanisms of action. Microcins target close relatives of the producing strain to promote its survival. Their narrow spectrum of antibacterial activity makes them a promising alternative to conventional antibiotics, as it should decrease the probability of resistance dissemination and collateral damage to the host's microbiota. To assess the therapeutic potential of microcins, there is a need to understand the mechanisms of resistance to these molecules. In this study, we performed genomic analyses of the resistance to four microcins [microcin C, a nucleotide peptide; microcin J25, a lasso peptide; microcin B17, a linear azol(in)e-containing peptide; and microcin E492, a siderophore peptide] on a collection of 54 Enterobacteriaceae from three species: Escherichia coli, Salmonella enterica and Klebsiella pneumoniae. A gene-targeted analysis revealed that about half of the microcin-resistant strains presented mutations of genes involved in the microcin mechanism of action, especially those involved in their uptake (fhuA, fepA, cirA and ompF). A genome-wide association study did not reveal any significant correlations, yet relevant genetic elements were associated with microcin resistance. These were involved in stress responses, biofilm formation, transport systems and acquisition of immunity genes. Additionally, microcin-resistant strains exhibited several mutations within genes involved in specific metabolic pathways, especially for S. enterica and K. pneumoniae.

Epidemiology of Bartonella henselae infection in pet and stray cats in Croatia with risk factors analysis

BACKGROUND

Cats are the primary reservoirs of the bacterium Bartonella henselae, the main cause of cat-scratch disease in humans. The main vector of the bacterium is the cat flea, Ctenocephalides felis. In southeastern Europe, data are lacking on the prevalence of B. henselae infection in cats, the strains of B. henselae involved and the risk factors associated with the infection.

METHODS

Blood samples collected in ethylenediaminetetraacetic acid-containing tubes from 189 domestic cats (156 pet cats and 33 stray cats) from Zagreb, the capital city of Croatia, and 10 counties throughout Croatia were cultured for Bartonella spp. Following culture, bacterial isolates were genotyped at eight loci after using PCR to amplify 16S ribosomal RNA (rRNA) and the internal transcribed spacer region between the 16S and 23S rRNA sequences. Univariate and multivariate logistic regression were used to identify risk factors for B. henselae infection in cats.

RESULTS

Bartonella spp. was detected in 31 cats (16.4%), and subsequent genotyping at the eight loci revealed B. henselae in all cases. Thirty complete multilocus sequence typing profiles were obtained, and the strains were identified as four sequence types that had been previously reported, namely ST5 (56.7%), ST6 (23.3%), ST1 (13.3%) and ST24 (3.3%), as well as a novel sequence type, ST33 (3.3%). The univariate analysis revealed a significantly higher risk of B. henselae infection in cats residing in coastal areas of Croatia (odds ratio [OR] 2.592, 95% confidence interval [CI] 1.150-5.838; P = 0.0191) and in cats with intestinal parasites (OR 3.207, 95% CI 1.088-9.457; P = 0.0279); a significantly lower risk was identified in cats aged > 1 year (OR 0.356, 95% CI 0.161-0.787; P = 0.0247) and in cats sampled between April and September (OR 0.325, 95% CI 0.147-0.715; P = 0.005). The multivariate analysis that controlled for age showed a positive association with the presence of intestinal parasites (OR 4.241, 95% CI 1.243-14.470; P = 0.0119) and coastal residence (OR 2.567, 95% CI 1.114-5.915; P = 0.0216) implying increased risk of infection, and a negative association with sampling between April and September (OR 0.379, 95% CI 0.169-0.848; P = 0.018) implying a decreased risk of infection. After controlling for the season, an increased risk of infection remained for the coastal region (OR 2.725, 95% CI 1.200-6.186; P = 0.012).

CONCLUSIONS

Bartonella henselae is prevalent throughout Croatia and is a public health threat. Environmental and host factors can significantly affect the risk of infection, and these should be explored in more detail. The presence of intestinal parasites highlights the need to eliminate the flea vector, Ctenocephalides felis, as the most effective approach to control infections in cats and humans.

Characterizing the diversity and commensal origins of penA mosaicism in the genus Neisseria

Mosaic penA alleles formed through horizontal gene transfer (HGT) have been instrumental to the rising incidence of ceftriaxone-resistant gonococcal infections. Although interspecies HGT of regions of the penA gene between Neisseria gonorrhoeae and commensal Neisseria species has been described, knowledge concerning which species are the most common contributors to mosaic penA alleles is limited, with most studies examining only a small number of alleles. Here, we investigated the origins of recombinant penA alleles through in silico analyses that incorporated 1700 penA alleles from 35 513 Neisseria isolates, comprising 15 different Neisseria species. We identified Neisseria subflava and Neisseria cinerea as the most common source of recombinant sequences in N. gonorrhoeae penA. This contrasted with Neisseria meningitidis penA, for which the primary source of recombinant DNA was other meningococci, followed by Neisseria lactamica. Additionally, we described the distribution of polymorphisms implicated in antimicrobial resistance in penA, and found that these are present across the genus. These results provide insight into resistance-related changes in the penA gene across human-associated Neisseria species, illustrating the importance of genomic surveillance of not only the pathogenic Neisseria, but also of the oral niche-associated commensals from which these pathogens are sourcing key genetic variation.

Comparative genomics of Brucella abortus and Brucella melitensis unravels the gene sharing, virulence factors and SNP diversity among the standard, vaccine and field strains

Brucella abortus and Brucella melitensis are the primary etiological agents of brucellosis in large and small ruminants, respectively. There are limited comparative genomic studies involving Brucella strains that explore the relatedness among both species. In this study, we involved strains (n=44) representing standard, vaccine and Indian field origin for pangenome, single nucleotide polymorphism (SNP) and phylogenetic analysis. Both species shared a common gene pool representing 2884 genes out of a total 3244 genes. SNP-based phylogenetic analysis indicated higher SNP diversity among B. melitensis (3824) strains in comparison to B. abortus (540) strains, and a clear demarcation was identified between standard/vaccine and field strains. The analysis for virulence genes revealed that virB3, virB7, ricA, virB5, ipx5, wbkC, wbkB, and acpXL genes were highly conserved in most of the Brucella strains. Interestingly, virB10 gene was found to have high variability among the B. abortus strains. The cgMLST analysis revealed distinct sequence types for the standard/vaccine and field strains. B. abortus strains from north-eastern India fall within similar sequence type differing from other strains. In conclusion, the analysis revealed a highly shared core genome among two Brucella species. SNP analysis revealed B. melitensis strains exhibit high diversity as compared to B. abortus strains. Strains with absence or high polymorphism of virulence genes can be exploited for the development of novel vaccine candidates effective against both B. abortus and B. melitensis.

Average Nucleotide Identity based Staphylococcus aureus strain grouping allows identification of strain-specific genes in the pangenome

Staphylococcus aureus causes both hospital and community acquired infections in humans worldwide. Due to the high incidence of infection S. aureus is also one of the most sampled and sequenced pathogens today, providing an outstanding resource to understand variation at the bacterial subspecies level. We processed and downsampled 83,383 public S. aureus Illumina whole genome shotgun sequences and 1,263 complete genomes to produce 7,954 representative substrains. Pairwise comparison of core gene Average Nucleotide Identity (ANI) revealed a natural boundary of 99.5% that could be used to define 145 distinct strains within the species. We found that intermediate frequency genes in the pangenome (present in 10-95% of genomes) could be divided into those closely linked to strain background ("strain-concentrated") and those highly variable within strains ("strain-diffuse"). Non-core genes had different patterns of chromosome location; notably, strain-diffuse associated with prophages, strain-concentrated with the vSaβ genome island and rare genes (<10% frequency) concentrated near the origin of replication. Antibiotic genes were enriched in the strain-diffuse class, while virulence genes were distributed between strain-diffuse, strain-concentrated, core and rare classes. This study shows how different patterns of gene movement help create strains as distinct subspecies entities and provide insight into the diverse histories of important S. aureus functions.

Phenotypic and genotypic characterization of Neisseria gonorrhoeae isolates among individuals at high risk for sexually transmitted diseases in Zurich, Switzerland

BACKGROUND

While ceftriaxone resistance remains scarce in Switzerland, global Neisseria gonorrhoeae (NG) antimicrobial resistance poses an urgent threat. This study describes clinical characteristics in MSM (men who have sex with men) diagnosed with NG infection and analyses NG resistance by phenotypic and genotypic means.

METHODS

Data of MSM enrolled in three clinical cohorts with a positive polymerase chain reaction test (PCR) for NG were analysed between January 2019 and December 2021 and linked with antibiotic susceptibility testing. Bacterial isolates were subjected to whole genome sequencing (WGS).

RESULTS

Of 142 participants, 141 (99%) were MSM and 118 (84%) living with HIV. Participants were treated with ceftriaxone (N = 79), azithromycin (N = 2), or a combination of both (N = 61). No clinical or microbiological failures were observed. From 182 positive PCR samples taken, 23 were available for detailed analysis. Based on minimal inhibitory concentrations (MICs), all isolates were susceptible to ceftriaxone, gentamicin, cefixime, cefpodoxime, ertapenem, zoliflodacin, and spectinomycin. Resistance to azithromycin, tetracyclines and ciprofloxacin was observed in 10 (43%), 23 (100%) and 11 (48%) of the cases, respectively. Analysis of WGS data revealed combinations of resistance determinants that matched with the corresponding phenotypic resistance pattern of each isolate.

CONCLUSION

Among the MSM diagnosed with NG mainly acquired in Switzerland, ceftriaxone MICs were low for a subset of bacterial isolates studied and no treatment failures were observed. For azithromycin, high occurrences of in vitro resistance were found. Gentamicin, cefixime, cefpodoxime, ertapenem, spectinomycin, and zoliflodacin displayed excellent in vitro activity against the 23 isolates underscoring their potential as alternative agents to ceftriaxone.

Role of transcriptomic and genomic analyses in improving the comprehension of cefiderocol activity in Acinetobacter baumannii

The mechanisms of action and resistance of cefiderocol (FDC) in Acinetobacter baumannii are still not fully elucidated, but iron transport systems have been evoked in its entry into the cell to reach the penicillin-binding proteins (PBPs). To capture the dynamics of gene expression related to FDC action in various conditions, we report on the genomic and transcriptomic features of seven A. baumannii strains with different FDC susceptibility, focusing on the variants in genes associated with β-lactam resistance and the expression of the siderophore biosynthesis and transport systems acinetobactin and baumannoferrin. We also investigated the expression of the TonB energy transduction system (ETS) and siderophore receptors piuA and pirA. The four clinical samples belonged to the same clonal complex (CC2), and the two strains with the highest FDC MICs showed peculiar variants in PBP2 and ampC. Similarly, the two clinical strains with the lowest MICs shared variants in an outer membrane protein as well as ampC. Gene expression analyses highlighted the up-regulation of the acinetobactin and baumannoferrin genes in response to iron depletion and a down-regulation in the presence of high iron concentrations. In response to FDC, gene expression seemed strain-dependent, probably due to the different metabolic features of each strain. Overall, FDC activates the ETS, confirming the active import of the drug; baumannoferrin, more than acinetobactin, appeared stimulated by FDC in an iron-depleted medium. In conclusion, iron transport systems play a clear role in the FDC uptake, and their expression likely contributes to MIC variation together with β-lactam resistance determinants.IMPORTANCEAcinetobacter baumannii poses a threat to healthcare due to its ability to give difficult-to-treat infections as a consequence of our shortage of antibiotic molecules active on this multidrug-resistant bacterium. Cefiderocol (FDC) represents one of the few drugs active on A. baumannii, and to preserve its activity, this study explored the transcriptomic and genomic features of seven strains with varying susceptibility to FDC. Transcriptomic analyses revealed the different effects of FDC on iron transport systems, promoting mainly baumannoferrin expression-thus more likely related to FDC entry-and the energy transduction systems. These findings suggest that not all iron transport systems are equally involved in FDC entry into A. baumannii cells. Finally, mutations in PBPs and β-lactamases may contribute to the resistance onset. Overall, the study sheds light on the importance of iron availability and metabolic differences in FDC resistance, offering insights into understanding the evolution of resistance in A. baumannii strains.

Comparative genomic analysis of uropathogenic Escherichia coli strains from women with recurrent urinary tract infection

Introduction

Recurrent urinary tract infections (RUTIs) caused by uropathogenic Escherichia coli are costly public health problems impacting patients' quality of life.

Aim

In this work, a comparative genomics analysis of three clinical RUTI strains isolated from bladder biopsy specimens was performed.

Materials and methods

One hundred seventy-two whole genomes of urinary tract E. coli strains were selected from the NCBI database. The search for virulence factors, fitness genes, regions of interest, and genetic elements associated with resistance was manually carried out. The phenotypic characterization of antibiotic resistance, haemolysis, motility, and biofilm formation was performed. Moreover, adherence and invasion assays with human bladder HTB-5 cells, and transmission electron microscopy (TEM) were performed.

Results

The UTI-1_774U and UTI-3_455U/ST1193 strains were associated with the extraintestinal pathotypes, and the UTI-2_245U/ST295 strain was associated with the intestinal pathotype, according to a phylogenetic analysis of 172 E. coli urinary strains. The three RUTI strains were of clinical, epidemiological, and zoonotic relevance. Several resistance genes were found within the plasmids of these strains, and a multidrug resistance phenotype was revealed. Other virulence genes associated with CFT073 were not identified in the three RUTI strains (genes for type 1 and P fimbriae, haemolysin hlyA, and sat toxin). Quantitative adherence analysis showed that UTI-1_774U was significantly (p < 0.0001) more adherent to human bladder HTB-5 cells. Quantitative invasion analysis showed that UTI-2_245U was significantly more invasive than the control strains. No haemolysis or biofilm activity was detected in the three RUTI strains. The TEM micrographs showed the presence of short and thin fimbriae only in the UTI-2_245U strain.

Conclusion

The high variability and genetic diversity of the RUTI strains indicate that are a mosaic of virulence, resistance, and fitness genes that could promote recurrence in susceptible patients.

Human Keratinocyte Entry of Noninvasive Streptococcus dysgalactiae Subsp. equisimilis from Humans and Companion Animals: Relatedness with Lancefield Group, Source, Virulence-Associated Genes, and Antimicrobial Resistance Phenotype

We evaluated the cell invasion ability (CIA) of non-invasive Streptococcus dysgalactiae subsp. equisimilis using human keratinocytes and determined the association of CIA populations with their hosts and microbiological traits. Forty-two isolates from humans and companion animals were selected with host information. In addition to CIA, virulence-associated gene (VAG, spegg-ska-scpA-inlA-sicG-brpA-prtF1-prtF2-lmb-cbp-srtp1-srtp2) profiling, emm genotyping, multilocus sequence typing, and antimicrobial resistance (AMR) phenotyping/genotyping were performed. We designated CIA values higher than the mean of all isolates as high-frequency and those lower than the mean as low-frequency. Differences in the CIA between the different sources and Lancefield groups were assessed. We analyzed the association between high- and low-frequency CIA and VAG, emm genotype, sequence type/clonal complex, and AMR phenotype/genotype. Based on the mean (19.368 colony-forming units/100 cells) of 42 isolates, eight isolates had high-frequency CIA, whereas 34 had low-frequency CIA. We found an association between low-frequency CIA population and group G isolates, as well as a link between high-frequency CIA population and group C isolates. We also observed associations between low-frequency CIA population and oral/respiratory tract origin, ska, scpA, and lmb detection, and the AMR phenotype. Our observations suggest potential associations between high-/low-frequency CIA and the group, source, VAG, and AMR phenotypes.

Identification of Burkholderia gladioli pv. cocovenenans in Black Fungus and Efficient Recognition of Bongkrekic Acid and Toxoflavin Producing Phenotype by Back Propagation Neural Network

Burkholderia gladioli pv. cocovenenans is a serious safety issue in black fungus due to the deadly toxin, bongkrekic acid. This has triggered the demand for an efficient toxigenic phenotype recognition method. The objective of this study is to develop an efficient method for the recognition of toxin-producing B. gladioli strains. The potential of multilocus sequence typing and a back propagation neural network for the recognition of toxigenic B. cocovenenans was explored for the first time. The virulent strains were isolated from a black fungus cultivation environment in Qinba Mountain area, Shaanxi, China. A comprehensive evaluation of toxigenic capability of 26 isolates were conducted using Ultra Performance Liquid Chromatography for determination of bongkrekic acid and toxoflavin production in different culturing conditions and foods. The isolates produced bongkrekic acid in the range of 0.05-6.24 mg/L in black fungus and a highly toxin-producing strain generated 201.86 mg/L bongkrekic acid and 45.26 mg/L toxoflavin in co-cultivation with Rhizopus oryzae on PDA medium. Multilocus sequence typing phylogeny (MLST) analysis showed that housekeeping gene sequences have a certain relationship with a strain toxigenic phenotype. We developed a well-trained, back-propagation neutral network for prediction of toxigenic phenotype in B. gladioli based on MLST sequences with an accuracy of 100% in the training set and an accuracy of 86.7% in external test set strains. The BP neutral network offers a highly efficient approach to predict toxigenic phenotype of strains and contributes to hazard detection and safety surveillance.

Genomic analysis of Citrobacter from Australian wastewater and silver gulls reveals novel sequence types carrying critically important antibiotic resistance genes

Antimicrobial resistance (AMR) is a major public health concern, and environmental bacteria have been recognized as important reservoirs of antimicrobial resistance genes (ARGs). Citrobacter, a common environmental bacterium and opportunistic pathogen in humans and other animals, has been largely understudied in terms of its diversity and AMR potential. Whole-genome (short-read) sequencing on a total of 77 Citrobacter isolates obtained from Australian silver gull (Chroicocephalus novaehollandiae) (n = 17) and influent wastewater samples (n = 60) was performed, revealing a diverse Citrobacter population, with seven different species and 33 sequence types, 17 of which were novel. From silver gull using non-selective media we isolated a broader range of species with little to no mobilised ARG carriage. Wastewater isolates (selected using Carbapenem- Resistant Enterobacterales (CRE) selective media) carried a heavy burden of ARGs (up to 21 ARGs, conferring resistance to nine classes of antibiotics), with several novel multidrug-resistant (MDR) lineages identified, including C. braakii ST1110, which carried ARGs conferring resistance to eight to nine classes of antibiotics, and C. freundii ST1105, which carried two carbapenemase genes, blaIMP-4 in class 1 integron structure, and blaKPC-2. Additionally, we identified an MDR C. portucalensis isolate carrying blaNDM-1, blaSHV-12, and mcr-9. We identified IncC, IncM2, and IncP6 plasmids as the likely vectors for many of the critically important mobilised ARGs. Phylogenetic analyses were performed to assess any epidemiological linkages between isolation sources, demonstrating low relatedness across sources beyond the ST level. However, these analyses did reveal some closer relationships between strains from disparate wastewater sources despite their collection some 13,000 km apart. These findings support the need for future surveillance of Citrobacter populations in wastewater and wildlife populations to monitor for potential opportunistic human pathogens.

Integrated analysis of patient networks and plasmid genomes reveals a regional, multi-species outbreak of carbapenemase-producing Enterobacterales carrying both blaIMP and mcr-9 genes

BACKGROUND

Carbapenemase-producing Enterobacterales (CPE) are challenging in healthcare, with resistance to multiple classes of antibiotics. This study describes the emergence of IMP-encoding CPE amongst diverse Enterobacterales species between 2016 and 2019 across a London regional network.

METHODS

We performed a network analysis of patient pathways, using electronic health records, to identify contacts between IMP-encoding CPE positive patients. Genomes of IMP-encoding CPE isolates were overlayed with patient contacts to imply potential transmission events.

RESULTS

Genomic analysis of 84 Enterobacterales isolates revealed diverse species (predominantly Klebsiella spp, Enterobacter spp, E. coli); 86% (72/84) harboured an IncHI2 plasmid carrying blaIMP and colistin resistance gene mcr-9 (68/72). Phylogenetic analysis of IncHI2 plasmids identified three lineages showing significant association with patient contacts and movements between four hospital sites and across medical specialities, which was missed on initial investigations.

CONCLUSIONS

Combined, our patient network and plasmid analyses demonstrate an interspecies, plasmid-mediated outbreak of blaIMPCPE, which remained unidentified during standard investigations. With DNA sequencing and multi-modal data incorporation, the outbreak investigation approach proposed here provides a framework for real-time identification of key factors causing pathogen spread. Plasmid-level outbreak analysis reveals that resistance spread may be wider than suspected, allowing more interventions to stop transmission within hospital networks.

High-risk lineages of extended spectrum cephalosporinase producing Escherichia coli from Eurasian griffon vultures (Gyps fulvus) foraging in landfills in north-eastern Spain

Extended-spectrum cephalosporinase producing (ESC) E. coli are regarded as key indicator microorganisms of antimicrobial resistance (AMR), calling for a One Health integrated global surveillance strategy. Wildlife is exposed to antibiotic contaminants and/or resistant bacteria that have been released into the environment, potentially acting as reservoirs and spreaders of resistance genes as well as sentinels of anthropogenic pressure. Monitoring AMR in wildlife has become crucial in determining anthropogenic environmental impacts as well as transmission routes. In this study, we determined the occurrence and potential sources of ESC E. coli in 218 Eurasian griffon vultures (Gyps fulvus) foraging regularly on human waste disposed at a dumpsite in north-eastern Spain. Minimal inhibitory concentration for 14 different antimicrobials was performed to evaluate the phenotype of the isolates, and whole genome sequencing was carried out to investigate lineages and plasmids harbouring ESC genes. Our sequences were compared to previously published Spanish sequences of human, animal, and wildlife origin. We report a high prevalence of CTX-M-15, as well as the presence of other resistance genes such as OXA-10, CTX-M-27, and CTX-M-65 which are rarely described in European livestock, suggesting a human origin. The isolates also carried a diverse range of additional AMR genes for a broad spectrum of drug families, with the majority being multi-drug resistant. The phylogenomic analyses suggests the transmission of high-risk lineages from humans to vultures, with 49 % of our isolates matching the most common extraintestinal pathogenic E. coli (ExPEC) lineages described in humans worldwide, including ST131, ST10 and ST58. We conclude that anthropogenically altered habitats, such as landfills, are hotspots for the acquisition and spread of high-risk ESC E. coli lineages associated with hospital infections. Measures must be implemented to limit their spread into natural environments.

Streptococcus suis Research Update: Serotype Prevalence and Antimicrobial Resistance Distribution in Swine Isolates Recovered in Spain from 2020 to 2022

This study aimed to update the Streptococcus suis serotype distribution in Spain by analysing 302 clinical isolates recovered from diseased pigs between 2020 and 2022. The main objectives were to identify prevalent serotypes, differentiate specific serotypes 1, 14, 2, and 1/2, investigate specific genotypic and phenotypic antimicrobial resistance features, and explore associations between resistance genes and phenotypic resistances. Serotypes 9 (21.2%), 1 (16.2%), 2 (15.6%), 3 (6%), and 7 (5.6%) were the most prevalent, whereas serotypes 14 and 1/2 corresponded with 4.3% and 0.7% of all isolates. Antimicrobial resistance genes, including tet(O), erm(B), lnu(B), lsa(E), tet(M), and mef(A/E), were analysed, which were present in 85.8%, 65.2%, 7%, 7%, 6.3%, and 1% of the samples, respectively. Susceptibility testing for 18 antimicrobials revealed high resistance levels, particularly for clindamycin (88.4%), chlortetracycline (89.4%), and sulfadimethoxine (94.4%). Notably, seven significant associations (p < 0.0001) were detected, correlating specific antimicrobial resistance genes to the observed phenotypic resistance. These findings contribute to understanding the S. suis serotype distribution and its antibiotic resistance profiles in Spain, offering valuable insights for veterinary and public health efforts in managing S. suis-associated infections.

Serological and molecular detection as well as typing of Leptospira spp. in foxes, raccoons, and other wild carnivores in North-Eastern Germany, 2021-2022

Leptospirosis is a worldwide zoonosis caused by pathogenic Leptospira spp. While the latter are reported from various mammal hosts such as humans, dogs, or rodents, less is known about their presence in wild carnivores. We therefore investigated the presence of Leptospira spp. in foxes, raccoons, badgers, raccoon dogs, and martens in North-Eastern Germany. Kidney, urine, and blood specimens obtained from legally hunted or road-killed animals were tested by real-time PCR and by serogroup specific antibody detection for the presence of Leptospira spp. Additionally, kidney and urine specimens were tested by real-time PCR for the presence of Brucella spp. and Francisella tularensis, with all being negative for these two zoonotic pathogens. Leptospira spp. were detected by PCR in 12.6 % (n = 21/166) and serologically in 26.2 % (n = 53/202) of tissue and serum samples, respectively. Antibodies to 15 different serogroups were identified with Javanica (n = 25) and Bataviae (n = 12) being predominant. A high sero-prevalence of 34.0 % and 18.6 % in foxes and raccoons, respectively, and the presence of ST17 associated with human and animal leptospirosis indicates a reservoir and the zoonotic potential of these wild animals.

Streptococcus pneumoniae serotype 33G: genetic, serological, and structural analysis of a new capsule type

IMPORTANCE

Streptococcus pneumoniae (the pneumococcus) is a bacterial pathogen with the greatest burden of disease in Asia and Africa. The pneumococcal capsular polysaccharide has biological relevance as a major virulence factor as well as public health importance as it is the target for currently licensed vaccines. These vaccines have limited valency, covering up to 23 of the >100 known capsular types (serotypes) with higher valency vaccines in development. Here, we have characterized a new pneumococcal serotype, which we have named 33G. We detected serotype 33G in nasopharyngeal swabs (n = 20) from children and adults hospitalized with pneumonia, as well as healthy children in Mongolia. We show that the genetic, serological, and biochemical properties of 33G differ from existing serotypes, satisfying the criteria to be designated as a new serotype. Future studies should focus on the geographical distribution of 33G and any changes in prevalence following vaccine introduction.

Successful treatment of streptococcal toxic shock syndrome complicated by primary peritonitis and bilateral empyema in a healthy young woman: Identification of uncommon clone emm103 and novel sequence type 1363

Streptococcal toxic shock syndrome (STSS) has a dramatic clinical course and high mortality rate. Here, we report a case of STSS complicated by primary peritonitis and bilateral empyema. A previously healthy young woman was diagnosed with STSS complicated by primary peritonitis and bilateral empyema. Blood culture results on admission were negative. Sever shock, respiratory failure, systemic inflammation, thrombocytopenia, renal failure, ascites, and pleural effusion occurred, mimicking thrombocytopenia, anasarca, fever, reticulin fibrosis/renal failure and organomegaly (TAFRO) syndrome. Retesting blood cultures identified Streptococcus pyogenes. Gram staining of ascites and pleural fluid indicated gram-positive cocci in chains. Antibiotics, immunoglobulins, and surgical intervention led to recovery without complications. Ex-post genotypic analyses showed uncommon emm103.0 (cluster E3) of emm long sequence (784 base) and novel sequence type 1363. STSS diagnosis can be difficult as it mimics other systemic inflammatory diseases. Therefore, it is crucial for clinicians to perform microbiological examinations from infection foci, even if the initial culture is negative.

Antibiotic resistance, plasmids, and virulence-associated markers in human strains of Campylobacter jejuni and Campylobacter coli isolated in Italy

Campylobacteriosis, a prevalent foodborne gastrointestinal infection in Europe, is primarily caused by Campylobacter jejuni and Campylobacter coli, with rising global concerns over antimicrobial resistance in these species. This study comprehensively investigates 133 human-origin Campylobacter spp. strains (102 C. jejuni and 31 C. coli) collected in Italy from 2013 to 2021. The predominant Multilocus Sequence Typing Clonal complexes (CCs) were ST-21 CC and ST-206 CC in C. jejuni and ST-828 CC in C. coli. Ciprofloxacin and tetracycline resistance, mainly attributed to GyrA (T86I) mutation and tet(O) presence, were prevalent, while erythromycin resistance was associated with 23S rRNA gene mutation (A2075G), particularly in C. coli exhibiting multidrug-resistant pattern CipTE. Notable disparities in virulence factors among strains were observed, with C. jejuni exhibiting a higher abundance compared to C. coli. Notably, specific C. jejuni sequence types, including ST-21, ST-5018, and ST-1263, demonstrated significantly elevated counts of virulence genes. This finding underscores the significance of considering both the species and strain-level variations in virulence factor profiles, shedding light on potential differences in the pathogenicity and clinical outcomes associated with distinct C. jejuni lineages. Campylobacter spp. plasmids were classified into three groups comprising pVir-like and pTet-like plasmids families, exhibiting diversity among Campylobacter spp. The study underscores the importance of early detection through Whole Genome Sequencing to identify potential emergent virulence, resistance/virulence plasmids, and new antimicrobial resistance markers. This approach provides actionable public health data, supporting the development of robust surveillance programs in Italy.

Investigating environmental transmission to resolve a Bacillus cereus group outbreak in a neonatal intensive care unit using core genome multilocus sequence typing

BACKGROUND

We analyzed an outbreak of Bacillus cereus group (Bcg) at a single-center neonatal intensive care unit level IV by conducting comprehensive sampling of both patients and the environment.

METHODS

Between 06/2020 and 10/2021, all Bcg isolates identified by both regular colonization screening and additional sampling of the environment were subjected to whole-genome sequencing, followed by in vitro extraction of MLST ST, resistance genes and virulence factors. Using publicly available genome sequences, we defined an ad hoc core genome multilocus sequence typing (cgMLST) scheme comprising 2759 target genes for Bcg typing, which we applied to the detected isolates. We have compared the results with a stable cgMLST that was published in the meantime and completed the investigation with a SNP analysis.

RESULTS

We analyzed 28 Bcg isolates from patient and environmental samples using MLST and cgMLST. This revealed multiple sequence types, with ST127 being the most common (n = 13). Both cgMLST schemes grouped ten of the 13 ST127 isolates into a cluster, including two invasive isolates from two different patients and several environmental samples. SNP analysis postulated a screen from a ventilation machine as a possible reservoir.

CONCLUSION

In sensitive settings such as neonatal intensive care units, considering the environment in outbreak analyses is crucial, especially when investigating potential transmission routes through shared devices. When dealing with widespread bacteria such as Bcg, high-resolution typing techniques are necessary. In this study, we successfully resolved an outbreak of Bcg infections using a custom cgMLST scheme combined with a SNP analysis.

gcPathogen: a comprehensive genomic resource of human pathogens for public health

Here, we present the manually curated Global Catalogue of Pathogens (gcPathogen), an extensive genomic resource designed to facilitate rapid and accurate pathogen analysis, epidemiological exploration and monitoring of antibiotic resistance features and virulence factors. The catalogue seamlessly integrates and analyzes genomic data and associated metadata for human pathogens isolated from infected patients, animal hosts, food and the environment. The pathogen list is supported by evidence from medical or government pathogenic lists and publications. The current version of gcPathogen boasts an impressive collection of 1 164 974 assemblies comprising 986 044 strains from 497 bacterial taxa, 4794 assemblies encompassing 4319 strains from 265 fungal taxa, 89 965 assemblies featuring 13 687 strains from 222 viral taxa, and 646 assemblies including 387 strains from 159 parasitic taxa. Through this database, researchers gain access to a comprehensive 'one-stop shop' that facilitates global, long-term public health surveillance while enabling in-depth analysis of genomes, sequence types, antibiotic resistance genes, virulence factors and mobile genetic elements across different countries, diseases and hosts. To access and explore the data and statistics, an interactive web interface has been developed, which can be accessed at https://nmdc.cn/gcpathogen/. This user-friendly platform allows seamless querying and exploration of the extensive information housed within the gcPathogen database.

Convergence of resistance and evolutionary responses in Escherichia coli and Salmonella enterica co-inhabiting chicken farms in China

Sharing of genetic elements among different pathogens and commensals inhabiting same hosts and environments has significant implications for antimicrobial resistance (AMR), especially in settings with high antimicrobial exposure. We analysed 661 Escherichia coli and Salmonella enterica isolates collected within and across hosts and environments, in 10 Chinese chicken farms over 2.5 years using data-mining methods. Most isolates within same hosts possessed the same clinically relevant AMR-carrying mobile genetic elements (plasmids: 70.6%, transposons: 78%), which also showed recent common evolution. Supervised machine learning classifiers revealed known and novel AMR-associated mutations and genes underlying resistance to 28 antimicrobials, primarily associated with resistance in E. coli and susceptibility in S. enterica. Many were essential and affected same metabolic processes in both species, albeit with varying degrees of phylogenetic penetration. Multi-modal strategies are crucial to investigate the interplay of mobilome, resistance and metabolism in cohabiting bacteria, especially in ecological settings where community-driven resistance selection occurs.

Molecular Characterization of Chimeric Staphylococcus aureus Strains from Waterfowl

Staphylococcus aureus is a versatile pathogen that does not only occur in humans but also in various wild and domestic animals, including several avian species. When characterizing S. aureus isolates from waterfowl, isolates were identified as atypical CC133 by DNA microarray analysis. They differed from previously sequenced CC133 strains in the presence of the collagen adhesin gene cna; some also showed a different capsule type and a deviant spa type. Thus, they were subjected to whole-genome sequencing. This revealed multiple insertions of large regions of DNA from other S. aureus lineages into a CC133-derived backbone genome. Three distinct strains were identified based on the size and extent of these inserts. One strain comprised two small inserts of foreign DNA up- and downstream of oriC; one of about 7000 nt or 0.25% originated from CC692 and the other, at ca. 38,000 nt or 1.3% slightly larger one was of CC522 provenance. The second strain carried a larger CC692 insert (nearly 257,000 nt or 10% of the strain's genome), and its CC522-derived insert was also larger, at about 53,500 nt or 2% of the genome). The third strain carried an identical CC692-derived region (in which the same mutations were observed as in the second strain), but it had a considerably larger CC522-like insertion of about 167,000 nt or 5.9% of the genome. Both isolates of the first, and two out of four isolates of the second strain also harbored a hemolysin-beta-integrating prophage carrying "bird-specific" virulence factors, ornithine cyclodeaminase D0K6J8 and a putative protease D0K6J9. Furthermore, isolates had two different variants of SCC elements that lacked mecA/mecC genes. These findings highlight the role of horizontal gene transfer in the evolution of S. aureus facilitated by SCC elements, by phages, and by a yet undescribed mechanism for large-scale exchange of core genomic DNA.

Genomic analysis reveals the population structure and antimicrobial resistance of avian Pasteurella multocida in China

OBJECTIVES

To investigate the population structure and antimicrobial resistance (AMR) of avian Pasteurella multocida in China.

METHODS

Utilizing WGS analysis, we explored the phylogeny using a dataset of 546 genomes, comprising avian P. multocida isolates from China (n = 121), the USA (n = 165), Australia(n = 153), Bangladesh (n = 3) and isolates of other hosts from China (n = 104). We examined the integrative and conjugative element (ICE) structures and the distribution of their components carrying resistance genes, and reconstructed the evolutionary history of A:L1:ST129 (n = 110).

RESULTS

The population structure of avian P. multocida in China was dominated by the A:L1:ST129 clone with limited genetic diversity. A:L1:ST129 isolates possessed a broader spectrum of resistance genes at comparatively higher frequencies than those from other hosts and countries. The novel putative ICEs harboured complex resistant clusters that were prevalent in A:L1:ST129. Bayesian analysis predicted that the A:L1:ST129 clone emerged around 1923, and evolved slowly.

CONCLUSIONS

A:L1:ST129 appears to possess a host predilection towards avian species in China, posing a potential health threat to other animals. The complex AMR determinants coupled with high frequencies may strengthen the population dominance of A:L1:ST129. The extensive antimicrobial utilization in poultry farming and the mixed rearing practices could have accelerated AMR accumulation in A:L1:ST129. ICEs, together with their resistant clusters, significantly contribute to resistance gene transfer and facilitate the adaptation of A:L1:ST129 to ecological niches. Despite the genetic stability and slow evolution rate, A:L1:ST129 deserves continued monitoring due to its propensity to retain resistance genes, warranting global attention to preclude substantial economic losses.

Prevalence, antibiotic resistance, and genomic characterisation of Campylobacter spp. in retail chicken in Hanoi, Vietnam

Campylobacter spp. are a leading cause of bacterial foodborne zoonosis worldwide, with poultry meat and products recognised as a significant source of human infection. In Vietnam there are few data regarding the occurrence, antimicrobial resistance, and genomic diversity of Campylobacter in poultry and poultry meat. The aim of this study was to estimate the prevalence of Campylobacter in chicken meat at retail in Hanoi, determine antimicrobial sensitivities of the Campylobacter isolated, and assess their genetic diversity. A total of 120 chicken meat samples were collected from eight traditional retail markets (n=80) and four supermarkets (n=40). Campylobacter was isolated following ISO 10272-1 : 2017 and identification verified by PCR. The prevalence of Campylobacter was 38.3 % (46/120) and C. coli was the most prevalent species in both retail markets (74 %) and supermarkets (88 %). The minimum inhibitory concentrations for ciprofloxacin, erythromycin, gentamicin, nalidixic acid, streptomycin, and tetracycline were determined by broth microdilution for 32 isolates. All characterised Campylobacter were resistant to ciprofloxacin, nalidixic acid, and tetracycline, with corresponding resistance determinants detected in the sequenced genomes. Most C. coli were multidrug resistant (24/28) and two harboured the erythromycin resistance gene ermB on a multiple drug-resistance genomic island, a potential mechanism for dissemination of resistance. The 32 isolates belonged to clonal complexes associated with both poultry and people, such as CC828 for C. coli. These results contribute to the One Health approach for addressing Campylobacter in Vietnam by providing detailed new insights into a main source of human infection and can inform the design of future surveillance approaches.

Wolbachia populations across organs of individual Culex pipiens: highly conserved intra-individual core pangenome with inter-individual polymorphisms

Wolbachia is a maternally inherited intracellular bacterium that infects a wide range of arthropods including mosquitoes. The endosymbiont is widely used in biocontrol strategies due to its capacity to modulate arthropod reproduction and limit pathogen transmission. Wolbachia infections in Culex spp. are generally assumed to be monoclonal but the potential presence of genetically distinct Wolbachia subpopulations within and between individual organs has not been investigated using whole genome sequencing. Here we reconstructed Wolbachia genomes from ovary and midgut metagenomes of single naturally infected Culex pipiens mosquitoes from Southern France to investigate patterns of intra- and inter-individual differences across mosquito organs. Our analyses revealed a remarkable degree of intra-individual conservancy among Wolbachia genomes from distinct organs of the same mosquito both at the level of gene presence-absence signal and single-nucleotide polymorphisms (SNPs). Yet, we identified several synonymous and non-synonymous substitutions between individuals, demonstrating the presence of some level of genomic heterogeneity among Wolbachia that infect the same C. pipiens field population. Overall, the absence of genetic heterogeneity within Wolbachia populations in a single individual confirms the presence of a dominant Wolbachia that is maintained under strong purifying forces of evolution.

Analysis of Whole Genome Sequencing Data for Detection of Antimicrobial Resistance Determinants

Genomic sequencing has revolutionized microbial typing methods and transformed high-throughput methods in reference, clinical, and research laboratories. The detection of antimicrobial-resistant (AMR) determinants using genomic methods can provide valuable information on the emergence of resistance. Here we describe an approach to detecting AMR determinants using an open access and freely available platform which does not require bioinformatic expertise.

Step-by-Step Bacterial Genome Comparison

Thanks to advancements in genome sequencing and bioinformatics, thousands of bacterial genome sequences are available in public databases. This presents an opportunity to study bacterial diversity in unprecedented detail. This chapter describes a complete bioinformatics workflow for comparative genomics of bacterial genomes, including genome annotation, pangenome reconstruction and visualization, phylogenetic analysis, and identification of sequences of interest such as antimicrobial-resistance genes, virulence factors, and phage sequences. The workflow uses state-of-the-art, open-source tools. The workflow is presented by means of a comparative analysis of Salmonella enterica serovar Typhimurium genomes. The workflow is based on Linux commands and scripts, and result visualization relies on the R environment. The chapter provides a step-by-step protocol that researchers with basic expertise in bioinformatics can easily follow to conduct investigations on their own genome datasets.

Computational structure-based approach to study chimeric antigens using a new protein scaffold displaying foreign epitopes

The identification and recombinant production of functional antigens and/or epitopes of pathogens represent a crucial step for the development of an effective protein-based vaccine. Many vaccine targets are outer membrane proteins anchored into the lipidic bilayer through an extended hydrophobic portion making their recombinant production challenging. Moreover, only the extracellular loops, and not the hydrophobic regions, are naturally exposed to the immune system. In this work, the Domain 3 (D3) from Group B Streptococcus (GBS) pilus 2a backbone protein has been identified and engineered to be used as a scaffold for the display of extracellular loops of two Neisseria gonorrhoeae membrane proteins (PorB.1b and OpaB). A computational structure-based approach has been applied to the design of both the scaffold and the model antigens. Once identified the best D3 engineerable site, several different chimeric D3 displaying PorB.1b and OpaB extracellular loops were produced as soluble proteins. Each molecule has been characterized in terms of solubility, stability, and ability to correctly display the foreign epitope. This antigen dissection strategy allowed the identification of most immunogenic extracellular loops of both PorB.1b and OpaB gonococcal antigens. The crystal structure of chimeric D3 displaying PorB.1b immunodominant loop has been obtained confirming that the engineerization did not alter the predicted native structure of this epitope. Taken together, the reported data suggest that D3 is a novel protein scaffold for epitope insertion and display, and a valid alternative to the production of whole membrane protein antigens. Finally, this work describes a generalized computational structure-based approach for the identification, design, and dissection of epitopes in target antigens through chimeric proteins.

Genomic and clinical characteristics of campylobacteriosis in Australia

Campylobacter spp. are a common cause of bacterial gastroenteritis in Australia, primarily acquired from contaminated meat. We investigated the relationship between genomic virulence characteristics and the severity of campylobacteriosis, hospitalisation, and other host factors.We recruited 571 campylobacteriosis cases from three Australian states and territories (2018-2019). We collected demographic, health status, risk factors, and self-reported disease data. We whole genome sequenced 422 C. jejuni and 84 C. coli case isolates along with 616 retail meat isolates. We classified case illness severity using a modified Vesikari scoring system, performed phylogenomic analysis, and explored risk factors for hospitalisation and illness severity.On average, cases experienced a 7.5 day diarrhoeal illness with additional symptoms including stomach cramps (87.1 %), fever (75.6 %), and nausea (72.0 %). Cases aged ≥75 years had milder symptoms, lower Vesikari scores, and higher odds of hospitalisation compared to younger cases. Chronic gastrointestinal illnesses also increased odds of hospitalisation. We observed significant diversity among isolates, with 65 C. jejuni and 21 C. coli sequence types. Antimicrobial resistance genes were detected in 20.4 % of isolates, but multidrug resistance was rare (0.04 %). Key virulence genes such as cdtABC (C. jejuni) and cadF were prevalent (>90 % presence) but did not correlate with disease severity or hospitalisation. However, certain genes (e.g. fliK, Cj1136, and Cj1138) appeared to distinguish human C. jejuni cases from food source isolates.Campylobacteriosis generally presents similarly across cases, though some are more severe. Genotypic virulence factors identified in the literature to-date do not predict disease severity but may differentiate human C. jejuni cases from food source isolates. Host factors like age and comorbidities have a greater influence on health outcomes than virulence factors.

Broad Genome Sequencing of Environmental and Clinical Strains and Genotyping

The Next-Generation Sequencing revolution had a great impact on the genomics of Pseudomonas aeruginosa. Since the first release of the P. aeruginosa PAO1 genome, there are more than 5700 genomes published. This wealth of information has been accompanied by the development of bioinformatic tools for handling genomic and phenotypic data. Bioinformatics, indeed, become de facto a big data science. In this chapter, we give a brief historical overview of the knowledge gained from P. aeruginosa genome sequencing, then we describe the wet-lab procedure to extract the DNA and prepare the library for broad genome sequencing using Illumina MiSeq technology. As last, we describe three user-friendly bioinformatics procedures to infer the P. aeruginosa genotype, starting from NGS data, with the Multi-Locus Sequence Typing method, and visualize it as a minimum spanning tree.

Antimicrobial susceptibility of U.S. porcine Brachyspira isolates and genetic diversity of B. hyodysenteriae by multilocus sequence typing

Swine dysentery, caused by Brachyspira hyodysenteriae and the newly recognized Brachyspira hampsonii in grower-finisher pigs, is a substantial economic burden in many swine-rearing countries. Antimicrobial therapy is the only commercially available measure to control and prevent Brachyspira-related colitis. However, data on antimicrobial susceptibility trends and genetic diversity of Brachyspira species from North America is limited. We evaluated the antimicrobial susceptibility profiles of U.S. Brachyspira isolates recovered between 2013 and 2022 to tiamulin, tylvalosin, lincomycin, doxycycline, bacitracin, and tylosin. In addition, we performed multilocus sequence typing (MLST) on 64 B. hyodysenteriae isolates. Overall, no distinct alterations in the susceptibility patterns over time were observed among Brachyspira species. However, resistance to the commonly used antimicrobials was seen sporadically with a higher resistance frequency to tylosin compared to other tested drugs. B. hampsonii was more susceptible to the tested drugs than B. hyodysenteriae and B. pilosicoli. MLST revealed 16 different sequence types (STs) among the 64 B. hyodysenteriae isolates tested, of which 5 STs were previously known, whereas 11 were novel. Most isolates belonged to the known STs: ST93 (n = 32) and ST107 (n = 13). Our findings indicate an overall low prevalence of resistance to clinically important antimicrobials other than tylosin and bacitracin, and high genetic diversity among the clinical Brachyspira isolates from pigs in the United States during the past decade. Further molecular, epidemiologic, and surveillance studies are needed to better understand the infection dynamics of Brachyspira on swine farms and to help develop effective control measures.

pyngoST: fast, simultaneous and accurate multiple sequence typing of Neisseria gonorrhoeae genome collections

Extensive gonococcal surveillance has been performed using molecular typing at global, regional, national and local levels. The three main genotyping schemes for this pathogen, multi-locus sequence typing (MLST), Neisseria gonorrhoeae multi-antigen sequence typing (NG-MAST) and N. gonorrhoeae sequence typing for antimicrobial resistance (NG-STAR), allow inter-laboratory and inter-study comparability and reproducibility and provide an approximation to the gonococcal population structure. With whole-genome sequencing (WGS), we obtain a substantially higher and more accurate discrimination between strains compared to previous molecular typing schemes. However, WGS remains unavailable or not affordable in many laboratories, and thus bioinformatic tools that allow the integration of data among laboratories with and without access to WGS are imperative for a joint effort to increase our understanding of global pathogen threats. Here, we present pyngoST, a command-line Python tool for fast, simultaneous and accurate sequence typing of N. gonorrhoeae from WGS assemblies. pyngoST integrates MLST, NG-MAST and NG-STAR, and can also designate NG-STAR clonal complexes, NG-MAST genogroups and penA mosaicism, facilitating multiple sequence typing from large WGS assembly collections. Exact and closest matches for existing alleles and sequence types are reported. The implementation of a fast multi-pattern searching algorithm allows pyngoST to be rapid and report results on 500 WGS assemblies in under 1 min. The mapping of typing results on a core genome tree of 2375 gonococcal genomes revealed that NG-STAR is the scheme that best represents the population structure of this pathogen, emphasizing the role of antimicrobial use and antimicrobial resistance as a driver of gonococcal evolution. This article contains data hosted by Microreact.

The important lessons lurking in the history of meningococcal epidemiology

INTRODUCTION

The epidemiology of invasive meningococcal disease (IMD), a rare but potentially fatal illness, is typically described as unpredictable and subject to sporadic outbreaks.

AREAS COVERED

Meningococcal epidemiology and vaccine use during the last ~ 200 years are examined within the context of meningococcal characterization and classification to guide future IMD prevention efforts.

EXPERT OPINION

Historical and contemporary data highlight the dynamic nature of meningococcal epidemiology, with continued emergence of hyperinvasive clones and affected regions. Recent shifts include global increases in serogroup W disease, meningococcal antimicrobial resistance (AMR), and meningococcal urethritis; additionally, unvaccinated populations have experienced disease resurgences following lifting of COVID-19 restrictions. Despite these changes, a close analysis of meningococcal epidemiology indicates consistent dominance of serogroups A, B, C, W, and Y and elevated IMD rates among infants and young children, adolescents/young adults, and older adults. Demonstrably effective vaccines against all 5 major disease-causing serogroups are available, and their prophylactic use represents a powerful weapon against IMD, including AMR. The World Health Organization's goal of defeating meningitis by the year 2030 demands broad protection against IMD, which in turn indicates an urgent need to expand meningococcal vaccination programs across major disease-causing serogroups and age-related risk groups.

Beyond the Wild MRSA: Genetic Features and Phylogenomic Review of mecC-Mediated Methicillin Resistance in Non-aureus Staphylococci and Mammaliicocci

Methicillin resistance, mediated by the mecA gene in staphylococci and mammaliicocci, has caused tremendous setbacks in the use of antibiotics in human and veterinary medicine due to its high potential of presenting the multidrug resistance (MDR) phenotype. Three other mec analogs exist, of which the mecC has evolutionary been associated with methicillin-resistant Staphylococcus aureus (MRSA) in wild animals, thus loosely referred to as the wild MRSA. In this study, we present an epidemiological review and genomic analysis of non-aureus staphylococci and mammaliicocci that carry the mecC-mediated methicillin resistance trait and determine whether this trait has any relevant link with the One Health niches. All previous studies (2007 till 2023) that described the mecC gene in non-aureus staphylococci and mammaliicocci were obtained from bibliometric databases, reviewed, and systematically analyzed to obtain the antimicrobial resistance (AMR) and virulence determinants, mobilome, and other genetic contents. Moreover, core genome single-nucleotide polymorphism analysis was used to assess the relatedness of these strains. Of the 533 articles analyzed, only 16 studies (on livestock, environmental samples, milk bulk tanks, and wild animals) were eligible for inclusion, of which 17 genomes from 6 studies were used for various in silico genetic analyses. Findings from this systematic review show that all mecC-carrying non-aureus staphylococci were resistant to only beta-lactam antibiotics and associated with the classical SCCmec XI of S. aureusLGA251. Similarly, two studies on wild animals reported mecC-carrying Mammaliicoccus stepanovicii associated with SCCmec XI. Nevertheless, most of the mecC-carrying Mammaliicoccus species presented an MDR phenotype (including linezolid) and carried the SCCmec-mecC hybrid associated with mecA. The phylogenetic analysis of the 17 genomes revealed close relatedness (<20 SNPs) and potential transmission of M. sciuri and M. lentus strains in livestock farms in Algeria, Tunisia, and Brazil. Furthermore, closely related M. sciuri strains from Austria, Brazil, and Tunisia (<40 SNPs) were identified. This systematic review enhances our comprehension of the epidemiology and genetic organization of mecC within the non-aureus staphylococci and mammaliicocci. It could be hypothesized that the mecC-carrying non-aureus staphylococci are evolutionarily related to the wild MRSA-mecC. The potential implications of clonal development of a lineage of mecA/mecC carrying strains across multiple dairy farms in a vast geographical region with the dissemination of MDR phenotype is envisaged. It was observed that most mecC-carrying non-aureus staphylococci and mammaliicocci were reported in mastitis cases. Therefore, veterinarians and veterinary microbiology laboratories must remain vigilant regarding the potential existence of mecA/mecC strains originating from mastitis as a potential niche for this resistance trait.

Clonal Diversity, Antibiotic Resistance, and Virulence Factor Prevalence of Community Associated Staphylococcus aureus in Southeastern Virginia

Staphylococcus aureus is a significant human pathogen with a formidable propensity for antibiotic resistance. Worldwide, it is the leading cause of skin and soft tissue infections (SSTI), septic arthritis, osteomyelitis, and infective endocarditis originating from both community- and healthcare-associated settings. Although often grouped by methicillin resistance, both methicillin-resistant (MRSA) and methicillin-sensitive (MSSA) strains are known to cause significant pathologies and injuries. Virulence factors and growing resistance to antibiotics play major roles in the pathogenicity of community-associated strains. In our study, we examined the genetic variability and acquired antibiograms of 122 S. aureus clinical isolates from SSTI, blood, and urinary tract infections originating from pediatric patients within the southeast region of Virginia, USA. We identified a suite of clinically relevant virulence factors and evaluated their prevalence within these isolates. Five genes (clfA, spA, sbi, scpA, and vwb) with immune-evasive functions were identified in all isolates. MRSA isolates had a greater propensity to be resistant to more antibiotics as well as significantly more likely to carry several virulence factors compared to MSSA strains. Further, the carriage of various genes was found to vary significantly based on the infection type (SSTI, blood, urine).

Genomics insights into a Mycobacterium pinnipedii isolate causing tuberculosis in a captive South American sea lion (Otaria flavescens) from Italy

Tuberculosis (TB) affects humans and other animals, and it is caused by bacteria within the Mycobacterium tuberculosis complex (MTBC). In this study, we report the characterisation of Mycobacterium pinnipedii that caused a TB case in a sea lion (Otaria flavescens) kept in an Italian zoo. The animal died due to severe, progressive disorders involving the respiratory and gastro-enteric systems and the skin. At necropsy, typical gross lesions referable to a TB generalised form were found. In particular, nodular granulomatous lesions were detected in the lungs and several lymph nodes, and colonies referable to Mycobacterium spp. were isolated from lung, mesenteric, and mediastinal lymph nodes. The isolate was identified by PCR as a MTBC, had a spoligotype SB 1480 ("seal lineage"), and was characterised and characterised by whole-genome sequencing analysis confirming that the MTBC involved was M. pinnipedii. The analysis of the resistome and virulome indicated the presence of macrolide and aminoglycoside resistance genes intrinsic in M. tuberculosis [erm-37 and aac(2')-Ic] and confirmed the presence of the region of difference 1 (RD1), harbouring the esxA and esxB virulence genes, differently from its closest taxon, M. microti. As for other MTCB members, M. pinnipedii infection can spill over into non-pinniped mammalian species; therefore, zoological gardens, veterinary practitioners, and public health officers should be aware of the hazard posed by tuberculosis from marine mammals. Since the isolate under study, as well as all available genomes of M. pinnipedii investigated in this study retains almost all the M. tuberculosis virulence genes, it could indeed cause infection, lesions, and disease in other animal species, including humans.

Seroprevalence and Molecular Characterization of Brucella abortus from the Himalayan Marmot in Qinghai, China

Objective

Brucellosis is a serious public health issue in Qinghai (QH), China. Surveying the seroprevalence and isolation of B. abortus strains from marmots is key to understanding the role of wildlife in the maintenance and spread of brucellosis.

Methods

In this study, a set of methods, including a serology survey, bacteriology, antibiotic susceptibility, molecular genotyping (MLST and MLVA), and genome sequencing, were employed to characterize the two B. abortus strains.

Results

The seroprevalence of brucellosis in marmots was 7.0% (80/1146) by serum tube agglutination test (SAT); one Brucella strain was recovered from these positive samples, and another Brucella strain from a human. Two strains were identified as B. abortus bv. 1 and were susceptible to all eight drugs examined. The distribution patterns of the accessory genes, virulence associated genes, and resistance genes of the two strains were consistent, and there was excellent collinearity between the two strains on chromosome I, but they had significant SVs in chromosome II, including inversions and translocations. MLST genotyping identified two B. abortus strains as ST2, and MLVA-16 analysis showed that the two strains clustered with strains from northern China. WGS-SNP phylogenetic analysis showed that the strains were genetically homogeneous with strains from the northern region, implying that strains from a common lineage were spread continuously in different regions and hosts.

Conclusion

Seroprevalence and molecular clues demonstrated frequent direct or indirect contact between sheep/goats, cattle, and marmots, implying that wildlife plays a vital role in the maintenance and spread of B. abortus in the Qinghai-Tibet Plateau.

Molecular characterization of Listeria monocytogenes strains isolated from imported food in China from 14 countries/regions, 2003-2018

Listeria monocytogenes (Lm) is associated with severe foodborne infections and ubiquitous in the nature. Identification of characteristics of Lm transmission through trading of food products is essential for rapidly tracking Lm sources and controlling dissemination of listeriosis. In this study, a total of 44 Lm strains were isolated from food products originating from 14 countries/regions during 2003-2018 at the Shanghai port. The genomes of these Lm strains were sequenced by high-throughput sequencing. Multilocus sequence typing (MLST) analysis showed that 43 isolates were divided into 17 sequence types (STs). The distribution of STs was decentralized, with the dominant ST2 accounting for only 18.18% of the strains. The LM63 strain did not match with any of the existing STs. Core-genome MLST (cgMLST) analysis based on 1748 core genes categorized the 44 strains into 30 cgMLST types (CTs), with CT10153 and CT7892 as the most predominant CTs. Notably, LM63 and LM67 shared the same CT in the cgMLST analysis. The phylogenetic analysis based on single-copy homologous genes revealed that the 44 Lm strains were primarily classified into two lineages. The SNP analysis also indicated that these strains were roughly divided into two clades, with strains in the first clade mainly collected earlier than those in the second clade, which were predominantly collected from 2010 onwards. The analysis using the virulence factor database (VFDB) indicated that the virulence gene inlJ was the most prevalent among these 44 strains. Notably, ddrA, msbA, and sugC were enriched in this dataset, requiring further clarification of their roles in Listeria through future studies. These results might provide a clue for understanding of the global epidemiology and surveillance of Lm and present insights for implementing effective measures to reduce or prevent Listeria contamination outbreaks in imported food products.

Evolution and exchange of plasmids in pathogenic Neisseria

IMPORTANCE

Horizontal gene transfer (HGT) is a major influence in driving the spread of antimicrobial resistance (AMR) in many bacteria. A conjugative plasmid which is widespread in Neisseria gonorrhoeae, pConj, prevented the use of tetracycline/doxycycline for treating gonococcal infection. Here, we show that pConj evolved in the related pathogen, Neisseria meningitidis, and has been acquired by the gonococcus from the meningococcus on multiple occasions. Following its initial acquisition, pConj spread to different gonococcal lineages; changes in the plasmid's conjugation machinery associated with another transfer event limit spread in the gonococcal populations. Our findings have important implications for the use of doxycycline to prevent bacterial sexually transmitted disease which is likely to exacerbate the spread of AMR through HGT in pathogenic bacteria.

Bactericidal killing of meningococcal W strains isolated in Argentina by the sera of adolescents and infants immunized with 4-component meningococcal serogroup B vaccine (4CMenB)

Invasive meningococcal disease (IMD) is a life-threatening disease caused by meningococcal serogroups A, B, C, W, X, and Y, of which B and W are most common in Argentina. The 4-component meningococcal serogroup B (4CMenB) vaccine contains three purified recombinant protein antigens (Neisseria adhesin A [NadA], factor H binding protein [fHbp], and Neisserial Heparin Binding Antigen [NHBA]) and outer membrane vesicles (OMV), which is derived from the New Zealand epidemic strain and contains Porin A 1.4. These antigens are present and conserved in strains that belong to other serogroups. In this study, we show that 10/11 (91%) meningococcal serogroup W (MenW) strains selected to be representative of MenW isolates that caused IMD in Argentina during 2010-2011 were killed in bactericidal assays by the sera of adolescents and infants who had been immunized with the 4CMenB vaccine. We also show that MenW strains that caused IMD in Argentina during 2018-2021 were genetically similar to the earlier strains, indicating that the 4CMenB vaccine would likely still provide protection against current MenW strains. These data highlight the potential of 4CMenB vaccination to protect adolescents and infants against MenW strains that are endemic in Argentina.

Complete genome sequence of a Clostridioides difficile cryptic C-III strain isolated from horse feces

We provide the complete genome of a non-toxigenic Clostridioides difficile strain isolated from horse feces. The strain represents a sub-cluster in the cryptic clade C-III. The genome consists of one chromosome (4,144,784 bp) and one plasmid (10,144 bp) and encodes 3,798 putative genes.

Plasmid genomic epidemiology of blaKPC carbapenemase-producing Enterobacterales in Canada, 2010-2021

Carbapenems are considered last-resort antibiotics for the treatment of infections caused by multidrug-resistant Enterobacterales, but carbapenem resistance due to acquisition of carbapenemase genes is a growing threat that has been reported worldwide. Klebsiella pneumoniae carbapenemase (blaKPC) is the most common type of carbapenemase in Canada and elsewhere; it can hydrolyze penicillins, cephalosporins, aztreonam, and carbapenems and is frequently found on mobile plasmids in the Tn4401 transposon. This means that alongside clonal expansion, blaKPC can disseminate through plasmid- and transposon-mediated horizontal gene transfer. We applied whole genome sequencing to characterize the molecular epidemiology of 829 blaKPC carbapenemase-producing isolates collected by the Canadian Nosocomial Infection Surveillance Program from 2010 to 2021. Using a combination of short-read and long-read sequencing, we obtained 202 complete and circular blaKPC-encoding plasmids. Using MOB-suite, 10 major plasmid clusters were identified from this data set which represented 87% (175/202) of the Canadian blaKPC-encoding plasmids. We further estimated the genomic location of incomplete blaKPC-encoding contigs and predicted a plasmid cluster for 95% (603/635) of these. We identified different patterns of carbapenemase mobilization across Canada related to different plasmid clusters, including clonal transmission of IncF-type plasmids (108/829, 13%) in K. pneumoniae clonal complex 258 and novel repE(pEh60-7) plasmids (44/829, 5%) in Enterobacter hormaechei ST316, and horizontal transmission of IncL/M (142/829, 17%) and IncN-type plasmids (149/829, 18%) across multiple genera. Our findings highlight the diversity of blaKPC genomic loci and indicate that multiple, distinct plasmid clusters have contributed to blaKPC spread and persistence in Canada.

Complete genomes of Pasteurella multocida isolates from porcine, bovine, and cervine samples found in New Zealand farms

Here, we present complete genome assemblies of Pasteurella multocida strains isolated from porcine, bovine, and cervine farms as part of bacteriology incursion investigations to identify pathogens that might present a sanitary risk to New Zealand.

The population-level impact of Enterococcus faecalis genetics on intestinal colonization and extraintestinal infection

IMPORTANCE

Enterococcus faecalis causes life-threatening invasive hospital- and community-associated infections that are usually associated with multidrug resistance globally. Although E. faecalis infections cause opportunistic infections typically associated with antibiotic use, immunocompromised immune status, and other factors, they also possess an arsenal of virulence factors crucial for their pathogenicity. Despite this, the relative contribution of these virulence factors and other genetic changes to the pathogenicity of E. faecalis strains remain poorly understood. Here, we investigated whether specific genomic changes in the genome of E. faecalis isolates influence its pathogenicity-infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection and intestinal colonization. Our findings indicate that E. faecalis genetics partially influence the infection of hospitalized and nonhospitalized individuals and the propensity to cause extraintestinal infection, possibly due to gut-to-bloodstream translocation, highlighting the potential substantial role of host and environmental factors, including gut microbiota, on the opportunistic pathogenic lifestyle of this bacterium.

Emergence of novel ST1299 vanA lineages as possible cause for the striking rise of vancomycin resistance among invasive strains of Enterococcus faecium at a German university hospital

IMPORTANCE

The proportion of VREfm among all Enterococcus faecium isolated from blood cultures in German hospitals has increased in the period 2015-2020 from 11.9% to 22.3% with a country-wide spread of the clonal lineage ST117/CT71 vanB. In this study, we provided useful information about the genetic diversity of invasive strains of E. faecium. Moreover, our findings confirm the nosocomial spread of novel ST1299 vanA lineages, which recently had a rapid expansion in Austria and the south-eastern part of Germany.

NDM-1- and OXA-23-producing Acinetobacter baumannii in wastewater of a Nigerian hospital

IMPORTANCE

Acinetobacter baumannii is a leading cause of hospital-associated infections globally. A. baumannii reservoirs outside hospital settings are still unknown, and their occurrence in the environment is linked to clinical and anthropogenic activities. Although the risk of transmission of A. baumannii from environmental sources to humans is not fully understood, these sources pose significant risks for the continued dissemination of A. baumannii and their resistance traits. This study provides evidence that diverse and clinically relevant A. baumannii strains, many of which are resistant to carbapenems, are constantly being discharged into the environment through inadequately treated hospital wastewater. We further elucidate potential transmission routes between the environment and clinical infections and demonstrate the high prevalence of carbapenem resistance genes on highly mobile transposons among these strains. Our findings highlight the pressing need to address hospital wastewater as a crucial factor in curtailing the spread of carbapenem-resistant A. baumannii.

Genomic analysis of penicillin-binding proteins and recombination events in an emerging amoxicillin- and meropenem-resistant PMEN3 (Spain9V-3, ST156) variant in Taiwan and comparison with global descendants of this lineage

From 2008 to 2020, the Taiwan National Notifiable Disease Surveillance System database demonstrated that the incidence of non-vaccine serotype 23A invasive pneumococcal disease (IPD) approximately doubled. In this study, 276 non-repetitive pneumococcal clinical isolates were collected from two medical centers in Taiwan between 2019 and 2021. Of these 267 pneumococci, 60 were serotype 23A. Among them, 50 (83%) of serotype 23A isolates belonged to the sequence type (ST) 166 variant of the Spain9V-3 clone. Pneumococcal 23A-ST166 isolates were collected to assess their evolutionary relationships using whole-genome sequencing. All 23A-ST166 isolates were resistant to amoxicillin and meropenem, and 96% harbored a novel combination of penicillin-binding proteins (PBPs) (1a:2b:2x):15:11:299, the newly identified PBP2x-299 in Taiwan. Transformation of the pbp1a, pbp2b, and pbp2x alleles into the β-lactam-susceptible R6 strain revealed that PBP2x-299 and PBP2b-11 increased the MIC of ceftriaxone and meropenem by 16-fold, respectively. Prediction analysis of recombination sites in PMEN3 descendants (23A-ST166 in Taiwan, 35B-ST156 in the United States, and 11A-ST838/ST6521 in Europe) showed that adaptive evolution involved repeated, selectively favored convergent recombination in the capsular polysaccharide synthesis region, PBPs, murM, and folP genome sites. In the late 13-valent pneumococcal conjugate vaccine era, PMEN3 continuously displayed an evolutionary capacity for global dissemination and persistence, increasing IPD incidence, leading to an offset in the decrease of pneumococcal conjugate vaccine serotype-related diseases, and contributing to high antibiotic resistance. A clonal shift with a highly β-lactam-resistant non-vaccine serotype 23A, from ST338 to ST166, increased in Taiwan. ST166 is a single-locus variant of the Spain9V-3 clone, which is also called the PMEN3 lineage. All 23A-ST166 isolates, in this study, were resistant to amoxicillin and meropenem, and 96% harbored a novel combination of penicillin-binding proteins (PBPs) (1a:2b:2x):15:11:299. PBP2x-299 and PBP2b-11 contributed to the increasing MIC of ceftriaxone and meropenem, respectively. Prediction analysis of recombination sites in PMEN3 descendants showed that adaptive evolution involved repeated, selectively favored convergent recombination in the capsular polysaccharide synthesis region, PBPs, murM, and folP genome sites. In the late 13-valent pneumococcal conjugate vaccine era, PMEN3 continuously displays the evolutionary capacity for dissemination, leading to an offset in the decrease of pneumococcal conjugate vaccine serotype-related diseases and contributing to high antibiotic resistance.

Tracing the Evolutionary Pathways of Serogroup O78 Avian Pathogenic Escherichia coli

Avian pathogenic E. coli (APEC) causes severe economic losses in the poultry industry, and O78 serogroup APEC strains are prevalent in chickens. In this study, we aimed to understand the evolutionary pathways and relationships between O78 APEC and other E. coli strains. To trace these evolutionary pathways, we classified 3101 E. coli strains into 306 subgenotypes according to the numbers and types of single nucleotide polymorphisms (RST0 to RST63-1) relative to the consensus sequence (RST0) of the RNA polymerase beta subunit gene and performed network analysis. The E. coli strains showed four apparently different evolutionary pathways (I-1, I-2, I-3, and II). The thirty-two Korean O78 APEC strains tested in this study were classified into RST4-4 (45.2%), RST3-1 (32.3%), RST21-1 (12.9%), RST4-5 (3.2%), RST5-1 (3.2%), and RST12-6 (3.2%), and all RSTs except RST21-1 (I-2) may have evolved through the same evolutionary pathway (I-1). A comparative genomic study revealed the highest relatedness between O78 strains of the same RST in terms of genome sequence coverage/identity and the spacer sequences of CRISPRs. The early-appearing RST3-1 and RST4-4 prevalence among O78 APEC strains may reflect the early settlement of O78 E. coli in chickens, after which these bacteria accumulated virulence and antibiotic resistance genes to become APEC strains. The zoonotic risk of the conventional O78 APEC strains is low at present, but the appearance of genetically distinct and multiple virulence gene-bearing RST21-1 O78 APEC strains may alert us to a need to evaluate their virulence in chickens as well as their zoonotic risk.