CarrierSeq: a sequence analysis workflow for low-input nanopore sequencing
BMC Bioinformatics volume 19, Article number: 108 (2018)
Long-read nanopore sequencing technology is of particular significance for taxonomic identification at or below the species level. For many environmental samples, the total extractable DNA is far below the current input requirements of nanopore sequencing, preventing “sample to sequence” metagenomics from low-biomass or recalcitrant samples.
Here we address this problem by employing carrier sequencing, a method to sequence low-input DNA by preparing the target DNA with a genomic carrier to achieve ideal library preparation and sequencing stoichiometry without amplification. We then use CarrierSeq, a sequence analysis workflow to identify the low-input target reads from the genomic carrier. We tested CarrierSeq experimentally by sequencing from a combination of 0.2 ng Bacillus subtilis ATCC 6633 DNA in a background of 1000 ng Enterobacteria phage λ DNA. After filtering of carrier, low quality, and low complexity reads, we detected target reads (B. subtilis), contamination reads, and “high quality noise reads” (HQNRs) not mapping to the carrier, target or known lab contaminants. These reads appear to be artifacts of the nanopore sequencing process as they are associated with specific channels (pores).
By treating sequencing as a Poisson arrival process, we implement a statistical test to reject data from channels dominated by HQNRs while retaining low-input target reads.
Environmental metagenomic sequencing poses a number of challenges. First, complex soil matrices and tough-to-lyse organisms can frustrate the extraction of deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) . Second, low-biomass samples require further extraction and concentration steps which increase the likelihood of contamination . Third, whole genome amplification may bias population results  while targeted amplification (e.g., 16S rRNA amplicon) may decrease taxonomic resolution . To address these challenges, we have developed extraction protocols compatible with low-biomass recalcitrant samples and difficult to lyse organisms . These protocols, developed using tough-to-lyse spores of Bacillus subtilis, allow us to achieve at least 5% extraction yield from a 50 mg sample containing 2 × 105 cells/g of soil without centrifugation . Furthermore, in order to avoid possible amplification biases and additional points of contamination, we have experimented with utilizing a genomic carrier (Enterobacteria phage λ) to shuttle low-input amounts of target DNA (B. subtilis) through library preparation and sequencing with ideal stoichiometry . This approach has allowed us to detect down to 0.2 ng of B. subtilis DNA prepared with 1000 ng of Lambda DNA using the Oxford Nanopore Technologies (ONT) MinION sequencer. Here we present CarrierSeq, a sequence analysis workflow developed to identify target reads from a low-input sequencing run employing a genomic carrier.
CarrierSeq implements bwa-mem  to first map all reads to the genomic carrier then extracts unmapped reads by using samtools  and seqtk . Thereafter, the user can define a quality score threshold and CarrierSeq proceeds to discard low-complexity reads  with fqtrim . This set of unmapped and filtered reads are labeled “reads of interest” (ROI) and should theoretically comprise target reads and likely contamination. However, ROIs also include “high-quality noise reads” (HQNRs), defined as reads that satisfy quality score and complexity filters yet do not match to any database and disproportionately originate from specific channels. By treating reads as a Poisson arrival process, CarrierSeq models the expected ROIs channel distribution and rejects data from channels exceeding a reads/channels threshold (xcrit).
Quality score filter
The default per-read quality score threshold (Q9) was determined through receiver operating characteristic curve (ROC) analysis  of carrier sequencing runs of B. subtilis and Lambda DNA (Fig. 1). This threshold is best suited for Lambda carriers that are 99% library by mass and essentially function as a pseudo “lambda burn-in”. Therefore, the user is encouraged to define their own threshold based on their libraries’ quality control metrics (e.g., quality distribution, sequencing accuracy achieved, and basecaller confidence).
Assuming that sequencing is a stochastic process, CarrierSeq is able to identify channels producing spurious reads by calculating the expected Poisson distribution of reads/channel. Given total ROIs and number of active sequencing channels, CarrierSeq will determine the arrival rate (λ = reads of interest/active channels). CarrierSeq then calculates an xcrit threshold (xcrit = poisson.ppf (1 – p-value), λ)) and sorts ROIs into target reads (reads/channel ≤ xcrit) or HQNRs (reads/channel > xcrit).
Here we test CarrierSeq by analyzing carrier sequencing data from a library containing 0.2 ng of B. subtilis DNA prepared with 1000 ng of Lambda DNA using the Oxford Nanopore Technologies (ONT) ligation sequencing kit (LSK-SQK108). Following the standard Nanopore Lambda calibration or “burn in” protocol recommended for every new Nanopore user, B. subtilis DNA was used in place of the 3.6 kb positive control DNA. The library was then sequenced on a MinION Mark-1B sequencer and R9.4 flowcell for 48 h and basecalled using ONT’s Albacore (v1.10) offline basecaller.
From the resulting 48 h of sequencing, we detected a total of 718,432 reads or 6.4 gigabases. Exactly 676,086 reads mapped to Lambda, 777 reads mapped to B. subtilis, and 41,569 reads mapped to neither.
ROIs and sorting
Applying the parameters p = 0.05 and q = 9, CarrierSeq identified 1811 ROIs and determined xcrit = 7. Therefore, channels producing greater than 7 reads were identified as HQNR-associated while channels producing less than or equal to 7 reads were identified as “good” channels (Fig. 2). CarrierSeq then sorted 1179 reads, including 1162 true negative reads (real HQNRs) and 17 false negative reads (B. subtilis), as likely HQNRs. The final 632 target reads consisted of 574 true positive reads (574 B. subtilis and 4 Homo sapiens) and 54 false positive reads (HQNRs). Overall, CarrierSeq identified 74% of all B. subtilis reads present. Moreover, from the discarded 203 B. subtilis reads, 186 were below Q9 while 17 originated from discarded HQNR-associated channels.
From experimenting with low-input carrier sequencing and CarrierSeq we observed that the abundance of HQNRs may vary per run, perhaps due to sub-optimal library preparation, delays in initializing sequencing, or other sequencing conditions. In addition, target DNA purity and lysis carryover (e.g., proteins) may conceivably contribute to HQNR abundance possibly due to pore blockages from unknown macromolecules that result in erroneous reads. While the cause or significance of HQNRs have yet to be determined, future work will focus on developing a method to identify HQNRs on a per-read basis. In contrast, the current approach discards entire HQNR-associated channels at the risk of discarding target reads. Moreover, some reads in non-HQNR-associated channels may also be artifacts. The ability to identify HQNRs on a per-read basis is especially important for metagenomic studies of novel microbial communities where HQNRs may complicate the identification of an unknown organism, or in a life detection application  where artefactual reads not mapping to known life could represent a false-positive.
CarrierSeq was developed to analyze low-input carrier sequencing data and identify target reads. We have since deployed CarrierSeq to test the limits of detection of ONT’s MinION sequencer from 0.2 ng down to 2 pg of low-input carrier sequencing. CarrierSeq may be a particularly valuable tool for in-situ metagenomic studies where limited sample availability (e.g., low biomass environmental samples) and laboratory resources (i.e., field deployments) may benefit from sequencing with a genomic carrier.
Availability and requirements
Project name: CarrierSeq.
Project home page: https://github.com/amojarro/carrierseq
Operating system(s): macOS and Linux.
Programming language: BASH and Python.
Other requirements: bwa, seqtk, samtools, fqtrim, Biopython, Docker (optional).
Any restrictions to use by non-academics: None.
High-quality noise reads
Oxford Nanopore Technologies
Reads of interest
Lever MA, Torti A, Eickenbusch P, Michaud AB, Šantl-Temkiv T, Jørgensen BB. A modular method for the extraction of DNA and RNA, and the separation of DNA pools from diverse environmental sample types. Front Microbiol. 2015;6:1–25.
Barton HA, Taylor NM, Lubbers BR, Pemberton AC. DNA extraction from low-biomass carbonate rock: an improved method with reduced contamination and the low-biomass contaminant database. J Microbiol Methods. 2006;66:21–31.
Sabina J, Leamon JH. Bias in whole genome amplification: causes and considerations. Whole genome amplification. New York: Springer New York; 2015. p. 15–41.
Poretsky R, LM R-R, Luo C, Tsementzi D, Konstantinidis KT. Strengths and limitations of 16S rRNA gene amplicon sequencing in revealing temporal microbial community dynamics. PLoS One. 2014;9:e93827.
Mojarro A, Ruvkun G, Zuber MT, Carr CE. Nucleic acid extraction from synthetic Mars analog soils for in situ life detection. Astrobiology. 2017;17:747–60.
Carr CE, Mojarro A, Hachey J, Saboda K, Tani J, Bhattaru SA, et al. Towards in situ sequencing for life detection: IEEE Aerospace Conference. IEEE: Big Sky; 2017. pp. 1–18. https://doi.org/10.1109/AERO.2017.7943896.
Mojarro A, Hachey J, Bailey R, Brown M, Doebler R, Ruvkun G, et al. Nucleic acid extraction and sequencing from low-biomass synthetic Mars analog soils. In: 48th Lunar and Planetary Science Conference, vol. 48; 2017. p. 1585.
Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv. 2013:arXiv:1303.3997.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Li H. Seqtk: a toolkit for processing sequences in FASTA/Q formats. 2012. https://github.com/lh3/seqtk.
Morgulis A, Gertz EM, Schäffer AA, Agarwala R. A fast and symmetric DUST implementation to mask low-complexity DNA sequences. J Comput Biol. 2006;13:1028–40.
Pertea G. Fqtrim: v0. 9.4 release. 2015. https://doi.org/10.5281/zenodo.20552.
Fawcett T. An introduction to ROC analysis. Pattern Recogn Lett. 2006;27:861–74.
The authors would like to thank Michael Micorescu at Oxford Nanopore Technologies for providing and granting us permission to utilize his fastq quality filter script.
This work has been supported by NASA MatISSE award NNX15AF85G.
Availability of data and materials
The dataset analyzed during the current study is available from Figshare, https://doi.org/10.6084/m9.figshare.5868825.v1
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Mojarro, A., Hachey, J., Ruvkun, G. et al. CarrierSeq: a sequence analysis workflow for low-input nanopore sequencing. BMC Bioinformatics 19, 108 (2018). https://doi.org/10.1186/s12859-018-2124-3