Integrated genes from bulk RNA seq and ATAC seq
1
0
Entering edit mode
13 months ago
Chris ▴ 290

Hi all, I would like to do integrated the genes from bulk-RNA and genes from ATAC-seq as this paper did but don't know the code to do this. For bulk RNA seq, I got around 300 up and down-regulated genes between 2 conditions but ATAC seq got around 10k genes that fold change > 2.5.

Would you please have a suggestion? Thank you so much.

https://www.frontiersin.org/articles/10.3389/fnut.2021.742672/full

enter image description here

ATAC-seq bulk-RNA • 2.0k views
ADD COMMENT
1
Entering edit mode

hi Chris - what is the p-value that you're using for the cutoff when looking at your ATAC data? It looks like they used the following: different peaks were filtered with P < 0.00001 |log2(fold change)| >= 1.

ADD REPLY
0
Entering edit mode

Hi Sasha, if I use log2 fold > 2.5 and p < 0.00001, my object from Diffbind has around 20k observations.

ADD REPLY
1
Entering edit mode

Hmm - maybe there is some confusion in the phrasing of what they are communicating and how we are interpreting it. It may be that they also had 20K observable peaks but the peaks happened in the area of specific genes. What is that number when you look at it from a gene point of view.

ADD REPLY
0
Entering edit mode

With the filtering values above, there are around 8k unique genes in my data. My question is the code to get the overlap genes as they did.

ADD REPLY
1
Entering edit mode

One thing that you could potentially try is bedtools intersect to find the overlapping regions between the ATAC-seq peaks and the RNA-seq gene annotations. Here's an example command:

bedtools intersect -a atac_peaks.bed -b rna_genes.gtf -wa -wb > overlap_regions.bed

You can read more about it here: https://bedtools.readthedocs.io/en/latest/content/tools/intersect.html

DISCLAIMER: I'm using my chatbot here (https://tinybio.cloud) to help generate this answer. This answer has not been tested and may be incorrect. You can download it from the website. If this answer does not work - please let me know in this thread!

ADD REPLY
0
Entering edit mode

Thanks Sasha for the suggestion. Interesting chatbot. Don't know why I didn't get notified about your reply. Maybe -b rna.bam because rna file output doesn't have a gtf file.

ADD REPLY
0
Entering edit mode

I have 4 bam files for 2 conditions, 2 replicates for each condition. Do you think I should merge them into one bam file to run the bedtools command?

ADD REPLY
1
Entering edit mode

Hi Chris - if you could post more information describing the ATAC and RNA separately that would be great.

But, given what I understand; maybe try to merge the gff files and the bed files first and then do an intersect. Other people may have better ideas.

ADD REPLY
0
Entering edit mode

Seem I don't have gff file in both RNA-seq and ATAC seq analysis. Would you please tell me what information about ATAC and RNA-seq you want me to add?

ADD REPLY
1
Entering edit mode

To do the intersect using bedtools intersect on ATAC and RNA data, you would need BED files representing the genomic regions of interest for both data types. For ATAC-seq data, you would typically have a BED file containing the genomic coordinates of open chromatin regions, while for RNA-seq data, you might have a BED file containing the genomic coordinates of expressed genes or exons.

See if that works for ya.

Here's a very basic example:

atac_bed_file="atac_regions.bed"

# RNA-seq data in a BED file
rna_bed_file="rna_expressed_genes.bed"

# Perform the intersection and save the results to a new BED file
bedtools intersect -a $atac_bed_file -b $rna_bed_file > overlapping_regions.bed
ADD REPLY
1
Entering edit mode

Thank you for your help. I found this tutorial mentioning the bed file from ATAC seq. Before this, I haven't seen the bed file from ATAC-seq analysis yet, only bigwig, bam, broadpeak files. https://cambiotraining.github.io/chipseq/Practicals/ATACseq/ATACseq_tutorial.html

ADD REPLY
1
Entering edit mode

You might get a way with just using the bams. You can see the options listed here for bedtools intersect: https://bedtools.readthedocs.io/en/latest/content/tools/intersect.html

Look at the -bed flag.

Hope this helps!

ADD REPLY

Login before adding your answer.

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