Hierarchical and spatially explicit clustering of DNA sequences with BAPS software

The advantage of intergenic regions as genomic features for machine-learning-based host attribution of Salmonella Typhimurium from the USA

Salmonella enterica is a taxonomically diverse pathogen with over 2600 serovars associated with a wide variety of animal hosts including humans, other mammals, birds and reptiles. Some serovars are host-specific or host-restricted and cause disease in distinct host species, while others, such as serovar S. Typhimurium (STm), are generalists and have the potential to colonize a wide variety of species. However, even within generalist serovars such as STm it is becoming clear that pathovariants exist that differ in tropism and virulence. Identifying the genetic factors underlying host specificity is complex, but the availability of thousands of genome sequences and advances in machine learning have made it possible to build specific host prediction models to aid outbreak control and predict the human pathogenic potential of isolates from animals and other reservoirs. We have advanced this area by building host-association prediction models trained on a wide range of genomic features and compared them with predictions based on nearest-neighbour phylogeny. SNPs, protein variants (PVs), antimicrobial resistance (AMR) profiles and intergenic regions (IGRs) were extracted from 3883 high-quality STm assemblies collected from humans, swine, bovine and poultry in the USA, and used to construct Random Forest (RF) machine learning models. An additional 244 recent STm assemblies from farm animals were used as a test set for further validation. The models based on PVs and IGRs had the best performance in terms of predicting the host of origin of isolates and outperformed nearest-neighbour phylogenetic host prediction as well as models based on SNPs or AMR data. However, the models did not yield reliable predictions when tested with isolates that were phylogenetically distinct from the training set. The IGR and PV models were often able to differentiate human isolates in clusters where the majority of isolates were from a single animal source. Notably, IGRs were the feature with the best performance across multiple models which may be due to IGRs acting as both a representation of their flanking genes, equivalent to PVs, while also capturing genomic regulatory variation, such as altered promoter regions. The IGR and PV models predict that ~45 % of the human infections with STm in the USA originate from bovine, ~40 % from poultry and ~14.5 % from swine, although sequences of isolates from other sources were not used for training. In summary, the research demonstrates a significant gain in accuracy for models with IGRs and PVs as features compared to SNP-based and core genome phylogeny predictions when applied within the existing population structure. This article contains data hosted by Microreact.

I Like This New Me: Unravelling Population Structure of Mediterranean Electric Rays and Taxonomic Uncertainties within Torpediniformes

The present study focused on the three species of electric rays known to occur in the Mediterranean Sea: Torpedo torpedo, Torpedo marmorata and Tetronarce nobiliana. Correct identification of specimens is needed to properly assess the impact of fisheries on populations and species. Unfortunately, torpedoes share high morphological similarities, boosting episodes of field misidentification. In this context, genetic data was used (1) to identify specimens caught during fishing operations, (2) to measure the diversity among and within these species, and (3) to shed light on the possible occurrence of additional hidden species in the investigated area. New and already published sequences of COI and NADH2 mitochondrial genes were analyzed, both at a small scale along the Sardinian coasts (Western Mediterranean) and at a large scale in the whole Mediterranean Sea. High levels of genetic diversity were found in Sardinian populations, being significantly different from other areas of the Eastern Mediterranean Sea due to the biotic and abiotic factors here discussed. Sardinian torpedoes can hence be indicated as priority populations/areas to be protected within the Mediterranean Sea. Moreover, sequence data confirmed that only the three species occur in the investigated area. The application of several 'species-delimitation' methods found evidence of cryptic species in the three species outside the Mediterranean Sea, as well as in other genera/families, suggesting the urgent need for future studies and a comprehensive revision of the order Torpediniformes for its effective conservation.

An improved framework for detecting discrete epidemiologically meaningful partitions in hierarchically clustered genetic data

Motivation

Hierarchical clustering of microbial genotypes has the limitation that hierarchical clusters are nested, where smaller groups of related isolates exist within larger groups that get progressively larger as relationships become increasingly distant. In an epidemiologic context, investigators must dissect hierarchical trees into discrete groupings that are epidemiologically meaningful. We recently described a statistical framework (Method A) for dissecting hierarchical trees that attempts to minimize investigator bias. Here, we apply a modified version of that framework (Method B) to a hierarchical tree constructed from 2111 genotypes of the foodborne parasite Cyclospora, including 639 genotypes linked to epidemiologically defined outbreaks. To evaluate Method B's performance, we examined the concordance between these epidemiologically defined groupings and the genetic partitions identified. We also used the same epidemiologic clusters to evaluate the performance of Method A, plus two tree-dissection methods (cutreeHybrid and cutreeDynamic) available within the Dynamic Tree Cut R package, in addition to the TreeCluster method and PARNAS.

Results

Compared to the other methods, Method B, TreeCluster, and PARNAS were the most accurate (99.4%) in identifying genetic groups that reflected the epidemiologic groupings, noting that TreeCluster and PARNAS performed identically on our dataset. CutreeHybrid identified groups reflecting patterns in the wider Cyclospora population structure but lacked finer, strain-level discrimination (Simpson's D: cutreeHybrid=0.785). CutreeDynamic displayed good strain discrimination (Simpson's D = 0.933), though lacked sensitivity (77%). At two different threshold/radius settings TreeCluster/PARNAS displayed similar utility to Method B. However, Method B computes a tree-dissection threshold automatically, and the threshold/radius settings used when executing TreeCluster/PARNAS here were computed using Method B. Using a TreeCluster threshold of 0.045 as recommended in the TreeCluster documentation, epidemiologic utility dropped markedly below that of Method B.

Availability and implementation

Relevant code and data are publicly available. Source code (Method B) and instructions for its use are available here: https://github.com/Joel-Barratt/Hierarchical-tree-dissection-framework.

Deciphering mechanisms of bla NDM gene transmission between human and animals: a genomics study of bacterial isolates from various sources in China, 2015 to 2017

BackgroundIn China, the bla NDM gene has been recovered from human bacterial isolates since 2011. After 2014, detections of this gene in animal and food bacterial isolates have increasingly been reported.AimWe aimed to understand how bla NDM-bearing bacteria could spread between humans, animals, and animal-derived food.MethodsA total of 288 non-duplicate Escherichia coli strains, including 130 bla NDM-carrying and 158 bla NDM-negative strains were collected from clinical (humans), food-producing animals (pigs) and food (retail pork) sources between 2015 and 2017. The strains were whole genome sequenced. Core-genome-multilocus-sequence-typing was conducted. To investigate if sequence types (STs) found in human, animal or food samples could have a prior origin in a clinical, animal or food-borne animal reservoir, discriminant analysis of principal components (DAPC) was used. Plasmids bearing bla NDM were characterised.ResultsThe 130 bla NDM-carrying E. coli strains comprised a total of 60 STs, with ST167 (10/51), ST77 (6/33) and ST48 (6/46) being most prevalent in clinical, animal and food sources, respectively. Some ST10 and ST167 strains were respectively found among all three sources sampled, suggesting they might enable transfer of bla NDM between sources. DAPC analysis indicated possible transmissions of ST167 from humans to animals and ST10 from animals to human. In 114 of 130 bla NDM-carrying isolates, bla NDM was located on an IncX3 plasmid.ConclusionThis study in a Chinese context suggests that cross-species transmission of certain STs of E. coli harbouring bla NDM on mobile elements, may facilitate the spread of carbapenem-resistant Enterobacteriaceae. Stringent monitoring of bla NDM-bearing E. coli in ecosystems is important.

Predictive modeling of antibiotic eradication therapy success for new-onset Pseudomonas aeruginosa pulmonary infections in children with cystic fibrosis

Chronic Pseudomonas aeruginosa (Pa) lung infections are the leading cause of mortality among cystic fibrosis (CF) patients; therefore, the eradication of new-onset Pa lung infections is an important therapeutic goal that can have long-term health benefits. The use of early antibiotic eradication therapy (AET) has been shown to clear the majority of new-onset Pa infections, and it is hoped that identifying the underlying basis for AET failure will further improve treatment outcomes. Here we generated machine learning models to predict AET outcomes based on pathogen genomic data. We used a nested cross validation design, population structure control, and recursive feature selection to improve model performance and showed that incorporating population structure control was crucial for improving model interpretation and generalizability. Our best model, controlling for population structure and using only 30 recursively selected features, had an area under the curve of 0.87 for a holdout test dataset. The top-ranked features were generally associated with motility, adhesion, and biofilm formation.

Molecular Surveillance of MRSA in Raw Milk Provides Insight into MRSA Cross Species Evolution

Methicillin-resistant Staphylococcus aureus (MRSA) in foods has been associated with severe infections in humans and animals worldwide. In the present study, the molecular characteristics of livestock-associated MRSA (LA-MRSA) and human-associated MRSA (hMRSA) isolates obtained in China, as well as MRSA isolates obtained from raw milk in 2018, were investigated. In total, 343 (20.38%; 343/1,683) S. aureus isolates were obtained from 1,683 raw milk samples from 100 dairy farms in 11 provinces across China. Among these, 49 (2.91%; 49/1,683) were mecA-positive MRSA. All LA-MRSA isolates were resistant to penicillin and highly resistant to erythromycin, sulfisoxazole, and clindamycin. Bioinformatic analysis the 49 genomes of LA-MRSA and 71 previously published hMRSA genomes isolated from Chinese individuals in 2018 indicated that blaZ, erm, ant(6)-Ia, aph(3')-III, tet(K), cat, and aph(2″)-Ia were more prevalent in MRSA from raw milk (P < 0.05) compared to hMRSA. Additionally, hMRSA isolates were more significantly associated with ST5 (P < 0.01) compared to LA-MRSA; in contrast, ST338 was more prevalent among LA-MRSA isolates (P < 0.01). Likewise, the SCCmec type II was only detected in hMRSA isolates, whereas SCCmec type V and IV were more prevalent among LA-MRSA (P < 0.01). Furthermore, core-genome phylogenetic analysis showed the endemic characteristics of LA-MRSA in local provinces, as well as the close evolutionary relationships between MRSA from cattle and humans. Finally, homology analysis of mecA and blaZ genetic contexts revealed a high possibility of horizontal transmission of MRSA resistance genes among raw milk-associated and hMRSA strains, which increases the risk for public health. IMPORTANCE Methicillin-resistant Staphylococcus aureus (MRSA) is considered a public health concern as it is resistant to multiple antibiotics, thus being in zoonotic transmission of antibiotic resistance genes. MRSA causes serious public health issues and leads to hard-to-treat infections in humans and animals; therefore, it was meaningful to determine the prevalence of MRSA in raw milk samples and investigate phenotype and genotype of antimicrobial resistance and molecular characteristics in livestock-associated MRSA (LA-MRSA) and human-associated MRSA (hMRSA) in China, which could provide a theoretical basis for preventing and controlling the spread of MRSA between livestock and humans.

Mesomycoplasma ovipneumoniae from goats with respiratory infection: pathogenic characteristics, population structure, and genomic features

BACKGROUND

Mycoplasma ovipneumoniae is a critical pathogen that causes respiratory diseases that threaten Caprini health and cause economic damage. A genome-wide study of M. ovipneumoniae will help understand the pathogenic characteristics of this microorganism.

RESULTS

Toxicological pathology and whole-genome sequencing of nine M. ovipneumoniae strains isolated from goats were performed using an epidemiological survey. These strains exhibited anterior ventral lung consolidation, typical of bronchopneumonia in goats. Average nucleotide identity and phylogenetic analysis based on whole-genome sequences showed that all M. ovipneumoniae strains clustered into two clades, largely in accordance with their geographical origins. The pan-genome of the 23 M. ovipneumoniae strains contained 5,596 genes, including 385 core, 210 soft core, and 5,001 accessory genes. Among these genes, two protein-coding genes were annotated as cilium adhesion and eight as paralog surface adhesins when annotated to VFDB, and no antibiotic resistance-related genes were predicted. Additionally, 23 strains carried glucosidase-related genes (ycjT and group_1595) and glucosidase-related genes (atpD_2), indicating that M. ovipneumoniae possesses a wide range of glycoside hydrolase activities.

CONCLUSIONS

The population structure and genomic features identified in this study will facilitate further investigations into the pathogenesis of M. ovipneumoniae and lay the foundation for the development of preventive and therapeutic methods.

The persistence of time: the lifespan of Bacillus anthracis spores in environmental reservoirs

Anthrax is a lethal bacterial zoonosis primarily affecting herbivorous wildlife and livestock. Upon host death Bacillus anthracis vegetative cells form spores capable of surviving for years in soil. Anthrax transmission requires host exposure to large spore doses. Thus, conditions that facilitate higher spore concentrations or promote spore survival will increase the probability that a pathogen reservoir infects future hosts. We investigated abiotic and pathogen genomic variation in relation to spore concentrations in surface soils (0-1 cm depth) at 40 plains zebra (Equus quagga) anthrax carcass sites in Namibia. Specifically, how initial spore concentrations and spore survival were affected by seasonality associated with the timing of host mortality, local soil characteristics, and pathogen genomic variation. Zebras dying of anthrax in wet seasons-the peak season for anthrax in Etosha National Park-had soil spore concentrations 1.36 orders of magnitude higher than those that died in dry seasons. No other variables considered affected spore concentrations, and spore survival rates did not differ among sites. Surface soils at these pathogen reservoirs remained culture positive for a range of 3.8-10.4 years after host death. Future research could evaluate if seasonal patterns in spore concentrations are driven by differences in sporulation success or levels of terminal bacteremia.

Barriers to chimpanzee gene flow at the south-east edge of their distribution

Populations on the edge of a species' distribution may represent an important source of adaptive diversity, yet these populations tend to be more fragmented and are more likely to be geographically isolated. Lack of genetic exchanges between such populations, due to barriers to animal movement, can not only compromise adaptive potential but also lead to the fixation of deleterious alleles. The south-eastern edge of chimpanzee distribution is particularly fragmented, and conflicting hypotheses have been proposed about population connectivity and viability. To address this uncertainty, we generated both mitochondrial and MiSeq-based microsatellite genotypes for 290 individuals ranging across western Tanzania. While shared mitochondrial haplotypes confirmed historical gene flow, our microsatellite analyses revealed two distinct clusters, suggesting two populations currently isolated from one another. However, we found evidence of high levels of gene flow maintained within each of these clusters, one of which covers an 18,000 km2 ecosystem. Landscape genetic analyses confirmed the presence of barriers to gene flow with rivers and bare habitats highly restricting chimpanzee movement. Our study demonstrates how advances in sequencing technologies, combined with the development of landscape genetics approaches, can resolve ambiguities in the genetic history of critical populations and better inform conservation efforts of endangered species.

Whole-genome sequencing of Chlamydia psittaci from Australasian avian hosts: A genomics approach to a pathogen that still ruffles feathers

Chlamydia psittaci is a globally distributed veterinary pathogen with zoonotic potential. Although C. psittaci infections have been reported in various hosts, isolation and culture of Chlamydia is challenging, hampering efforts to produce contemporary global C. psittaci genomes. This is particularly evident in the lack of avian C. psittaci genomes from Australia and New Zealand. In this study, we used culture-independent probe-based whole-genome sequencing to expand the global C. psittaci genome catalogue. Here, we provide new C. psittaci genomes from two pigeons, six psittacines, and novel hosts such as the Australian bustard (Ardeotis australis) and sooty shearwater (Ardenna grisea) from Australia and New Zealand. We also evaluated C. psittaci genetic diversity using multilocus sequence typing (MLST) and major outer membrane protein (ompA) genotyping on additional C. psittaci-positive samples from various captive avian hosts and field isolates from Australasia. We showed that the first C. psittaci genomes sequenced from New Zealand parrots and pigeons belong to the clonal sequence type (ST)24 and diverse 'pigeon-type' ST27 clade, respectively. Australian parrot-derived strains also clustered in the ST24 group, whereas the novel ST332 strain from the Australian bustard clustered in a genetically diverse clade of strains from a fulmar, parrot, and livestock. MLST and ompA genotyping revealed ST24/ompA genotype A in wild and captive parrots and a sooty shearwater, whilst 'pigeon-types' (ST27/35 and ompA genotypes B/E) were found in pigeons and other atypical hosts, such as captive parrots, a little blue penguin/Kororā (Eudyptula minor) and a zebra finch (Taeniopygia guttata castanotis) from Australia and New Zealand. This study provides new insights into the global phylogenomic diversity of C. psittaci and further demonstrates the multi-host generalist capacity of this pathogen.

Plasmid permissiveness of wastewater microbiomes can be predicted from 16S rRNA sequences by machine learning

MOTIVATION

Wastewater treatment plants (WWTPs) harbor a dense and diverse microbial community. They constantly receive antimicrobial residues and resistant strains, and therefore provide conditions for horizontal gene transfer (HGT) of antimicrobial resistance (AMR) determinants. This facilitates the transmission of clinically important genes between, e.g. enteric and environmental bacteria, and vice versa. Despite the clinical importance, tools for predicting HGT remain underdeveloped.

RESULTS

In this study, we examined to which extent water cycle microbial community composition, as inferred by partial 16S rRNA gene sequences, can predict plasmid permissiveness, i.e. the ability of cells to receive a plasmid through conjugation, based on data from standardized filter mating assays using fluorescent bio-reporter plasmids. We leveraged a range of machine learning models for predicting the permissiveness for each taxon in the community, representing the range of hosts a plasmid is able to transfer to, for three broad host-range resistance IncP plasmids (pKJK5, pB10, and RP4). Our results indicate that the predicted permissiveness from the best performing model (random forest) showed a moderate-to-strong average correlation of 0.49 for pB10 [95% confidence interval (CI): 0.44-0.55], 0.43 for pKJK5 (0.95% CI: 0.41-0.49), and 0.53 for RP4 (0.95% CI: 0.48-0.57) with the experimental permissiveness in the unseen test dataset. Predictive phylogenetic signals occurred despite the broad host-range nature of these plasmids. Our results provide a framework that contributes to the assessment of the risk of AMR pollution in wastewater systems.

AVAILABILITY AND IMPLEMENTATION

The predictive tool is available as an application at https://github.com/DaneshMoradigaravand/PlasmidPerm.

To Be, or Not to Be: That Is the Hamletic Question of Cryptic Evolution in the Eastern Atlantic and Mediterranean Raja miraletus Species Complex

Despite a high species diversity, skates (Rajiformes) exhibit remarkably conservative morphology and ecology. Limited trait variations occur within and between species, and cryptic species have been reported among sister and non-sister taxa, suggesting that species complexes may be subject to stabilising selection. Three sibling species are currently recognised in the Raja miraletus complex: (i) R. miraletus occurring along the Portuguese and Mediterranean coasts, (ii) R. parva in the Central-Eastern Atlantic off West Africa and (iii) R. ocellifera in the Western Indian Ocean off South Africa. In the present study, the genetic variation at mitochondrial and nuclear markers was estimated in the species complex by analysing 323 individuals sampled across most of its geographical distribution area to test the hypothesis that restricted gene flow and genetic divergence within species reflect known climate and bio-oceanographic discontinuities. Our results support previous morphological studies and confirm the known taxonomic boundaries of the three recognised species. In addition, we identified multiple weakly differentiated clades in the Northeastern Atlantic Ocean and Mediterranean, at least two additional cryptic taxa off Senegal and Angola, a pronounced differentiation of ancient South African clades. The hidden genetic structure presented here may represent a valuable support to species' conservation action plans.

Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa

Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.

Preliminary population studies of the grassland swallowtail butterfly Euryades corethrus (Lepidoptera, Papilionidae)

Euryades corethrus is a Troidini butterfly (Papilionidae, Papilioninae), endemic to grasslands in southern Brazil, Uruguay, Argentina and Paraguay. Formerly abundant, nowadays it is in the Red list of endangered species for those areas. During its larval stage, it feeds on Aristolochia spp, commonly found in southern grasslands. These native grassland areas are diminishing, being converted to crops and pastures, causing habitat loss for Aristolochia and E. corethrus. This study aimed to assess the genetic diversity, population structure and demographic history of E. corethrus. We sampled eight populations from Rio Grande do Sul, Brazil and based on Cytochrome Oxidase subunit I (COI) molecular marker, our results suggest a low genetic variability between populations, presence of gene flow and, consequently, lack of population structure. A single maternally inherited-genetic marker is insufficient for population-level decisions, but barcoding is a useful tool during early stages of population investigation, bringing out genomic diversity patterns within the target species. Those populations likely faced a bottleneck followed by a rapid expansion during the last glaciation and subsequent stabilization in effective population size. Habitat loss is a threat, which might cause isolation, loss of genetic variability and, ultimately, extinction of E. corethrus if no habitat conservation policy is adopted.

Highly drug resistant clone of Salmonella Kentucky ST198 in clinical infections and poultry in Zimbabwe

A highly multidrug-resistant strain of Salmonella enterica serotype Kentucky (S. Kentucky) of sequence type (ST)198 emerged in North Africa and has since spread widely. To investigate the genetic diversity and phylogenetic relationship of S. Kentucky in Zimbabwe and identify potential sources of infection, the whole-genome sequence of 37 S. Kentucky strains isolated from human clinical infections and from poultry farms between 2017 and 2020 was determined. Of 37 S. Kentucky isolates, 36 were ST198 and one was ST152. All ST198 isolates had between six and fifteen antimicrobial resistance (AMR) genes, and 92% carried at least ten AMRs. All ST198 isolates harbored the Salmonella genomic island K-Israel variant (SGI1-KIV) integrated into the chromosome with aac(3)-ld, aac(6)-laa, aadA7, blaTEM-1, sul1, and tetA genes, with occasional sporadic loss of one or more genes noted from five isolates. All ST198 isolates also had mutations in the quinolone resistance-determining region of the gyrA and parC genes. The blaCTX-M-14.1 and fosA3 genes were present in 92% of the ST198 isolates, conferring resistance to extended-spectrum cephalosporins and fosfomycin, respectively, were present on an IncHI2 plasmid with the aadA2b, aadA1, aph(3')-Ib, aph(6')-Id, cmlA1 and sul3 AMR genes. S. Kentucky ST198 isolates from Zimbabwe formed a closely related phylogenetic clade that emerged from a previously reported global epidemic population. The close genetic relationship and population structure of the human clinical and poultry isolates of ST198 in Zimbabwe are consistent with poultry being an important source of highly resistant strains of S. Kentucky in Zimbabwe.

A comparative study of intestinal Pseudomonas aeruginosa in healthy individuals and ICU inpatients

The human intestinal tract is considered the most important reservoir of the opportunistic pathogens, including Pseudomonas aeruginosa, which is often overlooked but critical due to its antimicrobial resistance and virulence. Public health interventions to control this pathogen require a comprehensive understanding of its epidemiology and genomics.  In the current study, we identified P. aeruginosa strains from 2,605 fecal samples collected between 2021 to 2022. Among these samples, 574 were from ICU inpatients in Zhejiang province, while 2,031 were obtained from healthy individuals residing in ten different provinces in China. The prevalence of P. aeruginosa intestinal carriage was found to be higher in ICU inpatients (10.28%, 95% CI: 7.79%-12.76%) than that in healthy individuals (3.99%, 81/2,031, 95% CI: 3.14%-4.84%). Similarly, the prevalence of carbapenem-resistant P. aeruginosa (CRPA) was higher in ICU inpatients (32.2%) compared to healthy individuals (7.41%). The population structure analysis of our isolates revealed a predominantly non-clonal distribution, with 41 distinct sequence types identified among 59 P. aeruginosa isolates from ICU inpatients and 38 different STs among 81 P. aeruginosa isolates from healthy individuals. These findings suggest that the individual acquisition of P. aeruginosa is more frequent than patient-to-patient transmission, as evidenced by the polyclonal population structure. Antimicrobial susceptibility testing and genome analysis indicated that P. aeruginosa strains from ICU inpatients exhibited significantly higher resistance rates to most antimicrobials and harbored a greater number of acquired resistance genes compared to strains from healthy individuals. Notably, in ICU inpatients, we identified three isolates of ST463, all of which shared the conserved Tn3-TnpR-ISKpn8-blaKPC-ISKpn6 genetic context. Additionally, five isolates carrying the qacE gene were also identified, these findings suggest that small-scale transmission events may still occur within the ICU setting, posing significant challenges for clinical management. With regard to virulence factors, we observed similar profiles between the two groups, except for phzA2, phzB2, and pilA, which were statistically higher in isolates from healthy individuals. This may be because the accumulating resistance mutations in ICU-derived P. aeruginosa are linked to a decrease in virulence.

Population genomics and pathotypic evaluation of the bacterial leaf blight pathogen of rice reveals rapid evolutionary dynamics of a plant pathogen

The Xanthomonas oryzae pv. oryzae (Xoo) is a bacterial pathogen causing bacterial blight disease in rice, resulting in significant yield reductions of up to 50% in rice production. Despite its serious threat to food production globally, knowledge of its population structure and virulence evolution is relatively limited. In this study, we employed whole-genome sequencing to explore the diversity and evolution of Xoo in the main rice-growing areas of China over the past 30 years. Using phylogenomic analysis, we revealed six lineages. CX-1 and CX-2 primarily contained Xoo isolates from South China, while CX-3 represented Xoo isolates from North China. Xoo isolates belonging to CX-5 and CX-6 were the most prevalent across all studied areas, persisting as dominant lineages for several decades. Recent sporadic disease outbreaks were primarily caused by Xoo isolates derived from the two major lineages, CX-5 and CX-6, although Xoo isolates from other lineages also contributed to these outbreaks. The lineage and sub-lineage distributions of Xoo isolates were strongly correlated with their geographical origin, which was found to be mainly determined by the planting of the two major rice subspecies, indica and japonica. Moreover, large-scale virulence testing was conducted to evaluate the diversity of pathogenicity for Xoo. We found rapid virulence evolution against rice, and its determinant factors included the genetic background of Xoo, rice resistance genes, and planting environment of rice. This study provides an excellent model for understanding the evolution and dynamics of plant pathogens in the context of their interactions with their hosts, which are shaped by a combination of geographical conditions and farming practices. The findings of this study may have important implications for the development of effective strategies for disease management and crop protection in rice production systems.

Genomic and Phylogenetic Analysis of a Multidrug-Resistant blaNDM-carrying Klebsiella michiganensis in China

Objective

Klebsiella michiganensis is an emerging hospital-acquired bacterial pathogen. However, there is a dearth of knowledge on the antimicrobial resistance and transmission of K. michiganensis. Here, we characterized the microbiological and genomic features of a carbapenem-resistant K. michiganensis strain harboring the bla_NDM-1 gene in China.

Methods

K. michiganensis strain 2563 was recovered from the sputum sample of a hospitalized patient with pulmonary infection. Whole-genome sequencing of K. michiganensis strain 2563 was conducted using both the short-read Illumina and long-read MinION platforms to thoroughly characterize the genetic context of bla_NDM-carrying plasmid in K. michiganensis 2563. Furthermore, BacWGSTdb server was utilized to perform in silico multilocus sequence typing (MLST), identify antimicrobial resistance genes, and conduct genomic epidemiological analyses of the closely related isolates deposited in the public database.

Results

K. michiganensis 2563 was resistant to piperacillin, aztreonam, meropenem, imipenem, amoxicillin-clavulanic acid, ampicillin, cefotaxime, cefazolin, ampicillin/sulbactam, cefepime, piperacillin-tazobactam, and ceftazidime. It belonged to sequence type (ST) 43, and the bla_NDM-1 gene was found to be located on the plasmid p2563_NDM (54,035 bp). This plasmid showed remarkable similarity to other bla_NDM-1-encoding plasmids found in various Enterobacterium species in the public database. The occurrence of global ST43 K. michiganensis was primarily sporadic, and the closest relative of K. michiganensis 2563 was another ST43 isolate 12,084 recovered from China in 2013, which differed by 171 SNPs.

Conclusion

Our study reports the genome characteristics of a carbapenem-resistant K. michiganensis strain carrying the bla_NDM-1 gene in China, highlighting the need for ongoing surveillance of this pathogen in clinical settings.

Out of East Asia: Early Warning of the Possible Invasion of the Important Bean Pest Stalk-Eyed Seed Bug Chauliops fallax (Heteroptera: Malcidae: Chauliopinae)

The short stay at the beginning of the invasion process is a critical time for invasive species identification and preventing invasive species from developing a wider distribution and significant economic impact. The stalk-eyed seed bug Chauliops fallax is an important agricultural pest of soybean and was first reported to occur outside East Asia. Here, we reported the native evolutionary history, recent invasion history, and potential invasion threats of C. fallax for the first time based on population genetic methods and ecological niche modelling. The results showed that four native East Asian genetic groups (EA, WE, TL, and XZ) were well supported, showing an east-west differentiation pattern consistent with the geographical characteristics of three-step landforms in China. Two main haplotypes existed: Hap1 might have experienced a rapid northwards expansion process after the LGM period, and Hap5 reflected local adaptation to the environment in southeastern China. The Kashmir sample was found to come from the recent invasion of populations in the coastal areas of southern China. Ecological niche modelling results suggested that North America has a high risk of invasion, which might pose a serious threat to local soybean production. In addition, with future global warming, the suitable habitat in Asia will move towards the higher latitude region and gradually deviate from the soybean planting area, which indicates the threat of C. fallax to soybean production in Asia will decrease in the future. The results could provide new insights into the monitoring and management of this agricultural pest in the early invasion stage.

Prevalence and molecular characteristics of mcr-1-positive Escherichia coli isolated from duck farms and the surrounding environments in coastal China

The emergence of colistin-resistance is considered a threat to public health and colistin-resistant bacteria have recently been reported in animal, environmental and human sources. Whereas, the epidemic and dissemination of colistin-resistant bacteria in duck farms have not been surveyed, especially the surrounding environmental contamination from duck farms. We investigated the prevalence and molecular characteristics of mcr-1-positive E. coli from duck farms in coastal China. 360 mcr-1-positive E. coli isolates were collected from 1112 samples from duck farms and surrounding environments. The prevalence of mcr-1-positive E. coli in Guangdong province was higher than other two provinces we examined. PFGE analysis indicated clonal spread of mcr-1-positive E. coli between duck farms and surrounding environments, including water and soil. MLST analysis demonstrated that ST10 was more common than ST1011, ST117, and ST48. Phylogenomic analysis also suggested mcr-1-positive E. coli collected from distinct cities were assigned to the same lineage and mcr-1 was primarily located on IncI2 and IncHI2 plasmids. Genomic environment analysis showed mobile gene elements ISApl1 most likely plays a key role in the horizontal transmission of mcr-1. WGS further revealed that mcr-1 was found associated with 27 different ARGs. Our findings emphasize the urgent need for effective colistin resistance surveillance in humans, animals and the environment.

Genome-wide association reveals host-specific genomic traits in Escherichia coli

BACKGROUND

Escherichia coli is an opportunistic pathogen which colonizes various host species. However, to what extent genetic lineages of E. coli are adapted or restricted to specific hosts and the genomic determinants of such adaptation or restriction is poorly understood.

RESULTS

We randomly sampled E. coli isolates from four countries (Germany, UK, Spain, and Vietnam), obtained from five host species (human, pig, cattle, chicken, and wild boar) over 16 years, from both healthy and diseased hosts, to construct a collection of 1198 whole-genome sequenced E. coli isolates. We identified associations between specific E. coli lineages and the host from which they were isolated. A genome-wide association study (GWAS) identified several E. coli genes that were associated with human, cattle, or chicken hosts, whereas no genes associated with the pig host could be found. In silico characterization of nine contiguous genes (collectively designated as nan-9) associated with the human host indicated that these genes are involved in the metabolism of sialic acids (Sia). In contrast, the previously described sialic acid regulon known as sialoregulon (i.e. nanRATEK-yhcH, nanXY, and nanCMS) was not associated with any host species. In vitro growth experiments with a Δnan-9 E. coli mutant strain, using the sialic acids 5-N-acetylneuraminic acid (Neu5Ac) and N-glycolylneuraminic acid (Neu5Gc) as sole carbon source, showed impaired growth behaviour compared to the wild-type.

CONCLUSIONS

This study provides an extensive analysis of genetic determinants which may contribute to host specificity in E. coli. Our findings should inform risk analysis and epidemiological monitoring of (antimicrobial resistant) E. coli.

The evolution and international spread of extensively drug resistant Shigella sonnei

Shigella sonnei causes shigellosis, a severe gastrointestinal illness that is sexually transmissible among men who have sex with men (MSM). Multidrug resistance in S. sonnei is common including against World Health Organisation recommended treatment options, azithromycin, and ciprofloxacin. Recently, an MSM-associated outbreak of extended-spectrum β-lactamase producing, extensively drug resistant S. sonnei was reported in the United Kingdom. Here, we aimed to identify the genetic basis, evolutionary history, and international dissemination of the outbreak strain. Our genomic epidemiological analyses of 3,304 isolates from the United Kingdom, Australia, Belgium, France, and the United States of America revealed an internationally connected outbreak with a most recent common ancestor in 2018 carrying a low-fitness cost resistance plasmid, previously observed in travel associated sublineages of S. flexneri. Our results highlight the persistent threat of horizontally transmitted antimicrobial resistance and the value of continuing to work towards early and open international sharing of genomic surveillance data.

Small population of the largest water strider after the late Pleistocene and the implications for its conservation

Climate oscillation and its synergistic impacts on habitat fragmentation have been identified as threatening the survival of some extant species. However, the mechanisms by which semi-aquatic insects impacted by such events remain poorly understood. Herein, we studied the largest water strider in the world, Gigantometra gigas, to explore the effect of these two factors on its evolutionary history. The sequences of mitogenomic and nrDNA cluster were utilized to reconstruct phylogenetic relationship among G. gigas populations and its demographic history. Mitochondrial genes were separately reconstructed topologies of that populations and detected remarkable differences. We found that G. gigas populations conform to the isolation-by-distance model, and decline occurred at about 120 ka, which was probably influenced by the climate change during the late Pleistocene and eventually maintained a small effective population size (Ne) around 85,717. The populations in Guangdong Province of China are worthy of note in that they exhibit low genetic diversity, a small Ne around 18,899 individuals, and occupy an area with little suitable future habitat for G. gigas. This work recommends that conservation efforts are implemented to ensure the long-term survival of small G. gigas populations, and notes that further evaluation of their extinction risk under the impacts of human activities is required.

Climatic oscillation promoted diversification of spinous assassin bugs during Pleistocene glaciation

Insect speciation is among the most fascinating topics in evolutionary biology; however, its underlying mechanisms remain unclear. Allopatric speciation represents one of the major types of speciation and is believed to have frequently occurred during glaciation periods, when climatic oscillation may have caused suitable habitats to be fragmented repeatedly, creating geographical isolation among populations. However, supporting evidence for allopatric speciation of insects in East Asia during the Pleistocene glaciation remains lacking. We aim to investigate the effect of climatic oscillation during the Pleistocene glaciation on the diversification pattern and evolutionary history of hemipteran insects and to test the hypothesis of Pleistocene species stability using spinous assassin bugs Sclomina (Hemiptera: Reduviidae), a small genus widely distributed in southern China but was later found to have cryptic species diversity. Here, using the whole mitochondrial genome (mitogenome) and nuclear ribosomal RNA genes, we investigated both interspecific and intraspecific diversification patterns of spinous assassin bugs. Approximate Bayesian computation, ecological niche modeling, and demographic history analyses were also applied to understand the diversification process and driven factors. Our data suggest that the five species of Sclomina are highly diverged, despite three of them currently being cryptic. Speciation occurred during the Pleistocene when suitable distribution areas were possibly fragmented. Six phylogeographic groups in the type species S. erinacea were identified, among which two groups underwent expansion during the early Last Glacial Period and after Last Glacier Maximum. Our analyses suggest that this genus may have experienced climate-driven habitat fragmentation and postglacial expansion in the Pleistocene, promoting allopatric speciation and intraspecific diversification. Our results reveal underestimated species diversity in a small insect group and illustrate a remarkable example of allopatric speciation of insects in East Asia promoted by Pleistocene climatic oscillations. These findings provide important insights into the speciation processes and aid the conservation of insect species diversity.

Weak genetic structure of flower thrips Frankliniella intonsa in China revealed by mitochondrial genomes

Human activities facilitate long-distance dispersal of insects beyond their native range. In particular, the transportation of live plants offers diffusion opportunities for some insects with weak flight abilities. The increase in urban afforestation also help insect reside in urban habitats. The flower thrips, Frankliniella intonsa, is a widespread pest that causes serious damage to many economically important plants. Human activities are likely to facilitate the dispersal of this pest, however, the population genetic structure of this pest remains unclear. Herein, high-throughput sequencing was used to obtain 149 whole mitochondrial genomes of flower thrips from 28 geographic populations in China. Population genetic analyses, phylogenetic reconstruction, and inference of demographic history were then performed. A weak genetic structure was found among all populations across large geographic distance in China, in which five mitochondrial haplotype lineages were resolved. One of the lineages was identified to be shared among most samples collected from central city areas, which may be derived from the surrounding areas. Demographic history analyses suggested a recent population expansion of F. intonsa. Overall, the present population genetic structure of flower thrips in China may be promoted by human-mediated urban afforestation across the country.

A phylogenomic analysis of Limosilactobacillus reuteri reveals ancient and stable evolutionary relationships with rodents and birds and zoonotic transmission to humans

BACKGROUND

Gut microbes play crucial roles in the development and health of their animal hosts. However, the evolutionary relationships of gut microbes with vertebrate hosts, and the consequences that arise for the ecology and lifestyle of the microbes are still insufficiently understood. Specifically, the mechanisms by which strain-level diversity evolved, the degree by which lineages remain stably associated with hosts, and how their evolutionary history influences their ecological performance remain a critical gap in our understanding of vertebrate-microbe symbiosis.

RESULTS

This study presents the characterization of an extended collection of strains of Limosilactobacillus reuteri and closely related species from a wide variety of hosts by phylogenomic and comparative genomic analyses combined with colonization experiments in mice to gain insight into the long-term evolutionary relationship of a bacterial symbiont with vertebrates. The phylogenetic analysis of L. reuteri revealed early-branching lineages that primarily consist of isolates from rodents (four lineages) and birds (one lineage), while lineages dominated by strains from herbivores, humans, pigs, and primates arose more recently and were less host specific. Strains from rodent lineages, despite their phylogenetic divergence, showed tight clustering in gene-content-based analyses. These L. reuteri strains but not those ones from non-rodent lineages efficiently colonize the forestomach epithelium of germ-free mice. The findings support a long-term evolutionary relationships of L. reuteri lineages with rodents and a stable host switch to birds. Associations of L. reuteri with other host species are likely more dynamic and transient. Interestingly, human isolates of L. reuteri cluster phylogenetically closely with strains from domesticated animals, such as chickens and herbivores, suggesting zoonotic transmissions.

CONCLUSIONS

Overall, this study demonstrates that the evolutionary relationship of a vertebrate gut symbiont can be stable in particular hosts over time scales that allow major adaptations and specialization, but also emphasizes the diversity of symbiont lifestyles even within a single bacterial species. For L. reuteri, symbiont lifestyles ranged from autochthonous, likely based on vertical transmission and stably aligned to rodents and birds over evolutionary time, to allochthonous possibly reliant on zoonotic transmission in humans. Such information contributes to our ability to use these microbes in microbial-based therapeutics.

Genetic variability highlights the invasion route of the Lutzomyia longipalpis complex, the main vector of Visceral Leishmaniasis in Uruguay

In the Americas, the sandfly Lutzomyia longipalpis is the main vector of the parasitic protozoa Leishmania infantum, the etiological agent of visceral leishmaniasis (VL). The Lu. longipalpis species complex is currently discontinuously distributed across the Neotropical region, from Mexico to the north of Argentina and Uruguay. During its continental spreading, it must have adapted to several biomes and temperature amplitudes, when founder events should have contributed to the high genetic divergence and geographical structure currently observed, reinforcing the speciation process. The first report of Lu. longipalpis in Uruguay was in 2010, calling the attention of Public Health authorities. Five years later, the parasite Le. infantum was recorded and in 2015 the first case of VL in canids was reported. Hitherto seven human cases by VL have been reported in Uruguay. Here, we publish the first DNA sequences from the mitochondrial genes ND4 and CYTB of Lu. longipalpis collected in Uruguay, and we used these molecular markers to investigate their genetic variability and population structure. We described four new ND4 haplotypes in a total of 98 (4/98) and one CYTB in a total of 77 (1/77). As expected, we were able to establish that the Lu. longipalpis collected in two localities (i.e. Salto and Bella Unión) from the north of Uruguay are closely related to the populations from neighbouring countries. We also propose that the possible route for the vector arrival to the region may have been through vegetation and forest corridors of the Uruguay River system, as well as it may have benefited from landscape modifications generated by commercial forestation. The ecological-scale processes shaping Lu. longipalpis populations, the identification of genetically homogeneous groups and the gene flow among them must be carefully investigated by using highly sensible molecular markers (i.e. genome wide SNPs) since it will help to the understanding of VL transmission and contribute to the planification of public policies on its control.

Aeromonas spp. from hospital sewage act as a reservoir of genes resistant to last-line antibiotics

BACKGROUND

Aeromonas species are opportunistic pathogens distributed widely in the ecosystem. They are known to be capable of acquiring antibiotic resistance genes, including those encoding proteins against last-line antibiotics, such as the tmexCD-toprJ, mcr and carbapenemase genes. We investigated the genomic and phenotypic characteristics of tmexCD-toprJ-positive Aeromonas strains collected from human, animals, and water samples, particularly those from hospital wastewater in China.

METHODS

Samples were collected from living animals, meat, water and human. Aeromonas strains in these samples were isolated in selective media. Antimicrobial resistance profiles of all Aeromonas strains were tested by the broth microdilution method. The presence of tmexCD-toprJ was verified by polymerase chain reaction (PCR). All tmexCD-toprJ-positive (n = 36) and selected tmexCD-toprJ-negative (n = 18) Aeromonas strains were subjected to whole genome sequencing. Carriage of antimicrobial resistance genes, the genetic environment of tmexCD-toprJ and genetic diversity of tmexCD-toprJ-positive Aeromonas strains were determined by bioinformatics analysis. Phylogenetic tree of the Aeromonas strains was built by using the Harvest Suite.

FINDINGS

Among the 636 Aeromonas strains isolated from different sources, 36 were positive for tmexCD-toprJ, with the highest prevalence of tmexCD-toprJ being found in fishes (8.8%, 95 CI% 3.6-17.2%), followed by hospital wastewater (6.5%, 95 CI% 4.3-9.3%), river water (2.0%, 0.1-10.9) and duck (1.2%, 95 CI% 3.6-17.2%). All tmexCD-toprJ-positive Aeromonas strains carried multiple antimicrobial resistance genes and exhibited resistance to different classes of antibiotics. Co-existence of tmexCD-toprJ, mcr and bla_KPC-2 were identified in 21 strains. The tmexCD-toprJ-positive Aeromonas strains were genetically diverse and found to belong to four different species that could be clustered into three major lineages. The tmexCD-toprJ gene clusters were predominantly located in the chromosome (35/36) of Aeromonas spp., with only one strain carrying the plasmid-borne tmexCD-toprJ cluster. The tmexCD-toprJ genes were associated with seven different types of genetic environments, each of which carried distinct types of mobile elements that may be responsible for mediating transmission of this gene cluster.

Multiple domestication events explain the origin of Gossypium hirsutum landraces in Mexico

Several Mesoamerican crops constitute wild-to-domesticated complexes generated by multiple initial domestication events, and continuous gene flow among crop populations and between these populations and their wild relatives. It has been suggested that the domestication of cotton (Gossypium hirsutum) started in the northwest of the Yucatán Peninsula, from where it spread to other regions inside and outside of Mexico. We tested this hypothesis by assembling chloroplast genomes of 23 wild, landraces, and breeding lines (transgene-introgressed and conventional). The phylogenetic analysis showed that the evolutionary history of cotton in Mexico involves multiple events of introgression and genetic divergence. From this, we conclude that Mexican landraces arose from multiple wild populations. Our results also revealed that their structural and functional chloroplast organizations had been preserved. However, genetic diversity decreases as a consequence of domestication, mainly in transgene-introgressed (TI) individuals (π = 0.00020, 0.00001, 0.00016, 0, and 0, of wild, TI-wild, landraces, TI-landraces, and breeding lines, respectively). We identified homologous regions that differentiate wild from domesticated plants and indicate a relationship among the samples. A decrease in genetic diversity associated with transgene introgression in cotton was identified for the first time, and our outcomes are therefore relevant to both biosecurity and agrobiodiversity conservation.

Carriage and Transmission of mcr-1 in Salmonella Typhimurium and Its Monophasic 1,4,[5],12:i:- Variants from Diarrheal Outpatients: a 10-Year Genomic Epidemiology in Guangdong, Southern China

The banning of colistin as a feed additive for food-producing animals in mainland China in 2017 caused the decline in the prevalence of Escherichia coli-mobilized colistin resistance (mcr-1) in China. Salmonella Typhimurium and its monophasic 1,4,[5],12:i:- variants are also the main species associated with the spread of mcr-1; however, the evidence of the prevalence and transmission of mcr-1 among Salmonella is lacking. Herein, the 5,354 Salmonella isolates recovered from fecal samples of diarrheal patients in Guangdong, Southern China, from 2009 to 2019 were screened for colistin resistance and mcr-1, and mcr-1-positive isolates were characterized based on whole-genome sequencing (WGS) data. Relatively high prevalence rates of colistin resistance and mcr-1 (4.05%/4.50%) were identified, and more importantly, the prevalence trends of colistin-resistant and mcr-1-positive Salmonella isolates had a similar dynamic profile, i.e., both were first detected in 2012 and rapidly increased during 2013 to 2016, followed by a sharp decrease since 2017. WGS and phylogenetic analysis indicate that, whether before or after the ban, the persistence and cross-hospital transmission of mcr-1 are primarily determined by IncHI2 plasmids with similar backbones and sequence type 34 (ST34) Salmonella in specific clades that are associated with a high prevalence of IncHI2 plasmids and clinically important antimicrobial resistance genes, including bla_CTX-M-14-fosA3-oqxAB-floR genotypes. Our work reveals the difference in the prevalence rate of mcr-1 in clinical Salmonella before and after the Chinese colistin ban, whereas mcr-1 transmission was closely linked to multidrug-resistant IncHI2 plasmid and ST34 Salmonella across diverse hospitals over 10 years. Continued surveillance is required to explore the factors related to a sharp decrease in mcr-1 after the recent ban and determine whether the ban has affected the carriage of mcr-1 in Salmonella circulating in the health care system. IMPORTANCE Colistin is one of the last-line antibiotics for the clinical treatment of Enterobacteriaceae. However, the emergence of the mobilized colistin resistance (mcr-1) gene has spread throughout the entire human health system and largely threatens the usage of colistin in the clinical setting. In this study, we investigated the existence of mcr-1 in clinical Salmonella from a 10-year continuous surveillance and genomic study. Overall, the colistin resistance rate and mcr-1 carriage of Salmonella in tertiary hospitals in Guangdong (2009 to 2019) were relatively high and, importantly, rapidly increased from 2013 to 2016 and significantly decreased after the Chinese colistin withdrawal. However, before or after the ban, the MDR IncHI2 plasmid with a similar backbone and ST34 Salmonella were the main vectors involved in the spread of mcr-1. Interestingly, these Chinese mcr-1-carrying Salmonella obtain phylogenetically and phylogeographically distinct patterns compared with those from other continents and are frequently associated with clinically important ARGs including the extended-spectrum β-lactamases. Our data confirmed that the national stewardship intervention seems to be successful in blocking antibiotic resistance determinants and that continued surveillance of colistin resistance in clinical settings, farm animals, and related products is necessary.

Novel Multilocus Sequence Typing and Global Sequence Clustering Schemes for Characterizing the Population Diversity of Streptococcus mitis

Streptococcus mitis is a common oral commensal and an opportunistic pathogen that causes bacteremia and infective endocarditis; however, the species has received little attention compared to other pathogenic streptococcal species. Effective and easy-to-use molecular typing tools are essential for understanding bacterial population diversity and biology, but schemes specific for S. mitis are not currently available. We therefore developed a multilocus sequence typing (MLST) scheme and defined sequence clusters or lineages of S. mitis using a comprehensive global data set of 322 genomes (148 publicly available and 174 newly sequenced). We used internal 450-bp sequence fragments of seven housekeeping genes (accA, gki, hom, oppC, patB, rlmN, and tsf) to define the MLST scheme and derived the global S. mitis sequence clusters using the PopPUNK clustering algorithm. We identified an initial set of 259 sequence types (STs) and 258 global sequence clusters. The schemes showed high concordance (100%), capturing extensive S. mitis diversity with strains assigned to multiple unique STs and global sequence clusters. The tools also identified extensive within- and between-host S. mitis genetic diversity among isolates sampled from a cohort of healthy individuals, together with potential transmission events, supported by both phylogeny and pairwise single nucleotide polymorphism (SNP) distances. Our novel molecular typing and strain clustering schemes for S. mitis allow for the integration of new strain data, are electronically portable at the PubMLST database (https://pubmlst.org/smitis), and offer a standardized approach to understanding the population structure of S. mitis. These robust tools will enable new insights into the epidemiology of S. mitis colonization, disease and transmission.

Genomic characterization of oxacillin-susceptible mecA-positive Staphylococcus aureus ST59

Objective. To characterize the genomes of oxacillin-susceptible mecA-positive Staphylococcus aureus ST59 isolated in St. Petersburg. Materials and Methods. Nine oxacillin-susceptible mecA-positive of S. aureus isolates (OS-MRSA) of ST59 were included in the study. The isolates were obtained from children who showed no clinical signs of staphylococcal infections during nasal screening of S. aureus in St. Petersburg in 2018–2019. One isolate was obtained from an adult patient with skin and soft tissue infection (SSTI). The susceptibility to antibiotics and whole genome sequencing were performed. The analysis included 242 genomes of S. aureus ST59 from open access databases. Results. By employing the broth serial dilution and VITEK, the isolates’ phenotypic susceptibility to oxacillin was determined. The cefoxitin inhibition zones ranged from 17 to 22 mm. All isolates showed a penicillinclavulanate MIC ≤ 0.5 µg/mL. Isolates obtained from carriers belonged to the ST59-t1950-SCCmec Vb (seb+) genotype whereas the isolate obtained from SSTI belonged to the ST59-t437-SCCmec Vb (seb/ lukF/lukS+) genotype. Nucleotide position -33 (C/T) of mecA promoter and mutations in PBP2a (S225R + E246G) were present in all isolates. Based on phylogenetic analysis and Bayesian clustering the ST59 genomes were divided into four clusters and all Russian genomes belonged to the East Asian ST59 sublineage. The PVL toxin was present in the genomes of the first cluster of the East Asian ST59 sublineage. Pairwise comparisons of nucleotide substitutions among the genomes of Russian isolates showed a high similarity: median 13, interquartile range 8–18. All ST59 clusters were characterized by the presence of enterotoxin B, as well as mutations in PBP2a (S225R and E246G) and the promoter regions of the mecA gene (-7 G/A or -33 C/T). The genomes of the Russian isolates differed from the globally spread ST59 by specific mutations at the following loci (relative to the reference genome of S. aureus M013TW): lactose catabolism regulator RS03495 (N168D), ribosomal protein L28 (V47A), putative glyoxalase RS07825 (V42A), and the hypothetical protein RS13235 (K32E). Conclusions. Russian MRSA-ST59 isolates belong to the East Asian sublineage and are characterized by the presence of the enterotoxin B gene. Oxacillin susceptibility and borderline resistance to cefoxitin are specific characteristics of MRSA-ST59. OS-MRSA phenotypes have a risk of improper sensitivity testing leading to ineffective antibiotic treatment. Detection of mecA gene is the most accurate method for differentiating between MSSA and MRSA.

The anthropogenic effect of land use on population genetics of Malcus inconspicuus

Since the beginning of the Holocene era, human activities have seriously impacted animal habitats and vegetative environments. Species that are dependent on natural habitats or with narrow niches might be more severely affected by habitat changes. Malcus inconspicuus is distributed in subtropical China and highly dependent on the mountain environment. Our study investigated the role of the mountainous landscape in the historical evolution of M. inconspicuus and the impact of Holocene human activities on it. A phylogeographical approach was implemented with integrative datasets including double-digest restriction site-associated DNA (ddRAD), mitochondrial data, and distribution data. Three obvious clades and an east-west phylogeographical pattern were found in subtropical China. Mountainous landscape has "multifaceted" effects on the evolutionary history of M. inconspicuus, it has contributed to population differentiation, provided glacial refuges, and provided population expansion corridors during the postglacial period. The effective population size (Ne) of M. inconspicuus showed a sharp decline during the Holocene era, which revealed a significantly negative correlation with the development of cropland in a hilly area at the same time and space. It supported that the species which are highly dependent on natural habitats might undergo greater impact when the habitat was damaged by agricultural activities and we should pay more attention to them, especially in the land development of their distribution areas.

An advanced sequence clustering and designation workflow reveals the enzootic maintenance of a dominant West Nile virus subclade in Germany

West Nile virus (WNV) is the most widespread arthropod-borne (arbo) virus and the primary cause of arboviral encephalitis globally. Members of WNV species genetically diverged and are classified into different hierarchical groups below species rank. However, the demarcation criteria for allocating WNV sequences into these groups remain individual and inconsistent, and the use of names for different levels of the hierarchical levels is unstructured. In order to have an objective and comprehensible grouping of WNV sequences, we developed an advanced grouping workflow using the 'affinity propagation clustering' algorithm and newly included the 'agglomerative hierarchical clustering' algorithm for the allocation of WNV sequences into different groups below species rank. In addition, we propose to use a fixed set of terms for the hierarchical naming of WNV below species level and a clear decimal numbering system to label the determined groups. For validation, we applied the refined workflow to WNV sequences that have been previously grouped into various lineages, clades, and clusters in other studies. Although our workflow regrouped some WNV sequences, overall, it generally corresponds with previous groupings. We employed our novel approach to the sequences from the WNV circulation in Germany 2020, primarily from WNV-infected birds and horses. Besides two newly defined minor (sub)clusters comprising only three sequences each, Subcluster 2.5.3.4.3c was the predominant WNV sequence group detected in Germany from 2018 to 2020. This predominant subcluster was also associated with at least five human WNV infections in 2019-20. In summary, our analyses imply that the genetic diversity of the WNV population in Germany is shaped by enzootic maintenance of the dominant WNV subcluster accompanied by sporadic incursions of other rare clusters and subclusters. Moreover, we show that our refined approach for sequence grouping yields meaningful results. Although we primarily aimed at a more detailed WNV classification, the presented workflow can also be applied to the objective genotyping of other virus species.

Prevalence and Persistence of Ceftiofur-Resistant Escherichia coli in A Chicken Layer Breeding Program

We determined the longitudinal persistence of ceftiofur-resistant Escherichia coli from a chicken breeding farm in China. A total of 150 samples were collected from 5 breeding periods in a flock of layer hens, and the prevalence of ceftiofur-resistant E. coli fluctuated across the 5 chicken breeding stages: eggs, 3.33%; growing period, 100%; early laying period, 36.7%; peak laying period, 66.7% and late laying period, 80%. The most prevalent ceftiofur resistance genes were bla_CTX-M-55, bla_CMY and bla_NDM, and ST101 was the most prevalent and persistent sequence type across the breeding periods. Our results indicated that this breeder flock was heavily contaminated by ST101 ceftiofur-resistant E. coli and that its presence should be intensively monitored on chicken farms.

Phenotypic and Genotypic Characterization of Carbapenem-Resistant Enterobacteriaceae Recovered from a Single Hospital in China, 2013 to 2017

Objective

Carbapenem-resistant Enterobacteriaceae (CRE) have become an increasingly common cause of healthcare-related infections and present a serious challenge to clinical treatment. This study examined the phenotypic, genotypic characterization, clinical, and microbiological data of CRE in the Huizhou Municipal Central Hospital.

Methods

We conducted a phenotypic susceptibility evaluation and whole genome sequence analysis for 52 CRE strains isolated from 37 patients and 2 medical device-related samples during 2013-2017 to characterize risk factors, antimicrobial resistance profiles, dominant clones and hospital transmission.

Results

Long-term hospitalization, treatment time with antibiotics and use of invasive devices were linked to the risk of CRE infection. The carbapenem resistance genes (CRGs) we found included bla_NDM (82.7%), bla_IMP (19.2%) and bla_KPC (3.8%), Escherichia coli (44.2%) and Klebsiella pneumoniae (44.2%) were the dominant species we identified, and the type of CRG carried by isolates was highly correlated with species. The coexistence of CRGs with a variety of other antibiotic resistance genes leads to an increased prevalence of high resistance levels for CRE to β-lactams and other antibiotic classes such as aminoglycosides and fluoroquinolones. These isolates were sensitive only to colistin and tigecycline. In addition to this, we observed significantly genomic diversity of CRE isolates in this hospital. Importantly, we found that long-term transmission of multiple CRE clones had occurred at this hospital between various wards.

Conclusion

Evaluating and improving the current infection control strategies may be necessary, and reducing nosocomial transmission remains the primary control element for CRE infections in China.

Impact of international travel and diarrhea on gut microbiome and resistome dynamics

International travel contributes to the global spread of antimicrobial resistance. Travelers' diarrhea exacerbates the risk of acquiring multidrug-resistant organisms and can lead to persistent gastrointestinal disturbance post-travel. However, little is known about the impact of diarrhea on travelers' gut microbiomes, and the dynamics of these changes throughout travel. Here, we assembled a cohort of 159 international students visiting the Andean city of Cusco, Peru and applied next-generation sequencing techniques to 718 longitudinally-collected stool samples. We find that gut microbiome composition changed significantly throughout travel, but taxonomic diversity remained stable. However, diarrhea disrupted this stability and resulted in an increased abundance of antimicrobial resistance genes that can remain high for weeks. We also identified taxa differentially abundant between diarrheal and non-diarrheal samples, which were used to develop a classification model that distinguishes between these disease states. Additionally, we sequenced the genomes of 212 diarrheagenic Escherichia coli isolates and found those from travelers who experienced diarrhea encoded more antimicrobial resistance genes than those who did not. In this work, we find the gut microbiomes of international travelers' are resilient to dysbiosis; however, they are also susceptible to colonization by multidrug-resistant bacteria, a risk that is more pronounced in travelers with diarrhea.

Genomic Shift in Population Dynamics of mcr-1-positive Escherichia coli in Human Carriage

Emergence of the colistin resistance gene, mcr-1, has attracted worldwide attention. Despite the prevalence of mcr-1-positive Escherichia coli (MCRPEC) strains in human carriage showing a significant decrease between 2016 and 2019, genetic differences in MCRPEC strains remain largely unknown. We therefore conducted a comparative genomic study on MCRPEC strains from fecal samples of healthy human subjects in 2016 and 2019. We identified three major differences in MCRPEC strains between these two time points. First, the insertion sequence ISApl1 was often deleted and the percentage of mcr-1-carrying IncI2 plasmids was increased in MCRPEC strains in 2019. Second, the antibiotic resistance genes (ARGs), aac(3)-IVa and bla_CTX-M-1, emerged and coexisted with mcr-1 in 2019. Third, MCRPEC strains in 2019 contained more virulence genes, resulting in an increased proportion of extraintestinal pathogenic E. coli (ExPEC) strains (36.1%) in MCRPEC strains in 2019 compared to that in 2016 (10.5%), which implies that these strains could occupy intestinal ecological niches by competing with other commensal bacteria. Our results suggest that despite the significant reduction in the prevalence of MCRPEC strains in humans, mcr-1 is now associated with more stable genetic structures as well as the widespread IncI2 plasmid exhibits increased coexistence with other clinically important ARGs, and is increasingly associated with ExPEC strains, thus posing a potential public health threat.

Comparative genomics and pangenomics of vancomycin-resistant and susceptible Enterococcus faecium from Irish hospitals

Introduction. Enterococcus faecium has emerged as an important nosocomial pathogen, which is increasingly difficult to treat due to the genetic acquisition of vancomycin resistance. Ireland has a recalcitrant vancomycin-resistant bloodstream infection rate compared to other developed countries.

Hypothesis/Gap statement. Vancomycin resistance rates persist amongst E. faecium isolates from Irish hospitals. The evolutionary genomics governing these trends have not been fully elucidated.

Methodology. A set of 28 vancomycin-resistant isolates was sequenced to construct a dataset alongside 61 other publicly available Irish genomes. This dataset was extensively analysed using in silico methodologies (comparative genomics, pangenomics, phylogenetics, genotypics and comparative functional analyses) to uncover distinct evolutionary, coevolutionary and clinically relevant population trends.

Results. These results suggest that a stable (in terms of genome size, GC% and number of genes), yet genetically diverse population (in terms of gene content) of E. faecium persists in Ireland with acquired resistance arising via plasmid acquisition (vanA) or, to a lesser extent, chromosomal recombination (vanB). Population analysis revealed five clusters with one cluster partitioned into four clades which transcend isolation dates. Pangenomic and recombination analyses revealed an open (whole genome and chromosomal specific) pangenome illustrating a rampant evolutionary pattern. Comparative resistomics and virulomics uncovered distinct chromosomal and mobilomal propensity for multidrug resistance, widespread chromosomal point-mutation-mediated resistance and chromosomally harboured arsenals of virulence factors. Interestingly, a potential difference in biofilm formation strategies was highlighted by coevolutionary analysis, suggesting differential biofilm genotypes between vanA and vanB isolates.

Conclusions. These results highlight the evolutionary history of Irish E. faecium isolates and may provide insight into underlying infection dynamics in a clinical setting. Due to the apparent ease of vancomycin resistance acquisition over time, susceptible E. faecium should be concurrently reduced in Irish hospitals to mitigate potential resistant infections.

Molecular Epidemiology of Staphylococcus aureus in China Reveals the Key Gene Features Involved in Epidemic Transmission and Adaptive Evolution

Staphylococcus aureus is a Gram-positive pathogen that causes various infections in humans and domestic animals. In China, S. aureus is the most common Gram-positive pathogen that causes clinical infections. However, there are few comprehensive genome-based molecular epidemiology studies to investigate the genotypic background of the major S. aureus clones that are epidemic in China. Here, four S. aureus isolates that were recovered from hospital personnel were sequenced. In combination with whole-genome sequencing (WGS) data of 328 S. aureus strains as references, we performed a comprehensive molecular epidemiology study to reveal the molecular epidemic characterization of S. aureus that is epidemic in China. It was found that 332 S. aureus isolates were phylogenetically categorized into 4 major epidemic groups with different epidemiology phenotypes. Each group has exclusive features in virulence genotypic profiles, antimicrobial resistance genotypic profiles, core and pangenome features representing the differences involved in genetic features, evolutionary processes, and potential future evolutionary directions. Moreover, a comparative core genome analysis of 332 S. aureus isolates indicated several key genes that contributed to differences in molecular epidemic characterization and promoted the adaptive evolutionary process of each group. This study provides a comprehensive understanding of molecular epidemiological characteristics and adaptive evolutionary directions of major S. aureus clones that are epidemic in China. IMPORTANCE Staphylococcus aureus is an important Gram-positive pathogen that is epidemic worldwide and causes various infections in humans and domestic animals. However, there has been relatively little research on comprehensive molecular epidemiology in China. In this research, we reconstructed the phylogenetic relationship based on whole-genome data of strains almost all over China, screened for resistance and virulence genes, and took core and pan genome analysis to perform a comprehensive molecular epidemiology study of S. aureus that is epidemic in China. Our results highlight that there are 4 major epidemic groups with different epidemiology phenotypes after phylogenetic categorization with exclusive genetic features in virulence genotypic profiles, antimicrobial-resistance genotypic profiles, and core and pangenome features, and we found key gene features involved in epidemic transmission and adaptive evolution. Our findings are critical in describing molecular characteristic profiles of S. aureus infection, which could update existing preventive measures and take appropriate strategies.

Molecular Epidemiology of Plasmid-Mediated Types 1 and 3 Fimbriae Associated with Biofilm Formation in Multidrug Resistant Escherichia coli from Diseased Food Animals in Guangdong, China

Types 1 and 3 fimbriae in Enterobacteriaceae play versatile roles in bacterial physiology including attachment, invasion, cell motility as well as with biofilm formation and urinary tract infections. Herein, we investigated the prevalence and transmission of plasmid-mediated types 1 and 3 fimbriae from 1753 non-duplicate Enterobacteriaceae from diseased food Animals. We identified 123 (7.01%) strong biofilm producers and all was identified as E. coli. WGS analysis of 43 selected strong biofilm producers revealed that they harbored multiple ARGs, including ESBLs, PMQR and mcr-1. The gene clusters mrkABCDF and fimACDH encoding types 1 and 3 fimbriae, respectively, were identified among 43 (34.96%) and 7 (5.7%) of 123 strong biofilm isolates, respectively. These two operons were able to confer strong biofilm-forming ability to an E. coli weak-biofilm forming laboratory strain. Plasmid analysis revealed that mrk and fim operons were found to co-exist with ARGs and were primarily located on IncX1 and IncFII plasmids with similar backbones, respectively. mrkABCDF operons was present in all of 9457 Klebsiella pneumoniae using archived WGS data, and shared high homology to those on plasmids of 8 replicon types and chromosomes from 6 Enterobacteriaceae species from various origins and countries. In contrast, fimACDH operons was present in most of Enterobacter cloacae (62.15%), and shared high homology to those with only a small group of plasmids and Enterobacteriaceae species. This is the first comprehensive report of the prevalence, transmission and homology of plasmid-encoded type 1 and 3 fimbriae among the Enterobacteriaceae. Our findings indicated that plasmid-encoded mrkABCDF and fimACDH were major contributors to enhanced biofilm formation among E. coli and these two operons, in particular mrk could be as a potential anti-biofilm target. IMPORTANCE Biofilms allow bacteria to tolerate disinfectants and antimicrobials, as well as mammalian host defenses, and are therefore difficult to treat clinically. Most research concerning biofilm-related infections is typically focused on chromosomal biofilm-associated factors, including types 1 and 3 fimbriae of biofilm-forming Enterobacterium. However, the transmission and homology of the mobile types 1 and 3 fimbriae among Enterobacteriaceae is largely unknown. The findings revealed that the plasmid-encoded type 3 fimbriae encoded by mrkABCDF and type 1 fimbriae encoded by fimACDH were major contributors to enhancing biofilm formation among strong biofilm E. coli from diseased food producing animals. Additionally, mrk operon with high homology at an amino acid sequence was present both on plasmids of various replicon types and on chromosomes from diverse Enterobacteriaceae species from numerous origins and countries. These findings provide important information on the transmission of the mobile types 1 and 3 fimbriae among Enterobacteriaceae, indicating a potential antibiofilm target.

A unified global genotyping framework of dengue virus serotype-1 for a stratified coordinated surveillance strategy of dengue epidemics

BACKGROUND

Dengue is the fastest spreading arboviral disease, posing great challenges on global public health. A reproduceable and comparable global genotyping framework for contextualizing spatiotemporal epidemiological data of dengue virus (DENV) is essential for research studies and collaborative surveillance.

METHODS

Targeting DENV-1 spreading prominently in recent decades, by reconciling all qualified complete E gene sequences of 5003 DENV-1 strains with epidemiological information from 78 epidemic countries/areas ranging from 1944 to 2018, we established and characterized a unified global high-resolution genotyping framework using phylogenetics, population genetics, phylogeography, and phylodynamics.

RESULTS

The defined framework was discriminated with three hierarchical layers of genotype, subgenotype and clade with respective mean pairwise distances 2-6%, 0.8-2%, and ≤ 0.8%. The global epidemic patterns of DENV-1 showed strong geographic constraints representing stratified spatial-genetic epidemic pairs of Continent-Genotype, Region-Subgenotype and Nation-Clade, thereby identifying 12 epidemic regions which prospectively facilitates the region-based coordination. The increasing cross-transmission trends were also demonstrated. The traditional endemic countries such as Thailand, Vietnam and Indonesia displayed as persisting dominant source centers, while the emerging epidemic countries such as China, Australia, and the USA, where dengue outbreaks were frequently triggered by importation, showed a growing trend of DENV-1 diffusion. The probably hidden epidemics were found especially in Africa and India. Then, our framework can be utilized in an accurate stratified coordinated surveillance based on the defined viral population compositions. Thereby it is prospectively valuable for further hampering the ongoing transition process of epidemic to endemic, addressing the issue of inadequate monitoring, and warning us to be concerned about the cross-national, cross-regional, and cross-continental diffusions of dengue, which can potentially trigger large epidemics.

CONCLUSIONS

The framework and its utilization in quantitatively assessing DENV-1 epidemics has laid a foundation and re-unveiled the urgency for establishing a stratified coordinated surveillance platform for blocking global spreading of dengue. This framework is also expected to bridge classical DENV-1 genotyping with genomic epidemiology and risk modeling. We will promote it to the public and update it periodically.

Mandrake: visualizing microbial population structure by embedding millions of genomes into a low-dimensional representation

In less than a decade, population genomics of microbes has progressed from the effort of sequencing dozens of strains to thousands, or even tens of thousands of strains in a single study. There are now hundreds of thousands of genomes available even for a single bacterial species, and the number of genomes is expected to continue to increase at an accelerated pace given the advances in sequencing technology and widespread genomic surveillance initiatives. This explosion of data calls for innovative methods to enable rapid exploration of the structure of a population based on different data modalities, such as multiple sequence alignments, assemblies and estimates of gene content across different genomes. Here, we present Mandrake, an efficient implementation of a dimensional reduction method tailored for the needs of large-scale population genomics. Mandrake is capable of visualizing population structure from millions of whole genomes, and we illustrate its usefulness with several datasets representing major pathogens. Our method is freely available both as an analysis pipeline (https://github.com/johnlees/mandrake) and as a browser-based interactive application (https://gtonkinhill.github.io/mandrake-web/).

This article is part of a discussion meeting issue ‘Genomic population structures of microbial pathogens’.

A scalable analytical approach from bacterial genomes to epidemiology

Recent years have seen a remarkable increase in the practicality of sequencing whole genomes from large numbers of bacterial isolates. The availability of this data has huge potential to deliver new insights into the evolution and epidemiology of bacterial pathogens, but the scalability of the analytical methodology has been lagging behind that of the sequencing technology. Here we present a step-by-step approach for such large-scale genomic epidemiology analyses, from bacterial genomes to epidemiological interpretations. A central component of this approach is the dated phylogeny, which is a phylogenetic tree with branch lengths measured in units of time. The construction of dated phylogenies from bacterial genomic data needs to account for the disruptive effect of recombination on phylogenetic relationships, and we describe how this can be achieved. Dated phylogenies can then be used to perform fine-scale or large-scale epidemiological analyses, depending on the proportion of cases for which genomes are available. A key feature of this approach is computational scalability and in particular the ability to process hundreds or thousands of genomes within a matter of hours. This is a clear advantage of the step-by-step approach described here. We discuss other advantages and disadvantages of the approach, as well as potential improvements and avenues for future research. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Phylogenomic analysis of the bowfin (Amia calva) reveals unrecognized species diversity in a living fossil lineage

The Bowfin (Amia calva), as currently recognized, represents the sole living member of the family Amiidae, which dates back to approximately 150 Ma. Prior to 1896, 13 species of extant Bowfins had been described, but these were all placed into a single species with no rationale or analysis given. This situation has persisted until the present day, with little attention given to re-evaluation of those previously described nominal forms. Here, we present a phylogenomic analysis based on over 21,000 single nucleotide polymorphisms (SNPs) from 94 individuals that unambiguously demonstrates the presence of at least two independent evolutionary lineages within extant Amia populations that merit species-level standing, as well as the possibility of two more. These findings not only expand the recognizable species diversity in an iconic, ancient lineage, but also demonstrate the utility of such methods in addressing previously intractable questions of molecular systematics and phylogeography in slowly evolving groups of ancient fishes.

One clone to rule them all: Culture-independent genomics of Chlamydia psittaci from equine and avian hosts in Australia

Chlamydia psittaci is an avian pathogen with zoonotic potential. In Australia, C. psittaci has been well reported as a cause of reproductive loss in mares which subsequently have been the source of infection and illness in some in-contact humans. To date, molecular typing studies describe the predominant and clonal C. psittaci sequence type (ST)24 strains in horse, psittacine, and human infections. We sought to assess the clonality between ST24 strains and the emergence of equine ST24 with a comprehensive genomics approach. We used culture-independent probe-based and metagenomic whole-genome sequencing to investigate 13 C. psittaci genomes from horses, psittacines, and a pigeon from Australia. Published genomes of 36 C. psittaci strains were also used to contextualise our Australian dataset and investigate lineage diversity. We utilised a single-nucleotide polymorphism (SNP) based clustering and multi-locus sequence typing (MLST) approach. C. psittaci has four major phylogenetic groups (PG1-4) based on core-genome SNP-based phylogeny. PG1 contained clonal global and Australian equine, psittacine, and human ST24 genomes, with a median pairwise SNP distance of 68 SNPs. PG2, PG3, and PG4 had greater genomic diversity, including diverse STs collected from birds, livestock, human, and horse hosts from Europe and North America and a racing pigeon from Australia. We show that the clustering of C. psittaci by MLST was congruent with SNP-based phylogeny. The monophyletic ST24 clade has four major sub-lineages. The genomes of 17 Australian human, equine, and psittacine strains collected between 2008 and 2021 formed the predominant ST24 sub-lineage 1 (emerged circa 1979). Despite a temporal distribution of 13 years, the genomes within sub-lineage 1 had a median pairwise SNP distance of 32 SNPs, suggesting a recent population expansion or potential cross-host transmission. However, two C. psittaci genomes collected in 2015 from Victorian parrots clustered into distinct ST24 sub-lineage 4 (emerged circa 1965) with ovine strain C19/98 from Germany. This work describes a comprehensive phylogenomic characterisation of ST24 and identifies a timeline of potential bird-to-equine spillover events.

Phylogeography and colonization pattern of subendemic round-leaved oxeye daisy from the Dinarides to the Carpathians

The Carpathians are an important biodiversity hotspot and a link between mountain ranges on the European continent. This study investigated the phylogeography of one the Carpathian subendemics, Leucanthemum rotundifolium, which is distributed throughout the range and in one isolated population outside it. Range-wide sampling was used to examine phylogeographic patterns by sequencing uniparentally inherited chloroplast markers that exemplify seed dispersal. Reconstruct Ancestral State in Phylogenies (RASP) software, Bayesian binary Markov Chain Monte Carlo (BBM) analysis, and ecological niche modeling based on concatenated results of five algorithms were used to infer migration routes and examine links with other species through phylogeny. The round-leaved oxeye daisy is an example of organisms that reached the Carpathians through a southern "Dacian" migration route, most probably through long-distance dispersal. Dating placed the events in the Pleistocene and supported migrations during cooler periods and stasis/isolation followed by separation in the interglacials. Haplotype diversification indicated that after L. rotundifolium reached the area around the Fagaras Mountains, several migration events occurred leading to colonization of the Southern Carpathians followed by migration to the Apuseni Mountains, the Eastern Carpathians, and finally the Western Carpathians. The results are consistent with previous phylogeographic studies in this region and indicate several novel patterns.

Comparative genomics of Staphylococcus capitis reveals species determinants

Staphylococcus capitis is primarily described as a human skin commensal but is now emergent as an opportunistic pathogen isolated from the bloodstream and prosthetic joint infections, and neonatal intensive care unit (NICU)-associated sepsis. We used comparative genomic analyses of S. capitis to provide new insights into commensal scalp isolates from varying skin states (healthy, dandruff lesional, and non-lesional), and to expand our current knowledge of the species populations (scalp isolates, n = 59; other skin isolates, n = 7; publicly available isolates, n = 120). A highly recombinogenic population structure was revealed, with genomes including the presence of a range of previously described staphylococcal virulence factors, cell wall-associated proteins, and two-component systems. Genomic differences between the two described S. capitis subspecies were explored, which revealed the determinants associated exclusively with each subspecies. The subspecies ureolyticus was distinguished from subspecies capitis based on the differences in antimicrobial resistance genes, β-lactam resistance genes, and β-class phenol soluble modulins and gene clusters linked to biofilm formation and survival on skin. This study will aid further research into the classification of S. capitis and virulence-linked phylogroups to monitor the spread and evolution of S. capitis.

Phylogenomic analysis of the genus Delftia reveals distinct major lineages with ecological specializations

Delftia is a diverse betaproteobacterial genus with many strains having agricultural and industrial relevance, including plant-growth promotion, bioremediation of hydrocarbon-contaminated soils, and heavy metal immobilization. Delftia spp. are broadly distributed in the environment, and have been isolated from plant hosts as well as healthy and diseased animal hosts, yet the genetic basis of this ecological versatility has not been characterized. Here, we present a phylogenomic comparison of published Delftia genomes and show that the genus is divided into two well-supported clades: one 'Delftia acidovorans' clade with isolates from soils and plant rhizospheres, and a second 'Delftia lacustris and Delftia tsuruhatensis' clade with isolates from humans and sludge. The pan-genome inferred from 61 Delftia genomes contained over 28 000 genes, of which only 884 were found in all genomes. Analysis of industrially relevant functions highlighted the ecological versatility of Delftia and supported their role as generalists.

Contrasting long-term dynamics of antimicrobial resistance and virulence plasmids in Salmonella Typhimurium from animals

Plasmids are mobile elements that can carry genes encoding traits of clinical concern, including antimicrobial resistance (AMR) and virulence. Population-level studies of Enterobacterales, including Escherichia coli, Shigella and Klebsiella, indicate that plasmids are important drivers of lineage expansions and dissemination of AMR genes. Salmonella Typhimurium is the second most common cause of salmonellosis in humans and livestock in the UK and Europe. The long-term dynamics of plasmids between S. Typhimurium were investigated using isolates collected through national surveillance of animals in England and Wales over a 25-year period. The population structure of S. Typhimurium and its virulence plasmid (where present) were inferred through phylogenetic analyses using whole-genome sequence data for 496 isolates. Antimicrobial resistance genes and plasmid markers were detected in silico. Phenotypic plasmid characterization, using the Kado and Liu method, was used to confirm the number and size of plasmids. The differences in AMR and plasmids between clades were striking, with livestock clades more likely to carry one or more AMR plasmid and be multi-drug-resistant compared to clades associated with wildlife and companion animals. Multiple small non-AMR plasmids were distributed across clades. However, all hybrid AMR–virulence plasmids and most AMR plasmids were highly clade-associated and persisted over decades, with minimal evidence of horizontal transfer between clades. This contrasts with the role of plasmids in the short-term dissemination of AMR between diverse strains in other Enterobacterales in high-antimicrobial-use settings, with implications for predicting plasmid dissemination amongst S. Typhimurium.