The BRCA2 R2645G variant increases DNA binding and induces hyper-recombination.

BRCA2 tumor suppressor protein ensures genome integrity by mediating DNA repair via homologous recombination (HR). This function is executed in part by its canonical DNA binding domain located at the C-terminus (BRCA2CTD), the only folded domain of the protein. Most germline pathogenic missense variants are located in this highly conserved region which binds to single-stranded DNA (ssDNA) and to the acidic protein DSS1. These interactions are essential for the HR function of BRCA2. Here, we report that the variant R2645G, identified in breast cancer and located at the DSS1 interface, unexpectedly increases the ssDNA binding activity of BRCA2CTDin vitro. Human cells expressing this variant display a hyper-recombination phenotype, chromosomal instability in the form of chromatid gaps when exposed to DNA damage, and increased PARP inhibitor sensitivity. In mouse embryonic stem cells (mES), this variant alters viability and confers sensitivity to cisplatin and Mitomycin C. These results suggest that BRCA2 interaction with ssDNA needs to be tightly regulated to limit HR and prevent chromosomal instability and we propose that this control mechanism involves DSS1. Given that several missense variants located within this region have been identified in breast cancer patients, these findings might have clinical implications for carriers.

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.

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.

Germline HPF1 retrogene insertion in RB1 gene involved in cancer predisposition.

About half of the human genome is composed of repeated sequences derived from mobile elements, mainly retrotransposons, generally without pathogenic effect. Familial forms of retinoblastoma are caused by germline pathogenic variants in RB1 gene. Here, we describe a family with retinoblastoma affecting a father and his son. No pathogenic variant was identified after DNA analysis of RB1 gene coding sequence and exon-intron junctions. However, RB1 mRNA analysis showed a chimeric transcript with insertion of 114 nucleotides from HPF1 gene inside RB1 gene. This chimeric transcript led to an insertion of 38 amino acids in functional domain of retinoblastoma protein. Subsequent DNA analysis in RB1 intron 17 revealed the presence of a full-length HPF1 retrogene insertion in opposite orientation. Functional assay shows that this insertion has a deleterious impact on retinoblastoma protein function. This is the first report of a full-length retrogene insertion involved in human Mendelian disease leading to a chimeric transcript and a non-functional chimeric protein. Some retrogene insertions may be missed by standard diagnostic genetic testing, so contribution of retrogene insertions to human disease may be underestimated. The increasing use of whole genome sequencing in diagnostic settings will help to get a more comprehensive view of retrogenes.

Adaptive nanopore sequencing to determine pathogenicity of BRCA1 exonic duplication.

BRCA1 and BRCA2 are tumour suppressor genes that have been characterised as predisposition genes for the development of hereditary breast and ovarian cancers among other malignancies. The molecular diagnosis of this predisposition syndrome is based on the detection of inactivating variants of any type in those genes. But in the case of structural variants, functional consequences can be difficult to assess using standard molecular methods, as the precise resolution of their sequence is often impossible with short-read next generation sequencing techniques. It has been recently demonstrated that Oxford Nanopore long-read sequencing technology can accurately and rapidly provide genetic diagnoses of Mendelian diseases, including those linked to pathogenic structural variants. Here, we report the accurate resolution of a germline duplication event of exons 18-20 of BRCA1 using Nanopore sequencing with adaptive sampling target enrichment. This allowed us to classify this variant as pathogenic within a short timeframe of 10 days. This study provides a proof-of-concept that nanopore adaptive sampling is a highly efficient technique for the investigation of structural variants of tumour suppressor genes in a clinical context.

A likelihood ratio approach for utilizing case-control data in the clinical classification of rare sequence variants: application to BRCA1 and BRCA2.

A large number of variants identified through clinical genetic testing in disease susceptibility genes, are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion), can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analyses of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC), and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity - findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared to classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and pre-formatted excel calculators for implementation of the method for rare variants in BRCA1, BRCA2 and other high-risk genes with known penetrance.

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.

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.

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

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.

A new hybrid record linkage process to make epidemiological databases interoperable: application to the GEMO and GENEPSO studies involving BRCA1 and BRCA2 mutation carriers.

BACKGROUND: Linking independent sources of data describing the same individuals enable innovative epidemiological and health studies but require a robust record linkage approach. We describe a hybrid record linkage process to link databases from two independent ongoing French national studies, GEMO (Genetic Modifiers of BRCA1 and BRCA2), which focuses on the identification of genetic factors modifying cancer risk of BRCA1 and BRCA2 mutation carriers, and GENEPSO (prospective cohort of BRCAx mutation carriers), which focuses on environmental and lifestyle risk factors. METHODS: To identify as many as possible of the individuals participating in the two studies but not registered by a shared identifier, we combined probabilistic record linkage (PRL) and supervised machine learning (ML). This approach (named "PRL + ML") combined together the candidate matches identified by both approaches. We built the ML model using the gold standard on a first version of the two databases as a training dataset. This gold standard was obtained from PRL-derived matches verified by an exhaustive manual review. Results The Random Forest (RF) algorithm showed a highest recall (0.985) among six widely used ML algorithms: RF, Bagged trees, AdaBoost, Support Vector Machine, Neural Network. Therefore, RF was selected to build the ML model since our goal was to identify the maximum number of true matches. Our combined linkage PRL + ML showed a higher recall (range 0.988-0.992) than either PRL (range 0.916-0.991) or ML (0.981) alone. It identified 1995 individuals participating in both GEMO (6375 participants) and GENEPSO (4925 participants). CONCLUSIONS: Our hybrid linkage process represents an efficient tool for linking GEMO and GENEPSO. It may be generalizable to other epidemiological studies involving other databases and registries.

Is BRCA2 involved in early onset colorectal cancer risk?

A, Closed symbols indicate patients affected with cancer. Open symbols indicate healthy individuals. The type of cancer and age at presentation are given in brackets. Blue circle represents c.4471_4474del variant and red circle represents the c.9648 + 1G > A. B, RNA was extracted from blood of patient III-3 and his sisters III-1 and III-4. RT-PCR analysis was performed with primers mapping to exons 25 and 27, and PCR products were separated by Bioanalyzer electrophoresis. The sizes of the DNA marker (M) are indicated to the left. LM, lower marker; UM, upper marker. C, Each RT-PCR product from patient III-3 was gel-purified and analyzed by Sanger sequencing. The 297-bp band corresponds to the reference BRCA2 transcript and the 150-bp band corresponds to a BRCA2 transcript lacking exon 26.

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

HRness in Breast and Ovarian Cancers.

Ovarian and breast cancers are currently defined by the main pathways involved in the tumorigenesis. The majority are carcinomas, originating from epithelial cells that are in constant division and subjected to cyclical variations of the estrogen stimulus during the female hormonal cycle, therefore being vulnerable to DNA damage. A portion of breast and ovarian carcinomas arises in the context of DNA repair defects, in which genetic instability is the backdrop for cancer initiation and progression. For these tumors, DNA repair deficiency is now increasingly recognized as a target for therapeutics. In hereditary breast/ovarian cancers (HBOC), tumors with BRCA1/2 mutations present an impairment of DNA repair by homologous recombination (HR). For many years, BRCA1/2 mutations were only screened on germline DNA, but now they are also searched at the tumor level to personalize treatment. The reason of the inactivation of this pathway remains uncertain for most cases, even in the presence of a HR-deficient signature. Evidence indicates that identifying the mechanism of HR inactivation should improve both genetic counseling and therapeutic response, since they can be useful as new biomarkers of response.

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.

The BRCA1 c. 5096G>A p.Arg1699Gln (R1699Q) intermediate risk variant: breast and ovarian cancer risk estimation and recommendations for clinical management from the ENIGMA consortium.

BACKGROUND: We previously showed that the BRCA1 variant c.5096G>A p.Arg1699Gln (R1699Q) was associated with an intermediate risk of breast cancer (BC) and ovarian cancer (OC). This study aimed to assess these cancer risks for R1699Q carriers in a larger cohort, including follow-up of previously studied families, to further define cancer risks and to propose adjusted clinical management of female BRCA1*R1699Q carriers. METHODS: Data were collected from 129 BRCA1*R1699Q families ascertained internationally by ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) consortium members. A modified segregation analysis was used to calculate BC and OC risks. Relative risks were calculated under both monogenic model and major gene plus polygenic model assumptions. RESULTS: In this cohort the cumulative risk of BC and OC by age 70 years was 20% and 6%, respectively. The relative risk for developing cancer was higher when using a model that included the effects of both the R1699Q variant and a residual polygenic component compared with monogenic model (for BC 3.67 vs 2.83, and for OC 6.41 vs 5.83). CONCLUSION: Our results confirm that BRCA1*R1699Q confers an intermediate risk for BC and OC. Breast surveillance for female carriers based on mammogram annually from age 40 is advised. Bilateral salpingo-oophorectomy should be considered based on family history.

BRCA1 and BRCA2 5' noncoding region variants identified in breast cancer patients alter promoter activity and protein binding.

The widespread use of next generation sequencing for clinical testing is detecting an escalating number of variants in noncoding regions of the genome. The clinical significance of the majority of these variants is currently unknown, which presents a significant clinical challenge. We have screened over 6,000 early-onset and/or familial breast cancer (BC) cases collected by the ENIGMA consortium for sequence variants in the 5' noncoding regions of BC susceptibility genes BRCA1 and BRCA2, and identified 141 rare variants with global minor allele frequency < 0.01, 76 of which have not been reported previously. Bioinformatic analysis identified a set of 21 variants most likely to impact transcriptional regulation, and luciferase reporter assays detected altered promoter activity for four of these variants. Electrophoretic mobility shift assays demonstrated that three of these altered the binding of proteins to the respective BRCA1 or BRCA2 promoter regions, including NFYA binding to BRCA1:c.-287C>T and PAX5 binding to BRCA2:c.-296C>T. Clinical classification of variants affecting promoter activity, using existing prediction models, found no evidence to suggest that these variants confer a high risk of disease. Further studies are required to determine if such variation may be associated with a moderate or low risk of BC.

Combined genetic and splicing analysis of BRCA1 c.[594-2A>C; 641A>G] highlights the relevance of naturally occurring in-frame transcripts for developing disease gene variant classification algorithms.

A recent analysis using family history weighting and co-observation classification modeling indicated that BRCA1 c.594-2A > C (IVS9-2A > C), previously described to cause exon 10 skipping (a truncating alteration), displays characteristics inconsistent with those of a high risk pathogenic BRCA1 variant. We used large-scale genetic and clinical resources from the ENIGMA, CIMBA and BCAC consortia to assess pathogenicity of c.594-2A > C. The combined odds for causality considering case-control, segregation and breast tumor pathology information was 3.23 × 10-8 Our data indicate that c.594-2A > C is always in cis with c.641A > G. The spliceogenic effect of c.[594-2A > C;641A > G] was characterized using RNA analysis of human samples and splicing minigenes. As expected, c.[594-2A > C; 641A > G] caused exon 10 skipping, albeit not due to c.594-2A > C impairing the acceptor site but rather by c.641A > G modifying exon 10 splicing regulatory element(s). Multiple blood-based RNA assays indicated that the variant allele did not produce detectable levels of full-length transcripts, with a per allele BRCA1 expression profile composed of ≈70-80% truncating transcripts, and ≈20-30% of in-frame Δ9,10 transcripts predicted to encode a BRCA1 protein with tumor suppression function.We confirm that BRCA1c.[594-2A > C;641A > G] should not be considered a high-risk pathogenic variant. Importantly, results from our detailed mRNA analysis suggest that BRCA-associated cancer risk is likely not markedly increased for individuals who carry a truncating variant in BRCA1 exons 9 or 10, or any other BRCA1 allele that permits 20-30% of tumor suppressor function. More generally, our findings highlight the importance of assessing naturally occurring alternative splicing for clinical evaluation of variants in disease-causing genes.