Abstract
Transcriptome-wide association studies (TWAS) integrate gene expression prediction models and genome-wide association studies (GWAS) to identify gene–trait associations. The power of TWAS is determined by the sample size of GWAS and the accuracy of the expression prediction model. Here, we present a new method, the Summary-level Unified Method for Modeling Integrated Transcriptome using Functional Annotations (SUMMIT-FA), that improves the accuracy of gene expression prediction by leveraging functional annotation resources and a large expression quantitative trait loci (eQTL) summary-level dataset. We build gene expression prediction models using SUMMIT-FA with a comprehensive functional database MACIE and the eQTL summary-level data from the eQTLGen consortium. By applying the resulting models to GWASs for 24 complex traits and exploring it through a simulation study, we show that SUMMIT-FA improves the accuracy of gene expression prediction models in whole blood, identifies significantly more gene-trait associations, and improves predictive power for identifying “silver standard” genes compared to several benchmark methods.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This study is supported by National Institutes of Health (R03 AG070669) and has been conducted using the UK Biobank recourse (application numbers 48240).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
The GWAS summary data used in this study are summarized in Supplementary Data 1 (with the download link). The eQTL summary data are available at https://www.eqtlgen.org/cis-eqtls.html The UK Biobank is an open-access resource requiring registration, available at https://www.ukbiobank.ac.uk/researchers/. The genotype and RNA sequencing data for the GTEx project may be found at the database of Genotypes and Phenotypes (accession number phs000424.v8.p2, https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000424.v8.p2). The processed gene expression data from the GTEx project is available from the GTEx portal (https://gtexportal.org). The MR-JTI, PrediXcan, and UTMOST models may be downloaded from https://doi.org/10.5281/zenodo.3842289. The TWAS-FUSION models may be down-loaded from http://gusevlab.org/projects/fusion/. The 1000 Genomes Project data may be downloaded from https://www.internationalgenome.org/data. The genetic distance data for 1000 Genomes Project may be downloaded from https://github.com/joepickrell/1000-genomes-genetic-maps.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
The GWAS summary data used in this study are summarized in Supplementary Data 1 (with the download link). The eQTL summary data are available at https://www.eqtlgen.org/cis-eqtls. html The UK Biobank is an open-access resource requiring registration, available at https://www.ukbiobank.ac.uk/researchers/. The genotype and RNA sequencing data for the GTEx project may be found at the database of Genotypes and Phenotypes (accession number phs000424.v8.p2, https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000424.v8.p2). The processed gene expression data from the GTEx project is available from the GTEx portal (https://gtexportal.org). The MR-JTI, PrediXcan, and UTMOST models may be downloaded from https://doi.org/10.5281/zenodo.3842289. The TWAS-FUSION models may be down-loaded from http://gusevlab.org/projects/fusion/. The 1000 Genomes Project data may be downloaded from https://www.internationalgenome.org/data. The genetic distance data for 1000 Genomes Project may be downloaded from https://github.com/joepickrell/1000-genomes-genetic-maps. The SUMMIT-FA models generated for this study will be available at OSF.IO.
https://www.eqtlgen.org/cis-eqtls.html
https://www.ukbiobank.ac.uk/researchers/
https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000424.v8.p2
https://doi.org/10.5281/zenodo.3842289
http://gusevlab.org/projects/fusion/