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.

Genomic diversity and CRISPR-Cas systems in the cyanobacterium Nostoc in the High Arctic.

Nostoc (Nostocales, Cyanobacteria) has a global distribution in the Polar Regions. However, the genomic diversity of Nostoc is little known and there are no genomes available for polar Nostoc. Here we carried out the first genomic analysis of the Nostoc commune morphotype with a recent sample from the High Arctic and a herbarium specimen collected during the British Arctic Expedition (1875-76). Comparisons of the polar genomes with 26 present-day non-polar members of the Nostocales family highlighted that there are pronounced genetic variations among Nostoc strains and species. Osmoprotection and other stress genes were found in all Nostoc strains, but the two Arctic strains had markedly higher numbers of biosynthetic gene clusters for uncharacterised non-ribosomal peptide synthetases, suggesting a high diversity of secondary metabolites. Since viral-host interactions contribute to microbial diversity, we analysed the CRISPR-Cas systems in the Arctic and two temperate Nostoc species. There were a large number of unique repeat-spacer arrays in each genome, indicating diverse histories of viral attack. All Nostoc strains had a subtype I-D system, but the polar specimens also showed evidence of a subtype I-B system that has not been previously reported in cyanobacteria, suggesting diverse cyanobacteria-virus interactions in the Arctic.

Machine learning analysis identifies genes differentiating triple negative breast cancers.

Triple negative breast cancer (TNBC) is one of the most aggressive form of breast cancer (BC) with the highest mortality due to high rate of relapse, resistance, and lack of an effective treatment. Various molecular approaches have been used to target TNBC but with little success. Here, using machine learning algorithms, we analyzed the available BC data from the Cancer Genome Atlas Network (TCGA) and have identified two potential genes, TBC1D9 (TBC1 domain family member 9) and MFGE8 (Milk Fat Globule-EGF Factor 8 Protein), that could successfully differentiate TNBC from non-TNBC, irrespective of their heterogeneity. TBC1D9 is under-expressed in TNBC as compared to non-TNBC patients, while MFGE8 is over-expressed. Overexpression of TBC1D9 has a better prognosis whereas overexpression of MFGE8 correlates with a poor prognosis. Protein-protein interaction analysis by affinity purification mass spectrometry (AP-MS) and proximity biotinylation (BioID) experiments identified a role for TBC1D9 in maintaining cellular integrity, whereas MFGE8 would be involved in various tumor survival processes. These promising genes could serve as biomarkers for TNBC and deserve further investigation as they have the potential to be developed as therapeutic targets for TNBC.

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.

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.

Predictive computational phenotyping and biomarker discovery using reference-free genome comparisons.

BACKGROUND: The identification of genomic biomarkers is a key step towards improving diagnostic tests and therapies. We present a reference-free method for this task that relies on a k-mer representation of genomes and a machine learning algorithm that produces intelligible models. The method is computationally scalable and well-suited for whole genome sequencing studies. RESULTS: The method was validated by generating models that predict the antibiotic resistance of C. difficile, M. tuberculosis, P. aeruginosa, and S. pneumoniae for 17 antibiotics. The obtained models are accurate, faithful to the biological pathways targeted by the antibiotics, and they provide insight into the process of resistance acquisition. Moreover, a theoretical analysis of the method revealed tight statistical guarantees on the accuracy of the obtained models, supporting its relevance for genomic biomarker discovery. CONCLUSIONS: Our method allows the generation of accurate and interpretable predictive models of phenotypes, which rely on a small set of genomic variations. The method is not limited to predicting antibiotic resistance in bacteria and is applicable to a variety of organisms and phenotypes. Kover, an efficient implementation of our method, is open-source and should guide biological efforts to understand a plethora of phenotypes ( http://github.com/aldro61/kover/ ).

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.

Partial recovery of microbiomes after antibiotic treatment.

Antibiotics profoundly affect the gut microbiome and modulate microbial communities. We recently observed that antimicrobial drugs also impact the abundance and distribution of antibiotic resistance genes. In this addendum, we reanalyze our ∼1 trillion nucleotide shotgun metagenomic dataset to quantify comprehensive genomic differences at the sequence level before and after antibiotic treatment. We show that 7 day exposure to cefprozil leads to a statistically significant loss of metagenome sequences. Recovery of gut microbiomes 3 months after antibiotherapy was characterized by the emergence of new genome sequences not observed prior to antibiotic exposure. Participants with low initial gut microbiome diversity had an increased amount of sequences related to antibiotic resistance. Therefore, we suggest that while the taxonomical composition of microbiomes is partially affected by the antibiotic, the genomic content and population structure of bacterial communities is noticeably impacted.

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.

Complete sequence of pOZ176, a 500-kilobase IncP-2 plasmid encoding IMP-9-mediated carbapenem resistance, from outbreak isolate Pseudomonas aeruginosa 96.

Pseudomonas aeruginosa 96 (PA96) was isolated during a multicenter surveillance study in Guangzhou, China, in 2000. Whole-genome sequencing of this outbreak strain facilitated analysis of its IncP-2 carbapenem-resistant plasmid, pOZ176. The plasmid had a length of 500,839 bp and an average percent G+C content of 57%. Of the 618 predicted open reading frames, 65% encode hypothetical proteins. The pOZ176 backbone is not closely related to any plasmids thus far sequenced, but some similarity to pQBR103 of Pseudomonas fluorescens SBW25 was observed. Two multiresistant class 1 integrons and several insertion sequences were identified. The blaIMP-9-carrying integron contained aacA4 → bla(IMP-9) → aacA4, flanked upstream by Tn21 tnpMRA and downstream by a complete tni operon of Tn402 and a mer module, named Tn6016. The second integron carried aacA4 → catB8a → bla(OXA-10) and was flanked by Tn1403-like tnpRA and a sul1-type 3' conserved sequence (3'-CS), named Tn6217. Other features include three resistance genes similar to those of Tn5, a tellurite resistance operon, and two pil operons. The replication and maintenance systems exhibit similarity to a genomic island of Ralstonia solanacearum GM1000. Codon usage analysis suggests the recent acquisition of bla(IMP-9). The origins of the integrons on pOZ176 indicated separate horizontal gene transfer events driven by antibiotic selection. The novel mosaic structure of pOZ176 suggests that it is derived from environmental bacteria.