Icarust, a real-time simulator for Oxford Nanopore adaptive sampling.

MOTIVATION: Oxford Nanopore Technologies (ONT) sequencers enable real-time generation of sequence data, which allows for concurrent analysis during a run. Adaptive sampling leverages this real-time capability in extremis, rejecting or accepting reads for sequencing based on assessment of the sequence from the start of each read. This functionality is provided by ONT's software, MinKNOW (Oxford Nanopore Technologies). Designing and developing software to take advantage of adaptive sampling can be costly in terms of sequencing consumables, using precious samples and preparing sequencing libraries. MinKNOW addresses this in part by allowing the replay of previously sequenced runs for testing. However, as we show, the sequencing output only partially changes in response to adaptive sampling instructions. Here we present Icarust, a tool enabling more accurate approximations of sequencing runs. Icarust recreates all the required endpoints of MinKNOW to perform adaptive sampling and writes output compatible with current base-callers and analysis pipelines. Icarust serves nanopore signal simulating a MinION or PromethION flow cell experiment from any reference genome using either R9 or R10 pore models. We show that simulating sequencing runs with Icarust provides a realistic testing and development environment for software exploiting the real-time nature of Nanopore sequencing. AVAILABILITY AND IMPLEMENTATION: All code is open source and freely available here-https://github.com/LooseLab/Icarust. Icarust is implemented in Rust, with a docker container also available. The data underlying this article will be shared on reasonable request to the corresponding author.

minoTour, real-time monitoring and analysis for nanopore sequencers.

SUMMARY: minoTour offers a Laboratory Informations Management System (LIMS) system for Oxford Nanopore Technology sequencers, with real-time metrics and analysis available permanently for review. Integration of unique real-time automated analysis can reduce the time required to answer biological questions, including mapping and classification of sequence while a run is in progress. Real-time sequence data require new methods of analysis which do not wait for the completion of a run and minoTour provides a framework to allow users to exploit these features. AVAILABILITY AND IMPLEMENTATION: Source code and documentation are available at https://github.com/LooseLab/minotourcli and https://github.com/LooseLab/minotourapp. Docker images are available from https://hub.docker.com/r/adoni5/, and can be installed using a preconfigured docker-compose script at https://github.com/LooseLab/minotour-docker. An example server is available at http://137.44.59.170. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The Stabilization of S100A9 Structure by Calcium Inhibits the Formation of Amyloid Fibrils.

The calcium-binding protein S100A9 is recognized as an important component of the brain neuroinflammatory response to the onset and development of neurodegenerative disease. S100A9 is intrinsically amyloidogenic and in vivo co-aggregates with amyloid-ß peptide and α-synuclein in Alzheimer's and Parkinson's diseases, respectively. It is widely accepted that calcium dyshomeostasis plays an important role in the onset and development of these diseases, and studies have shown that elevated levels of calcium limit the potential for S100A9 to adopt a fibrillar structure. The exact mechanism by which calcium exerts its influence on the aggregation process remains unclear. Here we demonstrate that despite S100A9 exhibiting α-helical secondary structure in the absence of calcium, the protein exhibits significant plasticity with interconversion between different conformational states occurring on the micro- to milli-second timescale. This plasticity allows the population of conformational states that favour the onset of fibril formation. Magic-angle spinning solid-state NMR studies of the resulting S100A9 fibrils reveal that the S100A9 adopts a single structurally well-defined rigid fibrillar core surrounded by a shell of approximately 15-20 mobile residues, a structure that persists even when fibrils are produced in the presence of calcium ions. These studies highlight how the dysregulation of metal ion concentrations can influence the conformational equilibria of this important neuroinflammatory protein to influence the rate and nature of the amyloid deposits formed.

Oligomerization processes limit photoactivation and recovery of the orange carotenoid protein.

The orange carotenoid protein (OCP) is a photoactive protein involved in cyanobacterial photoprotection by quenching of the excess of light-harvested energy. The photoactivation mechanism remains elusive, in part due to absence of data pertaining to the timescales over which protein structural changes take place. It also remains unclear whether or not oligomerization of the dark-adapted and light-adapted OCP could play a role in the regulation of its energy-quenching activity. Here, we probed photoinduced structural changes in OCP by a combination of static and time-resolved X-ray scattering and steady-state and transient optical spectroscopy in the visible range. Our results suggest that oligomerization partakes in regulation of the OCP photocycle, with different oligomers slowing down the overall thermal recovery of the dark-adapted state of OCP. They furthermore reveal that upon non-photoproductive excitation a numbed state forms, which remains in a non-photoexcitable structural state for at least ≈0.5 µs after absorption of a first photon.

Dynamic, adaptive sampling during nanopore sequencing using Bayesian experimental design.

Nanopore sequencers can select which DNA molecules to sequence, rejecting a molecule after analysis of a small initial part. Currently, selection is based on predetermined regions of interest that remain constant throughout an experiment. Sequencing efforts, thus, cannot be re-focused on molecules likely contributing most to experimental success. Here we present BOSS-RUNS, an algorithmic framework and software to generate dynamically updated decision strategies. We quantify uncertainty at each genome position with real-time updates from data already observed. For each DNA fragment, we decide whether the expected decrease in uncertainty that it would provide warrants fully sequencing it, thus optimizing information gain. BOSS-RUNS mitigates coverage bias between and within members of a microbial community, leading to improved variant calling; for example, low-coverage sites of a species at 1% abundance were reduced by 87.5%, with 12.5% more single-nucleotide polymorphisms detected. Such data-driven updates to molecule selection are applicable to many sequencing scenarios, such as enriching for regions with increased divergence or low coverage, reducing time-to-answer.

A framework for real-time monitoring, analysis and adaptive sampling of viral amplicon nanopore sequencing.

The ongoing SARS-CoV-2 pandemic demonstrates the utility of real-time sequence analysis in monitoring and surveillance of pathogens. However, cost-effective sequencing requires that samples be PCR amplified and multiplexed via barcoding onto a single flow cell, resulting in challenges with maximising and balancing coverage for each sample. To address this, we developed a real-time analysis pipeline to maximise flow cell performance and optimise sequencing time and costs for any amplicon based sequencing. We extended our nanopore analysis platform MinoTour to incorporate ARTIC network bioinformatics analysis pipelines. MinoTour predicts which samples will reach sufficient coverage for downstream analysis and runs the ARTIC networks Medaka pipeline once sufficient coverage has been reached. We show that stopping a viral sequencing run earlier, at the point that sufficient data has become available, has no negative effect on subsequent down-stream analysis. A separate tool, SwordFish, is used to automate adaptive sampling on Nanopore sequencers during the sequencing run. This enables normalisation of coverage both within (amplicons) and between samples (barcodes) on barcoded sequencing runs. We show that this process enriches under-represented samples and amplicons in a library as well as reducing the time taken to obtain complete genomes without affecting the consensus sequence.

Readfish enables targeted nanopore sequencing of gigabase-sized genomes.

Nanopore sequencers can be used to selectively sequence certain DNA molecules in a pool by reversing the voltage across individual nanopores to reject specific sequences, enabling enrichment and depletion to address biological questions. Previously, we achieved this using dynamic time warping to map the signal to a reference genome, but the method required substantial computational resources and did not scale to gigabase-sized references. Here we overcome this limitation by using graphical processing unit (GPU) base-calling. We show enrichment of specific chromosomes from the human genome and of low-abundance organisms in mixed populations without a priori knowledge of sample composition. Finally, we enrich targeted panels comprising 25,600 exons from 10,000 human genes and 717 genes implicated in cancer, identifying PML-RARA fusions in the NB4 cell line in <15 h sequencing. These methods can be used to efficiently screen any target panel of genes without specialized sample preparation using any computer and a suitable GPU. Our toolkit, readfish, is available at https://www.github.com/looselab/readfish .