Producing polished prokaryotic pangenomes with the Panaroo pipeline

FepR as a Central Genetic Target in the Adaptation to Quaternary Ammonium Compounds and Cross-Resistance to Ciprofloxacin in Listeria monocytogenes

The foodborne pathogen, Listeria monocytogenes, (Lm), frequently undergoes selection pressure associated with the extensive use of disinfectants, such as quaternary ammonium compounds, which are widely used in food processing plants. The repeated exposure to sub-inhibitory biocide concentrations can induce increased tolerance to these compounds, but can also trigger the development of antibiotic resistance, and both increase the risk of food contamination and persistence in food production environments. Although the acquisition of genes can explain biocide tolerance, the genetic mechanisms underlying the adaptive cross-resistance to antibiotics remain unclear. We previously showed that repeated exposure to benzalkonium chloride (BC) and didecyldimethyl ammonium chloride (DDAC) led to reduced susceptibility to ciprofloxacin in Lm strains from diverse sources. Here, we compared the genomes of 16 biocide-adapted and 10 parental strains to identify the molecular mechanisms of fluoroquinolone cross-resistance. A core genome SNP analysis identified various mutations in the transcriptional regulator fepR (lmo2088) for 94% of the adapted strains and mutations in other effectors at a lower frequency. FepR is a local repressor of the MATE fluoroquinolone efflux pump FepA. The impact of the mutations on the structure and function of the protein was assessed by performing in silico prediction and protein homology modeling. Our results show that 75% of the missense mutations observed in fepR are located in the HTH domain of the protein, within the DNA interaction site. These mutations are predicted to reduce the activity of the regulator, leading to the overexpression of the efflux pump responsible for the ciprofloxacin-enhanced resistance.

A curated collection of Klebsiella metabolic models reveals variable substrate usage and gene essentiality

The Klebsiella pneumoniae species complex (KpSC) is a set of seven Klebsiella taxa that are found in a variety of niches and are an important cause of opportunistic health care-associated infections in humans. Because of increasing rates of multi-drug resistance within the KpSC, there is a growing interest in better understanding the biology and metabolism of these organisms to inform novel control strategies. We collated 37 sequenced KpSC isolates isolated from a variety of niches, representing all seven taxa. We generated strain-specific genome-scale metabolic models (GEMs) for all 37 isolates and simulated growth phenotypes on 511 distinct carbon, nitrogen, sulfur, and phosphorus substrates. Models were curated and their accuracy was assessed using matched phenotypic growth data for 94 substrates (median accuracy of 96%). We explored species-specific growth capabilities and examined the impact of all possible single gene deletions using growth simulations in 145 core carbon substrates. These analyses revealed multiple strain-specific differences, within and between species, and highlight the importance of selecting a diverse range of strains when exploring KpSC metabolism. This diverse set of highly accurate GEMs could be used to inform novel drug design, enhance genomic analyses, and identify novel virulence and resistance determinants. We envisage that these 37 curated strain-specific GEMs, covering all seven taxa of the KpSC, provide a valuable resource to the Klebsiella research community.

Accessory Genomes Drive Independent Spread of Carbapenem-Resistant Klebsiella pneumoniae Clonal Groups 258 and 307 in Houston, TX

Carbapenem-resistant Klebsiella pneumoniae (CRKp) is an urgent public health threat. Worldwide dissemination of CRKp has been largely attributed to clonal group (CG) 258. However, recent evidence indicates the global emergence of a CRKp CG307 lineage. Houston, TX, is the first large city in the United States with detected cocirculation of both CRKp CG307 and CG258. We sought to characterize the genomic and clinical factors contributing to the parallel endemic spread of CG258 and CG307. CRKp isolates were collected as part of the prospective, Consortium on Resistance against Carbapenems in Klebsiella and other Enterobacterales 2 (CRACKLE-2) study. Hybrid short-read and long-read genome assemblies were generated from 119 CRKp isolates (95 originated from Houston hospitals). A comprehensive characterization of phylogenies, gene transfer, and plasmid content with pan-genome analysis was performed on all CRKp isolates. Plasmid mating experiments were performed with CG307 and CG258 isolates of interest. Dissection of the accessory genomes suggested independent evolution and limited horizontal gene transfer between CG307 and CG258 lineages. CG307 contained a diverse repertoire of mobile genetic elements, which were shared with other non-CG258 K. pneumoniae isolates. Three unique clades of Houston CG307 isolates clustered distinctly from other global CG307 isolates, indicating potential selective adaptation of particular CG307 lineages to their respective geographical niches. CG307 strains were often isolated from the urine of hospitalized patients, likely serving as important reservoirs for genes encoding carbapenemases and extended-spectrum β-lactamases. Our findings suggest parallel cocirculation of high-risk lineages with potentially divergent evolution. IMPORTANCE The prevalence of carbapenem-resistant Klebsiella pneumoniae (CRKp) infections in nosocomial settings remains a public health challenge. High-risk clones such as clonal group 258 (CG258) are particularly concerning due to their association with blaKPC carriage, which can severely complicate antimicrobial treatments. There is a recent emergence of clonal group 307 (CG307) worldwide with little understanding of how this successful clone has been able to adapt while cocirculating with CG258. We provide the first evidence of potentially divergent evolution between CG258 and CG307 with limited sharing of adaptive genes. Houston, TX, is home to the largest medical center in the world, with a large influx of domestic and international patients. Thus, our unique geographical setting, where two pandemic strains of CRKp are circulating, provides an indication of how differential accessory genome content can drive stable, endemic populations of CRKp. Pan-genomic analyses such as these can reveal unique signatures of successful CRKp dissemination, such as the CG307-associated plasmid (pCG307_HTX), and provide invaluable insights into the surveillance of local carbapenem-resistant Enterobacterales (CRE) epidemiology.

Linear plasmids in Klebsiella and other Enterobacteriaceae

Linear plasmids are extrachromosomal DNA elements that have been found in a small number of bacterial species. To date, the only linear plasmids described in the family Enterobacteriaceae belong to Salmonella, first found in Salmonella enterica Typhi. Here, we describe a collection of 12 isolates of the Klebsiella pneumoniae species complex in which we identified linear plasmids. Screening of assembly graphs assembled from public read sets identified linear plasmid structures in a further 13 K. pneumoniae species complex genomes. We used these 25 linear plasmid sequences to query all bacterial genome assemblies in the National Center for Biotechnology Information database, and discovered an additional 61 linear plasmid sequences in a variety of Enterobacteriaceae species. Gene content analysis divided these plasmids into five distinct phylogroups, with very few genes shared across more than two phylogroups. The majority of linear plasmid-encoded genes are of unknown function; however, each phylogroup carried its own unique toxin–antitoxin system and genes with homology to those encoding the ParAB plasmid stability system. Passage in vitro of the 12 linear plasmid-carrying Klebsiella isolates in our collection (which include representatives of all five phylogroups) indicated that these linear plasmids can be stably maintained, and our data suggest they can transmit between K. pneumoniae strains (including members of globally disseminated multidrug-resistant clones) and also between diverse Enterobacteriaceae species. The linear plasmid sequences, and representative isolates harbouring them, are made available as a resource to facilitate future studies on the evolution and function of these novel plasmids.

Pangenome Evolution Reconciles Robustness and Instability of Rhizobial Symbiosis

Root nodulating rhizobia are nearly ubiquitous in soils and provide the critical service of nitrogen fixation to thousands of legume species, including staple crops. However, the magnitude of fixed nitrogen provided to hosts varies markedly among rhizobia strains, despite host legumes having mechanisms to selectively reward beneficial strains and to punish ones that do not fix sufficient nitrogen. Variation in the services of microbial mutualists is considered paradoxical given host mechanisms to select beneficial genotypes. Moreover, the recurrent evolution of non-fixing symbiont genotypes is predicted to destabilize symbiosis, but breakdown has rarely been observed. Here, we deconstructed hundreds of genome sequences from genotypically and phenotypically diverse Bradyrhizobium strains and revealed mechanisms that generate variation in symbiotic nitrogen fixation. We show that this trait is conferred by a modular system consisting of many extremely large integrative conjugative elements and few conjugative plasmids. Their transmissibility and propensity to reshuffle genes generate new combinations that lead to uncooperative genotypes and make individual partnerships unstable. We also demonstrate that these same properties extend beneficial associations to diverse host species and transfer symbiotic capacity among diverse strains. Hence, symbiotic nitrogen fixation is underpinned by modularity, which engenders flexibility, a feature that reconciles evolutionary robustness and instability. These results provide new insights into mechanisms driving the evolution of mobile genetic elements. Moreover, they yield a new predictive model on the evolution of rhizobial symbioses, one that informs on the health of organisms and ecosystems that are hosts to symbionts and that helps resolve the long-standing paradox.

Comparative Genomics of Disease and Carriage Serotype 1 Pneumococci

The isolation of Streptococcus pneumoniae serotypes in systemic tissues of patients with invasive disease versus the nasopharynx of healthy individuals with asymptomatic carriage varies widely. Some serotypes are hyper-invasive, particularly serotype 1, but the underlying genetics remain poorly understood due to the rarity of carriage isolates, reducing the power of comparison with invasive isolates. Here, we use a well-controlled genome-wide association study to search for genetic variation associated with invasiveness of serotype 1 pneumococci from a serotype 1 endemic setting in Africa. We found no consensus evidence that certain genomic variation is overrepresented among isolates from patients with invasive disease than asymptomatic carriage. Overall, the genomic variation explained negligible phenotypic variability, suggesting a minimal effect on the disease status. Furthermore, changes in lineage distribution were seen with lineages replacing each other over time, highlighting the importance of continued pathogen surveillance. Our findings suggest that the hyper-invasiveness is an intrinsic property of the serotype 1 strains, not specific for a "disease-associated" subpopulation disproportionately harboring unique genomic variation.

Panakeia - a universal tool for bacterial pangenome analysis

Background

Development of new pan-genome analysis tools is important, as the pangenome of a microbial species has become an important method to define the diversity of a selected taxon, most commonly a species, in the last years. This enables comparison of strains from different ecological niches and can be used to define the functional potential in a bacterial population. It gives us a much better view of microbial genomics than can be gained from singular genomes which after all are just single representatives of a much more varied population.

Results

We present Panakeia, a tool which strives to be easy to use and providing a detailed view of the pangenome structure which can efficiently be utilised for discovery, or further in-depth analysis, of features of interest. It analyses synteny and multiple structural patterns of the pangenome, giving insights into the biological diversity and evolution of the studied taxon. Panakeia hence provides both broad and detailed information on the structure of a pangenome, for diverse and highly clonal populations of bacteria.

Conclusions

Previously published pangenome tools often reduce the information to a presence/absence matrix of unconnected genes or generate massive hard to interpret output graphs. However, Panakeia includes synteny and structural information and presents it in a way that can readily be used for further analysis. Panakeia can be downloaded at https://github.com/BioSina/Panakeiatogether with a detailed User Guide.

Horizontal transmission enables flexible associations with locally adapted symbiont strains in deep-sea hydrothermal vent symbioses

SignificanceIn marine ecosystems, transmission of microbial symbionts between host generations occurs predominantly through the environment. Yet, it remains largely unknown how host genetics, symbiont competition, environmental conditions, and geography shape the composition of symbionts acquired by individual hosts. To address this question, we applied population genomic approaches to four species of deep-sea hydrothermal vent snails that live in association with chemosynthetic bacteria. Our analyses show that environment is more important to strain-level symbiont composition than host genetics and that symbiont strains show genetic variation indicative of adaptation to the distinct geochemical conditions at each vent site. This corroborates a long-standing hypothesis that hydrothermal vent invertebrates affiliate with locally adapted symbiont strains to cope with the variable conditions characterizing their habitats.

PowerBacGWAS: a computational pipeline to perform power calculations for bacterial genome-wide association studies

Genome-wide association studies (GWAS) are increasingly being applied to investigate the genetic basis of bacterial traits. However, approaches to perform power calculations for bacterial GWAS are limited. Here we implemented two alternative approaches to conduct power calculations using existing collections of bacterial genomes. First, a sub-sampling approach was undertaken to reduce the allele frequency and effect size of a known and detectable genotype-phenotype relationship by modifying phenotype labels. Second, a phenotype-simulation approach was conducted to simulate phenotypes from existing genetic variants. We implemented both approaches into a computational pipeline (PowerBacGWAS) that supports power calculations for burden testing, pan-genome and variant GWAS; and applied it to collections of Enterococcus faecium, Klebsiella pneumoniae and Mycobacterium tuberculosis. We used this pipeline to determine sample sizes required to detect causal variants of different minor allele frequencies (MAF), effect sizes and phenotype heritability, and studied the effect of homoplasy and population diversity on the power to detect causal variants. Our pipeline and user documentation are made available and can be applied to other bacterial populations. PowerBacGWAS can be used to determine sample sizes required to find statistically significant associations, or the associations detectable with a given sample size. We recommend to perform power calculations using existing genomes of the bacterial species and population of study.

Global evolutionary dynamics and resistome analysis of Clostridioides difficile ribotype 017

Clostridioides difficile PCR ribotype (RT) 017 ranks among the most successful strains of C. difficile in the world. In the past three decades, it has caused outbreaks on four continents, more than other ‘epidemic’ strains, but our understanding of the genomic epidemiology underpinning the spread of C. difficile RT 017 is limited. Here, we performed high-resolution phylogenomic and Bayesian evolutionary analyses on an updated and more representative dataset of 282 non-clonal C. difficile RT 017 isolates collected worldwide between 1981 and 2019. These analyses place an estimated time of global dissemination between 1953 and 1983 and identified the acquisition of the ermB-positive transposon Tn6194 as a key factor behind global emergence. This coincided with the introduction of clindamycin, a key inciter of C. difficile infection, into clinical practice in the 1960s. Based on the genomic data alone, the origin of C. difficile RT 017 could not be determined; however, geographical data and records of population movement suggest that C. difficile RT 017 had been moving between Asia and Europe since the Middle Ages and was later transported to North America around 1860 (95 % confidence interval: 1622–1954). A focused epidemiological study of 45 clinical C. difficile RT 017 genomes from a cluster in a tertiary hospital in Thailand revealed that the population consisted of two groups of multidrug-resistant (MDR) C. difficile RT 017 and a group of early, non-MDR C. difficile RT 017. The significant genomic diversity within each MDR group suggests that although they were all isolated from hospitalized patients, there was probably a reservoir of C. difficile RT 017 in the community that contributed to the spread of this pathogen.

Accessory Genome Dynamics of Local and Global Staphylococcus pseudintermedius Populations

Staphylococcus pseudintermedius is a major bacterial colonizer and opportunistic pathogen in dogs. Methicillin-resistant S. pseudintermedius (MRSP) continues to emerge as a significant challenge to maintaining canine health. We sought to determine the phylogenetic relationships of S. pseudintermedius across five states in the New England region of the United States and place them in a global context. The New England dataset consisted of 125 previously published S. pseudintermedius genomes supplemented with 45 newly sequenced isolates. The core genome phylogenetic tree revealed many deep branching lineages consisting of 142 multi-locus sequence types (STs). In silico detection of the mecA gene revealed 40 MRSP and 130 methicillin-susceptible S. pseudintermedius (MSSP) isolates. MRSP were derived from five structural types of SCCmec, the mobile genetic element that carries the mecA gene conferring methicillin resistance. Although many genomes were MSSP, they nevertheless harbored genes conferring resistance to many other antibiotic classes, including aminoglycosides, macrolides, tetracyclines and penams. We compared the New England genomes to 297 previously published genomes sampled from five other states in the United States and 13 other countries. Despite the prevalence of the clonally expanding ST71 found worldwide and in other parts of the United States, we did not detect it in New England. We next sought to interrogate the combined New England and global datasets for the presence of coincident gene pairs linked to antibiotic resistance. Analysis revealed a large co-circulating accessory gene cluster, which included mecA as well as eight other resistance genes [aac (6′)-Ie-aph (2″)-Ia, aad (6), aph (3′)-IIIa, sat4, ermB, cat, blaZ, and tetM]. Furthermore, MRSP isolates carried significantly more accessory genes than their MSSP counterparts. Our results provide important insights to the evolution and geographic spread of high-risk clones that can threaten the health of our canine companions.

Genome wide association study of Escherichia coli bloodstream infection isolates identifies genetic determinants for the portal of entry but not fatal outcome

Escherichia coli is an important cause of bloodstream infections (BSI), which is of concern given its high mortality and increasing worldwide prevalence. Finding bacterial genetic variants that might contribute to patient death is of interest to better understand infection progression and implement diagnostic methods that specifically look for those factors. E. coli samples isolated from patients with BSI are an ideal dataset to systematically search for those variants, as long as the influence of host factors such as comorbidities are taken into account. Here we performed a genome-wide association study (GWAS) using data from 912 patients with E. coli BSI from hospitals in Paris, France. We looked for associations between bacterial genetic variants and three patient outcomes (death at 28 days, septic shock and admission to intensive care unit), as well as two portals of entry (urinary and digestive tract), using various clinical variables from each patient to account for host factors. We did not find any association between genetic variants and patient outcomes, potentially confirming the strong influence of host factors in influencing the course of BSI; we however found a strong association between the papGII operon and entrance of E. coli through the urinary tract, which demonstrates the power of bacterial GWAS when applied to actual clinical data. Despite the lack of associations between E. coli genetic variants and patient outcomes, we estimate that increasing the sample size by one order of magnitude could lead to the discovery of some putative causal variants. Given the wide adoption of bacterial genome sequencing of clinical isolates, such sample sizes may be soon available.

Pangenomic and functional investigations for dormancy and biodegradation features of an organic pollutant-degrading bacterium Rhodococcus biphenylivorans TG9

Environmental bacteria contain a wealth of untapped potential in the form of biodegradative genes. Leveraging this potential can often be confounded by a lack of understanding of fundamental survival strategies, like dormancy, for environmental stress. Investigating bacterial dormancy-to-degradation relationships enables improvement of bioremediation. Here, we couple genomic and functional assessment to provide context for key attributes of the organic pollutant-degrading strain Rhodococcus biphenylivorans TG9. Whole genome sequencing, pangenome analysis and functional characterization were performed to elucidate important genes and gene products, including antimicrobial resistance, dormancy, and degradation. Rhodococcus as a genus has strong potential for degradation and dormancy, which we demonstrate using R. biphenylivorans TG9 as a model. We identified four Resuscitation-promoting factor (Rpf) encoding genes in TG9 involved in dormancy and resuscitation. We demonstrate that R. biphenylivorans TG9 grows on fourteen typical organic pollutants, and exhibits a robust ability to degrade biphenyl and several congeners of polychlorinated biphenyls. We further induced TG9 into a dormant state and demonstrated pronounced differences in morphology and activity. Together, these results expand our understanding of the genus Rhodococcus and the relationship between dormancy and biodegradation in the presence of environmental stressors.

Oxytetracycline and Streptomycin Resistance Genes in Xanthomonas arboricola pv. pruni, the Causal Agent of Bacterial Spot in Peach

Xanthomonas arboricola pv. pruni (Xap) causes bacterial spot, a major worldwide disease of Prunus species. Very few chemical management options are available for this disease and frequent applications of oxytetracycline (OTC) in the United States peach orchards have raised concerns about resistance development. During 2017-2020, 430 Xap strains were collected from ten peach orchards in South Carolina. Seven OTC-resistant (OTC R ) Xap strains were found in 2017 and 2020 from four orchards about 20-270 km apart. Interestingly, the seven strains were also resistant to streptomycin (STR). Six strains grew on media amended with ≤100 μg/mL OTC, while one strain, R1, grew on ≤250 μg/mL OTC. Genome sequence analysis of four representative OTC R strains revealed a 14-20 kb plasmid carrying tetC, tetR, and strAB in each strain. These three genes were transferable to Xanthomonas perforans via conjugation, and they were PCR confirmed in all seven OTC R Xap strains. When tetC and tetR were cloned and expressed together in a sensitive strain, the transconjugants showed resistance to ≤100 μg/mL OTC. When tetC was cloned and expressed alone in a sensitive strain, the transconjugants showed resistance to ≤250 μg/mL OTC. TetC and tetR expression was inducible by OTC in all six wild-type strains resistant to ≤100 μg/mL OTC. However, in the R1 strain resistant to ≤250 μg/mL OTC, tetR was not expressed, possibly due to the presence of Tn3 in the tetR gene, and in this case tetC was constitutively expressed. These data suggest that tetC confers OTC resistance in Xap strains, and tetR regulates the level of OTC resistance conferred by tetC. To our knowledge, this is the first report of OTC resistance in plant pathogenic xanthomonads.

Genome-wide association, prediction and heritability in bacteria with application to Streptococcus pneumoniae

Whole-genome sequencing has facilitated genome-wide analyses of association, prediction and heritability in many organisms. However, such analyses in bacteria are still in their infancy, being limited by difficulties including genome plasticity and strong population structure. Here we propose a suite of methods including linear mixed models, elastic net and LD-score regression, adapted to bacterial traits using innovations such as frequency-based allele coding, both insertion/deletion and nucleotide testing and heritability partitioning. We compare and validate our methods against the current state-of-art using simulations, and analyse three phenotypes of the major human pathogen Streptococcus pneumoniae, including the first analyses of minimum inhibitory concentrations (MIC) for penicillin and ceftriaxone. We show that the MIC traits are highly heritable with high prediction accuracy, explained by many genetic associations under good population structure control. In ceftriaxone MIC, this is surprising because none of the isolates are resistant as per the inhibition zone criteria. We estimate that half of the heritability of penicillin MIC is explained by a known drug-resistance region, which also contributes a quarter of the ceftriaxone MIC heritability. For the within-host carriage duration phenotype, no associations were observed, but the moderate heritability and prediction accuracy indicate a moderately polygenic trait.

Identification of novel, cryptic Clostridioides species isolates from environmental samples collected from diverse geographical locations

Clostridioides difficile is a pathogen often associated with hospital-acquired infection or antimicrobial-induced disease; however, increasing evidence indicates infections can result from community or environmental sources. Most genomic sequencing of C. difficile has focused on clinical strains, although evidence is growing that C. difficile spores are widespread in soil and water in the environment. In this study, we sequenced 38 genomes collected from soil and water isolates in Flagstaff (AZ, USA) and Slovenia in an effort targeted towards environmental surveillance of C. difficile. At the average nucleotide identity (ANI) level, the genomes were divergent to C. difficile at a threshold consistent with different species. A phylogenetic analysis of these divergent genomes together with Clostridioides genomes available in public repositories confirmed the presence of three previously described, cryptic Clostridioides species and added two additional clades. One of the cryptic species (C-III) was almost entirely composed of Arizona and Slovenia genomes, and contained distinct sub-groups from each region (evidenced by SNP and gene-content differences). A comparative genomics analysis identified multiple unique coding sequences per clade, which can serve as markers for subsequent environmental surveys of these cryptic species. Homologues to the C. difficile toxin genes, tcdA and tcdB, were found in cryptic species genomes, although they were not part of the typical pathogenicity locus observed in C. difficile, and in silico PCR suggested that some would not amplify with widely used PCR diagnostic tests. We also identified gene homologues in the binary toxin cluster, including some present on phage and, for what is believed to be the first time, on a plasmid. All isolates were obtained from environmental samples, so the function and disease potential of these toxin homologues is currently unknown. Enzymatic profiles of a subset of cryptic isolates (n=5) demonstrated differences, suggesting that these isolates contain substantial metabolic diversity. Antimicrobial resistance (AMR) was observed across a subset of isolates (n=4), suggesting that AMR mechanisms are intrinsic to the genus, perhaps originating from a shared environmental origin. This study greatly expands our understanding of the genomic diversity of Clostridioides. These results have implications for C. difficile One Health research, for more sensitive C. difficile diagnostics, as well as for understanding the evolutionary history of C. difficile and the development of pathogenesis.

Streptococcus dysgalactiae subspecies dysgalactiae in Norwegian bovine dairy herds: Risk factors, sources, and genomic diversity

Despite the importance of Streptococcus dysgalactiae ssp. dysgalactiae (SDSD) as an udder pathogen, the reservoir and epidemiological characteristics of this bacterium are largely unexplored. The aims of this study were to investigate risk factors for SDSD intramammary infections (SDSD-IMI) in Norwegian bovine dairy herds, identify sources of SDSD on animals and in the environment, and elucidate the genetic diversity of SDSD isolates. Data from herd recordings and a questionnaire were used to investigate herd-level risk factors for SDSD-IMI in 359 freestall dairy herds. Seven herds with a suspected high prevalence of SDSD-IMI were visited to sample extramammary sources (e.g., skin, wounds, mucous membranes, and freestall environment). Bacterial isolates were whole-genome sequenced to investigate the distribution of SDSD genotypes within herds and to assess the phylogenetic relationship between SDSD isolates from 27 herds across Norway. Risk factors for high incidence of SDSD-IMI in freestall dairy herds were related to housing, including closed flooring in alleys and rubber mats in cubicle bases. Parlor milking was also a risk factor compared with automatic milking systems. From herd visits, a considerable proportion of extramammary samples were SDSD positive, particularly from wounds and skin of the animals and the cubicle bases. Samples from mucous surfaces (nostrils, rectum, and vagina) and water troughs were least frequently positive. Eight multilocus sequence types (ST) were identified among the sequenced isolates from 27 herds, and phylogenetic analyses revealed 8 clades corresponding to ST. No significant association was identified between sampling site (milk, body sites, and environment) and ST. In 4 of 6 herds from which 5 or more isolates were available, one ST dominated and was found in milk and extramammary samples. One ST (ST453) was found in 15 of 27 herds, which implies that this is a widely distributed and possibly a bovine-adapted strain. Findings in this study suggest that SDSD is a cow-adapted opportunist with potential for contagious transmission, and that the freestall environment is likely to play a role in transmission between cows.

International links between Streptococcus pneumoniae vaccine serotype 4 sequence type (ST) 801 in Northern European shipyard outbreaks of invasive pneumococcal disease

BACKGROUND

Pneumococcal disease outbreaks of vaccine preventable serotype 4 sequence type (ST)801 in shipyards have been reported in several countries. We aimed to use genomics to establish any international links between them.

METHODS

Sequence data from ST801-related outbreak isolates from Norway (n = 17), Finland (n = 11) and Northern Ireland (n = 2) were combined with invasive pneumococcal disease surveillance from the respective countries, and ST801-related genomes from an international collection (n = 41 of > 40,000), totalling 106 genomes. Raw data were mapped and recombination excluded before phylogenetic dating.

RESULTS

Outbreak isolates were relatively diverse, with up to 100 SNPs (single nucleotide polymorphisms) and a common ancestor estimated around the year 2000. However, 19 Norwegian and Finnish isolates were nearly indistinguishable (0-2 SNPs) with the common ancestor dated around 2017.

CONCLUSION

The total diversity of ST801 within the outbreaks could not be explained by recent transmission alone, suggesting that harsh environmental and associated living conditions reported in the shipyards may facilitate invasion of colonising pneumococci. However, near identical strains in the Norwegian and Finnish outbreaks does suggest that transmission between international shipyards also contributed to those outbreaks. This indicates the need for improved preventative measures in this working population including pneumococcal vaccination.

Curtobacterium spp. and Curtobacterium flaccumfaciens: Phylogeny, Genomics-Based Taxonomy, Pathogenicity, and Diagnostics

The genus of Curtobacterium, belonging to the Microbacteriaceae family of the Actinomycetales order, includes economically significant pathogenic bacteria of soybeans and other agricultural crops. Thorough phylogenetic and full-genome analysis using the latest genomic data has demonstrated a complex and contradictory taxonomic picture within the group of organisms classified as the Curtobacterium species. Based on these data, it is possible to delineate about 50 new species and to reclassify a substantial part of the Curtobacterium strains. It is suggested that 53 strains, including most of the Curtobacterium flaccumfaciens pathovars, can compose a monophyletic group classified as C. flaccumfaciens. A genomic analysis using the most recent inventory of bacterial chromosomal and plasmid genomes deposited to GenBank confirmed the possible role of Microbacteriaceae plasmids in pathogenicity and demonstrated the existence of a group of related plasmids carrying virulence factors and possessing a gene distantly related to DNA polymerase found in bacteriophages and archaeal and eukaryotic viruses. A PCR diagnostic assay specific to the genus Curtobacterium was developed and tested. The presented results assist in the understanding of the evolutionary relations within the genus and can lay the foundation for further taxonomic updates.

Changes in the skin microbiome associated with squamous cell carcinoma in transplant recipients

Actinic keratoses (AK) arise in severely photo-damaged skin and can progress to squamous cell carcinomas (SCC). AK and SCC are common in Caucasian populations, and immunosuppressed individuals have a markedly higher risk of developing SCC. An overabundance of Staphylococcus aureus has been reported in AK and SCC lesions of immunocompetent individuals, however, the AK/SCC microbiome in immunosuppressed cohorts has not been investigated. Here, the microbial profile and bacterial load of AK, SCC and control skin swabs from 32 immunosuppressed organ transplant recipients were characterised via SSU rRNA gene sequencing and qPCR, and compared to a previously described immunocompetent cohort. Although the taxonomic composition of skin swab samples was mostly subject-specific, significant differences were observed between control skin, AK, and SCC in both cohorts. Surface bacterial load was increased and alpha diversity decreased in AK and SCC compared to control skin due to an increased abundance of Staphylococcus species and relative decrease of skin commensals. Staphylococcus epidermidis predominated on SCC from transplant recipients in contrast to SCC of immunocompetent subjects dominated by S. aureus. In conclusion, AK and SCC of immunosuppressed and immunocompetent subjects present with distinctive microbial dysbioses, which may be relevant to SCC pathogenesis and progression.

Multiple phylogenetically-diverse, differentially-virulent Burkholderia pseudomallei isolated from a single soil sample collected in Thailand

Burkholderia pseudomallei is a soil-dwelling bacterium endemic to Southeast Asia and northern Australia that causes the disease, melioidosis. Although the global genomic diversity of clinical B. pseudomallei isolates has been investigated, there is limited understanding of its genomic diversity across small geographic scales, especially in soil. In this study, we obtained 288 B. pseudomallei isolates from a single soil sample (~100g; intensive site 2, INT2) collected at a depth of 30cm from a site in Ubon Ratchathani Province, Thailand. We sequenced the genomes of 169 of these isolates that represent 7 distinct sequence types (STs), including a new ST (ST1820), based on multi-locus sequence typing (MLST) analysis. A core genome SNP phylogeny demonstrated that all identified STs share a recent common ancestor that diverged an estimated 796-1260 years ago. A pan-genomics analysis demonstrated recombination between clades and intra-MLST phylogenetic and gene differences. To identify potential differential virulence between STs, groups of BALB/c mice (5 mice/isolate) were challenged via subcutaneous injection (500 CFUs) with 30 INT2 isolates representing 5 different STs; over the 21-day experiment, eight isolates killed all mice, 2 isolates killed an intermediate number of mice (1-2), and 20 isolates killed no mice. Although the virulence results were largely stratified by ST, one virulent isolate and six attenuated isolates were from the same ST (ST1005), suggesting that variably conserved genomic regions may contribute to virulence. Genomes from the animal-challenged isolates were subjected to a bacterial genome-wide association study to identify genomic regions associated with differential virulence. One associated region is a unique variant of Hcp1, a component of the type VI secretion system, which may result in attenuation. The results of this study have implications for comprehensive sampling strategies, environmental exposure risk assessment, and understanding recombination and differential virulence in B. pseudomallei.

Functional strain redundancy and persistent phage infection in Swiss hard cheese starter cultures

Undefined starter cultures are poorly characterized bacterial communities from environmental origin used in cheese making. They are phenotypically stable and have evolved through domestication by repeated propagation in closed and highly controlled environments over centuries. This makes them interesting for understanding eco-evolutionary dynamics governing microbial communities. While cheese starter cultures are known to be dominated by a few bacterial species, little is known about the composition, functional relevance, and temporal dynamics of strain-level diversity. Here, we applied shotgun metagenomics to an important Swiss cheese starter culture and analyzed historical and experimental samples reflecting 82 years of starter culture propagation. We found that the bacterial community is highly stable and dominated by only a few coexisting strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. lactis. Genome sequencing, metabolomics analysis, and co-culturing experiments of 43 isolates show that these strains are functionally redundant, but differ tremendously in their phage resistance potential. Moreover, we identified two highly abundant Streptococcus phages that seem to stably coexist in the community without any negative impact on bacterial growth or strain persistence, and despite the presence of a large and diverse repertoire of matching CRISPR spacers. Our findings show that functionally equivalent strains can coexist in domesticated microbial communities and highlight an important role of bacteria-phage interactions that are different from kill-the-winner dynamics.

Functional strain redundancy and persistent phage infection in Swiss hard cheese starter cultures

Undefined starter cultures are poorly characterized bacterial communities from environmental origin used in cheese making. They are phenotypically stable and have evolved through domestication by repeated propagation in closed and highly controlled environments over centuries. This makes them interesting for understanding eco-evolutionary dynamics governing microbial communities. While cheese starter cultures are known to be dominated by a few bacterial species, little is known about the composition, functional relevance, and temporal dynamics of strain-level diversity. Here, we applied shotgun metagenomics to an important Swiss cheese starter culture and analyzed historical and experimental samples reflecting 82 years of starter culture propagation. We found that the bacterial community is highly stable and dominated by only a few coexisting strains of Streptococcus thermophilus and Lactobacillus delbrueckii subsp. lactis. Genome sequencing, metabolomics analysis, and co-culturing experiments of 43 isolates show that these strains are functionally redundant, but differ tremendously in their phage resistance potential. Moreover, we identified two highly abundant Streptococcus phages that seem to stably coexist in the community without any negative impact on bacterial growth or strain persistence, and despite the presence of a large and diverse repertoire of matching CRISPR spacers. Our findings show that functionally equivalent strains can coexist in domesticated microbial communities and highlight an important role of bacteria-phage interactions that are different from kill-the-winner dynamics.

Beyond to the Stable: Role of the Insertion Sequences as Epidemiological Descriptors in Corynebacterium striatum

In recent years, epidemiological studies of infectious agents have focused mainly on the pathogen and stable components of its genome. The use of these stable components makes it possible to know the evolutionary or epidemiological relationships of the isolates of a particular pathogen. Under this approach, focused on the pathogen, the identification of resistance genes is a complementary stage of a bacterial characterization process or an appendix of its epidemiological characterization, neglecting its genetic components’ acquisition or dispersal mechanisms. Today we know that a large part of antibiotic resistance is associated with mobile elements. Corynebacterium striatum, a bacterium from the normal skin microbiota, is also an opportunistic pathogen. In recent years, reports of infections and nosocomial outbreaks caused by antimicrobial multidrug-resistant C. striatum strains have been increasing worldwide. Despite the different existing mobile genomic elements, there is evidence that acquired resistance genes are coupled to insertion sequences in C. striatum. This perspective article reviews the insertion sequences linked to resistance genes, their relationship to evolutionary lineages, epidemiological characteristics, and the niches the strains inhabit. Finally, we evaluate the potential of the insertion sequences for their application as a descriptor of epidemiological scenarios, allowing us to anticipate the emergence of multidrug-resistant lineages.

Expanding the Burkholderia pseudomallei Complex with the Addition of Two Novel Species: Burkholderia mayonis sp. nov. and Burkholderia savannae sp. nov

Distinct Burkholderia strains were isolated from soil samples collected in tropical northern Australia (Northern Territory and the Torres Strait Islands, Queensland). Phylogenetic analysis of 16S rRNA and whole genome sequences revealed these strains were distinct from previously described Burkholderia species and assigned them to two novel clades within the B. pseudomallei complex (Bpc). Because average nucleotide identity and digital DNA-DNA hybridization calculations are consistent with these clades representing distinct species, we propose the names Burkholderia mayonis sp. nov. and Burkholderia savannae sp. nov. Strains assigned to B. mayonis sp. nov. include type strain BDU6T (=TSD-80; LMG 29941; ASM152374v2) and BDU8. Strains assigned to B. savannae sp. nov. include type strain MSMB266T (=TSD-82; LMG 29940; ASM152444v2), MSMB852, BDU18, and BDU19. Comparative genomics revealed unique coding regions for both putative species, including clusters of orthologous genes associated with phage. Type strains of both B. mayonis sp. nov. and B. savannae sp. nov. yielded biochemical profiles distinct from each other and from other species in the Bpc, and profiles also varied among strains within B. mayonis sp. nov. and B. savannae sp. nov. Matrix-assisted laser desorption ionization time-of-flight (MLST) analysis revealed a B. savannae sp. nov. cluster separate from other species, whereas B. mayonis sp. nov. strains did not form a distinct cluster. Neither B. mayonis sp. nov. nor B. savannae sp. nov. caused mortality in mice when delivered via the subcutaneous route. The addition of B. mayonis sp. nov. and B. savannae sp. nov. results in a total of eight species currently within the Bpc. IMPORTANCEBurkholderia species can be important sources of novel natural products, and new species are of interest to diverse scientific disciplines. Although many Burkholderia species are saprophytic, Burkholderia pseudomallei is the causative agent of the disease melioidosis. Understanding the genomics and virulence of the closest relatives to B. pseudomallei, i.e., the other species within the B. pseudomallei complex (Bpc), is important for identifying robust diagnostic targets specific to B. pseudomallei and for understanding the evolution of virulence in B. pseudomallei. Two proposed novel species, B. mayonis sp. nov. and B. savannae sp. nov., were isolated from soil samples collected from multiple locations in northern Australia. The two proposed species belong to the Bpc but are phylogenetically distinct from all other members of this complex. The addition of B. mayonis sp. nov. and B. savannae sp. nov. results in a total of eight species within this significant complex of bacteria that are available for future studies.

Development of a Genomics-Based Approach To Identify Putative Hypervirulent Nontyphoidal Salmonella Isolates: Salmonella enterica Serovar Saintpaul as a Model

While differences in human virulence have been reported across nontyphoidal Salmonella (NTS) serovars and associated subtypes, a rational and scalable approach to identify Salmonella subtypes with differential ability to cause human diseases is not available. Here, we used NTS serovar Saintpaul (S. Saintpaul) as a model to determine if metadata and associated whole-genome sequence (WGS) data in the NCBI Pathogen Detection (PD) database can be used to identify (i) subtypes with differential likelihoods of causing human diseases and (ii) genes and single nucleotide polymorphisms (SNPs) potentially responsible for such differences. S. Saintpaul SNP clusters (n = 211) were assigned different epidemiology types (epi-types) based on statistically significant over- or underrepresentation of human clinical isolates, including human associated (HA; n = 29), non-human associated (NHA; n = 23), and other (n = 159). Comparative genomic analyses identified 384 and 619 genes overrepresented among isolates in 5 HA and 4 NHA SNP clusters most significantly associated with the respective isolation source. These genes included 5 HA-associated virulence genes previously reported to be present on Gifsy-1/Gifsy-2 prophages. Additionally, premature stop codons in 3 and 7 genes were overrepresented among the selected HA and NHA SNP clusters, respectively. Tissue culture experiments with strains representing 4 HA and 3 NHA SNP clusters did not reveal evidence for enhanced invasion or intracellular survival for HA strains. However, the presence of sodCI (encoding a superoxide dismutase), found in 4 HA and 1 NHA SNP clusters, was positively correlated with intracellular survival in macrophage-like cells. Post hoc analyses also suggested a possible difference in intracellular survival among S. Saintpaul lineages. IMPORTANCE Not all Salmonella isolates are equally likely to cause human disease, and Salmonella control strategies may unintentionally focus on serovars and subtypes with high prevalence in source populations but are rarely associated with human clinical illness. We describe a framework leveraging WGS data in the NCBI PD database to identify Salmonella subtypes over- and underrepresented among human clinical cases. While we identified genomic signatures associated with HA/NHA SNP clusters, tissue culture experiments failed to identify consistent phenotypic characteristics indicative of enhanced human virulence of HA strains. Our findings illustrate the challenges of defining hypo- and hypervirulent S. Saintpaul and potential limitations of phenotypic assays when evaluating human virulence, for which in vivo experiments are essential. Identification of sodCI, an HA-associated virulence gene associated with enhanced intracellular survival, however, illustrates the potential of the framework and is consistent with prior work identifying specific genomic features responsible for enhanced or reduced virulence of nontyphoidal Salmonella.

The NagY regulator: A member of the BglG/SacY antiterminator family conserved in Enterococcus faecalis and involved in virulence

Enterococcus faecalis is a commensal bacterium of the gastrointestinal tract but also a major nosocomial pathogen. This bacterium uses regulators like BglG/SacY family of transcriptional antiterminators to adapt its metabolism during host colonization. In this report, we investigated the role of the BglG/SacY family antiterminator NagY in the regulation of the nagY-nagE operon in presence of N-acetylglucosamine, with nagE encoding a transporter of this carbohydrate, as well as the expression of the virulence factor HylA. We showed that this last protein is involved in biofilm formation and glycosaminoglycans degradation that are important features in bacterial infection, confirmed in the Galleria mellonella model. In order to elucidate the evolution of these actors, we performed phylogenomic analyses on E. faecalis and Enterococcaceae genomes, identified orthologous sequences of NagY, NagE, and HylA, and we report their taxonomic distribution. The study of the conservation of the upstream region of nagY and hylA genes showed that the molecular mechanism of NagY regulation involves ribonucleic antiterminator sequence overlapping a rho-independent terminator, suggesting a regulation conforming to the canonical model of BglG/SacY family antiterminators. In the perspective of opportunism understanding, we offer new insights into the mechanism of host sensing thanks to the NagY antiterminator and its targets expression.

Pneumococcal within-host diversity during colonization, transmission and treatment

Characterizing the genetic diversity of pathogens within the host promises to greatly improve surveillance and reconstruction of transmission chains. For bacteria, it also informs our understanding of inter-strain competition and how this shapes the distribution of resistant and sensitive bacteria. Here we study the genetic diversity of Streptococcus pneumoniae within 468 infants and 145 of their mothers by deep sequencing whole pneumococcal populations from 3,761 longitudinal nasopharyngeal samples. We demonstrate that deep sequencing has unsurpassed sensitivity for detecting multiple colonization, doubling the rate at which highly invasive serotype 1 bacteria were detected in carriage compared with gold-standard methods. The greater resolution identified an elevated rate of transmission from mothers to their children in the first year of the child's life. Comprehensive treatment data demonstrated that infants were at an elevated risk of both the acquisition and persistent colonization of a multidrug-resistant bacterium following antimicrobial treatment. Some alleles were enriched after antimicrobial treatment, suggesting that they aided persistence, but generally purifying selection dominated within-host evolution. Rates of co-colonization imply that in the absence of treatment, susceptible lineages outcompeted resistant lineages within the host. These results demonstrate the many benefits of deep sequencing for the genomic surveillance of bacterial pathogens.

Genomic diversity of Helicobacter pylori populations from different regions of the human stomach

Individuals infected with Helicobacter pylori harbor unique and diverse populations of quasispecies, but diversity between and within different regions of the human stomach and the process of bacterial adaptation to each location are not yet well understood. We applied whole-genome deep sequencing to characterize the within- and between-stomach region genetic diversity of H. pylori populations from paired antrum and corpus biopsies of 15 patients, along with single biopsies from one region of an additional 3 patients, by scanning allelic diversity. We combined population deep sequencing with more conventional sequencing of multiple H. pylori single colony isolates from individual biopsies to generate a unique dataset. Single colony isolates were used to validate the scanning allelic diversity pipelines. We detected extensive population allelic diversity within the different regions of each patient's stomach. Diversity was most commonly found within non-coding, hypothetical, outer membrane, restriction modification system, virulence, lipopolysaccharide biosynthesis, efflux systems, and chemotaxis-associated genes. Antrum and corpus populations from the same patient grouped together phylogenetically, indicating that most patients were initially infected with a single strain, which then diversified. Single colonies from the antrum and corpus of the same patients grouped into distinct clades, suggesting mechanisms for within-location adaptation across multiple H. pylori isolates from different patients. The comparisons made available by combined sequencing and analysis of isolates and populations enabled comprehensive analysis of the genetic changes associated with H. pylori diversification and stomach region adaptation.

Integrating cultivation and metagenomics for a multi-kingdom view of skin microbiome diversity and functions

Human skin functions as a physical barrier to foreign pathogen invasion and houses numerous commensals. Shifts in the human skin microbiome have been associated with conditions ranging from acne to atopic dermatitis. Previous metagenomic investigations into the role of the skin microbiome in health or disease have found that much of the sequenced data do not match reference genomes, making it difficult to interpret metagenomic datasets. We combined bacterial cultivation and metagenomic sequencing to assemble the Skin Microbial Genome Collection (SMGC), which comprises 622 prokaryotic species derived from 7,535 metagenome-assembled genomes and 251 isolate genomes. The metagenomic datasets that we generated were combined with publicly available skin metagenomic datasets to identify members and functions of the human skin microbiome. The SMGC collection includes 174 newly identified bacterial species and 12 newly identified bacterial genera, including the abundant genus 'Candidatus Pellibacterium', which has been newly associated with the skin. The SMGC increases the characterized set of known skin bacteria by 26%. We validated the SMGC metagenome-assembled genomes by comparing them with sequenced isolates obtained from the same samples. We also recovered 12 eukaryotic species and assembled thousands of viral sequences, including newly identified clades of jumbo phages. The SMGC enables classification of a median of 85% of skin metagenomic sequences and provides a comprehensive view of skin microbiome diversity, derived primarily from samples obtained in North America.

Applied Bioinformatics and Public Health Microbiology: challenges, discoveries and innovations during a pandemic

The eighth Applied Bioinformatics and Public Health Microbiology (ABPHM) conference showcased the recent acceleration of bioinformatic approaches used in public health settings. This included approaches for the surveillance of infectious diseases, understanding microbial evolution and diversity and pathogen interactions. Overall, the meeting highlighted the importance of data-driven approaches used by scientists during the COVID-19 pandemic.

A catalogue of 1,167 genomes from the human gut archaeome

The human gut microbiome plays an important role in health, but its archaeal diversity remains largely unexplored. In the present study, we report the analysis of 1,167 nonredundant archaeal genomes (608 high-quality genomes) recovered from human gastrointestinal tract, sampled across 24 countries and rural and urban populations. We identified previously undescribed taxa including 3 genera, 15 species and 52 strains. Based on distinct genomic features, we justify the split of the Methanobrevibacter smithii clade into two separate species, with one represented by the previously undescribed 'Candidatus Methanobrevibacter intestini'. Patterns derived from 28,581 protein clusters showed significant associations with sociodemographic characteristics such as age groups and lifestyle. We additionally show that archaea are characterized by specific genomic and functional adaptations to the host and carry a complex virome. Our work expands our current understanding of the human archaeome and provides a large genome catalogue for future analyses to decipher its impact on human physiology.

Characterization of the Type I Restriction Modification System Broadly Conserved among Group A Streptococci

Although prokaryotic DNA methylation investigations have long focused on immunity against exogenous DNA, it has been recently recognized that DNA methylation impacts gene expression and phase variation in Streptococcus pneumoniae and Streptococcus suis. A comprehensive analysis of DNA methylation is lacking for beta-hemolytic streptococci, and thus we sought to examine DNA methylation in the major human pathogen group A Streptococcus (GAS). Using a database of 224 GAS genomes encompassing 80 emm types, we found that nearly all GAS strains encode a type I restriction modification (RM) system that lacks the hsdS′ alleles responsible for impacting gene expression in S. pneumoniae and S. suis. The GAS type I system is located on the core chromosome, while sporadically present type II orphan methyltransferases were identified on prophages. By combining single-molecule real-time (SMRT) analyses of 10 distinct emm types along with phylogenomics of 224 strains, we were able to assign 13 methylation patterns to the GAS population. Inactivation of the type I RM system, occurring either naturally through phage insertion or through laboratory-induced gene deletion, abrogated DNA methylation detectable via either SMRT or MinION sequencing. Contrary to a previous report, inactivation of the type I system did not impact transcript levels of the gene (mga) encoding the key multigene activator protein (Mga) or Mga-regulated genes. Inactivation of the type I system significantly increased plasmid transformation rates. These data delineate the breadth of the core chromosomal type I RM system in the GAS population and clarify its role in immunity rather than impacting Mga regulon expression.

IMPORTANCE The advent of whole-genome approaches capable of detecting DNA methylation has markedly expanded appreciation of the diverse roles of epigenetic modification in prokaryotic physiology. For example, recent studies have suggested that DNA methylation impacts gene expression in some streptococci. The data described herein are from the first systematic analysis of DNA methylation in a beta-hemolytic streptococcus and one of the few analyses to comprehensively characterize DNA methylation across hundreds of strains of the same bacterial species. We clarify that DNA methylation in group A Streptococcus (GAS) is primarily due to a type I restriction modification (RM) system present in the core genome and does not impact mga-regulated virulence gene expression, but does impact immunity against exogenous DNA. The identification of the DNA motifs recognized by each type I RM system may assist with optimizing methods for GAS genetic manipulation and help us understand how bacterial pathogens acquire exogenous DNA elements.

The Mouse Gastrointestinal Bacteria Catalogue enables translation between the mouse and human gut microbiotas via functional mapping

Human health and disease have increasingly been shown to be impacted by the gut microbiota, and mouse models are essential for investigating these effects. However, the compositions of human and mouse gut microbiotas are distinct, limiting translation of microbiota research between these hosts. To address this, we constructed the Mouse Gastrointestinal Bacteria Catalogue (MGBC), a repository of 26,640 high-quality mouse microbiota-derived bacterial genomes. This catalog enables species-level analyses for mapping functions of interest and identifying functionally equivalent taxa between the microbiotas of humans and mice. We have complemented this with a publicly deposited collection of 223 bacterial isolates, including 62 previously uncultured species, to facilitate experimental investigation of individual commensal bacteria functions in vitro and in vivo. Together, these resources provide the ability to identify and test functionally equivalent members of the host-specific gut microbiotas of humans and mice and support the informed use of mouse models in human microbiota research.

Protocol for post-processing of bacterial pangenome data using Pagoo pipeline

Multiple downstream analyses are necessary to interpret the output of bacterial pangenome reconstruction software. This requires integrating diverse kinds of genetic and phenotypic data, which to date are left to each user's criterion. To fill this gap, we created Pagoo, a pangenome post-processing tool that leverages a standardized but flexible and extensible framework for data integration, analysis, and storage. Here, we provide the protocol for running Pagoo and performing from simple to more complex comparative analyses on bacterial pangenome data. For complete details on the use and execution of this protocol, please refer to Ferrés and Iraola (2021).

Introduction and adaptation of an emerging pathogen to olive trees in Italy

The invasive plant pathogen Xylella fastidiosa currently threatens European flora through the loss of economically and culturally important host plants. This emerging vector-borne bacterium, native to the Americas, causes several important diseases in a wide range of plants including crops, ornamentals, and trees. Previously absent from Europe, and considered a quarantine pathogen, X. fastidiosa was first detected in Apulia, Italy in 2013 associated with a devastating disease of olive trees (Olive Quick Decline Syndrome, OQDS). OQDS has led to significant economic, environmental, cultural, as well as political crises. Although the biology of X. fastidiosa diseases have been studied for over a century, there is still no information on the determinants of specificity between bacterial genotypes and host plant species, which is particularly relevant today as X. fastidiosa is expanding in the naive European landscape. We analysed the genomes of 79 X. fastidiosa samples from diseased olive trees across the affected area in Italy as well as genomes of the most genetically closely related strains from Central America. We provided insights into the ecological and evolutionary emergence of this pathogen in Italy. We first showed that the outbreak in Apulia is due to a single introduction from Central America that we estimated to have occurred in 2008 [95 % HPD: 1930–2016]. By using a combination of population genomic approaches and evolutionary genomics methods, we further identified a short list of genes that could play a major role in the adaptation of X. fastidiosa to this new environment. We finally provided experimental evidence for the adaptation of the strain to this new environment.

A high-throughput multiplexing and selection strategy to complete bacterial genomes

BACKGROUND

Bacterial whole-genome sequencing based on short-read technologies often results in a draft assembly formed by contiguous sequences. The introduction of long-read sequencing technologies permits those contiguous sequences to be unambiguously bridged into complete genomes. However, the elevated costs associated with long-read sequencing frequently limit the number of bacterial isolates that can be long-read sequenced. Here we evaluated the recently released 96 barcoding kit from Oxford Nanopore Technologies (ONT) to generate complete genomes on a high-throughput basis. In addition, we propose an isolate selection strategy that optimizes a representative selection of isolates for long-read sequencing considering as input large-scale bacterial collections.

RESULTS

Despite an uneven distribution of long reads per barcode, near-complete chromosomal sequences (assembly contiguity = 0.89) were generated for 96 Escherichia coli isolates with associated short-read sequencing data. The assembly contiguity of the plasmid replicons was even higher (0.98), which indicated the suitability of the multiplexing strategy for studies focused on resolving plasmid sequences. We benchmarked hybrid and ONT-only assemblies and showed that the combination of ONT sequencing data with short-read sequencing data is still highly desirable (i) to perform an unbiased selection of isolates for long-read sequencing, (ii) to achieve an optimal genome accuracy and completeness, and (iii) to include small plasmids underrepresented in the ONT library.

CONCLUSIONS

The proposed long-read isolate selection ensures the completion of bacterial genomes that span the genome diversity inherent in large collections of bacterial isolates. We show the potential of using this multiplexing approach to close bacterial genomes on a high-throughput basis.

First Comparative Analysis of Clostridium septicum Genomes Provides Insights Into the Taxonomy, Species Genetic Diversity, and Virulence Related to Gas Gangrene

Clostridium septicum is a Gram-positive, toxin-producing, and spore-forming bacterium that is recognized, together with C. perfringens, as the most important etiologic agent of progressive gas gangrene. Clostridium septicum infections are almost always fatal in humans and animals. Despite its clinical and agricultural relevance, there is currently limited knowledge of the diversity and genome structure of C. septicum. This study presents the complete genome sequence of C. septicum DSM 7534T type strain as well as the first comparative analysis of five C. septicum genomes. The taxonomy of C. septicum, as revealed by 16S rRNA analysis as well as by genomic wide indices such as protein-based phylogeny, average nucleotide identity, and digital DNA–DNA hybridization indicates a stable clade. The composition and presence of prophages, CRISPR elements and accessory genetic material was variable in the investigated genomes. This is in contrast to the limited genetic variability described for the phylogenetically and phenotypically related species Clostridium chauvoei. The restriction-modification (RM) systems between two C. septicum genomes were heterogeneous for the RM types they encoded. C. septicum has an open pangenome with 2,311 genes representing the core genes and 1,429 accessory genes. The core genome SNP divergence between genome pairs varied up to 4,886 pairwise SNPs. A vast arsenal of potential virulence genes was detected in the genomes studied. Sequence analysis of these genes revealed that sialidase, hemolysin, and collagenase genes are conserved compared to the α-toxin and hyaluronidase genes. In addition, a conserved gene found in all C. septicum genomes was predicted to encode a leucocidin homolog (beta-channel forming cytolysin) similar (71.10% protein identity) to Clostridium chauvoei toxin A (CctA), which is a potent toxin. In conclusion, our results provide first, valuable insights into strain relatedness and genomic plasticity of C. septicum and contribute to our understanding of the virulence mechanisms of this important human and animal pathogen.

Phylogenetic Relatedness and Genome Structure of Yersinia ruckeri Revealed by Whole Genome Sequencing and a Comparative Analysis

Yersinia ruckeri is the causative agent of enteric redmouth disease (ERM), a serious infection that affects global aquaculture with high economic impact. The present study used whole genome sequences to perform a comparative analysis on 10 Y. ruckeri strains and to explore their genetic relatedness to other members of the genus. Y. ruckeri, Yersinia entomophaga, and Yersinia nurmii formed a species complex that constitutes the most basal lineage of the genus. The results showed that the taxonomy of Y. ruckeri strains is better defined by using a core genome alignment and phylogenetic analysis. The distribution of accessory genes in all Yersinia species revealed the presence of 303 distinctive genes in Y. ruckeri. Of these, 169 genes were distributed in 17 genomic islands potentially involved in the pathogenesis of ERM via (1) encoding virulence factors such as Afp18, Yrp1, phage proteins and (2) improving the metabolic capabilities by enhancing utilization and metabolism of iron, amino acids (specifically, arginine and histidine), and carbohydrates. The genome of Y. ruckeri is highly conserved regarding gene structure, gene layout and functional categorization of genes. It contains various components of mobile genetic elements but lacks the CRISPR-Cas system and possesses a stable set of virulence genes possibly playing a critical role in pathogenicity. Distinct virulence plasmids were exclusively restricted to a specific clonal group of Y. ruckeri (CG4), possibly indicating a selective advantage. Phylogenetic analysis of Y. ruckeri genomes revealed the co-presence of multiple genetically distant lineages of Y. ruckeri strains circulating in Germany. Our results also suggest a possible dissemination of a specific group of strains in the United States, Peru, Germany, and Denmark. In conclusion, this study provides new insights into the taxonomy and evolution of Y. ruckeri and contributes to a better understanding of the pathogenicity of ERM in aquaculture. The genomic analysis presented here offers a framework for the development of more efficient control strategies for this pathogen.

Genome-wide association study of gastric cancer- and duodenal ulcer-derived Helicobacter pylori strains reveals discriminatory genetic variations and novel oncoprotein candidates

Genome-wide association studies (GWASs) can reveal genetic variations associated with a phenotype in the absence of any hypothesis of candidate genes. The problem of false-positive sites linked with the responsible site might be bypassed in bacteria with a high homologous recombination rate, such as Helicobacter pylori, which causes gastric cancer. We conducted a small-sample GWAS (125 gastric cancer cases and 115 controls) followed by prediction of gastric cancer and control (duodenal ulcer) H. pylori strains. We identified 11 single nucleotide polymorphisms (eight amino acid changes) and three DNA motifs that, combined, allowed effective disease discrimination. They were often informative of the underlying molecular mechanisms, such as electric charge alteration at the ligand-binding pocket, alteration in subunit interaction, and mode-switching of DNA methylation. We also identified three novel virulence factors/oncoprotein candidates. These results provide both defined targets for further informatic and experimental analyses to gain insights into gastric cancer pathogenesis and a basis for identifying a set of biomarkers for distinguishing these H. pylori-related diseases.

Context-aware genomic surveillance reveals hidden transmission of a carbapenemase-producing Klebsiella pneumoniae

Genomic surveillance can inform effective public health responses to pathogen outbreaks. However, integration of non-local data is rarely done. We investigate two large hospital outbreaks of a carbapenemase-carrying Klebsiella pneumoniae strain in Germany and show the value of contextual data. By screening about 10 000 genomes, over 400 000 metagenomes and two culture collections using in silico and in vitro methods, we identify a total of 415 closely related genomes reported in 28 studies. We identify the relationship between the two outbreaks through time-dated phylogeny, including their respective origin. One of the outbreaks presents extensive hidden transmission, with descendant isolates only identified in other studies. We then leverage the genome collection from this meta-analysis to identify genes under positive selection. We thereby identify an inner membrane transporter (ynjC) with a putative role in colistin resistance. Contextual data from other sources can thus enhance local genomic surveillance at multiple levels and should be integrated by default when available.

Phylogenetic Analysis and Genome-Wide Association Study Applied to an Italian Listeria monocytogenes Outbreak

From May 2015 to March 2016, a severe outbreak due to Listeria monocytogenes ST7 strain occurred in Central Italy and caused 24 confirmed clinical cases. The epidemic strain was deeply investigated using whole-genome sequencing (WGS) analysis. In the interested area, the foodborne outbreak investigation identified a meat food-producing plant contaminated by the outbreak strain, carried by pork-ready-to-eat products. In the same region, in March 2018, the epidemic strain reemerged causing one listeriosis case in a 10-month-old child. The aim of this study was to investigate the phylogeny of the epidemic and reemergent strains over time and to compare them with a closer ST7 clone, detected during the outbreak and with different pulsed-field gel electrophoresis (PFGE) profiles, in order to identify genomic features linked to the persistence and the reemergence of the outbreak. An approach combining phylogenetic analysis and genome-wide association study (GWAS) revealed that the epidemic and reemergent clones were genetically closer to the ST7 clone with different PFGE profiles and strictly associated with the pork production chain. The repeated detection of both clones was probably correlated with (i) the presence of truly persistent clones and the repeated introduction of new ones and (ii) the contribution of prophage genes in promoting the persistence of the epidemic clones. Despite that no significant genomic differences were detected between the outbreak and the reemergent strain, the two related clones detected during the outbreak can be differentiated by transcriptional factor and phage genes associated with the phage LP-114.

High Genomic Identity between Clinical and Environmental Strains of Herbaspirillum frisingense Suggests Pre-Adaptation to Different Hosts and Intrinsic Resistance to Multiple Drugs

The genus Herbaspirillum is widely studied for its ability to associate with grasses and to perform biological nitrogen fixation. However, the bacteria of the Herbaspirillum genus have frequently been isolated from clinical samples. Understanding the genomic characteristics that allow these bacteria to switch environments and become able to colonize human hosts is essential for monitoring emerging pathogens and predicting outbreaks. In this work, we describe the sequencing, assembly, and annotation of the genome of H. frisingense AU14559 isolated from the sputum of patients with cystic fibrosis, and its comparison with the genomes of the uropathogenic strain VT-16-41 and the environmental strains GSF30, BH-1, IAC152, and SG826. The genes responsible for biological nitrogen fixation were absent from all strains except for GSF30. On the other hand, genes encoding virulence and host interaction factors were mostly shared with environmental strains. We also identified a large set of intrinsic antibiotic resistance genes that were shared across all strains. Unlike other strains, in addition to unique genomic islands, AU14559 has a mutation that renders the biosynthesis of rhamnose and its incorporation into the exopolysaccharide unfeasible. These data suggest that H. frisingense has characteristics that provide it with the metabolic diversity needed to infect and colonize human hosts.

Genomic Epidemiology and Antimicrobial Resistance Mechanisms of Imported Typhoid in Australia

Typhoid fever is an invasive bacterial disease of humans that disproportionately affects low- and middle-income countries. Antimicrobial resistance (AMR) has been increasingly prevalent in recent decades in Salmonella enterica serovar Typhi (S. Typhi), the causative agent of typhoid fever, limiting treatment options. In Australia, most cases of typhoid fever are imported due to travel to regions where typhoid fever is endemic. Here, all 116 isolates of S. Typhi isolated in Victoria, Australia, between 1 July 2018 and 30 June 2020, underwent whole-genome sequencing and antimicrobial susceptibility testing. Genomic data were linked to international travel data collected from routine case interviews. Travel to South Asia accounted for most cases, with 92.2% imported from seven primary countries (the top two were India, n = 87, and Pakistan, n = 12). A total of 17 S. Typhi genotypes were detected in the 2-year cohort, with 48.2% genotyped as part of global AMR lineages. Ciprofloxacin resistance was detected in two lineages, 3.3 and 4.3.1.2, all from cases with reported travel to India. Nearly all multidrug and extensively drug resistant isolates (90%) were from cases with reported travel to Pakistan in genotypes 4.3.1.1 and 4.3.1.1.P1. Extended spectrum beta-lactamases, blaCTX-M-15 and blaSHV-12, were detected in cases with travel to Pakistan and India, respectively. Linking epidemiological data with genomic studies of S. Typhi provides an opportunity to improve understanding of the emergence, spread and risk of drug-resistant S. Typhi infections and to better inform empirical treatment guidelines in returned travelers.

Genomic Epidemiology and Antimicrobial Susceptibility Profile of Enterotoxigenic Escherichia coli From Outpatients With Diarrhea in Shenzhen, China, 2015-2020

Enterotoxigenic Escherichia coli (ETEC) is the leading cause of severe diarrhea in children and the most common cause of diarrhea in travelers. However, most ETEC infections in Shenzhen, China were from indigenous adults. In this study, we characterized 106 ETEC isolates from indigenous outpatients with diarrhea (77% were adults aged >20 years) in Shenzhen between 2015 and 2020 by whole-genome sequencing and antimicrobial susceptibility testing. Shenzhen ETEC isolates showed a remarkable high diversity, which belonged to four E. coli phylogroups (A: 71%, B1: 13%, E: 10%, and D: 6%) and 15 ETEC lineages, with L11 (25%, O159:H34/O159:H43, ST218/ST3153), novel L2/4 (21%, O6:H16, ST48), and L4 (15%, O25:H16, ST1491) being major lineages. Heat-stable toxin (ST) was most prevalent (76%, STh: 60% STp: 16%), followed by heat-labile toxin (LT, 17%) and ST + LT (7%). One or multiple colonization factors (CFs) were identified in 68 (64%) isolates, with the common CFs being CS21 (48%) and CS6 (34%). Antimicrobial resistance mutation/gene profiles of genomes were concordant with the phenotype testing results of 52 representative isolates, which revealed high resistance rate to nalidixic acid (71%), ampicillin (69%), and ampicillin/sulbactam (46%), and demonstrated that the novel L2/4 was a multidrug-resistant lineage. This study provides novel insight into the genomic epidemiology and antimicrobial susceptibility profile of ETEC infections in indigenous adults for the first time, which further improves our understanding on ETEC epidemiology and has implications for the development of vaccine and future surveillance and prevention of ETEC infections.

Isolation, Characterization, and Genomic Investigation of a Phytopathogenic Strain of Stenotrophomonas maltophilia

Stenotrophomonas maltophilia is ubiquitous in diverse environmental habitats. It merits significant concern because of its increasing incidence of nosocomial and community-acquired infection in immunocompromised patients and multiple drug resistance. It is rarely reported as a phytopathogen except in causing white stripe disease of rice in India and postharvest fruit rot of Lanzhou lily. For this study, Dickeya zeae and S. maltophilia strains were simultaneously isolated from soft rot leaves of Clivia miniata in Guangzhou, China, and were both demonstrated to be pathogenic to the host. Compared with the D. zeae strains, S. maltophilia strains propagated faster for greater growth in lysogeny broth medium and produced no cellulases or polygalacturonases, but did produce more proteases and fewer extracellular polysaccharides. Furthermore, S. maltophilia strains swam and swarmed dramatically less on semisolid media, but formed a great many more biofilms. Both D. zeae and S. maltophilia strains isolated from clivia caused rot symptoms on other monocot hosts, but not on dicots. Similar to previously reported S. maltophilia strains isolated from other sources, the strain JZL8 survived under many antibiotic stresses. The complete genome sequence of S. maltophilia strain JZL8 consists of a chromosome of 4,635,432 bp without a plasmid. Pan-genome analysis of JZL8 and 180 other S. maltophilia strains identified 50 genes that are unique to JZL8, seven of which implicate JZL8 as the potential pathogen contributor in plants. JZL8 also contains three copies of Type I Secretion System machinery; this is likely responsible for its greater production of proteases. Findings from this study extend our knowledge on the host range of S. maltophilia and provide insight into the phenotypic and genetic features underlying the plant pathogenicity of JZL8.

Mutation rate dynamics reflect ecological change in an emerging zoonotic pathogen

Mutation rates vary both within and between bacterial species, and understanding what drives this variation is essential for understanding the evolutionary dynamics of bacterial populations. In this study, we investigate two factors that are predicted to influence the mutation rate: ecology and genome size. We conducted mutation accumulation experiments on eight strains of the emerging zoonotic pathogen Streptococcus suis. Natural variation within this species allows us to compare tonsil carriage and invasive disease isolates, from both more and less pathogenic populations, with a wide range of genome sizes. We find that invasive disease isolates have repeatedly evolved mutation rates that are higher than those of closely related carriage isolates, regardless of variation in genome size. Independent of this variation in overall rate, we also observe a stronger bias towards G/C to A/T mutations in isolates from more pathogenic populations, whose genomes tend to be smaller and more AT-rich. Our results suggest that ecology is a stronger correlate of mutation rate than genome size over these timescales, and that transitions to invasive disease are consistently accompanied by rapid increases in mutation rate. These results shed light on the impact that ecology can have on the adaptive potential of bacterial pathogens.

Drivers of methicillin-resistant Staphylococcus aureus (MRSA) lineage replacement in China

BACKGROUND

Methicillin-resistant Staphylococcus aureus (MRSA) is a major nosocomial pathogen subdivided into lineages termed sequence types (STs). Since the 1950s, successive waves of STs have appeared and replaced previously dominant lineages. One such event has been occurring in China since 2013, with community-associated (CA-MRSA) strains including ST59 largely replacing the previously dominant healthcare-associated (HA-MRSA) ST239. We previously showed that ST59 isolates tend to have a competitive advantage in growth experiments against ST239. However, the underlying genomic and phenotypic drivers of this replacement event are unclear.

METHODS

Here, we investigated the replacement of ST239 using whole-genome sequencing data from 204 ST239 and ST59 isolates collected in Chinese hospitals between 1994 and 2016. We reconstructed the evolutionary history of each ST and considered two non-mutually exclusive hypotheses for ST59 replacing ST239: antimicrobial resistance (AMR) profile and/or ability to colonise and persist in the environment through biofilm formation. We also investigated the differences in cytolytic activity, linked to higher virulence, between STs. We performed an association study using the presence and absence of accessory virulence genes.

RESULTS

ST59 isolates carried fewer AMR genes than ST239 and showed no evidence of evolving towards higher AMR. Biofilm production was marginally higher in ST59 overall, though this effect was not consistent across sub-lineages so is unlikely to be a sole driver of replacement. Consistent with previous observations of higher virulence in CA-MRSA STs, we observed that ST59 isolates exhibit significantly higher cytolytic activity than ST239 isolates, despite carrying on average fewer putative virulence genes. Our association study identified the chemotaxis inhibitory protein (chp) as a strong candidate for involvement in the increased virulence potential of ST59. We experimentally validated the role of chp in increasing the virulence potential of ST59 by creating Δchp knockout mutants, confirming that ST59 can carry chp without a measurable impact on fitness.

CONCLUSIONS

Our results suggest that the ongoing replacement of ST239 by ST59 in China is not associated to higher AMR carriage or biofilm production. However, the increase in ST59 prevalence is concerning since it is linked to a higher potential for virulence, aided by the carriage of the chp gene.

Phylogeographic Clustering Suggests that Distinct Clades of Salmonella enterica Serovar Mississippi Are Endemic in Australia, the United Kingdom, and the United States

Salmonella enterica serovar Mississippi is the 2nd and 14th leading cause of human clinical salmonellosis in the Australian island state of Tasmania and the United States, respectively. Despite its public health relevance, relatively little is known about this serovar. Comparison of whole-genome sequence (WGS) data of S. Mississippi isolates with WGS data for 317 additional S. enterica serovars placed one clade of S. Mississippi within S. enterica clade B (“clade B Mississippi”) and the other within section Typhi in S. enterica clade A (“clade A Mississippi”), suggesting that these clades evolved from different ancestors. Phylogenetic analysis of 364 S. Mississippi isolates from Australia, the United Kingdom, and the United States suggested that the isolates cluster geographically, with U.S. and Australian isolates representing different subclades (Ai and Aii, respectively) within clade A Mississippi and clade B isolates representing the predominant S. Mississippi isolates in the United Kingdom. Intraclade comparisons suggested that different mobile elements, some of which encode virulence factors, are responsible for the observed differences in gene content among isolates within these clades. Specifically, genetic differences among clade A isolates reflect differences in prophage contents, while differences among clade B isolates are due to the acquisition of a 47.1-kb integrative conjugative element (ICE). Phylogenies inferred from antigenic components (fliC, fljB, and O-antigen-processing genes) support that clade A and B Mississippi isolates acquired these loci from different ancestral serovars. Overall, these data support that different S. Mississippi phylogenetic clades are endemic in Australia, the United Kingdom, and the United States.

IMPORTANCE The number of known so-called “polyphyletic” serovars (i.e., phylogenetically distinct clades with the same O and H antigenic formulas) continues to increase as additional Salmonella isolates are sequenced. While serotyping remains a valuable tool for reporting and monitoring Salmonella, more discriminatory analyses for classifying polyphyletic serovars may improve surveillance efforts for these serovars, as we found that for S. Mississippi, distinct genotypes predominate at different geographic locations. Our results suggest that the acquisition of genes encoding O and H antigens from different ancestors led to the emergence of two Mississippi clades. Furthermore, our results suggest that different mobile elements contribute to the microevolution and diversification of isolates within these two clades, which has implications for the acquisition of novel adaptations, such as virulence factors.

Pan-Genome of the Genus Streptomyces and Prioritization of Biosynthetic Gene Clusters With Potential to Produce Antibiotic Compounds

Species of the genus Streptomyces are known for their ability to produce multiple secondary metabolites; their genomes have been extensively explored to discover new bioactive compounds. The richness of genomic data currently available allows filtering for high quality genomes, which in turn permits reliable comparative genomics studies and an improved prediction of biosynthetic gene clusters (BGCs) through genome mining approaches. In this work, we used 121 genome sequences of the genus Streptomyces in a comparative genomics study with the aim of estimating the genomic diversity by protein domains content, sequence similarity of proteins and conservation of Intergenic Regions (IGRs). We also searched for BGCs but prioritizing those with potential antibiotic activity. Our analysis revealed that the pan-genome of the genus Streptomyces is clearly open, with a high quantity of unique gene families across the different species and that the IGRs are rarely conserved. We also described the phylogenetic relationships of the analyzed genomes using multiple markers, obtaining a trustworthy tree whose relationships were further validated by Average Nucleotide Identity (ANI) calculations. Finally, 33 biosynthetic gene clusters were detected to have potential antibiotic activity and a predicted mode of action, which might serve up as a guide to formulation of related experimental studies.