pileUp and dmr showing significant results for m6A even with very low amount of mapped reads from Dorado's modBAM
まだ誰も着手していません。
評価
- 難易度
- 4/5
- 見積もり時間
- 3〜5日
- 初心者へのやさしさ
- 35/100
- issue の種類
- バグ
- 明瞭さ
- 説明が足りない
- 活発さ
- 静か
- 技術スタック
- rust
調査の方向性
Start with the modkit pileup and dmr entry points, using the reported modBAM inputs and the 0h, 24h, and 48h mapped-read counts. Reproduce the pass-threshold and missing-in-one-condition warnings, then compare the dmr bedMethyl results across read depths. Done means establishing whether the low-input 0h results are statistically valid and documenting or correcting the behavior.
索引モデルが issue の本文から書いたものです。
説明
Dear modkit developers,
I've been facing this curious issue with m6A pileup) and differential methylation (dmr), and would like to hear your opinion on whether or not to trust on results:
I'm analysing m6A on a viral RNAs in a timeseries designh (0h, 24h, and 48h, of infection). On the 0h timepoint, expectedly, I've got very few reads aligned to the virus ref genome ( ~3k mapped reads only).
I know that modkit pileup recommends an N sampling of 10-50k reads for performing that "pass" threshold calculation (and I use -n 50000 on my runs), but it does still seem to run well even for these 0h samples like 0h_rep1 (~2.1k reads mapped) and 0h_rep3 (~4.8k reads mapped), finding a confidence 'pass' threshold over 80%.
Then, modkit dmr also seems to run well (on both 0h-24h and 0h-vs-48h), even though warning on a great amount of failed sites due to "missing-in-one-condition", which is the 0h one.
After filtering for pval < 0.05 (column 15 from the dmr bedMethyl), I get tens of differentially modified m6A sites, with an effect size (column) bias towards the 0h sample group. All sites are highly differentially modified in the 0h group.
Shall I trust these results? Or it's just a statistical bias from the little input size (%mapped reads) of the 0h sample.
JFYI: 24h and 48h have several hundreds thousands (200-700k) mapped reds against the same genome. And it's dmr bedMethyl file displays a good balance of effect sizes (some more modified in the 24h, and some more in the 48h).
Thanks,
Best,
Elton
- 主要言語
- Rust
- スター
- 276
- フォーク
- 33
- PR マージ指標
- 30日以内にマージされた PR はありません
環境構築
このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
nanoporetech/modkit のほかの issue
-
bug
難易度 2/5 1〜3時間 初心者へのやさしさ 68/100
nanoporetech/modkit#520 · コメント 2 件 ·
-
documentation
難易度 1/5 1時間未満 初心者へのやさしさ 65/100
nanoporetech/modkit#336 · コメント 1 件 ·
-
難易度 5/5 1週間以上 初心者へのやさしさ 25/100
nanoporetech/modkit#723 ·
-
難易度 3/5 1〜2日 初心者へのやさしさ 66/100
nanoporetech/modkit#721 · コメント 1 件 ·
-
難易度 4/5 3〜5日 初心者へのやさしさ 55/100
nanoporetech/modkit#719 ·
nanoporetech/modkit の issue をすべて見る
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 84/100
vercel-labs/agent-browser#2017 ·
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
tursodatabase/turso#9405 ·
メンテナーはふだん 1 日以内に返信
-
bug
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
PolyMeilex/Neothesia#447 ·
メンテナーはふだん 1 日以内に返信
-
backend::vllm diffusion multimodal
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
メンテナーはふだん 1 日以内に返信
-
難易度 2/5 1〜3時間 初心者へのやさしさ 68/100
trezor/trezor-firmware#7985 ·
メンテナーはふだん 2 日以内に返信