Circos Plot to Visualize Chromatin Accessibility Changes from ATAC-seq Data

Question

[biology_inform]

Hi everyone,

I have analyzed ATAC-seq data using the nf-core ATAC-seq pipeline, and I have obtained peak files from MACS2. My goal is to create a Circos plot to visualize global chromatin accessibility changes across different chromosomes between two sample groups (e.g., A549-WT vs A549-RES).

Despite exploring various tools, I haven't been able to figure out how to generate this plot. Could anyone provide a step-by-step guide or recommend an easy way to accomplish this?

Any help or pointers would be greatly appreciated!

Thank you!

[Pierre Lindenbaum]

Don't forget to follow up on your threads. If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one answer if they all work. If an answer was not really helpful or did not work, provide detailed feedback so others know not to use that answer.

• RNA-Seq analysis with Python
• FastQC command line usage
• clustering changes in complex heatmap everytime
• complex heatmap with CPM
• rna seq design matrix
• GSEA with duplicate sample
• GSEA phenotype label file and expression data relation
• Microarray platform download in R
• gene library composition
• wgcna without clinical traits
• ggplot and DESeq2 PCA plot
• finding coverage of bam file
• trimming a paired end sra
• igv and vcf difference
• featurecounts index bai file
• RNA-seq analysis without replicate
• edgeR comparison of 2 samples
• filtering dge data with cpm
• Salmon alignment based quantification
• Salmon installation problem
• igv viewer
• STAR aligment generating genome index
• Rna-seq data
• Transcription factors and gene expression profile

[Pierre Lindenbaum]

I haven't been able to figure out how to generate this plot.

show us what you tried and the error messages

[biology_inform]

I've been exploring different visualization tools like Circos and shinyCircos, and I've noticed that they require an input file similar to a BED file. However, I'm having trouble creating this file. I have two groups, each with two replicates, and the input file also requires an additional column called "value" (image below). I'm not sure how to generate this file correctly. Could you provide guidance or an example on how to structure this file?

[Pierre Lindenbaum]

Please do not paste screenshots of plain text content, it is counterproductive. You can copy paste the content directly here (using the code formatting option shown below), or use a GitHub Gist if the content volume exceeds allowed length here.

Answer 1

[Papyrus]

If you want to plot global accessibility, maybe you would like to plot the ATAC "signal", instead of specific peaks. In a very general manner: one strategy you could do is to generate bedGraph files from your BAM alignment files (e.g. using deepTools) and then get that information into R (with packages such as plyranges or rtracklayer). This will get you your "signal" tracks into GenomicRanges format. Then, you could use the R package circlize, for which there is ample documentation, to plot circular visualizations of these tracks.