Megahertz Pangenome Documentation
website·Indexed 30 days ago
https://megahertz-docs.netlify.appTechnical documentation for pangenome sequence-to-graph (S2G) acceleration, featuring detailed notes on hardware architectures and algorithms including Locus (Reconfigurable PE Array, Graph Path Finder, Traceback Engine), PanGaea, TALCO (WFAA and X-Drop accelerators), SeGraM (BitAlign), and GenASM-DC.
What's inside Megahertz
- Banding is a strategy used in sequence alignment to avoid filling the entire $m \times n$ DP matrix. Instead, it calculates only the cells within a 'band' near the main diagonal, based on the assumption that the optimal alignment typically does not deviate significantly from the main diagonal. This reduces time complexity from $O(mn)$ to $O(nw)$, where $w$ is the bandwidth.
Overview of SeGraM Pre-processing Workflow
SeGraM requires two primary pre-processing steps before it can execute read mapping. These steps are performed once and the resulting static data structures are pre-loaded into main memory to be reused across multiple mapping executions, amortizing the overhead.
- Graph-Based Reference Generation: Creating a genomic graph from linear references and variations.
- Hash-Table-Based Index Generation: Creating an index for the reference graph to facilitate fast seeding.
Overview of the Northeast Semiconductor R&D Corridor
The 'Northeast Tech Corridor' spanning New York, Massachusetts, and New Jersey is the highest density research area for semiconductors in the US. It is characterized by a strong synergy between top-tier academic institutions and industrial R&D centers. Key focus areas include analog, RF, and connectivity chips, as well as cutting-edge semiconductor innovation and materials research.Overview of the California Semiconductor Cluster
California (CA) serves as the global hub for semiconductor design, software, IP, and equipment. The ecosystem is primarily centered around the axis of San Jose, Santa Clara, and San Diego, extending to Irvine, Los Angeles, and Seattle. This region forms a "computing innovation corridor" that dominates the development of GPUs, AI, mobile communications, and server SoCs, hosting the majority of US Fabless companies.Overview of Partial Order Alignment (POA) for genome polishing
Partial Order Alignment (POA) is a dynamic programming algorithm used during the genome assembly polishing phase to correct errors in de novo assembly. Unlike traditional sequence alignment that treats the reference as a linear sequence, POA represents the reference as a Partial Order Graph (a Directed Acyclic Graph or DAG). New reads are aligned to this graph to update node and edge weights, eventually producing a more accurate consensus sequence based on the highest-weight path.Overview of the Texas Semiconductor Ecosystem
Texas serves as a critical hub for Integrated Device Manufacturers (IDM), Microcontroller Units (MCU), and automotive electronics. Unlike California's fabless-centric model, Texas's strength is driven by a combination of industrial demand and manufacturing roots, specifically in automotive-grade power chips, smart energy, and embedded applications for industrial, oil, and power sectors. Austin has evolved into a primary innovation center for AI SoC, MCU, DSP, and automotive semiconductors.Overview of the Graph Path Finder (GPF) Engine
The Graph Path Finder (GPF) engine is a hardware module in the Locus accelerator dedicated to the Graph Chaining (GC) phase of S2G mapping. Its primary purpose is to handle shortest-path searches when linear chains need to extend across segment boundaries. It determines if a source segment can reach a target segment, calculates the shortest graph distance for scoring transitions, and identifies intermediate vertices for the 'bridging' phase to generate a complete vertex-by-vertex path for subsequent alignment.Overview of the TALCO tiling strategy
TALCO (Tiling Genome Sequence Alignment using Convergence of Traceback Pointers) is a tiling strategy designed for long genomic alignments. Unlike traditional tiling heuristics, TALCO guarantees optimality under banding constraints, meaning alignment scores are not degraded by the tiling process. It leverages the insight that traceback pointers in alignments have weak long-range dependencies, leading to the convergence of multiple alignment paths.Overview of the Sequence-to-Graph Mapping Pipeline
The sequence-to-graph mapping pipeline follows a seed-and-extend strategy consisting of two offline pre-processing steps and three online execution steps:
Offline Pre-processing:
- Graph Construction: Builds the genome graph using a linear reference genome and known variations.
- Indexing: Indexes graph nodes and generates a hash-table-based index for fast lookup.
Online Execution:
- Seeding: Fragments reads into sub-strings (seeds) and identifies candidate mapping locations within the graph nodes using the index.
- Filtering/Chaining/Clustering (Optional): Reduces the number of required alignments for the final step.
- Alignment: Performs alignment between the query read and the remaining candidate subgraphs to find the optimal alignment.
Overview of Banded Smith-Waterman (BSW) for Read Alignment
Banded Smith-Waterman (BSW) is a variation of the classic Smith-Waterman (SW) local sequence alignment algorithm. It is primarily used in short-read sequencing to align a read (typically 100–150 bp) to a reference genome. BSW identifies the most similar local region between two sequences while accounting for mismatches, insertions, and deletions. To reduce computational complexity, BSW restricts calculations to a 'band' around the main diagonal of the DP table, assuming that reasonable alignments will not deviate significantly from the diagonal.Overview of Genomics Analysis Pipelines
Genome sequencing analysis typically follows one of three primary pipelines after basecalling (the process of deriving reads from raw sequencer signals):
- Reference-guided Assembly: Reconstructs a sample genome by aligning reads to an existing reference genome using dynamic programming to estimate pairwise similarity. This is followed by variant calling (using HMMs or ML models) to detect mutations.
- De novo Assembly: Constructs a genome from scratch without a reference. Reads with overlapping regions are chained into an overlap graph, extended into contiguous regions, and refined via a graph-based dynamic programming polishing step.
- Metagenomics Classification: Used for real-time pathogen detection and microbial abundance estimation by aligning microbial reads to a reference pan-genome consisting of different species.
Overview of the GenASM acceleration framework
GenASM is a high-performance, low-power acceleration framework for Approximate String Matching (ASM) used in genome sequence analysis. It enhances the Bitap algorithm to support both short and long reads and eliminates loop-carried data dependencies to enable parallelization of single string matching operations. The framework consists of two primary hardware components:
- GenASM-DC (Distance Calculation): Executes the modified Bitap algorithm to generate bitvectors (representing matches, insertions, deletions, or substitutions) and calculates the minimum number of errors (edit distance) between a read and a reference segment.
- GenASM-TB (Traceback): Executes a novel Bitap-compatible traceback algorithm using the bitvectors from GenASM-DC to identify the optimal alignment of a read.