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.

The BRCA1-Δ11q Alternative Splice Isoform Bypasses Germline Mutations and Promotes Therapeutic Resistance to PARP Inhibition and Cisplatin

Breast and ovarian cancer patients harboring BRCA1/2 germline mutations have clinically benefitted from therapy with PARP inhibitor (PARPi) or platinum compounds, but acquired resistance limits clinical impact. In this study, we investigated the impact of mutations on BRCA1 isoform expression and therapeutic response. Cancer cell lines and tumors harboring mutations in exon 11 of BRCA1 express a BRCA1-Δ11q splice variant lacking the majority of exon 11. The introduction of frameshift mutations to exon 11 resulted in nonsense-mediated mRNA decay of full-length, but not the BRCA1-Δ11q isoform. CRISPR/Cas9 gene editing as well as overexpression experiments revealed that the BRCA1-Δ11q protein was capable of promoting partial PARPi and cisplatin resistance relative to full-length BRCA1, both in vitro and in vivo Furthermore, spliceosome inhibitors reduced BRCA1-Δ11q levels and sensitized cells carrying exon 11 mutations to PARPi treatment. Taken together, our results provided evidence that cancer cells employ a strategy to remove deleterious germline BRCA1 mutations through alternative mRNA splicing, giving rise to isoforms that retain residual activity and contribute to therapeutic resistance. Cancer Res; 76(9); 2778-90. ©2016 AACR.

Clinical Variant Classification: A Comparison of Public Databases and a Commercial Testing Laboratory

BACKGROUND

There is a growing move to consult public databases following receipt of a genetic test result from a clinical laboratory; however, the well-documented limitations of these databases call into question how often clinicians will encounter discordant variant classifications that may introduce uncertainty into patient management. Here, we evaluate discordance in BRCA1 and BRCA2 variant classifications between a single commercial testing laboratory and a public database commonly consulted in clinical practice.

MATERIALS AND METHODS

BRCA1 and BRCA2 variant classifications were obtained from ClinVar and compared with the classifications from a reference laboratory. Full concordance and discordance were determined for variants whose ClinVar entries were of the same pathogenicity (pathogenic, benign, or uncertain). Variants with conflicting ClinVar classifications were considered partially concordant if ≥1 of the listed classifications agreed with the reference laboratory classification.

RESULTS

Four thousand two hundred and fifty unique BRCA1 and BRCA2 variants were available for analysis. Overall, 73.2% of classifications were fully concordant and 12.3% were partially concordant. The remaining 14.5% of variants had discordant classifications, most of which had a definitive classification (pathogenic or benign) from the reference laboratory compared with an uncertain classification in ClinVar (14.0%).

CONCLUSION

Here, we show that discrepant classifications between a public database and single reference laboratory potentially account for 26.7% of variants in BRCA1 and BRCA2. The time and expertise required of clinicians to research these discordant classifications call into question the practicality of checking all test results against a database and suggest that discordant classifications should be interpreted with these limitations in mind.

IMPLICATIONS FOR PRACTICE

With the increasing use of clinical genetic testing for hereditary cancer risk, accurate variant classification is vital to ensuring appropriate medical management. There is a growing move to consult public databases following receipt of a genetic test result from a clinical laboratory; however, we show that up to 26.7% of variants in BRCA1 and BRCA2 have discordant classifications between ClinVar and a reference laboratory. The findings presented in this paper serve as a note of caution regarding the utility of database consultation.

Classification of the spliceogenic BRCA1 c.4096+3A>G variant as likely benign based on cosegregation data and identification of a healthy homozygous carrier

BRCA1, c.4096+3A>G was identified in a consanguineous Danish family with several cases of breast/ovarian cancer. In silico analysis and splicing assays indicated that the variant caused aberrant splicing. However, based on segregation data and the finding of a healthy homozygous carrier, we classify the BRCA1 c.4096+3A>G variant as likely benign.

Gene and Variant Annotation for Mendelian Disorders in the Era of Advanced Sequencing Technologies

Comprehensive annotations of genetic and noncoding regions and corresponding accurate variant classification for Mendelian diseases are the next big challenge in the new genomic era of personalized medicine. Progress in the development of faster and more accurate pipelines for genome annotation and variant classification will lead to the discovery of more novel disease associations and candidate therapeutic targets. This ultimately will facilitate better patient recruitment in clinical trials. In this review, we describe the trends in research at the intersection of basic and clinical genomics that aims to increase understanding of overall genomic complexity, complex inheritance patterns of disease, and patient-phenotype-specific genomic associations. We describe the emerging field of translational functional genomics, which integrates other functional "-omics" approaches that support next-generation sequencing genomic data in order to facilitate personalized diagnostics, disease management, biomarker discovery, and medicine. We also discuss the utility of this integrated approach for diagnostic clinics and medical databases and its role in the future of personalized medicine.

Targeting oncogenic vulnerabilities in triple negative breast cancer: biological bases and ongoing clinical studies

Triple negative breast cancer (TNBC) is still an incurable disease despite the great scientific effort performed during the last years. The huge heterogeneity of this disease has motivated the evaluation of a great number of therapies against different molecular alterations. In this article, we review the biological bases of this entity and how the known molecular evidence supports the current preclinical and clinical development of new therapies. Special attention will be given to ongoing clinical studies and potential options for future drug combinations.

Functional classification of DNA variants by hybrid minigenes: Identification of 30 spliceogenic variants of BRCA2 exons 17 and 18

Mutation screening of the breast cancer genes BRCA1 and BRCA2 identifies a large fraction of variants of uncertain clinical significance (VUS) whose functional and clinical interpretations pose a challenge for genomic medicine. Likewise, an increasing amount of evidence indicates that genetic variants can have deleterious effects on pre-mRNA splicing. Our goal was to investigate the impact on splicing of a set of reported variants of BRCA2 exons 17 and 18 to assess their role in hereditary breast cancer and to identify critical regulatory elements that may constitute hotspots for spliceogenic variants. A splicing reporter minigene with BRCA2 exons 14 to-20 (MGBR2_ex14-20) was constructed in the pSAD vector. Fifty-two candidate variants were selected with splicing prediction programs, introduced in MGBR2_ex14-20 by site-directed mutagenesis and assayed in triplicate in MCF-7 cells. Wild type MGBR2_ex14-20 produced a stable transcript of the expected size (1,806 nucleotides) and structure (V1-[BRCA2_exons_14–20]–V2). Functional mapping by microdeletions revealed essential sequences for exon recognition on the 3’ end of exon 17 (c.7944-7973) and the 5’ end of exon 18 (c.7979-7988, c.7999-8013). Thirty out of the 52 selected variants induced anomalous splicing in minigene assays with >16 different aberrant transcripts, where exon skipping was the most common event. A wide range of splicing motifs were affected including the canonical splice sites (15 variants), novel alternative sites (3 variants), the polypyrimidine tract (3 variants) and enhancers/silencers (9 variants). According to the guidelines of the American College of Medical Genetics and Genomics (ACMG), 20 variants could be classified as pathogenic (c.7806-2A>G, c.7806-1G>A, c.7806-1G>T, c.7806-1_7806-2dup, c.7976+1G>A, c.7977-3_7978del, c.7977-2A>T, c.7977-1G>T, c.7977-1G>C, c.8009C>A, c.8331+1G>T and c.8331+2T>C) or likely pathogenic (c.7806-9T>G, c.7976G>C, c.7976G>A, c.7977-7C>G, c.7985C>G, c.8023A>G, c.8035G>T and c.8331G>A), accounting for 30.8% of all pathogenic/likely pathogenic variants of exons 17–18 at the BRCA Share database. The remaining 8 variants (c.7975A>G, c.7977-6T>G, c.7988A>T, c.7992T>A, c.8007A>G, c.8009C>T, c.8009C>G, and c.8072C>T) induced partial splicing anomalies with important ratios of the full-length transcript (≥70%), so that they remained classified as VUS. Aberrant splicing is therefore especially prevalent in BRCA2 exons 17 and 18 due to the presence of active ESEs involved in exon recognition. Splicing functional assays with minigenes are a valuable strategy for the initial characterization of the splicing outcomes and the subsequent clinical interpretation of variants of any disease-gene, although these results should be checked, whenever possible, against patient RNA.

A significant proportion of disease-causing mutations of inherited disorders impair splicing. Massive sequencing projects of genetic diseases generate thousands of sequence variations that require functional and clinical interpretations. We have shown that splicing reporter minigenes of the breast cancer genes BRCA1 and BRCA2 are useful tools to functionally test DNA variants. In this work, we have constructed a 7-exon BRCA2 minigene (exons 14 to 20) where we mapped critical splicing regulatory sequences and tested 52 selected variants of exons 17 and 18 detected in breast cancer patients. We finely located three DNA segments on both exons that presumably contain splicing enhancer sequences. We observed that a total of 30 variants of any type disrupted the splicing patterns and, given the severity of their outcomes, we classified 20 of them as pathogenic or likely pathogenic. We also showed that a wide range of splicing elements were affected including canonical and novel 5’ and 3’ splice sites, the polypyrimidine tract and enhancer and silencer sequences. We concluded that splicing aberrations are frequent in Hereditary Breast and Ovarian Cancer and that minigenes are valuable tools to functionally classify DNA variants of any human disease gene under the splicing viewpoint.

Novel splice-switching oligonucleotide promotes BRCA1 aberrant splicing and susceptibility to PARP inhibitor action

Tumors carrying hereditary mutations in BRCA1, which attenuate the BRCA1 DNA damage repair pathway, are more susceptible to dual treatment with PARP inhibitors and DNA damaging therapeutics. Conversely, breast cancer tumors with nonmutated functional BRCA1 are less sensitive to PARP inhibition. We describe a method that triggers susceptibility to PARP inhibition in BRCA1-functional tumor cells. BRCA1 exon 11 is a key for the function of BRCA1 in DNA damage repair. Analysis of the BRCA1 exon 11 splicing mechanism identified a key region within this exon which, when deleted, induced exon 11 skipping. An RNA splice-switching oligonucleotide (SSO) developed to target this region was shown to artificially stimulate skipping of exon 11 in endogenous BRCA1 pre-mRNA. SSO transfection rendered wild-type BRCA1 expressing cell lines more susceptible to PARP inhibitor treatment, as demonstrated by a reduction in cell survival at all SSO concentrations tested. Combined SSO and PARP inhibitor treatment increased γH2AX expression indicating that SSO-dependent skipping of BRCA1 exon 11 was able to promote DSBs and therefore synthetic lethality. In conclusion, this SSO provides a new potential therapeutic strategy for targeting BRCA1-functional breast cancer by enhancing the effect of PARP inhibitors.

Hereditary breast and ovarian cancer: new genes in confined pathways

Genetic abnormalities in the DNA repair genes BRCA1 and BRCA2 predispose to hereditary breast and ovarian cancer (HBOC). However, only approximately 25% of cases of HBOC can be ascribed to BRCA1 and BRCA2 mutations. Recently, exome sequencing has uncovered substantial locus heterogeneity among affected families without BRCA1 or BRCA2 mutations. The new pathogenic variants are rare, posing challenges to estimation of risk attribution through patient cohorts. In this Review article, we examine HBOC genes, focusing on their role in genome maintenance, the possibilities for functional testing of putative causal variants and the clinical application of new HBOC genes in cancer risk management and treatment decision-making.

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.

Adding In Silico Assessment of Potential Splice Aberration to the Integrated Evaluation of BRCA Gene Unclassified Variants

Clinical mutation screening of the cancer susceptibility genes BRCA1 and BRCA2 generates many unclassified variants (UVs). Most of these UVs are either rare missense substitutions or nucleotide substitutions near the splice junctions of the protein coding exons. Previously, we developed a quantitative method for evaluation of BRCA gene UVs—the “integrated evaluation”—that combines a sequence analysis‐based prior probability of pathogenicity with patient and/or tumor observational data to arrive at a posterior probability of pathogenicity. One limitation of the sequence analysis‐based prior has been that it evaluates UVs from the perspective of missense substitution severity but not probability to disrupt normal mRNA splicing. Here, we calibrated output from the splice‐site fitness program MaxEntScan to generate spliceogenicity‐based prior probabilities of pathogenicity for BRCA gene variants; these range from 0.97 for variants with high probability to damage a donor or acceptor to 0.02 for exonic variants that do not impact a splice junction and are unlikely to create a de novo donor. We created a database http://priors.hci.utah.edu/PRIORS/ that provides the combined missense substitution severity and spliceogenicity‐based probability of pathogenicity for BRCA gene single‐nucleotide substitutions. We also updated the BRCA gene Ex‐UV LOVD, available at http://hci‐exlovd.hci.utah.edu, with 77 re‐evaluable variants.

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.

Unravelling the transcriptomic landscape of the major phase II UDP-glucuronosyltransferase drug metabolizing pathway using targeted RNA sequencing

A comprehensive view of the human UDP-glucuronosyltransferase (UGT) transcriptome is a prerequisite to the establishment of an individual's UGT metabolic glucuronidation signature. Here, we uncover the transcriptome landscape of the 10 human UGT gene loci in normal and tumoral metabolic tissues by targeted RNA next-generation sequencing. Alignment on the human hg19 reference genome identifies 234 novel exon-exon junctions. We recover all previously known UGT1 and UGT2 enzyme-coding transcripts and identify over 130 structurally and functionally diverse novel UGT variants. We further expose a revised genomic structure of UGT loci and provide a comprehensive repertoire of transcripts for each UGT gene. Data also uncover a remodelling of the UGT transcriptome occurring in a tissue- and tumor-specific manner. The complex alternative splicing program regulating UGT expression and protein functions is likely critical in determining detoxification capacity of an organ and stress-related responses, with significant impact on drug responses and diseases.

BRCA1 Alternative splicing landscape in breast tissue samples

Background

BRCA1 is a key protein in cell network, involved in DNA repair pathways and cell cycle. Recently, the ENIGMA consortium has reported a high number of alternative splicing (AS) events at this locus in blood-derived samples. However, BRCA1 splicing pattern in breast tissue samples is unknown. Here, we provide an accurate description of BRCA1 splicing events distribution in breast tissue samples.

Methods

BRCA1 splicing events were scanned in 70 breast tumor samples, 4 breast samples from healthy individuals and in 72 blood-derived samples by capillary electrophoresis (capillary EP). Molecular subtype was identified in all tumor samples. Splicing events were considered predominant if their relative expression level was at least the 10% of the full-length reference signal.

Results

54 BRCA1 AS events were identified, 27 of them were annotated as predominant in at least one sample. Δ5q, Δ13, Δ9, Δ5 and ▼1aA were significantly more frequently annotated as predominant in breast tumor samples than in blood-derived samples. Predominant splicing events were, on average, more frequent in tumor samples than in normal breast tissue samples (P = 0.010). Similarly, likely inactivating splicing events (PTC-NMDs, Non-Coding, Δ5 and Δ18) were more frequently annotated as predominant in tumor than in normal breast samples (P = 0.020), whereas there were no significant differences for other splicing events (No-Fs) frequency distribution between tumor and normal breast samples (P = 0.689).

Conclusions

Our results complement recent findings by the ENIGMA consortium, demonstrating that BRCA1 AS, despite its tremendous complexity, is similar in breast and blood samples, with no evidences for tissue specific AS events. Further on, we conclude that somatic inactivation of BRCA1 through spliciogenic mutations is, at best, a rare mechanism in breast carcinogenesis, albeit our data detects an excess of likely inactivating AS events in breast tumor samples.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-015-1145-9) contains supplementary material, which is available to authorized users.

Splicing analysis of 14 BRCA1 missense variants classifies nine variants as pathogenic

Pathogenic germline mutations in the BRCA1 gene predispose carriers to early onset breast and ovarian cancer. Clinical genetic screening of BRCA1 often reveals variants with uncertain clinical significance, complicating patient and family management. Therefore, functional examinations are urgently needed to classify whether these uncertain variants are pathogenic or benign. In this study, we investigated 14 BRCA1 variants by in silico splicing analysis and mini-gene splicing assay. All 14 alterations were missense variants located within the BRCT domain of BRCA1 and had previously been examined by functional analysis at the protein level. Results from a validated mini-gene splicing assay indicated that nine BRCA1 variants resulted in splicing aberrations leading to truncated transcripts and thus can be considered pathogenic (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5072C>T/p.Thr1691Ile, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala), whereas five BRCA1 variants had no effect on splicing (c.4985T>C/p.Phe1662Ser, c.5072C>A/p.Thr1691Lys, c.5153G>C/p.Trp1718Ser, c.5154G>T/p.Trp1718Cys, and c.5333A>G/p.Asp1778Gly). Eight of the variants having an effect on splicing (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala) were previously determined to have no or an uncertain effect on the protein level, whereas one variant (c.5072C>T/p.Thr1691Ile) were shown to have a strong effect on the protein level as well. In conclusion, our study emphasizes that in silico splicing prediction and mini-gene splicing analysis are important for the classification of BRCA1 missense variants located close to exon/intron boundaries.

Characterization of three alternative transcripts of the BRCA1 gene in patients with breast cancer and a family history of breast and/or ovarian cancer who tested negative for pathogenic mutations

The study of BRCA1 and BRCA2 genes and their alterations has been essential to the understanding of the development of familial breast and ovarian cancers. Many of the variants identified have an unknown pathogenic significance. These include variants which determine alternative mRNA splicing, identified in the intronic regions and those are capable of destroying the splicing ability. The aim of this study was to detect BRCA1/BRCA2 aberrant transcripts resulting from alternative splicing, in women with a known family history and/or early onset of breast and/or ovarian cancer, tested wild-type for BRCA1 and BRCA2. The identification and characterization of aberrant transcripts through the analysis of mRNA levels in blood lymphocytes may help us to recognize families otherwise misclassified as wild-type BRCA1 and BRCA2. Blood samples were collected from 13 women that had a family history of breast and/or ovarian cancer and tested negative for pathogenic mutations in the BRCA1 and BRCA2 genes. Total RNA was analyzed for the presence of BRCA1 and BRCA2 naturally occuring and pathological transcripts using RT-PCR. In 2 out of the 13 samples, 2 alternative transcripts of the BRCA1 gene were identified. These were probably pathogenic as they lacked exon 17 and exon 15, respectively, giving rise to a truncated protein. In addition to these, we identified the Δ17–19 transcript in 1 patient, which gives rise to a protein with an in-frame deletion of 69 amino acids. In conclusion, this study on alternative transcripts of the BRCA1 and BRCA2 genes revealed the presence of isoforms (prevalence of 15%) in blood samples from women with breast and ovarian cancer that were probably pathogenic, that were not detected by conventional methods of mutation screening based on direct sequencing of all coding regions, intron-exons junctions and MLPA analysis.