Module 2: Wet-Lab Quality Control, Library Metrics & Oncology Panel Design
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 2 of 12
0% Course Completed
1. Library Preparation Chemistries: Amplicon vs. Hybrid Capture
Clinical NGS workflow success is governed by the library preparation method chosen for targeted enrichment:
| Enrichment Technology | Biochemical Principle | Input DNA Tolerance | Strengths | Diagnostic Limitations |
|---|---|---|---|---|
| Multiplex Amplicon PCR | Dozens to hundreds of target-specific primer pairs amplify regions of interest simultaneously. | Ultra-low: $1-10\text{ ng}$ degraded FFPE DNA | Rapid protocol (3–5 hours), high on-target percentage (> 90%), simple workflow. | High PCR duplication artifacts, vulnerable to primer drop-out caused by SNPs under primer binding sites, poor performance on copy number variation (CNV). |
| Hybridization Capture | Biotinylated oligonucleotide probes (RNA/DNA baits) hybridize in solution to genomic shotgun libraries; captured using streptavidin magnetic beads. | Standard: $50-200\text{ ng}$ (down to 20 ng with enzymatic shearing) | Superior target uniformity, robust CNV and fusion detection, captures novel breakpoints, compatible with Unique Molecular Identifiers (UMIs). | Longer hybridization time (16–24 hours), requires higher sequencing depth to overcome non-specific off-target capture. |
2. FFPE Tissue Processing & Chemical Artifacts
Formalin-Fixed Paraffin-Embedded (FFPE) tissue remains the global standard for anatomical pathology preservation. However, formalin fixation introduces severe biochemical damage that distorts sequencing data:
- Cytosine Deamination: Formalin accelerates the hydrolytic deamination of cytosine into uracil ($C \to U$). During PCR amplification, DNA polymerase reads uracil as thymine, introducing artificial $C:G \to T:A$ transition errors. In low-VAF somatic calls, these deamination artifacts mimic true somatic mutations. Pre-treatment with Uracil-DNA Glycosylase (UDG) cleaves artificial uracil bases to suppress false positives.
- Formaldehyde Cross-Linking: Protein-DNA and DNA-DNA methylene cross-links stall polymerase enzymes, causing chimeric read artifacts and artificial short fragment lengths.
- DNA Fragmentation: Severe hydrolysis degrades FFPE DNA into fragments < 150 bp, requiring quantification via the DIN (DNA Integrity Number) or $\Delta C_q$ functional qPCR assays.
3. Clinical Run & Library QC Summary Interpretation
Prior to clinical bioinformatics processing, a molecular pathologist must evaluate the vendor sequencing metrics against pre-defined Clinical Laboratory Improvement Amendments (CLIA) and CAP acceptance criteria:
| Quality Metric | Clinical Passing Threshold | Root Cause if Failed | Recommended Action |
|---|---|---|---|
| Mean Target Depth | $\ge 500\times$ (FFPE) / $\ge 2,500\times$ (cfDNA) | Under-clustering on flowcell, library over-dilution, poor hybridization capture. | If coverage is insufficient, low-VAF subclonal mutations (5% VAF) will be missed; re-sequence or re-pool library. |
| Uniformity of Coverage (% > 0.2× mean) | $\ge 90\% - 95\%$ of targeted bases | Extreme GC-rich or AT-rich target regions, probe drop-out, poor PCR amplification kinetics. | Panel gap: report affected un-covered exons with explicit diagnostic caveats. |
| Percent On-Target Reads | $\ge 65\% - 80\%$ (hybrid capture) | Insufficient washing stringency, low hybridization temperature, non-specific bait binding. | Acceptable if total depth threshold is satisfied; otherwise repeat capture with optimized wash buffers. |
| Duplicate Read Rate | $< 30\%$ for $50\text{ ng}$ input | Low molecular library complexity, excessive PCR amplification cycles ($> 18$ cycles). | High duplication artificially inflates depth; deduct duplicates bioinformatically using Picard or UMIs. |
| Phred Quality Score ($Q_{30}$) | $\ge 80\% - 85\%$ of bases $\ge Q_{30}$ | Sequencing reagent degradation, flowcell bubble, optical laser calibration drift. | Base-call error probability $> 0.1\%$; re-run flowcell if $Q_{30}$ drops below clinical threshold. |
4. Hands-On Exercise: Designing a 10-Gene Targeted Lung Cancer Panel
Using NCCN Non-Small Cell Lung Cancer guidelines and OncoKB Level 1/2 therapeutic evidence, compile a minimal 10-gene diagnostic lung panel:
Core Actionable Genes for NSCLC Panel:
1. EGFR -> Exons 18-21 (Sensitizing L858R, Exon 19 del, Exon 20 ins; Resistance T790M, C797S)
2. ALK -> Exons 19-29 (Fusions with EML4; Resistance mutations L1196M, G1202R)
3. ROS1 -> Exons 32-41 (Fusions with CD74, SDC4; Crizotinib/Entrectinib sensitivity)
4. BRAF -> Exon 15 (V600E point mutation; Dabrafenib + Trametinib combination)
5. KRAS -> Exons 2-3 (G12C mutation for Sotorasib/Adagrasib; non-G12C drivers)
6. MET -> Exon 14 and flanking splice sites (Exon 14 skipping; MET amplification)
7. RET -> Exons 11-13 (Fusions with KIF5B, CCDC6; Selpercatinib/Pralsetinib)
8. ERBB2 -> Exon 20 insertions (Trastuzumab deruxtecan eligibility)
9. NTRK1-3 -> Kinase domain fusions (Pan-cancer Larotrectinib/Entrectinib sensitivity)
10. TP53 -> Exons 2-11 (Prognostic co-mutation indicator of genomic instability)
⚠️
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
chr1or1; divergent naming causes silent alignment drops.
Was this guide helpful to your research?
Help us maintain 100% accurate, peer-reviewed bioinformatics standards.