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.

Computational Tools for Splicing Defect Prediction in Breast/Ovarian Cancer Genes: How Efficient Are They at Predicting RNA Alterations?

In silico tools for splicing defect prediction have a key role to assess the impact of variants of uncertain significance. Our aim was to evaluate the performance of a set of commonly used splicing in silico tools comparing the predictions against RNA in vitro results. This was done for natural splice sites of clinically relevant genes in hereditary breast/ovarian cancer (HBOC) and Lynch syndrome. A study divided into two stages was used to evaluate SSF-like, MaxEntScan, NNSplice, HSF, SPANR, and dbscSNV tools. A discovery dataset of 99 variants with unequivocal results of RNA in vitro studies, located in the 10 exonic and 20 intronic nucleotides adjacent to exon-intron boundaries of BRCA1, BRCA2, MLH1, MSH2, MSH6, PMS2, ATM, BRIP1, CDH1, PALB2, PTEN, RAD51D, STK11, and TP53, was collected from four Spanish cancer genetic laboratories. The best stand-alone predictors or combinations were validated with a set of 346 variants in the same genes with clear splicing outcomes reported in the literature. Sensitivity, specificity, accuracy, negative predictive value (NPV) and Mathews Coefficient Correlation (MCC) scores were used to measure the performance. The discovery stage showed that HSF and SSF-like were the most accurate for variants at the donor and acceptor region, respectively. The further combination analysis revealed that HSF, HSF+SSF-like or HSF+SSF-like+MES achieved a high performance for predicting the disruption of donor sites, and SSF-like or a sequential combination of MES and SSF-like for predicting disruption of acceptor sites. The performance confirmation of these last results with the validation dataset, indicated that the highest sensitivity, accuracy, and NPV (99.44%, 99.44%, and 96.88, respectively) were attained with HSF+SSF-like or HSF+SSF-like+MES for donor sites and SSF-like (92.63%, 92.65%, and 84.44, respectively) for acceptor sites. We provide recommendations for combining algorithms to conduct in silico splicing analysis that achieved a high performance. The high NPV obtained allows to select the variants in which the study by in vitro RNA analysis is mandatory against those with a negligible probability of being spliceogenic. Our study also shows that the performance of each specific predictor varies depending on whether the natural splicing sites are donors or acceptors.

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.

Thorough in silico and in vitro cDNA analysis of 21 putative BRCA1 and BRCA2 splice variants and a complex tandem duplication in BRCA2 allowing the identification of activated cryptic splice donor sites in BRCA2 exon 11.

For 21 putative BRCA1 and BRCA2 splice site variants, the concordance between mRNA analysis and predictions by in silico programs was evaluated. Aberrant splicing was confirmed for 12 alterations. In silico prediction tools were helpful to determine for which variants cDNA analysis is warranted, however, predictions for variants in the Cartegni consensus region but outside the canonical sites, were less reliable. Learning algorithms like Adaboost and Random Forest outperformed the classical tools. Further validations are warranted prior to implementation of these novel tools in clinical settings. Additionally, we report here for the first time activated cryptic donor sites in the large exon 11 of BRCA2 by evaluating the effect at the cDNA level of a novel tandem duplication (5' breakpoint in intron 4; 3' breakpoint in exon 11) and of a variant disrupting the splice donor site of exon 11 (c.6841+1G > C). Additional sites were predicted, but not activated. These sites warrant further research to increase our knowledge on cis and trans acting factors involved in the conservation of correct transcription of this large exon. This may contribute to adequate design of ASOs (antisense oligonucleotides), an emerging therapy to render cancer cells sensitive to PARP inhibitor and platinum therapies.

Naturally occurring BRCA2 alternative mRNA splicing events in clinically relevant samples.

BACKGROUND: BRCA1 and BRCA2 are the two principal tumour suppressor genes associated with inherited high risk of breast and ovarian cancer. Genetic testing of BRCA1/2 will often reveal one or more sequence variants of uncertain clinical significance, some of which may affect normal splicing patterns and thereby disrupt gene function. mRNA analyses are therefore among the tests used to interpret the clinical significance of some genetic variants. However, these could be confounded by the appearance of naturally occurring alternative transcripts unrelated to germline sequence variation or defects in gene function. To understand which novel splicing events are associated with splicing mutations and which are part of the normal BRCA2 splicing repertoire, a study was undertaken by members of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium to characterise the spectrum of naturally occurring BRCA2 mRNA alternate-splicing events. METHODS: mRNA was prepared from several blood and breast tissue-derived cells and cell lines by contributing ENIGMA laboratories. cDNA representing BRCA2 alternate splice sites was amplified and visualised using capillary or agarose gel electrophoresis, followed by sequencing. RESULTS: We demonstrate the existence of 24 different BRCA2 mRNA alternate-splicing events in lymphoblastoid cell lines and both breast cancer and non-cancerous breast cell lines. CONCLUSIONS: These naturally occurring alternate-splicing events contribute to the array of cDNA fragments that may be seen in assays for mutation-associated splicing defects. Caution must be observed in assigning alternate-splicing events to potential splicing mutations.

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.

RAD51C germline mutations found in Spanish site-specific breast cancer and breast-ovarian cancer families.

BRCA1 and BRCA2 are the most well-known breast and ovarian cancer susceptibility genes. Additional genes involved in DNA repair have been identified as predisposing to breast cancer. Recently, RAD51C, a new Fanconi Anemia gene, essential for homologous recombination repair, has been reported to be a rare hereditary breast and ovarian cancer susceptibility gene. Indeed, several pathogenic mutations have been identified in BRCA1/BRCA2-negative hereditary breast and ovarian cancer families. Here, we present the results of the screening of RAD51C mutations in a large series of 516 BRCA1/BRCA2-negative Spanish patients from breast and/or ovarian cancer families, and the evaluation of these results in the context of all RAD51C carriers. RAD51C mutation screening was performed by DNA analysis for all index cases. All the genetic variants identified were analyzed in silico for splicing and protein predictions. cDNA analysis was performed for three selected variants. All previous RAD51C mutation studies on breast and/or ovarian cancer were reviewed. We identified three inactivating RAD51C mutations. Two mutations were found in breast and ovarian cancer families and one mutation in a site-specific breast cancer family. Based on the mean age of ovarian cancer diagnosis in RAD51C carriers, we would recommend prophylactic bilateral salpingo-ophorectomy in premenopausal RAD51C mutation carriers. Our results support that RAD51C is a rare breast and ovarian cancer susceptibility gene and may contribute to a small fraction of families including breast and ovarian cancer cases and families with only breast cancer. Thus, RAD51C testing should be offered to hereditary breast and/or ovarian cancer families without selecting for specific cancer origin.

Large genomic rearrangements of BRCA1 and BRCA2 among patients referred for genetic analysis in Galicia (NW Spain): delimitation and mechanism of three novel BRCA1 rearrangements.

In the Iberian Peninsula, which includes mainly Spain and Portugal, large genomic rearrangements (LGRs) of BRCA1 and BRCA2 have respectively been found in up to 2.33% and 8.4% of families with hereditary breast and/or ovarian cancer (HBOC) that lack point mutations and small indels. In Galicia (Northwest Spain), the spectrum and frequency of BRCA1/BRCA2 point mutations differs from the rest of the Iberian populations. However, to date there are no Galician frequency reports of BRCA1/BRCA2 LGRs. Here we used multiplex ligation-dependent probe amplification (MLPA) to screen 651 Galician index cases (out of the 830 individuals referred for genetic analysis) without point mutations or small indels. We identified three different BRCA1 LGRs in four families. Two of them have been previously classified as pathogenic LGRs: the complete deletion of BRCA1 (identified in two unrelated families) and the deletion of exons 1 to 13. We also identified the duplication of exons 1 and 2 that is a LGR with unknown pathogenicity. Determination of the breakpoints of the BRCA1 LGRs using CNV/SNP arrays and sequencing identified them as NG_005905.2:g.70536_180359del, NG_005905.2:g.90012_97270dup, and NC_000017.10:g.41230935_41399840delinsAluSx1, respectively; previous observations of BRCA1 exon1-24del, exon1-2dup, and exon1-13del LGRs have not characterized them in such detail. All the BRCA1 LGRs arose from unequal homologous recombination events involving Alu elements. We also detected, by sequencing, one BRCA2 LGR, the Portuguese founder mutation c.156_157insAluYa5. The low frequency of BRCA1 LGRs within BRCA1 mutation carriers in Galicia (2.34%, 95% CI: 0.61-7.22) seems to differ from the Spanish population (9.93%, 95% CI: 6.76-14.27, P-value = 0.013) and from the rest of the Iberian population (9.76%, 95% CI: 6.69-13.94, P-value = 0.014).

Comprehensive annotation of splice junctions supports pervasive alternative splicing at the BRCA1 locus: a report from the ENIGMA consortium.

Loss-of-function germline mutations in BRCA1 (MIM #113705) confer markedly increased risk of breast and ovarian cancer. The full-length transcript codifies for a protein involved in DNA repair pathways and cell-cycle checkpoints. Several BRCA1 splicing isoforms have been described in public domain databases, but the physiological role (if any) of BRCA1 alternative splicing remains to be established. An accurate description of 'naturally occurring' alternative splicing at this locus is a prerequisite to understand its biological significance. However, a systematic analysis of alternative splicing at the BRCA1 locus is yet to be conducted. Here, the Evidence-Based Network for the Interpretation of Germ-Line Mutant Alleles consortium combines RT-PCR, exon scanning, cloning, sequencing and relative semi-quantification to describe naturally occurring BRCA1 alternative splicing with unprecedented resolution. The study has been conducted in blood-related RNA sources, commonly used for clinical splicing assays, as well as in one healthy breast tissue. We have characterized a total of 63 BRCA1 alternative splicing events, including 35 novel findings. A minimum of 10 splicing events (Δ1Aq, Δ5, Δ5q, Δ8p, Δ9, Δ(9,10), Δ9_11, Δ11q, Δ13p and Δ14p) represent a substantial fraction of the full-length expression level (ranging from 5 to 100%). Remarkably, our data indicate that BRCA1 alternative splicing is similar in blood and breast, a finding supporting the clinical relevance of blood-based in vitro splicing assays. Overall, our data suggest an alternative splicing model in which most non-mutually exclusive alternative splicing events are randomly combined into individual mRNA molecules to produce hundreds of different BRCA1 isoforms.

Comparison of mRNA splicing assay protocols across multiple laboratories: recommendations for best practice in standardized clinical testing.

BACKGROUND: Accurate evaluation of unclassified sequence variants in cancer predisposition genes is essential for clinical management and depends on a multifactorial analysis of clinical, genetic, pathologic, and bioinformatic variables and assays of transcript length and abundance. The integrity of assay data in turn relies on appropriate assay design, interpretation, and reporting. METHODS: We conducted a multicenter investigation to compare mRNA splicing assay protocols used by members of the ENIGMA (Evidence-Based Network for the Interpretation of Germline Mutant Alleles) consortium. We compared similarities and differences in results derived from analysis of a panel of breast cancer 1, early onset (BRCA1) and breast cancer 2, early onset (BRCA2) gene variants known to alter splicing (BRCA1: c.135-1G>T, c.591C>T, c.594-2A>C, c.671-2A>G, and c.5467+5G>C and BRCA2: c.426-12_8delGTTTT, c.7988A>T, c.8632+1G>A, and c.9501+3A>T). Differences in protocols were then assessed to determine which elements were critical in reliable assay design. RESULTS: PCR primer design strategies, PCR conditions, and product detection methods, combined with a prior knowledge of expected alternative transcripts, were the key factors for accurate splicing assay results. For example, because of the position of primers and PCR extension times, several isoforms associated with BRCA1, c.594-2A>C and c.671-2A>G, were not detected by many sites. Variation was most evident for the detection of low-abundance transcripts (e.g., BRCA2 c.8632+1G>A Δ19,20 and BRCA1 c.135-1G>T Δ5q and Δ3). Detection of low-abundance transcripts was sometimes addressed by using more analytically sensitive detection methods (e.g., BRCA2 c.426-12_8delGTTTT ins18bp). CONCLUSIONS: We provide recommendations for best practice and raise key issues to consider when designing mRNA assays for evaluation of unclassified sequence variants.

About 1% of the breast and ovarian Spanish families testing negative for BRCA1 and BRCA2 are carriers of RAD51D pathogenic variants.

RAD51D mutations have been recently identified in breast (BC) and ovarian cancer (OC) families. Although an etiological role in OC appears to be present, the association of RAD51D mutations and BC risk is more unclear. We aimed to determine the prevalence of germline RAD51D mutations in Spanish BC/OC families negative for BRCA1/BRCA2 mutations. We analyzed 842 index patients: 491 from BC/OC families, 171 BC families, 51 OC families and 129 patients without family history but with early-onset BC or OC or metachronous BC and OC. Mutation detection was performed with high-resolution melting, denaturing high-performance liquid chromatography or Sanger sequencing. Three mutations were found in four families with BC and OC cases (0.82%). Two were novel: c.1A>T (p.Met1?) and c.667+2_667+23del, leading to the exon 7 skipping and one previously described: c.674C>T (p.Arg232*). All were present in BC/OC families with only one OC. The c.667+2_667+23del cosegregated in the family with one early-onset BC and two bilateral BC cases. We also identified the c.629C>T (p.Ala210Val) variant, which was predicted in silico to be potentially pathogenic. About 1% of the BC and OC Spanish families negative for BRCA1/BRCA2 are carriers of RAD51D mutations. The presence of several BC mutation carriers, albeit in the context of familial OC, suggests an increased risk for BC, which should be taken into account in the follow-up and early detection measures. RAD51D testing should be considered in clinical setting for families with BC and OC, irrespective of the number of OC cases in the family.

Annular elastolytic giant cell granuloma associated to late-onset X-linked dominant protoporphyria.

X-linked dominant protoporphyria (XLDPP) is a genetic disorder that affects the synthesis of the heme group due to an increase in delta-aminolaevulinate synthase 2 (ALAS2) enzyme activity. Moreover, annular elastolytic giant-cell granuloma (AEGCG) is a rare reactive granulomatous dermatosis, usually associated with actinic damage. An 86-year-old man presented with edematous-erythematous lesions in photoexposed areas of the face and on the dorsum of both hands. Protoporphyrin levels in serum and feces were significantly elevated and a heterozygous frameshift mutation in the exon 11 of the ALAS2 gene: c.1706-1709del (p.Glu569GlyfsX24) was identified. Concomitantly, we observed an annular plaque with raised borders on the back of his right hand, clinically and histologically compatible with a diagnosis of AEGCG. Skin lesions disappeared only upon use of a physical sunscreen. We report two rare photodermatoses in an elderly patient and discuss the significance of dermal elastic fiber damage induced by the XLDPP as a main triggering factor of AEGCG.

Evaluation of rare variants in the new fanconi anemia gene ERCC4 (FANCQ) as familial breast/ovarian cancer susceptibility alleles.

Recently, it has been reported that biallelic mutations in the ERCC4 (FANCQ) gene cause Fanconi anemia (FA) subtype FA-Q. To investigate the possible role of ERCC4 in breast and ovarian cancer susceptibility, as occurs with other FA genes, we screened the 11 coding exons and exon-intron boundaries of ERCC4 in 1573 index cases from high-risk Spanish familial breast and ovarian cancer pedigrees that had been tested negative for BRCA1 and BRCA2 mutations and 854 controls. The frequency of ERCC4 mutation carriers does not differ between cases and controls, suggesting that ERCC4 is not a cancer susceptibility gene. Interestingly, the prevalence of ERCC4 mutation carriers (one in 288) is similar to that reported for FANCA, whereas there are approximately 100-fold more FA-A than FA-Q patients, indicating that most biallelic combinations of ERCC4 mutations are embryo lethal. Finally, we identified additional bone-fide FA ERCC4 mutations specifically disrupting interstrand cross-link repair.

Evaluation of a 5-tier scheme proposed for classification of sequence variants using bioinformatic and splicing assay data: inter-reviewer variability and promotion of minimum reporting guidelines.

Splicing assays are commonly undertaken in the clinical setting to assess the clinical relevance of sequence variants in disease predisposition genes. A 5-tier classification system incorporating both bioinformatic and splicing assay information was previously proposed as a method to provide consistent clinical classification of such variants. Members of the ENIGMA Consortium Splicing Working Group undertook a study to assess the applicability of the scheme to published assay results, and the consistency of classifications across multiple reviewers. Splicing assay data were identified for 235 BRCA1 and 176 BRCA2 unique variants, from 77 publications. At least six independent reviewers from research and/or clinical settings comprehensively examined splicing assay methods and data reported for 22 variant assays of 21 variants in four publications, and classified the variants using the 5-tier classification scheme. Inconsistencies in variant classification occurred between reviewers for 17 of the variant assays. These could be attributed to a combination of ambiguity in presentation of the classification criteria, differences in interpretation of the data provided, nonstandardized reporting of results, and the lack of quantitative data for the aberrant transcripts. We propose suggestions for minimum reporting guidelines for splicing assays, and improvements to the 5-tier splicing classification system to allow future evaluation of its performance as a clinical tool.

Characterization of BRCA1 and BRCA2 splicing variants: a collaborative report by ENIGMA consortium members.

Mutations in BRCA1 and BRCA2 predispose carriers to early onset breast and ovarian cancer. A common problem in clinical genetic testing is interpretation of variants with unknown clinical significance. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium was initiated to evaluate and implement strategies to characterize the clinical significance of BRCA1 and BRCA2 variants. As an initial project of the ENIGMA Splicing Working Group, we report splicing and multifactorial likelihood analysis of 25 BRCA1 and BRCA2 variants from seven different laboratories. Splicing analysis was performed by reverse transcriptase PCR or mini gene assay, and sequencing to identify aberrant transcripts. The findings were compared to bioinformatic predictions using four programs. The posterior probability of pathogenicity was estimated using multifactorial likelihood analysis, including co-occurrence with a deleterious mutation, segregation and/or report of family history. Abnormal splicing patterns expected to lead to a non-functional protein were observed for 7 variants (BRCA1 c.441+2T>A, c.4184_4185+2del, c.4357+1G>A, c.4987-2A>G, c.5074G>C, BRCA2 c.316+5G>A, and c.8754+3G>C). Combined interpretation of splicing and multifactorial analysis classified an initiation codon variant (BRCA2 c.3G>A) as likely pathogenic, uncertain clinical significance for 7 variants, and indicated low clinical significance or unlikely pathogenicity for another 10 variants. Bioinformatic tools predicted disruption of consensus donor or acceptor sites with high sensitivity, but cryptic site usage was predicted with low specificity, supporting the value of RNA-based assays. The findings also provide further evidence that clinical RNA-based assays should be extended from analysis of invariant dinucleotides to routinely include all variants located within the donor and acceptor consensus splicing sites. Importantly, this study demonstrates the added value of collaboration between laboratories, and across disciplines, to collate and interpret information from clinical testing laboratories to consolidate patient management.

CDK redundancy guarantees cell cycle progression in Rb-negative tumor cells independently of their p16 status.

A high percentage of tumor cells bear mutations in the Rb tumor suppressor gene. They have high levels of the cdk inhibitor p16(ink4a) and no cyclin D/cdk4-6 complexes. Although p16 is known not to arrest the proliferation of Rb-negative cells, it is not known whether its presence affects how their cycle progresses, how E2F-dependent transcription is modulated, and whether and how Rb-related proteins are inactivated. We have assessed the relevance of p16(ink4a) for cell cycle progression of these cells. Using SaOS2 osteosarcoma cells as a model, we find that downregulation of p16(ink4a) by RNAi and reconstitution of active cyclin D/cdk4 complexes does not affect progression of the cycle of these cells or expression of E2F-dependent genes. Rb-negative tumor cells can functionally inactivate Rb-related proteins in G(1)/S transition independently of their p16(ink4a) status. Furthermore, Rb-negative tumor cells do not arrest when cdk1, cdk2 or cdk3 are inhibited by RNAi, independently of their p16(ink4a) status, and combined inhibition of these cdks is also not enough to arrest their cell cycle. However, cell cycle progression of Rb-negative tumor cells is sensitive to complete cdk inhibition, as it is arrested by the chemical cdk inhibitor roscovitine and the biological cdk inhibitor p27. These results suggest that, despite their lack of cyclin D-containing complexes, Rb-negative tumor cells, like normal untransformed cells, take advantage of the high degree of redundancy of cdks for their cell cycle progression.