Flexible protein database based on amino acid k-mers.

Identification of proteins is one of the most computationally intensive steps in genomics studies. It usually relies on aligners that do not accommodate rich information on proteins and require additional pipelining steps for protein identification. We introduce kAAmer, a protein database engine based on amino-acid k-mers that provides efficient identification of proteins while supporting the incorporation of flexible annotations on these proteins. Moreover, the database is built to be used as a microservice, to be hosted and queried remotely.

Complete Genome Sequences of Klebsiella michiganensis and Citrobacter farmeri, KPC-2-Producers Serially Isolated from a Single Patient.

Carbapenemase-producing Enterobacterales, including KPC-2 producers, have become a major clinical problem. During an outbreak in Quebec City, Canada, KPC-2-producing Klebsiella michiganensis and Citrobacter farmeri were isolated from a patient six weeks apart. We determined their complete genome sequences. Both isolates carried nearly identical IncN2 plasmids with blaKPC-2 on a Tn4401b element. Both strains also carried IncP1 plasmids, but that of C. farmeri did not carry a Beta-lactamase gene, whereas that of K. michiganensis carried a second copy of blaKPC-2 on Tn4401b. These results suggest recent plasmid transfer between the two species and a recent transposition event.

Serratia marcescens SCH909 as reservoir and source of genetic elements related to wide dissemination of antimicrobial resistance mechanisms.

Serratia marcescens SCH909 is a multidrug resistant strain isolated in 1988 harboring three class 1 integrons. We wondered if these integrons were retained over time and if there were other antimicrobial resistant determinants contributing to its multidrug resistant profile. Genomic analysis showed a fourth multidrug resistance integron, a Tn7 transposon with dfrA1-sat2-ybeA-ybfA-ybfB-ybgA gene cassettes in the variable region. Insertion sequences were involved in the genesis of novel composite transposons in the L4 subtype plasmid pSCH909, such as Tn6824 carrying an arsenic regulon and two head to head class 1 integrons surrounded by two complete IS1. Remarkably, a novel chromosomal genomic island, SmaR, was identified, closely related to Multiple Antimicrobial Resistance Regions (MARR), usually found in AbaR0-type and AbGRI2-0 from global clones of Acinetobacter baumannii, and in M-type plasmids circulating in Enterobacteriaceae. Maintenance studies showed that the three class 1 integrons were maintained over 1 month without antimicrobial pressure. Since S. marcescens is considered a relevant nosocomial pathogen that can have a wide range of niches - human, plant, animal, soil and inanimate surfaces, our findings support the ability of this species to capture, maintain and spread a broad variety of antimicrobial resistance elements.

Comment on "Conserved phylogenetic distribution and limited antibiotic resistance of class 1 integrons revealed by assessing the bacterial genome and plasmid collection" by A.N. Zhang et al.

An article published in Microbiome in July 2018 uses incorrect definitions of integron integrase IntI1 and of class 1 integrons that affect the interpretation of the data.

Genome Sequence of a Klebsiella pneumoniae NDM-1 Producer Isolated in Quebec City.

We report here the complete genome sequence of Klebsiella pneumoniae CCRI-22199, isolated from a patient from India treated in Quebec City, Canada. Genes encoding beta-lactamases NDM-1 and CTX-M-15 were identified on two distinct plasmids. While the chromosome is similar to that of strain BAA-2146, CCRI-22199 provides a further example of rearrangements in plasmids.

Culture-enriched human gut microbiomes reveal core and accessory resistance genes.

BACKGROUND: Low-abundance microorganisms of the gut microbiome are often referred to as a reservoir for antibiotic resistance genes. Unfortunately, these less-abundant bacteria can be overlooked by deep shotgun sequencing. In addition, it is a challenge to associate the presence of resistance genes with their risk of acquisition by pathogens. In this study, we used liquid culture enrichment of stools to assemble the genome of lower-abundance bacteria from fecal samples. We then investigated the gene content recovered from these culture-enriched and culture-independent metagenomes in relation with their taxonomic origin, specifically antibiotic resistance genes. We finally used a pangenome approach to associate resistance genes with the core or accessory genome of Enterobacteriaceae and inferred their propensity to horizontal gene transfer. RESULTS: Using culture-enrichment approaches with stools allowed assembly of 187 bacterial species with an assembly size greater than 1 million nucleotides. Of these, 67 were found only in culture-enriched conditions, and 22 only in culture-independent microbiomes. These assembled metagenomes allowed the evaluation of the gene content of specific subcommunities of the gut microbiome. We observed that differentially distributed metabolic enzymes were associated with specific culture conditions and, for the most part, with specific taxa. Gene content differences between microbiomes, for example, antibiotic resistance, were for the most part not associated with metabolic enzymes, but with other functions. We used a pangenome approach to determine if the resistance genes found in Enterobacteriaceae, specifically E. cloacae or E. coli, were part of the core genome or of the accessory genome of this species. In our healthy volunteer cohort, we found that E. cloacae contigs harbored resistance genes that were part of the core genome of the species, while E. coli had a large accessory resistome proximal to mobile elements. CONCLUSION: Liquid culture of stools contributed to an improved functional and comparative genomics study of less-abundant gut bacteria, specifically those associated with antibiotic resistance. Defining whether a gene is part of the core genome of a species helped in interpreting the genomes recovered from culture-independent or culture-enriched microbiomes.

Criibacterium bergeronii gen. nov., sp. nov., a new member of the family Peptostreptococcaceae, isolated from human clinical samples.

A rod-shaped, motile anaerobic bacterium, designated CCRI-22567T, was isolated from a vaginal sample of a woman diagnosed with bacterial vaginosis and subjected to a polyphasic taxonomic study. The novel strain was capable of growth at 30-42 °C (optimum, 42 °C), at pH 5.5-8.5 (optimum, pH 7.0-7.5) and in the presence of 0-1.5 % (w/v) NaCl (optimally at 0.5 % NaCl). The phylogenetic trees based on 16S rRNA gene sequences showed that strain CCRI-22567T forms a distinct evolutionary lineage independent of other taxa in the family Peptostreptococcaceae. Strain CCRI-22567T exhibited 90.1 % 16S rRNA gene sequence similarity to Peptoanaerobacter stomatis ACC19aT and 89.7 % to Eubacterium yurii subsp. schtitka ATCC 43716. The three closest organisms with an available whole genome were compared to strain CCRI-22567T for genomic relatedness assessment. The genomic average nucleotide identities (OrthoANIu) obtained with Peptoanaerobacter stomatis ACC19aT, Eubacterium yurii subsp. margaretiae ATCC 43715 and Filifactor alocis ATCC 35896T were 71.8, 70.3 and 69.6 %, respectively. Strain CCRI-22567T contained C18 : 1 ω9c and C18 : 1 ω9c DMA as the major fatty acids. The DNA G+C content of strain CCRI-22567T based on its genome sequence was 33.8 mol%. On the basis of the phylogenetic, chemotaxonomic and other phenotypic properties, strain CCRI-22567T is considered to represent a new genus and species within the family Peptostreptococcaceae, for which the name Criibacterium bergeronii gen. nov., sp. nov., is proposed. The type strain of Criibacterium bergeronii is CCRI-22567T (=LMG 31278T=DSM 107614T=CCUG 72594T).

Phenetic Comparison of Prokaryotic Genomes Using k-mers.

Bacterial genomics studies are getting more extensive and complex, requiring new ways to envision analyses. Using the Ray Surveyor software, we demonstrate that comparison of genomes based on their k-mer content allows reconstruction of phenetic trees without the need of prior data curation, such as core genome alignment of a species. We validated the methodology using simulated genomes and previously published phylogenomic studies of Streptococcus pneumoniae and Pseudomonas aeruginosa. We also investigated the relationship of specific genetic determinants with bacterial population structures. By comparing clusters from the complete genomic content of a genome population with clusters from specific functional categories of genes, we can determine how the population structures are correlated. Indeed, the strain clustering based on a subset of k-mers allows determination of its similarity with the whole genome clusters. We also applied this methodology on 42 species of bacteria to determine the correlational significance of five important bacterial genomic characteristics. For example, intrinsic resistance is more important in P. aeruginosa than in S. pneumoniae, and the former has increased correlation of its population structure with antibiotic resistance genes. The global view of the pangenome of bacteria also demonstrated the taxa-dependent interaction of population structure with antibiotic resistance, bacteriophage, plasmid, and mobile element k-mer data sets.

Complete Genome of a Panresistant Pseudomonas aeruginosa Strain, Isolated from a Patient with Respiratory Failure in a Canadian Community Hospital.

We report here the complete genome sequence of a panresistant Pseudomonas aeruginosa strain, isolated from a patient with respiratory failure in Canada. No carbapenemase genes were identified. Carbapenem resistance is attributable to a frameshift in the oprD gene; the basis for colistin resistance remains undetermined.

Use of phylogenetical analysis to predict susceptibility of pathogenic Candida spp. to antifungal drugs.

Successful treatment of a Candida infection relies on 1) an accurate identification of the pathogenic fungus and 2) on its susceptibility to antifungal drugs. In the present study we investigated the level of correlation between phylogenetical evolution and susceptibility of pathogenic Candida spp. to antifungal drugs. For this, we compared a phylogenetic tree, assembled with the concatenated sequences (2475-bp) of the ATP2, TEF1, and TUF1 genes from 20 representative Candida species, with published minimal inhibitory concentrations (MIC) of the four principal antifungal drug classes commonly used in the treatment of candidiasis: polyenes, triazoles, nucleoside analogues, and echinocandins. The phylogenetic tree revealed three distinct phylogenetic clusters among Candida species. Species within a given phylogenetic cluster have generally similar susceptibility profiles to antifungal drugs and species within Clusters II and III were less sensitive to antifungal drugs than Cluster I species. These results showed that phylogenetical relationship between clusters and susceptibility to several antifungal drugs could be used to guide therapy when only species identification is available prior to information pertaining to its resistance profile. An extended study comprising a large panel of clinical samples should be conducted to confirm the efficiency of this approach in the treatment of candidiasis.

Inhibition of Helicobacter pylori Glu-tRNAGln amidotransferase by novel analogues of the putative transamidation intermediate.

Glutaminyl-tRNAGln in Helicobacter pylori is formed by an indirect route requiring a noncanonical glutamyl-tRNA synthetase and a tRNA-dependent heterotrimeric amidotransferase (AdT) GatCAB. Widespread use of this pathway among prominent human pathogens, and its absence in the mammalian cytoplasm, identify AdT as a target for the development of antimicrobial agents. We present here the inhibitory properties of three dipeptide-like sulfone-containing compounds analogous to the transamidation intermediates, which are competitive inhibitors of AdT with respect to Glu-tRNAGln . Molecular docking revealed that AdT inhibition by these compounds depends on π-π stacking interactions between their aromatic groups and Tyr81 of the GatB subunit. The properties of these inhibitors indicate that the 3'-terminal adenine of Glu-tRNAGln plays a major role in binding to the AdT transamidation active site.

Draft Genome Sequence of Criibacterium bergeronii gen. nov., sp. nov., Strain CCRI-22567T, Isolated from a Vaginal Sample from a Woman with Bacterial Vaginosis.

Criibacterium bergeronii gen. nov., sp. nov., CCRI-22567 is the type strain of the new genus Criibacterium The strain was isolated from a woman with bacterial vaginosis. The genome assembly comprised 2,384,460 bp, with 34.4% G+C content. This is the first genome announcement of a strain belonging to the genus Criibacterium.

Genome and Plasmid Analysis of blaIMP-4-Carrying Citrobacter freundii B38.

Sequencing of the blaIMP-4-carrying C. freundii B38 using the PacBio SMRT technique revealed that the genome contained a chromosome of 5,134,500 bp and three plasmids, pOZ172 (127,005 bp), pOZ181 (277,592 bp), and pOZ182 (18,467 bp). Plasmid pOZ172 was identified as IncFIIY, like pP10164-NDM and pNDM-EcGN174. It carries a class 1 integron with four cassettes (blaIMP-4-qacG2-aacA4-aphA15) and a complete hybrid tni module (tniR-tniQ-tniB-tniA). The recombination of tniR from Tn402 (identical) with tniQBA from Tn5053 (99%) occurred within the res site of Tn402/5053 The Tn402/5053-like integron, named Tn6017, was inserted into Tn1722 at the res II site. The replication, partitioning, and transfer systems of pOZ181 were similar to those of IncHI2 plasmids (e.g., R478) and contained a sul1-type class 1 integron with the cassette array orf-dfrA1-orf-gcu37-aadA5 linked to an upstream Tn1696 tnpA-tnpR and to a downstream 3' conserved sequence (3'-CS) and ISCR1 A Tn2 transposon encoding a blaTEM-1 ß-lactamase was identified on pOZ182. Other interesting resistance determinants encoded on the B38 chromosome included multidrug resistance (MDR) efflux pumps, an AmpC ß-lactamase, and resistances to Cu, Ag, As, and Zn. This is the first report of a complete tni module linked to a blaIMP-4-carrying class 1 integron, which, together with other recently reported non-sul1 integrons, represents the emergence of a distinct evolutionary lineage of class 1 integrons lacking a 3'-CS (qacEΔ1-sul1). The unique cassette array, complete tni module of Tn6017, and incompatibility group of pOZ172 suggest a blaIMP-4 evolutionary pathway in C. freundii B38 different from that for other blaIMP-4 genes found in Gram-negative bacteria in the Western Pacific region.

The initial state of the human gut microbiome determines its reshaping by antibiotics.

Microbiome studies have demonstrated the high inter-individual diversity of the gut microbiota. However, how the initial composition of the microbiome affects the impact of antibiotics on microbial communities is relatively unexplored. To specifically address this question, we administered a second-generation cephalosporin, cefprozil, to healthy volunteers. Stool samples gathered before antibiotic exposure, at the end of the treatment and 3 months later were analysed using shotgun metagenomic sequencing. On average, 15 billion nucleotides were sequenced for each sample. We show that standard antibiotic treatment can alter the gut microbiome in a specific, reproducible and predictable manner. The most consistent effect of the antibiotic was the increase of Lachnoclostridium bolteae in 16 out of the 18 cefprozil-exposed participants. Strikingly, we identified a subgroup of participants who were enriched in the opportunistic pathogen Enterobacter cloacae after exposure to the antibiotic, an effect linked to lower initial microbiome diversity and to a Bacteroides enterotype. Although the resistance gene content of participants' microbiomes was altered by the antibiotic, the impact of cefprozil remained specific to individual participants. Resistance genes that were not detectable prior to treatment were observed after a 7-day course of antibiotic administration. Specifically, point mutations in beta-lactamase blaCfxA-6 were enriched after antibiotic treatment in several participants. This suggests that monitoring the initial composition of the microbiome before treatment could assist in the prevention of some of the adverse effects associated with antibiotics or other treatments.

The resistome of Pseudomonas aeruginosa in relationship to phenotypic susceptibility.

Many clinical isolates of Pseudomonas aeruginosa cause infections that are difficult to eradicate due to their resistance to a wide variety of antibiotics. Key genetic determinants of resistance were identified through genome sequences of 390 clinical isolates of P. aeruginosa, obtained from diverse geographic locations collected between 2003 and 2012 and were related to microbiological susceptibility data for meropenem, levofloxacin, and amikacin. ß-Lactamases and integron cassette arrangements were enriched in the established multidrug-resistant lineages of sequence types ST111 (predominantly O12) and ST235 (O11). This study demonstrates the utility of next-generation sequencing (NGS) in defining relevant resistance elements and highlights the diversity of resistance determinants within P. aeruginosa. This information is valuable in furthering the design of diagnostics and therapeutics for the treatment of P. aeruginosa infections.

Draft Genome Sequence of an International Clonal Lineage 1 Acinetobacter baumannii Strain from Argentina.

In the last few years Acinetobacter baumannii has emerged worldwide as an important nosocomial pathogen in medical institutions. Here, we present the draft genome sequence of the international clonal lineage 1 (ICL1) A. baumannii strain A144 that was isolated in a hospital in Buenos Aires City in the year 1997. The strain is susceptible to carbapenems and resistant to trimethoprim and gentamicin.

Genomic analysis of Pseudomonas aeruginosa PA96, the host of carbapenem resistance plasmid pOZ176.

Pseudomonas aeruginosa PA96 is a clinical isolate from Guangzhou, China, that is multiresistant to antibiotics. We previously described the 500-kb IncP-2 plasmid, pOZ176 that encodes many resistance genes including the IMP-9 carbapenemase. Whole-genome sequencing of PA96 enabled characterization of its genomic islands, virulence factors, and chromosomal resistance genes. We filled gaps using PCR and used optical mapping to confirm the correct contig order. We automatically annotated the core genome and manually annotated the genomic islands. The genome is 6 444 091 bp and encodes 5853 ORFs. From the whole-genome sequence, we constructed a physical map and constructed a phylogenetic tree for comparison with sequenced P. aeruginosa strains. Analysis of known core genome virulence factors and resistance genes revealed few differences with other strains, but the major virulence island is closer to that of DK2 than to PA14. PA96 most closely resembles the environmental strain M18, and notably shares a common serotype, pyoverdin type, flagellar operon, type IV pilin, and several genomic islands with M18.

"Distribution and functional identification of complex class 1 integrons".

The emergence of extended-spectrum ß-lactamases and plasmid-mediated resistance to quinolones has been previously found to be associated with the dissemination of complex class 1 integrons in Argentina. In this study, we analyzed their distribution through time and evaluated the functionality of the Orf513 protein, which is the putative recombinase of the ISCR1 mobile element. We investigated the presence of the orf513, blaCTX-M-2, dfrA3b, qnrB10 and blaDHA-1 genes by PCR and DNA sequencing as well as their linkage to class 1 integrons in 451 non-epidemiologically related nosocomial strains resistant to at least one expanded-spectrum cephalosporin and to one aminoglycoside, isolated between 1989 and 2010 from 7 hospitals from Buenos Aires City. The epidemiology of complex class 1 integrons was found to be notably different among fermenting (94/171) and non-fermenting clinical bacilli isolates (1/280). The ISCR1::qnrB10 positive isolates were found since 1993, confirming its presence in clinical isolates more than a decade before its first description. As expected, In35::ISCR1::blaCTX-M-2 was the most common complex class 1 integron among Enterobacteriaceae isolates, particularly in Proteus mirabilis. Experimental analysis corroborated the activity of the Orf513 protein, which was found to bind specific DNA sequences containing the previously suggested oriIS region. These findings showed the high dispersion and maintenance of complex class 1 integrons across time in our nosocomial isolates. The contribution of the ISCR1 mobile element to multidrug resistant phenotypes is significant due to its sustained association to class 1 integrons.