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/.

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.

Mutations in Citron Kinase Cause Recessive Microlissencephaly with Multinucleated Neurons.

Primary microcephaly is a neurodevelopmental disorder that is caused by a reduction in brain size as a result of defects in the proliferation of neural progenitor cells during development. Mutations in genes encoding proteins that localize to the mitotic spindle and centrosomes have been implicated in the pathogenicity of primary microcephaly. In contrast, the contractile ring and midbody required for cytokinesis, the final stage of mitosis, have not previously been implicated by human genetics in the molecular mechanisms of this phenotype. Citron kinase (CIT) is a multi-domain protein that localizes to the cleavage furrow and midbody of mitotic cells, where it is required for the completion of cytokinesis. Rodent models of Cit deficiency highlighted the role of this gene in neurogenesis and microcephaly over a decade ago. Here, we identify recessively inherited pathogenic variants in CIT as the genetic basis of severe microcephaly and neonatal death. We present postmortem data showing that CIT is critical to building a normally sized human brain. Consistent with cytokinesis defects attributed to CIT, multinucleated neurons were observed throughout the cerebral cortex and cerebellum of an affected proband, expanding our understanding of mechanisms attributed to primary microcephaly.

Identification of genetic variants for clinical management of familial colorectal tumors.

BACKGROUND: The genetic mechanisms for families who meet the clinical criteria for Lynch syndrome (LS) but do not carry pathogenic variants in the mismatch repair (MMR) genes are still undetermined. We aimed to study the potential contribution of genes other than MMR genes to the biological and clinical characteristics of Norwegian families fulfilling Amsterdam (AMS) criteria or revised Bethesda guidelines. METHODS: The Hereditary Cancer Biobank of the Norwegian Radium Hospital was interrogated to identify individuals with a high risk of developing colorectal cancer (CRC) for whom no pathogenic variants in MMR genes had been found in routine diagnostic DNA sequencing. Forty-four cancer susceptibility genes were selected and analyzed by using our in-house designed TruSeq amplicon-based assay for targeted sequencing. RNA splicing- and protein-dedicated in silico analyses were performed for all variants of unknown significance (VUS). Variants predicted as likely to affect splicing were experimentally analyzed by resorting to minigene assays. RESULTS: We identified a patient who met the revised Bethesda guidelines and carried a likely pathogenic variant in CHEK2 (c.470 T > C, p.I157T). In addition, 25 unique VUS were identified in 18 individuals, of which 2 exonic variants (MAP3K1 c.764A > G and NOTCH3 c.5854G >A) were analyzed in the minigene splicing assay and found not to have an effect on RNA splicing. CONCLUSIONS: Among high-risk CRC patients that fulfill the AMS criteria or revised Bethesda guidelines, targeted gene sequencing identified likely pathogenic variant and VUS in other genes than the MMR genes (CHEK2, NOTCH3 and MAP3K1). Our study suggests that the analysis of genes currently excluded from routine molecular diagnostic screens may confer cancer susceptibility.

Genetic variants of prospectively demonstrated phenocopies in BRCA1/2 kindreds.

BACKGROUND: In kindreds carrying path_BRCA1/2 variants, some women in these families will develop cancer despite testing negative for the family's pathogenic variant. These families may have additional genetic variants, which not only may increase the susceptibility of the families' path_BRCA1/2, but also be capable of causing cancer in the absence of the path_BRCA1/2 variants. We aimed to identify novel genetic variants in prospectively detected breast cancer (BC) or gynecological cancer cases tested negative for their families' pathogenic BRCA1/2 variant (path_BRCA1 or path_BRCA2). METHODS: Women with BC or gynecological cancer who had tested negative for path_BRCA1 or path_BRCA2 variants were included. Forty-four cancer susceptibility genes were screened for genetic variation through a targeted amplicon-based sequencing assay. Protein- and RNA splicing-dedicated in silico analyses were performed for all variants of unknown significance (VUS). Variants predicted as the ones most likely affecting pre-mRNA splicing were experimentally analyzed in a minigene assay. RESULTS: We identified 48 women who were tested negative for their family's path_BRCA1 (n = 13) or path_BRCA2 (n = 35) variants. Pathogenic variants in the ATM, BRCA2, MSH6 and MUTYH genes were found in 10% (5/48) of the cases, of whom 15% (2/13) were from path_BRCA1 and 9% (3/35) from path_BRCA2 families. Out of the 26 unique VUS, 3 (12%) were predicted to affect RNA splicing (APC c.721G > A, MAP3K1 c.764A > G and MSH2 c.815C > T). However, by using a minigene, assay we here show that APC c.721G > A does not cause a splicing defect, similarly to what has been recently reported for the MAP3K1 c.764A > G. The MSH2 c.815C > T was previously described as causing partial exon skipping and it was identified in this work together with the path_BRCA2 c.9382C > T (p.R3128X). CONCLUSION: All women in breast or breast/ovarian cancer kindreds would benefit from being offered genetic testing irrespective of which causative genetic variants have been demonstrated in their relatives.

A survey of the clinicopathological and molecular characteristics of patients with suspected Lynch syndrome in Latin America.

BACKGROUND: Genetic counselling and testing for Lynch syndrome (LS) have recently been introduced in several Latin America countries. We aimed to characterize the clinical, molecular and mismatch repair (MMR) variants spectrum of patients with suspected LS in Latin America. METHODS: Eleven LS hereditary cancer registries and 34 published LS databases were used to identify unrelated families that fulfilled the Amsterdam II (AMSII) criteria and/or the Bethesda guidelines or suggestive of a dominant colorectal (CRC) inheritance syndrome. RESULTS: We performed a thorough investigation of 15 countries and identified 6 countries where germline genetic testing for LS is available and 3 countries where tumor testing is used in the LS diagnosis. The spectrum of pathogenic MMR variants included MLH1 up to 54%, MSH2 up to 43%, MSH6 up to 10%, PMS2 up to 3% and EPCAM up to 0.8%. The Latin America MMR spectrum is broad with a total of 220 different variants which 80% were private and 20% were recurrent. Frequent regions included exons 11 of MLH1 (15%), exon 3 and 7 of MSH2 (17 and 15%, respectively), exon 4 of MSH6 (65%), exons 11 and 13 of PMS2 (31% and 23%, respectively). Sixteen international founder variants in MLH1, MSH2 and MSH6 were identified and 41 (19%) variants have not previously been reported, thus representing novel genetic variants in the MMR genes. The AMSII criteria was the most used clinical criteria to identify pathogenic MMR carriers although microsatellite instability, immunohistochemistry and family history are still the primary methods in several countries where no genetic testing for LS is available yet. CONCLUSION: The Latin America LS pathogenic MMR variants spectrum included new variants, frequently altered genetic regions and potential founder effects, emphasizing the relevance implementing Lynch syndrome genetic testing and counseling in all of Latin America countries.

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.

Skipping Nonsense to Maintain Function: The Paradigm of BRCA2 Exon 12.

Germline nonsense and canonical splice site variants identified in disease-causing genes are generally considered as loss-of-function (LoF) alleles and classified as pathogenic. However, a fraction of such variants could maintain function through their impact on RNA splicing. To test this hypothesis, we used the alternatively spliced BRCA2 exon 12 (E12) as a model system because its in-frame skipping leads to a potentially functional protein. All E12 variants corresponding to putative LoF variants or predicted to alter splicing (n = 40) were selected from human variation databases and characterized for their impact on splicing in minigene assays and, when available, in patient lymphoblastoid cell lines. Moreover, a selection of variants was analyzed in a mouse embryonic stem cell-based functional assay. Using these complementary approaches, we demonstrate that a subset of variants, including nonsense variants, induced in-frame E12 skipping through the modification of splice sites or regulatory elements and, consequently, led to an internally deleted but partially functional protein. These data provide evidence, for the first time in a cancer-predisposition gene, that certain presumed null variants can retain function due to their impact on splicing. Further studies are required to estimate cancer risk associated with these hypomorphic variants. More generally, our findings highlight the need to exercise caution in the interpretation of putative LoF variants susceptible to induce in-frame splicing modifications. SIGNIFICANCE: This study presents evidence that certain presumed loss-of-function variants in a cancer predisposition gene can retain function due to their direct impact on RNA splicing.

Calibration of Pathogenicity Due to Variant-Induced Leaky Splicing Defects by Using BRCA2 Exon 3 as a Model System.

BRCA2 is a clinically actionable gene implicated in breast and ovarian cancer predisposition that has become a high priority target for improving the classification of variants of unknown significance (VUS). Among all BRCA2 VUS, those causing partial/leaky splicing defects are the most challenging to classify because the minimal level of full-length (FL) transcripts required for normal function remains to be established. Here, we explored BRCA2 exon 3 (BRCA2e3) as a model for calibrating variant-induced spliceogenicity and estimating thresholds for BRCA2 haploinsufficiency. In silico predictions, minigene splicing assays, patients' RNA analyses, a mouse embryonic stem cell (mESC) complementation assay and retrieval of patient-related information were combined to determine the minimal requirement of FL BRCA2 transcripts. Of 100 BRCA2e3 variants tested in the minigene assay, 64 were found to be spliceogenic, causing mild to severe RNA defects. Splicing defects were also confirmed in patients' RNA when available. Analysis of a neutral leaky variant (c.231T>G) showed that a reduction of approximately 60% of FL BRCA2 transcripts from a mutant allele does not cause any increase in cancer risk. Moreover, data obtained from mESCs suggest that variants causing a decline in FL BRCA2 with approximately 30% of wild-type are not pathogenic, given that mESCs are fully viable and resistant to DNA-damaging agents in those conditions. In contrast, mESCs producing lower relative amounts of FL BRCA2 exhibited either null or hypomorphic phenotypes. Overall, our findings are likely to have broader implications on the interpretation of BRCA2 variants affecting the splicing pattern of other essential exons. SIGNIFICANCE: These findings demonstrate that BRCA2 tumor suppressor function tolerates substantial reduction in full-length transcripts, helping to determine the pathogenicity of BRCA2 leaky splicing variants, some of which may not increase cancer risk.

Results of multigene panel testing in familial cancer cases without genetic cause demonstrated by single gene testing.

We have surveyed 191 prospectively sampled familial cancer patients with no previously detected pathogenic variant in the BRCA1/2, PTEN, TP53 or DNA mismatch repair genes. In all, 138 breast cancer (BC) cases, 34 colorectal cancer (CRC) and 19 multiple early-onset cancers were included. A panel of 44 cancer-predisposing genes identified 5% (9/191) pathogenic or likely pathogenic variants and 87 variants of uncertain significance (VUS). Pathogenic or likely pathogenic variants were identified mostly in familial BC individuals (7/9) and were located in 5 genes: ATM (3), BRCA2 (1), CHEK2 (1), MSH6 (1) and MUTYH (1), followed by multiple early-onset (2/9) individuals, affecting the CHEK2 and ATM genes. Eleven of the 87 VUS were tested, and 4/11 were found to have an impact on splicing by using a minigene splicing assay. We here report for the first time the splicing anomalies using this assay for the variants ATM c.3806A > G and BUB1 c.677C > T, whereas CHEK1 c.61G > A did not result in any detectable splicing anomaly. Our study confirms the presence of pathogenic or likely pathogenic variants in genes that are not routinely tested in the context of the above-mentioned clinical phenotypes. Interestingly, more than half of the pathogenic germline variants were found in the moderately penetrant ATM and CHEK2 genes, where only truncating variants from these genes are recommended to be reported in clinical genetic testing practice.

Potentially pathogenic germline CHEK2 c.319+2T>A among multiple early-onset cancer families.

To study the potential contribution of genes other than BRCA1/2, PTEN, and TP53 to the biological and clinical characteristics of multiple early-onset cancers in Norwegian families, including early-onset breast cancer, Cowden-like and Li-Fraumeni-like syndromes (BC, CSL and LFL, respectively). The Hereditary Cancer Biobank from the Norwegian Radium Hospital was used to identify early-onset BC, CSL or LFL for whom no pathogenic variants in BRCA1/2, PTEN, or TP53 had been found in routine diagnostic DNA sequencing. Forty-four cancer susceptibility genes were selected and analyzed by our in-house designed TruSeq amplicon-based assay for targeted sequencing. Protein- and RNA splicing-dedicated in silico analyses were performed for all variants of unknown significance (VUS). Variants predicted as the more likely to affect splicing were experimentally analyzed by minigene assay. We identified a CSL individual carrying a variant in CHEK2 (c.319+2T>A, IVS2), here considered as likely pathogenic. Out of the five VUS (BRCA2, CDH1, CHEK2, MAP3K1, NOTCH3) tested in the minigene splicing assay, only NOTCH3 c.14090C>T (p.Ser497Leu) showed a significant effect on RNA splicing, notably by inducing partial skipping of exon 9. Among 13 early-onset BC, CSL and LFL patients, gene panel sequencing identified a potentially pathogenic variant in CHEK2 that affects a canonical RNA splicing signal. Our study provides new information on genetic loci that may affect the risk of developing cancer in these patients and their families, demonstrating that genes presently not routinely tested in molecular diagnostic settings may be important for capturing cancer predisposition in these families.

Biallelic Loss of Function of SORL1 in an Early Onset Alzheimer's Disease Patient.

Heterozygous SORL1 protein truncating variants (PTV) are a strong risk factor for early-onset Alzheimer's disease (EOAD). In case control studies performed at the genome-wide level, PTV definition is usually straightforward. Regarding splice site variants, only those affecting canonical sites are typically included. Some other variants, not annotated as PTV, could, however, affect splicing and hence result in a loss of SORL1 function. We took advantage of the whole exome sequencing data from the 9/484 patients with a previously reported SORL1 PTV in the French EOAD series and searched for a second variant which may affect splicing and eventually result in more than 50% loss of function overall. We found that one patient, known to carry a variant predicted to disrupt the canonical 5' splice site of exon 8, also carried a second novel intronic variant predicted to affect SORL1 splicing of exon 29. Segregation analysis showed that the second variant was located in trans from the known PTV. We performed ex vivo minigene splicing assays and showed that both variants led to the generation of transcripts containing a premature stop codon. This is therefore the first evidence of a human carrying biallelic SORL1 PTV. This patient had a family history of dementia in both maternal and paternal lineages with later ages of onset than the proband himself. However, his 55 years age at onset was in the same ranges as previously published SORL1 heterozygous PTV carriers. This suggests that biallelic loss of SORL1 function is an extremely rare event that was not associated with a dramatically earlier age at onset than heterozygous SORL1 loss-of-function variant carriers, in this single patient.

Full in-frame exon 3 skipping of BRCA2 confers high risk of breast and/or ovarian cancer.

Germline pathogenic variants in the BRCA2 gene are associated with a cumulative high risk of breast/ovarian cancer. Several BRCA2 variants result in complete loss of the exon-3 at the transcript level. The pathogenicity of these variants and the functional impact of loss of exon 3 have yet to be established. As a collaboration of the COVAR clinical trial group (France), and the ENIGMA consortium for investigating breast cancer gene variants, this study evaluated 8 BRCA2 variants resulting in complete deletion of exon 3. Clinical information for 39 families was gathered from Portugal, France, Denmark and Sweden. Multifactorial likelihood analyses were conducted using information from 293 patients, for 7 out of the 8 variants (including 6 intronic). For all variants combined the likelihood ratio in favor of causality was 4.39*1025. These results provide convincing evidence for the pathogenicity of all examined variants that lead to a total exon 3 skipping, and suggest that other variants that result in complete loss of exon 3 at the molecular level could be associated with a high risk of cancer comparable to that associated with classical pathogenic variants in BRCA1 or BRCA2 gene. In addition, our functional study shows, for the first time, that deletion of exon 3 impairs the ability of cells to survive upon Mitomycin-C treatment, supporting lack of function for the altered BRCA2 protein in these cells. Finally, this study demonstrates that any variant leading to expression of only BRCA2 delta-exon 3 will be associated with an increased risk of breast and ovarian cancer.

Clinical and molecular characterization of cystinuria in a French cohort: relevance of assessing large-scale rearrangements and splicing variants.

BACKGROUND: Cystinuria is an autosomal recessive disorder of dibasic amino acid transport in the kidney and the intestine leading to increased urinary cystine excretion and nephrolithiasis. Two genes, SLC3A1 and SLC7A9, coding respectively for rBAT and b0,+AT, account for the genetic basis of cystinuria. METHODS: This study reports the clinical and molecular characterization of a French cohort including 112 cystinuria patients and 25 relatives from 99 families. Molecular screening was performed using sequencing and Quantitative Multiplex PCR of Short Fluorescent Fragments analyses. Functional minigene-based assays have been used to characterize splicing variants. RESULTS: Eighty-eight pathogenic nucleotide changes were identified in SLC3A1 (63) and SLC7A9 (25) genes, of which 42 were novel. Interestingly, 17% (15/88) and 11% (10/88) of the total number of variants correspond, respectively, to large-scale rearrangements and splicing mutations. Functional minigene-based assays were performed for six variants located outside the most conserved sequences of the splice sites; three variants affect splice sites, while three others modify exonic splicing regulatory elements (ESR), in good agreement with a new in silico prediction based on ΔtESRseq values. CONCLUSION: This report expands the spectrum of SLC3A1 and SLC7A9 variants and supports that digenic inheritance is unlikely. Furthermore, it highlights the relevance of assessing large-scale rearrangements and splicing mutations to fully characterize cystinuria patients at the molecular level.