How to add annotation information to a plant genome in the treatment of single-cell ATAC with Signac?
0
0
Entering edit mode
3.6 years ago
qzLi0031 • 0

I tried to use Signac to process plant single-cell ATAC data, but I encountered problems when adding gene annotation information to the Seurat object. Here is my code.

gtf <- rtracklayer::import("F:\\scATAC_scSeq\\Arabidopsis_thaliana.TAIR10.50.gtf")

gene.coords <- gtf[gtf$type == 'gene']

seqlevelsStyle(gene.coords) <- 'Ensembl'

gene.coords <- keepStandardChromosomes(gene.coords,pruning.mode = 'coarse')

Annotation(ATAC1) <- gene.coords

ATAC1 <- NucleosomeSignal(object = ATAC1)

ATAC1 <- TSSEnrichment(object = ATAC1,fast = FALSE)

enter image description here

annotation Signac plant • 678 views
ADD COMMENT

Login before adding your answer.

Traffic: 2815 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