The c.5242C>A BRCA1 missense variant induces exon skipping by increasing splicing repressors binding

Understanding and Rescuing the Splicing Defect Caused by the Frequent ABCA4 Variant c.4253+43G>A Underlying Stargardt Disease

Pathogenic variants in ABCA4 are the underlying molecular cause of Stargardt disease (STGD1), an autosomal recessive macular dystrophy characterized by a progressive loss of central vision. Among intronic ABCA4 variants, c.4253+43G>A is frequently detected in STGD1 cases and is classified as a hypomorphic allele, generally associated with late-onset cases. This variant was previously reported to alter splicing regulatory sequences, but the splicing outcome is not fully understood yet. In this study, we attempted to better understand its effect on splicing and to rescue the aberrant splicing via antisense oligonucleotides (AONs). Wild-type and c.4253+43G>A variant-harboring maxigene vectors revealed additional skipping events, which were not previously detected upon transfection in HEK293T cells. To restore exon inclusion, we designed a set of 27 AONs targeting either splicing silencer motifs or the variant region and screened these in maxigene-transfected HEK293T cells. Candidate AONs able to promote exon inclusion were selected for further testing in patient-derived photoreceptor precursor cells. Surprisingly, no robust splicing modulation was observed in this model system. Overall, this research helped to adequately characterize the splicing alteration caused by the c.4253+43G>A variant, although future development of AON-mediated exon inclusion therapy for ABCA4 is needed.

Systematic characterization of short intronic splicing-regulatory elements in SMN2 pre-mRNA

Intronic splicing enhancers and silencers (ISEs and ISSs) are two groups of splicing-regulatory elements (SREs) that play critical roles in determining splice-site selection, particularly for alternatively spliced introns or exons. SREs are often short motifs; their mutation or dysregulation of their cognate proteins frequently causes aberrant splicing and results in disease. To date, however, knowledge about SRE sequences and how they regulate splicing remains limited. Here, using an SMN2 minigene, we generated a complete pentamer-sequence library that comprises all possible combinations of 5 nucleotides in intron 7, at a fixed site downstream of the 5′ splice site. We systematically analyzed the effects of all 1023 mutant pentamers on exon 7 splicing, in comparison to the wild-type minigene, in HEK293 cells. Our data show that the majority of pentamers significantly affect exon 7 splicing: 584 of them are stimulatory and 230 are inhibitory. To identify actual SREs, we utilized a motif set enrichment analysis (MSEA), from which we identified groups of stimulatory and inhibitory SRE motifs. We experimentally validated several strong SREs in SMN1/2 and other minigene settings. Our results provide a valuable resource for understanding how short RNA sequences regulate splicing. Many novel SREs can be explored further to elucidate their mechanism of action.

Knockin mouse model of the human CFL2 p.A35T mutation results in a unique splicing defect and severe myopathy phenotype

Cofilin-2 is an actin-binding protein that is predominantly expressed in skeletal and cardiac muscles and belongs to the AC group of proteins, which includes cofilin-1 and destrin. In humans, cofilin-2 (CFL2) mutations have been associated with congenital myopathies that include nemaline and myofibrillar myopathy. To understand the pathogenicity of the human CFL2 mutation, p.A35T, that first linked cofilin-2 with the human disease, we created a knock-in mouse model. The Cfl2A35T/A35T (KI) mice were indistinguishable from their wild-type littermates at birth, but they rapidly worsened and died by postnatal day 9. The phenotypic, histopathologic and molecular findings mimicked the constitutive Cfl2-knockout (KO) mice described previously, including sarcomeric disruption and actin accumulations in skeletal muscles and negligible amounts of cofilin-2 protein. In addition, KI mice demonstrated a marked reduction in Cfl2 mRNA levels in various tissues including skeletal muscles. Further investigation revealed evidence of alternative splicing with the presence of two alternate transcripts of smaller size. These alternate transcripts were expressed at very low levels in the wild-type mice and were significantly upregulated in the mutant mice, indicating that pre-translational splicing defects may be a critical component of the disease mechanism associated with the mutation. Evidence of reduced expression of the full-length CFL2 transcript was also observed in the muscle biopsy sample of the patient with p.A35T mutation.

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.

Investigating the effect of 28 BRCA1 and BRCA2 mutations on their related transcribed mRNA

Germline inactivating mutations in the BRCA1 and BRCA2 genes are responsible for hereditary breast and ovarian cancer syndrome (HBOCS). Genetic testing of these genes identifies a significant proportion of variants of uncertain significance (VUS). Elucidation of the clinical impact of these variants is an important challenge in genetic diagnostics and counseling. In this study, we assess the RNA effect of 28 BRCA1 and BRCA2 VUS identified in our set of HBOCS families with the aim of gaining insight into their clinical relevance. mRNA was extracted from VUS carriers and controls lymphocytes cultured for 5-6 days and treated with puromycin. RNA was reverse transcribed to perform transcriptional analysis for the study of splicing aberrations. In silico prediction tools were used to select those variants most likely to affect the RNA splicing process. Six out of the 28 variants analyzed showed an aberrant splicing pattern and could therefore be classified as probably pathogenic mutations. Reclassification of VUS improves the genetic counseling and clinical surveillance of carriers of these mutations and highlights the importance of RNA studies in routine diagnostic laboratories.

CAPG and GIPC1: Breast Cancer Biomarkers for Bone Metastasis Development and Treatment

Background:

Bone is the predominant site of metastasis from breast cancer, and recent trials have demonstrated that adjuvant bisphosphonate therapy can reduce bone metastasis development and improve survival. There is an unmet need for prognostic and predictive biomarkers so that therapy can be appropriately targeted.

Methods:

Potential biomarkers for bone metastasis were identified using proteomic comparison of bone-metastatic, lung-metastatic, and nonmetastatic variants of human breast cancer MDA-MB-231 cells. Clinical validation was performed using immunohistochemical staining of tumor tissue microarrays from patients in a large randomized trial of adjuvant zoledronic acid (zoledronate) (AZURE-ISRCTN79831382). We used Cox proportional hazards regression, the Kaplan-Meier estimate of the survival function, and the log-rank test to investigate associations between protein expression, clinical variables, and time to distant recurrence events. All statistical tests were two-sided.

Results:

Two novel biomarker candidates, macrophage-capping protein (CAPG) and PDZ domain–containing protein GIPC1 (GIPC1), were identified for clinical validation. Cox regression analysis of AZURE training and validation sets showed that control patients (no zoledronate) were more likely to develop first distant recurrence in bone (hazard ratio [HR] = 4.5, 95% confidence interval [CI] = 2.1 to 9.8, P < .001) and die (HR for overall survival = 1.8, 95% CI = 1.01 to 3.24, P = .045) if both proteins were highly expressed in the primary tumor. In patients with high expression of both proteins, zoledronate had a substantial effect, leading to 10-fold hazard ratio reduction (compared with control) for first distant recurrence in bone (P = .008).

Conclusions:

The composite biomarker, CAPG and GIPC1 in primary breast tumors, predicted disease outcomes and benefit from zoledronate and may facilitate patient selection for adjuvant bisphosphonate treatment.

Functional classification of BRCA2 DNA variants by splicing assays in a large minigene with 9 exons

Numerous pathogenic DNA variants impair the splicing mechanism in human genetic diseases. Minigenes are optimal approaches to test variants under the splicing viewpoint without the need of patient samples. We aimed to design a robust minigene construct of the breast cancer gene BRCA2 in order to investigate the impact of variants on splicing. BRCA2 exons 19-27 (MGBR2_ex19-27) were cloned in the new vector pSAD. It produced a large transcript of the expected size (2,174 nucleotides) and exon structure (V1-ex19-27-V2). Splicing assays showed that 18 (17 splice-site and 1 silencer variants) out of 40 candidate DNA variants induced aberrant patterns. Twenty-four anomalous transcripts were accurately detected by fluorescent-RT-PCR that were generated by exon-skipping, alternative site usage, and intron-retention events. Fourteen variants induced major anomalies and were predicted to disrupt protein function so they could be classified as pathogenic. Furthermore, minigene mimicked previously reported patient RNA outcomes of seven variants supporting the reproducibility of minigene assays. Therefore, a relevant fraction of variants are involved in breast cancer through splicing alterations. MGBR2_ex19-27 is the largest reported BRCA2 minigene and constitutes a valuable tool for the functional and clinical classification of sequence variations.

Identifying splicing regulatory elements with de Bruijn graphs

Splicing regulatory elements (SREs) are short, degenerate sequences on pre-mRNA molecules that enhance or inhibit the splicing process via the binding of splicing factors, proteins that regulate the functioning of the spliceosome. Existing methods for identifying SREs in a genome are either experimental or computational. Here, we propose a formalism based on de Bruijn graphs that combines genomic structure, word count enrichment analysis, and experimental evidence to identify SREs found in exons. In our approach, SREs are not restricted to a fixed length (i.e., k-mers, for a fixed k). As a result, we identify 2001 putative exonic enhancers and 3080 putative exonic silencers for human genes, with lengths varying from 6 to 15 nucleotides. Many of the predicted SREs overlap with experimentally verified binding sites. Our model provides a novel method to predict variable length putative regulatory elements computationally for further experimental investigation.

DNA repair capacity is impaired in healthy BRCA1 heterozygous mutation carriers

BRCA1 germline mutations increase the lifetime risk of developing breast and ovarian cancers. However, taking into account the differences in disease manifestation among mutation carriers, it is probable that different BRCA1 mutations have distinct haploinsufficiency effects and lead to the formation of different phenotypes. Using lymphoblastoid cell lines derived from heterozygous BRCA1 mutation carriers and non-carriers, we investigated the haploinsufficiency effects of various mutation types using qPCR, immunofluorescence, and microarray technology. Lymphoblastoid cell lines carrying a truncating mutation showed significantly lower BRCA1 mRNA and protein levels and higher levels of gamma-H2AX than control cells or those harboring a missense mutation, indicating greater spontaneous DNA damage. Cells carrying either BRCA1 mutation type showed impaired RAD51 foci formation, suggesting defective repair in mutated cells. Moreover, compared to controls, cell lines carrying missense mutations displayed a more distinct expression profile than cells with truncating mutations, which is consistent with different mutations giving rise to distinct phenotypes. Alterations in the immune response pathway in cells harboring missense mutations point to possible mechanisms of breast cancer initiation in carriers of these mutations. Our findings offer insight into how various heterozygous mutations in BRCA1 could lead to impairment of BRCA1 function and provide strong evidence of haploinsufficiency in BRCA1 mutation carriers.

Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes

Introduction

The underlying pathogenic mechanism of a large fraction of DNA variants of disease-causing genes is the disruption of the splicing process. We aimed to investigate the effect on splicing of the BRCA2 variants c.8488-1G > A (exon 20) and c.9026_9030del (exon 23), as well as 41 BRCA2 variants reported in the Breast Cancer Information Core (BIC) mutation database.

Methods

DNA variants were analyzed with the splicing prediction programs NNSPLICE and Human Splicing Finder. Functional analyses of candidate variants were performed by lymphocyte RT-PCR and/or hybrid minigene assays. Forty-one BIC variants of exons 19, 20, 23 and 24 were bioinformatically selected and generated by PCR-mutagenesis of the wild type minigenes.

Results

Lymphocyte RT-PCR of c.8488-1G > A showed intron 19 retention and a 12-nucleotide deletion in exon 20, whereas c.9026_9030del did not show any splicing anomaly. Minigene analysis of c.8488-1G > A displayed the aforementioned aberrant isoforms but also exon 20 skipping. We further evaluated the splicing outcomes of 41 variants of four BRCA2 exons by minigene analysis. Eighteen variants presented splicing aberrations. Most variants (78.9%) disrupted the natural splice sites, whereas four altered putative enhancers/silencers and had a weak effect. Fluorescent RT-PCR of minigenes accurately detected 14 RNA isoforms generated by cryptic site usage, exon skipping and intron retention events. Fourteen variants showed total splicing disruptions and were predicted to truncate or eliminate essential domains of BRCA2.

Conclusions

A relevant proportion of BRCA2 variants are correlated with splicing disruptions, indicating that RNA analysis is a valuable tool to assess the pathogenicity of a particular DNA change. The minigene system is a straightforward and robust approach to detect variants with an impact on splicing and contributes to a better knowledge of this gene expression step.

Characterisation of unclassified variants in the BRCA1/2 genes with a putative effect on splicing

A subset of the unclassified variants (UVs) identified during genetic screening of BRCA1/2 genes may affect splicing. We assessed at RNA level the effect of four BRCA1 and ten BRCA2 UVs with a putative splice effect, as predicted in silico. The variants selected for this study were beyond the positions -1, -2 or +1, +2 from the exon, and were not previously described (n = 8) or their effect on splicing was not assessed previously (n = 6). Lymphocytes from UV carriers and healthy controls were cultured and treated with puromycin to prevent nonsense-mediated mRNA decay. The relative contribution of each allele to the various transcripts was assessed using combinations of allele-specific and transcript-specific primers. BRCA2 c.425G>T, c.7976+3_7976+4del and c.8754+3G>C give rise to aberrant transcripts BRCA2Δ4, BRCA2Δ17 and retention of 46nt of intron 21, respectively, and were considered pathogenic. BRCA1 c.4987-3C>G gives rise to BRCA1Δ17 that is likely pathogenic; however, residual expression of the full-length transcript from the variant allele could not be excluded. BRCA1 c.692C>T, c.693G>A and BRCA2 c.6935A>T, besides expressing the full-length transcript, increased expression of BRCA1Δ11 and BRCA2Δ12, respectively. As these are naturally occurring isoforms, also observed in controls, the clinical relevance is unclear. The seven remaining UVs did not affect splicing and three intronic variants were therefore classified as neutral. In conclusion, the RNA analysis results clarified the clinical relevance of 6 of the 14 studied UVs and thereby greatly improve the genetic counselling of high-risk breast/ovarian cancer patients carrying these classified variants.

Prediction of single-nucleotide substitutions that result in exon skipping: identification of a splicing silencer in BRCA1 exon 6

Missense, nonsense, and translationally silent mutations can inactivate genes by altering the inclusion of mutant exons in mRNA, but their overall frequency among disease-causing exonic substitutions is unknown. Here, we have tested missense and silent mutations deposited in the BRCA1 mutation databases of unclassified variants for their effects on exon inclusion. Analysis of 21 BRCA1 variants using minigene assays revealed a single exon-skipping mutation c.231G>T. Comprehensive mutagenesis of an adjacent 12-nt segment showed that this silent mutation resulted in a higher level of exon skipping than the 35 other single-nucleotide substitutions. Exon inclusion levels of mutant constructs correlated significantly with predicted splicing enhancers/silencers, prompting the development of two online utilities freely available at http://www.dbass.org.uk. EX-SKIP quickly estimates which allele is more susceptible to exon skipping, whereas HOT-SKIP examines all possible mutations at each exon position and identifies candidate exon-skipping positions/substitutions. We demonstrate that the distribution of exon-skipping and disease-associated substitutions previously identified in coding regions was biased toward top-ranking HOT-SKIP mutations. Finally, we show that proteins 9G8, SC35, SF2/ASF, Tra2, and hnRNP A1 were associated with significant alterations of BRCA1 exon 6 inclusion in the mRNA. Together, these results facilitate prediction of exonic substitutions that reduce exon inclusion in mature transcripts.

Effect of BRCA2 sequence variants predicted to disrupt exonic splice enhancers on BRCA2 transcripts

Background

Genetic screening of breast cancer patients and their families have identified a number of variants of unknown clinical significance in the breast cancer susceptibility genes, BRCA1 and BRCA2. Evaluation of such unclassified variants may be assisted by web-based bioinformatic prediction tools, although accurate prediction of aberrant splicing by unclassified variants affecting exonic splice enhancers (ESEs) remains a challenge.

Methods

This study used a combination of RT-PCR analysis and splicing reporter minigene assays to assess five unclassified variants in the BRCA2 gene that we had previously predicted to disrupt an ESE using bioinformatic approaches.

Results

Analysis of BRCA2 c.8308 G > A (p.Ala2770Thr) by mRNA analysis, and BRCA2 c.8962A > G (p.Ser2988Gly), BRCA2 c.8972G > A (p.Arg2991His), BRCA2 c.9172A > G (p.Ser3058Gly), and BRCA2 c.9213G > T (p.Glu3071Asp) by a minigene assay, revealed no evidence for aberrant splicing.

Conclusions

These results illustrate the need for improved methods for predicting functional ESEs and the potential consequences of sequence variants contained therein.

A glycine-rich domain of hnRNP H/F promotes nucleocytoplasmic shuttling and nuclear import through an interaction with transportin 1

Heterogeneous nuclear ribonucleoprotein (hnRNP) H and F are members of a closely related subfamily of hnRNP proteins that are implicated in many aspects of RNA processing. hnRNP H and F are alternative splicing factors for numerous U2- and U12-dependent introns. The proteins have three RNA binding domains and two glycine-rich domains and localize to both the nucleus and cytoplasm, but little is known about which domains govern subcellular localization or splicing activity. We show here that the central glycine-tyrosine-arginine-rich (GYR) domain is responsible for nuclear localization, and a nonclassical nuclear localization signal (NLS) was mapped to a short, highly conserved sequence whose activity was compromised by point mutations. Glutathione S-transferase (GST) pulldown assays demonstrated that the hnRNP H NLS interacts with the import receptor transportin 1. Finally, we show that hnRNP H/F are transcription-dependent shuttling proteins. Collectively, the results suggest that hnRNP H and F are GYR domain-dependent shuttling proteins whose posttranslational modifications may alter nuclear localization and hence function.

A high proportion of DNA variants of BRCA1 and BRCA2 is associated with aberrant splicing in breast/ovarian cancer patients

PURPOSE

Most BRCA1/2 mutations are of unknown clinical relevance. An increasing amount of evidence indicates that there can be deleterious effects through the disruption of the splicing process. We have investigated the effect of aberrant splicing of BRCA1/2 on hereditary breast/ovarian cancer (HBOC).

EXPERIMENTAL DESIGN

DNA variants were analyzed with splicing prediction programs to select putative splicing mutations. Splicing assays of 57 genetic variants were done by lymphocyte reverse transcription-PCR and/or hybrid minigenes in HeLa and nontumor breast epithelial cells.

RESULTS

Twenty-four BRCA1/2 variants of Spanish HBOC patients were bioinformatically preselected. Functional assays showed that 12 variants induced anomalous splicing patterns, 6 of which accounted for 58.5% of BRCA1 families. To further evaluate the defective splicing of BRCA1/2, we analyzed 31 Breast Cancer Information Core Database (BIC) and two artificial variants that were generated by mutagenesis. Sixteen variants induced different degrees of aberrant splicing. Altogether, anomalous splicing was caused by 28 BRCA1/2 variants of all types, indicating that any DNA change can disrupt pre-mRNA processing. We show that a wide range of regulatory elements can be involved, including the canonical and cryptic splice sites, the polypyrimidine tract, and splicing enhancers/silencers. Twenty mutations were predicted to truncate the BRCA proteins and/or to delete essential domains, thus supporting a role in HBOC.

CONCLUSIONS

An important fraction of DNA variants of BRCA1/2 presents splicing aberrations that may represent a relevant disease-causing mechanism in HBOC. The identification of splicing disruptions by functional assays is a valuable tool to discriminate between benign polymorphisms and pathogenic mutations.