Abstract
Background Transcriptome-wide association study (TWAS) is an influential tool for identifying novel genes associated with complex diseases, where their genetic effects may be mediated through transcriptome. TWAS utilizes reference genetic and transcriptomic data to estimate genetic effect sizes on expression quantitative traits of target genes (i.e., effect sizes of a broad sense of expression quantitative trait loci, eQTL). These estimated effect sizes are then employed as variant weights in burden gene-based association test statistics, facilitating the mapping of risk genes for complex diseases with genome-wide association study (GWAS) data. However, most existing TWAS of Alzheimer’s disease (AD) dementia have primarily focused on cis-eQTL, disregarding potential trans-eQTL. To overcome this limitation, we applied the Bayesian Genome-wide TWAS (BGW-TWAS) method which incorporated both cis- and trans-eQTL of brain and blood tissues to enhance mapping risk genes for AD dementia.
Methods We first applied BGW-TWAS to the Genotype-Tissue Expression (GTEx) V8 dataset to estimate cis- and trans-eQTL effect sizes of the prefrontal cortex, cortex, and whole blood tissues. Subsequently, estimated eQTL effect sizes were integrated with the summary data of the most recent GWAS of AD dementia to obtain BGW-TWAS (i.e., gene-based association test) p-values of AD dementia per tissue type. Finally, we used the aggregated Cauchy association test to combine TWAS p-values across three tissues to obtain omnibus TWAS p-values per gene.
Results We identified 37 genes in prefrontal cortex, 55 in cortex, and 51 in whole blood that were significantly associated with AD dementia. By combining BGW-TWAS p-values across these three tissues, we obtained 93 significant risk genes including 29 genes primarily due to trans-eQTL and 50 novel genes. Utilizing protein-protein interaction network and phenotype enrichment analyses with these 93 significant risk genes, we detected 5 functional clusters comprised of both known and novel AD risk genes and 7 enriched phenotypes.
Conclusion We applied BGW-TWAS and aggregated Cauchy test methods to integrate both cis- and trans-eQTL data of brain and blood tissues with GWAS summary data to identify risk genes of AD dementia. The risk genes we identified provide novel insights into the underlying biological pathways implicated in AD dementia.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
S.G. and J.Y. are supported by grant R35GM138313 provided by the National Institute of General Medical Sciences, National Institute of Health.
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:
All genomic data analyzed during the course of this study are either de-identified or summary statistics, which are considered as human data according to the guidelines of National Institute of Health.
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, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Data Availability
The RNAseq transcriptomic and WGS genetic data of GTEx V8 are publicly available from dbGAP (phs000424.v8.p2) and the GTEx portal website (https://gtexportal.org/home/). The summary-level GWAS data of AD dementia were generated by Wightman D.P. et. al. by meta-analysis excluding 23andMe samples, which are publicly available at https://ctg.cncr.nl/software/summary_statistics.
List of Abbreviations
- AD
- Alzheimer’s disease
- GWAS
- Genome-wide association studies
- TWAS
- Transcriptome-wide association study
- BGW-TWAS
- Bayesian Genome-wide TWAS
- eQTL
- Expression quantitative trait loci
- GTEx
- Genotype-Tissue Expression
- ACAT-O
- Omnibus aggregated Cauchy association test
- GReX
- Genetically regulated gene expression
- BVSR
- Bayesian variable selection regression
- EM-MCMC
- Expectation-maximization Markov Chain Monte Carlo
- PCP
- Posterior causal probabilities
- WGS
- Whole genome sequencing
- TPM
- Transcripts Per Million
- PPI
- Protein-protein interaction
- MS4
- Membrane-spanning 4-domains
- ApoB
- Apolipoprotein b
- LDL-C
- Low-density lipoprotein cholesterol
- CRP
- C-reactive protein
