Mosaic BRCA1 promoter methylation contribution in hereditary breast/ovarian cancer pedigrees.

PURPOSE: Mosaic BRCA1 promoter methylation (BRCA1meth) increases the risk of early-onset breast cancer, triple-negative breast cancer and ovarian cancer. As mosaic BRCA1meth are believed to occur de novo, their role in family breast/ovarian cancer has not been assessed. PATIENTS: Blood-derived DNA from 20 unrelated affected cases from families with aggregation of breast/ovarian cancer, but with no germline pathogenic variants in BRCA1/2, PALB2 or RAD51C/D, were screened by methylation-sensitive high-resolution melting. CpG analysis was performed by pyrosequencing on blood and buccal swab. Two probands carried a pathogenic variant in a moderate-penetrance gene (ATM and BARD1), and 8 of 18 others (44%) carried BRCA1meth (vs none of the 20 age-matched controls). Involvement of BRCA1 in tumourigenesis in methylated probands was demonstrated in most tested cases by detection of a loss of heterozygosity and a homologous recombination deficiency signature. Among the eight methylated probands, two had relatives with breast cancer with detectable BRCA1meth in blood, including one with high methylation levels in two non-tumour tissues. CONCLUSIONS: The high prevalence of mosaic BRCA1meth in patients with breast/ovarian cancer with affected relatives, as well as this first description of a family aggregation of mosaic BRCA1meth, shows how this de novo event can contribute to hereditary breast/ovarian cancer pedigrees.

APC germline pathogenic variants and epithelial ovarian cancer: causal or coincidental findings?

APC germline pathogenic variants result in predisposition to familial adenomatous polyposis and extraintestinal tumours such as desmoid fibromatosis, medulloblastomas and thyroid cancers. They have also been recently involved in ovarian microcystic stromal tumours. APC inactivation has been described at the tumour level in epithelial ovarian cancers (EOCs). Here, we report the identification of APC germline pathogenic variants in two patients diagnosed with premenopausal EOC in early 30s, with no other pathogenic variant detected in the known ovarian cancer predisposing genes. Subsequent tumour analysis showed neither a second hit of APC inactivation nor ß-catenin activation. Both tumours did not have a homologous recombination (HR) deficiency, pointing towards the implication of other genes than those involved in HR. APC may contribute to the carcinogenesis of EOC in a multifactorial context. Further studies are required to clarify the role of APC in predisposition to EOC.

Hereditary cancer predispositions: Comparison of multigene panel sequencing on fresh-frozen breast/ovarian tumor versus blood.

In breast or ovarian cancer (BC/OC) patients with evocative personal and/or family history, multigene panel sequencing is performed on blood to diagnose hereditary predispositions. Additionally, BRCA1/BRCA2 testing can be performed on tumor sample for therapeutic purpose. The accuracy of multigene panel tumor analysis on BC/OC to detect predisposing germline pathogenic variants (gPV) has not been precisely assessed. By comparing sequencing data from blood and fresh-frozen tumor we show that tumor genomic instability causes pitfalls to consider when performing tumor testing to detect gPV. Even if loss of heterozygosity increases germline signal in most cases, somatic copy number variants (CNV) can mask germline CNV and collapse point gPV variant allele frequency (VAF). Moreover, VAF does not allow an accurate distinction between germline and somatic pathogenic variants.

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/.

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.

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/).

Functional classification of ATM variants in ataxia-telangiectasia patients.

Ataxia-telangiectasia (A-T) is a recessive disorder caused by biallelic pathogenic variants of ataxia-telangiectasia mutated (ATM). This disease is characterized by progressive ataxia, telangiectasia, immune deficiency, predisposition to malignancies, and radiosensitivity. However, hypomorphic variants may be discovered associated with very atypical phenotypes, raising the importance of evaluating their pathogenic effects. In this study, multiple functional analyses were performed on lymphoblastoid cell lines from 36 patients, comprising 49 ATM variants, 24 being of uncertain significance. Thirteen patients with atypical phenotype and presumably hypomorphic variants were of particular interest to test strength of functional analyses and to highlight discrepancies with typical patients. Western-blot combined with transcript analyses allowed the identification of one missing variant, confirmed suspected splice defects and revealed unsuspected minor transcripts. Subcellular localization analyses confirmed the low level and abnormal cytoplasmic localization of ATM for most A-T cell lines. Interestingly, atypical patients had lower kinase defect and less altered cell-cycle distribution after genotoxic stress than typical patients. In conclusion, this study demonstrated the pathogenic effects of the 49 variants, highlighted the strength of KAP1 phosphorylation test for pathogenicity assessment and allowed the establishment of the Ataxia-TeLangiectasia Atypical Score to predict atypical phenotype. Altogether, we propose strategies for ATM variant detection and classification.

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.

Genomic hallmarks of homologous recombination deficiency in invasive breast carcinomas.

Therapeutic strategies targeting Homologous Recombination Deficiency (HRD) in breast cancer requires patient stratification. The LST (Large-scale State Transitions) genomic signature previously validated for triple-negative breast carcinomas (TNBC) was evaluated as biomarker of HRD in luminal (hormone receptor positive) and HER2-overexpressing (HER2+) tumors. The LST genomic signature related to the number of large-scale chromosomal breakpoints in SNP-array tumor profile was applied to identify HRD in in-house and TCGA sets of breast tumors, in which the status of BRCA1/2 and other genes was also investigated. In the in-house dataset, HRD was predicted in 5% (20/385) of sporadic tumors luminal or HER2+ by the LST genomic signature and the inactivation of BRCA1, BRCA2 or RAD51C confirmed this prediction in 75% (12/16) of the tested cases. In 14% (6/43) of tumors occurring in BRCA1/2 mutant carriers, the corresponding wild-type allele was retained emphasizing the importance of determining the tumor status. In the TCGA luminal and HER2+ subtypes HRD incidence was estimated at 5% (18/329, 95%CI: 5-8%) and 2% (1/59, 95%CI: 2-9%), respectively. In TNBC cisplatin-based neo-adjuvant clinical trials, HRD is shown to be a necessary condition for cisplatin sensitivity. This analysis demonstrates the high performance of the LST genomic signature for HRD detection in breast cancers, which suggests its potential as a biomarker for genetic testing and patient stratification for clinical trials evaluating platinum salts and PARP inhibitors.

Germline BAP1 mutations predispose to renal cell carcinomas.

The genetic cause of some familial nonsyndromic renal cell carcinomas (RCC) defined by at least two affected first-degree relatives is unknown. By combining whole-exome sequencing and tumor profiling in a family prone to cases of RCC, we identified a germline BAP1 mutation c.277A>G (p.Thr93Ala) as the probable genetic basis of RCC predisposition. This mutation segregated with all four RCC-affected relatives. Furthermore, BAP1 was found to be inactivated in RCC-affected individuals from this family. No BAP1 mutations were identified in 32 familial cases presenting with only RCC. We then screened for germline BAP1 deleterious mutations in familial aggregations of cancers within the spectrum of the recently described BAP1-associated tumor predisposition syndrome, including uveal melanoma, malignant pleural mesothelioma, and cutaneous melanoma. Among the 11 families that included individuals identified as carrying germline deleterious BAP1 mutations, 6 families presented with 9 RCC-affected individuals, demonstrating a significantly increased risk for RCC. This strongly argues that RCC belongs to the BAP1 syndrome and that BAP1 is a RCC-predisposition gene.

Male breast cancer: No evidence for mosaic BRCA1 promoter methylation involvement.

Breast cancers (BC) are rare in men and are often caused by constitutional predisposing factors. In women, mosaic BRCA1 promoter methylations (MBPM) are frequent events, detected in 4-8% of healthy subjects. This constitutional epimutation increases risk of early-onset and triple-negative BC. However, the role of MBPM in male BC predisposition has never been assessed. We screened 40 blood samples from men affected by BC, and performed extensive tumour analysis on MBPM-positive patients. We detected two patients carrying MBPM. Surprisingly, tumour analysis revealed that neither of these two male BCs were caused by the constitutional BRCA1 epimutations carried by the patients.

Loss of heterozygosity at 13q13 and 14q32 predicts BRCA2 inactivation in luminal breast carcinomas.

BRCA2 is the major high-penetrance predisposition gene for luminal (estrogen receptor [ER] positive) breast cancers. However, many BRCA2 mutant carriers lack family history of breast/ovarian cancers and do not benefit from genetic testing. Specific genomic features associated with BRCA2 inactivation in tumors could help identify patients for whom a genetic test for BRCA2 may be proposed. A series of ER-positive invasive ductal carcinomas (IDCs) including 30 carriers of BRCA2 mutations and 215 control cases was studied by single-nucleotide polymorphism (SNP) arrays. Cases and controls were stratified by grade and HER2 status. Independently, 7 BRCA2 and 51 control cases were used for validation. Absolute copy number and Loss of heterozygosity (LOH) profiles were obtained from SNP arrays by the genome alteration print (GAP) method. BRCA2 tumors were observed to display a discriminatively greater number of chromosomal breaks calculated after filtering out and smoothing <3 Mb variations. This argues for a BRCA2-associated genomic instability responsible for long-segment aberrations. Co-occurrence of two genomic features-LOH of 13q13 and 14q32-was found to predict BRCA2 status with 90% of sensitivity and 87% of specificity in discovery series of high-grade HER2-negative IDCs and 100% of sensitivity and 88% of specificity in an independent series of 58 IDCs. Estimated positive predictive value was 17.2% (confidence interval: 6.7-33.5) in the whole series. In conclusion, the simplified BRCA2 classifier based on the co-occurrence of LOH at 13q13 and 14q32 could provide an indication to test for BRCA2 mutation in patients with ER-positive IDC.

Common variants of the BRCA1 wild-type allele modify the risk of breast cancer in BRCA1 mutation carriers.

Mutations in the BRCA1 gene substantially increase a woman's lifetime risk of breast cancer. However, there is great variation in this increase in risk with several genetic and non-genetic modifiers identified. The BRCA1 protein plays a central role in DNA repair, a mechanism that is particularly instrumental in safeguarding cells against tumorigenesis. We hypothesized that polymorphisms that alter the expression and/or function of BRCA1 carried on the wild-type (non-mutated) copy of the BRCA1 gene would modify the risk of breast cancer in carriers of BRCA1 mutations. A total of 9874 BRCA1 mutation carriers were available in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) for haplotype analyses of BRCA1. Women carrying the rare allele of single nucleotide polymorphism rs16942 on the wild-type copy of BRCA1 were at decreased risk of breast cancer (hazard ratio 0.86, 95% confidence interval 0.77-0.95, P = 0.003). Promoter in vitro assays of the major BRCA1 haplotypes showed that common polymorphisms in the regulatory region alter its activity and that this effect may be attributed to the differential binding affinity of nuclear proteins. In conclusion, variants on the wild-type copy of BRCA1 modify risk of breast cancer among carriers of BRCA1 mutations, possibly by altering the efficiency of BRCA1 transcription.

Germline mutation in the RAD51B gene confers predisposition to breast cancer.

BACKGROUND: Most currently known breast cancer predisposition genes play a role in DNA repair by homologous recombination. Recent studies conducted on RAD51 paralogs, involved in the same DNA repair pathway, have identified rare germline mutations conferring breast and/or ovarian cancer predisposition in the RAD51C, RAD51D and XRCC2 genes. The present study analysed the five RAD51 paralogs (RAD51B, RAD51C, RAD51D, XRCC2, XRCC3) to estimate their contribution to breast and ovarian cancer predisposition. METHODS: The study was conducted on 142 unrelated patients with breast and/or ovarian cancer either with early onset or with a breast/ovarian cancer family history. Patients were referred to a French family cancer clinic and had been previously tested negative for a BRCA1/2 mutation. Coding sequences of the five genes were analysed by EMMA (Enhanced Mismatch Mutation Analysis). Detected variants were characterized by Sanger sequencing analysis. RESULTS: Three splicing mutations and two likely deleterious missense variants were identified: RAD51B c.452 + 3A > G, RAD51C c.706-2A > G, RAD51C c.1026 + 5_1026 + 7del, RAD51B c.475C > T/p.Arg159Cys and XRCC3 c.448C > T/p.Arg150Cys. No RAD51D and XRCC2 gene mutations were detected. These mutations and variants were detected in families with both breast and ovarian cancers, except for the RAD51B c.475C > T/p.Arg159Cys variant that occurred in a family with 3 breast cancer cases. CONCLUSIONS: This study identified the first RAD51B mutation in a breast and ovarian cancer family and is the first report of XRCC3 mutation analysis in breast and ovarian cancer. It confirms that RAD51 paralog mutations confer breast and ovarian cancer predisposition and are rare events. In view of the low frequency of RAD51 paralog mutations, international collaboration of family cancer clinics will be required to more accurately estimate their penetrance and establish clinical guidelines in carrier individuals.

Guidelines for splicing analysis in molecular diagnosis derived from a set of 327 combined in silico/in vitro studies on BRCA1 and BRCA2 variants.

Assessing the impact of variants of unknown significance (VUS) on splicing is a key issue in molecular diagnosis. This impact can be predicted by in silico tools, but proper evaluation and user guidelines are lacking. To fill this gap, we embarked upon the largest BRCA1 and BRCA2 splice study to date by testing 272 VUSs (327 analyses) within the BRCA splice network of Unicancer. All these VUSs were analyzed by using six tools (splice site prediction by neural network, splice site finder (SSF), MaxEntScan (MES), ESE finder, relative enhancer and silencer classification by unanimous enrichment, and human splicing finder) and the predictions obtained were compared with transcript analysis results. Combining MES and SSF gave 96% sensitivity and 83% specificity for VUSs occurring in the vicinity of consensus splice sites, that is, the surrounding 11 and 14 bases for the 5' and 3' sites, respectively. This study was also an opportunity to define guidelines for transcript analysis along with a tentative classification of splice variants. The guidelines drawn from this large series should be useful for the whole community, particularly in the context of growing sequencing capacities that require robust pipelines for variant interpretation.

EMMA, a cost- and time-effective diagnostic method for simultaneous detection of point mutations and large-scale genomic rearrangements: application to BRCA1 and BRCA2 in 1,525 patients.

The detection of unknown mutations remains a serious challenge and, despite the expected benefits for the patient's health, a large number of genes are not screened on a routine basis. We present the diagnostic application of EMMA (Enhanced Mismatch Mutation Analysis(®) , Fluigent, Paris, France), a novel method based on heteroduplex analysis by capillary electrophoresis using innovative matrices. BRCA1 and BRCA2 were screened for point mutations and large rearrangements in 1,525 unrelated patients (372 for the validation step and 1,153 in routine diagnosis) using a single analytical condition. Seven working days were needed for complete BRCA1/2 screening in 30 patients by one technician (excluding DNA extraction and sequencing). A total of 137 mutations were found, including a BRCA2 duplication of exons 19 and 20, previously missed by Comprehensive BRACAnalysis(®) . The mutation detection rate was 11.9%, which is consistent with patient inclusions. This study therefore suggests that EMMA represents a valuable short-term and midterm option for many diagnostic laboratories looking for an easy, reliable, and affordable strategy, enabling fast and sensitive analysis for a large number of genes.

5' Region Large Genomic Rearrangements in the BRCA1 Gene in French Families: Identification of a Tandem Triplication and Nine Distinct Deletions with Five Recurrent Breakpoints.

BACKGROUND: Large genomic rearrangements (LGR) in BRCA1 consisting of deletions/duplications of one or several exons have been found throughout the gene with a large proportion occurring in the 5' region from the promoter to exon 2. The aim of this study was to better characterize those LGR in French high-risk breast/ovarian cancer families. METHODS: DNA from 20 families with one apparent duplication and nine deletions was analyzed with a dedicated comparative genomic hybridization (CGH) array, high-resolution BRCA1 Genomic Morse Codes analysis and Sanger sequencing. RESULTS: The apparent duplication was in fact a tandem triplication of exons 1 and 2 and part of intron 2 of BRCA1, fully characterized here for the first time. We calculated a causality score with the multifactorial model from data obtained from six families, classifying this variant as benign. Among the nine deletions detected in this region, eight have never been identified. The breakpoints fell in six recurrent regions and could confirm some specific conformation of the chromatin. CONCLUSIONS: Taken together, our results firmly establish that the BRCA1 5' region is a frequent site of different LGRs and highlight the importance of the segmental duplication and Alu sequences, particularly the very high homologous region, in the mechanism of a recombination event. This also confirmed that those events are not systematically deleterious.

The Spectrum of FANCM Protein Truncating Variants in European Breast Cancer Cases.

Germline protein truncating variants (PTVs) in the FANCM gene have been associated with a 2-4-fold increased breast cancer risk in case-control studies conducted in different European populations. However, the distribution and the frequency of FANCM PTVs in Europe have never been investigated. In the present study, we collected the data of 114 European female breast cancer cases with FANCM PTVs ascertained in 20 centers from 13 European countries. We identified 27 different FANCM PTVs. The p.Gln1701* PTV is the most common PTV in Northern Europe with a maximum frequency in Finland and a lower relative frequency in Southern Europe. On the contrary, p.Arg1931* seems to be the most common PTV in Southern Europe. We also showed that p.Arg658*, the third most common PTV, is more frequent in Central Europe, and p.Gln498Thrfs*7 is probably a founder variant from Lithuania. Of the 23 rare or unique FANCM PTVs, 15 have not been previously reported. We provide here the initial spectrum of FANCM PTVs in European breast cancer cases.