Combined Tumor-Based BRCA1/2 and TP53 Mutation Testing in Ovarian Cancer.

Somatic/germline BRCA1/2 mutations (m)/(likely) pathogenic variants (PV) (s/gBRCAm) remain the best predictive biomarker for PARP inhibitor efficacy. As >95% of high-grade serous ovarian cancers (HGSOC) have a somatic TP53m, combined tumor-based BRCA1/2 (tBRCA) and TP53 mutation testing (tBRCA/TP53m) may improve the quality of results in somatic BRCAm identification and interpretation of the 'second hit' event, i.e., loss of heterozygosity (LOH). A total of 237 patients with HGSOC underwent tBRCA/TP53m testing. The ratio of allelic fractions (AFs) for tBRCA/TP53m was calculated to estimate the proportion of cells carrying BRCAm and to infer LOH. Among the 142/237 gBRCA results, 16.2% demonstrated a pathogenic/deleterious variant (DEL) gBRCA1/2m. Among the 195 contributive tumor samples, 43 DEL of tBRCAm (22.1%) were identified (23 gBRCAm and 20 sBRCAm) with LOH identified in 37/41 conclusive samples. The median AF of TP53m was 0.52 (0.01-0.93), confirming huge variability in tumor cellularity. Initially, three samples were considered as wild type with <10% cellularity. However, additional testing detected a very low AF (<0.05) in both BRCA1/2m and TP53m, thus reidentifying them as sBRCA1/2m. Combined tBRCA/TP53m testing is rapid, sensitive, and identifies somatic and germline BRCA1/2m. AF TP53m is essential for interpreting sBRCA1/2m in low-cellularity samples and provides indirect evidence for LOH as the 'second hit' of BRCA1/2-related tumorigenesis.

SPiP: Splicing Prediction Pipeline, a machine learning tool for massive detection of exonic and intronic variant effects on mRNA splicing.

Modeling splicing is essential for tackling the challenge of variant interpretation as each nucleotide variation can be pathogenic by affecting pre-mRNA splicing via disruption/creation of splicing motifs such as 5'/3' splice sites, branch sites, or splicing regulatory elements. Unfortunately, most in silico tools focus on a specific type of splicing motif, which is why we developed the Splicing Prediction Pipeline (SPiP) to perform, in one single bioinformatic analysis based on a machine learning approach, a comprehensive assessment of the variant effect on different splicing motifs. We gathered a curated set of 4616 variants scattered all along the sequence of 227 genes, with their corresponding splicing studies. The Bayesian analysis provided us with the number of control variants, that is, variants without impact on splicing, to mimic the deluge of variants from high-throughput sequencing data. Results show that SPiP can deal with the diversity of splicing alterations, with 83.13% sensitivity and 99% specificity to detect spliceogenic variants. Overall performance as measured by area under the receiving operator curve was 0.986, better than SpliceAI and SQUIRLS (0.965 and 0.766) for the same data set. SPiP lends itself to a unique suite for comprehensive prediction of spliceogenicity in the genomic medicine era. SPiP is available at: https://sourceforge.net/projects/splicing-prediction-pipeline/.

The mitochondrial seryl-tRNA synthetase SARS2 modifies onset in spastic paraplegia type 4.

PURPOSE: Hereditary spastic paraplegia type 4 is extremely variable in age at onset; the same variant can cause onset at birth or in the eighth decade. We recently discovered that missense variants in SPAST, which influences microtubule dynamics, are associated with earlier onset and more severe disease than truncating variants, but even within the early and late-onset groups there remained significant differences in onset. Given the rarity of the condition, we adapted an extreme phenotype approach to identify genetic modifiers of onset. METHODS: We performed a genome-wide association study on 134 patients bearing truncating pathogenic variants in SPAST, divided into early- and late-onset groups (aged ≤15 and ≥45 years, respectively). A replication cohort of 419 included patients carrying either truncating or missense variants. Finally, age at onset was analyzed in the merged cohort (N = 553). RESULTS: We found 1 signal associated with earlier age at onset (rs10775533, P = 8.73E-6) in 2 independent cohorts and in the merged cohort (N = 553, Mantel-Cox test, P < .0001). Western blotting in lymphocytes of 20 patients showed that this locus tends to upregulate SARS2 expression in earlier-onset patients. CONCLUSION: SARS2 overexpression lowers the age of onset in hereditary spastic paraplegia type 4. Lowering SARS2 or improving mitochondrial function could thus present viable approaches to therapy.

Clinical and genetic spectra of 1550 index patients with hereditary spastic paraplegia.

Hereditary spastic paraplegia refers to rare genetic neurodevelopmental and/or neurodegenerative disorders in which spasticity due to length-dependent damage to the upper motor neuron is a core sign. Their high clinical and genetic heterogeneity makes their diagnosis challenging. Multigene panels allow a high-throughput targeted analysis of the increasing number of genes involved using next-generation sequencing. We report here the clinical and genetic results of 1550 index cases tested for variants in a panel of hereditary spastic paraplegia related genes analysed in routine diagnosis. A causative variant was found in 475 patients (30.7%) in 35/65 screened genes. SPAST and SPG7 were the most frequently mutated genes, representing 142 (9.2%) and 75 (4.8%) index cases of the whole series, respectively. KIF1A, ATL1, SPG11, KIF5A and REEP1 represented more than 1% (>17 cases) each. There were 661 causative variants (382 different ones) and 30 of them were structural variants. This large cohort allowed us to obtain an overview of the clinical and genetic spectrum of hereditary spastic paraplegia in clinical practice. Because of the wide phenotypic variability, there was no very specific sign that could predict the causative gene, but there were some constellations of symptoms that were found often related to specific subtypes. Finally, we confirmed the diagnostic effectiveness of a targeted sequencing panel as a first-line genetic test in hereditary spastic paraplegia. This is a pertinent strategy because of the relative frequency of several known genes (i.e. SPAST, KIF1A) and it allows identification of variants in the rarest involved genes and detection of structural rearrangements via coverage analysis, which is less efficient in exome datasets. It is crucial because these structural variants represent a significant proportion of the pathogenic hereditary spastic paraplegia variants (∼6% of patients), notably for SPAST and REEP1. In a subset of 42 index cases negative for the targeted multigene panel, subsequent whole-exome sequencing allowed a theoretical diagnosis yield of ∼50% to be reached. We then propose a two-step strategy combining the use of a panel of genes followed by whole-exome sequencing in negative cases.

Novel germline MET pathogenic variants in French patients with papillary renal cell carcinomas type I.

Hereditary papillary renal cell carcinoma (HPRC) is a rare inherited renal cancer syndrome characterized by bilateral and multifocal papillary type 1 renal tumors (PRCC1). Activating germline pathogenic variants of the MET gene were identified in HPRC families. We reviewed the medical and molecular records of a large French series of 158 patients screened for MET oncogenic variants. MET pathogenic and likely pathogenic variants rate was 12.4% with 40.6% among patients with familial PRCC1 and 5% among patients with sporadic PRCC1. The phenotype in cases with MET pathogenic and likely pathogenic variants was characteristic: PRCC1 tumors were mainly bilateral (84.3%) and multifocal (87.5%). Histologically, six out of seven patients with MET pathogenic variant harbored biphasic squamoid alveolar PRCC. Genetic screening identified one novel pathogenic variant MET c.3389T>C, p.(Leu1130Ser) and three novel likely pathogenic variants: MET c.3257A>T, p.(His1086Leu); MET c.3305T>C, p.(Ile1102Thr) and MET c.3373T>G, p.(Cys1125Gly). Functional assay confirmed their oncogenic effect as they induced an abnormal focus formation. The genotype-phenotype correlation between MET pathogenic variants and PRCC1 presentation should encourage to widen the screening, especially toward nonfamilial PRCC1. This precise phenotype also constitutes a strong argument for the classification of novel missense variants within the tyrosine kinase domain when functional assays are not accessible.

Phenoconversion from Spastic Paraplegia to ALS/FTD Associated with CYP7B1 Compound Heterozygous Mutations.

Biallelic mutations in the CYP7B1 gene lead to spastic paraplegia-5 (SPG5). We report herein the case of a patient whose clinical symptoms began with progressive lower limb spasticity during childhood, and who secondly developed amyotrophic lateral sclerosis/frontotemporal dementia (ALS/FTD) at the age of 67 years. Hereditary spastic paraplegia (HSP) gene analysis identified the compound heterozygous mutations c.825T>A (pTyr275*) and c.1193C>T (pPro398Leu) in CYP7B1 gene. No other pathogenic variant in frequent ALS/FTD causative genes was found. The CYP7B1 gene seems, therefore, to be the third gene associated with the phenoconversion from HSP to ALS, after the recently described UBQLN2 and ERLIN2 genes. We therefore expand the phenotype associated with CYP7B1 biallelic mutations and make an assumption about a link between cholesterol dyshomeostasis and ALS/FTD.

Classification of 101 BRCA1 and BRCA2 variants of uncertain significance by cosegregation study: A powerful approach.

Up to 80% of BRCA1 and BRCA2 genetic variants remain of uncertain clinical significance (VUSs). Only variants classified as pathogenic or likely pathogenic can guide breast and ovarian cancer prevention measures and treatment by PARP inhibitors. We report the first results of the ongoing French national COVAR (cosegregation variant) study, the aim of which is to classify BRCA1/2 VUSs. The classification method was a multifactorial model combining different associations between VUSs and cancer, including cosegregation data. At this time, among the 653 variants selected, 101 (15%) distinct variants shared by 1,624 families were classified as pathogenic/likely pathogenic or benign/likely benign by the COVAR study. Sixty-six of the 101 (65%) variants classified by COVAR would have remained VUSs without cosegregation data. Of note, among the 34 variants classified as pathogenic by COVAR, 16 remained VUSs or likely pathogenic when following the ACMG/AMP variant classification guidelines. Although the initiation and organization of cosegregation analyses require a considerable effort, the growing number of available genetic tests results in an increasing number of families sharing a particular variant, and thereby increases the power of such analyses. Here we demonstrate that variant cosegregation analyses are a powerful tool for the classification of variants in the BRCA1/2 breast-ovarian cancer predisposition genes.

Evidence of mosaicism in SPAST variant carriers in four French families.

Hereditary spastic paraplegias (HSP) are heterogeneous disorders, with more than 70 causative genes. Variants in SPAST are the most frequent genetic etiology and are responsible for spastic paraplegia type 4 (SPG4). Age at onset can vary, even between patients from the same family, and incomplete penetrance is described. Somatic mosaicism is extremely rare with only three patients reported in the literature. We report here SPAST mosaic variants in four unrelated patients. We confirm that mosaicism in SPAST is a very rare event with only four identified cases on more than 300 patients with a SPAST variant previously described by our clinical diagnostic laboratory.

Response to systemic therapy in fumarate hydratase-deficient renal cell carcinoma.

PURPOSE: Fumarate hydratase-deficient (FHdef) renal cell carcinoma (RCC) is a rare entity associated with the hereditary leiomyomatosis and RCC syndrome with no standard therapy approved. The aim of this retrospective study was to evaluate the efficacy of different systemic treatments in this population. METHODS: We performed a multicentre retrospective analysis of Fhdef RCC patients to determine the response to systemic treatments. The endpoints were objective response rate (ORR), time-to-treatment failure (TTF), and overall survival (OS). The two latter were estimated using the Kaplan-Meier method. RESULTS: Twenty-four Fhdef RCC patients were identified, and 21 under systemic therapy were included in the analysis: ten received cabozantinib, 14 received sunitinib, nine received "other antiangiogenics" (sorafenib, pazopanib, and axitinib), three received erlotinib-bevacizumab (E-B), three received mTOR inhibitors, and 11 received immune checkpoint blockers (ICBs). ORR for treatments were 50% for cabozantinib, 43% for sunitinib, 63% for "other antiangiogenics," and 30% for E-B, whereas ORR was 0% for mTOR inhibitors and 18% for ICBs. The median TTF (mTTF) was significantly higher with antiangiogenics (11.6 months) than with mTOR inhibitors (4.4 months) or ICBs (2.7 months). In the first-line setting, antiangiogenics presented a higher ORR compared with nivolumab-ipilimumab (64% versus 25%) and a significantly superior mTTF (11.0 months vs 2.5 months; p = 0.0027). The median OS from the start of the first systemic treatment was 44.0 months (95% confidence interval: 13.0-95.0). CONCLUSIONS: We report the first European retrospective study of Fhdef RCC patients treated with systemic therapy with a remarkably long median OS of 44.0 months. Our results suggest that antiangiogenics may be superior to ICB/mTOR inhibitors in this population.

Calibration of Pathogenicity Due to Variant-Induced Leaky Splicing Defects by Using BRCA2 Exon 3 as a Model System.

BRCA2 is a clinically actionable gene implicated in breast and ovarian cancer predisposition that has become a high priority target for improving the classification of variants of unknown significance (VUS). Among all BRCA2 VUS, those causing partial/leaky splicing defects are the most challenging to classify because the minimal level of full-length (FL) transcripts required for normal function remains to be established. Here, we explored BRCA2 exon 3 (BRCA2e3) as a model for calibrating variant-induced spliceogenicity and estimating thresholds for BRCA2 haploinsufficiency. In silico predictions, minigene splicing assays, patients' RNA analyses, a mouse embryonic stem cell (mESC) complementation assay and retrieval of patient-related information were combined to determine the minimal requirement of FL BRCA2 transcripts. Of 100 BRCA2e3 variants tested in the minigene assay, 64 were found to be spliceogenic, causing mild to severe RNA defects. Splicing defects were also confirmed in patients' RNA when available. Analysis of a neutral leaky variant (c.231T>G) showed that a reduction of approximately 60% of FL BRCA2 transcripts from a mutant allele does not cause any increase in cancer risk. Moreover, data obtained from mESCs suggest that variants causing a decline in FL BRCA2 with approximately 30% of wild-type are not pathogenic, given that mESCs are fully viable and resistant to DNA-damaging agents in those conditions. In contrast, mESCs producing lower relative amounts of FL BRCA2 exhibited either null or hypomorphic phenotypes. Overall, our findings are likely to have broader implications on the interpretation of BRCA2 variants affecting the splicing pattern of other essential exons. SIGNIFICANCE: These findings demonstrate that BRCA2 tumor suppressor function tolerates substantial reduction in full-length transcripts, helping to determine the pathogenicity of BRCA2 leaky splicing variants, some of which may not increase cancer risk.

Skipping Nonsense to Maintain Function: The Paradigm of BRCA2 Exon 12.

Germline nonsense and canonical splice site variants identified in disease-causing genes are generally considered as loss-of-function (LoF) alleles and classified as pathogenic. However, a fraction of such variants could maintain function through their impact on RNA splicing. To test this hypothesis, we used the alternatively spliced BRCA2 exon 12 (E12) as a model system because its in-frame skipping leads to a potentially functional protein. All E12 variants corresponding to putative LoF variants or predicted to alter splicing (n = 40) were selected from human variation databases and characterized for their impact on splicing in minigene assays and, when available, in patient lymphoblastoid cell lines. Moreover, a selection of variants was analyzed in a mouse embryonic stem cell-based functional assay. Using these complementary approaches, we demonstrate that a subset of variants, including nonsense variants, induced in-frame E12 skipping through the modification of splice sites or regulatory elements and, consequently, led to an internally deleted but partially functional protein. These data provide evidence, for the first time in a cancer-predisposition gene, that certain presumed null variants can retain function due to their impact on splicing. Further studies are required to estimate cancer risk associated with these hypomorphic variants. More generally, our findings highlight the need to exercise caution in the interpretation of putative LoF variants susceptible to induce in-frame splicing modifications. SIGNIFICANCE: This study presents evidence that certain presumed loss-of-function variants in a cancer predisposition gene can retain function due to their direct impact on RNA splicing.

Assessment of branch point prediction tools to predict physiological branch points and their alteration by variants.

BACKGROUND: Branch points (BPs) map within short motifs upstream of acceptor splice sites (3'ss) and are essential for splicing of pre-mature mRNA. Several BP-dedicated bioinformatics tools, including HSF, SVM-BPfinder, BPP, Branchpointer, LaBranchoR and RNABPS were developed during the last decade. Here, we evaluated their capability to detect the position of BPs, and also to predict the impact on splicing of variants occurring upstream of 3'ss. RESULTS: We used a large set of constitutive and alternative human 3'ss collected from Ensembl (n = 264,787 3'ss) and from in-house RNAseq experiments (n = 51,986 3'ss). We also gathered an unprecedented collection of functional splicing data for 120 variants (62 unpublished) occurring in BP areas of disease-causing genes. Branchpointer showed the best performance to detect the relevant BPs upstream of constitutive and alternative 3'ss (99.48 and 65.84% accuracies, respectively). For variants occurring in a BP area, BPP emerged as having the best performance to predict effects on mRNA splicing, with an accuracy of 89.17%. CONCLUSIONS: Our investigations revealed that Branchpointer was optimal to detect BPs upstream of 3'ss, and that BPP was most relevant to predict splicing alteration due to variants in the BP area.

SpliceLauncher: a tool for detection, annotation and relative quantification of alternative junctions from RNAseq data.

SUMMARY: Alternative splicing is an important biological process widely analyzed in molecular diagnostic settings. Indeed, a variant can be pathogenic by splicing alteration and a suspected pathogenic variant (e.g. truncating variant) can be rescued by splicing. In this context, detecting and quantifying alternative splicing is challenging. We developed SpliceLauncher, a fast and easy to use open source tool that aims at detecting, annotating and quantifying alternative splice junctions at high resolution. AVAILABILITY AND IMPLEMENTATION: SpliceLauncher is available at https://github.com/raphaelleman/SpliceLauncher. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Full in-frame exon 3 skipping of BRCA2 confers high risk of breast and/or ovarian cancer.

Germline pathogenic variants in the BRCA2 gene are associated with a cumulative high risk of breast/ovarian cancer. Several BRCA2 variants result in complete loss of the exon-3 at the transcript level. The pathogenicity of these variants and the functional impact of loss of exon 3 have yet to be established. As a collaboration of the COVAR clinical trial group (France), and the ENIGMA consortium for investigating breast cancer gene variants, this study evaluated 8 BRCA2 variants resulting in complete deletion of exon 3. Clinical information for 39 families was gathered from Portugal, France, Denmark and Sweden. Multifactorial likelihood analyses were conducted using information from 293 patients, for 7 out of the 8 variants (including 6 intronic). For all variants combined the likelihood ratio in favor of causality was 4.39*1025. These results provide convincing evidence for the pathogenicity of all examined variants that lead to a total exon 3 skipping, and suggest that other variants that result in complete loss of exon 3 at the molecular level could be associated with a high risk of cancer comparable to that associated with classical pathogenic variants in BRCA1 or BRCA2 gene. In addition, our functional study shows, for the first time, that deletion of exon 3 impairs the ability of cells to survive upon Mitomycin-C treatment, supporting lack of function for the altered BRCA2 protein in these cells. Finally, this study demonstrates that any variant leading to expression of only BRCA2 delta-exon 3 will be associated with an increased risk of breast and ovarian cancer.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort.

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5' and 3' consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Pattern multiplicity and fumarate hydratase (FH)/S-(2-succino)-cysteine (2SC) staining but not eosinophilic nucleoli with perinucleolar halos differentiate hereditary leiomyomatosis and renal cell carcinoma-associated renal cell carcinomas from kidney tumors without FH gene alteration.

Hereditary leiomyomatosis and renal cell carcinoma syndrome is characterized by an increased risk of agressive renal cell carcinoma, often of type 2 papillary histology, and is caused by FH germline mutations. A prominent eosinophilic macronucleolus with a perinucleolar clear halo is distinctive of hereditary leiomyomatosis and renal cell carcinoma syndrome-associated renal cell carcinoma according to the 2012 ISUP and 2016 WHO kidney tumor classification. From an immunohistochemistry perspective, tumors are often FH-negative and S-(2-succino)-cysteine (2SC) positive. We performed a pathology review of 24 renal tumors in 23 FH mutation carriers, and compared them to 12 type 2 papillary renal cell carcinomas from FH wild-type patients. Prominent eosinophilic nucleoli with perinucleolar halos were present in almost all FH-deficient renal cell carcinomas (23/24). Unexpectedly, they were also present in 58% of type 2 papillary renal cell carcinomas from wild-type patients. Renal cell carcinoma in mutation carriers displayed a complex architecture with multiple patterns, typically papillary, tubulopapillary, and tubulocystic, but also sarcomatoid and rhabdoid. Such pattern diversity was not seen in non-carriers. FH/2SC immunohistochemistry was informative as all hereditary leiomyomatosis and renal cell carcinoma-associated renal cell carcinomas were either FH- or 2SC+. For FH and 2SC immunohistochemistries taken separately, sensitivity of negative anti-FH immunohistochemistry was 87.5% and specificity was 100%. For positive anti-2SC immunohistochemistry, sensitivity, and specificity were 91.7% and 91.7%, respectively. All FH wild-type renal cell carcinoma were FH-positive, and all but one were 2SC-negative. In conclusion, multiplicity of architectural patterns, rhabdoid/sarcomatoid components and combined FH/2SC staining, but not prominent eosinophilic nucleoli with perinucleolar halos, differentiate hereditary leiomyomatosis and renal cell carcinoma-associated renal cell carcinoma from type 2 papillary renal cell carcinoma with efficient FH gene. Our findings are crucial in identifying who should be referred to Cancer Genetics clinics for genetic counseling and testing.

Germline SUFU mutation carriers and medulloblastoma: clinical characteristics, cancer risk, and prognosis.

Background: Germline mutations of suppressor of fused homolog (SUFU) predispose to sonic hedgehog (SHH) medulloblastoma. Germline SUFU mutations have been reported in nevoid basal cell carcinoma syndrome (NBCCS), but little is known about the cancer risk and clinical spectrum. Methods: We performed a retrospective review of all patients with medulloblastoma and a germline SUFU mutation in France. Results: Twenty-two patients from 17 families were identified with medulloblastoma and a germline SUFU mutation (median age at diagnosis: 16.5 mo). Macrocrania was present in 20 patients, but only 5 met the diagnostic criteria for NBCCS. Despite treatment with surgery and chemotherapy, to avoid radiotherapy in all patients except one, the outcome was worse than expected for SHH medulloblastoma, due to the high incidence of local relapses (8/22 patients) and second malignancies (n = 6 in 4/22 patients). The 5-year progression-free survival and overall survival rates were 42% and 66%. Mutations were inherited in 79% of patients, and 34 additional SUFU mutation carriers were identified within 14 families. Medulloblastoma penetrance was incomplete, but higher than in Patched 1 (PTCH1) mutation carriers. Besides medulloblastoma, 19 other tumors were recorded among the 56 SUFU mutation carriers, including basal cell carcinoma (BCC) in 2 patients and meningioma in 3 patients. Conclusion: Germline SUFU mutations strongly predispose to medulloblastoma in the first years of life, with worse prognosis than usually observed for SHH medulloblastoma. The clinical spectrum differs between SUFU and PTCH1 mutation carriers, and BCC incidence is much lower in SUFU mutation carriers. The optimal treatment of SUFU mutation-associated medulloblastoma has not been defined.

Detecting splicing patterns in genes involved in hereditary breast and ovarian cancer.

Interpretation of variants of unknown significance (VUS) is a major challenge for laboratories performing molecular diagnosis of hereditary breast and ovarian cancer (HBOC), especially considering that many genes are now known to be involved in this syndrome. One important way these VUS can have a functional impact is through their effects on RNA splicing. Here we present a custom RNA-Seq assay plus bioinformatics and biostatistics pipeline to analyse specifically alternative and abnormal splicing junctions in 11 targeted HBOC genes. Our pipeline identified 14 new alternative splices in BRCA1 and BRCA2 in addition to detecting the majority of known alternative spliced transcripts therein. We provide here the first global splicing pattern analysis for the other nine genes, which will enable a comprehensive interpretation of splicing defects caused by VUS in HBOC. Previously known splicing alterations were consistently detected, occasionally with a more complex splicing pattern than expected. We also found that splicing in the 11 genes is similar in blood and breast tissue, supporting the utility and simplicity of blood splicing assays. Our pipeline is ready to be integrated into standard molecular diagnosis for HBOC, but it could equally be adapted for an integrative analysis of any multigene disorder.