EnteroBase: hierarchical clustering of 100 000s of bacterial genomes into species/subspecies and populations

Modulation of Salmonella virulence by a novel SPI-2 injectisome effector that interacts with the dystrophin-associated protein complex

The injectisome encoded by Salmonella pathogenicity island 2 (SPI-2) had been thought to translocate 28 effectors. Here, we used a proteomic approach to characterize the secretome of a clinical strain of invasive non-typhoidal Salmonella enterica serovar Enteritidis that had been mutated to cause hyper-secretion of the SPI-2 injectisome effectors. Along with many known effectors, we discovered the novel SseM protein. sseM is widely distributed among the five subspecies of Salmonella enterica, is found in many clinically relevant serovars, and is co-transcribed with pipB2, a SPI-2 effector gene. The translocation of SseM required a functional SPI-2 injectisome. Following expression in human cells, SseM interacted with five components of the dystrophin-associated protein complex (DAPC), namely, β-2-syntrophin, utrophin/dystrophin, α-catulin, α-dystrobrevin, and β-dystrobrevin. The interaction between SseM and β-2-syntrophin and α-dystrobrevin was verified in Salmonella Typhimurium-infected cells and relied on the postsynaptic density-95/discs large/zonula occludens-1 (PDZ) domain of β-2-syntrophin and a sequence corresponding to a PDZ-binding motif (PBM) in SseM. A ΔsseM mutant strain had a small competitive advantage over the wild-type strain in the S. Typhimurium/mouse model of systemic disease. This phenotype was complemented by a plasmid expressing wild-type SseM from S. Typhimurium or S. Enteritidis and was dependent on the PBM of SseM. Therefore, a PBM within a Salmonella effector mediates interactions with the DAPC and modulates the systemic growth of bacteria in mice. Furthermore, the ΔsseM mutant strain displayed enhanced replication in bone marrow-derived macrophages, demonstrating that SseM restrains intracellular bacterial growth to modulate Salmonella virulence.

IMPORTANCE

In Salmonella enterica, the injectisome machinery encoded by Salmonella pathogenicity island 2 (SPI-2) is conserved among the five subspecies and delivers proteins (effectors) into host cells, which are required for Salmonella virulence. The identification and functional characterization of SPI-2 injectisome effectors advance our understanding of the interplay between Salmonella and its host(s). Using an optimized method for preparing secreted proteins and a clinical isolate of the invasive non-typhoidal Salmonella enterica serovar Enteritidis strain D24359, we identified 22 known SPI-2 injectisome effectors and one new effector-SseM. SseM modulates bacterial growth during murine infection and has a sequence corresponding to a postsynaptic density-95/discs large/zonula occludens-1 (PDZ)-binding motif that is essential for interaction with the PDZ-containing host protein β-2-syntrophin and other components of the dystrophin-associated protein complex (DAPC). To our knowledge, SseM is unique among Salmonella effectors in containing a functional PDZ-binding motif and is the first bacterial protein to target the DAPC.

Antimicrobial resistance and phylogenetic relatedness of extended-spectrum ß-lactamase (ESBL)-producing Escherichia coli in peridomestic rats (Rattus norvegicus and Rattus tanezumi) linked to city areas and animal farms in Hong Kong

Extended-spectrum β-lactamase-producing Escherichia (E.) coli (ESBL-EC) in the clinical setting have emerged as a major threat to public and animal health. Wildlife, including Rattus spp. may serve as reservoirs and spreaders of ESBL-EC in the environment. Peridomestic rats are well adapted to living in proximity to humans and animals in a variety of urban and agricultural environments and may serve as sentinels to identify variations of ESBL-EC within their different habitats. In this study, a set of 221 rats (Rattus norvegicus, R. tanezumi, R. andamanensis, and Niviventer huang) consisting of 104 rats from city areas, 44 from chicken farms, 52 from pig farms, and 21 from stables of horse-riding schools were screened for ESBL-EC. Overall, a total of 134 ESBL-EC were isolated from the caecal samples of 130 (59%) rats. The predominant blaESBL genes were blaCTX-M-14, blaCTX-M-15, blaCTX-M-55, and blaCTX-M-65. Phylogenetic analysis revealed a total of 62 sequence types (STs) and 17 SNP clusters. E. coli ST10 and ST155 were common to ESBL-EC from city areas and chicken farms, and ST44 were found among ESBL-EC from city areas and pig farms. Extra-intestinal pathogenic E. coli (ExPEC) ST69, ST131 and ST1193 were found exclusively among rats from city areas, and avian pathogenic E. coli (APEC) ST177 was restricted to ESBL-EC originating from chicken farms. Phylogenetic analysis showed that the populations of rodent ESBL-EC from city areas, chicken farms and pig farms were genetically different, suggesting a certain degree of partitioning between the human and animal locations. This study contributes to current understanding of ESBL-EC occurring in rats in ecologically diverse locations.

Homology of Escherichia coli isolated from urine and vagina and their antimicrobial susceptibility in postmenopausal women with recurrent cystitis

Uropathogenic Escherichia coli (UPEC) is a typical cystitis-causing organism that can migrate from the vagina to the bladder and cause recurrent cystitis (RC). Few reports have compared the characteristics of urinary and vaginal UPEC in patients with RC. We carried out molecular biological analyses of Escherichia coli (E. coli) strains and their antimicrobial susceptibility to assess the association between urinary and vaginally UPEC. We included E. coli isolated from urinary and vaginal samples at the onset of cystitis in postmenopausal women with RC between 2014 and 2019 in our hospital. Pulsed-field gel electrophoresis (PFGE) was performed using a restriction enzyme (Xba I). These sequences were compared with 17 antimicrobial susceptibilities determined by a micro-liquid dilution method. Multilocus sequence typing (MLST) and classification of extended-spectrum β-lactamase (ESBL) genotypes by multiplex polymerase chain reaction (PCR) were performed on ESBL-producing E. coli. We analyzed 14 specimens (each seven urine and vaginal) from seven patients in total. On PFGE, the similarity of urinary and vaginal E. coli per patient ranged from 89.5 to 100 %, including four patients with 100 % matches. MLST demonstrated that 29 % (4/14 specimens) were strain sequence type 131. Two specimens contained ESBL-producing strains and identified the CTX-M-27 genotype for each specimen. For each patient, antimicrobial susceptibilities between urinary and vaginal E. coli were mostly identical. Thus, urinary- and vaginally-derived E. coli were identical in postmenopausal women with RC. Management targeting both urinary and vaginal UPEC is essential for RC, indicating the importance of a vagina-targeted approach.

Centralized industrialization of pork in Europe and America contributes to the global spread of Salmonella enterica

Salmonella enterica causes severe food-borne infections through contamination of the food supply chain. Its evolution has been associated with human activities, especially animal husbandry. Advances in intensive farming and global transportation have substantially reshaped the pig industry, but their impact on the evolution of associated zoonotic pathogens such as S. enterica remains unresolved. Here we investigated the population fluctuation, accumulation of antimicrobial resistance genes and international serovar Choleraesuis transmission of nine pig-enriched S. enterica populations comprising more than 9,000 genomes. Most changes were found to be attributable to the developments of the modern pig industry. All pig-enriched salmonellae experienced host transfers in pigs and/or population expansions over the past century, with pigs and pork having become the main sources of S. enterica transmissions to other hosts. Overall, our analysis revealed strong associations between the transmission of pig-enriched salmonellae and the global pork trade.

SaLTy: a novel Staphylococcus aureus Lineage Typer

Staphylococcus aureus asymptomatically colonises 30 % of humans but can also cause a range of diseases, which can be fatal. In 2017 S. aureus was associated with 20 000 deaths in the USA alone. Dividing S. aureus isolates into smaller sub-groups can reveal the emergence of distinct sub-populations with varying potential to cause infections. Despite multiple molecular typing methods categorising such sub-groups, they do not take full advantage of S. aureus genome sequences when describing the fundamental population structure of the species. In this study, we developed Staphylococcus aureus Lineage Typing (SaLTy), which rapidly divides the species into 61 phylogenetically congruent lineages. Alleles of three core genes were identified that uniquely define the 61 lineages and were used for SaLTy typing. SaLTy was validated on 5000 genomes and 99.12 % (4956/5000) of isolates were assigned the correct lineage. We compared SaLTy lineages to previously calculated clonal complexes (CCs) from BIGSdb (n=21 173). SALTy improves on CCs by grouping isolates congruently with phylogenetic structure. SaLTy lineages were further used to describe the carriage of Staphylococcal chromosomal cassette containing mecA (SCCmec) which is carried by methicillin-resistant S. aureus (MRSA). Most lineages had isolates lacking SCCmec and the four largest lineages varied in SCCmec over time. Classifying isolates into SaLTy lineages, which were further SCCmec typed, allowed SaLTy to describe high-level MRSA epidemiology. We provide SaLTy as a simple typing method that defines phylogenetic lineages (https://github.com/LanLab/SaLTy). SaLTy is highly accurate and can quickly analyse large amounts of S. aureus genome data. SaLTy will aid the characterisation of S. aureus populations and ongoing surveillance of sub-groups that threaten human health.

Survival of highly related ESBL- and pAmpC- producing Escherichia coli in broiler farms identified before and after cleaning and disinfection using cgMLST

BACKGROUND

Broiler chickens are frequently colonized with Extended-Spectrum Beta-Lactamase- (ESBL-) and plasmid mediated AmpC Beta-Lactamase- (pAmpC-) producing Enterobacterales, and we are confronted with the potential spread of these resistant bacteria in the food chain, in the environment, and to humans. Research focused on identifying of transmission routes and investigating potential intervention measures against ESBL- and pAmpC- producing bacteria in the broiler production chain. However, few data are available on the effects of cleaning and disinfection (C&D) procedures in broiler stables on ESBL- and pAmpC- producing bacteria.

RESULTS

We systematically investigated five broiler stables before and after C&D and identified potential ESBL- and pAmpC- colonization sites after C&D in the broiler stables, including the anteroom and the nearby surrounding environment of the broiler stables. Phenotypically resistant E. coli isolates grown on MacConkey agar with cefotaxime were further analyzed for their beta-lactam resistance genes and phylogenetic groups, as well as the relation of isolates from the investigated stables before and after C&D by whole genome sequencing. Survival of ESBL- and pAmpC- producing E. coli is highly likely at sites where C&D was not performed or where insufficient cleaning was performed prior to disinfection. For the first time, we showed highly related ESBL-/pAmpC- producing E. coli isolates detected before and after C&D in four of five broiler stables examined with cgMLST. Survival of resistant isolates in investigated broiler stables as well as transmission of resistant isolates from broiler stables to the anteroom and surrounding environment and between broiler farms was shown. In addition, enterococci (frequently utilized to detect fecal contamination and for C&D control) can be used as an indicator bacterium for the detection of ESBL-/pAmpC- E. coli after C&D.

CONCLUSION

We conclude that C&D can reduce ESBL-/pAmpC- producing E. coli in conventional broiler stables, but complete ESBL- and pAmpC- elimination does not seem to be possible in practice as several factors influence the C&D outcome (e.g. broiler stable condition, ESBL-/pAmpC- status prior to C&D, C&D procedures used, and biosecurity measures on the farm). A multifactorial approach, combining various hygiene- and management measures, is needed to reduce ESBL-/pAmpC- E. coli in broiler farms.

A new Salmonella enterica serovar that was isolated from a wild sparrow presents a distinct genetic, metabolic and virulence profile

Salmonella enterica is a ubiquitous and clinically-important bacterial pathogen, able to infect and cause different diseases in a wide range of hosts. Here, we report the isolation and characterization of a new S. enterica serovar (13,23:i:-; S. Tirat-Zvi), belonging to the Havana supper-lineage that was isolated from a wild house sparrow (Passer domesticus) in Israel. Whole genome sequencing and complete assembly of its genome indicated a plasmid-free, 4.7 Mb genome that carries the Salmonella pathogenicity islands 1-6, 9, 19 and an integrative and conjugative element (ICE), encoding arsenic resistance genes. Phenotypically, S. Tirat-Zvi isolate TZ282 was motile, readily formed biofilm, more versatile in carbon source utilization than S. Typhimurium and highly tolerant to arsenic, but impaired in host cell invasion. In-vivo infection studies indicated that while S. Tirat-Zvi was able to infect and cause an acute inflammatory enterocolitis in young chicks, it was compromised in mice colonization and did not cause an inflammatory colitis in mice compared to S. Typhimurium. We suggest that these phenotypes reflect the distinctive ecological niche of this new serovar and its evolutionary adaptation to passerine birds, as a permissive host. Moreover, these results further illuminate the genetic, phenotypic and ecological diversity of S. enterica pathovars.

An imputed ancestral reference genome for the Mycobacterium tuberculosis complex better captures structural genomic diversity for reference-based alignment workflows

Reference-based alignment of short-reads is a widely used technique in genomic analysis of the Mycobacterium tuberculosis complex (MTBC) and the choice of reference sequence impacts the interpretation of analyses. The most widely used reference genomes include the ATCC type strain (H37Rv) and the putative MTBC ancestral sequence of Comas et al. both of which are based on a lineage 4 sequence. As such, these reference sequences do not capture all of the structural variation known to be present in the ancestor of the MTBC. To better represent the base of the MTBC, we generated an imputed ancestral genomic sequence, termed MTBC0 from reference-free alignments of closed MTBC genomes. When used as a reference sequence in alignment workflows, MTBC0 mapped more short sequencing reads and called more pairwise SNPs relative to the Comas et al. sequence while exhibiting minimal impact on the overall phylogeny of MTBC. The results also show that MTBC0 provides greater fidelity in capturing genomic variation and allows for the inclusion of regions absent from H37Rv in standard MTBC workflows without additional steps. The use of MTBC0 as an ancestral reference sequence in standard workflows modestly improved read mapping, SNP calling and intuitively facilitates the study of structural variation and evolution in MTBC.

Epidemiological trends and antimicrobial resistance in Salmonella enterica serovar Typhimurium clones in Taiwan between 2004 and 2019

OBJECTIVES

We investigated the temporal trends of Salmonella enterica serovar Typhimurium (S. Typhimurium) clones in Taiwan from 2004 to 2019, focusing on antimicrobial resistance (AMR), resistance genetic determinants, and plasmid types.

METHODS

Salmonella isolates were characterized using pulsed-field gel electrophoresis (PFGE), whole-genome sequencing, and antimicrobial susceptibility testing. Clones were defined using PFGE clustering and the hierarchical cgMLST clustering (HierCC) assignments.

RESULTS

Seven major S. Typhimurium clones, HC100_2, 13, 41, 305, 310, 501, and 46261, accounted for 97.6% (8079/8275) of human isolates in Taiwan. Each clone displayed a unique AMR profile, resistance genetic determinants, and plasmid types. Four highly resistant clones (HC100_2, 41, 305, and 310) exhibited multiple resistance in 86.5% to 96.1% of isolates. HC100_305 and HC100_2 were pandemic multidrug-resistant clones, characterized by resistance to ampicillin, chloramphenicol, streptomycin, sulfonamides, and tetracycline (ACSSuT) and ASSuT, respectively. The prevalence of the ACSSuT clone decreased from 68.7% of S. Typhimurium isolates in 2004 to 1.7% in 2019, while the ASSuT clone emerged in 2007 and became the largest clone after 2010. Several plasmids, including IncHI2-IncHI2A, IncC, IncFIB(K), and IncI1-1(α), carried multiple resistance genes or were associated with the carriage of mph(A), blaCMY-2, and blaDHA-1.

CONCLUSIONS

Between 2004 and 2019, Taiwan experienced the emergence, prevalence, and subsequent decline of several highly resistant S. Typhimurium clones. The clones defined using the HierCC approach have global comparability. The increasing resistance to third-generation cephalosporins, cephamycins, ciprofloxacin, and azithromycin in recent years poses a significant medical concern.

Genomic analysis of Shigella isolates from Lebanon reveals marked genetic diversity and antimicrobial resistance

In this study, we characterized 54 clinical isolates of Shigella collected in North Lebanon between 2009 and 2017 through phenotypic and genomic analyses. The most prevalent serogroup was S. sonnei, accounting for 46.3 % (25/54) of the isolates, followed by S. flexneri (27.8 %, 15/54), S. boydii (18.5 %, 10/54) and S. dysenteriae (7.4 %, 4/54). Only three isolates were pan-susceptible, and 87 % (47/54) of the isolates had multidrug resistance phenotypes. Notably, 27.8 % (15/54) of the isolates were resistant to third-generation cephalosporins (3GCs) and 77.8 % (42/54) were resistant to nalidixic acid. 3GC resistance was mediated by the extended-spectrum beta-lactamase genes bla CTX-M-15 and bla CTX-M-3, which were present on various plasmids. Quinolone resistance was conferred by single point mutations in the gyrA DNA gyrase gene, leading to GyrA S83L, GyrA D87Y or GyrA S83A amino acid substitutions. This is the first study, to our knowledge, to provide genomic insights into the serotypes of Shigella circulating in Lebanon and the various antimicrobial resistance determinants carried by these strains.

Distributed genotyping and clustering of Neisseria strains reveal continual emergence of epidemic meningococcus over a century

Core genome multilocus sequence typing (cgMLST) is commonly used to classify bacterial strains into different types, for taxonomical and epidemiological applications. However, cgMLST schemes require central databases for the nomenclature of new alleles and sequence types, which must be synchronized worldwide and involve increasingly intensive calculation and storage demands. Here, we describe a distributed cgMLST (dcgMLST) scheme that does not require a central database of allelic sequences and apply it to study evolutionary patterns of epidemic and endemic strains of the genus Neisseria. We classify 69,994 worldwide Neisseria strains into multi-level clusters that assign species, lineages, and local disease outbreaks. We divide Neisseria meningitidis into 168 endemic lineages and three epidemic lineages responsible for at least 9 epidemics in the past century. According to our analyses, the epidemic and endemic lineages experienced very different population dynamics in the past 100 years. Epidemic lineages repetitively emerged from endemic lineages, disseminated worldwide, and apparently disappeared rapidly afterward. We propose a stepwise model for the evolutionary trajectory of epidemic lineages in Neisseria, and expect that the development of similar dcgMLST schemes will facilitate epidemiological studies of other bacterial pathogens.

Exchange of Carbapenem-Resistant Escherichia coli Sequence Type 38 Intercontinentally and among Wild Bird, Human, and Environmental Niches

Carbapenem-resistant Enterobacteriaceae (CRE) are a global threat to human health and are increasingly being isolated from nonclinical settings. OXA-48-producing Escherichia coli sequence type 38 (ST38) is the most frequently reported CRE type in wild birds and has been detected in gulls or storks in North America, Europe, Asia, and Africa. The epidemiology and evolution of CRE in wildlife and human niches, however, remains unclear. We compared wild bird origin E. coli ST38 genome sequences generated by our research group and publicly available genomic data derived from other hosts and environments to (i) understand the frequency of intercontinental dispersal of E. coli ST38 clones isolated from wild birds, (ii) more thoroughly measure the genomic relatedness of carbapenem-resistant isolates from gulls sampled in Turkey and Alaska, USA, using long-read whole-genome sequencing and assess the spatial dissemination of this clone among different hosts, and (iii) determine whether ST38 isolates from humans, environmental water, and wild birds have different core or accessory genomes (e.g., antimicrobial resistance genes, virulence genes, plasmids) which might elucidate bacterial or gene exchange among niches. Our results suggest that E. coli ST38 strains, including those resistant to carbapenems, are exchanged between humans and wild birds, rather than separately maintained populations within each niche. Furthermore, despite close genetic similarity among OXA-48-producing E. coli ST38 clones from gulls in Alaska and Turkey, intercontinental dispersal of ST38 clones among wild birds is uncommon. Interventions to mitigate the dissemination of antimicrobial resistance throughout the environment (e.g., as exemplified by the acquisition of carbapenem resistance by birds) may be warranted. IMPORTANCE Carbapenem-resistant bacteria are a threat to public health globally and have been found in the environment as well as the clinic. Some bacterial clones are associated with carbapenem resistance genes, such as Escherichia coli sequence type 38 (ST38) and the carbapenemase gene blaOXA-48. This is the most frequently reported carbapenem-resistant clone in wild birds, though it was unclear if it circulated within wild bird populations or was exchanged among other niches. The results from this study suggest that E. coli ST38 strains, including those resistant to carbapenems, are frequently exchanged among wild birds, humans, and the environment. Carbapenem-resistant E. coli ST38 clones in wild birds are likely acquired from the local environment and do not constitute an independent dissemination pathway within wild bird populations. Management actions aimed at preventing the environmental dissemination and acquisition of antimicrobial resistance by wild birds may be warranted.

Bacterial genome sequencing tracks the housefly-associated dispersal of fluoroquinolone- and cephalosporin-resistant Escherichia coli from a pig farm

The regular use of antimicrobials in livestock production selects for antimicrobial resistance. The potential impact of this practice on human health needs to be studied in more detail, including the role of the environment for the persistence and transmission of antimicrobial-resistant bacteria. During an investigation of a pig farm and its surroundings in Brandenburg, Germany, we detected abundant cephalosporin- and fluoroquinolone-resistant Escherichia coli in pig faeces, sedimented dust, and house flies (Musca domestica). Genome sequencing of E. coli isolates revealed large phylogenetic diversity and plasmid-borne extended-spectrum beta lactamase (ESBL) genes CTX-M-1 in multiple strains. [Correction added on 28 February 2023, after first online publication: In the preceding sentence, 'and TEM-1' was previously included but has been deleted in this version.] Close genomic relationships indicated frequent transmission of antimicrobial-resistant E. coli between pigs from different herds and across buildings of the farm and suggested dust and flies as vectors for dissemination of faecal pathogens. Strikingly, we repeatedly recovered E. coli from flies collected up to 2 km away from the source, whose genome sequences were identical or closely related to those from pig faeces isolates, indicating the fly-associated transport of diverse ESBL-producing E. coli from the pig farm into urban habitation areas. The observed proximity of contaminated flies to human households poses a risk of transmission of antimicrobial-resistant enteric pathogens from livestock to man.

Recent progress on methods for estimating and updating large phylogenies

With the increased availability of sequence data and even of fully sequenced and assembled genomes, phylogeny estimation of very large trees (even of hundreds of thousands of sequences) is now a goal for some biologists. Yet, the construction of these phylogenies is a complex pipeline presenting analytical and computational challenges, especially when the number of sequences is very large. In the past few years, new methods have been developed that aim to enable highly accurate phylogeny estimations on these large datasets, including divide-and-conquer techniques for multiple sequence alignment and/or tree estimation, methods that can estimate species trees from multi-locus datasets while addressing heterogeneity due to biological processes (e.g. incomplete lineage sorting and gene duplication and loss), and methods to add sequences into large gene trees or species trees. Here we present some of these recent advances and discuss opportunities for future improvements. This article is part of a discussion meeting issue 'Genomic population structures of microbial pathogens'.