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

5' Region Large Genomic Rearrangements in the BRCA1 Gene in French Families: Identification of a Tandem Triplication and Nine Distinct Deletions with Five Recurrent Breakpoints.

BACKGROUND: Large genomic rearrangements (LGR) in BRCA1 consisting of deletions/duplications of one or several exons have been found throughout the gene with a large proportion occurring in the 5' region from the promoter to exon 2. The aim of this study was to better characterize those LGR in French high-risk breast/ovarian cancer families. METHODS: DNA from 20 families with one apparent duplication and nine deletions was analyzed with a dedicated comparative genomic hybridization (CGH) array, high-resolution BRCA1 Genomic Morse Codes analysis and Sanger sequencing. RESULTS: The apparent duplication was in fact a tandem triplication of exons 1 and 2 and part of intron 2 of BRCA1, fully characterized here for the first time. We calculated a causality score with the multifactorial model from data obtained from six families, classifying this variant as benign. Among the nine deletions detected in this region, eight have never been identified. The breakpoints fell in six recurrent regions and could confirm some specific conformation of the chromatin. CONCLUSIONS: Taken together, our results firmly establish that the BRCA1 5' region is a frequent site of different LGRs and highlight the importance of the segmental duplication and Alu sequences, particularly the very high homologous region, in the mechanism of a recombination event. This also confirmed that those events are not systematically deleterious.

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.

Mutation analysis of PALB2 gene in French breast cancer families.

Several population-based and family-based studies have demonstrated that germline mutations of the PALB2 gene (Partner and Localizer of BRCA2) are associated with an increased risk of breast cancer. Distinct mutation frequencies and spectrums have been described depending on the population studied. Here we describe the first complete PALB2 coding sequence screening in the French population. We screened the complete coding sequence and intron-exon boundaries of PALB2, using the EMMA technique, to assess the contribution of pathogenic mutations in a set of 835 familial breast cancer cases and 662 unrelated controls from the French national study GENESIS and the Paul Strauss Cancer Centre, all previously tested negative for BRCA1 and BRCA2 pathogenic mutations. Our analysis revealed the presence of four novel deleterious mutations: c.1186insT, c.1857delT and c.2850delC in three cases, c.3418dupT in one control. In addition, we identified two in-frame insertion/deletion, 19 missense substitutions (two of them predicted as pathogenic), 9 synonymous variants, 28 variants located in introns and 2 in UTRs, as well as frequent variants. Truncating PALB2 mutations were found in 0.36% of familial breast cancer cases, a frequency lower than the one detected in comparable studies in other populations (0.73-3.40%). This suggests a small but significant contribution of PALB2 mutations to the breast cancer susceptibility in the French population.

An entire exon 3 germ-line rearrangement in the BRCA2 gene: pathogenic relevance of exon 3 deletion in breast cancer predisposition.

BACKGROUND: Germ-line mutations in the BRCA1 and BRCA2 genes are major contributors to hereditary breast/ovarian cancer. Large rearrangements are less frequent in the BRCA2 gene than in BRCA1. We report, here, the first total deletion of exon 3 in the BRCA2 gene that was detected during screening of 2058 index cases from breast/ovarian cancer families for BRCA2 large rearrangements. Deletion of exon 3, which is in phase, does not alter the reading frame. Low levels of alternative transcripts lacking exon 3 (Δ3 delta3 transcript) have been reported in normal tissues, which raises the question whether deletion of exon 3 is pathogenic. METHODS: Large BRCA2 rearrangements were analysed by QMPSF (Quantitative Multiplex PCR of Short Fluorescent Fragments) or MLPA (Multiplex Ligation-Dependent Probe Amplification). The exon 3 deletion was characterized with a "zoom-in" dedicated CGH array to the BRCA2 gene and sequencing. To determine the effect of exon 3 deletion and assess its pathogenic effect, three methods of transcript quantification were used: fragment analysis of FAM-labelled PCR products, specific allelic expression using an intron 2 polymorphism and competitive quantitative RT-PCR. RESULTS: Large rearrangements of BRCA2 were detected in six index cases out of 2058 tested (3% of all deleterious BRCA2 mutations). This study reports the first large rearrangement of the BRCA2 gene that includes all of exon 3 and leads to an in frame deletion of exon 3 at the transcriptional level. Thirty five variants in exon 3 and junction regions of BRCA2 are also reported, that contribute to the interpretation of the pathogenicity of the deletion. The quantitative approaches showed that there are three classes of delta3 BRCA2 transcripts (low, moderate and exclusive). Exclusive expression of the delta3 transcript by the mutant allele and segregation data provide evidence for a causal effect of the exon 3 deletion. CONCLUSION: This paper highlights that large rearrangements and total deletion of exon 3 in the BRCA2 gene could contribute to hereditary breast and/or ovarian cancer. In addition, our findings suggest that, to interpret the pathogenic effect of any variants of exon 3, both accurate transcript quantification and co-segregation analysis are required.

Functional interaction between poly(ADP-Ribose) polymerase 2 (PARP-2) and TRF2: PARP activity negatively regulates TRF2.

The DNA damage-dependent poly(ADP-ribose) polymerase-2 (PARP-2) is, together with PARP-1, an active player of the base excision repair process, thus defining its key role in genome surveillance and protection. Telomeres are specialized DNA-protein structures that protect chromosome ends from being recognized and processed as DNA strand breaks. In mammals, telomere protection depends on the T(2)AG(3) repeat binding protein TRF2, which has been shown to remodel telomeres into large duplex loops (t-loops). In this work we show that PARP-2 physically binds to TRF2 with high affinity. The association of both proteins requires the N-terminal domain of PARP-2 and the myb domain of TRF2. Both partners colocalize at promyelocytic leukemia bodies in immortalized telomerase-negative cells. In addition, our data show that PARP activity regulates the DNA binding activity of TRF2 via both a covalent heteromodification of the dimerization domain of TRF2 and a noncovalent binding of poly(ADP-ribose) to the myb domain of TRF2. PARP-2(-/-) primary cells show normal telomere length as well as normal telomerase activity compared to wild-type cells but display a spontaneously increased frequency of chromosome and chromatid breaks and of ends lacking detectable T(2)AG(3) repeats. Altogether, these results suggest a functional role of PARP-2 activity in the maintenance of telomere integrity.

Poly(ADP-ribose) polymerase-2 (PARP-2) is required for efficient base excision DNA repair in association with PARP-1 and XRCC1.

The DNA damage dependence of poly(ADP-ribose) polymerase-2 (PARP-2) activity is suggestive of its implication in genome surveillance and protection. Here we show that the PARP-2 gene, mainly expressed in actively dividing tissues follows, but to a smaller extent, that of PARP-1 during mouse development. We found that PARP-2 and PARP-1 homo- and heterodimerize; the interacting interfaces, sites of reciprocal modification, have been mapped. PARP-2 was also found to interact with three other proteins involved in the base excision repair pathway: x-ray cross complementing factor 1 (XRCC1), DNA polymerase beta, and DNA ligase III, already known as partners of PARP-1. XRCC1 negatively regulates PARP-2 activity, as it does for PARP-1, while being a polymer acceptor for both PARP-1 and PARP-2. To gain insight into the physiological role of PARP-2 in response to genotoxic stress, we developed by gene disruption mice deficient in PARP-2. Following treatment by the alkylating agent N-nitroso-N-methylurea (MNU), PARP-2-deficient cells displayed an important delay in DNA strand breaks resealing, similar to that observed in PARP-1 deficient cells, thus confirming that PARP-2 is also an active player in base excision repair despite its low capacity to synthesize ADP-ribose polymers.