DIAlignR provides precise retention time alignment across distant runs in DIA and targeted proteomics
Bibliographic record
Abstract
Abstract SWATH-MS has been widely used for proteomics analysis given its high-throughput and reproducibility but ensuring consistent quantification of analytes across large-scale studies of heterogeneous samples such as human-plasma remains challenging. Heterogeneity in large-scale studies can be caused by large time intervals between data-acquisition, acquisition by different operators or instruments, intermittent repair or replacement of parts, such as the liquid chromatography column, all of which affect retention time (RT) reproducibility and successively performance of SWATH-MS data analysis. Here, we present a novel algorithm for retention time alignment of SWATH-MS data based on direct alignment of raw MS2 chromatograms using a hybrid dynamic programming approach. The algorithm does not impose a chronological order of elution and allows for alignment of elution-order swapped peaks. Furthermore, allowing RT-mapping in a certain window around coarse global fit makes it robust against noise. On a manually validated dataset, this strategy outperforms the current state-of-the-art approaches. In addition, on a real-world clinical data, our approach outperforms global alignment methods by mapping 98% of peaks compared to 67% cumulatively and DIAlignR can reduce alignment error up to 30-fold for extremely distant runs. The robustness of technical parameters used in this pairwise alignment strategy has also been demonstrated. The source code is released under the BSD license at https://github.com/Roestlab/DIAlignR . Abbreviations: AUC Area Under the Curve DIA Data-independent acquisition LC Liquid chromatography LOESS Local weighted regression RSE Residual Standard Error RT Retention time XIC Extracted ion chromatograms Data Availability: Raw chromatograms and features extracted by OpenSWATH are available on PeptideAtlas. Servername: ftp.peptideatlas.org Username: PASS01280 Password: KQ2592b
Fetched live from OpenAlex and de-inverted. Abstracts are not stored in this database: the inverted indexes are 8.6 GB of the frame’s 9.3 GB of text, and the host has 13 GB free.
How this classification was reachedexpand
Full frame machine prediction
Teacher imitationNot calibrated prevalence, not ground truth. Human validation pending. The Gemma side is a direct model label for every work in the frame, read from the title-only record. The Codex side is a classifier learned from the 10,348 direct Codex labels and calibrated to design-weighted sample rates; fields without enough sample support carry no Codex call. Candidate is the union of the two sides; consensus is their intersection. These outputs are machine_predicted_unvalidated and are not human labels.
Distilled classifier scores by category (both heads)
| Category | Codex | Gemma |
|---|---|---|
| Metaresearch | 0.004 | 0.006 |
| Meta-epidemiology (narrow) | 0.001 | 0.001 |
| Meta-epidemiology (broad) | 0.001 | 0.001 |
| Bibliometrics | 0.002 | 0.002 |
| Science and technology studies | 0.001 | 0.001 |
| Scholarly communication | 0.002 | 0.001 |
| Open science | 0.001 | 0.002 |
| Research integrity | 0.001 | 0.002 |
| Insufficient payload (model declined to judge) | 0.003 | 0.004 |
Machine scores (provisional)
The two teacher heads of the student model, read on this work. A score orders the frame for review; it never asserts a category, and the validation status ships verbatim with every row.
Baseline scores from an immature model (maturity gate not passed, 7 training rounds). Scores rank; they never assert a category.
score_only:v0-immature-baseline · verbatim from the scoring run: score_only means the number may rank works, and no category label ships from itClassification
machine, unvalidatedMachine predicted; a candidate call from one source (direct Gemma or distilled Codex), not a consensus.
How this classification was reached, model by model and score by score, is at the end of the page under "How this classification was reached".