LMAS: evaluating metagenomic short de novo assembly methods through defined communities

BACKGROUND

The de novo assembly of raw sequence data is key in metagenomic analysis. It allows recovering draft genomes from a pool of mixed raw reads, yielding longer sequences that offer contextual information and provide a more complete picture of the microbial community.

FINDINGS

To better compare de novo assemblers for metagenomic analysis, LMAS (Last Metagenomic Assembler Standing) was developed as a flexible platform allowing users to evaluate assembler performance given known standard communities. Overall, in our test datasets, k-mer De Bruijn graph assemblers outperformed the alternative approaches but came with a greater computational cost. Furthermore, assemblers branded as metagenomic specific did not consistently outperform other genomic assemblers in metagenomic samples. Some assemblers still in use, such as ABySS, MetaHipmer2, minia, and VelvetOptimiser, perform relatively poorly and should be used with caution when assembling complex samples. Meaningful strain resolution at the single-nucleotide polymorphism level was not achieved, even by the best assemblers tested.

CONCLUSIONS

The choice of a de novo assembler depends on the computational resources available, the replicon of interest, and the major goals of the analysis. No single assembler appeared an ideal choice for short-read metagenomic prokaryote replicon assembly, each showing specific strengths. The choice of metagenomic assembler should be guided by user requirements and characteristics of the sample of interest, and LMAS provides an interactive evaluation platform for this purpose. LMAS is open source, and the workflow and its documentation are available at https://github.com/B-UMMI/LMAS and https://lmas.readthedocs.io/, respectively.

Software testing in microbial bioinformatics: a call to action

Computational algorithms have become an essential component of research, with great efforts by the scientific community to raise standards on development and distribution of code. Despite these efforts, sustainability and reproducibility are major issues since continued validation through software testing is still not a widely adopted practice. Here, we report seven recommendations that help researchers implement software testing in microbial bioinformatics. We have developed these recommendations based on our experience from a collaborative hackathon organised prior to the American Society for Microbiology Next Generation Sequencing (ASM NGS) 2020 conference. We also present a repository hosting examples and guidelines for testing, available from https://github.com/microbinfie-hackathon2020/CSIS.

Heterogeneity of penicillin-non-susceptible group B streptococci isolated from a single patient in Germany

OBJECTIVES

Streptococcus agalactiae [group B streptococci (GBS)] have been considered uniformly susceptible to penicillin. However, increasing reports from Asia and North America are documenting penicillin-non-susceptible GBS (PRGBS) with mutations in pbp genes. Here we report, to the best of our knowledge, the first two PRGBS isolates recovered in Europe (AC-13238-1 and AC-13238-2), isolated from the same patient.

METHODS

Two different colony morphologies of GBS were noted from a surgical abscess drainage sample. Both were serotyped and antimicrobial susceptibility testing was performed by different methodologies. High-throughput sequencing was done to compare the isolates at the genomic level, to identify their capsular type and ST, to evaluate mutations in the pbp genes and to compare the isolates with the genomes of other PRGBS isolates sharing the same serotype and ST.

RESULTS

Isolates AC-13238-1 and AC-13238-2 presented MICs above the EUCAST and CLSI breakpoints for penicillin susceptibility. Both shared the capsular type Ia operon and ST23. Genomic analysis uncovered differences between the two isolates in seven genes, including altered pbp genes. Deduced amino acid sequences revealed critical substitutions in PBP2X in both isolates. Comparison with serotype Ia clonal complex 23 PRGBS from the USA reinforced the similarity between AC-13238-1 and AC-13238-2, and their divergence from the US strains.

CONCLUSIONS

Our results support the in-host evolution of β-lactam-resistant GBS, with two PRGBS variants being isolated from one patient.

Chewie Nomenclature Server (chewie-NS): a deployable nomenclature server for easy sharing of core and whole genome MLST schemas

Chewie Nomenclature Server (chewie-NS, https://chewbbaca.online/) allows users to share genome-based gene-by-gene typing schemas and to maintain a common nomenclature, simplifying the comparison of results. The combination between local analyses and a public repository of allelic data strikes a balance between potential confidentiality issues and the need to compare results. The possibility of deploying private instances of chewie-NS facilitates the creation of nomenclature servers with a restricted user base to allow compliance with the strictest data policies. Chewie-NS allows users to easily share their own schemas and to explore publicly available schemas, including informative statistics on schemas and loci presented in interactive charts and tables. Users can retrieve all the information necessary to run a schema locally or all the alleles identified at a particular locus. The integration with the chewBBACA suite enables users to directly upload new schemas to chewie-NS, download existing schemas and synchronize local and remote schemas from chewBBACA command line version, allowing an easier integration into high-throughput analysis pipelines. The same REST API linking chewie-NS and the chewBBACA suite supports the interaction of other interfaces or pipelines with the databases available at chewie-NS, facilitating the reusability of the stored data.