Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions.

BACKGROUND: The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. METHODS: We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. RESULTS: We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. CONCLUSIONS: Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Large-scale cross-cancer fine-mapping of the 5p15.33 region reveals multiple independent signals.

Genome-wide association studies (GWASs) have identified thousands of cancer risk loci revealing many risk regions shared across multiple cancers. Characterizing the cross-cancer shared genetic basis can increase our understanding of global mechanisms of cancer development. In this study, we collected GWAS summary statistics based on up to 375,468 cancer cases and 530,521 controls for fourteen types of cancer, including breast (overall, estrogen receptor [ER]-positive, and ER-negative), colorectal, endometrial, esophageal, glioma, head/neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancer, to characterize the shared genetic basis of cancer risk. We identified thirteen pairs of cancers with statistically significant local genetic correlations across eight distinct genomic regions. Specifically, the 5p15.33 region, harboring the TERT and CLPTM1L genes, showed statistically significant local genetic correlations for multiple cancer pairs. We conducted a cross-cancer fine-mapping of the 5p15.33 region based on eight cancers that showed genome-wide significant associations in this region (ER-negative breast, colorectal, glioma, lung, melanoma, ovarian, pancreatic, and prostate cancer). We used an iterative analysis pipeline implementing a subset-based meta-analysis approach based on cancer-specific conditional analyses and identified ten independent cross-cancer associations within this region. For each signal, we conducted cross-cancer fine-mapping to prioritize the most plausible causal variants. Our findings provide a more in-depth understanding of the shared inherited basis across human cancers and expand our knowledge of the 5p15.33 region in carcinogenesis.

Multitrait transcriptome-wide association study (TWAS) tests.

Multitrait tests can improve power to detect associations between individual single-nucleotide polymorphisms (SNPs) and several related traits. Here, we develop methods for multi-SNP transcriptome-wide association (TWAS) tests to test the association between predicted gene expression levels and multiple phenotypes. We show that the correlation in TWAS test statistics for multiple phenotypes has the same form as multitrait statistics for the single-SNP setting. Thus, established methods for combining single-SNP test statistics across multiple traits can be extended directly to the TWAS setting. We performed an extensive evaluation across eight multitrait methods in simulations that varied gene-phenotype effect sizes in addition to the underlying covariance structure among the phenotypes. We found that all multitrait TWAS tests have well-calibrated Type I error (except ASSET, which can have a slightly elevated or depressed Type I error rate). Our results show that multitrait TWAS can improve statistical power compared with multiple single-trait TWAS followed by Bonferroni correction. To illustrate our approach to real data, we conducted a multitrait TWAS of four circulating lipid traits from the Global Lipids Genetics Consortium. We found that our multitrait Wald TWAS approach identified 506 genes associated with lipid levels compared with 87 identified through Bonferroni-corrected single-trait TWAS. Overall, we find that our proposed multitrait TWAS framework outperforms single-trait approaches to identify new genetic associations, especially for functionally correlated phenotypes and phenotypes with overlapping genome-wide association studies samples, leading to insights into the genetic architecture of multiple phenotypes.

AAV-Txnip prolongs cone survival and vision in mouse models of retinitis pigmentosa.

Retinitis pigmentosa (RP) is an inherited retinal disease affecting >20 million people worldwide. Loss of daylight vision typically occurs due to the dysfunction/loss of cone photoreceptors, the cell type that initiates our color and high-acuity vision. Currently, there is no effective treatment for RP, other than gene therapy for a limited number of specific disease genes. To develop a disease gene-agnostic therapy, we screened 20 genes for their ability to prolong cone photoreceptor survival in vivo. Here, we report an adeno-associated virus vector expressing Txnip, which prolongs the survival of cone photoreceptors and improves visual acuity in RP mouse models. A Txnip allele, C247S, which blocks the association of Txnip with thioredoxin, provides an even greater benefit. Additionally, the rescue effect of Txnip depends on lactate dehydrogenase b (Ldhb) and correlates with the presence of healthier mitochondria, suggesting that Txnip saves RP cones by enhancing their lactate catabolism.

Retinitis pigmentosa is an inherited eye disease affecting around one in every 4,000 people. It results from genetic defects in light sensitive cells of the retina, called photoreceptor cells, which line the back of the eye. Though vision loss can occur from birth, retinitis pigmentosa usually involves a gradual loss of vision, sometimes leading to blindness. Rod photoreceptors, which are responsible for vision in low light, are impacted first. The disease then affects cone photoreceptors, the cells that detect light during the day, providing both color and sharp vision. Around 100 mutated genes associated with retinitis pigmentosa have been identified, but only a handful of families with one of these mutant genes have been treated with a gene therapy specific for their mutated gene. There are currently no therapies available to treat the vast number of people with this disease. The mutations that cause retinitis pigmentosa directly affect the rod cells that detect dim light, leading to loss of night vision. There is also an indirect effect that causes cone photoreceptors to stop working and die. One theory to explain this two-step disease process relates to the fact that cone photoreceptors are very active cells, requiring a high level of energy, nutrients and oxygen. If surrounding rod cells die, cone photoreceptors may be deprived of some essential supplies, leading to cone cell death and daylight vision loss. To examine this theory, Xue et al. tested a new gene therapy designed to alleviate the potential shortfall in nutrients. The experiments used three different strains of mice that had the same genetic mutations as humans with retinitis pigmentosa. The gene therapy used a virus, called adeno-associated virus (AAV), to deliver 20 different genes to cone cells. Each of the 20 genes tested plays a different role in cells' processing of nutrients to provide energy. After administering the treatment, Xue et al. monitored the mice to see whether or not their vision was affected, and how cone cells responded. Only one of the 20 genes, Txnip, delivered using gene therapy, had a beneficial effect, prolonging cone cell survival in all three mouse strains. The mice that received Txnip also retained their ability to discern moving stripes on vision tests. Further investigations demonstrated that activating Txnip forced the cones to start using a molecule called lactate as an energy source, which could be more available to them than glucose, their usual fuel. These cells also had healthier mitochondria ­ the compartments inside cells that produce and manage energy supplies. This dual effect on fuel use and mitochondrial health is thought to be the basis for the extended cone survival and function. These experiments by Xue et al. have identified a good gene therapy candidate for treating retinitis pigmentosa independently of which genes are causing the disease. Further research will be required to test the safety of the gene therapy, and whether its beneficial effects translate to humans with retinitis pigmentosa, and potentially other diseases with unhealthy photoreceptors.

Leveraging expression from multiple tissues using sparse canonical correlation analysis and aggregate tests improves the power of transcriptome-wide association studies.

Transcriptome-wide association studies (TWAS) test the association between traits and genetically predicted gene expression levels. The power of a TWAS depends in part on the strength of the correlation between a genetic predictor of gene expression and the causally relevant gene expression values. Consequently, TWAS power can be low when expression quantitative trait locus (eQTL) data used to train the genetic predictors have small sample sizes, or when data from causally relevant tissues are not available. Here, we propose to address these issues by integrating multiple tissues in the TWAS using sparse canonical correlation analysis (sCCA). We show that sCCA-TWAS combined with single-tissue TWAS using an aggregate Cauchy association test (ACAT) outperforms traditional single-tissue TWAS. In empirically motivated simulations, the sCCA+ACAT approach yielded the highest power to detect a gene associated with phenotype, even when expression in the causal tissue was not directly measured, while controlling the Type I error when there is no association between gene expression and phenotype. For example, when gene expression explains 2% of the variability in outcome, and the GWAS sample size is 20,000, the average power difference between the ACAT combined test of sCCA features and single-tissue, versus single-tissue combined with Generalized Berk-Jones (GBJ) method, single-tissue combined with S-MultiXcan, UTMOST, or summarizing cross-tissue expression patterns using Principal Component Analysis (PCA) approaches was 5%, 8%, 5% and 38%, respectively. The gain in power is likely due to sCCA cross-tissue features being more likely to be detectably heritable. When applied to publicly available summary statistics from 10 complex traits, the sCCA+ACAT test was able to increase the number of testable genes and identify on average an additional 400 additional gene-trait associations that single-trait TWAS missed. Our results suggest that aggregating eQTL data across multiple tissues using sCCA can improve the sensitivity of TWAS while controlling for the false positive rate.

Transcriptome-wide association study of breast cancer risk by estrogen-receptor status.

Previous transcriptome-wide association studies (TWAS) have identified breast cancer risk genes by integrating data from expression quantitative loci and genome-wide association studies (GWAS), but analyses of breast cancer subtype-specific associations have been limited. In this study, we conducted a TWAS using gene expression data from GTEx and summary statistics from the hitherto largest GWAS meta-analysis conducted for breast cancer overall, and by estrogen receptor subtypes (ER+ and ER-). We further compared associations with ER+ and ER- subtypes, using a case-only TWAS approach. We also conducted multigene conditional analyses in regions with multiple TWAS associations. Two genes, STXBP4 and HIST2H2BA, were specifically associated with ER+ but not with ER- breast cancer. We further identified 30 TWAS-significant genes associated with overall breast cancer risk, including four that were not identified in previous studies. Conditional analyses identified single independent breast-cancer gene in three of six regions harboring multiple TWAS-significant genes. Our study provides new information on breast cancer genetics and biology, particularly about genomic differences between ER+ and ER- breast cancer.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer.

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Transcriptome-wide association studies accounting for colocalization using Egger regression.

Integrating genome-wide association (GWAS) and expression quantitative trait locus (eQTL) data into transcriptome-wide association studies (TWAS) based on predicted expression can boost power to detect novel disease loci or pinpoint the susceptibility gene at a known disease locus. However, it is often the case that multiple eQTL genes colocalize at disease loci, making the identification of the true susceptibility gene challenging, due to confounding through linkage disequilibrium (LD). To distinguish between true susceptibility genes (where the genetic effect on phenotype is mediated through expression) and colocalization due to LD, we examine an extension of the Mendelian randomization (MR) egger regression method that allows for LD while only requiring summary association data for both GWAS and eQTL. We derive the standard TWAS approach in the context of MR and show in simulations that the standard TWAS does not control type I error for causal gene identification when eQTLs have pleiotropic or LD-confounded effects on disease. In contrast, LD-aware MR-Egger (LDA MR-Egger) regression can control type I error in this case while attaining similar power as other methods in situations where these provide valid tests. However, when the direct effects of genetic variants on traits are correlated with the eQTL associations, all of the methods we examined including LDA MR-Egger regression can have inflated type I error. We illustrate these methods by integrating gene expression within a recent large-scale breast cancer GWAS to provide guidance on susceptibility gene identification.