Efficient and Robust Search of Microbial Genomes via Phylogenetic Compression.

Comprehensive collections approaching millions of sequenced genomes have become central information sources in the life sciences. However, the rapid growth of these collections makes it effectively impossible to search these data using tools such as BLAST and its successors. Here, we present a technique called phylogenetic compression, which uses evolutionary history to guide compression and efficiently search large collections of microbial genomes using existing algorithms and data structures. We show that, when applied to modern diverse collections approaching millions of genomes, lossless phylogenetic compression improves the compression ratios of assemblies, de Bruijn graphs, and k-mer indexes by one to two orders of magnitude. Additionally, we develop a pipeline for a BLAST-like search over these phylogeny-compressed reference data, and demonstrate it can align genes, plasmids, or entire sequencing experiments against all sequenced bacteria until 2019 on ordinary desktop computers within a few hours. Phylogenetic compression has broad applications in computational biology and may provide a fundamental design principle for future genomics infrastructure.

High-resolution genomic surveillance elucidates a multilayered hierarchical transfer of resistance between WWTP- and human/animal-associated bacteria.

BACKGROUND: Our interconnected world and the ability of bacteria to quickly swap antibiotic resistance genes (ARGs) make it particularly important to establish the epidemiological links of multidrug resistance (MDR) transfer between wastewater treatment plant (WWTP)- and human/animal-associated bacteria, under the One Health framework. However, evidence of ARGs exchange and potential factors that contribute to this transfer remain limited. RESULTS: Here, by combining culture-based population genomics and genetic comparisons with publicly available datasets, we reconstructed the complete genomes of 82 multidrug-resistant isolates from WWTPs and found that most WWTP-associated isolates were genetically distinct from their closest human/animal-associated relatives currently available in the public database. Even in the minority of lineages that were closely related, WWTP-associated isolates were characterized by quite different plasmid compositions. We identified a high diversity of circular plasmids (264 in total, of which 141 were potentially novel), which served as the main source of resistance, and showed potential horizontal transfer of ARG-bearing plasmids between WWTP- and humans/animal-associated bacteria. Notably, the potentially transferred ARGs and virulence factors (VFs) with different genetic backgrounds were closely associated with flanking insertion sequences (ISs), suggesting the importance of synergy between plasmids and ISs in mediating a multilayered hierarchical transfer of MDR and potentiating the emergence of MDR-hypervirulent clones. CONCLUSION: Our findings advance the current efforts to establish potential epidemiological links of MDR transmission between WWTP- and human/animal-associated bacteria. Plasmids play an important role in mediating the transfer of ARGs and the IS-associated ARGs that are carried by conjugative plasmids should be prioritized to tackle the spread of resistance. Video Abstract.

SARS-CoV-2 epidemic after social and economic reopening in three U.S. states reveals shifts in age structure and clinical characteristics.

State-level reopenings in late spring 2020 facilitated the resurgence of severe acute respiratory syndrome coronavirus 2 transmission. Here, we analyze age-structured case, hospitalization, and death time series from three states-Rhode Island, Massachusetts, and Pennsylvania-that had successful reopenings in May 2020 without summer waves of infection. Using 11 daily data streams, we show that from spring to summer, the epidemic shifted from an older to a younger age profile and that elderly individuals were less able to reduce contacts during the lockdown period when compared to younger individuals. Clinical case management improved from spring to summer, resulting in fewer critical care admissions and lower infection fatality rate. Attack rate estimates through 31 August 2020 are 6.2% [95% credible interval (CI), 5.7 to 6.8%] of the total population infected for Rhode Island, 6.7% (95% CI, 5.4 to 7.6%) in Massachusetts, and 2.7% (95% CI, 2.5 to 3.1%) in Pennsylvania.

SARS-CoV-2 epidemic after social and economic reopening in three US states reveals shifts in age structure and clinical characteristics.

In the United States, state-level re-openings in spring 2020 presented an opportunity for the resurgence of SARS-CoV-2 transmission. One important question during this time was whether human contact and mixing patterns could increase gradually without increasing viral transmission, the rationale being that new mixing patterns would likely be associated with improved distancing, masking, and hygiene practices. A second key question to follow during this time was whether clinical characteristics of the epidemic would improve after the initial surge of cases. Here, we analyze age-structured case, hospitalization, and death time series from three states - Rhode Island, Massachusetts, and Pennsylvania - that had successful re-openings in May 2020 without summer waves of infection. Using a Bayesian inference framework on eleven daily data streams and flexible daily population contact parameters, we show that population-average mixing rates dropped by >50% during the lockdown period in March/April, and that the correlation between overall population mobility and transmission-capable mobility was broken in May as these states partially re-opened. We estimate the reporting rates (fraction of symptomatic cases reporting to health system) at 96.0% (RI), 72.1% (MA), and 75.5% (PA); in Rhode Island, when accounting for cases caught through general-population screening programs, the reporting rate estimate is 94.5%. We show that elderly individuals were less able to reduce contacts during the lockdown period when compared to younger individuals. Attack rate estimates through August 31 2020 are 6.4% (95% CI: 5.8% ‒ 7.3%) of the total population infected for Rhode Island, 5.7% (95% CI: 5.0% ‒ 6.8%) in Massachusetts, and 3.7% (95% CI: 3.1% ‒ 4.5%) in Pennsylvania, with some validation available through published seroprevalence studies. Infection fatality rates (IFR) estimates for the spring epidemic are higher in our analysis (>2%) than previously reported values, likely resulting from the epidemics in these three states affecting the most vulnerable sub-populations, especially the most vulnerable of the ≥80 age group.

Optimal SARS-CoV-2 vaccine allocation using real-time attack-rate estimates in Rhode Island and Massachusetts.

BACKGROUND: When three SARS-CoV-2 vaccines came to market in Europe and North America in the winter of 2020-2021, distribution networks were in a race against a major epidemiological wave of SARS-CoV-2 that began in autumn 2020. Rapid and optimized vaccine allocation was critical during this time. With 95% efficacy reported for two of the vaccines, near-term public health needs likely require that distribution is prioritized to the elderly, health care workers, teachers, essential workers, and individuals with comorbidities putting them at risk of severe clinical progression. METHODS: We evaluate various age-based vaccine distributions using a validated mathematical model based on current epidemic trends in Rhode Island and Massachusetts. We allow for varying waning efficacy of vaccine-induced immunity, as this has not yet been measured. We account for the fact that known COVID-positive cases may not have been included in the first round of vaccination. And, we account for age-specific immune patterns in both states at the time of the start of the vaccination program. Our analysis assumes that health systems during winter 2020-2021 had equal staffing and capacity to previous phases of the SARS-CoV-2 epidemic; we do not consider the effects of understaffed hospitals or unvaccinated medical staff. RESULTS: We find that allocating a substantial proportion (>75%) of vaccine supply to individuals over the age of 70 is optimal in terms of reducing total cumulative deaths through mid-2021. This result is robust to different profiles of waning vaccine efficacy and several different assumptions on age mixing during and after lockdown periods. As we do not explicitly model other high-mortality groups, our results on vaccine allocation apply to all groups at high risk of mortality if infected. A median of 327 to 340 deaths can be avoided in Rhode Island (3444 to 3647 in Massachusetts) by optimizing vaccine allocation and vaccinating the elderly first. The vaccination campaigns are expected to save a median of 639 to 664 lives in Rhode Island and 6278 to 6618 lives in Massachusetts in the first half of 2021 when compared to a scenario with no vaccine. A policy of vaccinating only seronegative individuals avoids redundancy in vaccine use on individuals that may already be immune, and would result in 0.5% to 1% reductions in cumulative hospitalizations and deaths by mid-2021. CONCLUSIONS: Assuming high vaccination coverage (>28%) and no major changes in distancing, masking, gathering size, hygiene guidelines, and virus transmissibility between 1 January 2021 and 1 July 2021 a combination of vaccination and population immunity may lead to low or near-zero transmission levels by the second quarter of 2021.

Simplitigs as an efficient and scalable representation of de Bruijn graphs.

de Bruijn graphs play an essential role in bioinformatics, yet they lack a universal scalable representation. Here, we introduce simplitigs as a compact, efficient, and scalable representation, and ProphAsm, a fast algorithm for their computation. For the example of assemblies of model organisms and two bacterial pan-genomes, we compare simplitigs to unitigs, the best existing representation, and demonstrate that simplitigs provide a substantial improvement in the cumulative sequence length and their number. When combined with the commonly used Burrows-Wheeler Transform index, simplitigs reduce memory, and index loading and query times, as demonstrated with large-scale examples of GenBank bacterial pan-genomes.

Conjugative plasmids interact with insertion sequences to shape the horizontal transfer of antimicrobial resistance genes.

It is well established that plasmids play an important role in the dissemination of antimicrobial resistance (AMR) genes; however, little is known about the role of the underlying interactions between different plasmid categories and other mobile genetic elements (MGEs) in shaping the promiscuous spread of AMR genes. Here, we developed a tool designed for plasmid classification, AMR gene annotation, and plasmid visualization and found that most plasmid-borne AMR genes, including those localized on class 1 integrons, are enriched in conjugative plasmids. Notably, we report the discovery and characterization of a massive insertion sequence (IS)-associated AMR gene transfer network (245 combinations covering 59 AMR gene subtypes and 53 ISs) linking conjugative plasmids and phylogenetically distant pathogens, suggesting a general evolutionary mechanism for the horizontal transfer of AMR genes mediated by the interaction between conjugative plasmids and ISs. Moreover, our experimental results confirmed the importance of the observed interactions in aiding the horizontal transfer and expanding the genetic range of AMR genes within complex microbial communities.

Optimal SARS-CoV-2 vaccine allocation using real-time seroprevalence estimates in Rhode Island and Massachusetts.

As three SARS-CoV-2 vaccines come to market in Europe and North America in the winter of 2020-2021, distribution networks will be in a race against a major epidemiological wave of SARS-CoV-2 that began in autumn 2020. Rapid and optimized vaccine allocation is critical during this time. With 95% efficacy reported for two of the vaccines, near-term public health needs require that distribution is prioritized to the elderly, health-care workers, teachers, essential workers, and individuals with co-morbidities putting them at risk of severe clinical progression. Here, we evaluate various age-based vaccine distributions using a validated mathematical model based on current epidemic trends in Rhode Island and Massachusetts. We allow for varying waning efficacy of vaccine-induced immunity, as this has not yet been measured. We account for the fact that known COVID-positive cases may not be included in the first round of vaccination. And, we account for current age-specific immune patterns in both states. We find that allocating a substantial proportion ( > 75%) of vaccine supply to individuals over the age of 70 is optimal in terms of reducing total cumulative deaths through mid-2021. As we do not explicitly model other high mortality groups, this result on vaccine allocation applies to all groups at high risk of mortality if infected. Our analysis confirms that for an easily transmissible respiratory virus, allocating a large majority of vaccinations to groups with the highest mortality risk is optimal. Our analysis assumes that health systems during winter 2020-2021 have equal staffing and capacity to previous phases of the SARS-CoV-2 epidemic; we do not consider the effects of understaffed hospitals or unvaccinated medical staff. Vaccinating only seronegative individuals avoids redundancy in vaccine use on individuals that may already be immune, and will result in 1% to 2% reductions in cumulative hospitalizations and deaths by mid-2021. Assuming high vaccination coverage ( > 28%) and no major relaxations in distancing, masking, gathering size, or hygiene guidelines between now and spring 2021, our model predicts that a combination of vaccination and population immunity will lead to low or near-zero transmission levels by the second quarter of 2021.

Using Genetic Distance from Archived Samples for the Prediction of Antibiotic Resistance in Escherichia coli.

The rising rates of antibiotic resistance increasingly compromise empirical treatment. Knowing the antibiotic susceptibility of a pathogen's close genetic relative(s) may improve empirical antibiotic selection. Using genomic and phenotypic data for Escherichia coli isolates from three separate clinically derived databases, we evaluated multiple genomic methods and statistical models for predicting antibiotic susceptibility, focusing on potentially rapidly available information, such as lineage or genetic distance from archived isolates. We applied these methods to derive and validate the prediction of antibiotic susceptibility to common antibiotics. We evaluated 968 separate episodes of suspected and confirmed infection with Escherichia coli from three geographically and temporally separated databases in Ontario, Canada, from 2010 to 2018. Across all approaches, model performance (area under the curve [AUC]) ranges for predicting antibiotic susceptibility were the greatest for ciprofloxacin (AUC, 0.76 to 0.97) and the lowest for trimethoprim-sulfamethoxazole (AUC, 0.51 to 0.80). When a model predicted that an isolate was susceptible, the resulting (posttest) probabilities of susceptibility were sufficient to warrant empirical therapy for most antibiotics (mean, 92%). An approach combining multiple models could permit the use of narrower-spectrum oral agents in 2 out of every 3 patients while maintaining high treatment adequacy (∼90%). Methods based on genetic relatedness to archived samples of E. coli could be used to predict antibiotic resistance and improve antibiotic selection.

Rapid inference of antibiotic resistance and susceptibility by genomic neighbour typing.

Surveillance of drug-resistant bacteria is essential for healthcare providers to deliver effective empirical antibiotic therapy. However, traditional molecular epidemiology does not typically occur on a timescale that could affect patient treatment and outcomes. Here, we present a method called 'genomic neighbour typing' for inferring the phenotype of a bacterial sample by identifying its closest relatives in a database of genomes with metadata. We show that this technique can infer antibiotic susceptibility and resistance for both Streptococcus pneumoniae and Neisseria gonorrhoeae. We implemented this with rapid k-mer matching, which, when used on Oxford Nanopore MinION data, can run in real time. This resulted in the determination of resistance within 10 min (91% sensitivity and 100% specificity for S. pneumoniae and 81% sensitivity and 100% specificity for N. gonorrhoeae from isolates with a representative database) of starting sequencing, and within 4 h of sample collection (75% sensitivity and 100% specificity for S. pneumoniae) for clinical metagenomic sputum samples. This flexible approach has wide application for pathogen surveillance and may be used to greatly accelerate appropriate empirical antibiotic treatment.

RNF: a general framework to evaluate NGS read mappers.

MOTIVATION: Read simulators combined with alignment evaluation tools provide the most straightforward way to evaluate and compare mappers. Simulation of reads is accompanied by information about their positions in the source genome. This information is then used to evaluate alignments produced by the mapper. Finally, reports containing statistics of successful read alignments are created.In default of standards for encoding read origins, every evaluation tool has to be made explicitly compatible with the simulator used to generate reads. RESULTS: To solve this obstacle, we have created a generic format Read Naming Format (Rnf) for assigning read names with encoded information about original positions. Futhermore, we have developed an associated software package RnfTools containing two principal components. MIShmash applies one of popular read simulating tools (among DwgSim, Art, Mason, CuReSim, etc.) and transforms the generated reads into Rnf format. LAVEnder evaluates then a given read mapper using simulated reads in Rnf format. A special attention is payed to mapping qualities that serve for parametrization of Roc curves, and to evaluation of the effect of read sample contamination. AVAILABILITY AND IMPLEMENTATION: RnfTools: http://karel-brinda.github.io/rnftools Spec. of Rnf: http://karel-brinda.github.io/rnf-spec CONTACT: karel.brinda@univ-mlv.fr.

Spaced seeds improve k-mer-based metagenomic classification.

MOTIVATION: Metagenomics is a powerful approach to study genetic content of environmental samples, which has been strongly promoted by next-generation sequencing technologies. To cope with massive data involved in modern metagenomic projects, recent tools rely on the analysis of k-mers shared between the read to be classified and sampled reference genomes. RESULTS: Within this general framework, we show that spaced seeds provide a significant improvement of classification accuracy, as opposed to traditional contiguous k-mers. We support this thesis through a series of different computational experiments, including simulations of large-scale metagenomic projects.Availability and implementation, Supplementary information: Scripts and programs used in this study, as well as supplementary material, are available from http://github.com/gregorykucherov/spaced-seeds-for-metagenomics. CONTACT: gregory.kucherov@univ-mlv.fr.