Hierarchical and spatially explicit clustering of DNA sequences with BAPS software

Comparative genomics of Salmonella enterica serovars Paratyphi A, Typhi and Typhimurium reveals distinct profiles of their pangenome, mobile genetic elements, antimicrobial resistance and defense systems repertoire

Salmonella enterica (S. enterica) is a highly ubiquitous and diverse animal and human pathogen. Distinct S. enterica serovars may present varying host-specificity and cause different diseases. While the human-restricted serovars S. Typhi (STY) and S. Paratyphi A (SPA) cause in humans a systemic life-threatening enteric fever, the host-generalist serovar, S. Typhimurium (STM) causes in immunocompetent individuals a self-limited gastroenteritis. Here, we have performed whole-genome sequencing and hybrid assembly of new SPA and STY typhoidal strains and took a comparative genomics approach to examine their phylogeny, pangenome structure and accessory genome content in comparison to the reference non-typhoidal serovar, STM. Our results identified previously uncharacterized lineages of SPA and refined the presence and distribution of core pseudogenes in typhoidal serovars. Pangenome analysis showed that while these serovars have a relatively similar core-genome size, the accessory genome of STM is more than four times larger than those of typhoidal Salmonellae and that STY and SPA display a more closed pangenome than STM. Unexpectedly, we demonstrate that STY and SPA present distinct differences in their pangenome composition, with a noticeable lower number of prophages, conjugative elements and antimicrobial genes per genome in SPA vs. STY. These results suggest that although SPA and STY are closely related at the DNA level, share a similar lifestyle and cause a symptomatic-indistinguishable disease, their genomic evolution and accessory genomes are markedly different. Moreover, these results may provide genomic explanation to phenotypic and epidemiological differences in antimicrobial resistance profiles associated with these serovars globally.

Integrative and conjugative elements associated with antimicrobial resistance in multidrug resistant Pasteurella multocida isolates from bovine respiratory disease (BRD)-affected animals in Spanish feedlots

The emergence of multidrug-resistance (MDR) in Pasteurella multocida, a major contributor to bovine respiratory disease (BRD) is being increasingly reported, often linked to the carriage of antimicrobial resistance genes (ARGs) on integrative and conjugative elements (ICEs). The resistance phenotype for 19 antimicrobials was determined using broth microdilution in 75 Pasteurella multocida isolates from healthy and BRD-affected cattle from five feedlots. The genomes of 32 isolates were sequenced to identify ARG) and mobile genetic elements (MGEs) and assess their genetic diversity. MDR isolates (with phenotypic resistance to aminoglycosides, macrolides, fluoroquinolones and/or tetracyclines) were primarily found among BRD-affected compared to healthy animals. Non-susceptible isolates, belonging to ST79 and ST13, harbored point mutations and four to nine ARGs, including rarely reported mechanisms in Europe (mph(E), msr(E) and aadA31 ARGs and newly described mutations in the gyrA/parC genes). All ARGs were linked to the presence of MGEs including two ICEs, Tn7407 and the novel Tn7809, a prophage and a putative composite transposon. Clonally related isolates were found in different batches from the same feedlot, suggesting maintenance of MDR strains. Our findings demonstrate the diverse genetic basis of AMR in P. multocida from BRD-affected cattle in Spain, emphasizing the role of MGEs in the ARG dissemination.

Genetic landscape and evolution of Acinetobacter pittii, an underestimated emerging nosocomial pathogen

As a member of Acinetobacter calcoaceticus-baumannii complex, Acinetobacter pittii has been an emerging concern in nosocomial infection due to its increasing prevalence and multidrug resistance (MDR). However, its population structure remains broadly unknown, hampering efficient tracing of its transmission and evolution. In this study, we developed a distributed core genome multilocus sequence typing (dcgMLST) for A. pittii based on 750 genomes and employed it to map the genetic landscape and evolution of A. pittii. The results demonstrated that two hierarchical clustering (HC) levels effectively correspond to genetic diversity from species (HC1100) to natural populations (HC450), as well as that a predominant lineage, HC1100_4, accounts for 33.9% of A. pittii strains. Subsequent analysis revealed that specific gene gain and loss events within HC1100_4 are linked to adaptations to environmental stress. Moreover, we identified a cluster of multidrug-resistant plasmids PT_712 responsible for the dissemination of blaNDM-1 genes within the genus of Acinetobacter. This study provides a framework for characterizing genetic diversity, evolutionary dynamics, molecular population distribution, and tracing of A. pittii, which has the potential to improve infection control strategies and public health policy.

Molecular epidemiology of New Delhi metallo-β-lactamase-producing Escherichia coli in retail market chickens, Shandong, China

Background

The global spread of carbapenem-resistant Escherichia coli is a major public health concern. An investigation of their presence in the human and food chain products would facilitate the elucidation of the route of their food-borne transmission. Thus, the aim of this study was to investigate the prevalence of NDM-positive E. coli isolates in chicken at retail markets in Shandong, China.

Methods

A total of 60 NDM-positive isolates were recovered from 531 E. coli isolates obtained from chickens at the retail market in Shandong. Antimicrobial susceptibility testing and polymerase chain reaction screening were performed to investigate the phenotype and genotype of carbapenemase resistance. Genomic characteristics of the -producing isolates were determined by WGS and bioinformatic analysis.

Results

All of these isolates were multidrug-resistant (MDR), with a majority exhibiting resistance to meropenem, ampicillin, ceftazidime, cefotaxime, florfenicol, sulfamethoxazole/trimethoprim, and tetracycline. Whole genome sequencing (WGS) analysis indicated that these isolates were belonged to 18 distinct sequence types (STs), with the most prevalent STs being ST515 (17/60) and ST69 (11/60). Additionally, WGS analysis revealed that clonal spread of NDM-positive ST69 and ST515 E. coli isolates at markets in different cities in Shandong. Phylogenomic analysis showed that NDM-positive E. coli isolates from chickens were closely related to those of human origin.

Conclusion

This study provides a new insight into the spread of NDM-positive E. coli isolates from retail chicken, and offers essential data for public health management.

Molecular epidemiology of Salmonella Enteritidis in humans and animals in Spain

Salmonella Enteritidis, the most prevalent serovar-causing human gastroenteritis, has been traditionally linked to poultry sources. Although antimicrobial resistance (AMR) is not common in this serovar, increasing levels of resistance to fluoroquinolones and ampicillin have been reported in the last few years. Here, 298 isolates retrieved from different sources (human, livestock, wildlife, food, and environment) and years (2002-2021) in Spain were analyzed to evaluate their diversity, the distribution of AMR-conferring genes (ARGs), and mutations and reconstruct the epidemiology of infection due to this serovar. Isolates were clustered in two major clades (I and II), with strains in clade I (including 61.5% of all human isolates) displaying a pan-susceptible phenotype and not carrying AMR determinants. In contrast, clade II included 80.7% of isolates from animal/food/environmental sources, with the majority (69.8%) harboring mutations in the quinolone resistance determinant regions (QRDR). ARGs, although rare, were mostly found in clade II strains that also carried plasmid replicons, among which IncX1 was the most common. Although higher levels of phenotypic resistance were found in animal isolates, extended-spectrum beta-lactamase, plasmid-mediated AmpC, and carbapenemase-encoding genes were only found among human isolates. In summary, the majority of human and animal isolates from Spanish sources in our collection were classified in different phylogenetic branches, suggesting that additional sources are contributing to the occurrence of foodborne infections in Spain. Furthermore, the different distributions of virulence factors and ARGs in isolates from different sources and their association with specific plasmids suggest the presence of different dynamics contributing to the selection of resistant strains.

Diversification of an emerging bacterial plant pathogen; insights into the global spread of Xanthomonas euvesicatoria pv. perforans

Emerging and re-emerging plant diseases continue to present multifarious threats to global food security. Considerable recent efforts are therefore being channeled towards understanding the nature of pathogen emergence, their spread and evolution. Xanthomonas euvesicatoria pv. perforans (Xep), one of the causal agents of bacterial spot of tomato, rapidly emerged and displaced other bacterial spot xanthomonads in many tomato production regions around the world. In less than three decades, it has become a dominant xanthomonad pathogen in tomato production systems across the world and presents a compelling example for understanding diversification of recently emerged bacterial plant pathogens. Although Xep has been continuously monitored in Florida since its discovery, the global population structure and evolution at the genome-scale is yet to be fully explored. The objectives of this work were to determine genetic diversity globally to ascertain if different tomato production regions contain genetically distinct Xep populations, to examine genetic relatedness of strains collected in tomato seed production areas in East Asia and other production regions, and to evaluate variation in type III secretion effectors, which are critical pathogenicity and virulence factors, in relationship to population structure. We used genome data from 270 strains from 13 countries for phylogenetic analysis and characterization of type III effector gene diversity among strains. Our results showed notable genetic diversity in the pathogen. We found genetically similar strains in distant tomato production regions, including seed production regions, and diversification over the past 100 years, which is consistent with intercontinental dissemination of the pathogen in hybrid tomato production chains. Evolution of the Xep pangenome, including the acquisition and loss of type III secreted effectors, is apparent within and among phylogenetic lineages. The apparent long-distance movement of the pathogen, together with variants that may not yet be widely distributed, poses risks of emergence of new variants in tomato production.

Whole-genome-based characterization of Escherichia albertii strains isolated from paediatric diarrhoeal cases in Kolkata, India

Escherichia albertii is a Gram-negative facultative anaerobic bacterium that causes diarrhoea in humans. This study shows the isolation of E. albertii from hospitalized paediatric diarrhoeal cases and genome-based characteristics with putative virulence factors and antimicrobial resistance. E. albertii isolates were identified by species-specific PCR, targeting the gene encoding cytolethal distending toxin (Ea-cdt). The genome of E. albertii was sequenced to identify (i) genes encoding virulence factors (ii) antibiotic resistance-encoding genes, including the mobile genetic elements and (iii) core gene-based phylogenetic relationships and pan-genome features. A total of 10 (1.2%) E. albertii isolates were isolated from 854 faecal samples, of which 6 (60%) were found as the sole pathogen and the remaining 4 (40%) were identified along with other pathogens, such as enteroaggregative Escherichia coli, rotavirus and adenovirus. Patients from whom E. albertii was isolated presented cholera-like diarrhoea, i.e. with watery stool (60%) with moderate dehydration (100%), fever (20%) and abdominal pain (20%). The antimicrobial susceptibility testing of E. albertii showed that most of the isolates were susceptible or reduced susceptible to most of the antibiotics except resistance to erythromycin (80%), tetracycline (50%), nalidixic acid (40%), ampicillin (40%), doxycycline (30%) and ceftriaxone (20%). In the whole-genome sequence, E. albertii isolates revealed several virulence-encoding genes, namely the intimin (eae, E. coli attaching and effacing), the cytolethal distending toxin type II subunit A (cdt-IIA), adhesion (paa, porcine attaching- and effacing-associated), non-LEE (locus of enterocyte effacement) encoded effector A (nleA) and antimicrobial resistance genes (ARGs) conferring resistance to tetracycline (tetA, tetR), sulphonamides (sul2), fluoroquinolones (qnrS) and beta-lactamases (bla CTX-M, blaTEM). The SNP-based phylogenetic analysis of 647 whole genomes of E. albertii isolates from the National Center for Biotechnology Information databases did not reveal any comparable clustering pattern based on the biological source and place of isolation. The genome of some of the E. albertii was closely related to those of the isolates from China and the United Kingdom. The PFGE patterns revealed that most of the E. albertii isolates were distinct clones. This study reports on the extensive genome analysis of diarrhoea-associated E. albertii harbouring multiple virulence and ARGs.

skDER & CiDDER: two scalable approaches for microbial genome dereplication

An abundance of microbial genomes have been sequenced in the past two decades. For fundamental comparative genomic investigations, where the goal is to determine the major gain and loss events shaping the pangenome of a species, it is often unnecessary and computationally onerous to include all available genomes in studies. In addition, over-representation of specific lineages due to sampling and sequencing bias can have undesired effects on evolutionary analyses. To assist users with genomic dereplication, selecting a subset of representative genomes, for downstream comparative genomic investigations, we developed skDER & CiDDER (https://github.com/raufs/skDER). skDER combines recent advances to efficiently estimate average nucleotide identity (ANI) between thousands of microbial genomes with two efficient algorithms for genomic dereplication. Further, CiDDER implements an approach whereby protein clusters are determined across all genomes and genomes are iteratively selected as representatives until a user-defined saturation of the total protein space is achieved. To support ease of use, several auxiliary functionalities are implemented within the two programs, including arguments to: (i) test the number of representative genomes resulting from a variety of clustering parameters, (ii) automate downloading of genomes belonging to a bacterial species or genus by name, (iii) cluster non-representative genomes to their closest representative genomes, and (iv) automatically filter predicted plasmids and phages prior to dereplication. We further assess the effects of filtering mobile genetic elements (MGEs) on ANI and alignment fraction (AF) estimates between pairs of genomes and find that MGEs tend to slightly deflate both metrics in one species.

Avian-specific Salmonella transition to endemicity is accompanied by localized resistome and mobilome interaction

Bacterial regional demonstration after global dissemination is an essential pathway for selecting distinct finesses. However, the evolution of the resistome during the transition to endemicity remains unaddressed. Using the most comprehensive whole-genome sequencing dataset of Salmonella enterica serovar Gallinarum (S. Gallinarum) collected from 15 countries, including 45 newly recovered samples from two related local regions, we established the relationship among avian-specific pathogen genetic profiles and localization patterns. Initially, we revealed the international transmission and evolutionary history of S. Gallinarum to recent endemicity through phylogenetic analysis conducted using a spatiotemporal Bayesian framework. Our findings indicate that the independent acquisition of the resistome via the mobilome, primarily through plasmids and transposons, shapes a unique antimicrobial resistance profile among different lineages. Notably, the mobilome-resistome combination among distinct lineages exhibits a geographical-specific manner, further supporting a localized endemic mobilome-driven process. Collectively, this study elucidates resistome adaptation in the endemic transition of an avian-specific pathogen, likely driven by the localized farming style, and provides valuable insights for targeted interventions.

Salmonella Infantis outbreak on six broiler units in Great Britain: investigation, epidemiology, and control

AIMS

To describe the analysis, epidemiology, and control of six contemporaneous and linked outbreaks of Salmonella enterica subsp. enterica serovar Infantis on British broiler farms. Salmonella Infantis is a potentially multidrug-resistant foodborne zoonosis and can persistently colonize poultry flocks and farms.

METHODS AND RESULTS

Routine monitoring initially identified the organism, which was tracked to six farms associated with a single company. Extensive, repeat sampling identified widespread and, in some cases, persistent contamination. Salmonella Infantis was also isolated from three associated processing factories and catching crew equipment, but not from associated hatcheries and feed mills. Whole genome sequencing and resistance phenotyping revealed one strain was present in the processing plants and on five farms. However, on one of those farms, several highly genetically distinct strains were also detected, including one also found in one of the processing plants. The sixth farm had a strain that was genetically unrelated to strains collected from the other premises and which exhibited an extended spectrum beta-lactamase phenotype. Cleaning and disinfection were enhanced, and the organism was eventually cleared from all farms.

CONCLUSIONS

There were multiple incursions of varied strains, with a possible link to processing factories. Elimination of S. Infantis from premises can be challenging but achievable.

A genomic and phenotypic investigation of pigeon-adaptive Salmonella

Salmonella, a significant threat to public safety, inflicts substantial economic losses on the poultry industry. The unique "parental feeding" breeding model of pigeon farms, against the "all-in & all-out" biosecurity strategy, makes them susceptible to Salmonella infections and subsequent outbreaks of pigeon paratyphoid. This study initially studied three pigeon paratyphoid outbreak incidents in Henan, China, in which 53 strains of pigeon-origin Salmonella Typhimurium (STM) were identified. Whole-genome sequencing (WGS) and antimicrobial-resistant profile analysis revealed that the three outbreaks were caused by distinct STM clones (ST128-DT2, ST19-DT99). Global phylogenetic analysis suggested that the United States is a possible origin, indicating a risk of intercontinental transmission via pigeon eggs. Further bacterial virulence and invasion assays, including in vitro and in vivo assays, revealed that pigeon-host-adaptive STM, compared to broad-host-range STM, carried fewer resistance genes, exhibited higher invasion indices and pseudogene levels, displayed a non-rdar (red dry and rough) phenotype, and had strong biofilm formation capability. Additionally, they showed reduced virulence and invasiveness in mice but a pigeon-adaptive feature in cogent models. The collective results support the host adaptation for pigeons among DT2 and DT99 phage-type isolates.

An MDR Salmonella Enteritidis sublineage associated with gastroenteritis outbreaks and invasive disease in China

OBJECTIVES

Salmonella enterica serovar Enteritidis (S. Enteritidis) is a commonly reported pathogen which adapts to multiple hosts and causes critical disease burden at a global level. Here, we investigated a recently derived epidemic sublineage with multidrug resistance (MDR), which have caused extended time-period and cross-regional gastroenteritis outbreaks and even invasive nontyphoidal Salmonella disease (iNTS) in China.

METHODS

Whole-genome sequencing and antimicrobial resistance (AMR) testing were applied to 729 Chinese S. Enteritidis isolates in relation to gastroenteritis outbreaks, gastrointestinal-sporadic and iNTS infections, spanning 28 years (1994-2021) in China. Phylogenomic analysis was performed to explore the population structure and evolutionary history of the Chinese isolates within a global context. Molecular investigations of AMR genes, virulence factors, mobile genetic elements and pan-genomes were also performed.

RESULTS

The Chinese S. Enteritidis collections exhibited a high level of multidrug resistance (MDR), including high resistance to nalidixic acid (97.67%). Notably, the multidrug resistance rate of iNTS strains has significantly increased over the past decade. Phylogenomic analysis showed that the majority of the Chinese isolates (98.63%) were distributed in the global pandemic lineage L1, while the other lineages were highly continent-specific. Particularly, the Chinese isolates were predominantly distributed in sublineages L1.2 (37.45%) and L1.3 (59.26%), forming two main Chinese clades (MCC1&2). The most recent common ancestor of MCC1&2 dated back to 1944 and 2004, respectively. The lineage L1, especially MCC1&2, harbored the most amount of AMR determinants and virulence genes, which was mainly due to the presence of a hybrid virulence-resistance plasmid and coexistence of different types of AMR plasmids in S. Enteritidis.

CONCLUSIONS

S. Enteritidis has evolved unique clonal clusters, MCC1&2, with critical MDR in China, which phylogenetically constitute an extension of the globally epidemic lineage and were characterized by distinct genetic traits. These clades have induced extensive outbreaks of gastroenteritis and serious cases of iNTS in China, underscoring the pressing nature and severity of this public health crisis. Implementing the One-Health strategy, longstanding routine surveillance and further genomic epidemiological studies are urgently required to capture epidemics, monitor changes in bacterial populations and determine the consequent risk to global public health.

Genomic survey of multidrug resistant Salmonella enterica serovar Minnesota clones in chicken products

Salmonella enterica serovar Minnesota (S. Minnesota) is an emerging serovar that persists within poultry supply chains, potentially causing outbreaks in humans. Understanding its population genomics is crucial for designing preventive measures. We performed a genomic surveillance study of S. Minnesota by analyzing 259 isolates from poultry in Saudi Arabia. Whole-genome sequencing data for these isolates were analyzed to characterize emerging clones and the genetic factors underlying antimicrobial resistance and virulence. We compared the isolates to all available global genomes of S. Minnesota. Our results revealed the emergence of four clones, three of which were mixed with global strains. These clones exhibited higher levels of antimicrobial resistance and virulence due to the acquisition of multiple plasmids, particularly IncC plasmids, carrying resistance and virulence genes. IncC plasmids underwent genomic rearrangements, presenting diverse configurations of resistance genes. Our findings demonstrate the emergence and persistence of pathogenic and multidrug-resistant S. Minnesota clones.

The evolutionary history of the common bean (Phaseolus vulgaris) revealed by chloroplast and nuclear genomes analysis

KEY MESSAGE

The origin of common bean was investigated throughout chloroplast and nuclear WGS data considering recombination events. Our results support the Mesoamerican origin of common bean. The remarkable evolutionary history of the common bean (Phaseolus vulgaris L.) has led to the emergence of three wild main gene pools corresponding to three different eco-geographical areas: Mesoamerica, the Andes and northern Peru/Ecuador. Recent works proposed novel scenarios, and the northern Peru/Ecuador population has been described as a new species called P. debouckii, rekindling the debate about the origin of P. vulgaris. Here we shed light on the origin of P. vulgaris by analyzing the chloroplast and nuclear genomes of a large varietal collection representing the entire geographical distribution of wild forms including a large collection of Mesoamerican and Andean individuals. We assembled 37 chloroplast genomes de novo and used them to construct a time frame for the divergence of the genotypes under investigation, revealing that the separation of the Mesoamerican and northern Peru/Ecuador gene pools occurred ~ 0.15 Mya. Our results clearly support a Mesoamerican origin of the common bean and reject the recent P. deboukii hypothesis. These results also imply two independent migratory events from Mesoamerica to the North and South Andes, probably facilitated by birds. Our work represents a paradigmatic example of the importance of taking into account the genetic rearrangements produced by recombination when investigating phylogeny and of the analysis of wild forms when studying the evolutionary history of a crop species.

Genomic diversity and evolutionary patterns of Edwardsiella ictaluri affecting farmed striped catfish (Pangasianodon hypophthalmus) in Vietnam over 20 years

Edwardsiella ictaluri continues to pose a significant risk to the health and production of striped catfish (Pangasianodon hypophthalmus) in Vietnam. Whilst recent advances in genomic sequencing provide an insight into the global genomic diversity of this important fish pathogen, genome-wide analysis of Vietnamese isolates recovered over time is lacking. In this study, we used a whole-genome sequencing approach to compare the genomes of 31 E. ictaluri isolates recovered over a 20-year period (2001-2021) and performed comparative genomic analysis to explore temporal changes in genome diversity, population structure and mechanisms driving pathogenesis and antimicrobial resistance. Our findings revealed an open pan-genome with 4148 genes and a core genome (3 060 genes) accounting for over two-thirds of the genome. Moreover, we found the genomes sequenced to classify into two distinct lineages and estimated the ancestral origin of these lineages within Vietnam to date back to the 1950s. Plasmids were highly prevalent in Vietnamese E. ictaluri, with isolates harbouring up to four plasmids within their genome. Further, a diverse mobilome was observed with nine different plasmid types detected across the genome collection. Exploration of putative plasmids revealed a diverse set of antimicrobial resistance genes (ARGs) against key antibiotics used in Vietnamese aquaculture and virulence genes associated with protein secretion systems. Correlation analysis revealed the total number of ARGs detected in genomes to increase with isolate recovery time. Whilst the number of virulence genes remained relatively stable, temporal variation was noted in several virulence factors related to motility and immune system modulation. Findings from this study highlight the need for continued genomic surveillance to monitor changes in antimicrobial resistance and pathogenesis, to help inform the development of disease control and management strategies.

Genomic Epidemiology of Pseudomonas aeruginosa Sequence Type 111

PURPOSE

Pseudomonas aeruginosa ST111 is a global multidrug resistant (MDR) high-risk clone and comprehensive data about its molecular epidemiology is limited in Canada. Comprehensive data about the evolution of ST111 clades is limited. We characterized a Canadian collection of ST111 causing bloodstream infections and investigated the genomic relationship between Canadian and global ST111.

MATERIAL AND METHODS

We used long and short read WGS to characterize Canadian ST111 (n = 10 from 2010-18). We performed phylogenetic analysis on a global collection of ST111 (n = 969) and investigated the evolutionary history of clades using BEAST.

RESULTS

ST111 belonged to 3 clades (A, B, C) and two subclades (C1, C2). ST111-A was the ancestral clade while clades B, C1 and C2 emerged during the 1700s and 1800s. ST111-C2 dominated the global ST111 population. Serotype switching from O4 to O12 and the acquisition of Tn21, gyrA_T83I, parC_S87L, In59 with blaVIM-2 and aacA29 over time, were important in the evolution of ST111-C2. The Calgary ST111 strains consisted of a diverse collection that belonged to ST111-A (O4), ST111-C1 (O4) and ST111-C2 (O12) with different transposon structures.

CONCLUSIONS

We provided details on the emergence and evolution of different ST111 clades over time and highlighted the roles of serotype switching and the acquisition of certain AMR determinants and transposon structures in the evolution of ST111-C2.

Global epidemiology of Mycobacterium tuberculosis lineage 4 insights from Ecuadorian genomic data

Tuberculosis is a global public health concern, and understanding Mycobacterium tuberculosis transmission routes and genetic diversity of M. tuberculosis is crucial for outbreak control. This study aimed to explore the genomic epidemiology and genetic diversity of M. tuberculosis in Ecuador by analyzing 88 local isolates and 415 public genomes from 19 countries within the Euro-American lineage (L4). Our results revealed significant genomic diversity among the isolates, particularly in the genes related to protein processing, carbohydrate metabolism, lipid metabolism, and xenobiotic biodegradation and metabolism. The population structure analysis showed that sub-lineages 4.3.2/3 (35.4%), 4.1.2.1 (22.7%), 4.4.1 (12.7%), and 4.1.1. (10.7%) were the most prevalent. Phylogenetic and transmission network analyses suggest that these isolates circulating within Ecuador share genetic ties with isolates from other continents, implying historical and ongoing intercontinental transmission events. Our findings underscore the importance of integrating genomic data into public health strategies for tuberculosis control and suggest that enhanced genomic surveillance is essential for understanding and mitigating the global spread of M. tuberculosis. This study provides a comprehensive genomic framework for future epidemiological investigations and control measures targeting M. tuberculosis L4 in Ecuador.

Global Invasion History and Genomic Signatures of Adaptation of the Highly Invasive Sycamore Lace Bug

Invasive species cause massive economic and ecological damages. Climate change has resulted in an unprecedented increase in the number and impact of invasive species; however, the mechanisms underlying these invasions are unclear. The sycamore lace bug, Corythucha ciliata, is a highly invasive species originating from North America and has expanded across the Northern Hemisphere since the 1960s. In this study, we assembled the C. ciliata genome using high-coverage Pacific Biosciences (PacBio), Illumina, and high-throughput chromosome conformation capture (Hi-C) sequencing. A total of 15,278 protein-coding genes were identified, and expansions of gene families with oxidoreductase and metabolic activities were observed. In-depth resequencing of 402 samples from native and nine invaded countries across three continents revealed 2.74 million single nucleotide polymorphisms. Two major invasion routes of C. ciliata were identified from North America to Europe and Japan, with a contact zone forming in East Asia. Genomic signatures of selection associated with invasion and long-term balancing selection in native ranges were identified. These genomic signatures overlapped with each other as well as with expanded genes, suggesting improvements in the oxidative stress and thermal tolerance of C. ciliata. These findings offer valuable insights into the genomic architecture and adaptive evolution underlying the invasive capabilities of species during rapid environmental changes.

Phylogeography of the Mesa Silverside fish Chirostoma jordani (Woolman, 1894) throughout the Mexican Plateau

Background

Understanding the processes that influence distribution of organisms is a major goal in evolutionary biology. Speciation in freshwater fishes is mainly associated with the "island-like" model of evolution, in which the formation of land barriers between different hydrographic basins interrupts gene flow and promotes isolation. Freshwater fish therefore provide an excellent model system for macro- and micro-evolutionary studies. The Mesa Silverside, Chirostoma jordani, is one of the most widespread freshwater fish species in the Mexican Plateau, a geologically complex physiographic region with a long history of genesis, destruction and compartmentalization of hydrographic basins that has promoted the dispersal and isolation of freshwater fishes.

Methods

We used mitochondrial (Cytb and D-loop) and nuclear (first intron of the ribosomal protein S7) data and used phylogeographic and coalescent based methods to elucidate the evolutionary history of C. jordani throughout its distributional range on the Mexican Plateau.

Results

The results obtained in the present study revealed that C. jordani consists of two main genetic groups with geographical correspondence. Clade I occur exclusively in north-western basin and shows population structure. Clade II is widely distributed across the west, central and eastern basins without population structure. The split between these two main clades was estimated at 1.4 Mya. This cladogenetic event may be associated with the allopatric process promoted by the fragmentation and compartmentalization of hydrographic basins induced by the geological and climatic history of the Mexican Plateau.

Comprehensive lung microbial gene and genome catalogs assist the mechanism survey of Mesomycoplasma hyopneumoniae strains causing pig lung lesions

Understanding the community structure of the lower respiratory tract microbiome is crucial for elucidating its roles in respiratory tract diseases. However, there are few studies about this topic due to the difficulty in obtaining microbial samples from both healthy and disease individuals. Here, using 744 high-depth metagenomic sequencing data of lower respiratory tract microbial samples from 675 well-phenotyped pigs, we constructed a lung microbial gene catalog containing the largest scale of 10,031,593 nonredundant genes to date, 44.8% of which are novel. We obtained 356 metagenome-assembled genomes (MAGs) which were further clustered into 256 species-level genome bins with 41.8% being first reported in the current databases. Based on these data sets and through integrated analysis of the isolation of the related bacterial strains, in vitro infection, and RNA sequencing, we identified and confirmed that Mesomycoplasma hyopneumoniae (M. hyopneumoniae) MAG_47 and its adhesion-related virulence factors (VFs) were associated with lung lesions in pigs. Differential expression levels of adhesion- and immunomodulation-related VFs likely determined the heterogenicity of adhesion and pathogenicity among M. hyopneumoniae strains. M. hyopneumoniae adhesion activated several pathways, including nuclear factor kappa-light-chain-enhancer of activated B, mitogen-activated protein kinase, cell apoptosis, T helper 1 and T helper 2 cell differentiation, tumor necrosis factor signaling, interleukin-6/janus kinase 2/signal transducer and activator of transcription signaling, and response to reactive oxygen species, leading to cilium loss, epithelial cell‒cell barrier disruption, and lung tissue lesions. Finally, we observed the similar phylogenetic compositions of the lung microbiome between humans with Mycoplasma pneumoniae and pigs infected with M. hyopneumoniae. The results provided important insights into pig lower respiratory tract microbiome and its relationship with lung health.

Genomic epidemiology and phenotypic characterisation of Salmonella enterica serovar Panama in Victoria, Australia

Salmonella enterica serovar Panama, a causative agent of non-typhoidal salmonellosis (NTS), is one of several serovars that causes invasive NTS disease (iNTS) in humans. S. Panama is an understudied pathogen, with its pathobiology poorly understood. It is a predominant iNTS serovar in Australia, a high-income country with high rates of salmonellosis, where S. Panama has been documented to have a high odds ratio (13.9-15.26) for causing iNTS. This study investigates the genomic epidemiology and antimicrobial resistance profiles of all S. Panama isolates recovered in Victoria, Australia, between 2000 and 2021. We examined the infection dynamics of S. Panama in seven isolates, representing the genetic diversity of the study population. Two sub-lineages, encompassed within a previously described Asian lineage, were identified. Multi-drug resistance (resistance to ≥3 drug classes) was detected in 46 (51.7%) Australian isolates. The plasmid-mediated colistin resistance gene, mcr1.1, was detected in one Australian S. Panama isolate, carried by an IncI plasmid previously reported in Salmonella and Escherichia coli isolates collected from poultry in South-East Asia. Examination of the intracellular replication dynamics of S. Panama isolates demonstrated diverse phenotypes. In THP-1 derived macrophages, despite low host cell uptake, S. Panama showed higher replication rates over time compared to S. enterica serovar Typhimurium. However, a causative genotype could not be identified to explain this observed phenotype. This study provides insights into the S. Panama isolates circulating in Australia over two-decades, finding that 78% were linked to international travel suggesting importation in Australia. It shows MDR was common in this iNTS serovar, and colistin resistance reported for the first time. It provides the first data on the host-pathogen interactions of S. Panama in Australia, which will aid our collective understanding of the pathobiology of S. Panama and iNTS serovars more broadly.

In-host intra- and inter-species transfer of blaKPC-2 and blaNDM-1 in Serratia marcescens and its local and global epidemiology

OBJECTIVES

The aim of this study was to investigate interspecies transfer of resistance gene blaNDM-1 and intraspecies transfer of resistance gene blaKPC-2 in Serratia marcescens, and explore the epidemical and evolutionary characteristics of carbapenemase-producing S. marcescens (CPSM) regionally and globally.

METHODS

Interspecies and intraspecies transfer of blaKPC-2- or blaNDM-1 were identified by antimicrobial susceptibility testing, plasmid conjugation and curing, discovery of transposable units (TUs), outer membrane vesicles (OMVs), qPCR, whole-genome sequencing (WGS) and bioinformatic analysis. The genomic evolution of CPSM strains was explored by cgSNP and maximum-likelihood phylogenetic tree.

RESULTS

CPSM S50079 strain, co-carrying blaKPC-2 and blaNDM-1 on one plasmid, was isolated from the blood of a patient with acute pancreatitis and could generate TUs carrying either blaKPC-2 or blaNDM-1. The interspecies transfer of blaNDM-1-carrying plasmid from Providencia rettgeri P50213, producing the identical blaNDM-1-carrying TUs, to S. marcescens S50079K, an S50079 variant via plasmid curing, was identified through blaNDM-1-harbouring plasmid conjugation and OMVs transfer. Moreover, the intraspecies transfer of blaKPC-2, mediated by IS26 from plasmid to chromosome in S50079, was also identified. In another patient, who underwent lung transplantation, interspecies transfer of blaNDM-1 carried by IncX3 plasmid was identified among S. marcescens and Citrobacter freundii as well as Enterobacter hormaechei via plasmid transfer. Furthermore, 11 CPSM from 349 non-repetitive S. marcescens strains were identified in the same hospital, and clonal dissemination, with carbapenemase evolution from blaKPC-2 to both blaKPC-2 and blaNDM-1, was found in the 8 CPSM across 4 years. Finally, the analysis of 236 global CPSM from 835 non-repetitive S. marcescens genomes, retrieved from the NCBI database, revealed long-term spread and evolution worldwide, and would cause the convergence of more carbapenemase genes.

CONCLUSIONS

Interspecies transfer of resistance gene blaNDM-1 and intraspecies transfer of resistance gene blaKPC-2 in CPSM were identified. Nosocomial and global dissemination of CPSM were revealed and more urgent surveillance was acquired.

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.

Genomic insight into the diversity of Glaesserella parasuis isolates from 19 countries

Glaesserella parasuis is a commensal bacterial organism found in the upper respiratory tract of healthy pigs and the etiological agent of Glässer's disease, which causes severe economic losses in the swine industry. This study aimed to better understand the epidemiological characteristics of this opportunistic pathogen. We investigated the prevalence and distribution of sequence types (STs), serovars, antimicrobial resistance genes (ARGs), and potential virulence factors (VFs) in 764 G. parasuis isolates collected from diseased and healthy pigs from 19 countries, including China. Multilocus sequence typing showed a high degree of variation with 334 STs, of which 93 were not previously recognized. Phylogenetic analysis revealed two major clades distinguished by isolation year, source, country, and serovar. The dominant serovars of G. parasuis were serovars 4 (19.50%), 7 (15.97%), 5/12 (13.87%), and 13 (12.30%). Serovar 7 gradually became one of the dominant serovars in G. parasuis with more VFs and fewer ARGs. Serovars 4 and 5/12 were the most frequent serovars in diseased pigs, whereas serovars 2, 8, and 11 were predominant in healthy pigs. Serovars 7 and 13 possessed more VFs than the other serovars. This study provides novel insights into the global prevalence and epidemiology of G. parasuis and valuable clues for further investigation into the pathogenicity of G. parasuis, which will facilitate the development of effective vaccines.IMPORTANCEGlaesserella parasuis is a clinically important gram-negative opportunistic pathogen, which causes serious financial losses in swine industry on a global scale. No vaccine is known that provides cross-protection against all 15 serovars; furthermore, the correlation between serovar and virulence is largely unknown. This study provides a large number of sequenced strains in 19 countries and compares the genomic diversity of G. parasuis between diseased and healthy pigs. We found a slight change in the dominant serovar of G. parasuis in the world, with serovar 7 gradually emerging as one of the predominant serovars. The observed higher average number of VFs in this particular serovar strain challenges the previously held notion that serovar 7 is non-virulent, indicating a more complex virulence landscape than previously understood. Our analysis indicating that six ARGs [tet(B), sul2, aph(3')-Ia, aph (6)-Id, blaROB-1, and aph(3'')-Ib] are likely to be transmitted horizontally in their entirety. By analyzing VFs, we provided an improved understanding of the virulence of G. parasuis, and these key findings suggest that vaccine development will be challenging.

Synergistic effects of Pleistocene geological and climatic events on complex phylogeographic history of widespread sympatric species of Megaloptera in East Asia

Unraveling the phylogeographic histories of species remains a key endeavor for comprehending the evolutionary processes contributing to the rich biodiversity and high endemism found in East Asia. In this study, we explored the phylogeographic patterns and demographic histories of three endemic fishfly and dobsonfly species ( Neochauliodes formosanus, Protohermes costalis, and Neoneuromus orientalis) belonging to the holometabolan order Megaloptera. These species, which share a broad and largely overlapping distribution, were analyzed using comprehensive mitogenomic data. Our findings revealed a consistent influence of vicariance on the population isolation of Neoc. formosanus and P. costalis between Hainan, Taiwan, and the East Asian mainland during the early Pleistocene, potentially hindering subsequent colonization of the later diverged Neon. orientalis to these islands. Additionally, we unveiled the dual function of the major mountain ranges in East Asia, serving both as barriers and conduits, in shaping the population structure of all three species. Notably, we demonstrated that these co-distributed species originated from Southwest, Southern, and eastern Central China, respectively, then subsequently migrated along multi-directional routes, leading to their sympatric distribution on the East Asian mainland. Furthermore, our results highlighted the significance of Pleistocene land bridges along the eastern coast of East Asia in facilitating the dispersal of mountain-dwelling insects with low dispersal ability. Overall, this study provides novel insight into the synergistic impact of Pleistocene geological and climatic events in shaping the diversity and distribution of aquatic insects in East Asia.

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.

Phenotypic and Genetic Diversity of Xanthomonads Isolated from Pepper (Capsicum spp.) in Taiwan from 1989 to 2019

Bacterial spot caused by Xanthomonas spp. is an economically important disease of pepper causing significant yield losses in Taiwan. Monitoring the pathogen population on a continuous basis is necessary for developing disease management strategies. We analyzed a collection of xanthomonad strains isolated from pepper in Taiwan between 1989 and 2019. Among the sequenced genomes, 65 were identified as Xanthomonas euvesicatoria, and 10 were X. perforans. Thirty-five X. euvesicatoria and 10 X. perforans strains were copper tolerant, whereas only four X. euvesicatoria and none of the X. perforans strains were tolerant to streptomycin. Nine X. euvesicatoria strains were amylolytic, which is considered an unusual characteristic for X. euvesicatoria. Bayesian analysis of the population structure based on core gene single-nucleotide polymorphisms clustered the strains into five clusters for X. euvesicatoria and three clusters for X. perforans. One X. perforans cluster, designated as TP-2019, appears to be a novel genetic cluster based on core genes, accessory gene content, and effector profile. This knowledge of pathogen diversity with whole genomic information will be useful in future comparative studies and in improving breeding programs to develop disease-resistant cultivars and other disease management options.

Can genomics and meteorology predict outbreaks of legionellosis in urban settings?

Legionella pneumophila is ubiquitous and sporadically infects humans causing Legionnaire's disease (LD). Globally, reported cases of LD have risen fourfold from 2000 to 2014. In 2016, Sydney, Australia was the epicenter of an outbreak caused by L. pneumophila serogroup 1 (Lpsg1). Whole-genome sequencing was instrumental in identifying the causal clone which was found in multiple locations across the city. This study examined the epidemiology of Lpsg1 in an urban environment, assessed typing schemes to classify resident clones, and investigated the association between local climate variables and LD outbreaks. Of 223 local Lpsg1 isolates, we identified dominant clones with one clone isolated from patients in high frequency during outbreak investigations. The core genome multi-locus sequence typing scheme was the most reliable in identifying this Lpsg1 clone. While an increase in humidity and rainfall was found to coincide with a rise in LD cases, the incidence of the major L. pneumophila outbreak clone did not link to weather phenomena. These findings demonstrated the role of high-resolution typing and weather context assessment in determining source attribution for LD outbreaks in urban settings, particularly when clinical isolates remain scarce.IMPORTANCEWe investigated the genomic and meteorological influences of infections caused by Legionella pneumophila in Sydney, Australia. Our study contributes to a knowledge gap of factors that drive outbreaks of legionellosis compared to sporadic infections in urban settings. In such cases, clinical isolates can be rare, and thus, other data are needed to inform decision-making around control measures. The study revealed that core genome multi-locus sequence typing is a reliable and adaptable technique when investigating Lpsg1 outbreaks. In Sydney, the genomic profile of Lpsg1 was dominated by a single clone, which was linked to numerous community cases over a period of 40 years. Interestingly, the peak in legionellosis cases during Autumn was not associated with this prevalent outbreak clone. Incorporating meteorological data with Lpsg1 genomics can support risk assessment strategies for legionellosis in urban environments, and this approach may be relevant for other densely populated regions globally.

Clonal spread of blaNDM-1-carrying Salmonella enterica serovar Typhimurium clone ST34 and wide spread of IncHI2/ST3-blaNDM-5 plasmid in China

OBJECTIVES

To characterize blaNDM-carrying Salmonella recovered from a pig slaughterhouse.

METHODS

In this study, 46 environment samples were collected from a slaughterhouse in China, and screened for carbapenem-resistant Enterobacterales. WGS, antimicrobial susceptibility testing and conjugation experiments were carried out to identify the isolates' resistance phenotypes and genetic characteristics. The phylogenetic relatedness of the Salmonella isolates obtained in this study and Salmonella (ST34 and ST29) in GenBank was determined.

RESULTS

Two ST34 Salmonella Typhimurium and one ST29 Salmonella Stanley, recovered from three environmental samples (6.52%), were positive for blaNDM-1 and blaNDM-5, respectively. The two ST34 S. Typhimurium strains exhibited a close relationship (10-36 SNPs) with two human-derived blaNDM-1-bearing isolates from China (Hong Kong and Guangxi Province) and two blaNDM-negative ST34 Salmonella strains from the UK. The blaNDM-1 genes were located on IncHI2/ST3 plasmids. The capture of blaNDM-1 by the IncHI2/ST3 plasmid seems to be due to homologous recombination mediated by circular structures, as the genetic arrangements of the blaNDM-1 gene contain two IS26 elements of the same orientation. The blaNDM-5 gene was also carried by the IncHI2/ST3 plasmid, which shares highly similar structures with other blaNDM-5-bearing IncHI2/ST3 plasmids from other sources (fish, chicken, duck, human).

CONCLUSIONS

This is the first report of a blaNDM-5-carrying IncHI2/ST3 plasmid in Salmonella. The clonal spread of NDM-1-producing ST34 S. Typhimurium across human and animal-associated environments, and the widespread dissemination of epidemic blaNDM-5-carrying IncHI2/ST3 plasmids among Enterobacteriaceae in China indicate the potential of further dissemination of blaNDM among Salmonella, which poses a threat to public health.

Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group

Oceanic islands have been recognized as natural laboratories in which to study a great variety of evolutionary processes. One such process is evolutionary radiations, the diversification of a single ancestor into a number of species that inhabit different environments and differ in the traits that allow them to exploit those environments. The factors that drive evolutionary radiations have been studied for decades in charismatic organisms such as birds or lizards, but are lacking in lichen-forming fungi, despite recent reports of some lineages showing diversification patterns congruent with radiation. Here we propose the Ramalina decipiens group as a model system in which to carry out such studies. This group is currently thought to be comprised of five saxicolous species, all of them endemic to the Macaronesian region (the Azores, Madeira, Selvagens, Canary and Cape Verde islands). Three species are single-island endemics (a rare geographic distribution pattern in lichens), whereas two are widespread and show extreme morphological variation. The latter are suspected to harbor unrecognized species-level lineages. In order to use the Ramalina decipiens group as a model system it is necessary to resolve the group's phylogeny and to clarify its species boundaries. In this study we attempt to do so following an integrative taxonomy approach. We constructed a phylogenetic tree based on six molecular markers, four of which are newly developed and generated competing species hypotheses based on molecular (species discovery strategies based on both single locus and multilocus datasets) and phenotypic data (unsupervised clustering algorithms based on morphology, secondary chemistry and geographic origin). We found that taxonomic diversity in the Ramalina decipiens group has been highly underestimated in previous studies. In consequence, we describe six new species, most of them single-island endemics and provide a key to the group. Phylogenetic relationships among species have been reconstructed with almost full support which, coupled with the endemic character of the group, makes it an excellent system for the study of island radiations in lichen-forming fungi. Citation: Blázquez M, Pérez-Vargas I, Garrido-Benavent I, et al. 2024. Endless forms most frustrating: disentangling species boundaries in the Ramalina decipiens group (Lecanoromycetes, Ascomycota), with the description of six new species and a key to the group. Persoonia 52: 44-93. https://doi.org/10.3767/persoonia.2024.52.03 .

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.

Characterization of a Salmonella enterica serovar Typhimurium lineage with rough colony morphology and multidrug resistance

Salmonella enterica serovar Typhimurium (S. Typhimurium) is a major cause of salmonellosis, and the emergence of multidrug-resistant pathovariants has become a growing concern. Here, we investigate a distinct rough colony variant exhibiting a strong biofilm-forming ability isolated in China. Whole-genome sequencing on 2,212 Chinese isolates and 1,739 publicly available genomes reveals the population structure and evolutionary history of the rough colony variants. Characterized by macro, red, dry, and rough (mrdar) colonies, these variants demonstrate enhanced biofilm formation at 28 °C and 37 °C compared to typical rdar colonies. The mrdar variants exhibit extensive multidrug resistance, with significantly higher resistance to at least five classes of antimicrobial agents compared to non-mrdar variants. This resistance is primarily conferred by an IncHI2 plasmid harboring 19 antimicrobial resistance genes. Phylogenomic analysis divides the global collections into six lineages. The majority of mrdar variants belong to sublineage L6.5, which originated from Chinese smooth colony strains and possibly emerged circa 1977. Among the mrdar variants, upregulation of the csgDEFG operons is observed, probably due to a distinct point mutation (-44G > T) in the csgD gene promoter. Pangenome and genome-wide association analyses identify 87 specific accessory genes and 72 distinct single nucleotide polymorphisms associated with the mrdar morphotype.

Rapid identification and subsequent contextualization of an outbreak of methicillin-resistant Staphylococcus aureus in a neonatal intensive care unit using nanopore sequencing

Outbreaks of methicillin-resistant Staphylococcus aureus (MRSA) are well described in the neonatal intensive care unit (NICU) setting. Genomics has revolutionized the investigation of such outbreaks; however, to date, this has largely been completed retrospectively and has typically relied on short-read platforms. In 2022, our laboratory established a prospective genomic surveillance system using Oxford Nanopore Technologies sequencing for rapid outbreak detection. Herein, using this system, we describe the detection and control of an outbreak of sequence-type (ST)97 MRSA in our NICU. The outbreak was identified 13 days after the first MRSA-positive culture and at a point where there were only two known cases. Ward screening rapidly defined the extent of the outbreak, with six other infants found to be colonized. There was minimal transmission once the outbreak had been detected and appropriate infection control measures had been instituted; only two further ST97 cases were detected, along with three unrelated non-ST97 MRSA cases. To contextualize the outbreak, core-genome single-nucleotide variants were identified for phylogenetic analysis after de novo assembly of nanopore data. Comparisons with global (n=45) and national surveillance (n=35) ST97 genomes revealed the stepwise evolution of methicillin resistance within this ST97 subset. A distinct cluster comprising nine of the ten ST97-IVa genomes from the NICU was identified, with strains from 2020 to 2022 national surveillance serving as outgroups to this cluster. One ST97-IVa genome presumed to be part of the outbreak formed an outgroup and was retrospectively excluded. A second phylogeny was created using Illumina sequencing, which considerably reduced the branch lengths of the NICU isolates on the phylogenetic tree. However, the overall tree topology and conclusions were unchanged, with the exception of the NICU outbreak cluster, where differences in branch lengths were observed. This analysis demonstrated the ability of a nanopore-only prospective genomic surveillance system to rapidly identify and contextualize an outbreak of MRSA in a NICU.

Epidemiological and genetic characterization of multidrug-resistant non-O1 and non-O139 Vibrio cholerae from food in southern China

This study reports a comprehensive epidemiological and genetic analysis of V. cholerae strains, specifically non-O1/non-O139 serogroups, isolated from animal-derived food samples in Guangdong province from 2015 to 2019. A total of 21 V. cholerae strains were obtained, which exhibited high resistance rates for nalidixic acid (57.14 %, 12/21), ampicillin (33.33 %, 7/21), and ciprofloxacin (19.05 %, 4/21). The quinolone resistance-related gene, qnrVC, was prevalent in 80.95 % (17/21) of the isolates. Additionally, chromosomally mediated quinolone-resistance mutations, including mutations in GyrA at position 83 (S83I) and ParC at position 85 (S85L), were detected in 47.62 % of the isolates. The combination of target mutation and qnrVC genes was shown to mediate resistance or intermediate resistance to ciprofloxacin in V. cholerae. Furthermore, an IncC-type conjugative plasmid carrying thirteen antibiotic resistance genes, including genes conferring resistance to two clinically important antibiotics, cephalosporins and fluoroquinolones, was identified in the shrimp-derived strain Vc516. While none of our food isolates harbored the toxigenic CTX- and TCP-encoding genes, they did possess genes encoding toxins such as HlyA and Autoinducer-2. Notably, some V. cholerae strains from this study exhibited a close genetic relationship with clinical strains, suggesting their potential to cause human infections. Taken together, this study provides a comprehensive view of the epidemiological features and genetic basis of antimicrobial resistance and virulence potential of V. cholerae strains isolated from food in southern China, thereby advancing our understanding of this important pathogen.

Genetic analysis of harvest samples reveals population structure in a highly mobile generalist carnivore

Delineating wildlife population boundaries is important for effective population monitoring and management. The bobcat (Lynx rufus) is a highly mobile generalist carnivore that is ecologically and economically important. We sampled 1225 bobcats harvested in South Dakota, USA (2014-2019), of which 878 were retained to assess genetic diversity and infer population genetic structure using 17 microsatellite loci. We assigned individuals to genetic clusters (K) using spatial and nonspatial Bayesian clustering algorithms and quantified differentiation (F ST and G ST ″ ) among clusters. We found support for population genetic structure at K = 2 and K = 4, with pairwise F ST and G ST ″ values indicating weak to moderate differentiation, respectively, among clusters. For K = 2, eastern and western clusters aligned closely with historical bobcat management units and were consistent with a longitudinal suture zone for bobcats previously identified in the Great Plains. We did not observe patterns of population genetic structure aligning with major rivers or highways. Genetic divergence observed at K = 4 aligned roughly with ecoregion breaks and may be associated with environmental gradients, but additional sampling with more precise locational data may be necessary to validate these patterns. Our findings reveal that cryptic population structure may occur in highly mobile and broadly distributed generalist carnivores, highlighting the importance of considering population structure when establishing population monitoring programs or harvest regulations. Our study further demonstrates that for elusive furbearers, harvest can provide an efficient, broad-scale sampling approach for genetic population assessments.

Phylogeographic Pattern of the Assassin Bug Sycanus bifidus Inferred from Mitochondrial Genomes and Nuclear Genes

The assassin bug Sycanus bifidus has a wide distribution across southern China. This study explored its distribution and evolution by analyzing mitochondrial and nuclear ribosomal RNA genes, revealing how Pleistocene climate and geological changes shaped its phylogeography. We identified two main clades, A and B, that diverged in the Middle Pleistocene. Hainan Island's populations form a unique group within Clade A, suggesting that the Qiongzhou Strait served as a dispersal corridor during glaciation. Rising sea levels likely separated the Hainan population afterward. Ecological niche modeling showed that both populations have been viable since the last interglacial period, with demographic analyses indicating possible expansions during the Middle and Late Pleistocene, driven by favorable climates. This study highlights the significant effects of Pleistocene sea-level and climatic changes on the distribution and evolution of S. bifidus in China.

The emergence of highly resistant and hypervirulent Klebsiella pneumoniae CC14 clone in a tertiary hospital over 8 years

BACKGROUND

Klebsiella pneumoniae is a major bacterial and opportunistic human pathogen, increasingly recognized as a healthcare burden globally. The convergence of resistance and virulence in K. pneumoniae strains has led to the formation of hypervirulent and multidrug-resistant strains with dual risk, limiting treatment options. K. pneumoniae clones are known to emerge locally and spread globally. Therefore, an understanding of the dynamics and evolution of the emerging strains in hospitals is warranted to prevent future outbreaks.

METHODS

In this study, we conducted an in-depth genomic analysis on a large-scale collection of 328 multidrug-resistant (MDR) K. pneumoniae strains recovered from 239 patients from a single major hospital in the western coastal city of Jeddah in Saudi Arabia from 2014 through 2022. We employed a broad range of phylogenetic and phylodynamic methods to understand the evolution of the predominant clones on epidemiological time scales, virulence and resistance determinants, and their dynamics. We also integrated the genomic data with detailed electronic health record (EHR) data for the patients to understand the clinical implications of the resistance and virulence of different strains.

RESULTS

We discovered a diverse population underlying the infections, with most strains belonging to Clonal Complex 14 (CC14) exhibiting dominance. Specifically, we observed the emergence and continuous expansion of strains belonging to the dominant ST2096 in the CC14 clade across hospital wards in recent years. These strains acquired resistance mutations against colistin and extended spectrum β-lactamase (ESBL) and carbapenemase genes, namely blaOXA-48 and blaOXA-232, located on three distinct plasmids, on epidemiological time scales. Strains of ST2096 exhibited a high virulence level with the presence of the siderophore aerobactin (iuc) locus situated on the same mosaic plasmid as the ESBL gene. Integration of ST2096 with EHR data confirmed the significant link between colonization by ST2096 and the diagnosis of sepsis and elevated in-hospital mortality (p-value < 0.05).

CONCLUSIONS

Overall, these results demonstrate the clinical significance of ST2096 clones and illustrate the rapid evolution of an emerging hypervirulent and MDR K. pneumoniae in a clinical setting.

Unravelling the complex interplay between antibiotic consumption and adaptive changes in methicillin-resistant Staphylococcus aureus

OBJECTIVES

This study aims to elucidate the genomic dynamics driving the emergence of antimicrobial resistance (AMR), with a specific focus on the interplay between AMR and antimicrobial usage.

METHODS

We conducted a comprehensive analysis using a ST239 methicillin-resistant Staphylococcus aureus (MRSA) dataset over a continuous 12-year period from a single hospital. Genomic analyses were performed tracking the changes in MRSA populations, particularly the emergence of reduced vancomycin susceptibility, and assessing the impact of glycopeptide use on these emergence events.

RESULTS

Our findings reveal a significant correlation between hospital glycopeptide usage and the selection of MRSA strains with reduced vancomycin susceptibility. Genomic analyses provided insights into the molecular mechanisms driving resistance emergence, including the slowing of the molecular clock rate in response to heightened antimicrobial consumption.

CONCLUSIONS

In conclusion, this study the highlights the complex dynamics between AMR and antimicrobial use at the hospital level. The observed correlation between antimicrobial consumption and the development of less susceptible MRSA strains underscores the importance of antimicrobial stewardship programmes and the establishment of optimal consumption thresholds for mitigating AMR effectively.

Recovery and genetic characterization of clinically-relevant ST2 carbapenem-resistant Acinetobacter baumannii isolates from untreated hospital sewage in Zhejiang Province, China

The global transmission of carbapenem-resistant Acinetobacter baumannii (CRAB) poses a significant and grave threat to human health. To investigate the potential relationship between hospital sewage and the transmission of CRAB within healthcare facilities, isolates of Acinetobacter spp. obtained from untreated hospital sewage samples were subjected to antimicrobial susceptibility tests, genome sequencing, and bioinformatic and phylogenetic tree analysis, and that data were matched with those of the clinical isolates. Among the 70 Acinetobacter spp. sewage isolates tested, A. baumannii was the most prevalent and detectable in 5 hospitals, followed by A. nosocomialis and A. gerneri. Worryingly, 57.14 % (40/70) of the isolates were MDR, with 25.71 % (18/70) being resistant to carbapenem. When utilizing the Pasteur scheme, ST2 was the predominant type among these CRAB isolates, with Tn2006 (ΔISAba1-blaOXA-23-ATPase-yeeB-yeeA-ΔISAba1) and Tn2009 (ΔISAba1-blaOXA-23-ATPase-hp-parA-yeeC-hp-yeeB-ΔISAba1) being the key mobile genetic elements that encode carbapenem resistance. Seven A. gerneri isolates which harbored Tn2008 (ISAba1-blaOXA-23 -ATPase) and the blaPER-1 gene were also identified. Besides, an A. soil isolate was found to exhibit high-level of meropenem resistance (MIC ≥128 mg/L) and harbor a blaNDM-1 gene located in a core genetic structure of ISAba125-blaNDM-1-ble-trpF-dsbC-cutA. To investigate the genetic relatedness between isolates recovered from hospital sewage and those collected from ICUs, a phylogenetic tree was constructed for 242 clinical isolates and 9 sewage isolates. The results revealed the presence of two evolutionary clades, each containing isolates from both ICU and sewage water, suggesting that CRAB isolates in untreated sewage water were also the transmission clones or closely related evolutionary isolates recoverable in hospital settings. Findings in this work confirm that hospital sewage is a potential reservoir of CRAB.

Development and validation of a random forest algorithm for source attribution of animal and human Salmonella Typhimurium and monophasic variants of S. Typhimurium isolates in England and Wales utilising whole genome sequencing data

Source attribution has traditionally involved combining epidemiological data with different pathogen characterisation methods, including 7-gene multi locus sequence typing (MLST) or serotyping, however, these approaches have limited resolution. In contrast, whole genome sequencing data provide an overview of the whole genome that can be used by attribution algorithms. Here, we applied a random forest (RF) algorithm to predict the primary sources of human clinical Salmonella Typhimurium (S. Typhimurium) and monophasic variants (monophasic S. Typhimurium) isolates. To this end, we utilised single nucleotide polymorphism diversity in the core genome MLST alleles obtained from 1,061 laboratory-confirmed human and animal S. Typhimurium and monophasic S. Typhimurium isolates as inputs into a RF model. The algorithm was used for supervised learning to classify 399 animal S. Typhimurium and monophasic S. Typhimurium isolates into one of eight distinct primary source classes comprising common livestock and pet animal species: cattle, pigs, sheep, other mammals (pets: mostly dogs and horses), broilers, layers, turkeys, and game birds (pheasants, quail, and pigeons). When applied to the training set animal isolates, model accuracy was 0.929 and kappa 0.905, whereas for the test set animal isolates, for which the primary source class information was withheld from the model, the accuracy was 0.779 and kappa 0.700. Subsequently, the model was applied to assign 662 human clinical cases to the eight primary source classes. In the dataset, 60/399 (15.0%) of the animal and 141/662 (21.3%) of the human isolates were associated with a known outbreak of S. Typhimurium definitive type (DT) 104. All but two of the 141 DT104 outbreak linked human isolates were correctly attributed by the model to the primary source classes identified as the origin of the DT104 outbreak. A model that was run without the clonal DT104 animal isolates produced largely congruent outputs (training set accuracy 0.989 and kappa 0.985; test set accuracy 0.781 and kappa 0.663). Overall, our results show that RF offers considerable promise as a suitable methodology for epidemiological tracking and source attribution for foodborne pathogens.

Impact of Homologous Recombination on Core Genome Evolution and Host Adaptation of Pectobacterium parmentieri

Homologous recombination is a major force mechanism driving bacterial evolution, host adaptability, and acquisition of novel virulence traits. Pectobacterium parmentieri is a plant bacterial pathogen distributed worldwide, primarily affecting potatoes, by causing soft rot and blackleg diseases. The goal of this investigation was to understand the impact of homologous recombination on the genomic evolution of P. parmentieri. Analysis of P. parmentieri genomes using Roary revealed a dynamic pan-genome with 3,742 core genes and over 55% accessory genome variability. Bayesian population structure analysis identified 7 lineages, indicating species heterogeneity. ClonalFrameML analysis displayed 5,125 recombination events, with the lineage 4 exhibiting the highest events. fastGEAR analysis identified 486 ancestral and 941 recent recombination events ranging from 43 bp to 119 kb and 36 bp to 13.96 kb, respectively, suggesting ongoing adaptation. Notably, 11% (412 genes) of the core genome underwent recent recombination, with lineage 1 as the main donor. The prevalence of recent recombination (double compared to ancient) events implies continuous adaptation, possibly driven by global potato trade. Recombination events were found in genes involved in vital cellular processes (DNA replication, DNA repair, RNA processing, homeostasis, and metabolism), pathogenicity determinants (type secretion systems, cell-wall degrading enzymes, iron scavengers, lipopolysaccharides (LPS), flagellum, etc.), antimicrobial compounds (phenazine and colicin) and even CRISPR-Cas genes. Overall, these results emphasize the potential role of homologous recombination in P. parmentieri's evolutionary dynamics, influencing host colonization, pathogenicity, adaptive immunity, and ecological fitness.

Population-based genomic surveillance of Pseudomonas aeruginosa causing bloodstream infections in a large Canadian health region

PURPOSE

Population-based surveillance was undertaken to determine clinical factors, susceptibility patterns, and incidence rates (IR) of Pseudomonas aeruginosa causing bloodstream infections (BSIs) in a Canadian region (2010-2018).

METHODS

We combined clinical data with genomics to characterize P. aeruginosa (BSIs) (n = 167) in a well-defined Canadian (Calgary) human population over a 9-year period (2010-2018).

RESULTS

The annual population IR per 100,000 patient years increased from 3.4/100,000 in 2010 to 5.9/100,000 in 2018, with the highest IRs in elderly males from the hospital setting. Over a quarter of patients presented with febrile neutropenia, followed by urinary tract infections and pneumonia. Antimicrobial resistance (AMR) rates and determinants were rare. The P. aeruginosa population was polyclonal consisting of three dominant sequence types (STs), namely ST244, ST111, and ST17. Antimicrobial-susceptible ST244 was the most common clone and belonged to three clades (A, B, C). The ST244 IR/100,000 increased over time due to the expansion of clade C. Multidrug-resistant ST111 was the second most common clone and IR/100,000 decreased over time. ST111 belonged to three clades (A, B, C) with clade C containing blaVIM-2. Different serotypes were linked to various STs. The IR/100,000 of P. aeruginosa that belonged to serotypes O6 increased significantly over time.

CONCLUSION

An effective multivalent vaccine consisting of five serotypes (O1, O3, O5, O6, O11) would confer protection to > 70% of Calgary residents with P. aeruginosa BSIs. This study has provided a unique perspective of the population dynamics over time of P. aeruginosa STs, clades, and serotypes responsible for BSIs.

Natural selection and genetic diversity maintenance in a parasitic wasp during continuous biological control application

Aphidius gifuensis is a parasitoid wasp and primary endoparasitoid enemy of the peach potato aphid, Myzus persicae. Artificially reared, captive wasps of this species have been extensively and effectively used to control populations of aphids and limit crop loss. However, the consequences of large-scale releasing of captive A. gifuensis, such as genetic erosion and reduced fitness in wild populations of this species, remains unclear. Here, we sequence the genomes of 542 A. gifuensis individuals collected across China, including 265 wild and 277 human-intervened samples. Population genetic analyses on wild individuals recovered Yunnan populations as the ancestral group with the most complex genetic structure. We also find genetic signature of environmental adaptation during the dispersal of wild populations from Yunnan to other regions. While comparative genomic analyses of captive wasps revealed a decrease in genetic diversity during long-term rearing, population genomic analyses revealed signatures of natural selection by several biotic (host plants) or abiotic (climate) factors, which support maintenance of the gene pool of wild populations in spite of the introduction of captive wasps. Therefore, the impact of large-scale release is reduced. Our study suggests that A. gifuensis is a good system for exploring the genetic and evolutionary effects of mass rearing and release on species commonly used as biocontrol agents.

Geographical and temporal distribution of multidrug-resistant Salmonella Infantis in Europe and the Americas

Recently emerged S. Infantis strains carrying resistance to several commonly used antimicrobials have been reported from different parts of the globe, causing human cases of salmonellosis and with occurrence reported predominantly in broiler chickens. Here, we performed phylogenetic and genetic clustering analyses to describe the population structure of 417 S. Infantis originating from multiple European countries and the Americas collected between 1985 and 2019. Of these, 171 were collected from 56 distinct premises located in England and Wales (E/W) between 2009 and 2019, including isolates linked to incursions of multidrug-resistant (MDR) strains from Europe associated with imported poultry meat. The analysis facilitated the comparison of isolates from different E/W sources with isolates originating from other countries. There was a high degree of congruency between the outputs of different types of population structure analyses revealing that the E/W and central European (Germany, Hungary, and Poland) isolates formed several disparate groups, which were distinct from the cluster relating to the United States (USA) and Ecuador/Peru, but that isolates from Brazil were closely related to the E/W and the central European isolates. Nearly half of the analysed strains/genomes (194/417) harboured the IncFIB(pN55391) replicon typical of the "parasitic" pESI-like megaplasmid found in diverse strains of S. Infantis. The isolates that contained the IncFIB(pN55391) replicon clustered together, despite originating from different parts of the globe. This outcome was corroborated by the time-measured phylogeny, which indicated that the initial acquisition of IncFIB(pN55391) likely occurred in Europe in the late 1980s, with a single introduction of IncFIB(pN55391)-carrying S. Infantis to the Americas several years later. Most of the antimicrobial resistance (AMR) genes were identified in isolates that harboured one or more different plasmids, but based on the short-read assemblies, only a minority of the resistance genes found in these isolates were identified as being associated with the detected plasmids, whereas the hybrid assemblies comprising the short and long reads demonstrated that the majority of the identified AMR genes were associated with IncFIB(pN55391) and other detected plasmid replicon types. This finding underlies the importance of applying appropriate methodologies to investigate associations of AMR genes with bacterial plasmids.

Genomic analysis of Ralstonia pickettii reveals the genetic features for potential pathogenicity and adaptive evolution in drinking water

Ralstonia pickettii, the most critical clinical pathogen of the genus Ralstonia, has been identified as a causative agent of numerous harmful infections. Additionally, Ralstonia pickettii demonstrates adaptability to extreme environmental conditions, such as those found in drinking water. In this study, we conducted a comprehensive genomic analysis to investigate the genomic characteristics related to potential pathogenicity and adaptive evolution in drinking water environments of Ralstonia pickettii. Through phylogenetic analysis and population genetic analysis, we divided Ralstonia pickettii into five Groups, two of which were associated with drinking water environments. The open pan-genome with a large and flexible gene repertoire indicated a high genetic plasticity. Significant differences in functional enrichment were observed between the core- and pan-genome of different groups. Diverse mobile genetic elements (MGEs), extensive genomic rearrangements, and horizontal gene transfer (HGT) events played a crucial role in generating genetic diversity. In drinking water environments, Ralstonia pickettii exhibited strong adaptability, and the acquisition of specific adaptive genes was potentially facilitated by genomic islands (GIs) and HGT. Furthermore, environmental pressures drove the adaptive evolution of Ralstonia pickettii, leading to the accumulation of unique mutations in key genes. These mutations may have a significant impact on various physiological functions, particularly carbon metabolism and energy metabolism. The presence of virulence-related elements associated with macromolecular secretion systems, virulence factors, and antimicrobial resistance indicated the potential pathogenicity of Ralstonia pickettii, making it capable of causing multiple nosocomial infections. This study provides comprehensive insights into the potential pathogenicity and adaptive evolution of Ralstonia pickettii in drinking water environments from a genomic perspective.

Comparative genomics reveals the correlations of stress response genes and bacteriophages in developing antibiotic resistance of Staphylococcus saprophyticus

IMPORTANCE

Staphylococcus saprophyticus is the second most common bacteria associated with urinary tract infections (UTIs) in women. The antimicrobial treatment regimen for uncomplicated UTI is normally nitrofurantoin, trimethoprim-sulfamethoxazole (TMP-SMX), or a fluoroquinolone without routine susceptibility testing of S. saprophyticus recovered from urine specimens. However, TMP-SMX-resistant S. saprophyticus has been detected recently in UTI patients, as well as in our cohort. Herein, we investigated the understudied resistance patterns of this pathogenic species by linking genomic antibiotic resistance gene (ARG) content to susceptibility phenotypes. We describe ARG associations with known and novel SCCmec configurations as well as phage elements in S. saprophyticus, which may serve as intervention or diagnostic targets to limit resistance transmission. Our analyses yielded a comprehensive database of phenotypic data associated with the ARG sequence in clinical S. saprophyticus isolates, which will be crucial for resistance surveillance and prediction to enable precise diagnosis and effective treatment of S. saprophyticus UTIs.

Genomics analysis of KPC-2 and NDM-5-producing Enterobacteriaceae in migratory birds from Qinghai Lake, China

The study examined the epidemiological characteristics of carbapenem-resistant Enterobacteriaceae (CRE) isolated from migratory birds and surroundings in Qinghai Lake, China. We identified 69 (15.7%) CRE isolates from a total of 439 samples including 29 (6.6%) blaNDM-5 Escherichia coli and 40 (9.1%) blaKPC-2 Klebsiella pneumoniae. WGS analysis indicated that ST746, ST48, ST1011, and ST167 were the primary sequence types (ST) for blaNDM-5 E. coli, while all blaKPC-2 K. pneumoniae were ST11 and harbored numerous antibiotic resistance gene types including blaCTX-M, qnrS, and rmtB. A phylogenetic tree based on core genomes revealed that blaNDM-5 E. coli was highly heterogeneous while the blaKPC-2 K. pneumoniae was highly genetically similar within the group and to human Chinese isolates. IncX3, IncHI2, and IncFIB-HI2 plasmid replicon types were associated with blaNDM-5 spread, while IncFII-R and IncFII plasmids mediated blaKPC-2 spread. We also identified IncFII-R hybrid plasmids most likely formed by IS26-mediated integration of IncFII into IncR plasmid backbones. This also facilitated the persistence of IncFII-R plasmids and antibiotic resistance genes including blaKPC-2. In addition, all of the blaKPC-2 K. pneumoniae isolates harbored a pLVKP-like virulence plasmid carrying a combination of two or more hypervirulence markers that included peg-344, iroB, iucA, rmpA, and rmpA2. This is the first description of ST11 K. pneumoniae that co-carried blaKPC-2- and pLVKP-like virulence plasmids from migratory birds. The blaKPC-2 K. pneumoniae carried by migratory birds displayed high genetic relatedness to human isolates highlighting a high risk of transmission of these K. pneumoniae. KEY POINTS: • Multidrug resistance plasmids (blaKPC-2, bla436NDM-5, bla CTX-M, qnrS, and rmtB). • Co-occurrence of plasmid-mediated resistance and virulence genes. • High similarity between migratory bird genomes and humans.

The Genetic Diversity of White-Backed Planthoppers (Sogatella furcifera) between Myanmar and Yunnan Province of China

In order to clarify the migration route and the source of white-backed planthopper (WBPH) (Sogatella furcifera) between Myanmar and Yunnan Province, China, we collected six populations throughout Myanmar and five populations around the border areas in Yunnan Province, China. A total of 790 base pairs in the mtDNA COI genes from 416 individuals were obtained. A total of 43 haplotypes were identified, among which 37 were unique haplotypes, and the remaining 6 were shared among different populations. Two common shared haplotypes (H_1 and H_2) had a widespread distribution in all populations and accounted for 88.8% of the total haplotype frequency, suggesting a high-level gene flow among the Myanmar and Yunnan populations. Bayesian skyline plot (BSP) analysis results indicated that the effective population size of WBPH expanded between about 10,000 and 7000 years ago, and S. furcifera might follow the post-LGM (Last Glacial Maximum) expansion pattern. Based on the total migrant (Nem) value, it can be deduced that north and northeast Myanmar were the primary migration sources for WBPH populations in the southwest and south Yunnan regions. This study aims to contribute to the sustainable regional management of this important rice pest and provide new insights into the genetic diversity of WBPH in Southeast Asia.

Insights on Pinna nobilis population genetic structure in the Aegean and Ionian Sea

The fan mussel Pinna nobilis Linnaeus, 1758 is an endemic species of the Mediterranean Sea, protected by international agreements. It is one of the largest bivalves in the world, playing an important role in the benthic communities; yet it has been recently characterized as Critically Endangered by the IUCN, due to mass mortality events. In this context, the assessment of the genetic variation of the remaining P. nobilis populations and the evaluation of connectivity among them are crucial elements for the conservation of the species. For this purpose, samples were collected from six regions of the Eastern Mediterranean Sea; the Islands of Karpathos, Lesvos and Crete; the Chalkidiki and Attica Peninsulas; and the Amvrakikos Gulf. Sampling was performed either by collecting tissue from the individuals or by using a non-invasive method, i.e., by scraping the inside of their shells aiming to collect their mucus and thus avoid stress induction to them. Conventional molecular techniques with the use of the COI and 16S rRNA mitochondrial markers were selected for the depiction of the intra-population genetic variability. The analyses included 104 samples from the present study and publicly available sequences of individuals across the whole Mediterranean Sea. The results of this work (a) suggest the use of eDNA as an efficient sampling method for protected bivalves and (b) shed light to the genetic structure of P. nobilis population in the Eastern Mediterranean; this latter knowledge might prove to be fundamental for the species conservation and hence the ecosystem resilience. The haplotype analyses reinforced the evidence that there is a certain degree of connectivity among the distinct regions of the Mediterranean; yet there is evidence of population distinction within the basin, namely between the Western and the Eastern basins. The combination of both genetic markers in the same analysis along with the inclusion of a large number of individuals produced more robust results, revealing a group of haplotypes being present only in the Eastern Mediterranean and providing insights for the species' most suitable conservation management.

An integrated nationwide genomics study reveals transmission modes of typhoid fever in China

IMPORTANCE

Typhoid fever is a life-threatening disease caused by Salmonella enterica serovar Typhi, resulting in a significant disease burden across developing countries. Historically, China was very much close to the global epicenter of typhoid, but the role of typhoid transmission within China and among epicenter remains overlooked in previous investigations. By using newly produced genomics on a national scale, we clarify the complex local and global transmission history of such a notorious disease agent in China spanning the most recent five decades, which largely undermines the global public health network.