Fine mapping of RNA isoform diversity using an innovative targeted long-read RNA sequencing protocol with novel dedicated bioinformatics pipeline

BACKGROUND

Solving the structure of mRNA transcripts is a major challenge for both research and molecular diagnostic purposes. Current approaches based on short-read RNA sequencing and RT-PCR techniques cannot fully explore the complexity of transcript structure. The emergence of third-generation long-read sequencing addresses this problem by solving this sequence directly. However, genes with low expression levels are difficult to study with the whole transcriptome sequencing approach. To fix this technical limitation, we propose a novel method to capture transcripts of a gene panel using a targeted enrichment approach suitable for Pacific Biosciences and Oxford Nanopore Technologies platforms.

RESULTS

We designed a set of probes to capture transcripts of a panel of genes involved in hereditary breast and ovarian cancer syndrome. We present SOSTAR (iSofOrmS annoTAtoR), a versatile pipeline to assemble, quantify and annotate isoforms from long read sequencing using a new tool specially designed for this application. The significant enrichment of transcripts by our capture protocol, together with the SOSTAR annotation, allowed the identification of 1,231 unique transcripts within the gene panel from the eight patients sequenced. The structure of these transcripts was annotated with a resolution of one base relative to a reference transcript. All major alternative splicing events of the BRCA1 and BRCA2 genes described in the literature were found. Complex splicing events such as pseudoexons were correctly annotated. SOSTAR enabled the identification of abnormal transcripts in the positive controls. In addition, a case of unexplained inheritance in a family with a history of breast and ovarian cancer was solved by identifying an SVA retrotransposon in intron 13 of the BRCA1 gene.

CONCLUSIONS

We have validated a new protocol for the enrichment of transcripts of interest using probes adapted to the ONT and PacBio platforms. This protocol allows a complete description of the alternative structures of transcripts, the estimation of their expression and the identification of aberrant transcripts in a single experiment. This proof-of-concept opens new possibilities for RNA structure exploration in both research and molecular diagnostics.

Functional Characterization of the Human BRCA1 ∆11 Splicing Isoforms in Yeast

BRCA1, a crucial tumor suppressor gene, has several splicing isoforms, including Δ9-11, Δ11, and Δ11q, which lack exon 11, coding for significant portions of the protein. These isoforms are naturally present in both normal and cancerous cells, exhibiting altered activity compared to the full-length BRCA1. Despite this, the impact on cancer risk of the germline intronic variants promoting the exclusive expression of these Δ11 isoforms remains uncertain. Consequently, they are classified as variants of uncertain significance (VUS), posing challenges for traditional genetic classification methods due to their rarity and complexity. Our research utilizes a yeast-based functional assay, previously validated for assessing missense BRCA1 variants, to compare the activity of the Δ11 splicing isoforms with known pathogenic missense variants. This approach allows us to elucidate the functional implications of these isoforms and determine whether their exclusive expression could contribute to increased cancer risk. By doing so, we aim to provide insights into the pathogenic potential of intronic VUS-generating BRCA1 splicing isoforms and improve the classification of BRCA1 variants.

Reclassifying BRCA1 c.4358-2A > G and BRCA2 c.475 + 5G > C variants from "Uncertain Significance" to "Pathogenic" based on minigene assays and clinical evidence

BACKGROUND

Pathogenic variants in BRCA genes play a crucial role in the pathogenesis of ovarian cancer. Intronic variants of uncertain significance (VUS) may contribute to pathogenicity by affecting splicing. Currently, the significance of many intronic variants in BRCA has not been clarified, impacting patient treatment strategies and the management of familial cases.

METHOD

A retrospective study was conducted to analyze BRCA intronic VUS in a cohort of 707 unrelated ovarian cancer patients at a single institution from 2018 to 2023. Three splicing predictors were employed to analyze detected intronic VUS. Variants predicted to have splicing alterations were selected for further validation through minigene assays. Patient and familial investigations were conducted to comprehend cancer incidence within pedigrees and the application of poly (ADP-ribose) polymerase inhibitors (PARPi) by the patients. In accordance with the guidelines of the American College of Medical Genetics and Genomics (ACMG), the intronic VUS were reclassified based on minigene assay results and clinical evidence.

RESULT

Approximately 9.8% (69/707) of patients were identified as carriers of 67 different VUS in BRCA1/2, with four intronic variants accounting for 6% (4/67) of all VUS. Splicing predictors indicated potential splicing alterations in splicing for BRCA1 c.4358-2A>G and BRCA2 c.475+5G>C variants. Minigene assays utilizing the pSPL3 exon trapping vector revealed that these variants induced changes in splicing sites and frameshift, resulting in premature termination of translation (p. Ala1453Glyfs and p. Pro143Glyfs). According to ACMG guidelines, BRCA1 c.4358-2A>G and BRCA2 c.475+5G>C were reclassified as pathogenic variants. Pedigree investigations were conducted on patients with BRCA1 c.4358-2A>G variant, and the detailed utilization of PARPi provided valuable insights into research on PARPi resistance.

CONCLUSION

Two intronic VUS were reclassified as pathogenic variants. A precise classification of variants is crucial for the effective treatment and management of both patients and healthy carriers.

Parallel DNA/RNA NGS Using an Identical Target Enrichment Panel in the Analysis of Hereditary Cancer Predisposition

Germline DNA testing using the next-gene-ration sequencing (NGS) technology has become the analytical standard for the diagnostics of hereditary diseases, including cancer. Its increasing use places high demands on correct sample identification, independent confirmation of prioritized variants, and their functional and clinical interpretation. To streamline these processes, we introduced parallel DNA and RNA capture-based NGS using identical capture panel CZECANCA, which is routinely used for DNA analysis of hereditary cancer predisposition. Here, we present the analytical workflow for RNA sample processing and its analytical and diagnostic performance. Parallel DNA/RNA analysis allowed credible sample identification by calculating the kinship coefficient. The RNA capture-based approach enriched transcriptional targets for the majority of clinically relevant cancer predisposition genes to a degree that allowed analysis of the effect of identified DNA variants on mRNA processing. By comparing the panel and whole-exome RNA enrichment, we demonstrated that the tissue-specific gene expression pattern is independent of the capture panel. Moreover, technical replicates confirmed high reproducibility of the tested RNA analysis. We concluded that parallel DNA/RNA NGS using the identical gene panel is a robust and cost-effective diagnostic strategy. In our setting, it allows routine analysis of 48 DNA/RNA pairs using NextSeq 500/550 Mid Output Kit v2.5 (150 cycles) in a single run with sufficient coverage to analyse 226 cancer predisposition and candidate ge-nes. This approach can replace laborious Sanger confirmatory sequencing, increase testing turnaround, reduce analysis costs, and improve interpretation of the impact of variants by analysing their effect on mRNA processing.

The emergence of Fanconi anaemia type S: a phenotypic spectrum of biallelic BRCA1 mutations

BRCA1 is involved in the Fanconi anaemia (FA) pathway, which coordinates repair of DNA interstrand cross-links. FA is a rare genetic disorder characterised by bone marrow failure, cancer predisposition and congenital abnormalities, caused by biallelic mutations affecting proteins in the FA pathway. Germline monoallelic pathogenic BRCA1 mutations are known to be associated with hereditary breast/ovarian cancer, however biallelic mutations of BRCA1 were long predicted to be incompatible with embryonic viability, hence BRCA1 was not considered to be a canonical FA gene. Despite this, several patients with biallelic pathogenic BRCA1 mutations and FA-like phenotypes have been identified - defining a new FA type (FA-S) and designating BRCA1 as an FA gene. This report presents a scoping review of the cases of biallelic BRCA1 mutations identified to date, discusses the functional effects of the mutations identified, and proposes a phenotypic spectrum of BRCA1 mutations based upon available clinical and genetic data. We report that this FA-S cohort phenotype includes short stature, microcephaly, facial dysmorphisms, hypo/hyperpigmented lesions, intellectual disability, chromosomal sensitivity to crosslinking agents and predisposition to breast/ovarian cancer and/or childhood cancers, with some patients exhibiting sensitivity to chemotherapy. Unlike most other types of FA, FA-S patients lack bone marrow failure.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

Analysis of BRCA1 and BRCA2 alternative splicing in predisposition to ovarian cancer

BACKGROUND

The mRNA splicing is regulated on multiple levels, resulting in the proper distribution of genes' transcripts in each cell and maintaining cell homeostasis. At the same time, the expression of alternative transcripts can change in response to underlying genetic variants, often missed during routine diagnostics.

AIM

The main aim of this study was to define the frequency of aberrant splicing in BRCA1 and BRCA2 genes in blood RNA extracted from ovarian cancer patients who were previously found negative for the presence of pathogenic alterations in the 25 most commonly analysed ovarian cancer genes, including BRCA1 and BRCA2.

MATERIAL AND METHODS

Frequency and spectrum of splicing alterations in BRCA1 and BRCA2 genes were analysed in blood RNA from 101 ovarian cancer patients and healthy controls (80 healthy women) using PCR followed by gel electrophoresis and Sanger sequencing. The expression of splicing events was examined using RT-qPCR.

RESULTS

We did not identify any novel, potentially pathogenic splicing alterations. Nevertheless, we detected six naturally occurring transcripts, named BRCA1ΔE9-10, BRCA1ΔE11, BRCA1ΔE11q, and BRCA2ΔE3, BRCA2ΔE12 and BRCA2ΔE17-18 of which three (BRCA1ΔE11q, BRCA1ΔE11 and BRCA2ΔE3) were significantly higher expressed in the ovarian cancer cohort than in healthy controls (p ≤ 0.0001).

CONCLUSIONS

This observation indicates that the upregulation of selected naturally occurring transcripts can be stimulated by non-genetic mechanisms and be a potential systemic response to disease progression and/or treatment. However, this hypothesis requires further examination.

A novel BRCA1 splicing variant detected in an early onset triple-negative breast cancer patient additionally carrying a pathogenic variant in ATM: A case report

The widespread adoption of gene panel testing for cancer predisposition is leading to the identification of an increasing number of individuals with clinically relevant allelic variants in two or more genes. The potential combined effect of these variants on cancer risks is mostly unknown, posing a serious problem for genetic counseling in these individuals and their relatives, in whom the variants may segregate singly or in combination. We report a female patient who developed triple-negative high grade carcinoma in the right breast at the age of 36 years. The patient underwent bilateral mastectomy followed by combined immunotherapy and chemotherapy (IMpassion030 clinical trial). Two years later she developed a skin recurrence on the right anterior chest wall. Despite intensive treatment, the patient died at 40-year-old due to disease progression. Gene panel testing of patient's DNA revealed the presence of a protein truncating variant in ATM [c.1672G>T; p.(Gly558Ter)] and of a not previously reported variant in the BRCA1 exon 22 donor splice site [c.5406+6T>C], whose clinical significance was unknown. The analysis of patient's RNA revealed the up-regulation of two alternative BRCA1 mRNA isoforms derived from skipping of exon 22 and of exons 22-23. The corresponding predicted protein products, p.(Asp1778GlyfsTer27) and p.(Asp1778_His1822del) are both expected to affect the BRCA1 C Terminus (BRCT) domain. The two variants were observed to co-occur also in the proband's brother who, in addition, was heterozygous for a common variant (c.4837A>G) mapped to BRCA1 exon 16. This allowed to ascertain, by transcript-specific amplification, the lack of functional mRNA isoforms expressed by the c.5406+6T>C allele and provided evidence to classify the BRCA1 variant as pathogenic, according to the guidelines of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium. To our knowledge, excluding two cases detected following the screening of population specific recurrent variants, only one ATM/BRCA1 double heterozygote has been reported in the literature, being the case here described the one with the youngest age at cancer onset. The systematic collection of cases with pathogenic variants in more than one cancer predisposition gene is needed to verify if they deserve ad hoc counseling and clinical management.

Assessment of small in-frame indels and C-terminal nonsense variants of BRCA1 using a validated functional assay

BRCA1 (Breast Cancer 1, early onset) is linked to breast and ovarian cancer predisposition. Still, the risks conferred by a significant portion of BRCA1 variants identified in the population remains unknown. Most of these variants of uncertain significance are missense alterations. However, the functional implications of small in-frame deletions and/or insertions (indels) are also difficult to predict. Our group has previously evaluated the functional impact of 347 missense variants using an extensively validated transcriptional activity assay. Here we show a systematic assessment of 30 naturally occurring in-frame indels located at the C-terminal region of BRCA1. We identified positions sensitive and tolerant to alterations, expanding the knowledge of structural determinants of BRCA1 function. We further designed and assessed the impact of four single codon deletions in the tBRCT linker region and six nonsense variants at the C-terminus end of BRCA1. Amino acid substitutions, deletions or insertions in the disordered region do not significantly impact activity and are not likely to constitute pathogenic alleles. On the other hand, a sizeable fraction of in-frame indels at the BRCT domain significantly impact function. We then use a Bayesian integrative statistical model to derive the probability of pathogenicity for each variant. Our data highlights the importance of assessing the impact of small in-frame indels in BRCA1 to improve risk assessment and clinical decisions for carriers.

Comprehensive RNA and protein functional assessments contribute to the clinical interpretation of MSH2 variants causing in-frame splicing alterations

Background

Spliceogenic variants in disease-causing genes are often presumed pathogenic since most induce frameshifts resulting in loss of function. However, it was recently shown in cancer predisposition genes that some may trigger in-frame anomalies that preserve function. Here, we addressed this question by using MSH2, a DNA mismatch repair gene implicated in Lynch syndrome, as a model system.

Methods

Eighteen MSH2 variants, mostly localised within canonical splice sites, were analysed by using minigene splicing assays. The impact of the resulting protein alterations was assessed in a methylation tolerance-based assay. Clinicopathological characteristics of variant carriers were collected.

Results

Three in-frame RNA biotypes were identified based on variant-induced spliceogenic outcomes: exon skipping (E3, E4, E5 and E12), segmental exonic deletions (E7 and E15) and intronic retentions (I3, I6, I12 and I13). The 10 corresponding protein isoforms exhibit either large deletions (49–93 amino acids (aa)), small deletions (12 or 16 aa) or insertions (3–10 aa) within different functional domains. We showed that all these modifications abrogate MSH2 function, in agreement with the clinicopathological features of variant carriers.

Conclusion

Altogether, these data demonstrate that MSH2 function is intolerant to in-frame indels caused by the spliceogenic variants analysed in this study, supporting their pathogenic nature. This work stresses the importance of combining complementary RNA and protein approaches to ensure accurate clinical interpretation of in-frame spliceogenic variants.

Consistency of variant interpretations among bioinformaticians and clinical geneticists in hereditary cancer panels

Next-generation sequencing (NGS) is used increasingly in hereditary cancer patients' (HCP) management. While enabling evaluation of multiple genes simultaneously, the technology brings to light the dilemma of variant interpretation. Here, we aimed to reveal the underlying reasons for the discrepancy in the evidence titles used during variant classification according to ACMG guidelines by two different bioinformatic specialists (BIs) and two different clinical geneticists (CGs). We evaluated final reports of 1920 cancer patients and 189 different variants from 285 HCP were enrolled to the study. A total of 173 of these variants were classified as pathogenic (n = 132) and likely pathogenic (n = 41) by the BI and an additional 16 variants, that were classified as VUS by at least one interpreter and their classification would change the clinical management, were compared for their evidence titles between different specialists. The attributed evidence titles and the final classification of the variants among BIs and CGs were compared. The discrepancy between P/LP final reports was 22.5%. The discordance between CGs was 30% whereas the discordance between two BIs was almost 75%. The use of PVS1, PS3, PP3, PP5, PM1, PM2, BP1, BP4 criteria markedly varied from one expert to another. This difference was particularly noticeable in PP3, PP5, and PM1 evidence and mostly in the variants affecting splice sites like BRCA1(NM_007294.4) c.4096 + 1 G > A and CHEK2(NM_007194.4) c.592 + 3 A > T. With recent advancements in precision medicine, the importance of variant interpretations is emerging. Our study shows that variant interpretation is subjective process that is in need of concrete definitions for accurate and standard interpretation.

Systematic analysis of the relationship between ovarian cancer prognosis and alternative splicing

Background

Ovarian cancer(OC) is the gynecological tumor with the highest mortality rate, effective biomarkers are of great significance in improving its prognosis. In recent years, there have been many studies on alternative splicing (AS) events, and the role of AS events in tumor has become a focus of attention.

Methods

Data were downloaded from the TCGA database and Univariate Cox regression analysis was performed to determine AS events associated with OC prognosis.Eight prognostic models of OC were constructed in R package, and the accuracy of the models were evaluated by the time-dependent receiver operating characteristic (ROC) curves.Eight types of survival curves were drawn to evaluate the differences between the high and low risk groups.Independent prognostic factors of OC were analyzed by single factor independent analysis and multi-factor independent prognostic analysis.Again, Univariate Cox regression analysis was used to analyze the relationship between splicing factors(SF) and AS events, and Gene Ontology(GO) and Kyoto Encyclopedia of Genes and Genomes(KEGG) enrichment analysis were performed on OS-related SFs to understand the pathways.

Results

Univariate Cox regression analysis showed that among the 15,278 genes, there were 31,286 overall survival (OS) related AS events, among which 1524 AS events were significantly correlated with OS. The area under the time-dependent receiver operating characteristic curve (AUC) of AT and ME were the largest and the RI was the smallest,which were 0.757 and 0.68 respectively. The constructed models have good value for the prognosis assessment of OC patients. Among the eight survival curves, AP was the most significant difference between the high and low risk groups, with a P value of 1.61e − 1.The results of single factor independent analysis and multi-factor independent prognostic analysis showed that risk score calculated by the model and age could be used as independent risk factors.According to univariate COX regression analysis,109 SFs were correlated with AS events and adjusted in two ways: positive and negative.

Conclusions

SFs and AS events can directly or indirectly affect the prognosis of OC patients. It is very important to find effective prognostic markers to improve the survival rate of OC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13048-021-00866-1.

Altered regulation of BRCA1 exon 11 splicing is associated with breast cancer risk in carriers of BRCA1 pathogenic variants

Germline pathogenic variants in BRCA1 confer a high risk of developing breast and ovarian cancer. The BRCA1 exon 11 (formally exon 10) is one of the largest exons and codes for the nuclear localization signals of the corresponding gene product. This exon can be partially or entirely skipped during pre-mRNA splicing, leading to three major in-frame isoforms that are detectable in most cell types and tissue, and in normal and cancer settings. However, it is unclear whether the splicing imbalance of this exon is associated with cancer risk. Here we identify a common genetic variant in intron 10, rs5820483 (NC_000017.11:g.43095106_43095108dup), which is associated with exon 11 isoform expression and alternative splicing, and with the risk of breast cancer, but not ovarian cancer, in BRCA1 pathogenic variant carriers. The identification of this genetic effect was confirmed by analogous observations in mouse cells and tissue in which a loxP sequence was inserted in the syntenic intronic region. The prediction that the rs5820483 minor allele variant would create a binding site for the splicing silencer hnRNP A1 was confirmed by pull-down assays. Our data suggest that perturbation of BRCA1 exon 11 splicing modifies the breast cancer risk conferred by pathogenic variants of this gene.

New Approach for Detection of Normal Alternative Splicing Events and Aberrant Spliceogenic Transcripts with Long-Range PCR and Deep RNA Sequencing

RNA sequencing is a promising technique for detecting normal and aberrant RNA isoforms. Here, we present a new single-gene, straightforward 1-day hands-on protocol for detection of splicing alterations with deep RNA sequencing from blood. We have validated our method’s accuracy by detecting previously published normal splicing isoforms of STK11 gene. Additionally, the same technique was used to provide the first comprehensive catalogue of naturally occurring alternative splicing events of the NBN gene in blood. Furthermore, we demonstrate that our approach can be used for detection of splicing impairment caused by genetic variants. Therefore, we were able to reclassify three variants of uncertain significance: NBN:c.584G>A, STK11:c.863-5_863-3delCTC and STK11:c.615G>A. Due to the simplicity of our approach, it can be incorporated into any molecular diagnostics laboratory for determination of variant’s impact on splicing.

BRCA1 and BRCA2 whole cDNA analysis in unsolved hereditary breast/ovarian cancer patients

Germline pathogenic variants in BRCA1 and BRCA2 genes (BRCA1/2) explain an important fraction of hereditary breast/ovarian cancer (HBOC) cases. Genetic testing generally involves examining coding regions and exon/intron boundaries, thus the frequency of deleterious variants in non-coding regions is unknown. Here we analysed BRCA1/2 whole cDNA in a large cohort of 320 unsolved high-risk HBOC cases in order to identify potential splicing alterations explained by variants in BRCA1/2 deep intronic regions. Whole RNA splicing profiles were analysed by RT-PCR using Sanger sequencing or high-resolution electrophoresis in a QIAxcel instrument. Known predominant BRCA1/2 alternative splicing events were detected, together with two novel events BRCA1 ▼21 and BRCA2 Δ18q_27p. BRCA2 exon 3 skipping was detected in one patient (male) affected with breast cancer, caused by a known Portuguese founder mutation (c.156_157insAluYa5). An altered BRCA2 splicing pattern was detected in three patients, consisting in the up-regulation of ▼20A, Δ22 and ▼20A+Δ22 transcripts. In silico analysis and semi-quantitative data identified the polymorphism BRCA2 c.8755-66T>C as a potential modifier of Δ22 levels. Our findings suggest that mRNA alterations in BRCA1/2 caused by deep intronic variants are rare in Spanish population. However, RNA analysis complements DNA-based strategies allowing the identification of alterations that could go undetected by conventional testing.

Spectrum of splicing variants in disease genes and the ability of RNA analysis to reduce uncertainty in clinical interpretation

The complexities of gene expression pose challenges for the clinical interpretation of splicing variants. To better understand splicing variants and their contribution to hereditary disease, we evaluated their prevalence, clinical classifications, and associations with diseases, inheritance, and functional characteristics in a 689,321-person clinical cohort and two large public datasets. In the clinical cohort, splicing variants represented 13% of all variants classified as pathogenic (P), likely pathogenic (LP), or variants of uncertain significance (VUSs). Most splicing variants were outside essential splice sites and were classified as VUSs. Among all individuals tested, 5.4% had a splicing VUS. If RNA analysis were to contribute supporting evidence to variant interpretation, we estimated that splicing VUSs would be reclassified in 1.7% of individuals in our cohort. This would result in a clinically significant result (i.e., P/LP) in 0.1% of individuals overall because most reclassifications would change VUSs to likely benign. In ClinVar, splicing VUSs were 4.8% of reported variants and could benefit from RNA analysis. In the Genome Aggregation Database (gnomAD), splicing variants comprised 9.4% of variants in protein-coding genes; most were rare, precluding unambiguous classification as benign. Splicing variants were depleted in genes associated with dominant inheritance and haploinsufficiency, although some genes had rare variants at essential splice sites or had common splicing variants that were most likely compatible with normal gene function. Overall, we describe the contribution of splicing variants to hereditary disease, the potential utility of RNA analysis for reclassifying splicing VUSs, and how natural variation may confound clinical interpretation of splicing variants.

The BRCA1 c.788G > T (NM_007294.4) variant in a high grade serous ovarian cancer (HGSOC) patient: foods for thought

In this report we described the case of a BRCA1/2 (BRCA) molecular testing performed on tumor sample in a High Grade Serous Ovarian Cancer (HGSOC) patient with two different Next Generation Tumor Sequencing (NGTS) pipelines. The two clinical reports leaded to apparently different BRCA status, providing important foods for thought. After NGTS, the gene sequencing information (i.e., reads) are aligned to the reference gene sequences obtained from public databases, in order to provide an uniform nomenclature for unambiguous variant designation. However, the criteria adopted for variant reporting in tissue test are not always univocal. Particularly, this is the case of rare and unclassified BRCA variants for which the molecular evaluation may be a relevant challenge. Here we described a BRCA1 unclassified variant that may be re-evaluated in the context of alternative BRCA1 transcripts due to its different biological effect. We underlined that an in-depth knowledge of BRCA testing is mandatory for its appropriate use.

RNA splicing alteration in the response to platinum chemotherapy in ovarian cancer: A possible biomarker and therapeutic target

Since its discovery, alternative splicing has been recognized as a powerful way for a cell to amplify the genetic information and for a living organism to adapt, evolve, and survive. We now know that a very high number of genes are regulated by alternative splicing and that alterations of splicing have been observed in different types of human diseases, including cancer. Here, we review the accumulating knowledge that links the regulation of alternative splicing to the response to chemotherapy, focusing our attention on ovarian cancer and platinum-based treatments. Moreover, we discuss how expanding information could be exploited to identify new possible biomarkers of platinum response, to better select patients, and/or to design new therapies able to overcome platinum resistance.

Implementing NGS-based BRCA tumour tissue testing in FFPE ovarian carcinoma specimens: hints from a real-life experience within the framework of expert recommendations

AIMS

Next Generation Sequencing (NGS)-based BRCA tumour tissue testing poses several challenges. As a first step of its implementation within a regional health service network, an in-house validation study was compared with published recommendations.

METHODS

Epithelial ovarian cancer (EOC) formalin-fixed paraffin-embedded specimens stored in the archives of the eight regional pathology units were selected from a consecutive series of patients with known BRCA germline status. Two expert pathologists evaluated tumour cell content for manual macrodissection. DNA extraction, library preparation and NGS analyses were performed blinded to the germinal status. Parameters used in the study were confronted with guidelines for the validation of NGS-based oncology panels and for BRCA tumour tissue testing.

RESULTS

NGS analyses were successful in 66 of 67 EOC specimens, with good quality metrics and high reproducibility among different runs. In all, 19 BRCA pathogenic variants were identified: 12 were germline and 7 were somatic. A 100% concordance with blood tests was detected for germline variants. A BRCA1 variant showed a controversial classification. In different areas of two early stage EOCs showing somatic variants, intratumour heterogeneity not relevant for test results (variant allele frequency >5%) was observed. Compared with expert recommendations, main limitations of the study were absence of controls with known somatic BRCA status and exclusion from the validation of BRCA copy number variations (CNV).

CONCLUSIONS

A close collaboration between pathology and genetics units provides advantages in the implementation of BRCA tumour tissue testing. The development of tools for designing and interpreting complex testing in-house validation could improve process quality.

Variant effect on splicing regulatory elements, branchpoint usage, and pseudoexonization: Strategies to enhance bioinformatic prediction using hereditary cancer genes as exemplars

It is possible to estimate the prior probability of pathogenicity for germline disease gene variants based on bioinformatic prediction of variant effect/s. However, routinely used approaches have likely led to the underestimation and underreporting of variants located outside donor and acceptor splice site motifs that affect messenger RNA (mRNA) processing. This review presents information about hereditary cancer gene germline variants, outside native splice sites, with experimentally validated splicing effects. We list 95 exonic variants that impact splicing regulatory elements (SREs) in BRCA1, BRCA2, MLH1, MSH2, MSH6, and PMS2. We utilized a pre-existing large-scale BRCA1 functional data set to map functional SREs, and assess the relative performance of different tools to predict effects of 283 variants on such elements. We also describe rare examples of intronic variants that impact branchpoint (BP) sites and create pseudoexons. We discuss the challenges in predicting variant effect on BP site usage and pseudoexonization, and suggest strategies to improve the bioinformatic prioritization of such variants for experimental validation. Importantly, our review and analysis highlights the value of considering impact of variants outside donor and acceptor motifs on mRNA splicing and disease causation.

Criteria of the German Consortium for Hereditary Breast and Ovarian Cancer for the Classification of Germline Sequence Variants in Risk Genes for Hereditary Breast and Ovarian Cancer

More than ten years ago, the German Consortium for Hereditary Breast and Ovarian Cancer (GC-HBOC) set up a panel of experts (VUS Task Force) which was tasked with reviewing the classifications of genetic variants reported by individual centres of the GC-HBOC to the central database in Leipzig and reclassifying them, where necessary, based on the most recent data. When it evaluates variants, the VUS Task Force must arrive at a consensus. The resulting classifications are recorded in a central database where they serve as a basis for ensuring the consistent evaluation of previously known and newly identified variants in the different centres of the GC-HBOC. The standardised VUS evaluation by the VUS Task Force is a key element of the recall system which has also been set up by the GC-HBOC. The system will be used to pass on information to families monitored and managed by GC-HBOC centres in the event that previously classified variants are reclassified based on new information. The evaluation algorithm of the VUS Task Force was compiled using internationally established assessment methods (IARC, ACMG, ENIGMA) and is presented here together with the underlying evaluation criteria used to arrive at the classification decision using a flow chart. In addition, the characteristics and special features of specific individual risk genes associated with breast and/or ovarian cancer are discussed in separate subsections. The URLs of relevant databases have also been included together with extensive literature references to provide additional information and cover the scope and dynamism of the current state of knowledge on the evaluation of genetic variants. In future, if criteria are updated based on new information, the update will be published on the website of the GC-HBOC ( https://www.konsortium-familiaerer-brustkrebs.de/ ).

The intronic BRCA1 c.5407-25T>A variant causing partly skipping of exon 23-a likely pathogenic variant with reduced penetrance?

Rare sequence variants in the non-coding part of the BRCA genes are often reported as variants of uncertain significance (VUS), which leave patients and doctors in a challenging position. The aim of this study was to determine the pathogenicity of the BRCA1 c.5407-25T>A variant found in 20 families from Norway, France and United States with suspected hereditary breast and ovarian cancer. This was done by combining clinical and family information with allele frequency data, and assessment of the variant’s effect on mRNA splicing. Mean age at breast (n = 12) and ovarian (n = 11) cancer diagnosis in female carriers was 49.9 and 60.4 years, respectively. The mean Manchester score in the 20 families was 16.4. The allele frequency of BRCA1 c.5407-25T>A was 1/64,566 in non-Finnish Europeans (gnomAD database v2.1.1). We found the variant in 1/400 anonymous Norwegian blood donors and 0/784 in-house exomes. Sequencing of patient-derived cDNA from blood, normal breast and ovarian tissue showed that BRCA1 c.5407-25T>A leads to skipping of exon 23, resulting in frameshift and protein truncation: p.(Gly1803GlnfsTer11). Western blot analysis of transiently expressed BRCA1 proteins in HeLa cells showed a reduced amount of the truncated protein compared with wild type. Noteworthily, we found that a small amount of full-length transcript was also generated from the c.5407-25T>A allele, potentially explaining the intermediate cancer burden in families carrying this variant. In summary, our results show that BRCA1 c.5407-25T>A leads to partial skipping of exon 23, and could represent a likely pathogenic variant with reduced penetrance.

Spliceogenic analysis of BRCA1 c.439T>C (rs794727800) variant by High Resolution Melting Analysis

Correct classification of genomic variants causing potentially aberrant splicing is of utmost importance for patient management, especially in clinically actionable genes such as BRCA1/2. In this article, we report molecular evaluation of the BRCA1 c.439T>C (rs794727800, p.Leu147=) variant based on RNA of a patient suffering with high-grade serous ovarian cancer syndrome, to add new evidence to the only in silico data available for this variant. High Resolution Melting Analysis (HRMA) was used for the first time to investigate the spliceogenicity of a BRCA1 variant. HRMA with Sanger sequencing provided evidence that the c.439C allele does not cause aberrant splicing of the BRCA1 exon 7. In addition, HRMA with Sanger highlighted a different expression of the naturally occurring BRCA1 r.442_444del (c.442_444delCAG, p.Gln148del, at DNA level) isoform between blood and tumor, in this patient. HRMA is an alternative molecular approach to analyze spliceogenic properties of the c.439T>C variant and potentially for all those BRCA1/2 variants affecting splicing sites. These new evidences allowed to classify definitively the c.439T>C variant as benign. Furthermore, the different BRCA1 r.442_444del expression opens the discussion to consider a wider classification criteria for the splicing variants, including molecular evaluation at tissue level, which is an aspect currently scarcely considered in BRCA1/2 variant classification recommendations.

The BRCA1 c.4096+3A>G Variant Displays Classical Characteristics of Pathogenic BRCA1 Mutations in Hereditary Breast and Ovarian Cancers, But Still Allows Homozygous Viability

Mutations in BRCA1 result in predisposal to breast and ovarian cancers, but many variants exist with unknown clinical significance (VUS). One is BRCA1 c.4096+3A>G, which affects production of the full-length BRCA1 transcript, while augmenting transcripts lacking most or all of exon 11. Nonetheless, homozygosity of this variant has been reported in a healthy woman. We saw this variant cosegregate with breast and ovarian cancer in several family branches of four Icelandic pedigrees, with instances of phenocopies and a homozygous woman with lung cancer. We found eight heterozygous carriers (0.44%) in 1820 unselected breast cancer cases, and three (0.15%) in 1968 controls (p = 0.13). Seeking conclusive evidence, we studied tumors from carriers in the pedigrees for wild-type-loss of heterozygosity (wtLOH) and BRCA1-characteristic prevalence of estrogen receptor (ER) negativity. Of 15 breast and six ovarian tumors, wtLOH occurred in nine breast and all six ovarian tumours, and six of the nine breast tumors with wtLOH were ER-negative. These data accord with a pathogenic BRCA1-mutation. Our findings add to the current knowledge of BRCA1, and the role of its exon 11 in cancer pathogenicity, and will be of use in clinical genetic counselling.

ZC3H18 specifically binds and activates the BRCA1 promoter to facilitate homologous recombination in ovarian cancer

Reduced BRCA1 expression causes homologous recombination (HR) repair defects in high-grade serous ovarian cancers (HGSOCs). Here, we demonstrate that BRCA1 is transcriptionally activated by a previously unknown function of ZC3H18. We show that ZC3H18 is a DNA-binding protein that interacts with an E2F site in the BRCA1 promoter where it facilitates recruitment of E2F4 to an adjacent E2F site to promote BRCA1 transcription. Consistent with ZC3H18 role in activating BRCA1 expression, ZC3H18 depletion induces BRCA1 promoter methylation, reduces BRCA1 expression, disrupts HR, and sensitizes cells to DNA crosslinkers and poly(ADP-ribose) polymerase inhibitors. Moreover, in patient-derived xenografts and primary HGSOC tumors, ZC3H18 and E2F4 mRNA levels are positively correlated with BRCA1 mRNA levels, further supporting ZC3H18 role in regulating BRCA1. Given that ZC3H18 lies within 16q24.2, a region with frequent copy number loss in HGSOC, these findings suggest that ZC3H18 copy number losses could contribute to HR defects in HGSOC.

Incorporation of semi-quantitative analysis of splicing alterations for the clinical interpretation of variants in BRCA1 and BRCA2 genes

BRCA1 and BRCA2 (BRCA1/2) genetic variants that disrupt messenger RNA splicing are commonly associated with increased risks of developing breast/ovarian cancer. The majority of splicing studies published to date rely on qualitative methodologies (i.e., Sanger sequencing), but it is necessary to incorporate semi-quantitative or quantitative approaches to accurately interpret the clinical significance of spliceogenic variants. Here, we characterize the splicing impact of 31 BRCA1/2 variants using semi-quantitative capillary electrophoresis of fluorescent amplicons (CE), Sanger sequencing and allele-specific assays. A total of 14 variants were found to disrupt splicing. Allelic-specific assays could be performed for BRCA1 c.302-1G>A and BRCA2 c.516+2T>A, c.1909+1G>A, c.8332-13T>G, c.8332-2A>G, c.8954-2A>T variants, showing a monoallelic contribution to full-length transcript expression that was concordant with semi-quantitative data. The splicing fraction of alternative and aberrant transcripts was also measured by CE, facilitating variant interpretation. Following Evidence-based Network for the Interpretation of Germline Mutant Alleles criteria, we successfully classified eight variants as pathogenic (Class 5), five variants as likely pathogenic (Class 4), and 14 variants as benign (Class 1). We also provide splicing data for four variants classified as uncertain (Class 3), which produced a "leaky" splicing effect or introduced a missense change in the protein sequence, that will require further assessment to determine their clinical significance.

Diagnostic mRNA splicing assay for variants in BRCA1 and BRCA2 identified two novel pathogenic splicing aberrations

Background

Pathogenic variants in BRCA1 and BRCA2 cause hereditary breast and ovarian cancer. Screening of these genes has become easily accessible in diagnostic laboratories. Sequencing and copy number analyses are used to detect pathogenic variants, but also lead to identification of variants of unknown clinical significance (VUS). If the effect of a VUS can be clarified, it has direct consequence for the clinical management of the patient and family members. A splicing assay is one of several tools that might help in the classification of VUS. We therefore established mRNA analyses for BRCA1 and BRCA2 in the diagnostic laboratory in 2015. We hereby report the results of mRNA analysis variants in BRCA1 and BRCA2 after three years.

Methods

Variants predicted to alter splicing and variants within the canonical splice sites were selected for splicing analyses. Splicing assays were performed by reverse transcription-PCR of patient RNA. A biallalic expression analysis was carried out whenever possible.

Results

Twenty-five variants in BRCA1 and BRCA2 were analyzed by splicing assays; nine showed altered transcripts and 16 showed normal splicing patterns. The two novel pathogenic variants in BRCA1 c.4484 + 3 A > C and c.5407–10G > A were characterized.

Conclusions

We conclude that mRNA analyses are useful in characterization of variants that may affect splicing. The results can guide classification of variants from unknown clinical significance to pathogenic or benign in a diagnostic laboratory, and thus be of direct clinical importance.

BRCA1-No Matter How You Splice It

BRCA1 (breast cancer 1, early onset), a well-known breast cancer susceptibility gene, is a highly alternatively spliced gene. BRCA1 alternative splicing may serve as an alternative regulatory mechanism for the inactivation of the BRCA1 gene in both hereditary and sporadic breast cancers, and other BRCA1-associated cancers. The alternative transcripts of BRCA1 can mimic known functions, possess unique functions compared with the full-length BRCA1 transcript, and in some cases, appear to function in opposition to full-length BRCA1 In this review, we will summarize the functional "naturally occurring" alternative splicing transcripts of BRCA1 and then discuss the latest next-generation sequencing-based detection methods and techniques to detect alternative BRCA1 splicing patterns and their potential use in cancer diagnosis, prognosis, and therapy.

Mis-splicing in breast cancer: identification of pathogenic BRCA2 variants by systematic minigene assays

Splicing disruption is a common mechanism of gene inactivation associated with germline variants of susceptibility genes. To study the role of BRCA2 mis-splicing in hereditary breast/ovarian cancer (HBOC), we performed a comprehensive analysis of variants from BRCA2 exons 2-9, as well as the initial characterization of the regulatory mechanisms of such exons. A pSAD-based minigene with exons 2-9 was constructed and validated in MCF-7 cells, producing the expected transcript (1016-nt/V1-BRCA2_exons_2-9-V2). DNA variants from mutational databases were analyzed by NNSplice and Human Splicing Finder softwares. To refine ESE-variant prediction, we mapped the regulatory regions through a functional strategy whereby 26 exonic microdeletions were introduced into the minigene and tested in MCF-7 cells. Thus, we identified nine spliceogenic ESE-rich intervals where ESE-variants may be located. Combining bioinformatics and microdeletion assays, 83 variants were selected and genetically engineered in the minigene. Fifty-three changes impaired splicing: 28 variants disrupted the canonical sites, four created new ones, 10 abrogated enhancers, eight created silencers and three caused a double-effect. Notably, nine spliceogenic-ESE variants were located within ESE-containing intervals. Capillary electrophoresis and sequencing revealed more than 23 aberrant transcripts, where exon skipping was the most common event. Interestingly, variant c.67G>A triggered the usage of a noncanonical GC-donor 4-nt upstream. Thirty-six variants that induced severe anomalies (>60% aberrant transcripts) were analyzed according to the ACMG guidelines. Thus, 28 variants were classified as pathogenic, five as likely pathogenic and three as variants of uncertain significance. Interestingly, 13 VUS were reclassified as pathogenic or likely pathogenic variants. In conclusion, a large fraction of BRCA2 variants (∼64%) provoked splicing anomalies lending further support to the high prevalence of this disease-mechanism. The low accuracy of ESE-prediction algorithms may be circumvented by functional ESE-mapping that represents an optimal strategy to identify spliceogenic ESE-variants. Finally, systematic functional assays by minigenes depict a valuable tool for the initial characterization of splicing anomalies and the clinical interpretation of variants. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Prevalence of BRCA1 and BRCA2 gene mutations in Chinese patients with high-risk breast cancer

Background

Breast cancer is the most common cancer among women worldwide. Here, we report the prevalence of BRCA1/2 mutations in patients with high‐risk breast cancer from Inner Mongolia and Jilin, China, which was a part of a nationwide project on the detection of BRCA1/2 mutations in Chinese patients with hereditary breast cancer.

Methods

According to the criteria, index patients from a total of 245 independent families were initially recruited. All 49 exons of BRCA1 and BRCA2 and adjacent noncoding regions were screened for mutations based on next‐generation sequencing from collected saliva.

Results

We detected 17 BRCA1/2 variants in 18 of 216 (8.3%) index patients with high‐risk breast cancer. Among these, seven mutations were novel, including four BRCA1 mutations (c.123_124delCAinsAT, c.5093_5096delCTAA, c.5396‐2A>G, and c.2054delinsGAAGAGTAACAAGTAAGAAGAGTAACAAGAAG), and three BRCA2 mutations (c.304A>T, c.7552_7553insT, and c.9548_9549insA). The BRCA1/2 variants were identified in 14% (8/57) of the patients with triple‐negative breast cancer and in 6.3% (10/159) of the patients with non‐triple‐negative breast cancer. There was no significant difference between the two groups (p = 0.07). A higher frequency for BRCA1 mutations was observed in patients with triple‐negative breast cancer than in those with non‐triple‐negative breast cancer (12.3% vs. 2.5%, p = 0.004). The frequencies of the BRCA2 mutations were not significantly different between patients with triple‐negative breast cancer and those with non‐triple‐negative breast cancer (1.8% vs. 3.8%, p = 0.46).

Conclusion

We found that patients with triple‐negative breast cancer had a higher frequency of BRCA1 mutations than those with non‐triple‐negative breast cancer. In this study, no significant associations between the BRCA1/2 mutation status and age, family history of breast cancer, ovarian cancer, pancreatic cancer and prostate cancer, number of primary lesions, tumor size, or lymph node metastasis were observed.

Alternative splicing and ACMG-AMP-2015-based classification of PALB2 genetic variants: an ENIGMA report

BACKGROUND

PALB2 monoallelic loss-of-function germ-line variants confer a breast cancer risk comparable to the average BRCA2 pathogenic variant. Recommendations for risk reduction strategies in carriers are similar. Elaborating robust criteria to identify loss-of-function variants in PALB2-without incurring overprediction-is thus of paramount clinical relevance. Towards this aim, we have performed a comprehensive characterisation of alternative splicing in PALB2, analysing its relevance for the classification of truncating and splice site variants according to the 2015 American College of Medical Genetics and Genomics-Association for Molecular Pathology guidelines.

METHODS

Alternative splicing was characterised in RNAs extracted from blood, breast and fimbriae/ovary-related human specimens (n=112). RNAseq, RT-PCR/CE and CloneSeq experiments were performed by five contributing laboratories. Centralised revision/curation was performed to assure high-quality annotations. Additional splicing analyses were performed in PALB2 c.212-1G>A, c.1684+1G>A, c.2748+2T>G, c.3113+5G>A, c.3350+1G>A, c.3350+4A>C and c.3350+5G>A carriers. The impact of the findings on PVS1 status was evaluated for truncating and splice site variant.

RESULTS

We identified 88 naturally occurring alternative splicing events (81 newly described), including 4 in-frame events predicted relevant to evaluate PVS1 status of splice site variants. We did not identify tissue-specific alternate gene transcripts in breast or ovarian-related samples, supporting the clinical relevance of blood-based splicing studies.

CONCLUSIONS

PVS1 is not necessarily warranted for splice site variants targeting four PALB2 acceptor sites (exons 2, 5, 7 and 10). As a result, rare variants at these splice sites cannot be assumed pathogenic/likely pathogenic without further evidences. Our study puts a warning in up to five PALB2 genetic variants that are currently reported as pathogenic/likely pathogenic in ClinVar.

Usefulness and Limitations of Comprehensive Characterization of mRNA Splicing Profiles in the Definition of the Clinical Relevance of BRCA1/2 Variants of Uncertain Significance

Highly penetrant variants of BRCA1/2 genes are involved in hereditary predisposition to breast and ovarian cancer. The detection of pathogenic BRCA variants has a considerable clinical impact, allowing appropriate cancer-risk management. However, a major drawback is represented by the identification of variants of uncertain significance (VUS). Many VUS potentially affect mRNA splicing, making transcript analysis an essential step for the definition of their pathogenicity. Here, we characterize the impact on splicing of ten BRCA1/2 variants. Aberrant splicing patterns were demonstrated for eight variants whose alternative transcripts were fully characterized. Different events were observed, including exon skipping, intron retention, and usage of de novo and cryptic splice sites. Transcripts with premature stop codons or in-frame loss of functionally important residues were generated. Partial/complete splicing effect and quantitative contribution of different isoforms were assessed, leading to variant classification according to Evidence-based Network for the Interpretation of Mutant Alleles (ENIGMA) consortium guidelines. Two variants could be classified as pathogenic and two as likely benign, while due to a partial splicing effect, six variants remained of uncertain significance. The association with an undefined tumor risk justifies caution in recommending aggressive risk-reduction treatments, but prevents the possibility of receiving personalized therapies with potential beneficial effect. This indicates the need for applying additional approaches for the analysis of variants resistant to classification by gene transcript analyses.

Targeted RNA-seq successfully identifies normal and pathogenic splicing events in breast/ovarian cancer susceptibility and Lynch syndrome genes

A subset of genetic variants found through screening of patients with hereditary breast and ovarian cancer syndrome (HBOC) and Lynch syndrome impact RNA splicing. Through target enrichment of the transcriptome, it is possible to perform deep-sequencing and to identify the different and even rare mRNA isoforms. A targeted RNA-seq approach was used to analyse the naturally-occurring splicing events for a panel of 8 breast and/or ovarian cancer susceptibility genes (BRCA1, BRCA2, RAD51C, RAD51D, PTEN, STK11, CDH1, TP53), 3 Lynch syndrome genes (MLH1, MSH2, MSH6) and the fanconi anaemia SLX4 gene, in which monoallelic mutations were found in non-BRCA families. For BRCA1, BRCA2, RAD51C and RAD51D the results were validated by capillary electrophoresis and were compared to a non-targeted RNA-seq approach. We also compared splicing events from lymphoblastoid cell-lines with those from breast and ovarian fimbriae tissues. The potential of targeted RNA-seq to detect pathogenic changes in RNA-splicing was validated by the inclusion of samples with previously well characterized BRCA1/2 genetic variants. In our study, we update the catalogue of normal splicing events for BRCA1/2, provide an extensive catalogue of normal RAD51C and RAD51D alternative splicing, and list splicing events found for eight other genes. Additionally, we show that our approach allowed the identification of aberrant splicing events due to the presence of BRCA1/2 genetic variants and distinguished between complete and partial splicing events. In conclusion, targeted-RNA-seq can be very useful to classify variants based on their putative pathogenic impact on splicing.

Quantifying BRCA1 and BRCA2 mRNA Isoform Expression Levels in Single Cells

BRCA1 and BRCA2 spliceogenic variants are often associated with an elevated risk of breast and ovarian cancers. Analyses of BRCA1 and BRCA2 splicing patterns have traditionally used technologies that sample a population of cells but do not account for the variation that may be present between individual cells. This novel proof of concept study utilises RNA in situ hybridisation to measure the absolute expression of BRCA1 and BRCA2 mRNA splicing events in single lymphoblastoid cells containing known spliceogenic variants (BRCA1c.671-2 A>G or BRCA2c.7988 A>T). We observed a large proportion of cells (>42%) in each sample that did not express mRNA for the targeted gene. Increased levels (average mRNA molecules per cell) of BRCA2 ∆17_18 were observed in the cells containing the known spliceogenic variant BRCA2c.7988 A>T, but cells containing BRCA1c.671-2 A>G were not found to express significantly increased levels of BRCA1 ∆11, as had been shown previously. Instead, we show for each variant carrier sample that a higher proportion of cells expressed the targeted splicing event compared to control cells. These results indicate that BRCA1/2 mRNA is expressed stochastically, suggesting that previously reported results using RT-PCR may have been influenced by the number of cells with BRCA1/2 mRNA expression and may not represent an elevation of constitutive mRNA expression. Detection of mRNA expression in single cells allows for a more comprehensive understanding of how spliceogenic variants influence the expression of mRNA isoforms. However, further research is required to assess the utility of this technology to measure the expression of predicted spliceogenic BRCA1 and BRCA2 variants in a diagnostic setting.

Screening of BRCA1/2 deep intronic regions by targeted gene sequencing identifies the first germline BRCA1 variant causing pseudoexon activation in a patient with breast/ovarian cancer

BACKGROUND

Genetic analysis of BRCA1 and BRCA2 for the diagnosis of hereditary breast and ovarian cancer (HBOC) is commonly restricted to coding regions and exon-intron boundaries. Although germline pathogenic variants in these regions explain about ~20% of HBOC cases, there is still an important fraction that remains undiagnosed. We have screened BRCA1/2 deep intronic regions to identify potential spliceogenic variants that could explain part of the missing HBOC susceptibility.

METHODS

We analysed BRCA1/2 deep intronic regions by targeted gene sequencing in 192 high-risk HBOC families testing negative for BRCA1/2 during conventional analysis. Rare variants (MAF <0.005) predicted to create/activate splice sites were selected for further characterisation in patient RNA. The splicing outcome was analysed by RT-PCR and Sanger sequencing, and allelic imbalance was also determined when heterozygous exonic loci were present.

RESULTS

A novel transcript was detected in BRCA1 c.4185+4105C>T variant carrier. This variant promotes the inclusion of a pseudoexon in mature mRNA, generating an aberrant transcript predicted to encode for a non-functional protein. Quantitative and allele-specific assays determined haploinsufficiency in the variant carrier, supporting a pathogenic effect for this variant. Genotyping of 1030 HBOC cases and 327 controls did not identify additional carriers in Spanish population.

CONCLUSION

Screening of BRCA1/2 intronic regions has identified the first BRCA1 deep intronic variant associated with HBOC by pseudoexon activation. Although the frequency of deleterious variants in these regions appears to be low, our study highlights the importance of studying non-coding regions and performing comprehensive RNA assays to complement genetic diagnosis.

GFP-Fragment Reassembly Screens for the Functional Characterization of Variants of Uncertain Significance in Protein Interaction Domains of the BRCA1 and BRCA2 Genes

Genetic testing for BRCA1 and BRCA2 genes has led to the identification of many unique variants of uncertain significance (VUS). Multifactorial likelihood models that predict the odds ratio for VUS in favor or against cancer causality, have been developed, but their use is conditioned by the amount of necessary data, which are difficult to obtain if a variant is rare. As an alternative, variants mapping to the coding regions can be examined using in vitro functional assays. BRCA1 and BRCA2 proteins promote genome protection by interacting with different proteins. In this study, we assessed the functional effect of two sets of variants in BRCA genes by exploiting the green fluorescent protein (GFP)-reassembly in vitro assay, which was set-up to test the BRCA1/BARD1, BRCA1/UbcH5a, and BRCA2/DSS1 interactions. Based on the findings observed for the validation panels of previously classified variants, BRCA1/UbcH5a and BRCA2/DSS1 binding assays showed 100% sensitivity and specificity in identifying pathogenic and non-pathogenic variants. While the actual efficiency of these assays in assessing the clinical significance of BRCA VUS has to be verified using larger validation panels, our results suggest that the GFP-reassembly assay is a robust method to identify variants affecting normal protein functioning and contributes to the classification of VUS.

BRCA1 and Breast Cancer: a Review of the Underlying Mechanisms Resulting in the Tissue-Specific Tumorigenesis in Mutation Carriers

Since the first cloning of BRCA1 in 1994, many of its cellular interactions have been elucidated. However, its highly specific role in tumorigenesis in the breast tissue—carriers of BRCA1 mutations are predisposed to life-time risks of up to 80%—relative to many other tissues that remain unaffected, has not yet been fully enlightened. In this article, we have applied a universal model of tissue-specificity of cancer genes to BRCA1 and present a systematic review of proposed concepts classified into 4 categories. Firstly, tissue-specific differences in levels of BRCA1 expression and secondly differences in expression of proteins with redundant functions are outlined. Thirdly, cell-type specific interactions of BRCA1 are presented: its regulation of aromatase, its interaction with Progesterone- and receptor activator of nuclear factor-κB ligand-signaling that controls proliferation of luminal progenitor cells, and its influence on cell differentiation via modulation of the key regulators jagged 1-NOTCH and snail family transcriptional repressor 2. Fourthly, factors specific to the cell-type as well as the environment of the breast tissue are elucidated: distinct frequency of losses of heterozygosity, interaction with X inactivation specific transcript RNA, estrogen-dependent induction of genotoxic metabolites and nuclear factor (erythroid-derived 2)-like 2, and regulation of sirtuin 1. In conclusion, the impact of these concepts on the formation of hormone-sensitive and -insensitive breast tumors is outlined.

Additional molecular and clinical evidence open the way to definitive IARC classification of the BRCA1 c.5332G > A (p.Asp1778Asn) variant

OBJECTIVES

In silico splicing analysis, a mini-gene assay and splicing data, obtained using RNA from blood samples, have shown that the BRCA1 c.5332G > A variant induces exon 21 skipping. However, despite these evidences, up to date, this variant is unclassified. The aim of this study is to provide further molecular and clinical evidence for the BRCA1 c.5332G > A variant in a patient with high grade serous ovarian carcinoma (HGSOC) to allow a definitive classification of this variant.

DESIGN AND METHOD

The effect of the BRCA1 c.5332G > A variant on RNA splicing was evaluated by amplifying regions of BRCA1 from the cDNA of the patient. Loss of heterozygosity (LOH) in tumor tissue was also investigated.

RESULTS

The c.5332G > A allele causes significantly aberrant splicing of the BRCA1 exon 21. Evaluation of the c.5332A allele in tumor tissue highlights a possible loss of heterozygosity, supporting her pathogenic effect.

CONCLUSIONS

Our results regarding the c.5332G > A variant confirm that it contributed to predisposition and onset of ovarian carcinoma in the patient. We propose to classify this variant as 'likely-pathogenic' (class IV).

Accurate classification of BRCA1 variants with saturation genome editing

Variants of uncertain significance fundamentally limit the clinical utility of genetic information. The challenge they pose is epitomized by BRCA1, a tumour suppressor gene in which germline loss-of-function variants predispose women to breast and ovarian cancer. Although BRCA1 has been sequenced in millions of women, the risk associated with most newly observed variants cannot be definitively assigned. Here we use saturation genome editing to assay 96.5% of all possible single-nucleotide variants (SNVs) in 13 exons that encode functionally critical domains of BRCA1. Functional effects for nearly 4,000 SNVs are bimodally distributed and almost perfectly concordant with established assessments of pathogenicity. Over 400 non-functional missense SNVs are identified, as well as around 300 SNVs that disrupt expression. We predict that these results will be immediately useful for the clinical interpretation of BRCA1 variants, and that this approach can be extended to overcome the challenge of variants of uncertain significance in additional clinically actionable genes.

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.

Quantitative Analysis of BRCA1 and BRCA2 Germline Splicing Variants Using a Novel RNA-Massively Parallel Sequencing Assay

Clinical genetic testing for hereditary breast and ovarian cancer (HBOC) is becoming widespread. However, the interpretation of variants of unknown significance (VUS) in HBOC genes, such as the clinically actionable genes BRCA1 and BRCA2, remain a challenge. Among the variants that are frequently classified as VUS are those with unclear effects on splicing. In order to address this issue we developed a high-throughput RNA-massively parallel sequencing assay—CloneSeq—capable to perform quantitative and qualitative analysis of transcripts in cell lines and HBOC patients. This assay is based on cloning of RT-PCR products followed by massive parallel sequencing of the cloned transcripts. To validate this assay we compared it to the RNA splicing assays recommended by members of the ENIGMA (Evidence-based Network for the Interpretation of Germline Mutant Alleles) consortium. This comparison was performed using well-characterized lymphoblastoid cell lines (LCLs) generated from carriers of the BRCA1 or BRCA2 germline variants that have been previously described to be associated with splicing defects. CloneSeq was able to replicate the ENIGMA results, in addition to providing quantitative characterization of BRCA1 and BRCA2 germline splicing alterations in a high-throughput fashion. Furthermore, CloneSeq was used to analyze blood samples obtained from carriers of BRCA1 or BRCA2 germline sequence variants, including the novel uncharacterized alteration BRCA1 c.5152+5G>T, which was identified in a HBOC family. CloneSeq provided a high-resolution picture of all the transcripts induced by BRCA1 c.5152+5G>T, indicating it results in significant levels of exon skipping. This analysis proved to be important for the classification of BRCA1 c.5152+5G>T as a clinically actionable likely pathogenic variant. Reclassifications such as these are fundamental in order to offer preventive measures, targeted treatment, and pre-symptomatic screening to the correct individuals.

Identification of Eight Spliceogenic Variants in BRCA2 Exon 16 by Minigene Assays

Genetic testing of BRCA1 and BRCA2 identifies a large number of variants of uncertain clinical significance whose functional and clinical interpretations pose a challenge for genetic counseling. Interestingly, a relevant fraction of DNA variants can disrupt the splicing process in cancer susceptibility genes. We have tested more than 200 variants throughout 19 BRCA2 exons mostly by minigene assays, 54% of which displayed aberrant splicing, thus confirming the utility of this assay to check genetic variants in the absence of patient RNA. Our goal was to investigate BRCA2 exon 16 with a view to characterizing spliceogenic variants recorded at the mutational databases. Seventy-two different BIC and UMD variants were analyzed with NNSplice and Human Splicing Finder, 12 of which were selected because they were predicted to disrupt essential splice motifs: canonical splice sites (ss; eight variants) and exonic/intronic splicing enhancers (four variants). These 12 candidate variants were introduced into the BRCA2 minigene with seven exons (14–20) by site-directed mutagenesis and then transfected into MCF-7 cells. Seven variants (six intronic and one missense) induced complete abnormal splicing patterns: c.7618-2A>T, c.7618-2A>G, c.7618-1G>C, c.7618-1G>A, c.7805G>C, c.7805+1G>A, and c.7805+3A>C, as well as a partial anomalous outcome by c.7802A>G. They generated at least 10 different transcripts: Δ16p₄₄ (alternative 3’ss 44-nt downstream; acceptor variants), Δ16 (exon 16-skipping; donor variants), Δ16p₅₅ (alternative 3’ss 55-nt downstream), Δ16q₄ (alternative 5’ss 4-nt upstream), Δ16q₁₀₀ (alternative 5’ss 4-nt upstream), ▾16q²⁰ (alternative 5’ss 20-nt downstream), as well as minor (Δ16p₉₃ and Δ16,17p₆₉) and uncharacterized transcripts of 893 and 954 nucleotides. Isoforms Δ16p₄₄, Δ16, Δ16p₅₅, Δ16q₄, Δ16q₁₀₀, and ▾16q²⁰ introduced premature termination codons which presumably inactivate BRCA2. According to the guidelines the American College of Medical Genetics and Genomics these eight variants could be classified as pathogenic or likely pathogenic whereas the Evidence-based Network for the Interpretation of Germline Mutant Alleles rules suggested seven class 4 and one class 3 variants. In conclusion, our study highlights the relevance of splicing functional assays by hybrid minigenes for the clinical classification of genetic variations. Hence, we provide new data about spliceogenic variants of BRCA2 exon 16 that are directly correlated with breast cancer susceptibility.

Investigation of Experimental Factors That Underlie BRCA1/2 mRNA Isoform Expression Variation: Recommendations for Utilizing Targeted RNA Sequencing to Evaluate Potential Spliceogenic Variants

PCR-based RNA splicing assays are commonly used in diagnostic and research settings to assess the potential effects of variants of uncertain clinical significance in BRCA1 and BRCA2. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium completed a multicentre investigation to evaluate differences in assay design and the integrity of published data, raising a number of methodological questions associated with cell culture conditions and PCR-based protocols. We utilized targeted RNA-seq to re-assess BRCA1 and BRCA2 mRNA isoform expression patterns in lymphoblastoid cell lines (LCLs) previously used in the multicentre ENIGMA study. Capture of the targeted cDNA sequences was carried out using 34 BRCA1 and 28 BRCA2 oligonucleotides from the Illumina Truseq Targeted RNA Expression platform. Our results show that targeted RNA-seq analysis of LCLs overcomes many of the methodology limitations associated with PCR-based assays leading us to make the following observations and recommendations: (1) technical replicates (n > 2) of variant carriers to capture methodology induced variability associated with RNA-seq assays, (2) LCLs can undergo multiple freeze/thaw cycles and can be cultured up to 2 weeks without noticeably influencing isoform expression levels, (3) nonsense-mediated decay inhibitors are essential prior to splicing assays for comprehensive mRNA isoform detection, (4) quantitative assessment of exon:exon junction levels across BRCA1 and BRCA2 can help distinguish between normal and aberrant isoform expression patterns. Experimentally derived recommendations from this study will facilitate the application of targeted RNA-seq platforms for the quantitation of BRCA1 and BRCA2 mRNA aberrations associated with sequence variants of uncertain clinical significance.

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.

Break Breast Cancer Addiction by CRISPR/Cas9 Genome Editing

Breast cancer is the leading diagnosed cancer for women globally. Evolution of breast cancer in tumorigenesis, metastasis and treatment resistance appears to be driven by the aberrant gene expression and protein degradation encoded by the cancer genomes. The uncontrolled cancer growth relies on these cellular events, thus constituting the cancerous programs and rendering the addiction towards them. These programs are likely the potential anticancer biomarkers for Personalized Medicine of breast cancer. This review intends to delineate the impact of the CRSPR/Cas-mediated genome editing in identification and validation of these anticancer biomarkers. It reviews the progress in three aspects of CRISPR/Cas9-mediated editing of the breast cancer genomes: Somatic genome editing, transcription and protein degradation addictions.

Nanopore sequencing of full-length BRCA1 mRNA transcripts reveals co-occurrence of known exon skipping events

Background

Laboratory assays evaluating the effect of DNA sequence variants on BRCA1 mRNA splicing may contribute to classification by providing molecular evidence. However, our knowledge of normal and aberrant BRCA1 splicing events to date has been limited to data derived from assays targeting partial transcript sequences. This study explored the utility of nanopore sequencing to examine whole BRCA1 mRNA transcripts and to provide accurate categorisation of in-frame and out-of-frame splicing events.

Methods

The exon structure of BRCA1 transcripts from a previously studied control lymphoblastoid cell line were assessed using MinION nanopore sequencing of long-range reverse transcriptase-PCR amplicons.

Results

Our study identified and characterised 32 complete BRCA1 isoforms, including 18 novel isoforms which showed skipping of multiple contiguous and/or non-contiguous exons. Furthermore, we show that known BRCA1 exon skipping events, such as Δ(9,10) and Δ21, can co-occur in a single transcript, with some isoforms containing four or more alternative splice junctions. Fourteen novel isoforms were formed entirely from a combination of previously identified alternative splice junctions, suggesting that the total number of BRCA1 isoforms might be lower than the number of splicing events reported previously.

Conclusions

Our results highlight complexity in BRCA1 transcript structure that has not been described previously. This finding has key implications for predicting the translation frame of splicing transcripts, important for interpreting the clinical significance of spliceogenic variants. Future research is warranted to quantitatively assess full-length BRCA1 transcript levels, and to assess the application of nanopore sequencing for routine evaluation of potential spliceogenic variants.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-017-0919-1) contains supplementary material, which is available to authorized users.

Characterization of a novel germline BRCA1 splice variant, c.5332+4delA

PURPOSE

Germline mutations in BRCA1 and BRCA2 confer a significant increase in risk for cancer, and determining pathogenicity of a BRCA variant can guide the clinical management of the disease. About 1/3 of BRCA1 variants reported in the public databases have uncertain clinical significance due to lack of conclusive evidence. This study aims to characterize a novel BRCA1 deletion affecting the + 4 splice donor site identified in an individual with early-onset breast cancer.

METHODS

The effect of BRCA1 c.5332+4delA variant on RNA splicing was evaluated by amplifying regions of BRCA1 from cDNA derived from the patient. The proportion of abnormal transcript in the total transcripts was quantified. Loss of heterozygosity (LOH) in tumor tissue was investigated using Sanger sequencing and fragment analysis.

RESULTS

BRCA1 c.5332+4delA caused skipping of exon 21 in patient-derived samples. Semi-quantitative analysis indicated that this aberrant RT-PCR product accounts for about 40% of the total transcript levels. Loss of heterozygosity (LOH) was observed in patient's tumor tissue.

CONCLUSIONS

Our results indicate that the BRCA1 c.5332+4delA variant contributes to cancer predisposition through disruption of normal mRNA splicing. We classify this variant as likely pathogenic.

Alternative Splicing in Breast Cancer and the Potential Development of Therapeutic Tools

Alternative splicing is a key molecular mechanism now considered as a hallmark of cancer that has been associated with the expression of distinct isoforms during the onset and progression of the disease. The leading cause of cancer-related deaths in women worldwide is breast cancer, and even when the role of alternative splicing in this type of cancer has been established, the function of this mechanism in breast cancer biology is not completely decoded. In order to gain a comprehensive view of the role of alternative splicing in breast cancer biology and development, we summarize here recent findings regarding alternative splicing events that have been well documented for breast cancer evolution, considering its prognostic and therapeutic value. Moreover, we analyze how the response to endocrine and chemical therapies could be affected due to alternative splicing and differential expression of variant isoforms. With all this knowledge, it becomes clear that targeting alternative splicing represents an innovative approach for breast cancer therapeutics and the information derived from current studies could guide clinical decisions with a direct impact in the clinical advances for breast cancer patients nowadays.

Characterization of BRCA1 and BRCA2 variants found in a Norwegian breast or ovarian cancer cohort

Germline mutations in BRCA1 and BRCA2 cause hereditary breast and ovarian cancer. Molecular screening of these two genes in patients with a family history of breast or ovarian cancer has revealed pathogenic variants as well as genetic variants of unknown significance (VUS). These VUS may cause a challenge in the genetic counseling process regarding clinical management of the patient and the family. We investigated 32 variants previously detected in 33 samples from patients with a family history of breast or ovarian cancer. cDNA was analyzed for alternative transcripts and selected missense variants located in the BRCT domains of BRCA1 were assessed for their trans-activation ability. Although an extensive cDNA analysis was done, only three of the 32 variants appeared to affect the splice-process (BRCA1 c.213-5T>A, BRCA1 c.5434C>G and BRCA2 c.68-7T>A). In addition, two variants located in the BRCT domains of BRCA1 (c.5075A>C p.Asp1692Ala and c.5513T>G p.Val1838Gly) were shown to abolish the BRCT domain trans-activation ability, whereas BRCA1 c.5125G>A p.Gly1709Arg exhibited equal trans-activation capability as the WT domain. These functional studies may offer further insights into the pathogenicity of certain identified variants; however, this assay is only applicable for a subset of missense variants.

Multiplex PCR and NGS-based identification of mRNA splicing variants: Analysis of BRCA1 splicing pattern as a model

Alternative pre-mRNA splicing increases transcriptome plasticity by forming naturally-occurring alternative splicing variants (ASVs). Alterations of splicing processes, caused by DNA mutations, result in aberrant splicing and the formation of aberrant mRNA isoforms. Analyses of hereditary cancer predisposition genes reveal many DNA variants with unknown clinical significance (VUS) that potentially affect pre-mRNA splicing. Therefore, a comprehensive description of ASVs is an essential prerequisite for the interpretation of germline VUS in high-risk individuals. To identify ASVs in a gene of interest, we have proposed an approach based on multiplex PCR (mPCR) amplification of all theoretically possible exon-exon junctions and subsequent characterization of size-selected and pooled mPCR products by next-generation sequencing (NGS). The efficiency of this method is illustrated by a comprehensive analysis of BRCA1 ASVs in human leukocytes, normal mammary, and adipose tissues and stable cell lines. We revealed 94 BRCA1 ASVs, including 29 variants present in all tested samples. While differences in the qualitative expression of BRCA1 ASVs among the analyzed human tissues were minor, larger differences were detected between tissue and cell line samples. Compared with other ASV analysis methods, this approach represents a highly sensitive and rapid alternative for the identification of ASVs in any gene of interest.