The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes

Demographic, clinical and molecular epidemiology of extended-spectrum beta-lactamase-producing Escherichia coli bloodstream infections in Central Australia

We describe the demographics, clinical and molecular epidemiology of extended-spectrum β-lactamase (ESBL) Escherichia coli bloodstream infections (BSI) in Central Australia. All ESBL-producing E. coli bloodstream isolates from January 2018 to December 2020 were retrospectively identified. Demographic and clinical information was extracted by chart review. Whole-genome sequencing was performed for multi-locus sequence typing, antibiotic-resistance genes, and phylogenetic relationships. We identified 41 non-duplicate episodes of ESBL E. coli BSI. Median age was 55 years (IQR 47-63), 78% were female, 93% were Aboriginal, and half came from a remote community. Infections were predominantly urinary (68%, 28/41). In the 12 months prior, 70% (26/37) of identified patients had been hospitalised and 81% (30/37) prescribed antibiotics. Meropenem and piperacillin-tazobactam susceptibility was maintained in 100% and 95% of isolates, respectively. Co-resistance to non-β-lactam antibiotics was 32% to gentamicin, 61% to trimethoprim/sulfamethoxazole, and 68% to ciprofloxacin. For sequenced isolates, 41% (16/35) were sequence type 131 (ST131). Mean acquired antibiotic-resistance genes for each isolate was 12.3 (SD 3.1). Four isolates carried an OXA-1 gene. Only non-ST131 isolates carried AmpC and acquired quinolone-resistance genes. There was some evidence of clustering of closely related strains, but no evidence of community or healthcare admission overlap. ESBL rates are rapidly rising in Central Australia, which is a conducive environment for antibiotic resistance development (e.g., overcrowding, socioeconomic disadvantages, high healthcare exposure and high antibiotic use). Future research is required to explore resistance-transmission dynamics in this unique setting.

Emergence of Staphylococcus argenteus in pediatrics: Molecular insights from a hospital in East China

Staphylococcus argenteus is a novel species within the Staphylococcus aureus complex and can cause serious bloodstream infections (BSIs) in humans, which have been mainly reported in adults, especially the elderly. In this study, we analyzed the molecular characterization of a strain of S. argenteus (22WJ8192) isolated from the peripheral vein blood sample of a seven-month-old female infant in Eastern China. The 22WJ8192 belonged to sequence type (ST)2250 and harbored six antibiotic-resistance genes and 53 virulence genes and was resistant to penicillin. Additionally, we conducted a comparative analysis of the molecular characteristics of S. argenteus sourced from various origins within the dataset, predominantly from the National Center for Biotechnology Information Collection (NCBI) genome database. Antibiotic-resistance genes blaR1, blaI_of_Z, blaZ, fosB-Saur, tet(L), aph(3")-IIIa, mecA, and dfrG were more prevalent among the strains of human origin. Virulence genes lukF-PV, sak, sdrE, scn, sdrC, and sdrD were more prevalent among strains of human origin. The presence of antibiotic-resistance genes blaR1, blaI_of_Z, blaZ, fosB-Saur, and aph(3")-IIIa in strain 22WJ8192 was also more common among strains of human origin in the dataset. Conversely, the antibiotic-resistance genes tet(L), mecA, and dfrG, typically found in strains of human origin, were not detected in 22WJ8192. Additionally, virulence genes lukF-PV, sak, sdrE, scn, sdrC, and sdrD present in 22WJ8192 exhibited a higher prevalence among strains of human origin in the dataset. In conclusion, this study emphasizes the potential of S. argenteus ST2250 to induce severe bloodstream infections in infants, shedding light on the molecular characteristics of this strain.

Genomic characterization and clinical evaluation of prosthetic joint infections caused by Cutibacterium acnes

Cutibacterium acnes is a major skin commensal that may act as an opportunistic pathogen. It is difficult to interpret findings of C. acnes in tissue cultures obtained during arthroplasty revision surgery, since they may represent true infection or contamination. This study investigated whether C. acnes obtained from prosthetic joint infections (PJIs) were related and shared common genomic traits that might correlate with clinical courses and patient outcomes. C. acnes isolates from revision surgery of patients with PJIs of the hip, shoulder, and knee were characterized using molecular methods to determine the sequence type (ST) and the presence of possible virulence determinants (Christie-Atkins-Munch-Peterson factors, dermatan sulfate-binding adhesion 1, hyaluronidase lyase, and linear plasmid). A standardized review of the patients' medical charts was performed. The study included 37 patients with C. acnes culture-positive tissue samples where multiple isolates of C. acnes belonged to the same ST. Most of the isolates belonged to phylotype IA1. Phylogenetic analysis of virulence determinants revealed no shared pattern among PJI isolates. Seven patients had a polymicrobial infection. Exchange revision was performed in 70% of the patients, and >50% of all patients received antibiotic treatment for ≥3 months. Failure was noted in seven patients. No specific ST or any identifiable unique feature among virulence determinants were found among C. acnes isolated from PJIs of hips and shoulders. The majority of patients had low inflammatory markers and were treated successfully, even polymicrobial infections. However, failure was more common among shoulder infections compared with hip infections.

IMPORTANCE

Prosthetic joint infection (PJI) is a rare complication after arthroplasty surgery. The infection seldom resolves without a combination of both surgical and antibiotic treatment and can cause significant suffering among affected patients. Cutibacterium acnes is a common skin bacterium that is most often found in shoulder PJIs but can also infect other prostheses. In this study, we conducted a review of patients with previously verified PJIs involving C. acnes in hip or shoulder prostheses, along with a genomic analysis of the bacteria causing the infections. The majority of patients had successful outcomes. We did not identify any specific phylogenetic lineage or specific molecular signature of virulence factors among these PJI-associated C. acnes isolates that seemed to be associated with increased potential to cause infection among this species. This indicates that C. acnes isolated from PJIs originates from the patients' own skin microbiome and is inoculated during the arthroplasty surgery.

First detection of VEB-1 extended-spectrum β-lactamase-producing Escherichia coli clinical isolate in Japan

The extended-spectrum β-lactamase (ESBL) gene, blaVEB-1, was identified for the first time in an Escherichia coli clinical isolate, JARB-RN-0061, from blood cultures in a Japanese general hospital in 2021. The isolate exhibited high resistance to broad-spectrum cephalosporins, including ceftazidime (MIC >128 mg/L) and cefepime (MIC = 16 mg/L). blaVEB-1 was identified during whole-genome sequencing and characterization of the isolate. JARB-RN-0061 belonged to the B2-O2:K1:H7-ST95-fimH41 lineage and was classified as presumptive extraintestinal pathogenic E. coli (ExPEC) and uropathogenic E. coli (UPEC). Moreover, the strain harbored multiple virulence genes on the chromosome. The Col156/IncFIB(AP001918)/IncFII(29)-type plasmid (114,216 bp), with clbB and tcpC genes involved in bacteremia, was unique to the fimH41 subclone. The blaVEB-1 gene was located on a non-typeable and non-conjugative plasmid, pJARB-RN-0061_VEB-1 (17,093 bp). It was embedded in the class 1 integron In1883-like, with multidrug resistance gene cassettes for aacA4, aadB, cmlA5, qnrVC4, and dfrA14. Notably, comparative analysis of the complete sequence of plasmid pJARB-RN-0061_VEB-1 revealed that it was highly homologous to the blaVEB-1-harboring plasmid, pMS2H7VEB-1 (100% coverage and 99.99% identity), except for the Tn3 family transposon (4,931 bp) and the plasmid pRHBSTW-00138_5 (97% coverage and 100% identity) harbored by Klebsiella quasipneumoniae subsp. similipneumoniae strains from hospital sewage in Japan and wastewater influent in the United Kingdom, respectively. The emergence of a human pathogenic E. coli clinical isolate with the blaVEB-1-carrying plasmid in the B2-ST95 worldwide pandemic lineage, characterized by the virulence potential of ExPEC or UPEC but a low prevalence of antimicrobial resistance, would raise public health concerns.

IMPORTANCE

ESBLs are plasmid-mediated enzymes that confer resistance to clinically significant antimicrobial agents, such as broad-spectrum cephalosporins. Recently, the rapid spread of CTX-M-type ESBL-producing E. coli has become a global issue, including in Japan, where ESBL production in human pathogenic E. coli, such as the ExPEC and UPEC lineages, which typically harbor several virulence genes, is a severe public health concern. To date, VEB (Vietnamese extended-spectrum β-lactamase) producers have been found only in hospital wastewater and rivers in Japan. Thus, we describe the first detection of a very rare human-derived blaVEB-1 gene in the E. coli B2-ST95 pandemic clonal lineage that is highly associated with ExPEC and UPEC in a Japanese clinical setting. Furthermore, we characterized the genomic features of plasmids harboring the class 1 integron-borne blaVEB-1. Our findings highlight the significance of closely monitoring ESBL-producing E. coli isolates to prevent the potential dissemination of this resistance determinant in Japan.

Assessment and genomic analysis of Salmonella and Campylobacter from different stages of an integrated no-antibiotics-ever (NAE) broiler complex: a longitudinal study

The objective of this study was to determine prevalence and perform genomic analysis of Salmonella spp. and Campylobacter spp. isolated from different stages of an integrated NAE broiler complex. Environmental samples were screened with 3M-Molecular Detection System (MDS) and MDS positive samples were further processed for confirmation of results and identification. Core genome-based phylogenies were built for both bacteria isolated from this study along with selected NCBI genomes. The odds ratios and 95% confidence limits were compared among stages and sample types (α < 0.05) using multivariable model. Based on MDS results, 4% and 18% of total samples were positive for Salmonella spp. and Campylobacter spp. respectively. The odds of Salmonella detection in hatchery samples were 2.58 times as likely as compared to its detection in production farms' samples (P = 0.151) while the odds of Campylobacter detection in production farms' samples were 32.19 times as likely as its detection in hatchery (P = 0.0015). Similarly, the odds of Campylobacter detection in boot swabs, soil, water, and miscellaneous samples were statistically significant (P < 0.05) as compared with fly paper as reference group. The serovars identified for Salmonella were Typhimurium, Barranquilla, Liverpool, Kentucky, Enteritidis, Luciana, and Rough_O:r:1,5. For Campylobacter, the species identified were Campylobacter jejuni and Campylobacter coli. Phylogeny results show close genetic relatedness among bacterial strains isolated from different locations within the same stage and between different stages. The results show possibility of multiple entry points of such bacteria entering broiler complex and can potentially contaminate the final raw product in the processing plant. It suggests the need for a comprehensive control strategy with strict biosecurity measures and best management practices to minimize or eliminate such pathogens from the poultry food chain.

Prevalence and genetic characterization of linezolid resistance gene reservoirs in hospital sewage from Zhejiang Province, China

Hospital sewage represented important hotspots for the aggregation and dissemination of clinically relevant pathogens and antimicrobial resistance genes. To investigate the prevalence and molecular epidemiology of linezolid resistance genes in hospital sewage, both influent and effluent samples from 11 hospitals in Zhejiang Province, China, were collected and analyzed for linezolid resistance gene carriers. Thirty colonies of putative isolates that grew on the selective media with 10 mg/L florfenicol were randomly picked per sample. A total of 420 Gram-positive isolates, including 330 from 11 influent samples and 90 from three effluent samples, were obtained. Each isolate carried at least one of the linezolid resistance genes, including optrA, poxtA, cfr, and cfr(D), and the optrA gene was highly dominant (388/420). Enterococci displayed predominance among the linezolid resistance gene carriers in the hospital sewage, exhibiting a resistance rate to linezolid of 77.8 %. The wild-type OptrA and OptrA variants KLDP, RDK, and KLDK, all associated with high linezolid MICs, were most frequently detected. Phylogenetic analysis revealed the multispecies and polyclonal distribution of linezolid-resistant bacteria in hospital sewage, while Enterococcus faecalis sequence types (STs) 16 and 179 demonstrated the widest dissemination across different hospitals. Despite generally high genetic diversity, phylogenetic analysis showed that 87 isolates, assigned to ten STs from both sewage and other sources, were genetically related. Moreover, the genetic environment of linezolid resistance genes in isolates from sewage was similar to that from animals, humans, or the environment, with "Tn554-fexA-optrA" as the most common structure. These findings revealed the potential risk of the transmission of linezolid resistance genes through hospital sewage to other environments.

Diversity of Campylobacter spp. circulating in a rhesus macaque (Macaca mulatta) breeding colony using culture and molecular methods

Campylobacter jejuni and Campylobacter coli represent the leading causes of bacterial gastroenteritis in humans, and infections can produce post-infectious irritable bowel syndrome (PI-IBS). Rhesus macaques (Macaca mulatta) (RM) are similarly susceptible to acute campylobacteriosis and represent a potential model of PI-IBS. We characterized the Campylobacter species circulating in an RM breeding colony using culture, qPCR, and whole genome sequencing (WGS). We also compared the C. jejuni and C. coli prevalence in RM as detected with qPCR versus culture and identified risk factors for bacteria presence and intestinal disease. Culture of 275 samples yielded C. coli (103) and C. jejuni (8), of which 21.6% were resistant to quinolones and 3.6% were resistant to macrolides. Multidrug-resistant isolates were obtained exclusively from animals exhibiting diarrhea or with histologically confirmed chronic enterocolitis. WGS revealed a non-clonal population of Campylobacter spp. Genotypic predictions of resistance were excellent except for aminoglycosides. All sequenced isolates contained genes for all subunits of cytolethal distending toxin. qPCR detected a prevalence of 45.9% for C. coli and 29.6% for C. jejuni. The quantity of either bacteria was significantly higher (P < 0.05) in animals with intestinal disease compared to healthy animals, though only young age was significantly associated with the presence of Campylobacter sp. or intestinal disease. Significantly more C. jejuni positive animals were detected with qPCR than with culture. These results provide a comprehensive characterization of Campylobacter spp. circulating in a breeding colony of RM in the United States and suggest that qPCR is superior for the detection of C. jejuni in RM.

IMPORTANCE

Gastrointestinal disease is one of the most common reasons for hospitalization in non-human primate colonies and accounts for over one-third of non-research related euthanasia. In rhesus macaques, this manifests as both acute diarrhea and chronic enterocolitis (CE), a syndrome of chronic diarrhea resulting in poor weight gain or weight loss which is minimally responsive to treatment. Campylobacter spp. are major causes of acute enterocolitis in rhesus macaques and may predispose individuals to the development of CE, similar to post-infectious irritable bowel syndrome in humans. Despite these concerns, there are few studies characterizing Campylobacter in rhesus macaque colonies, in particular utilizing whole genome sequencing and assessing findings with respect to the health status of the host. Our findings provide insight into Campylobacter strains circulating in rhesus macaque colonies, which can improve clinical monitoring, assist in treatment decisions, and provide new avenues of investigation into campylobacteriosis as a catalyst for CE.

Global prevalence of Mycobacterium massiliense strains with recombinant rpoB genes (Rec-Mas) horizontally transferred from Mycobacterium abscessus: two major types, dominant circulating clone 7 and MLST ST46 sequence type

Mycobacterium abscessus is a group of emerging antimicrobial-resistant nontuberculous mycobacteria that causes severe lung disease in infected patients globally. Recently, molecular epidemiology studies have indicated that horizontal gene transfer (HGT) events in the rpoB gene are prevalent between M. abscessus subspecies. To determine the global prevalence of M. abscessus strains subjected to rpoB HGT, we performed phylogenetic inference using a 711-bp rpoB sequence extracted from 1,786 M. abscessus isolates for which the whole-genome sequence was publicly available. Our data showed that a total of 74 isolates (4.1%) from 1,786 strains are subject to rpoB HGT, which is more prevalent than strains with hsp65 HGT (19 isolates from 1,786, 1.1%). Most of these (69 isolates) belong to two major groups of Mycobacterium massiliense, of which the rpoB gene is horizontally transferred from M. abscessus (Rec-mas), dominant circulating clone 7 (DCC7) (44 isolates) and ST46 type by multilocus sequence typing (25 isolates). The Rec-mas strains of the two groups have distinct geographical patient distributions, of which the former is mainly distributed in the United States, while the latter is prevalent in Asia. Our further genome-based analysis indicated that the ST46 type is a novel DCC candidate of M. massiliense that is responsible for dissemination between noncystic fibrosis patients in Asia. In conclusion, our global phylogenetic analysis revealed two major Rec-mas clones with distinct geographical distributions, namely, DCC7 and ST46. This study provides insights into the genetic clustering and person-to-person transmission of globally dominant and area-specific strains harboring the HGT rpoB gene.

IMPORTANCE

Horizontal gene transfer (HGT) events play a pivotal role in the evolution of Mycobacterium abscessus into dominant circulating clones (DCCs), which is capable of causing patient-to-patient transmission. In particular, HGT of the rpoB gene between strains of different subspecies of M. abscessus could also compromise differentiation between strains of M. abscessus. Here, for the first time, using 1,786 M. abscessus genome sequences, we evaluated the global prevalence of M. abscessus strains subjected to rpoB HGT. We found a greater prevalence of M. abscessus subjected to rpoB HGT than to those subjected to hsp65 HGT, which is mainly due to two Rec-mas clones, dominant circulating clone 7 and ST46, which are responsible for dissemination between non-CF patients in Asia. Our data highlight the importance of rpoB HGT in the evolution of M. abscessus, particularly Mycobacterium massiliense, into virulent DCC clones.

Emerging blaNDM-positive Salmonella enterica in Chinese pediatric infections

Non-typhoidal Salmonella (NTS) is a common zoonotic foodborne pathogen, whose rising antimicrobial resistance has been an urgent threat to global public health. Here, we reported two carbapenem-resistant Salmonella enterica (CRSE) strains (NBFE-049 and NBFE-164) carrying the blaNDM genes, which were recovered from two Chinese children, belonged to Salmonella enterica serovar Typhimurium (S. Typhimurium) monophasic variant (S. 1,4,[5],12:i:-) ST34 (a sequence type) and S. Typhimurium ST19, respectively. Genes blaNDM-5 and blaNDM-13 were detected in NBFE-049 and NBFE-164, respectively. The blaNDM-5 in NBFE-049 was located in an IncHI2-type plasmid, named pNBFE-049. In NBFE-164, the blaNDM-13 was located in an IncI1-type plasmid, named pNBFE-164. The plasmid pNBFE-164 successfully transferred its resistance phenotype into the recipient strain Escherichia coli J53 with a high efficiency of 1.1 × 10-2, while no transconjugants were obtained in pNBFE-049 conjugation assays. We further elucidated the genetic relationships of globally occurring New Delhi Metallo-β-lactamase (NDM)-positive strains and locally distributed clinical strains within the same serovar. The closest relative of NBFE-049 was clinical Salmonella strain 1722, which was recovered in 2020 and differed by only three Single Nucleotide Polymorphisms (SNPs). No NDM-positive ST19 could be found in the National Center for Biotechnology Information (NCBI) database, and NBFE-164 showed a close genetic relationship with the other ST19 in this area. To sum up, we suggested the potential contributions of clonal spread and plasmid-mediated blaNDM transfer in CRSE dissemination. This study reported the complete genome of two blaNDM-carrying S. Typhimurium isolates, shedding new insights into the antimicrobial resistance mechanisms and dissemination patterns of the emerging CRSE.IMPORTANCENTS is one of the most common zoonotic pathogens that causes foodborne illnesses, while S. Typhimurium is one of the most common serovars. With the rising prevalence of multi-resistant Salmonella worldwide, carbapenems have emerged as the last-line antibiotics for treating severe bacterial infections. In this study, we reported the genomic characteristics of two carbapenem-resistant S. Typhimurium strains, which were recovered from two pediatric patients, carrying blaNDM-5 and blaNDM-13, providing new insights into the antimicrobial resistance deteriminants and transmission risk of blaNDM-positive NTS in China. We suggested the potential contributions of clonal spread and plasmid-mediated blaNDM transfer in CRSE dissemination. Future enhanced surveillance policy should mitigate CRSE spreading, and more importantly, clinical antimicrobial therapeutic regimens should be adjusted accordingly.

Development of a High-Resolution Melting Method for the Detection of Clarithromycin-Resistant Helicobacter pylori in the Gastric Microbiome

Background: The issue of Helicobacter pylori (H. pylori) resistance to clarithromycin (CLR) has consistently posed challenges for clinical treatment. Hence, a rapid susceptibility testing (AST) method urgently needs to be developed. Methods: In the present study, 35 isolates of H. pylori were isolated from 203 gastritis patients of the Guangzhou cohort, and the antimicrobial resistance phenotypes were associated with their genomes to analyze the relevant mutations. Based on these mutations, a rapid detection system utilizing high-resolution melting (HRM) curve analysis was designed and verified by the Shenzhen cohort, which consisted of 38 H. pylori strains. Results: Genomic analysis identified the mutation of the 2143 allele from A to G (A2143G) of 23S rRNA as the most relevant mutation with CLR resistance (p < 0.01). In the HRM system, the wild-type H. pylori showed a melting temperature (Tm) of 79.28 ± 0.01 °C, while the mutant type exhibited a Tm of 79.96 ± 0.01 °C. These differences enabled a rapid distinction between two types of H. pylori (p < 0.01). Verification examinations showed that this system could detect target DNA as low as 0.005 ng/μL in samples without being affected by other gastric microorganisms. The method also showed a good performance in the Shenzhen validation cohort, with 81.58% accuracy, and 100% specificity. Conclusions: We have developed an HRM system that can accurately and quickly detect CLR resistance in H. pylori. This method can be directly used for the detection of gastric microbiota samples and provides a new benchmark for the simple detection of H. pylori resistance.

Serratia marcescens outbreak in a neonatal intensive care unit associated with contaminated handwashing sinks

PURPOSE

This article describes the origin of a S. marcescens outbreak in a neonatal intensive care unit (NICU).

MATERIALS AND METHODS

A retrospective case-control study including 12 S. marcescens-positive and 22 S. marcescens-negative neonates in the NICU was performed to identify the source of the outbreak. S. marcescens isolates were collected during the outbreak and analyzed using whole-genome sequencing (WGS). IQ-Tree software, BEAST2 software package and SCOTTI software were used to construct a phylogenetic tree and a propagation path map.

RESULTS

The index case occurred on February 21st and outbreak ended on March 9th, 2021, affecting a total of 12 neonates (2 with S. marcescens infection and 10 with S. marcescens colonization). Multivariate logistic regression identified that the distance of <0.8 m between the bed unit and the sink (odds ratio [OR], 20.50; 95 % confidence interval [CI], 1.09-384.86), a large number of rotating nurses within a week (OR 2.58, 95 % CI, 1.09-6.11) and use of humidification water in the incubator (OR 189.70, 95 % CI, 2.76-13027.31) were significant increased risk factors for S. marcescens infection or colonization in the outbreak. WGS sifted out a predominant clone between contaminated handwashing sinks and patients, suggesting that cross-transmission was involved in the dissemination of S. marcescens.

CONCLUSION

Contaminated handwashing sinks can be a communication intermediary of S. marcescens infection or colonization of neonates in the NICU. A distance of <0.8 m between the bed unit and the sink, and a large number of rotating nurses might play important roles in this outbreak. Attention should be paid to sinks contamination and contact transmission to prevent outbreaks.

Expanded genome and proteome reallocation in a novel, robust Bacillus coagulans strain capable of utilizing pentose and hexose sugars

Bacillus coagulans, a Gram-positive thermophilic bacterium, is recognized for its probiotic properties and recent development as a microbial cell factory. Despite its importance for biotechnological applications, the current understanding of B. coagulans' robustness is limited, especially for undomesticated strains. To fill this knowledge gap, we characterized the metabolic capability and performed functional genomics and systems analysis of a novel, robust strain, B. coagulans B-768. Genome sequencing revealed that B-768 has the largest B. coagulans genome known to date (3.94 Mbp), about 0.63 Mbp larger than the average genome of sequenced B. coagulans strains, with expanded carbohydrate metabolism and mobilome. Functional genomics identified a well-equipped genetic portfolio for utilizing a wide range of C5 (xylose, arabinose), C6 (glucose, mannose, galactose), and C12 (cellobiose) sugars present in biomass hydrolysates, which was validated experimentally. For growth on individual xylose and glucose, the dominant sugars in biomass hydrolysates, B-768 exhibited distinct phenotypes and proteome profiles. Faster growth and glucose uptake rates resulted in lactate overflow metabolism, which makes B. coagulans a lactate overproducer; however, slower growth and xylose uptake diminished overflow metabolism due to the high energy demand for sugar assimilation. Carbohydrate Transport and Metabolism (COG-G), Translation (COG-J), and Energy Conversion and Production (COG-C) made up 60%-65% of the measured proteomes but were allocated differently when growing on xylose and glucose. The trade-off in proteome reallocation, with high investment in COG-C over COG-G, explains the xylose growth phenotype with significant upregulation of xylose metabolism, pyruvate metabolism, and tricarboxylic acid (TCA) cycle. Strain B-768 tolerates and effectively utilizes inhibitory biomass hydrolysates containing mixed sugars and exhibits hierarchical sugar utilization with glucose as the preferential substrate.IMPORTANCEThe robustness of B. coagulans makes it a valuable microorganism for biotechnology applications; yet, this phenotype is not well understood at the cellular level. Through phenotypic characterization and systems analysis, this study elucidates the functional genomics and robustness of a novel, undomesticated strain, B. coagulans B-768, capable of utilizing inhibitory switchgrass biomass hydrolysates. The genome of B-768, enriched with carbohydrate metabolism genes, demonstrates high regulatory capacity. The coordination of proteome reallocation in Carbohydrate Transport and Metabolism (COG-G), Translation (COG-J), and Energy Conversion and Production (COG-C) is critical for effective cell growth, sugar utilization, and lactate production via overflow metabolism. Overall, B-768 is a novel, robust, and promising B. coagulans strain that can be harnessed as a microbial biomanufacturing platform to produce chemicals and fuels from biomass hydrolysates.

Exploring genetic determinants of antimicrobial resistance in Brucella melitensis strains of human and animal origin from India

Introduction

Antimicrobial resistance (AMR) in Brucella melitensis, the causative agent of brucellosis, is of growing concern, particularly in low and middle-income countries. This study aimed to explore the genetic basis of AMR in B. melitensis strains from India.

Methods

Twenty-four isolates from humans and animals were subjected to antimicrobial susceptibility testing and whole-genome sequencing.

Results

Resistance to doxycycline (20.80%), ciprofloxacin (16.67%), cotrimoxazole (4.17%), and rifampicin (16.67%) was observed. Genome analysis revealed efflux-related genes like mprF, bepG, bepF, bepC, bepE, and bepD across all isolates, however, classical AMR genes were not detected. Mutations in key AMR-associated genes such as rpoB, gyrA, and folP were identified, intriguingly present in both resistant and susceptible isolates, suggesting a complex genotype-phenotype relationship in AMR among Brucella spp. Additionally, mutations in efflux genes were noted in resistant and some susceptible isolates, indicating their potential role in resistance mechanisms. However, mutations in AMR-associated genes did not consistently align with phenotypic resistance, suggesting a multifactorial basis for resistance.

Discussion

The study underscores the complexity of AMR in B. melitensis and advocates for a holistic multi-omics approach to fully understand resistance mechanisms. These findings offer valuable insights into genetic markers associated with AMR, guiding future research and treatment strategies.

Characterization of acquired β-lactamases in Pseudomonas aeruginosa and quantification of their contributions to resistance

Pseudomonas aeruginosa is a highly problematic opportunistic pathogen that causes a range of different infections. Infections are commonly treated with β-lactam antibiotics, including cephalosporins, monobactams, penicillins, and carbapenems, with carbapenems regarded as antibiotics of last resort. Isolates of P. aeruginosa can contain horizontally acquired bla genes encoding β-lactamase enzymes, but the extent to which these contribute to β-lactam resistance in this species has not been systematically quantified. The overall aim of this research was to address this knowledge gap by quantifying the frequency of β-lactamase-encoding genes in P. aeruginosa and by determining the effects of β-lactamases on susceptibility of P. aeruginosa to β-lactams. Genome analysis showed that β-lactamase-encoding genes are present in 3% of P. aeruginosa but are enriched in carbapenem-resistant isolates (35%). To determine the substrate antibiotics, 10 β-lactamases were expressed from an integrative plasmid in the chromosome of P. aeruginosa reference strain PAO1. The β-lactamases reduced susceptibility to a variety of clinically used antibiotics, including carbapenems (meropenem, imipenem), penicillins (ticarcillin, piperacillin), cephalosporins (ceftazidime, cefepime), and a monobactam (aztreonam). Different enzymes acted on different β-lactams. β-lactamases encoded by the genomes of P. aeruginosa clinical isolates had similar effects to the enzymes expressed in strain PAO1. Genome engineering was used to delete β-lactamase-encoding genes from three carbapenem-resistant clinical isolates and increased susceptibility to substrate β-lactams. Our findings demonstrate that acquired β-lactamases play an important role in β-lactam resistance in P. aeruginosa, identifying substrate antibiotics for a range of enzymes and quantifying their contributions to resistance.IMPORTANCEPseudomonas aeruginosa is an extremely problematic pathogen, with isolates that are resistant to the carbapenem class of β-lactam antibiotics being in critical need of new therapies. Genes encoding β-lactamase enzymes that degrade β-lactam antibiotics can be present in P. aeruginosa, including carbapenem-resistant isolates. Here, we show that β-lactamase genes are over-represented in carbapenem-resistant isolates, indicating their key role in resistance. We also show that different β-lactamases alter susceptibility of P. aeruginosa to different β-lactam antibiotics and quantify the effects of selected enzymes on β-lactam susceptibility. This research significantly advances the understanding of the contributions of acquired β-lactamases to antibiotic resistance, including carbapenem resistance, in P. aeruginosa and by implication in other species. It has potential to expedite development of methods that use whole genome sequencing of infecting bacteria to inform antibiotic treatment, allowing more effective use of antibiotics, and facilitate the development of new antibiotics.

WGS Analysis of Staphylococcus warneri Outbreak in a Neonatal Intensive Care Unit

Purpose

Staphylococcus warneri is an opportunistic pathogen responsible for hospital-acquired infections (HAIs). The aim of this study was to describe an outbreak caused by S. warneri infection in a neonatal intensive care unit (NICU) and provide investigation, prevention and control strategies for this outbreak.

Methods

We conducted an epidemiological investigation of the NICU S. warneri outbreak, involving seven neonates, staff, and environmental screening, to identify the source of infection. WGS analyses were performed on S. warneri isolates, including species identification, core genome single-nucleotide polymorphism (cgSNP) analysis, pan-genome analysis, and genetic characterization assessment of the prevalence of specific antibiotic resistance and virulence genes.

Results

Eight S. warneri strains were isolated from this outbreak, with seven from neonates and one from environment. Six clinical cases within three days in 2021 were linked to one strain isolated from environmental samples; isolates varied by 0-69 SNPs and were confirmed to be from an outbreak through WGS. Multiple infection prevention measures were implemented, including comprehensive environmental disinfection and stringent protocols, and all affected neonates were transferred to the isolation wards. Following these interventions, no further cases of S. warneri infections were observed. Furthermore, pan-genome analysis results suggested that in human S. warneri may exhibit host specificity.

Conclusion

The investigation has revealed that the outbreak was linked to the milk preparation workbench by the WGS. It is recommended that there be a stronger focus on environmental disinfection management in order to raise awareness, improve identification, and prevention of healthcare-associated infections that are associated with the hospital environment.

Genomic insights into epidemic plasmids carrying bla CTX-M and mcr-1 genes in Escherichia coli from Lebanese broiler production

Background

In a previous nationwide survey in the Lebanese broiler production, multidrug-resistant CTX-M-producing E. coli were found to carry the mobile colistin resistance gene mcr-1.

Objectives

To investigate the mobile genetic supports responsible for the spread of these resistance genes among E. coli in healthy broilers in Lebanon.

Methods

Thirty-three bla CTX-M and mcr-1 positive E. coli of various sequence types from 17 broilers farms were subjected to conjugation assays. Long-read sequencing (Oxford Nanopore Technologies) and hybrid assembly were performed to determine complete plasmid sequences and their phylogenetic diversity.

Results

Twenty-nine conjugative IncFII plasmids harboured the extended-spectrum β-lactamase genes bla CTX-M-3 (n = 25) or bla CTX-M-55 (n = 4). Highly related IncF2:A-:B-/bla CTX-M-3 plasmids differing only through IS-mediated genetic rearrangements in antibiotic resistance gene clusters were found in genetically diverse E. coli strains isolated from distant farms. The mobile colistin resistance genes mcr-1.1 and mcr-1.26 were carried by IncX4 and IncI2 plasmids. Worryingly, in one isolate, the ISEcp1-bla CTX-M-55 transposable unit was found integrated in a mcr-1.26-carrying IncX4 plasmid. Beside expanded cephalosporins and colistin resistances, all E. coli isolates were multidrug-resistant with different additional resistances against aminoglycosides, (fluoro)quinolones, fosfomycin, phenicols, sulphonamides, tetracycline and trimethoprim.

Conclusions

Closely related blaCTX-M-3/55-borne IncF2:A-:B- plasmids harbouring variable MDR regions and mcr-1 carrying IncX4 plasmids are widely disseminated in the E. coli population of healthy broilers in Lebanon. Further surveillance programmes of antimicrobial resistance and interventions to reduce the abusive use of medically important antibiotics are necessary to limit the spread of resistances in food-producing animals in Lebanon.

Plasmid-mediated azithromycin resistance in non-typhoidal Salmonella recovered from human infections

OBJECTIVES

Mechanisms of non-typhoidal Salmonella (NTS) resistance to azithromycin have rarely been reported. Here we investigate the epidemiology and genetic features of 10 azithromycin-resistant NTS isolates.

METHODS

A total of 457 NTS isolates were collected from a tertiary hospital in Guangzhou. We performed antimicrobial susceptibility tests, conjugation experiments, efflux pump expression tests, whole-genome sequencing and bioinformatics analysis to conduct the study.

RESULTS

The results showed that 10 NTS isolates (2.8%) were resistant to azithromycin with minimum inhibitory concentration values ranging from 128 to 512 mg/L and exhibited multidrug resistance. The phylogenetic tree revealed that 5 S. London isolates (AR1-AR5) recognized at different times and departments were closely related [3-74 single-nucleotide polymorphisms (SNPs)] and 2 S. Typhimurium isolates (AR7 and AR8) were clones (<3 SNPs) at 3-month intervals. The azithromycin resistance was conferred by mph(A) gene found on different plasmids, including IncFIB, IncHI2, InFII, IncC and IncI plasmids. Among them, IncFIB, InFII and IncHI2 plasmids carried different IS26-class 1 integron (intI1) arrangement patterns that mediated multidrug resistance transmission. Conjugative IncC plasmid encoded resistance to ciprofloxacin, ceftriaxone and azithromycin. Furthermore, phylogenetic analysis demonstrated that mph(A)-positive plasmids closely related to 10 plasmids in this study were mainly discovered from NTS, Escherichia coli, Klebsiella pneumonia and Enterobacter hormaechei. The genetic environment of mph(A) in 10 NTS isolates was IS26-mph(A)-mrx(A)-mphR(A)-IS6100/IS26 that co-arranged with intI1 harbour multidrug-resistant (MDR) gene cassettes on diverse plasmids.

CONCLUSIONS

These findings highlighted that the dissemination of these plasmids carrying mph(A) and various intI1 MDR gene cassettes would seriously restrict the availability of essential antimicrobial agents for treating NTS infections.

Genome Analysis of Potential Probiotic Levilactobacillus brevis AcCh91 Isolated from Indian Home-Made Fermented Milk Product (Chhurpi)

Consumption of naturally fermented milk (NFM) products is the dietary culture in India. The mountainous people of Arunachal Pradesh in India prepare the assorted artisanal home-made NFM products from cow and yak milk. Previously, we isolated and identified 76 strains of lactic acid bacteria (LAB) from NFM products of Arunachal Pradesh, viz. mar, chhurpi and churkam. We hypothesized that some of these LAB strains may possess probiotic potentials; hence, we investigated the probiotic potentials of these strains. On the basis of in vitro and genetic screening for probiotic attributes including haemolytic ability, 20 LAB strains were selected out of 76 strains, for further analysis. Using in silico analysis, viz. multivariate heatmap and PCA (principal component analysis) biplot, Levilactobacillus brevis AcCh91 was selected as the most promising probiotic strain, which was further characterized by the whole-genome analysis. Lev. brevis AcCh91 showed the highest survival rate of 93.38% in low pH and 86.68 ± 2.69% in low bile and the highest hydrophobicity average of 86.34 ± 5.53%. This strain also showed auto-aggregation and co-aggregation with antimicrobial properties against the pathogens, showed ability to produce beta-galactosidase and cholesterol reduction property and, most importantly, produced GABA, an important psychobiotic element. Genomic analysis of Lev. brevis AcCh91 showed the presence of genes corresponding to GABA, vitamins, amino acids, cholesterol reduction, immunomodulation, bioactive peptides and antioxidant activity. The absence of antimicrobial-resistant genes and virulence factors was observed. Hence, genome analysis supports the probiotic potentials of Lev. brevis AcCh91, which may be further investigated to understand its health-promoting properties.

Contribution of the Mobilome to the Configuration of the Resistome of Corynebacterium striatum

Corynebacterium striatum, present in the microbiota of human skin and nasal mucosa, has recently emerged as a causative agent of hospital-acquired infections, notable for its resistance to multiple antimicrobials. Its mobilome comprises several mobile genetic elements, such as plasmids, transposons, insertion sequences and integrons, which contribute to the acquisition of antimicrobial resistance genes. This study analyzes the contribution of the C. striatum mobilome in the transfer and dissemination of resistance genes. In addition, integrative and conjugative elements (ICEs), essential in the dissemination of resistance genes between bacterial populations, whose role in C. striatum has not yet been studied, are examined. This study examined 365 C. striatum genomes obtained from the NCBI Pathogen Detection database. Phylogenetic and pangenome analyses were performed, the resistance profile of the bacterium was recognized, and mobile elements, including putative ICE, were detected. Bioinformatic analyses identified 20 antimicrobial resistance genes in this species, with the Ermx gene being the most predominant. Resistance genes were mainly associated with plasmid sequence regions and class 1 integrons. Although an ICE was detected, no resistance genes linked to this element were found. This study provided valuable information on the geographic spread and prevalence of outbreaks observed through phylogenetic and pangenome analyses, along with identifying antimicrobial resistance genes and mobile genetic elements that carry many of the resistance genes and may be the subject of future research and therapeutic approaches.

First Case of Candida Auris Sepsis in Southern Italy: Antifungal Susceptibility and Genomic Characterisation of a Difficult-to-Treat Emerging Yeast

Candida auris is an emerging yeast considered a serious threat to global health. We report the first case of C. auris candidemia in Southern Italy, characterized using whole genome sequencing (WGS), and compared with a second strain isolated from a patient who presented as C. auris-colonized following screening. The C. auris strain was isolated from clinical samples, identified via MALDI-TOF, and subjected to WGS. Antifungal susceptibility testing was performed using commercial broth microdilution plates, and resistance protein sequences were evaluated with TBLASTN-2.15.0. Following the initial C. auris isolation from patient A, active surveillance and environmental investigations were implemented for all ICU patients. Of the 26 ICU surfaces sampled, 46.1% tested positive for C. auris via real-time PCR. Screening identified a second patient (patient B) as C. auris-colonized. The phylogenetic characterization of strains from patients A and B, based on the D1/D2 region of the 28s rDNA and the internal transcribed spacer (ITS) region, showed high similarity with strains from Lebanon. SNP analysis revealed high clonality, assigning both strains to clade I, indicating a significant similarity with Lebanese strains. This case confirms the alarming spread of C. auris infections and highlights the need for stringent infection control measures to manage outbreaks.

Association analysis of antibiotic and disinfectant resistome in human and foodborne E. coli in Beijing, China

The widespread use of chemical disinfectants and antibiotics poses a major threat to food safety and human health. However, the mechanisms of co-transmission of antimicrobial resistance genes (ARGs) and biocides and metal resistance genes (BMRGs) of foodborne pathogens in the food chain is still unclear. This study isolated 343 E. coli strains from animal-derived foods in Beijing and incorporated online data of human-derived E. coli strains from NCBI. Our results demonstrated a relatively uniform distribution of strains from various regions in Beijing, indicating a lack of region-specific clustering. Additionally, predominant sequence types varied between food- and human-derived strains, suggesting a preference for different hosts and environments. Phenotypic association analysis showed that the chlorine disinfectants peroxides had a significant positive correlation with tetracyclines. Many more ARGs and BMRGs were enriched in human-associated E. coli compared with those in chicken- and pork-origin. The quaternary ammonium compounds (QACs) resistance gene qacEΔ1 had a strong correlation with aminoglycoside resistance gene aadA5, folate pathway antagonist resistance gene dfrA17, sul1 and macrolide resistance gene mph(A). The correlation results indicated a significant association between the copper resistance gene cluster pco and the silver resistance gene cluster sil. Coexistence of many resistance genes was observed within the qacEΔ1 gene structure, with qacEΔ1-sul1 being the most common combination. Our findings demonstrated that the epidemiological spread of resistance is affected by a combination of heavy metals, disinfectants and antibiotic use, suggesting that the prevention and control strategies of antimicrobial resistance need to be multifaceted and comprehensive.

Sequence Types and Antimicrobial Resistance Profiles of Salmonella Typhimurium in the Food Chain in Singapore

Salmonella remains a significant foodborne pathogen globally with S. Typhimurium presenting as a frequently occurring serovar. This study aimed to characterize 67 S. Typhimurium isolates from humans, food, farms, and slaughterhouses collected in Singapore from 2016 to 2017. Using whole-genome sequencing analysis, the isolates were found to belong to either ST19 (n = 33) or ST36 (n = 34). ST36 predominated in human intestinal and chicken isolates, while human extra-intestinal and non-chicken food isolates belonged to ST19. Plasmids were predicted in 88.1% (n = 59) of the isolates with the most common incompatibility group profiles being IncFIB(S), IncFII(S) and IncQ1. IncFIB(S) (adjusted p-value < 0.05) and IncFII(S) (adjusted p-value < 0.05) were significantly more prevalent in ST19 isolates, while Col156 (adjusted p-value < 0.05) was more significantly found in ST36 isolates. ST36 isolates exhibited higher resistance to multiple antibiotic classes such as penicillins, phenicols, folate pathway inhibitors, aminoglycosides, β-lactam/β-lactamase inhibitor combinations, tetracyclines, and fluoroquinolones. Phylogenetics analysis suggested potential shared routes of transmission among human, chicken, farm and slaughterhouse environments. Taken together, this study offers a cross-sectional epidemiological insight into the genomic epidemiology and antimicrobial landscape of S. Typhimurium isolates in Singapore, informing strategies for future public health and food safety surveillance.

Phenotypic and genomic analysis of the hypervirulent methicillin-resistant Staphylococcus aureus ST630 clone in China

Methicillin-resistant Staphylococcus aureus (MRSA) sequence type 630 (ST630) is a rarely reported lineage worldwide. This study aimed to trace the dissemination of the emerging MRSA ST630 clones in China and investigate their virulence potential. We collected 22 ST630-MRSA isolates from across China and performed whole-genome sequencing analysis and virulence characterization on these isolates. Epidemiological results showed that MRSA ST630 isolates were primarily isolated from pus/wound secretions, mainly originating from Jiangxi province, and carried diverse virulence and drug resistance genes. Staphylococcal cassette chromosome mec type V (SCCmec V) predominated (11/22, 50.0%) among the MRSA ST630 isolates. Interestingly, nearly half (45.5%) of the 22 ST630-MRSA isolates tested lacked intact SCCmec elements. Phylogenetic analysis demonstrated that ST630-MRSA could be divided into two distinct clades, with widespread dissemination mainly in Chinese regions. Five representative isolates were selected for phenotypic assays, including hemolysin activity, real-time fluorescence quantitative PCR, western blot analysis, hydrogen peroxide killing assay, blood killing assay, cell adhesion and invasion assay, and mouse skin abscess model. The results showed that, compared to the USA300-LAC strain, ST630 isolates exhibited particularly strong invasiveness and virulence in the aforementioned phenotypic assays. This study described the emergence of a highly virulent ST630-MRSA lineage and improved our insight into the molecular epidemiology of ST630 clones in China.IMPORTANCEMethicillin-resistant Staphylococcus aureus (MRSA) sequence type 630 (ST630) is an emerging clone with an increasing isolation rate in China. This study raises awareness of the hypervirulent MRSA ST630 clones in China and alerts people to their widespread dissemination. ST630-staphylococcal cassette chromosome mec V is a noteworthy clone in China, and we present the first comprehensive genetic and phenotypic analysis of this lineage. Our findings provide valuable insights for the prevention and control of infections caused by this emerging MRSA clone.

Whole-genome comparison using complete genomes from Campylobacter fetus strains revealed single nucleotide polymorphisms on non-genomic islands for subspecies differentiation

Introduction

Bovine Genital Campylobacteriosis (BGC), caused by Campylobacter fetus subsp. venerealis, is a sexually transmitted bacterium that significantly impacts cattle reproductive performance. However, current detection methods lack consistency and reliability due to the close genetic similarity between C. fetus subsp. venerealis and C. fetus subsp. fetus. Therefore, this study aimed to utilize complete genome analysis to distinguish genetic features between C. fetus subsp. venerealis and other subspecies, thereby enhancing BGC detection for routine screening and epidemiological studies.

Methods and results

This study reported the complete genomes of four C. fetus subsp. fetus and five C. fetus subsp. venerealis, sequenced using long-read sequencing technologies. Comparative whole-genome analyses (n = 25) were conducted, incorporating an additional 16 complete C. fetus genomes from the NCBI database, to investigate the genomic differences between these two closely related C. fetus subspecies. Pan-genomic analyses revealed a core genome consisting of 1,561 genes and an accessory pangenome of 1,064 genes between the two C. fetus subspecies. However, no unique predicted genes were identified in either subspecies. Nonetheless, whole-genome single nucleotide polymorphisms (SNPs) analysis identified 289 SNPs unique to one or the C. fetus subspecies. After the removal of SNPs located on putative genomic islands, recombination sites, and those causing synonymous amino acid changes, the remaining 184 SNPs were functionally annotated. Candidate SNPs that were annotated with the KEGG "Peptidoglycan Biosynthesis" pathway were recruited for further analysis due to their potential association with the glycine intolerance characteristic of C. fetus subsp. venerealis and its biovar variant. Verification with 58 annotated C. fetus genomes, both complete and incomplete, from RefSeq, successfully classified these seven SNPs into two groups, aligning with their phenotypic identification as CFF (Campylobacter fetus subsp. fetus) or CFV/CFVi (Campylobacter fetus subsp. venerealis and its biovar variant). Furthermore, we demonstrated the application of mraY SNPs for detecting C. fetus subspecies using a quantitative PCR assay.

Discussion

Our results highlighted the high genetic stability of C. fetus subspecies. Nevertheless, Campylobacter fetus subsp. venerealis and its biovar variants encoded common SNPs in genes related to glycine intolerance, which differentiates them from C. fetus subsp. fetus. This discovery highlights the potential of employing a multiple-SNP assay for the precise differentiation of C. fetus subspecies.

Evolutionary dynamics of the successful expansion of pandemic Vibrio parahaemolyticus ST3 in Latin America

The underlying evolutionary mechanisms driving global expansions of pathogen strains are poorly understood. Vibrio parahaemolyticus is one of only two marine pathogens where variants have emerged in distinct climates globally. The success of a Vibrio parahaemolyticus clone (VpST3) in Latin America- the first spread identified outside its endemic region of tropical Asia- provided an invaluable opportunity to investigate mechanisms of VpST3 expansion into a distinct marine climate. A global collection of VpST3 isolates and novel Latin American isolates were used for evolutionary population genomics, pangenome analysis and combined with oceanic climate data. We found a VpST3 population (LatAm-VpST3) introduced in Latin America well before the emergence of this clone in India, previously considered the onset of the VpST3 epidemic. LatAm-VpST3 underwent successful adaptation to local conditions over its evolutionary divergence from Asian VpST3 isolates, to become dominant in Latin America. Selection signatures were found in genes providing resilience to the distinct marine climate. Core genome mutations and accessory gene presences that promoted survival over long dispersals or increased environmental fitness were associated with environmental conditions. These results provide novel insights into the global expansion of this successful V. parahaemolyticus clone into regions with different climate scenarios.

Transmission dynamics of ESBL/AmpC and carbapenemase-producing Enterobacterales between companion animals and humans

Antimicrobial resistance mediated by extended-spectrum beta-lactamase (ESBL)- and plasmid-mediated cephalosporinase (AmpC)-producing Enterobacterales, as well as carbapenemase-producing Enterobacterales have globally increased among companion animals, posing a potential health risk to humans in contact with them. This prospective longitudinal study investigates the transfer of ESBL/AmpC- and carbapenemase-producing Enterobacterales between companion animals and their cohabitant humans in Portugal (PT) and the United Kingdom (UK) during animal infection. Fecal samples and nasal swabs collected from dogs and cats with urinary tract infection (UTI) or skin and soft tissue infection (SSTI), and their cohabitant humans were screened for resistant strains. Relatedness between animal and human strains was established by whole-genome sequencing (WGS). ESBL/AmpC-producing Enterobacterales were detected in companion animals (PT = 55.8%; UK = 36.4%) and humans (PT = 35.9%; UK = 12.5%). Carbapenemase-producing Enterobacterales carriage was observed in one dog from Portugal (2.6%) and another dog from the UK (4.5%). Transmission of index clinical ESBL-producing Escherichia coli and Klebsiella pneumoniae strains to cohabitant humans was observed in three Portuguese households (6.9%, n = 43), with repeated isolation of the index strains on fecal samples from the animals and their cohabiting humans. In addition, longitudinal sharing of E. coli strains carried by companion animals and their owners was observed in other two Portuguese households and two households from the UK. Furthermore, a multidrug-resistant ACT-24-producing Enterobacter hormaechei subsp. hoffmannii strains were also shared within another Portuguese household. These results highlight the importance of the household as an epidemiological unit in the efforts to mitigate the spread of antimicrobial resistance, further emphasizing the need for antimicrobial surveillance in this context, capable of producing data that can inform and evaluate public health actions.

A microbiological and genomic perspective of globally collected Escherichia coli from adults hospitalized with invasive E. coli disease

OBJECTIVES

Escherichia coli can cause infections in the urinary tract and in normally sterile body sites leading to invasive E. coli disease (IED), including bacteraemia and sepsis, with older populations at increased risk. We aimed to estimate the theoretical coverage rate by the ExPEC4V and 9V vaccine candidates. In addition, we aimed at better understanding the diversity of E. coli isolates, including their genetic and phenotypic antimicrobial resistance (AMR), sequence types (STs), O-serotypes and the bacterial population structure.

METHODS

Blood and urine culture E. coli isolates (n = 304) were collected from hospitalized patients ≥60 years (n = 238) with IED during a multicentric, observational study across three continents. All isolates were tested for antimicrobial susceptibility, O-serotyped, whole-genome sequenced and bioinformatically analysed.

RESULTS

A large diversity of STs and of O-serotypes were identified across all centres, with O25b-ST131, O6-ST73 and O1-ST95 being the most prevalent types. A total of 45.4% and 64.7% of all isolates were found to have an O-serotype covered by the ExPEC4V and ExPEC9V vaccine candidates, respectively. The overall frequency of MDR was 37.4% and ST131 was predominant among MDR isolates. Low in-patient genetic variability was observed in cases where multiple isolates were collected from the same patient.

CONCLUSIONS

Our results highlight the predominance of MDR O25b-ST131 E. coli isolates across diverse geographic areas. These findings provide further baseline data on the theoretical coverage of novel vaccines targeting E. coli associated with IED in older adults and their associated AMR levels.

Penicillin-binding protein 3 sequence variations reduce susceptibility of Pseudomonas aeruginosa to β-lactams but inhibit cell division

BACKGROUND

β-lactam antibiotics, which inhibit penicillin-binding protein 3 (PBP3) that is required for cell division, play a key role in treating P. aeruginosa infections. Some sequence variations in PBP3 have been associated with β-lactam resistance but the effects of variations on antibiotic susceptibility and on cell division have not been quantified. Antibiotic efflux can also reduce susceptibility.

OBJECTIVES

To quantify the effects of PBP3 variations on β-lactam susceptibility and cell morphology in P. aeruginosa.

METHODS

Nineteen PBP3 variants were expressed from a plasmid in the reference strain P. aeruginosa PAO1 and genome engineering was used to construct five mutants expressing PBP3 variants from the chromosome. The effects of the variations on β-lactam minimum inhibitory concentration (MIC) and cell morphology were measured.

RESULTS

Some PBP3 variations reduced susceptibility to a variety of β-lactam antibiotics including meropenem, ceftazidime, cefepime and ticarcillin with different variations affecting different antibiotics. None of the tested variations reduced susceptibility to imipenem or piperacillin. Antibiotic susceptibility was further reduced when PBP3 variants were expressed in mutant bacteria overexpressing the MexAB-OprM efflux pump, with some variations conferring clinical levels of resistance. Some PBP3 variations, and sub-MIC levels of β-lactams, reduced bacterial growth rates and inhibited cell division, causing elongated cells.

CONCLUSIONS

PBP3 variations in P. aeruginosa can increase the MIC of multiple β-lactam antibiotics, although not imipenem or piperacillin. PBP3 variations, or the presence of sub-lethal levels of β-lactams, result in elongated cells indicating that variations reduce the activity of PBP3 and may reduce bacterial fitness.

Drug resistance and genomic variations among Mycobacterium tuberculosis isolates from The Nile Delta, Egypt

Tuberculosis is a global public health concern. Earlier reports suggested the emergence of high rates of drug resistant tuberculosis in Egypt. This study included 102 isolates of Mycobacterium tuberculosis collected from two reference laboratories in Cairo and Alexandria. All clinical isolates were sub-cultured on Löwenstein-Jensen medium and analyzed using both BD BACTEC MGIT 960 SIRE Kit and standard diffusion disk assays to identify the antibiotic sensitivity profile. Extracted genomic DNA was subjected to whole genome sequencing (WGS) using Illumina platform. Isolates that belong to lineage 4 represented > 80%, while lineage 3 represented only 11% of the isolates. The percentage of drug resistance for the streptomycin, isoniazid, rifampicin and ethambutol were 31.0, 17.2, 19.5 and 20.7, respectively. Nearly 47.1% of the isolates were sensitive to the four anti-tuberculous drugs, while only one isolate was resistant to all four drugs. In addition, several new and known mutations were identified by WGS. High rates of drug resistance and new mutations were identified in our isolates. Tuberculosis control measures should focus on the spread of mono (S, I, R, E)- and double (S, E)-drug resistant strains present at higher rates throughout the whole Nile Delta, Egypt.

Prevalence and molecular characteristics of intestinal pathogenic Escherichia coli isolated from diarrheal pigs in Southern China

Intestinal pathogenic Escherichia coli (InPEC) is one of the most common causes of bacterial diarrhea in farm animals, including profuse neonatal diarrhea and post weaning diarrhea (PWD) in piglets. In this study, we investigated the prevalence of InPEC and associated primary virulence factors among 543 non-duplicate E. coli isolates from diarrheal pigs from 15 swine farms in southern China. Six major virulence genes associated with InPEC were identified among 69 (12.71 %) E. coli isolates and included est (6.62 %), K88 (4.79 %), elt (3.68 %), eae (1.47 %), stx2 (0.92 %) and F18 (0.55 %). Three pathotypes of InPEC were identified including ETEC (8.10 %), EPEC (1.29 %) and STEC/ETEC (0.92 %). In particular, K88 was only found in ETEC from breeding farms, whereas F18 was only present in STEC/ETEC hybrid from finishing farms. Whole genome sequence analysis of 37 E. coli isolates revealed that InPEC strains frequently co-carried multiple antibiotic resistance gene (ARG). est, elt and F18 were also found to co-locate with ARGs on a single IncFIB/IncFII plasmid. InPEC isolates from different pathotypes also possessed different profiles of virulence genes and antimicrobial resistance genes. Population structure analysis demonstrated that InPEC isolates from different pathotypes were highly heterogeneous whereas those of the same pathotype were extremely similar. Plasmid analysis revealed that K88 and/or est/elt were found on pGX18-2-like/pGX203-2-like and pGX203-1-like IncFII plasmids, while F18 and elt/est, as well as diverse ARGs were found to co-locate on IncFII/IncFIB plasmids with a non-typical backbone. Moreover, these key virulence genes were flanked by or adjacent to IS elements. Our findings indicated that both clonal expansion and horizontal spread of epidemic IncFII plasmids contributed to the prevalence of InPEC and the specific virulence genes (F4, F18, elt and est) in the tested swine farms.

Genomic insights into in-ICU emergence of last-resort antimicrobial resistance in a rare, carbapenem resistant, ST16 Klebsiella pneumoniae strain from Jodhpur, India

OBJECTIVES

To investigate the genomic differences between two extensively drug resistant, ST16 strains of Klebsiella pneumoniae recovered from patients in the same ICU, one of which was colistin resistant.

METHODS

Antimicrobial susceptibilities of the isolates were determined using VITEK-2. Hybrid assemblies for both strains were generated using Oxford Nanopore and Illumina technologies. The sequence type, capsule type, O-locus type, antimicrobial resistance determinants and plasmids carried by the isolates were inferred from the genome sequence. The phylogenetic placement, antimicrobial resistance, and virulence determinants of the isolates relative to a collection (n = 871) of ST16 isolates were assessed.

RESULTS

Both BC16, a colistin-resistant blood stream isolate and U23, a colistin-sensitive urinary isolate displayed near-identical antimicrobial resistance profiles and genome sequences with varying plasmid profiles. The BC16 genome only had 21 SNPs relative to U23 and belonged to the same capsule, O-antigen locus and multi-locus sequence types. The mgrB locus in BC16 was disrupted by an IS5 element. Phylogenetically, U23 and BC16 were placed on a clade with 4 strains belonging to K-type K48 and O-type O2a as opposed to majority (n = 807) of the strains (K-type K51 and O-type O3b).

CONCLUSIONS

BC16 was a colistin resistant derivative of U23, which evolved colistin resistance by an IS5-mediated disruption of the mgrB locus, likely during treatment of the index patient with colistin in the ICU. The strains belong to a rare subtype of ST16 with unique capsular and O-antigen types underscoring the utility of genomic surveillance networks and open-access genomic surveillance data in tracking problem clones.

Draft genome sequences of clinical mastitis-associated Enterococcus faecalis and Enterococcus faecium carrying multiple antimicrobial resistance genes isolated from dairy cows

OBJECTIVES

The emergence of antimicrobial-resistant and mastitis-associated Enterococcus faecalis and Enterococcus faecium is of great concern due to the huge economic losses associated with enterococcal infections. Here we report the draft genome sequences of E. faecalis and E. faecium strains that were isolated from raw milk samples obtained from mastitis-infected cows in Bangladesh.

METHODS

The two strains were isolated, identified, and genomic DNA was sequenced using the Illumina NextSeq 550 platform. The assembled contigs were analysed for virulence, antimicrobial resistance genes, and multilocus sequence type. The genomes were compared to previously reported E. faecalis and E. faecium genomes to generate core genome phylogenetic trees.

RESULTS

E. faecalis strain BR-MHR218Efa and E. faecium strain BR-MHR268Efe belonged to multilocus sequence types ST-190 and ST-22, respectively, both of which appear to represent relatively rare sequence types. BR-MHR268Efe harboured only one antibiotic resistance gene encoding resistance towards macrolides (lsa(A)), while BR-MHR218Efa harboured ten different antibiotic resistance genes encoding resistance to aminoglycosides (ant[6]-Ia, aph(3')-III), sulphonamides (aac(6')-II), lincosamides (lnu(B)), macrolides (erm(B)), MLSB antibiotics (msr(C)), tetracyclines (tet(M), tet(L)), trimethoprim (dfrG), and pleuromutilin-lincosamide-streptogramin A (lsa(E)). Virulence gene composition was different between the two isolates. BR-MHR218Efa harboured only two virulence genes involved in adherence (acm and scm). BR-MHR268Efe harboured eight complete virulence operons including three operons involved in adherence (Ace, Ebp pili, and EfaA), two operons involved in biofilm formation (BopD and Fsr), and three exoenzymes (gelatinase, hyaluronidase, SprE).

CONCLUSIONS

The genome sequences of the strains BR-MHR268Efe and BR-MHR218Efa will serve as a reference point for molecular epidemiological studies of mastitis-associated E. faecalis and E. faecium. Additionally, the findings will help understand the complex antimicrobial-resistance in livestock-assoiated Enterococci.

Occurrence and characteristics of blaOXA-181-carrying Klebsiella aerogenes from swine in China

OBJECTIVES

Klebsiella aerogenes is a largely understudied opportunistic pathogen that can cause sepsis and lead to high mortality rates. In this study, we reported the occurrence of carbapenem-resistant blaOXA-181-carrying Klebsiella aerogenes from swine in China and elucidate their genomic characteristics.

METHODS

A total of 126 samples, including 109 swine fecal swabs, 14 environmental samples, and three feed samples were collected from a pig farm in China. The samples were enriched with LB broth culture and then inoculated into MacConkey agar plates for bacterial isolation. After PCR detection of carbapenemases genes, the blaOXA-181-carrying isolates were subjected to antimicrobial susceptibility testing, and whole-genome sequence analysis.

RESULTS

Four Klebsiella aerogenes isolates carrying the blaOXA-181 gene were obtained from swine faecal samples. All the 4 strains were belonged to ST438. The blaOXA-181 genes were located in IncX3-ColKP3 hybrid plasmids with the core genetic structure of IS26-ΔIS3000-ΔISEcp1-blaOXA-181-ΔlysR-ΔereA-ΔrepA-ISKpn19-tinR-qnrS1-ΔIS2-IS26, which suggests the potential for horizontal transfer and further dissemination of this resistance gene among Enterobacteriaceae and other sources.

CONCLUSIONS

This study represents the first instance of OXA-181-producing K. aerogenes being identified from swine faeces in China. It is crucial to maintain continuous monitoring and ongoing attention to the detection of K. aerogenes carrying blaOXA-181 and other resistance genes in pigs.

Detection and characterisation of a sixth Candida auris clade in Singapore: a genomic and phenotypic study

BACKGROUND

The emerging fungal pathogen Candida auris poses a serious threat to global public health due to its worldwide distribution, multidrug resistance, high transmissibility, propensity to cause outbreaks, and high mortality. We aimed to characterise three unusual C auris isolates detected in Singapore, and to determine whether they constitute a novel clade distinct from all previously known C auris clades (I-V).

METHODS

In this genotypic and phenotypic study, we characterised three C auris clinical isolates, which were cultured from epidemiologically unlinked inpatients at a large tertiary hospital in Singapore. The index isolate was detected in April, 2023. We performed whole-genome sequencing (WGS) and obtained hybrid assemblies of these C auris isolates. The complete genomes were compared with representative genomes of all known C auris clades. To provide a global context, 3651 international WGS data from the National Center for Biotechnology Information (NCBI) database were included in a high-resolution single nucleotide polymorphism (SNP) analysis. Antifungal susceptibility testing was done and antifungal resistance genes, mating-type locus, and chromosomal rearrangements were characterised from the WGS data of the three investigated isolates. We further implemented Bayesian logistic regression models to classify isolates into known clades and simulate the automatic detection of isolates belonging to novel clades as their WGS data became available.

FINDINGS

The three investigated isolates were separated by at least 37 000 SNPs (range 37 000-236 900) from all existing C auris clades. These isolates had opposite mating-type allele and different chromosomal rearrangements when compared with their closest clade IV relatives. The isolates were susceptible to all tested antifungals. Therefore, we propose that these isolates represent a new clade of C auris, clade VI. Furthermore, an independent WGS dataset from Bangladesh, accessed via the NCBI Sequence Read Archive, was found to belong to this new clade. As a proof-of-concept, our Bayesian logistic regression model was able to flag these outlier genomes as a potential new clade.

INTERPRETATION

The discovery of a new C auris clade in Singapore and Bangladesh in the Indomalayan zone, showing a close relationship to clade IV members most commonly found in South America, highlights the unknown genetic diversity and origin of C auris, particularly in under-resourced regions. Active surveillance in clinical settings, along with effective sequencing strategies and downstream analysis, will be essential in the identification of novel strains, tracking of transmission, and containment of adverse clinical effects of C auris infections.

FUNDING

Duke-NUS Academic Medical Center Nurturing Clinician Researcher Scheme, and the Genedant-GIS Innovation Program.

Characteristics and comparative genome analysis of Yersinia enterocolitica and related species associated with human infections in Switzerland 2019-2023

PURPOSE

We aimed to characterise Yersinia enterocolitica from human clinical specimens in Switzerland using epidemiological, microbiological and whole-genome sequencing (WGS) data.

METHODS

Isolates (n = 149) were collected between January 2019 and December 2023. Epidemiological data was noted and strains were characterized by biochemical and serological typing, antimicrobial susceptibility testing (AST), and WGS-based analysis.

RESULTS

Most of the isolates (86%) were from stool specimens and 52% were from male patients. The patients' median age was 28 years (range < 1-94 years). Typing assigned the isolates to bioserotype 4/O:3 (44%), biotype 1A (34%), bioserotype 2/O:9 (21%), and bioserotype 3/O:3 (1%). WGS identified Y. enterocolitica (n = 147), Y. alsatica (n = 1) and Y. proxima (n = 1). Seven isolates were multidrug resistant (MDR) and harboured plasmid pAB829 carrying aph(3″)-Ib, aph(6)-Id, and tet(Y) (n = 1), pAC120 carrying aph(6)-Id and tet(A) (n = 2), or a 12.6 kb Tn2670-like transposon containing catA1, aadA12, sul1, and qacEΔ1 (n = 4). Virulence factors (VFs) included ail (n = 99), invB, (n = 145), ystA (n = 99), ystB (n = 48) and pYV-associated VFs (n = 93). MLST and cgMLST analysis showed that BT 1A strains consisted of several STs and were highly diverse, whereas BT 2/O:9 strains were all ST12 and clustered closely, and BT 4/O:3 strains mostly belonged to ST18 but were more diverse. SNP analysis revealed two highly clonal BT 4/O:3 subpopulations with wide spatio-temporal distribution.

CONCLUSIONS

Y. enterocolitica BT 1A, BT 2/O:9 and BT 4/O:3 are frequently associated with human yersiniosis in Switzerland. WGS-based subtyping of Y. enterocolitica is a powerful tool to explore the genetic diversity and the pathogenic potential of human isolates.

Evolutionary trajectory and characteristics of Mpox virus in 2023 based on a large-scale genomic surveillance in Shenzhen, China

The global epidemic of Mpox virus (MPXV) continues, and a local outbreak has occurred in Shenzhen city since June 2023. Herein, the evolutionary trajectory and characteristics of MPXV in 2023 were analyzed using 92 MPXV sequences from the Shenzhen outbreak and the available genomes from GISAID and GenBank databases. Phylogenetic tracing of the 92 MPXVs suggests that MPXVs in Shenzhen may have multiple sources of importation, and two main transmission chains have been established. The combination of phylogenetic relationships, epidemiological features, and mutation characteristics supports the emergence of a new lineage C.1.1. Together with the B.1 lineage diverging from the A.1 lineage, C.1.1 lineage diverging from the C.1 lineage may serve as another significant evolutionary events of MPXV. Moreover, increasing apolipoprotein B mRNA-editing catalytic polypeptide-like 3 (APOBEC3) related mutations, higher rate of missense mutations, and less mutations in the non-coding regions have been shown during MPXV evolution. Host regulation proteins of MPXV have accumulated considerable amino acid mutations since the B.1 lineage, and a lineage-defining APOBEC3-related mutation that disrupts the N2L gene encoding a viral innate immune modulator has been identified in the C.1.1 lineage. In summary, our study provides compelling evidence for the ongoing evolution of MPXV with specific features.

Adaptive evolution of carbapenem-resistant hypervirulent Klebsiella pneumoniae in the urinary tract of a single patient

The emergence of carbapenem-resistant hypervirulent Klebsiella pneumoniae (CR-hvKp) is a growing concern due to its high mortality and limited treatment options. Although hypermucoviscosity is crucial for CR-hvKp infection, the role of changes in bacterial mucoviscosity in the host colonization and persistence of CR-hvKp is not clearly defined. Herein, we observed a phenotypic switch of CR-hvKp from a hypermucoviscous to a hypomucoviscous state in a patient with scrotal abscess and urinary tract infection (UTI). This switch was attributed to decreased expression of rmpADC, the regulator of mucoid phenotype, caused by deletion of the upstream insertion sequence ISKpn26. Postswitching, the hypomucoid variant showed a 9.0-fold decrease in mice sepsis mortality, a >170.0-fold reduction in the ability to evade macrophage phagocytosis in vitro, and an 11.2- to 40.9-fold drop in growth rate in normal mouse serum. Conversely, it exhibited an increased residence time in the mouse urinary tract (21 vs. 6 d), as well as a 216.4-fold boost in adhesion to bladder epithelial cells and a 48.7% enhancement in biofilm production. Notably, the CR-hvKp mucoid switch was reproduced in an antibiotic-free mouse UTI model. The in vivo generation of hypomucoid variants was primarily associated with defective or low expression of rmpADC or capsule synthesis gene wcaJ, mediated by ISKpn26 insertion/deletion or base-pair insertion. The spontaneous hypomucoid variants also outcompeted hypermucoid bacteria in the mouse urinary tract. Collectively, the ISKpn26-associated mucoid switch in CR-hvKp signifies the antibiotic-independent host adaptive evolution, providing insights into the role of mucoid switch in the persistence of CR-hvKp.

The integrated genomic surveillance system of Andalusia (SIEGA) provides a One Health regional resource connected with the clinic

The One Health approach, recognizing the interconnectedness of human, animal, and environmental health, has gained significance amid emerging zoonotic diseases and antibiotic resistance concerns. This paper aims to demonstrate the utility of a collaborative tool, the SIEGA, for monitoring infectious diseases across domains, fostering a comprehensive understanding of disease dynamics and risk factors, highlighting the pivotal role of One Health surveillance systems. Raw whole-genome sequencing is processed through different species-specific open software that additionally reports the presence of genes associated to anti-microbial resistances and virulence. The SIEGA application is a Laboratory Information Management System, that allows customizing reports, detect transmission chains, and promptly alert on alarming genetic similarities. The SIEGA initiative has successfully accumulated a comprehensive collection of more than 1900 bacterial genomes, including Salmonella enterica, Listeria monocytogenes, Campylobacter jejuni, Escherichia coli, Yersinia enterocolitica and Legionella pneumophila, showcasing its potential in monitoring pathogen transmission, resistance patterns, and virulence factors. SIEGA enables customizable reports and prompt detection of transmission chains, highlighting its contribution to enhancing vigilance and response capabilities. Here we show the potential of genomics in One Health surveillance when supported by an appropriate bioinformatic tool. By facilitating precise disease control strategies and antimicrobial resistance management, SIEGA enhances global health security and reduces the burden of infectious diseases. The integration of health data from humans, animals, and the environment, coupled with advanced genomics, underscores the importance of a holistic One Health approach in mitigating health threats.

Bacterial homologs of innate eukaryotic antiviral defenses with anti-phage activity highlight shared evolutionary roots of viral defenses

Prokaryotes have evolved a multitude of defense systems to protect against phage predation. Some of these resemble eukaryotic genes involved in antiviral responses. Here, we set out to systematically project the current knowledge of eukaryotic-like antiviral defense systems onto prokaryotic genomes, using Pseudomonas aeruginosa as a model organism. Searching for phage defense systems related to innate antiviral genes from vertebrates and plants, we uncovered over 450 candidates. We validated six of these phage defense systems, including factors preventing viral attachment, R-loop-acting enzymes, the inflammasome, ubiquitin pathway, and pathogen recognition signaling. Collectively, these defense systems support the concept of deep evolutionary links and shared antiviral mechanisms between prokaryotes and eukaryotes.

Genomic diversity of Salmonella enterica serovar Typhimurium isolated from chicken processing facilities in New South Wales, Australia

Contamination of poultry products by Salmonella enterica serovar Typhimurium (STm) is a major cause of foodborne infections and outbreaks. This study aimed to assess the diversity and antimicrobial resistance (AMR) carriage of STm in three chicken processing plants using genomic sequencing. It also aimed to investigate whether any particular strain types were associated with cases of human illness. Multilevel genome typing (MGT) was used to analyze 379 STm isolates from processed chicken carcasses. The diversity of chicken STm sequence types (STs) increased from MGT1 (2 STs) to MGT9 (257 STs). STs at MGT5 to MGT9 levels that were unique to one processing plant and shared among the processing plants were identified, likely reflecting the diversity of STm at their farm source. Fifteen medium resolution MGT5 STs matched those from human infections in Australia and globally. However, no STs matched between the chicken and human isolates at high resolution levels (MGT8 or MGT9), indicating the two STm populations were phylogenetically related but were unlikely to be directly epidemiologically linked. AMR genes were rare, with only a bla TEM-1 gene carried by a 95 kb IncI1 Alpha plasmid being identified in 20 isolates. In conclusion, subpopulations that were widespread in processing plants and had caused human infections were described using MGT5 STs. In this STM population, AMR was rare with only sporadic resistance to a single drug class observed. The genomic analysis of STm from chicken processing plants in this study provided insights into STm that contaminate meat chickens early in the food production chain.

OrthoPhyl-streamlining large-scale, orthology-based phylogenomic studies of bacteria at broad evolutionary scales

There are a staggering number of publicly available bacterial genome sequences (at writing, 2.0 million assemblies in NCBI's GenBank alone), and the deposition rate continues to increase. This wealth of data begs for phylogenetic analyses to place these sequences within an evolutionary context. A phylogenetic placement not only aids in taxonomic classification but informs the evolution of novel phenotypes, targets of selection, and horizontal gene transfer. Building trees from multi-gene codon alignments is a laborious task that requires bioinformatic expertise, rigorous curation of orthologs, and heavy computation. Compounding the problem is the lack of tools that can streamline these processes for building trees from large-scale genomic data. Here we present OrthoPhyl, which takes bacterial genome assemblies and reconstructs trees from whole genome codon alignments. The analysis pipeline can analyze an arbitrarily large number of input genomes (>1200 tested here) by identifying a diversity-spanning subset of assemblies and using these genomes to build gene models to infer orthologs in the full dataset. To illustrate the versatility of OrthoPhyl, we show three use cases: E. coli/Shigella, Brucella/Ochrobactrum and the order Rickettsiales. We compare trees generated with OrthoPhyl to trees generated with kSNP3 and GToTree along with published trees using alternative methods. We show that OrthoPhyl trees are consistent with other methods while incorporating more data, allowing for greater numbers of input genomes, and more flexibility of analysis.

FosA3 emerging in clinical carbapenemase-producing C. freundii

Fosfomycin (FOS) is an effective antibiotic against multidrug-resistant Enterobacterales, but its effectiveness is reducing. Little is known on the current prevalence of FosA enzymes in low-risk pathogens, such as Citrobacter freundii. The aim of the study was the molecular characterization of a carbapenemase- and FosA-producing C. freundii collected in Italy. AK867, collected in 2023, showed an XDR profile, retaining susceptibility only to colistin. AK867 showed a FOS MIC >128 mg/L by ADM. Based on WGS, AK867 belonged to ST116 and owned a wide resistome, including fosA3, blaKPC-2, and blaVIM-1. fosA3 was carried by a conjugative pKPC-CAV1312 plasmid of 320,480 bp, on a novel composite transposon (12,907 bp). FosA3 transposon shared similarities with other fosA3-harboring pKPC-CAV1312 plasmids among Citrobacter spp. We report the first case of FosA3 production in clinical carbapenemase-producing C. freundii ST116. The incidence of FosA3 enzymes is increasing among Enterobacterales, affecting even low-virulence pathogens, as C. freundii.

Methicillin-resistant Staphylococcus aureus outbreak in a Dutch equine referral clinic

In 2020 and 2022, nine cases of surgical site infections with a methicillin-resistant Staphylococcus aureus (MRSA) were diagnosed in horses in an equine referral clinic. Sixteen isolates (horses, n=9; environment, n=3; and staff members, n=4) were analysed retrospectively using Nanopore whole-genome sequencing to investigate the relatedness of two suspected MRSA outbreaks (2020 and 2022). The MRSA isolates belonged to ST398 and ST612. ST398 genomes from 2020 and 2022 formed three phylogenetic clusters. The first ST398 cluster from 2020 consisted of isolates from five horses and one staff member, and we suspected within clinic transmission. The second cluster of ST398 isolates from 2022 originated from two horses and two staff members but showed higher single nucleotide polymorphism (SNP) distances. One ST398 isolate from an individual staff member was not related to the other two clusters. The ST612 isolates were isolated in 2022 from two horses and three environmental samples and showed very low SNP distances (<7 SNPs), indicating the transmission of MRSA ST612 in this clinic in 2022. Molecular characterization revealed an abundant set of virulence genes and plasmids in the ST612 isolates in comparison to ST398 isolates. Phenotypic antimicrobial susceptibility showed that differences between the two sequence types were consistent with the genetic characteristics. MRSA ST612 has not been reported in Europe before, but it is a dominant clone in African hospitals and has been described in horses and people working with horses in Australia, indicating the importance of surveillance.

PanVA: Pangenomic Variant Analysis

Genomics researchers increasingly use multiple reference genomes to comprehensively explore genetic variants underlying differences in detectable characteristics between organisms. Pangenomes allow for an efficient data representation of multiple related genomes and their associated metadata. However, current visual analysis approaches for exploring these complex genotype-phenotype relationships are often based on single reference approaches or lack adequate support for interpreting the variants in the genomic context with heterogeneous (meta)data. This design study introduces PanVA, a visual analytics design for pangenomic variant analysis developed with the active participation of genomics researchers. The design uniquely combines tailored visual representations with interactions such as sorting, grouping, and aggregation, allowing users to navigate and explore different perspectives on complex genotype-phenotype relations. Through evaluation in the context of plants and pathogen research, we show that PanVA helps researchers explore variants in genes and generate hypotheses about their role in phenotypic variation.

Genome sequence data of Saccharomyces cerevisiae CBS 493.94

Whole genome sequencing (WGS) and data concerning identity and safety for Saccharomyces cerevisiae CBS 493.94 are reported. This strain was isolated from a British brewery in 1958 and deposited at the CBS culture collection Westerdijk Fungal Biodiversity Institute under the accession number CBS 493.94. The long-reads sequencing data, obtained via PacBio Sequel, and short-reads data, via Illumina NovaSeq 6000, were deposited at NCBI under accession number PRJNA1044661. The hybrid assembly was made publicly available via Zenodo and NCBI. For strain identification, data from 18S rRNA, ANI dendrogram and Core Genome single nucleotide polymorphism (SNP) Tree showed that the present isolate belongs to the genus Saccharomyces, species cerevisiae. The potential genes of concern, e.g. antimycotic resestance genes, were not detected. This strain is commonly used as a feed additive for animal health improvement and the present data summarise the unambiguous identity and strain's FKS1 gene does not code for any amino acid variants of concern.

Phylogenomic Analyses of Three Distinct Lineages Uniting Staphylococcus cohnii and Staphylococcus urealyticus from Diverse Hosts

We sequenced and assembled genomes for 17 isolates of Staphylococcus cohnii isolated from osteomyelitis lesions in young broilers from two separate experiments where we induced lameness using a hybrid wire-litter flooring system. Whole genome comparisons using three different methods support a close relationship of genomes from both S. cohnii and Staphylococcus urealyticus. The data support three different lineages, which we designated as Lineage 1, Lineage 2, and Lineage 3, uniting these two species within an evolving complex. We present evidence for horizontal transfer between lineages of genomic regions from 50-440 kbp. The transfer of a 186 kbp region from Lineage 1 to Lineage 2 appears to have generated Lineage 3. Human-associated isolates appear to be limited to Lineages 2 and 3 but Lineage 2 appears to contain a higher number of human pathogenic isolates. The chicken isolates from our lameness trials included genomically diverse isolates from both Lineage 1 and 2, and isolates from both lineages were obtained from osteomyelitis lesions of individual birds. Our results expand the diversity of Staphylococci associated with osteomyelitis in poultry and suggest a high diversity in the microbiome of day-old chicks. Our data also support a reevaluation and unification of the taxonomic classifications of S. cohnii and S. urealyticus.

Evidence of Reduced Virulence and Increased Colonization Among Pneumococcal Isolates of Serotype 3 Clade II Lineage in Mice

Recent phylogenetic profiling of pneumococcal serotype 3 (Pn3) isolates revealed a dynamic interplay among major lineages with the emergence and global spread of a variant termed clade II. The cause of Pn3 clade II dissemination along with epidemiological and clinical ramifications are currently unknown. Here, we sought to explore biological characteristics of dominant Pn3 clades in a mouse model of pneumococcal invasive disease and carriage. Carriage and virulence potential were strain dependent with marked differences among clades. We found that clinical isolates from Pn3 clade II are less virulent and less invasive in mice compared to clade I isolates. We also observed that clade II isolates are carried for longer and at higher bacterial densities in mice compared to clade I isolates. Taken together, our data suggest that the epidemiological success of Pn3 clade II could be related to alterations in the pathogen's ability to cause invasive disease and to establish a robust carriage episode.

Pangenome and genomic signatures linked to the dominance of the lineage-4 of Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

BACKGROUND

The lineage 4 (L4) of Mycobacterium tuberculosis (MTB) is not only globally prevalent but also locally dominant, surpassing other lineages, with lineage 2 (L2) following in prevalence. Despite its widespread occurrence, factors influencing the expansion of L4 and its sub-lineages remain poorly understood both at local and global levels. Therefore, this study aimed to conduct a pan-genome and identify genomic signatures linked to the elevated prevalence of L4 sublineages among extrapulmonary TB (EPTB) patients in western Ethiopia.

METHODS

A cross-sectional study was conducted at an institutional level involving confirmed cases of extrapulmonary tuberculosis (EPTB) patients from August 5, 2018, to December 30, 2019. A total of 75 MTB genomes, classified under lineage 4 (L4), were used for conducting pan-genome and genome-wide association study (GWAS) analyses. After a quality check, variants were identified using MTBseq, and genomes were de novo assembled using SPAdes. Gene prediction and annotation were performed using Prokka. The pan-genome was constructed using GET_HOMOLOGUES, and its functional analysis was carried out with the Bacterial Pan-Genome Analysis tool (BPGA). For GWAS analysis, Scoary was employed with Benjamini-Hochberg correction, with a significance threshold set at p-value ≤ 0.05.

RESULTS

The analysis revealed a total of 3,270 core genes, predominantly associated with orthologous groups (COG) functions, notably in the categories of '[R] General function prediction only' and '[I] Lipid transport and metabolism'. Conversely, functions related to '[N] Cell motility' and '[Q] Secondary metabolites biosynthesis, transport, and catabolism' were primarily linked to unique and accessory genes. The pan-genome of MTB L4 was found to be open. Furthermore, the GWAS study identified genomic signatures linked to the prevalence of sublineages L4.6.3 and L4.2.2.2.

CONCLUSIONS

Apart from host and environmental factors, the sublineage of L4 employs distinct virulence factors for successful dissemination in western Ethiopia. Given that the functions of these newly identified genes are not well understood, it is advisable to experimentally validate their roles, particularly in the successful transmission of specific L4 sublineages over others.

Tracing the possible evolutionary trends of Morganella morganii: insights from molecular epidemiology and phylogenetic analysis

Morganella morganii, encompassing two subspecies, subsp. morganii and subsp. sibonii, is a common opportunistic pathogen, notable for intrinsic resistance to multiple antimicrobial agents. Despite its clinical significance, research into the potential evolutionary dynamics of M. morganii remains limited. This study involved the analysis of genome sequences from 431 M. morganii isolates, comprising 206 isolates that cause host infections, obtained from this study and 225 from the NCBI genome data sets. A diverse array of antimicrobial resistance genes (ARGs) was identified in M. morganii isolates, including mcr-1, tet(X4), tmexCD-toprJ, and various carbapenemase genes. In addition, a novel blaKPC-2-bearing plasmid with demonstrated conjugative capability was discovered in M. morganii. The majority of virulence-related genes (VRGs), except for the hlyCABD gene cluster, were found in almost all M. morganii. Three novel genospecies of M. morganii were identified, designated as M. chanii, M. variant1, and M. variant2. Compared to M. sibonii, M. chanii genospecies possessed a greater number of flagellar-related genes, typically located within mobile genetic elements (MGEs), suggesting potential for better environmental adaptability. Phylogenetic analysis further disclosed that M. morganii was divided into 12 sequence clusters (SCs). Particularly, SC9 harbored an elevated abundance of ARGs and VRGs, mainly toxin-related genes, and was associated with a higher presence of MGEs compared to non-SC9 strains. The collective findings suggest that M. morganii undergoes evolution driven by the influence of MGEs, thereby significantly enhancing its adaptability to selective pressures of environmental changes and clinical antimicrobial agents.IMPORTANCEThe growing clinical significance of Morganella morganii arises from its abundant virulence factors and antimicrobial resistance genes, resulting in elevated infection rates and increased clinical scrutiny. However, research on the molecular epidemiology and evolutionary trends of M. morganii has been scarce. Our study established a list of virulence-related genes (VRGs) for M. morganii and conducted a large-scale epidemiological investigation into these VRGs. Based on genomic classification, three novel genotypes of M. morganii were identified, representing evolutionary adaptations and responses to environmental challenges. Furthermore, we discovered the emergence of a sequence cluster enriched with antimicrobial resistance genes, VRGs, and mobile genetic elements, attributed to the selective pressure of antimicrobial agents. In addition, we identified a novel conjugative plasmid harboring the blaKPC-2 gene. These findings hold significance in monitoring and comprehending the epidemiology of M. morganii.

Geographic EBV variants confound disease-specific variant interpretation and predict variable immune therapy responses

ABSTRACT

Epstein-Barr virus (EBV) is a potent carcinogen linked to hematologic and solid malignancies and causes significant global morbidity and mortality. Therapy using allogeneic EBV-specific lymphocytes shows promise in certain populations, but the impact of EBV genome variation on these strategies remains unexplored. To address this, we sequenced 217 EBV genomes, including hematologic malignancies from Guatemala, Peru, Malawi, and Taiwan, and analyzed them alongside 1307 publicly available EBV genomes from cancer, nonmalignant diseases, and healthy individuals across Africa, Asia, Europe, North America, and South America. These included, to our knowledge, the first natural killer (NK)/T-cell lymphoma (NKTCL) EBV genomes reported outside of East Asia. Our findings indicate that previously proposed EBV genome variants specific to certain cancer types are more closely tied to geographic origin than to cancer histology. This included variants previously reported to be specific to NKTCL but were prevalent in EBV genomes from other cancer types and healthy individuals in East Asia. After controlling for geographic region, we did identify multiple NKTCL-specific variants associated with a 7.8-fold to 21.9-fold increased risk. We also observed frequent variations in EBV genomes that affected peptide sequences previously reported to bind common major histocompatibility complex alleles. Finally, we found several nonsynonymous variants spanning the coding sequences of current vaccine targets BALF4, BKRF2, BLLF1, BXLF2, BZLF1, and BZLF2. These results highlight the need to consider geographic variation in EBV genomes when devising strategies for exploiting adaptive immune responses against EBV-related cancers, ensuring greater global effectiveness and equity in prevention and treatment.