Gene Tracks from 10x Data
1
0
Entering edit mode
9 weeks ago
ntsopoul ▴ 10

hi, I am interested in the 3’ UTR of genes and would like to see the gene tracks of certain populations. How can I extract from 10x data e.g. bam files the specific cell population and make bigwig files out of this subset? Is there a tool where I can inspect a bam file, select a population from reduced dim plots and visualize the gene tracks? I could not find such a function with the loupe browser for rna-seq. However the loupe browser can do it for atac-seq …

scRNA-seq polyadenylation 10x • 319 views
ADD COMMENT
0
Entering edit mode
9 weeks ago

The bams that 10xGenomics makes are tagged with cell barcode. So figure out which call barcodes you want, then you can filter the bam accordingly.

I don't know off the top of my head if you can get cell barcodes out of Loupe. You might need to use something like Seurat.

ADD COMMENT
0
Entering edit mode

Thanks a lot! I thought already that something similar might be the case... So by filtering you mean "grep" the cell ID of interest from the .bam file into a new .bam file

It is just sad that there is no tool to inspect the tracks of cells. It will be a pain to do this ...

ADD REPLY
0
Entering edit mode

Yeah, pretty much. Possibly Picardtools has something for filtering by tags?

ADD REPLY
0
Entering edit mode

IGV

right click reads -> Experiment Type -> RNA

this will split reads out by the "CB" tag for you.

This won't generate .wig files for you; but filtering by a read tag (instead of a read FLAG) would be a great feature request for deeptools (https://github.com/deeptools/deepTools)

ADD REPLY

Login before adding your answer.

Traffic: 2528 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6