Question: Explaining PCA plots made from BigWig coverage files in DeepTools
0
gravatar for a.rex
16 months ago by
a.rex210
a.rex210 wrote:

I have 3 ATAC-seq conditions, with 2 replicates in each. I want to compare how similar they are to each other. I mapped with BWA and the output bam file was converted to bigwig format using bamcoverage, and normalising using RPKM.

I then used multibigwigsummary as following:

./multiBigwigSummary bins -b cellpop1_rep1 cellpop1_rep2 cellpop2_rep1 cellpop2_rep2 cellpop3_1 cellpop3_2 -out cell123.npz --binSize=100

The PCA plot of the resultant count summary looks like this: enter image description here

Am I right in thinking that my replicates cluster well together, meaning they are quite similar. However, there are differences in PC2. But the eigenvalues of this second component are minimal? Meaning that the differences are small between the samples, but they do exist?

atac-seq deeptools • 797 views
ADD COMMENTlink modified 16 months ago by RamRS26k • written 16 months ago by a.rex210

Rerun plotPCA with --transpose.

ADD REPLYlink written 16 months ago by Devon Ryan94k

what does this do? I don't understand what transpose does to the information?

ADD REPLYlink modified 16 months ago • written 16 months ago by a.rex210

It ... transposes the data, turning the rows into columns and columns into rows. In processes where the columns and rows are significantly different in meaning, such as when working with data.frames and PCA, this can make a significant difference in the meaning of the output, as opposed to when you're working with data structures such as 2D matrices of numbers.

ADD REPLYlink written 16 months ago by RamRS26k

Basically it does what Ram indicated. The gist is that at the moment PC1 is genomic position level changes, which will tend to be huge for ATAC data. That's going to end up masking what you're actually interested in, namely how well your samples actually cluster together. By transposing the matrix you look more at that (this is the standard in things like RNA-seq).

ADD REPLYlink written 16 months ago by Devon Ryan94k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1250 users visited in the last hour