Population genomics of bacterial host adaptation

Heterogeneity in recombination rates and accessory gene co-occurrence distinguish Pseudomonas aeruginosa phylogroups

Pseudomonas aeruginosa (class Gammaproteobacteria) is a ubiquitous, ecologically widespread, and metabolically versatile species. It is also an opportunistic pathogen that causes a variety of chronic and acute infections in humans. Its ability to thrive in diverse environments and exhibit a wide range of phenotypes lies in part on its large gene pool, but the processes that govern inter-strain genomic variation remain unclear. Here, we aim to characterize the recombination features and accessory genome structure of P. aeruginosa using 840 globally distributed genome sequences. The species can be subdivided into five phylogenetic sequence clusters (corresponding to known phylogroups), two of which are most prominent. Notable epidemic clones are found in the two phylogroups: ST17, ST111, ST146, ST274, and ST395 in phylogroup 1, and ST235 and ST253 in phylogroup 2. The two phylogroups differ in the frequency and characteristics of homologous recombination in their core genomes, including the specific genes that most frequently recombine and the impact of recombination on sequence diversity. Each phylogroup's accessory genome is characterized by a unique gene pool, co-occurrence networks of shared genes, and anti-phage defense systems. Different pools of antimicrobial resistance and virulence genes exist in the two phylogroups and display dissimilar patterns of co-occurrence. Altogether, our results indicate that each phylogroup displays distinct histories and patterns of acquiring exogenous DNA, which may contribute in part to their predominance in the global population. Our study has important implications for understanding the genome dynamics, within-species heterogeneity, and clinically relevant traits of P. aeruginosa.

IMPORTANCE

The consummate opportunist Pseudomonas aeruginosa inhabits many nosocomial and non-clinical environments, posing a major health burden worldwide. Our study reveals phylogroup-specific differences in recombination features and co-occurrence networks of accessory genes within the species. This genomic variation partly explains its remarkable ability to exhibit diverse ecological and phenotypic traits, and thus contribute to circumventing clinical and public health intervention strategies to contain it. Our results may help inform efforts to control and prevent P. aeruginosa diseases, including managing transmission, therapeutic efforts, and pathogen circulation in non-clinical environmental reservoirs.

Comparative Genomic Analysis of Two Vibrio harveyi Strains from Larimichthys crocea with Divergent Virulence Profiles

Vibrio harveyi is a significant pathogen in marine aquaculture, causing vibriosis in various marine species. This study presents a comparative genomic analysis of two V. harveyi strains, N8T11 and 45T2, which exhibit differing virulence profiles. Virulence assays revealed that N8T11 caused 92% mortality in infected fish, while 45T2 resulted in 0% mortality. Whole-genome sequencing revealed that strain N8T11 harbors five plasmids (pN8T11a, pN8T11b, pN8T11c, pN8T11d and pN8T11e) absent in 45T2, encoding genes potentially linked to virulence, such as siderophore-mediated iron acquisition and stress response mechanisms. Pan-genome analysis highlighted substantial genomic plasticity within V. harveyi, with mobile genetic elements, including plasmids and prophages, contributing to horizontal gene transfer. Conjugation experiments demonstrated that all five N8T11 plasmids can transfer to 45T2 with efficiencies up to 87%, with pN8T11b remaining stable across multiple subcultures, enabling the dissemination of virulence-associated genes. These findings suggest that plasmid-mediated gene transfer plays a role in the virulence variability observed between V. harveyi strains. This study contributes to understanding the genomic factors underlying pathogenicity in V. harveyi and provides insights for future research aimed at controlling vibriosis in aquaculture.

Biofilm-disrupting DNA nanomedicines for targeted elimination of resistant wound microbiota

Biofilms are complex bacterial communities that significantly hinder the treatment of chronic and recurrent infections by enhancing bacterial virulence and conferring resistance to antimicrobial therapies. To address this challenge, an intelligent DNA nanomedicine has been engineered to dismantle biofilms and target resistant bacteria, offering an innovative solution for chronic wound infections. These nanomedicines initiate biofilm degradation through in situ generation of potent oxidative radicals, enabling deep biofilm penetration and precise bacterial targeting. Utilizing aptamers for specific bacterial identification, the nanomedicines concentrate therapeutic agents directly at infection sites. The combined effect of severe oxidative stress and sustained silver ion release ensures a continuous, focused assault on pathogens, effectively eradicating resistant bacteria. This strategy demonstrated broad-spectrum efficacy against both Gram-positive and Gram-negative bacteria, significantly enhancing wound healing in a diabetic infection model. By integrating intelligent bacterial eradication with modulation of the wound microenvironment, this approach presents a promising solution for overcoming biofilm-associated resistance and advancing chronic wound infection treatment.

Emerging genetic diversity of Mycoplasma gallisepticum in Mexican house finches: Evidence of possible independent spillover events

In a previous study, we reported for the first time the detection of Mycoplasma gallisepticum (MG) in Mexican house finches (Haemorhous mexicanus). Building on this finding, we explored the genetic diversity of MG, addressing the potential independent spillover events. Samples from 247 wild finches across Mexico revealed MG infection in 72 % of choanal swabs and 24 % of conjunctival swabs, with no clinical signs observed. Phylogenetic analysis identified two novel MG clades distinct from U.S. house finch-associated and poultry-associated MG strains, suggesting independent evolution within Mexico. Coinfections with diverse haplotypes were common, raising concerns about recombination and shifts in virulence. This research highlights the asymptomatic carrier state of Mexican house finches, their potential as reservoirs, and the ecological implications of pathogen spread and adaptation. These findings underscore the need for enhanced surveillance and further study on MG's dynamics in Mexican avifauna.

Genomic evidence of symbiotic adaptations in fungus-associated bacteria

Fungi harbor diverse bacteria that engage in various relationships. While these relationships potentially influence fungal functioning, their underlying genetic mechanisms remain unexplored. Here, we aimed to elucidate the key genomic features of fungus-associated bacteria (FaB) by comparing 163 FaB genomes to 1,048 bacterial genomes from other hosts and habitats. Our analyses revealed several distinctive genomic features of FaB. We found that FaB are enriched in carbohydrate transport/metabolism- and motility-related genes, suggesting an adaptation for utilizing complex fungal carbon sources. They are also enriched in genes targeting fungal biomass, likely reflecting their role in recycling and rebuilding fungal structures. Additionally, FaB associated with plant-mutualistic fungi possess a wider array of carbon-acquisition enzymes specific to fungal and plant substrates compared to those residing with saprotrophic fungi. These unique genomic features highlight FaB' potential as key players in fungal nutrient acquisition and decomposition, ultimately influencing plant-fungal symbiosis and ecosystem functioning.

Methicillin-sensitive Staphylococcus aureus lineages contribute towards poor patient outcomes in orthopaedic device-related infections

Staphylococci are the most common cause of orthopaedic device-related infections (ODRIs), with Staphylococcus aureus responsible for a third or more of cases. This prospective clinical and laboratory study investigated the association of genomic and phenotypic variation with treatment outcomes in ODRI isolates. Eighty-six invasive S. aureus isolates were collected from patients with ODRI, and clinical outcome was assessed after a follow-up examination of 24 months. Each patient was then considered to have been 'cured' or 'not cured' based on predefined clinical criteria. Whole-genome sequencing and molecular characterization identified isolates belonging to globally circulating community- and hospital-acquired lineages. Most isolates were phenotypically susceptible to methicillin and lacked the staphylococcal cassette chromosome mec cassette [methicillin-susceptible S. aureus (MSSA); 94%] but contained several virulence genes, including toxins and biofilm genes. Whilst recognizing the role of the host immune response, we identified genetic variance, which could be associated with the infection severity or clinical outcome. Whilst this and several other studies reinforce the role antibiotic resistance [e.g. methicillin-resistant S. aureus (MRSA) infection] has on treatment failure, it is important not to overlook MSSA that can cause equally destructive infections and lead to poor patient outcomes.

Streamlining whole genome sequencing for clinical diagnostics with ONT technology

Recent advances in whole-genome sequencing (WGS) have increased the accessibility of this tool, offering substantial potential for pathogen surveillance, outbreak response, and diagnostics. However, the routine clinical adoption of WGS is hindered by factors such as high costs, technical complexity, and the requirement for bioinformatics expertise for data analysis. To address these challenges, we propose RapidONT, a workflow designed for cost-effective and accessible WGS-based pathogen analysis. RapidONT employs a mechanical shearing-based DNA extraction protocol, followed by library construction by using a multiplexing Oxford nanopore technologies (ONT) rapid barcoding kit. Flye software is used for de novo assembly without manual intervention, followed by basic assembly polishing using Medaka and Homopolish. The polished assemblies are then analyzed using the user-friendly web-based platform Pathogenwatch, which facilitates species identification, molecular typing, and antimicrobial resistance (AMR) prediction, all while requiring minimal bioinformatics expertise. The efficacy of RapidONT was evaluated using nine clinically relevant pathogens, encompassing a total of 90 gram-positive and gram-negative bacterial strains. The workflow demonstrated high accuracy in critical tasks such as multilocus sequence typing (MLST) and AMR identification, using only ONT R9.4.1 flowcell data. Notably, limitations were observed with Salmonella spp. and Neisseria gonorrhoeae. Furthermore, RapidONT enabled the generation of genomic information for 48 bacterial isolates by using a single flow cell, significantly reducing sequencing costs. This approach eliminates the need for extensive experimentation in obtaining crucial genomic information. This workflow facilitates broader WGS implementation in clinical pathogen analysis and diagnostics.

A One Health approach for the genomic characterization of antibiotic-resistant Campylobacter isolates using Nanopore whole-genome sequencing

In response to the growing threat posed by the spread of antimicrobial resistance in zoonotic Campylobacter, a One Health approach was used to examine the genomic diversity, phylogenomic relationships, and the distribution of genetic determinants of resistance (GDR) in C. jejuni and C. coli isolates from humans, animals (ruminants, swine, and chickens), and avian food products collected during a regionally (Basque Country, Spain) and temporally (mostly 2021-2022) restricted sampling. Eighty-three C. jejuni and seventy-one C. coli isolates, most exhibiting resistance to ciprofloxacin and/or erythromycin, were whole-genome sequenced using Oxford Nanopore Technologies long-fragment sequencing (ONT). Multilocus sequence typing (MLST) analysis identified a high genomic diversity among isolates. Phylogenomic analysis showed that clustering based on the core genome was aligned with MLST profiles, regardless of the sample source. In contrast, accessory genome content sometimes discriminated isolates within the same STs and occasionally differentiated isolates from different sources. The majority of the identified GDRs were present in isolates from different sources, and a good correlation was observed between GDR distribution and phenotypic susceptibility profiles (based on minimum inhibitory concentrations interpreted according to the EUCAST epidemiological cutoff values). Genotypic resistance profiles were independent of genotypes, indicating no apparent association between resistance and phylogenetic origin. This study demonstrates that ONT sequencing is a powerful tool for molecular surveillance of bacterial pathogens in the One Health framework.

CspZ variant-specific interaction with factor H incorporates a metal site to support Lyme borreliae complement evasion

Polymorphic microbial immune evasion proteins dictate the pathogen species- or strain-specific virulence. Metals can impact how microbial proteins confer host-pathogen interactions, but whether this activity can be allelically variable is unclear. Here, we investigate the polymorphic CspZ protein of Lyme disease spirochete bacteria to assess the role of metals in protein-protein interaction. CspZ facilitates evasion of the complement system, the first line of immune defense through binding to the complement regulator factor H (FH). By obtaining a high-resolution cocrystal CspZ-FH structure, we identified a zinc coordinating the binding of FH SCR6-7 domains to a Glu65 on a loop from CspZ of Borrelia burgdorferi B31. However, zinc is dispensable for human FH binding for CspZ orthologs with a different loop orientation and/or lacking this glutamate. Phylogenetic analysis of all known human FH-binding CspZ variants further grouped the proteins into three unique lineages correlating with loop sequences. This suggests multiple FH-binding mechanisms evolved through Lyme disease spirochete-host interactions. Overall, this multidisciplinary work elucidates how the allelically specific immune evasion role of metals is impacted by microbial protein polymorphisms.

Evolution and spread of Xanthomonas citri subsp. citri in the São Paulo, Brazil, citrus belt inferred from 758 novel genomes

The São Paulo state citrus belt in Brazil is a major citrus production region. Since at least 1957, citrus plantations in this region have been affected by citrus canker, an economically damaging disease caused by Xanthomonas citri subsp. citri (Xcc). For about 50 years, until 2017, a citrus canker eradication programme was carried out in this region. In this work, our aim was to investigate the effects of the eradication programme on genetic variability and evolution of Xcc. To this end, we sequenced and analysed 758 Xcc genomes sampled in the São Paulo citrus belt, together with 730 publicly available Xcc genomes from around the world. Our phylogenomic analyses show that these genomes can be grouped into seven major lineages and that in São Paulo, lineage L7 is dominant. Our time estimate for its appearance closely matches the date when citrus production expanded. L7 can be subdivided into lineages L7.1 and L7.2. In our samples, L7.2, which we estimate to have emerged around 1964, is by far the most abundant, showing that the eradication programme had little impact on strain diversification. On the other hand, oscillations in the estimated effective population size of L7.2 strains over time closely match the shifts in the eradication programme. In sum, we present a detailed view of the genomic diversity of Xcc in the world and in São Paulo, the largest such effort in terms of a number of genomes for a crop pathogen undertaken so far. The methods employed here can form the basis for active genomic surveillance of Xcc in major citrus production areas.

Using GWAS and Machine Learning to Identify and Predict Genetic Variants Associated with Foodborne Bacteria Phenotypic Traits

One of the main challenges in food microbiology is to prevent the risk of outbreaks by avoiding the distribution of food contaminated by bacteria. This requires constant monitoring of the circulating strains throughout the food production chain. Bacterial genomes contain signatures of natural evolution and adaptive markers that can be exploited to better understand the behavior of pathogen in the food industry. The monitoring of foodborne strains can therefore be facilitated by the use of these genomic markers capable of rapidly providing essential information on isolated strains, such as the source of contamination, risk of illness, potential for biofilm formation, and tolerance or resistance to biocides. The increasing availability of large genome datasets is enhancing the understanding of the genetic basis of complex traits such as host adaptation, virulence, and persistence. Genome-wide association studies have shown very promising results in the discovery of genomic markers that can be integrated into rapid detection tools. In addition, machine learning has successfully predicted phenotypes and classified important traits. Genome-wide association and machine learning tools have therefore the potential to support decision-making circuits intending at reducing the burden of foodborne diseases. The aim of this chapter review is to provide knowledge on the use of these two methods in food microbiology and to recommend their use in the field.

Simple and accurate genomic classification model for distinguishing between human and pig Staphylococcus aureus

Staphylococcus aureus (S. aureus) can cause various infections in humans and animals, contributing to high morbidity and mortality. To prevent and control cross-species transmission of S. aureus, it is necessary to understand the host-associated genetic variants. We performed a two-stage genome-wide association study (GWAS) including initial screening and further validation to compare genomic differences between human and pig S. aureus, aiming to identify host-associated determinants. Our multiple GWAS analyses found six consensus significant k-mers associated with host species, providing novel genetic evidence for distinguishing human from pig S. aureus. The best k-mer predictor achieved a high classification accuracy of 98.12% on its own and had extremely high resolution similar to the SNPs-based phylogeny, offering a very simple target for predicting the cross-species transmission risk of S. aureus. The final k-mer model revealed that 90% of S. aureus isolates from farm workers were predicted as livestock origin, suggesting a high risk of cross-species transmission. Bayesian inference revealed different cross-species transmission directions, with the human-to-pig transmission for ST5 and the pig-to-human transmission for ST398. Our findings provide a simple and accurate k-mer model for identifying and predicting the cross-species transmission risk of S. aureus.

Proximity to humans is associated with antimicrobial-resistant enteric pathogens in wild bird microbiomes

Humans are radically altering global ecology, and one of the most apparent human-induced effects is urbanization, where high-density human habitats disrupt long-established ecotones. Changes to these transitional areas between organisms, especially enhanced contact among humans and wild animals, provide new opportunities for the spread of zoonotic pathogens. This poses a serious threat to global public health, but little is known about how habitat disruption impacts cross-species pathogen spread. Here, we investigated variation in the zoonotic enteric pathogen Campylobacter jejuni. The ubiquity of C. jejuni in wild bird gut microbiomes makes it an ideal organism for understanding how host behavior and ecology influence pathogen transition and spread. We analyzed 700 C. jejuni isolate genomes from 30 bird species in eight countries using a scalable generalized linear model approach. Comparing multiple behavioral and ecological traits showed that proximity to human habitation promotes lineage diversity and is associated with antimicrobial-resistant (AMR) strains in natural populations. Specifically, wild birds from urban areas harbored up to three times more C. jejuni genotypes and AMR genes. This study provides novel methodology and much-needed quantitative evidence linking urbanization to gene pool spread and zoonoses.

Host species shapes genotype, antimicrobial resistance, and virulence profiles of enterotoxigenic Escherichia coli (ETEC) from livestock in the United States

Enterotoxigenic Escherichia coli (ETEC) are significant pathogen in both cattle and pigs, causing diarrhea in these animals and leading to economic losses in the livestock industry. Understanding the dissimilarity in genotype, antimicrobial resistance (AMR), and virulence between bovine and swine ETEC is crucial for development of targeted preventive and therapeutic approaches for livestock. However, a comprehensive study on this area remains lacking. Here, we performed whole-genome sequencing-based analyses of bovine (n = 554) and swine (n = 623) ETEC collected in the United States over a 53-year period. We identified distinct ETEC genotypes (fimH type, O antigen, H antigen, sequence type) in cattle and pigs. Furthermore, specific AMR and virulence profiles were associated with bovine and swine ETEC. Compared to swine ETEC, bovine ETEC were less diverse in genotypes and had a significantly (P < 0.001) lower number of AMR genes per isolate but higher co-occurrence of Shiga toxin and enterotoxin genes. Our results provide an overview of the key genomic differences between bovine and swine ETEC in the United States, which might be attributed to host adaptation and antibiotic usage practice. Ongoing surveillance and research are essential to monitor the genetic diversity and AMR patterns of ETEC in different host species.

IMPORTANCE

Enterotoxigenic Escherichia coli (ETEC)-associated diarrhea represent one of the most economically important diseases in the livestock industry. By analyzing over a thousand livestock-derived ETEC samples in the United States, our study unveiled a clear distinction in ETEC's genetic traits (i.e., genotypes, antimicrobial resistance [AMR], and virulence profiles) that might be tied to the different use of antibiotics in cattle and pigs, and the bacteria's adaptation to their specific animal hosts. This understanding is crucial for tailoring preventive and therapeutic strategies. It also highlights the significance of ongoing surveillance and research into the evolution of bacterial pathogens like ETEC in livestock by using advanced techniques such as whole-genome sequencing.

Integrating research on bacterial pathogens and commensals to fight infections-an ecological perspective

The incidence of antibiotic-resistant bacterial infections is increasing, and development of new antibiotics has been deprioritised by the pharmaceutical industry. Interdisciplinary research approaches, based on the ecological principles of bacterial fitness, competition, and transmission, could open new avenues to combat antibiotic-resistant infections. Many facultative bacterial pathogens use human mucosal surfaces as their major reservoirs and induce infectious diseases to aid their lateral transmission to new host organisms under some pathological states of the microbiome and host. Beneficial bacterial commensals can outcompete specific pathogens, thereby lowering the capacity of the pathogens to spread and cause serious infections. Despite the clinical relevance, however, the understanding of commensal-pathogen interactions in their natural habitats remains poor. In this Personal View, we highlight directions to intensify research on the interactions between bacterial pathogens and commensals in the context of human microbiomes and host biology that can lead to the development of innovative and sustainable ways of preventing and treating infectious diseases.

Comparative genomic analysis identifies potential adaptive variation in Mycoplasma ovipneumoniae

Mycoplasma ovipneumoniae is associated with respiratory disease in wild and domestic Caprinae globally, with wide variation in disease outcomes within and between host species. To gain insight into phylogenetic structure and mechanisms of pathogenicity for this bacterial species, we compared M. ovipneumoniae genomes for 99 samples from 6 countries (Australia, Bosnia and Herzegovina, Brazil, China, France and USA) and 4 host species (domestic sheep, domestic goats, bighorn sheep and caribou). Core genome sequences of M. ovipneumoniae assemblies from domestic sheep and goats fell into two well-supported phylogenetic clades that are divergent enough to be considered different bacterial species, consistent with each of these two clades having an evolutionary origin in separate host species. Genome assemblies from bighorn sheep and caribou also fell within these two clades, indicating multiple spillover events, most commonly from domestic sheep. Pangenome analysis indicated a high percentage (91.4 %) of accessory genes (i.e. genes found only in a subset of assemblies) compared to core genes (i.e. genes found in all assemblies), potentially indicating a propensity for this pathogen to adapt to within-host conditions. In addition, many genes related to carbon metabolism, which is a virulence factor for Mycoplasmas, showed evidence for homologous recombination, a potential signature of adaptation. The presence or absence of annotated genes was very similar between sheep and goat clades, with only two annotated genes significantly clade-associated. However, three M. ovipneumoniae genome assemblies from asymptomatic caribou in Alaska formed a highly divergent subclade within the sheep clade that lacked 23 annotated genes compared to other assemblies, and many of these genes had functions related to carbon metabolism. Overall, our results suggest that adaptation of M. ovipneumoniae has involved evolution of carbon metabolism pathways and virulence mechanisms related to those pathways. The genes involved in these pathways, along with other genes identified as potentially involved in virulence in this study, are potential targets for future investigation into a possible genomic basis for the high variation observed in disease outcomes within and between wild and domestic host species.

Genomic and functional determinants of host spectrum in Group B Streptococcus

Group B Streptococcus (GBS) is a major human and animal pathogen that threatens public health and food security. Spill-over and spill-back between host species is possible due to adaptation and amplification of GBS in new niches but the evolutionary and functional mechanisms underpinning those phenomena are poorly known. Based on analysis of 1,254 curated genomes from all major GBS host species and six continents, we found that the global GBS population comprises host-generalist, host-adapted and host-restricted sublineages, which are found across host groups, preferentially within one host group, or exclusively within one host group, respectively, and show distinct levels of recombination. Strikingly, the association of GBS genomes with the three major host groups (humans, cattle, fish) is driven by a single accessory gene cluster per host, regardless of sublineage or the breadth of host spectrum. Moreover, those gene clusters are shared with other streptococcal species occupying the same niche and are functionally relevant for host tropism. Our findings demonstrate (1) the heterogeneity of genome plasticity within a bacterial species of public health importance, enabling the identification of high-risk clones; (2) the contribution of inter-species gene transmission to the evolution of GBS; and (3) the importance of considering the role of animal hosts, and the accessory gene pool associated with their microbiota, in the evolution of multi-host bacterial pathogens. Collectively, these phenomena may explain the adaptation and clonal expansion of GBS in animal reservoirs and the risk of spill-over and spill-back between animals and humans.

Mechanisms of host adaptation by bacterial pathogens

The emergence of new infectious diseases poses a major threat to humans, animals, and broader ecosystems. Defining factors that govern the ability of pathogens to adapt to new host species is therefore a crucial research imperative. Pathogenic bacteria are of particular concern, given dwindling treatment options amid the continued expansion of antimicrobial resistance. In this review, we summarize recent advancements in the understanding of bacterial host species adaptation, with an emphasis on pathogens of humans and related mammals. We focus particularly on molecular mechanisms underlying key steps of bacterial host adaptation including colonization, nutrient acquisition, and immune evasion, as well as suggest key areas for future investigation. By developing a greater understanding of the mechanisms of host adaptation in pathogenic bacteria, we may uncover new strategies to target these microbes for the treatment and prevention of infectious diseases in humans, animals, and the broader environment.

Epistasis, core-genome disharmony, and adaptation in recombining bacteria

Recombination of short DNA fragments via horizontal gene transfer (HGT) can introduce beneficial alleles, create genomic disharmony through negative epistasis, and create adaptive gene combinations through positive epistasis. For non-core (accessory) genes, the negative epistatic cost is likely to be minimal because the incoming genes have not co-evolved with the recipient genome and are frequently observed as tightly linked cassettes with major effects. By contrast, interspecific recombination in the core genome is expected to be rare because disruptive allelic replacement is likely to introduce negative epistasis. Why then is homologous recombination common in the core of bacterial genomes? To understand this enigma, we take advantage of an exceptional model system, the common enteric pathogens Campylobacter jejuni and C. coli that are known for very high magnitude interspecies gene flow in the core genome. As expected, HGT does indeed disrupt co-adapted allele pairings, indirect evidence of negative epistasis. However, multiple HGT events enable recovery of the genome's co-adaption between introgressing alleles, even in core metabolism genes (e.g., formate dehydrogenase). These findings demonstrate that, even for complex traits, genetic coalitions can be decoupled, transferred, and independently reinstated in a new genetic background-facilitating transition between fitness peaks. In this example, the two-step recombinational process is associated with C. coli that are adapted to the agricultural niche.IMPORTANCEGenetic exchange among bacteria shapes the microbial world. From the acquisition of antimicrobial resistance genes to fundamental questions about the nature of bacterial species, this powerful evolutionary force has preoccupied scientists for decades. However, the mixing of genes between species rests on a paradox: 0n one hand, promoting adaptation by conferring novel functionality; on the other, potentially introducing disharmonious gene combinations (negative epistasis) that will be selected against. Taking an interdisciplinary approach to analyze natural populations of the enteric bacteria Campylobacter, an ideal example of long-range admixture, we demonstrate that genes can independently transfer across species boundaries and rejoin in functional networks in a recipient genome. The positive impact of two-gene interactions appears to be adaptive by expanding metabolic capacity and facilitating niche shifts through interspecific hybridization. This challenges conventional ideas and highlights the possibility of multiple-step evolution of multi-gene traits by interspecific introgression.

Genomic tailoring of autogenous poultry vaccines to reduce Campylobacter from farm to fork

Campylobacter is a leading cause of food-borne gastroenteritis worldwide, linked to the consumption of contaminated poultry meat. Targeting this pathogen at source, vaccines for poultry can provide short-term caecal reductions in Campylobacter numbers in the chicken intestine. However, this approach is unlikely to reduce Campylobacter in the food chain or human incidence. This is likely as vaccines typically target only a subset of the high genomic strain diversity circulating among chicken flocks, and rapid evolution diminishes vaccine efficacy over time. To address this, we used a genomic approach to develop a whole-cell autogenous vaccine targeting isolates harbouring genes linked to survival outside of the host. We hyper-immunised a whole major UK breeder farm to passively target offspring colonisation using maternally-derived antibody. Monitoring progeny, broiler flocks revealed a near-complete shift in the post-vaccination Campylobacter population with an ~50% reduction in isolates harbouring extra-intestinal survival genes and a significant reduction of Campylobacter cells surviving on the surface of meat. Based on these findings, we developed a logistic regression model that predicted that vaccine efficacy could be extended to target 65% of a population of clinically relevant strains. Immuno-manipulation of poultry microbiomes towards less harmful commensal isolates by competitive exclusion, has major potential for reducing pathogens in the food production chain.

Analysis of twelve genomes of the bacterium Kerstersia gyiorum from brown-throated sloths (Bradypus variegatus), the first from a non-human host

Kerstersia gyiorum is a Gram-negative bacterium found in various animals, including humans, where it has been associated with various infections. Knowledge of the basic biology of K. gyiorum is essential to understand the evolutionary strategies of niche adaptation and how this organism contributes to infectious diseases; however, genomic data about K. gyiorum is very limited, especially from non-human hosts. In this work, we sequenced 12 K. gyiorum genomes isolated from healthy free-living brown-throated sloths (Bradypus variegatus) in the Parque Estadual das Fontes do Ipiranga (São Paulo, Brazil), and compared them with genomes from isolates of human origin, in order to gain insights into genomic diversity, phylogeny, and host specialization of this species. Phylogenetic analysis revealed that these K. gyiorum strains are structured according to host. Despite the fact that sloth isolates were sampled from a single geographic location, the intra-sloth K. gyiorum diversity was divided into three clusters, with differences of more than 1,000 single nucleotide polymorphisms between them, suggesting the circulation of various K. gyiorum lineages in sloths. Genes involved in mobilome and defense mechanisms against mobile genetic elements were the main source of gene content variation between isolates from different hosts. Sloth-specific K. gyiorum genome features include an IncN2 plasmid, a phage sequence, and a CRISPR-Cas system. The broad diversity of defense elements in K. gyiorum (14 systems) may prevent further mobile element flow and explain the low amount of mobile genetic elements in K. gyiorum genomes. Gene content variation may be important for the adaptation of K. gyiorum to different host niches. This study furthers our understanding of diversity, host adaptation, and evolution of K. gyiorum, by presenting and analyzing the first genomes of non-human isolates.

Resequencing and characterization of the first Corynebacterium pseudotuberculosis genome isolated from camel

Background

Corynebacterium pseudotuberculosis is a zoonotic Gram-positive bacterial pathogen known to cause different diseases in many mammals, including lymph node abscesses in camels. Strains from biovars equi and ovis of C. pseudotuberculosis can infect camels. Comparative genomics could help to identify features related to host adaptation, and currently strain Cp162 from biovar equi is the only one from camel with a sequenced genome.

Methods

In this work, we compared the quality of three genome assemblies of strain Cp162 that used data from the DNA sequencing platforms SOLiD v3 Plus, IonTorrent PGM, and Illumina HiSeq 2500 with an optical map and investigate the unique features of this strain. For this purpose, we applied comparative genomic analysis on the different Cp162 genome assembly versions and included other 129 genomes from the same species.

Results

Since the first version of the genome, there was an increase of 88 Kbp and 121 protein-coding sequences, a decrease of pseudogenes from 139 to 53, and two inversions and one rearrangement corrected. We identified 30 virulence genes, none associated to the camel host, and the genes rpob2 and rbpA predicted to confer resistance to rifampin. In comparison to 129 genomes of the same species, strain Cp162 has four genes exclusively present, two of them code transposases and two truncated proteins, and the three exclusively absent genes lysG, NUDIX domain protein, and Hypothetical protein. All 130 genomes had the rifampin resistance genes rpob2 and rbpA. Our results found no unique gene that could be associated with tropism to camel host, and further studies should include more genomes and genome-wide association studies testing for genes and SNPs.

Purifying selection drives distinctive arsenic metabolism pathways in prokaryotic and eukaryotic microbes

Microbes play a crucial role in the arsenic biogeochemical cycle through specific metabolic pathways to adapt to arsenic toxicity. However, the different arsenic-detoxification strategies between prokaryotic and eukaryotic microbes are poorly understood. This hampers our comprehension of how microbe-arsenic interactions drive the arsenic cycle and the development of microbial methods for remediation. In this study, we utilized conserved protein domains from 16 arsenic biotransformation genes (ABGs) to search for homologous proteins in 670 microbial genomes. Prokaryotes exhibited a wider species distribution of arsenic reduction- and arsenic efflux-related genes than fungi, whereas arsenic oxidation-related genes were more prevalent in fungi than in prokaryotes. This was supported by significantly higher acr3 (arsenite efflux permease) expression in bacteria (upregulated 3.72-fold) than in fungi (upregulated 1.54-fold) and higher aoxA (arsenite oxidase) expression in fungi (upregulated 5.11-fold) than in bacteria (upregulated 2.05-fold) under arsenite stress. The average values of nonsynonymous substitutions per nonsynonymous site to synonymous substitutions per synonymous site (dN/dS) of homologous ABGs were higher in archaea (0.098) and bacteria (0.124) than in fungi (0.051). Significant negative correlations between the dN/dS of ABGs and species distribution breadth and gene expression levels in archaea, bacteria, and fungi indicated that microbes establish the distinct strength of purifying selection for homologous ABGs. These differences contribute to the distinct arsenic metabolism pathways in prokaryotic and eukaryotic microbes. These observations facilitate a significant shift from studying individual or several ABGs to characterizing the comprehensive microbial strategies of arsenic detoxification.

Draft genome sequencing of Tilletia caries inciting common bunt of wheat provides pathogenicity-related genes

Common bunt of wheat caused by Tilletia caries is an important disease worldwide. The T. caries TC1_MSG genome was sequenced using the Illumina HiSeq 2500 and Nanopore ONT platforms. The Nanopore library was prepared using the ligation sequencing kit SQK-LSK110 to generate approximately 24 GB for sequencing. The assembly size of 38.18 Mb was generated with a GC content of 56.10%. The whole genome shotgun project was deposited at DDBJ/ENA/GenBank under the accession number JALUTQ000000000. Forty-six contigs were obtained with N50 of 1,798,756 bp. In total, 10,698 genes were predicted in the assembled genome. Out of 10,698 genes, 10,255 genes were predicted significantly in the genome. The repeat sequences made up approximately 1.57% of the genome. Molecular function, cellular components, and biological processes for predicted genes were mapped into the genome. In addition, repeat elements in the genome were assessed. In all, 0.89% of retroelements were observed, followed by long terminal repeat elements (0.86%) in the genome. In simple sequence repeat (SSR) analysis, 8,582 SSRs were found in the genome assembly. The trinucleotide SSR type (3,703) was the most abundant. Few putative secretory signal peptides and pathogenicity-related genes were predicted. The genomic information of T. caries will be valuable in understanding the pathogenesis mechanism as well as developing new methods for the management of the common bunt disease of wheat.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

Contrasting genes conferring short- and long-term biofilm adaptation in Listeria

Listeria monocytogenes is an opportunistic food-borne bacterium that is capable of infecting humans with high rates of hospitalization and mortality. Natural populations are genotypically and phenotypically variable, with some lineages being responsible for most human infections. The success of L. monocytogenes is linked to its capacity to persist on food and in the environment. Biofilms are an important feature that allow these bacteria to persist and infect humans, so understanding the genetic basis of biofilm formation is key to understanding transmission. We sought to investigate the biofilm-forming ability of L. monocytogenes by identifying genetic variation that underlies biofilm formation in natural populations using genome-wide association studies (GWAS). Changes in gene expression of specific strains during biofilm formation were then investigated using RNA sequencing (RNA-seq). Genetic variation associated with enhanced biofilm formation was identified in 273 genes by GWAS and differential expression in 220 genes by RNA-seq. Statistical analyses show that the number of overlapping genes flagged by either type of experiment is less than expected by random sampling. This novel finding is consistent with an evolutionary scenario where rapid adaptation is driven by variation in gene expression of pioneer genes, and this is followed by slower adaptation driven by nucleotide changes within the core genome.

Evidence of potential Campylobacter jejuni zooanthroponosis in captive macaque populations

Non-human primates share recent common ancestry with humans and exhibit comparable disease symptoms. Here, we explored the transmission potential of enteric bacterial pathogens in monkeys exhibiting symptoms of recurrent diarrhoea in a biomedical research facility in China. The common zoonotic bacterium Campylobacter jejuni was isolated from macaques (Macaca mulatta and Macaca fascicularis) and compared to isolates from humans and agricultural animals in Asia. Among the monkeys sampled, 5 % (44/973) tested positive for C. jejuni, 11 % (5/44) of which displayed diarrhoeal symptoms. Genomic analysis of monkey isolates, and 1254 genomes from various sources in Asia, were used to identify the most likely source of human infection. Monkey and human isolates shared high average nucleotide identity, common MLST clonal complexes and clustered together on a phylogeny. Furthermore, the profiles of putative antimicrobial resistance genes were similar between monkeys and humans. Taken together these findings suggest that housed macaques became infected with C. jejuni either directly from humans or via a common contamination source.

What are the missing pieces needed to stop antibiotic resistance?

As recognized by several international agencies, antibiotic resistance is nowadays one of the most relevant problems for human health. While this problem was alleviated with the introduction of new antibiotics into the market in the golden age of antimicrobial discovery, nowadays few antibiotics are in the pipeline. Under these circumstances, a deep understanding on the mechanisms of emergence, evolution and transmission of antibiotic resistance, as well as on the consequences for the bacterial physiology of acquiring resistance is needed to implement novel strategies, beyond the development of new antibiotics or the restriction in the use of current ones, to more efficiently treat infections. There are still several aspects in the field of antibiotic resistance that are not fully understood. In the current article, we make a non-exhaustive critical review of some of them that we consider of special relevance, in the aim of presenting a snapshot of the studies that still need to be done to tackle antibiotic resistance.

Multi-host infection and phylogenetically diverse lineages shape the recombination and gene pool dynamics of Staphylococcus aureus

BACKGROUND

Staphylococcus aureus can infect and adapt to multiple host species. However, our understanding of the genetic and evolutionary drivers of its generalist lifestyle remains inadequate. This is particularly important when considering local populations of S. aureus, where close physical proximity between bacterial lineages and between host species may facilitate frequent and repeated interactions between them. Here, we aim to elucidate the genomic differences between human- and animal-derived S. aureus from 437 isolates sampled from disease cases in the northeast region of the United States.

RESULTS

Multi-locus sequence typing revealed the existence of 75 previously recognized sequence types (ST). Our population genomic analyses revealed heterogeneity in the accessory genome content of three dominant S. aureus lineages (ST5, ST8, ST30). Genes related to antimicrobial resistance, virulence, and plasmid types were differentially distributed among isolates according to host (human versus non-human) and among the three major STs. Across the entire population, we identified a total of 1,912 recombination events that occurred in 765 genes. The frequency and impact of homologous recombination were comparable between human- and animal-derived isolates. Low-frequency STs were major donors of recombined DNA, regardless of the identity of their host. The most frequently recombined genes (clfB, aroA, sraP) function in host infection and virulence, which were also frequently shared between the rare lineages.

CONCLUSIONS

Taken together, these results show that frequent but variable patterns of recombination among co-circulating S. aureus lineages, including the low-frequency lineages, that traverse host barriers shape the structure of local gene pool and the reservoir of host-associated genetic variants. Our study provides important insights to the genetic and evolutionary factors that contribute to the ability of S. aureus to colonize and cause disease in multiple host species. Our study highlights the importance of continuous surveillance of S. aureus circulating in different ecological host niches and the need to systematically sample from them. These findings will inform development of effective measures to control S. aureus colonization, infection, and transmission across the One Health continuum.

Facultative symbiont virulence determines horizontal transmission rate without host strain specificity

In facultative symbioses, only a fraction of hosts are associated with a symbiont. Understanding why specific host and symbiont strains are associated can inform us of the evolutionary forces affecting facultative symbioses. Possibilities include ongoing host-symbiont coevolution driven by reciprocal selection, or priority effects that are neutral in respect to the host-symbiont interaction. We hypothesized that ongoing host-symbiont coevolution would lead to higher fitness estimates for naturally co-occurring (native) host and symbiont combinations compared to nonnative combinations. We used the Dictyostelium discoideum - Paraburkholderia bonniea system to test this hypothesis. P. bonniea features a reduced genome size relative to another Paraburkholderia symbiont of D. discoideum, indicating a significant history of coevolution with its host. Facultative symbionts may experience continued genome reduction if coevolution is ongoing, or their genome size may have reached a stable state if the symbiosis has also stabilized. Our work demonstrates that ongoing coevolution is unlikely for D. discoideum and P. bonniea. The system instead represents a stable facultative symbiosis. Specifically associated host and symbiont strains in this system are the result of priority effects, and presently unassociated hosts are simply uncolonized. We find evidence for a virulence-transmission trade-off without host strain specificity, and identify candidate virulence factors in the genomes of P. bonniea strains that may contribute to variation in benevolence.

Whole-Genome Subtyping Reveals Population Structure and Host Adaptation of Salmonella Typhimurium from Wild Birds

Within-host evolution of bacterial pathogens can lead to host-associated variants of the same species or serovar. Identification and characterization of closely related variants from diverse host species are crucial to public health and host-pathogen adaptation research. However, the work remained largely underexplored at a strain level until the advent of whole-genome sequencing (WGS). Here, we performed WGS-based subtyping and analyses of Salmonella enterica serovar Typhimurium (n = 787) from different wild birds across 18 countries over a 75-year period. We revealed seven avian host-associated S. Typhimurium variants/lineages. These lineages emerged globally over short timescales and presented genetic features distinct from S. Typhimurium lineages circulating among humans and domestic animals. We further showed that, in terms of virulence, host adaptation of these variants was driven by genome degradation. Our results provide a snapshot of the population structure and genetic diversity of S. Typhimurium within avian hosts. We also demonstrate the value of WGS-based subtyping and analyses in unravelling closely related variants at the strain level.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km² region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Comparative genomics and phylogenomics of Campylobacter unveil potential novel species and provide insights into niche segregation

Campylobacter is a bacterial genus associated with community outbreaks and gastrointestinal symptoms. Studies on Campylobacter generally focus on specific pathogenic species such as C. coli and C. jejuni. Currently, there are thousands of publicly available Campylobacter genomes, allowing a more complete assessment of the genus diversity. In this work, we report a network-based analysis of all available Campylobacter genomes to explore the genus structure and diversity, revealing potentially new species and elucidating genus features. We also hypothesize that the previously established Clade III of C. coli is in fact a novel species (referred here as Campylobacter spp12). Finally, we found a negative correlation between pangenome fluidity and saturation coefficient, with potential implications to the lifestyles of distinct Campylobacter species. Since pangenome analysis depends on the number of available genomes, this correlation could help estimate pangenome metrics of Campylobacter species with less sequenced genomes, helping understand their lifestyle and niche adaptation. Together, our results indicate that the Campylobacter genus should be re-evaluated, with particular attention to the interplay between genome structure and niche segregation.

Genomic epidemiology and transmission characteristics of mcr1-positive colistin-resistant Escherichia coli strains circulating at natural environment

MCR-positive Escherichia coli (MCRPEC) have been reported in humans worldwide. The high prevalence of mcr-1 poses clinical and environmental risks due to its diverse genetic mechanisms. Given the vital role of animals and the environment in the spread of antibiotic resistance, a "One Health" perspective should be taken when addressing antimicrobial resistance issues. This study conducted a prospective study in six farms (located in Jiaxing City, Zhejiang province, China) in 2019. MCRPEC strains were screened from samples of different sources. The molecular epidemiological surveys and transmission potential were investigated by whole-genome sequencing and phylogenetic analysis. MCRPEC were detected in different farms with various sources. Sequence type complex 10 was dominant and distributed widely in multiple sources. Core-genome multilocus sequence type (cgMLST) analysis indicated that clonal transmission could occur within and between farms. In addition, mcr-1 genes with different locations showed different transmission tendencies. The study indicated that interspecies and cross-regional transmission of MCRPEC could occur between different sectors in farms. Further surveillance and research of non-clinical MCRPEC strains are necessary to reduce the threat of MCRPEC.

Antimicrobial Resistance in the Global Health Network: Known Unknowns and Challenges for Efficient Responses in the 21st Century

Antimicrobial resistance (AMR) is one of the Global Health challenges of the 21st century. The inclusion of AMR on the global map parallels the scientific, technological, and organizational progress of the healthcare system and the socioeconomic changes of the last 100 years. Available knowledge about AMR has mostly come from large healthcare institutions in high-income countries and is scattered in studies across various fields, focused on patient safety (infectious diseases), transmission pathways and pathogen reservoirs (molecular epidemiology), the extent of the problem at a population level (public health), their management and cost (health economics), cultural issues (community psychology), and events associated with historical periods (history of science). However, there is little dialogue between the aspects that facilitate the development, spread, and evolution of AMR and various stakeholders (patients, clinicians, public health professionals, scientists, economic sectors, and funding agencies). This study consists of four complementary sections. The first reviews the socioeconomic factors that have contributed to building the current Global Healthcare system, the scientific framework in which AMR has traditionally been approached in such a system, and the novel scientific and organizational challenges of approaching AMR in the fourth globalization scenario. The second discusses the need to reframe AMR in the current public health and global health contexts. Given that the implementation of policies and guidelines are greatly influenced by AMR information from surveillance systems, in the third section, we review the unit of analysis ("the what" and "the who") and the indicators (the "operational units of surveillance") used in AMR and discuss the factors that affect the validity, reliability, and comparability of the information to be applied in various healthcare (primary, secondary, and tertiary), demographic, and economic contexts (local, regional, global, and inter-sectorial levels). Finally, we discuss the disparities and similarities between distinct stakeholders' objectives and the gaps and challenges of combatting AMR at various levels. In summary, this is a comprehensive but not exhaustive revision of the known unknowns about how to analyze the heterogeneities of hosts, microbes, and hospital patches, the role of surrounding ecosystems, and the challenges they represent for surveillance, antimicrobial stewardship, and infection control programs, which are the traditional cornerstones for controlling AMR in human health.

Assessment of plasmids for relating the 2020 Salmonella enterica serovar Newport onion outbreak to farms implicated by the outbreak investigation

BACKGROUND

The Salmonella enterica serovar Newport red onion outbreak of 2020 was the largest foodborne outbreak of Salmonella in over a decade. The epidemiological investigation suggested two farms as the likely source of contamination. However, single nucleotide polymorphism (SNP) analysis of the whole genome sequencing data showed that none of the Salmonella isolates collected from the farm regions were linked to the clinical isolates-preventing the use of phylogenetics in source identification. Here, we explored an alternative method for analyzing the whole genome sequencing data driven by the hypothesis that if the outbreak strain had come from the farm regions, then the clinical isolates would disproportionately contain plasmids found in isolates from the farm regions due to horizontal transfer.

RESULTS

SNP analysis confirmed that the clinical isolates formed a single, nearly-clonal clade with evidence for ancestry in California going back a decade. The clinical clade had a large core genome (4,399 genes) and a large and sparsely distributed accessory genome (2,577 genes, at least 64% on plasmids). At least 20 plasmid types occurred in the clinical clade, more than were found in the literature for Salmonella Newport. A small number of plasmids, 14 from 13 clinical isolates and 17 from 8 farm isolates, were found to be highly similar (> 95% identical)-indicating they might be related by horizontal transfer. Phylogenetic analysis was unable to determine the geographic origin, isolation source, or time of transfer of the plasmids, likely due to their promiscuous and transient nature. However, our resampling analysis suggested that observing a similar number and combination of highly similar plasmids in random samples of environmental Salmonella enterica within the NCBI Pathogen Detection database was unlikely, supporting a connection between the outbreak strain and the farms implicated by the epidemiological investigation.

CONCLUSION

Horizontally transferred plasmids provided evidence for a connection between clinical isolates and the farms implicated as the source of the outbreak. Our case study suggests that such analyses might add a new dimension to source tracking investigations, but highlights the need for detailed and accurate metadata, more extensive environmental sampling, and a better understanding of plasmid molecular evolution.

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.

Staphylococcus aureus host interactions and adaptation

Invasive Staphylococcus aureus infections are common, causing high mortality, compounded by the propensity of the bacterium to develop drug resistance. S. aureus is an excellent case study of the potential for a bacterium to be commensal, colonizing, latent or disease-causing; these states defined by the interplay between S. aureus and host. This interplay is multidimensional and evolving, exemplified by the spread of S. aureus between humans and other animal reservoirs and the lack of success in vaccine development. In this Review, we examine recent advances in understanding the S. aureus-host interactions that lead to infections. We revisit the primary role of neutrophils in controlling infection, summarizing the discovery of new immune evasion molecules and the discovery of new functions ascribed to well-known virulence factors. We explore the intriguing intersection of bacterial and host metabolism, where crosstalk in both directions can influence immune responses and infection outcomes. This Review also assesses the surprising genomic plasticity of S. aureus, its dualism as a multi-mammalian species commensal and opportunistic pathogen and our developing understanding of the roles of other bacteria in shaping S. aureus colonization.

Integrated omic techniques and their genomic features for invasive weeds

Many emerging invasive weeds display rapid adaptation against different stressful environments compared to their natives. Rapid adaptation and dispersal habits helped invasive populations have strong diversity within the population compared to their natives. Advances in molecular marker techniques may lead to an in-depth understanding of the genetic diversity of invasive weeds. The use of molecular techniques is rapidly growing, and their implications in invasive weed studies are considered powerful tools for genome purposes. Here, we review different approach used multi-omics by invasive weed studies to understand the functional structural and genomic changes in these species under different environmental fluctuations, particularly, to check the accessibility of advance-sequencing techniques used by researchers in genome sequence projects. In this review-based study, we also examine the importance and efficiency of different molecular techniques in identifying and characterizing different genes, associated markers, proteins, metabolites, and key metabolic pathways in invasive and native weeds. Use of these techniques could help weed scientists to further reduce the knowledge gaps in understanding invasive weeds traits. Although these techniques can provide robust insights about the molecular functioning, employing a single omics platform can rarely elucidate the gene-level regulation and the associated real-time expression of weedy traits due to the complex and overlapping nature of biological interactions. We conclude that different multi-omic techniques will provide long-term benefits in launching new genome projects to enhance the understanding of invasive weeds' invasion process.

PRAWNS: compact pan-genomic features for whole-genome population genomics

MOTIVATION

Scientists seeking to understand the genomic basis of bacterial phenotypes, such as antibiotic resistance, today have access to an unprecedented number of complete and nearly complete genomes. Making sense of these data requires computational tools able to perform multiple-genome comparisons efficiently, yet currently available tools cannot scale beyond several tens of genomes.

RESULTS

We describe PRAWNS, an efficient and scalable tool for multiple-genome analysis. PRAWNS defines a concise set of genomic features (metablocks), as well as pairwise relationships between them, which can be used as a basis for large-scale genotype-phenotype association studies. We demonstrate the effectiveness of PRAWNS by identifying genomic regions associated with antibiotic resistance in Acinetobacter baumannii.

AVAILABILITY AND IMPLEMENTATION

PRAWNS is implemented in C++ and Python3, licensed under the GPLv3 license, and freely downloadable from GitHub (https://github.com/KiranJavkar/PRAWNS.git).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparative Genomics of the Genus Pseudomonas Reveals Host- and Environment-Specific Evolution

Each Earth ecosystem has unique microbial communities. Pseudomonas bacteria have evolved to occupy a plethora of different ecological niches, including living hosts, such as animals and plants. Many genes necessary for the Pseudomonas-niche interaction and their encoded functions remain unknown. Here, we describe a comparative genomic study of 3,274 genomes with 19,056,667 protein-coding sequences from Pseudomonas strains isolated from diverse environments. We detected functional divergence of Pseudomonas that depends on the niche. Each group of strains from a certain environment harbored a distinctive set of metabolic pathways or functions. The horizontal transfer of genes, which mainly proceeded between closely related taxa, was dependent on the isolation source. Finally, we detected thousands of undescribed proteins and functions associated with each Pseudomonas lifestyle. This research represents an effort to reveal the mechanisms underlying the ecology, pathogenicity, and evolution of Pseudomonas, and it will enable clinical, ecological, and biotechnological advances. IMPORTANCE Microbes play important roles in the health of living beings and in the environment. The knowledge of these functions may be useful for the development of new clinical and biotechnological applications and the restoration and preservation of natural ecosystems. However, most mechanisms implicated in the interaction of microbes with the environment remain poorly understood; thus, this field of research is very important. Here, we try to understand the mechanisms that facilitate the differential adaptation of Pseudomonas-a large and ubiquitous bacterial genus-to the environment. We analyzed more than 3,000 Pseudomonas genomes and searched for genetic patterns that can be related with their coevolution with different hosts (animals, plants, or fungi) and environments. Our results revealed that thousands of genes and genetic features are associated with each niche. Our data may be useful to develop new technical and theoretical advances in the fields of ecology, health, and industry.

How to survive pig farming: Mechanism of SCCmec element deletion and metabolic stress adaptation in livestock-associated MRSA

Previous research on methicillin susceptible Staphylococcus aureus (MSSA) belonging to livestock-associated (LA-) sequence type (ST) 398, isolated from pigs and their local surroundings, indicated that differences between these MSSA and their methicillin resistant predecessors (MRSA) are often limited to the absence of the staphylococcal cassette chromosome mec (SCCmec) and few single nucleotide polymorphisms. So far, our understanding on how LA-MRSA endure the environmental conditions associated with pig-farming as well as the putative impact of this particular environment on the mobilisation of SCCmec elements is limited. Thus, we performed in-depth genomic and transcriptomic analyses using the LA-MRSA ST398 strain IMT38951 and its methicillin susceptible descendant. We identified a mosaic-structured SCCmec region including a putative replicative SCCmecVc which is absent from the MSSA chromosome through homologous recombination. Based on our data, such events occur between short repetitive sequences identified within and adjacent to two distinct alleles of the large cassette recombinase genes C (ccrC). We further evaluated the global transcriptomic response of MRSA ST398 to particular pig-farm associated conditions, i.e., contact with host proteins (porcine serum) and a high ammonia concentration. Differential expression of global regulators involved in stress response control were identified, i.e., ammonia-induced alternative sigma factor B-depending activation of genes for the alkaline shock protein 23, the heat shock response and the accessory gene regulator (agr)-controlled transcription of virulence factors. Exposure to serum transiently induced the transcription of distinct virulence factor encoding genes. Transcription of genes reported for mediating the loss of methicillin resistance, especially ccrC, was not significantly different compared to the unchallenged controls. We concluded that, from an evolutionary perspective, bacteria may save energy by incidentally dismissing a fully replicative SCCmec element in contrast to the induction of ccr genes on a population scale. Since the genomic SCCmec integration site is a hot-spot of recombination, occasional losses of elements of 16 kb size may restore capacities for the uptake of foreign genetic material. Subsequent spread of resistance, on the other hand, might depend on the autonomous replication machinery of the deleted SCCmec elements that probably enhance chances for reintegration of SCCmec into susceptible genomes by mere multiplication.

Genetic diversity and variation in antimicrobial-resistance determinants of non-serotype 2 Streptococcus suis isolates from healthy pigs

Streptococcus suis is a leading cause of bacterial meningitis in South-East Asia, with frequent zoonotic transfer to humans associated with close contact with pigs. A small number of invasive lineages are responsible for endemic infection in the swine industry, causing considerable global economic losses. A lack of surveillance and a rising trend in clinical treatment failure has raised concerns of growing antimicrobial resistance (AMR) among invasive S. suis. Gene flow between healthy and disease isolates is poorly understood and, in this study, we sample and sequence a collection of isolates predominantly from healthy pigs in Chiang Mai province, Northern Thailand. Pangenome characterization identified extensive genetic diversity and frequent AMR carriage in isolates from healthy pigs. Multiple AMR genes were identified, conferring resistance to aminoglycosides, lincosamides, tetracycline and macrolides. All isolates were non-susceptible to three or more different antimicrobial classes, and 75 % of non-serotype 2 isolates were non-susceptible to six or more classes (compared to 37.5 % of serotype 2 isolates). AMR genes were found on integrative and conjugative elements previously observed in other species, suggesting a mobile gene pool that can be accessed by invasive disease isolates. This article contains data hosted by Microreact.

Functional analysis of intergenic regulatory regions of genes encoding surface adhesins in Staphylococcus aureus isolates from periprosthetic joint infections

Staphylococcus aureus is a leading cause of prosthetic joint infections (PJI). Surface adhesins play an important role in the primary attachment to plasma proteins that coat the surface of prosthetic devices after implantation. Previous efforts to identify a genetic component of the bacterium that confers an enhanced capacity to cause PJI have focused on gene content, kmers, or single-nucleotide polymorphisms (SNPs) in coding sequences. Here, using a collection of S. aureus strains isolated from PJI and wounds, we investigated whether genetic variations in the regulatory region of genes encoding surface adhesins lead to differences in their expression levels and modulate the capacity of S. aureus to colonize implanted prosthetic devices. The data revealed that S. aureus isolates from the same clonal complex (CC) contain a specific pattern of SNPs in the regulatory region of genes encoding surface adhesins. As a consequence, each clonal lineage shows a specific profile of surface proteins expression. Co-infection experiments with representative isolates of the most prevalent CCs demonstrated that some lineages have a higher capacity to colonize implanted catheters in a murine infection model, which correlated with a greater ability to form a biofilm on coated surfaces with plasma proteins. Together, results indicate that differences in the expression level of surface adhesins may modulate the propensity of S. aureus strains to cause PJI. Given the high conservation of surface proteins among staphylococci, our work lays the framework for investigating how diversification at intergenic regulatory regions affects the capacity of S. aureus to colonize the surface of medical implants.

Distinct evolutionary trajectories in the Escherichia coli pangenome occur within sequence types

The Escherichia coli species contains a diverse set of sequence types and there remain important questions regarding differences in genetic content within this population that need to be addressed. Pangenomes are useful vehicles for studying gene content within sequence types. Here, we analyse 21 E. coli sequence type pangenomes using comparative pangenomics to identify variance in both pangenome structure and content. We present functional breakdowns of sequence type core genomes and identify sequence types that are enriched in metabolism, transcription and cell membrane biogenesis genes. We also uncover metabolism genes that have variable core classification, depending on which allele is present. Our comparative pangenomics approach allows for detailed exploration of sequence type pangenomes within the context of the species. We show that ongoing gene gain and loss in the E. coli pangenome is sequence type-specific, which may be a consequence of distinct sequence type-specific evolutionary drivers.

Horizontal transfer and phylogenetic distribution of the immune evasion factor tarP

Methicillin-resistant Staphylococcus aureus (MRSA), a major human pathogen, uses the prophage-encoded tarP gene as an important immune evasion factor. TarP glycosylates wall teichoic acid (WTA) polymers, major S. aureus surface antigens, to impair WTA immunogenicity and impede host defence. However, tarP phages appear to be restricted to only a few MRSA clonal lineages, including clonal complexes (CC) 5 and 398, for unknown reasons. We demonstrate here that tarP-encoding prophages can be mobilized to lysogenize other S. aureus strains. However, transfer is largely restricted to closely related clones. Most of the non-transducible clones encode tarM, which generates a WTA glycosylation pattern distinct from that mediated by TarP. However, tarM does not interfere with infection by tarP phages. Clonal complex-specific Type I restriction-modification systems were the major reasons for resistance to tarP phage infection. Nevertheless, tarP phages were found also in unrelated S. aureus clones indicating that tarP has the potential to spread to distant clonal lineages and contribute to the evolution of new MRSA clones.

Lactic Acid Resistance and Population Structure of Escherichia coli from Meat Processing Environment

To explore the effect of beef processing on Escherichia coli populations in relation to lactic acid resistance, this study investigated the links among acid response, phylogenetic structure, genome diversity, and genotypes associated with acid resistance of meat plant E. coli. Generic E. coli isolates (n = 700) were from carcasses, fabrication equipment, and beef products. Acid treatment was carried out in Luria-Bertani broth containing 5.5% lactic acid (pH 2.9). Log reductions of E. coli ranged from <0.5 to >5 log CFU/mL (median: 1.37 log). No difference in lactic acid resistance was observed between E. coli populations recovered before and after a processing step or antimicrobial interventions. E. coli from the preintervention carcasses were slightly more resistant than E. coli isolated from equipment, differing by <0.5 log unit. Acid-resistant E. coli (log reduction <1, n = 45) had a higher prevalence of genes related to energy metabolism (ydj, xap, ato) and oxidative stress (fec, ymjC) than the less resistant E. coli (log reduction >1, n = 133). The ydj and ato operons were abundant in E. coli from preintervention carcasses. In contrast, fec genes were abundant in E. coli from equipment surfaces. The preintervention E. coli contained phylogroups A and B1 in relatively equal proportions. Phylogroup B1 predominated (95%) in the population from equipment. Of note, E. coli collected after sanitation shared either the antigens of O8 or H21. Additionally, genome diversity decreased after chilling and equipment sanitation. Overall, beef processing did not select for E. coli resistant to lactic acid but shaped the population structure. IMPORTANCE Antimicrobial interventions have significantly reduced the microbial loads on carcasses/meat products; however, the wide use of chemical and physical biocides has raised concerns over their potential for selecting resistant populations in the beef processing environment. Phenotyping of acid resistance and whole-genome analysis described in this study demonstrated beef processing practices led to differences in acid resistance, genotype, and population structure between carcass- and equipment-associated E. coli but did not select for the acid-resistant population. Results indicate that genes coding for the metabolism of long-chain sugar acids (ydj) and short-chain fatty acids (ato) were more prevalent in carcass-associated than equipment-associated E. coli. These results suggest E. coli from carcasses and equipment surfaces have been exposed to different selective pressures. The findings improve our understanding of the microbial ecology of E. coli in food processing environments and in general.

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