Software for plotting gene presence absence matrix

Question

[natasha]

Hi

I have a gene presence-absence matrix of accessory genes, for multiple strains of bacteria. I would like to produce a good figure with this data, and maybe include a tree as well . I was wondering what the best software programmes are available to do this? As I have >4000 genes, I don't each gene to be labelled!

Thanks

[Nicolas Rosewick]

how many strains do you have ?

[natasha]

50 strains. I make a mistake, I have 1800 genes to plot not 4000.

Answer 1

[Nicolas Rosewick]

Maybe try a PCA on your data. 4000 observations is quiet big to plot in a figure..

or you could do both (a. PCA ; b. heatmap)

Answer 2

[Jean-Karim Heriche]

Perhaps a heatmap would do ?

[WouterDeCoster]

For 4000 observations? I doubt whether that would send a clear message. Although a heatmap would be a nice combination with a hierarchical clustering tree.

You want to produce a good figure, but you should think first which message you want that figure to contain.

[natasha]

Because I am looking at bacteria of the the same ST, the accessory genes are what make each isolate unique from one another. Sorry I made a mistake -- I only have 1800 genes to plot (4000 is the core genome)

[WouterDeCoster]

That's still a lot. What is the message your figure should send? You could try as suggested a PCA with an additional hierarchical clustering (tree) and a PCA and see how that looks. There are many R packages that can do that for you.

Answer 3

[Philipp Bayer]

How about a gene presence/absence network, or project the presence/absence onto a known phylogeny, like Figures 2 and 3 in the CoPAP paper? http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3692100/