Toward a clinical diagnostic pipeline for SPINK1 intronic variants

Heterozygous Spink1 Deficiency Promotes Trypsin-dependent Chronic Pancreatitis in Mice

BACKGROUND & AIMS

Heterozygous SPINK1 mutations are strong risk factors for chronic pancreatitis in humans, yet heterozygous disruption of mouse Spink1 yielded no pancreatic phenotype. To resolve this contradiction, we used CRISPR/Cas9-mediated genome editing to generate heterozygous Spink1-deleted mice (Spink1-KOhet) in the C57BL/6N strain and studied the effect of this allele in trypsin-independent and trypsin-dependent pancreatitis models.

METHODS

We investigated severity of acute pancreatitis and progression to chronic pancreatitis in Spink1-KOhet mice after transient (10 injections) and prolonged (2 × 8 injections) cerulein hyperstimulation. We crossed Spink1-KOhet mice with T7D23A and T7D22N,K24R mice that carry strongly autoactivating trypsinogen mutants and exhibit spontaneous chronic pancreatitis.

RESULTS

Prolonged but not transient cerulein stimulation resulted in increased intrapancreatic trypsin activity and more severe acute pancreatitis in Spink1-KOhet mice relative to the C57BL/6N control strain. After the acute episode, Spink1-KOhet mice developed progressive disease with chronic pancreatitis-like features, whereas C57BL/6N mice recovered rapidly. Trypsinogen mutant mice carrying the Spink1-KOhet allele exhibited strikingly more severe chronic pancreatitis than the respective parent strains.

CONCLUSIONS

Heterozygous Spink1 deficiency caused more severe acute pancreatitis after prolonged cerulein stimulation and promoted chronic pancreatitis after the cerulein-induced acute episode, in two strains of trypsinogen mutant mice with spontaneous disease. In contrast, acute pancreatitis induced with limited cerulein hyperstimulation was unaffected by heterozygous Spink1 deletion, in agreement with recent observations that trypsin activity does not mediate pathologic responses in this model. Taken together, the findings strongly support the notion that loss-of-function SPINK1 mutations in humans increase chronic pancreatitis risk in a trypsin-dependent manner.

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.

Splicing Outcomes of 5' Splice Site GT>GC Variants That Generate Wild-Type Transcripts Differ Significantly Between Full-Length and Minigene Splicing Assays

Combining data derived from a meta-analysis of human disease-associated 5' splice site GT>GC (i.e., +2T>C) variants and a cell culture-based full-length gene splicing assay (FLGSA) of forward engineered +2T>C substitutions, we recently estimated that ∼15-18% of +2T>C variants can generate up to 84% wild-type transcripts relative to their wild-type counterparts. Herein, we analyzed the splicing outcomes of 20 +2T>C variants that generate some wild-type transcripts in two minigene assays. We found a high discordance rate in terms of the generation of wild-type transcripts, not only between FLGSA and the minigene assays but also between the different minigene assays. In the pET01 context, all 20 wild-type minigene constructs generated the expected wild-type transcripts; of the 20 corresponding variant minigene constructs, 14 (70%) generated wild-type transcripts. In the pSPL3 context, only 18 of the 20 wild-type minigene constructs generated the expected wild-type transcripts whereas 8 of the 18 (44%) corresponding variant minigene constructs generated wild-type transcripts. Thus, in the context of a particular type of variant, we raise awareness of the limitations of minigene splicing assays and emphasize the importance of sequence context in regulating splicing. Whether or not our findings apply to other types of splice-altering variant remains to be investigated.

Functional Roles of SPINK1 in Cancers

Serine Peptidase Inhibitor Kazal Type 1 (SPINK1) is a secreted protein known as a protease inhibitor of trypsin in the pancreas. However, emerging evidence shows its function in promoting cancer progression in various types of cancer. SPINK1 modulated tumor malignancies and induced the activation of the downstream signaling of epidermal growth factor receptor (EGFR) in cancer cells, due to the structural similarity with epidermal growth factor (EGF). The discoverable SPINK1 somatic mutations, expressional signatures, and prognostic significances in various types of cancer have attracted attention as a cancer biomarker in clinical applications. Emerging findings further clarify the direct and indirect biological effects of SPINK1 in regulating cancer proliferation, metastasis, drug resistance, transdifferentiation, and cancer stemness, warranting the exploration of the SPINK1-mediated molecular mechanism to identify a therapeutic strategy. In this review article, we first integrate the transcriptomic data of different types of cancer with clinical information and recent findings of SPINK1-mediated malignant phenotypes. In addition, a comprehensive summary of SPINK1 expression in a pan-cancer panel and individual cell types of specific organs at the single-cell level is presented to indicate the potential sites of tumorigenesis, which has not yet been reported. This review aims to shed light on the roles of SPINK1 in cancer and provide guidance and potential directions for scientists in this field.

Clinical interpretation of SPINK1 and CTRC variants in pancreatitis

Since the description of the SPINK1 gene encoding the serine protease inhibitor Kazal type 1 and the CTRC gene encoding the Chymotrypsin C as being involved in chronic pancreatitis, more than 56 SPINK1 and 87 CTRC variants have been reported. Assessing the clinical relevance of SPINK1 and CTRC variants is often complicated in the absence of functional evidence and interpretation of rare variants is not very easy in clinical practice. The aim of this study was to review the different variants identified in these two genes and to classify them according to their degree of damaging effect. This classification was based on the results of in vitro experiments, in silico analysis using different prediction tools, and on population data, in comparing the allelic frequency of each variant in patients with pancreatitis and in unaffected control individuals. This review should help geneticists and clinicians in charge of patient's care and genetic counseling to interpret the results of genetic studies.

5' splice site GC>GT and GT>GC variants differ markedly in terms of their functionality and pathogenicity

In the human genome, most 5' splice sites (~99%) employ the canonical GT dinucleotide whereas a small minority (~1%) use the noncanonical GC dinucleotide. The functionality and pathogenicity of 5' splice site GT>GC (+2T>C) variants have been extensively studied but we know very little about 5' splice site GC>GT (+2C>T) variants. Herein, we have addressed this deficiency by performing a meta-analysis of reported +2C>T "pathogenic" variants together with a functional analysis of engineered +2C>T substitutions using a cell culture-based full-length gene splicing assay. Our results establish proof of concept that +2C>T variants are qualitatively different from +2T>C variants in terms of their functionality and suggest that, in sharp contrast to +2T>C variants, most if not all +2C>T variants have no pathological relevance. Our findings have important implications for interpreting the clinical relevance of +2C>T variants and understanding the evolutionary switching between GT and GC 5' splice sites in mammalian genomes.

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.