plasmidSPAdes: assembling plasmids from whole genome sequencing data

Genetic Configuration of Genomic Resistance Islands in Acinetobacter baumannii Clinical Isolates From Egypt

In Acinetobacter baumannii (A. baumannii), a wide repertoire of resistance genes is often carried within genomic resistance islands (RIs), particularly in high-risk global clones (GCs). As the first in Egypt, the current study aimed at exploring the diversity and genetic configuration of RIs in the clinical isolates of A. baumannii. For this purpose, draft genomes of 18 isolates were generated by Illumina sequencing. Disk diffusion susceptibility profiling revealed multidrug resistance (MDR) and extensive drug resistance (XDR) phenotypes in 27.7 and 72.2%, respectively. The highest susceptibility was noted for tigecycline (100.0%) followed by colistin (94.4%), for which an MIC₅₀ of 0.25 μg/ml was recorded by the broth microdilution assay. Sequence typing (ST) showed that the majority of the isolates belonged to high-risk global clones (GC1, GC2, and GC9). A novel Oxford sequence type (ST2329) that also formed a novel clonal complex was submitted to the PubMLST database. A novel bla_ADC variant (bla_ADC−258) was also identified in strain M18 (ST85^Pas/1089^Oxf). In addition to a wide array of resistance determinants, whole-genome sequencing (WGS) disclosed at least nine configurations of genomic RIs distributed over 16/18 isolates. GC2 isolates accumulated the largest number of RIs (three RIs/isolate) followed by those that belong to GC1 (two RIs/isolate). In addition to Tn6022 (44.4%), the comM gene was interrupted by AbaR4 (5.5%) and three variants of A. baumanniigenomic resistance island 1(AbGRI)-type RIs (44.4%), including AbaR4b (16.6%) and two novel configurations of AbGRI1-like RIs (22.2%). Three of which (AbaR4, AbaR4b, and AbGRI1-like-2) carried bla_OXA−23 within Tn2006. With less abundance (38.8%), IS26-bound RIs were detected exclusively in GC2 isolates. These included a short version of AbGRI2 (AbGRI2-15) carrying the genes bla_TEM−1 and aphA1 and two variants of AbGRI3 RIs carrying up to seven resistance genes [mphE-msrE-armA-sul1-aadA1-catB8-aacA4]. Confined to GC1 (22.2%), sulfonamide resistance was acquired by an ISAba1 bracketed GIsul2 RI. An additional RI (RI-PER-7) was also identified on a plasmid carried by strain M03. Among others, RI-PER-7 carried the resistance genes armA and bla_PER−7. Here, we provided a closer view of the diversity and genetic organization of RIs carried by a previously unexplored population of A. baumannii.

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.

Whole-Genome Sequence of Multidrug-Resistant Methicillin-Resistant Staphylococcus epidermidis Carrying Biofilm-Associated Genes and a Unique Composite of SCCmec

Staphylococcus epidermidis is part of the normal human flora that has recently become an important opportunistic pathogen causing nosocomial infections and tends to be multidrug-resistant. In this investigation, we aimed to study the genomic characteristics of methicillin-resistant S. epidermidis isolated from clinical specimens. Three isolates were identified using biochemical tests and evaluated for drug susceptibility. Genomic DNA sequences were obtained using Illumina, and were processed for analysis using various bioinformatics tools. The isolates showed multidrug resistance to most of the antibiotics tested in this study, and were identified with three types (III(3A), IV(2B&5), and VI(4B)) of the mobile genetic element SCCmec that carries the methicillin resistance gene (mecA) and its regulators (mecI and mecR1). A total of 11 antimicrobial resistance genes (ARGs) was identified as chromosomally mediated or in plasmids; these genes encode for proteins causing decreased susceptibility to methicillin (mecA), penicillin (blaZ), fusidic acid (fusB), fosfomycin (fosB), tetracycline (tet(K)), aminoglycosides (aadD, aac(6′)-aph(2′’)), fluoroquinolone (MFS antibiotic efflux pump), trimethoprim (dfrG), macrolide (msr(A)), and chlorhexidine (qacA)). Additionally, the 9SE strain belongs to the globally disseminated ST2, and harbors biofilm-formation genes (icaA, icaB, icaC, icaD, and IS256) with phenotypic biofilm production capability. It also harbors the fusidic acid resistance gene (fusB), which could increase the risk of device-associated healthcare infections, and 9SE has been identified as having a unique extra SCC gene (ccrB4); this new composite element of the ccr type needs more focus to better understand its role in the drug resistance mechanism.

Comparative genomics of the Western Hemisphere soft tick-borne relapsing fever borreliae highlights extensive plasmid diversity

BACKGROUND

Tick-borne relapsing fever (TBRF) is a globally prevalent, yet under-studied vector-borne disease transmitted by soft and hard bodied ticks. While soft TBRF (sTBRF) spirochetes have been described for over a century, our understanding of the molecular mechanisms facilitating vector and host adaptation is poorly understood. This is due to the complexity of their small (~ 1.5 Mb) but fragmented genomes that typically consist of a linear chromosome and both linear and circular plasmids. A majority of sTBRF spirochete genomes' plasmid sequences are either missing or are deposited as unassembled sequences. Consequently, our goal was to generate complete, plasmid-resolved genomes for a comparative analysis of sTBRF species of the Western Hemisphere.

RESULTS

Utilizing a Borrelia specific pipeline, genomes of sTBRF spirochetes from the Western Hemisphere were sequenced and assembled using a combination of short- and long-read sequencing technologies. Included in the analysis were the two recently isolated species from Central and South America, Borrelia puertoricensis n. sp. and Borrelia venezuelensis, respectively. Plasmid analyses identified diverse sequences that clustered plasmids into 30 families; however, only three families were conserved and syntenic across all species. We also compared two species, B. venezuelensis and Borrelia turicatae, which were isolated ~ 6,800 km apart and from different tick vector species but were previously reported to be genetically similar.

CONCLUSIONS

To truly understand the biological differences observed between species of TBRF spirochetes, complete chromosome and plasmid sequences are needed. This comparative genomic analysis highlights high chromosomal synteny across the species yet diverse plasmid composition. This was particularly true for B. turicatae and B. venezuelensis, which had high average nucleotide identity yet extensive plasmid diversity. These findings are foundational for future endeavors to evaluate the role of plasmids in vector and host adaptation.

Whole-genome sequence analysis of Clostridium chauvoei isolated from clinical case of black quarter (BQ) from India

Black quarter (BQ) is an infectious disease affecting cattle and small ruminants worldwide caused by Gram-positive anaerobic bacterium Clostridium chauvoei. In this study, a draft genome sequence of C. chauvoei NIVEDIBQ1 strain isolated from clinical case of black quarter was analyzed. Sequence analysis indicated that genome had 2653 predicted coding DNA sequences, harbored numerous genes, mobile genetic elements for pathogenesis, and virulence factors. Computational analysis revealed that strain contained 30 virulence-associated genes. An intact genomic region highly similar to the Clostridium phage was present in the genome. Presence of CRISPR systems and the transposon components likely contribute to the genome plasticity. Strain encode diverse spectrum of degradative carbohydrate-active enzymes (CAZymes). Comparative SNP analysis revealed that the genomes of the C. chauvoei strains analyzed were highly conserved. Phylogenetic analysis of strains and available genome (n = 21) based on whole-genome multi-locus sequence typing (wgMLST) and core orthologous genes showed the clustering of strains into two different clusters suggesting geographical links.

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.

Profiling the plasmid conjugation potential of urinary Escherichia coli

Escherichia coli is often associated with urinary tract infection (UTI). Antibiotic resistance in E. coli is an ongoing challenge in managing UTI. Extrachromosomal elements - plasmids - are vectors for clinically relevant traits, such as antibiotic resistance, with conjugation being one of the main methods for horizontal propagation of plasmids in bacterial populations. Targeting of conjugation components has been proposed as a strategy to curb the spread of plasmid-borne antibiotic resistance. Understanding the types of conjugative systems present in urinary E. coli isolates is fundamental to assessing the viability of this strategy. In this study, we profile two well-studied conjugation systems (F-type and P-type) in the draft genomes of 65 urinary isolates of E. coli obtained from the bladder urine of adult women with and without UTI-like symptoms. Most of these isolates contained plasmids and we found that conjugation genes were abundant/ubiquitous, diverse and often associated with IncF plasmids. To validate conjugation of these urinary plasmids, the plasmids from two urinary isolates, UMB1223 (predicted to have F-type genes) and UMB1284 (predicted to have P-type genes), were transferred by conjugation into the K-12 E. coli strain MG1655. Overall, the findings of this study support the notion that care should be taken in targeting any individual component of a urinary E. coli isolate's conjugation system, given the inherent mechanistic redundancy, gene diversity and different types of conjugation systems in this population.

Genomic Analysis of Multidrug-Resistant Hypervirulent (Hypermucoviscous) Klebsiella pneumoniae Strain Lacking the Hypermucoviscous Regulators (rmpA/rmpA2)

Hypervirulent K. pneumoniae (hvKP) strains possess distinct characteristics such as hypermucoviscosity, unique serotypes, and virulence factors associated with high pathogenicity. To better understand the genomic characteristics and virulence profile of the isolated hvKP strain, genomic data were compared to the genomes of the hypervirulent and typical K. pneumoniae strains. The K. pneumoniae strain was isolated from a patient with a recurrent urinary tract infection, and then the string test was used for the detection of the hypermucoviscosity phenotype. Whole-genome sequencing was conducted using Illumina, and bioinformatics analysis was performed for the prediction of the isolate resistome, virulome, and phylogenetic analysis. The isolate was identified as hypermucoviscous, type 2 (K2) capsular polysaccharide, ST14, and multidrug-resistant (MDR), showing resistance to ciprofloxacin, ceftazidime, cefotaxime, trimethoprim-sulfamethoxazole, cephalexin, and nitrofurantoin. The isolate possessed four antimicrobial resistance plasmids (pKPN3-307_type B, pECW602, pMDR, and p3K157) that carried antimicrobial resistance genes (ARGs) (bla_OXA-1,bla_CTX-M-15, sul2, APH(3″)-Ib, APH(6)-Id, and AAC(6′)-Ib-cr6). Moreover, two chromosomally mediated ARGs (fosA6 and SHV-28) were identified. Virulome prediction revealed the presence of 19 fimbrial proteins, one aerobactin (iutA) and two salmochelin (iroE and iroN). Four secretion systems (T6SS-I (13), T6SS-II (9), T6SS-III (12), and Sci-I T6SS (1)) were identified. Interestingly, the isolate lacked the known hypermucoviscous regulators (rmpA/rmpA2) but showed the presence of other RcsAB capsule regulators (rcsA and rcsB). This study documented the presence of a rare MDR hvKP with hypermucoviscous regulators and lacking the common capsule regulators, which needs more focus to highlight their epidemiological role.

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.

Genomics of human and chicken Salmonella isolates in Senegal: Broilers as a source of antimicrobial resistance and potentially invasive nontyphoidal salmonellosis infections

Salmonella enterica is the most common foodborne pathogen worldwide. It causes two types of diseases, a self-limiting gastroenteritis and an invasive, more threatening, infection. Salmonella gastroenteritis is caused by several serotypes and is common worldwide. In contrast, invasive salmonellosis is rare in high-income countries (HIC) while frequent in low- and middle-income countries (LMIC), especially in sub-Saharan Africa (sSA). Invasive Nontyphoidal Salmonella (iNTS), corresponding to serotypes other than Typhi and Paratyphi, have emerged in sSA and pose a significant risk to public health. We conducted a whole-genome sequence (WGS) analysis of 72 strains of Salmonella isolated from diarrheic human patients and chicken meat sold in multipurpose markets in Dakar, Senegal. Antimicrobial susceptibility testing combined with WGS data analysis revealed frequent resistance to fluoroquinolones and the sulfamethoxazole-trimethoprim combination that are among the most used treatments for invasive Salmonella. In contrast, resistance to the historical first-line drugs chloramphenicol and ampicillin, and to cephalosporins was rare. Antimicrobial resistance (AMR) was lower in clinical isolates compared to chicken strains pointing to the concern posed by the excessive use of antimicrobials in farming. Phylogenetic analysis suggested possible transmission of the emerging multidrug resistant (MDR) Kentucky ST198 and serotype Schwarzengrund from chicken to human. These results stress the need for active surveillance of Salmonella and AMR in order to address invasive salmonellosis caused by nontyphoidal Salmonella strains and other important bacterial diseases in sSA.

Comparative genomics of 26 complete circular genomes of 18 different serotypes of Actinobacillus pleuropneumoniae

Actinobacillus pleuropneumoniae is a Gram-negative, rod-shaped bacterium of the family Pasteurellaceae causing pig pleuropneumonia associated with great economic losses worldwide. Nineteen serotypes with distinctive lipopolysaccharide (LPS) and capsular (CPS) compositions have been described so far, yet complete circular genomes are publicly available only for the reference strains of serotypes 1, 4 and 5b, and for field strains of serotypes 1, 3, 7 and 8. We aimed to complete this picture by sequencing the reference strains of 17 different serotypes with the MinION sequencer (Oxford Nanopore Technologies, ONT) and on an Illumina HiSeq (Illumina) platform. We also included two field isolates of serotypes 2 and 3 that were PacBio- and MinION-sequenced, respectively. Genome assemblies were performed following two different strategies, i.e. PacBio- or ONT-only de novo assemblies polished with Illumina reads or a hybrid assembly by directly combining ONT and Illumina reads. Both methods proved successful in obtaining accurate circular genomes with comparable qualities. blast-based genome comparisons and core-genome phylogeny based on core genes, SNP typing and multi-locus sequence typing (cgMLST) of the 26 circular genomes indicated well-conserved genomes across the 18 different serotypes, differing mainly in phage insertions, and CPS, LPS and RTX-toxin clusters, which, consistently, encode serotype-specific antigens. We also identified small antibiotic resistance plasmids, and complete subtype I-F and subtype II-C CRISPR-Cas systems. Of note, highly similar clusters encoding all those serotype-specific traits were also found in other pathogenic and commensal Actinobacillus species. Taken together with the presence of transposable elements surrounding these loci, we speculate a dynamic intra- and interspecies exchange of such virulence-related factors by horizontal gene transfer. In conclusion, our comprehensive genomics analysis provides useful information for diagnostic test and vaccine development, but also for whole-genome-based epidemiological studies, as well as for the surveillance of the evolution of antibiotic resistance and virulence genes in A. pleuropneumoniae.

Species-Level Analysis of the Human Gut Microbiome Shows Antibiotic Resistance Genes Associated With Colorectal Cancer

Colorectal cancer (CRC) is the second leading cause of cancer deaths and continuously increases new cancer cases globally. Accumulating evidence links risks of CRC to antibiotic use. Long-term use and abuse of antibiotics increase the resistance of the gut microbiota; however, whether CRC is associated with antibiotic resistance in gut microbiota is still unclear. In this study, we performed a de novo assembly to metagenomic sequences in 382 CRC patients and 387 healthy controls to obtain representative species-level genome bins (rSGBs) and plasmids and analyzed the abundance variation of species and antibiotic resistance genes (ARGs). Twenty-five species and 65 ARGs were significantly enriched in the CRC patients, and among these ARGs, 12 were multidrug-resistant genes (MRGs), which mainly included acrB, TolC, marA, H-NS, Escherichia coli acrR mutation, and AcrS. These MRGs could confer resistance to fluoroquinolones, tetracyclines, cephalosporins, and rifamycin antibiotics by antibiotic efflux and inactivation. A classification model was built using the abundance of species and ARGs and achieved areas under the curve of 0.831 and 0.715, respectively. Our investigation has identified the antibiotic resistance types of ARGs and suggested that E. coli is the primary antibiotic resistance reservoir of ARGs in CRC patients, providing valuable evidence for selecting appropriate antibiotics in the CRC treatment.

Ligand Selection for Affinity Chromatography Using Phage Display

Phage display coupled with in vitro affinity selection to mimic evolutionary principles has propelled the discovery of specific binding peptides and proteins for diverse applications, including affinity chromatography. By tailoring screening conditions, ligands with desired predefined properties, such as pH- or ion strength-responsive binding, can be identified from phage-displayed combinatorial peptide libraries. Initial hit peptides can be further optimized through directed evolution by focused mutagenesis and rescreening. Quantitative analysis of eluted binders with next-generation sequencing (NGS) assists in reducing enrichment bias and simplifies picking the most promising ligand candidate(s) through enrichment ranking. We describe, in detail, procedures of ligand selection for affinity chromatography using peptide phage display library screening, focused mutagenesis, and NGS. Furthermore, we outline the subsequent workflow for ligand characterization and affinity column construction.

Reconstruction of Archaeal Genomes from Short-Read Metagenomes

As the majority of biological diversity remains unexplored and uncultured, investigating it requires culture-independent approaches. Archaea in particular suffer from a multitude of issues that make their culturing problematic, from them being frequently members of the rare biosphere, to low growth rates, to them thriving under very specific and often extreme environmental and community conditions that are difficult to replicate. OMICs techniques are state of the art approaches that allow direct high-throughput investigations of environmental samples at all levels from nucleic acids to proteins, lipids, and secondary metabolites. Metagenomics, as the foundation for other OMICs techniques, facilitates the identification and functional characterization of the microbial community members and can be combined with other methods to provide insights into the microbial activities, both on the RNA and protein levels. In this chapter, we provide a step-by-step workflow for the recovery of archaeal genomes from metagenomes, starting from raw short-read sequences. This workflow can be applied to recover bacterial genomes as well.

Whole-Genome Sequencing of Gram-Negative Bacteria Isolated From Bovine Mastitis and Raw Milk: The First Emergence of Colistin mcr-10 and Fosfomycin fosA5 Resistance Genes in Klebsiella pneumoniae in Middle East

Antimicrobial resistance is a major concern in the dairy industry. This study investigated the prevalence, antimicrobial resistance phenotypes, and genome sequencing of Gram-negative bacteria isolated from clinical (n = 350) and subclinical (n = 95) bovine mastitis, and raw unpasteurized milk (n = 125). Klebsiella pneumoniae, Aeromonas hydrophila, Enterobacter cloacae (100% each), Escherichia coli (87.78%), and Proteus mirabilis (69.7%) were the most prevalent multidrug-resistant (MDR) species. Extensive drug-resistance (XDR) phenotype was found in P. mirabilis (30.30%) and E. coli (3.33%) isolates. Ten isolates (four E. coli, three Klebsiella species and three P. mirabilis) that displayed the highest multiple antibiotic resistance (MAR) indices (0.54–0.83), were exposed to whole-genome sequencing (WGS). Two multilocus sequence types (MLST): ST2165 and ST7624 were identified among the sequenced E. coli isolates. Three E. coli isolates (two from clinical mastitis and one from raw milk) belonging to ST2165 showed similar profile of plasmid replicon types: IncFIA, IncFIB, IncFII, and IncQ1 with an exception to an isolate that contained IncR, whereas E. coli ST7624 showed a different plasmid profile including IncHI2, IncHI2A, IncI1α, and IncFII replicon types. ResFinder findings revealed the presence of plasmid-mediated colistin mcr-10 and fosfomycin fosA5 resistance genes in a K. pneumoniae (K1) isolate from bovine milk. Sequence analysis of the reconstructed mcr-10 plasmid from WGS of K1 isolate, showed that mcr-10 gene was bracketed by xerC and insertion sequence IS26 on an IncFIB plasmid. Phylogenetic analysis revealed that K1 isolate existed in a clade including mcr-10-harboring isolates from human and environment with different STs and countries [United Kingdom (ST788), Australia (ST323), Malawi (ST2144), Myanmar (ST705), and Laos (ST2355)]. This study reports the first emergence of K. pneumoniae co-harboring mcr-10 and fosA5 genes from bovine milk in the Middle East, which constitutes a public health threat and heralds the penetration of the last-resort antibiotics. Hence, prudent use of antibiotics in both humans and animals and antimicrobial surveillance plans are urgently required.

Characterizing Plasmids in Bacteria Species Relevant to Urinary Health

The urinary tract has a microbial community (the urinary microbiota or urobiota) that has been associated with human health. Whole genome sequencing of bacteria is a powerful tool, allowing investigation of the genomic content of the urobiota, also called the urinary microbiome (urobiome). Bacterial plasmids are a significant component of the urobiome yet are understudied. Because plasmids can be vectors and reservoirs for clinically relevant traits, they are important for urobiota dynamics and thus may have relevance to urinary health. In this project, we sought plasmids in 11 clinically relevant urinary species: Aerococcus urinae, Corynebacterium amycolatum, Enterococcus faecalis, Escherichia coli, Gardnerella vaginalis, Klebsiella pneumoniae, Lactobacillus gasseri, Lactobacillus jensenii, Staphylococcus epidermidis, Streptococcus anginosus, and Streptococcus mitis. We found evidence of plasmids in E. faecalis, E. coli, K. pneumoniae, S. epidermidis, and S. anginosus but insufficient evidence in other species sequenced thus far. Some identified plasmidic assemblies were predicted to have putative virulence and/or antibiotic resistance genes, although the majority of their annotated coding regions were of unknown predicted function. In this study, we report on plasmids from urinary species as a first step to understanding the role of plasmids in the bacterial urobiota. IMPORTANCE The microbial community of the urinary tract (urobiota) has been associated with human health. Whole genome sequencing of bacteria permits examination of urobiota genomes, including plasmids. Because plasmids are vectors and reservoirs for clinically relevant traits, they are important for urobiota dynamics and thus may have relevance to urinary health. Currently, urobiota plasmids are understudied. Here, we sought plasmids in 11 clinically relevant urinary species. We found evidence of plasmids in E. faecalis, E. coli, K. pneumoniae, S. epidermidis, and S. anginosus but insufficient evidence in the other 6 species. We identified putative virulence and/or antibiotic resistance genes in some of the plasmidic assemblies, but most of their annotated coding regions were of unknown function. This is a first step to understanding the role of plasmids in the bacterial urobiota.

Genomic Characterization of Extensively Drug-Resistant NDM-Producing Acinetobacter baumannii Clinical Isolates With the Emergence of Novel bla ADC-257

Acinetobacter baumannii has become a major challenge to clinicians worldwide due to its high epidemic potential and acquisition of antimicrobial resistance. This work aimed at investigating antimicrobial resistance determinants and their context in four extensively drug-resistant (XDR) NDM-producing A. baumannii clinical isolates collected between July and October 2020 from Kasr Al-Ainy Hospital, Cairo, Egypt. A total of 20 A. baumannii were collected and screened for acquired carbapenemases (bla_NDM, bla_VIM and bla_IMP) using PCR. Four NDM producer A. baumannii isolates were identified and selected for whole-genome sequencing, in silico multilocus sequence typing, and resistome analysis. Antimicrobial susceptibility profiles were determined using disk diffusion and broth microdilution tests. All bla_NDM-positive A. baumannii isolates were XDR. Three isolates belonged to high-risk international clones (IC), namely, IC2 corresponding to ST570^Pas/1701^Oxf (M20) and IC9 corresponding to ST85^Pas/ST1089^Oxf (M02 and M11). For the first time, we report bla_NDM-1 gene on the chromosome of an A. baumannii strain that belongs to sequence type ST164^Pas/ST1418^Oxf. Together with AphA6, bla_NDM-1 was bracketed by two copies of ISAba14 in ST85^Pas isolates possibly facilitating co-transfer of amikacin and carbapenem resistance. A novel bla_ADC allele (bla_ADC-257) with an upstream ISAba1 element was identified in M19 (ST/CC164^Pas and ST1418^Oxf/CC234^Oxf). bla_ADC genes harbored by M02 and M11 were uniquely interrupted by IS1008. Tn2006-associated bla_OXA-23 was carried by M20. bla_OXA-94 genes were preceded by ISAba1 element in M02 and M11. AbGRI3 was carried by M20 hosting the resistance genes aph(3`)-Ia, aac(6`)-Ib`, catB8, ant(3``)-Ia, sul1, armA, msr(E), and mph(E). Nonsynonymous mutations were identified in the quinolone resistance determining regions (gyrA and parC) of all isolates. Resistance to colistin in M19 was accompanied by missense mutations in lpxACD and pmrABC genes. The current study provided an insight into the genomic background of XDR phenotype in A. baumannii recovered from patients in Egypt. WGS revealed strong association between resistance genes and diverse mobile genetic elements with novel insertion sites and genetic organizations.

Long-Read Sequencing Reveals Evolution and Acquisition of Antimicrobial Resistance and Virulence Genes in Salmonella enterica

Salmonella enterica is a significant and phylogenetically diverse zoonotic pathogen. To understand its genomic heterogeneity and antimicrobial resistance, we performed long-read sequencing on Salmonella isolated from retail meats and food animals. A collection of 134 multidrug-resistant isolates belonging to 33 serotypes were subjected to PacBio sequencing. One major locus of diversity among these isolates was the presence and orientation of Salmonella pathogenic islands (SPI), which varied across different serotypes but were largely conserved within individual serotypes. We also identified insertion of an IncQ resistance plasmid into the chromosome of fourteen strains of serotype I 4,[5],12:i:- and the Salmonella genomic island 1 (SGI-1) in five serotypes. The presence of various SPIs, SGI-1 and integrated plasmids contributed significantly to the genomic variability and resulted in chromosomal resistance in 55.2% (74/134) of the study isolates. A total of 93.3% (125/134) of isolates carried at least one plasmid, with isolates carrying up to seven plasmids. We closed 233 plasmid sequences of thirteen replicon types, along with twelve hybrid plasmids. Some associations between Salmonella isolate source, serotype, and plasmid type were seen. For instance, IncX plasmids were more common in serotype Kentucky from retail chicken. Plasmids IncC and IncHI had on average more than five antimicrobial resistance genes, whereas in IncX, it was less than one per plasmid. Overall, 60% of multidrug resistance (MDR) strains that carried >3 AMR genes also carried >3 heavy metal resistance genes, raising the possibility of co-selection of antimicrobial resistance in the presence of heavy metals. We also found nine isolates representing four serotypes that carried virulence plasmids with the spv operon. Together, these data demonstrate the power of long-read sequencing to reveal genomic arrangements and integrated plasmids with a high level of resolution for tracking and comparing resistant strains from different sources. Additionally, the findings from this study will help expand the reference set of closed Salmonella genomes that can be used to improve genome assembly from short-read data commonly used in One Health antimicrobial resistance surveillance.

Genetic and Phenotypic Diversity of Morganella morganii Isolated From Cheese

The bacterium Morganella morganii can produce the biogenic amines (BA) cadaverine, putrescine, and histamine in vitro and is responsible for high histamine concentrations in fish products. These BA can have toxic effects upon ingestion and are undesired in food. The purpose of this study was to characterize the phenotype and genotype of 11 M. morganii isolated from cheese in regard to the BA formation. In addition, we investigated the phylogeny, trehalose fermentation ability, and antibiotic resistance of the cheese isolates. To do so, we sequenced their genomes using both long and short read technologies. Due to the presence of the trehalose operon and the ability to ferment trehalose, the cheese isolates can be assigned to the subsp. sibonii. Comparative genomics with public available M. morganii genomes shows that the genomes of the cheese isolates cluster together with other subsp. sibonii genomes. All genomes between subsp. morganii and subsp. sibonii are separated by an average nucleotide identity (ANI) of less than 95.0%. Therefore, the subspecies could represent two distinct species. Nine of the strains decarboxylated lysine yielding cadaverine in vitro. This metabolic activity is linked to a previously unknown gene cluster comprising genes encoding a lysine-tRNA ligase (lysS), an HTH-transcriptional regulator (argP), a cadaverine-lysine antiporter (cadB), and a lysine decarboxylase (cadA). The formation of putrescine is linked to the speF gene encoding an ornithine decarboxylase. The gene is disrupted in five strains by an insertion sequence, and these strains only exhibit a weak putrescine production. Antimicrobial susceptibility profiling revealed that all cheese strains are resistant to tetracycline, chloramphenicol, tigecycline, colistin, and ampicillin. These phenotypes, except for colistin which is intrinsic, could be linked to antimicrobial resistance genes located on the chromosome.

Draft Genome Sequences of Cefepime-Resistant and -Susceptible Escherichia coli Strains and Imipenem-Resistant and -Susceptible Pseudomonas aeruginosa Strains

Draft genome sequences of Escherichia coli and Pseudomonas aeruginosa strains collected from clinical infections were used to determine the prevalence of newly emerging antibiotic resistance genes in Maine. Comparisons between cefepime-resistant and -susceptible E. coli strains and imipenem-resistant and -susceptible P. aeruginosa strains are being conducted.

RFPlasmid: predicting plasmid sequences from short-read assembly data using machine learning

Antimicrobial-resistance (AMR) genes in bacteria are often carried on plasmids and these plasmids can transfer AMR genes between bacteria. For molecular epidemiology purposes and risk assessment, it is important to know whether the genes are located on highly transferable plasmids or in the more stable chromosomes. However, draft whole-genome sequences are fragmented, making it difficult to discriminate plasmid and chromosomal contigs. Current methods that predict plasmid sequences from draft genome sequences rely on single features, like k-mer composition, circularity of the DNA molecule, copy number or sequence identity to plasmid replication genes, all of which have their drawbacks, especially when faced with large single-copy plasmids, which often carry resistance genes. With our newly developed prediction tool RFPlasmid, we use a combination of multiple features, including k-mer composition and databases with plasmid and chromosomal marker proteins, to predict whether the likely source of a contig is plasmid or chromosomal. The tool RFPlasmid supports models for 17 different bacterial taxa, including Campylobacter, Escherichia coli and Salmonella, and has a taxon agnostic model for metagenomic assemblies or unsupported organisms. RFPlasmid is available both as a standalone tool and via a web interface.

Evolution of ColE1-like plasmids across γ-Proteobacteria: From bacteriocin production to antimicrobial resistance

Antimicrobial resistance is one of the major threats to Public Health worldwide. Understanding the transfer and maintenance of antimicrobial resistance genes mediated by mobile genetic elements is thus urgent. In this work, we focus on the ColE1-like plasmid family, whose distinctive replication and multicopy nature has given rise to key discoveries and tools in molecular biology. Despite being massively used, the hosts, functions, and evolutionary history of these plasmids remain poorly known. Here, we built specific Hidden Markov Model (HMM) profiles to search ColE1 replicons within genomes. We identified 1,035 ColE1 plasmids in five Orders of γ-Proteobacteria, several of which are described here for the first time. The phylogenetic analysis of these replicons and their characteristic MOB_P5/HEN relaxases suggest that ColE1 plasmids have diverged apart, with little transfer across orders, but frequent transfer across families. Additionally, ColE1 plasmids show a functional shift over the last decades, losing their characteristic bacteriocin production while gaining several antimicrobial resistance genes, mainly enzymatic determinants and including several extended-spectrum betalactamases and carbapenemases. Furthermore, ColE1 plasmids facilitate the intragenomic mobilization of these determinants, as various replicons were identified co-integrated with large non-ColE1 plasmids, mostly via transposases. These results illustrate how families of plasmids evolve and adapt their gene repertoires to bacterial adaptive requirements.

The extraordinary adaptability of bacteria and the massive prevalence of mobile genetic elements within populations has turned antimicrobial resistance into a growing threat to Public Health. Among all the mobile genetic elements, plasmids have been the focus of attention as these extrachromosomal molecules of DNA are able to mobilize several antimicrobial resistance genes at once through conjugation. However, although small mobilizable and non-conjugative replicons have been traditionally overlooked when analyzing plasmid-mediated antimicrobial resistance, they have recently been described as important carriers of AMR genes. In this work, we have analyzed the ColE1-like plasmid family, whose study has been neglected even if they are one of the main groups of small plasmids in natural populations of Proteobacteria. We observed that these plasmids have evolved for a long time within γ-Proteobacteria acquiring different genetic features in specific hosts, being major players in the spread of antimicrobial resistance determinants.

Comparative genomics of Photobacterium species from terrestrial and marine habitats

Photobacterium (P.) is a genus widely studied in regards to its association with and ubiquitous presence in marine environments. However, certain species (P. phosphoreum, P. carnosum, P. iliopiscarium) have been recently described to colonize and spoil raw meats without a marine link. We have studied 27 strains from meat as well as 26 strains from marine environments in order to probe for intraspecies marine/terrestrial subpopulations and identify distinct genomic features acquired by environmental adaptation. We have conducted phylogenetic analysis (MLSA, ANI, fur, codon usage), search of plasmids (plasmidSPADES), phages (PHASTER), CRISPR-cas operons (CRISPR-finder) and secondary metabolites gene clusters (antiSMASH, BAGEL), in addition to a targeted gene search for specific pathways (e.g. TCA cycle, pentose phosphate, respiratory chain) and elements relevant for growth, adaptation and competition (substrate utilization, motility, bioluminescence, sodium and iron transport). P. carnosum appears as a conserved single clade, with one isolate from MAP fish clustering apart that doesn't, however, show distinct features that could indicate different adaptation. The species harbors genes for a wide carbon source utilization (glycogen/starch, maltose, pullulan, fucose) for colonization of diverse niches in its genome. P. phosphoreum is represented by two different clades on the phylogenetic analyses not correlating to their origin or distribution of other features analyzed that can be divided into two novel subspecies based on genome-wide values. A more diverse antimicrobial activity (sactipeptides, microcins), production of secondary metabolites (siderophores and arylpolyenes), stress response and adaptation (bioluminescence, sodium transporters, catalase, high affinity for oxygen cytochrome cbb3 oxidase, DMSO reductase and proton translocating NADH dehydrogenase) is predicted compared to the other species. P. iliopiscarium was divided into two clades based on source of isolation correlating with phylogeny and distribution of several traits. The species shows traits common to the other two species, similar carbon utilization/transport gene conservation as P. carnosum for the meat-isolated strains, and predicted utilization of marine-common DMSO and flagellar cluster for the sea-isolated strains. Results additionally suggest that photobacteria are highly prone to horizontal acquisition/loss of genetic material and genetic transduction, and that it might be a strategy for increasing the frequency of strain- or species-specific features that offers a growth/competition advantage.

Genomic diversity of antimicrobial resistance in non-typhoidal Salmonella in Victoria, Australia

Non-typhoidal Salmonella (NTS) is the second most common cause of foodborne bacterial gastroenteritis in Australia with antimicrobial resistance (AMR) increasing in recent years. Whole-genome sequencing (WGS) provides opportunities for in silico detection of AMR determinants. The objectives of this study were two-fold: (1) establish the utility of WGS analyses for inferring phenotypic resistance in NTS, and (2) explore clinically relevant genotypic AMR profiles to third generation cephalosporins (3GC) in NTS lineages. The concordance of 2490 NTS isolates with matched WGS and phenotypic susceptibility data against 13 clinically relevant antimicrobials was explored. In silico serovar prediction and typing was performed on assembled reads and interrogated for known AMR determinants. The surrounding genomic context, plasmid determinants and co-occurring AMR patterns were further investigated for multidrug resistant serovars harbouring bla CMY-2, bla CTX-M-55 or bla CTX-M-65. Our data demonstrated a high correlation between WGS and phenotypic susceptibility testing. Phenotypic-genotypic concordance was observed between 2440/2490 (98.0 %) isolates, with overall sensitivity and specificity rates >98 % and positive and negative predictive values >97 %. The most common AMR determinants were bla TEM-1, sul2 , tet (A), strA-strB and floR . Phenotypic resistance to cefotaxime and azithromycin was low and observed in 6.2 % (151/2486) and 0.9 % (16/1834) of the isolates, respectively. Several multi-drug resistant NTS lineages were resistant to 3GC due to different genetic mechanisms including bla CMY-2, bla CTX-M-55 or bla CTX-M-65. This study shows WGS can enhance existing AMR surveillance in NTS datasets routinely produced in public health laboratories to identify emerging AMR in NTS. These approaches will be critical for developing capacity to detect emerging public health threats such as resistance to 3GC.

pLS20 is the archetype of a new family of conjugative plasmids harboured by Bacillus species

Conjugation plays important roles in genome plasticity, adaptation and evolution but is also the major horizontal gene-transfer route responsible for spreading toxin, virulence and antibiotic resistance genes. A better understanding of the conjugation process is required for developing drugs and strategies to impede the conjugation-mediated spread of these genes. So far, only a limited number of conjugative elements have been studied. For most of them, it is not known whether they represent a group of conjugative elements, nor about their distribution patterns. Here we show that pLS20 from the Gram-positive bacterium Bacillus subtilis is the prototype conjugative plasmid of a family of at least 35 members that can be divided into four clades, and which are harboured by different Bacillus species found in different global locations and environmental niches. Analyses of their phylogenetic relationship and their conjugation operons have expanded our understanding of a family of conjugative plasmids of Gram-positive origin.

Comparative pangenomic analyses and biotechnological potential of cocoa-related Acetobacter senegalensis strains

Acetobacter senegalensis belongs to the group of acetic acid bacteria (AAB) that present potential biotechnological applications, for production of D-gluconate, cellulose and acetic acid. AAB can overcome heat and acid stresses by using strategies involving the overexpression of heat-shock proteins and enzymes from the complex pyrroquinoline-ADH, besides alcohol dehydrogenases (ADH). Nonetheless, the isolation of A. senegalensis and other AAB from food may be challenging due to presence of viable but non-culturable (VBNC) cells and due to uncertainties about nutritional requirements. To contribute for a better understanding of the ecology of AAB, this paper reports on the pangenome analysis of five strains of A. senegalensis recently isolated from a Brazilian spontaneous cocoa fermentation. The results showed biosynthetic clusters exclusively found in some cocoa-related AAB, such as those related to terpene pathways, which are important for flavour development. Genes related to oxidative stress were conserved in all the genomes, with multiple clusters. Moreover, there were genes coding for ADH and putative ABC transporters distributed in core, shell and cloud genomes, while chaperonin-encoding genes were present only in the core and soft-core genomes. Regarding quorum sensing, a response regulator gene was in the shell genome, and the gene encoding for acyl-homoserine lactone efflux protein was in the soft-core genome. There were quorum quenching-related genes, mainly encoding for lactonases, but also for acylases. Moreover, A. senegalensis did not have determinants of virulence or antibiotic resistance, which are good traits for strains intended to be applied in food fermentation.

Extensive Genome Exploration of Clostridium botulinum Group III Field Strains

In animals, botulism is commonly sustained by botulinum neurotoxin C, D or their mosaic variants, which are produced by anaerobic bacteria included in Clostridium botulinum group III. In this study, a WGS has been applied to a large collection of C. botulinum group III field strains in order to expand the knowledge on these BoNT-producing Clostridia and to evaluate the potentiality of this method for epidemiological investigations. Sixty field strains were submitted to WGS, and the results were analyzed with respect to epidemiological information and compared to published sequences. The strains were isolated from biological or environmental samples collected in animal botulism outbreaks which occurred in Italy from 2007 to 2016. The new sequenced strains belonged to subspecific groups, some of which were already defined, while others were newly characterized, peculiar to Italian strains and contained genomic features not yet observed. This included, in particular, two new flicC types (VI and VII) and new plasmids which widen the known plasmidome of the species. The extensive genome exploration shown in this study improves the C. botulinum and related species classification scheme, enriching it with new strains of rare genotypes and permitting the highest grade of discrimination among strains for forensic and epidemiological applications.

Characterization of Komagataeibacter Isolate Reveals New Prospects in Waste Stream Valorization for Bacterial Cellulose Production

Komagataeibacter spp. has been used for the bioconversion of industrial wastes and lignocellulosic hydrolysates to bacterial cellulose (BC). Recently, studies have demonstrated the capacity of Komagataeibacter spp. in the biotransformation of inhibitors found in lignocellulosic hydrolysates, aromatic lignin-derived monomers (LDMs) and acetate. In general, detoxification and BC synthesis from lignocellulosic inhibitors requires a carbon flow from acetyl-coA towards tricarboxylic acid and gluconeogenesis, respectively. However, the related molecular aspects have not yet been identified in Komagataeibacter spp. In this study, we isolated a cellulose-producing bacterium capable of synthesizing BC in a minimal medium containing crude glycerol, a by-product from the biodiesel production process. The isolate, affiliated to Komagataeibacter genus, synthesized cellulose in a minimal medium containing glucose (3.3 ± 0.3 g/L), pure glycerol (2.2 ± 0.1 g/L) and crude glycerol (2.1 ± 0.1 g/L). Genome assembly and annotation identified four copies of bacterial cellulose synthase operon and genes for redirecting the carbon from the central metabolic pathway to gluconeogenesis. According to the genome annotations, a BC production route from acetyl-CoA, a central metabolic intermediate, was hypothesized and was validated using acetate. We identified that when K. rhaeticus ENS9b was grown in a minimal medium supplemented with acetate, BC production was not observed. However, in the presence of readily utilizable substrates, such as spent yeast hydrolysate, acetate supplementation improved BC synthesis.

Clinical Molecular and Genomic Epidemiology of Morganella morganii in China

Objectives: Ongoing acquisition of antimicrobial resistance genes has made Morganella morganii a new clinical treatment challenge. Understanding the molecular epidemiology of M. morganii will contribute to clinical treatment and prevention.

Methods: We undertook a 6-year clinical molecular epidemiological investigation of M. morganii from three tertiary hospitals in China since 2014. Antimicrobial susceptibility testing was performed using a VITEK-2 system. All isolates were screened for β-lactam and plasmid-mediated quinolone resistance genes by PCR. Isolates carrying carbapenem-resistant genes were subjected to whole-genome sequencing (WGS). The variation and evolution of these mobile genetic elements (MGEs) were then systematically analyzed.

Results: Among all M. morganii isolates (n = 335), forty (11.9%) were recognized as multidrug resistant strains. qnrD1, aac(6′)-Ib-cr, bla_TEM–104, and bla_{CTX–M–162} were the top four most prevalent resistance genes. Notably, phylogenomic and population structure analysis suggested clade 1 (rhierBAPS SC3 and SC5) associated with multiple resistance genes seemed to be widely spread. WGS showed a bla_OXA–181-carrying IncX3 plasmid and a Proteus genomic island 2 variant carrying bla_CTX–M–3, aac(6′)-Ib-cr coexisted in the same multidrug resistant strain zy_m28. Additionally, a bla_IMP–1-carrying IncP-1β type plasmid was found in the strain nx_m63.

Conclusion: This study indicates a clade of M. morganii is prone to acquire resistance genes, and multidrug resistant M. morganii are increasing by harboring a variety of MGEs including two newly discovered ones in the species. We should be vigilant that M. morganii may bring more extensive and challenging antimicrobial resistance issue.

Genomic Characterization of a Nalidixic Acid-Resistant Salmonella Enteritidis Strain Causing Persistent Infections in Broiler Chickens

Virulent strains of Salmonella enterica subsp. enterica serovar Enteritidis (SE) harbored by poultry can cause disease in poultry flocks and potentially result in human foodborne illness. Two broiler flocks grown a year apart on the same premises experienced mortality throughout the growing period due to septicemic disease caused by SE. Gross lesions predominantly consisted of polyserositis followed by yolk sacculitis, arthritis, osteomyelitis, and spondylitis. Tissues with lesions were cultured yielding 59 SE isolates. These were genotyped by Rep-PCR followed by whole-genome sequencing (WGS) of 15 isolates which were clonal. The strain, SE_TAU19, was further characterized for antimicrobial susceptibility and virulence in a broiler embryo lethality assay. SE_TAU19 was resistant to nalidixic acid and sulfadimethoxine and was virulent to embryos with 100% mortality of all challenged broiler embryos within 3.5 days. Screening the SE_TAU19 whole-genome sequence revealed seven antimicrobial resistance (AMR) genes, 120 virulence genes, and two IncF plasmid replicons corresponding to a single, serovar-specific pSEV virulence plasmid. The pef, spv, and rck virulence genes localized to the plasmid sequence assembly. We report phenotypic and genomic features of a virulent SE strain from persistently infected broiler flocks and present a workflow for SE characterization from isolate collection to genome assembly and sequence analysis. Further SE surveillance and investigation of SE virulence in broiler chickens is warranted.

Comparative Genomic Analysis of Novel Bifidobacterium longum subsp. longum Strains Reveals Functional Divergence in the Human Gut Microbiota

Bifidobacterium longum subsp. longum is a prevalent group in the human gut microbiome. Its persistence in the intestinal microbial community suggests a close host-microbe relationship according to age. The subspecies adaptations are related to metabolic capabilities and genomic and functional diversity. In this study, 154 genomes from public databases and four new Chilean isolates were genomically compared through an in silico approach to identify genomic divergence in genes associated with carbohydrate consumption and their possible adaptations to different human intestinal niches. The pangenome of the subspecies was open, which correlates with its remarkable ability to colonize several niches. The new genomes homogenously clustered within subspecies longum, as observed in phylogenetic analysis. B. longum SC664 was different at the sequence level but not in its functions. COG analysis revealed that carbohydrate use is variable among longum subspecies. Glycosyl hydrolases participating in human milk oligosaccharide use were found in certain infant and adult genomes. Predictive genomic analysis revealed that B. longum M12 contained an HMO cluster associated with the use of fucosylated HMOs but only endowed with a GH95, being able to grow in 2-fucosyllactose as the sole carbon source. This study identifies novel genomes with distinct adaptations to HMOs and highlights the plasticity of B. longum subsp. longum to colonize the human gut microbiota.

Colistin and Carbapenem-Resistant Acinetobacter baumannii Aci46 in Thailand: Genome Analysis and Antibiotic Resistance Profiling

Resistance to the last-line antibiotics against invasive Gram-negative bacterial infection is a rising concern in public health. Multidrug resistant (MDR) Acinetobacter baumannii Aci46 can resist colistin and carbapenems with a minimum inhibitory concentration of 512 µg/mL as determined by microdilution method and shows no zone of inhibition by disk diffusion method. These phenotypic characteristics prompted us to further investigate the genotypic characteristics of Aci46. Next generation sequencing was applied in this study to obtain whole genome data. We determined that Aci46 belongs to Pasture ST2 and is phylogenetically clustered with international clone (IC) II as the predominant strain in Thailand. Interestingly, Aci46 is identical to Oxford ST1962 that previously has never been isolated in Thailand. Two plasmids were identified (pAci46a and pAci46b), neither of which harbors any antibiotic resistance genes but pAci46a carries a conjugational system (type 4 secretion system or T4SS). Comparative genomics with other polymyxin and carbapenem-resistant A. baumannii strains (AC30 and R14) identified shared features such as CzcCBA, encoding a cobalt/zinc/cadmium efflux RND transporter, as well as a drug transporter with a possible role in colistin and/or carbapenem resistance in A. baumannii. Single nucleotide polymorphism (SNP) analyses against MDR ACICU strain showed three novel mutations i.e., Glu229Asp, Pro200Leu, and Ala138Thr, in the polymyxin resistance component, PmrB. Overall, this study focused on Aci46 whole genome data analysis, its correlation with antibiotic resistance phenotypes, and the presence of potential virulence associated factors.

Characterization of Campylobacter spp. Strains Isolated From Wild Birds in Turkey

Turkey is an important stopover site for many migrating birds between Europe, Asia and Africa. Campylobacter spp. are frequently found in wildlife, in particular waterfowl, and distinct strains are disseminated within this reservoir. In this study, 183 wild birds of hunting areas in Turkey were collected and thermophilic Campylobacter spp. from cloacal swabs were isolated at a prevalence of 5.2% from song thrushes (6/116) and 93% from Eurasian coots (41/44). After PCR species differentiation and flaA restriction profiles determination, C. jejuni and C. coli strains were further investigated by whole genome sequencing. PCR target amplification of the ceuE gene, commonly used for C. coli species-identification was inefficient and even hampered in one isolate. A close look on the ceuE sequence revealed that various mismatches in the ceuE oligo annealing sites caused less efficient diagnostic detection. All C. coli isolates belonged to the environmental clade II and clade III, for which thirty-six novel MLST types were identified. Further single nucleotide polymorphism (SNP) analysis showed a high genomic divergence between the C. coli isolates. High variability was also implicated for putative plasmid-located genes detected in 51% of the C. coli isolates. Distinct gene variants in clades II and III C. coli were identified by a k-mer analysis. After substracting k-mers in common with C. coli clade I database, 11 and 35 distinct genes were identified in clades II and III isolates, mainly involved in surface structures and modifications as well as signal transduction, suggesting niche adaptation of C. coli strains in wild birds. All strains were susceptible against (fluoro-)quinolones, erythromycin, tetracycline, gentamicin and only one isolate was resistant against streptomycin, suggesting that the sensitive phenotype was due to absence of selective pressure and niche separation in wild birds in Turkey. We conclude that Campylobacter spp. isolates from wildlife and environmental sources are still scarce in the databases and that there is a need for more studies on thermophilic Campylobacter spp. from different places all over the world in order to complement our understanding on dissemination and adaptation to distinct niches of this global food-borne pathogen.

Functional Genomic Insights into Probiotic Bacillus siamensis Strain B28 from Traditional Korean Fermented Kimchi

Bacillus siamensis strain B28 was previously isolated from traditional Korean fermented kimchi and inhibited expression of the microphthalmia-associated transcription factor and β-catenin in human embryonic kidney 293 cells. Here, we determined the complete genome sequence of strain B28 and compared it with other strains to elucidate its potential probiotic properties. Strain B28 does not contain antibiotic resistance-, hemolysin- or enterotoxin-encoding genes. The genome includes genes related to survival in extreme conditions, adhesion in the gut, and synthesis of the bacteriocin. Considering the potential for enhancement of human health, the strain B28 genome encodes genes related to production of eight essential amino acids, γ-aminobutyric acid, branched-chain fatty acids, γ-glutamyltransferase, and subtilisin. There are genes for the synthesis of uracil, lipoteichoic acid, glutathione, and several reactive oxygen species-scavenging enzymes. Experimentally, strain B28 exhibited sensitivity to eight antibiotics and antibacterial activity against seven foodborne pathogens. B. siamensis B28 is a safe strain with potential for development as a probiotic.

CONSULT: accurate contamination removal using locality-sensitive hashing

A fundamental question appears in many bioinformatics applications: Does a sequencing read belong to a large dataset of genomes from some broad taxonomic group, even when the closest match in the set is evolutionarily divergent from the query? For example, low-coverage genome sequencing (skimming) projects either assemble the organelle genome or compute genomic distances directly from unassembled reads. Using unassembled reads needs contamination detection because samples often include reads from unintended groups of species. Similarly, assembling the organelle genome needs distinguishing organelle and nuclear reads. While k-mer-based methods have shown promise in read-matching, prior studies have shown that existing methods are insufficiently sensitive for contamination detection. Here, we introduce a new read-matching tool called CONSULT that tests whether k-mers from a query fall within a user-specified distance of the reference dataset using locality-sensitive hashing. Taking advantage of large memory machines available nowadays, CONSULT libraries accommodate tens of thousands of microbial species. Our results show that CONSULT has higher true-positive and lower false-positive rates of contamination detection than leading methods such as Kraken-II and improves distance calculation from genome skims. We also demonstrate that CONSULT can distinguish organelle reads from nuclear reads, leading to dramatic improvements in skim-based mitochondrial assemblies.

Comparative Genomic Analysis of Extended-spectrum β-lactamase and mcr-1 Positive Escherichia coli from Gut Microbiota of Healthy Singaporeans

Multidrug resistant (MDR) Escherichia coli strains that carry extended-spectrum β-lactamases (ESBLs) or colistin resistance gene mcr-1 have been identified in the human gut at an increasing incidence worldwide. In this study, we isolated and characterized MDR Enterobacteriaceae from the gut microbiota of healthy Singaporeans and show that the detection rates for ESBL-producing and mcr-positive Enterobacteriaceae are 25.7% (28/109) and 7.3% (8/109), respectively. Whole-genome sequencing analysis of the 37 E. coli isolates assigned them into 25 sequence types and six different phylogroups, suggesting that the MDR E. coli gut colonizers are highly diverse. We then analysed the genetic context of the resistance genes and found that composite transposons played important roles in the co-transfer of bla_CTX-M-15/55 and qnrS1, as well as the acquisition of mcr-1. Furthermore, comparative genomic analysis showed that 12 of the 37 MDR E. coli isolates showed high similarity to ESBL-producing E. coli isolates from raw meat products in local markets. By analyzing the core genome SNPs shared by these isolates, we identified possible clonal transmission of a MDR E. coli clone between human and raw meat, as well as a group of highly similar IncI2 (Delta) plasmids that might be responsible for the dissemination of mcr-1 in a much wider geographic region. Together, these results suggest that antibiotic resistance may be transmitted between different environmental settings by the expansion of MDR E. coli clones, as well as by the dissemination of resistance plasmids. Importance The human gut can harbor both antibiotic resistant and virulent E. coli which may subsequently cause infections. In this study, we found that MDR E. coli isolates from the gut of healthy Singaporeans carry a diverse range of antibiotic resistance mechanisms and virulence factor genes, and are highly diverse to each other. By comparing their genomes with the ESBL-producing E. coli isolates from raw meat products that were sampled at a similar time from local markets, we detected a MDR E. coli clone that was possibly transmitted between humans and raw meat products. Furthermore, we also found that a group of resistance plasmids might be responsible for the dissemination of colistin resistance gene mcr-1 in Singapore, Malaysia and Europe. Our findings call for better countermeasures to block the transmission of antibiotic resistance.

Recovering Escherichia coli Plasmids in the Absence of Long-Read Sequencing Data

The incidence of infections caused by multidrug-resistant E. coli strains has risen in the past years. Antibiotic resistance in E. coli is often mediated by acquisition and maintenance of plasmids. The study of E. coli plasmid epidemiology and genomics often requires long-read sequencing information, but recently a number of tools that allow plasmid prediction from short-read data have been developed. Here, we reviewed 25 available plasmid prediction tools and categorized them into binary plasmid/chromosome classification tools and plasmid reconstruction tools. We benchmarked six tools (MOB-suite, plasmidSPAdes, gplas, FishingForPlasmids, HyAsP and SCAPP) that aim to reliably reconstruct distinct plasmids, with a special focus on plasmids carrying antibiotic resistance genes (ARGs) such as extended-spectrum beta-lactamase genes. We found that two thirds (n = 425, 66.3%) of all plasmids were correctly reconstructed by at least one of the six tools, with a range of 92 (14.58%) to 317 (50.23%) correctly predicted plasmids. However, the majority of plasmids that carried antibiotic resistance genes (n = 85, 57.8%) could not be completely recovered as distinct plasmids by any of the tools. MOB-suite was the only tool that was able to correctly reconstruct the majority of plasmids (n = 317, 50.23%), and performed best at reconstructing large plasmids (n = 166, 46.37%) and ARG-plasmids (n = 41, 27.9%), but predictions frequently contained chromosome contamination (40%). In contrast, plasmidSPAdes reconstructed the highest fraction of plasmids smaller than 18 kbp (n = 168, 61.54%). Large ARG-plasmids, however, were frequently merged with sequences derived from distinct replicons. Available bioinformatic tools can provide valuable insight into E. coli plasmids, but also have important limitations. This work will serve as a guideline for selecting the most appropriate plasmid reconstruction tool for studies focusing on E. coli plasmids in the absence of long-read sequencing data.

SAUTE: sequence assembly using target enrichment

Background

Illumina is the dominant sequencing technology at this time. Short length, short insert size, some systematic biases, and low-level carryover contamination in Illumina reads continue to make assembly of repeated regions a challenging problem. Some applications also require finding multiple well supported variants for assembled regions.

Results

To facilitate assembly of repeat regions and to report multiple well supported variants when a user can provide target sequences to assist the assembly, we propose SAUTE and SAUTE_PROT assemblers. Both assemblers use de Bruijn graph on reads. Targets can be transcripts or proteins for RNA-seq reads and transcripts, proteins, or genomic regions for genomic reads. Target sequences are nucleotide and protein sequences for SAUTE and SAUTE_PROT, respectively.

Conclusions

For RNA-seq, comparisons with Trinity, rnaSPAdes, SPAligner, and SPAdes assembly of reads aligned to target proteins by DIAMOND show that SAUTE_PROT finds more coding sequences that translate to benchmark proteins. Using AMRFinderPlus calls, we find SAUTE has higher sensitivity and precision than SPAdes, plasmidSPAdes, SPAligner, and SPAdes assembly of reads aligned to target regions by HISAT2. It also has better sensitivity than SKESA but worse precision.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04174-9.

Directed Recovery and Molecular Characterization of Antibiotic Resistance Plasmids from Cheese Bacteria

Resistance to antimicrobials is a growing problem of worldwide concern. Plasmids are thought to be major drivers of antibiotic resistance spread. The present work reports a simple way to recover replicative plasmids conferring antibiotic resistance from the bacteria in cheese. Purified plasmid DNA from colonies grown in the presence of tetracycline and erythromycin was introduced into plasmid-free strains of Lactococcus lactis, Lactiplantibacillus plantarum and Lacticaseibacillus casei. Following antibiotic selection, the plasmids from resistant transformants were isolated, analyzed by restriction enzyme digestion, and sequenced. Seven patterns were obtained for the tetracycline-resistant colonies, five from L. lactis, and one each from the lactobacilli strains, as well as a single digestion profile for the erythromycin-resistant transformants obtained in L. lactis. Sequence analysis respectively identified tet(S) and ermB in the tetracycline- and erythromycin-resistance plasmids from L. lactis. No dedicated resistance genes were detected in plasmids conferring tetracycline resistance to L. casei and L. plantarum. The present results highlight the usefulness of the proposed methodology for isolating functional plasmids that confer antibiotic resistance to LAB species, widen our knowledge of antibiotic resistance in the bacteria that inhabit cheese, and emphasize the leading role of plasmids in the spread of resistance genes via the food chain.

Diversity and Evolution of Clostridium beijerinckii and Complete Genome of the Type Strain DSM 791T

Clostridium beijerinckii is a relatively widely studied, yet non-model, bacterium. While 246 genome assemblies of its various strains are available currently, the diversity of the whole species has not been studied, and it has only been analyzed in part for a missing genome of the type strain. Here, we sequenced and assembled the complete genome of the type strain Clostridium beijerinckii DSM 791T, composed of a circular chromosome and a circular megaplasmid, and used it for a comparison with other genomes to evaluate diversity and capture the evolution of the whole species. We found that strains WB53 and HUN142 were misidentified and did not belong to the Clostridium beijerinckii species. Additionally, we filtered possibly misassembled genomes, and we used the remaining 237 high-quality genomes to define the pangenome of the whole species. By its functional annotation, we showed that the core genome contains genes responsible for basic metabolism, while the accessory genome has genes affecting final phenotype that may vary among different strains. We used the core genome to reconstruct the phylogeny of the species and showed its great diversity, which complicates the identification of particular strains, yet hides possibilities to reveal hitherto unreported phenotypic features and processes utilizable in biotechnology.

ORFograph: search for novel insecticidal protein genes in genomic and metagenomic assembly graphs

BACKGROUND

Since the prolonged use of insecticidal proteins has led to toxin resistance, it is important to search for novel insecticidal protein genes (IPGs) that are effective in controlling resistant insect populations. IPGs are usually encoded in the genomes of entomopathogenic bacteria, especially in large plasmids in strains of the ubiquitous soil bacteria, Bacillus thuringiensis (Bt). Since there are often multiple similar IPGs encoded by such plasmids, their assemblies are typically fragmented and many IPGs are scattered through multiple contigs. As a result, existing gene prediction tools (that analyze individual contigs) typically predict partial rather than complete IPGs, making it difficult to conduct downstream IPG engineering efforts in agricultural genomics.

METHODS

Although it is difficult to assemble IPGs in a single contig, the structure of the genome assembly graph often provides clues on how to combine multiple contigs into segments encoding a single IPG.

RESULTS

We describe ORFograph, a pipeline for predicting IPGs in assembly graphs, benchmark it on (meta)genomic datasets, and discover nearly a hundred novel IPGs. This work shows that graph-aware gene prediction tools enable the discovery of greater diversity of IPGs from (meta)genomes.

CONCLUSIONS

We demonstrated that analysis of the assembly graphs reveals novel candidate IPGs. ORFograph identified both already known genes "hidden" in assembly graphs and potential novel IPGs that evaded existing tools for IPG identification. As ORFograph is fast, one could imagine a pipeline that processes many (meta)genomic assembly graphs to identify even more novel IPGs for phenotypic testing than would previously be inaccessible by traditional gene-finding methods. While here we demonstrated the results of ORFograph only for IPGs, the proposed approach can be generalized to any class of genes. Video abstract.

Klebsiella pneumoniae ST307 with OXA-181: threat of a high-risk clone and promiscuous plasmid in a resource-constrained healthcare setting

INTRODUCTION

Klebsiella pneumoniae with OXA-48-like enzymes were introduced into Tshwane Tertiary Hospital (TTH) (Pretoria, South Africa) during September 2015, causing nosocomial outbreaks.

METHODS

PCR methodologies and WGS were used to characterize K. pneumoniae with carbapenemases (n = 124) from TTH (July 2015-December 2016).

RESULTS

PCR was used to track K. pneumoniae ST307 with OXA-181 among 60% of carbapenemase-positive isolates in different wards/units over time and showed the transmission of IncX3 plasmids to other K. pneumoniae clones. WGS identified different ST307 clades: 307_OXA181 (consisting of two lineages, A and B) with OXA-181 on IncX3 plasmids (named p72_X3_OXA181) and clade 307_VIM with VIM-1 on IncFII plasmids. Clade 307_OXA181 lineage B was responsible for the rapid increase and transmission of OXA-181 K. pneumoniae in various wards/units throughout TTH, while the numbers of clade 307_OXA181 lineage A and clade 307_VIM remained low. Separate outbreaks were due to K. pneumoniae ST17 and ST29 with p72_X3_OXA181 plasmids.

CONCLUSIONS

The high-risk clone K. pneumoniae ST307 with OXA-181 rapidly spread to different wards/units despite infection and prevention measures. ST307 clades and lineages seemingly acted differently in outbreak situations. This study also highlighted the threat of promiscuous plasmids such as p72_X3_OXA181.

Phylogenomic Reconstruction and Metabolic Potential of the Genus Aminobacter

Bacteria belonging to the genus Aminobacter are metabolically versatile organisms thriving in both natural and anthropized terrestrial environments. To date, the taxonomy of this genus is poorly defined due to the unavailability of the genomic sequence of A. anthyllidis LMG 26462T and the presence of unclassified Aminobacter strains. Here, we determined the genome sequence of A. anthyllidis LMG 26462T and performed phylogenomic, average nucleotide identity and digital DNA-DNA hybridization analyses of 17 members of genus Aminobacter. Our results indicate that 16S rRNA-based phylogeny does not provide sufficient species-level discrimination, since most of the unclassified Aminobacter strains belong to valid Aminobacter species or are putative new species. Since some members of the genus Aminobacter can utilize certain C1 compounds, such as methylamines and methyl halides, a comparative genomic analysis was performed to characterize the genetic basis of some degradative/assimilative pathways in the whole genus. Our findings suggest that all Aminobacter species are heterotrophic methylotrophs able to generate the methylene tetrahydrofolate intermediate through multiple oxidative pathways of C1 compounds and convey it in the serine cycle. Moreover, all Aminobacter species carry genes implicated in the degradation of phosphonates via the C-P lyase pathway, whereas only A. anthyllidis LMG 26462T contains a symbiosis island implicated in nodulation and nitrogen fixation.

Comparative genomics and in silico gene evaluation involved in the probiotic potential of Bifidobacterium longum 51A

The Bifidobacterium longum 5^1A strain of isolated from feces of a healthy child, has demonstrated probiotic properties by in vivo and in vitro studies, which may be assigned to its production of metabolites such as acetate. Thus, through the study of comparative genomics, the present work sought to identify unique genes that might be related to the production of acetate. To perform the study, the DNA strain was sequenced using Illumina HiSeq technology, followed by assembly and manual curation of coding sequences. Comparative analysis was performed including 19 complete B. longum genomes available in Genbank/NCBI. In the phylogenetic analysis, the CECT 7210 and 157F strains of B. longum subsp. infantis aggregated within the subsp. longum cluster, suggesting that their taxonomic classification should be reviewed. The strain 5^1A of B. longum has 26 unique genes, six of which are possibly related to carbohydrate metabolism and acetate production. The phosphoketolase pathway from B. longum 5^1A showed a difference in acetyl-phosphate production. This result seems to corroborate the analysis of their unique genes, whose presence suggests the strain may use different sources of carbohydrates that allow a greater production of acetate and consequently offer benefits to the host health.

Polyphenism of a Novel Trait Integrated Rapidly Evolving Genes into Ancestrally Plastic Networks

Developmental polyphenism, the ability to switch between phenotypes in response to environmental variation, involves the alternating activation of environmentally sensitive genes. Consequently, to understand how a polyphenic response evolves requires a comparative analysis of the components that make up environmentally sensitive networks. Here, we inferred coexpression networks for a morphological polyphenism, the feeding-structure dimorphism of the nematode Pristionchus pacificus. In this species, individuals produce alternative forms of a novel trait—moveable teeth, which in one morph enable predatory feeding—in response to environmental cues. To identify the origins of polyphenism network components, we independently inferred coexpression modules for more conserved transcriptional responses, including in an ancestrally nonpolyphenic nematode species. Further, through genome-wide analyses of these components across the nematode family (Diplogastridae) in which the polyphenism arose, we reconstructed how network components have changed. To achieve this, we assembled and resolved the phylogenetic context for five genomes of species representing the breadth of Diplogastridae and a hypothesized outgroup. We found that gene networks instructing alternative forms arose from ancestral plastic responses to environment, specifically starvation-induced metabolism and the formation of a conserved diapause (dauer) stage. Moreover, loci from rapidly evolving gene families were integrated into these networks with higher connectivity than throughout the rest of the P. pacificus transcriptome. In summary, we show that the modular regulatory outputs of a polyphenic response evolved through the integration of conserved plastic responses into networks with genes of high evolutionary turnover.

Genomics of Atlantic Forest Mycobacteriaceae strains unravels a mobilome diversity with a novel integrative conjugative element and plasmids harbouring T7SS

Mobile genetic elements (MGEs) are agents of bacterial evolution and adaptation. Genome sequencing provides an unbiased approach that has revealed an abundance of MGEs in prokaryotes, mainly plasmids and integrative conjugative elements. Nevertheless, many mobilomes, particularly those from environmental bacteria, remain underexplored despite their representing a reservoir of genes that can later emerge in the clinic. Here, we explored the mobilome of the Mycobacteriaceae family, focusing on strains from Brazilian Atlantic Forest soil. Novel Mycolicibacterium and Mycobacteroides strains were identified, with the former ones harbouring linear and circular plasmids encoding the specialized type-VII secretion system (T7SS) and mobility-associated genes. In addition, we also identified a T4SS-mediated integrative conjugative element (ICEMyc226) encoding two T7SSs and a number of xenobiotic degrading genes. Our study uncovers the diversity of the Mycobacteriaceae mobilome, providing the evidence of an ICE in this bacterial family. Moreover, the presence of T7SS genes in an ICE, as well as plasmids, highlights the role of these mobile genetic elements in the dispersion of T7SS.

Investigation of Plasmids Among Clinical Staphylococcus aureus and Staphylococcus haemolyticus Isolates From Egypt

Staphylococci can cause a wide array of infections that can be life threatening. These infections become more deadly when the isolates are antibiotic resistant and thus harder to treat. Many resistance determinants are plasmid-mediated; however, staphylococcal plasmids have not yet been fully characterized. In particular, plasmids and their contributions to antibiotic resistance have not been investigated within the Arab states, where antibiotic use is not universally regulated. Here, we characterized the putative plasmid content among 56 Staphylococcus aureus and 10 Staphylococcus haemolyticus clinical isolates from Alexandria, Egypt. Putative plasmid sequences were detected in over half of our collection. In total, we identified 72 putative plasmid sequences in 27 S. aureus and 1 S. haemolyticus isolates. While these isolates typically carried one or two plasmids, we identified one isolate-S. aureus AA53-with 11 putative plasmids. The plasmid sequences most frequently encoded a Rep_1, RepL, or PriCT_1 type replication protein. As expected, antibiotic resistance genes were widespread among the identified plasmid sequences. Related plasmids were identified amongst our clinical isolates; homologous plasmids present in multiple isolates clustered into 11 groups based upon sequence similarity. Plasmids from the same cluster often shared antibiotic resistance genes, including blaZ, which is associated with β-lactam resistance. Our analyses suggest that plasmids are a key factor in the pathology and epidemiology of S. aureus in Egypt. A better characterization of plasmids and the role they contribute to the success of Staphylococci as pathogens will guide the design of effective control strategies to limit their spread.

A Peek into the Plasmidome of Global Sewage

Plasmids can provide a selective advantage for microorganisms to survive and adapt to new environmental conditions. Plasmid-encoded traits, such as antimicrobial resistance (AMR) or virulence, impact the ecology and evolution of bacteria and can significantly influence the burden of infectious diseases. Insight about the identity and functions encoded on plasmids on the global scale are largely lacking. Here, we investigate the plasmidome of 24 samples (22 countries, 5 continents) from the global sewage surveillance project. We obtained 105-Gbp Oxford Nanopore and 167-Gbp Illumina NextSeq DNA sequences from plasmid DNA preparations and assembled 165,302 contigs (159,322 circular). Of these, 58,429 carried genes encoding for plasmid-related and 11,222 for virus/phage-related proteins. About 90% of the circular DNA elements did not have any similarity to known plasmids. Those that exhibited similarity had similarity to plasmids whose hosts were previously detected in these sewage samples (e.g., Acinetobacter, Escherichia, Moraxella, Enterobacter, Bacteroides, and Klebsiella). Some AMR classes were detected at a higher abundance in plasmidomes (e.g., macrolide-lincosamide-streptogramin B, macrolide, and quinolone) compared to the respective complex sewage samples. In addition to AMR genes, a range of functions were encoded on the candidate plasmids, including plasmid replication and maintenance, mobilization, and conjugation. In summary, we describe a laboratory and bioinformatics workflow for the recovery of plasmids and other potential extrachromosomal DNA elements from complex microbiomes. Moreover, the obtained data could provide further valuable insight into the ecology and evolution of microbiomes, knowledge about AMR transmission, and the discovery of novel functions.

IMPORTANCE This is, to the best of our knowledge, the first study to investigate plasmidomes at a global scale using long read sequencing from complex untreated domestic sewage. Previous metagenomic surveys have detected AMR genes in a variety of environments, including sewage. However, it is unknown whether the AMR genes were present on the microbial chromosome or located on extrachromosomal elements, such as plasmids. Using our approach, we recovered a large number of plasmids, of which most appear novel. We identified distinct AMR genes that were preferentially located on plasmids, potentially contributing to their transmissibility. Overall, plasmids are of great importance for the biology of microorganisms in their natural environments (free-living and host-associated), as well as for molecular biology and biotechnology. Plasmidome collections may therefore be valuable resources for the discovery of fundamental biological mechanisms and novel functions useful in a variety of contexts.

Antibiotic resistance: Time of synthesis in a post-genomic age

Antibiotic resistance has been highlighted by international organizations, including World Health Organization, World Bank and United Nations, as one of the most relevant global health problems. Classical approaches to study this problem have focused in infected humans, mainly at hospitals. Nevertheless, antibiotic resistance can expand through different ecosystems and geographical allocations, hence constituting a One-Health, Global-Health problem, requiring specific integrative analytic tools. Antibiotic resistance evolution and transmission are multilayer, hierarchically organized processes with several elements (from genes to the whole microbiome) involved. However, their study has been traditionally gene-centric, each element independently studied. The development of robust-economically affordable whole genome sequencing approaches, as well as other -omic techniques as transcriptomics and proteomics, is changing this panorama. These technologies allow the description of a system, either a cell or a microbiome as a whole, overcoming the problems associated with gene-centric approaches. We are currently at the time of combining the information derived from -omic studies to have a more holistic view of the evolution and spread of antibiotic resistance. This synthesis process requires the accurate integration of -omic information into computational models that serve to analyse the causes and the consequences of acquiring AR, fed by curated databases capable of identifying the elements involved in the acquisition of resistance. In this review, we analyse the capacities and drawbacks of the tools that are currently in use for the global analysis of AR, aiming to identify the more useful targets for effective corrective interventions.

Molecular Epidemiology of Third-Generation-Cephalosporin-Resistant Enterobacteriaceae in Southeast Queensland, Australia

Third-generation cephalosporin-resistant (3GC-R) Enterobacteriaceae represent a major threat to human health. Here, we captured 288 3GC-R Enterobacteriaceae clinical isolates from 264 patients presenting at a regional Australian hospital over a 14-month period. In addition to routine mass spectrometry and antibiotic sensitivity testing, isolates were examined using rapid (∼40-min) real-time PCR assays targeting the most common extended-spectrum β-lactamases (ESBLs; bla _CTX-M-1 and bla _CTX-M-9 groups, plus bla _TEM, bla _SHV, and an internal 16S rRNA gene control). AmpC CMY β-lactamase (bla _CMY) prevalence was also examined. Escherichia coli (80.2%) and Klebsiella pneumoniae (17.0%) were dominant, with Klebsiella oxytoca, Klebsiella aerogenes, and Enterobacter cloacae infrequently identified. Ceftriaxone and cefoxitin resistance were identified in 97.0% and 24.5% of E. coli and K. pneumoniae isolates, respectively. Consistent with global findings in Enterobacteriaceae, most (98.3%) isolates harbored at least one β-lactamase gene, with 144 (50%) harboring bla _CTX-M-1 group, 92 (31.9%) harboring bla _CTX-M-9 group, 48 (16.7%) harboring bla _SHV, 133 (46.2%) harboring bla _TEM, and 34 (11.8%) harboring bla _CMY genes. A subset of isolates (n = 98) were subjected to whole-genome sequencing (WGS) to identify the presence of cryptic resistance determinants and to verify genotyping accuracy. WGS of β-lactamase-negative or carbapenem-resistant isolates identified uncommon ESBL and carbapenemase genes, including bla _NDM and bla _IMP, and confirmed all PCR-positive genotypes. We demonstrate that our PCR assays enable the rapid and cost-effective identification of ESBLs in the hospital setting, which has important infection control and therapeutic implications.