Comprehensive splicing analysis of the alternatively spliced CHEK2 exons 8 and 10 reveals three enhancer/silencer-rich regions and 38 spliceogenic variants.

Splicing is controlled by a large set of regulatory elements (SREs) including splicing enhancers and silencers, which are involved in exon recognition. Variants at these motifs may dysregulate splicing and trigger loss-of-function transcripts associated with disease. Our goal here was to study the alternatively spliced exons 8 and 10 of the breast cancer susceptibility gene CHEK2. For this purpose, we used a previously published minigene with exons 6-10 that produced the expected minigene full-length transcript and replicated the naturally occurring events of exon 8 [Δ(E8)] and exon 10 [Δ(E10)] skipping. We then introduced 12 internal microdeletions of exons 8 and 10 by mutagenesis in order to map SRE-rich intervals by splicing assays in MCF-7 cells. We identified three minimal (10-, 11-, 15-nt) regions essential for exon recognition: c.863_877del [ex8, Δ(E8): 75%] and c.1073_1083del and c.1083_1092del [ex10, Δ(E10): 97% and 62%, respectively]. Then 87 variants found within these intervals were introduced into the wild-type minigene and tested functionally. Thirty-eight of them (44%) impaired splicing, four of which (c.883G>A, c.883G>T, c.884A>T, and c.1080G>T) induced negligible amounts (<5%) of the minigene full-length transcript. Another six variants (c.886G>A, c.886G>T, c.1075G>A, c.1075G>T, c.1076A>T, and c.1078G>T) showed significantly strong impacts (20-50% of the minigene full-length transcript). Thirty-three of the 38 spliceogenic variants were annotated as missense, three as nonsense, and two as synonymous, underlying the fact that any exonic change is capable of disrupting splicing. Moreover, c.883G>A, c.883G>T, and c.884A>T were classified as pathogenic/likely pathogenic variants according to ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based criteria. © 2024 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

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.

Systematic Minigene-Based Splicing Analysis and Tentative Clinical Classification of 52 CHEK2 Splice-Site Variants.

BACKGROUND: Disrupted pre-mRNA splicing is a frequent deleterious mechanism in hereditary cancer. We aimed to functionally analyze candidate spliceogenic variants of the breast cancer susceptibility gene CHEK2 by splicing reporter minigenes. METHODS: A total of 128 CHEK2 splice-site variants identified in the Breast Cancer After Diagnostic Gene Sequencing (BRIDGES) project (https://cordis.europa.eu/project/id/634935) were analyzed with MaxEntScan and subsetted to 52 variants predicted to impact splicing. Three CHEK2 minigenes, which span all 15 exons, were constructed and validated. The 52 selected variants were then genetically engineered into the minigenes and assayed in MCF-7 (human breast adenocarcinoma) cells. RESULTS: Of 52 variants, 46 (88.5%) impaired splicing. Some of them led to complex splicing patterns with up to 11 different transcripts. Thirty-four variants induced splicing anomalies without any trace or negligible amounts of the full-length transcript. A total of 89 different transcripts were annotated, which derived from different events: single- or multi-exon skipping, alternative site-usage, mutually exclusive exon inclusion, intron retention or combinations of the abovementioned events. Fifty-nine transcripts were predicted to introduce premature termination codons, 7 kept the original open-reading frame, 5 removed the translation start codon, 6 affected the 5'UTR (Untranslated Region), and 2 included missense variations. Analysis of variant c.684-2A > G revealed the activation of a non-canonical TG-acceptor site and exon 6 sequences critical for its recognition. CONCLUSIONS: Incorporation of minigene read-outs into an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme allowed us to classify 32 CHEK2 variants (27 pathogenic/likely pathogenic and 5 likely benign). However, 20 variants (38%) remained of uncertain significance, reflecting in part the complex splicing patterns of this gene.

Associations of height, body mass index, and weight gain with breast cancer risk in carriers of a pathogenic variant in BRCA1 or BRCA2: the BRCA1 and BRCA2 Cohort Consortium.

INTRODUCTION: Height, body mass index (BMI), and weight gain are associated with breast cancer risk in the general population. It is unclear whether these associations also exist for carriers of pathogenic variants in the BRCA1 or BRCA2 genes. PATIENTS AND METHODS: An international pooled cohort of 8091 BRCA1/2 variant carriers was used for retrospective and prospective analyses separately for premenopausal and postmenopausal women. Cox regression was used to estimate breast cancer risk associations with height, BMI, and weight change. RESULTS: In the retrospective analysis, taller height was associated with risk of premenopausal breast cancer for BRCA2 variant carriers (HR 1.20 per 10 cm increase, 95% CI 1.04-1.38). Higher young-adult BMI was associated with lower premenopausal breast cancer risk for both BRCA1 (HR 0.75 per 5 kg/m2, 95% CI 0.66-0.84) and BRCA2 (HR 0.76, 95% CI 0.65-0.89) variant carriers in the retrospective analysis, with consistent, though not statistically significant, findings from the prospective analysis. In the prospective analysis, higher BMI and adult weight gain were associated with higher postmenopausal breast cancer risk for BRCA1 carriers (HR 1.20 per 5 kg/m2, 95% CI 1.02-1.42; and HR 1.10 per 5 kg weight gain, 95% CI 1.01-1.19, respectively). CONCLUSION: Anthropometric measures are associated with breast cancer risk for BRCA1 and BRCA2 variant carriers, with relative risk estimates that are generally consistent with those for women from the general population.

Ovarian cancer pathology characteristics as predictors of variant pathogenicity in BRCA1 and BRCA2.

BACKGROUND: The distribution of ovarian tumour characteristics differs between germline BRCA1 and BRCA2 pathogenic variant carriers and non-carriers. In this study, we assessed the utility of ovarian tumour characteristics as predictors of BRCA1 and BRCA2 variant pathogenicity, for application using the American College of Medical Genetics and the Association for Molecular Pathology (ACMG/AMP) variant classification system. METHODS: Data for 10,373 ovarian cancer cases, including carriers and non-carriers of BRCA1 or BRCA2 pathogenic variants, were collected from unpublished international cohorts and consortia and published studies. Likelihood ratios (LR) were calculated for the association of ovarian cancer histology and other characteristics, with BRCA1 and BRCA2 variant pathogenicity. Estimates were aligned to ACMG/AMP code strengths (supporting, moderate, strong). RESULTS: No histological subtype provided informative ACMG/AMP evidence in favour of BRCA1 and BRCA2 variant pathogenicity. Evidence against variant pathogenicity was estimated for the mucinous and clear cell histologies (supporting) and borderline cases (moderate). Refined associations are provided according to tumour grade, invasion and age at diagnosis. CONCLUSIONS: We provide detailed estimates for predicting BRCA1 and BRCA2 variant pathogenicity based on ovarian tumour characteristics. This evidence can be combined with other variant information under the ACMG/AMP classification system, to improve classification and carrier clinical management.

Minigene-based splicing analysis and ACMG/AMP-based tentative classification of 56 ATM variants.

The ataxia telangiectasia-mutated (ATM) protein is a major coordinator of the DNA damage response pathway. ATM loss-of-function variants are associated with 2-fold increased breast cancer risk. We aimed at identifying and classifying spliceogenic ATM variants detected in subjects of the large-scale sequencing project BRIDGES. A total of 381 variants at the intron-exon boundaries were identified, 128 of which were predicted to be spliceogenic. After further filtering, we ended up selecting 56 variants for splicing analysis. Four functional minigenes (mgATM) spanning exons 4-9, 11-17, 25-29, and 49-52 were constructed in the splicing plasmid pSAD. Selected variants were genetically engineered into the four constructs and assayed in MCF-7/HeLa cells. Forty-eight variants (85.7%) impaired splicing, 32 of which did not show any trace of the full-length (FL) transcript. A total of 43 transcripts were identified where the most prevalent event was exon/multi-exon skipping. Twenty-seven transcripts were predicted to truncate the ATM protein. A tentative ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme that integrates mgATM data allowed us to classify 29 ATM variants as pathogenic/likely pathogenic and seven variants as likely benign. Interestingly, the likely pathogenic variant c.1898+2T>G generated 13% of the minigene FL-transcript due to the use of a noncanonical GG-5'-splice-site (0.014% of human donor sites). Circumstantial evidence in three ATM variants (leakiness uncovered by our mgATM analysis together with clinical data) provides some support for a dosage-sensitive expression model in which variants producing ≥30% of FL-transcripts would be predicted benign, while variants producing ≤13% of FL-transcripts might be pathogenic. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Splicing predictions, minigene analyses, and ACMG-AMP clinical classification of 42 germline PALB2 splice-site variants.

PALB2 loss-of-function variants confer high risk of developing breast cancer. Here we present a systematic functional analysis of PALB2 splice-site variants detected in approximately 113,000 women in the large-scale sequencing project Breast Cancer After Diagnostic Gene Sequencing (BRIDGES; https://bridges-research.eu/). Eighty-two PALB2 variants at the intron-exon boundaries were analyzed with MaxEntScan. Forty-two variants were selected for the subsequent splicing functional assays. For this purpose, three splicing reporter minigenes comprising exons 1-12 were constructed. The 42 potential spliceogenic variants were introduced into the minigenes by site-directed mutagenesis and assayed in MCF-7/MDA-MB-231 cells. Splicing anomalies were observed in 35 variants, 23 of which showed no traces or minimal amounts of the expected full-length transcripts of each minigene. More than 30 different variant-induced transcripts were characterized, 23 of which were predicted to truncate the PALB2 protein. The pathogenicity of all variants was interpreted according to an in-house adaptation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) variant classification scheme. Up to 23 variants were classified as pathogenic/likely pathogenic. Remarkably, three ±1,2 variants (c.49-2A>T, c.108+2T>C, and c.211+1G>A) were classified as variants of unknown significance, as they produced significant amounts of either in-frame transcripts of unknown impact on the PALB2 protein function or the minigene full-length transcripts. In conclusion, we have significantly contributed to the ongoing effort of identifying spliceogenic variants in the clinically relevant PALB2 cancer susceptibility gene. Moreover, we suggest some approaches to classify the findings in accordance with the ACMG-AMP rationale. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Clinical, splicing, and functional analysis to classify BRCA2 exon 3 variants: Application of a points-based ACMG/AMP approach.

Skipping of BRCA2 exon 3 (∆E3) is a naturally occurring splicing event, complicating clinical classification of variants that may alter ∆E3 expression. This study used multiple evidence types to assess pathogenicity of 85 variants in/near BRCA2 exon 3. Bioinformatically predicted spliceogenic variants underwent mRNA splicing analysis using minigenes and/or patient samples. ∆E3 was measured using quantitative analysis. A mouse embryonic stem cell (mESC) based assay was used to determine the impact of 18 variants on mRNA splicing and protein function. For each variant, population frequency, bioinformatic predictions, clinical data, and existing mRNA splicing and functional results were collated. Variant class was assigned using a gene-specific adaptation of ACMG/AMP guidelines, following a recently proposed points-based system. mRNA and mESC analysis combined identified six variants with transcript and/or functional profiles interpreted as loss of function. Cryptic splice site use for acceptor site variants generated a transcript encoding a shorter protein that retains activity. Overall, 69/85 (81%) variants were classified using the points-based approach. Our analysis shows the value of applying gene-specific ACMG/AMP guidelines using a points-based approach and highlights the consideration of cryptic splice site usage to appropriately assign PVS1 code strength.

Risks of breast and ovarian cancer for women harboring pathogenic missense variants in BRCA1 and BRCA2 compared with those harboring protein truncating variants.

PURPOSE: Germline genetic testing for BRCA1 and BRCA2 variants has been a part of clinical practice for >2 decades. However, no studies have compared the cancer risks associated with missense pathogenic variants (PVs) with those associated with protein truncating (PTC) variants. METHODS: We collected 582 informative pedigrees segregating 1 of 28 missense PVs in BRCA1 and 153 pedigrees segregating 1 of 12 missense PVs in BRCA2. We analyzed 324 pedigrees with PTC variants in BRCA1 and 214 pedigrees with PTC variants in BRCA2. Cancer risks were estimated using modified segregation analysis. RESULTS: Estimated breast cancer risks were markedly lower for women aged >50 years carrying BRCA1 missense PVs than for the women carrying BRCA1 PTC variants (hazard ratio [HR] = 3.9 [2.4-6.2] for PVs vs 12.8 [5.7-28.7] for PTC variants; P = .01), particularly for missense PVs in the BRCA1 C-terminal domain (HR = 2.8 [1.4-5.6]; P = .005). In case of BRCA2, for women aged >50 years, the HR was 3.9 (2.0-7.2) for those heterozygous for missense PVs compared with 7.0 (3.3-14.7) for those harboring PTC variants. BRCA1 p.[Cys64Arg] and BRCA2 p.[Trp2626Cys] were associated with particularly low risks of breast cancer compared with other PVs. CONCLUSION: These results have important implications for the counseling of at-risk women who harbor missense PVs in the BRCA1/2 genes.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants.

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

A Collaborative Effort to Define Classification Criteria for ATM Variants in Hereditary Cancer Patients.

BACKGROUND: Gene panel testing by massive parallel sequencing has increased the diagnostic yield but also the number of variants of uncertain significance. Clinical interpretation of genomic data requires expertise for each gene and disease. Heterozygous ATM pathogenic variants increase the risk of cancer, particularly breast cancer. For this reason, ATM is included in most hereditary cancer panels. It is a large gene, showing a high number of variants, most of them of uncertain significance. Hence, we initiated a collaborative effort to improve and standardize variant classification for the ATM gene. METHODS: Six independent laboratories collected information from 766 ATM variant carriers harboring 283 different variants. Data were submitted in a consensus template form, variant nomenclature and clinical information were curated, and monthly team conferences were established to review and adapt American College of Medical Genetics and Genomics/Association for Molecular Pathology (ACMG/AMP) criteria to ATM, which were used to classify 50 representative variants. RESULTS: Amid 283 different variants, 99 appeared more than once, 35 had differences in classification among laboratories. Refinement of ACMG/AMP criteria to ATM involved specification for twenty-one criteria and adjustment of strength for fourteen others. Afterwards, 50 variants carried by 254 index cases were classified with the established framework resulting in a consensus classification for all of them and a reduction in the number of variants of uncertain significance from 58% to 42%. CONCLUSIONS: Our results highlight the relevance of data sharing and data curation by multidisciplinary experts to achieve improved variant classification that will eventually improve clinical management.

Haplotype analysis of the internationally distributed BRCA1 c.3331_3334delCAAG founder mutation reveals a common ancestral origin in Iberia.

BACKGROUND: The BRCA1 c.3331_3334delCAAG founder mutation has been reported in hereditary breast and ovarian cancer families from multiple Hispanic groups. We aimed to evaluate BRCA1 c.3331_3334delCAAG haplotype diversity in cases of European, African, and Latin American ancestry. METHODS: BC mutation carrier cases from Colombia (n = 32), Spain (n = 13), Portugal (n = 2), Chile (n = 10), Africa (n = 1), and Brazil (n = 2) were genotyped with the genome-wide single nucleotide polymorphism (SNP) arrays to evaluate haplotype diversity around BRCA1 c.3331_3334delCAAG. Additional Portuguese (n = 13) and Brazilian (n = 18) BC mutation carriers were genotyped for 15 informative SNPs surrounding BRCA1. Data were phased using SHAPEIT2, and identical by descent regions were determined using BEAGLE and GERMLINE. DMLE+ was used to date the mutation in Colombia and Iberia. RESULTS: The haplotype reconstruction revealed a shared 264.4-kb region among carriers from all six countries. The estimated mutation age was ~ 100 generations in Iberia and that it was introduced to South America early during the European colonization period. CONCLUSIONS: Our results suggest that this mutation originated in Iberia and later introduced to Colombia and South America at the time of Spanish colonization during the early 1500s. We also found that the Colombian mutation carriers had higher European ancestry, at the BRCA1 gene harboring chromosome 17, than controls, which further supported the European origin of the mutation. Understanding founder mutations in diverse populations has implications in implementing cost-effective, ancestry-informed screening.

Alternative mRNA splicing can attenuate the pathogenicity of presumed loss-of-function variants in BRCA2.

PURPOSE: Current interpretation guidelines for germline variants in high-risk cancer susceptibility genes consider predicted loss-of-function (LoF) variants, such as nonsense variants and variants in the canonical splice site sequences ofBRCA2, to be associated with high cancer risk. However, some variant alleles produce alternative transcripts that encode (partially) functional protein isoforms leading to possible incorrect risk estimations. For accurate classification of variants it is therefore essential that alternative transcripts are identified and functionally characterized. METHODS: We systematically evaluated a large panel of human BRCA2 variants for the production of alternative transcripts and assessed their capacity to exert BRCA2 protein functionality. Evaluated variants included all single-exon deletions, various multiple-exon deletions, intronic variants at the canonical splice donor and acceptor sequences, and variants that previously have been shown to affect messenger RNA (mRNA) splicing in carriers. RESULTS: Multiple alternative transcripts encoding (partially) functional protein isoforms were identified (e.g., ∆[E4-E7], ∆[E6-E7], ∆E[6q39_E8], ∆[E10], ∆[E12], ∆E[12-14]). Expression of these transcripts did attenuate the impact of predicted LoF variants such as the canonical splice site variants c.631+2T>G, c.517-2A>G, c.6842-2A>G, c.6937+1G>A, and nonsense variants c.491T>A, c.581G>A, and c.6901G>T. CONCLUSION: These results allow refinement of variant interpretation guidelines for BRCA2 by providing insight into the functional consequences of naturally occurring and variant-related alternative splicing events.

Correction: Alternative mRNA splicing can attenuate the pathogenicity of presumed loss-of-function variants in BRCA2.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Polygenic risk scores and breast and epithelial ovarian cancer risks for carriers of BRCA1 and BRCA2 pathogenic variants.

PURPOSE: We assessed the associations between population-based polygenic risk scores (PRS) for breast (BC) or epithelial ovarian cancer (EOC) with cancer risks for BRCA1 and BRCA2 pathogenic variant carriers. METHODS: Retrospective cohort data on 18,935 BRCA1 and 12,339 BRCA2 female pathogenic variant carriers of European ancestry were available. Three versions of a 313 single-nucleotide polymorphism (SNP) BC PRS were evaluated based on whether they predict overall, estrogen receptor (ER)-negative, or ER-positive BC, and two PRS for overall or high-grade serous EOC. Associations were validated in a prospective cohort. RESULTS: The ER-negative PRS showed the strongest association with BC risk for BRCA1 carriers (hazard ratio [HR] per standard deviation = 1.29 [95% CI 1.25-1.33], P = 3×10-72). For BRCA2, the strongest association was with overall BC PRS (HR = 1.31 [95% CI 1.27-1.36], P = 7×10-50). HR estimates decreased significantly with age and there was evidence for differences in associations by predicted variant effects on protein expression. The HR estimates were smaller than general population estimates. The high-grade serous PRS yielded the strongest associations with EOC risk for BRCA1 (HR = 1.32 [95% CI 1.25-1.40], P = 3×10-22) and BRCA2 (HR = 1.44 [95% CI 1.30-1.60], P = 4×10-12) carriers. The associations in the prospective cohort were similar. CONCLUSION: Population-based PRS are strongly associated with BC and EOC risks for BRCA1/2 carriers and predict substantial absolute risk differences for women at PRS distribution extremes.

Towards controlled terminology for reporting germline cancer susceptibility variants: an ENIGMA report.

The vocabulary currently used to describe genetic variants and their consequences reflects many years of studying and discovering monogenic disease with high penetrance. With the recent rapid expansion of genetic testing brought about by wide availability of high-throughput massively parallel sequencing platforms, accurate variant interpretation has become a major issue. The vocabulary used to describe single genetic variants in silico, in vitro, in vivo and as a contributor to human disease uses terms in common, but the meaning is not necessarily shared across all these contexts. In the setting of cancer genetic tests, the added dimension of using data from genetic sequencing of tumour DNA to direct treatment is an additional source of confusion to those who are not experienced in cancer genetics. The language used to describe variants identified in cancer susceptibility genetic testing typically still reflects an outdated paradigm of Mendelian inheritance with dichotomous outcomes. Cancer is a common disease with complex genetic architecture; an improved lexicon is required to better communicate among scientists, clinicians and patients, the risks and implications of genetic variants detected. This review arises from a recognition of, and discussion about, inconsistencies in vocabulary usage by members of the ENIGMA international multidisciplinary consortium focused on variant classification in breast-ovarian cancer susceptibility genes. It sets out the vocabulary commonly used in genetic variant interpretation and reporting, and suggests a framework for a common vocabulary that may facilitate understanding and clarity in clinical reporting of germline genetic tests for cancer susceptibility.

Alternative splicing and ACMG-AMP-2015-based classification of PALB2 genetic variants: an ENIGMA report.

BACKGROUND: PALB2 monoallelic loss-of-function germ-line variants confer a breast cancer risk comparable to the average BRCA2 pathogenic variant. Recommendations for risk reduction strategies in carriers are similar. Elaborating robust criteria to identify loss-of-function variants in PALB2-without incurring overprediction-is thus of paramount clinical relevance. Towards this aim, we have performed a comprehensive characterisation of alternative splicing in PALB2, analysing its relevance for the classification of truncating and splice site variants according to the 2015 American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines. METHODS: Alternative splicing was characterised in RNAs extracted from blood, breast and fimbriae/ovary-related human specimens (n=112). RNAseq, RT-PCR/CE and CloneSeq experiments were performed by five contributing laboratories. Centralised revision/curation was performed to assure high-quality annotations. Additional splicing analyses were performed in PALB2 c.212-1G>A, c.1684+1G>A, c.2748+2T>G, c.3113+5G>A, c.3350+1G>A, c.3350+4A>C and c.3350+5G>A carriers. The impact of the findings on PVS1 status was evaluated for truncating and splice site variant. RESULTS: We identified 88 naturally occurring alternative splicing events (81 newly described), including 4 in-frame events predicted relevant to evaluate PVS1 status of splice site variants. We did not identify tissue-specific alternate gene transcripts in breast or ovarian-related samples, supporting the clinical relevance of blood-based splicing studies. CONCLUSIONS: PVS1 is not necessarily warranted for splice site variants targeting four PALB2 acceptor sites (exons 2, 5, 7 and 10). As a result, rare variants at these splice sites cannot be assumed pathogenic/likely pathogenic without further evidences. Our study puts a warning in up to five PALB2 genetic variants that are currently reported as pathogenic/likely pathogenic in ClinVar.

RECQL5: Another DNA helicase potentially involved in hereditary breast cancer susceptibility.

There is still around 50% of the familial breast cancer (BC) cases with an undefined genetic cause, here we have used next-generation sequencing (NGS) technology to identify new BC susceptibility genes. This approach has led to the identification of RECQL5, a member of RECQL-helicases family, as a new BC susceptibility candidate, which deserves further study. We have used a combination of whole exome sequencing in a family negative for mutations in BRCA1/2 throughout (BRCAX), in which we found a probably deleterious variant in RECQL5, and targeted NGS of the complete coding regions and exon-intron boundaries of the candidate gene in 699 BC Spanish BRCAX families and 665 controls. Functional characterization and in silico inference of pathogenicity were performed to evaluate the deleterious effect of detected variants. We found at least seven deleterious or likely deleterious variants among the cases and only one in controls. These results prompt us to propose RECQL5 as a gene that would be worth to analyze in larger studies to explore its possible implication in BC susceptibility.

Large scale multifactorial likelihood quantitative analysis of BRCA1 and BRCA2 variants: An ENIGMA resource to support clinical variant classification.

The multifactorial likelihood analysis method has demonstrated utility for quantitative assessment of variant pathogenicity for multiple cancer syndrome genes. Independent data types currently incorporated in the model for assessing BRCA1 and BRCA2 variants include clinically calibrated prior probability of pathogenicity based on variant location and bioinformatic prediction of variant effect, co-segregation, family cancer history profile, co-occurrence with a pathogenic variant in the same gene, breast tumor pathology, and case-control information. Research and clinical data for multifactorial likelihood analysis were collated for 1,395 BRCA1/2 predominantly intronic and missense variants, enabling classification based on posterior probability of pathogenicity for 734 variants: 447 variants were classified as (likely) benign, and 94 as (likely) pathogenic; and 248 classifications were new or considerably altered relative to ClinVar submissions. Classifications were compared with information not yet included in the likelihood model, and evidence strengths aligned to those recommended for ACMG/AMP classification codes. Altered mRNA splicing or function relative to known nonpathogenic variant controls were moderately to strongly predictive of variant pathogenicity. Variant absence in population datasets provided supporting evidence for variant pathogenicity. These findings have direct relevance for BRCA1 and BRCA2 variant evaluation, and justify the need for gene-specific calibration of evidence types used for variant classification.

Targeted RNA-seq successfully identifies normal and pathogenic splicing events in breast/ovarian cancer susceptibility and Lynch syndrome genes.

A subset of genetic variants found through screening of patients with hereditary breast and ovarian cancer syndrome (HBOC) and Lynch syndrome impact RNA splicing. Through target enrichment of the transcriptome, it is possible to perform deep-sequencing and to identify the different and even rare mRNA isoforms. A targeted RNA-seq approach was used to analyse the naturally-occurring splicing events for a panel of 8 breast and/or ovarian cancer susceptibility genes (BRCA1, BRCA2, RAD51C, RAD51D, PTEN, STK11, CDH1, TP53), 3 Lynch syndrome genes (MLH1, MSH2, MSH6) and the fanconi anaemia SLX4 gene, in which monoallelic mutations were found in non-BRCA families. For BRCA1, BRCA2, RAD51C and RAD51D the results were validated by capillary electrophoresis and were compared to a non-targeted RNA-seq approach. We also compared splicing events from lymphoblastoid cell-lines with those from breast and ovarian fimbriae tissues. The potential of targeted RNA-seq to detect pathogenic changes in RNA-splicing was validated by the inclusion of samples with previously well characterized BRCA1/2 genetic variants. In our study, we update the catalogue of normal splicing events for BRCA1/2, provide an extensive catalogue of normal RAD51C and RAD51D alternative splicing, and list splicing events found for eight other genes. Additionally, we show that our approach allowed the identification of aberrant splicing events due to the presence of BRCA1/2 genetic variants and distinguished between complete and partial splicing events. In conclusion, targeted-RNA-seq can be very useful to classify variants based on their putative pathogenic impact on splicing.