Mash Screen: high-throughput sequence containment estimation for genome discovery

A survey of k-mer methods and applications in bioinformatics

The rapid progression of genomics and proteomics has been driven by the advent of advanced sequencing technologies, large, diverse, and readily available omics datasets, and the evolution of computational data processing capabilities. The vast amount of data generated by these advancements necessitates efficient algorithms to extract meaningful information. K-mers serve as a valuable tool when working with large sequencing datasets, offering several advantages in computational speed and memory efficiency and carrying the potential for intrinsic biological functionality. This review provides an overview of the methods, applications, and significance of k-mers in genomic and proteomic data analyses, as well as the utility of absent sequences, including nullomers and nullpeptides, in disease detection, vaccine development, therapeutics, and forensic science. Therefore, the review highlights the pivotal role of k-mers in addressing current genomic and proteomic problems and underscores their potential for future breakthroughs in research.

Metagenomic functional profiling: to sketch or not to sketch?

MOTIVATION

Functional profiling of metagenomic samples is essential to decipher the functional capabilities of microbial communities. Traditional and more widely used functional profilers in the context of metagenomics rely on aligning reads against a known reference database. However, aligning sequencing reads against a large and fast-growing database is computationally expensive. In general, k-mer-based sketching techniques have been successfully used in metagenomics to address this bottleneck, notably in taxonomic profiling. In this work, we describe leveraging FracMinHash (implemented in sourmash, a publicly available software), a k-mer-sketching algorithm, to obtain functional profiles of metagenome samples.

RESULTS

We show how pieces of the sourmash software (and the resulting FracMinHash sketches) can be put together in a pipeline to functionally profile a metagenomic sample. We named our pipeline fmh-funprofiler. We report that the functional profiles obtained using this pipeline demonstrate comparable completeness and better purity compared to the profiles obtained using other alignment-based methods when applied to simulated metagenomic data. We also report that fmh-funprofiler is 39-99× faster in wall-clock time, and consumes up to 40-55× less memory. Coupled with the KEGG database, this method not only replicates fundamental biological insights but also highlights novel signals from the Human Microbiome Project datasets.

AVAILABILITY AND IMPLEMENTATION

This fast and lightweight metagenomic functional profiler is freely available and can be accessed here: https://github.com/KoslickiLab/fmh-funprofiler. All scripts of the analyses we present in this manuscript can be found on GitHub.

Genome sequence of Bacillus pumilus RI06-95 isolated during a Microcystis bloom in lake Champlain, USA

Here, we present the complete genome sequence of Bacillus pumilus RI06-95 isolated during a Microcystis bloom in Lake Champlain. The assembled genome comprises a 3.8 Mbp chromosome with a GC content of 42%, and two plasmids, pPZZ84 6.4 Kbp, GC content 37% and pSHB9 97Kbp, GC content 36%.

CRISPR-AMRtracker: A novel toolkit to monitor the antimicrobial resistance gene transfer in fecal microbiota

The spread of antibiotic resistance genes (ARGs), particularly those carried on plasmids, poses a major risk to global health. However, the extent and frequency of ARGs transfer in microbial communities among human, animal, and environmental sectors is not well understood due to a lack of effective tracking tools. We have developed a novel fluorescent tracing tool, CRISPR-AMRtracker, to study ARG transfer. It combines CRISPR/Cas9 fluorescence tagging, fluorescence-activated cell sorting, 16S rRNA gene sequencing, and microbial community analysis. CRISPR-AMRtracker integrates a fluorescent tag immediately downstream of ARGs, enabling the tracking of ARG transfer without compromising the host cell's antibiotic susceptibility, fitness, conjugation, and transposition. Notably, our experiments demonstrate that sfGFP-tagged plasmid-borne mcr-1 can transfer across diverse bacterial species within fecal samples. This innovative approach holds the potential to illuminate the dynamics of ARG dissemination and provide valuable insights to shape effective strategies in mitigating the escalating threat of antibiotic resistance.

The integrated genomic surveillance system of Andalusia (SIEGA) provides a One Health regional resource connected with the clinic

The One Health approach, recognizing the interconnectedness of human, animal, and environmental health, has gained significance amid emerging zoonotic diseases and antibiotic resistance concerns. This paper aims to demonstrate the utility of a collaborative tool, the SIEGA, for monitoring infectious diseases across domains, fostering a comprehensive understanding of disease dynamics and risk factors, highlighting the pivotal role of One Health surveillance systems. Raw whole-genome sequencing is processed through different species-specific open software that additionally reports the presence of genes associated to anti-microbial resistances and virulence. The SIEGA application is a Laboratory Information Management System, that allows customizing reports, detect transmission chains, and promptly alert on alarming genetic similarities. The SIEGA initiative has successfully accumulated a comprehensive collection of more than 1900 bacterial genomes, including Salmonella enterica, Listeria monocytogenes, Campylobacter jejuni, Escherichia coli, Yersinia enterocolitica and Legionella pneumophila, showcasing its potential in monitoring pathogen transmission, resistance patterns, and virulence factors. SIEGA enables customizable reports and prompt detection of transmission chains, highlighting its contribution to enhancing vigilance and response capabilities. Here we show the potential of genomics in One Health surveillance when supported by an appropriate bioinformatic tool. By facilitating precise disease control strategies and antimicrobial resistance management, SIEGA enhances global health security and reduces the burden of infectious diseases. The integration of health data from humans, animals, and the environment, coupled with advanced genomics, underscores the importance of a holistic One Health approach in mitigating health threats.

ModDotPlot-rapid and interactive visualization of tandem repeats

MOTIVATION

A common method for analyzing genomic repeats is to produce a sequence similarity matrix visualized via a dot plot. Innovative approaches such as StainedGlass have improved upon this classic visualization by rendering dot plots as a heatmap of sequence identity, enabling researchers to better visualize multi-megabase tandem repeat arrays within centromeres and other heterochromatic regions of the genome. However, computing the similarity estimates for heatmaps requires high computational overhead and can suffer from decreasing accuracy.

RESULTS

In this work, we introduce ModDotPlot, an interactive and alignment-free dot plot viewer. By approximating average nucleotide identity via a k-mer-based containment index, ModDotPlot produces accurate plots orders of magnitude faster than StainedGlass. We accomplish this through the use of a hierarchical modimizer scheme that can visualize the full 128 Mb genome of Arabidopsis thaliana in under 5 min on a laptop. ModDotPlot is bundled with a graphical user interface supporting real-time interactive navigation of entire chromosomes.

AVAILABILITY AND IMPLEMENTATION

ModDotPlot is available at https://github.com/marbl/ModDotPlot.

Comprehensive blueprint of Salmonella genomic plasticity identifies hotspots for pathogenicity genes

Understanding the dynamic evolution of Salmonella is vital for effective bacterial infection management. This study explores the role of the flexible genome, organised in regions of genomic plasticity (RGP), in shaping the pathogenicity of Salmonella lineages. Through comprehensive genomic analysis of 12,244 Salmonella spp. genomes covering 2 species, 6 subspecies, and 46 serovars, we uncover distinct integration patterns of pathogenicity-related gene clusters into RGP, challenging traditional views of gene distribution. These RGP exhibit distinct preferences for specific genomic spots, and the presence or absence of such spots across Salmonella lineages profoundly shapes strain pathogenicity. RGP preferences are guided by conserved flanking genes surrounding integration spots, implicating their involvement in regulatory networks and functional synergies with integrated gene clusters. Additionally, we emphasise the multifaceted contributions of plasmids and prophages to the pathogenicity of diverse Salmonella lineages. Overall, this study provides a comprehensive blueprint of the pathogenicity potential of Salmonella. This unique insight identifies genomic spots in nonpathogenic lineages that hold the potential for harbouring pathogenicity genes, providing a foundation for predicting future adaptations and developing targeted strategies against emerging human pathogenic strains.

Fluoroquinolone-resistant Campylobacter in backyard and commercial broiler production systems in the United States

Objectives

Campylobacter spp. are one of the leading foodborne pathogens in the world, and chickens are a known reservoir. This is significant considering broiler chicken is the top consumed meat worldwide. In the USA, backyard poultry production is increasing, but little research has been done to investigate prevalence and antimicrobial resistance associated with Campylobacter in these environments.

Methods

Our study encompasses a farm-to-genome approach to identify Campylobacter and investigate its antimicrobial resistance phenotypically and genotypically. We travelled to 10 backyard and 10 integrated commercial broiler farms to follow a flock throughout production. We sampled at days 10, 31 and 52 for backyard and 10, 24 and 38 for commercial farms. Bird faecal (n = 10) and various environmental samples (soil n = 5, litter/compost n = 5, and feeder and waterer swabs n = 6) were collected at each visit and processed for Campylobacter.

Results

Our results show a higher prevalence of Campylobacter in samples from backyard farms (21.9%) compared to commercial (12.2%). Most of our isolates were identified as C. jejuni (70.8%) and the remainder as C. coli (29.2%). Antimicrobial susceptibility testing reveals phenotypic resistance to ciprofloxacin (40.2%), an important treatment drug for Campylobacter infection, and tetracycline (46.6%). A higher proportion of resistance was found in C. jejuni isolates and commercial farms. Whole-genome sequencing revealed resistance genes, such as tet(O) and gyrA_T86I point mutation, that may confer resistance.

Conclusion

Overall, our research emphasizes the need for interventions to curb prevalence of resistant Campylobacter spp. on broiler production systems.

Genome assembly, Full-length transcriptome, and isoform diversity of Red Snapper, Lutjanus argentimaculatus

The mangrove red snapper, Lutjanus argentimaculatus, is a marine food fish of economic and aquaculture importance. The application of genomic selection-based breeding programs for this species is limited by the absence of a reference genome and transcriptome profiles. The current study attempted to fill this void by generating genomic and transcriptomic resources for red snapper. Using PacBio long reads, and Arima Hi-C linked reads, a scaffold-level genome assembly was generated for L. argentimaculatus. The assembly is of 1.03 Gb comprising of 400 scaffolds with N50 of 33.8 Mb and was assessed to be 97.2% complete upon benchmarking with BUSCO. Full-length transcriptome generated with PacBio Iso-Sequencing strategy using six tissues (muscle, gills, liver, kidney, stomach, and gonad) contained 56,515 isoforms belonging to 18,108 unique genes with N50 length of 3,973 bp. The resources generated will have potential applications in the functional studies, conservation, broodstock management and selective breeding programmes of L. argentimaculatus.

HyperGen: Compact and Efficient Genome Sketching using Hyperdimensional Vectors

MOTIVATION

Genomic distance estimation is a critical workload since exact computation for whole-genome similarity metrics such as Average Nucleotide Identity (ANI) incurs prohibitive runtime overhead. Genome sketching is a fast and memory-efficient solution to estimate ANI similarity by distilling representative k-mers from the original sequences. In this work, we present HyperGen that improves accuracy, runtime performance, and memory efficiency for large-scale ANI estimation. Unlike existing genome sketching algorithms that convert large genome files into discrete k-mer hashes, HyperGen leverages the emerging hyperdimensional computing (HDC) to encode genomes into quasi-orthogonal vectors (Hypervector, HV) in high-dimensional space. HV is compact and can preserve more information, allowing for accurate ANI estimation while reducing required sketch sizes. In particular, the HV sketch representation in HyperGen allows efficient ANI estimation using vector multiplication, which naturally benefits from highly optimized general matrix multiply (GEMM) routines. As a result, HyperGen enables the efficient sketching and ANI estimation for massive genome collections.

RESULTS

We evaluate HyperGen 's sketching and database search performance using several genome datasets at various scales. HyperGen is able to achieve comparable or superior ANI estimation error and linearity compared to other sketch-based counterparts. The measurement results show that HyperGen is one of the fastest tools for both genome sketching and database search. Meanwhile, HyperGen produces memory-efficient sketch files while ensuring high ANI estimation accuracy.

AVAILABILITY

A Rust implementation of HyperGen is freely available under the MIT license as an open-source software project at https://github.com/wh-xu/Hyper-Gen. The scripts to reproduce the experimental results can be accessed at https://github.com/wh-xu/experiment-hyper-gen.

Comparative genomics of quinolone-resistant Escherichia coli from broilers and humans in Norway

BACKGROUND

The usage of fluoroquinolones in Norwegian livestock production is very low, including in broiler production. Historically, quinolone-resistant Escherichia coli (QREC) isolated from Norwegian production animals rarely occur. However, with the introduction of a selective screening method for QREC in the Norwegian monitoring programme for antimicrobial resistance in the veterinary sector in 2014; 89.5% of broiler caecal samples and 70.7% of broiler meat samples were positive. This triggered the concern if there could be possible links between broiler and human reservoirs of QREC. We are addressing this by characterizing genomes of QREC from humans (healthy carriers and patients) and broiler isolates (meat and caecum).

RESULTS

The most frequent mechanism for quinolone resistance in both broiler and human E. coli isolates were mutations in the chromosomally located gyrA and parC genes, although plasmid mediated quinolone resistance (PMQR) was also identified. There was some relatedness of the isolates within human and broiler groups, but little between these two groups. Further, some overlap was seen for isolates with the same sequence type isolated from broiler and humans, but overall, the SNP distance was high.

CONCLUSION

Based on data from this study, QREC from broiler makes a limited contribution to the incidence of QREC in humans in Norway.

Core genome multilocus sequence typing (cgMLST) applicable to the monophyletic Klebsiella oxytoca species complex

The environmental bacterium Klebsiella oxytoca displays an alarming increase of antibiotic-resistant strains that frequently cause outbreaks in intensive care units. Due to its prevalence in the environment and opportunistic presence in humans, molecular surveillance (including resistance marker screening) and high-resolution cluster analysis are of high relevance. Furthermore, K. oxytoca previously described in studies is rather a species complex (KoSC) than a single species comprising at least six closely related species that are not easily differentiated by standard typing methods. To reach a discriminatory power high enough to identify and resolve clusters within these species, whole genome sequencing is necessary. The resolution is achievable with core genome multilocus sequence typing (cgMLST) extending typing of a few housekeeping genes to thousands of core genome genes. CgMLST is highly standardized and provides a nomenclature enabling cross laboratory reproducibility and data exchange for routine diagnostics. Here, we established a cgMLST scheme not only capable of resolving the KoSC species but also producing reliable and consistent results for published outbreaks. Our cgMLST scheme consists of 2,536 core genome and 2,693 accessory genome targets, with a percentage of good cgMLST targets of 98.31% in 880 KoSC genomes downloaded from the National Center for Biotechnology Information (NCBI). We also validated resistance markers against known resistance gene patterns and successfully linked genetic results to phenotypically confirmed toxic strains carrying the til gene cluster. In conclusion, our novel cgMLST enables highly reproducible typing of four different clinically relevant species of the KoSC and thus facilitates molecular surveillance and cluster investigations.

Exploration of mobile genetic elements in the ruminal microbiome of Nellore cattle

Metagenomics has made it feasible to elucidate the intricacies of the ruminal microbiome and its role in the differentiation of animal production phenotypes of significance. The search for mobile genetic elements (MGEs) has taken on great importance, as they play a critical role in the transfer of genetic material between organisms. Furthermore, these elements serve a dual purpose by controlling populations through lytic bacteriophages, thereby maintaining ecological equilibrium and driving the evolutionary progress of host microorganisms. In this study, we aimed to identify the association between ruminal bacteria and their MGEs in Nellore cattle using physical chromosomal links through the Hi-C method. Shotgun metagenomic sequencing and the proximity ligation method ProxiMeta were used to analyze DNA, getting 1,713,111,307 bp, which gave rise to 107 metagenome-assembled genomes from rumen samples of four Nellore cows maintained on pasture. Taxonomic analysis revealed that most of the bacterial genomes belonged to the families Lachnospiraceae, Bacteroidaceae, Ruminococcaceae, Saccharofermentanaceae, and Treponemataceae and mostly encoded pathways for central carbon and other carbohydrate metabolisms. A total of 31 associations between host bacteria and MGE were identified, including 17 links to viruses and 14 links to plasmids. Additionally, we found 12 antibiotic resistance genes. To our knowledge, this is the first study in Brazilian cattle that connect MGEs with their microbial hosts. It identifies MGEs present in the rumen of pasture-raised Nellore cattle, offering insights that could advance biotechnology for food digestion and improve ruminant performance in production systems.

ModDotPlot-Rapid and interactive visualization of complex repeats

Motivation

A common method for analyzing genomic repeats is to produce a sequence similarity matrix visualized via a dot plot. Innovative approaches such as StainedGlass have improved upon this classic visualization by rendering dot plots as a heatmap of sequence identity, enabling researchers to better visualize multi-megabase tandem repeat arrays within centromeres and other heterochromatic regions of the genome. However, computing the similarity estimates for heatmaps requires high computational overhead and can suffer from decreasing accuracy.

Results

In this work we introduce ModDotPlot, an interactive and alignment-free dot plot viewer. By approximating average nucleotide identity via a k-mer-based containment index, ModDotPlot produces accurate plots orders of magnitude faster than StainedGlass. We accomplish this through the use of a hierarchical modimizer scheme that can visualize the full 128 Mbp genome of Arabidopsis thaliana in under 5 minutes on a laptop. ModDotPlot is bundled with a graphical user interface supporting real-time interactive navigation of entire chromosomes.

Availability and Implementation

ModDotPlot is available at https://github.com/marbl/ModDotPlot.

Metagenomics untangles potential adaptations of Antarctic endolithic bacteria at the fringe of habitability

Survival and growth strategies of Antarctic endolithic microbes residing in Earth's driest and coldest desert remain virtually unknown. From 109 endolithic microbiomes, 4539 metagenome-assembled genomes were generated, 49.3 % of which were novel candidate bacterial species. We present evidence that trace gas oxidation and atmospheric chemosynthesis may be the prevalent strategies supporting metabolic activity and persistence of these ecosystems at the fringe of life and the limits of habitability.

Prolonged carriage of OXA-244-carbapenemase-producing Escherichia coli complicates epidemiological investigations

The rapid increase of OXA-244-producing Escherichia coli, predominantly driven by genetically clustered isolates of sequence type (ST)38, has been observed in at least nine European countries, including Germany. However, the reasons for the spread of OXA-244-producing E. coli remain unclear. Here, we aim to evaluate the possibility of prolonged carriage. We identified a total of six different patients with repeated detection of OXA-244-producing E. coli isolates, which were subjected to both short and long-read whole-genome sequencing (WGS). Besides allelic differences using core genome multilocus sequence typing (cgMLST) analyses, we obtained numbers of single-nucleotide polymorphisms (SNPs) to calculate individual base-pair substitution (BPS) rates. To assess possible re-exposure and risk factors for prolonged carriage, case interviews were conducted. The time between detections ranged from eleven months to more than three years. Initial isolates originated in three+ out of six cases from clinical samples, whereas remaining samples were from screening, mostly in the inpatient setting. As expected, cgMLST analyses showed low numbers of allelic differences between isolates of each case ranging from 1 to 4, whereas numbers of SNPs were between 2 and 99 (mean = 36), thus clearly highlighting the discrepancy between these different bacterial typing approaches. For five out of six cases, observed BPS rates suggest that patients can be colonized with OXA-244-producing E. coli, including ST38 cluster isolates, for extensively long times. Thus, we may have previously missed the epidemiological link between cases because exposure to OXA-244-producing E. coli could have occurred in a time frame, which has not been evaluated in previous investigations. Our results may help to guide future epidemiological investigations as well as to support the interpretation of genetic diversity of OXA-244-producing E. coli, particularly among ST38 cluster isolates.

A Chromosome-Level Genome Assembly and Annotation for the Clouded Apollo Butterfly (Parnassius mnemosyne): A Species of Global Conservation Concern

The clouded apollo (Parnassius mnemosyne) is a palearctic butterfly distributed over a large part of western Eurasia, but population declines and fragmentation have been observed in many parts of the range. The development of genomic tools can help to shed light on the genetic consequences of the decline and to make informed decisions about direct conservation actions. Here, we present a high-contiguity, chromosome-level genome assembly of a female clouded apollo butterfly and provide detailed annotations of genes and transposable elements. We find that the large genome (1.5 Gb) of the clouded apollo is extraordinarily repeat rich (73%). Despite that, the combination of sequencing techniques allowed us to assemble all chromosomes (nc = 29) to a high degree of completeness. The annotation resulted in a relatively high number of protein-coding genes (22,854) compared with other Lepidoptera, of which a large proportion (21,635) could be assigned functions based on homology with other species. A comparative analysis indicates that overall genome structure has been largely conserved, both within the genus and compared with the ancestral lepidopteran karyotype. The high-quality genome assembly and detailed annotation presented here will constitute an important tool for forthcoming efforts aimed at understanding the genetic consequences of fragmentation and decline, as well as for assessments of genetic diversity, population structure, inbreeding, and genetic load in the clouded apollo butterfly.

ICEberg 3.0: functional categorization and analysis of the integrative and conjugative elements in bacteria

ICEberg 3.0 (https://tool2-mml.sjtu.edu.cn/ICEberg3/) is an upgraded database that provides comprehensive insights into bacterial integrative and conjugative elements (ICEs). In comparison to the previous version, three key enhancements were introduced: First, through text mining and manual curation, it now encompasses details of 2065 ICEs, 607 IMEs and 275 CIMEs, including 430 with experimental support. Secondly, ICEberg 3.0 systematically categorizes cargo gene functions of ICEs into six groups based on literature curation and predictive analysis, providing a profound understanding of ICEs'diverse biological traits. The cargo gene prediction pipeline is integrated into the online tool ICEfinder 2.0. Finally, ICEberg 3.0 aids the analysis and exploration of ICEs from the human microbiome. Extracted and manually curated from 2405 distinct human microbiome samples, the database comprises 1386 putative ICEs, offering insights into the complex dynamics of Bacteria-ICE-Cargo networks within the human microbiome. With the recent updates, ICEberg 3.0 enhances its capability to unravel the intricacies of ICE biology, particularly in the characterization and understanding of cargo gene functions and ICE interactions within the microbiome. This enhancement may facilitate the investigation of the dynamic landscape of ICE biology and its implications for microbial communities.

A two-sequence motif-based method for the inventory of gene families in fragmented and poorly annotated genome sequences

Databases of genome sequences are growing exponentially, but, in some cases, assembly is incomplete and genes are poorly annotated. For evolutionary studies, it is important to identify all members of a given gene family in a genome. We developed a method for identifying most, if not all, members of a gene family from raw genomes in which assembly is of low quality, using the P-type ATPase superfamily as an example. The method is based on the translation of an entire genome in all six reading frames and the co-occurrence of two family-specific sequence motifs that are in close proximity to each other. To test the method's usability, we first used it to identify P-type ATPase members in the high-quality annotated genome of barley (Hordeum vulgare). Subsequently, after successfully identifying plasma membrane H+-ATPase family members (P3A ATPases) in various plant genomes of varying quality, we tested the hypothesis that the number of P3A ATPases correlates with the ability of the plant to tolerate saline conditions. In 19 genomes of glycophytes and halophytes, the total number of P3A ATPase genes was found to vary from 7 to 22, but no significant difference was found between the two groups. The method successfully identified P-type ATPase family members in raw genomes that are poorly assembled.

Novel strains of Culex flavivirus and Hubei chryso-like virus 1 from the Anopheles mosquito in western Kenya

Surveillance of mosquito vectors is critical for early detection, prevention and control of vector borne diseases. In this study we used advanced molecular tools, such as DNA barcoding in combination with novel sequencing technologies to discover new and already known viruses in genetically identified mosquito species. Mosquitoes were captured using BG sentinel traps in Western Kenya during May and July 2019, and homogenized individually before pooled into groups of ten mosquitoes. The pools and individual samples were then used for molecular analysis and to infect cell cultures. Of a total of fifty-four (54) 10-pools, thirteen (13) showed cytopathic effect (CPE) on VeroB4 cells, eighteen (18) showed CPE on C6/36 cells. Eight (8) 10-pools out of the 31 CPE positive pools showed CPE on both VeroB4 and C6/36 cells. When using reverse transcriptase polymerase chain reaction (RT-PCR), Sanger sequencing and Twist Comprehensive Viral Research Panel (CVRP) (Twist Biosciences), all pools were found negative by RT-PCR when using genus specific primers targeting alphaviruses, orthobunyaviruses and virus specific primers towards o'nyong-nyong virus, chikungunya virus and Sindbis virus (previously reported to circulate in the region). Interestingly, five pools were RT-PCR positive for flavivirus. Two of the RT-PCR positive pools showed CPE on both VeroB4 and C6/36 cells, two pools showed CPE on C6/36 cells alone and one pool on VeroB4 cells only. Fifty individual mosquito homogenates from the five RT-PCR positive 10-pools were analyzed further for flavivirus RNA. Of these, 19 out of the 50 individual mosquito homogenates indicated the presence of flavivirus RNA. Barcoding of the flavivirus positive mosquitoes revealed the mosquito species as Aedes aegypti (1), Mansonia uniformis (6), Anopheles spp (3), Culex pipiens (5), Culex spp (1), Coquilletidia metallica (2) and Culex quinquefasciatus (1). Of the 19 flavivirus positive individual mosquitoes, five (5) virus positive homogenates were sequenced. Genome sequences of two viruses were completed. One was identified as the single-stranded RNA Culex flavivirus and the other as the double-stranded RNA Hubei chryso-like virus 1. Both viruses were found in the same Anopheles spp. homogenate extracted from a sample that showed CPE on both VeroB4 and C6/36 cells. The detection of both viruses in a single mosquito homogenate indicated coinfection. Phylogenetic analyses suggested that the Culex flavivirus sequence detected was closely related to a Culex flavivirus isolated from Uganda in 2008. All four Hubei chryso-like virus 1 segments clusters closely to Hubei chryso-like virus 1 strains isolated in Australia, China and USA. Two novel strains of insect-specific viruses in Anopheles mosquitoes were detected and characterized.

The protective role of commensal gut microbes and their metabolites against bacterial pathogens

Multidrug-resistant microorganisms have become a major public health concern around the world. The gut microbiome is a gold mine for bioactive compounds that protect the human body from pathogens. We used a multi-omics approach that integrated whole-genome sequencing (WGS) of 74 commensal gut microbiome isolates with metabolome analysis to discover their metabolic interaction with Salmonella and other antibiotic-resistant pathogens. We evaluated differences in the functional potential of these selected isolates based on WGS annotation profiles. Furthermore, the top altered metabolites in co-culture supernatants of selected commensal gut microbiome isolates were identified including a series of dipeptides and examined for their ability to prevent the growth of various antibiotic-resistant bacteria. Our results provide compelling evidence that the gut microbiome produces metabolites, including the compound class of dipeptides that can potentially be applied for anti-infection medication, especially against antibiotic-resistant pathogens. Our established pipeline for the discovery and validation of bioactive metabolites from the gut microbiome as novel candidates for multidrug-resistant infections represents a new avenue for the discovery of antimicrobial lead structures.

Sketched reference databases for genome-based taxonomy and comparative genomics

Abstract The analysis of curated genomic, metagenomic and proteomic data is of paramount importance in the fields of biology, medicine, education, and bioinformatics. Although this type of data is usually hosted in raw format on free international repositories, the full access requires lots of computing power and large storage disk space for the domestic user. The purpose of the study is to offer a comprehensive set of microbial genomic and proteomic reference databases in an accessible and easy-to-use form to the scientific community and demonstrate its advantages and usefulness. Also, we present a case study on the applicability of the sketched data, for the determination of overall genomic coherence between two members of the Brucellacea family, which suggests they belong to the same genomospecies that remain as discrete ecotypes. A representative set of genomes, proteomes (from type material), and metagenomes were directly collected from the NCBI Assembly database and Genome Taxonomy Database (GTDB), associated with the major groups of Bacteria, Archaea, Virus, and Fungi. Sketched databases were subsequently created and stored on handy reduced representations by using the MinHash algorithm implemented in Mash software. The obtained dataset contains more than 133 GB of space disk reduced to 883.25 MB and represents 125,110 genomics/proteomic records from eight informative contexts, which have been prefiltered to make them accessible, usable, and user-friendly with limited computational resources. Potential uses of these sketched databases are discussed, including but not limited to microbial species delimitation, estimation of genomic distances and genomic novelties, paired comparisons between proteomes, genomes, and metagenomes; phylogenetic neighbor’s exploration and selection, among others.

Molecular and Pathological Characterization of Classical Swine Fever Virus Genotype 2 Strains Responsible for the 2013-2018 Outbreak in Colombia

Classical swine fever (CSF) is a highly contagious transboundary viral disease of domestic and wild pigs. Despite mass vaccination and continuous eradication programs, CSF remains endemic in Asia, some countries in Europe, the Caribbean and South America. Since June 2013, Northern Colombia has reported 137 CSF outbreaks, mostly in backyard production systems with low vaccination coverage. The purpose of this study was to characterize the virus responsible for the outbreak. Phylogenetic analysis based on the full-length E2 sequence shows that the virus is closely related to CSF virus (CSFV) genotype 2.6 strains circulating in Southeast Asia. The pathotyping experiment suggests that the virus responsible is a moderately virulent strain. The 190 nucleotide stretch of the E2 hypervariable region of these isolates also shows high similarity to the CSFV isolates from Colombia in 2005 and 2006, suggesting a common origin for the CSF outbreaks caused by genotype 2.6 strains. The emergence of genotype 2.6 in Colombia suggests a potential transboundary spread of CSFV from Asia to the Americas, complicating the ongoing CSF eradication efforts in the Americas, and emphasizes the need for continuous surveillance in the region.

Evolutionary Insight into the Association between New Jersey Polyomavirus and Humans

Advances in viral discovery techniques have led to the identification of numerous novel viruses in human samples. However, the low prevalence of certain viruses in humans raises doubts about their association with our species. To ascertain the authenticity of a virus as a genuine human-infecting agent, it can be useful to investigate the diversification of its lineage within hominines, the group encompassing humans and African great apes. Building upon this rationale, we examined the case of the New Jersey polyomavirus (NJPyV; Alphapolyomavirus terdecihominis), which has only been detected in a single patient thus far. In this study, we obtained and analyzed sequences from closely related viruses infecting all African great ape species. We show that NJPyV nests within the diversity of these viruses and that its lineage placement is compatible with an ancient origin in humans, despite its apparent rarity in human populations.

A genomic catalogue of soil microbiomes boosts mining of biodiversity and genetic resources

Soil harbors a vast expanse of unidentified microbes, termed as microbial dark matter, presenting an untapped reservo)ir of microbial biodiversity and genetic resources, but has yet to be fully explored. In this study, we conduct a large-scale excavation of soil microbial dark matter by reconstructing 40,039 metagenome-assembled genome bins (the SMAG catalogue) from 3304 soil metagenomes. We identify 16,530 of 21,077 species-level genome bins (SGBs) as unknown SGBs (uSGBs), which expand archaeal and bacterial diversity across the tree of life. We also illustrate the pivotal role of uSGBs in augmenting soil microbiome's functional landscape and intra-species genome diversity, providing large proportions of the 43,169 biosynthetic gene clusters and 8545 CRISPR-Cas genes. Additionally, we determine that uSGBs contributed 84.6% of previously unexplored viral-host associations from the SMAG catalogue. The SMAG catalogue provides an useful genomic resource for further studies investigating soil microbial biodiversity and genetic resources.

Phosphorylation of Human Polyomavirus Large and Small T Antigens: An Ignored Research Field

Protein phosphorylation and dephosphorylation are the most common post-translational modifications mediated by protein kinases and protein phosphatases, respectively. These reversible processes can modulate the function of the target protein, such as its activity, subcellular localization, stability, and interaction with other proteins. Phosphorylation of viral proteins plays an important role in the life cycle of a virus. In this review, we highlight biological implications of the phosphorylation of the monkey polyomavirus SV40 large T and small t antigens, summarize our current knowledge of the phosphorylation of these proteins of human polyomaviruses, and conclude with gaps in the knowledge and a proposal for future research directions.

Outbreak of Pandoraea commovens Infections among Non-Cystic Fibrosis Intensive Care Patients, Germany, 2019-2021

Pandoraea spp. are gram-negative, nonfermenting rods mainly known to infect patients with cystic fibrosis (CF). Outbreaks have been reported from several CF centers. We report a Pandoraea spp. outbreak comprising 24 non-CF patients at a large university hospital and a neighboring heart center in Germany during July 2019-December 2021. Common features in the patients were critical illness, invasive ventilation, antimicrobial pretreatment, and preceding surgery. Complicated and relapsing clinical courses were observed in cases with intraabdominal infections but not those with lower respiratory tract infections. Genomic analysis of 15 isolates identified Pandoraea commovens as the genetically most similar species and confirmed the clonality of the outbreak strain, designated P. commovens strain LB-19-202-79. The strain exhibited resistance to most antimicrobial drugs except ampicillin/sulbactam, imipenem, and trimethoprim/sulfamethoxazole. Our findings suggest Pandoraea spp. can spread among non-CF patients and underscore that clinicians and microbiologists should be vigilant in detecting and assessing unusual pathogens.

Cyanobacteria as a critical reservoir of the environmental antimicrobial resistome

Antimicrobial resistance (AMR) is predicted to cause a worldwide annual toll of 10 million deaths by 2050. This looming public health threat has been linked to antibiotic overuse and pollution, which places selective pressures on AMR maintenance and transfer in and between microbial populations. We examined the distribution, diversity and potential mobility of AMR genes in cyanobacteria. While cyanobacteria are not pathogenic, we hypothesised that they could be a major environmental reservoir for AMR genes. Genes encoding AMR to seven antimicrobial drug classes were found in 10% of cyanobacterial genomes. AMR genes were found in 13% of freshwater, 19% of terrestrial, 34% of symbiotic, 2% of thermal spring, and 3% of marine genomes. AMR genes were found in five cyanobacterial orders with 23% of Nostocales and 8% of Oscillatoriales strains containing AMR genes. The most frequently observed alleles were ansamycin resistance genes, which were present in 7% of strains. AMR genes responsible for resistance to broad-spectrum β-lactams, chloramphenicols, tetracyclines, macrolides, and aminoglycosides were associated with mobile genetic elements or plasmid replicons or both. These results suggest that cyanobacteria are an extensive reservoir, and potential vector, for AMR genes in diverse terrestrial and aquatic habitats.

Nuclear exchange generates population diversity in the wheat leaf rust pathogen Puccinia triticina

In clonally reproducing dikaryotic rust fungi, non-sexual processes such as somatic nuclear exchange are postulated to play a role in diversity but have been difficult to detect due to the lack of genome resolution between the two haploid nuclei. We examined three nuclear-phased genome assemblies of Puccinia triticina, which causes wheat leaf rust disease. We found that the most recently emerged Australian lineage was derived by nuclear exchange between two pre-existing lineages, which originated in Europe and North America. Haplotype-specific phylogenetic analysis reveals that repeated somatic exchange events have shuffled haploid nuclei between long-term clonal lineages, leading to a global P. triticina population representing different combinations of a limited number of haploid genomes. Thus, nuclear exchange seems to be the predominant mechanism generating diversity and the emergence of new strains in this otherwise clonal pathogen. Such genomics-accelerated surveillance of pathogen evolution paves the way for more accurate global disease monitoring.

RabbitKSSD: accelerating genome distance estimation on modern multi-core architectures

SUMMARY

We propose RabbitKSSD, a high-speed genome distance estimation tool. Specifically, we leverage load-balanced task partitioning, fast I/O, efficient intermediate result accesses, and high-performance data structures to improve overall efficiency. Our performance evaluation demonstrates that RabbitKSSD achieves speedups ranging from 5.7× to 19.8× over Kssd for the time-consuming sketch generation and distance computation on commonly used workstations. In addition, it significantly outperforms Mash, BinDash, and Dashing2. Moreover, RabbitKSSD can efficiently perform all-vs-all distance computation for all RefSeq complete bacterial genomes (455 GB in FASTA format) in just 2 min on a 64-core workstation.

AVAILABILITY AND IMPLEMENTATION

RabbitKSSD is available at https://github.com/RabbitBio/RabbitKSSD.

Pangenome analysis reveals the genetic basis for taxonomic classification of the Lactobacillaceae family

Lactobacillaceae represent a large family of important microbes that are foundational to the food industry. Many genome sequences of Lactobacillaceae strains are now available, enabling us to conduct a comprehensive pangenome analysis of this family. We collected 3591 high-quality genomes from public sources and found that: 1) they contained enough genomes for 26 species to perform a pangenomic analysis, 2) the normalized Heap's coefficient λ (a measure of pangenome openness) was found to have an average value of 0.27 (ranging from 0.07 to 0.37), 3) the pangenome openness was correlated with the abundance and genomic location of transposons and mobilomes, 4) the pangenome for each species was divided into core, accessory, and rare genomes, that highlight the species-specific properties (such as motility and restriction-modification systems), 5) the pangenome of Lactiplantibacillus plantarum (which contained the highest number of genomes found amongst the 26 species studied) contained nine distinct phylogroups, and 6) genome mining revealed a richness of detected biosynthetic gene clusters, with functions ranging from antimicrobial and probiotic to food preservation, but ∼93% were of unknown function. This study provides the first in-depth comparative pangenomics analysis of the Lactobacillaceae family.

Reference genomes for BALB/c Nude and NOD/SCID mouse models

Mouse xenograft models play a vital role in tumor studies for research as well as for screening of drugs for the pharmaceutical industry. In particular, models with compromised immunity are favorable to increase the success of transplantation, such as, e.g. NOD/SCID and BALB/c Nude strains. The genomic sequence and alterations of many of these models still remain elusive and might hamper a model's further optimization or proper adapted usage. This can be in respect to treatments (e.g. NOD/SCID sensitivity to radiation), experiments or analysis of derived sequencing data of such models. Here we present the genome assemblies for the NOD/SCID and BALB/c Nude strains to overcome this short-coming for the future and improve our understanding of these models in the process. We highlight as well first insights into observed genomic differences for these models compared to the C57BL/6 reference genome. Genome assemblies for both are close to full-chromosome representations and provided with liftover annotations from the GRCm39 reference genome.

The phageome in normal and inflamed human skin

Dysbiosis of skin microbiota drives the progression of atopic dermatitis (AD). The contribution of bacteriophages to bacterial community compositions in normal and inflamed skin is unknown. Using shotgun metagenomics from skin swabs of healthy individuals and patients with AD, we found 13,586 potential viral contiguous DNA sequences, which could be combined into 164 putative viral genomes including 133 putative phages. The Shannon diversity index for the viral metagenome-assembled genomes (vMAGs) did not correlate with AD. In total, we identified 28 vMAGs that differed significantly between normal and AD skin. Quantitative polymerase chain reaction validation of three complete vMAGs revealed their independence from host bacterium abundance. Our data indicate that normal and inflamed skin harbor distinct phageomes and suggest a causative relationship between changing viral and bacterial communities as a driver of skin pathology.

Diverse Durham collection phages demonstrate complex BREX defense responses

Bacteriophages (phages) outnumber bacteria ten-to-one and cause infections at a rate of 1025 per second. The ability of phages to reduce bacterial populations makes them attractive alternative antibacterials for use in combating the rise in antimicrobial resistance. This effort may be hindered due to bacterial defenses such as Bacteriophage Exclusion (BREX) that have arisen from the constant evolutionary battle between bacteria and phages. For phages to be widely accepted as therapeutics in Western medicine, more must be understood about bacteria-phage interactions and the outcomes of bacterial phage defense. Here, we present the annotated genomes of 12 novel bacteriophage species isolated from water sources in Durham, UK, during undergraduate practical classes. The collection includes diverse species from across known phylogenetic groups. Comparative analyses of two novel phages from the collection suggest they may be founding members of a new genus. Using this Durham phage collection, we determined that particular BREX defense systems were likely to confer a varied degree of resistance against an invading phage. We concluded that the number of BREX target motifs encoded in the phage genome was not proportional to the degree of susceptibility. IMPORTANCE Bacteriophages have long been the source of tools for biotechnology that are in everyday use in molecular biology research laboratories worldwide. Phages make attractive new targets for the development of novel antimicrobials. While the number of phage genome depositions has increased in recent years, the expected bacteriophage diversity remains underrepresented. Here we demonstrate how undergraduates can contribute to the identification of novel phages and that a single City in England can provide ample phage diversity and the opportunity to find novel technologies. Moreover, we demonstrate that the interactions and intricacies of the interplay between bacterial phage defense systems such as Bacteriophage Exclusion (BREX) and phages are more complex than originally thought. Further work will be required in the field before the dynamic interactions between phages and bacterial defense systems are fully understood and integrated with novel phage therapies.

Minmers are a generalization of minimizers that enable unbiased local Jaccard estimation

MOTIVATION

The Jaccard similarity on k-mer sets has shown to be a convenient proxy for sequence identity. By avoiding expensive base-level alignments and comparing reduced sequence representations, tools such as MashMap can scale to massive numbers of pairwise comparisons while still providing useful similarity estimates. However, due to their reliance on minimizer winnowing, previous versions of MashMap were shown to be biased and inconsistent estimators of Jaccard similarity. This directly impacts downstream tools that rely on the accuracy of these estimates.

RESULTS

To address this, we propose the minmer winnowing scheme, which generalizes the minimizer scheme by use of a rolling minhash with multiple sampled k-mers per window. We show both theoretically and empirically that minmers yield an unbiased estimator of local Jaccard similarity, and we implement this scheme in an updated version of MashMap. The minmer-based implementation is over 10 times faster than the minimizer-based version under the default ANI threshold, making it well-suited for large-scale comparative genomics applications.

AVAILABILITY AND IMPLEMENTATION

MashMap3 is available at https://github.com/marbl/MashMap.

Intensive care unit sinks are persistently colonized with multidrug resistant bacteria and mobilizable, resistance-conferring plasmids

Contamination of hospital sinks with microbial pathogens presents a serious potential threat to patients, but our understanding of sink colonization dynamics is largely based on infection outbreaks. Here, we investigate the colonization patterns of multidrug-resistant organisms (MDROs) in intensive care unit sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. Using culture-based methods, we recovered 822 bacterial isolates representing 104 unique species and genomospecies. Genomic analyses revealed long-term colonization by Pseudomonas spp. and Serratia marcescens strains across multiple rooms. Nanopore sequencing uncovered examples of long-term persistence of resistance-conferring plasmids in unrelated hosts. These data indicate that antibiotic resistance (AR) in Pseudomonas spp. is maintained both by strain colonization and horizontal gene transfer (HGT), while HGT maintains AR within Acinetobacter spp. and Enterobacterales, independent of colonization. These results emphasize the importance of proactive, genomic-focused surveillance of built environments to mitigate MDRO spread. IMPORTANCE Hospital sinks are frequently linked to outbreaks of antibiotic-resistant bacteria. Here, we used whole-genome sequencing to track the long-term colonization patterns in intensive care unit (ICU) sinks and water from two hospitals in the USA and Pakistan collected over 27 months of prospective sampling. We analyzed 822 bacterial genomes, representing over 100 different species. We identified long-term contamination by opportunistic pathogens, as well as transient appearance of other common pathogens. We found that bacteria recovered from the ICU had more antibiotic resistance genes (ARGs) in their genomes compared to matched community spaces. We also found that many of these ARGs are harbored on mobilizable plasmids, which were found shared in the genomes of unrelated bacteria. Overall, this study provides an in-depth view of contamination patterns for common nosocomial pathogens and identifies specific targets for surveillance.

Widespread Detection of Yersiniabactin Gene Cluster and Its Encoding Integrative Conjugative Elements (ICEKp) among Nonoutbreak OXA-48-Producing Klebsiella pneumoniae Clinical Isolates from Spain and the Netherlands

In this study, we determined the presence of virulence factors in nonoutbreak, high-risk clones and other isolates belonging to less common sequence types associated with the spread of OXA-48-producing Klebsiella pneumoniae clinical isolates from The Netherlands (n = 61) and Spain (n = 53). Most isolates shared a chromosomally encoded core of virulence factors, including the enterobactin gene cluster, fimbrial fim and mrk gene clusters, and urea metabolism genes (ureAD). We observed a high diversity of K-Locus and K/O loci combinations, KL17 and KL24 (both 16%), and the O1/O2v1 locus (51%) being the most prevalent in our study. The most prevalent accessory virulence factor was the yersiniabactin gene cluster (66.7%). We found seven yersiniabactin lineages-ybt 9, ybt 10, ybt 13, ybt 14, ybt 16, ybt 17, and ybt 27-which were chromosomally embedded in seven integrative conjugative elements (ICEKp): ICEKp3, ICEKp4, ICEKp2, ICEKp5, ICEKp12, ICEKp10, and ICEKp22, respectively. Multidrug-resistant lineages-ST11, ST101, and ST405-were associated with ybt 10/ICEKp4, ybt 9/ICEKp3, and ybt 27/ICEKp22, respectively. The fimbrial adhesin kpi operon (kpiABCDEFG) was predominant among ST14, ST15, and ST405 isolates, as well as the ferric uptake system kfuABC, which was also predominant among ST101 isolates. No convergence of hypervirulence and resistance was observed in this collection of OXA-48-producing K. pneumoniae clinical isolates. Nevertheless, two isolates, ST133 and ST792, were positive for the genotoxin colibactin gene cluster (ICEKp10). In this study, the integrative conjugative element, ICEKp, was the major vehicle for yersiniabactin and colibactin gene clusters spreading. IMPORTANCE Convergence of multidrug resistance and hypervirulence in Klebsiella pneumoniae isolates has been reported mostly related to sporadic cases or small outbreaks. Nevertheless, little is known about the real prevalence of carbapenem-resistant hypervirulent K. pneumoniae since these two phenomena are often separately studied. In this study, we gathered information on the virulent content of nonoutbreak, high-risk clones (i.e., ST11, ST15, and ST405) and other less common STs associated with the spread of OXA-48-producing K. pneumoniae clinical isolates. The study of virulence content in nonoutbreak isolates can help us to expand information on the genomic landscape of virulence factors in K. pneumoniae population by identifying virulence markers and their mechanisms of spread. Surveillance should focus not only on antimicrobial resistance but also on virulence characteristics to avoid the spread of multidrug and (hyper)virulent K. pneumoniae that may cause untreatable and more severe infections.

Short- and Long-Read Sequencing Reveals the Presence and Evolution of an IncF Plasmid Harboring blaCTX-M-15 and blaCTX-M-27 Genes in Escherichia coli ST131

Escherichia coli sequence type 131 (ST131) has contributed to the spread of extended-spectrum beta-lactamase (ESBL) and has emerged as the dominant cause of hospital- and community-acquired urinary tract infections. Here, we report for the first time an in-depth analysis of whole-genome sequencing (WGS) of 4 ESBL-producing E. coli ST131 isolates recovered from patients in two hospitals in Armenia using Illumina short-read sequencing for accurate base calling to determine their genotype and to infer their phylogeny and using Oxford Nanopore Technologies long-read sequencing to resolve plasmid and chromosomal genetic elements. Genotypically, the four Armenian isolates were identified as part of the H30Rx/clade C2 (n = 2) and H41/clade A (n = 2) lineages and were phylogenetically closely related to isolates from the European Nucleotide Archive (ENA) database previously recovered from patients in the United States, Australia, and New Zealand. The Armenian isolates recovered in this study had chromosomal integration of the blaCTX-M-15 gene in the H30Rx isolates and a high number of virulence genes found in the H41 isolates associated with the carriage of a rare genomic island (in the context of E. coli ST131) containing the S fimbrial, salmochelin siderophore, and microcin H47 virulence genes. Furthermore, our data show the evolution of the IncF[2:A2:B20] plasmid harboring both blaCTX-M-15 and blaCTX-M-27 genes, derived from the recombination of genes from an IncF[F2:A-:B-] blaCTX-M-15-associated plasmid into the IncF[F1:A2:B20] blaCTX-M-27-associated plasmid backbone seen in two genetically closely related H41 Armenian isolates. IMPORTANCE Combining short and long reads from whole-genome sequencing analysis provided a genetic context for uncommon genes of clinical importance to better understand transmission and evolutionary features of ESBL-producing uropathogenic E. coli (UPEC) ST131 isolates recovered in Armenia. Using hybrid genome assembly in countries lacking genomic surveillance studies can inform us about new lineages not seen in other countries with genes encoding high virulence and antibiotic resistance harbored on mobile genetic elements.

bettercallsal: better calling of Salmonella serotypes from enrichment cultures using shotgun metagenomic profiling and its application in an outbreak setting

Most current Salmonella subtyping analyses rely on whole genome sequencing (WGS), which focuses on the high-resolution analysis of single genomes or multiple single genomes from the isolated colonies on microbiological agar plates. In this study, we introduce bioinformatics innovations for a metagenomic outbreak response workflow that accurately identifies multiple Salmonella serovars at the same time. bettercallsal is one of the first analysis tools to identify multiple Salmonella enterica serotypes from metagenomic or quasi-metagenomic datasets with high accuracy, allowing these isolate-independent methods to be incorporated into surveillance and root cause investigations. It was tested on an in silico benchmark dataset comprising 29 unique Salmonella serovars, 46 non-Salmonella bacterial genomes, and 10 viral genomes at varying read depths and on previously well-characterized and sequenced non-selective primary and selective enrichments of papaya and peach samples from separate outbreak investigations that resulted in the identification of multiple Salmonella serovars using traditional isolate culturing and WGS as well as nucleic acid assays. Analyses were also conducted on these datasets using a custom-built k-mer tool, SeqSero2, and Kallisto to compare serotype calling to bettercallsal. The in silico dataset analyzed with bettercallsal achieved the maximum precision, recall, and accuracy of 100, 83, and 94%, respectively. In the papaya outbreak samples, bettercallsal identified the presence of multiple serovars in agreement with the Luminex® xMAP assay results and also identified more serovars per sample, as evidenced by NCBI SNP clustering. In peach outbreak samples, bettercallsal identified two serovars in concordance with k-mer analysis and the Luminex xMAP assay. The genome hit reported by bettercallsal clustered with the chicken isolate genome, as reported by the FDA peach outbreak investigation from sequenced isolates (WGS). Overall, bettercallsal outperformed k-mer, Seqsero2, and Kallisto in identifying multiple serovars from enrichment cultures using shotgun metagenomic sequencing.

Molecular epidemiology of rotavirus among children in Western Canada: Dynamic changes in genotype prevalence in four consecutive seasons

Rotavirus molecular surveillance remains important in the postvaccine era to monitor the changes in transmission patterns, identify vaccine-induced antigenic changes and discover potentially pathogenic vaccine-related strains. The Canadian province of Alberta introduced rotavirus vaccination into its provincial vaccination schedule in June 2015. To evaluate the impact of this program on stool rotavirus positivity rate, strain diversity, and seasonal trends, we analyzed a prospective cohort of children with acute gastroenteritis recruited between December 2014 and August 2018. We identified dynamic changes in rotavirus positivity and genotype trends during pre- and post-rotavirus vaccine introduction periods. Genotypes G9P[8], G1P[8], G2P[4], and G12P[8] predominated consecutively each season with overall lower rotavirus incidence rates in 2016 and 2017. The demographic and clinical features of rotavirus gastroenteritis were comparable among wild-type rotaviruses; however, children with G12P[8] infections were older (p < 0.001). Continued efforts to monitor changes in the molecular epidemiology of rotavirus using whole genome sequence characterization are needed to further understand the impact of the selection pressure of vaccination on rotavirus evolution.

Genomic diversity and epidemiological significance of non-typhoidal Salmonella found in retail food collected in Norfolk, UK

Non-typhoidal Salmonella (NTS) is a major cause of bacterial gastroenteritis. Although many countries have implemented whole genome sequencing (WGS) of NTS, there is limited knowledge on NTS diversity on food and its contribution to human disease. In this study, the aim was to characterise the NTS genomes from retail foods in a particular region of the UK and assess the contribution to human NTS infections. Raw food samples were collected at retail in a repeated cross-sectional design in Norfolk, UK, including chicken (n=311), leafy green (n=311), pork (n=311), prawn (n=279) and salmon (n=157) samples. Up to eight presumptive NTS isolates per positive sample underwent WGS and were compared to publicly available NTS genomes from UK human cases. NTS was isolated from chicken (9.6 %), prawn (2.9 %) and pork (1.3 %) samples and included 14 serovars, of which Salmonella Infantis and Salmonella Enteritidis were the most common. The S. Enteritidis isolates were only isolated from imported chicken. No antimicrobial resistance determinants were found in prawn isolates, whilst 5.1 % of chicken and 0.64 % of pork samples contained multi-drug resistant NTS. The maximum number of pairwise core non-recombinant single nucleotide polymorphisms (SNPs) amongst isolates from the same sample was used to measure diversity and most samples had a median of two SNPs (range: 0-251). NTS isolates that were within five SNPs to clinical UK isolates belonged to specific serovars: S. Enteritidis and S. Infantis (chicken), and S. I 4,[5],12:i- (pork and chicken). Most NTS isolates that were closely related to human-derived isolates were obtained from imported chicken, but further epidemiological data are required to assess definitively the probable source of the human cases. Continued WGS surveillance of Salmonella on retail food involving multiple isolates from each sample is necessary to capture the diversity of Salmonella and determine the relative importance of different sources of human disease.

Genomic sketching with multiplicities and locality-sensitive hashing using Dashing 2

A genomic sketch is a small, probabilistic representation of the set of k-mers in a sequencing data set. Sketches are building blocks for large-scale analyses that consider similarities between many pairs of sequences or sequence collections. Although existing tools can easily compare tens of thousands of genomes, data sets can reach millions of sequences and beyond. Popular tools also fail to consider k-mer multiplicities, making them less applicable in quantitative settings. Here, we describe a method called Dashing 2 that builds on the SetSketch data structure. SetSketch is related to HyperLogLog (HLL) but discards use of leading zero count in favor of a truncated logarithm of adjustable base. Unlike HLL, SetSketch can perform multiplicity-aware sketching when combined with the ProbMinHash method. Dashing 2 integrates locality-sensitive hashing to scale all-pairs comparisons to millions of sequences. It achieves superior similarity estimates for the Jaccard coefficient and average nucleotide identity compared with the original Dashing, but in much less time while using the same-sized sketch. Dashing 2 is a free, open source software.

Whole genome sequencing of multidrug resistant Enterobacterales identified in children and their household members within Siem Reap, Cambodia

Objectives

To explore the association of recent hospitalization and asymptomatic carriage of multidrug-resistant Enterobacterales (MDRE) and determine the prevailing strains and antibiotic resistance genes in Siem Reap, Cambodia using WGS.

Methods

In this cross-sectional study, faecal samples were collected from two arms: a hospital-associated arm consisted of recently hospitalized children (2-14 years), with their family members; and a community-associated arm comprising children in the matching age group and their family members with no recent hospitalization. Forty-two families in each study arm were recruited, with 376 enrolled participants (169 adults and 207 children) and 290 stool specimens collected from participants. The DNA of ESBL- and carbapenemase-producing Enterobacterales cultured from the faecal samples was subject to WGS on the Illumina NovaSeq platform.

Results

Of the 290 stool specimens, 277 Escherichia coli isolates and 130 Klebsiella spp. were identified on CHROMagar ESBL and KPC plates. The DNA of 276 E. coli (one isolate failed quality control test), 89 Klebsiella pneumoniae, 40 Klebsiella quasipneumoniae and 1 Klebsiella variicola was sequenced. CTX-M-15 was the most common ESBL gene found in E. coli (n = 104, 38%), K. pneumoniae (n = 50, 56%) and K. quasipneumoniae (n = 16, 40%). The prevalence of bacterial lineages and ESBL genes was not associated with any specific arm.

Conclusions

Our results demonstrate that MDRE is likely to be endemic within the Siem Reap community. ESBL genes, specifically bla_CTX-M, can be found in almost all E. coli commensals, indicating that these genes are continuously propagated in the community through various unknown channels at present.

RabbitTClust: enabling fast clustering analysis of millions of bacteria genomes with MinHash sketches

We present RabbitTClust, a fast and memory-efficient genome clustering tool based on sketch-based distance estimation. Our approach enables efficient processing of large-scale datasets by combining dimensionality reduction techniques with streaming and parallelization on modern multi-core platforms. 113,674 complete bacterial genome sequences from RefSeq, 455 GB in FASTA format, can be clustered within less than 6 min and 1,009,738 GenBank assembled bacterial genomes, 4.0 TB in FASTA format, within only 34 min on a 128-core workstation. Our results further identify 1269 redundant genomes, with identical nucleotide content, in the RefSeq bacterial genomes database.

Arresting microbiome development limits immune system maturation and resistance to infection in mice

Disruptions to the intestinal microbiome during weaning lead to negative effects on host immune function. However, the critical host-microbe interactions during weaning that are required for immune system development remain poorly understood. We find that restricting microbiome maturation during weaning stunts immune system development and increases susceptibility to enteric infection. We developed a gnotobiotic mouse model of the early-life microbiome Pediatric Community (PedsCom). These mice develop fewer peripheral regulatory T cells and less IgA, hallmarks of microbiota-driven immune system development. Furthermore, adult PedsCom mice retain high susceptibility to Salmonella infection, which is characteristic of young mice and children. Altogether, our work illustrates how the post-weaning transition in microbiome composition contributes to normal immune maturation and protection from infection. Accurate modeling of the pre-weaning microbiome provides a window into the microbial requirements for healthy development and suggests an opportunity to design microbial interventions at weaning to improve immune development in human infants.

Navigating bottlenecks and trade-offs in genomic data analysis

Genome sequencing and analysis allow researchers to decode the functional information hidden in DNA sequences as well as to study cell to cell variation within a cell population. Traditionally, the primary bottleneck in genomic analysis pipelines has been the sequencing itself, which has been much more expensive than the computational analyses that follow. However, an important consequence of the continued drive to expand the throughput of sequencing platforms at lower cost is that often the analytical pipelines are struggling to keep up with the sheer amount of raw data produced. Computational cost and efficiency have thus become of ever increasing importance. Recent methodological advances, such as data sketching, accelerators and domain-specific libraries/languages, promise to address these modern computational challenges. However, despite being more efficient, these innovations come with a new set of trade-offs, both expected, such as accuracy versus memory and expense versus time, and more subtle, including the human expertise needed to use non-standard programming interfaces and set up complex infrastructure. In this Review, we discuss how to navigate these new methodological advances and their trade-offs.

Aspergillus luchuensis, an Endophyte Fungus from the Metal Hyperaccumulator Plant Prosopis laevigata, Promotes Its Growth and Increases Metal Translocation

Heavy metal pollution is a worldwide environmental and human health problem. Prosopis laevigata is a hyperaccumulator legume that bioaccumulates Pb, Cu and Zn. With interest in designing phytoremediation strategies for sites contaminated with heavy metals, we isolated and characterized endophytic fungi from the roots of P. laevigata growing on mine tailings located in Morelos, Mexico. Ten endophytic isolates were selected by morphological discrimination and a preliminary minimum inhibitory concentration was determined for zinc, lead and copper. A novel strain of Aspergillus closest to Aspergillus luchuensis was determined to be a metallophile and presented a marked tolerance to high concentrations of Cu, Zn and Pb, so it was further investigated for removal of metals and promotion of plant growth under greenhouse conditions. The control substrate with fungi promoted larger size characters in P. laevigata individuals in comparison with the other treatments, demonstrating that A. luchuensis strain C7 is a growth-promoting agent for P. laevigata individuals. The fungus favors the translocation of metals from roots to leaves in P. laevigata, promoting an increased Cu translocation. This new A. luchuensis strain showed endophytic character and plant growth-promotion activity, high metal tolerance, and an ability to increase copper translocation. We propose it as a novel, effective and sustainable bioremediation strategy for copper-polluted soils.

Molecular surveillance reveals the emergence and dissemination of NDM-5-producing Escherichia coli high-risk clones in Germany, 2013 to 2019

BackgroundCarbapenemase-producing Enterobacterales (CPE) are rapidly increasing worldwide, also in Europe. Although prevalence of CPE in Germany is comparatively low, the National Reference Centre for Multidrug-resistant Gram-negative Bacteria noted annually increasing numbers of NDM-5-producing Escherichia coli isolates.AimAs part of our ongoing surveillance programme, we characterised NDM-5-producing E. coli isolates received between 2013 and 2019 using whole genome sequencing (WGS).MethodsFrom 329 identified NDM-5-producing E. coli, 224 isolates from known geographical locations were subjected to Illumina WGS. Analyses of 222 sequenced isolates included multilocus sequence typing (MLST), core genome (cg)MLST and single-nucleotide polymorphism (SNP)-based analyses.ResultsResults of cgMLST revealed genetically distinct clusters for many of the 43 detected sequence types (ST), of which ST167, ST410, ST405 and ST361 predominated. The SNP-based phylogenetic analyses combined with geographical information identified sporadic cases of nosocomial transmission on a small spatial scale. However, we identified large clusters corresponding to clonal dissemination of ST167, ST410, ST405 and ST361 strains in consecutive years in different regions in Germany.ConclusionOccurrence of NDM-5-producing E. coli rose in Germany, which was to a large extent due to the increased prevalence of isolates belonging to the international high-risk clones ST167, ST410, ST405 and ST361. Of particular concern is the supra-regional dissemination of these epidemic clones. Available information suggest community spread of NDM-5-producing E. coli in Germany, highlighting the importance of epidemiological investigation and an integrated surveillance system in the One Health framework.

JC polyomavirus: a short review of its biology, its association with progressive multifocal leukoencephalopathy, and the diagnostic value of different methods to manifest its activity or presence

INTRODUCTION

JC polyomavirus is the causative agent of progressive multifocal leukoencephalopathy (PML), a demyelinating disease resulting from the lytic infection of oligodendrocytes that may develop in immunosuppressed individuals: HIV1 infected or individuals under immunosuppressive therapies. Understanding the biology of JCPyV is necessary for a proper patient management, the development of diagnostic tests, and risk stratification.

AREAS COVERED

The review covers different areas of expertise including the genomic characterization of JCPyV strains detected in different body compartments (urine, plasma, and cerebrospinal fluid) of PML patients, viral mutations, molecular diagnostics, viral miRNAs, and disease.

EXPERT OPINION

The implementation of molecular biology techniques improved our understanding of JCPyV biology. Deep sequencing analysis of viral genomes revealed the presence of viral quasispecies in the cerebrospinal fluid of PML patients characterized by noncoding control region rearrangements and VP1 mutations. These neurotropic JCPyV variants present enhanced replication and an altered cell tropism that contribute to PML development. Monitoring these variants may be relevant for the identification of patients at risk of PML. Multiplex realtime PCR targeting both the LTAg and the archetype NCCR could be used to identify them. Failure to amplify NCCR should indicate the presence of a JCPyV prototype speeding up the diagnostic process.

Pre-epidemic evolution of the MRSA USA300 clade and a molecular key for classification

Introduction

USA300 has remained the dominant community and healthcare associated methicillin-resistant Staphylococcus aureus (MRSA) clone in the United States and in northern South America for at least the past 20 years. In this time, it has experienced epidemic spread in both of these locations. However, its pre-epidemic evolutionary history and origins are incompletely understood. Large sequencing databases, such as NCBI, PATRIC, and Staphopia, contain clues to the early evolution of USA300 in the form of sequenced genomes of USA300 isolates that are representative of lineages that diverged prior to the establishment of the South American epidemic (SAE) clade and North American epidemic (NAE) clade. In addition, historical isolates collected prior to the emergence of epidemics can help reconstruct early events in the history of this lineage.

Methods

Here, we take advantage of the accrued, publicly available data, as well as two newly sequenced pre-epidemic historical isolates from 1996, and a very early diverging ACME-negative NAE genome, to understand the pre-epidemic evolution of USA300. We use database mining techniques to emphasize genomes similar to pre-epidemic isolates, with the goal of reconstructing the early molecular evolution of the USA300 lineage.

Results

Phylogenetic analysis with these genomes confirms that the NAE and SAE USA300 lineages diverged from a most recent common ancestor around 1970 with high confidence, and it also pinpoints the independent acquisition events of the of the ACME and COMER loci with greater precision than in previous studies. We provide evidence for a North American origin of the USA300 lineage and identify multiple introductions of USA300 into South and North America. Notably, we describe a third major USA300 clade (the pre-epidemic branching clade; PEB1) consisting of both MSSA and MRSA isolates circulating around the world that diverged from the USA300 lineage prior to the establishment of the South and North American epidemics. We present a detailed analysis of specific sequence characteristics of each of the major clades, and present diagnostic positions that can be used to classify new genomes.