omics.co.in logo omics.co.in

Module 1: Long-Read Sequencing Chemistries: Oxford Nanopore (ONT) vs. PacBio SMRT

Tested Stacks: GATK 4.6.1 Samtools 1.21 Nextflow 24.10 Seurat v5.1 Python 3.12
Prerequisites: Basic Linux CLI & Sequencing Formats
📚 Chapter 1 of 12 0% Course Completed

1. Biophysical Principles of Third-Generation Sequencing

First-generation capillary Sanger sequencing and second-generation short-read sequencing-by-synthesis (Illumina, MGI) rely on ensemble fluorescent measurements across clonal DNA clusters. While achieving extraordinary throughput and low per-base error rates ($< 0.1\%$), short reads ($100 - 300\text{ bp}$) cannot span complex repetitive elements, segmental duplications, or large structural variants.

Third-generation single-molecule sequencing resolves native, unamplified DNA and RNA strands. Two dominant platforms define this paradigm:

Biophysical ParameterOxford Nanopore Technologies (ONT)Pacific Biosciences (PacBio SMRT)
Core Detection MechanismResistive pulse sensing across an engineered protein pore (CsgG/R10.4.2) embedded in a synthetic electro-resistant polymer membrane.Optical zero-mode waveguide (ZMW) nanostructure monitoring processive incorporation of fluorescently labeled hexaphosphate dNTPs by $\Phi29$ DNA polymerase.
Native Signal MeasuredIonic current disruptions ($\Delta I$ in picoamperes, $\text{pA}$) caused by the steric and electrostatic occlusion of a 5-6 nucleotide k-mer in the pore barrel.Fluorescent light pulses emitted in zeptoliter detection volumes ($20 \times 10^{-21}\text{ L}$); base identity determined by emission wavelength, kinetics by inter-pulse duration (IPD).
Raw Read LengthYields reads proportional to intact input DNA fragment length. Standard runs: $10 - 50\text{ kb}$; "Ultra-long" protocols: up to $4.2\text{ Mb}$.Single continuous template length $10 - 25\text{ kb}$. Reads pass multiple times over a closed circular dumbbell template (SMRTbell).
Accuracy ParadigmRaw single-pass accuracy: $98.5 - 99.2\%$ with R10.4.2 chemistry and Super-Accurate (SUP) basecalling models; Duplex mode achieves $\text{Q}30+$ ($>99.9\%$).Circular Consensus Sequencing (CCS) passes over the same template $\ge 10-15$ times, yielding High-Fidelity (HiFi) reads exceeding $\text{Q}30$ ($> 99.9\%$).
Epigenetic ModificationsDirect chemical detection of $5\text{mC}$, $5\text{hmC}$, $4\text{mC}$, and $6\text{mA}$ from native current disruptions without chemical bisulfite conversion or enzyme deamination.Kinetic signature detection: Polymerase pausing or altered Inter-Pulse Duration (IPD) indicates methylated or damaged bases.

2. The Long-Read Imperative: Clinical & Genomic Scenarios

Long-read sequencing has transformed human genomics by resolving genomic regions previously deemed "dark" or inaccessible to short-read NGS:

  • Structural Variation & Balanced Rearrangements: Over $70\%$ of human structural variations ($> 50\text{ bp}$) cannot be discovered by short reads because the breakpoints fall inside repetitive or low-complexity sequence. Long reads span the entirety of deletions, insertions, inversions, and translocations in single continuous reads.
  • Segmental Duplications & Pseudogenes: High-sequence-identity duplications (e.g. CYP2D6 vs. CYP2D7, SMN1 vs. SMN2) cause mapping ambiguities in short reads. Long reads anchor uniquely into flanking divergence sites.
  • Tandem Repeat Expansions: Pathogenic trinucleotide and hexanucleotide expansions (e.g. C9orf72 in ALS/FTD, HTT in Huntington disease, FMR1 in Fragile X) exceed $1-10\text{ kb}$ in length and cannot be sized or phased with short-read sequencing.
  • Chromosome-Scale Haplotype Phasing: Long reads connect heterozygous single-nucleotide polymorphisms (SNPs) across tens of kilobases, allowing complete reconstruction of maternal and paternal haplotypes without parental sequencing.

3. Hands-On Benchmark Exercises

Calculate the theoretical physical length of a $100\text{ kb}$ Oxford Nanopore read and determine the translocation time required at standard motor enzyme speed ($400\text{ bases/second}$):

# Calculate translocation duration in seconds and minutes
python3 -c "
read_length_bp = 100000
translocation_speed = 400 # bp/sec
duration_sec = read_length_bp / translocation_speed
duration_min = duration_sec / 60
physical_length_um = read_length_bp * 0.34e-3 # 0.34 nm per base pair in B-DNA

print(f'Physical length: {physical_length_um:.2f} micrometers')
print(f'Translocation duration: {duration_sec:.1f} seconds ({duration_min:.2f} minutes)')
"
⚠️

Common Pitfalls & Troubleshooting Checklist

  • Java Heap Exhaustion (GATK / Picard): If processes terminate with java.lang.OutOfMemoryError, set your environment variable explicitly: export _JAVA_OPTIONS="-Xmx16g -Djava.io.tmpdir=/tmp".
  • Index Timestamp Mismatch: If tools complain that a BAM index is missing or out of date, regenerate with samtools index -@ 4 sample.bam.
  • Conda Channel Conflicts: Avoid ABI incompatibilities between bioconda and conda-forge by setting strict priority: conda config --set channel_priority strict.
  • Chromosome Naming Contamination: Always verify whether reference FASTA headers use chr1 or 1; divergent naming causes silent alignment drops.
Was this guide helpful to your research?
Help us maintain 100% accurate, peer-reviewed bioinformatics standards.
📝 Propose Edit on GitHub ↗