Large-scale transcriptome-wide association study identifies new prostate cancer risk regions

Mancuso, Nicholas; Gayther, Simon; Gusev, Alexander; Zheng, Wei; Penney, Kathryn L.; Kote-Jarai, Zsofia; Eeles, Rosalind; Freedman, Matthew; Haiman, Christopher; Pasaniuc, Bogdan

doi:10.1038/s41467-018-06302-1

Download PDF

Article
Open access
Published: 04 October 2018

Large-scale transcriptome-wide association study identifies new prostate cancer risk regions

Nature Communications volume 9, Article number: 4079 (2018) Cite this article

13k Accesses
87 Citations
13 Altmetric
Metrics details

Subjects

An Author Correction to this article was published on 08 January 2019

This article has been updated

Abstract

Although genome-wide association studies (GWAS) for prostate cancer (PrCa) have identified more than 100 risk regions, most of the risk genes at these regions remain largely unknown. Here we integrate the largest PrCa GWAS (N = 142,392) with gene expression measured in 45 tissues (N = 4458), including normal and tumor prostate, to perform a multi-tissue transcriptome-wide association study (TWAS) for PrCa. We identify 217 genes at 84 independent 1 Mb regions associated with PrCa risk, 9 of which are regions with no genome-wide significant SNP within 2 Mb. 23 genes are significant in TWAS only for alternative splicing models in prostate tumor thus supporting the hypothesis of splicing driving risk for continued oncogenesis. Finally, we use a Bayesian probabilistic approach to estimate credible sets of genes containing the causal gene at a pre-defined level; this reduced the list of 217 associations to 109 genes in the 90% credible set. Overall, our findings highlight the power of integrating expression with PrCa GWAS to identify novel risk loci and prioritize putative causal genes at known risk loci.

Association of imputed prostate cancer transcriptome with disease risk reveals novel mechanisms

Article Open access 15 July 2019

Nima C. Emami, Linda Kachuri, … John S. Witte

Integrating transcription factor occupancy with transcriptome-wide association analysis identifies susceptibility genes in human cancers

Article Open access 19 November 2022

Jingni He, Wanqing Wen, … Xingyi Guo

An integrative multi-omics analysis to identify candidate DNA methylation biomarkers related to prostate cancer risk

Article Open access 06 August 2020

Lang Wu, Yaohua Yang, … PEGASUS Consortium

Introduction

Prostate cancer (PrCa) affects ~1 in 7 men during their lifetime and is one of the most common cancers worldwide, with up to 58% of risk due to genetic factors^1,2. Genome-wide association studies (GWAS) have identified over 100 genomic regions harboring risk variants for PrCa which explain roughly one-third of familial risk^3,4,5,6,7. With few exceptions⁸, the causal variants and target susceptibility genes at most GWAS risk loci have yet to be identified. Multiple studies have shown that PrCa- and other disease-associated variants are enriched near variants that correlate with gene expression levels^{9,10,11,12,13}. In fact, recent approaches have integrated expression quantitative trait loci (eQTLs) with GWAS to implicate several plausible genes for PrCa risk (e.g., IRX4, MSMB, NCOA4, NUDT11, and SLC22A3)^{5,14,15,16,17,18,19,20,21}. While overlapping eQTLs and GWAS is powerful, the high prevalence of eQTLs²² coupled with linkage disequilibrium (LD) renders it difficult to distinguish the true susceptibility gene from spurious co-localization at the same locus²³. Therefore, disentangling LD is critical for prioritization and causal gene identification at risk loci.

Gene expression imputation followed by a transcriptome-wide association study^24,25,26 (TWAS) has been recently proposed as a powerful approach to prioritize candidate risk genes underlying complex traits. By taking LD into account across SNPs, the resulting association statistics reflect the underlying effect of steady-state gene or alternative splicing expression levels on disease risk^25,27, which can be used to identify new regions or to rank genes for functional validation at known risk regions^{24,25,26,27,28}. Here we perform a multi-tissue transcriptome-wide association study^24,25,26 to identify new risk regions and to prioritize genes at known risk regions for PrCa. Specifically, we integrate gene expression data from 48 panels measured in 45 tissues across 4448 individuals with GWAS of prostate cancer from the OncoArray in 142,392 men²⁹. Notably, we include alternatively spliced and total gene expression data measured in tumor prostate to identify genes contributing to prostate cancer risk or to continued oncogenesis. We identify 217 gene-trait associations for PrCa with 23 (11) genes identified uniquely using models of alternative spliced (total) expression in tumor. Significant genes were found in 84 independent 1 Mb regions, of which 9 regions are located more than 2 Mb away from any OncoArray GWAS significant variants, thus identifying new candidate risk regions. Second, we use TWAS to investigate genes previously reported as susceptibility genes for prostate cancer identified by eQTL-based analyses. We find a significant overlap with 56 out of 102 previously reported genes assayed in our study also significant in TWAS. Third, we use a novel Bayesian prioritization approach to compute credible sets of genes and prioritize 109 genes that explain at least 90% of the posterior density for association signal at TWAS risk regions. One notable example, IRX4, had 97% posterior probability to explain the association signal at its region with the remaining 3% explained by 9 neighboring genes. Overall, our findings highlight the power of integrating gene expression data with GWAS and provide testable hypotheses for future functional validation of prostate cancer risk.

Results

Overview of methods

To identify genes associated with PrCa risk, we performed a TWAS using 48 gene expression panels measured in 45 tissues^{22,30,31,32,33,34,35,36} integrated with summary data from the OncoArray PrCa GWAS of 142,392 individuals of European ancestry (81,318/61,074 cases/controls; Methods)²⁹. We performed the summary-based TWAS approach as described in ref. ²⁵ using the FUSION software (Methods). Briefly, this approach uses reference linkage disequilibrium (LD) and reference gene expression panels with GWAS summary statistics to estimate the association between the cis-genetic component of gene expression, or alternative splicing events, and PrCa risk²⁵. First, for each panel, FUSION estimated the heritability of steady-state gene and alternative splicing expression levels explained by SNPs local to each gene (i.e., 1 Mb flanking window) using the mixed-linear model (see Methods). Genes with nominally significant (P < 0.05) estimates of SNP-heritability (cis-$h_g^2$), are then put forward for training predictive models. Genes with non-significant estimates of heritability are pruned, as they are unlikely to be accurately predicted. Next, FUSION fits predictive linear models (e.g., Elastic Net, LASSO, GBLUP³⁷, and BSLMM³⁸) for every gene using local SNPs. The model with the best cross-validation prediction accuracy (significant out-of-sample R²; nominal P < 0.05) was used for prediction into the GWAS cohort. This was repeated for all expression datasets, resulting in 109,170 tissue-specific models spanning 15,383 unique genes using total expression and 4990 using alternatively spliced introns for a combined 16,389 unique genes. The average number of models per expression panel was 2228 (Supplementary Data 1). Gene expression measured in normal prostate tissue from GTEx²² resulted in only 710 gene models, which can be explained due to smaller sample size (N = 87) compared with the average (N = 234; Supplementary Data 1). Indeed, the number of gene models per panel was highly correlated with sample size, which implies that statistical power to detect genes with cis-regulatory control is limited by sample size (Supplementary Figure 1). Focusing only on models capturing total gene expression, genes on average had heritable levels of expression in 6.1 different panels (median 3) with 10,628/15,383 genes having heritable expression in at least two panels (Fig. 1). We found R² for predictive models was largely consistent across genomic locations, and predominantly affected by the number of non-zero weights used for prediction (Supplementary Figure 2). Predictive power of linear gene expression models is upper-bounded by heritability; thus, we use a normalized R² to measure in-sample prediction accuracy (R²/ cis-$h_g^2$). We found the average R²/ cis-$h_g^2$ across all tissue-specific models was 65%, which indicates that most of the signal in cis-regulated total expression and alternative splicing levels is captured by the fitted models (Fig. 1). To assess the predictive stability for models of normal prostate gene expression, we compared measured and predicted gene expression for TCGA^36,39 normal prostate samples using models fitted in GTEx²² normal prostate. We found a highly significant replication (mean R² = 0.07; P = 1.5 × 10⁻²⁹), explaining 41% of in-sample cross-validation R² (Supplementary Figure 3), which is consistent with previous out-of-sample estimates^24,25. We performed a cross-tissue analysis within TCGA and found tumor prostate gene expression models replicated in normal prostate (total expression R² = 0.06; splicing R² = 0.05; Supplementary Table 1). Given the large number of genes having evidence of genetic control across multiple tissues, we next aimed to measure the similarity of different tissue models (Methods). Across all reference panels for each gene we observed an average R² = 0.64 (Supplementary Figure 4). Similarly, when averaging across genes, reference panels displayed an average cross-tissue R² = 0.52 (Supplementary Figure 5). Together, these results suggest that trained models predict similar levels of cis-regulated expression on average, despite reference panels measuring expression in different tissues, with varying QC, and differing capture technologies. Next, we performed simulations to measure the statistical power of TWAS under a variety of trait architectures (Supplementary Note 1). Consistent with previous work, we found TWAS to be well-powered at various effect-sizes and heritability levels for gene expression. Importantly, we found no inflation under the null when cis-regulated gene expression has no effect on downstream trait (Supplementary Figure 6).

TWAS identifies 217 genes associated with PrCa status

In total, we tested 109,170 tissue-specific gene models of expression for association with PrCa status and observed 892 reaching transcriptome-wide significance (P_TWAS < 4.58 × 10⁻⁷; two-tailed Z-test), resulting in 217 unique genes, of which 114 were significant in more than one panel (Supplementary Data 2; Fig. 2). On average, we found 18.2 tissue-specific models associated with PrCa per reference expression panel (Supplementary Data 1). In 1 Mb regions with at least 1 transcriptome-wide significant gene, we observed 10.6 tissue-specific associated models on average, and 2.6 associated genes on average, indicating that further refinement of association signal at TWAS risk loci is necessary. To quantify the overlap between non-HLA, autosomal risk loci in the OncoArray PrCa GWAS and our TWAS results, we partitioned GWAS summary data into 1 Mb regions and observed 131 harboring at least one genome-wide significant SNP. Of these, 127/131 overlapped at least one gene model in our data and 68/131 overlapped at least one transcriptome-wide significant gene (Supplementary Figure 7). Associated genes were the closest gene to the top GWAS SNP 20% of the time when using 26,292 RefSeq genes. This result is consistent with previous reports^9,25,26 and suggests that prioritizing genes based on distance to index SNPs is suboptimal. We found gene model associations were largely consistent, further supporting the predictive stability of models using cis-SNPs (Supplementary Figure 8; Supplementary Note 1). We observed little evidence of prediction accuracy introducing biased results (Supplementary Figure 9; Supplementary Note 1). As a partial control, we compared TWAS results with S-PrediXcan, a related method for predicting gene expression into GWAS summary statistics, using independently trained models and observed a strong correlation (R = 0.90; see Supplementary Figure 10; Supplementary Note 1), further supporting the validity of the TWAS approach.

Most of the gene models captured total expression levels in normal tissues, however as a positive control we included models for total expression in tumor prostate tissue (Methods). Predicted expression using tumor prostate models accounted only for 43/217 significant genes compared with 6/217 in normal prostate which is likely due to the large difference in sample size between the original reference panels (Supplementary Data 1). Given this, we found no significant increase in proportion of tumor prostate associated models compared with normal prostate (Fisher’s exact P = 0.22). Of the 309 genes with models trained in both reference panels a single shared gene, MLPH (OMIM: 606526, a gene whose function is related to melanosome transport⁴⁰), was associated with PrCa risk. In all, 11/43 genes were significant only in tumor prostate models of total expression. We found, 7/11 genes were modeled in other panels but did not reach transcriptome-wide significance while the other 4/11 were not significantly heritable, and thus not testable, in other panels. We also tested models of alternatively spliced introns for association to PrCa risk. We identified predicted expression of alternatively spliced introns in tumor prostate accounted for 68/217 genes, with an average of 2.5 (median 1) alternatively spliced intron associations per significant gene. We next quantified the amount of overlap between results driven from models of alternative splicing events versus models of total gene expression. In all, 23/68 genes were found only in alternatively spliced introns, and 14/23 genes had models of total gene expression but did not reach transcriptome-wide significance. The remaining 9/23 were tested solely in alternatively spliced introns, due to heritability of total gene expression not reaching significance. Together these results emphasize earlier work demonstrating that sQTLs for a gene commonly capture signal independent of eQTLs⁴¹.

TWAS analysis increases power to find PrCa associations

Most of the power in the TWAS approach can be attributed to large GWAS sample size. However, two other factors can increase power over GWAS. First, TWAS carries a reduced testing burden compared with that of GWAS, due to TWAS having many fewer genes compared with SNPs. In all, 9/217 genes were located at nine novel independent 1 Mb regions (i.e., no overlapping GWAS SNP), all of which remained significant under a summary-based permutation test (P < 0.05/9; Table 1; Supplementary Data 2; Methods). We found this result was stable to increasing region sizes (Supplementary Data 3) and unlikely to be the result of long-range tagging with known GWAS risk (Supplementary Data 4; Supplementary Note 1). We observed increased association signal for SNPs at these regions compared to the genome-wide background after accounting for similar MAF and LD patterns (Supplementary Figure 11), which, together with observed TWAS associations, suggests that GWAS sample size is still a limiting factor in identifying PrCa risk SNPs. As a partially independent check, we performed a multi-tissue TWAS using summary data from an earlier PrCa GWAS (N = 49,346)⁷ and found 2 novel regions. We found both regions to overlap a genome-wide significant SNP within 1 Mb in this data further supporting the robustness of TWAS (Supplementary Table 2). Second, we expect to observe increased association signal when expression of a risk gene is regulated by multiple local SNPs²⁵. We observed 88/892 instances across 28 genes where TWAS association statistics were stronger than the respective top overlapping GWAS SNP statistics (one-sided Fisher’s exact P < 2.2 × 10⁻¹⁶; 6.5% higher χ² statistics on average). For example, GRHL3 (OMIM:608317; a gene associated with suppression of squamous cell carcinoma tumors⁴²) exhibited stronger signal in TWAS using expression in prostate tumor (P_TWAS = 9.38 × 10⁻¹⁰) compared with the lead SNP signal (P_GWAS = 1.49 × 10⁻⁵). Similarly, POLI (OMIM:605252, a DNA repair gene associated with mutagenesis of cancer cells^43,44) resulted in larger TWAS associations (P_TWAS = 2.29 × 10⁻⁸) compared with the best proximal SNP (P_GWAS = 5.44 × 10⁻⁷).

Table 1 Novel risk loci

Full size table

TWAS replicates previously reported genes

We next sought to quantify the extent of overlapping results between TWAS and previous studies that integrated eQTL data measured in normal and tumor prostate tissues at PrCa risk regions (Methods; Supplementary Table 3)^{5,14,15,16,17,18,19,20}. We considered only autosomal, non-HLA genes which resulted in 130 previously reported genes. We found a significant overlap between reported genes, with 102/130 assayed in our study and 56/102 reaching transcriptome-wide significance in at least one of our panels (Fisher’s exact P < 2.2 × 10⁻¹⁶; Supplementary Table 3, Supplementary Data 5). For example, MLPH was reported in 4/8 studies. We found significant associations suggesting that decreased expression of MLPH in normal and tumor prostate tissue increases risk for PrCa (e.g., GTEx prostate MLPH Z_TWAS = −5.80; P_TWAS = 6.69 × 10⁻⁹; TCGA prostate Z_TWAS = −6.77; P_TWAS = 1.25 × 10⁻¹¹). Predicted MLPH in tumor prostate remained significant under permutation, which suggests that chance co-localization with GWAS risk is unlikely (Supplementary Data 2). To assess the amount of residual association signal due to genetic variation in the GWAS risk region after accounting for predicted expression of MLPH, we performed a summary-based conditional analysis (Methods). We found MLPH to explain most of the signal at its region (lead SNP P_GWAS = 4.03 × 10⁻¹¹; conditioned on MLPH lead SNP P_GWAS = 1.13 × 10⁻³; Fig. 3). Our findings are consistent with recent work that found decreased expression levels of MLPH to be associated with increased PrCa risk⁴⁵. Despite previous eQTL data focusing on normal and tumor prostate tissue, we observed associations in 45 expression panels overlapping the 56 observed genes in total, underscoring earlier works demonstrating the consistency of cross-tissue cis-regulatory effects⁴⁶.

Prioritization pinpoints a single gene for most risk regions

TWAS genes are indicative of association and do not necessarily reflect causality (e.g., due to co-regulation at the same region). To prioritize genes at regions with multiple TWAS signals (Fig. 2), we used a Bayesian formulation to estimate 90%-credible gene sets (Methods). We found 109 unique genes across 84 non-overlapping 1 Mb regions comprising our 90% credible sets (Supplementary Data 6, 7). In all, 68/84 credible sets contained either a single gene or the same gene in multiple tissues. The average number of unique genes per credible set was 1.29 (median 1). We observed that 28/109 prioritized genes were previously reported in eQTL analyses^{5,14,15,16,17,18,19,20}, which supports the hypothesis that TWAS followed by Bayesian prioritization refines associations to relevant disease genes. For example, MLPH was the sole gene defining its region’s 90% credible set with a posterior probability of 94%. Similarly, SLC22A3 (OMIM: 604842; a gene involved in poly-specific organic cation transporters⁴⁷ and previously implicated in PrCa risk¹⁸) exhibited >94% posterior probability to be causal.

Prostate tissue genes have largest average effect

Given the large number of significant associations observed for non-prostate tissues in our data, we wanted to quantify which tissue is most relevant for PrCa risk. We first grouped TWAS PrCa associations into prostate/non-prostate and tested for enrichment in normal and tumor prostate expression models. Predicted expression and splicing events in normal and tumor prostate made up 221/892 associations with PrCa (Supplementary Data 2) which was highly significant compared to the grouping of all other tissues (Fisher’s exact P = 7.3 × 10⁻⁹). This measure only quantifies the total amount of observed associations and neglects average association strength. Next, we computed the mean TWAS association statistic using all genes predicted from each expression reference panel (Fig. 4). We observed the largest average TWAS associations in genes predicted from normal and tumor prostate tissue, which reaffirms our intuition of expression and splice events in prostate being the most relevant for PrCa risk. We re-ranked mean associations using only genes found to be transcriptome-wide significant and observed a similar ordering with total expression in normal prostate ranked highest (average χ² = 176.2; Supplementary Figure 12).

Discussion

Prostate cancer is a common male cancer that is expected to affect more than 180,000 men in the United States in 2017 alone⁴⁸. While GWAS has been successful in localizing risk for PrCa due to genetic variation, the underlying susceptibility genes remain elusive. Here we have presented results of a transcriptome-wide association study using the OncoArray PrCa GWAS summary statistics for over 142,000 case/control samples. This approach utilizes imputed expression levels and splicing events in the GWAS samples to identify and prioritize putative susceptibility genes. We identified 217 genes whose expression is associated with PrCa risk. These genes localized at 84 genomic regions, of which nine regions do not overlap with a genome-wide significant SNP in the OncoArray GWAS. We found 23 genes using predictive models for alternatively spliced introns in tumor prostate, which supports the its role in continued risk for tumor oncogenesis. A large fraction of identified genes was confirmed in earlier work, with 56 genes previously reported in eQTL/PrCa GWAS overlap studies. We used a novel Bayesian prioritization approach to refine our associations to credible sets of 109 genes with statistical evidence of causality under standard assumptions. Our results provide a functional map for PrCa risk which can be explored for follow-up and validation.

In this study, we compared our reported TWAS results with genes identified in previous works focusing on expression measured in normal and tumor prostate tissue. Several of these studies considered an eQTL and GWAS risk SNP to overlap if they are in linkage at a specified threshold. While these approaches are sound, they may be limited in statistical power for several reasons. First, if multiple local SNPs independently contribute to risk, overlap studies relying only on the top risk SNP will lose power. Second, earlier overlap studies used thresholds for association signal (i.e., GWAS P < 5 × 10⁻⁸) and linkage strength (i.e., LD > 0.5) to consider pairs of SNPs for evidence of expression influencing risk of PrCa. TWAS is largely agnostic to both issues as it jointly considers all SNPs in the region, regardless of reported GWAS association strength. However, when expression of a risk gene is regulated by a single causal SNP, we expect TWAS and earlier overlap approaches to have similar levels in power²⁵.

Previous works have strongly implicated expression of certain genes in PrCa risk that were not assayed in our study (e.g., MSMB^18,49) due to non-significant heritability estimates. TWAS operates by fitting predictive linear models of gene expression based on local genotype data, followed by prediction into large cohorts and subsequent association testing. Expression of genes that are not significantly heritable at current sample sizes are not included in the pipeline. This is the consequence of heritability providing an upper bound on the predictive accuracy under a linear model for genotype; therefore, if a gene has undetectable heritability at a given sample size, it will be difficult to predict using linear combinations of SNPs. To compute TWAS weights for normal prostate tissue, we used samples collected in the GTEx v6 panel (n = 87). Thus, our inability to detect heritable levels of gene expression can be explained due to the relatively small number of samples compared with other tissues. Indeed, previous work has shown a strong correlation between sample size in expression panels and the number of identified eGenes²⁷; therefore, as sample size increases for relevant tissues, we expect the number of genes included in the TWAS framework to increase. TWAS will lose power in situations where gene expression is a nonlinear function of local SNPs, or when trans (or distal) regulation is a major component in modulating expression levels.

We conclude with several caveats and possible future directions. First, while TWAS associations are consistent with models of steady-state gene expression levels altering risk for PrCa, they may be the result of confounding^25,26. Imputed gene expression levels are the result of weighted linear combinations of SNPs, many of which may tag non-regulatory mechanisms driving risk and result in inflated association statistics. Second, our results relied on validating prediction models using multiple approaches: within-reference methods (i.e., cross-validation), cross-reference methods (e.g., GTEx into TCGA), and external-reference methods (i.e., 1000 Genomes predictive stability). While results from these approaches support our models generalizing out-of-sample, we still lack within-GWAS replication of predictive models. Third, since genes with eQTLs are common, associations may be the result of chance co-localization between eQTLs and PrCa risk. Finally, we note recent work has extended TWAS-like methods to expose regulatory mechanisms for susceptibility genes by incorporating chromatin information⁵⁰. An extension to our work would be to pinpoint chromatin variation regulating expression levels at identified risk genes, thus describing a richer landscape of the molecular cascade where SNP → chromatin → expression → PrCa risk.

Methods

OncoArray GWAS summary statistics

Genome-wide association summary statistics for the OncoArray PrCa study were obtained from ref. ²⁹. Summary statistics were computed using a fixed-effect meta-analysis for 142,392 total samples of European ancestry from the OncoArray (81,318/61,074 cases/controls), UK stage 1 (1854/1894) and UK stage 2 (3706/3884), CaPS 1 (474/482) and CaPS 2 (1458/512), BPC3 (2068/3011), NCI PEGASUS (4600/2941) and iCOGS (20,219/20,440). The initial summary data contained association statistics for 19,726,430 variants. We filtered out summary statistics for SNPs with MAF <0.01 and any SNPs with ambiguous alternative alleles (e.g., A → T; C → G; or vice-versa). Finally, we kept only SNPs with rsIDs defined by dbSNP144. Our QC pipeline resulted in association statistics at 10,516,237 SNPs for downstream TWAS analyses.

Previous prostate-tissue eQTL studies

We collected previous studies that investigated the overlap of eQTLs in normal and tumor prostate tissue at known PrCa risk loci^{5,14,15,16,17,18,19,20}. We compared TWAS statistics versus reported eQTL overlap results as aggregated in refs. ^14,15. Across these studies, overlap of eQTLs and PrCa risk loci are computed by one of two possible methods. The first method tests known PrCa risk SNPs for association with expression levels of nearby genes/transcripts. The second method takes a two-step approach. First, genes nearby PrCa risk loci are tested for harboring eQTLs at some significance level. Next, genes with identified eQTL SNPs are tested to be in LD with known PrCa risk variants at some level (e.g., r² > 0.5).

Reference gene expression data and predictive models of expression

We downloaded the FUSION software (see URLs) along with its prepackaged weights for gene expression data. FUSION is an R package that implements the TWAS scheme described in ref. ²⁵. Weights for gene expression measured using RNA sequencing data were obtained from the CommonMind Consortium³⁰ (dorsolateral prefrontal cortex, n = 452), the Genotype-Tissue Expression Project²² (GTEx; 44 tissues; n = 449), the Metabolic Syndrome in Men study^32,33 (adipose, n = 563), and The Cancer Genome Atlas (TCGA; prostate adenocarcinoma, n = 483)³⁹. Expression microarray data were obtained from the Netherlands Twins Registry³⁵ (NTR; blood, n = 1247), and the Young Finns Study^31,34 (YFS; blood, n = 1264). All non-TCGA expression panel individuals were PrCa controls. Detailed description of quality control procedures on measured gene expression and genotype information for all non-TCGA reference panels are described in refs. ^25,27. TCGA genotype, gene expression, and exon-junction data for 525 samples were downloaded using the Broad GDAC FireHose version 2016_1_28 (see URLs). Genotypes were imputed to the Haplotype Reference Consortium⁵¹ and restricted to well-imputed (INFO > 0.9) HapMap3⁵² sites. Genes (exon junctions) missing in more than half of samples were removed. RPKM and log-adjusted gene expression levels were estimated in a generalized linear model controlling for three gene expression PCs. The estimated log-abundances were quantile-normalized and inverse-normal rank-normalized. We estimated alternatively spliced introns using the software MapSplice version 2 (see URLs). A total of 482 samples passed quality control procedures in both genotype and gene expression data. We note that batch effects from measurement biases (e.g., RNA-degradation) should be uncorrelated with SNPs local to a gene body definition, and therefore not impact prediction accuracy. By maximizing the sample size, predictive power when using cis-SNPs should increase and be largely unbiased. This is evidenced by the fact that models are largely stable across and within TCGA PRAD datasets (Supplementary Table 1).

We filtered genes that did not exhibit cis-genetic regulation at current samples sizes by keeping only genes with nominally significant (P < 0.05) estimates of cis-SNP heritability (cis-$h_g^2$), which resulted in 117,459 total tissue-gene pairs from 17,023 unique genes. We refrain from reporting genes from the HLA region due to complicated LD patterns.

To train predictive models, FUSION defines gene expression for n samples (y_GE) as a linear function of p SNPs (X) in a 1 Mb region flaking the gene as

$${\mathbf{y}}_{{\mathrm{GE}}} = {\mathbf{C\beta}} + {\mathbf{Xw}}_{{\mathrm{GE}}} + {\bf{\epsilon }},$$

where w_GE are the p SNP weights, Cβ are covariates (e.g., sex, age, genotype principal components, genotyping platform, and PEER factors) and their effects, and $\epsilon$ is random environmental noise. FUSION estimated weights for expression of a gene in a tissue using multiple penalized linear models. Generally, FUSION optimizes for

$$\left[ {\begin{array}{*{20}{c}} {{\hat{\mathbf w}}_{{\mathrm{GE}}}} \\ {{\hat{\mathbf \beta }}} \end{array}} \right] = \arg \mathop {{\min }}\limits_{{\mathbf{w}}_{{\mathrm{GE}}},{\mathbf{\eta }}} \|{\mathbf{y}}_{{\mathrm{GE}}} - {\mathbf{Xw}}_{{\mathrm{GE}}} - {\mathbf{C\beta }}\|_2^2 + f\left( {{\mathbf{w}}_{{\mathrm{GE}}}} \right),$$

where f(w_GE) is a parameterized penalty function specific to each model (e.g., GBLUP³⁷, LASSO, the Elastic Net). The exception to this optimization criterion is the Bayesian sparse linear mixed model (i.e., BSLMM)³⁸ which fits the posterior mean for w_GE using MCMC in the GEMMA v 0.94 software (see URLs) to obtain weights. To determine which model has the best prediction accuracy for a given gene-tissue pair, FUSION computes out-of-sample R² by performing fivefold cross-validation for each model. We compute the normalized prediction accuracy for a gene as ${\mathrm{min}}\left( {\frac{{R^2}}{{h_g^2}},1} \right)$. Weights from the model with the largest R² that was also nominally non-zero (P < 0.05) were used to compute TWAS association statistics. This resulted in a final tally of 109,170 tissue-specific models at 16,389 unique genes.

Cis-heritability of gene expression

FUSION reports the estimated SNP-heritability (i.e., $h_g^2$) for measured gene expression levels explained by SNPs in the cis-region (1 Mb region surrounding the TSS). This is modeled under a mixed-linear model as

$${\mathrm{var}}\left( {{\mathbf{y}}\prime _{{\mathrm{GE}}}} \right) = {\mathbf{A}}{{\sigma }}_{{g}}^2 + {\mathbf{I}}{{\sigma }}_{{e}}^2,$$

where ${\mathbf{y}}\prime _{{\rm GE}}$ is the residual gene expression after regressing out fixed-effect covariates C, A is the estimated kinship matrix from SNPs in the cis-region and ${\mathrm{\sigma }}_g^2$ (${\mathrm{\sigma }}_e^2$) is the variance explained by the cis-SNPs (environment). SNP-heritability is then defined to be ratio of genotypic variance and total trait variance as, $h_g^2 = \frac{{{\mathrm{\sigma }}_g^2}}{{{\mathrm{\sigma }}_g^2 + {\mathrm{\sigma }}_e^2}}$. Variance parameters are estimated using the AI-REML algorithm implemented in GCTA v1.26 (see URLs) with the top 3 genotypic principal components, sex, age, genotyping platform, and PEER factors as covariates.

Measuring cross-tissue similarity in predicted expression

We took an unbiased approach to identify susceptibility genes for PrCa by using gene expression panels measured in various tissues. To quantify how similar predicted expression levels are for the same gene across different tissues, we measured the squared Pearson correlation (R²). This value represents how well predicted expression from one tissue predicts expression in another tissue. To dissect similarities and differences of tissue-specific models, the ideal scenario would be to inspect effects at individual SNPs defining the models. In practice this is not possible due to predictive models not including the same set of SNPs due to QC and technological differences in the original studies. Therefore, as a proxy we predict gene expression into the 489 samples of European ancestry from 1000 Genomes⁵³ and compute R² across shared genes for pairs of tissues (Supplementary Note 1).

Transcriptome-wide association study using GWAS summary statistics

FUSION estimates the strength of association between predicted expression of a gene and PrCa (z_TWAS) as function of the vector of GWAS summary Z-scores at a given cis locus z_GWAS (i.e., vector of SNP association Wald statistics) and the LD-adjusted weights vector learned from the gene expression data w_GE as

$${\rm z}_{{\rm TWAS}} = \frac{{{\mathbf{{w}}}_{{\rm GE}}^\prime {\mathbf{{z}}}_{{\rm GWAS}}}}{{\sqrt {{\mathop{\rm{var}}} ({\mathbf{{w}}}_{{\rm GE}}^\prime {\mathbf{{z}}}_{{\rm GWAS}})} }} = \frac{{{\mathbf{{w}}}_{{\rm GE}}^\prime {\mathbf{{z}}}_{{\rm GWAS}}}}{{\sqrt {{\mathbf{{w}}}_{{\rm GE}}^\prime {\mathbf{{Vw}}}_{{\rm GE}}} }},$$

where V is a correlation matrix across SNPs at the locus (i.e., LD) and “‘” indicates transpose. A P-value for z_TWAS is obtained using a two-tailed test under N(0,1). In this work, we estimated V using 489 samples of European ancestry in 1000 Genomes⁵³. To account for the large number of hypotheses tested, we perform Bonferroni correction at α = 0.05/M, where M = 109,170 is the number of predictive models, which is conservative as many gene models are correlated. As reported by ref. ²⁵, there may be inflation at GWAS risk loci, due to chance co-varying of SNP effects between expression and PrCa. The same work described a permutation procedure that assesses likelihood of observing association by chance conditioned on GWAS signal. The algorithm works by permuting the eQTL weights w_GE while keeping z_GWAS fixed and computing z_TWAS,_perm. FUSION implements an adaptive procedure that stops once enough scores (i.e. |z_TWAS,_perm|≥|z_TWAS|) have been observed such that the empirical null cannot be rejected at a specified level. We define novel risk regions as a flanking region around a transcriptome-wide significant gene (splicing event; P_TWAS < 4.58 × 10⁻⁷; two-tailed Z-test) that does not harbor a genome-wide significant SNP (P_GWAS < 5 × 10⁻⁸; two-tailed Z-test). We consider 2 Mb windows by default (i.e. TSS ± 1 Mb) and show that the results are robust to the choice of window size (Supplementary Data 3).

GWAS analyses conditional on predicted expression

To assess the extent of residual association of SNP with PrCa risk after accounting for predicted gene expression levels, FUSION estimates conditional SNP association scores using GWAS summary statistics. Namely, define V as LD for SNPs in the region, V_GE as the correlation between predicted expression levels, and C as the correlation between SNPs and predicted expression. The least-squares estimates of z_GWAS|z_TWAS are determined by,

$${\mathbf{z}}_{{\mathrm{GWAS}}}|{\mathbf{z}}_{{\mathrm{TWAS}}} = {\mathbf{z}}_{{\mathrm{GWAS}}} - {\mathbf{CV}}_{{\mathbf{GE}}}^{ - 1}{\mathbf{z}}_{{\mathrm{TWAS}}}.$$

The variance of the residual association strength is given by,

$${\mathrm{var}}\left[ {{\mathbf{z}}_{{\mathrm{GWAS}}}|{\mathbf{z}}_{{\mathrm{TWAS}}}} \right] = {\mathrm{var}}\left[ {{\mathbf{z}}_{{\mathrm{GWAS}}}} \right] - {\mathrm{var}}\left[ {{\mathbf{CV}}_{{\mathbf{GE}}}^{ - 1}{\mathbf{z}}_{{\mathrm{TWAS}}}} \right] = {\mathbf{V}} - {\mathbf{CV}}_{{\mathbf{GE}}}^{ - 1}{\mathbf{C}}^{\prime}.$$

This results in the final conditional association score for the ith SNP as,

$${{z}}_{{i}} = \left[ {{\mathbf{z}}_{{\mathrm{GWAS}}} - {\mathbf{CV}}_{{\mathbf{GE}}}^{ - 1}{\mathbf{z}}_{{\mathrm{TWAS}}}} \right]_{{i}}/\surd {\mathrm{diag}}\left[ {{\mathbf{V}} - {\mathbf{CV}}_{{\mathbf{GE}}}^{ - 1}{\mathbf{C}}^\prime } \right]_{{{ii}}}.$$

Bayes factors and posterior inference of causal genes

Complex correlations between predicted expression levels at a given region can yield multiple associated genes in TWAS (Fig. 2). Thus, for the vast majority of risk regions it remains unclear which gene is causally influencing PrCa risk. Here, modeling under the assumption of a single causal gene per risk region and relying on the central limit theorem for normality, we can compute the Bayes Factor that the ith gene in a region is causal as,

$${\rm BF}_i = \frac{{N\left( {{\rm {z}}_{{\mathrm{TWAS}},{{i}}}|0,1 + n\sigma _\alpha ^2} \right)}}{{N\left( {{\rm {z}}_{{\mathrm{TWAS}},{{i}}}|0,1} \right)}} = \left( {1 + n\sigma _\alpha ^2} \right)^{ - 1/2}\exp \left( {\frac{{{\rm {z}}_{{\mathrm{TWAS}},{{i}}}^2}}{2}\frac{{n\sigma _\alpha ^2}}{{1 + n\sigma _\alpha ^2}}} \right),$$

where ${\rm z}_{{\mathrm{TWAS}},{{i}}}^2$ is the squared TWAS association statistic for the ith gene, n is the GWAS sample size, and $\sigma _\alpha ^2$ is prior effect-size variance for gene expression on PrCa risk (Supplementary Note 1). This model is structurally similar in form to earlier works^54,55,56 describing Bayes Factors for fine mapping SNPs at GWAS risk regions. The important distinction is that here, we formulate a Bayes Factor for genes at TWAS risk regions. The Bayes Factor for each gene quantifies the amount of evidence in favor of the causal model (ith gene drives risk) versus the null (ith gene has no causal effect). We extend individual Bayes Factors for k genes at a PrCa risk region to compute the posterior probability that a gene is causal as,

$$\Pr \left( {{\mathrm{gene}}\,i\,{\mathrm{is}}\,{\mathrm{causal}}|{\mathbf{{z}}}_{{\mathrm{TWAS}}},n\sigma _\alpha ^2} \right) = \frac{{{\rm BF}_{{i}}}}{{\mathop {\sum }\nolimits_{\mathrm{k}} {\rm {BF}}_{\mathrm{k}}}}.$$

Equipped with our definition of posterior probability for each gene being causal, we define ρ-credible gene sets for a PrCa risk region. Formally, a set of indices $i \in I$ defines a ρ-credible gene set if

$\rho = \mathop {\sum}\limits_{i \in I} {\Pr \left( {{\mathrm{gene}}\,i\,{\mathrm{is}}\,{\mathrm{causal|}}{\mathbf{{z}}}_{{\mathrm{TWAS}}},\,n\sigma _\alpha ^2} \right).}$

For a fixed ρ we optimize over k genes at a region by greedily adding genes until the total density is at least ρ.

To ensure that our ρ-credible sets are well-calibrated we performed simulations by predicting expression levels into 489 samples of European ancestry from 1000 Genomes⁵³ and estimating the local correlation structure to sample TWAS Z-scores directly (Supplementary Note 1). Under the assumption of a single causal gene at a risk region, we sampled TWAS Z-scores for 1000 independent regions. We then performed Bayesian prioritization at each region and computed ρ-credible sets for various levels of ρ while counting the proportion of causal genes identified across all simulations.

Pathway analyses

To determine which pathways may be enriched with genes identified from our Bayesian prioritization approach, we used the R package GOseq⁵⁷ which internally links gene identifiers to GO terms (GO db: 2017-09-02). We categorized all 16,389 genes into prioritized/not-prioritized and ran the analysis using custom R scripts linking GOseq. GOseq obtains P-values for overrepresented genes using the Wallenius approximation to the non-central hypergeometric distribution. We limited analysis to Gene Ontology Biological Pathways (GO:BP). GOSeq drops genes without GO annotations from analysis. We observed 4711 genes dropped from analyses resulting in 11,678 genes put forward for enrichment tests (Supplementary Data 8; Supplementary Note 1).

URLs

1000 Genomes Phase3: http://www.internationalgenome.org/

Fire Hose v2016_1_28: http://gdac.broadinstitute.org/

FUSION: http://gusevlab.org/projects/fusion/

GCTA v1.26: http://cnsgenomics.com/software/gcta/

GEMMA v0.94: http://www.xzlab.org/software.html

GOseq v1.26: http://bioinf.wehi.edu.au/software/goseq/

MapSplice v2: http://www.netlab.uky.edu/p/bioinfo/MapSplice2

PLINK v1.9: https://www.cog-genomics.org/plink2/

OncoArray: https://epi.grants.cancer.gov/oncoarray/

Data availability

Complete TWAS and fine-mapping results are available at http://github.com/bogdanlab/prca_twas/. OncoArray PrCa GWAS summary data used in this study are available at http://practical.icr.ac.uk/blog/. Relevant TCGA data are available from Broad Firehouse at http://gdac.broadinstitute.org. FUSION software, weights/models, and reference LD are available at http://gusevlab.org/projects/fusion/.

Change history

08 January 2019
The original version of this Article contained an error in the spelling of a member of the PRACTICAL Consortium, Manuela Gago-Dominguez, which was incorrectly given as Manuela Gago Dominguez. This has now been corrected in both the PDF and HTML versions of the Article. Furthermore, In the original HTML version of this Article, the order of authors within the author list was incorrect. The consortium PRACTICAL consortium was incorrectly listed after Bogdan Pasaniuc and should have been listed after Kathryn L. Penney. This error has been corrected in the HTML version of the Article; the PDF version was correct at the time of publication.

References

Hjelmborg, J. B. et al. The heritability of prostate cancer in the nordic twin study of cancer. Cancer Epidemiol. Biomark. 23, 2303 (2014).
Article Google Scholar
Mucci, L. A. et al. Familial risk and heritability of cancer among twins in nordic countries. JAMA 315, 68–76 (2016).
Article CAS Google Scholar
Eeles, R. A. et al. Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping array. Nat. Genet. 45, 385–391 (2013).
Article CAS Google Scholar
Amin Al Olama et al. Multiple novel prostate cancer susceptibility signals identified by fine-mapping of known risk loci among Europeans. Human Mol. Genet. 24, 5589–5602 (2015).
Article Google Scholar
Al Olama, A. A. et al. A meta-analysis of 87,040 individuals identifies 23 new susceptibility loci for prostate cancer. Nat. Genet. 46, 1103–1109 (2014).
Article CAS Google Scholar
Al Olama, A. A. et al. Multiple loci on 8q24 associated with prostate cancer susceptibility. Nat. Genet. 41, 1058–1060 (2009).
Article CAS Google Scholar
Eeles, R. A. et al. Multiple newly identified loci associated with prostate cancer susceptibility. Nat. Genet. 40, 316–321 (2008).
Article CAS Google Scholar
Spisak, S. et al. CAUSEL: an epigenome- and genome-editing pipeline for establishing function of noncoding GWAS variants. Nat. Med. 21, 1357–1363 (2015).
Article CAS Google Scholar
Nicolae, D. L. et al. Trait-associated SNPs are more likely to be eQTLs: annotation to enhance discovery from GWAS. PLoS Genet. 6, e1000888 (2010).
Article Google Scholar
Maurano, M. T. et al. Systematic localization of common disease-associated variation in regulatory DNA. Science 337, 1190–1195 (2012).
Article ADS CAS Google Scholar
Roadmap Epigenomics, C. et al. Integrative analysis of 111 reference human epigenomes. Nature 518, 317–330 (2015).
Article Google Scholar
Hazelett, D. J. et al. Comprehensive functional annotation of 77 prostate cancer risk loci. PLoS Genet. 10, e1004102 (2014).
Article Google Scholar
Gusev, A. et al. Atlas of prostate cancer heritability in European and African-American men pinpoints tissue-specific regulation. Nat. Commun. 7, 10979 (2016).
Article ADS CAS Google Scholar
Thibodeau, S. N. et al. Identification of candidate genes for prostate cancer-risk SNPs utilizing a normal prostate tissue eQTL data set. Nat. Commun. 6, 8653 (2015).
Article CAS Google Scholar
Whitington, T. et al. Gene regulatory mechanisms underpinning prostate cancer susceptibility. Nat. Genet. 48, 387–397 (2016).
Article CAS Google Scholar
Penney, K. L. et al. Association of prostate cancer risk variants with gene expression in normal and tumor tissue. Cancer Epidemiol. Biomark. 24, 255 (2015).
Article CAS Google Scholar
Li, Q. et al. Expression QTL-based analyses reveal candidate causal genes and loci across five tumor types. Human Mol. Genet. 23, 5294–5302 (2014).
Article CAS Google Scholar
Grisanzio, C. et al. Genetic and functional analyses implicate the NUDT11, HNF1B, and SLC22A3 genes in prostate cancer pathogenesis. Proc. Natl Acad. Sci. USA 109, 11252–11257 (2012).
Article ADS CAS Google Scholar
Xu, X. et al. Variants at IRX4 as prostate cancer expression quantitative trait loci. Eur. J. Hum. Genet. 22, 558–563 (2014).
Article CAS Google Scholar
Huang, Q. et al. A prostate cancer susceptibility allele at 6q22 increases RFX6 expression by modulating HOXB13 chromatin binding. Nat. Genet. 46, 126–135 (2014).
Article CAS Google Scholar
Pomerantz, M. M. et al. Analysis of the 10q11 cancer risk locus implicates MSMB and NCOA4 in human prostate tumorigenesis. PLoS Genet. 6, e1001204 (2010).
Article Google Scholar
Lonsdale, J. et al. The genotype-tissue expression (GTEx) project. Nat. Genet. 45, 580–585 (2013).
Article CAS Google Scholar
Chun, S. et al. Limited statistical evidence for shared genetic effects of eQTLs and autoimmune-disease-associated loci in three major immune-cell types. Nat. Genet. 49, 600 (2017).
Article CAS Google Scholar
Gamazon, E. R. et al. A gene-based association method for mapping traits using reference transcriptome data. Nat. Genet. 47, 1091–1098 (2015).
Article CAS Google Scholar
Gusev, A. et al. Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet. 48, 245–52 (2016).
Article CAS Google Scholar
Zhu, Z. et al. Integration of summary data from GWAS and eQTL studies predicts complex trait gene targets. Nat. Genet. 48, 481 (2016).
Article CAS Google Scholar
Mancuso, N. et al. Integrating gene expression with summary association statistics to identify genes associated with 30 complex traits. Am. J. Human Genet. 100, 473–487 (2017).
Article CAS Google Scholar
Pavlides, J. M. W. et al. Predicting gene targets from integrative analyses of summary data from GWAS and eQTL studies for 28 human complex traits. Genome Med. 8, 1–6 (2016).
Article Google Scholar
Schumacher, F. R. et al. Prostate cancer meta-analysis of more than 140,000 men identifies 63 novel prostate cancer susceptibility loci. Nat. Genet. 50, 928–936 (2017).
Article Google Scholar
Fromer, M. et al. Gene expression elucidates functional impact of polygenic risk for schizophrenia. Nat. Neurosci. 19, 1442–1453 (2016).
Article CAS Google Scholar
Raitakari, O. T. et al. Cohort profile: the cardiovascular risk in young finns study. Int. J. Epidemiol. 37, 1220–1226 (2008).
Article Google Scholar
Stančáková, A. et al. Hyperglycemia and a common variant of GCKR are associated with the levels of eight amino acids in 9,369 finnish men. Diabetes 61, 1895–1902 (2012).
Article Google Scholar
Stančáková, A. et al. Changes in insulin sensitivity and insulin release in relation to glycemia and glucose tolerance in 6,414 finnish men. Diabetes 58, 1212–1221 (2009).
Article Google Scholar
Nuotio, J. et al. Cardiovascular risk factors in 2011 and secular trends since 2007: the Cardiovascular Risk in Young Finns Study. Scand. J. Public Health 42, 563–571 (2014).
Article Google Scholar
Wright, F. A. et al. Heritability and genomics of gene expression in peripheral blood. Nat. Genet. 46, 430–437 (2014).
Article CAS Google Scholar
The Cancer Genome Atlas Research Network et al. The cancer genome atlas pan-cancer analysis project. Nat. Genet. 45, 1113–1120 (2013).
Article Google Scholar
de los Campos, G., Vazquez, A. I., Fernando, R., Klimentidis, Y. C. & Sorensen, D. Prediction of complex human traits using the genomic best linear unbiased predictor. PLoS Genet. 9, e1003608 (2013).
Article Google Scholar
Zhou, X., Carbonetto, P. & Stephens, M. Polygenic modeling with bayesian sparse linear mixed models. PLoS Genet. 9, e1003264 (2013).
Article CAS Google Scholar
Abeshouse, A. et al. The molecular taxonomy of primary prostate cancer. Cell 163, 1011–1025 (2015).
Article CAS Google Scholar
Matesic, L. E. et al. Mutations in Mlph, encoding a member of the Rab effector family, cause the melanosome transport defects observed in leaden mice. Proc. Natl Acad. Sci. USA 98, 10238–10243 (2001).
Article ADS CAS Google Scholar
Li, Y. I. et al. RNA splicing is a primary link between genetic variation and disease. Science 352, 600 (2016).
Article ADS CAS Google Scholar
Darido, C. et al. Targeting of the tumor suppressor GRHL3 by a miR-21-dependent proto-oncogenic network results in PTEN loss and tumorigenesis. Cancer Cell 20, 635–648 (2011).
Article CAS Google Scholar
Yang, J. et al. Altered DNA polymerase iota expression in breast cancer cells leads to a reduction in DNA replication fidelity and a higher rate of mutagenesis. Cancer Res. 64, 5597–607 (2004).
Article CAS Google Scholar
Yuan, F. et al. Overexpressed DNA polymerase iota regulated by JNK/c-Jun contributes to hypermutagenesis in bladder cancer. PLoS ONE 8, e69317 (2013).
Article ADS CAS Google Scholar
Bu, H. et al. Putative prostate cancer risk SNP in an androgen receptor-binding site of the melanophilin gene illustrates enrichment of risk SNPs in androgen receptor target sites. Human Mutat. 37, 52–64 (2016).
Article CAS Google Scholar
Gutierrez-Arcelus, M. et al. Tissue-specific effects of genetic and epigenetic variation on gene regulation and splicing. PLoS Genet. 11, e1004958 (2015).
Article Google Scholar
Verhaagh, S., Schweifer, N., Barlow, D. P. & Zwart, R. Cloning of the mouse and human solute carrier 22a3 (Slc22a3/SLC22A3) identifies a conserved cluster of three organic cation transporters on mouse chromosome 17 and human 6q26–q27. Genomics 55, 209–218 (1999).
Article CAS Google Scholar
Siegel, R. L., Miller, K. D. & Jemal, A. Cancer statistics, 2016. Cancer J. Clin. 66, 7–30 (2016).
Article Google Scholar
Sutcliffe, S., De Marzo, A. M., Sfanos, K. S. & Laurence, M. MSMB variation and prostate cancer risk: clues towards a possible fungal etiology. Prostate 74, 569–578 (2014).
Article CAS Google Scholar
Gusev, A. et al. Transcriptome-wide association study of schizophrenia and chromatin activity yields mechanistic disease insights. Nat. Genet. 50, 538–548 (2018).
Article CAS Google Scholar
McCarthy, S. et al. A reference panel of 64,976 haplotypes for genotype imputation. Nat. Genet. 48, 1279–1283 (2016).
Article CAS Google Scholar
Consortium, T.I.H. Integrating common and rare genetic variation in diverse human populations. Nature 467, 52–58 (2010).
Article ADS Google Scholar
The Genomes Project, C. A global reference for human genetic variation. Nature 526, 68–74 (2015).
Article ADS Google Scholar
Hormozdiari, F., Kichaev, G., Yang, W.-Y., Pasaniuc, B. & Eskin, E. Identification of causal genes for complex traits. Bioinformatics 31, i206–i213 (2015).
Article CAS Google Scholar
Maller, J. B. et al. Bayesian refinement of association signals for 14 loci in 3 common diseases. Nat. Genet. 44, 1294–1301 (2012).
Article CAS Google Scholar
Chen, W. et al. Fine mapping causal variants with an approximate bayesian method using marginal test statistics. Genetics 200, 719 (2015).
Article Google Scholar
Young, M. D., Wakefield, M. J., Smyth, G. K. & Oshlack, A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11, R14 (2010).
Article Google Scholar

Download references

Acknowledgements

We wish to thank all GWAS study groups contributing to the meta-analysis dataset from which the transcriptome-wide association analyses were conducted: BPC3 (Breast and Prostate Cancer Cohort Consortium); CAPS (Cancer of the Prostate in Sweden); PEGASUS (Prostate Cancer Genome-wide Association Study of Uncommon Susceptibility Loci); APCB BioResource (Australian Prostate Cancer BioResource); and The PRACTICAL (Prostate Cancer Association Group to Investigate Cancer-Associated Alterations in the Genome) Consortium. This work was supported by NIH grants R01-HG009120, R01-HG006399, and U01-CA194393. This work was supported by the Canadian Institutes of Health Research, European Commission's Seventh Framework Programme grant agreement n° 223175 (HEALTH-F2-2009-223175), Cancer Research UK Grants C5047/A7357, C1287/A10118, C1287/A16563, C5047/A3354, C5047/A10692, C16913/A6135, and The National Institute of Health (NIH) Cancer Post-Cancer GWAS initiative grant: No. 1U19 CA 148537-01 (the GAME-ON initiative). We thank the following for funding support: The Institute of Cancer Research and The Everyman Campaign, The Prostate Cancer Research Foundation, Prostate Research Campaign UK (now Prostate Action), The Orchid Cancer Appeal, The National Cancer Research Network UK, The National Cancer Research Institute (NCRI) UK. We are grateful for support of NIHR funding to the NIHR Biomedical Research Centre at The Institute of Cancer Research and The Royal Marsden NHS Foundation Trust and the NIHR Biomedical Research Centre at the University of Cambridge. The Prostate Cancer Program of Cancer Council Victoria also acknowledge grant support from The National Health and Medical Research Council, Australia (126402, 209057, 251533, 396414, 450104, 504700, 504702, 504715, 623204, 940394, and 614296), VicHealth, Cancer Council Victoria, The Prostate Cancer Foundation of Australia, The Whitten Foundation, Price Waterhouse Coopers, and Tattersall’s. Genotyping of the OncoArray was funded by the US National Institutes of Health (NIH) [U19 CA 148537 for ELucidating Loci Involved in Prostate cancer SuscEptibility (ELLIPSE) project and X01HG007492 to the Center for Inherited Disease Research (CIDR) under contract number HHSN268201200008I]. Additional analytic support was provided by NIH NCI U01 CA188392. Funding for the iCOGS infrastructure came from: the European Community's Seventh Framework Programme under grant agreement n° 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A 10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692, and C8197/A16565), the National Institutes of Health (CA128978) and Post-Cancer GWAS initiative (1U19 CA148537, 1U19 CA148065 and 1U19 CA148112—the GAME-ON initiative), the Department of Defence (W81XWH-10-1-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, Komen Foundation for the Cure, the Breast Cancer Research Foundation, and the Ovarian Cancer Research Fund. The BPC3 was supported by the US National Institutes of Health, National Cancer Institute (cooperative agreements U01-CA98233 to D.J.H., U01-CA98710 to S.M.G., U01-CA98216 to E.R., and U01-CA98758 to B.E.H., and Intramural Research Program of NIH/National Cancer Institute, Division of Cancer Epidemiology and Genetics). CAPS GWAS study was supported by the Swedish Cancer Foundation (grant no 09-0677, 11-484, 12-823), the Cancer Risk Prediction Center (CRisP; www.crispcenter.org), a Linneus Centre (Contract ID 70867902) financed by the Swedish Research Council, Swedish Research Council (grant no K2010-70X-20430-04-3, 2014-2269). PEGASUS was supported by the Intramural Research Program, Division of Cancer Epidemiology and Genetics, National Cancer Institute, National Institutes of Health.

Author information

Authors and Affiliations

Department of Pathology and Laboratory Medicine, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, 90095, CA, USA
Nicholas Mancuso & Bogdan Pasaniuc
The Center for Bioinformatics and Functional Genomics, Cedars-Sinai Medical Center, Los Angeles, 90048, CA, USA
Simon Gayther
Dana Farber Cancer Institute, Boston, 02215, MA, USA
Alexander Gusev
Division of Epidemiology, Department of Medicine, Vanderbilt Epidemiology Center, Vanderbilt-Ingram Cancer Center, Vanderbilt University School of Medicine, Nashville, 37232, TN, USA
Wei Zheng
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, 02115, MA, USA
Kathryn L. Penney & Lorelei Mucci
Channing Division of Network Medicine, Department of Medicine, Brigham and Women’s Hospital/Harvard Medical School, Boston, 02115, MA, USA
Kathryn L. Penney
Division of Genetics and Epidemiology, The Institute of Cancer Research, London, SW7 3RP, UK
Zsofia Kote-Jarai, Rosalind Eeles & Sara Benlloch
Royal Marsden NHS Foundation Trust, London, SW3 6JJ, UK
Zsofia Kote-Jarai & Rosalind Eeles
Department of Medical Oncology, Dana-Farber Cancer Institute and Harvard Medical School, Boston, 02215, MA, USA
Matthew Freedman
Department of Preventive Medicine, Norris Comprehensive Cancer Center, Keck School of Medicine, University of Southern California, Los Angeles, 90015, CA, USA
Christopher Haiman, Brian E. Henderson, David V. Conti & Sue Ann Ingles
Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, 90095, CA, USA
Bogdan Pasaniuc
Bioinformatics Interdepartmental Program, University of California, Los Angeles, Los Angeles, 90095, CA, USA
Bogdan Pasaniuc
Centre for Cancer Genetic Epidemiology, Department of Public Health and Primary Care, University of Cambridge, Strangeways Research Laboratory, Cambridge, CB1 8RN, UK
Sara Benlloch & Ali Amin Al Olama
Department of Epidemiology and Biostatistics, Case Western Reserve University, Cleveland, 44106-7219, OH, USA
Fredrick R. Schumacher
Seidman Cancer Center, University Hospitals, Cleveland, 44106, OH, USA
Fredrick R. Schumacher
University of Cambridge, Department of Clinical Neurosciences, Cambridge, CB2 0QQ, UK
Ali Amin Al Olama
Institute of Population Health, University of Manchester, Manchester, M13 9PL, UK
Kenneth Muir
Warwick Medical School, University of Warwick, Coventry, CV4 7AL, UK
Kenneth Muir
Division of Cancer Epidemiology and Genetics, National Cancer Institute, NIH, Bethesda, 21701, MD, USA
Sonja I. Berndt, Stephen Chanock, Demetrius Albanes, Stephanie Weinstein & Stella Koutros
Department of Medical Epidemiology and Biostatistics, Karolinska Institute, Stockholm, SE-171 77, Sweden
Fredrik Wiklund & Henrik Gronberg
Epidemiology Research Program, American Cancer Society, 250 Williams Street, Atlanta, 30303, GA, USA
Victoria L. Stevens
SWOG Statistical Center, Fred Hutchinson Cancer Research Center, Seattle, 98109-1024, WA, USA
Catherine M. Tangen
Australian Prostate Cancer Research Centre-Qld, Institute of Health and Biomedical Innovation and School of Biomedical Science, Queensland University of Technology, Brisbane, 4059, Queensland, Australia
Jyotsna Batra & Judith Clements
Translational Research Institute, Brisbane, 4102, Queensland, Australia
Jyotsna Batra & Judith Clements
University College London, Department of Applied Health Research, London, WC1E 7HB, UK
Nora Pashayan
Centre for Cancer Genetic Epidemiology, Department of Oncology, University of Cambridge, Strangeways Laboratory, Cambridge, WC1E 7HB, UK
Nora Pashayan
Department of Medical Biochemistry and Genetics, Institute of Biomedicine, University of Turku, Turku, FI-20014, Finland
Johanna Schleutker
Tyks Microbiology and Genetics, Department of Medical Genetics, Turku University Hospital, Hospital, 20521, Finland
Johanna Schleutker
BioMediTech, University of Tampere, Tampere, FI-33014, Finland
Johanna Schleutker
Division of Nutritional Epidemiology, Institute of Environmental Medicine, Karolinska Institutet, SE-171 77, Sweden
Alicja Wolk
Institute of Cancer Sciences, University of Manchester, Manchester Academic Health Science Centre, Radiotherapy Related Research, The Christie Hospital NHS Foundation Trust, Manchester, M13 9PL, UK
Catharine West
CeRePP, Pitie-Salpetriere Hospital, Paris, F-75020, France
Géraldine Cancel-Tassin
UPMC Univ Paris 06, GRC N°5 ONCOTYPE-URO, CeRePP, Tenon Hospital, Paris, F-75020, France
Géraldine Cancel-Tassin
Department of Molecular Medicine, Aarhus University Hospital, Aarhus N, 8200, Denmark
Karina Dalsgaard Sorensen
Department of Clinical Medicine, Aarhus University, Aarhus N, 8200, Denmark
Karina Dalsgaard Sorensen
Department of Medical Genetics, Oslo University Hospital, Oslo, 0424, Norway
Lovise Maehle
University of Cambridge, Department of Oncology, Addenbrooke’s Hospital, Cambridge, CB2 0QQ, UK
David E. Neal
Cancer Research UK Cambridge Research Institute, Li Ka Shing Centre, Cambridge, CB2 0RE, UK
David E. Neal
Nuffield Department of Surgical Sciences, University of Oxford, Oxford, OX1 2JD, UK
Freddie C. Hamdy
Faculty of Medical Science, University of Oxford, John Radcliffe Hospital, Oxford, OX1 2JD, UK
Freddie C. Hamdy
School of Social and Community Medicine, University of Bristol, Bristol, BS8 2PS, UK
Jenny L. Donovan
Cancer Epidemiology, Nuffield Department of Population Health University of Oxford, Oxford, OX3 7LF, UK
Ruth C. Travis
Department of Surgical Oncology, Princess Margaret Cancer Centre, Toronto, M5G 2M9, Canada
Robert J. Hamilton
Department of Radiation Oncology, Icahn School of Medicine at Mount Sinai, New York, 10029, NY, USA
Barry Rosenstein
Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, 10029-5674, NY, USA
Barry Rosenstein
Centre for Molecular Oncology, Barts Cancer Institute, Queen Mary University of London, John Vane Science Centre, London, EC1M 6BQ, UK
Yong-Jie Lu
Cancer Epidemiology and Intelligence Division, The Cancer Council Victoria, Melbourne, Victoria, 3004, Australia
Graham G. Giles
Centre for Epidemiology and Biostatistics, Melbourne School of Population and Global Health, The University of Melbourne, Melbourne, VIC 3010, Australia
Graham G. Giles
Division of Urologic Surgery, Brigham and Womens Hospital, Boston, 02115, MA, USA
Adam S. Kibel
Fundacion Publica Galega de Medicina Xenomica-SERGAS, Grupo de Medicina Xenomica, CIBERER, IDIS, Santiago de Compostela, 15706, Spain
Ana Vega
Centre for Research in Environmental Epidemiology (CREAL), Barcelona Institute for Global Health (ISGlobal), Barcelona, 08003, Spain
Manolis Kogevinas
CIBER Epidemiologia y Salud Publica (CIBERESP), Madrid, 28029, Spain
Manolis Kogevinas
IMIM (Hospital del Mar Research Institute), Barcelona, 08003, Spain
Manolis Kogevinas
Universitat Pompeu Fabra (UPF), Barcelona, 08002, Spain
Manolis Kogevinas
Department of Cancer Epidemiology, Moffitt Cancer Center, Tampa, 33612, USA
Jong Y. Park
Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, Washington, 98109-1024, USA
Janet L. Stanford & Lisa F. Newcomb
Department of Epidemiology, School of Public Health, University of Washington, Seattle, Washington, 98195, USA
Janet L. Stanford
International Hereditary Cancer Center, Department of Genetics and Pathology, Pomeranian Medical University, Szczecin, 70-115, Poland
Cezary Cybulski
Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, 2200, Denmark
Børge G. Nordestgaard
Department of Clinical Biochemistry, Herlev and Gentofte Hospital, Copenhagen University Hospital, Herlev, 2200, Denmark
Børge G. Nordestgaard
Division of Clinical Epidemiology and Aging Research, German Cancer Research Center (DKFZ), Heidelberg, D-69120, Germany
Hermann Brenner
German Cancer Consortium (DKTK), German Cancer Research Center (DKFZ), Heidelberg, D-69120, Germany
Hermann Brenner
Division of Preventive Oncology, German Cancer Research Center (DKFZ) and National Center for Tumor Diseases (NCT), Heidelberg, D-69120, Germany
Hermann Brenner
Institute for Human Genetics, University Hospital Ulm, Ulm, 89075, Germany
Christiane Maier
The University of Texas M. D. Anderson Cancer Center, Department of Genitourinary Medical Oncology, Houston, 77030, TX, USA
Jeri Kim
Cancer Prevention Institute of California, Fremont, 94538, CA, USA
Esther M. John
Department of Health Research and Policy (Epidemiology) and Stanford Cancer Institute, Stanford University School of Medicine, Stanford, 94305-5101, CA, USA
Esther M. John
Department of Genetics, Portuguese Oncology Institute of Porto, Porto, 4200-072, Portugal
Manuel R. Teixeira
Biomedical Sciences Institute (ICBAS), University of Porto, Porto, 4050-313, Portugal
Manuel R. Teixeira
Department of Population Sciences, Beckman Research Institute of the City of Hope, Duarte, 91010, CA, USA
Susan L. Neuhausen
Ghent University, Faculty of Medicine and Health Sciences, Basic Medical Sciences, Gent, B-9000, Belgium
Kim De Ruyck
Department of Surgery, Faculty of Medicine, University of Malaya, Kuala Lumpur, 50603, Malaysia
Azad Razack
Department of Urology, University of Washington, Seattle, 98195, WA, USA
Lisa F. Newcomb
Institute of Human Genetics, University Medical Center Hamburg-Eppendorf, Hamburg, D-20246, Germany
Davor Lessel
Molecular Medicine Center, Department of Medical Chemistry and Biochemistry, Medical University, Sofia, 1431, Bulgaria
Radka Kaneva
Department of Oncology, Cross Cancer Institute, University of Alberta, Edmonton, AB T6G 1Z2, Alberta, Canada
Nawaid Usmani
Division of Radiation Oncology, Cross Cancer Institute, Edmonton, AB T6G 1Z2, Alberta, Canada
Nawaid Usmani
Molecular Endocrinology Laboratory, Department of Cellular and Molecular Medicine, KU Leuven, BE-3000, Leuven, Belgium
Frank Claessens
Institute of Cancer Sciences, Manchester Cancer Research Centre, University of Manchester, Manchester Academic Health Science Centre, St Mary’s Hospital, Manchester, M13 9WL, UK
Paul A. Townsend
Genomic Medicine Group, Galician Foundation of Genomic Medicine, Instituto de Investigacion Sanitaria de Santiago de Compostela (IDIS), Complejo Hospitalario Universitario de Santiago, Servicio Galego de Saude, SERGAS, Santiago De Compostela, 15706, Spain
Manuela Gago-Dominguez
University of California San Diego, Moores Cancer Center, La Jolla, 92037, CA, USA
Manuela Gago-Dominguez
Department of Urology, Erasmus University Medical Center, Rotterdam, 3015 CE, The Netherlands
Monique J. Roobol
Cancer and Environment Group, Center for Research in Epidemiology and Population Health (CESP), INSERM, University Paris-Sud, University Paris-Saclay, Villejuif, 94807, France
Florence Menegaux
Clinical Gerontology Unit, University of Cambridge, Cambridge, CB2 2QQ, UK
Kay-Tee Khaw
Division of Genetic Epidemiology, Department of Medicine, University of Utah School of Medicine, Salt Lake City, 84112, Utah, USA
Lisa Cannon-Albright
George E. Wahlen Department of Veterans Affairs Medical Center, Salt Lake City, 84148, UT, USA
Lisa Cannon-Albright
The University of Surrey, Guildford, GU2 7XH, Surrey, UK
Hardev Pandha
Department of Laboratory Medicine and Pathology, Mayo Clinic, Rochester, 55905, MN, USA
Stephen N. Thibodeau
Program in Genetic Epidemiology and Statistical Genetics, Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, 02115, MA, USA
David J. Hunter & Peter Kraft

Authors

Nicholas Mancuso
View author publications
You can also search for this author in PubMed Google Scholar
Simon Gayther
View author publications
You can also search for this author in PubMed Google Scholar
Alexander Gusev
View author publications
You can also search for this author in PubMed Google Scholar
Wei Zheng
View author publications
You can also search for this author in PubMed Google Scholar
Kathryn L. Penney
View author publications
You can also search for this author in PubMed Google Scholar
Zsofia Kote-Jarai
View author publications
You can also search for this author in PubMed Google Scholar
Rosalind Eeles
View author publications
You can also search for this author in PubMed Google Scholar
Matthew Freedman
View author publications
You can also search for this author in PubMed Google Scholar
Christopher Haiman
View author publications
You can also search for this author in PubMed Google Scholar
Bogdan Pasaniuc
View author publications
You can also search for this author in PubMed Google Scholar

Consortia

The PRACTICAL consortium

Brian E. Henderson
, Sara Benlloch
, Fredrick R. Schumacher
, Ali Amin Al Olama
, Kenneth Muir
, Sonja I. Berndt
, David V. Conti
, Fredrik Wiklund
, Stephen Chanock
, Victoria L. Stevens
, Catherine M. Tangen
, Jyotsna Batra
, Judith Clements
, Henrik Gronberg
, Nora Pashayan
, Johanna Schleutker
, Demetrius Albanes
, Stephanie Weinstein
, Alicja Wolk
, Catharine West
, Lorelei Mucci
, Géraldine Cancel-Tassin
, Stella Koutros
, Karina Dalsgaard Sorensen
, Lovise Maehle
, David E. Neal
, Freddie C. Hamdy
, Jenny L. Donovan
, Ruth C. Travis
, Robert J. Hamilton
, Sue Ann Ingles
, Barry Rosenstein
, Yong-Jie Lu
, Graham G. Giles
, Adam S. Kibel
, Ana Vega
, Manolis Kogevinas
, Jong Y. Park
, Janet L. Stanford
, Cezary Cybulski
, Børge G. Nordestgaard
, Hermann Brenner
, Christiane Maier
, Jeri Kim
, Esther M. John
, Manuel R. Teixeira
, Susan L. Neuhausen
, Kim De Ruyck
, Azad Razack
, Lisa F. Newcomb
, Davor Lessel
, Radka Kaneva
, Nawaid Usmani
, Frank Claessens
, Paul A. Townsend
, Manuela Gago-Dominguez
, Monique J. Roobol
, Florence Menegaux
, Kay-Tee Khaw
, Lisa Cannon-Albright
, Hardev Pandha
, Stephen N. Thibodeau
, David J. Hunter
& Peter Kraft

Contributions

N.M., S.G., A.G., W.Z., K.L.P., Z.K., R.E., M.F., C.H., and B.P. planned the study. N.M. and A.G. performed primary analyses. PRACTICAL collected data. N.M. and B.P. wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Nicholas Mancuso.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

A full list of consortium members appears at the end of the paper.

Electronic supplementary material

Supplementary Information

Description of Additional Supplementary Files

Supplementary Data 1

Supplementary Data 2

Supplementary Data 3

Supplementary Data 4

Supplementary Data 5

Supplementary Data 6

Supplementary Data 7

Supplementary Data 8

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Mancuso, N., Gayther, S., Gusev, A. et al. Large-scale transcriptome-wide association study identifies new prostate cancer risk regions. Nat Commun 9, 4079 (2018). https://doi.org/10.1038/s41467-018-06302-1

Download citation

Received: 11 June 2018
Accepted: 28 August 2018
Published: 04 October 2018
DOI: https://doi.org/10.1038/s41467-018-06302-1

This article is cited by

Large-scale integrative analysis of juvenile idiopathic arthritis for new insight into its pathogenesis
- Daeun Kim
- Jaeseung Song
- Wonhee Jang
Arthritis Research & Therapy (2024)
Prioritizing susceptibility genes for the prognosis of male-pattern baldness with transcriptome-wide association study
- Eunyoung Choi
- Jaeseung Song
- Wonhee Jang
Human Genomics (2024)
Association between night work and prostate cancer: a systematic review and meta-analysis
- Alexander Jahn
- Mathilde Lumbye Nielsen
- Annett Dalbøge
International Archives of Occupational and Environmental Health (2024)
Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets
- Xiaoguang Xu
- Chachrit Khunsriraksakul
- Maciej Tomaszewski
Nature Communications (2024)
Transcriptome-wide association study identifies new susceptibility genes and pathways for spondyloarthritis
- Xiaochen Su
- Anfa Chen
- Yingang Zhang
Journal of Orthopaedic Surgery and Research (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.