locateVariants returns genes from both Forward and Reverse strands in PRECEDEID and FOLLOWID
Nadie ha tomado este issue todavía.
Evaluación
- Dificultad
- 3/5
- Tiempo estimado
- 1-2 días
- Aptitud para principiantes
- 25/100
- Tipo de issue
- Documentación
- Claridad
- Necesita aclaración
- Estado de actividad
- Estancado
- Stack tecnológico
- r
- Área
- bioinformatics
Línea de trabajo
Empieza por la documentación y la implementación de locateVariants() e IntergenicVariants(upstream=1000000, downstream=1000000); después, inspecciona cómo se rellenan PRECEDEID y FOLLOWID y cómo se gestiona la hebra. Se considera terminado cuando la interpretación consciente de la hebra está documentada con claridad, o cuando el comportamiento informado tiene una prueba de regresión específica y existe una corrección acordada.
Escrito por el modelo de indexación a partir del texto del issue.
Descripción
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?
- Lenguaje dominante
- R
- Estrellas
- 32
- Forks
- 21
- Métricas de merge de PR
- Sin PR fusionados en 30 d
Guía de contribución
No hay ninguna guía de contribución indexada para este repositorio
Primeros pasos
- Lee el issue completo y luego la guía de contribución del proyecto.
- Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
- Haz un fork del repositorio y trabaja en una rama.
- Abre un pull request que haga referencia al número del issue.
Más de Bioconductor/VariantAnnotation
-
refactor the package Abierto
Dificultad 5/5 Más de una semana Aptitud para principiantes 25/100
Bioconductor/VariantAnnotation#115 · 4 comentarios ·
-
Bioconductor/VariantAnnotation#114 · 1 reacción · 2 asignados ·
-
Dificultad 5/5 Más de una semana Aptitud para principiantes 35/100
Bioconductor/VariantAnnotation#113 · 2 comentarios ·
-
Dificultad 3/5 1-2 días Aptitud para principiantes 35/100
-
ensemblVEP::parseCSQToGRanges Abierto
Dificultad 5/5 Más de una semana Aptitud para principiantes 30/100
Bioconductor/VariantAnnotation#88 · 2 comentarios ·
Todos los issues de Bioconductor/VariantAnnotation
Issues similares
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 75/100
robjhyndman/forecast#1220 ·
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 65/100
JamesHWade/deputy#192 ·
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 75/100
-
bug triage_needed
Dificultad 2/5 1-3 horas Aptitud para principiantes 75/100
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 72/100
pharmaverse/rtables#1123 · 1 comentario · 1 reacción ·