Dataset Open Access

GWAS and GTEx QTL integration

Barbeira, Alvaro Numa; Bonazzola, Rodrigo; Gamazon, Eric R; Liang, Yanyu; Park, YoSon; Kim-Hellmuth, Sarah; Wang, Gao; Jiang, Zhuoxun; Zhou, Dan; Hormozdiari, Farhad; Liu, Boxiang; Rao, Abhiram; Hamel, Andrew R; Pividori, Milton D; Aguet, François; Bastarache, Lisa; Jordan, Daniel M; Verbanck, Marie; Do, Ron; Stephens, Matthew; Ardlie, Kristin; McCarthy, Mark; Montgomery, Stephen B; Segré, Ayellet V; Brown, Christopher D; Lappalainen, Tuuli; Wen, Xiaoquan; Im, Hae Kyung


Citation Style Language JSON Export

{
  "publisher": "Zenodo", 
  "DOI": "10.5281/zenodo.3518299", 
  "container_title": "Widespread dose-dependent effects of RNA expression and splicing on complex diseases and traits", 
  "title": "GWAS and GTEx QTL integration", 
  "issued": {
    "date-parts": [
      [
        2019, 
        10, 
        24
      ]
    ]
  }, 
  "abstract": "<p># Data usage policy</p>\n\n<p>When using this data, you must acknowledge the source by citing the publication &quot;Widespread dose-dependent effects of RNA expression and splicing on complex diseases and traits&quot; (https://doi.org/10.1101/814350).</p>\n\n<p>&nbsp;</p>\n\n<pre><em># GTEx GWAS integration\n</em>\nThis package contains the application of several GWAS-QTL integration methods.\nThe results were analyzed in [this preprint](<em>https://www.biorxiv.org/content/10.1101/814350v1</em>)\nabout GTEx v8 application to several GWAS traits.\n \n<em>``` \n</em><em>.\n</em><em>|-- colocalization\n</em><em>|   |-- coloc\n</em><em>|   |   `-- coloc_enloc_priors_eqtl.tar.gz\n</em><em>|   |-- enloc\n</em><em>|   |   |-- enloc_eqtl_eur.tar.gz\n</em><em>|   |   `-- enloc_sqtl_eur.tar.gz\n</em><em>|   `-- eur_ld.bed.gz\n</em><em>|-- prediction_models\n</em><em>|   |-- gtex_v8_expression_mashr_snp_smultixcan_covariance.txt.gz\n</em><em>|   |-- gtex_v8_splicing_mashr_snp_smultixcan_covariance.txt.gz\n</em><em>|   |-- mashr_eqtl.tar\n</em><em>|   `-- mashr_sqtl.tar\n</em><em>|-- smr\n</em><em>|   |-- SMR_gtex_v8_README.txt\n</em><em>|   `-- SMRresults_GTEx_v8_peQTL5e-08.tar.gz\n</em><em>|-- smultixcan\n</em><em>|   |-- smultixcan_eqtl.tar.gz\n</em><em>|   `-- smultixcan_sqtl.tar.gz\n</em><em>`-- spredixcan\n</em><em>    |-- spredixcan_eqtl.tar.gz\n</em><em>    `-- spredixcan_sqtl.tar.gz\n</em>\n<em> ```\n</em><em> \n</em>You can uncompress gzipped tarball packages <em>`*.tar.gz` </em>in a UNIX command line with an instruction such as:\n<em>```bash\n</em><em>tar -xzvpf smultixcan_eqtl.tar.gz\n</em><em>```\n</em>, and the tar packages (<em>`*.tar`</em>) with an analogous instruction:\n<em>```bash\n</em><em>tar -xvpf mashr_eqtl.tar\n</em><em>```\n</em>\n\n<em>## Preliminaries\n</em>\n<strong>**</strong>Finemapping<strong>** </strong>results are contained in a separate release due to size constraints.\n\nGWAS summary statistics for 114 traits were harmonized and imputed to GTEx v8 variants with MAF&gt;0.01 using only european samples.\n(summary imputation software [here](<em>https://github.com/hakyimlab/summary-gwas-imputation</em>)). \nSome of the following analyses used the full set of 114 traits,\nwhile some focused only on 87 traits whose imputed associations showed no deflation\n(the imputation algorithm is conservative, and studies with too few available variants have a depleted distribution of association p-values after imputation).\n\nThe harmonized and imputed GWAS summary statistics are contained in a separate release due to size constraints. \nFor completeness&#39; sake, the imputed summary statistics look like:\n<em>```\n</em><em>variant_id panel_variant_id   chromosome position   effect_allele  non_effect_allele  current_build  frequency  sample_size    zscore pvalue effect_size    standard_error imputation_status  n_cases\n</em><em>rs554008981    chr1_13550_G_A_b38 chr1   13550  A  G  hg38   0.017316017316017316   336474 -2.2919929353647097    0.021906050841240293   NA NA imputed    NA\n</em><em>rs201055865    chr1_14671_G_C_b38 chr1   14671  C  G  hg38   0.012987012987012988   336474 -0.9559192804440632    0.33911301727494103    NA NA imputed    NA\n</em><em>...\n</em><em>```\n</em>\nThe GWAS were split in approximately independent LD regions (Berisa-Pickrell)/\nGWAS regions are defined in <em>`eur_ld.bed.gz` </em>(note that a few of them are ill-defined in hg38 and where ignored; only completely defined regions were used). \n\n<em>## Colocalization\n</em>\n<em>### Enloc\n</em>\nENLOC ([see fotware here](<em>https://github.com/xqwen/integrative</em>))\nwas run for sQTLs and eQTLs using individuals of european ancestry and DAP-G QTL enrichment results on 87 traits.\nResult files are included in <em>`enloc_eqtl_eur.tar.gz` </em>and <em>`enloc_sqtl_eur.tar.gz`\n</em>Each file contains a particular tissue-trait combination.\nEach row details colocalization between a GWAS region (Berisa-Pickrell) and gene&#39;s or intron&#39;s cis-window.\n\nA region might overlap multiple genes/introns or viceversa.\nEach ENLOC file contains the following columns:\n\n<strong>* </strong>gwas_locus: GWAS LD region\n<strong>* </strong>molecular_qtl_trait: gene or intron\n<strong>* </strong>locus_gwas_pip: posterior inclusion probability of variants in the GWAS LD region\n<strong>* </strong>locus_rcp: regional colocalization probability (main colocalization measure)\n<strong>* </strong>lead_coloc_SNP: snp with highest RCP\n<strong>* </strong>lead_snp_rcp: rcp of the lead coloc snp\n\n\n<em>### Coloc\n</em>\nColoc ([see software here](<em>https://cran.r-project.org/web/packages/coloc/index.html</em>))\nwas run using prior probabilities estimated from QTL enrichment of GWAS variants (computed via ENLOC).\nResults for eQTL are available in <em>`coloc_enloc_priors_eqtl.tar.gz`</em>. \nEach file contains results for a trait-tissue combination. Columns are:\n<strong>* </strong>gene_id: gene or intron id\n<strong>* </strong>p0: probability that neither QTL nor GWAS contain a causal variant\n<strong>* </strong>p1: probability that only GWAS contains a causal variant\n<strong>* </strong>p2: probability that only QTL has a causal variant\n<strong>* </strong>p3: probability that GWAS and QTL have a causal variant and it&#39;s distinct\n<strong>* </strong>p4: probability that GWAS and QTL have a causal variant and it&#39;s the same (main colocalization measure)\n\n<em>## PrediXcan\n</em>\n<em>`mashr_eqtl.tar` </em>and <em>`mashr_sqtl.tar` </em>contain prediction models \n(trained on expression or splicing data respectively, for 49 GTEx tissues) and LD compilations \nto be used with PrediXcan, S-PrediXcan, MultiXcan and S-MultiXcan.\n\nFor every tissue, the <em>`mashr_{tissue}.db` </em>file is a SQLite file with the prediction model definitions.\n<em>`mashr_{tissue}.txt.gz` </em>is a gzipped-text file with the upper triangular matrices of covariance between snps\nwithin a gene/intron prediction model.\n\nMany variants in these models don&#39;t have an rsid. To fully leverage the information in these models, \nit is advised to at least harmonize to GTEx variants, and if possible impute as we did [here](<em>https://github.com/hakyimlab/summary-gwas-imputation</em>).\n\n<em>### S-PrediXcan\n</em>\nS-PrediXcan was run for the 114 harmonized and imputed traits, on eQTL and sQTL mashr prediction models.\nAll of the GWAS traits had the same format, so that the following format parameters were used with S-PrediXcan:\n\n<em>```\n</em><em>--snp_column panel_variant_id --effect_allele_column effect_allele --non_effect_allele_column non_effect_allele --zscore_column zscore \\\n</em><em>--keep_non_rsid --additional_output --model_db_snp_key varID \\\n</em><em>```\n</em>\nEach file is a CSV, with each row containing a gene/intron association at a given trait-tissue combination:\n<strong>* </strong>gene: ENSEMBLE ID or intron id\n<strong>* </strong>gene_name: HUGO name or intron id\n<strong>* </strong>zscore: predicted association z-score\n<strong>* </strong>effect_size: estimated effect size\n<strong>* </strong>pvalue: association p-value\n<strong>* </strong>var_g: estimated variance of predicted expression or splicing\n<strong>* </strong>pred_perf_r2: prediction model cross-validated performance\n<strong>* </strong>pred_perf_pval: prediction model cross-validated performance\n<strong>* </strong>pred_perf_qval: deprecated, empty field left for compatibility\n<strong>* </strong>n_snps_used: number of snps in the intersection of GWAS and model\n<strong>* </strong>n_snps_in_cov: number of snps in the LD compilation\n<strong>* </strong>n_snps_in_model: number of snps in the model\n<strong>* </strong>best_gwas_p: smallest p-value acros GWAS snps used in this model\n<strong>* </strong>largest_weight: largest prediction model weight\n\n<em>### S-Multixcan\n</em>\nS-MultiXcan results were generated from the above S-PrediXcan results. Each fiel contains multi-tissue associations for a given trait:\n\n\n<strong>* </strong>gene: ENSEMBLE ID or intron id\n<strong>* </strong>gene_name: HUGO name or intron id\n<strong>* </strong>pvalue: multi-tissue association p-value\n<strong>* </strong>n: number of models avialble for this gene/intron\n<strong>* </strong>n_indep: number of independent components of variation in predicted expression/splicing (surviving principal components) \n<strong>* </strong>p_i_best: highest single-tissue p-value (S-PrediXcan) \n<strong>* </strong>t_i_best: tissue of highest p-value\n<strong>* </strong>p_i_worst: lowest single-tissue p-value (S-PrediXcan)\n<strong>* </strong>t_i_worst: tissue of lowest p-value\n<strong>* </strong>eigen_max: maximum eigenvalue of SVD\n<strong>* </strong>eigen_min: minimum eigenvalue of SVD\n<strong>* </strong>eigen_min_kept: smallest eigenvalue retained after discarding smallest variations\n<strong>* </strong>z_min: minimum single-tissue z-score\n<strong>* </strong>z_max: maximum single-tissue z-score\n<strong>* </strong>z_mean: mean single-tissue zscre\n<strong>* </strong>z_sd: standard deviation of the single-tissue z-scores\n<strong>* </strong>tmi: trace of M * M_i where M is predicted expression/splicing covariance across tissues for a gene, and M_i is its SVD pseudo-inverse\n<strong>* </strong>status: computation status, 0 if no errors\n\n<em>## SMR\n</em>\nSee <em>`SMR_gtex_v8_README.txt` </em>for details.</pre>\n\n<p>&nbsp;</p>\n\n<p>&nbsp;</p>\n\n<p># Disclaimer</p>\n\n<p>The data is provided &quot;as is&quot;, and the authors assume no responsibility for errors or omissions. &nbsp;<br>\nThe User assumes the entire risk associated with its use of these data. &nbsp;<br>\nThe authors shall not be held liable for any use or misuse of the data described and/or contained herein. &nbsp;<br>\nThe User bears all responsibility in determining whether these data are fit for the User&#39;s intended use. &nbsp;</p>\n\n<p>The information contained in these data is not better than the original sources from which they were derived,<br>\nand both scale and accuracy may vary across the data set. &nbsp;<br>\nThese data may not have the accuracy, resolution, completeness, timeliness, or other characteristics<br>\nappropriate for applications that potential users of the data may contemplate. &nbsp;<br>\n&nbsp;<br>\nThe user is responsible to comply with any data usage policy from the original GWAS studies;<br>\nrefer to the list of traits described [here](https://www.biorxiv.org/content/10.1101/814350v1)<br>\nto identify their respective Consortia&#39;s requirements.</p>\n\n<p><br>\nTHE DATA IS PROVIDED WITHOUT WARRANTY OF ANY KIND,<br>\nEXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,<br>\nFITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.<br>\nIN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY,<br>\nWHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,<br>\nOUT OF OR IN CONNECTION WITH THE DATA OR THE USE OR OTHER DEALINGS IN THE DATA.</p>", 
  "author": [
    {
      "family": "Barbeira, Alvaro Numa"
    }, 
    {
      "family": "Bonazzola, Rodrigo"
    }, 
    {
      "family": "Gamazon, Eric R"
    }, 
    {
      "family": "Liang, Yanyu"
    }, 
    {
      "family": "Park, YoSon"
    }, 
    {
      "family": "Kim-Hellmuth, Sarah"
    }, 
    {
      "family": "Wang, Gao"
    }, 
    {
      "family": "Jiang, Zhuoxun"
    }, 
    {
      "family": "Zhou, Dan"
    }, 
    {
      "family": "Hormozdiari, Farhad"
    }, 
    {
      "family": "Liu, Boxiang"
    }, 
    {
      "family": "Rao, Abhiram"
    }, 
    {
      "family": "Hamel, Andrew R"
    }, 
    {
      "family": "Pividori, Milton D"
    }, 
    {
      "family": "Aguet, Fran\u00e7ois"
    }, 
    {
      "family": "Bastarache, Lisa"
    }, 
    {
      "family": "Jordan, Daniel M"
    }, 
    {
      "family": "Verbanck, Marie"
    }, 
    {
      "family": "Do, Ron"
    }, 
    {
      "family": "Stephens, Matthew"
    }, 
    {
      "family": "Ardlie, Kristin"
    }, 
    {
      "family": "McCarthy, Mark"
    }, 
    {
      "family": "Montgomery, Stephen B"
    }, 
    {
      "family": "Segr\u00e9, Ayellet V"
    }, 
    {
      "family": "Brown, Christopher D"
    }, 
    {
      "family": "Lappalainen, Tuuli"
    }, 
    {
      "family": "Wen, Xiaoquan"
    }, 
    {
      "family": "Im, Hae Kyung"
    }
  ], 
  "note": "This data was analyzed in:\nhttps://www.biorxiv.org/content/early/2019/10/22/814350.full.pdf", 
  "type": "dataset", 
  "id": "3518299"
}
1,681
2,768
views
downloads
All versions This version
Views 1,6811,680
Downloads 2,7682,768
Data volume 4.7 TB4.7 TB
Unique views 1,4261,425
Unique downloads 1,1661,166

Share

Cite as