When fold-enrichment is more informative than p-values: #TopAnat analysis of autism genes from GWAS

In a recent post on this blog we saw how to analyze results from a breast cancer GWAS. In that case, we did not have very strong expectations of tissue-specificity for the genes; it was more of an exploratory analysis.

This time, let’s do the same but searching GWAS Catalog for the term “autism”. Here we have a clear expectation of finding genes expressed in the brain.

Using the same methodology as for breast cancer, we find 87 genes with a significant SNP inside the gene, and we obtain the following in TopAnat:

gwas_topanat_autism1

We notice that the top terms, ranked by FDR (the default), are not necessarily brain related. Other the other hand, looking down the list we notice some higher fold-changes. Ranking by “Fold enrichment” (just click on the head of column), we get:

gwas_topanat_autism2

Now all the top terms are parts of the brain! In case you are wondering what the paraflocculus is, you can click on the Uberon ID and get to the definition: it’s a cerebellar
tonsil.

What this illustrates, and is well known in other contexts, is that the p-value (and its close friend the FDR) only takes you so far, since it is so dependent on sample size. The top terms by FDR have very large numbers of genes called expressed in them (≈20k!), whereas the more specific brain parts have a few thousand “only” genes expressed.

Thus we have a strong biological signal in the patterns of gene expression, but we have to be weary of relying on p-values. As always in statistics. 😉

 

Update: see also this more recent post on autism and epilepsy genes.

Advertisements

About bgeedb

This is the blog of the database of gene expression evolution Bgee, at University of Lausanne and Swiss Institute of Bioinformatics. http://www.bgee.org/
This entry was posted in topanat and tagged , , . Bookmark the permalink.

One Response to When fold-enrichment is more informative than p-values: #TopAnat analysis of autism genes from GWAS

  1. Pingback: Confirming that autism and epilepsy genes are expressed in specific brain areas | bgeedb

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s