PLSDB: a resource of complete bacterial plasmids

Using a combination of short- and long-read sequencing to investigate the diversity in plasmid- and chromosomally encoded extended-spectrum beta-lactamases (ESBLs) in clinical Shigella and Salmonella isolates in Belgium

For antimicrobial resistance (AMR) surveillance, it is important not only to detect AMR genes, but also to determine their plasmidic or chromosomal location, as this will impact their spread differently. Whole-genome sequencing (WGS) is increasingly used for AMR surveillance. However, determining the genetic context of AMR genes using only short-read sequencing is complicated. The combination with long-read sequencing offers a potential solution, as it allows hybrid assemblies. Nevertheless, its use in surveillance has so far been limited. This study aimed to demonstrate its added value for AMR surveillance based on a case study of extended-spectrum beta-lactamases (ESBLs). ESBL genes have been reported to occur also on plasmids. To gain insight into the diversity and genetic context of ESBL genes detected in clinical isolates received by the Belgian National Reference Center between 2013 and 2018, 100 ESBL-producing Shigella and 31 ESBL-producing Salmonella were sequenced with MiSeq and a representative selection of 20 Shigella and six Salmonella isolates additionally with MinION technology, allowing hybrid assembly. The bla CTX-M-15 gene was found to be responsible for a rapid rise in the ESBL Shigella phenotype from 2017. This gene was mostly detected on multi-resistance-carrying IncFII plasmids. Based on clustering, these plasmids were determined to be distinct from the circulating plasmids before 2017. They were spread to different Shigella species and within Shigella sonnei between multiple genotypes. Another similar IncFII plasmid was detected after 2017 containing bla CTX-M-27 for which only clonal expansion occurred. Matches of up to 99 % to plasmids of various bacterial hosts from all over the world were found, but global alignments indicated that direct or recent ESBL-plasmid transfers did not occur. It is most likely that travellers introduced these in Belgium and subsequently spread them domestically. However, a clear link to a specific country could not be made. Moreover, integration of bla CTX-M in the chromosome of two Shigella isolates was determined for the first time, and shown to be related to ISEcp1. In contrast, in Salmonella, ESBL genes were only found on plasmids, of which bla CTX-M-55 and IncHI2 were the most prevalent, respectively. No matching ESBL plasmids or cassettes were detected between clinical Shigella and Salmonella isolates. The hybrid assembly data allowed us to check the accuracy of plasmid prediction tools. MOB-suite showed the highest accuracy. However, these tools cannot replace the accuracy of long-read and hybrid assemblies. This study illustrates the added value of hybrid assemblies for AMR surveillance and shows that a strategy where even just representative isolates of a collection used for hybrid assemblies could improve international AMR surveillance as it allows plasmid tracking.

Characterization of antibiotic resistomes by reprogrammed bacteriophage-enabled functional metagenomics in clinical strains

Functional metagenomics is a powerful experimental tool to identify antibiotic resistance genes (ARGs) in the environment, but the range of suitable host bacterial species is limited. This limitation affects both the scope of the identified ARGs and the interpretation of their clinical relevance. Here we present a functional metagenomics pipeline called Reprogrammed Bacteriophage Particle Assisted Multi-species Functional Metagenomics (DEEPMINE). This approach combines and improves the use of T7 bacteriophage with exchanged tail fibres and targeted mutagenesis to expand phage host-specificity and efficiency for functional metagenomics. These modified phage particles were used to introduce large metagenomic plasmid libraries into clinically relevant bacterial pathogens. By screening for ARGs in soil and gut microbiomes and clinical genomes against 13 antibiotics, we demonstrate that this approach substantially expands the list of identified ARGs. Many ARGs have species-specific effects on resistance; they provide a high level of resistance in one bacterial species but yield very limited resistance in a related species. Finally, we identified mobile ARGs against antibiotics that are currently under clinical development or have recently been approved. Overall, DEEPMINE expands the functional metagenomics toolbox for studying microbial communities.

Characterization and whole-genome sequencing of an extreme arsenic-tolerant Citrobacter freundii SRS1 strain isolated from Savar area in Bangladesh

Citrobacter freundii SRS1, gram-negative bacteria, were isolated from Savar, Bangladesh. The strain could tolerate up to 80 mmol L-1 sodium arsenite, 400 mmol L-1 sodium arsenate, 5 mmol L-1 manganese sulfate, 3 mmol L-1 lead nitrate, 2.5 mmol L-1 cobalt chloride, 2.5  mmol L-1 cadmium acetate, and 2.5 mmol L-1 chromium chloride. The whole-genome sequencing revealed that the genome size of C. freundii SRS1 is estimated to be 5.4 Mb long, and the G + C content is 51.7%. The genome of C. freundii SRS1 contains arsA, arsB, arsC, arsD, arsH, arsR, and acr3 genes for arsenic resistance; czcA, czcD, cbiN, and cbiM genes for cobalt resistance; chrA and chrB genes for chromium resistance; mntH, sitA, sitB, sitC, and sitD genes for manganese resistance; and zntA gene for lead and cadmium resistance. This novel acr3 gene has never previously been reported in any C. freundii strain except SRS1. A set of 130 completely sequenced strains of C. freundii was selected for phylogenomic analysis. The phylogenetic tree showed that the SRS1 strain is closely related to the C. freundii 62 strain. Further analyses of the genes involved in metal and metalloid resistance might facilitate identifying the mechanisms and pathways involved in high metal resistance in the C. freundii SRS1 strain.

The Complete Genome Resource of Xanthomonas oryzae pv. oryzae CIX2779 Includes the First Sequence of a Plasmid for an African Representative of This Rice Pathogen

The bacterial plant pathogen Xanthomonas oryzae pv. oryzae is responsible for the foliar rice bacterial blight disease. Genetically contrasted, continent-specific, sublineages of this species can cause important damages to rice production both in Asia and Africa. We report on the genome of the CIX2779 strain of this pathogen, previously named NAI1 and originating from Niger. Oxford Nanopore long reads assembly and Illumina short reads polishing produced a genome sequence composed of a 4,725,792-bp circular chromosome and a 39,798-bp-long circular plasmid designated pCIX2779_1. The chromosome structure and base-level sequence are highly related to reference strains of African X. oryzae pv. oryzae and encode identical transcription activator-like effectors for virulence. Importantly, our in silico analysis strongly indicates that pCIX2779_1 is a genuine conjugative plasmid, the first indigenous one sequenced from an African strain of the X. oryzae species. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Sharing of Antimicrobial Resistance Genes between Humans and Food Animals

The prevalence and propagation of antimicrobial resistance (AMR) are serious global public health concerns. The large and the ever-increasing use of antibiotics in livestock is also considered a great concern. The extent of the similarity of acquired antibiotic resistance genes (ARGs) between humans and food animals and the driving factors underlying AMR transfer between them are not clear, although a link between ARGs in both hosts was proposed. To address this question, with swine and chicken as examples of food animals, we analyzed over 1,000 gut metagenomes of humans and food animals from over the world. A relatively high abundance and diversity of ARGs were observed in swine compared with those in humans as a whole. Commensal bacteria, particularly species from Clostridiales, contribute the most ARGs associated with mobile genetic elements (MGEs) and were found in both humans and food animals. Further studies demonstrate that overrepresented MGEs, namely, Tn4451/Tn4453 and TnAs3, are attributed mainly to the sharing between humans and food animals. A member of large resolvase family site-specific recombinases, TnpX, is found in Tn4451/Tn4453 which facilitates the insertions of the transient circular molecule. Although the variance in the transferability of ARGs in humans is higher than that in swine, a higher average transferability was observed in swine than that in humans. In conclusion, the potential antibiotic resistance hot spots with higher transferability in food animals observed in the present study emphasize the importance of surveillance for emerging resistance threats before they spread. IMPORTANCE Antimicrobial resistance (AMR) has proven to be a global public health concern. To conquer this increasingly worrying trend, an overarching, One Health approach has been used that brings together different sectors, but the fundamental knowledge of the relationship between humans, food animals, and their environments is not mature yet or is lacking in some aspect. With swine and chicken as examples of food animals, a large global data set of over 1,000 human and food animal gut metagenomes was analyzed with a focus on acquired antibiotic resistance genes (ARGs) associated with mobile genetic elements (MGEs) to answer this question. Outputs from this work open a new avenue to further our understanding of ARG transferability in food animals. It is a necessary milestone to better equip governmental agencies to monitor and pre-empt antibiotic resistance hot spots. This work will assist and give guidance on how to decipher other links within any One Health initiatives with expected positive feedback to human health.

Classification of bacterial plasmid and chromosome derived sequences using machine learning

Plasmids are important genetic elements that facilitate horizonal gene transfer between bacteria and contribute to the spread of virulence and antimicrobial resistance. Most bacterial genome sequences in the public archives exist in draft form with many contigs, making it difficult to determine if a contig is of chromosomal or plasmid origin. Using a training set of contigs comprising 10,584 chromosomes and 10,654 plasmids from the PATRIC database, we evaluated several machine learning models including random forest, logistic regression, XGBoost, and a neural network for their ability to classify chromosomal and plasmid sequences using nucleotide k-mers as features. Based on the methods tested, a neural network model that used nucleotide 6-mers as features that was trained on randomly selected chromosomal and plasmid subsequences 5kb in length achieved the best performance, outperforming existing out-of-the-box methods, with an average accuracy of 89.38% ± 2.16% over a 10-fold cross validation. The model accuracy can be improved to 92.08% by using a voting strategy when classifying holdout sequences. In both plasmids and chromosomes, subsequences encoding functions involved in horizontal gene transfer-including hypothetical proteins, transporters, phage, mobile elements, and CRISPR elements-were most likely to be misclassified by the model. This study provides a straightforward approach for identifying plasmid-encoding sequences in short read assemblies without the need for sequence alignment-based tools.

Identification of Genetic Markers for the Detection of Bacillus thuringiensis Strains of Interest for Food Safety

Bacillus thuringiensis (Bt), belonging to the Bacillus cereus (Bc) group, is commonly used as a biopesticide worldwide due to its ability to produce insecticidal crystals during sporulation. The use of Bt, especially subspecies aizawai and kurstaki, to control pests such as Lepidoptera, generally involves spraying mixtures containing spores and crystals on crops intended for human consumption. Recent studies have suggested that the consumption of commercial Bt strains may be responsible for foodborne outbreaks (FBOs). However, its genetic proximity to Bc strains has hindered the development of routine tests to discriminate Bt from other Bc, especially Bacillus cereus sensu stricto (Bc ss), well known for its involvement in FBOs. Here, to develop tools for the detection and the discrimination of Bt in food, we carried out a genome-wide association study (GWAS) on 286 complete genomes of Bc group strains to identify and validate in silico new molecular markers specific to different Bt subtypes. The analyses led to the determination and the in silico validation of 128 molecular markers specific to Bt, its subspecies aizawai, kurstaki and four previously described proximity clusters associated with these subspecies. We developed a command line tool based on a 14-marker workflow, to carry out a computational search for Bt-related markers from a putative Bc genome, thereby facilitating the detection of Bt of interest for food safety, especially in the context of FBOs.

Comparative Genomics of Lentilactobacillus parabuchneri isolated from dairy, KEM complex, Makgeolli, and Saliva Microbiomes

BACKGROUND

Lentilactobacillus parabuchneri is of particular concern in fermented food bioprocessing due to causing unwanted gas formation, cracks, and off-flavor in fermented dairy foods. This species is also a known culprit of histamine poisonings because of decarboxylating histidine to histamine in ripening cheese. Twenty-eight genomes in NCBI GenBank were evaluated via comparative analysis to determine genomic diversity within this species and identify potential avenues for reducing health associated risks and economic losses in the food industry caused by these organisms.

RESULT

Core genome-based phylogenetic analysis revealed four distinct major clades. Eight dairy isolates, two strains from an unknown source, and a saliva isolate formed the first clade. Three out of five strains clustered on clade 2 belonged to dairy, and the remaining two strains were isolated from the makgeolli and Korean effective microorganisms (KEM) complex. The third and fourth clade members were isolated from Tete de Moine and dairy-associated niches, respectively. Whole genome analysis on twenty-eight genomes showed ~ 40% of all CDS were conserved across entire strains proposing a considerable diversity among L. parabuchneri strains analyzed. After assigning CDS to their corresponding function, ~ 79% of all strains were predicted to carry putative intact prophages, and ~ 43% of the strains harbored at least one plasmid; however, all the strains were predicted to encode genomic island, insertion sequence, and CRISPR-Cas system. A type I-E CRISPR-Cas subgroup was identified in all the strains, with the exception of DSM15352, which carried a type II-A CRISPR-Cas system. Twenty strains were predicted to encode histidine decarboxylase gene cluster that belongs to not only dairy but also saliva, KEM complex, and unknown source. No bacteriocin-encoding gene(s) or antibiotic resistome was found in any of the L. parabuchneri strains screened.

CONCLUSION

The findings of the present work provide in-depth knowledge of the genomics of L. parabuchneri by comparing twenty-eight genomes available to date. For example, the hdc gene cluster was generally reported in cheese isolates; however, our findings in the current work indicated that it could also be encoded in those strains isolated from saliva, KEM complex, and unknown source. We think prophages are critical mobile elements of L. parabuchneri genomes that could pave the way for developing novel tools to reduce the occurrence of this unwanted species in the food industry.

IncP-type plasmids carrying genes for antibiotic resistance or for aromatic compound degradation are prevalent in sequenced Aromatoleum and Thauera strains

Self-transferable plasmids of the incompatibility group P-1 (IncP-1) are considered important carriers of genes for antibiotic resistance and other adaptive functions. In the laboratory, these plasmids have a broad host range; however, little is known about their in situ host profile. In this study, we discovered that Thauera aromatica K172^T , a facultative denitrifying microorganism capable of degrading various aromatic compounds, contains a plasmid highly similar to the IncP-1 ε archetype pKJK5. The plasmid harbours multiple antibiotic resistance genes and is maintained in strain K172^T for at least 1000 generations without selection pressure from antibiotics. In a subsequent search, we found additional nine IncP-type plasmids in a total of 40 sequenced genomes of the closely related genera Aromatoleum and Thauera. Six of these plasmids form a novel IncP-1 subgroup designated θ, four of which carry genes for anaerobic or aerobic degradation of aromatic compounds. Pentanucleotide sequence analyses (k-mer profiling) indicated that Aromatoleum spp. and Thauera spp. are among the most suitable hosts for the θ plasmids. Our results highlight the importance of IncP-1 plasmids for the genetic adaptation of these common facultative denitrifying bacteria and provide novel insights into the in situ host profile of these plasmids.

Environmental surveillance of ESBL and carbapenemase-producing gram-negative bacteria in a Ghanaian Tertiary Hospital

Background

The burden of antibiotic resistant infection is mainly felt in low-to-middle income countries, where the rate of antimicrobial resistance is largely under-surveyed and under huge pressure from unregulated, disparate and often self-guided access to antimicrobials. Nosocomial infections from hospital environments have been shown to be a particularly prevalent source of multi-drug resistant strains, yet surveillance of hospital environmental contamination is often not investigated.

Methods

The study was prospective, observational and cross-sectional, sampling 231 high and low touch surfaces from 15th March to 13th April 2021, from five wards in the Cape Coast Teaching Hospital, Ghana. Microbial growth in the presence of vancomycin and either meropenem or cefotaxime was examined and bacterial species were identified by MALDI-TOF. The presence of common extended-spectrum β-lactamases (ESBL) and carbapenemase antimicrobial resistance genes (ARG) were identified through PCR screening, which were confirmed by phenotypic antimicrobial susceptibility determination. Isolates positive for carbapenem resistance genes were sequenced using a multi-platform approach.

Results

We recovered microbial growth from 99% of swabs (n = 229/231) plated on agar in the absence of antimicrobials. Multiple sites were found to be colonised with resistant bacteria throughout the hospital setting. Bacteria with multi-drug resistance and ARG of concern were isolated from high and low touch points with evidence of strain dissemination throughout the environment. A total of 21 differing species of bacteria carrying ARG were isolated. The high prevalence of Acinetobacter baumannii carrying bla_NDM-1 observed was further characterised by whole genome sequencing and phylogenetic analysis to determine the relationship between resistant strains found in different wards.

Conclusion

Evidence of multiple clonal incursions of MDR bacteria of high sepsis risk were found in two separate wards for a regional hospital in Ghana. The prevalence of multiple bla_NDM carrying species in combination with combinations of ESBLs was particularly concerning and unexpected in Africa. We also identify strains carrying tet(X3), bla_VIM-5 or bla_DIM-1 showing a high diversity of carbapenamases present as a reservoir in a hospital setting. Findings of multi-drug resistant bacteria from multiple environmental sites throughout the hospital will inform future IPC practices and aid research prioritisation for AMR in Ghana.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13756-022-01090-2.

Finding of extended-spectrum beta-lactamase (ESBL)-producing Enterobacterales in wild game meat originating from several European countries: predominance of Moellerella wisconsensis producing CTX-M-1, November 2021

IntroductionMeat can be a vehicle for food-borne transmission of antimicrobial resistant bacteria and antimicrobial resistance genes. The occurrence of extended-spectrum beta-lactamase (ESBL) producing Enterobacterales has been observed in meat from livestock production but has not been well studied in meat from wild game.AimWe aimed to investigate, particularly in central Europe, to what extent ESBL-producing Enterobacterales may be present in wild game meat.MethodsA total of 111 samples of different types of game meat supplied by butchers, hunters, retail stores and a large game-processing establishment in Europe were screened for ESBL-producing Enterobacterales using a selective culture medium. Isolates were genotypically and phenotypically characterised.ResultsThirty-nine samples (35% of the total) yielded ESBL-producing Enterobacterales, with most (35/39) supplied by the game-processing establishment. Isolates included 32 Moellerella wisconsensis, 18 Escherichia coli and one Escherichia marmotae. PCR screening identified bla CTX-M-1 (n = 31), bla CTX-M-32 (n = 8), bla CTX-M-65 (n = 4), bla CTX-M-15 (n = 3), bla CTX-M-8 (n = 1), bla CTX-M-14 (n = 1), bla CTX-M-55 (n = 1), and bla SHV-12 (n = 2). Most E. coli belonged to phylogenetic group A (n = 7) or B1 (n = 9), but several isolates belonged to extraintestinal pathogenic E. coli (ExPEC) sequence types (ST)58 (n = 4), ST68 (n = 1) and ST540 (n = 1). Whole genome sequencing of six selected isolates localised bla CTX-M-1 on megaplasmids in four M. wisconsensis and bla CTX-M-32 on IncN_1 plasmids in one M. wisconsensis and one E. marmotae. Forty-eight isolates (94%) exhibited a multidrug-resistance phenotype.ConclusionWe found a high occurrence of ESBL-producing Enterobacterales in wild game meat, suggesting wildlife habitat pollution and possible microbial contamination events occurring during skinning or cutting carcasses.

Genomic analysis of sewage from 101 countries reveals global landscape of antimicrobial resistance

Antimicrobial resistance (AMR) is a major threat to global health. Understanding the emergence, evolution, and transmission of individual antibiotic resistance genes (ARGs) is essential to develop sustainable strategies combatting this threat. Here, we use metagenomic sequencing to analyse ARGs in 757 sewage samples from 243 cities in 101 countries, collected from 2016 to 2019. We find regional patterns in resistomes, and these differ between subsets corresponding to drug classes and are partly driven by taxonomic variation. The genetic environments of 49 common ARGs are highly diverse, with most common ARGs carried by multiple distinct genomic contexts globally and sometimes on plasmids. Analysis of flanking sequence revealed ARG-specific patterns of dispersal limitation and global transmission. Our data furthermore suggest certain geographies are more prone to transmission events and should receive additional attention.

Evidence of pediatric sepsis caused by a drug resistant Lactococcus garvieae contaminated platelet concentrate

Owing to an increasing number of infections in adults, Lactococcus (L.) garvieae has gained recognition as an emerging human pathogen, causing bacteraemia and septicaemia. In September 2020, four paediatric onco-hematologic patients received a platelet concentrate from the same adult donor at Bambino Gesù Children’s Hospital IRCCS, Rome. Three of four patients experienced L. garvieae sepsis one day after transfusion. The L. garvieae pediatric isolates and the donor’s platelet concentrates were retrospectively collected for whole-genome sequencing and shot-gun metagenomics, respectively (Illumina HiSeq). By de novo assembly of the L. garvieae genomes, we found that all three pediatric isolates shared a 99.9% identity and were characterized by 440 common SNPs. Plasmid pUC11C (conferring virulence properties) and the temperate prophage Plg-Tb25 were detected in all three strains. Core SNP genome-based maximum likelihood and Bayesian trees confirmed their phylogenetic common origin and revealed their relationship with L. garvieae strains affecting cows and humans (bootstrap values >100 and posterior probabilities = 1.00). Bacterial reads obtained by the donor’s platelet concentrate have been profiled with MetaPhlAn2 (v.2.7.5); among these, 29.9% belonged to Firmicutes, and 5.16% to Streptococcaceae (>97% identity with L. garvieae), confirming the presence of L. garvieae in the platelet concentrate transfusion. These data showed three episodes of sepsis for the first time due to a transfusion-associated transmission of L. garvieae in three pediatric hospitalized hematology patients. This highlights the importance to implement the screening of platelet components with new human-defined pathogens for ensuring the safety of blood supply, and more broadly, for the surveillance of emerging pathogens.

Plasmids of the urinary microbiota

Studies of the last decade have identified a phylogenetically diverse community of bacteria within the urinary tract of individuals with and without urinary symptoms. Mobile genetic elements (MGEs), including plasmids and phages, within this niche have only recently begun to be explored. These MGEs can expand metabolic capacity and increase virulence, as well as confer antibiotic resistance. As such, they have the potential to contribute to urinary symptoms. While plasmids for some of the bacterial taxa found within the urinary microbiota (urobiome) have been well characterized, many urinary species are under-studied with few genomes sequenced to date. Using a two-pronged bioinformatic approach, we have conducted a comprehensive investigation of the plasmid content of urinary isolates representative of 102 species. The bioinformatic tools plasmidSPAdes and Recycler were used in tandem to identify plasmid sequences from raw short-read sequence data followed by manual curation. In total, we identified 603 high-confidence plasmid sequences in 20 different genera of the urobiome. In total, 70 % of these high-confidence plasmids exhibit sequence similarity to plasmid sequences from the gut. This observation is primarily driven by plasmids from E. coli , which is found in both anatomical niches. To confirm our bioinformatic predictions, long-read sequencing was performed for 23 of the E. coli isolates in addition to two E. coli strains that were sequenced as part of a prior study. Overall, 66.95 % of these predictions were confirmed highlighting the strengths and weaknesses of current bioinformatic tools. Future studies of the urobiome, especially concerning under-studied species in the urobiome, should employ long-read sequencing to expand the catalogue of plasmids for this niche.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Nanopore-based long-read metagenomics uncover the resistome intrusion by antibiotic resistant bacteria from treated wastewater in receiving water body

Wastewater treatment plant (WWTP) effluent discharge could induce the resistome enrichment in the receiving water environments. However, because of the general lack of a robust antibiotic-resistant bacteria (ARB) identification method, the driving mechanism for resistome accumulation in receiving environment is unclear. Here, we took advantage of the enhanced ARBs recognition by nanopore long reads to distinguish the indigenous ARBs and the accumulation of WWTP-borne ARBs in the receiving water body of a domestic WWTP. A bioinformatic framework (named ARGpore2: https://github.com/sustc-xylab/ARGpore2) was constructed and evaluate to facilitate antibiotic resistance genes (ARGs) and ARBs identification in nanopore reads. ARGs identification by ARGpore2 showed comparable precision and recall to that of the commonly adopt BLASTP-based method, whereas the spectrum of ARBs doubled that of the assembled Illumina dataset. Totally, we identified 33 ARBs genera carrying 65 ARG subtypes in the receiving seawater, whose concentration was in general 10 times higher than clean seawater's. Notably we report a primary resistome intrusion caused by the revival of residual microbes survived from disinfection treatment. These WWTP-borne ARBs, including several animal/human enteric pathogens, contributed up to 85% of the receiving water resistome. Plasmids and class 1 integrons were reckoned as major vehicles facilitating the persistence and dissemination of ARGs. Moreover, our work demonstrated the importance of extensive carrier identification in determining the driving force of multifactor coupled resistome booming in complicated environmental conditions, thereby paving the way for establishing priority for effective ARGs mitigation strategies.

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

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

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

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

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

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

Global Distribution and Diversity of Prevalent Sewage Water Plasmidomes

Sewage water from around the world contains an abundance of short plasmids, several of which harbor antimicrobial resistance genes (ARGs). The global dynamics of plasmid-derived antimicrobial resistance and functions are only starting to be unveiled. Here, we utilized a previously created data set of 159,332 assumed small plasmids from 24 different global sewage samples. The detailed phylogeny, as well as the interplay between their protein domains, ARGs, and predicted bacterial host genera, were investigated to understand sewage plasmidome dynamics globally. A total of 58,429 circular elements carried genes encoding plasmid-related features, and MASH distance analyses showed a high degree of diversity. A single (yet diverse) cluster of 520 predicted Acinetobacter plasmids was predominant among the European sewage water. Our results suggested a prevalence of plasmid-backbone gene combinations over others. This could be related to selected bacterial genera that act as bacterial hosts. These combinations also mirrored the geographical locations of the sewage samples. Our functional domain network analysis identified three groups of plasmids. However, these backbone domains were not exclusive to any given group, and Acinetobacter was the dominant host genus among the theta-replicating plasmids, which contained a reservoir of the macrolide resistance gene pair msr(E) and mph(E). Macrolide resistance genes were the most common in the sewage plasmidomes and were found in the largest number of unique plasmids. While msr(E) and mph(E) were limited to Acinetobacter, erm(B) was disseminated among a range of Firmicutes plasmids, including Staphylococcus and Streptococcus, highlighting a potential reservoir of antibiotic resistance for these pathogens from around the globe. IMPORTANCE Antimicrobial resistance is a global threat to human health, as it inhibits our ability to treat infectious diseases. This study utilizes sewage water plasmidomes to identify plasmid-derived features and highlights antimicrobial resistance genes, particularly macrolide resistance genes, as abundant in sewage water plasmidomes in Firmicutes and Acinetobacter hosts. The emergence of macrolide resistance in these bacteria suggests that macrolide selective pressure exists in sewage water and that the resident bacteria can readily acquire macrolide resistance via small plasmids.

Two New Rhizobiales Species Isolated from Root Nodules of Common Sainfoin (Onobrychis viciifolia) Show Different Plant Colonization Strategies

Root nodules of legume plants are primarily inhabited by rhizobial nitrogen-fixing bacteria. Here, we propose two new Rhizobiales species isolated from root nodules of common sainfoin (Onobrychis viciifolia), as shown by core-gene phylogeny, overall genome relatedness indices, and pan-genome analysis. Mesorhizobium onobrychidis sp. nov. actively induces nodules and achieves atmospheric nitrogen and carbon dioxide fixation. This species appears to be depleted in motility genes and is enriched in genes for direct effects on plant growth performance. Its genome reveals functional and plant growth-promoting signatures, like a large unique chromosomal genomic island with high density of symbiotic genetic traits. Onobrychidicola muellerharveyae gen. nov. sp. nov. is described as a type species of the new genus Onobrychidicola in Rhizobiaceae. This species comprises unique genetic features and plant growth-promoting traits (PGPTs), which strongly indicate its function in biotic stress reduction and motility. We applied a newly developed bioinformatics approach for in silico prediction of PGPTs (PGPT-Pred), which supports the different lifestyles of the two new species and the plant growth-promoting performance of M. onobrychidis in the greenhouse trial. IMPORTANCE The intensive use of chemical fertilizers has a variety of negative effects on the environment. Increased utilization of biological nitrogen fixation (BNF) is one way to mitigate those negative impacts. In order to optimize BNF, suitable candidates for different legume species are required. Despite intensive search for new rhizobial bacteria associated with legumes, no new rhizobia have recently been identified from sainfoin (Onobrychis viciifolia). Here, we report on the discovery of two new rhizobial species associated with sainfoin, which are of high importance for the host and may help to increase sustainability in agricultural practices. We employed the combination of in silico prediction and in planta experiments, which is an effective way to detect promising plant growth-promoting bacteria.

Comparative Genomics of Listeria Species Recovered from Meat and Food Processing Facilities

Listeria species (spp.) are contaminants that can survive in food, on equipment, and on food processing premises if appropriate hygiene measures are not used. Homologous stress tolerance genes, virulence gene clusters such as the prfA cluster, and clusters of internalin genes that contribute to the pathogenic potential of the strains can be carried by both pathogenic and nonpathogenic Listeria spp. To enhance understanding of the genome evolution of virulence and virulence-associated properties, a comparative genome approach was used to analyze 41 genome sequences belonging to L. innocua and L. welshimeri isolated from food and food processing facilities. Genetic determinants responsible for disinfectant and stress tolerance were identified, including the efflux cassette bcrABC and Tn6188_qac_1 disinfectant resistance determinant, and stress survival islets. These disinfectant-resistant genes were more frequently found in L. innocua (12%) than in L. welshimeri (2%). Several isolates representing the presumed nonpathogenic L. innocua still carried virulence-associated genes, including LGI2, LGI3, LIPI-3, and LIPI-4 which were absent in all L. welshimeri isolates. The mobile genetic elements identified were plasmids (pLGUG1 and J1776) and prophages (PHAGE_Lister_vB_LmoS_188, PHAGE_Lister_LP_030_3, PHAGE_Lister_A118, PHAGE_Lister_B054, and PHAGE_Lister_vB_LmoS_293). The results suggest that the presumed nonpathogenic isolates especially L. innocua can carry genes relevant to the strain's virulence and stress tolerance in the food and food processing facilities. IMPORTANCE This study provides genomic insights into the recently expanded genus in order to gain valuable information about the evolution of the virulence and stress tolerance properties of the genus Listeria and the distribution of these genetic elements pertinent to the pathogenic potential across Listeria spp. and clonal lineages in South Africa (SA).

Pathogenic strains of Shewanella putrefaciens contain plasmids that are absent in the probiotic strain Pdp11

Shewanella putrefaciens Pdp11 is a strain described as a probiotic for use in aquaculture. However, S. putrefaciens includes strains reported to be pathogenic or saprophytic to fish. Although the probiotic trait has been related to the presence of a group of genes in its genome, the existence of plasmids that could determine the probiotic or pathogenic character of this bacterium is unknown. In the present work, we searched for plasmids in several strains of S. putrefaciens that differ in their pathogenic and probiotic character. Under the different conditions tested, plasmids were only found in two of the five pathogenic strains, but not in the probiotic strain nor in the two saprophytic strains tested. Using a workflow integrating Sanger and Illumina reads, the complete consensus sequences of the plasmids were obtained. Plasmids differed in one ORF and encoded a putative replication initiator protein of the repB family, as well as proteins related to plasmid stability and a toxin-antitoxin system. Phylogenetic analysis showed some similarity to functional repB proteins of other Shewanella species. The implication of these plasmids in the probiotic or pathogenic nature of S. putrefaciens is discussed.

Comparative genome analysis of mycobacteria focusing on tRNA and non-coding RNA

BACKGROUND

The Mycobacterium genus encompasses at least 192 named species, many of which cause severe diseases such as tuberculosis. Non-tuberculosis mycobacteria (NTM) can also infect humans and animals. Some are of emerging concern because they show high resistance to commonly used antibiotics while others are used and evaluated in bioremediation or included in anticancer vaccines.

RESULTS

We provide the genome sequences for 114 mycobacterial type strains and together with 130 available mycobacterial genomes we generated a phylogenetic tree based on 387 core genes and supported by average nucleotide identity (ANI) data. The 244 genome sequences cover most of the species constituting the Mycobacterium genus. The genome sizes ranged from 3.2 to 8.1 Mb with an average of 5.7 Mb, and we identified 14 new plasmids. Moreover, mycobacterial genomes consisted of phage-like sequences ranging between 0 and 4.64% dependent on mycobacteria while the number of IS elements varied between 1 and 290. Our data also revealed that, depending on the mycobacteria, the number of tRNA and non-coding (nc) RNA genes differ and that their positions on the chromosome varied. We identified a conserved core set of 12 ncRNAs, 43 tRNAs and 18 aminoacyl-tRNA synthetases among mycobacteria.

CONCLUSIONS

Phages, IS elements, tRNA and ncRNAs appear to have contributed to the evolution of the Mycobacterium genus where several tRNA and ncRNA genes have been horizontally transferred. On the basis of our phylogenetic analysis, we identified several isolates of unnamed species as new mycobacterial species or strains of known mycobacteria. The predicted number of coding sequences correlates with genome size while the number of tRNA, rRNA and ncRNA genes does not. Together these findings expand our insight into the evolution of the Mycobacterium genus and as such they establish a platform to understand mycobacterial pathogenicity, their evolution, antibiotic resistance/tolerance as well as the function and evolution of ncRNA among mycobacteria.

A low-cost genomics workflow enables isolate screening and strain-level analyses within microbiomes

Earth's environments harbor complex consortia of microbes that affect processes ranging from host health to biogeochemical cycles. Understanding their evolution and function is limited by an inability to isolate genomes in a high-throughput manner. Here, we present a workflow for bacterial whole-genome sequencing using open-source labware and the OpenTrons robotics platform, reducing costs to approximately $10 per genome. We assess genomic diversity within 45 gut bacterial species from wild-living chimpanzees and bonobos. We quantify intraspecific genomic diversity and reveal divergence of homologous plasmids between hosts. This enables population genetic analyses of bacterial strains not currently possible with metagenomic data alone.

Effects of Nutrient Level and Growth Rate on the Conjugation Process That Transfers Mobile Antibiotic Resistance Genes in Continuous Cultures

Bacteria in the effluent of wastewater treatment plants (WWTPs) can transfer antibiotic resistance genes (ARGs) to the bacteria in receiving water through conjugation; however, there is a lack of quantitative assessment of this phenomenon in continuous cultures. Our objective was to determine the effects of background nutrient levels in river water column and growth rates of bacteria on the conjugation frequency of ARGs from effluent bacteria to river bacteria, as well as on the resulting resistance level (i.e., MICs) of the river bacteria. Chemostats were employed to simulate the discharge points of WWTPs into rivers, where effluent bacteria (donor cells) meet river bacteria (recipient cells). Both donor and recipient cells were Escherichia coli cells, and the donor cells were constructed by filter mating with bacteria in the effluent of a local WWTP. Results showed that higher bacterial growth rate (0.45 h^-1 versus 0.15 h^-1) led to higher conjugation frequencies (10^-4 versus 10^-6 transconjugant per recipient). The nutrient level also significantly affected the conjugation frequency, albeit to a lesser extent than the growth rate. The MIC against tetracycline increased from 2 mg/L in the recipient to 64 to 128 mg/L in transconjugants. In comparison, the MIC only increased to as high as 8 mg/L in mutants. Whole-genome sequencing showed that the tet-containing plasmid in both the donor and the transconjugant cells also occur in other fecal bacterial genera. The quantitative information obtained from this study can inform hazard identification related to the proliferation of wastewater-associated ARGs in surface water. IMPORTANCE WWTPs have been regarded as an important hot spot of ARGs. The discharge point of WWTP effluent, where ARGs may be horizontally transferred from bacteria of treated wastewater to bacteria of receiving water, is an important interface between the human-dominated ecosystem and the natural environment. The use of batch cultures in previous studies cannot adequately simulate the nutrient conditions and growth rates in receiving water. In this study, chemostats were employed to simulate the continuous growth of bacteria in receiving water. Furthermore, the experimental setup allowed for separate investigations on the effects of nutrient levels (i.e., simulating background nutrients in river water) and bacterial growth rates on conjugation frequencies and resulting resistance levels. The study generates statistically sound ecological data that can be used to estimate the risk of wastewater-originated ARGs as part of the One Health framework.

Revealing bacteria-phage interactions in human microbiome through the CRISPR-Cas immune systems

The human gut microbiome is composed of a diverse consortium of microorganisms. Relatively little is known about the diversity of the bacteriophage population and their interactions with microbial organisms in the human microbiome. Due to the persistent rivalry between microbial organisms (hosts) and phages (invaders), genetic traces of phages are found in the hosts' CRISPR-Cas adaptive immune system. Mobile genetic elements (MGEs) found in bacteria include genetic material from phage and plasmids, often resultant from invasion events. We developed a computational pipeline (BacMGEnet), which can be used for inference and exploratory analysis of putative interactions between microbial organisms and MGEs (phages and plasmids) and their interaction network. Given a collection of genomes as the input, BacMGEnet utilizes computational tools we have previously developed to characterize CRISPR-Cas systems in the genomes, which are then used to identify putative invaders from publicly available collections of phage/prophage sequences. In addition, BacMGEnet uses a greedy algorithm to summarize identified putative interactions to produce a bacteria-MGE network in a standard network format. Inferred networks can be utilized to assist further examination of the putative interactions and for discovery of interaction patterns. Here we apply the BacMGEnet pipeline to a few collections of genomic/metagenomic datasets to demonstrate its utilities. BacMGEnet revealed a complex interaction network of the Phocaeicola vulgatus pangenome with its phage invaders, and the modularity analysis of the resulted network suggested differential activities of the different P. vulgatus' CRISPR-Cas systems (Type I-C and Type II-C) against some phages. Analysis of the phage-bacteria interaction network of human gut microbiome revealed a mixture of phages with a broad host range (resulting in large modules with many bacteria and phages), and phages with narrow host range. We also showed that BacMGEnet can be used to infer phages that invade bacteria and their interactions in wound microbiome. We anticipate that BacMGEnet will become an important tool for studying the interactions between bacteria and their invaders for microbiome research.

Nodules of wild legumes as unique natural hotspots of antibiotic resistance genes

Root nodules (RN) of legumes have distinct microenvironment from their symbiotic roots and surrounding soils. The rhizobia can withstand the host-produced phytoalexins and antimicrobial compounds. We thus hypothesize that the wild legume RN may develop unique natural resistome and be antibiotic resistance gene (ARG) hotspots. In this study, in comparison with rhizosphere soil (RS) and bulk soil (BS), we characterized the feature of antibiotic resistance in the RN of two wild legumes, Medicago polymorpha and Astragalus sinicus, by metagenomics. It was shown that the total relative abundance of ARGs followed the order of RN > RS > BS for both legumes. ARGs encoding antibiotic efflux pump predominated in all samples with increased proportion from BS to RN samples for both legumes. Totally 275 ARG subtypes were detected, and diversity of ARGs in RN was significantly lower than in BS samples for both legumes. 32 and 25 unique ARGs subtypes were detected in RN of both legumes. Bacterial community played a key role in shaping nodule-associated resistome because both ARG profiles and bacterial community differed greatly among BS, RS and RN. Rhizobia potentially hosted 10 and 15 ARGs subtypes for both legumes. The number and proportion of plasmid- and ARG-carrying contigs (ACCs) were higher in RN than in BS. Host tracking analysis of plasmid-ACCs suggests that proportion of rhizobial bacteria identified as their hosts decreased from BS to RN samples. No plasmid-ACCs with multiple ARGs were observed in BS samples, whereas they were detected in RN samples of both legumes. Our study showed that even wild legume nodules are unique natural ARG hotspots and enough attention should be paid to the dissemination risk of ARGs posed by globally produced legume crops.

Long-term ecological and evolutionary dynamics in the gut microbiomes of carbapenemase-producing Enterobacteriaceae colonized subjects

Long-term colonization of the gut microbiome by carbapenemase-producing Enterobacteriaceae (CPE) is a growing area of public health concern as it can lead to community transmission and rapid increase in cases of life-threatening CPE infections. Here, leveraging the observation that many subjects are decolonized without interventions within a year, we used longitudinal shotgun metagenomics (up to 12 timepoints) for detailed characterization of ecological and evolutionary dynamics in the gut microbiome of a cohort of CPE-colonized subjects and family members (n = 46; 361 samples). Subjects who underwent decolonization exhibited a distinct ecological shift marked by recovery of microbial diversity, key commensals and anti-inflammatory pathways. In addition, colonization was marked by elevated but unstable Enterobacteriaceae abundances, which exhibited distinct strain-level dynamics for different species (Escherichia coli and Klebsiella pneumoniae). Finally, comparative analysis with whole-genome sequencing data from CPE isolates (n = 159) helped identify substrain variation in key functional genes and the presence of highly similar E. coli and K. pneumoniae strains with variable resistance profiles and plasmid sharing. These results provide an enhanced view into how colonization by multi-drug-resistant bacteria associates with altered gut ecology and can enable transfer of resistance genes, even in the absence of overt infection and antibiotic usage.

Assessing the genomic composition, putative ecological relevance and biotechnological potential of plasmids from sponge bacterial symbionts

Plasmid-mediated transfer of genes can have direct consequences in several biological processes within sponge microbial communities. However, very few studies have attempted genomic and functional characterization of plasmids from marine host-associated microbial communities in general and those of sponges in particular. In the present study, we used an endogenous plasmid isolation method to obtain plasmids from bacterial symbionts of the marine sponges Stylissa carteri and Paratetilla sp. and investigated the genomic composition, putative ecological relevance and biotechnological potential of these plasmids. In total, we isolated and characterized three complete plasmids, three plasmid prophages and one incomplete plasmid. Our results highlight the importance of plasmids to transfer relevant genetic traits putatively involved in microbial symbiont adaptation and host-microbe and microbe-microbe interactions. For example, putative genes involved in bacterial response to chemical stress, competition, metabolic versatility and mediation of bacterial colonization and pathogenicity were detected. Genes coding for enzymes and toxins of biotechnological potential were also detected. Most plasmid prophage coding sequences were, however, hypothetical proteins with unknown functions. Overall, this study highlights the ecological relevance of plasmids in the marine sponge microbiome and provides evidence that plasmids of sponge bacterial symbionts may represent an untapped resource of genes of biotechnological interest.

Methicillin Resistant Staphylococci Isolated from Goats and Their Farm Environments in Saudi Arabia Genotypically Linked to Known Human Clinical Isolates: a Pilot Study

We conducted a pilot whole genome sequencing (WGS) study to characterize the genotypes of nine methicillin resistant staphylococci (MRS) isolates recovered from goats and their farm environments in Eastern Province, Saudi Arabia, between November 2019 to August 2020. Seven out of nine isolates were methicillin resistant Staphylococcus aureus (MRSA), and two were methicillin resistant Staphylococcus epidermidis (MRSE). All MRSA isolates possessed genotypes previously identified to infect humans, including isolates harboring ST6-SCCmec IV-t304 (n = 4), ST5-SCCmec VI- t688 (n = 2) and ST5-SCCmec V-t311 (n = 1). 2 MRSA isolates possessed plasmids that were genetically similar to those identified in S. aureus isolates recovered from humans and poultry. In contrast, plasmids found in three MRSA isolates and one MRSE isolate were genetically similar to those recovered from humans. All MRSA isolates harbored the host innate modulate genes sak and scn previously associated with human infections. The genotypes of MRSE isolates were determined as ST35, a well-known zoonotic sequence type and ST153, which has been associated with humans. However, the MRSE isolates were untypeable due to extra ccr complexes identified in their SCCmec elements. Moreover, we identified in ST153 isolate SCCmec element also harbored the Arginine Catabolic Mobile Element (ACME) IV. All MRS isolates were phenotypically resistant to trimethoprim-sulfamethoxazole, an antibiotic for the decolonization of MRS. Three isolates carried antibiotic resistance genes in their SCCmec elements that were not previously described, including those encoding fusidic acid resistance (fusC) and trimethoprim resistance (dfrC) incorporated in the MRSA SCCmec VI. IMPORTANCE Our findings demonstrate a possible cross-transmission of methicillin resistant staphylococci between goats and their local environments and between goats and humans. Due to ever increasing resistance to multiple antibiotics, the burden of MRS has a significant impact on livestock farming, public health, and the economy worldwide. This study highlights that implementing a holistic approach to whole genome sequencing surveillance in livestock and farm environments would aid our understanding of the transmission of methicillin resistant staphylococci and, most importantly, allow us to implement appropriate infection control and hygiene practices.

Extensive diversity and rapid turnover of phage defense repertoires in cheese-associated bacterial communities

BACKGROUND

Phages are key drivers of genomic diversity in bacterial populations as they impose strong selective pressure on the evolution of bacterial defense mechanisms across closely related strains. The pan-immunity model suggests that such diversity is maintained because the effective immune system of a bacterial species is the one distributed across all strains present in the community. However, only few studies have analyzed the distribution of bacterial defense systems at the community-level, mostly focusing on CRISPR and comparing samples from complex environments. Here, we studied 2778 bacterial genomes and 188 metagenomes from cheese-associated communities, which are dominated by a few bacterial taxa and occur in relatively stable environments.

RESULTS

We corroborate previous laboratory findings that in cheese-associated communities nearly identical strains contain diverse and highly variable arsenals of innate and adaptive (i.e., CRISPR-Cas) immunity systems suggesting rapid turnover. CRISPR spacer abundance correlated with the abundance of matching target sequences across the metagenomes providing evidence that the identified defense repertoires are functional and under selection. While these characteristics align with the pan-immunity model, the detected CRISPR spacers only covered a subset of the phages previously identified in cheese, providing evidence that CRISPR does not enable complete immunity against all phages, and that the innate immune mechanisms may have complementary roles.

CONCLUSIONS

Our findings show that the evolution of bacterial defense mechanisms is a highly dynamic process and highlight that experimentally tractable, low complexity communities such as those found in cheese, can help to understand ecological and molecular processes underlying phage-defense system relationships. These findings can have implications for the design of robust synthetic communities used in biotechnology and the food industry. Video Abstract.

Diversity and dynamics of the CRISPR-Cas systems associated with Bacteroides fragilis in human population

Background

CRISPR-Cas (clustered regularly interspaced short palindromic repeats—CRISPR-associated proteins) systems are adaptive immune systems commonly found in prokaryotes that provide sequence-specific defense against invading mobile genetic elements (MGEs). The memory of these immunological encounters are stored in CRISPR arrays, where spacer sequences record the identity and history of past invaders. Analyzing such CRISPR arrays provide insights into the dynamics of CRISPR-Cas systems and the adaptation of their host bacteria to rapidly changing environments such as the human gut.

Results

In this study, we utilized 601 publicly available Bacteroides fragilis genome isolates from 12 healthy individuals, 6 of which include longitudinal observations, and 222 available B. fragilis reference genomes to update the understanding of B. fragilis CRISPR-Cas dynamics and their differential activities. Analysis of longitudinal genomic data showed that some CRISPR array structures remained relatively stable over time whereas others involved radical spacer acquisition during some periods, and diverse CRISPR arrays (associated with multiple isolates) co-existed in the same individuals with some persisted over time. Furthermore, features of CRISPR adaptation, evolution, and microdynamics were highlighted through an analysis of host-MGE network, such as modules of multiple MGEs and hosts, reflecting complex interactions between B. fragilis and its invaders mediated through the CRISPR-Cas systems.

Conclusions

We made available of all annotated CRISPR-Cas systems and their target MGEs, and their interaction network as a web resource at https://omics.informatics.indiana.edu/CRISPRone/Bfragilis. We anticipate it will become an important resource for studying of B. fragilis, its CRISPR-Cas systems, and its interaction with mobile genetic elements providing insights into evolutionary dynamics that may shape the species virulence and lead to its pathogenicity.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08770-8).

Antibiotic resistance genes in the gut microbiota of mothers and linked neonates with or without sepsis from low- and middle-income countries

Early development of the microbiome has been shown to affect general health and physical development of the infant and, although some studies have been undertaken in high-income countries, there are few studies from low- and middle-income countries. As part of the BARNARDS study, we examined the rectal microbiota of 2,931 neonates (term used up to 60 d) with clinical signs of sepsis and of 15,217 mothers screening for bla_CTX-M-15, bla_NDM, bla_KPC and bla_OXA-48-like genes, which were detected in 56.1%, 18.5%, 0% and 4.1% of neonates’ rectal swabs and 47.1%, 4.6%, 0% and 1.6% of mothers’ rectal swabs, respectively. Carbapenemase-positive bacteria were identified by MALDI-TOF MS and showed a high diversity of bacterial species (57 distinct species/genera) which exhibited resistance to most of the antibiotics tested. Escherichia coli, Klebsiella pneumoniae and Enterobacter cloacae/E. cloacae complex, the most commonly found isolates, were subjected to whole-genome sequencing analysis and revealed close relationships between isolates from different samples, suggesting transmission of bacteria between neonates, and between neonates and mothers. Associations between the carriage of antimicrobial resistance genes (ARGs) and healthcare/environmental factors were identified, and the presence of ARGs was a predictor of neonatal sepsis and adverse birth outcomes.

Analysis of gut microbiota of mothers and its neonates—as part of the BARNARDS study—reveals associations between β-lactamase gene carriage and neonatal sepsis risk in low-income settings.

Shooting hoops: globetrotting plasmids spreading more than just antimicrobial resistance genes across One Health

Our study provides novel insights into the global nature of antimicrobial resistance (AMR) plasmids across the food chain. We provide compelling evidence of the globetrotting nature of AMR plasmids and the need for surveillance to sequence plasmids with a template of analyses for others to expand these data. The AMR plasmids analysed were detected in 63 countries and in samples from humans, animals and the environment. They contained a combination of known and novel AMR genes, metal resistance genes, virulence factors, phage and replicon types.

Draft-genome analysis provides insights into the virulence properties and genome plasticity of Vibrio fluvialis organisms isolated from shrimp farms and Turag river in Bangladesh

Vibrio fluvialis is an opportunistic waterborne and seafood-borne enteric pathogen capable of causing severe diarrhea leading to death. This pathogen is endemic to Bangladesh, a country which is a major producer of cultured shrimp and wild-caught prawns. In this study, we carried out whole-genome sequencing of three V. fluvialis organisms isolated from shrimp farm and river sediment showing strong pathogenic characteristics in vivo and in vitro and compared their genomes against other V. fluvialis and related pathogenic species to glean insights into their potential as pathogens. Numerous virulence-associated genes including hemolysins, cytolysins, three separate Type IV pili, Types II and VI secretion systems, biofilm, and the V. cholerae pathogenesis regulating gene, toxR, were identified. Moreover, we found strain S-10 to have the propensity to acquire antibiotic resistance genes through horizontal gene transfer. These findings indicate that shrimp farms and rivers could be potential sources of V. fluvialis organisms which are an infection threat of public health concern.

Revealing Genomic Insights of the Unexplored Porcine Pathogen Actinobacillus pleuropneumoniae Using Whole Genome Sequencing

Actinobacillus pleuropneumoniae (APP) is the causative agent of pleuropneumonia in pigs, one of the most relevant bacterial respiratory diseases in the swine industry. To date, 19 serotypes have been described based on capsular polysaccharide typing with significant virulence dissimilarities. In this study, 16 APP isolates from Spanish origin were selected to perform antimicrobial susceptibility tests and comparative genomic analysis using whole genome sequencing (WGS). To obtain a more comprehensive worldwide molecular epidemiologic analyses, all APP whole genome assemblies available at the National Center for Biotechnology Information (NCBI) at the time of the study were also included. An in-house in silico PCR approach enabled the correct serotyping of unserotyped or incorrectly serotyped isolates and allowed for the discrimination between serotypes 9 and 11. A pangenome analysis identified the presence or absence of gene clusters to be serotype specific, as well as virulence profile analyses targeting the apx operons. Antimicrobial resistance genes were correlated to the presence of specific plasmids. Altogether, this study provides new insights into the genetic variability within APP serotypes, correlates phenotypic tests with bioinformatic analyses and manifests the benefits of populated databases for a better assessment of diversity and variability of relatively unknown pathogens. Overall, genomic comparative analysis enhances the understanding of transmission and epidemiological patterns of this species and suggests vertical transmission of the pathogen, including the resistance genes, within the Spanish integrated systems.

IMPORTANCE Pleuropneumonia is one of the most relevant respiratory infections in the swine industry. Despite Actinobacillus pleuropneumoniae (APP) being one of the most important pathogens in the pig production, this is the first comparative study including all available whole genome sequencing data from NCBI. Moreover, this study also includes 16 APP isolates of Spanish origin with known epidemiological relationships through vertical integrated systems. Genomic comparisons provided a deeper understanding of molecular and epidemiological knowledge between different APP serotypes. Furthermore, determination of resistance and toxin profiles allowed correlation with the presence of mobile genetic elements and specific serotype, respectively.

Genomic Analysis of an I1 Plasmid Hosting a sul3-Class 1 Integron and blaSHV-12 within an Unusual Escherichia coli ST297 from Urban Wildlife

Wild birds, particularly silver gulls (Chroicocephalus novaehollandiae) that nest near anthropogenic sites, often harbour bacteria resistant to multiple antibiotics, including those considered of clinical importance. Here, we describe the whole genome sequence of Escherichia coli isolate CE1867 from a silver gull chick sampled in 2012 that hosted an I1 pST25 plasmid with bla_SHV-12, a β-lactamase gene that encodes the ability to hydrolyze oxyimino β-lactams, and other antibiotic resistance genes. Isolate CE1867 is an ST297 isolate, a phylogroup B1 lineage, and clustered with a large ST297 O130:H11 clade, which carry Shiga toxin genes. The I1 plasmid belongs to plasmid sequence type 25 and is notable for its carriage of an atypical sul3-class 1 integron with mefB_∆260, a structure most frequently reported in Australia from swine. This integron is a typical example of a Tn21-derived element that captured sul3 in place of the standard sul1 structure. Interestingly, the mercury resistance (mer) module of Tn21 is missing and has been replaced with Tn2-bla_TEM-1 and a bla_SHV-12 encoding module flanked by direct copies of IS26. Comparisons to similar plasmids, however, demonstrate a closely related family of ARG-carrying plasmids that all host variants of the sul3-associated integron with conserved Tn21 insertion points and a variable presence of both mer and mefB truncations, but predominantly mefB_∆260.

UG/Abi: a highly diverse family of prokaryotic reverse transcriptases associated with defense functions

Reverse transcriptases (RTs) are enzymes capable of synthesizing DNA using RNA as a template. Within the last few years, a burst of research has led to the discovery of novel prokaryotic RTs with diverse antiviral properties, such as DRTs (Defense-associated RTs), which belong to the so-called group of unknown RTs (UG) and are closely related to the Abortive Infection system (Abi) RTs. In this work, we performed a systematic analysis of UG and Abi RTs, increasing the number of UG/Abi members up to 42 highly diverse groups, most of which are predicted to be functionally associated with other gene(s) or domain(s). Based on this information, we classified these systems into three major classes. In addition, we reveal that most of these groups are associated with defense functions and/or mobile genetic elements, and demonstrate the antiphage role of four novel groups. Besides, we highlight the presence of one of these systems in novel families of human gut viruses infecting members of the Bacteroidetes and Firmicutes phyla. This work lays the foundation for a comprehensive and unified understanding of these highly diverse RTs with enormous biotechnological potential.

High Prevalence of bla CTXM-1/IncI1-Iγ/ST3 Plasmids in Extended-Spectrum β-Lactamase-Producing Escherichia coli Isolates Collected From Domestic Animals in Guadeloupe (French West Indies)

Extended-spectrum β-lactamase-producing Enterobacteriaceae (ESBL-E) have been classified in the group of resistant bacteria of highest priority. We determined the prevalence of ESBL-E collected in feces from household and shelter pets in Guadeloupe (French West Indies). A single rectal swab was taken from 125 dogs and 60 cats between June and September 2019. The prevalence of fecal carriage of ESBL-E was 7.6% (14/185, 95% CI: 4.2-12.4), within the range observed worldwide. The only risk factor associated with a higher prevalence of ESBL-E rectal carriage was a stay in a shelter, suggesting that refuges could be hotspots for their acquisition. All but one (Klebsiella pneumoniae from a cat) were Escherichia coli. We noted the presence of a bla _CTX-M-1/IncI1-Iγ/sequence type (ST3) plasmid in 11 ESBL-producing E. coli isolates belonging to ST328 (n = 6), ST155 (n = 4) and ST953 (n = 1). A bla _CTX-M-15 gene was identified in the three remaining ESBL-E isolates. The bla _CTX-M-1 and most of the antimicrobial resistance genes were present in a well-conserved large conjugative IncI1-Iγ/ST3 plasmid characterized by two accessory regions containing antibiotic resistance genes. The plasmid has been detected worldwide in E. coli isolates from humans and several animal species, such as food-producing animals, wild birds and pets, and from the environment. This study shows the potential role of pets as a reservoir of antimicrobial-resistant bacteria or genes for humans and underlines the importance of basic hygiene measures by owners of companion animals.

Virulence plasmid pINV as a genetic signature for Shigella flexneri phylogeny

Shigella flexneri is a major health burden in low- and middle-income countries, where it is a leading cause of mortality associated with diarrhoea in children, and shows an increasing incidence among travellers and men having sex with men. Like all Shigella spp., S. flexneri has evolved from commensal Escherichia coli following the acquisition of a large plasmid pINV, which contains genes essential for virulence. Current sequence typing schemes of Shigella are based on combinations of chromosomal genetic loci, since pINV-encoded virulence genes are often lost during growth in the laboratory, making these elements inappropriate for sequence typing. By performing comparative analysis of pINVs from S. flexneri strains isolated from different geographical regions and belonging to different serotypes, we found that in contrast to plasmid-encoded virulence genes, plasmid maintenance genes are highly stable pINV-encoded elements. For the first time, to our knowledge, we have developed a S. flexneri plasmid multilocus sequence typing (pMLST) method based on different combinations of alleles of the vapBC and yacAB toxin–antitoxin (TA) systems, and the parAB partitioning system. This enables typing of S. flexneri pINV plasmids into distinct ‘virulence sequence types’ (vSTs). Furthermore, the phylogenies of vST alleles and bacterial host core genomes suggests an intimate co-evolution of pINV with the chromosome of its bacterial host, consistent with previous findings. This work demonstrates the potential of plasmid maintenance loci as genetic characteristics to study as well as to trace the molecular phylogenesis of S. flexneri pINV and the phylogenetic relationship of this plasmid with its bacterial host.

Phenotypic characterization and analysis of complete genomes of two distinct strains of the proposed species "L. swaminathanii"

Recently, a new Listeria species, “Listeria swaminathanii”, was proposed. Here, we phenotypically and genotypically characterize two additional strains that were previously obtained from soil samples and compare the results to the type strain. Complete genomes for both strains were assembled from hybrid Illumina and Nanopore sequencing reads and annotated. Further genomic analysis including average nucleotide identity (ANI) and detection of mobile genetic elements and genes of interest (e.g., virulence-associated) were conducted. The strains showed 98.7–98.8% ANI with the type strain. The UTK C1-0015 genome contained a partial monocin locus and a plasmid, while the UTK C1-0024 genome contained a full monocin locus and a prophage. Phenotypic characterization consistent with those performed on the proposed type strain was conducted to assess consistency of phenotypes across a greater diversity of the proposed species (n = 3 instead of n = 1). Only a few findings were notably different from those of the type strain, such as catalase activity, glycerol metabolism, starch metabolism, and growth at 41 °C. This study further expands our understanding of this newly proposed sensu stricto Listeria species.

CRISPR-Cas systems are widespread accessory elements across bacterial and archaeal plasmids

Many prokaryotes encode CRISPR-Cas systems as immune protection against mobile genetic elements (MGEs), yet a number of MGEs also harbor CRISPR-Cas components. With a few exceptions, CRISPR-Cas loci encoded on MGEs are uncharted and a comprehensive analysis of their distribution, prevalence, diversity, and function is lacking. Here, we systematically investigated CRISPR-Cas loci across the largest curated collection of natural bacterial and archaeal plasmids. CRISPR-Cas loci are widely but heterogeneously distributed across plasmids and, in comparison to host chromosomes, their mean prevalence per Mbp is higher and their distribution is distinct. Furthermore, the spacer content of plasmid CRISPRs exhibits a strong targeting bias towards other plasmids, while chromosomal arrays are enriched with virus-targeting spacers. These contrasting targeting preferences highlight the genetic independence of plasmids and suggest a major role for mediating plasmid-plasmid conflicts. Altogether, CRISPR-Cas are frequent accessory components of many plasmids, which is an overlooked phenomenon that possibly facilitates their dissemination across microbiomes.

Whole genome sequencing of Klebsiella pneumoniae clinical isolates sequence type 627 isolated from Egyptian patients

Klebsiella pneumoniae is considered a threat to public health especially due to multidrug resistance emergence. It is largely oligoclonal based on multi-locus sequence typing (MLST); in Egypt, ST 627 was recently detected. Despites the global dissemination of this ST, there is still paucity of information about it. Herein, we used 4 K. pneumoniae ST627 for whole genome sequencing utilizing an Illumina MiSeq platform. Genome sequences were examined for resistance and virulence determinants, capsular types, plasmids, insertion sequences, phage regions, and Clustered Regularly Interspaced Palindromic Repeats (CRISPR) regions using bioinformatic analysis. The molecular characterization revealed 15 and 65 antimicrobial resistance and virulence genes, respectively. Resistance genes such as tet(D), aph(3’’)-Ib, aph(6)-Id, blaTEM-234, fosA, and fosA6; were mainly responsible for tetracycline, aminoglycoside, and fosfomycin resistance; respectively. The capsular typing revealed that the four strains are KL-24 and O1v1. One plasmid was found in all samples known as pC17KP0052-1 and another plasmid with accession no. NZ_CP032191.1 was found only in K90. IncFIB(K) and IncFII(K) are two replicons found in all samples, while ColRNAI replicon was found only in K90. Entero P88, Salmon SEN5, and Klebsi phiKO2 intact phage regions were identified. All samples harbored CRISPR arrays including CRISPR1 and CRISPR2. Our results shed light on critical tasks of mobile genetic elements in ST 627 in antibiotic resistance spreading.

Emergence of Neonatal Sepsis Caused by MCR-9- and NDM-1-Co-Producing Enterobacter hormaechei in China

Mobile colistin resistance (mcr) genes represent an emerging threat to public health. Reports on the prevalence, antimicrobial profiles, and clonality of MCR-9-producing Enterobacter cloacae complex (ECC) isolates on a national scale in China are limited. We screened 3,373 samples from humans, animals, and the environment and identified eleven MCR-9-positive ECC isolates. We further investigated their susceptibility, epidemiology, plasmid profiles, genetic features, and virulence potential. Ten strains were isolated from severe bloodstream infection cases, especially three of them were recovered from neonatal sepsis. Enterobacter hormaechei was the most predominant species among the MCR-9-producing ECC population. Moreover, the co-existence of MCR-9, CTX-M, and SHV-12 encoding genes in MCR-9-positive isolates was globally observed. Notably, mcr-9 was mainly carried by IncHI2 plasmids, and we found a novel ~187 kb IncFII plasmid harboring mcr-9, with low similarity with known plasmids. In summary, our study presented genomic insights into genetic characteristics of MCR-9-producing ECC isolates retrieved from human, animal, and environment samples with one health perspective. This study is the first to reveal NDM-1- and MCR-9-co-producing ECC from neonatal sepsis in China. Our data highlights the risk for the hidden spread of the mcr-9 colistin resistance gene.

Metagenomics Insights Into the Microbial Diversity and Microbiome Network Analysis on the Heterogeneity of Influent to Effluent Water

Sanitizing the water sources of local communities is important to control the spread of microbial resistance genes, especially those for water-borne illnesses. The activities of antibiotic resistance gene (ARG)-host pathogens pose a threat to public health, and it has been estimated that the infection will lead up to 10 million deaths globally by the year 2050. Hence, in this study, we aim to analyze the efficiency of our municipal wastewater treatment plant (WWTP) process in producing pathogen-free water by investigating the microbial composition between influent and effluent water sites. Shotgun metagenomics sequencing using the Illumina platform was performed on the influent and effluent samples of six different WWTP sites located in Johore, Malaysia. After raw data pre-processing, the non-redundant contigs library was then aligned against BLASTP for taxonomy profiling and the Comprehensive Antibiotic Resistance Database for ARG annotation. Interestingly, the alpha-diversity result reported that effluent site samples showed higher abundance and diverse heterogeneity compared to the influent site. The principal component analysis (PCA) and non-metric multidimensional scaling (NMDS) plots also suggested that effluent sites showed high variation in the genetic material due to loosely clustered sample plots, as compared to the tightly clustered influent samples. This study has successfully identified the top three abundant phyla in influent-Proteobacteria, Firmicutes, and Bacteroidetes-and effluent-Proteobacteria, Actinobacteria, and Bacteroidetes-water. Despite the overlap within the top three abundant phyla in influent and effluent sites (Proteobacteria and Bacteroidetes), the ARG composition heat map and drug class phenotype plot bar exhibits a general trend of a downward shift, showing the efficiency of WWTP in reducing opportunistic pathogens. Overall, it was demonstrated that our municipal WWTP efficiently eliminated pathogenic microbes from the influent water before its total discharge to the environment, though not with the total elimination of microorganisms. This metagenomics study allowed for an examination of our water source and showed the potential interaction of species and ARGs residing in the influent and effluent environment. Both microbial profile structure and co-occurrence network analysis provide integrated understanding regarding the diversity of microorganisms and interactions for future advanced water sanitation treatments.

In-depth genome analysis of Bacillus sp. BH32, a salt stress-tolerant endophyte obtained from a halophyte in a semiarid region

Endophytic strains belonging to the Bacillus cereus group were isolated from the halophytes Atriplex halimus L. (Amaranthaceae) and Tamarix aphylla L. (Tamaricaceae) from costal and continental regions in Algeria. Based on their salt tolerance (up to 5%), the strains were tested for their ability to alleviate salt stress in tomato and wheat. Bacillus sp. strain BH32 showed the highest potential to reduce salinity stress (up to + 50% and + 58% of dry weight improvement, in tomato and wheat, respectively, compared to the control). To determine putative mechanisms involved in salt tolerance and plant growth promotion, the whole genome of Bacillus sp. BH32 was sequenced, annotated, and used for comparative genomics against the genomes of closely related strains. The pangenome of Bacillus sp. BH32 and its closest relative was further analyzed. The phylogenomic analyses confirmed its taxonomic position, a member of the Bacillus cereus group, with intergenomic distances (GBDP analysis) pinpointing to a new taxon (digital DNA-DNA hybridization, dDDH < 70%). Genome mining unveiled several genes involved in stress tolerance, production of anti-oxidants and genes involved in plant growth promotion as well as in the production of secondary metabolites. KEY POINTS : • Bacillus sp. BH32 and other bacterial endophytes were isolated from halophytes, to be tested on tomato and wheat and to limit salt stress adverse effects. • The strain with the highest potential was then studied at the genomic level to highlight numerous genes linked to plant growth promotion and stress tolerance. • Pangenome approaches suggest that the strain belongs to a new taxon within the Bacillus cereus group.

Use of Nanopore Sequencing to Characterise the Genomic Architecture of Mobile Genetic Elements Encoding bla CTX-M-15 in Escherichia coli Causing Travellers' Diarrhoea

Increasing levels of antimicrobial resistance (AMR) have been documented in Escherichia coli causing travellers’ diarrhoea, particularly to the third-generation cephalosporins. Diarrhoeagenic E. coli (DEC) can act as a reservoir for the exchange of AMR genes between bacteria residing in the human gut, enabling them to survive and flourish through the selective pressures of antibiotic treatments. Using Oxford Nanopore Technology (ONT), we sequenced eight isolates of DEC from four patients’ specimens who had all recently returned to the United Kingdome from Pakistan. Sequencing yielded two DEC harbouring bla_CTX-M-15 per patient, all with different sequence types (ST) and belonging to five different pathotypes. The study aimed to determine whether bla_CTX-M-15 was located on the chromosome or plasmid and to characterise the drug-resistant regions to better understand the mechanisms of onward transmission of AMR determinants. Patients A and C both had one isolate where bla_CTX-M-15 was located on the plasmid (899037 & 623213, respectively) and one chromosomally encoded (899091 & 623214, respectively). In patient B, bla_CTX-M-15 was plasmid-encoded in both DEC isolates (786605 & 7883090), whereas in patient D, bla_CTX-M-15 was located on the chromosome in both DEC isolates (542093 & 542099). The two bla_CTX-M-15-encoding plasmids associated with patient B were different although the bla_CTX-M-15-encoding plasmid isolated from 788309 (IncFIB) exhibited high nucleotide similarity to the bla_CTX-M-15-encoding plasmid isolated from 899037 (patient A). In the four isolates where bla_CTX-M-15 was chromosomally encoded, two isolates (899091 & 542099) shared the same insertion site. The bla_CTX-M-15 insertion site in isolate 623214 was described previously, whereas that of isolate 542093 was unique to this study. Analysis of Nanopore sequencing data enables us to characterise the genomic architecture of mobile genetic elements encoding AMR determinants. These data may contribute to a better understanding of persistence and onward transmission of AMR determinants in multidrug-resistant (MDR) E. coli causing gastrointestinal and extra-intestinal infections.

Precise genotyping of circular mobile elements from metagenomic data uncovers human-associated plasmids with recent common ancestors

Mobile genetic elements with circular genomes play a key role in the evolution of microbial communities. Their circular genomes correspond to circular walks in metagenome graphs, and yet, assemblies derived from natural microbial communities produce graphs riddled with spurious cycles, complicating the accurate reconstruction of circular genomes. We present DomCycle, an algorithm that reconstructs likely circular genomes based on the identification of so-called 'dominant' graph cycles. In the implementation we leverage paired reads to bridge assembly gaps and scrutinize cycles through a nucleotide-level analysis, making the approach robust to misassembly artifacts. We validated the approach using simulated and real sequencing data. Application of DomCycle to 32 publicly available DNA shotgun sequence data sets from diverse natural environments led to the reconstruction of hundreds of circular mobile genomes. Clustering revealed 20 highly prevalent and cryptic plasmids that have clonal population structures with recent common ancestors. This method facilitates the study of microbial communities that evolve through horizontal gene transfer.