Comparative in silico analyses and experimental validation of novel splice site and missense mutations in the genes MLH1 and MSH2

Understanding the Genetic Landscape of Pancreatic Ductal Adenocarcinoma to Support Personalized Medicine: A Systematic Review

Pancreatic ductal adenocarcinoma (PDAC) is one of the most fatal malignancies worldwide. While population-wide screening recommendations for PDAC in asymptomatic individuals are not achievable due to its relatively low incidence, pancreatic cancer surveillance programs are recommended for patients with germline causative variants in PDAC susceptibility genes or a strong family history. In this study, we sought to determine the prevalence and significance of germline alterations in major genes (ATM, BRCA1, BRCA2, CDKN2A, EPCAM, MLH1, MSH2, MSH6, PALB2, PMS2, STK11, TP53) involved in PDAC susceptibility. We performed a systematic review of PubMed publications reporting germline variants identified in these genes in PDAC patients. Overall, the retrieved articles included 1493 PDAC patients. A high proportion of these patients (n = 1225/1493, 82%) were found to harbor alterations in genes (ATM, BRCA1, BRCA2, PALB2) involved in the homologous recombination repair (HRR) pathway. Specifically, the remaining PDAC patients were reported to carry alterations in genes playing a role in other cancer pathways (CDKN2A, STK11, TP53; n = 181/1493, 12.1%) or in the mismatch repair (MMR) pathway (MLH1, MSH2, MSH6, PMS2; n = 87/1493, 5.8%). Our findings highlight the importance of germline genetic characterization in PDAC patients for better personalized targeted therapies, clinical management, and surveillance.

Synonymous Variants: Necessary Nuance in our Understanding of Cancer Drivers and Treatment Outcomes

Once called "silent mutations" and assumed to have no effect on protein structure and function, synonymous variants are now recognized to be drivers for some cancers. There have been significant advances in our understanding of the numerous mechanisms by which synonymous single nucleotide variants (sSNVs) can affect protein structure and function by affecting pre-mRNA splicing, mRNA expression, stability, folding, miRNA binding, translation kinetics, and co-translational folding. This review highlights the need for considering sSNVs in cancer biology to gain a better understanding of the genetic determinants of human cancers and to improve their diagnosis and treatment. We surveyed the literature for reports of sSNVs in cancer and found numerous studies on the consequences of sSNVs on gene function with supporting in vitro evidence. We also found reports of sSNVs that have statistically significant associations with specific cancer types but for which in vitro studies are lacking to support the reported associations. Additionally, we found reports of germline and somatic sSNVs that were observed in numerous clinical studies and for which in silico analysis predicts possible effects on gene function. We provide a review of these investigations and discuss necessary future studies to elucidate the mechanisms by which sSNVs disrupt protein function and are play a role in tumorigeneses, cancer progression, and treatment efficacy. As splicing dysregulation is one of the most well recognized mechanisms by which sSNVs impact protein function, we also include our own in silico analysis for predicting which sSNVs may disrupt pre-mRNA splicing.

Comprehensive Constitutional Genetic and Epigenetic Characterization of Lynch-Like Individuals

The causal mechanism for cancer predisposition in Lynch-like syndrome (LLS) remains unknown. Our aim was to elucidate the constitutional basis of mismatch repair (MMR) deficiency in LLS patients throughout a comprehensive (epi)genetic analysis. One hundred and fifteen LLS patients harboring MMR-deficient tumors and no germline MMR mutations were included. Mutational analysis of 26 colorectal cancer (CRC)-associated genes was performed. Pathogenicity of MMR variants was assessed by splicing and multifactorial likelihood analyses. Genome-wide methylome analysis was performed by the Infinium Human Methylation 450K Bead Chip. The multigene panel analysis revealed the presence of two MMR gene truncating mutations not previously found. Of a total of 15 additional MMR variants identified, five -present in 6 unrelated individuals- were reclassified as pathogenic. In addition, 13 predicted deleterious variants in other CRC-predisposing genes were found in 12 probands. Methylome analysis detected one constitutional MLH1 epimutation, but no additional differentially methylated regions were identified in LLS compared to LS patients or cancer-free individuals. In conclusion, the use of an ad-hoc designed gene panel combined with pathogenicity assessment of variants allowed the identification of deleterious MMR mutations as well as new LLS candidate causal genes. Constitutional epimutations in non-LS-associated genes are not responsible for LLS.

TAMM41 is required for heart valve differentiation via regulation of PINK-PARK2 dependent mitophagy

TAMM41, located within the congenital heart diseases (CHD) sensitive region of 3p25 deletion syndrome, is a mitochondrial membrane maintenance protein critical for yeast survival, but its function in higher vertebrates remains unknown. Via in vivo zebrafish model, we found that tamm41 is highly expressed in the developing heart and deficiency of which led to heart valve abnormalities. Molecular mechanistic studies revealed that TAMM41 interacts and modulates the PINK1-PARK2 dependent mitophagy pathway, thereby implicating TAMM41 in heart valve development during zebrafish embryonic cardiogenesis. Furthermore, through screening of the congenital heart diseases (CHD) sensitive region of 3p25 deletion syndrome among 118 sporadic atrioventricular septal defect (AVSD) patients, we identified three cases carrying heterozygous pathogenic intronic variants of TAMM41. All three cases lacked normal full-length TAMM41 transcripts, most likely due to specific expression of the mutant allele. Collectively, our studies highlight essential roles for TAMM41-dependent mitophagy in development of the heart and provide novel insights into the etiology of AVSD.

In Silico Systems Biology Analysis of Variants of Uncertain Significance in Lynch Syndrome Supports the Prioritization of Functional Molecular Validation

Lynch syndrome (LS) is a genetic condition secondary to germline alterations in the DNA mismatch repair (MMR) genes with 30% of changes being variants of uncertain significance (VUS). Our aim was to perform an in silico reclassification of VUS from a large single institutional cohort that will help prioritizing functional validation. A total of 54 VUS were detected with 33 (61%) novel variants. We integrated family history, pathology, and genetic information along with supporting evidence from eight different in silico tools at the RNA and protein level. Our assessment allowed us to reclassify 54% (29/54) of the VUS as probably damaging, 13% (7/54) as possibly damaging, and 28% (15/54) as probably neutral. There are more than 1,000 VUS reported in MMR genes and our approach facilitates the prioritization of further functional efforts to assess the pathogenicity to those classified as probably damaging. Cancer Prev Res; 10(10); 580-7. ©2017 AACR.

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.

Exon-skipping and mRNA decay in human liver tissue: molecular consequences of pathogenic bile salt export pump mutations

The bile salt export pump BSEP mediates bile formation. Over 150 BSEP mutations are associated with progressive familial intrahepatic cholestasis type 2 (PFIC-2), with few characterised specifically. We examined liver tissues from two PFIC-2 patients compound heterozygous for the splice-site mutation c.150 + 3A > C and either c.2783_2787dup5 resulting in a frameshift with a premature termination codon (child 1) or p.R832C (child 2). Splicing was analysed with a minigene system and mRNA sequencing from patients’ livers. Protein expression was shown by immunofluorescence. Using the minigene, c.150 + 3A > C causes complete skipping of exon 3. In liver tissue of child 1, c.2783_2787dup5 was found on DNA but not on mRNA level, implying nonsense-mediated mRNA decay (NMD) when c.2783_2787dup5 is present. Still, BSEP protein as well as mRNA with and without exon 3 were detectable and can be assigned to the c.150 + 3A > C allele. Correctly spliced transcripts despite c.150 + 3A > C were also confirmed in liver of child 2. In conclusion, we provide evidence (1) for effective NMD due to a BSEP frameshift mutation and (2) partial exon-skipping due to c.150 + 3A > C. The results illustrate that the extent of exon-skipping depends on the genomic and cellular context and that regulation of splicing may have therapeutic potential.

RNA splicing. The human splicing code reveals new insights into the genetic determinants of disease

To facilitate precision medicine and whole-genome annotation, we developed a machine-learning technique that scores how strongly genetic variants affect RNA splicing, whose alteration contributes to many diseases. Analysis of more than 650,000 intronic and exonic variants revealed widespread patterns of mutation-driven aberrant splicing. Intronic disease mutations that are more than 30 nucleotides from any splice site alter splicing nine times as often as common variants, and missense exonic disease mutations that have the least impact on protein function are five times as likely as others to alter splicing. We detected tens of thousands of disease-causing mutations, including those involved in cancers and spinal muscular atrophy. Examination of intronic and exonic variants found using whole-genome sequencing of individuals with autism revealed misspliced genes with neurodevelopmental phenotypes. Our approach provides evidence for causal variants and should enable new discoveries in precision medicine.

Alternative splicing of DNA damage response genes and gastrointestinal cancers

Alternative splicing, which is a common phenomenon in mammalian genomes, is a fundamental process of gene regulation and contributes to great protein diversity. Alternative splicing events not only occur in the normal gene regulation process but are also closely related to certain diseases including cancer. In this review, we briefly demonstrate the concept of alternative splicing and DNA damage and describe the association of alternative splicing and cancer pathogenesis, focusing on the potential relationship of alternative splicing, DNA damage, and gastrointestinal cancers. We will also discuss whether alternative splicing leads to genetic instability, which is considered to be a driving force for tumorigenesis. Better understanding of the role and mechanism of alternative splicing in tumorigenesis may provide new directions for future cancer studies.

Functional testing strategy for coding genetic variants of unclear significance in MLH1 in Lynch syndrome diagnosis

Lynch syndrome is caused by inactivating mutations in the MLH1 gene, but genetic variants of unclear significance frequently preclude diagnosis. Functional testing can reveal variant-conferred defects in gene or protein function. Based on functional defect frequencies and clinical applicability of test systems, we developed a functional testing strategy aimed at efficiently detecting pathogenic defects in coding MLH1 variants. In this strategy, tests of repair activity and expression are prioritized over analyses of subcellular protein localization and messenger RNA (mRNA) formation. This strategy was used for four unclear coding MLH1 variants (p.Asp41His, p.Leu507Phe, p.Gln689Arg, p.Glu605del + p.Val716Met). Expression was analyzed using a transfection system, mismatch repair (MMR) activity by complementation in vitro, mRNA formation by reverse transcriptase-PCR in carrier lymphocyte mRNA, and subcellular localization with dye-labeled fusion constructs. All tests included clinically meaningful controls. The strategy enabled efficient identification of defects in two unclear variants: the p.Asp41His variant showed loss of MMR activity, whereas the compound variant p.Glu605del + p.Val716Met had a defect of expression. This expression defect was significantly stronger than the pathogenic expression reference variant analyzed in parallel, therefore the defect of the compound variant is also pathogenic. Interestingly, the expression defect was caused additively by both of the compound variants, at least one of which is non-pathogenic when occurring by itself. Tests were neutral for p.Leu507Phe and p.Gln689Arg, and the results were consistent with available clinical data. We finally discuss the improved sensitivity and efficiency of the applied strategy and its limitations in analyzing unclear coding MLH1 variants.

Genomic HEXploring allows landscaping of novel potential splicing regulatory elements

Effective splice site selection is critically controlled by flanking splicing regulatory elements (SREs) that can enhance or repress splice site use. Although several computational algorithms currently identify a multitude of potential SRE motifs, their predictive power with respect to mutation effects is limited. Following a RESCUE-type approach, we defined a hexamer-based ‘HEXplorer score’ as average Z-score of all six hexamers overlapping with a given nucleotide in an arbitrary genomic sequence. Plotted along genomic regions, HEXplorer score profiles varied slowly in the vicinity of splice sites. They reflected the respective splice enhancing and silencing properties of splice site neighborhoods beyond the identification of single dedicated SRE motifs. In particular, HEXplorer score differences between mutant and reference sequences faithfully represented exonic mutation effects on splice site usage. Using the HIV-1 pre-mRNA as a model system highly dependent on SREs, we found an excellent correlation in 29 mutations between splicing activity and HEXplorer score. We successfully predicted and confirmed five novel SREs and optimized mutations inactivating a known silencer. The HEXplorer score allowed landscaping of splicing regulatory regions, provided a quantitative measure of mutation effects on splice enhancing and silencing properties and permitted calculation of the mutationally most effective nucleotide.

Functional examination of MLH1, MSH2, and MSH6 intronic mutations identified in Danish colorectal cancer patients

Background

Germ-line mutations in the DNA mismatch repair genes MLH1, MSH2, and MSH6 predispose to the development of colorectal cancer (Lynch syndrome or hereditary nonpolyposis colorectal cancer). These mutations include disease-causing frame-shift, nonsense, and splicing mutations as well as large genomic rearrangements. However, a large number of mutations, including missense, silent, and intronic variants, are classified as variants of unknown clinical significance.

Methods

Intronic MLH1, MSH2, or MSH6 variants were investigated using in silico prediction tools and mini-gene assay to asses the effect on splicing.

Results

We describe in silico and in vitro characterization of nine intronic MLH1, MSH2, or MSH6 mutations identified in Danish colorectal cancer patients, of which four mutations are novel. The analysis revealed aberrant splicing of five mutations (MLH1 c.588 + 5G > A, MLH1 c.677 + 3A > T, MLH1 c.1732-2A > T, MSH2 c.1276 + 1G > T, and MSH2 c.1662-2A > C), while four mutations had no effect on splicing compared to wild type (MLH1 c.117-34A > T, MLH1 c.1039-8 T > A, MSH2 c.2459-18delT, and MSH6 c.3439-16C > T).

Conclusions

In conclusion, we classify five MLH1/MSH2 mutations as pathogenic, whereas four MLH1/MSH2/MSH6 mutations are classified as neutral. This study supports the notion that in silico prediction tools and mini-gene assays are important for the classification of intronic variants, and thereby crucial for the genetic counseling of patients and their family members.

RNA splicing: a new player in the DNA damage response

It is widely accepted that tumorigenesis is a multistep process characterized by the sequential accumulation of genetic alterations. However, the molecular basis of genomic instability in cancer is still partially understood. The observation that hereditary cancers are often characterized by mutations in DNA repair and checkpoint genes suggests that accumulation of DNA damage is a major contributor to the oncogenic transformation. It is therefore of great interest to identify all the cellular pathways that contribute to the response to DNA damage. Recently, RNA processing has emerged as a novel pathway that may contribute to the maintenance of genome stability. In this review, we illustrate several different mechanisms through which pre-mRNA splicing and genomic stability can influence each other. We specifically focus on the role of splicing factors in the DNA damage response and describe how, in turn, activation of the DDR can influence the activity of splicing factors.

Genetic screens to identify pathogenic gene variants in the common cancer predisposition Lynch syndrome

In many individuals suspected of the common cancer predisposition Lynch syndrome, variants of unclear significance (VUS), rather than an obviously pathogenic mutations, are identified in one of the DNA mismatch repair (MMR) genes. The uncertainty of whether such VUS inactivate MMR, and therefore are pathogenic, precludes targeted healthcare for both carriers and their relatives. To facilitate the identification of pathogenic VUS, we have developed an in cellulo genetic screen-based procedure for the large-scale mutagenization, identification, and cataloging of residues of MMR genes critical for MMR gene function. When a residue identified as mutated in an individual suspected of Lynch syndrome is listed as critical in such a reverse diagnosis catalog, there is a high probability that the corresponding human VUS is pathogenic. To investigate the applicability of this approach, we have generated and validated a prototypic reverse diagnosis catalog for the MMR gene MutS Homolog 2 (Msh2) by mutagenizing, identifying, and cataloging 26 deleterious mutations in 23 amino acids. Extensive in vivo and in vitro analysis of mutants listed in the catalog revealed both recessive and dominant-negative phenotypes. Nearly half of these critical residues match with VUS previously identified in individuals suspected of Lynch syndrome. This aids in the assignment of pathogenicity to these human VUS and validates the approach described here as a diagnostic tool. In a wider perspective, this work provides a model for the translation of personalized genomics into targeted healthcare.

Comprehensive functional assessment of MLH1 variants of unknown significance

Lynch syndrome is associated with germline mutations in DNA mismatch repair (MMR) genes. Up to 30% of DNA changes found are variants of unknown significance (VUS). Our aim was to assess the pathogenicity of eight MLH1 VUS identified in patients suspected of Lynch syndrome. All of them are novel or not previously characterized. For their classification, we followed a strategy that integrates family history, tumor pathology, and control frequency data with a variety of in silico and in vitro analyses at RNA and protein level, such as MMR assay, MLH1 and PMS2 expression, and subcellular localization. Five MLH1 VUS were classified as pathogenic: c.[248G>T(;)306G>C], c.[780C>G;788A>C], and c.791-7T>A affected mRNA processing, whereas c.218T>C (p.L73P) and c.244A>G [corrected] (p.T82A) impaired MMR activity. Two other VUS were considered likely neutral: the silent c.702G>A variant did not affect mRNA processing or stability, and c.974G>A (p.R325Q) did not influence MMR function. In contrast, variant c.25C>T (p.R9W) could not be classified, as it associated with intermediate levels of MMR activity. Comprehensive functional assessment of MLH1 variants was useful in their classification and became relevant in the diagnosis and genetic counseling of carrier families.

Classification of mismatch repair gene missense variants with PON-MMR

Numerous mismatch repair (MMR) gene variants have been identified in Lynch syndrome and other cancer patients, but knowledge about their pathogenicity is frequently missing. The diagnosis and treatment of patients would benefit from knowing which variants are disease related. Bioinformatic approaches are well suited to the problem and can handle large numbers of cases. Functional effects were revealed based on literature for 168 MMR missense variants. Performance of numerous prediction methods was tested with this dataset. Among the tested tools, only the results of tolerance prediction methods correlated to functional information, however, with poor performance. Therefore, a novel consensus-based predictor was developed. The novel prediction method, pathogenic-or-not mismatch repair (PON-MMR), achieved accuracy of 0.87 and Matthews correlation coefficient of 0.77 on the experimentally verified variants. When applied to 616 MMR cases with unknown effects, 81 missense variants were predicted to be pathogenic and 167 neutral. With PON-MMR, the number of MMR missense variants with unknown effect was reduced by classifying a large number of cases as likely pathogenic or benign. The results can be used, for example, to prioritize cases for experimental studies and assist in the classification of cases.

Complexities of 5'splice site definition: implications in clinical analyses

In higher eukaryotes, the 5' splice site (5'ss) is initially recognized through an RNA-RNA interaction by U1 small nuclear ribonucleoprotein (U1 snRNP). This event represents one of the key steps in initial spliceosomal assembly and many disease-associated mutations in humans often disrupt this process. Beside base pair complementarity, 5'ss recognition can also be modified by additional factors such as RNA secondary structures or the specific binding of other nuclear proteins. In this work, we have focused on investigating a few examples of changes detected within the 5'ss in patients, that would not be immediately considered "disease causing mutations". We show that the splicing outcome of very similar mutations can be very different due to variations in trans-acting factor(s) interactions and specific context influences. Using several NF1 donor sites and SELEX approaches as experimental models, we have examined the binding properties of particular sequence motifs such as GGGU found in donor sites, and how the sequence context can change their interaction with hnRNPs such as H/F and A1/A2. Our results clearly show that even minor differences in local nucleotide context can differentially affect the binding ability of these factors to the GGGU core. Finally, using a previously identified mutation in KCNH2 that resulted in intron retention we show how very similar 5'ss mutations found in patients can have a very different splicing outcome due to the neighbouring sequence context, thus highlighting the general need to approach splicing problems with suitable experimental approaches.

A novel splicing mutation alters DSPP transcription and leads to dentinogenesis imperfecta type II

Dentinogenesis imperfecta (DGI) type II is an autosomal dominant disease characterized by a serious disorders in teeth. Mutations of dentin sialophosphoprotein (DSPP) gene were revealed to be the causation of DGI type II (DGI-II). In this study, we identified a novel mutation (NG_011595.1:g.8662T>C, c.135+2T>C) lying in the splice donor site of intron 3 of DSPP gene in a Chinese Han DGI-II pedigree. It was found in all affected subjects but not in unaffected ones or other unrelated healthy controls. The function of the mutant DSPP gene, which was predicted online and subsequently confirmed by in vitro splicing analysis, was the loss of splicing of intron 3, leading to the extended length of DSPP mRNA. For the first time, the functional non-splicing of intron was revealed in a novel DSPP mutation and was considered as the causation of DGI-II. It was also indicated that splicing was of key importance to the function of DSPP and this splice donor site might be a sensitive mutation hot spot. Our findings combined with other reports would facilitate the genetic diagnosis of DGI-II, shed light on its gene therapy and help to finally conquer human diseases.

Cryptic splice sites and split genes

We describe a new program called cryptic splice finder (CSF) that can reliably identify cryptic splice sites (css), so providing a useful tool to help investigate splicing mutations in genetic disease. We report that many css are not entirely dormant and are often already active at low levels in normal genes prior to their enhancement in genetic disease. We also report a fascinating correlation between the positions of css and introns, whereby css within the exons of one species frequently match the exact position of introns in equivalent genes from another species. These results strongly indicate that many introns were inserted into css during evolution and they also imply that the splicing information that lies outside some introns can be independently recognized by the splicing machinery and was in place prior to intron insertion. This indicates that non-intronic splicing information had a key role in shaping the split structure of eukaryote genes.

Correct mRNA processing at a mutant TT splice donor in FANCC ameliorates the clinical phenotype in patients and is enhanced by delivery of suppressor U1 snRNAs

The U1 small nuclear RNA (U1 snRNA) as a component of the major U2-dependent spliceosome recognizes 5' splice sites (5'ss) containing GT as the canonical dinucleotide in the intronic positions +1 and +2. The c.165+1G>T germline mutation in the 5'ss of exon 2 of the Fanconi anemia C (FANCC) gene commonly predicted to prevent correct splicing was identified in nine FA patients from three pedigrees. RT-PCR analysis of the endogenous FANCC mRNA splicing pattern of patient-derived fibroblasts revealed aberrant mRNA processing, but surprisingly also correct splicing at the TT dinucleotide, albeit with lower efficiency. This consequently resulted in low levels of correctly spliced transcript and minute levels of normal posttranslationally processed FANCD2 protein, indicating that this naturally occurring TT splicing might contribute to the milder clinical manifestations of the disease in these patients. Functional analysis of this FANCC 5'ss within splicing reporters revealed that both the noncanonical TT dinucleotide and the genomic context of FANCC were required for the residual correct splicing at this mutant 5'ss. Finally, use of lentiviral vectors as a delivery system to introduce expression cassettes for TT-adapted U1 snRNAs into primary FANCC patient fibroblasts allowed the correction of the DNA-damage-induced G2 cell-cycle arrest in these cells, thus representing an alternative transcript-targeting approach for genetic therapy of inherited splice-site mutations.

Germline mutations in breast and ovarian cancer pedigrees establish RAD51C as a human cancer susceptibility gene

Germline mutations in a number of genes involved in the recombinational repair of DNA double-strand breaks are associated with predisposition to breast and ovarian cancer. RAD51C is essential for homologous recombination repair, and a biallelic missense mutation can cause a Fanconi anemia-like phenotype. In index cases from 1,100 German families with gynecological malignancies, we identified six monoallelic pathogenic mutations in RAD51C that confer an increased risk for breast and ovarian cancer. These include two frameshift-causing insertions, two splice-site mutations and two nonfunctional missense mutations. The mutations were found exclusively within 480 pedigrees with the occurrence of both breast and ovarian tumors (BC/OC; 1.3%) and not in 620 pedigrees with breast cancer only or in 2,912 healthy German controls. These results provide the first unambiguous evidence of highly penetrant mutations associated with human cancer in a RAD51 paralog and support the 'common disease, rare allele' hypothesis.