Analysis of the Impact of Known SPINK1 Missense Variants on Pre-mRNA Splicing and/or mRNA Stability in a Full-Length Gene Assay

Combining full-length gene assay and SpliceAI to interpret the splicing impact of all possible SPINK1 coding variants

BACKGROUND

Single-nucleotide variants (SNVs) within gene coding sequences can significantly impact pre-mRNA splicing, bearing profound implications for pathogenic mechanisms and precision medicine. In this study, we aim to harness the well-established full-length gene splicing assay (FLGSA) in conjunction with SpliceAI to prospectively interpret the splicing effects of all potential coding SNVs within the four-exon SPINK1 gene, a gene associated with chronic pancreatitis.

RESULTS

Our study began with a retrospective analysis of 27 SPINK1 coding SNVs previously assessed using FLGSA, proceeded with a prospective analysis of 35 new FLGSA-tested SPINK1 coding SNVs, followed by data extrapolation, and ended with further validation. In total, we analyzed 67 SPINK1 coding SNVs, which account for 9.3% of the 720 possible coding SNVs. Among these 67 FLGSA-analyzed SNVs, 12 were found to impact splicing. Through detailed comparison of FLGSA results and SpliceAI predictions, we inferred that the remaining 653 untested coding SNVs in the SPINK1 gene are unlikely to significantly affect splicing. Of the 12 splice-altering events, nine produced both normally spliced and aberrantly spliced transcripts, while the remaining three only generated aberrantly spliced transcripts. These splice-impacting SNVs were found solely in exons 1 and 2, notably at the first and/or last coding nucleotides of these exons. Among the 12 splice-altering events, 11 were missense variants (2.17% of 506 potential missense variants), and one was synonymous (0.61% of 164 potential synonymous variants). Notably, adjusting the SpliceAI cut-off to 0.30 instead of the conventional 0.20 would improve specificity without reducing sensitivity.

CONCLUSIONS

By integrating FLGSA with SpliceAI, we have determined that less than 2% (1.67%) of all possible coding SNVs in SPINK1 significantly influence splicing outcomes. Our findings emphasize the critical importance of conducting splicing analysis within the broader genomic sequence context of the study gene and highlight the inherent uncertainties associated with intermediate SpliceAI scores (0.20 to 0.80). This study contributes to the field by being the first to prospectively interpret all potential coding SNVs in a disease-associated gene with a high degree of accuracy, representing a meaningful attempt at shifting from retrospective to prospective variant analysis in the era of exome and genome sequencing.

Crosstalk between KIF1C and PRKAR1A in left atrial myxoma

Cardiac myxoma (CM) is the most common benign cardiac tumor, and most CMs are left atrial myxomas (LAMs). Six variations of KIF1C, c.899 A > T, c.772 T > G, c.352 A > T, c.2895 C > T, c.3049 G > A, and c.*442_*443dup in left atrial myxoma tissues are identified by whole-exome sequencing (WES) and Sanger sequencing. RNA-seq and function experiments show the reduction of the expression of KIF1C and PRKAR1A caused by rare variations of KIF1C. KIF1C is observed to be located in the nucleus, bind to the promoter region of PRKAR1A, and regulate its transcription. Reduction of KIF1C decreases PRKAR1A expression and activates the PKA, which causes an increase in ERK1/2 phosphorylation and SRC-mediated STAT3 activation, a reduction of CDH1, TP53, CDKN1A, and BAX, and eventually promotes tumor formation both in vitro and in vivo. The results suggest that inhibition of KIF1C promotes the pathogenesis of LAM through positive feedback formed by the crosstalk between KIF1C and PRKAR1A.

Analysis of the impact of DGAT1 p.M435L and p.K232A variants on pre-mRNA splicing in a full-length gene assay

DGAT1 is playing a major role in fat metabolism and triacylglyceride synthesis. Only two DGAT1 loss-of-function variants altering milk production traits in cattle have been reported to date, namely p.M435L and p.K232A. The p.M435L variant is a rare alteration and has been associated with skipping of exon 16 which results in a non-functional truncated protein, and the p.K232A-containing haplotype has been associated with modifications of the splicing rate of several DGAT1 introns. In particular, the direct causality of the p.K232A variant in decreasing the splicing rate of the intron 7 junction was validated using a minigene assay in MAC-T cells. As both these DGAT1 variants were shown to be spliceogenic, we developed a full-length gene assay (FLGA) to re-analyse p.M435L and p.K232A variants in HEK293T and MAC-T cells. Qualitative RT-PCR analysis of cells transfected with the full-length DGAT1 expression construct carrying the p.M435L variant highlighted complete skipping of exon 16. The same analysis performed using the construct carrying the p.K232A variant showed moderate differences compared to the wild-type construct, suggesting a possible effect of this variant on the splicing of intron 7. Finally, quantitative RT-PCR analyses of cells transfected with the p.K232A-carrying construct did not show any significant modification on the splicing rate of introns 1, 2 and 7. In conclusion, the DGAT1 FLGA confirmed the p.M435L impact previously observed in vivo, but invalidated the hypothesis whereby the p.K232A variant strongly decreased the splicing rate of intron 7.

Inhibition of mouse trypsin isoforms by SPINK1 and effect of human pancreatitis-associated mutations

Serine protease inhibitor Kazal type 1 (SPINK1) is a trypsin-selective inhibitor protein secreted by the exocrine pancreas. Loss-of-function SPINK1 mutations predispose to chronic pancreatitis through either reduced expression, secretion, or impaired trypsin inhibition. In this study, we aimed to characterize the inhibitory activity of mouse SPINK1 against cationic (T7) and anionic (T8, T9, T20) mouse trypsin isoforms. Kinetic measurements with a peptide substrate, and digestion experiments with β-casein indicated that the catalytic activity of all mouse trypsins is comparable. Human SPINK1 and its mouse ortholog inhibited mouse trypsins with comparable efficiency (K_D range 0.7-2.2 pM), with the sole exception of T7 trypsin, which was inhibited less effectively by the human inhibitor (K_D 21.9 pM). Characterization of four chronic pancreatitis-associated human SPINK1 mutations in the context of the mouse inhibitor revealed that the reactive-loop mutations R42N (human K41N) and I43M (human I42M) impaired SPINK1 binding to trypsin (K_D 60 nM and 47.5 pM, respectively), whereas mutations D35S (human N34S) and A56S (human P55S) had no impact on trypsin inhibition. Our results confirmed that high-affinity trypsin inhibition by SPINK1 is conserved in the mouse, and the functional consequences of human pancreatitis-associated SPINK1 mutations can be replicated in the mouse inhibitor.

Chronic Pancreatitis: The True Pathogenic Culprit within the SPINK1 N34S-Containing Haplotype Is No Longer at Large

A diverse range of loss-of-function variants in the SPINK1 gene (encoding pancreatic secretory trypsin inhibitor) has been identified in patients with chronic pancreatitis (CP). The haplotype harboring the SPINK1 c.101A>G (p.Asn34Ser or N34S) variant (rs17107315:T>C) is one of the most important heritable risk factors for CP as a consequence of its relatively high prevalence worldwide (population allele frequency ≈ 1%) and its considerable effect size (odds ratio ≈ 11). The causal variant responsible for this haplotype has been intensively investigated over the past two decades. The different hypotheses tested addressed whether the N34S missense variant has a direct impact on enzyme structure and function, whether c.101A>G could affect pre-mRNA splicing or mRNA stability, and whether another variant in linkage disequilibrium with c.101A>G might be responsible for the observed association with CP. Having reviewed the currently available genetic and experimental data, we conclude that c.-4141G>T (rs142703147:C>A), which disrupts a PTF1L-binding site within an evolutionarily conserved HNF1A-PTF1L cis-regulatory module located ∼4 kb upstream of the SPINK1 promoter, can be designated as the causal variant beyond reasonable doubt. This case illustrates the difficulties inherent in determining the identity of the causal variant underlying an initially identified disease association.

Defective binding of SPINK1 variants is an uncommon mechanism for impaired trypsin inhibition in chronic pancreatitis

The serine protease inhibitor Kazal type 1 (SPINK1) protects the pancreas from intrapancreatic trypsin activation that can lead to pancreatitis. Loss-of-function genetic variants of SPINK1 increase the risk for chronic pancreatitis, often by diminishing inhibitor expression or secretion. Variants that are secreted normally have been presumed to be pathogenic because of defective trypsin inhibition, but evidence has been lacking. Here, we report quantitative studies on the inhibition of human trypsins by wildtype SPINK1 and seven secreted missense variants. We found that tyrosine sulfation of human trypsins weakens binding of SPINK1 because of altered interactions with Tyr43 in the SPINK1 reactive loop. Using authentic sulfated human trypsins, we provide conclusive evidence that SPINK1 variants N34S, N37S, R65Q, and Q68R have unimpaired inhibitory activity, whereas variant P55S exhibits a small and clinically insignificant binding defect. In contrast, rare variants K41N and I42M that affect the reactive-site peptide bond of SPINK1 decrease inhibitor binding by 20,000- to 30,000-fold and three- to sevenfold, respectively. Taken together, the observations indicate that defective trypsin inhibition by SPINK1 variants is an uncommon mechanism in chronic pancreatitis. The results also strengthen the notion that a decline in inhibitor levels explains pancreatitis risk associated with the large majority of SPINK1 variants.

Splicing analysis of SLC40A1 missense variations and contribution to hemochromatosis type 4 phenotypes

Hemochromatosis type 4, or ferroportin disease, is considered as the second leading cause of primary iron overload after HFE-related hemochromatosis. The disease, which is predominantly associated with missense variations in the SLC40A1 gene, is characterized by wide clinical heterogeneity. We tested the possibility that some of the reported missense mutations, despite their positions within exons, cause splicing defects. Fifty-eight genetic variants were selected from the literature based on two criteria: a precise description of the nucleotide change and individual evidence of iron overload. The selected variants were investigated by different in silico prediction tools and prioritized for midigene splicing assays. Of the 15 variations tested in vitro, only two were associated with splicing changes. We confirm that the c.1402G>A transition (p.Gly468Ser) disrupts the exon 7 donor site, leading to the use of an exonic cryptic splicing site and the generation of a truncated reading frame. We observed, for the first time, that the p.Gly468Ser substitution has no effect on the ferroportin iron export function. We demonstrate alternative splicing of exon 5 in different cell lines and show that the c.430A>G (p.Asn144Asp) variant promotes exon 5 inclusion. This could be part of a gain-of-function mechanism. We conclude that splicing mutations rarely contribute to hemochromatosis type 4 phenotypes. An in-depth investigation of exon 5 auxiliary splicing sequences may help to elucidate the mechanism by which splicing regulatory proteins regulate the production of the full length SLC40A1 transcript and to clarify its physiological importance.

First estimate of the scale of canonical 5' splice site GT>GC variants capable of generating wild-type transcripts

It has long been known that canonical 5' splice site (5'SS) GT>GC variants may be compatible with normal splicing. However, to date, the actual scale of canonical 5'SSs capable of generating wild-type transcripts in the case of GT>GC substitutions remains unknown. Herein, combining data derived from a meta-analysis of 45 human disease-causing 5'SS GT>GC variants and a cell culture-based full-length gene splicing assay of 103 5'SS GT>GC substitutions, we estimate that ~15-18% of canonical GT 5'SSs retain their capacity to generate between 1% and 84% normal transcripts when GT is substituted by GC. We further demonstrate that the canonical 5'SSs in which substitution of GT by GC-generated normal transcripts exhibit stronger complementarity to the 5' end of U1 snRNA than those sites whose substitutions of GT by GC did not lead to the generation of normal transcripts. We also observed a correlation between the generation of wild-type transcripts and a milder than expected clinical phenotype but found that none of the available splicing prediction tools were capable of reliably distinguishing 5'SS GT>GC variants that generated wild-type transcripts from those that did not. Our findings imply that 5'SS GT>GC variants in human disease genes may not invariably be pathogenic.

The phenotypic spectrum of WWOX-related disorders: 20 additional cases of WOREE syndrome and review of the literature

Purpose

Germline WWOX pathogenic variants have been associated with disorder of sex differentiation (DSD), spinocerebellar ataxia (SCA), and WWOX-related epileptic encephalopathy (WOREE syndrome). We review clinical and molecular data on WWOX-related disorders, further describing WOREE syndrome and phenotype/genotype correlations.

Methods

We report clinical and molecular findings in 20 additional patients from 18 unrelated families with WOREE syndrome and biallelic pathogenic variants in the WWOX gene. Different molecular screening approaches were used (quantitative polymerase chain reaction/multiplex ligation-dependent probe amplification [qPCR/MLPA], array comparative genomic hybridization [array-CGH], Sanger sequencing, epilepsy gene panel, exome sequencing), genome sequencing.

Results

Two copy-number variations (CNVs) or two single-nucleotide variations (SNVs) were found respectively in four and nine families, with compound heterozygosity for one SNV and one CNV in five families. Eight novel missense pathogenic variants have been described. By aggregating our patients with all cases reported in the literature, 37 patients from 27 families with WOREE syndrome are known. This review suggests WOREE syndrome is a very severe epileptic encephalopathy characterized by absence of language development and acquisition of walking, early-onset drug-resistant seizures, ophthalmological involvement, and a high likelihood of premature death. The most severe clinical presentation seems to be associated with null genotypes.

Conclusion

Germline pathogenic variants in WWOX are clearly associated with a severe early-onset epileptic encephalopathy. We report here the largest cohort of individuals with WOREE syndrome.