Lastly, per-gene permutations were performed to identify genes with at least one significant eQTL after correcting for multiple marker testing13. 1000 permutations were performed per gene and FDR was estimated on the permutation P values using the qvalue R package (Dabney A and Storey JD. qvalue: Q-value estimation for false discovery rate control. R package version 1.43.0).
