SpliceAI-10k calculator for the prediction of pseudoexonization, intron retention, and exon deletion.

SUMMARY: SpliceAI is a widely used splicing prediction tool and its most common application relies on the maximum delta score to assign variant impact on splicing. We developed the SpliceAI-10k calculator (SAI-10k-calc) to extend use of this tool to predict: the splicing aberration type including pseudoexonization, intron retention, partial exon deletion, and (multi)exon skipping using a 10 kb analysis window; the size of inserted or deleted sequence; the effect on reading frame; and the altered amino acid sequence. SAI-10k-calc has 95% sensitivity and 96% specificity for predicting variants that impact splicing, computed from a control dataset of 1212 single-nucleotide variants (SNVs) with curated splicing assay results. Notably, it has high performance (≥84% accuracy) for predicting pseudoexon and partial intron retention. The automated amino acid sequence prediction allows for efficient identification of variants that are expected to result in mRNA nonsense-mediated decay or translation of truncated proteins. AVAILABILITY AND IMPLEMENTATION: SAI-10k-calc is implemented in R (https://github.com/adavi4/SAI-10k-calc) and also available as a Microsoft Excel spreadsheet. Users can adjust the default thresholds to suit their target performance values.

Association analysis identifies 65 new breast cancer risk loci.

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

The splicing effect of variants at branchpoint elements in cancer genes.

PURPOSE: Branchpoint elements are required for intron removal, and variants at these elements can result in aberrant splicing. We aimed to assess the value of branchpoint annotations generated from recent large-scale studies to select branchpoint-abrogating variants, using hereditary cancer genes as model. METHODS: We identified branchpoint elements in 119 genes associated with hereditary cancer from 3 genome-wide experimentally-inferred and 2 predicted branchpoint data sets. We then identified variants that occur within branchpoint elements from public databases. We compared conservation, unique variant observations, and population frequencies at different nucleotides within branchpoint motifs. Finally, selected minigene assays were performed to assess the splicing effect of variants at branchpoint elements within mismatch repair genes. RESULTS: There was poor overlap between predicted and experimentally-inferred branchpoints. Our analysis of cancer genes suggested that variants at -2 nucleotide, -1 nucleotide, and branchpoint positions in experimentally-inferred canonical motifs are more likely to be clinically relevant. Minigene assay data showed the -2 nucleotide to be more important to branchpoint motif integrity but also showed fluidity in branchpoint usage. CONCLUSION: Data from cancer gene analysis suggest that there are few high-risk alleles that severely impact function via branchpoint abrogation. Results of this study inform a general scheme to prioritize branchpoint motif variants for further study.

Genetic analyses of gynecological disease identify genetic relationships between uterine fibroids and endometrial cancer, and a novel endometrial cancer genetic risk region at the WNT4 1p36.12 locus.

Endometriosis, polycystic ovary syndrome (PCOS) and uterine fibroids have been proposed as endometrial cancer risk factors; however, disentangling their relationships with endometrial cancer is complicated due to shared risk factors and comorbidities. Using genome-wide association study (GWAS) data, we explored the relationships between these non-cancerous gynecological diseases and endometrial cancer risk by assessing genetic correlation, causal relationships and shared risk loci. We found significant genetic correlation between endometrial cancer and PCOS, and uterine fibroids. Adjustment for genetically predicted body mass index (a risk factor for PCOS, uterine fibroids and endometrial cancer) substantially attenuated the genetic correlation between endometrial cancer and PCOS but did not affect the correlation with uterine fibroids. Mendelian randomization analyses suggested a causal relationship between only uterine fibroids and endometrial cancer. Gene-based analyses revealed risk regions shared between endometrial cancer and endometriosis, and uterine fibroids. Multi-trait GWAS analysis of endometrial cancer and the genetically correlated gynecological diseases identified a novel genome-wide significant endometrial cancer risk locus at 1p36.12, which replicated in an independent endometrial cancer dataset. Interrogation of functional genomic data at 1p36.12 revealed biologically relevant genes, including WNT4 which is necessary for the development of the female reproductive system. In summary, our study provides genetic evidence for a causal relationship between uterine fibroids and endometrial cancer. It further provides evidence that the comorbidity of endometrial cancer, PCOS and uterine fibroids may partly be due to shared genetic architecture. Notably, this shared architecture has revealed a novel genome-wide risk locus for endometrial cancer.

Variant effect on splicing regulatory elements, branchpoint usage, and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars.

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.

Two integrated and highly predictive functional analysis-based procedures for the classification of MSH6 variants in Lynch syndrome.

PURPOSE: Variants in the DNA mismatch repair (MMR) gene MSH6, identified in individuals suspected of Lynch syndrome, are difficult to classify owing to the low cancer penetrance of defects in that gene. This not only obfuscates personalized health care but also the development of a rapid and reliable classification procedure that does not require clinical data. METHODS: The complete in vitro MMR activity (CIMRA) assay was calibrated against clinically classified MSH6 variants and, employing Bayes' rule, integrated with computational predictions of pathogenicity. To enable the validation of this two-component classification procedure we have employed a genetic screen to generate a large set of inactivating Msh6 variants, as proxies for pathogenic variants. RESULTS: The genetic screen-derived variants established that the two-component classification procedure displays high sensitivities and specificities. Moreover, these inactivating variants enabled the direct reclassification of human variants of uncertain significance (VUS) as (likely) pathogenic. CONCLUSION: The two-component classification procedure and the genetic screens provide complementary approaches to rapidly and cost-effectively classify the large majority of human MSH6 variants. The approach followed here provides a template for the classification of variants in other disease-predisposing genes, facilitating the translation of personalized genomics into personalized health care.

A functional assay-based procedure to classify mismatch repair gene variants in Lynch syndrome.

PURPOSE: To enhance classification of variants of uncertain significance (VUS) in the DNA mismatch repair (MMR) genes in the cancer predisposition Lynch syndrome, we developed the cell-free in vitro MMR activity (CIMRA) assay. Here, we calibrate and validate the assay, enabling its integration with in silico and clinical data. METHODS: Two sets of previously classified MLH1 and MSH2 variants were selected from a curated MMR gene database, and their biochemical activity determined by the CIMRA assay. The assay was calibrated by regression analysis followed by symmetric cross-validation and Bayesian integration with in silico predictions of pathogenicity. CIMRA assay reproducibility was assessed in four laboratories. RESULTS: Concordance between the training runs met our prespecified validation criterion. The CIMRA assay alone correctly classified 65% of variants, with only 3% discordant classification. Bayesian integration with in silico predictions of pathogenicity increased the proportion of correctly classified variants to 87%, without changing the discordance rate. Interlaboratory results were highly reproducible. CONCLUSION: The CIMRA assay accurately predicts pathogenic and benign MMR gene variants. Quantitative combination of assay results with in silico analysis correctly classified the majority of variants. Using this calibration, CIMRA assay results can be integrated into the diagnostic algorithm for MMR gene variants.

Long-Range Modulation of PAG1 Expression by 8q21 Allergy Risk Variants.

The gene(s) whose expression is regulated by allergy risk variants is unknown for many loci identified through genome-wide association studies. Addressing this knowledge gap might point to new therapeutic targets for allergic disease. The aim of this study was to identify the target gene(s) and the functional variant(s) underlying the association between rs7009110 on chromosome 8q21 and allergies. Eight genes are located within 1 Mb of rs7009110. Multivariate association analysis of publicly available exon expression levels from lymphoblastoid cell lines (LCLs) identified a significant association between rs7009110 and the expression of a single gene, PAG1 (p = 0.0017), 732 kb away. Analysis of histone modifications and DNase I hypersensitive sites in LCLs identified four putative regulatory elements (PREs) in the region. Chromosome conformation capture confirmed that two PREs interacted with the PAG1 promoter, one in allele-specific fashion. To determine whether these PREs were functional, LCLs were transfected with PAG1 promoter-driven luciferase reporter constructs. PRE3 acted as a transcriptional enhancer for PAG1 exclusively when it carried the rs2370615:C allergy predisposing allele, a variant in complete linkage disequilibrium with rs7009110. As such, rs2370615, which overlaps RelA transcription factor (TF) binding in LCLs and was found to disrupt Foxo3a binding to PRE3, represents the putative functional variant in this locus. Our studies suggest that the risk-associated allele of rs2370615 predisposes to allergic disease by increasing PAG1 expression, which might promote B cell activation and have a pro-inflammatory effect. Inhibition of PAG1 expression or function might have therapeutic potential for allergic diseases.

Fine-scale mapping of the 5q11.2 breast cancer locus reveals at least three independent risk variants regulating MAP3K1.

Genome-wide association studies (GWASs) have revealed SNP rs889312 on 5q11.2 to be associated with breast cancer risk in women of European ancestry. In an attempt to identify the biologically relevant variants, we analyzed 909 genetic variants across 5q11.2 in 103,991 breast cancer individuals and control individuals from 52 studies in the Breast Cancer Association Consortium. Multiple logistic regression analyses identified three independent risk signals: the strongest associations were with 15 correlated variants (iCHAV1), where the minor allele of the best candidate, rs62355902, associated with significantly increased risks of both estrogen-receptor-positive (ER(+): odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.21-1.27, ptrend = 5.7 × 10(-44)) and estrogen-receptor-negative (ER(-): OR = 1.10, 95% CI = 1.05-1.15, ptrend = 3.0 × 10(-4)) tumors. After adjustment for rs62355902, we found evidence of association of a further 173 variants (iCHAV2) containing three subsets with a range of effects (the strongest was rs113317823 [pcond = 1.61 × 10(-5)]) and five variants composing iCHAV3 (lead rs11949391; ER(+): OR = 0.90, 95% CI = 0.87-0.93, pcond = 1.4 × 10(-4)). Twenty-six percent of the prioritized candidate variants coincided with four putative regulatory elements that interact with the MAP3K1 promoter through chromatin looping and affect MAP3K1 promoter activity. Functional analysis indicated that the cancer risk alleles of four candidates (rs74345699 and rs62355900 [iCHAV1], rs16886397 [iCHAV2a], and rs17432750 [iCHAV3]) increased MAP3K1 transcriptional activity. Chromatin immunoprecipitation analysis revealed diminished GATA3 binding to the minor (cancer-protective) allele of rs17432750, indicating a mechanism for its action. We propose that the cancer risk alleles act to increase MAP3K1 expression in vivo and might promote breast cancer cell survival.

Discovery and functional assessment of gene variants in the vascular endothelial growth factor pathway.

Angiogenesis is a host-mediated mechanism in disease pathophysiology. The vascular endothelial growth factor (VEGF) pathway is a major determinant of angiogenesis, and a comprehensive annotation of the functional variation in this pathway is essential to understand the genetic basis of angiogenesis-related diseases. We assessed the allelic heterogeneity of gene expression, population specificity of cis expression quantitative trait loci (eQTLs), and eQTL function in luciferase assays in CEU and Yoruba people of Ibadan, Nigeria (YRI) HapMap lymphoblastoid cell lines in 23 resequenced genes. Among 356 cis-eQTLs, 155 and 174 were unique to CEU and YRI, respectively, and 27 were shared between CEU and YRI. Two cis-eQTLs provided mechanistic evidence for two genome-wide association study findings. Five eQTLs were tested for function in luciferase assays and the effect of two KRAS variants was concordant with the eQTL effect. Two eQTLs found in each of PRKCE, PIK3C2A, and MAP2K6 could predict 44%, 37%, and 45% of the variance in gene expression, respectively. This is the first analysis focusing on the pattern of functional genetic variation of the VEGF pathway genes in CEU and YRI populations and providing mechanistic evidence for genetic association studies of diseases for which angiogenesis plays a pathophysiologic role.

Dietary Factors and Endometrial Cancer Risk: A Mendelian Randomization Study.

Given the strong association between obesity and endometrial cancer risk, dietary factors may play an important role in the development of this cancer. However, observational studies of micro- and macronutrients and their role in endometrial cancer risk have been inconsistent. Clarifying these relationships are important to develop nutritional recommendations for cancer prevention. We performed two-sample Mendelian randomization (MR) to investigate the effects of circulating levels of 15 micronutrients (vitamin A (retinol), folate, vitamin B6, vitamin B12, vitamin C, vitamin D, vitamin E, ß-carotene, calcium, copper, iron, magnesium, phosphorus, selenium, and zinc) as well as corrected relative macronutrient intake (protein, carbohydrate, sugar and fat) on risks of endometrial cancer and its subtypes (endometrioid and non-endometrioid histologies). Genetically predicted vitamin C levels were found to be strongly associated with endometrial cancer risk. There was some evidence that genetically predicted relative intake of macronutrients (carbohydrate, sugar and fat) affects endometrial cancer risk. No other significant association were observed. Conclusions: In summary, these findings suggest that vitamin C and macronutrients influence endometrial cancer risk but further investigation is required.

Splicing annotation of endometrial cancer GWAS risk loci reveals potentially causal variants and supports a role for NF1 and SKAP1 as susceptibility genes.

Alternative splicing contributes to cancer development. Indeed, splicing analysis of cancer genome-wide association study (GWAS) risk variants has revealed likely causal variants. To systematically assess GWAS variants for splicing effects, we developed a prioritization workflow using a combination of splicing prediction tools, alternative transcript isoforms, and splicing quantitative trait locus (sQTL) annotations. Application of this workflow to candidate causal variants from 16 endometrial cancer GWAS risk loci highlighted single-nucleotide polymorphisms (SNPs) that were predicted to upregulate alternative transcripts. For two variants, sQTL data supported the predicted impact on splicing. At the 17q11.2 locus, the protective allele for rs7502834 was associated with increased splicing of an exon in a NF1 alternative transcript encoding a truncated protein in adipose tissue and is consistent with an endometrial cancer transcriptome-wide association study (TWAS) finding in adipose tissue. Notably, NF1 haploinsufficiency is protective for obesity, a well-established risk factor for endometrial cancer. At the 17q21.32 locus, the rs2278868 risk allele was predicted to upregulate a SKAP1 transcript that is subject to nonsense-mediated decay, concordant with a corresponding sQTL in lymphocytes. This is consistent with a TWAS finding that indicates decreased SKAP1 expression in blood increases endometrial cancer risk. As SKAP1 is involved in T cell immune responses, decreased SKAP1 expression may impact endometrial tumor immunosurveillance. In summary, our analysis has identified potentially causal endometrial cancer GWAS risk variants with plausible biological mechanisms and provides a splicing annotation workflow to aid interpretation of other GWAS datasets.

Multi-trait genome-wide association study identifies a novel endometrial cancer risk locus that associates with testosterone levels.

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

The clinical utility and costs of whole-genome sequencing to detect cancer susceptibility variants-a multi-site prospective cohort study.

BACKGROUND: Many families and individuals do not meet criteria for a known hereditary cancer syndrome but display unusual clusters of cancers. These families may carry pathogenic variants in cancer predisposition genes and be at higher risk for developing cancer. METHODS: This multi-centre prospective study recruited 195 cancer-affected participants suspected to have a hereditary cancer syndrome for whom previous clinical targeted genetic testing was either not informative or not available. To identify pathogenic disease-causing variants explaining participant presentation, germline whole-genome sequencing (WGS) and a comprehensive cancer virtual gene panel analysis were undertaken. RESULTS: Pathogenic variants consistent with the presenting cancer(s) were identified in 5.1% (10/195) of participants and pathogenic variants considered secondary findings with potential risk management implications were identified in another 9.7% (19/195) of participants. Health economic analysis estimated the marginal cost per case with an actionable variant was significantly lower for upfront WGS with virtual panel ($8744AUD) compared to standard testing followed by WGS ($24,894AUD). Financial analysis suggests that national adoption of diagnostic WGS testing would require a ninefold increase in government annual expenditure compared to conventional testing. CONCLUSIONS: These findings make a case for replacing conventional testing with WGS to deliver clinically important benefits for cancer patients and families. The uptake of such an approach will depend on the perspectives of different payers on affordability.

10 Years of GWAS discovery in endometrial cancer: Aetiology, function and translation.

Endometrial cancer is a common gynaecological cancer with increasing incidence and mortality. In the last decade, endometrial cancer genome-wide association studies (GWAS) have provided a resource to explore aetiology and for functional interpretation of heritable risk variation, informing endometrial cancer biology. Indeed, GWAS data have been used to assess relationships with other traits through correlation and Mendelian randomisation analyses, establishing genetic relationships and potential risk factors. Cross-trait GWAS analyses have increased statistical power and identified novel endometrial cancer risk variation related to other traits. Functional analysis of risk loci has helped prioritise candidate susceptibility genes, revealing molecular mechanisms and networks. Lastly, risk scores generated using endometrial cancer GWAS data may allow for clinical translation through identification of patients at high risk of disease. In the next decade, this knowledge base should enable substantial progress in our understanding of endometrial cancer and, potentially, new approaches for its screening and treatment.

Multi-tissue transcriptome-wide association study identifies eight candidate genes and tissue-specific gene expression underlying endometrial cancer susceptibility.

Genome-wide association studies (GWAS) have revealed sixteen risk loci for endoemtrial cancer but the identification of candidate susceptibility genes remains challenging. Here, we perform transcriptome-wide association study (TWAS) analyses using the largest endometrial cancer GWAS and gene expression from six relevant tissues, prioritizing eight candidate endometrial cancer susceptibility genes, one of which (EEFSEC) is located at a potentially novel endometrial cancer risk locus. We also show evidence of biologically relevant tissue-specific expression associations for CYP19A1 (adipose), HEY2 (ovary) and SKAP1 (whole blood). A phenome-wide association study demonstrates associations of candidate susceptibility genes with anthropometric, cardiovascular, diabetes, bone health and sex hormone traits that are related to endometrial cancer risk factors. Lastly, analysis of TWAS data highlights candidate compounds for endometrial cancer repurposing. In summary, this study reveals endometrial cancer susceptibility genes, including those with evidence of tissue specificity, providing insights into endometrial cancer aetiology and avenues for therapeutic development.

Associations between Genetically Predicted Circulating Protein Concentrations and Endometrial Cancer Risk.

Endometrial cancer (EC) is the leading female reproductive tract malignancy in developed countries. Currently, genome-wide association studies (GWAS) have identified 17 risk loci for EC. To identify novel EC-associated proteins, we used previously reported protein quantitative trait loci for 1434 plasma proteins as instruments to evaluate associations between genetically predicted circulating protein concentrations and EC risk. We studied 12,906 cases and 108,979 controls of European descent included in the Endometrial Cancer Association Consortium, the Epidemiology of Endometrial Cancer Consortium, and the UK Biobank. We observed associations between genetically predicted concentrations of nine proteins and EC risk at a false discovery rate of <0.05 (p-values range from 1.14 × 10-10 to 3.04 × 10-4). Except for vascular cell adhesion protein 1, all other identified proteins were independent from known EC risk variants identified in EC GWAS. The respective odds ratios (95% confidence intervals) per one standard deviation increase in genetically predicted circulating protein concentrations were 1.21 (1.13, 1.30) for DNA repair protein RAD51 homolog 4, 1.27 (1.14, 1.42) for desmoglein-2, 1.14 (1.07, 1.22) for MHC class I polypeptide-related sequence B, 1.05 (1.02, 1.08) for histo-blood group ABO system transferase, 0.77 (0.68, 0.89) for intestinal-type alkaline phosphatase, 0.82 (0.74, 0.91) for carbohydrate sulfotransferase 15, 1.07 (1.03, 1.11) for D-glucuronyl C5-epimerase, and 1.07 (1.03, 1.10) for CD209 antigen. In conclusion, we identified nine potential EC-associated proteins. If validated by additional studies, our findings may contribute to understanding the pathogenesis of endometrial tumor development and identifying women at high risk of EC along with other EC risk factors and biomarkers.

Identification of a Locus Near ULK1 Associated With Progression-Free Survival in Ovarian Cancer.

BACKGROUND: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. METHODS: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. RESULTS: We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10-8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10-8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. CONCLUSIONS: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. IMPACT: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options.See related commentary by Peres and Monteiro, p. 1604.