Genome-scale rates of evolutionary change in bacteria

Multidrug-resistant Shigella flexneri outbreak affecting humans and non-human primates in New Mexico, USA

Shigellosis is a gastrointestinal infection caused by species of Shigella. A large outbreak of Shigella flexneri serotype 2a occurred in Albuquerque, New Mexico between May 2021 and November 2023 that involved humans and non-human primates (NHP) from a local zoo. We analyzed the genomes of 202 New Mexican isolates as well as 15 closely related isolates from other states, and four from NHP. The outbreak was initially detected within men who have sex with men but then predominantly affected people experiencing homelessness. Nearly 70% of cases were hospitalized and there was one human death. The outbreak extended into Albuquerque's BioPark Zoo, causing high morbidity and six deaths in NHPs. All isolates were multidrug-resistant, including towards fluoroquinolones, a first line treatment option which led to treatment failures in human and NHP populations. We show the circulation of the same S. flexneri strain in humans and NHPs, causing fatalities in both populations. This study demonstrates the threat of antimicrobial resistant organisms to vulnerable human and NHP populations and emphasizes the value of genomic surveillance within a One Health framework.

How does date-rounding affect phylodynamic inference for public health?

Phylodynamic analyses infer epidemiological parameters from pathogen genome sequences for enhanced genomic surveillance in public health. Pathogen genome sequences and their associated sampling dates are the essential data in every analysis. However, sampling dates are usually associated with hospitalisation or testing and can sometimes be used to identify individual patients, posing a threat to patient confidentiality. To lower this risk, sampling dates are often given with reduced date-resolution to the month or year, which can potentially bias inference. Here, we introduce a practical guideline on when date-rounding biases the inference of epidemiologically important parameters across a diverse range of empirical and simulated datasets. We show that the direction of bias varies for different parameters, datasets, and tree priors, while compounding with lower date-resolution and higher substitution rates. We also find that bias decreases for datasets with longer sampling intervals, implying that our guideline is most applicable to emerging datasets. We conclude by discussing future solutions that prioritise patient confidentiality and propose a method for safer sharing of sampling dates that translates them them uniformly by a random number.

Integrated population clustering and genomic epidemiology with PopPIPE

Genetic distances between bacterial DNA sequences can be used to cluster populations into closely related subpopulations and as an additional source of information when detecting possible transmission events. Due to their variable gene content and order, reference-free methods offer more sensitive detection of genetic differences, especially among closely related samples found in outbreaks. However, across longer genetic distances, frequent recombination can make calculation and interpretation of these differences more challenging, requiring significant bioinformatic expertise and manual intervention during the analysis process. Here, we present a Population analysis PIPEline (PopPIPE) which combines rapid reference-free genome analysis methods to analyse bacterial genomes across these two scales, splitting whole populations into subclusters and detecting plausible transmission events within closely related clusters. We use k-mer sketching to split populations into strains, followed by split k-mer analysis and recombination removal to create alignments and subclusters within these strains. We first show that this approach creates high-quality subclusters on a population-wide dataset of Streptococcus pneumoniae. When applied to nosocomial vancomycin-resistant Enterococcus faecium samples, PopPIPE finds transmission clusters that are more epidemiologically plausible than core genome or multilocus sequence typing (MLST) approaches. Our pipeline is rapid and reproducible, creates interactive visualizations and can easily be reconfigured and re-run on new datasets. Therefore, PopPIPE provides a user-friendly pipeline for analyses spanning species-wide clustering to outbreak investigations.

Ecology, global diversity and evolutionary mechanisms in the Mycobacterium tuberculosis complex

With the COVID-19 pandemic receding, tuberculosis (TB) is again the number one cause of human death to a single infectious agent. TB is caused by bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Recent advances in genome sequencing have provided new insights into the ecology and evolution of the MTBC. This includes the discovery of new phylogenetic lineages within the MTBC, a deeper understanding of the host tropism among the various animal-adapted lineages, enhanced knowledge on the evolutionary dynamics of antimicrobial resistance and transmission, as well as a better grasp of the within-host MTBC diversity. Moreover, advances in long-read sequencing are increasingly highlighting the relevance of structural genomic variation in the MTBC. These findings not only shed new light on the biology and epidemiology of TB, but also give rise to new questions and research avenues. The purpose of this Review is to summarize these new insights and discuss their implications for global TB control.

Evolution of pathogenic Escherichia coli harboring the transmissible locus of stress tolerance: from food sources to clinical environments

Escherichia coli (E. coli) carrying the transmissible locus of stress tolerance (tLST) are able to overcome numerous environmental challenges. In our in-silico study, we aimed to characterize tLST in terms of its variants in 793 genomes of E. coli from Brazil originating from food, environmental and clinical (animal and human) sources, and to perform a temporal analysis in order to identify the historical moment of its emergence. We also analyzed the presence of two Yersinia high pathogenicity island (HPI) variants in E. coli genomes, describing other genes and accessory for resistance, persistence, mobile elements (plasmids) and sequence types. The prevalence of the tLST was 10% in E. coli from Brazil, predominantly observed in milk-originating genomes, within the prevalent tLSTCP010237 variant. In E. coli from other sources (clinical/environmental), only part of the tLST was present. Remarkably, our temporal analysis pinpointed the emergence of tLST back to around 1914, coinciding with major societal events. Regarding virulence genes, we found a prevalence of 38.5% for HPI of Y. pestis across genomes from all sources. Our global analysis also showed a high diversity of other virulence genes for milk E. coli (+ 100 genes). These genomes also stood out from the overall metadata for presenting a greater variety of resistance genes to other stresses, such as metals, biocides and acids, as well as persistence genes (biofilm formation). This study demonstrated the historical background of E. coli with tLST genes dating back more than 100 years, and the acquisition of a wide range of virulence and resistance genes that allow it to circulate in different environments: from food to clinic or from clinic to food, making this bacterium a pathogen that requires rigorous surveillance and strategic interventions to mitigate potential risks.

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

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

Persistent, Private, and Mobile Genes: A Model for Gene Dynamics in Evolving Pangenomes

The pangenome of a species is the set of all genes carried by at least one member of the species. In bacteria, pangenomes can be much larger than the set of genes carried by a single organism. Many questions remain unanswered regarding the evolutionary forces shaping the patterns of the presence/absence of genes in pangenomes of a given species. We introduce a new model for bacterial pangenome evolution along a species phylogeny that explicitly describes the timing of appearance of each gene in the species and accounts for three generic types of gene evolutionary dynamics: persistent genes that are present in the ancestral genome, private genes that are specific to a given clade, and mobile genes that are imported once into the gene pool and then undergo frequent horizontal gene transfers. We call this model the Persistent-Private-Mobile (PPM) model. We develop an algorithm fitting the PPM model and apply it to a dataset of 902 Salmonella enterica genomes. We show that the best fitting model is able to reproduce the global pattern of some multivariate statistics like the gene frequency spectrum and the parsimony vs. frequency plot. Moreover, the gene classification induced by the PPM model allows us to study the position of accessory genes on the chromosome depending on their category, as well as the gene functions that are most present in each category. This work paves the way for a mechanistic understanding of pangenome evolution, and the PPM model developed here could be used for dynamics-aware gene classification.

Proteomic evidence for aerobic methane production in groundwater by methylotrophic Methylotenera

Members of Methylotenera are signature denitrifiers and methylotrophs commonly found together with methanotrophic bacteria in lakes and freshwater sediments. Here, we show that three distinct Methylotenera ecotypes were abundant in methane-rich groundwaters recharged during the Pleistocene. Just like in surface water biomes, groundwater Methylotenera often co-occurred with methane-oxidizing bacteria, even though they were generally unable to denitrify. One abundant Methylotenera ecotype expressed a pathway for aerobic methane production from methylphosphonate. This phosphate-acquisition strategy was recently found to contribute to methane production in the oligotrophic, oxic upper ocean. Gene organization, phylogeny, and 3D protein structure of the key enzyme, carbon-phosphorus lyase subunit PhnJ, were consistent with a role in phosphate uptake. We conclude that phosphate may be a limiting nutrient in productive, methane-rich aquifers, and that methylphosphonate degradation might contribute to groundwater methane production.

Research advancements on the diversity and host interaction of gut microbiota in chickens

The maintenance of host health and immune function is heavily dependent on the gut microbiota. However, the precise contribution of individual microbial taxa to regulating the overall functionality of the gut microbiome remains inadequately investigated. Chickens are commonly used as models for studying poultry gut microbiota, with high-throughput 16S rRNA sequencing has emerged as a valuable tool for assessing both its composition and functionality. The interactions between the gut's microbial community and its host significantly influence health outcomes, disease susceptibility, and various mechanisms affecting gastrointestinal function. Despite substantial research efforts, the dynamic nature of this microbial ecosystem has led to inconsistencies in findings related to chicken gut microbiota, which is largely attributed to variations in rearing conditions. Consequently, the interaction between the chickens' gut microflora and its host remains inadequately explored. This review highlights recent advances in understanding these relationships, with a specific focus on microbial composition, diversity, functional mechanisms, and their potential implications for improving poultry production.

Assessing the effect of model specification and prior sensitivity on Bayesian tests of temporal signal

Our understanding of the evolution of many microbes has been revolutionised by the molecular clock, a statistical tool to infer evolutionary rates and timescales from analyses of biomolecular sequences. In all molecular clock models, evolutionary rates and times are jointly unidentifiable and 'calibration' information must therefore be used. For many organisms, sequences sampled at different time points can be employed for such calibration. Before attempting to do so, it is recommended to verify that the data carry sufficient information for molecular dating, a practice referred to as evaluation of temporal signal. Recently, a fully Bayesian approach, BETS (Bayesian Evaluation of Temporal Signal), was proposed to overcome known limitations of other commonly used techniques such as root-to-tip regression or date randomisation tests. BETS requires the specification of a full Bayesian phylogenetic model, posing several considerations for untangling the impact of model choice on the detection of temporal signal. Here, we aimed to (i) explore the effect of molecular clock model and tree prior specification on the results of BETS and (ii) provide guidelines for improving our confidence in molecular clock estimates. Using microbial molecular sequence data sets and simulation experiments, we assess the impact of the tree prior and its hyperparameters on the accuracy of temporal signal detection. In particular, highly informative priors that are inconsistent with the data can result in the incorrect detection of temporal signal. In consequence, we recommend: (i) using prior predictive simulations to determine whether the prior generates a reasonable expectation of parameters of interest, such as the evolutionary rate and age of the root node, (ii) conducting prior sensitivity analyses to assess the robustness of the posterior to the choice of prior, and (iii) selecting a molecular clock model that reasonably describes the evolutionary process.

Rifaximin prophylaxis causes resistance to the last-resort antibiotic daptomycin

Multidrug-resistant bacterial pathogens like vancomycin-resistant Enterococcus faecium (VREfm) are a critical threat to human health1. Daptomycin is a last-resort antibiotic for VREfm infections with a novel mode of action2, but for which resistance has been widely reported but is unexplained. Here we show that rifaximin, an unrelated antibiotic used prophylactically to prevent hepatic encephalopathy in patients with liver disease3, causes cross-resistance to daptomycin in VREfm. Amino acid changes arising within the bacterial RNA polymerase in response to rifaximin exposure cause upregulation of a previously uncharacterized operon (prdRAB) that leads to cell membrane remodelling and cross-resistance to daptomycin through reduced binding of the antibiotic. VREfm with these mutations are spread globally, making this a major mechanism of resistance. Rifaximin has been considered 'low risk' for the development of antibiotic resistance. Our study shows that this assumption is flawed and that widespread rifaximin use, particularly in patients with liver cirrhosis, may be compromising the clinical use of daptomycin, a major last-resort intervention for multidrug-resistant pathogens. These findings demonstrate how unanticipated antibiotic cross-resistance can undermine global strategies designed to preserve the clinical use of critical antibiotics.

Seamless, rapid, and accurate analyses of outbreak genomic data using split k-mer analysis

Sequence variation observed in populations of pathogens can be used for important public health and evolutionary genomic analyses, especially outbreak analysis and transmission reconstruction. Identifying this variation is typically achieved by aligning sequence reads to a reference genome, but this approach is susceptible to reference biases and requires careful filtering of called genotypes. There is a need for tools that can process this growing volume of bacterial genome data, providing rapid results, but that remain simple so they can be used without highly trained bioinformaticians, expensive data analysis, and long-term storage and processing of large files. Here we describe split k-mer analysis (SKA2), a method that supports both reference-free and reference-based mapping to quickly and accurately genotype populations of bacteria using sequencing reads or genome assemblies. SKA2 is highly accurate for closely related samples, and in outbreak simulations, we show superior variant recall compared with reference-based methods, with no false positives. SKA2 can also accurately map variants to a reference and be used with recombination detection methods to rapidly reconstruct vertical evolutionary history. SKA2 is many times faster than comparable methods and can be used to add new genomes to an existing call set, allowing sequential use without the need to reanalyze entire collections. With an inherent absence of reference bias, high accuracy, and a robust implementation, SKA2 has the potential to become the tool of choice for genotyping bacteria. SKA2 is implemented in Rust and is freely available as open-source software.

Genetic heterogeneity in the Salmonella Typhi Vi capsule locus: a population genomic study from Fiji

Typhoid fever is endemic in many parts of the world and remains a major public health concern in tropical and sub-tropical developing nations, including Fiji. To address high rates of typhoid fever, the Northern Division of Fiji implemented a mass vaccination with typhoid conjugate vaccine (Vi-polysaccharide conjugated to tetanus toxoid) as a public health control measure in 2023. In this study we define the genomic epidemiology of Salmonella Typhi in the Northern Division prior to island-wide vaccination, sequencing 85% (n=419) of the total cases from the Northern and Central Divisions of Fiji that occurred in the period 2017-2019. We found elevated rates of nucleotide polymorphisms in the tviD and tviE genes (responsible for Vi-polysaccharide synthesis) relative to core genome levels within the Fiji endemic S. Typhi genotype 4.2. Expansion of these findings within a globally representative database of 12 382 S. Typhi (86 genotyphi clusters) showed evidence of convergent evolution of the same tviE mutations across the S. Typhi population, indicating that tvi selection has occurred both independently and globally. The functional impact of tvi mutations on the Vi-capsular structure and other phenotypic characteristics are not fully elucidated, yet commonly occurring tviE polymorphisms localize adjacent to predicted active site residues when overlayed against the predicted TviE protein structure. Given the central role of the Vi-polysaccharide in S. Typhi biology and vaccination, further integrated epidemiological, genomic and phenotypic surveillance is required to determine the spread and functional implications of these mutations.

Multidrug-resistant Shigella flexneri outbreak associated with a high-mortality spillover event into nonhuman primates

Shigellosis is a gastrointestinal infection caused by species of Shigella. A large outbreak of Shigella flexneri serotype 2a occurred in Albuquerque, New Mexico (NM) between May 2021 and November 2023 that involved humans and nonhuman primates (NHP) from a local zoo. We analyzed the genomes of 202 New Mexico isolates as well as 15 closely related isolates from other states, and four from NHP. The outbreak was initially detected within men who have sex with men (MSM) but then predominantly affected people experiencing homelessness (PEH). Nearly 70% of cases were hospitalized and there was one human death. The outbreak extended into Albuquerque's BioPark Zoo, causing high morbidity and six deaths in NHPs. The NHP isolates were identical to those in the human outbreak. All isolates were multidrug-resistant, including towards fluoroquinolones, a first line treatment option which led to treatment failures in human and NHP populations. We demonstrate the transmission of this S. flexneri strain between humans and NHPs, causing fatalities in both populations. This study demonstrates the threat of antimicrobial resistant organisms to vulnerable human and primate populations and emphasizes the value of vigilant genomic surveillance within a One Health framework.

Factors Associated with Inappropriate Use of Antibiotics Among Animal Health Professionals in Selected Districts of Rwanda, 2021

BACKGROUND

Antibiotic resistance is a global health concern. Humans can acquire antibiotic resistance through human-to-human transmission, from the environment, via the food chain, and through the contact with animals. The National Action Plan on antimicrobial resistance 2020-2024 highlights the prudent use of antibiotics in veterinary activities as the key element in keeping antibiotics effective. We determined the factors associated with misuse of antibiotics among animal health professionals in Rwanda.

METHODS

This was a cross-sectional study that enrolled animal health field professionals from five districts, where stratified random sampling was used to select one district by each province of Rwanda. Structured questions were used during face-to-face interviews. The misuse of antibiotics was defined as the use of antibiotics for reasons other than treatment, the non-completion of required courses, or the use of a high dose (i.e., an overdose) of antibiotics. We collected socio-demographic data of respondents, as well as elementary knowledge and perceptions on veterinary antibiotics and antibiotic resistance. A backward stepwise logistic regression model was used to identify the factors that were predictive of the inappropriate use of antibiotics.

RESULTS

There were 256 respondents to the survey. Of those, 198 were male and 58 were female. Almost three quarters of respondents (n = 174/256; 68%) reported the misuse of antibiotics at least once in the previous 12 months. The final logistic regression analysis identified the following factors to be predictive of antibiotics misuse: aged ≤ 24 years (aOR 0.92; 95% CI [0.88, 0.96]; p < 0.001); low trust in veterinary antibiotics available in the local market (aOR 8.45; 95% CI [4.18, 17.07]; p < 0.01), insufficient knowledge about basic understanding of antibiotics and antibiotic resistance (aOR 2.78; 95% CI [1.38, 5.58], p < 0.01) and not acquiring any continuing education (aOR 1.97; 95% CI [1.02, 4.19]; p = 0.04).

CONCLUSIONS

This study identified inadequate perceptions of proper antibiotic use among animal health professionals. There is a need for continuous education on appropriate antibiotic use among animal health professionals to lessen the negative impact of antibiotic resistance on public health security.

Population genomic evidence of structured and connected Plasmodium vivax populations under host selection in Latin America

Pathogen genomic epidemiology has the potential to provide a deep understanding of population dynamics, facilitating strategic planning of interventions, monitoring their impact, and enabling timely responses, and thereby supporting control and elimination efforts of parasitic tropical diseases. Plasmodium vivax, responsible for most malaria cases outside Africa, shows high genetic diversity at the population level, driven by factors like sub-patent infections, a hidden reservoir of hypnozoites, and early transmission to mosquitoes. While Latin America has made significant progress in controlling Plasmodium falciparum, it faces challenges with residual P. vivax. To characterize genetic diversity and population structure and dynamics, we have analyzed the largest collection of P. vivax genomes to date, including 1474 high-quality genomes from 31 countries across Asia, Africa, Oceania, and America. While P. vivax shows high genetic diversity globally, Latin American isolates form a distinctive population, which is further divided into sub-populations and occasional clonal pockets. Genetic diversity within the continent was associated with the intensity of transmission. Population differentiation exists between Central America and the North Coast of South America, vs. the Amazon Basin, with significant gene flow within the Amazon Basin, but limited connectivity between the Northwest Coast and the Amazon Basin. Shared genomic regions in these parasite populations indicate adaptive evolution, particularly in genes related to DNA replication, RNA processing, invasion, and motility - crucial for the parasite's survival in diverse environments. Understanding these population-level adaptations is crucial for effective control efforts, offering insights into potential mechanisms behind drug resistance, immune evasion, and transmission dynamics.

Geographic population structure and distinct intra-population dynamics of globally abundant freshwater bacteria

Implications of geographic separation and temporal dynamics on the evolution of free-living bacterial species are widely unclear. However, the vast amount of metagenome sequencing data generated during the last decades from various habitats around the world provides an unprecedented opportunity for such investigations. Here, we exploited publicly available and new freshwater metagenomes in combination with the genomes of abundant freshwater bacteria to reveal geographic and temporal population structure. We focused on species that were detected across broad geographic ranges at high enough sequence coverage for meaningful population genomic analyses, associated with the predominant freshwater taxa acI, LD12, Polynucleobacter, and Candidatus Methylopumilus. Despite the broad geographic ranges, each species appeared as a sequence-discrete cluster, in contrast to abundant marine taxa, for which continuous diversity structures were reported on a global scale. Population differentiation increased significantly with spatial distance in all species, but notable dispersal barriers (e.g. oceanic) were not apparent. Yet, the different species showed contrasting rates of geographic divergence and strikingly different intra-population dynamics in time series within individual habitats. The change in an LD12 population over 7 years was minor (FST = 0.04) compared to differentiation between lakes, whereas a Polynucleobacter population displayed strong changes within merely 2 months (FST up to 0.54), similar in scale to differentiation between populations separated by thousands of kilometers. The slowly and steadily evolving LD12 population showed high strain diversity, whereas the dynamic Polynucleobacter population exhibited alternating clonal expansions of mostly two strains only. Based on the contrasting population structures, we propose distinct models of speciation.

Comparative Genomics of Clostridioides difficile

Clostridioides difficile, a Gram-positive spore-forming anaerobic bacterium, has rapidly emerged as the leading cause of nosocomial diarrhoea in hospitals. The availability of large numbers of genome sequences, mainly due to the use of next-generation sequencing methods, has undoubtedly shown their immense advantages in the determination of C. difficile population structure. The implementation of fine-scale comparative genomic approaches has paved the way for global transmission and recurrence studies, as well as more targeted studies, such as the PaLoc or CRISPR/Cas systems. In this chapter, we provide an overview of recent and significant findings on C. difficile using comparative genomic studies with implications for epidemiology, infection control and understanding of the evolution of C. difficile.

Infectious disease in the Pleistocene: Old friends or old foes?

The impact of endemic and epidemic disease on humans has traditionally been seen as a comparatively recent historical phenomenon associated with the Neolithisation of human groups, an increase in population size led by sedentarism, and increasing contact with domesticated animals as well as species occupying opportunistic symbiotic and ectosymbiotic relationships with humans. The orthodox approach is that Neolithisation created the conditions for increasing population size able to support a reservoir of infectious disease sufficient to act as selective pressure. This orthodoxy is the result of an overly simplistic reliance on skeletal data assuming that no skeletal lesions equated to a healthy individual, underpinned by the assumption that hunter-gatherer groups were inherently healthy while agricultural groups acted as infectious disease reservoirs. The work of van Blerkom, Am. J. Phys. Anthropol., vol. suppl 37 (2003), Wolfe et al., Nature, vol. 447 (2007) and Houldcroft and Underdown, Am. J. Phys. Anthropol., vol. 160, (2016) has changed this landscape by arguing that humans and pathogens have long been fellow travelers. The package of infectious diseases experienced by our ancient ancestors may not be as dissimilar to modern infectious diseases as was once believed. The importance of DNA, from ancient and modern sources, to the study of the antiquity of infectious disease, and its role as a selective pressure cannot be overstated. Here we consider evidence of ancient epidemic and endemic infectious diseases with inferences from modern and ancient human and hominin DNA, and from circulating and extinct pathogen genomes. We argue that the pandemics of the past are a vital tool to unlock the weapons needed to fight pandemics of the future.

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

IMPORTANCE

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

Novel multidrug-resistant sublineages of Staphylococcus aureus clonal complex 22 discovered in India

Staphylococcus aureus is a major pathogen in India causing community and nosocomial infections, but little is known about its molecular epidemiology and mechanisms of resistance in hospital settings. Here, we use whole-genome sequencing (WGS) to characterize 478 S. aureus clinical isolates (393 methicillin-resistant Staphylococcus aureus (MRSA) and 85 methicilin-sensitive Staphylococcus aureus (MSSA) collected from 17 sentinel sites across India between 2014 and 2019. Sequencing results confirmed that sequence type 22 (ST22) (142 isolates, 29.7%), ST239 (74 isolates, 15.48%), and ST772 (67 isolates, 14%) were the most common clones. An in-depth analysis of 175 clonal complex (CC) 22 Indian isolates identified two novel ST22 MRSA lineages, both Panton-Valentine leukocidin+, both resistant to fluoroquinolones and aminoglycosides, and one harboring the the gene for toxic shock syndrome toxin 1 (tst). A temporal analysis of 1797 CC22 global isolates from 14 different studies showed that the two Indian ST22 lineages shared a common ancestor in 1984 (95% highest posterior density [HPD]: 1982-1986), as well as evidence of transmission to other parts of the world. Moreover, the study also gives a comprehensive view of ST2371, a sublineage of CC22, as a new emerging lineage in India and describes it in relationship with the other Indian ST22 isolates. In addition, the retrospective identification of a putative outbreak of multidrug-resistant (MDR) ST239 from a single hospital in Bangalore that persisted over a period of 3 years highlights the need for the implementation of routine surveillance and simple infection prevention and control measures to reduce these outbreaks. To our knowledge, this is the first WGS study that characterized CC22 in India and showed that the Indian clones are distinct from the EMRSA-15 clone. Thus, with the improved resolution afforded by WGS, this study substantially contributed to our understanding of the global population of MRSA. IMPORTANCE The study conducted in India between 2014 and 2019 presents novel insights into the prevalence of MRSA in the region. Previous studies have characterized two dominant clones of MRSA in India, ST772 and ST239, using whole-genome sequencing. However, this study is the first to describe the third dominant clone, ST22, using the same approach. The ST22 Indian isolates were analyzed in-depth, leading to the discovery of two new sublineages of hospital-acquired Staphylococcus aureus in India, both carrying antimicrobial resistance genes and mutations, which limit treatment options for patients. One of the newly characterized sublineages, second Indian cluster, carries the tsst-1 virulence gene, increasing the risk of severe infections. The geographic spread of the two novel lineages, both within India and internationally, could pose a global public health threat. The study also sheds light on ST2371 in India, a single-locus variant of ST22. The identification of a putative outbreak of MDR ST239 in a single hospital in Bangalore emphasizes the need for routine surveillance and simple infection prevention and control measures to reduce these outbreaks. Overall, this study significantly contributes to our understanding of the global population of MRSA, thanks to the improved resolution afforded by WGS.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

Global diversity and antimicrobial resistance of typhoid fever pathogens: Insights from a meta-analysis of 13,000 Salmonella Typhi genomes

BACKGROUND

The Global Typhoid Genomics Consortium was established to bring together the typhoid research community to aggregate and analyse Salmonella enterica serovar Typhi (Typhi) genomic data to inform public health action. This analysis, which marks 22 years since the publication of the first Typhi genome, represents the largest Typhi genome sequence collection to date (n=13,000).

METHODS

This is a meta-analysis of global genotype and antimicrobial resistance (AMR) determinants extracted from previously sequenced genome data and analysed using consistent methods implemented in open analysis platforms GenoTyphi and Pathogenwatch.

RESULTS

Compared with previous global snapshots, the data highlight that genotype 4.3.1 (H58) has not spread beyond Asia and Eastern/Southern Africa; in other regions, distinct genotypes dominate and have independently evolved AMR. Data gaps remain in many parts of the world, and we show the potential of travel-associated sequences to provide informal 'sentinel' surveillance for such locations. The data indicate that ciprofloxacin non-susceptibility (>1 resistance determinant) is widespread across geographies and genotypes, with high-level ciprofloxacin resistance (≥3 determinants) reaching 20% prevalence in South Asia. Extensively drug-resistant (XDR) typhoid has become dominant in Pakistan (70% in 2020) but has not yet become established elsewhere. Ceftriaxone resistance has emerged in eight non-XDR genotypes, including a ciprofloxacin-resistant lineage (4.3.1.2.1) in India. Azithromycin resistance mutations were detected at low prevalence in South Asia, including in two common ciprofloxacin-resistant genotypes.

CONCLUSIONS

The consortium's aim is to encourage continued data sharing and collaboration to monitor the emergence and global spread of AMR Typhi, and to inform decision-making around the introduction of typhoid conjugate vaccines (TCVs) and other prevention and control strategies.

FUNDING

No specific funding was awarded for this meta-analysis. Coordinators were supported by fellowships from the European Union (ZAD received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 845681), the Wellcome Trust (SB, Wellcome Trust Senior Fellowship), and the National Health and Medical Research Council (DJI is supported by an NHMRC Investigator Grant [GNT1195210]).

Genome-wide phylodynamic approach reveals the epidemic dynamics of the main Mycoplasma bovis subtype circulating in France

Mycoplasma bovis is a major aetiological agent of bovine respiratory disease worldwide. Genome-based analyses are increasingly being used to monitor the genetic diversity and global distribution of M. bovis, complementing existing subtyping schemes based on locus sequencing. However, these analyses have so far provided limited information on the spatiotemporal and population dynamics of circulating subtypes. Here we applied a genome-wide phylodynamic approach to explore the epidemic dynamics of 88 French M. bovis strains collected between 2000 and 2019 in France and belonging to the currently dominant polC subtype 2 (st2). A strong molecular clock signal detected in the genomic data enabled robust phylodynamic inferences, which estimated that the M. bovis st2 population in France is composed of two lineages that successively emerged from independent introductions of international strains. The first lineage appeared around 2000 and supplanted the previously established antimicrobial-susceptible polC subtype 1. The second lineage, which is likely more transmissible, progressively replaced the first M. bovis st2 lineage population from 2005 onward and became predominant after 2010. Analyses also showed a brief decline in this second M. bovis st2 lineage population in around 2011, possibly due to the challenge from the concurrent emergence of M. bovis polC subtype 3 in France. Finally, we identified non-synonymous mutations in genes associated with lineages, which raises prospects for identifying new surveillance molecular markers. A genome-wide phylodynamic approach provides valuable resources for monitoring the evolution and epidemic dynamics of circulating M. bovis subtypes, and may prove critical for developing more effective surveillance systems and disease control strategies.

Phylogenomics of Globally Spread Clonal Groups 14 and 15 of Klebsiella pneumoniae

Klebsiella pneumoniae sequence type 14 (ST14) and ST15 caused outbreaks of CTX-M-15 and/or carbapenemase producers worldwide, but their phylogeny and global dynamics remain unclear. We clarified the evolution of K. pneumoniae clonal group 14 (CG14) and CG15 by analyzing the capsular locus (KL), resistome, virulome, and plasmidome of public genomes (n = 481) and de novo sequences (n = 9) representing main sublineages circulating in Portugal. CG14 and CG15 evolved independently within 6 main subclades defined according to the KL and the accessory genome. The CG14 (n = 65) clade was structured in two large monophyletic subclades, CG14-I (KL2, 86%) and CG14-II (KL16, 14%), whose emergences were dated to 1932 and 1911, respectively. Genes encoding extended-spectrum β-lactamase (ESBL), AmpC, and/or carbapenemases were mostly observed in CG14-I (71% versus 22%). CG15 clade (n = 170) was segregated into subclades CG15-IA (KL19/KL106, 9%), CG15-IB (variable KL types, 6%), CG15-IIA (KL24, 43%) and CG15-IIB (KL112, 37%). Most CG15 genomes carried specific GyrA and ParC mutations and emerged from a common ancestor in 1989. CTX-M-15 was especially prevalent in CG15 (68% CG15 versus 38% CG14) and in CG15-IIB (92%). Plasmidome analysis revealed 27 predominant plasmid groups (PG), including particularly pervasive and recombinant F-type (n = 10), Col (n = 10), and new plasmid types. While blaCTX-M-15 was acquired multiple times by a high diversity of F-type mosaic plasmids, other antibiotic resistance genes (ARGs) were dispersed by IncL (blaOXA-48) or IncC (blaCMY/TEM-24) plasmids. We first demonstrate an independent evolutionary trajectory for CG15 and CG14 and how the acquisition of specific KL, quinolone-resistance determining region (QRDR) mutations (CG15), and ARGs in highly recombinant plasmids could have shaped the expansion and diversification of particular subclades (CG14-I and CG15-IIA/IIB). IMPORTANCE Klebsiella pneumoniae represents a major threat in the burden of antibiotic resistance (ABR). Available studies to explain the origin, the diversity, and the evolution of certain ABR K. pneumoniae populations have mainly been focused on a few clonal groups (CGs) using phylogenetic analysis of the core genome, the accessory genome being overlooked. Here, we provide unique insights into the phylogenetic evolution of CG14 and CG15, two poorly characterized CGs which have contributed to the global dissemination of genes responsible for resistance to first-line antibiotics such as β-lactams. Our results point out an independent evolution of these two CGs and highlight the existence of different subclades structured by the capsular type and the accessory genome. Moreover, the contribution of a turbulent flux of plasmids (especially multireplicon F type and Col) and adaptive traits (antibiotic resistance and metal tolerance genes) to the pangenome reflect the exposure and adaptation of K. pneumoniae under different selective pressures.

Widespread extinctions of co-diversified primate gut bacterial symbionts from humans

Humans and other primates harbour complex gut bacterial communities that influence health and disease, but the evolutionary histories of these symbioses remain unclear. This is partly due to limited information about the microbiota of ancestral primates. Here, using phylogenetic analyses of metagenome-assembled genomes (MAGs), we show that hundreds of gut bacterial clades diversified in parallel (that is, co-diversified) with primate species over millions of years, but that humans have experienced widespread losses of these ancestral symbionts. Analyses of 9,460 human and non-human primate MAGs, including newly generated MAGs from chimpanzees and bonobos, revealed significant co-diversification within ten gut bacterial phyla, including Firmicutes, Actinobacteriota and Bacteroidota. Strikingly, ~44% of the co-diversifying clades detected in African apes were absent from available metagenomic data from humans and ~54% were absent from industrialized human populations. In contrast, only ~3% of non-co-diversifying clades detected in African apes were absent from humans. Co-diversifying clades present in both humans and chimpanzees displayed consistent genomic signatures of natural selection between the two host species but differed in functional content from co-diversifying clades lost from humans, consistent with selection against certain functions. This study discovers host-species-specific bacterial symbionts that predate hominid diversification, many of which have undergone accelerated extinctions from human populations.

Estimating the fitness cost and benefit of antimicrobial resistance from pathogen genomic data

Increasing levels of antibiotic resistance in many bacterial pathogen populations are a major threat to public health. Resistance to an antibiotic provides a fitness benefit when the bacteria are exposed to this antibiotic, but resistance also often comes at a cost to the resistant pathogen relative to susceptible counterparts. We lack a good understanding of these benefits and costs of resistance for many bacterial pathogens and antibiotics, but estimating them could lead to better use of antibiotics in a way that reduces or prevents the spread of resistance. Here, we propose a new model for the joint epidemiology of susceptible and resistant variants, which includes explicit parameters for the cost and benefit of resistance. We show how Bayesian inference can be performed under this model using phylogenetic data from susceptible and resistant lineages and that by combining data from both we are able to disentangle and estimate the resistance cost and benefit parameters separately. We applied our inferential methodology to several simulated datasets to demonstrate good scalability and accuracy. We analysed a dataset of Neisseria gonorrhoeae genomes collected between 2000 and 2013 in the USA. We found that two unrelated lineages resistant to fluoroquinolones shared similar epidemic dynamics and resistance parameters. Fluoroquinolones were abandoned for the treatment of gonorrhoea due to increasing levels of resistance, but our results suggest that they could be used to treat a minority of around 10% of cases without causing resistance to grow again.

Rapid adaptations of Legionella pneumophila to the human host

Legionella pneumophila are host-adapted bacteria that infect and reproduce primarily in amoeboid protists. Using similar infection mechanisms, they infect human macrophages, and cause Legionnaires' disease, an atypical pneumonia, and the milder Pontiac fever. We hypothesized that, despite the similarities in infection mechanisms, the hosts are different enough that there exist high-selective value mutations that would dramatically increase the fitness of Legionella inside the human host. By comparing a large number of isolates from independent infections, we identified two genes, mutated in three unrelated patients, despite the short duration of the incubation period (2-14 days). One is a gene coding for an outer membrane protein (OMP) belonging to the OmpP1/FadL family. The other is a gene coding for an EAL-domain-containing protein involved in cyclic-di-GMP regulation, which in turn modulates flagellar activity. The clinical strain, carrying the mutated EAL-domain-containing homologue, grows faster in macrophages than the wild-type strain, and thus appears to be better adapted to the human host. As human-to-human transmission is very rare, fixation of these mutations into the population and spread into the environment is unlikely. Therefore, parallel evolution - here mutations in the same genes observed in independent human infections - could point to adaptations to the accidental human host. These results suggest that despite the ability of L. pneumophila to infect, replicate in and exit from macrophages, its human-specific adaptations are unlikely to be fixed in the population.

Genomic and Evolutionary Insights into Australian Toxigenic Vibrio cholerae O1 Strains

Vibrio cholerae O1 is the causative agent of cholera, a severe diarrheal disease which can cause death if left untreated. In this study, a collection of clinical and environmental V. cholerae serogroup O1 isolates from Australia (1977 to 1987) (from local cases and cases acquired through international travel) and publicly available international isolates were characterized for genotypic features (virulence genes, mobile genetic elements [MGEs], and antimicrobial resistance gene profiles). Whole-genome sequencing (WGS) was used to investigate and compare the genetic relatedness between the 44 Australian and nine travel-associated isolates and the 60 publicly available international V. cholerae sequences representing pre-seventh-pandemic (pre-7PET) isolates and different waves of 7PET isolates. In this study, 36 (81%) Australian clinical and aquatic isolates harbored the cholera toxin-producing genes located in the CTX bacteriophage region. All the Australian environmental and clinical isolates lacked the seventh-pandemic virulence-associated genomic islands (VSP-I and -II). In silico multilocus sequence typing (MLST) classified all nine internationally acquired isolates as sequence type 69 (ST69), 36 clinical and aquatic isolates as ST70, and eight isolates from Australia as ST71. Most of the nontoxigenic clinical and aquatic isolates of ST71 had diverse genetic variations compared to ST70 Australian strains. The antimicrobial resistance-associated genes gyrA, parC, and parE had no mutations in all the environmental and clinical isolates from Australia. The SXT genetic element and class 1 integron gene sequences were not detected in Australian strains. Moreover, in this study, a Bayesian evolutionary study suggests that two distinct lineages of ST71 (new set of strains) and ST70 strains were prevalent around similar times in Australia, in ~1973 and 1969. IMPORTANCE Australia has its own indigenous V. cholerae strains, both toxigenic and nontoxigenic, that are associated with disease. Exotic strains are also detected in Australian patients returning from overseas travel. The clinical and aquatic V. cholerae O1 toxin gene-positive isolates from Australia responsible for cases in 1977 to 1987 were linked to acquisition from Queensland waterways but until now had not been characterized genetically. It is important to determine the genetic relatedness of Australian strains to international strains to assist in understanding their origin. This is the first extensive study to provide sequences and genomic analysis focused on toxigenic O1 V. cholerae clinical and environmental strains from Australia and its possible evolutionary relationship with other publicly available pre-7PET and 7PET V. cholerae strains. It is important to understand the population genetics of Australian V. cholerae from a public health perspective to assist in devising control measures and management plans for reducing V. cholerae exposure in Australia, given previous Australian disease clusters.

Whole-genome-based characterization of Campylobacter jejuni from human patients with gastroenteritis collected over an 18 year period reveals increasing prevalence of antimicrobial resistance

Campylobacteriosis is the most common cause of acute gastrointestinal bacterial infection in Europe, with most infections linked to the consumption of contaminated food. While previous studies found an increasing rate of antimicrobial resistance (AMR) in Campylobacter spp. over the past decades, the investigation of additional clinical isolates is likely to provide novel insights into the population structure and mechanisms of virulence and drug resistance of this important human pathogen. Therefore, we combined whole-genome sequencing and antimicrobial-susceptibility testing of 340 randomly selected Campylobacter jejuni isolates from humans with gastroenteritis, collected in Switzerland over an 18 year period. In our collection, the most common multilocus sequence types (STs) were ST-257 (n=44), ST-21 (n=36) and ST-50 (n=35); the most common clonal complexes (CCs) were CC-21 (n=102), CC-257 (n=49) and CC-48 (n=33). High heterogeneity was observed among STs, with the most abundant STs recurring over the entire study period, while others were observed only sporadically. Source attribution based on ST assigned more than half of the strains to the 'generalist' category (n=188), 25  % as 'poultry specialist' (n=83), and only a few to 'ruminant specialist' (n=11) or 'wild bird' origin (n=9). The isolates displayed an increased frequency of AMR from 2003 to 2020, with the highest rates of resistance observed for ciprofloxacin and nalidixic acid (49.8 %), followed by tetracycline (36.9 %). Quinolone-resistant isolates carried chromosomal gyrA mutations T86I (99.4 %) and T86A (0.6 %), whereas tetracycline-resistant isolates carried tet(O) (79.8 %) or mosaic tetO/32/O (20.2 %) genes. A novel chromosomal cassette carrying several resistance genes, including aph(3')-III, satA and aad(6), and flanked by insertion sequence elements was detected in one isolate. Collectively, our data revealed an increasing prevalence of resistance to quinolones and tetracycline in C. jejuni isolates from Swiss patients over time, linked to clonal expansion of gyrA mutants and acquisition of the tet(O) gene. Investigation of source attribution suggests that infections are most likely related to isolates from poultry or generalist backgrounds. These findings are relevant to guide future infection prevention and control strategies.

Persistence and In Vivo Evolution of Vaginal Bacterial Strains over a Multiyear Time Period

It is not clear whether the bacterial strains that comprise our microbiota are mostly long-term colonizers or transient residents. Studies have demonstrated decades-long persistence of bacterial strains within the gut, but persistence at other body sites has yet to be determined. The vaginal microbiota (VMB) is often dominated by Lactobacillus, although it is also commonly comprised of a more diverse set of other facultative and obligate anaerobes. Longitudinal studies have demonstrated that these communities can be stable over several menstrual cycles or can fluctuate temporally in species composition. We sought to determine whether the bacterial strains that comprise the VMB were capable of persisting over longer time periods. We performed shotgun metagenomics on paired samples from 10 participants collected 1 and 2 years apart. The resulting sequences were de novo assembled and binned into high-quality metagenome assembled genomes. Persistent strains were identified based on the sequence similarity between the genomes present at the two time points and were found in the VMB of six of the participants, three of which had multiple persistent strains. The VMB of the remaining four participants was similar in species composition at the two time points but was comprised of different strains. For the persistent strains, we were able to identify the mutations that were fixed in the populations over the observed time period, giving insight into the evolution of these bacteria. These results indicate that bacterial strains can persist in the vagina for extended periods of time, providing an opportunity for them to evolve in the host microenvironment. IMPORTANCE The stability of strains within the vaginal microbiota is largely uncharacterized. Should these strains be capable of persisting for extended periods of time, they could evolve within their host in response to selective pressures exerted by the host or by other members of the community. Here, we present preliminary findings demonstrating that bacterial strains can persist in the vagina for at least 1 year. We further characterized in vivo evolution of the persistent strains. Several participants were also found to not have persistent strains, despite having a vaginal microbiota (VMB) with similar species composition at the two time points. Our observations motivate future studies that collect samples from more participants, at more time points, and over even longer periods of time. Understanding which strains persist, what factors drive their persistence, and what selective pressures they face will inform the development and delivery of rationally designed live biotherapeutics for the vagina.

Competitive Transmission of Carbapenem-Resistant Klebsiella pneumoniae in a Newly Opened Intensive Care Unit

We conducted a 6-month prospective study in a newly opened ICU for high-resolution tracking of carbapenem-resistant Klebsiella pneumoniae (CRKP) through environmental surveillance, patient screening, and genome sequencing. Among all ICU patients (n = 348) screened, 3.5% carried CRKP on admission and 16.3% acquired CRKP thereafter. CRKP was not detected in the environment until 10 weeks and was then isolated from 98 of 2,989 environmental samples (3.3%). The first CRKP isolate from rectal swabs (n = 37) and the first clinical isolate (n = 8) of each patient as well as the 98 isolates from environmental were subjected to whole-genome sequencing. The 143 CRKP isolates from patients and environment samples were assigned to four sequence types, with ST11 dominating (95.8%) and further divided into 14 clones, suggesting introduction of multiple clones. Subsequent CRKP transmission was complex and dynamic with 10 clones found in multiple patients and seven also detected in the environment. Two particular ST11 clones caused extensive (≥5 rooms) and persistent (≥10 weeks) environmental contamination. Both clones were associated with patients who carried CRKP throughout their prolonged ICU stay. Such "super-contaminators" are a priority for isolation and environmental surveillance. IMPORTANCE Carbapenem-resistant Klebsiella pneumoniae (CRKP) is a global challenge for human health. In health care settings, patients have frequent interactions with other patients and the environment, rendering challenges for untangling the introduction and transmission of CRKP. We conducted a prospective surveillance study in a newly opened ICU for high-resolution tracking of CRKP. Our study demonstrated the dynamic, complicated transmission of CRKP and has important findings that may help to curb its spread in health care settings. First, compliance with basic measures such as routine environment cleaning and postdischarge terminal cleaning is needed to minimize the environmental contamination-driven spread. Second, active screening could demonstrate the scale of the problem, and room transfer of patients with CRKP should be prohibited whenever possible. Third, the priority for single-room isolation should be given to patients with prolonged carriage of CRKP, especially in resource-limited settings. Good infection control practice lays a foundation for tackling multidrug-resistant organisms like CRKP.

Multiclonal human origin and global expansion of an endemic bacterial pathogen of livestock

Most new pathogens of humans and animals arise via switching events from distinct host species. However, our understanding of the evolutionary and ecological drivers of successful host adaptation, expansion, and dissemination are limited. Staphylococcus aureus is a major bacterial pathogen of humans and a leading cause of mastitis in dairy cows worldwide. Here we trace the evolutionary history of bovine S. aureus using a global dataset of 10,254 S. aureus genomes including 1,896 bovine isolates from 32 countries in 6 continents. We identified 7 major contemporary endemic clones of S. aureus causing bovine mastitis around the world and traced them back to 4 independent host-jump events from humans that occurred up to 2,500 y ago. Individual clones emerged and underwent clonal expansion from the mid-19th to late 20th century coinciding with the commercialization and industrialization of dairy farming, and older lineages have become globally distributed via established cattle trade links. Importantly, we identified lineage-dependent differences in the frequency of host transmission events between humans and cows in both directions revealing high risk clones threatening veterinary and human health. Finally, pangenome network analysis revealed that some bovine S. aureus lineages contained distinct sets of bovine-associated genes, consistent with multiple trajectories to host adaptation via gene acquisition. Taken together, we have dissected the evolutionary history of a major endemic pathogen of livestock providing a comprehensive temporal, geographic, and gene-level perspective of its remarkable success.

Characterization of Pantoea ananatis from rice planthoppers reveals a clade of rice-associated P. ananatis undergoing genome reduction

Pantoea ananatis is a bacterium that is found in many agronomic crops and agricultural pests. Here, we isolated a P. ananatis strain (Lstr) from the rice planthopper Laodelphax striatellus, a notorious pest that feeds on rice plant sap and transmits rice viruses, in order to examine its genome and biology. P. ananatis Lstr is an insect symbiont that is pathogenic to the host insect and appears to mostly inhabit the gut. Its pathogenicity thus raises the possibility of using the Lstr strain as a biological agent. To this end, we analysed the genome of the Lstr strain and compared it with the genomes of other Pantoea species. Our analysis of these genomes shows that P. ananatis can be divided into two mono-phylogenetic clades (clades one and two). The Lstr strain belongs to clade two and is grouped with P. ananatis strains that were isolated from rice or rice-associated samples. A comparative genomic analysis shows that clade two differs from clade one in many genomic characteristics including genome structures, mobile elements, and categories of coding proteins. The genomes of clade two P. ananatis are significantly smaller, have much fewer coding sequences but more pseudogenes than those of clade one, suggesting that clade two species are at the early stage of genome reduction. On the other hand, P. ananatis has a type VI secretion system that is highly variable but cannot be separated by clades. These results clarify our understanding of P. ananatis' phylogenetic diversity and provide clues to the interactions between P. ananatis, host insect, and plant that may lead to advances in rice protection and pest control.

Low-level colonization of methicillin-resistant Staphylococcus aureus in pigs is maintained by slowly evolving, closely related strains in Finnish pig farms

BACKGROUND

Over the past two decades, livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) has become widely prevalent in pig production in Europe. The carriage status of LA-MRSA is known to vary among individual pigs, but bacterial load in pigs has rarely been studied. We assessed the quantity of LA-MRSA in nasal and skin samples of pigs and investigated the genetic diversity of the strains together with sequenced strains from national surveillance and pathology samples from the Finnish Food Authority. On two farms with assumed MRSA-positive status, farm 1 and farm 2, 10 healthy pigs were sampled three times during 2 weeks from the nares and skin (study A). On farm 1, 54 additional pigs were sampled and from confirmed MRSA-positive animals, 10 were randomly selected and transported to a clean, controlled environment for further sampling (study B). From the samples taken on farms 1 and 2 and in the controlled environment, MRSA was isolated both by direct plating and enrichment on selective media. spa types, multilocus sequence types, staphylococcal cassette chromosome mec types, resistance and virulence genes were determined. Core genome multilocus sequence typing (cgMLST) analysis was performed, including the sequences deriving from the surveillance/pathology samples from the Finnish Food Authority.

RESULTS

All pigs on farm 1 carried LA-MRSA in the nares at all three time points and five pigs on farm 2 at one time point. Nasal quantity varied between 10 and 10³ CFU/swab and quantity on the skin between 10 and 10² CFU/swab. In the controlled environment, MRSA was detected in at least one of the nasal samples from each animal. spa type t034 was predominant. cgMLST showed one cluster with minimum allele differences between 0 and 11.

CONCLUSIONS

The study shows predominantly low-level carriage (< 10³ CFU/swab) of LA-MRSA on farms. In the controlled environment we observed a decline in nasal carriage but constant skin carriage. cgMLST showed that strains of spa type t034 are closely related at the national level.

Phylodynamic signatures in the emergence of community-associated MRSA

Community-associated, methicillin-resistant <i>Staphylococcus aureus</i> (MRSA) lineages have emerged in many geographically distinct regions around the world during the past 30 y. Here, we apply consistent phylodynamic methods across multiple community-associated MRSA lineages to describe and contrast their patterns of emergence and dissemination. We generated whole-genome sequencing data for the Australian sequence type (ST) ST93-MRSA-IV from remote communities in Far North Queensland and Papua New Guinea, and the Bengal Bay ST772-MRSA-V clone from metropolitan communities in Pakistan. Increases in the effective reproduction number (R_e) and sustained transmission (R_e &gt; 1) coincided with spread of progenitor methicillin-susceptible <i>S. aureus</i> (MSSA) in remote northern Australian populations, dissemination of the ST93-MRSA-IV genotype into population centers on the Australian East Coast, and subsequent importation into the highlands of Papua New Guinea and Far North Queensland. Applying the same phylodynamic methods to existing lineage datasets, we identified common signatures of epidemic growth in the emergence and epidemiological trajectory of community-associated <i>S. aureus</i> lineages from America, Asia, Australasia, and Europe. Surges in R_e were observed at the divergence of antibiotic-resistant strains, coinciding with their establishment in regional population centers. Epidemic growth was also observed among drug-resistant MSSA clades in Africa and northern Australia. Our data suggest that the emergence of community-associated MRSA in the late 20th century was driven by a combination of antibiotic-resistant genotypes and host epidemiology, leading to abrupt changes in lineage-wide transmission dynamics and sustained transmission in regional population centers.

Identification of shared bacterial strains in the vaginal microbiota of related and unrelated reproductive-age mothers and daughters using genome-resolved metagenomics

It has been suggested that the human microbiome might be vertically transmitted from mother to offspring and that early colonizers may play a critical role in development of the immune system. Studies have shown limited support for the vertical transmission of the intestinal microbiota but the derivation of the vaginal microbiota remains largely unknown. Although the vaginal microbiota of children and reproductive age women differ in composition, the vaginal microbiota could be vertically transmitted. To determine whether there was any support for this hypothesis, we examined the vaginal microbiota of daughter-mother pairs from the Baltimore metropolitan area (ages 14-27, 32-51; n = 39). We assessed whether the daughter's microbiota was similar in composition to their mother's using metataxonomics. Permutation tests revealed that while some pairs did have similar vaginal microbiota, the degree of similarity did not exceed that expected by chance. Genome-resolved metagenomics was used to identify shared bacterial strains in a subset of the families (n = 22). We found a small number of bacterial strains that were shared between mother-daughter pairs but identified more shared strains between individuals from different families, indicating that vaginal bacteria may display biogeographic patterns. Earlier-in-life studies are needed to demonstrate vertical transmission of the vaginal microbiota.

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'.

Detecting bacterial adaptation within individual microbiomes

The human microbiome harbours a large capacity for within-person adaptive mutations. Commensal bacterial strains can stably colonize a person for decades, and billions of mutations are generated daily within each person's microbiome. Adaptive mutations emerging during health might be driven by selective forces that vary across individuals, vary within an individual, or are completely novel to the human population. Mutations emerging within individual microbiomes might impact the immune system, the metabolism of nutrients or drugs, and the stability of the community to perturbations. Despite this potential, relatively little attention has been paid to the possibility of adaptive evolution within complex human-associated microbiomes. This review discusses the promise of studying within-microbiome adaptation, the conceptual and technical limitations that may have contributed to an underappreciation of adaptive de novo mutations occurring within microbiomes to date, and methods for detecting recent adaptive evolution. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.

Emergence of Tn1999.7, a New Transposon in blaOXA-48-Harboring Plasmids Associated with Increased Plasmid Stability

OXA-48 is the most common carbapenemase in Enterobacterales in Germany and many other European countries. Depending on the genomic location of bla_OXA-48, OXA-48-producing isolates vary in phenotype and intra- and interspecies transferability of bla_OXA-48. In most bacterial isolates, bla_OXA-48 is located on one of seven variants of Tn1999 (Tn1999.1 to Tn1999.6 and invTn1999.2). Here, a novel Tn1999 variant, Tn1999.7, is described, which was identified in 11 clinical isolates from 2016 to 2020. Tn1999.7 differs from Tn1999.1 by the insertion of the 8,349-bp Tn3 family transposon Tn7442 between the lysR gene and bla_OXA-48 open reading frame. Tn7442 carries genes coding for a restriction endonuclease and a DNA methyltransferase as cargo, forming a type III restriction modification system. Tn1999.7 was carried on an ~71-kb IncL plasmid in 9/11 isolates. In one isolate, Tn1999.7 was situated on an ~76-kb plasmid, harboring an additional insertion sequence in the plasmid backbone. In one isolate, the plasmid size is only ~63 kb due to a deletion adjacent to Tn7442 that extends into the plasmid backbone. Mean conjugation rates of the Tn1999.7-harboring plasmids in J53 ranged from 4.47 × 10^-5 to 2.03 × 10^-2, similar to conjugation rates of other pOXA-48-type IncL plasmids. The stability of plasmids with Tn1999.7 was significantly higher than that of a Tn1999.2-harboring plasmid in vitro. This increase in stability could be related to the insertion of a restriction-modification system, which can promote postsegregational killing. The increased plasmid stability associated with Tn1999.7 could contribute to the further spread of OXA-48.

Comment on "Biodegradation of Polystyrene by Pseudomonas sp. Isolated from the Gut of Superworms (Larvae of Zophobas atratus)"

Kim et al. isolated a polystyrene-degrading bacterium from the gut of the superworm (Zophobas atratus) and identified it as Pseudomonas aeruginosa DSM 50071. However, the original P. aeruginosa DSM 50071 was isolated before 1946 from an unknown source by researchers from Merck Sharp & Dohme. To claim that the superworm bacterial isolate is the same P. aeruginosa DSM 50071 strain is misleading. The two strains are different, and the original DSM 50071 strain may or may not degrade polystyrene. The recommendation to Kim et al. is to label the strain isolated from the superworm with a unique code, to re-sequence its 16S rRNA gene to near completeness, and to deposit the culture in a public collection so that other researchers and society may benefit from its study.

Colonization dynamics of extended-spectrum beta-lactamase-producing Enterobacterales in the gut of Malawian adults

Drug-resistant bacteria of the order Enterobacterales which produce extended-spectrum beta-lactamase enzymes (ESBL-Enterobacterales, ESBL-E) are global priority pathogens. Antimicrobial stewardship interventions proposed to curb their spread include shorter courses of antimicrobials to reduce selection pressure but individual-level acquisition and selection dynamics are poorly understood. We sampled stool of 425 adults (aged 16-76 years) in Blantyre, Malawi, over 6 months and used multistate modelling and whole-genome sequencing to understand colonization dynamics of ESBL-E. Models suggest a prolonged effect of antimicrobials such that truncating an antimicrobial course at 2 days has a limited effect in reducing colonization. Genomic analysis shows largely indistinguishable diversity of healthcare-associated and community-acquired isolates, hence some apparent acquisition of ESBL-E during hospitalization may instead represent selection from a patient's microbiota by antimicrobial exposure. Our approach could help guide stewardship protocols; interventions that aim to review and truncate courses of unneeded antimicrobials may be of limited use in preventing ESBL-E colonization.

First isolation of ST398 methicillin-resistant Staphylococcus aureus carrying staphylococcal cassette chromosome mec type IVd from pig ears in Japan

The emergence and increasing prevalence of livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) are a global concern. To investigate the prevalence and characteristics of sequence type 398 (ST398) MRSA in pig ears, 102 pig's ears were collected from 102 animals shipped from 51 farms at an abattoir. Eight ST398 MRSA isolates were isolated from the ears of eight pigs shipped from seven farms. Of the eight ST398 isolates, seven had the staphylococcal cassette chromosome mec (SCCmec) type IVd and these were obtained from seven pigs shipped from six farms. Single nucleotide polymorphisms ranging from 13 to 26 were observed in the core-genome regions in the seven SCCmec type IVd isolates. We believe that this is the first report on the isolation of ST398 MRSA SCCmec type IVd in Japan.

Genomic epidemiology of Salmonella Typhi in Central Division, Fiji, 2012 to 2016

BACKGROUND

Typhoid fever is endemic in some Pacific Island Countries including Fiji and Samoa yet genomic surveillance is not routine in such settings. Previous studies suggested imports of the global H58 clade of Salmonella enterica var Typhi (Salmonella Typhi) contribute to disease in these countries which, given the MDR potential of H58, does not auger well for treatment. The objective of the study was to define the genomic epidemiology of Salmonella Typhi in Fiji.

METHODS

Genomic sequencing approaches were implemented to study the distribution of 255 Salmonella Typhi isolates from the Central Division of Fiji. We augmented epidemiological surveillance and Bayesian phylogenomic approaches with a multi-year typhoid case-control study to define geospatial patterns among typhoid cases.

FINDINGS

Genomic analyses showed Salmonella Typhi from Fiji resolved into 2 non-H58 genotypes with isolates from the two dominant ethnic groups, the Indigenous (iTaukei) and non-iTaukei genetically indistinguishable. Low rates of international importation of clones was observed and overall, there were very low levels an antibiotic resistance within the endemic Fijian typhoid genotypes. Genomic epidemiological investigations were able to identify previously unlinked case clusters. Bayesian phylodynamic analyses suggested that genomic variation within the larger endemic Salmonella Typhi genotype expanded at discreet times, then contracted.

INTERPRETATION

Cyclones and flooding drove 'waves' of typhoid outbreaks in Fiji which, through population aggregation, poor sanitation and water safety, and then mobility of the population, spread clones more widely. Minimal international importations of new typhoid clones suggest that targeted local intervention strategies may be useful in controlling endemic typhoid infection. These findings add to our understanding of typhoid transmission networks in an endemic island country with broad implications, particularly across Pacific Island Countries.

FUNDING

This work was supported by the Coalition Against Typhoid through the Bill and Melinda Gates Foundation [grant number OPP1017518], the Victorian Government, the National Health and Medical Research Council Australia, the Australian Research Council, and the Fiji Ministry of Health and Medical Services.

Evaluation of within-host evolution of methicillin-resistant Staphylococcus aureus (MRSA) by comparing cgMLST and SNP analysis approaches

Whole genome sequencing (WGS) of methicillin-resistant Staphylococcus aureus (MRSA) provides high-resolution typing, facilitating surveillance and outbreak investigations. The aim of this study was to evaluate the genomic variation rate in MRSA, by comparing commonly used core genome multilocus sequencing (cgMLST) against single nucleotide polymorphism (SNP) analyses. WGS was performed on 95 MRSA isolates, collected from 20 carriers during years 2003–2019. To assess variation and methodological-related differences, two different cgMLST schemes were obtained using Ridom SeqSphere+ and the cloud-based 1928 platform. In addition, two SNP methods, 1928 platform and Northern Arizona SNP Pipeline (NASP) were used. The cgMLST using Ridom SeqSphere+ and 1928 showed a median of 5.0 and 2.0 allele variants/year, respectively. In the SNP analysis, performed with two reference genomes COL and Newman, 1928 showed a median of 13 and 24 SNPs (including presumed recombination) and 3.8 respectively 4.0 SNPs (without recombination) per individual/year. Accordantly, NASP showed a median of 5.5 and 5.8 SNPs per individual/year. In conclusion, an estimated genomic variation rate of 2.0–5.8 genetic events per year (without recombination), is suggested as a general guideline to be used at clinical laboratories for surveillance and outbreak investigations independently of analysis approach used.

A Whole Genome Sequencing-Based Epidemiological Investigation of a Pregnancy-Related Invasive Listeriosis Case in Central Italy

Listeriosis is currently the fifth most common foodborne disease in Europe. Most cases are sporadic; however, outbreaks have also been reported. Compared to other foodborne infections, listeriosis has a modest incidence but can cause life-threatening complications, especially in elderly or immunocompromised people and pregnant women. In the latter case, the pathology can be the cause of premature birth or spontaneous abortion, especially if the fetus is affected during the first months of gestation. The causative agent of listeriosis, Listeria monocytogenes, is characterized by the innate ability to survive in the environment and in food, even in adverse conditions and for long periods. Ready-to-eat food represents the category most at risk for contracting listeriosis. This study presents the result of an investigation carried out on a case of maternal-fetal transmission of listeriosis which occurred in 2020 in central Italy and which was linked, with a retrospective approach, to other cases residing in the same city of the pregnant woman. Thanks to the use of next-generation sequencing methodologies, it was possible to identify an outbreak of infection, linked to the consumption of ready-to-eat sliced products sold in a supermarket in the investigated city.

The periodic table of fermented foods: limitations and opportunities

Fermentation is one of the oldest methods of food processing and accounts for a substantial proportion of human foods, including not only staple foods such as bread, cereal porridges or fermented legumes but also fermented vegetables, meats, fish and dairy, alcoholic beverages as well as coffee, cocoa and condiments such as vinegar, soy sauce and fish sauces. Adding the regional varieties to these diverse product categories makes for an almost immeasurable diversity of fermented foods. The periodic table of fermented foods aims to map this diversity on the 118 entries of the periodic table of chemical elements. While the table fails to represent the diversity of fermented foods, it represents major fermentation substrates, product categories, fermentation processes and fermentation organisms. This communication not only addresses limitations of the graphical display on a "periodic table of fermented foods", but also identifies opportunities that relate to questions that are facilitated by this graphical presentation: on the origin and purpose of food fermentation, which fermented foods represent "indigenous" foods, differences and similarities in the assembly of microbial communities in different fermentations, differences in the global preferences for food fermentation, the link between microbial diversity, fermentation time and product properties, and opportunities of using traditional food fermentations as template for development of new products. KEY POINTS: • Fermented foods are produced in an almost immeasurable diversity. • Fermented foods were mapped on a periodic table of fermented foods. • This table facilitates identification of communalities and differences of products.

Whole genome sequence analysis of Salmonella Typhi in Papua New Guinea reveals an established population of genotype 2.1.7 sensitive to antimicrobials

BACKGROUND

Typhoid fever, a systemic infection caused by Salmonella enterica serovar Typhi, remains a considerable public health threat in impoverished regions within many low- and middle-income settings. However, we still lack a detailed understanding of the emergence, population structure, molecular mechanisms of antimicrobial resistance (AMR), and transmission dynamics of S. Typhi across many settings, particularly throughout the Asia-Pacific islands. Here we present a comprehensive whole genome sequence (WGS) based overview of S. Typhi populations circulating in Papua New Guinea (PNG) over 30 years.

PRINCIPLE FINDINGS

Bioinformatic analysis of 86 S. Typhi isolates collected between 1980-2010 demonstrated that the population structure of PNG is dominated by a single genotype (2.1.7) that appears to have emerged in the Indonesian archipelago in the mid-twentieth century with minimal evidence of inter-country transmission. Genotypic and phenotypic data demonstrated that the PNG S. Typhi population appears to be susceptible to former first line drugs for treating typhoid fever (chloramphenicol, ampicillin and co-trimoxazole), as well as fluoroquinolones, third generation cephalosporins, and macrolides. PNG genotype 2.1.7 was genetically conserved, with very few deletions, and no evidence of plasmid or prophage acquisition. Genetic variation among this population was attributed to either single point mutations, or homologous recombination adjacent to repetitive ribosomal RNA operons.

SIGNIFICANCE

Antimicrobials remain an effective option for the treatment of typhoid fever in PNG, along with other intervention strategies including improvements to water, sanitation and hygiene (WaSH) related infrastructure and potentially the introduction of Vi-conjugate vaccines. However, continued genomic surveillance is warranted to monitor for the emergence of AMR within local populations, or the introduction of AMR associated genotypes of S. Typhi in this setting.

Considerations and best practices in animal science 16S ribosomal RNA gene sequencing microbiome studies

Microbiome studies in animal science using 16S rRNA gene sequencing have become increasingly common in recent years as sequencing costs continue to fall and bioinformatic tools become more powerful and user-friendly. The combination of molecular biology, microbiology, microbial ecology, computer science, and bioinformatics-in addition to the traditional considerations when conducting an animal science study-makes microbiome studies sometimes intimidating due to the intersection of different fields. The objective of this review is to serve as a jumping-off point for those animal scientists less familiar with 16S rRNA gene sequencing and analyses and to bring up common issues and concerns that arise when planning an animal microbiome study from design through analysis. This review includes an overview of 16S rRNA gene sequencing, its advantages, and its limitations; experimental design considerations such as study design, sample size, sample pooling, and sample locations; wet lab considerations such as field handing, microbial cell lysis, low biomass samples, library preparation, and sequencing controls; and computational considerations such as identification of contamination, accounting for uneven sequencing depth, constructing diversity metrics, assigning taxonomy, differential abundance testing, and, finally, data availability. In addition to general considerations, we highlight some special considerations by species and sample type.

regentrans: a framework and R package for using genomics to study regional pathogen transmission

Increasing evidence of regional pathogen transmission networks highlights the importance of investigating the dissemination of multidrug-resistant organisms (MDROs) across a region to identify where transmission is occurring and how pathogens move across regions. We developed a framework for investigating MDRO regional transmission dynamics using whole-genome sequencing data and created regentrans, an easy-to-use, open source R package that implements these methods (https://github.com/Snitkin-Lab-Umich/regentrans). Using a dataset of over 400 carbapenem-resistant isolates of Klebsiella pneumoniae collected from patients in 21 long-term acute care hospitals over a one-year period, we demonstrate how to use our framework to gain insights into differences in inter- and intra-facility transmission across different facilities and over time. This framework and corresponding R package will allow investigators to better understand the origins and transmission patterns of MDROs, which is the first step in understanding how to stop transmission at the regional level.