Site directed mutagenesis of hMLH1 exonic splicing enhancers does not correlate with splicing disruption

Splice-disrupt genomic variants in prostate cancer

Background

Splice-disrupt genomic variants are one of the causes of cancer-causing errors in gene expression. Little is known about splice-disrupt genomic variants.

Methods and results

Here, pattern of splice-disrupt variants was investigated using 21,842,764 genomic variants in different types of prostate cancer. A particular attention was paid to genomic locations of splice-disrupt variants on target genes. HLA-A in prostate cancer, MSR1 in familial prostate cancer, and EGFR in both castration-resistant prostate cancer and metastatic castration-resistant had the highest allele frequencies of splice-disrupt variations. Some splice-disrupt variants, located on coding sequences of NCOR2, PTPRC, and CRP, were solely present in the advanced metastatic castration-resistant prostate cancer. High-risk splice-disrupt variants were identified based on computationally calculated Polymorphism Phenotyping (PolyPhen), Sorting Intolerant From Tolerant (SIFT), and Genomic Evolutionary Rate Profiling (GERP) + + scores as well as the recorded clinical significance in dbSNP database of NCBI. Functional annotation of damaging splice-disrupt variants highlighted important cancer-associated functions, including endocrine resistance, lipid metabolic process, steroid metabolic process, regulation of mitotic cell cycle, and regulation of metabolic process. This is the first study that profiles the splice-disrupt genomic variants and their target genes in prostate cancer. Literature mining based variant analysis highlighted the importance of rs1800716 variant, located on the CYP2D6 gene, involved in a range of important functions, such as RNA spicing, drug interaction, death, and urotoxicity.

Conclusions

This is the first study that profiles the splice-disrupt genomic variants and their target genes in different types of prostate cancer. Unravelling alternative splicing opens a new avenue towards the establishment of new diagnostic and prognostic markers for prostate cancer progression and metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11033-022-07257-9.

PAX7 nucleotide variants and the risk of non-syndromic orofacial clefts in the Polish population

OBJECTIVE

The etiology of non-syndromic cleft lip with or without cleft palate (nsCL/P) is multifactorial, heterogeneous, and still not completely understood. The aim of the present study was to examine the associations between common and rare PAX7 nucleotide variants and the risk of this common congenital anomaly in a Polish population.

SUBJECTS AND METHODS

Eight top nsCL/P-associated PAX7 variants identified in our cleft genome-wide association study (GWAS) were selected for replication analysis in an independent group of patients and controls (n = 247 and n = 445, respectively). In addition, mutation screening of the PAX7 protein-coding region was conducted.

RESULTS

Analysis of the pooled data from the GWAS and replication study confirmed that common PAX7 nucleotide variants are significantly associated with the increased risk of nsCL/P. The strongest individual variant was rs1339062 (c.586 + 15617T > C) with a p-value = 2.47E-05 (OR = 1.4, 95%CI: 1.20-1.64). Sequencing analysis identified a novel synonymous PAX7 substitution (c.87G > A, p.Val29Val) in a single patient with nsCLP. This transition located in the early exonic position was predicted to disrupt potential splice enhancer elements.

CONCLUSION

Our study confirmed that PAX7 is a strong candidate gene for nsCL/P. Nucleotide variants of this gene contribute to the etiology of nsCL/P in the homogenous Polish population.

RNA splicing in human disease and in the clinic

Defects at the level of the pre-mRNA splicing process represent a major cause of human disease. Approximately 15-50% of all human disease mutations have been shown to alter functioning of basic and auxiliary splicing elements. These elements are required to ensure proper processing of pre-mRNA splicing molecules, with their disruption leading to misprocessing of the pre-mRNA molecule and disease. The splicing process is a complex process, with much still to be uncovered before we are able to accurately predict whether a reported genomic sequence variant (GV) represents a splicing-associated disease mutation or a harmless polymorphism. Furthermore, even when a mutation is correctly identified as affecting the splicing process, there still remains the difficulty of providing an exact evaluation of the potential impact on disease onset, severity and duration. In this review, we provide a brief overview of splicing diagnostic methodologies, from in silico bioinformatics approaches to wet lab in vitro and in vivo systems to evaluate splicing efficiencies. In particular, we provide an overview of how the latest developments in high-throughput sequencing can be applied to the clinic, and are already changing clinical approaches.

Biallelic DICER1 mutations occur in Wilms tumours

DICER1 is an endoribonuclease central to the generation of microRNAs (miRNAs) and short interfering RNAs (siRNAs). Germline mutations in DICER1 have been associated with a pleiotropic tumour predisposition syndrome and Wilms tumour (WT) is a rare manifestation of this syndrome. Three WTs, each in a child with a deleterious germline DICER1 mutation, were screened for somatic DICER1 mutations and were found to bear specific mutations in either the RNase IIIa (n = 1) or the RNase IIIb domain (n = 2). In the two latter cases, we demonstrate that the germline and somatic DICER1 mutations were in trans, suggesting that the two-hit hypothesis of tumour formation applies for these examples of WT. Among 191 apparently sporadic WTs, we identified five different missense or deletion somatic DICER1 mutations (2.6%) in four individual WTs; one tumour had two very likely deleterious somatic mutations in trans in the RNase IIIb domain (c.5438A>G and c.5452G>A). In vitro studies of two somatic single-base substitutions (c.5429A>G and c.5438A>G) demonstrated exon 25 skipping from the transcript, a phenomenon not previously reported in DICER1. Further we show that DICER1 transcripts lacking exon 25 can be translated in vitro. This study has demonstrated that a subset of WTs exhibits two 'hits' in DICER1, suggesting that these mutations could be key events in the pathogenesis of these tumours.

Aberrant splicing caused by a MLH1 splice donor site mutation found in a young Japanese patient with Lynch syndrome

Lynch syndrome, also known as hereditary non-polyposis colorectal cancer, characterized by predisposition to colorectal cancer and other associated cancers, is an autosomal-dominant disorder mainly caused by germline mutations in DNA mismatch repair (MMR) genes such as MLH1, MSH2, and MSH6. Some mutations that disrupt splice donor or acceptor sites cause aberrant mRNA splicing. These mutations are generally considered as pathogenic ones, however, it is sometimes uneasy to accurately predict their pathogenicity without functional assays, particularly when the mutation is a single nucleotide substitution. In this report, we describe a 25-year-old patient with Lynch syndrome who carries a germline variant in a splice donor site of the MLH1 gene (c.790 + 5 G > T), which was first detected among Asian populations. The immunohistochemical analysis revealed loss of MLH1 protein expression in the tumor. Our splicing assay confirmed that the intronic MLH1 variant actually caused aberrant splicing, supporting its pathogenic effect. Our data accumulate more information on the genotype-phenotype relationships in patients with Lynch syndrome.

Integrated analysis of unclassified variants in mismatch repair genes

PURPOSE

Lynch syndrome is a genetic disease that predisposes to colorectal tumors, caused by mutation in mismatch repair genes. The use of genetic tests to identify mutation carriers does not always give perfectly clear results, as happens when an unclassified variant is found. This study aimed to define the pathogenic role of 35 variants present in MSH2, MLH1, MSH6, and PMS2 genes identified in our 15-year case study.

METHODS

We collected clinical and molecular data of all carriers, and then we analyzed the variants pathogenic role with web tools and molecular analyses. Using a Bayesian approach, we derived a posterior probability of pathogenicity and classified each variant according to a standardized five-class system.

RESULTS

The MSH2 p.Pro349Arg, p.Met688Arg, the MLH1 p.Gly67Arg, p.Thr82Ala, p.Lys618Ala, the MSH6 p.Ala1236Pro, and the PMS2 p.Arg20Gln were classified as pathogenic, and the MSH2 p.Cys697Arg and the PMS2 p.Ser46Ile were classified as likely pathogenic. Seven variants were likely nonpathogenic, 3 were nonpathogenic, and 16 remained uncertain.

CONCLUSION

Quantitative assessment of several parameters and their integration in a multifactorial likelihood model is the method of choice for classifying the variants. As such classifications can be associated with surveillance and testing recommendations, the results and the method developed in our study can be useful for helping laboratory geneticists in evaluation of genetic tests and clinicians in the management of carriers.

Expression and regulation of a low-density lipoprotein receptor exon 12 splice variant

As low-density lipoprotein receptor (LDLR) contributes to cholesterol and amyloid beta homeostasis, insights into LDLR regulation may facilitate our understanding of cardiovascular disease and Alzheimer's disease. Previously, we identified LDLR isoforms that lacked exon 12 or exons 11-12 and that are predicted to encode soluble, dominant negative, LDLR. Moreover, these isoforms were associated with rs688, an exon 12 polymorphism that was associated with LDL-cholesterol and Alzheimer's disease risk. In this study, we present evidence that although the truncated LDLR isoforms are translated in vitro, they represent < 0.1% of CSF proteins. As these LDLR isoforms likely represent a loss of mRNA-encoding functional LDLR, we then focused upon identifying intron-exon boundary and exonic splicing enhancer elements critical to splicing. Exon 12 inclusion is enhanced by altering the 5' splice site in intron 12 towards a consensus splice donor sequence, consistent with its being a weak 5' splice site. Additionally, of the nine evolutionarily conserved putative splicing enhancer regions within exon 12, two regions that flank rs688 were critical to exon 12 inclusion. Overall, these results suggest that LDLR splice variants represent a loss of mRNA encoding functional LDLR and provide insights into the regulatory elements critical for LDLR exon 12 splicing.

Genomic features defining exonic variants that modulate splicing

A comparative analysis of SNPs and their exonic and intronic environments identifies the features predictive of splice affecting variants.

Background

Single point mutations at both synonymous and non-synonymous positions within exons can have severe effects on gene function through disruption of splicing. Predicting these mutations in silico purely from the genomic sequence is difficult due to an incomplete understanding of the multiple factors that may be responsible. In addition, little is known about which computational prediction approaches, such as those involving exonic splicing enhancers and exonic splicing silencers, are most informative.

Results

We assessed the features of single-nucleotide genomic variants verified to cause exon skipping and compared them to a large set of coding SNPs common in the human population, which are likely to have no effect on splicing. Our findings implicate a number of features important for their ability to discriminate splice-affecting variants, including the naturally occurring density of exonic splicing enhancers and exonic splicing silencers of the exon and intronic environment, extensive changes in the number of predicted exonic splicing enhancers and exonic splicing silencers, proximity to the splice junctions and evolutionary constraint of the region surrounding the variant. By extending this approach to additional datasets, we also identified relevant features of variants that cause increased exon inclusion and ectopic splice site activation.

Conclusions

We identified a number of features that have statistically significant representation among exonic variants that modulate splicing. These analyses highlight putative mechanisms responsible for splicing outcome and emphasize the role of features important for exon definition. We developed a web-tool, Skippy, to score coding variants for these relevant splice-modulating features.

Functional characterization of rare missense mutations in MLH1 and MSH2 identified in Danish colorectal cancer patients

Recently, we have performed a population based study to analyse the frequency of colorectal cancer related MLH1 and MSH2 missense mutations in the Danish population. Half of the analyzed mutations were rare and most likely only present in the families where they were identified originally. Some of the missense mutations were located in conserved regions in the MLH1 and MSH2 proteins indicating a relation to disease development. In the present study, we functionally characterized 10 rare missense mutations in MLH1 and MSH2 identified in 13 Danish CRC families. To elucidate the pathogenicity of the missense mutations, we carried out in vitro functional analyses. The missense mutations were analyzed for their effect on protein expression and repair efficiency. The results of the functional analysis were correlated with clinical data on the families carrying these mutations. Eight missense mutations resulted in proteins with expression and repair efficiency similar to the wild type. One missense mutation (MSH2 p.Met688Val) caused reduced protein expression and one (MSH2 p.Leu187Arg) caused both reduced protein expression and repair deficiency. The MSH2 p.Leu187Arg mutation was found in an Amsterdam II family presenting with high microsatellite instability and loss of MSH2 and MSH6 proteins in tumours. In conclusion, only 1/10 missense mutations displayed repair deficiency and could be classified as pathogenic. No final conclusion can be drawn on the MSH2 p.Met688Val mutation, which caused reduced protein expression. Although, no deficiencies have been identified in the proteins harbouring the other missense mutations, pathogenicity of these variants cannot be unambiguously excluded.

Determination of splice-site mutations in Lynch syndrome (hereditary non-polyposis colorectal cancer) patients using functional splicing assay

Lynch syndrome (hereditary non-polyposis colorectal cancer) is an inherited disease caused by germ-line mutation in mismatch repair genes such as MLH1, MSH2, and MSH6. The mutations include missense and nonsense mutations, small insertions and deletions, and gross genetic alterations including large deletions and duplications. In addition to these genetic changes, mutations in introns are also involved in the pathogenesis. However, it is sometimes difficult to interpret correctly the pathogenicity of variants in exons as well as introns. To evaluate the effect of splice-site mutations in two Lynch syndrome patients, we carried out a functional splicing assay using minigenes. Consequently, this assay showed that the mutation of c.1731+5G>A in MLH1 led to exon15 skipping, and that the mutation of c.211+1G>C in MSH2 created an activated cryptic splice-site 17-nucleotides upstream in exon1. These aberrant splicing patterns were not observed when wild type sequence was used for the assay. We also obtained concordant results by RT-PCR experiments with transcripts from the patients. Furthermore, additional functional splicing assays using two different intronic mutations described in earlier studies revealed splicing alterations that were in complete agreement with the reports. Therefore, functional splicing assay is helpful for evaluating the effects of genetic variants on splicing.

A large fraction of unclassified variants of the mismatch repair genes MLH1 and MSH2 is associated with splicing defects

Numerous unclassified variants (UVs) have been found in the mismatch repair genes MLH1 and MSH2 involved in hereditary nonpolyposis colorectal cancer (HNPCC or Lynch syndrome). Some of these variants may have an effect on pre-mRNA splicing, either by altering degenerate positions of splice site sequences or by affecting intronic or exonic splicing regulatory sequences such as exonic splicing enhancers (ESEs). In order to determine the consequences of UVs on splicing, we used a functional assay of exon inclusion. For each variant, mutant and wild-type exons to be tested were PCR-amplified from patient genomic DNA together with approximately 150 bp of flanking sequences and were inserted into a splicing reporter minigene. After transfection into HeLa cells, the effects on splicing were evaluated by RT-PCR analysis and systematic sequencing. A total of 22 UVs out of 85 different variant alleles examined in 82 families affected splicing, including four exonic variants that affected putative splicing regulatory elements. We analyzed short stretches spanning the latter variants by cloning them into the ESE-dependent central exon of a three-exon splicing minigene and we showed in cell transfection experiments that the wild-type sequences indeed contain functional ESEs. We then used this construct to query for ESE elements in the MLH1 or MSH2 regions affected by 14 previously reported exonic splicing mutations and showed that they also contain functional ESEs. These splicing assays represent a valuable tool for the interpretation of UVs and should contribute to the optimization of the molecular diagnosis of the Lynch syndrome and of other genetic diseases.

In silico and in vivo splicing analysis of MLH1 and MSH2 missense mutations shows exon- and tissue-specific effects

BACKGROUND

Abnormalities of pre-mRNA splicing are increasingly recognized as an important mechanism through which gene mutations cause disease. However, apart from the mutations in the donor and acceptor sites, the effects on splicing of other sequence variations are difficult to predict. Loosely defined exonic and intronic sequences have been shown to affect splicing efficiency by means of silencing and enhancement mechanisms. Thus, nucleotide substitutions in these sequences can induce aberrant splicing. Web-based resources have recently been developed to facilitate the identification of nucleotide changes that could alter splicing. However, computer predictions do not always correlate with in vivo splicing defects. The issue of unclassified variants in cancer predisposing genes is very important both for the correct ascertainment of cancer risk and for the understanding of the basic mechanisms of cancer gene function and regulation. Therefore we aimed to verify how predictions that can be drawn from in silico analysis correlate with results obtained in an in vivo splicing assay.

RESULTS

We analysed 99 hMLH1 and hMSH2 missense mutations with six different algorithms. Transfection of three different cell lines with 20 missense mutations, showed that a minority of them lead to defective splicing. Moreover, we observed that some exons and some mutations show cell-specific differences in the frequency of exon inclusion.

CONCLUSION

Our results suggest that the available algorithms, while potentially helpful in identifying splicing modulators especially when they are located in weakly defined exons, do not always correspond to an obvious modification of the splicing pattern. Thus caution must be used in assessing the pathogenicity of a missense or silent mutation with prediction programs. The variations observed in the splicing proficiency in three different cell lines suggest that nucleotide changes may dictate alternative splice site selection in a tissue-specific manner contributing to the widely observed phenotypic variability in inherited cancers.

Aberrant splicing in MLH1 and MSH2 due to exonic and intronic variants

Single base substitutions in DNA mismatch repair genes which are predicted to lead either to missense or silent mutations, or to intronic variants outside the highly conserved splicing region are often found in hereditary nonpolyposis colorectal cancer (HNPCC) families. In order to use the variants for predictive testing in persons at risk, their pathogenicity has to be evaluated. There is growing evidence that some substitutions have a detrimental influence on splicing. We examined 19 unclassified variants (UVs) detected in MSH2 or MLH1 genes in patients suspected of HNPCC for expression at RNA level. We demonstrate that 10 of the 19 UVs analyzed affect splicing. For example, the substitution MLH1,c.2103G > C in the last position of exon 18 does not result in a missense mutation as theoretically predicted (p.Gln701His), but leads to a complete loss of exon 18. The substitution MLH1,c.1038G > C (predicted effect p.Gln346His) leads to complete inactivation of the mutant allele by skipping of exons 10 and 11, and by activation of a cryptic intronic splice site. Similarly, the intronic variant MLH1,c.306+2dupT results in loss of exon 3 and a frameshift mutation due to a new splice donor site 5 bp upstream. Furthermore, we confirmed complete exon skipping for the mutations MLH1,c.1731G > A and MLH1,c.677G > A. Partial exon skipping was demonstrated for the mutations MSH2,c.1275A > G, MLH1,c.588+5G > A, MLH1,c.790+4A > G and MLH1,c.1984A > C. In contrast, five missense mutations (MSH2,c.4G > A, MSH2,c.2123T > A, MLH1,c.464T > G, MLH1,c.875T > C and MLH1,c.2210A > T) were found in similar proportions in the mRNA as in the genomic DNA. We conclude that the mRNA examination should precede functional tests at protein level.

Disruption of an exon splicing enhancer in exon 3 of MLH1 is the cause of HNPCC in a Quebec family

BACKGROUND

A 3 bp deletion located at the 5' end of exon 3 of MLH1, resulting in deletion of exon 3 from RNA, was recently identified.

HYPOTHESIS

That this mutation disrupts an exon splicing enhancer (ESE) because it occurs in a purine-rich sequence previously identified as an ESE in other genes, and ESEs are often found in exons with splice signals that deviate from the consensus signals, as does the 3' splice signal in exon 3 of MLH1.

DESIGN

The 3 bp deletion and several other mutations were created by polymerase chain reaction mutagenesis and tested using an in vitro splicing assay. Both mutant and wild type exon 3 sequences were cloned into an exon trapping vector and transiently expressed in Cos-1 cells.

RESULTS

Analysis of the RNA indicates that the 3 bp deletion c.213_215delAGA (gi:28559089, NM_000249.2), a silent mutation c.216T-->C, a missense mutation c.214G-->C, and a nonsense mutation c.214G-->T all cause varying degrees of exon skipping, suggesting the presence of an ESE at the 5' end of exon 3. These mutations are situated in a GAAGAT sequence 3 bp downstream from the start of exon 3.

CONCLUSIONS

The results of the splicing assay suggest that inclusion of exon 3 in the mRNA is ESE dependent. The exon 3 ESE is not recognised by all available motif scoring matrices, highlighting the importance of RNA analysis in the detection of ESE disrupting mutations.