Single-molecule real-time sequencing

(Redirected from SMRT sequencing)

Single-molecule real-time (SMRT) sequencing is a parallelized single molecule DNA sequencing method. Single-molecule real-time sequencing utilizes a zero-mode waveguide (ZMW).[1] A single DNA polymerase enzyme is affixed at the bottom of a ZMW with a single molecule of DNA as a template. The ZMW is a structure that creates an illuminated observation volume that is small enough to observe only a single nucleotide of DNA being incorporated by DNA polymerase. Each of the four DNA bases is attached to one of four different fluorescent dyes. When a nucleotide is incorporated by the DNA polymerase, the fluorescent tag is cleaved off and diffuses out of the observation area of the ZMW where its fluorescence is no longer observable. A detector detects the fluorescent signal of the nucleotide incorporation, and the base call is made according to the corresponding fluorescence of the dye.[2]

Technology

edit

The DNA sequencing is done on a chip that contains many ZMWs. Inside each ZMW, a single active DNA polymerase with a single molecule of single stranded DNA template is immobilized to the bottom through which light can penetrate and create a visualization chamber that allows monitoring of the activity of the DNA polymerase at a single molecule level. The signal from a phospho-linked nucleotide incorporated by the DNA polymerase is detected as the DNA synthesis proceeds which results in the DNA sequencing in real time.

Template preparation

edit

To prepare the library, DNA fragments are put into a circular form using hairpin adapter ligations.[3]

Phospholinked nucleotide

edit

For each of the nucleotide bases, there is a corresponding fluorescent dye molecule that enables the detector to identify the base being incorporated by the DNA polymerase as it performs the DNA synthesis. The fluorescent dye molecule is attached to the phosphate chain of the nucleotide. When the nucleotide is incorporated by the DNA polymerase, the fluorescent dye is cleaved off with the phosphate chain as a part of a natural DNA synthesis process during which a phosphodiester bond is created to elongate the DNA chain. The cleaved fluorescent dye molecule then diffuses out of the detection volume so that the fluorescent signal is no longer detected.[4]

Zero-Mode Waveguide

edit

The zero-mode waveguide (ZMW) is a nanophotonic confinement structure that consists of a circular hole in an aluminum cladding film deposited on a clear silica substrate.[5]

The ZMW holes are ~70 nm in diameter and ~100 nm in depth. Due to the behavior of light when it travels through a small aperture, the optical field decays exponentially inside the chamber.[6][7]

The observation volume within an illuminated ZMW is ~20 zeptoliters (20 X 10−21 liters). Within this volume, the activity of DNA polymerase incorporating a single nucleotide can be readily detected.[4][8]

Sequencing Performance

edit

Sequencing performance can be measured in read length, accuracy, and total throughput per experiment. PacBio sequencing systems using ZMWs have the advantage of long read lengths, although error rates are on the order of 5-15% and sample throughput is lower than Illumina sequencing platforms.[9]

On 19 Sep 2018, Pacific Biosciences [PacBio] released the Sequel 6.0 chemistry, synchronizing the chemistry version with the software version. Performance is contrasted for large-insert libraries with high molecular weight DNA versus shorter-insert libraries below ~15,000 bases in length. For larger templates average read lengths are up to 30,000 bases. For shorter-insert libraries, average read length are up to 100,000 bases while reading the same molecule in a circle several times. The latter shorter-insert libraries then yield up to 50 billion bases from a single SMRT Cell.[10]

History

edit

Pacific Biosciences (PacBio) commercialized SMRT sequencing in 2011,[11] after releasing a beta version of its RS instrument in late 2010.[12]

RS and RS II

edit
 
SMRT Cell for a RS or RS II Sequencer

At commercialization, read length had a normal distribution with a mean of about 1100 bases. A new chemistry kit released in early 2012 increased the sequencer's read length; an early customer of the chemistry cited mean read lengths of 2500 to 2900 bases.[13]

The XL chemistry kit released in late 2012 increased average read length to more than 4300 bases.[14][15]

On August 21, 2013, PacBio released a new DNA polymerase Binding Kit P4. This P4 enzyme has average read lengths of more than 4,300 bases when paired with the C2 sequencing chemistry and more than 5,000 bases when paired with the XL chemistry.[16] The enzyme’s accuracy is similar to C2, reaching QV50 between 30X and 40X coverage. The resulting P4 attributes provided higher-quality assemblies using fewer SMRT Cells and with improved variant calling.[16] When coupled with input DNA size selection (using an electrophoresis instrument such as BluePippin) yields average read length over 7 kilobases.[17]

On October 3, 2013, PacBio released new reagent combination for PacBio RS II, the P5 DNA polymerase with C3 chemistry (P5-C3). Together, they extend sequencing read lengths to an average of approximately 8,500 bases, with the longest reads exceeding 30,000 bases.[18] Throughput per SMRT cell is around 500 million bases demonstrated by sequencing results from the CHM1 cell line.[19]

On October 15, 2014, PacBio announced the release of new chemistry P6-C4 for the RS II system, which represents the company's 6th generation of polymerase and 4th generation chemistry--further extending the average read length to 10,000 - 15,000 bases, with the longest reads exceeding 40,000 bases. The throughput with the new chemistry was estimated between 500 million to 1 billion bases per SMRT Cell, depending on the sample being sequenced.[20][21] This was the final version of chemistry released for the RS instrument.

Throughput per experiment for the technology is both influenced by the read length of DNA molecules sequenced as well as total multiplex of a SMRT Cell. The prototype of the SMRT Cell contained about 3000 ZMW holes that allowed parallelized DNA sequencing. At commercialization, the SMRT Cells were each patterned with 150,000 ZMW holes that were read in two sets of 75,000.[22] In April 2013, the company released a new version of the sequencer called the "PacBio RS II" that uses all 150,000 ZMW holes concurrently, doubling the throughput per experiment.[23][24] The highest throughput mode in November 2013 used P5 binding, C3 chemistry, BluePippin size selection, and a PacBio RS II officially yielded 350 million bases per SMRT Cell though a human de novo data set released with the chemistry averaging 500 million bases per SMRT Cell. Throughput varies based on the type of sample being sequenced.[25] With the introduction of P6-C4 chemistry typical throughput per SMRT Cell increased to 500 million bases to 1 billion bases.

RS Performance
C1 C2 P4-XL P5-C3 P6-C4
Average read length bases 1100 2500 - 2900 4300 - 5000 8500 10,000 - 15,000
Throughput per SMRT Cell 30M - 40M 60M - 100M 250M - 300M 350M - 500M 500M - 1B

Sequel

edit
 
SMRT Cell for a Sequel Sequencer

In September 2015, the company announced the launch of a new sequencing instrument, the Sequel System, that increased capacity to 1 million ZMW holes.[26][27]

With the Sequel instrument initial read lengths were comparable to the RS, then later chemistry releases increased read length.

On January 23, 2017, the V2 chemistry was released. It increased average read lengths to between 10,000 and 18,000 bases.[28]

On March 8, 2018, the 2.1 chemistry was released. It increased average read length to 20,000 bases and half of all reads above 30,000 bases in length. Yield per SMRT Cell increased to 10 or 20 billion bases, for either large-insert libraries or shorter-insert (e.g. amplicon) libraries respectively.[29]

 
Pipette tip in an 8M SMRT Cell

On 19 September 2018, the company announced the Sequel 6.0 chemistry with average read lengths increased to 100,000 bases for shorter-insert libraries and 30,000 for longer-insert libraries. SMRT Cell yield increased up to 50 billion bases for shorter-insert libraries.[10]

Sequel Performance
V2 2.1 6.0
Average read length bases 10,000 - 18,000 20,000 - 30,000 30,000 - 100,000
Throughput per SMRT Cell 5B - 8B 10B - 20B 20B - 50B

8M Chip

edit

In April 2019 the company released a new SMRT Cell with eight million ZMWs,[30] increasing the expected throughput per SMRT Cell by a factor of eight.[31] Early access customers in March 2019 reported throughput over 58 customer run cells of 250 GB of raw yield per cell with templates about 15 kb in length, and 67.4 GB yield per cell with templates in higher weight molecules.[32] System performance is now reported in either high-molecular-weight continuous long reads or in pre-corrected HiFi (also known as Circular Consensus Sequence (CCS)) reads. For high-molecular-weight reads roughly half of all reads are longer than 50 kb in length.

Sequel II High-Molecular-Weight Performance
Early Access 1.0 2.0
Throughput per SMRT Cell ~67.4 GB Up to 160 GB Up to 200 GB

The HiFi performance includes corrected bases with quality above Phred score Q20, using repeated amplicon passes for correction. These take amplicons up to 20kb in length.

Sequel II HiFi Corrected Read Performance
Early Access 1.0 2.0
Raw reads per SMRT Cell ~250 GB Up to 360 GB Up to 500 GB
Corrected reads per SMRT Cell (>Q20) ~25 GB Up to 36 GB Up to 50 GB

Application

edit

Single-molecule real-time sequencing may be applicable for a broad range of genomics research.

For de novo genome sequencing, read lengths from the single-molecule real-time sequencing are comparable to or greater than that from the Sanger sequencing method based on dideoxynucleotide chain termination. The longer read length allows de novo genome sequencing and easier genome assemblies.[2][33][34] Scientists are also using single-molecule real-time sequencing in hybrid assemblies for de novo genomes to combine short-read sequence data with long-read sequence data.[35][36] In 2012, several peer-reviewed publications were released demonstrating the automated finishing of bacterial genomes,[37][38] including one paper that updated the Celera Assembler with a pipeline for genome finishing using long SMRT sequencing reads.[39] In 2013, scientists estimated that long-read sequencing could be used to fully assemble and finish the majority of bacterial and archaeal genomes.[40]

The same DNA molecule can be resequenced independently by creating the circular DNA template and utilizing a strand displacing enzyme that separates the newly synthesized DNA strand from the template.[41] In August 2012, scientists from the Broad Institute published an evaluation of SMRT sequencing for SNP calling.[42]

The dynamics of polymerase can indicate whether a base is methylated.[43] Scientists demonstrated the use of single-molecule real-time sequencing for detecting methylation and other base modifications.[44][45][46] In 2012 a team of scientists used SMRT sequencing to generate the full methylomes of six bacteria.[47] In November 2012, scientists published a report on genome-wide methylation of an outbreak strain of E. coli.[48]

Long reads make it possible to sequence full gene isoforms, including the 5' and 3' ends. This type of sequencing is useful to capture isoforms and splice variants.[49][50]

SMRT sequencing has several applications in reproductive medical genetics research when investigating families with suspected parental gonadal mosaicism. Long reads enable haplotype phasing in patients to investigate parent-of-origin of mutations. Deep sequencing enables determination of allele frequencies in sperm cells, of relevance for estimation of recurrence risk for future affected offspring.[51][52]

References

edit
  1. ^ Levene MJ, Korlach J, Turner SW, et al. (2003). "Zero-Mode Waveguides for Single-Molecule Analysis at High Concentrations". Science. 299 (5607): 682–6. Bibcode:2003Sci...299..682L. doi:10.1126/science.1079700. PMID 12560545. S2CID 6060239.
  2. ^ a b Eid J, Fehr A, Gray J, et al. (2009). "Real-Time DNA Sequencing from Single Polymerase Molecules". Science. 323 (5910): 133–8. Bibcode:2009Sci...323..133E. doi:10.1126/science.1162986. PMID 19023044. S2CID 54488479.
  3. ^ Friedmann, Theodore (2012). Advances in genetics (in Dutch). Oxford: Academic. ISBN 978-0-12-394395-8. OCLC 813987819.
  4. ^ a b "Pacific Biosciences Develops Transformative DNA Sequencing Technology" (PDF). Pacific Biosciences Technology Backgrounder. 2008.
  5. ^ Korlach J, Marks PJ, Cicero RL, et al. (2008). "Selective aluminum passivation for targeted immobilization of single DNA polymerase molecules in zero-mode waveguide nanostructures". PNAS. 105 (4): 1176–81. Bibcode:2008PNAS..105.1176K. doi:10.1073/pnas.0710982105. PMC 2234111. PMID 18216253.
  6. ^ Foquet M, Samiee KT, Kong X, et al. (2008). "Improved fabrication of zero-mode waveguides for single-molecule detection". J. Appl. Phys. 103 (3): 034301–034301–9. Bibcode:2008JAP...103c4301F. doi:10.1063/1.2831366. S2CID 38892226.
  7. ^ Zhu, Paul; Craighead, Harold G. (2012-06-09). "Zero-Mode Waveguides for Single-Molecule Analysis". Annual Review of Biophysics. 41 (1). Annual Reviews: 269–293. doi:10.1146/annurev-biophys-050511-102338. ISSN 1936-122X. PMID 22577821.
  8. ^ Baibakov, Mikhail; Barulin, Aleksandr; Roy, Prithu; Claude, Jean-Benoît; Patra, Satyajit; Wenger, Jérôme (1999-02-22). "Zero-mode waveguides can be made better: fluorescence enhancement with rectangular aluminum nanoapertures from the visible to the deep ultraviolet". Nanoscale Advances. 2 (9): 4153–4160. doi:10.1039/D0NA00366B. PMC 9417158. PMID 36132755.
  9. ^ Pollock, Jolinda; Glendinning, Laura; Wisedchanwet, Trong; Watson, Mick (2018). "The Madness of Microbiome: Attempting To Find Consensus "Best Practice" for 16S Microbiome Studies". Applied and Environmental Microbiology. 84 (7): e02627-17. Bibcode:2018ApEnM..84E2627P. doi:10.1128/AEM.02627-17. PMC 5861821. PMID 29427429.
  10. ^ a b "PacBio Post". Twitter. 19 Sep 2018.
  11. ^ Karow J (3 May 2011). "PacBio Ships First Two Commercial Systems; Order Backlog Grows to 44". GenomeWeb.
  12. ^ Karow J (7 Dec 2010). "PacBio Reveals Beta System Specs for RS; Says Commercial Release is on Track for First Half of 2011". GenomeWeb.
  13. ^ Karow J (10 Jan 2012). "After a Year of Testing, Two Early PacBio Customers Expect More Routine Use of RS Sequencer in 2012". GenomeWeb.
  14. ^ Heger M (13 Nov 2012). "PacBio's XL Chemistry Increases Read Lengths and Throughput; CSHL Tests the Tech on Rice Genome". GenomeWeb.
  15. ^ Heger M (5 Mar 2013). "PacBio Users Report Progress in Long Reads for Plant Genome Assembly, Tricky Regions of Human Genome". GenomeWeb.
  16. ^ a b "New DNA Polymerase P4 Delivers Higher-Quality Assemblies Using Fewer SMRT Cells". PacBio Blog. 21 Aug 2013.
  17. ^ lexnederbragt (19 Jun 2013). "Longing for the longest reads: PacBio and BluePippin". In between lines of code.
  18. ^ "New Chemistry for PacBio RS II Provides Average 8.5 kb Read Lengths for Complex Genome Studies". PacBio Blog. 3 Oct 2013.
  19. ^ Chaisson MJ, Huddleston J, Dennis MY, et al. (2014). "Resolving the complexity of the human genome using single-molecule sequencing". Nature. 517 (7536): 608–11. Bibcode:2015Natur.517..608C. doi:10.1038/nature13907. PMC 4317254. PMID 25383537.
  20. ^ "Pacific Biosciences Releases New DNA Sequencing Chemistry to Enhance Read Length and Accuracy for the Study of Human and Other Complex Genomes". Pacific Biosciences (Press Release). 15 Oct 2014.
  21. ^ "New Chemistry Boosts Average Read Length to 10 kb – 15 kb for PacBio RS II". PacBio Blog. 15 Oct 2014.
  22. ^ "SMRT Cells, sequencing reagent kits, and accessories for the PacBio RS II". Pacific Biosciences. 2020. Archived from the original on 2013-04-21. Retrieved 2012-04-28.
  23. ^ "PacBio Launches PacBio RS II Sequencer". Next Gen Seek. 11 Apr 2013. Archived from the original on 19 December 2019. Retrieved 18 April 2013.
  24. ^ "New Products: PacBio's RS II; Cufflinks". GenomeWeb. 16 Apr 2013.
  25. ^ "Duke Sequencing Post". Twitter. 30 Aug 2013.
  26. ^ "PacBio Announces Sequel Sequencing System". Bio-IT World. 30 Sep 2015. Archived from the original on 29 July 2020. Retrieved 16 November 2015.
  27. ^ Heger M (1 Oct 2015). "PacBio Launches Higher-Throughput, Lower-Cost Single-Molecule Sequencing System". GenomeWeb.
  28. ^ "New Chemistry and Software for Sequel System Improve Read Length, Lower Project Costs". PacBio Blog. 9 Jan 2017.
  29. ^ "New Software, Polymerase for Sequel System Boost Throughput and Affordability". PacBio Blog. 7 Mar 2018.
  30. ^ "PacBio Launches Sequel II System". Bio-IT World. 26 Apr 2019.
  31. ^ "Archived copy". Archived from the original on 2018-09-24. Retrieved 2018-09-24.{{cite web}}: CS1 maint: archived copy as title (link)
  32. ^ Heger M (7 Mar 2019). "PacBio Shares Early-Access Customer Experiences, New Applications for Sequel II". GenomeWeb.
  33. ^ Rasko DA, Webster DR, Sahl JW, et al. (2011). "Origins of the E. coli Strain Causing an Outbreak of Hemolytic–Uremic Syndrome in Germany". N. Engl. J. Med. 365 (8): 709–17. doi:10.1056/NEJMoa1106920. PMC 3168948. PMID 21793740.
  34. ^ Chin CS, Sorenson J, Harris JB, et al. (2011). "The Origin of the Haitian Cholera Outbreak Strain". N. Engl. J. Med. 364 (1): 33–42. doi:10.1056/NEJMoa1012928. PMC 3030187. PMID 21142692.
  35. ^ Gao H, Green SJ, Jafari N, et al. (2012). "Tech Tips: Next-Generation Sequencing". Genetic Engineering & Biotechnology News. 32 (8).
  36. ^ Schatz M (7 Sep 2011). "SMRT-assembly approaches" (PDF). schatzlab.cshl.edu (PacBio Users Meeting).
  37. ^ Ribeiro FJ, Przybylski D, Yin S, et al. (2012). "Finished bacterial genomes from shotgun sequence data". Genome Res. 22 (11): 2270–7. doi:10.1101/gr.141515.112. PMC 3483556. PMID 22829535.
  38. ^ Bashir A, Klammer A, Robins WP, et al. (2012). "A hybrid approach for the automated finishing of bacterial genomes". Nat. Biotechnol. 30 (7): 701–7. doi:10.1038/nbt.2288. PMC 3731737. PMID 22750883.
  39. ^ Koren S, Schatz MC, Walenz BP, et al. (2012). "Hybrid error correction and de novo assembly of single-molecule sequencing reads". Nat. Biotechnol. 30 (7): 693–700. doi:10.1038/nbt.2280. PMC 3707490. PMID 22750884.
  40. ^ Koren S, Harhay GP, Smith TP, et al. (2013). "Reducing assembly complexity of microbial genomes with single-molecule sequencing". Genome Biol. 14 (9): R101. arXiv:1304.3752. Bibcode:2013arXiv1304.3752K. doi:10.1186/gb-2013-14-9-r101. PMC 4053942. PMID 24034426.
  41. ^ Smith CC, Wang Q, Chin CS, et al. (2012). "Validation of ITD mutations in FLT3 as a therapeutic target in human acute myeloid leukaemia". Nature. 485 (7397): 260–3. Bibcode:2012Natur.485..260S. doi:10.1038/nature11016. PMC 3390926. PMID 22504184.
  42. ^ Carneiro MO, Russ C, Ross MG, et al. (2012). "Pacific Biosciences Sequencing Technology for Genotyping and Variation Discovery in Human Data". BMC Genom. 13 (1): 375. doi:10.1186/1471-2164-13-375. PMC 3443046. PMID 22863213.
  43. ^ Flusberg BA, Webster DR, Lee JH, et al. (2010). "Direct detection of DNA methylation during single-molecule, real-time sequencing". Nat. Methods. 7 (6): 461–5. doi:10.1038/nmeth.1459. PMC 2879396. PMID 20453866.
  44. ^ Clark TA, Murray IA, Morgan RD, et al. (2012). "Characterization of DNA Methyltransferase Specificities Using Single-Molecule, Real-Time DNA Sequencing". Nucleic Acids Res. 40 (4): e29. doi:10.1093/nar/gkr1146. PMC 3287169. PMID 22156058.
  45. ^ Song CX, Clark TA, Lu XY, et al. (2011). "Sensitive and Specific Single-Molecule Sequencing of 5-hydroxymethylcytosine". Nat Methods. 9 (1): 75–7. doi:10.1038/nmeth.1779. PMC 3646335. PMID 22101853.
  46. ^ Clark TA, Spittle KE, Turner SW, et al. (2011). "Direct Detection and Sequencing of Damaged DNA Bases". Genome Integr. 2 (1): 10. doi:10.1186/2041-9414-2-10. PMC 3264494. PMID 22185597.
  47. ^ Murray IA, Clark TA, Morgan RD, et al. (2012). "The Methylomes of Six Bacteria". Nucleic Acids Res. 40 (22): 11450–62. doi:10.1093/nar/gks891. PMC 3526280. PMID 23034806.
  48. ^ Fang G, Munera D, Friedman DI, et al. (2012). "Genome-wide Mapping of Methylated Adenine Residues in Pathogenic Escherichia Coli Using Single-Molecule Real-Time Sequencing". Nat. Biotechnol. 30 (12): 1232–9. doi:10.1038/nbt.2432. PMC 3879109. PMID 23138224.
  49. ^ Sharon D, Tilgner H, Grubert F, et al. (2013). "A Single-Molecule Long-Read Survey of the Human Transcriptome". Nat. Biotechnol. 31 (11): 1009–14. doi:10.1038/nbt.2705. PMC 4075632. PMID 24108091.
  50. ^ Au KF, Sebastiano V, Afshar PT, et al. (2013). "Characterization of the human ESC transcriptome by hybrid sequencing". PNAS. 110 (50): E4821–30. Bibcode:2013PNAS..110E4821A. doi:10.1073/pnas.1320101110. PMC 3864310. PMID 24282307.
  51. ^ Ardui S, Ameur A, Vermeesch JR, et al. (2018). "Single Molecule Real-Time (SMRT) Sequencing Comes of Age: Applications and Utilities for Medical Diagnostics". Nucleic Acids Res. 46 (5): 2159–68. doi:10.1093/nar/gky066. PMC 5861413. PMID 29401301.
  52. ^ Wilbe M, Gudmundsson S, Johansson J, et al. (2017). "A Novel Approach Using Long-Read Sequencing and ddPCR to Investigate Gonadal Mosaicism and Estimate Recurrence Risk in Two Families With Developmental Disorders". Prenatal Diagnosis. 37 (11): 1146–54. doi:10.1002/pd.5156. PMC 5725701. PMID 28921562.
edit