RNA splicing is a key mediator of tumour cell plasticity and a therapeutic vulnerability in colorectal cancer.

Tumour cell plasticity is a major barrier to the efficacy of targeted cancer therapies but the mechanisms that mediate it are poorly understood. Here, we identify dysregulated RNA splicing as a key driver of tumour cell dedifferentiation in colorectal cancer (CRC). We find that Apc-deficient CRC cells have dysregulated RNA splicing machinery and exhibit global rewiring of RNA splicing. We show that the splicing factor SRSF1 controls the plasticity of tumour cells by controlling Kras splicing and is required for CRC invasion in a mouse model of carcinogenesis. SRSF1 expression maintains stemness in human CRC organoids and correlates with cancer stem cell marker expression in human tumours. Crucially, partial genetic downregulation of Srsf1 does not detrimentally affect normal tissue homeostasis, demonstrating that tumour cell plasticity can be differentially targeted. Thus, our findings link dysregulation of the RNA splicing machinery and control of tumour cell plasticity.

A role for endothelial nitric oxide synthase in intestinal stem cell proliferation and mesenchymal colorectal cancer.

BACKGROUND: Nitric oxide (NO) has been highlighted as an important agent in cancer-related events. Although the inducible nitric oxide synthase (iNOS) isoform has received most attention, recent studies in the literature indicate that the endothelial isoenzyme (eNOS) can also modulate different tumor processes including resistance, angiogenesis, invasion, and metastasis. However, the role of eNOS in cancer stem cell (CSC) biology and mesenchymal tumors is unknown. RESULTS: Here, we show that eNOS was significantly upregulated in VilCre ERT2 Apc fl/+ and VilCre ERT2 Apc fl/fl mouse intestinal tissue, with intense immunostaining in hyperproliferative crypts. Similarly, the more invasive VilCre ERT2 Apc fl/+ Pten fl/+ mouse model showed an overexpression of eNOS in intestinal tumors whereas this isoform was not expressed in normal tissue. However, none of the three models showed iNOS expression. Notably, when 40 human colorectal tumors were classified into different clinically relevant molecular subtypes, high eNOS expression was found in the poor relapse-free and overall survival mesenchymal subtype, whereas iNOS was absent. Furthermore, Apc fl/fl organoids overexpressed eNOS compared with wild-type organoids and NO depletion with the scavenger carboxy-PTIO (c-PTIO) decreased the proliferation and the expression of stem-cell markers, such as Lgr5, Troy, Vav3, and Slc14a1, in these intestinal organoids. Moreover, specific NO depletion also decreased the expression of CSC-related proteins in human colorectal cancer cells such as ß-catenin and Bmi1, impairing the CSC phenotype. To rule out the contribution of iNOS in this effect, we established an iNOS-knockdown colorectal cancer cell line. NO-depleted cells showed a decreased capacity to form tumors and c-PTIO treatment in vivo showed an antitumoral effect in a xenograft mouse model. CONCLUSION: Our data support that eNOS upregulation occurs after Apc loss, emerging as an unexpected potential new target in poor-prognosis mesenchymal colorectal tumors, where NO scavenging could represent an interesting therapeutic alternative to targeting the CSC subpopulation.

A study of splicing mutations in disorders of sex development.

The presence of splicing sequence variants in genes responsible for sex development in humans may compromise correct biosynthesis of proteins involved in the normal development of gonads and external genitalia. In a cohort of Brazilian patients, we identified mutations in HSD17B3 and SRD5A2 which are both required for human sexual differentiation. A number of these mutations occurred within regions potentially critical for splicing regulation. Minigenes were used to validate the functional effect of mutations in both genes. We evaluated the c.277 + 2 T > G mutation in HSD17B3, and the c.544 G > A, c.548-44 T > G and c.278delG mutations in SRD5A2. We demonstrated that these mutations altered the splicing pattern of these genes. In a genomic era these results illustrate, and remind us, that sequence variants within exon-intron boundaries, which are primarily identified for diagnostic purposes and have unknown pathogenicity, need to be assessed with regards to their impact not only on protein expression, but also on mRNA splicing.

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

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

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

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

BRCA1 EXON 11, a CERES (composite regulatory element of splicing) element involved in splice regulation.

Unclassified variants (UV) of BRCA1 can affect normal pre-mRNA splicing. Here, we investigate the UV c.693G>A, a "silent" change in BRCA1 exon 11, which we have found induces aberrant splicing in patient carriers and in vitro. Using a minigene assay, we show that the UV c.693G>A has a strong effect on the splicing isoform ratio of BRCA1. Systematic site-directed mutagenesis of the area surrounding the nucleotide position c.693G>A induced variable changes in the level of exon 11 inclusion/exclusion in the mRNA, pointing to the presence of a complex regulatory element with overlapping enhancer and silencer functions. Accordingly, protein binding analysis in the region detected several splicing regulatory factors involved, including SRSF1, SRSF6 and SRSF9, suggesting that this sequence represents a composite regulatory element of splicing (CERES).

BRCA1 exon 11 a model of long exon splicing regulation.

BRCA1 exon 11 is one of the biggest human exons, spanning 3426 bases. This gene is potentially involved in DNA repair as well as cell growth and cell cycle control. Exon 11 is regulated at the splicing level producing three main different combinations of BRCA1 mature transcripts; one including the whole of exon 11 (full isoform), one skipping the entire exon (D11 isoform), and one including only 117 base pairs of exon 11 (D11q isoform). Using minigene and deletion analyses, we have previously described important splicing regulatory sequences located at the beginning of this exon (5' end). We have now found additional important sequences located at its 3' end. In particular, we describe the presence of a strong splicing enhancer adjacent to the downstream 5' splice site, which minimizes competition from an upstream 5' splice site and so ensures long exon inclusion. Analyses of the proteins binding these RNA sequences have revealed that Tra2beta and hnRNP L are involved in the regulation of BRCA1 exon 11 by influencing the recognition of donor sites. Interestingly, BRCA1 exon 11 carrying deletion of the regulatory sequences bound by these factors also showed unexpected responses to up- or downregulation of these regulatory proteins, suggesting that they can also bind elsewhere in this large exon and elicit different effects on its recognition.   The identification of sequences and proteins relevant for the regulation of BRCA1 exon 11 now provides better knowledge on how this exon is recognized and may represent an important step toward understanding how large exons are regulated.

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

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

BRCA1 exon 11 alternative splicing, multiple functions and the association with cancer.

BRCA1 (breast cancer early-onset 1) alternative splicing levels are regulated in a cell-cycle- and cell-type-specific manner, with splice variants being present in different proportions in tumour cell lines as well as in normal mammary epithelial cells. The importance of this difference in the pathogenesis of breast cancer has yet to be determined. Developing an understanding of the impact of BRCA1 isoform ratio changes on cell phenotype will be of value in breast cancer and may offer therapeutic options. In the present paper, we describe the splicing isoforms of BRCA1 exon 11, their possible role in cancer biology and the importance of maintaining a balanced ratio.

Evolutionary constraint helps unmask a splicing regulatory region in BRCA1 exon 11.

BACKGROUND: Alternative splicing across exon 11 produces several BRCA1 isoforms. Their proportion varies during the cell cycle, between tissues and in cancer suggesting functional importance of BRCA1 splicing regulation around this exon. Although the regulatory elements driving exon 11 splicing have never been identified, a selective constraint against synonymous substitutions (silent nucleotide variations that do not alter the amino acid residue sequence) in a critical region of BRCA1 exon 11 has been reported to be associated with the necessity to maintain regulatory sequences. METHODOLOGY/PRINCIPAL FINDINGS: Here we have designed a specific minigene to investigate the possibility that this bias in synonymous codon usage reflects the need to preserve the BRCA1 alternative splicing program. We report that in-frame deletions and translationally silent nucleotide substitutions in the critical region affect splicing regulation of BRCA1 exon 11. CONCLUSIONS/SIGNIFICANCE: Using a hybrid minigene approach, we have experimentally validated the hypothesis that the need to maintain correct alternative splicing is a selective pressure against translationally silent sequence variations in the critical region of BRCA1 exon 11. Identification of the trans-acting factors involved in regulating exon 11 alternative splicing will be important in understanding BRCA1-associated tumorigenesis.

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.

Low U1 snRNP dependence at the NF1 exon 29 donor splice site.

Many disease-causing splicing mutations described in the literature produce changes in splice sites (SS) or in exon-regulatory sequences. The delineation of these splice aberrations can provide important insights into novel regulation mechanisms. In this study, we evaluated the effect of patient variations in neurofibromatosis type 1 (NF1) exon 29 and its 5'SS surrounding area on its splicing process. Only two of all nonsense, missense, synonymous and intronic variations analyzed in this study clearly altered exon 29 inclusion/exclusion levels. In particular, the intronic mutation +5g>a had the strongest effect, resulting in total exon exclusion. This finding prompted us to evaluate the exon 29 5'SS in relation to its ability to bind U1 snRNP. This was performed by direct analysis of the ability of U1 to bind to wild-type and mutant donor sites, by engineering an in vitro splicing system to directly evaluate the functional importance of U1 snRNA base pairing with the exon 29 donor site, and by coexpression of mutant U1 snRNP molecules to try to rescue exon 29 inclusion in vivo. The results revealed a low dependency on the presence of U1 snRNP, and suggest that exon 29 donor site definition may depend on alternative mechanisms of 5'SS recognition.

Polypyrimidine tract binding protein regulates alternative splicing of an aberrant pseudoexon in NF1.

In disease-associated genes, understanding the functional significance of deep intronic nucleotide variants represents a difficult challenge. We previously reported that an NF1 intron 30 exonization event is triggered from a single correct nomenclature is 'c.293-279 A>G' mutation [Raponi M, Upadhyaya M & Baralle D (2006) Hum Mutat 27, 294-295]. In this paper, we investigate which characteristics play a role in regulating inclusion of the aberrant pseudoexon. Our investigation shows that pseudoexon inclusion levels are strongly downregulated by polypyrimidine tract binding protein and its homologue neuronal polypyrimidine tract binding protein. In particular, we provide evidence that the functional effect of polypyrimidine tract binding protein is proportional to its concentration, and map the cis-acting elements that are principally responsible for this negative regulation. These results highlight the importance of evaluating local sequence context for diagnostic purposes, and the utility of developing therapies to turn off activated pseudoexons.

Can donor splice site recognition occur without the involvement of U1 snRNP?

Many disease-causing mutations affecting donor splice site recognition are reported in the literature. One of the more frequently observed nucleotide changes causing aberrant splicing are due to mutations in the donor splice site which lower the strength of base pairing with U1 snRNA (small nuclear RNA). However, recent data have highlighted the possibility of a recognition mechanism for weak donor splice sites that are at least partially U1-independent. This is important as most of the donor splice site prediction programs currently in use are based on the U1 snRNA 5'-splice site base pairing and would not pick this up. We review these mechanisms and how an up-to-date donor splice site mutation repertoire indicates the heterogeneity of the molecular mechanism. We suggest that, in clinical molecular genetics, it is important to evaluate sequence variants for aberrant splicing even in those cases where the variant is not thought to alter the U1 snRNA interaction.

Reduced splicing efficiency induced by synonymous substitutions may generate a substrate for natural selection of new splicing isoforms: the case of CFTR exon 12.

Alternative splicing has been associated with increased evolutionary changes and with recent exon creation or loss. The addition of a new exon can be explained by its inclusion in only a fraction of the transcripts leaving the original form intact and giving to the new form the possibility to evolve independently but the exon loss phenomenon is less clear. To explore the mechanism that could be involved in CFTR exon 12 lower splicing efficiency in primates, we have analyzed the effect of multiple synonymous variations. Random patterns of synonymous variations were created in CFTR exon12 and the majority of them induced exon inclusion, suggesting a suboptimal splicing efficiency of the human gene. In addition, the effect of each single synonymous substitution on splicing is strongly dependent on the exonic context and does not correlate with available in silico exon splicing prediction programs. We propose that casual synonymous substitutions may lead to a reduced splicing efficiency that can result in a variable proportion of exon loss. If this phenomenon happens in in-frame exons and to an extent tolerated by the cells it can have an important evolutionary effect since it may generate a substrate for natural selection of new splicing isoforms.

Synonymous mutations in CFTR exon 12 affect splicing and are not neutral in evolution.

It is well established that exonic sequences contain regulatory elements of splicing that overlap with coding capacity. However, the conflict between ensuring splicing efficiency and preserving the coding capacity for an optimal protein during evolution has not been specifically analyzed. In fact, studies on genomic variability in fields as diverse as clinical genetics and molecular evolution mainly focus on the effect of mutations on protein function. Synonymous variations, in particular, are assumed to be functionally neutral both in clinical diagnosis and when measuring evolutionary distances between species. Using the cystic fibrosis transmembrane conductance regulator (CFTR) exon 12 splicing as a model, we have established that about one quarter of synonymous variations result in exon skipping and, hence, in an inactive CFTR protein. Furthermore, comparative splicing evaluation of mammalian sequence divergences showed that artificial combinations of CFTR exon 12 synonymous and nonsynonymous substitutions are incompatible with normal RNA processing. In particular, the combination of the mouse synonymous with the human missense variations causes exon skipping. It follows that there are two sequential levels at which evolutionary selection of genomic variants take place: splicing control and protein function optimization.

A centrosomal mechanism involving CDK5RAP2 and CENPJ controls brain size.

Autosomal recessive primary microcephaly is a potential model in which to research genes involved in human brain growth. We show that two forms of the disorder result from homozygous mutations in the genes CDK5RAP2 and CENPJ. We found neuroepithelial expression of the genes during prenatal neurogenesis and protein localization to the spindle poles of mitotic cells, suggesting that a centrosomal mechanism controls neuron number in the developing mammalian brain.