# Question about the practical relevance of chromosome position and p-value

Article of interest: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3023908/#SD1

Referring to Figure 2b of that article, if I understand correctly, the x-axis refers to the position of each variant in a chromosome while the y-axis is a log transform of the raw p-value of association of each variant with some variable of interest.

My questions is, is there practical importance in finding strong association (high log transform p-value) around one particular position versus finding it on some other position in the same chromosome?

My follow up question is, are these log transformed p-values publicly available?

Context: I am a statistics student studying multiple testing. My biology background is very limited.

In the specific case that you are talking about, the answer is very much yes that the specific position is important. Genes occur at (fairly) fixed positions along chromosomes.

The point of that figure, as indicated in the legend, is that the very low p-values are around the gene RPM1. That gene is known to have important natural genetic variation associated with the phenotype in question (the one for whichg they are looking for genetic associations). That figure is therefore an important validation of the methods of the authors.

In terms of finding such data, a fun account on twitter is the GWASbot, that just posts random manhattan plots (which is the informal name for the plot you are talking about), along with links to source data. I believe that that source data includes the p-values for each ascertained position, along with other summary statistics. One example is here.

If you are specifically interested in the topic of multiple test correction, an important related concept in GWAS is that of genome-wide significance.

Update:

To unzip and look at the file, use bgzip:

$bgzip -d ~/Downloads/continuous-20491-both_sexes.tsv.bgz$ head ~/Downloads/continuous-20491-both_sexes.tsv
# prints the following:
chr pos ref alt af_meta beta_meta   se_meta pval_meta   pval_heterogeneity  af_AFR  af_CSA  af_EUR  beta_AFR    beta_CSA    beta_EUR    se_AFR  se_CSA  se_EUR  pval_AFR    pval_CSA    pval_EUR    low_confidence_AFR  low_confidence_CSA  low_confidence_EUR
1   11063   T   G   2.088e-02   2.247e-02   2.072e-01   9.137e-01   0.000e+00   2.088e-02   NA  5.212e-05   2.247e-02   NA  1.287e-01   2.072e-01   NA  5.080e-01   9.137e-01   NA  8.000e-01   false   NA  true
1   13259   G   A   2.805e-04   9.840e-02   8.909e-02   2.694e-01   1.000e+00   NA  NA  2.805e-04   NA  NA  9.840e-02   NA  NA  8.909e-02   NA  NA  2.694e-01   NA  NA  false
1   17641   G   A   8.908e-04   -1.244e-02  4.922e-02   8.005e-01   1.000e+00   NA  NA  8.908e-04   NA  NA  -1.244e-02  NA  NA  4.922e-02   NA  NA  8.005e-01   NA  NA  false
1   30741   C   A   NA  NA  NA  NA  NA  9.400e-04   NA  NA  3.935e-03   NA  NA  1.297e+00   NA  NA  9.976e-01   NA  NA  true    NA  NA
1   51427   T   G   NA  NA  NA  NA  NA  1.127e-03   NA  NA  8.780e-01   NA  NA  9.366e-01   NA  NA  3.486e-01   NA  NA  true    NA  NA
1   57222   T   C   7.078e-04   -4.153e-03  5.656e-02   9.415e-01   1.000e+00   NA  NA  7.078e-04   NA  NA  -4.153e-03  NA  NA  5.656e-02   NA  NA  9.415e-01   NA  NA  false
1   58396   T   C   2.735e-04   3.011e-02   8.221e-02   7.142e-01   1.000e+00   2.622e-03   NA  2.735e-04   -1.069e-01  NA  3.011e-02   7.079e-01   NA  8.221e-02   8.800e-01   NA  7.142e-01   true    NA  false
1   62745   C   G   NA  NA  NA  NA  NA  2.701e-03   NA  NA  5.524e-01   NA  NA  5.529e-01   NA  NA  3.178e-01   NA  NA  true    NA  NA
1   63668   G   A   NA  NA  NA  NA  NA  NA  1.129e-03   3.673e-05   NA  2.205e-01   -4.462e-02  NA  9.738e-01   2.544e-01   NA  8.209e-01   8.608e-01   NA  true    true

• Thank you. I tried downloading an item from the GWAS bot. It was a .bgz file. May I ask what the best way is to extract data from this filetype? I tried to open it with Excel but it did not work. Aug 2 at 23:49
• @variancekills See update- you need to decompress the file, then it will be a TSV that you can in theory look at in excel, though it is pretty big for excel (several GB). I would suggest R, python or similar if you are familiar Aug 3 at 0:28
• R is fine for me. Thank you very much for your help! I've been working on a small pet project for multiple testing that requires a very specific data setup. I need thousands of hypotheses that I can line up meaningfully. This setup works very well for testing my method. Aug 3 at 13:08
• thank you very much. i was able to open the file. May I ask you one more question? Which column in the example you posted is the actual p-value used to compute the logtransformed p-values on the plot? Aug 3 at 23:38
• Thanks. Yes, I was able to reconstruct the map using pval_meta. This is very helpful. Thank you very much! Aug 4 at 1:01

Are these log transformed p-values publicly available?

For human population data, it's pretty common for the association statistics to be either included with the research paper, or uploaded to a public database (e.g. see curated compilations like GWAS Central or NHGRI GWAS Catalog).

For plant and other non-human data, there doesn't seem to have been as much drive to include GWAS data with the research paper, and it's often up to the authors to decide to do that. I find this a little odd, given that there should be fewer ethical barriers to releasing non-human public data.

If you're interested in a particular dataset, you can contact the authors, but I acknowledge that such an approach can be too high of a barrier for most interested researchers.

Is there practical importance in finding strong association (high log transform p-value) around one particular position versus finding it on some other position in the same chromosome?

It's common for statistical tests for association to assume independence of variant locations. With this in mind, two associated regions that are one base apart should not be any more interesting than two associated regions that are 50 megabases apart.

The problem with this assumption is that DNA variation is not statistically independent. Recombination of DNA (and associated reassortment of genetic variation) between two locations is more common when the locations are close, and this linkage should really be accounted for when doing association tests. One way to do this is to combine variants that share the same ancestral origin into a haplotype, then do association tests on the haplotype block rather than the single-base genotype. Even if the particular observed variants are not directly causative, there may be another untyped causative variant (or combination of variants) on the same haplotype block, meaning that the block becomes a proxy for the causative variant.

But this haplotype-based approach becomes quickly computationally intractable when looking at larger regions of chromosomes, and is very difficult to use for comparing populations with substantially different ancestral backgrounds. The workaround (if it could be called that) is to ignore the linkage aspect of genetic variation, and return to those same statistical tests, but declare regions that have peaks of close association to be statistically interesting, and zoom into them to look in more detail at the region.

That zoom in would be a great opportunity to use more complex / accurate haplotype-based tests for confirming an associative signal around a particular region. But for the most part, as with the paper you've indicated, the tests are identical.

However, even ignoring the haplotype issue... <begin rant>

The most common manhattan plots represent only p-value on the Y-axis. I think this is a mistake, because the p-value only represents the confidence in a statistic, and the important thing for association is the statistic itself, rather than the confidence in that statistic.

• there may be values that have less confidence, but a much higher associative statistic, and contribute more to the trait of interest (e.g. see here)
• values with a low association statistic are more sensitive to genotyping error, leading to spurious associations.
• the p-value doesn't indicate the direction of association
• p-values should not be used for ranking, and this ranking is implicitly done when using the p-value only manhattan plot

In addition, the generated p-values often don't make sense, and I haven't seen any commentary on this. Anything smaller than $$10^{-40}$$ doesn't represent any reasonable natural random occurrence, even on galactic timescales, but often these astronomically-small p-values can become insignificant through a few random mutations in one of the populations.

As a quick-fix using the existing statistics that are typically presented in GWAS results, where single points are desired, I recommend that people instead use the standard deviation to determine the minimum expected beta score:

see here for more details.

More rigorously, I think that people should consider doing their own internal validation of associations via bootstrap sub-sampling, and report a non-parametric statistic: the best observed rank after multiple bootstrap sub-samples. This requires access to genotype-level data, so cannot be done on existing GWAS results. I explore this in more detail in this preprint.

Yes, I'm aware this goes against the grain of how GWAS is typically done. I have been emboldened in my position on this after reading the ASA principles about the p-value, and seeing that many of the principles are in direct opposition to the typical manhattan plot:

https://www.tandfonline.com/doi/full/10.1080/00031305.2016.1154108#_i31

• I think that this is a fair, well-argued criticism of GWAS methodology, supported by references, and which I personally tend to agree with. I myself have serious problems with NHST, specifically in high-n cases like genomics. But I'm not sure that it really gets at the asker's question regarding getting a baseline understanding of what GWAS is and some of the fundamental ideas that motivate it. And yet, I do think that it's valuable to have this answer's viewpoint represented on this page. Aug 4 at 0:32
• Yes, fair comment. I've added a couple of direct answers to the specific questions asked. My additional rant was on a related issue which answers other questions which I think are more relevant for interpreting GWAS data.
– gringer
Aug 4 at 5:24