locateVariants returns genes from both Forward and Reverse strands in PRECEDEID and FOLLOWID
まだ誰も着手していません。
評価
- 難易度
- 3/5
- 見積もり時間
- 1〜2日
- 初心者へのやさしさ
- 25/100
- issue の種類
- ドキュメント
- 明瞭さ
- 説明が足りない
- 活発さ
- 停滞
- 技術スタック
- r
調査の方向性
locateVariants() と IntergenicVariants(upstream=1000000, downstream=1000000) のドキュメントと実装から始め、次に PRECEDEID と FOLLOWID がどのように設定されるか、また strand がどのように扱われるかを確認します。strand を考慮した解釈が明確に文書化されているか、報告された挙動に焦点を絞った回帰テストがあり、合意された修正があることが完了の条件です。
索引モデルが issue の本文から書いたものです。
説明
Hello,
I'm looking at an intergenic variant (in the bovine genome) and want to figure the genes relative to which it is downstream or upstream, i.e. relative to the gene's position and its strand. I used locateVariants() with region=IntergenicVariants(upstream=1000000, downstream=1000000) in order to do this. However, some results puzzle me.
Here's what my variant looks like in the gene annotation results:
> all_var_df[rownames(all_var_df) == "AX-106756303", ]
seqnames start end width strand LOCATION LOCSTART LOCEND QUERYID TXID CDSID GENEID PRECEDEID FOLLOWID
AX-106756303 1 34617002 34617002 1 * intergenic NA NA 1 <NA> <NA> ENSBTAG0.... ENSBTAG0....
Here's what PRECEDEID looks like:
lapply(all_var_df[rownames(all_var_df) == "AX-106756303", ]$PRECEDEID, function(X) {mapIds( org.Bt.eg.db, keys=X, column="SYMBOL", keytype="ENSEMBL", multiVals="first") } )
ENSBTAG00000019313 ENSBTAG00000016711 ENSBTAG00000003877 ENSBTAG00000044714 ENSBTAG00000020940 ENSBTAG00000020939 ENSBTAG00000001656 ENSBTAG00000045788 ENSBTAG00000006536
"ZMIZ1" "PPIF" "ZCCHC24" NA "ANXA11" "PLAC9" "TMEM254" NA "CL46"
When looking closer at these genes in Ensembl, I notice that they are all located "to the right" of the SNP location on the forward strand and that some of them are on the Forward strand (e.g. ZMIZ1 and PPIF), while others are on the Reverse strand (e.g. ZCCHC24 and PLAC9):
This seems a bit confusing. Shouldn't the variant be considered upstream relative to the two genes on the Forward strand, and downstream relative to the genes in the Reverse strand?
How should the PRECEDEID gene list be interpreted, more precisely?
- 主要言語
- R
- スター
- 32
- フォーク
- 21
- PR マージ指標
- 30日以内にマージされた PR はありません
コントリビューションガイド
このリポジトリのコントリビューションガイドは索引されていません
はじめの一歩
- issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
- 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
- リポジトリをフォークし、ブランチを切って変更します。
- issue 番号を参照したプルリクエストを送ります。
Bioconductor/VariantAnnotation のほかの issue
-
refactor the package オープン
難易度 5/5 1週間以上 初心者へのやさしさ 25/100
Bioconductor/VariantAnnotation#115 · コメント 4 件 ·
-
Bioconductor/VariantAnnotation#114 · リアクション 1 件 · 担当者 2 名 ·
-
難易度 5/5 1週間以上 初心者へのやさしさ 35/100
Bioconductor/VariantAnnotation#113 · コメント 2 件 ·
-
難易度 3/5 1〜2日 初心者へのやさしさ 35/100
-
難易度 5/5 1週間以上 初心者へのやさしさ 30/100
Bioconductor/VariantAnnotation#88 · コメント 2 件 ·
Bioconductor/VariantAnnotation の issue をすべて見る
似ている issue
-
難易度 2/5 1〜3時間 初心者へのやさしさ 75/100
robjhyndman/forecast#1220 ·
-
難易度 2/5 1〜3時間 初心者へのやさしさ 65/100
JamesHWade/deputy#192 ·
-
難易度 2/5 1〜3時間 初心者へのやさしさ 75/100
-
bug triage_needed
難易度 2/5 1〜3時間 初心者へのやさしさ 75/100
-
難易度 2/5 1〜3時間 初心者へのやさしさ 72/100
pharmaverse/rtables#1123 · コメント 1 件 · リアクション 1 件 ·