Skip to main content
Download PDF

Contribution and clinical relevance of germline variation to the cancer transcriptome

BMC Cancer volume 22, Article number: 675 (2022) Cite this article

Abstract

Background

Somatic alterations in the cancer genome, some of which are associated with changes in gene expression, have been characterized in multiple studies across diverse cancer types. However, less is known about germline variants that influence tumor biology by shaping the cancer transcriptome.

Methods

We performed expression quantitative trait loci (eQTL) analyses using multi-dimensional data from The Cancer Genome Atlas to explore the role of germline variation in mediating the cancer transcriptome. After accounting for associations between somatic alterations and gene expression, we determined the contribution of inherited variants to the cancer transcriptome relative to that of somatic variants. Finally, we performed an interaction analysis using estimates of tumor cellularity to identify cell type-restricted eQTLs.

Results

The proportion of genes with at least one eQTL varied between cancer types, ranging between 0.8% in melanoma to 28.5% in thyroid cancer and was correlated more strongly with intratumor heterogeneity than with somatic alteration rates. Although contributions to variance in gene expression was low for most genes, some eQTLs accounted for more than 30% of expression of proximal genes. We identified cell type-restricted eQTLs in genes known to be cancer drivers including LPP and EZH2 that were associated with disease-specific mortality in TCGA but not associated with disease risk in published GWAS. Together, our results highlight the need to consider germline variation in interpreting cancer biology beyond risk prediction.

Peer Review reports

Background

Deregulated gene expression is a defining feature of the cancer cell and often results in the disruption of key signaling pathways that control cellular growth and proliferation [1]. Analyses of multidimensional data from large-scale projects such as The Cancer Genome Atlas (TCGA) have demonstrated important associations between somatic genetic alterations and changes in gene expression, some of which represent oncogenic alterations in cancer driver genes [2, 3].

In contrast, the role of inherited polymorphisms in influencing the cancer transcriptome has been less well studied. Identification of associations between germline single nucleotide polymorphisms (SNP) and gene expression has been established as a strategy to understand the mechanisms by which inherited variants may influence defined phenotypes, including cancer incidence [4]. However, there have been few studies addressing the biological and clinical relevance of eQTLs in the specific context of the malignant transcriptome.

Tumor biopsies consist of numerous cell populations including immune and stromal cells. As a result, eQTLs identified in tumor samples may be derived from either malignant or non-malignant cell types. This was demonstrated in an analysis of 24 cancer types in which eQTL mapping revealed that genes in which an eQTL was detected (eGenes) were enriched for ontology terms related to immune function [5]. Approaches to identify cell type-restricted eQTLs include adapting the standard additive eQTL model by considering how an eQTL’s effect varies by the estimated fraction of a given cell type in a tissue biopsy [6].

Application of this approach to breast cancer datasets using tumor cellularity as an estimate of tumor cell fraction revealed that only a few of the eQTLs detected in these studies were likely to be tumor cell-specific [7]. More recently, analyses using in silico deconvolution methods have revealed the presence of cell type-restricted eQTLs in at least 3000 genes in the Genotype-Tissue Expression (GTEx) project [8].

In this study, we explored the role of inherited variants in the cancer transcriptome using genotype, somatic alteration and gene expression data from 24 cancer types from the TCGA project. We studied the tissue specificity of eQTLs in cancer, with a focus on characterizing putative functional enrichment in genes with at least one proximal eQTL (eGenes). We quantified the contribution of eQTLs to variation in gene expression relative to somatic alterations and determined that the role and potential relevance of cancer eQTLs likely varies by cancer type. In particular, thyroid carcinoma and prostate adenocarcinoma had the highest number of detected eQTLs across the study, likely driven by the low somatic alteration rates and low levels of intratumor heterogeneity within these cancer types. We also identified cell type-restricted eQTLs that are associated with patient prognosis, demonstrating the potential impact of inherited variants to mediating tumor biology and clinical outcome.

Methods

Data sources

We downloaded genotype, mutation and expression datasets from the TCGA from the Genomics Data Commons [4, 9, 10] (GDC). Data stored on the GDC have been processed using standardized pipelines with the hg38 reference genome build. We obtained the clinical data for subjects with samples in TCGA that had already been curated as previously described [11]. Raw genotype data were downloaded on November 28, 2017 and the version of the RNA-seq data used for the final analyses were accessed on October 18, 2019. All somatic alteration data used for the study were obtained on March 13, 2019.

Genotype processing

We processed genotype data previously run through BIRDSEED [12] and used PLINK [13] to remove low-quality SNPs and outlier samples. First, we remapped probe locations to hg38 using the remapping files provided by the GDC. We then removed genotype calls with BIRDSEED quality scores greater than 0.05. In addition, we removed SNPs with MAF < 0.01 and those with missing calls in more than 95% samples across the entire TCGA cohort from the analysis. We also filtered out individuals with more than 2.5% SNP calls missing. To account for technical artifacts including population stratification in the eQTL model, we performed principal component analysis (PCA) with PLINK to identify potential genetic covariates. By comparing TCGA genotypes with genotypes from Phase 3 of the 1000 Genomes project [14], we observed that the first five principal components were related to ancestry and accounted for 94.4% of the variance in genotypes (Supplementary Fig. 1). Thus, we included 5 principal components in our eQTL model.

Genotype phasing and imputation

We first phased genotypes with Eagle 2 [15], using data from Phase 3 of the 1000 Genomes project as a reference. All 1000 Genomes Project SNP coordinates were lifted over to reference genome hg38 prior to phasing. We performed phasing for all samples in the TCGA dataset together, and then separated samples by cancer type. We then imputed genotypes across the genome using Minimac 4 [16] and SNP calls from 1000 Genomes data with default settings. Finally, we filtered SNPs with MAF < 0.01 and R-squared less than 0.3, as recommended in the Minimac 4 guidelines. The final dataset consisted of approximately 10–11 million SNPs per sample.

RNA-seq data processing

We obtained raw HT-seq [17, 18] counts and RNA FPKMs for all samples from the GDC, which were generated with standardized alignment and read counting pipelines. We normalized data using the geometric mean method in DESeq2 [19] and further transformed the normalized counts using DESeq2’s variance stabilizing transformation for QC assessment and visualization. We performed PCA using the transformed data to identify sample outliers and removed individuals whose expression profiles lay more than three standard deviations away from the mean on any one of the six first principal components. In addition, we filtered lowly expressed genes, retaining only those with at least 0.1TPM and 6 reads in a minimum of 20% of the samples for each cancer type. These thresholds are similar to those used for GTEx analyses. To identify expression covariates for inclusion in the eQTL model, we used PEER [20], and included the first 15 PEER factors for each cancer type in the subsequent analyses.

eQTL mapping

We normalized expression data using a rank-based inverse normal transformation as described for GTEx [21]. We used a two-part regression for eQTL mapping, first regressing gene expression on genotype PCs, PEER factors, patient sex, copy number-changes and inactivating point mutations for each gene. We focused only on high-level amplifications and deep deletions and coded these + 1 and -1 respectively. We then normalized residuals from this regression using the rank-based inversed normal transformation, and used the resulting gene expression values as input to FastQTL [22]. We restricted analyses to variants for which the variant allele was observed at least five times in each cancer and then ran FastQTL with the parameters ‘—permute = 1000,10,000, –window = 1e6’. Mapping was performed using the imputed genotype dosages. The resulting p-values were corrected using Storey’s q-value method, and associations with FDR = 0.05 were retained. We used the same multiple testing correction strategy for somatic copy number alterations (CNAs) and inactivating mutations. We randomly sampled 200 patients from each cohort for downsampling analyses.

eQTL downstream analyses

We used the package variancePartition [23] to deconvolute gene expression variance into the relevant technical and genetic factors within each cancer type. In general, we deconvoluted variance for all covariates included in the additive eQTL model except for cases where data was missing or there were no somatic alterations in a gene. To evaluate the effects of individual variants, we used the allelic fold change method that has previously been described [24]. Functional enrichment of GO terms was analyzed using the R package goSeq [25], and the resulting p-values were corrected for using the Benjamini–Hochberg method. We defined a similarity index for two cancer types as the intersection divided by the union of the detected eGenes. This index was scaled between 0–1 for the purpose of visualization and interpretation. To assign gene type and gene boundaries, we used gene annotations from GENCODE v.22 (hg38), which was also used by the GDC for counting RNA-seq reads. We obtained ABSOLUTE-derived copy number instability estimates, genomic purity estimates and CIBERSORT [26] cellular fractions from previously published analyses [27]. We used the Consensus Purity Estimates (CPE) previously computed [28,29,30] as measurements of tumor purity, which considers IHC, gene expression, methylation and copy number data. In addition, we computed a version of the previously defined mutant allele tumor heterogeneity (MATH) statistic [31] to summarize intratumor heterogeneity for each sample. We calculated cancer cell fractions (CCF) [32] for all mutations in the dataset and defined the MATH score as the median absolute deviation in CCF divided by the median CCF for each patient.

Interaction and survival models

Following previously described strategies [33] we built cell type-restricted models by including an interaction term between genotype and tumor purity (expression = genotype + technical covariates + somatic alterations + purity + genotype*purity). We used a modified version of FastQTL to perform interaction analyses (https://github.com/francois-a/fastqt)l. and an FDR = 0.1 for identifying significant interactions thus accounting for the low power of detection in the interaction model. Where appropriate, we performed survival analyses using Cox proportional hazards models accounting for covariates as described in the text.

Results

The eQTL landscape in human cancers

We mapped proximal associations (within 1 Mb of gene boundaries as defined in GENCODE v. 22) between common polymorphisms (minor allele frequency; MAF > 1%) and gene expression in each of 24 cancer types from the TCGA project (estrogen receptor-positive (ER +) and ER- breast cancers were considered separately) (Table 1). In each cancer type, we employed an additive model that accounted for high-level somatic copy number events and inactivating point mutations (Methods). Overall, we identified 8,857 eGenes in at least one cancer type at a 5% genome-wide FDR, including 54 eGenes shared between all cancers and 4,346 eGenes present in only one cancer type (Supplementary Fig. 2A). The proportion of expressed genes (Supplementary Fig. 2B) that were also eGenes ranged between 0.9% (esophageal cancer) – 28.1% (thyroid cancer; median = 4.1%) with the highest proportion of eGenes/expressed genes in thyroid, prostate (22.0%), and ER + breast (15.0%) cancers (Fig. 1A). As the power to identify an eQTL is related to the number of patients in an indication (Spearman’s rho between sample size and eGene/expressed gene fraction = 0.83, p < 0.001; Supplementary Fig. 3), we repeated the analysis after downsampling each cancer type to 200 patients. In the 16 cancer datasets with at least 200 patients available, we identified eGenes in between 1.5% (ER- breast cancer) – 11.7% (thyroid cancer; median = 2.7%) of expressed genes using the downsampled data (Fig. 1B). Interestingly, previous eQTL analyses in normal tissue have also identified high number of eGenes in thyroid and prostate tissues [34, 35]. When using the downsampled data, there were 71 eGenes shared between all 16 cancer types and 2,184 eGenes unique to a single cancer type. Of these 2,184 unique eGenes, 1,572 (72.0%) were expressed in all 16 cancer types. This suggests that even genes that are expressed in multiple cancer types may be regulated by germline variants in only a subset of cancers.

Table 1 TCGA cohorts used in the study. Patient numbers include those who passed filtering criteria and had complete datasets (somatic alteration, gene expression, genotyping data) available. Only cancer types with at least 200 patients were used in downsampling analyses
Full size table
Fig. 1
figure 1

The eQTL landscape in 24 cancer types. A The proportions of expressed genes with evidence of germline regulation are shown for 24 cancer types from the TCGA. The number of cancer types in which eGenes are found are indicated by the colored bars. B eGenes identified after randomly selecting 200 patients from each cancer cohort. As before, bar colors represent tissue sharing between eGenes. Only cancers for which at least 200 patients were available are shown. C Similarities in eGene profiles between cancer types after downsampling. The similarity index was derived by considering the ratio of the intersection to the union of eGenes present in all pairs of cancers. On the adjusted scale, 0 = the lowest and 1 = the highest similarity observed between two different cancer types. D Overrepresentation analysis of GO enrichment terms when considering only cancer type-specific eGenes among all ubiquitously eGenes, E  Overrepresentation analysis of GO enrichemtn terms when considering shared eGenes (present in at least 10 cancer types) among all ubiquitously expressed eGenes. Enrichment analyses were performed considering all three GO categories (CC cellular compartment; BP biological process, MF molecular function) using goSeq. Adjusted p-values are shown for ten most significant GO terms

Full size image

eGene sharing between cancers

To identify similarities in eQTL profiles between cancers, we computed Jaccard similarity coefficients based on the overlap in eGenes between each pair of cancer types in the downsampled data (Fig. 1C; Methods). The coefficients ranged between 0.10 (thyroid and ER- breast cancer) – 0.43 (lung squamous cell carcinoma and lung adenomacarcinoma).

To compare relative similarities, we scaled the coefficients between 0–1, where 1 represented the highest similarity observed between two different cancers. As expected, given the relatively high proportion of unique eGenes in this cancer type, thyroid cancer was, on average, most dissimilar from the other cancer types. Nevertheless, thyroid cancer was most similar to prostate adenocarcinoma (scaled similarity coefficient = 0.73) and kidney renal cell cancer (scaled similarity coefficient = 0.58), supporting the notion that these three cancer types have substantially different eGene landscapes. Cancers with common histologies including lung adenocarcinoma and lung squamous cell cancer (scaled similarity coefficient = 0.94), and kidney renal cell cancer and kidney renal papilloma (scaled similarity coefficient = 0.95) also had similar eGene profiles. Interestingly, there was a lower overlap between ER + and ER- breast cancer (scaled similarity coefficient = 0.83), and ER- breast cancer was more similar in eGene composition to uterine cancer (0.98), bladder cancer (0.92) and cervical cancer (0.9). This may reflect similarities in somatic alteration profiles at the genomic level: all three cancer types, for example, are characterized by frequent mutation of TP53 and basal cell-like phenotypes. Together, these results demonstrate that although the overall eGene distribution resembles that found in normal tissues, differences in eGene profiles are sometimes apparent, potentially due to variation in tumor biology.

Next, to better understand the relevance of genes whose expression is under genetic control, we characterized the functional relevance of eGene sharing by performing Gene Ontology (GO) enrichment analysis. We first analyzed the 17,827 ubiquitously expressed genes to generate a background model. When comparing to this background set of genes, there was an overall enrichment for terms related to immune and MHC function, and metabolism within the 3,696 eGenes identified using downsampled data (Supplementary Fig. 4A). The eGene enrichment for immune function is consistent with the presence of eQTLs in infiltrating immune cells in heterogeneous tumor biopsies, and has been previously described in the cancer setting [36]. In the remaining 14,131 non-eGenes, there was significant overrepresentation of terms related to RNA biosynthesis and transcription (Supplementary Fig. 4B) indicating that the expression levels of genes involved in biosynthesis are less subject to germline regulation. As we were specifically interested in better understanding eGene sharing between cancers, we repeated the enrichment analysis this time using only the 3,696 eGenes in the background model. Comparison of unique eGenes (present in only one cancer type) with all identified eGenes included overrepresentation of functional terms related to developmental processes, although we did not find any statistically significant enrichment when considering these unique eGenes alongside the entire set of 17,287 genes in the background model (Fig. 1D). On the other hand, terms related to immune function were enriched for among shared eGenes (defined as eGenes detected in at least 10 cancer types) (Fig. 1E). Given the large number of thyroid cancer-specific eGenes detected and the intersection between developmental and cancer pathways, this result suggests that germline variation may be especially relevant for cancer biology in thyroid cancer and should therefore be considered alongside the spectrum of somatic alterations in this cancer type.

Contribution of eQTLs to the cancer transcriptome

Having characterized the global landscape of eQTLs in tumors, we next sought to understand the relative contribution of eQTLs to the cancer transcriptome and differences between cancer types. We first assessed the number of eQTLs and significant somatic variant-expression associations (somQTLs; high-level copy number amplifications, deep deletions, inactivating point mutations), as determined using the additive model, for all samples within each cancer type. The total number of variant-expression associations varied by the number of patients available for each cancer type, ranging between 527 (uveal melanoma) – 12,388 (ER + breast cancer; median = 4294) (Fig. 2A; Supplementary Fig. 5A). eQTL/somQTL ratios varied between 0.07 (ovarian cancer) – 62.3 (thyroid cancer; median = 0.47) (Supplementary Fig. 5B). After downsampling, the ratios of eQTL/somQTL were similar to the ratios observed in the complete dataset (Supplementary Fig. 5C). This reflects the fact that the relative power to detect eQTLs and somQTLs across cancer types is similar given that both eQTL and somQTL analyses are performed using the same datasets within each cancer type, and we subsequently used the complete datasets for our following analyses.

Fig. 2
figure 2

Contributions of eQTLs to the cancer transcriptome. (A) Numbers of total variant-expression associations identified across all cancer, where variants refer to inherited SNPs, somatic point mutations and somatic copy number changes. The proportions of these different classes of variants are also indicated. SNP  germline variants, SOMATIC  somatic variants, BOTH  germline and somatic variants. (B) Fraction of variance of individual genes attributable to different genetic factors. Gene experience variance was deconvoluted into all technical and biological factors included in the eQTL model (Methods) and the results for the genetic factors are shown here. Colors correspond to those in (A)

Full size image

We next used a linear mixed-effects model [5] to quantify the relative contributions of genetic factors to the variance in expression for each gene (Fig. 2B). Overall, median fractions of the proportion of variance attributable to either inherited variants or somatic alterations in a single gene were low. Nevertheless, somatic alterations accounted for the greatest proportion of variance in expression on a per-gene basis across all cancer types (mean proportions of variance for somQTLs = 0.02 (thyroid cancer) – 0.15 (ER + breast cancer); mean proportions of variance in expression accounted for by eQTLs were lower (0.01 (skin cutaneous melanoma) – 0.05 (uveal melanoma)). For example, in ER + breast cancer, somatic alterations accounted for at least 50% of variance in the expression of 4,158 genes but germline eQTLs accounted for at least 50% of variance in the expression of only 24 genes. Across all cancer types, germline variants accounted for at least 30% of expression variance in 522 genes, with 243/522 (46.6%) of these high-variance eGenes present in a single cancer type. In addition, 21 of these genes were eGenes almost universally across cancer types including five putative long non-coding RNAs, genes known to play a role in the immune response [24] (ICOSLG, ERAP2) and U2AF1, which is involved in spliceosome function. Our results indicate that somatic alterations dominate the cancer transcriptome, but that eQTLs can be important mediators of expression of some genes.

Finally, to understand the eQTL landscape in tumors relative to normal tissue, we compared eQTLs identified in 94 paired breast tumor and normal samples from the TCGA as breast was the indication with sufficient normal samples to facilitate a comparison. In the tumor samples, we identified 45 eGenes, of which 23/45 (51%) were unique to tumors. We identified 50 eGenes in the normal samples, including 28/50 (56.0%) eGenes detected only in normal samples.

Factors contributing to eGene detection

Given the large differences in the number of eGenes across indications, we examined additional genomic data to test for factors that may explain this difference (Supplementary Fig. 6). We first explored correlations between somatic alteration rates and eGene fractions in the downsampled datasets (Fig. 3A). The analysis revealed that indications with higher numbers of eGenes had lower somatic alteration rates than other cancer types in general. To further explore whether the effects of somatic alterations on gene expression was impacting our power to detect eQTLs, we focused on ‘quiet’ genes that are infrequently altered by copy number changes or point mutations (Fig. 3B). We observed that, even when removing genes with high somatic alteration rates from the analysis, cancer-specific patterns of eGene/expressed gene fractions remained similar to those from the dataset with all genes. This observation held true when defining quiet genes using different thresholds (genes with no somatic alterations, genes altered in no more than 5% of samples and genes altered in no more than 20% of samples). This result suggests that somatic alteration rates are not solely responsible for the differing numbers of detected eGenes across cancer types.

Fig. 3
figure 3

Detection of eQTLs in cancer tissue. A Correlation between somatic alteration rates and eGene fractions in the downsampled dataset. (Left) Correlation between tumor mutation burden and eGene fraction. (Right) Copy number alteration rate was measured by the fraction of genome altered by any type of copy number alteration. Spearman’s correlation statistics are shown. B Correlation between eGene fraction when considering all genes (x-axis) and when considering quiet genes only (those altered in subsets of samples) in the downsampled dataset. (Left) Quiet genes defined as those with no somatic alteration in any sample. (Middle) Quiet genes defined as those with somatic alterations in at most 5% of samples. (Right) Quiet genes defined as those with somatic alterations in at most 20% of samples. Only genes expressed in all samples were used in this analysis. Spearman’s correlation statistics are shown. C Correlations between tumor microenvironment scores as defined by ESTIMATE and eGene fraction. Standardized immune (left) and stromal (right) scores defined by ESTIMATE were obtained for all cancer types using the downsampled data, and the median score is shown for each cancer type. Spearman’s correlation statistics are shown. D Correlation between levels of intratumor heterogeneity and eGene fraction in downsampled data. The mutant allele tumor heterogeneity (MATH) score was computed for every sample (Methods) and the median MATH score was used to summarize intratumor heterogeneity for each cancer type. Higher MATH scores indicate higher levels of intratumor heterogeneity. Spearman’s correlation statistics are shown

Full size image

We next considered whether cellular heterogeneity within a biopsy influenced eQTL detection as the effect of an eQTL in a single cell type may be diluted in the presence of numerous other cell types. We first explored associations between stromal and immune infiltrate in a biopsy and eGene fraction in cancer using previously published estimates of immune and stromal cell tumor infiltration [5] (Fig. 3C; Supplementary Fig. 7). Although there were no significant correlations between eGene fraction and levels of either cell subpopulation, we observed that prostate adenocarcinoma and thyroid cancer had the lowest levels of immune infiltration across all cancer types. Expanding on this observation, we next reasoned that cellular diversity in cancer may also be mediated by varying levels of intratumor heterogeneity, which reflect the number of subclones present within a tumor biopsy. We adapted the previously published MATH score to measure global intratumor heterogeneity (Methods) [33] and observed that eGene fraction was positively correlated with the median MATH score for each cancer type (Fig. 3D). As higher MATH scores indicate higher levels of intratumor heterogeneity (i.e. the presence of more subclones within a tumor), this result suggests that more eGenes were detected in tumors with lower levels of intratumor heterogeneity. Thus, the strongest predictors of number of eGenes in a particular cancer type appear to be tumor intrinsic (subclonality) although we cannot rule out an additional role for tumor extrinsic factors in some cancer types.

Using an interaction model to identify cell type-restricted eQTLs

Given the association of cell-type heterogeneity with eQTL identification, we sought identify eQTLs specific to tumor cells or specific to the tumor microenvironment. We repeated our analyses, this time including an interaction term between genotype and tumor purity (Supplementary Fig. 6C) in the additive model. This updated model allowed us to detect eQTLs whose effects varied by tumor purity, allowing us to infer cell-type restricted eQTLs in heterogeneous tissue biopsies, although the exact cell types in which the eQTL is acting cannot be determined [31]. Using this approach in the entire cohorts for each cancer, we identified 2,271 interaction eGenes (ieGenes) at FDR = 10%, across all cancer types (Fig. 4A). These included 412 ieGenes in thyroid cancer and 360 ieGenes in prostate adenocarcinoma. A high number of ieGenes was also found in liver cancer (254) despite the relatively small number of eGenes identified in this cancer type using the base model. ieGenes are largely indication-specific with only 160 of the 2,271 (7%) ieGenes shared across indications and only 21 of 2,271 (0.9%) shared between more than 2 cancers. There was no association between the number of ieGenes observed and the median purity estimate obtained for a cancer type.

Fig. 4
figure 4

Detection of interaction eQTLs in heterogeneous biopsies. An interaction term was added to the eQTL model to identify ieQTLs, which are likely to be cell-type specific. A Total number of ieGenes identified using all patients for each cancer type. As for Fig. 1, the degree of ieGene sharing is also indicated. B GO enrichment terms related to the ieQTLs. The three GO categories were considered separately. C ieQTL associations in high and low-purity tumors as defined by the upper and lower tertiles of the Consensus Purity Estimate score (Methods). Germline variants are associated with CD1E expression in only low-purity but not in high-purity ER + breast tumors (upper panel), but are associated with AKR1C3 expression in high but not in low-purity ER + breast tumors (middle panel). In thyroid cancer, an ieQTL is associated with BCL7A expression in both high and low-purity tumors but in opposite directions (lower panel). Q-values indicate the adjusted p-values (Methods) associated with the interaction term in the eQTL model

Full size image

We next sought to identify the potential functional roles of ieGenes using GO enrichment analysis to study overrepresentation of functional terms (Fig. 4B). Terms related to immune function, stromal cells and cell adhesion were overrepresented over the full set of ieGenes, emphasizing that modelling expression based on tumor purity allows for identification of ieGenes in the non-tumor cells within a biopsy. Nevertheless, we did observe an enrichment of Molecular Function terms related to signaling and growth, mainly through modulation of the phosphatidylinositol 3-kinase pathway, which may be from tumor-restricted eQTL.

We hypothesized that we may be able to differentiate between tumor-specific and other ieQTLs based on the effects observed in when considering ‘high-purity’ and ‘low-purity’ tumors separately. We separated tumors based on whether their associated purity estimates were within the upper or low purity tertile for that cancer type. This allowed us to identify three classes of ieQTLs: those that were associated with gene expression in low-purity tumors but not in high-purity tumors, those that were associated with gene expression in high-purity tumors but not in low-purity tumors, and those were associated with gene expression in high- and low-purity tumors but in opposite directions (Fig. 4C). For example, an ieQTL was associated with CD1E expression in ER + breast cancer only in low-purity tumors. The CD1E gene encodes a protein involved in lipid antigen presentation and was expressed at higher levels in low-purity tumors suggesting that the gene is likely to be expressed in non-tumor cells. In contrast, an association between an ieQTL and AKR1C3 expression was detected in only high-purity tumors suggesting that this ieQTL may be tumor-restricted. Finally, germline variants in BCL7A had opposite effects in thyroid cancer when stratifying tumors by purity. In low-purity thyroid tumors, the variant ieQTL allele was associated with slightly lower BCL7A expression, but the same allele was associated with higher BCL7A expression in high-purity thyroid tumors. This observation suggests that the ieQTL function in BCL7A may differ between tumor and non-tumor cells within the context of thyroid cancer. These examples demonstrate that cell type-restricted eQTLs can be detected in tumor biopsies and that it may be possible to infer whether these act within or outside the tumor compartment.

Tumor-specific eGenes and patient outcome

Based on our observations that AKR1C3 and BCL7A have known functions in cancer and that ieQTLs in these genes are likely to be specific to tumor cells, we looked for other ieGenes across the datasets that may also play a role in cancer biology. Using the curated driver gene list from the COSMIC Cancer Gene Census [37]; genes are included on this list based on the presence of likely oncogenic somatic genetic variants. We identified 68 cancer-driver ieGenes (Fig. 5A). These included ABL1 and the master regulator NKX2-1 in thyroid cancer, ERBB3 in liver cancer and AKT3 in colorectal adenocarcinoma. Interestingly, expression levels of the ieGenes LPP in prostate adenocarcinoma and EZH2 in thyroid cancer have previously been associated with patient prognosis in cancer [38], and we observed a similar association within these two cancer types in the TCGA dataset (Supplementary Fig. 8).

Fig. 5
figure 5

Clinical relevance of interaction eQTLs. A ieGenes that were also identified in the COSMIC Cancer Gene Census. B Associations between genotypes in LPP and EZH2 and gene expression across all samples in PRAD and THCA. As in Fig. 4, q-values indicate the adjusted p-values associated with the interaction terms in the expanded model. C Kaplan–Meier curves showing the association between SNP genotypes and progression-free survival. Hazard ratios were computed when considering genotypes in an additive model

Full size image

Given the associations between patient genotype and gene expression for both LPP and EZH2 (Fig. 5B), we looked for associations between genotype and progression free- survival for these two genes. Our analysis revealed that the genotypes of the lead ieQTLs for both these genes were associated with patient outcomes (Fig. 5C). In thyroid cancer, the alternate allele for EZH2 was associated with lower expression of the gene and consequently longer progression-free survival (hazard ratio, HR, from univariate Cox proportional hazards model = 0.68; 95% confidence interval, CI = 0.47–1.00) consistent with the oncogenic role of EZH2 in cancer. Similarly, the alternate allele for LPP in prostate adenocarcinoma was associated with higher expression of the gene and better outcome (HR = 0.65, CI = 0.45–0.93). To test whether these observations were influenced by residual patient ancestry effects, we repeated the analysis focusing only on patients of European ancestry. In 301/448 (67.2%) thyroid cancer patients of European ancestry, the hazard ratio associated with the EZH2 eQTL modelled as an additive term was 0.64 (CI = 0.40 – 1.02), whereas the LPP eQTL was associated with a hazard ratio of 0.62 (CI = 0.45–0.93) in 356/443 (80.4%) prostate adenocarcinoma patients with European ancestry. These observations highlight the potential of germline variants acting through eQTLs to influence patient clinical trajectories.

Discussion

In this study, we have characterized the landscape of eQTLs in cancer, focusing on tumor-specific factors that may contribute to detection of eGenes, tissue specificity of eGenes and the possible roles ieGenes may play in mediating tumor biology and patient outcome.

We note that the overall landscape of eGenes across the cancer types reflects the landscape found in normal/non-cancerous tissues for GTEx, with the highest number of eGenes found in thyroid and prostate in GTEx, and in thyroid cancer and prostate adenocarcinoma in our analyses. Given the differences in sample numbers, tissue types and covariates used in the eQTL models, it is not possible to robustly compare eGene detection between the two datasets. While there were also insufficient samples in the TCGA dataset to perform a comprehensive analysis of eGenes in matched tumor/normal pairs, we observed differences in the detected eGenes between breast tumor and normal samples. Overall, we expect some differences in eQTL detection between tumor and normal samples to be due to the presence of somatic alterations in most cancers that dominate gene expression variance in the cancer transcriptome. Consequently, underlying tissue-specific regulatory mechanisms remain similar between normal and malignant tissue, but that there may be eQTLs that are specific to the cancer disease-state and future studies with larger numbers of paired tumor and normal samples would be needed to investigate the role of such tumor-specific eQTLs. The observation that gene expression variance attributable to somatic variants is generally higher than that attributable to germline genetic variants is in contrast to a previous study [5]. This discrepancy is most likely due to differences in coding copy number state in the additive models used in eQTL detection. We note that other strategies for modelling somatic copy number alterations may also influence eQTL detection and that the magnitude of change in expression is likely to be different between deep deletions and high-level amplifications. Factors contributing to eGene discovery both in GTEx and in TCGA include the cellular diversity present within a tissue sample, which most likely varies in an organ- or indication-specific manner [5]. In addition, the possible associations between eGene detection and cancer-specific features such as tumor mutation burden (TMB), chromosomal instability (CIN) and intratumor heterogeneity suggest additional factors that are relevant for eQTL detection in the malignant context. In particular, the association of higher levels of intratumor (sub-clonal) heterogeneity with lower numbers of eQTLs points to the complexity of the cancer transcriptome and the possibility of sub-clonal eQTLs, which would be difficult to identify in bulk gene expression data.

As has been described for normal tissues, genes without an eQTL (non-eGenes) are enriched for genes whose functions are related to cellular development. In tumors, non-eGenes were also enriched for those involved in transcription and RNA metabolism. We expect that the relative depletion of transcription-associated terms among cancer eGenes reflects increased transcriptional activity in cancer cells that are driven either directly or indirectly by somatic alterations and subsequently changes in pathway activity. These changes may disrupt the existing germline regulation of transcription in normal cells. Nevertheless, we observe significant overrepresentation of GO terms related to development when considering cancer type-specific eGenes relative to shared eGenes, which are enriched for terms related to immune and stromal function.

Previous studies have used eQTL analyses in the context of normal tissue to map the biological mechanisms by which SNPs identified through GWAS increase risk of cancer onset [39]. In contrast, the role of eQTLs in the context of the malignant transcriptome and the effects of germline SNPs on patient outcome post-diagnosis is less well understood. Here, we have shown that cell type-restricted eQTLs can be identified by modelling tumor purity as part of the standard eQTL analysis. We identified examples of germline SNPs near ieGenes that are associated with progression-free survival (LPP and EZH2), demonstrating the potential importance of germline SNPs acting through gene expression in cancer patient outcome post diagnosis. These results point to the importance of both considering all contributors to gene expression variance in cancer including germ line polymorphisms as well as the potential clinical importance of germline variation post cancer diagnosis.

Availability of data and materials

All data used in this study are available to the public. Somatic alteration data, gene expression data and clinical data from the TCGA used during the current study are available in the Genomics Data Commons repository (https://gdc.cancer.gov/). Germline data are hosted at dbGaP (https://dbgap.ncbi.nlm.nih.gov) and are available to the public following approval.

Abbreviations

Eqtl:

Expression quantitative trait loci

CIN:

Chromosomal instability

ER + /ER:

Estrogen-receptor positive/negative

FDR:

False discovery rate

GDC :

Genomics Data Commons

GO:

Gene Ontology

GTEx Project:

Genotype Expression Project

PCA:

Principal component analysis

SNP:

Single nucleotide polymorphism

somQTL:

Somatic quantitative trait loci

TCGA:

The Cancer Genome Atlas

TMB:

Tumor mutation burden

References

  1. L. Zhang, Gene Expression Profiles in Normal and Cancer Cells. Science (80-. ). 276, 1268–1272 (1997).

  2. Weinstein JN, Collisson EA, Mills GB, Shaw KRM, Ozenberger BA, Ellrott K, Sander C, Stuart JM, Chang K, Creighton CJ, Davis C, Donehower L, Drummond J, Wheeler D, Ally A, Balasundaram M, Birol I, Butterfield YSN, Chu A, Chuah E, Chun HJE, Dhalla N, Guin R, Hirst M, Hirst C, Holt RA, Jones SJM, Lee D, Li HI, Marra MA, Mayo M, Moore RA, Mungall AJ, Robertson AG, Schein JE, Sipahimalani P, Tam A, Thiessen N, Varhol RJ, Beroukhim R, Bhatt AS, Brooks AN, Cherniack AD, Freeman SS, Gabriel SB, Helman E, Jung J, Meyerson M, Ojesina AI, Pedamallu CS, Saksena G, Schumacher SE, Tabak B, Zack T, Lander ES, Bristow CA, Hadjipanayis A, Haseley P, Kucherlapati R, Lee S, Lee E, Luquette LJ, Mahadeshwar HS, Pantazi A, Parfenov M, Park PJ, Protopopov A, Ren X, Santoso N, Seidman J, Seth S, Song X, Tang J, Xi R, Xu AW, Yang L, Zeng D, Auman JT, Balu S, Buda E, Fan C, Hoadley KA, Jones CD, Meng S, Mieczkowski PA, Parker JS, Perou CM, Roach J, Shi Y, Silva GO, Tan D, Veluvolu U, Waring S, Wilkerson MD, Wu J, Zhao W, Bodenheimer T, Hayes DN, Hoyle AP, Jeffreys SR, Mose LE, Simons JV, Soloway MG, Baylin SB, Berman BP, Bootwalla MS, Danilova L, Herman JG, Hinoue T, Laird PW, Rhie SK, Shen H, Triche T, Weisenberger DJ, Carter SL, Cibulskis K, Chin L, Zhang J, Sougnez C, Wang M, Getz G, Dinh H, Doddapaneni HV, Gibbs R, Gunaratne P, Han Y, Kalra D, Kovar C, Lewis L, Morgan M, Morton D, Muzny D, Reid J, Xi L, Cho J, Dicara D, Frazer S, Gehlenborg N, Heiman DI, Kim J, Lawrence MS, Lin P, Liu Y, Noble MS, Stojanov P, Voet D, Zhang H, Zou L, Stewart C, Bernard B, Bressler R, Eakin A, Iype L, Knijnenburg T, Kramer R, Kreisberg R, Leinonen K, Lin J, Liu Y, Miller M, Reynolds SM, Rovira H, Shmulevich I, Thorsson V, Yang D, Zhang W, Amin S, Wu CJ, Wu CC, Akbani R, Aldape K, Baggerly KA, Broom B, Casasent TD, Cleland J, Dodda D, Edgerton M, Han L, Herbrich SM, Ju Z, Kim H, Lerner S, Li J, Liang H, Liu W, Lorenzi PL, Lu Y, Melott J, Nguyen L, Su X, Verhaak R, Wang W, Wong A, Yang Y, Yao J, Yao R, Yoshihara K, Yuan Y, Yung AK, Zhang N, Zheng S, Ryan M, Kane DW, Aksoy BA, Ciriello G, Dresdner G, Gao J, Gross B, Jacobsen A, Kahles A, Ladanyi M, Lee W, Van Lehmann K, Miller ML, Ramirez R, Rätsch G, Reva B, Schultz N, Senbabaoglu Y, Shen R, Sinha R, Sumer SO, Sun Y, Taylor BS, Weinhold N, Fei S, Spellman P, Benz C, Carlin D, Cline M, Craft B, Goldman M, Haussler D, Ma S, Ng S, Paull E, Radenbaugh A, Salama S, Sokolov A, Swatloski T, Uzunangelov V, Waltman P, Yau C, Zhu J, Hamilton SR, Abbott S, Abbott R, Dees ND, Delehaunty K, Ding L, Dooling DJ, Eldred JM, Fronick CC, Fulton R, Fulton LL, Kalicki-Veizer J, Kanchi KL, Kandoth C, Koboldt DC, Larson DE, Ley TJ, Lin L, Lu C, Magrini VJ, Mardis ER, McLellan MD, McMichael JF, Miller CA, O’Laughlin M, Pohl C, Schmidt H, Smith SM, Walker J, Wallis JW, Wendl MC, Wilson RK, Wylie T, Zhang Q, Burton R, Jensen MA, Kahn A, Pihl T, Pot D, Wan Y, Levine DA, Black AD, Bowen J, Frick J, Gastier-Foster JM, Harper HA, Helsel C, Leraas KM, Lichtenberg TM, McAllister C, Ramirez NC, Sharpe S, Wise L, Zmuda E, Chanock SJ, Davidsen T, Demchok JA, Eley G, Felau I, Sheth M, Sofia H, Staudt L, Tarnuzzer R, Wang Z, Yang L, Zhang J, Omberg L, Margolin A, Raphael BJ, Vandin F, Wu HT, Leiserson MDM, Benz SC, Vaske CJ, Noushmehr H, Wolf D, Veer LVT, Anastassiou D, Yang THO, Lopez-Bigas N, Gonzalez-Perez A, Tamborero D, Xia Z, Li W, Cho DY, Przytycka T, Hamilton M, McGuire S, Nelander S, Johansson P, Jörnsten R, Kling T. The cancer genome atlas pan-cancer analysis project. Nat Genet. 2013;45:1113–20.

    Article  Google Scholar 

  3. Calabrese C, Davidson NR, Demircioğlu D, Fonseca NA, He Y, Kahles A, Van Lehmann K, Liu F, Shiraishi Y, Soulette CM, Urban L, Calabrese C, Davidson NR, Demircioğlu D, Fonseca NA, He Y, Kahles A, Van Lehmann K, Liu F, Shiraishi Y, Soulette CM, Urban L, Greger L, Li S, Liu D, Perry MD, Xiang Q, Zhang F, Zhang J, Bailey P, Erkek S, Hoadley KA, Hou Y, Huska MR, Kilpinen H, Korbel JO, Marin MG, Markowski J, Nandi T, Pan-Hammarström Q, Pedamallu CS, Siebert R, Stark SG, Su H, Tan P, Waszak SM, Yung C, Zhu S, Awadalla P, Creighton CJ, Meyerson M, Ouellette BFF, Wu K, Yang H, Fonseca NA, Kahles A, Van Lehmann K, Urban L, Soulette CM, Shiraishi Y, Liu F, He Y, Demircioğlu D, Davidson NR, Calabrese C, Zhang J, Perry MD, Xiang Q, Greger L, Li S, Liu D, Stark SG, Zhang F, Amin SB, Bailey P, Chateigner A, Cortés-Ciriano I, Craft B, Erkek S, Frenkel-Morgenstern M, Goldman M, Hoadley KA, Hou Y, Huska MR, Khurana E, Kilpinen H, Korbel JO, Lamaze FC, Li C, Li X, Li X, Liu X, Marin MG, Markowski J, Nandi T, Nielsen MM, Ojesina AI, Pan-Hammarström Q, Park PJ, Pedamallu CS, Pedersen JS, Siebert R, Su H, Tan P, Teh BT, Wang J, Waszak SM, Xiong H, Yakneen S, Ye C, Yung C, Zhang X, Zheng L, Zhu J, Zhu S, Awadalla P, Creighton CJ, Meyerson M, Ouellette BFF, Wu K, Yang H, Göke J, Schwarz RF, Stegle O, Zhang Z, Brazma A, Rätsch G, Brooks AN, Brazma A, Brooks AN, Göke J, Rätsch G, Schwarz RF, Stegle O, Zhang Z. Genomic basis for RNA alterations in cancer. Nature. 2020;578:129–36.

    Article  CAS  Google Scholar 

  4. L. Fachal, H. Aschard, J. Beesley, D. R. Barnes, J. Allen, S. Kar, K. A. Pooley, J. Dennis, K. Michailidou, C. Turman, P. Soucy, A. Lemaçon, M. Lush, J. P. Tyrer, M. Ghoussaini, M. Moradi Marjaneh, X. Jiang, S. Agata, K. Aittomäki, M. R. Alonso, I. L. Andrulis, H. Anton-Culver, N. N. Antonenkova, A. Arason, V. Arndt, K. J. Aronson, B. K. Arun, B. Auber, P. L. Auer, J. Azzollini, J. Balmaña, R. B. Barkardottir, D. Barrowdale, A. Beeghly-Fadiel, J. Benitez, M. Bermisheva, K. Białkowska, A. M. Blanco, C. Blomqvist, W. Blot, N. V. Bogdanova, S. E. Bojesen, M. K. Bolla, B. Bonanni, A. Borg, K. Bosse, H. Brauch, H. Brenner, I. Briceno, I. W. Brock, A. Brooks-Wilson, T. Brüning, B. Burwinkel, S. S. Buys, Q. Cai, T. Caldés, M. A. Caligo, N. J. Camp, I. Campbell, F. Canzian, J. S. Carroll, B. D. Carter, J. E. Castelao, J. Chiquette, H. Christiansen, W. K. Chung, K. B. M. Claes, C. L. Clarke, J. M. Collée, S. Cornelissen, F. J. Couch, A. Cox, S. S. Cross, C. Cybulski, K. Czene, M. B. Daly, M. de la Hoya, P. Devilee, O. Diez, Y. C. Ding, G. S. Dite, S. M. Domchek, T. Dörk, I. Dos-Santos-Silva, A. Droit, S. Dubois, M. Dumont, M. Duran, L. Durcan, M. Dwek, D. M. Eccles, C. Engel, M. Eriksson, D. G. Evans, P. A. Fasching, O. Fletcher, G. Floris, H. Flyger, L. Foretova, W. D. Foulkes, E. Friedman, L. Fritschi, D. Frost, M. Gabrielson, M. Gago-Dominguez, G. Gambino, P. A. Ganz, S. M. Gapstur, J. Garber, J. A. García-Sáenz, M. M. Gaudet, V. Georgoulias, G. G. Giles, G. Glendon, A. K. Godwin, M. S. Goldberg, D. E. Goldgar, A. González-Neira, M. G. Tibiletti, M. H. Greene, M. Grip, J. Gronwald, A. Grundy, P. Guénel, E. Hahnen, C. A. Haiman, N. Håkansson, P. Hall, U. Hamann, P. A. Harrington, J. M. Hartikainen, M. Hartman, W. He, C. S. Healey, B. A. M. Heemskerk-Gerritsen, J. Heyworth, P. Hillemanns, F. B. L. Hogervorst, A. Hollestelle, M. J. Hooning, J. L. Hopper, A. Howell, G. Huang, P. J. Hulick, E. N. Imyanitov, C. Isaacs, M. Iwasaki, A. Jager, M. Jakimovska, A. Jakubowska, P. A. James, R. Janavicius, R. C. Jankowitz, E. M. John, N. Johnson, M. E. Jones, A. Jukkola-Vuorinen, A. Jung, R. Kaaks, D. Kang, P. M. Kapoor, B. Y. Karlan, R. Keeman, M. J. Kerin, E. Khusnutdinova, J. I. Kiiski, J. Kirk, C. M. Kitahara, Y.-D. Ko, I. Konstantopoulou, V.-M. Kosma, S. Koutros, K. Kubelka-Sabit, A. Kwong, K. Kyriacou, Y. Laitman, D. Lambrechts, E. Lee, G. Leslie, J. Lester, F. Lesueur, A. Lindblom, W.-Y. Lo, J. Long, A. Lophatananon, J. T. Loud, J. Lubiński, R. J. MacInnis, T. Maishman, E. Makalic, A. Mannermaa, M. Manoochehri, S. Manoukian, S. Margolin, M. E. Martinez, K. Matsuo, T. Maurer, D. Mavroudis, R. Mayes, L. McGuffog, C. McLean, N. Mebirouk, A. Meindl, A. Miller, N. Miller, M. Montagna, F. Moreno, K. Muir, A. M. Mulligan, V. M. Muñoz-Garzon, T. A. Muranen, S. A. Narod, R. Nassir, K. L. Nathanson, S. L. Neuhausen, H. Nevanlinna, P. Neven, F. C. Nielsen, L. Nikitina-Zake, A. Norman, K. Offit, E. Olah, O. I. Olopade, H. Olsson, N. Orr, A. Osorio, V. S. Pankratz, J. Papp, S. K. Park, T.-W. Park-Simon, M. T. Parsons, J. Paul, I. S. Pedersen, B. Peissel, B. Peshkin, P. Peterlongo, J. Peto, D. Plaseska-Karanfilska, K. Prajzendanc, R. Prentice, N. Presneau, D. Prokofyeva, M. A. Pujana, K. Pylkäs, P. Radice, S. J. Ramus, J. Rantala, R. Rau-Murthy, G. Rennert, H. A. Risch, M. Robson, A. Romero, M. Rossing, E. Saloustros, E. Sánchez-Herrero, D. P. Sandler, M. Santamariña, C. Saunders, E. J. Sawyer, M. T. Scheuner, D. F. Schmidt, R. K. Schmutzler, A. Schneeweiss, M. J. Schoemaker, B. Schöttker, P. Schürmann, C. Scott, R. J. Scott, L. Senter, C. M. Seynaeve, M. Shah, P. Sharma, C.-Y. Shen, X.-O. Shu, C. F. Singer, T. P. Slavin, S. Smichkoska, M. C. Southey, J. J. Spinelli, A. B. Spurdle, J. Stone, D. Stoppa-Lyonnet, C. Sutter, A. J. Swerdlow, R. M. Tamimi, Y. Y. Tan, W. J. Tapper, J. A. Taylor, M. R. Teixeira, M. Tengström, S. H. Teo, M. B. Terry, A. Teulé, M. Thomassen, D. L. Thull, M. Tischkowitz, A. E. Toland, R. A. E. M. Tollenaar, I. Tomlinson, D. Torres, G. Torres-Mejía, M. A. Troester, T. Truong, N. Tung, M. Tzardi, H.-U. Ulmer, C. M. Vachon, C. J. van Asperen, L. E. van der Kolk, E. J. van Rensburg, A. Vega, A. Viel, J. Vijai, M. J. Vogel, Q. Wang, B. Wappenschmidt, C. R. Weinberg, J. N. Weitzel, C. Wendt, H. Wildiers, R. Winqvist, A. Wolk, A. H. Wu, D. Yannoukakos, Y. Zhang, W. Zheng, D. Hunter, P. D. P. Pharoah, J. Chang-Claude, M. García-Closas, M. K. Schmidt, R. L. Milne, V. N. Kristensen, J. D. French, S. L. Edwards, A. C. Antoniou, G. Chenevix-Trench, J. Simard, D. F. Easton, P. Kraft, A. M. Dunning, Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes. Nat. Genet. 52, 56–73 (2020).

  5. Lim YW, Chen-Harris H, Mayba O, Lianoglou S, Wuster A, Bhangale T, Khan Z, Mariathasan S, Daemen A, Reeder J, Haverty PM, Forrest WF, Brauer M, Mellman I, Albert ML. Germline genetic polymorphisms influence tumor gene expression and immune cell infiltration. Proc Natl Acad Sci U S A. 2018;115:E11701–10.

    CAS  PubMed  PubMed Central  Google Scholar 

  6. Perola M, Wang DY, Makino S, Kettunen J, Meitinger T, Prokisch H, Veldink JH, Jansen RC, Milani L, Völker U, Schurmann C, Peters MJ, Yaghootkar H, Herder C, Grallert H, Esko T, Wijmenga C, Zolezzi F, Schramm K, Petersmann A, Fairfax BP, Strauch K, Fu J, van den Berg LH, Hernandez DG, Singleton AB, Melzer D, Homuth G, Melchiotti R, Roden M, Salomaa V, van Meurs JBJ, Kasela S, Ferrucci L, Andiappan AK, Ripatti S, Visschedijk M, Larbi A, Uitterlinden AG, Wood AR, Westra H-J, Franke L, Reinmaa E, Hofman A, Arends D, Peterson P, Poidinger M, Tserel L, Weersma RK, Rivadeneira F, Karjalainen J, Platteel M, Li Y, Rotzschke O, Knight JC, Lorbeer R, Lee B, Metspalu A. Cell Specific eQTL Analysis without Sorting Cells. PLOS Genet. 2015;11: e1005223.

    Article  Google Scholar 

  7. P. Geeleher, A. Nath, F. Wang, Z. Zhang, A. N. Barbeira, J. Fessler, R. L. Grossman, C. Seoighe, R. Stephanie Huang, Cancer expression quantitative trait loci (eQTLs) can be determined from heterogeneous tumor gene expression data by modeling variation in tumor purity. Genome Biol. 19, 130 (2018).

  8. S. Kim-Hellmuth, F. Aguet, M. Oliva, M. Muñoz-Aguirre, S. Kasela, V. Wucher, S. E. Castel, A. R. Hamel, A. Viñuela, A. L. Roberts, S. Mangul, X. Wen, G. Wang, A. N. Barbeira, D. Garrido-Martín, B. B. Nadel, Y. Zou, R. Bonazzola, J. Quan, A. Brown, A. Martinez-Perez, J. M. Soria, G. Getz, E. T. Dermitzakis, K. S. Small, M. Stephens, H. S. Xi, H. K. Im, R. Guigó, A. V. Segrè, B. E. Stranger, K. G. Ardlie, T. Lappalainen, Cell type–specific genetic regulation of gene expression across human tissues. Science (80-. ). 369, eaaz8528 (2020).

  9. Baxter JS, Leavy OC, Dryden NH, Maguire S, Johnson N, Fedele V, Simigdala N, Martin L, Andrews S, Wingett SW, Assiotis I, Fenwick K, Chauhan R, Rust AG, Orr N, Dudbridge F, Haider S, Fletcher O. Capture Hi-C identifies putative target genes at 33 breast cancer risk loci. Nat Commun. 2018;9:1028.

    Article  Google Scholar 

  10. Y. Qian, L. Zhang, M. Cai, H. Li, H. Xu, H. Yang, Z. Zhao, S. K. Rhie, P. J. Farnham, J. Shi, W. Lu, The prostate cancer risk variant rs55958994 regulates multiple gene expression through extreme long-range chromatin interaction to control tumor progression, 1–13 (2019).

  11. Grossman RL, Heath AP, Ferretti V, Varmus HE, Lowy DR, Kibbe WA, Staudt LM. Toward a Shared Vision for Cancer Genomic Data. N Engl J Med. 2016;375:1109–12.

    Article  Google Scholar 

  12. J. Liu, T. Lichtenberg, K. A. Hoadley, L. M. Poisson, A. J. Lazar, A. D. Cherniack, A. J. Kovatich, C. C. Benz, D. A. Levine, A. V. Lee, L. Omberg, D. M. Wolf, C. D. Shriver, V. Thorsson, S. J. Caesar-Johnson, J. A. Demchok, I. Felau, M. Kasapi, M. L. Ferguson, C. M. Hutter, H. J. Sofia, R. Tarnuzzer, Z. Wang, L. Yang, J. C. Zenklusen, J. (Julia) Zhang, S. Chudamani, J. Liu, L. Lolla, R. Naresh, T. Pihl, Q. Sun, Y. Wan, Y. Wu, J. Cho, T. DeFreitas, S. Frazer, N. Gehlenborg, G. Getz, D. I. Heiman, J. Kim, M. S. Lawrence, P. Lin, S. Meier, M. S. Noble, G. Saksena, D. Voet, H. Zhang, B. Bernard, N. Chambwe, V. Dhankani, T. Knijnenburg, R. Kramer, K. Leinonen, Y. Liu, M. Miller, S. Reynolds, I. Shmulevich, V. Thorsson, W. Zhang, R. Akbani, B. M. Broom, A. M. Hegde, Z. Ju, R. S. Kanchi, A. Korkut, J. Li, H. Liang, S. Ling, W. Liu, Y. Lu, G. B. Mills, K. S. Ng, A. Rao, M. Ryan, J. Wang, J. N. Weinstein, J. Zhang, A. Abeshouse, J. Armenia, D. Chakravarty, W. K. Chatila, I. de Bruijn, J. Gao, B. E. Gross, Z. J. Heins, R. Kundra, K. La, M. Ladanyi, A. Luna, M. G. Nissan, A. Ochoa, S. M. Phillips, E. Reznik, F. Sanchez-Vega, C. Sander, N. Schultz, R. Sheridan, S. O. Sumer, Y. Sun, B. S. Taylor, J. Wang, H. Zhang, P. Anur, M. Peto, P. Spellman, C. Benz, J. M. Stuart, C. K. Wong, C. Yau, D. N. Hayes, J. S. Parker, M. D. Wilkerson, A. Ally, M. Balasundaram, R. Bowlby, D. Brooks, R. Carlsen, E. Chuah, N. Dhalla, R. Holt, S. J. M. Jones, K. Kasaian, D. Lee, Y. Ma, M. A. Marra, M. Mayo, R. A. Moore, A. J. Mungall, K. Mungall, A. G. Robertson, S. Sadeghi, J. E. Schein, P. Sipahimalani, A. Tam, N. Thiessen, K. Tse, T. Wong, A. C. Berger, R. Beroukhim, A. D. Cherniack, C. Cibulskis, S. B. Gabriel, G. F. Gao, G. Ha, M. Meyerson, S. E. Schumacher, J. Shih, M. H. Kucherlapati, R. S. Kucherlapati, S. Baylin, L. Cope, L. Danilova, M. S. Bootwalla, P. H. Lai, D. T. Maglinte, D. J. Van Den Berg, D. J. Weisenberger, J. T. Auman, S. Balu, T. Bodenheimer, C. Fan, K. A. Hoadley, A. P. Hoyle, S. R. Jefferys, C. D. Jones, S. Meng, P. A. Mieczkowski, L. E. Mose, A. H. Perou, C. M. Perou, J. Roach, Y. Shi, J. V. Simons, T. Skelly, M. G. Soloway, D. Tan, U. Veluvolu, H. Fan, T. Hinoue, P. W. Laird, H. Shen, W. Zhou, M. Bellair, K. Chang, K. Covington, C. J. Creighton, H. Dinh, H. V. Doddapaneni, L. A. Donehower, J. Drummond, R. A. Gibbs, R. Glenn, W. Hale, Y. Han, J. Hu, V. Korchina, S. Lee, L. Lewis, W. Li, X. Liu, M. Morgan, D. Morton, D. Muzny, J. Santibanez, M. Sheth, E. Shinbro, L. Wang, M. Wang, D. A. Wheeler, L. Xi, F. Zhao, J. Hess, E. L. Appelbaum, M. Bailey, M. G. Cordes, L. Ding, C. C. Fronick, L. A. Fulton, R. S. Fulton, C. Kandoth, E. R. Mardis, M. D. McLellan, C. A. Miller, H. K. Schmidt, R. K. Wilson, D. Crain, E. Curley, J. Gardner, K. Lau, D. Mallery, S. Morris, J. Paulauskis, R. Penny, C. Shelton, T. Shelton, M. Sherman, E. Thompson, P. Yena, J. Bowen, J. M. Gastier-Foster, M. Gerken, K. M. Leraas, T. M. Lichtenberg, N. C. Ramirez, L. Wise, E. Zmuda, N. Corcoran, T. Costello, C. Hovens, A. L. Carvalho, A. C. de Carvalho, J. H. Fregnani, A. Longatto-Filho, R. M. Reis, C. Scapulatempo-Neto, H. C. S. Silveira, D. O. Vidal, A. Burnette, J. Eschbacher, B. Hermes, A. Noss, R. Singh, M. L. Anderson, P. D. Castro, M. Ittmann, D. Huntsman, B. Kohl, X. Le, R. Thorp, C. Andry, E. R. Duffy, V. Lyadov, O. Paklina, G. Setdikova, A. Shabunin, M. Tavobilov, C. McPherson, R. Warnick, R. Berkowitz, D. Cramer, C. Feltmate, N. Horowitz, A. Kibel, M. Muto, C. P. Raut, A. Malykh, J. S. Barnholtz-Sloan, W. Barrett, K. Devine, J. Fulop, Q. T. Ostrom, K. Shimmel, Y. Wolinsky, A. E. Sloan, A. De Rose, F. Giuliante, M. Goodman, B. Y. Karlan, C. H. Hagedorn, J. Eckman, J. Harr, J. Myers, K. Tucker, L. A. Zach, B. Deyarmin, H. Hu, L. Kvecher, C. Larson, R. J. Mural, S. Somiari, A. Vicha, T. Zelinka, J. Bennett, M. Iacocca, B. Rabeno, P. Swanson, M. Latour, L. Lacombe, B. Têtu, A. Bergeron, M. McGraw, S. M. Staugaitis, J. Chabot, H. Hibshoosh, A. Sepulveda, T. Su, T. Wang, O. Potapova, O. Voronina, L. Desjardins, O. Mariani, S. Roman-Roman, X. Sastre, M. H. Stern, F. Cheng, S. Signoretti, A. Berchuck, D. Bigner, E. Lipp, J. Marks, S. McCall, R. McLendon, A. Secord, A. Sharp, M. Behera, D. J. Brat, A. Chen, K. Delman, S. Force, F. Khuri, K. Magliocca, S. Maithel, J. J. Olson, T. Owonikoko, A. Pickens, S. Ramalingam, D. M. Shin, G. Sica, E. G. Van Meir, H. Zhang, W. Eijckenboom, A. Gillis, E. Korpershoek, L. Looijenga, W. Oosterhuis, H. Stoop, K. E. van Kessel, E. C. Zwarthoff, C. Calatozzolo, L. Cuppini, S. Cuzzubbo, F. DiMeco, G. Finocchiaro, L. Mattei, A. Perin, B. Pollo, C. Chen, J. Houck, P. Lohavanichbutr, A. Hartmann, C. Stoehr, R. Stoehr, H. Taubert, S. Wach, B. Wullich, W. Kycler, D. Murawa, M. Wiznerowicz, K. Chung, W. J. Edenfield, J. Martin, E. Baudin, G. Bubley, R. Bueno, A. De Rienzo, W. G. Richards, S. Kalkanis, T. Mikkelsen, H. Noushmehr, L. Scarpace, N. Girard, M. Aymerich, E. Campo, E. Giné, A. L. Guillermo, N. Van Bang, P. T. Hanh, B. D. Phu, Y. Tang, H. Colman, K. Evason, P. R. Dottino, J. A. Martignetti, H. Gabra, H. Juhl, T. Akeredolu, S. Stepa, D. Hoon, K. Ahn, K. J. Kang, F. Beuschlein, A. Breggia, M. Birrer, D. Bell, M. Borad, A. H. Bryce, E. Castle, V. Chandan, J. Cheville, J. A. Copland, M. Farnell, T. Flotte, N. Giama, T. Ho, M. Kendrick, J. P. Kocher, K. Kopp, C. Moser, D. Nagorney, D. O’Brien, B. P. O’Neill, T. Patel, G. Petersen, F. Que, M. Rivera, L. Roberts, R. Smallridge, T. Smyrk, M. Stanton, R. H. Thompson, M. Torbenson, J. D. Yang, L. Zhang, F. Brimo, J. A. Ajani, A. M. Angulo Gonzalez, C. Behrens, J. Bondaruk, R. Broaddus, B. Czerniak, B. Esmaeli, J. Fujimoto, J. Gershenwald, C. Guo, C. Logothetis, F. Meric-Bernstam, C. Moran, L. Ramondetta, D. Rice, A. Sood, P. Tamboli, T. Thompson, P. Troncoso, A. Tsao, I. Wistuba, C. Carter, L. Haydu, P. Hersey, V. Jakrot, H. Kakavand, R. Kefford, K. Lee, G. Long, G. Mann, M. Quinn, R. Saw, R. Scolyer, K. Shannon, A. Spillane, J. Stretch, M. Synott, J. Thompson, J. Wilmott, H. Al-Ahmadie, T. A. Chan, R. Ghossein, A. Gopalan, D. A. Levine, V. Reuter, S. Singer, B. Singh, N. V. Tien, T. Broudy, C. Mirsaidi, P. Nair, P. Drwiega, J. Miller, J. Smith, H. Zaren, J. W. Park, N. P. Hung, E. Kebebew, W. M. Linehan, A. R. Metwalli, K. Pacak, P. A. Pinto, M. Schiffman, L. S. Schmidt, C. D. Vocke, N. Wentzensen, R. Worrell, H. Yang, M. Moncrieff, C. Goparaju, J. Melamed, H. Pass, N. Botnariuc, I. Caraman, M. Cernat, I. Chemencedji, A. Clipca, S. Doruc, G. Gorincioi, S. Mura, M. Pirtac, I. Stancul, D. Tcaciuc, M. Albert, I. Alexopoulou, A. Arnaout, J. Bartlett, J. Engel, S. Gilbert, J. Parfitt, H. Sekhon, G. Thomas, D. M. Rassl, R. C. Rintoul, C. Bifulco, R. Tamakawa, W. Urba, N. Hayward, H. Timmers, A. Antenucci, F. Facciolo, G. Grazi, M. Marino, R. Merola, R. de Krijger, A. P. Gimenez-Roqueplo, A. Piché, S. Chevalier, G. McKercher, K. Birsoy, G. Barnett, C. Brewer, C. Farver, T. Naska, N. A. Pennell, D. Raymond, C. Schilero, K. Smolenski, F. Williams, C. Morrison, J. A. Borgia, M. J. Liptay, M. Pool, C. W. Seder, K. Junker, L. Omberg, M. Dinkin, G. Manikhas, D. Alvaro, M. C. Bragazzi, V. Cardinale, G. Carpino, E. Gaudio, D. Chesla, S. Cottingham, M. Dubina, F. Moiseenko, R. Dhanasekaran, K. F. Becker, K. P. Janssen, J. Slotta-Huspenina, M. H. Abdel-Rahman, D. Aziz, S. Bell, C. M. Cebulla, A. Davis, R. Duell, J. B. Elder, J. Hilty, B. Kumar, J. Lang, N. L. Lehman, R. Mandt, P. Nguyen, R. Pilarski, K. Rai, L. Schoenfield, K. Senecal, P. Wakely, P. Hansen, R. Lechan, J. Powers, A. Tischler, W. E. Grizzle, K. C. Sexton, A. Kastl, J. Henderson, S. Porten, J. Waldmann, M. Fassnacht, S. L. Asa, D. Schadendorf, M. Couce, M. Graefen, H. Huland, G. Sauter, T. Schlomm, R. Simon, P. Tennstedt, O. Olabode, M. Nelson, O. Bathe, P. R. Carroll, J. M. Chan, P. Disaia, P. Glenn, R. K. Kelley, C. N. Landen, J. Phillips, M. Prados, J. Simko, K. Smith-McCune, S. VandenBerg, K. Roggin, A. Fehrenbach, A. Kendler, S. Sifri, R. Steele, A. Jimeno, F. Carey, I. Forgie, M. Mannelli, M. Carney, B. Hernandez, B. Campos, C. Herold-Mende, C. Jungk, A. Unterberg, A. von Deimling, A. Bossler, J. Galbraith, L. Jacobus, M. Knudson, T. Knutson, D. Ma, M. Milhem, R. Sigmund, A. K. Godwin, R. Madan, H. G. Rosenthal, C. Adebamowo, S. N. Adebamowo, A. Boussioutas, D. Beer, T. Giordano, A. M. Mes-Masson, F. Saad, T. Bocklage, L. Landrum, R. Mannel, K. Moore, K. Moxley, R. Postier, J. Walker, R. Zuna, M. Feldman, F. Valdivieso, R. Dhir, J. Luketich, E. M. Mora Pinero, M. Quintero-Aguilo, C. G. Carlotti, J. S. Dos Santos, R. Kemp, A. Sankarankuty, D. Tirapelli, J. Catto, K. Agnew, E. Swisher, J. Creaney, B. Robinson, C. S. Shelley, E. M. Godwin, S. Kendall, C. Shipman, C. Bradford, T. Carey, A. Haddad, J. Moyer, L. Peterson, M. Prince, L. Rozek, G. Wolf, R. Bowman, K. M. Fong, I. Yang, R. Korst, W. K. Rathmell, J. L. Fantacone-Campbell, J. A. Hooke, A. J. Kovatich, C. D. Shriver, J. DiPersio, B. Drake, R. Govindan, S. Heath, T. Ley, B. Van Tine, P. Westervelt, M. A. Rubin, J. Il Lee, N. D. Aredes, A. Mariamidze, H. Hu, An Integrated TCGA Pan-Cancer Clinical Data Resource to Drive High-Quality Survival Outcome Analytics. Cell. 173, 400–416.e11 (2018).

  13. Korn JM, Kuruvilla FG, McCarroll SA, Wysoker A, Nemesh J, Cawley S, Hubbell E, Veitch J, Collins PJ, Darvishi K, Lee C, Nizzari MM, Gabriel SB, Purcell S, Daly MJ, Altshuler D. Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat Genet. 2008;40:1253–60.

    Article  CAS  Google Scholar 

  14. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC. PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage Analyses. Am J Hum Genet. 2007;81:559–75.

    Article  CAS  Google Scholar 

  15. The 1000 Genomes Consortium, A global reference for human genetic variation. Nature. 526, 68–74 (2015).

  16. Loh P, Danecek P, Palamara PF, Fuchsberger C, Reshef YA, Finucane HK, Schoenherr S, Forer L, McCarthy S, Abecasis GR, Durbin R, Price AL. Reference-based phasing using the Haplotype Reference Consortium panel. Nat Genet. 2016;48:1443–8.

    Article  CAS  Google Scholar 

  17. Howie B, Fuchsberger C, Stephens M, Marchini J, Abecasis GR. Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat Genet. 2012;44:955–9.

    Article  CAS  Google Scholar 

  18. S. Das, L. Forer, S. Schönherr, C. Sidore, A. E. Locke, A. Kwong, S. I. Vrieze, E. Y. Chew, S. Levy, M. Mcgue, D. Schlessinger, D. Stambolian, P. Loh, W. G. Iacono, A. Swaroop, L. J. Scott, F. Cucca, F. Kronenberg, M. Boehnke, G. R. Abecasis, C. Fuchsberger, Next-generation genotype imputation service and methods. 48 (2016), doi:https://doi.org/10.1038/ng.3656.

  19. Anders S, Pyl PT, Huber W. HTSeq–a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.

    Article  CAS  Google Scholar 

  20. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:1–21.

    Article  Google Scholar 

  21. Stegle O, Parts L, Piipari M, Winn J, Durbin R. Using probabilistic estimation of expression residuals (PEER) to obtain increased power and interpretability of gene expression analyses. Nat Protoc. 2012;7:500–7.

    Article  CAS  Google Scholar 

  22. Aguet F, Brown AA, Castel SE, Davis JR, He Y, Jo B, Mohammadi P, Park YS, Parsana P, Segrè AV, Strober BJ, Zappala Z, Cummings BB, Gelfand ET, Hadley K, Huang KH, Lek M, Li X, Nedzel JL, Nguyen DY, Noble MS, Sullivan TJ, Tukiainen T, MacArthur DG, Getz G, Addington A, Guan P, Koester S, Little AR, Lockhart NC, Moore HM, Rao A, Struewing JP, Volpi S, Brigham LE, Hasz R, Hunter M, Johns C, Johnson M, Kopen G, Leinweber WF, Lonsdale JT, McDonald A, Mestichelli B, Myer K, Roe B, Salvatore M, Shad S, Thomas JA, Walters G, Washington M, Wheeler J, Bridge J, Foster BA, Gillard BM, Karasik E, Kumar R, Miklos M, Moser MT, Jewell SD, Montroy RG, Rohrer DC, Valley D, Mash DC, Davis DA, Sobin L, Barcus ME, Branton PA, Abell NS, Balliu B, Delaneau O, Frésard L, Gamazon ER, Garrido-Martín D, Gewirtz ADH, Gliner G, Gloudemans MJ, Han B, He AZ, Hormozdiari F, Li X, Liu B, Kang EY, McDowell IC, Ongen H, Palowitch JJ, Peterson CB, Quon G, Ripke S, Saha A, Shabalin AA, Shimko TC, Sul JH, Teran NA, Tsang EK, Zhang H, Zhou YH, Bustamante CD, Cox NJ, Guigó R, Kellis M, McCarthy MI, Conrad DF, Eskin E, Li G, Nobel AB, Sabatti C, Stranger BE, Wen X, Wright FA, Ardlie KG, Dermitzakis ET, Lappalainen T, Battle A, Brown CD, Engelhardt BE, Montgomery SB, Handsaker RE, Kashin S, Karczewski KJ, Nguyen DT, Trowbridge CA, Barshir R, Basha O, Bogu GK, Chen LS, Chiang C, Damani FN, Ferreira PG, Hall IM, Howald C, Im HK, Kim Y, Kim-Hellmuth S, Mangul S, Monlong J, Muñoz-Aguirre M, Ndungu AW, Nicolae DL, Oliva M, Panousis N, Papasaikas P, Payne AJ, Quan J, Reverter F, Sammeth M, Scott AJ, Sodaei R, Stephens M, Urbut S, Van De Bunt M, Wang G, Xi HS, Yeger-Lotem E, Zaugg JB, Akey JM, Bates D, Chan J, Claussnitzer M, Demanelis K, Diegel M, Doherty JA, Feinberg AP, Fernando MS, Halow J, Hansen KD, Haugen E, Hickey PF, Hou L, Jasmine F, Jian R, Jiang L, Johnson A, Kaul R, Kibriya MG, Lee K, Li JB, Li Q, Lin J, Lin S, Linder S, Linke C, Liu Y, Maurano MT, Molinie B, Nelson J, Neri FJ, Park Y, Pierce BL, Rinaldi NJ, Rizzardi LF, Sandstrom R, Skol A, Smith KS, Snyder MP, Stamatoyannopoulos J, Tang H, Wang L, Wang M, Van Wittenberghe N, Wu F, Zhang R, Nierras CR, Carithers LJ, Vaught JB, Gould SE, Lockart NC, Martin C, Addington AM, Koester SE, Undale AH, Smith AM, Tabor DE, Roche NV, McLean JA, Vatanian N, Robinson KL, Valentino KM, Qi L, Hunter S, Hariharan P, Singh S, Um KS, Matose T, Tomaszewski MM, Barker LK, Mosavel M, Siminoff LA, Traino HM, Flicek P, Juettemann T, Ruffier M, Sheppard D, Taylor K, Trevanion SJ, Zerbino DR, Craft B, Goldman M, Haeussler M, Kent WJ, Lee CM, Paten B, Rosenbloom KR, Vivian J, Zhu J. Genetic effects on gene expression across human tissues. Nature. 2017;550:204–13.

    Article  Google Scholar 

  23. Ongen H, Buil A, Brown AA, Dermitzakis ET, Delaneau O. Fast and efficient QTL mapper for thousands of molecular phenotypes. Bioinformatics. 2016;32:1479–85.

    Article  CAS  Google Scholar 

  24. Hoffman GE, Schadt EE. variancePartition: interpreting drivers of variation in complex gene expression studies. BMC Bioinformatics. 2016;17:483.

    Article  Google Scholar 

  25. Mohammadi P, Castel SE, Brown AA, Lappalainen T. Quantifying the regulatory effect size of cis-acting genetic variation using allelic fold change. Genome Res. 2017;27:1872–84.

    Article  CAS  Google Scholar 

  26. Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.

    Article  Google Scholar 

  27. Newman AM, Steen CB, Liu CL, Gentles AJ, Chaudhuri AA, Scherer F, Khodadoust MS, Esfahani MS, Luca BA, Steiner D, Diehn M, Alizadeh AA. Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol. 2019;37:773–82.

    Article  CAS  Google Scholar 

  28. V. V. Thorsson, D. L. Gibbs, S. D. Brown, D. Wolf, D. S. Bortone, T. H. Ou Yang, E. Porta-Pardo, G. F. Gao, C. L. Plaisier, J. A. Eddy, E. Ziv, A. C. Culhane, E. O. Paull, I. K. A. Sivakumar, A. J. Gentles, R. Malhotra, F. Farshidfar, A. Colaprico, J. S. Parker, L. E. Mose, N. S. Vo, J. Liu, Y. Liu, J. Rader, V. Dhankani, S. M. Reynolds, R. Bowlby, A. Califano, A. D. Cherniack, D. Anastassiou, D. Bedognetti, A. Rao, K. Chen, A. Krasnitz, H. Hu, T. M. Malta, H. Noushmehr, C. S. Pedamallu, S. Bullman, A. I. Ojesina, A. Lamb, W. Zhou, H. Shen, T. K. Choueiri, J. N. Weinstein, J. Guinney, J. Saltz, R. A. Holt, C. E. Rabkin, S. J. Caesar-Johnson, J. A. Demchok, I. Felau, M. Kasapi, M. L. Ferguson, C. M. Hutter, H. J. Sofia, R. Tarnuzzer, Z. Wang, L. Yang, J. C. Zenklusen, J. J. (Julia) Zhang, S. Chudamani, J. Liu, L. Lolla, R. Naresh, T. Pihl, Q. Sun, Y. Wan, Y. Wu, J. Cho, T. DeFreitas, S. Frazer, N. Gehlenborg, G. Getz, D. I. Heiman, J. Kim, M. S. Lawrence, P. Lin, S. Meier, M. S. Noble, G. Saksena, D. Voet, H. H. H. Zhang, B. Bernard, N. Chambwe, V. Dhankani, T. Knijnenburg, R. Kramer, K. Leinonen, Y. Liu, M. Miller, S. M. Reynolds, I. Shmulevich, V. V. Thorsson, W. Zhang, R. Akbani, B. M. Broom, A. M. Hegde, Z. Ju, R. S. Kanchi, A. Korkut, J. Li, H. Liang, S. Ling, W. Liu, Y. Lu, G. B. Mills, K. S. Ng, A. Rao, M. Ryan, J. J. Wang, J. N. Weinstein, J. J. (Julia) Zhang, A. Abeshouse, J. Armenia, D. Chakravarty, W. K. Chatila, I. de Bruijn, J. Gao, B. E. Gross, Z. J. Heins, R. Kundra, K. La, M. Ladanyi, A. Luna, M. G. Nissan, A. Ochoa, S. M. Phillips, E. Reznik, F. Sanchez-Vega, C. Sander, N. Schultz, R. Sheridan, S. O. Sumer, Y. Sun, B. S. Taylor, J. J. Wang, H. H. H. Zhang, P. Anur, M. Peto, P. Spellman, C. Benz, J. M. Stuart, C. K. Wong, C. Yau, D. N. Hayes, J. S. Parker, M. D. Wilkerson, A. Ally, M. Balasundaram, R. Bowlby, D. Brooks, R. Carlsen, E. Chuah, N. Dhalla, R. A. Holt, S. J. M. Jones, K. Kasaian, D. Lee, Y. Ma, M. A. Marra, M. Mayo, R. A. Moore, A. J. Mungall, K. Mungall, A. G. Robertson, S. Sadeghi, J. E. Schein, P. Sipahimalani, A. Tam, N. Thiessen, K. Tse, T. Wong, A. C. Berger, R. Beroukhim, A. D. Cherniack, C. Cibulskis, S. B. Gabriel, G. F. Gao, G. Ha, M. Meyerson, S. E. Schumacher, J. Shih, M. H. Kucherlapati, R. S. Kucherlapati, S. Baylin, L. Cope, L. Danilova, M. S. Bootwalla, P. H. Lai, D. T. Maglinte, D. J. Van Den Berg, D. J. Weisenberger, J. T. Auman, S. Balu, T. Bodenheimer, C. Fan, K. A. Hoadley, A. P. Hoyle, S. R. Jefferys, C. D. Jones, S. Meng, P. A. Mieczkowski, L. E. Mose, A. H. Perou, C. M. Perou, J. Roach, Y. Shi, J. V. Simons, T. Skelly, M. G. Soloway, D. Tan, U. Veluvolu, H. Fan, T. Hinoue, P. W. Laird, H. Shen, W. Zhou, M. Bellair, K. Chang, K. Covington, C. J. Creighton, H. Dinh, H. V. Doddapaneni, L. A. Donehower, J. Drummond, R. A. Gibbs, R. Glenn, W. Hale, Y. Han, J. Hu, V. Korchina, S. Lee, L. Lewis, W. Li, X. Liu, M. Morgan, D. Morton, D. Muzny, J. Santibanez, M. Sheth, E. Shinbrot, L. Wang, M. Wang, D. A. Wheeler, L. Xi, F. Zhao, J. Hess, E. L. Appelbaum, M. Bailey, M. G. Cordes, L. Ding, C. C. Fronick, L. A. Fulton, R. S. Fulton, C. Kandoth, E. R. Mardis, M. D. McLellan, C. A. Miller, H. K. Schmidt, R. K. Wilson, D. Crain, E. Curley, J. Gardner, K. Lau, D. Mallery, S. Morris, J. Paulauskis, R. Penny, C. Shelton, T. Shelton, M. Sherman, E. Thompson, P. Yena, J. Bowen, J. M. Gastier-Foster, M. Gerken, K. M. Leraas, T. M. Lichtenberg, N. C. Ramirez, L. Wise, E. Zmuda, N. Corcoran, T. Costello, C. Hovens, A. L. Carvalho, A. C. de Carvalho, J. H. Fregnani, A. Longatto-Filho, R. M. Reis, C. Scapulatempo-Neto, H. C. S. Silveira, D. O. Vidal, A. Burnette, J. Eschbacher, B. Hermes, A. Noss, R. Singh, M. L. Anderson, P. D. Castro, M. Ittmann, D. Huntsman, B. Kohl, X. Le, R. Thorp, C. Andry, E. R. Duffy, V. Lyadov, O. Paklina, G. Setdikova, A. Shabunin, M. Tavobilov, C. McPherson, R. Warnick, R. Berkowitz, D. Cramer, C. Feltmate, N. Horowitz, A. Kibel, M. Muto, C. P. Raut, A. Malykh, J. S. Barnholtz-Sloan, W. Barrett, K. Devine, J. Fulop, Q. T. Ostrom, K. Shimmel, Y. Wolinsky, A. E. Sloan, A. De Rose, F. Giuliante, M. Goodman, B. Y. Karlan, C. H. Hagedorn, J. Eckman, J. Harr, J. Myers, K. Tucker, L. A. Zach, B. Deyarmin, H. Hu, L. Kvecher, C. Larson, R. J. Mural, S. Somiari, A. Vicha, T. Zelinka, J. Bennett, M. Iacocca, B. Rabeno, P. Swanson, M. Latour, L. Lacombe, B. Têtu, A. Bergeron, M. McGraw, S. M. Staugaitis, J. Chabot, H. Hibshoosh, A. Sepulveda, T. Su, T. Wang, O. Potapova, O. Voronina, L. Desjardins, O. Mariani, S. Roman-Roman, X. Sastre, M. H. Stern, F. Cheng, S. Signoretti, A. Berchuck, D. Bigner, E. Lipp, J. Marks, S. McCall, R. McLendon, A. Secord, A. Sharp, M. Behera, D. J. Brat, A. Chen, K. Delman, S. Force, F. Khuri, K. Magliocca, S. Maithel, J. J. Olson, T. Owonikoko, A. Pickens, S. Ramalingam, D. M. Shin, G. Sica, E. G. Van Meir, H. H. H. Zhang, W. Eijckenboom, A. Gillis, E. Korpershoek, L. Looijenga, W. Oosterhuis, H. Stoop, K. E. van Kessel, E. C. Zwarthoff, C. Calatozzolo, L. Cuppini, S. Cuzzubbo, F. DiMeco, G. Finocchiaro, L. Mattei, A. Perin, B. Pollo, C. Chen, J. Houck, P. Lohavanichbutr, A. Hartmann, C. Stoehr, R. Stoehr, H. Taubert, S. Wach, B. Wullich, W. Kycler, D. Murawa, M. Wiznerowicz, K. Chung, W. J. Edenfield, J. Martin, E. Baudin, G. Bubley, R. Bueno, A. De Rienzo, W. G. Richards, S. Kalkanis, T. Mikkelsen, H. Noushmehr, L. Scarpace, N. Girard, M. Aymerich, E. Campo, E. Giné, A. L. Guillermo, N. Van Bang, P. T. Hanh, B. D. Phu, Y. Tang, H. Colman, K. Evason, P. R. Dottino, J. A. Martignetti, H. Gabra, H. Juhl, T. Akeredolu, S. Stepa, D. Hoon, K. Ahn, K. J. Kang, F. Beuschlein, A. Breggia, M. Birrer, D. Bell, M. Borad, A. H. Bryce, E. Castle, V. Chandan, J. Cheville, J. A. Copland, M. Farnell, T. Flotte, N. Giama, T. Ho, M. Kendrick, J. P. Kocher, K. Kopp, C. Moser, D. Nagorney, D. O’Brien, B. P. O’Neill, T. Patel, G. Petersen, F. Que, M. Rivera, L. Roberts, R. Smallridge, T. Smyrk, M. Stanton, R. H. Thompson, M. Torbenson, J. D. Yang, L. Zhang, F. Brimo, J. A. Ajani, A. M. A. Gonzalez, C. Behrens, J. Bondaruk, R. Broaddus, B. Czerniak, B. Esmaeli, J. Fujimoto, J. Gershenwald, C. Guo, A. J. Lazar, C. Logothetis, F. Meric-Bernstam, C. Moran, L. Ramondetta, D. Rice, A. Sood, P. Tamboli, T. Thompson, P. Troncoso, A. Tsao, I. Wistuba, C. Carter, L. Haydu, P. Hersey, V. Jakrot, H. Kakavand, R. Kefford, K. Lee, G. Long, G. Mann, M. Quinn, R. Saw, R. Scolyer, K. Shannon, A. Spillane, onathan Stretch, M. Synott, J. Thompson, J. Wilmott, H. Al-Ahmadie, T. A. Chan, R. Ghossein, A. Gopalan, D. A. Levine, V. Reuter, S. Singer, B. Singh, N. V. Tien, T. Broudy, C. Mirsaidi, P. Nair, P. Drwiega, J. Miller, J. Smith, H. Zaren, J. W. Park, N. P. Hung, E. Kebebew, W. M. Linehan, A. R. Metwalli, K. Pacak, P. A. Pinto, M. Schiffman, L. S. Schmidt, C. D. Vocke, N. Wentzensen, R. Worrell, H. Yang, M. Moncrieff, C. Goparaju, J. Melamed, H. Pass, N. Botnariuc, I. Caraman, M. Cernat, I. Chemencedji, A. Clipca, S. Doruc, G. Gorincioi, S. Mura, M. Pirtac, I. Stancul, D. Tcaciuc, M. Albert, I. Alexopoulou, A. Arnaout, J. Bartlett, J. Engel, S. Gilbert, J. Parfitt, H. Sekhon, G. Thomas, D. M. Rassl, R. C. Rintoul, C. Bifulco, R. Tamakawa, W. Urba, N. Hayward, H. Timmers, A. Antenucci, F. Facciolo, G. Grazi, M. Marino, R. Merola, R. de Krijger, A. P. Gimenez-Roqueplo, A. Piché, S. Chevalier, G. McKercher, K. Birsoy, G. Barnett, C. Brewer, C. Farver, T. Naska, N. A. Pennell, D. Raymond, C. Schilero, K. Smolenski, F. Williams, C. Morrison, J. A. Borgia, M. J. Liptay, M. Pool, C. W. Seder, K. Junker, L. Omberg, M. Dinkin, G. Manikhas, D. Alvaro, M. C. Bragazzi, V. Cardinale, G. Carpino, E. Gaudio, D. Chesla, S. Cottingham, M. Dubina, F. Moiseenko, R. Dhanasekaran, K. F. Becker, K. P. Janssen, J. Slotta-Huspenina, M. H. Abdel-Rahman, D. Aziz, S. Bell, C. M. Cebulla, A. Davis, R. Duell, J. B. Elder, J. Hilty, B. Kumar, J. Lang, N. L. Lehman, R. Mandt, P. Nguyen, R. Pilarski, K. Rai, L. Schoenfield, K. Senecal, P. Wakely, P. Hansen, R. Lechan, J. Powers, A. Tischler, W. E. Grizzle, K. C. Sexton, A. Kastl, J. Henderson, S. Porten, J. Waldmann, M. Fassnacht, S. L. Asa, D. Schadendorf, M. Couce, M. Graefen, H. Huland, G. Sauter, T. Schlomm, R. Simon, P. Tennstedt, O. Olabode, M. Nelson, O. Bathe, P. R. Carroll, J. M. Chan, P. Disaia, P. Glenn, R. K. Kelley, C. N. Landen, J. Phillips, M. Prados, J. Simko, K. Smith-McCune, S. VandenBerg, K. Roggin, A. Fehrenbach, A. Kendler, S. Sifri, R. Steele, A. Jimeno, F. Carey, I. Forgie, M. Mannelli, M. Carney, B. Hernandez, B. Campos, C. Herold-Mende, C. Jungk, A. Unterberg, A. von Deimling, A. Bossler, J. Galbraith, L. Jacobus, M. Knudson, T. Knutson, D. Ma, M. Milhem, R. Sigmund, A. K. Godwin, R. Madan, H. G. Rosenthal, C. Adebamowo, S. N. Adebamowo, A. Boussioutas, D. Beer, T. Giordano, A. M. Mes-Masson, F. Saad, T. Bocklage, L. Landrum, R. Mannel, K. Moore, K. Moxley, R. Postier, J. Walker, R. Zuna, M. Feldman, F. Valdivieso, R. Dhir, J. Luketich, E. M. M. Pinero, M. Quintero-Aguilo, C. G. Carlotti, J. S. Dos Santos, R. Kemp, A. Sankarankuty, D. Tirapelli, J. Catto, K. Agnew, E. Swisher, J. Creaney, B. Robinson, C. S. Shelley, E. M. Godwin, S. Kendall, C. Shipman, C. Bradford, T. Carey, A. Haddad, J. Moyer, L. Peterson, M. Prince, L. Rozek, G. Wolf, R. Bowman, K. M. Fong, I. Yang, R. Korst, W. K. Rathmell, J. L. Fantacone-Campbell, J. A. Hooke, A. J. Kovatich, C. D. Shriver, J. DiPersio, B. Drake, R. Govindan, S. Heath, T. Ley, B. Van Tine, P. Westervelt, M. A. Rubin, J. Il Lee, N. D. Aredes, A. Mariamidze, A. J. Lazar, J. S. Serody, E. G. Demicco, M. L. Disis, B. G. Vincent, llya Shmulevich, The Immune Landscape of Cancer. Immunity. 48, 812–830.e14 (2018).

  29. Carter SL, Cibulskis K, Helman E, McKenna A, Shen H, Zack T, Laird PW, Onofrio RC, Winckler W, Weir BA, Beroukhim R, Pellman D, Levine DA, Lander ES, Meyerson M, Getz G. Absolute quantification of somatic DNA alterations in human cancer. Nat Biotechnol. 2012;30:413–21.

    Article  CAS  Google Scholar 

  30. Newman AM, Liu CL, Green MR, Gentles AJ, Feng W, Xu Y, Hoang CD, Diehn M, Alizadeh AA. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. 2015;12:453–7.

    Article  CAS  Google Scholar 

  31. Mroz EA, Rocco JW. MATH, a novel measure of intratumor genetic heterogeneity, is high in poor-outcome classes of head and neck squamous cell carcinoma. Oral Oncol. 2013;49:211–5.

    Article  CAS  Google Scholar 

  32. McGranahan N, Favero F, De Bruin EC, Birkbak NJ, Szallasi Z, Swanton C, De Bruin EC, Birkbak NJ, Szallasi Z, Swanton C, De Bruin EC, Birkbak NJ, Szallasi Z, Swanton C, De Bruin EC, Birkbak NJ, Szallasi Z, Swanton C, De Bruin EC, Birkbak NJ, Szallasi Z, Swanton C, De Bruin EC, Birkbak NJ, Szallasi Z, Swanton C. Clonal status of actionable driver events and the timing of mutational processes in cancer evolution. Sci Transl Med. 2015;7:1–12.

    Article  Google Scholar 

  33. Aran D, Sirota M, Butte AJ. Systematic pan-cancer analysis of tumour purity. Nat Commun. 2015;6:1–11.

    Article  Google Scholar 

  34. Kops GJPL, Weaver BAA, Cleveland DW. On the road to cancer: Aneuploidy and the mitotic checkpoint. Nat Rev Cancer. 2005;5:773–85.

    Article  CAS  Google Scholar 

  35. Genotype T, Expression T. The GTEx Consortium atlas of genetic regulatory effects across human tissues The Genotype Tissue Expression Consortium. 2019. https://doi.org/10.1101/787903.

    Article  Google Scholar 

  36. GTEx Consortium, Genetic effects on gene expression across human tissues. Nature. 550, 204–213 (2017).

  37. A. Carbone, L. Bernardini, F. Valenzano, I. Bottillo, C. De Simone, R. Capizzi, A. Capalbo, F. Romano, A. Novelli, B. Dallapiccola, P. Amerio, Array-based comparative genomic hybridization in early-stage mycosis fungoides: Recurrent deletion of tumor suppressor genes BCL7A, SMAC/DIABLO , and RHOF. Genes, Chromosom. Cancer. 47, 1067–1075 (2008).

  38. Sondka Z, Bamford S, Cole CG, Ward SA, Dunham I, Forbes SA. The COSMIC Cancer Gene Census: describing genetic dysfunction across all human cancers. Nat Rev Cancer. 2018;18:696–705.

    Article  CAS  Google Scholar 

  39. Peinado H, Zhang H, Matei IR, Costa-Silva B, Hoshino A, Rodrigues G, Psaila B, Kaplan RN, Bromberg JF, Kang Y, Bissell MJ, Cox TR, Giaccia AJ, Erler JT, Hiratsuka S, Ghajar CM, Lyden D. Pre-metastatic niches: organ-specific homes for metastases. Nat Rev Cancer. 2017;17:302–17.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are grateful for feedback from the Oncology Data Science group at Novartis Institutes for Biomedical Research.

Funding

The project was funded by Novartis using publicly available data.

Author information

Authors and Affiliations

  1. Novartis Institutes for Biomedical Research, 250 Massachusetts Avenue, Cambridge, MA, 02139, USA

    Bernard Pereira, Emma Labrot, Joshua M. Korn & Catarina D. Campbell

  2. Novartis Institutes for Biomedical Research, Novartis Campus, Fabrikstrasse 2, CH-4056, Basel, Switzerland

    Eric Durand & Audrey Kauffmann

Authors
  1. Bernard Pereira
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Emma Labrot
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Eric Durand
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Joshua M. Korn
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Audrey Kauffmann
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Catarina D. Campbell
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

B.P. performed analyses and interpreted results. B.P. and C.D.C. designed the study. E.L. downloaded and formatted data and interpreted results. E.D., J.M.K., A.K. and C.D.C. helped interpret results and develop hypotheses. All authors have read and contributed to the drafting and revision of the manuscript. “The author(s) read and approved the final manuscript.”

Corresponding author

Correspondence to Catarina D. Campbell.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

All authors are employees of Novartis.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pereira, B., Labrot, E., Durand, E. et al. Contribution and clinical relevance of germline variation to the cancer transcriptome. BMC Cancer 22, 675 (2022). https://doi.org/10.1186/s12885-022-09757-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12885-022-09757-0

Keywords

  • eQTL
  • TCGA
  • Cancer genomics
  • Germline variants

BMC Cancer

ISSN: 1471-2407

Contact us