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.

SPiP: Splicing Prediction Pipeline, a machine learning tool for massive detection of exonic and intronic variant effects on mRNA splicing.

Modeling splicing is essential for tackling the challenge of variant interpretation as each nucleotide variation can be pathogenic by affecting pre-mRNA splicing via disruption/creation of splicing motifs such as 5'/3' splice sites, branch sites, or splicing regulatory elements. Unfortunately, most in silico tools focus on a specific type of splicing motif, which is why we developed the Splicing Prediction Pipeline (SPiP) to perform, in one single bioinformatic analysis based on a machine learning approach, a comprehensive assessment of the variant effect on different splicing motifs. We gathered a curated set of 4616 variants scattered all along the sequence of 227 genes, with their corresponding splicing studies. The Bayesian analysis provided us with the number of control variants, that is, variants without impact on splicing, to mimic the deluge of variants from high-throughput sequencing data. Results show that SPiP can deal with the diversity of splicing alterations, with 83.13% sensitivity and 99% specificity to detect spliceogenic variants. Overall performance as measured by area under the receiving operator curve was 0.986, better than SpliceAI and SQUIRLS (0.965 and 0.766) for the same data set. SPiP lends itself to a unique suite for comprehensive prediction of spliceogenicity in the genomic medicine era. SPiP is available at: https://sourceforge.net/projects/splicing-prediction-pipeline/.

Functional characterization of ABCC8 variants of unknown significance based on bioinformatics predictions, splicing assays, and protein analyses: Benefits for the accurate diagnosis of congenital hyperinsulinism.

ABCC8 encodes the SUR1 subunit of the ß-cell ATP-sensitive potassium channel whose loss of function causes congenital hyperinsulinism (CHI). Molecular diagnosis is critical for optimal management of CHI patients. Unfortunately, assessing the impact of ABCC8 variants on RNA splicing remains very challenging as this gene is poorly expressed in leukocytes. Here, we performed bioinformatics analysis and cell-based minigene assays to assess the impact on splicing of 13 ABCC8 variants identified in 20 CHI patients. Next, channel properties of SUR1 proteins expected to originate from minigene-detected in-frame splicing defects were analyzed after ectopic expression in COSm6 cells. Out of the analyzed variants, seven induced out-of-frame splicing defects and were therefore classified as recessive pathogenic, whereas two led to skipping of in-frame exons. Channel functional analysis of the latter demonstrated their pathogenicity. Interestingly, the common rs757110 SNP increased exon skipping in our system suggesting that it may act as a disease modifier factor. Our strategy allowed determining the pathogenicity of all selected ABCC8 variants, and CHI-inheritance pattern for 16 out of the 20 patients. This study highlights the value of combining RNA and protein functional approaches in variant interpretation and reveals the minigene splicing assay as a new tool for CHI molecular diagnostics.

Large-scale comparative evaluation of user-friendly tools for predicting variant-induced alterations of splicing regulatory elements.

Discriminating which nucleotide variants cause disease or contribute to phenotypic traits remains a major challenge in human genetics. In theory, any intragenic variant can potentially affect RNA splicing by altering splicing regulatory elements (SREs). However, these alterations are often ignored mainly because pioneer SRE predictors have proved inefficient. Here, we report the first large-scale comparative evaluation of four user-friendly SRE-dedicated algorithms (QUEPASA, HEXplorer, SPANR, and HAL) tested both as standalone tools and in multiple combined ways based on two independent benchmark datasets adding up to >1,300 exonic variants studied at the messenger RNA level and mapping to 89 different disease-causing genes. These methods display good predictive power, based on decision thresholds derived from the receiver operating characteristics curve analyses, with QUEPASA and HAL having the best accuracies either as standalone or in combination. Still, overall there was a tight race between the four predictors, suggesting that all methods may be of use. Additionally, QUEPASA and HEXplorer may be beneficial as well for predicting variant-induced creation of pseudoexons deep within introns. Our study highlights the potential of SRE predictors as filtering tools for identifying disease-causing candidates among the plethora of variants detected by high-throughput DNA sequencing and provides guidance for their use in genomic medicine settings.

Assessment of branch point prediction tools to predict physiological branch points and their alteration by variants.

BACKGROUND: Branch points (BPs) map within short motifs upstream of acceptor splice sites (3'ss) and are essential for splicing of pre-mature mRNA. Several BP-dedicated bioinformatics tools, including HSF, SVM-BPfinder, BPP, Branchpointer, LaBranchoR and RNABPS were developed during the last decade. Here, we evaluated their capability to detect the position of BPs, and also to predict the impact on splicing of variants occurring upstream of 3'ss. RESULTS: We used a large set of constitutive and alternative human 3'ss collected from Ensembl (n = 264,787 3'ss) and from in-house RNAseq experiments (n = 51,986 3'ss). We also gathered an unprecedented collection of functional splicing data for 120 variants (62 unpublished) occurring in BP areas of disease-causing genes. Branchpointer showed the best performance to detect the relevant BPs upstream of constitutive and alternative 3'ss (99.48 and 65.84% accuracies, respectively). For variants occurring in a BP area, BPP emerged as having the best performance to predict effects on mRNA splicing, with an accuracy of 89.17%. CONCLUSIONS: Our investigations revealed that Branchpointer was optimal to detect BPs upstream of 3'ss, and that BPP was most relevant to predict splicing alteration due to variants in the BP area.

Novel diagnostic tool for prediction of variant spliceogenicity derived from a set of 395 combined in silico/in vitro studies: an international collaborative effort.

Variant interpretation is the key issue in molecular diagnosis. Spliceogenic variants exemplify this issue as each nucleotide variant can be deleterious via disruption or creation of splice site consensus sequences. Consequently, reliable in silico prediction of variant spliceogenicity would be a major improvement. Thanks to an international effort, a set of 395 variants studied at the mRNA level and occurring in 5' and 3' consensus regions (defined as the 11 and 14 bases surrounding the exon/intron junction, respectively) was collected for 11 different genes, including BRCA1, BRCA2, CFTR and RHD, and used to train and validate a new prediction protocol named Splicing Prediction in Consensus Elements (SPiCE). SPiCE combines in silico predictions from SpliceSiteFinder-like and MaxEntScan and uses logistic regression to define optimal decision thresholds. It revealed an unprecedented sensitivity and specificity of 99.5 and 95.2%, respectively, and the impact on splicing was correctly predicted for 98.8% of variants. We therefore propose SPiCE as the new tool for predicting variant spliceogenicity. It could be easily implemented in any diagnostic laboratory as a routine decision making tool to help geneticists to face the deluge of variants in the next-generation sequencing era. SPiCE is accessible at (https://sourceforge.net/projects/spicev2-1/).

Exonic Splicing Mutations Are More Prevalent than Currently Estimated and Can Be Predicted by Using In Silico Tools.

The identification of a causal mutation is essential for molecular diagnosis and clinical management of many genetic disorders. However, even if next-generation exome sequencing has greatly improved the detection of nucleotide changes, the biological interpretation of most exonic variants remains challenging. Moreover, particular attention is typically given to protein-coding changes often neglecting the potential impact of exonic variants on RNA splicing. Here, we used the exon 10 of MLH1, a gene implicated in hereditary cancer, as a model system to assess the prevalence of RNA splicing mutations among all single-nucleotide variants identified in a given exon. We performed comprehensive minigene assays and analyzed patient's RNA when available. Our study revealed a staggering number of splicing mutations in MLH1 exon 10 (77% of the 22 analyzed variants), including mutations directly affecting splice sites and, particularly, mutations altering potential splicing regulatory elements (ESRs). We then used this thoroughly characterized dataset, together with experimental data derived from previous studies on BRCA1, BRCA2, CFTR and NF1, to evaluate the predictive power of 3 in silico approaches recently described as promising tools for pinpointing ESR-mutations. Our results indicate that ΔtESRseq and ΔHZEI-based approaches not only discriminate which variants affect splicing, but also predict the direction and severity of the induced splicing defects. In contrast, the ΔΨ-based approach did not show a compelling predictive power. Our data indicates that exonic splicing mutations are more prevalent than currently appreciated and that they can now be predicted by using bioinformatics methods. These findings have implications for all genetically-caused diseases.

Assessment of the InSiGHT Interpretation Criteria for the Clinical Classification of 24 MLH1 and MSH2 Gene Variants.

Pathogenicity assessment of DNA variants in disease genes to explain their clinical consequences is an integral component of diagnostic molecular testing. The International Society for Gastrointestinal Hereditary Tumors (InSiGHT) has developed specific criteria for the interpretation of mismatch repair (MMR) gene variants. Here, we performed a systematic investigation of 24 MLH1 and MSH2 variants. The assessments were done by analyzing population frequency, segregation, tumor molecular characteristics, RNA effects, protein expression levels, and in vitro MMR activity. Classifications were confirmed for 15 variants and changed for three, and for the first time determined for six novel variants. Overall, based on our results, we propose the introduction of some refinements to the InSiGHT classification rules. The proposed changes have the advantage of homogenizing the InSIGHT interpretation criteria with those set out by the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium for the BRCA1/BRCA2 genes. We also observed that the addition of only few clinical data was sufficient to obtain a more stable classification for variants considered as "likely pathogenic" or "likely nonpathogenic." This shows the importance of obtaining as many as possible points of evidence for variant interpretation, especially from the clinical setting.

Mutational analysis of JAK2, CBL, RUNX1, and NPM1 genes in familial aggregation of hematological malignancies.

Familial aggregation of hematological malignancies has been reported highlighting inherited genetic predisposition. In this study, we targeted four candidate genes: JAK2 and RUNX1 genes assuring a prominent function in hematological process and CBL and NPM1 as proto-oncogenes. Their disruption was described in several sporadic hematological malignancies. The aim of this study is to determine whether JAK2, CBL, RUNX1, and NPM1 germline genes mutations are involved in familial hematological malignancies. Using direct sequencing, we analyzed JAK2 (exons 12 and 14); CBL (exons 7, 8 and 9); NPM1 (exon 12) and the entire RUNX1 in 88 independent families belonging to Tunisian and French populations. Twenty-one sporadic acute leukemias were included in this study. We reported a heterozygous intronic c.1641 + 6 T > C JAK2 variant (rs182123615) found in two independent familial cases diagnosed with gastric lymphoma and Hodgkin lymphoma. The in silico analysis suggested a potential impact on splicing, but the functional splicing minigene reporter assay on rs182123615 variant showed no aberrant transcripts. In one sporadic acute myeloblastic leukemia, we reported an insertion 846 in. TGTT in exon 12 of NPM1 gene that may impact the normal reading frame. The rs182123615 JAK2 variant was described in several contexts including myeloproliferative neoplasms and congenital erythrocytosis and was supposed to be pathogenic. Through this current study, we established the assessment of pathogenicity of rs182123615 and we classified it rather as rare polymorphism.

Identification of variants in the 4q35 gene FAT1 in patients with a facioscapulohumeral dystrophy-like phenotype.

Facioscapulohumeralmuscular dystrophy (FSHD) is linked to copy-number reduction (N < 10) of the 4q D4Z4 subtelomeric array, in association with DUX4-permissive haplotypes. This main form is indicated as FSHD1. FSHD-like phenotypes may also appear in the absence of D4Z4 copy-number reduction. Variants of the SMCHD1 gene have been reported to associate with D4Z4 hypomethylation in DUX4-compatible haplotypes, thus defining FSHD2. Recently, mice carrying a muscle-specific knock-out of the protocadherin gene Fat1 or its constitutive hypomorphic allele were shown to develop muscular and nonmuscular defects mimicking human FSHD. Here, we report FAT1 variants in a group of patients presenting with neuromuscular symptoms reminiscent of FSHD. The patients do not carry D4Z4 copy-number reduction, 4q hypomethylation, or SMCHD1 variants. However, abnormal splicing of the FAT1 transcript is predicted for all identified variants. To determine their pathogenicity, we elaborated a minigene approach coupled to an antisense oligonucleotide (AON) assay. In vitro, four out of five selected variants induced partial or complete alteration of splicing by creating new splice sites or modifying splicing regulators. AONs confirmed these effects. Altered transcripts may affect FAT1 protein interactions or stability. Altogether, our data suggest that defective FAT1 is associated with an FSHD-like phenotype.

Functional analysis of a large set of BRCA2 exon 7 variants highlights the predictive value of hexamer scores in detecting alterations of exonic splicing regulatory elements.

Exonic variants can alter pre-mRNA splicing either by changing splice sites or by modifying splicing regulatory elements. Often these effects are difficult to predict and are only detected by performing RNA analyses. Here, we analyzed, in a minigene assay, 26 variants identified in the exon 7 of BRCA2, a cancer predisposition gene. Our results revealed eight new exon skipping mutations in this exon: one directly altering the 5' splice site and seven affecting potential regulatory elements. This brings the number of splicing regulatory mutations detected in BRCA2 exon 7 to a total of 11, a remarkably high number considering the total number of variants reported in this exon (n = 36), all tested in our minigene assay. We then exploited this large set of splicing data to test the predictive value of splicing regulator hexamers' scores recently established by Ke et al. (). Comparisons of hexamer-based predictions with our experimental data revealed high sensitivity in detecting variants that increased exon skipping, an important feature for prescreening variants before RNA analysis. In conclusion, hexamer scores represent a promising tool for predicting the biological consequences of exonic variants and may have important applications for the interpretation of variants detected by high-throughput sequencing.

Multiple sequence variants of BRCA2 exon 7 alter splicing regulation.

BACKGROUND: Exonic variants of unknown biological significance (VUS) identified in patients can affect mRNA splicing, either by changing 5' or 3' splice sites or by modifying splicing regulatory elements. Bioinformatic predictions of these elements are still inaccurate and only few such elements have been functionally mapped in BRCA2. We studied the effect on splicing of eight exon 7 VUS, selected from the French UMD-BRCA2 mutation database. METHODS: We performed splicing minigene assays and analyses of patient RNA. We also developed a pyrosequencing-based quantitative assay, to measure, in patient RNA, the relative contribution of each allele to the production of exon 7-containing transcripts. Moreover, an exonic splicing enhancer (ESE)-dependent minigene assay was used to evaluate the splicing regulatory properties of wild-type and mutant segments. RESULTS: Six out of the eight variants induced splicing defects. In the minigene assay, c.517G>T and c.631G>A altered the natural splice sites, c.572A>G created a new 5' splice site, and c.520C>T, c.587G>A and c.617C>G induced exon 7 skipping (66%, 25% and 46%, respectively). Pyrosequencing of patient RNA confirmed these levels of exon skipping for c.520C>T and c.617C>G. Results from the ESE-dependent minigene assay indicated that c.520C>T and c.587G>A disturb splicing regulatory elements. CONCLUSIONS: BRCA2 exon 7 splicing is regulated by multiple exonic elements and is sensitive to disease-associated sequence variations. Measurements of allelic imbalance in patient-derived RNA and/or quantitative analyses using minigene assays provide valuable estimates of the extent of partial splicing defects. Assessment of pathogenicity of variants with partial splicing effect awaits additional evidence and especially the completion of segregation analyses.

Guidelines for splicing analysis in molecular diagnosis derived from a set of 327 combined in silico/in vitro studies on BRCA1 and BRCA2 variants.

Assessing the impact of variants of unknown significance (VUS) on splicing is a key issue in molecular diagnosis. This impact can be predicted by in silico tools, but proper evaluation and user guidelines are lacking. To fill this gap, we embarked upon the largest BRCA1 and BRCA2 splice study to date by testing 272 VUSs (327 analyses) within the BRCA splice network of Unicancer. All these VUSs were analyzed by using six tools (splice site prediction by neural network, splice site finder (SSF), MaxEntScan (MES), ESE finder, relative enhancer and silencer classification by unanimous enrichment, and human splicing finder) and the predictions obtained were compared with transcript analysis results. Combining MES and SSF gave 96% sensitivity and 83% specificity for VUSs occurring in the vicinity of consensus splice sites, that is, the surrounding 11 and 14 bases for the 5' and 3' sites, respectively. This study was also an opportunity to define guidelines for transcript analysis along with a tentative classification of splice variants. The guidelines drawn from this large series should be useful for the whole community, particularly in the context of growing sequencing capacities that require robust pipelines for variant interpretation.

Splicing analyses for variants in MMR genes: best practice recommendations from the European Mismatch Repair Working Group.

Over 20% of the DNA mismatch repair (MMR) germline variants in suspected Lynch syndrome patients are classified as variants of uncertain significance (VUS). Well-established functional assays are pivotal for assessing the biological impact of these variants and provide relevant evidence for clinical classification. In our collaborative European Mismatch Repair Working Group (EMMR-WG) we compared three different experimental approaches for evaluating the effect of seven variants on mRNA splicing in MMR genes: (i) RT-PCR of full-length transcripts (FLT), (ii) RT-PCR of targeted transcript sections (TTS), both from patient biological samples and (iii) minigene splicing assays. An overall good concordance was observed between splicing patterns in TTS, FLT and minigene analyses for all variants. The FLT analysis depicted a higher number of different isoforms and mitigated PCR-bias towards shorter isoforms. TTS analyses may miss aberrant isoforms and minigene assays may under/overestimate the severity of certain splicing defects. The interpretation of the experimental findings must be cautious to adequately discriminate abnormal events from physiological complex alternative splicing patterns. A consensus strategy for investigating the impact of MMR variants on splicing was defined. First, RNA should be obtained from patient's cell cultures (such as fresh lymphocyte cultures) incubated with/without a nonsense-mediated decay inhibitor. Second, FLT RT-PCR analysis is recommended to oversee all generated isoforms. Third, TTS analysis and minigene assays are useful independent approaches for verifying and clarifying FLT results. The use of several methodologies is likely to increase the strength of the experimental evidence which contributes to improve variant interpretation.

The BRCA1 c.5434C->G (p.Pro1812Ala) variant induces a deleterious exon 23 skipping by affecting exonic splicing regulatory elements.

BACKGROUND A large fraction of the sequence variants of unknown significance or unclassified variants (UVs) could be pathogenic by affecting mRNA splicing. The breast and ovarian cancer susceptibility gene BRCA1 exhibits a large spectrum of sequence variation but only two variants, both located in exon 18, have been shown experimentally to affect splicing regulatory elements. The present study investigated the impact on splicing of the variant BRCA1 c.5434C-->G (p.Pro1812Ala), identified in an ovarian cancer patient. This variant has previously been studied at the protein level with inconclusive results concerning its pathogenic role. METHODS Analysis of RNA from patient peripheral blood was performed by RT-PCR. The effect of the variant was tested by using splicing reporter hybrid minigene assays. RESULTS Using patient RNA analyses and hybrid minigene assays, we showed that this variant induces a major splicing defect, with skipping of exon 23, resulting in frameshift and predicted protein termination within the second BRCT domain. Moreover, we showed that the segment c.5420-5449 of BRCA1, in the centre of exon 23, exhibits splicing enhancer properties. This enhancement is abolished by the c.5434C-->G mutation, indicating that the nucleotide change, in this highly conserved region, affects a splicing regulatory element. Bioinformatics analyses predict that the mutation c.5434C-->G creates an hnRNPA1 dependent splicing silencer. CONCLUSION These data, together with segregation data, argue for the classification of BRCA1 c.5434C-->G as a pathogenic splicing mutation. These results also suggest that UVs in highly conserved nucleotide sequences of short exons may be good candidates for detecting functionally relevant splicing regulatory elements.

Germline APC mutation spectrum derived from 863 genomic variations identified through a 15-year medical genetics service to French patients with FAP.

Heterozygous APC germline alteration is responsible for familial adenomatous polyposis, a colon cancer predisposition with almost complete penetrance. Point mutations generally lead to truncated proteins or no protein at all. They mainly involve exon 3 to codon 1700 (exon 15). The work presented here delineates precisely the APC mutation spectrum from 15 years of systematic molecular screening which identified 863 independent alterations in the French population.