Identification of a functional enhancer variant within the chronic pancreatitis-associated SPINK1 c.101A>G (p.Asn34Ser)-containing haplotype

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.

The PRSS3P2 and TRY7 deletion copy number variant modifies risk for chronic pancreatitis

BACKGROUND

PRSS1 and PRSS2 constitute the only functional copies of a tandemly-arranged five-trypsinogen-gene cluster (i.e., PRSS1, PRSS3P1, PRSS3P2, TRY7 and PRSS2) on chromosome 7q35. Variants in PRSS1 and PRSS2, including missense and copy number variants (CNVs), have been reported to predispose to or protect against chronic pancreatitis (CP). We wondered whether a common trypsinogen pseudogene deletion CNV (that removes two of the three trypsinogen pseudogenes, PRSS3P2 and TRY7) might be associated with CP causation/predisposition.

METHODS

We analyzed the common PRSS3P2 and TRY7 deletion CNV in a total of 1536 CP patients and 3506 controls from France, Germany, India and Japan by means of quantitative fluorescent multiplex polymerase chain reaction.

RESULTS

We demonstrated that the deletion CNV variant was associated with a protective effect against CP in the French, German and Japanese cohorts whilst a trend toward the same association was noted in the Indian cohort. Meta-analysis under a dominant model yielded a pooled odds ratio (OR) of 0.68 (95% confidence interval (CI) 0.52-0.89; p = 0.005) whereas an allele-based meta-analysis yielded a pooled OR of 0.84 (95% CI 0.77-0.92; p = 0.0001). This protective effect is explicable by reference to the recent finding that the still functional PRSS3P2/TRY7 pseudogene enhancers upregulate pancreatic PRSS2 expression.

CONCLUSIONS

The common PRSS3P2 and TRY7 deletion CNV was associated with a reduced risk for CP. This finding provides additional support for the emerging view that dysregulated PRSS2 expression represents a discrete mechanism underlying CP predisposition or protection.

Expanding ACMG variant classification guidelines into a general framework

BACKGROUND

The American College of Medical Genetics and Genomics (ACMG)-recommended five variant classification categories (pathogenic, likely pathogenic, uncertain significance, likely benign, and benign) have been widely used in medical genetics. However, these guidelines are fundamentally constrained in practice owing to their focus upon Mendelian disease genes and their dichotomous classification of variants as being either causal or not. Herein, we attempt to expand the ACMG guidelines into a general variant classification framework that takes into account not only the continuum of clinical phenotypes, but also the continuum of the variants' genetic effects, and the different pathological roles of the implicated genes.

MAIN BODY

As a disease model, we employed chronic pancreatitis (CP), which manifests clinically as a spectrum from monogenic to multifactorial. Bearing in mind that any general conceptual proposal should be based upon sound data, we focused our analysis on the four most extensively studied CP genes, PRSS1, CFTR, SPINK1 and CTRC. Based upon several cross-gene and cross-variant comparisons, we first assigned the different genes to two distinct categories in terms of disease causation: CP-causing (PRSS1 and SPINK1) and CP-predisposing (CFTR and CTRC). We then employed two new classificatory categories, "predisposing" and "likely predisposing", to replace ACMG's "pathogenic" and "likely pathogenic" categories in the context of CP-predisposing genes, thereby classifying all pathologically relevant variants in these genes as "predisposing". In the case of CP-causing genes, the two new classificatory categories served to extend the five ACMG categories whilst two thresholds (allele frequency and functional) were introduced to discriminate "pathogenic" from "predisposing" variants.

CONCLUSION

Employing CP as a disease model, we expand ACMG guidelines into a five-category classification system (predisposing, likely predisposing, uncertain significance, likely benign, and benign) and a seven-category classification system (pathogenic, likely pathogenic, predisposing, likely predisposing, uncertain significance, likely benign, and benign) in the context of disease-predisposing and disease-causing genes, respectively. Taken together, the two systems constitute a general variant classification framework that, in principle, should span the entire spectrum of variants in any disease-related gene. The maximal compliance of our five-category and seven-category classification systems with the ACMG guidelines ought to facilitate their practical application.

Colocalization analysis of pancreas eQTLs with risk loci from alcoholic and novel non-alcoholic chronic pancreatitis GWAS suggests potential disease causing mechanisms

BACKGROUND

Previous genome-wide association studies (GWAS) identified genome-wide significant risk loci in chronic pancreatitis and investigated underlying disease causing mechanisms by simple overlaps with expression quantitative trait loci (eQTLs), a procedure which may often result in false positive conclusions.

METHODS

We conducted a GWAS in 584 non-alcoholic chronic pancreatitis (NACP) patients and 6040 healthy controls. Next, we applied Bayesian colocalization analysis of identified genome-wide significant risk loci from both, our recently published alcoholic chronic pancreatitis (ACP) and the novel NACP dataset, with pancreas eQTLs from the GTEx V8 European cohort to prioritize candidate causal genes and extracted credible sets of shared causal variants.

RESULTS

Variants at the CTRC (p = 1.22 × 10^-21) and SPINK1 (p = 6.59 × 10^-47) risk loci reached genome-wide significance in NACP. CTRC risk variants colocalized with CTRC eQTLs in ACP (PP4 = 0.99, PP4/PP3 = 95.51) and NACP (PP4 = 0.99, PP4/PP3 = 95.46). For both diseases, the 95% credible set of shared causal variants consisted of rs497078 and rs545634. CLDN2-MORC4 risk variants colocalized with CLDN2 eQTLs in ACP (PP4 = 0.98, PP4/PP3 = 42.20) and NACP (PP4 = 0.67, PP4/PP3 = 7.18), probably driven by the shared causal variant rs12688220.

CONCLUSIONS

A shared causal CTRC risk variant might unfold its pathogenic effect in ACP and NACP by reducing CTRC expression, while the CLDN2-MORC4 shared causal variant rs12688220 may modify ACP and NACP risk by increasing CLDN2 expression.

The landscape of GWAS validation; systematic review identifying 309 validated non-coding variants across 130 human diseases

Background

The remarkable growth of genome-wide association studies (GWAS) has created a critical need to experimentally validate the disease-associated variants, 90% of which involve non-coding variants.

Methods

To determine how the field is addressing this urgent need, we performed a comprehensive literature review identifying 36,676 articles. These were reduced to 1454 articles through a set of filters using natural language processing and ontology-based text-mining. This was followed by manual curation and cross-referencing against the GWAS catalog, yielding a final set of 286 articles.

Results

We identified 309 experimentally validated non-coding GWAS variants, regulating 252 genes across 130 human disease traits. These variants covered a variety of regulatory mechanisms. Interestingly, 70% (215/309) acted through cis-regulatory elements, with the remaining through promoters (22%, 70/309) or non-coding RNAs (8%, 24/309). Several validation approaches were utilized in these studies, including gene expression (n = 272), transcription factor binding (n = 175), reporter assays (n = 171), in vivo models (n = 104), genome editing (n = 96) and chromatin interaction (n = 33).

Conclusions

This review of the literature is the first to systematically evaluate the status and the landscape of experimentation being used to validate non-coding GWAS-identified variants. Our results clearly underscore the multifaceted approach needed for experimental validation, have practical implications on variant prioritization and considerations of target gene nomination. While the field has a long way to go to validate the thousands of GWAS associations, we show that progress is being made and provide exemplars of validation studies covering a wide variety of mechanisms, target genes, and disease areas.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01216-w.

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.

Regulatory SNPs: Altered Transcription Factor Binding Sites Implicated in Complex Traits and Diseases

The vast majority of the genetic variants (mainly SNPs) associated with various human traits and diseases map to a noncoding part of the genome and are enriched in its regulatory compartment, suggesting that many causal variants may affect gene expression. The leading mechanism of action of these SNPs consists in the alterations in the transcription factor binding via creation or disruption of transcription factor binding sites (TFBSs) or some change in the affinity of these regulatory proteins to their cognate sites. In this review, we first focus on the history of the discovery of regulatory SNPs (rSNPs) and systematized description of the existing methodical approaches to their study. Then, we brief the recent comprehensive examples of rSNPs studied from the discovery of the changes in the TFBS sequence as a result of a nucleotide substitution to identification of its effect on the target gene expression and, eventually, to phenotype. We also describe state-of-the-art genome-wide approaches to identification of regulatory variants, including both making molecular sense of genome-wide association studies (GWAS) and the alternative approaches the primary goal of which is to determine the functionality of genetic variants. Among these approaches, special attention is paid to expression quantitative trait loci (eQTLs) analysis and the search for allele-specific events in RNA-seq (ASE events) as well as in ChIP-seq, DNase-seq, and ATAC-seq (ASB events) data.

A Hypothesized Mechanism for Chronic Pancreatitis Caused by the N34S Mutation of Serine Protease Inhibitor Kazal-Type 1 Based on Conformational Studies

Purpose

Although strongly related, the pathophysiological effect of the N34S mutation in the serine protease inhibitor Kazal type 1 (SPINK1) in chronic pancreatitis is still unknown. In this study, we investigate the conformational space of the human cationic trypsin-serine protease inhibitor complex.

Methods

Simulations with molecular dynamics, replica exchange, and transition pathway methods are used.

Results

Two main binding states of the inhibitor to the complex were found, which explicitly relate the influence of the mutation site to conformational changes in the active site of trypsin.

Conclusion

Based on our result, a hypothesis is formulated that explains the development of chronic pancreatitis through accelerated digestion of the mutant by trypsin.

Scale and Scope of Gene-Alcohol Interactions in Chronic Pancreatitis: A Systematic Review

BACKGROUND

Excessive alcohol consumption has long been known to be the primary cause of chronic pancreatitis (CP) but genetic risk factors have been increasingly identified over the past 25 years. The scale and scope of gene-alcohol interactions in CP nevertheless remain unclear.

METHODS

All studies that had obtained genetic variant data concurrently on alcoholic CP (ACP) patients, non-ACP (NACP) patients and normal controls were collated. Employing normal controls as a common baseline, paired OR_ACP and OR_NACP (odds ratios associated with ACP and NACP, respectively) values were calculated and used to assess gene-alcohol interactions.

RESULTS

Thirteen variants involving PRSS1, SPINK1, CTRC, CLDN2, CPA1, CEL and CTRB1-CTRB2, and varying from very rare to common, were collated. Seven variants had an OR_ACP > OR_NACP, which was regarded as an immediate indicator of gene-alcohol interactions in CP. Variants with an OR_ACP < OR_NACP were also found to interact with alcohol consumption by virtue of their impact on age at first pancreatitis symptoms in ACP.

CONCLUSIONS

This study revealed evidence for extensive gene-alcohol interactions in CP. Our findings lend support to the hypothesis that alcohol affects the expression of genetically determined CP and highlight a predominant role of weak-effect variants in the development of ACP.

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.

Acute Pancreatitis: Genetic Risk and Clinical Implications

Acute pancreatitis (AP) is one of the most common gastroenterological indications for emergency admittance and hospitalization. Gallstones, alcohol consumption or the presence of additional initiating factors give rise to a disease with a diverse clinical appearance and a hard-to predict course of progression. One major challenge in the treatment of AP patients is the early identification of patients at risk for the development of systemic complications and organ failure. In addition, 20%-30% of patients with a first episode of AP later experience progress to recurrent or chronic disease. Complex gene-environment interactions have been identified to play a role in the pathogenesis of pancreatitis, but so far no predictive genetic biomarkers could be implemented into the routine clinical care of AP patients. The current review explains common and rare etiologies of acute pancreatitis with emphasis on underlying genetic aberrations and ensuing clinical management.

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.

Complex Genetics in Pancreatitis: Insights Gained From a New Candidate Locus Panel

OBJECTIVES

Chronic pancreatitis is the end stage of a pathologic inflammatory syndrome with multiple etiological factors, including genetic. We hypothesized that some pancreatitis etiology originates in pancreatic acinar or duct cells and requires both injury and compensatory mechanism failure.

METHODS

One hundred pancreatitis patients were assessed using a DNA sequencing panel for pancreatitis. Cooccurrence of variants within and between genes was measured. Gene coexpression was confirmed via published single-cell RNA sequencing.

RESULTS

One hundred and twenty-one variants were identified in 2 or more patients, 15 of which were enriched compared with reference populations. Single cell RNA-sequencing data verified coexpression of GGT1, CFTR, and PRSS1 in duct cells, PRSS1, CPA1, CEL, CTRC, and SPINK1 in acinar cells, and UBR1 in both. Multiple-risk variants with injury/stress effects (CEL, CFTR, CPA1, PRSS1) and impaired cell protection (CTRC, GGT1, SPINK1, UBR1) cooccur within duct cells, acinar cells, or both.

CONCLUSIONS

Pancreatitis is a complex disorder with genetic interactions across genes and cell types. These findings suggest a new, non-Mendelian genetic risk/etiology paradigm where a combination of nonpathogenic genetic risk variants in groups of susceptibility genes and injury/dysfunction response genes contribute to acquired pancreatic disease.

Toward a clinical diagnostic pipeline for SPINK1 intronic variants

Background

The clinical significance of SPINK1 intronic variants in chronic pancreatitis has been previously assessed by various approaches including a cell culture-based full-length gene assay. A close correlation between the results of this assay and in silico splicing prediction was apparent. However, until now, a clinical diagnostic pipeline specifically designed to classify SPINK1 intronic variants accurately and efficiently has been lacking. Herein, we present just such a pipeline and explore its efficacy and potential utility in potentiating the classification of newly described SPINK1 intronic variants.

Results

We confirm a close correlation between in silico splicing prediction and results from the cell culture-based full-length gene assay in the context of three recently reported pathogenic SPINK1 intronic variants. We then integrated in silico splicing prediction and the full-length gene assay into a stepwise approach and tested its utility in the classification of two novel datasets of SPINK1 intronic variants. The first dataset comprised 16 deep intronic variants identified in 52 genetically unexplained Chinese chronic pancreatitis patients by sequencing the entire intronic sequence of the SPINK1 gene. The second dataset comprised five novel rare proximal intronic variants identified through the routine analysis of the SPINK1 gene in French pancreatitis patients. Employing a minor allele frequency of > 5% as a population frequency filter, 6 of the 16 deep intronic variants were immediately classified as benign. In silico prediction of the remaining ten deep intronic variants and the five rare proximal intronic variants with respect to their likely impact on splice site selection suggested that only one proximal intronic variant, c.194 + 5G > A, was likely to be of functional significance. Employing the cell culture-based full-length gene assay, we functionally analyzed c.194 + 5G > A, together with seven predicted non-functional variants, thereby validating their predicted effects on splicing in all cases.

Conclusions

We demonstrated the accuracy and efficiency of in silico prediction in combination with the cell culture-based full-length gene assay for the classification of SPINK1 intronic variants. Based upon these findings, we propose an operational pipeline for classifying SPINK1 intronic variants in the clinical diagnostic setting.

Electronic supplementary material

The online version of this article (10.1186/s40246-019-0193-7) contains supplementary material, which is available to authorized users.

SPINK1, PRSS1, CTRC, and CFTR Genotypes Influence Disease Onset and Clinical Outcomes in Chronic Pancreatitis

Objectives

Rare pathogenic variants in the SPINK1, PRSS1, CTRC, and CFTR genes have been strongly associated with a risk of developing chronic pancreatitis (CP). However, their potential impact on the age of disease onset and clinical outcomes, as well as their potential interactions with environmental risk factors, remain unclear. These issues are addressed here in a large Chinese CP cohort.

Methods

We performed targeted next-generation sequencing of the four CP-associated genes in 1061 Han Chinese CP patients and 1196 controls. To evaluate gene–environment interactions, the patients were divided into three subgroups, idiopathic CP (ICP; n = 715), alcoholic CP (ACP; n = 206), and smoking-associated CP (SCP; n = 140). The potential impact of rare pathogenic variants on the age of onset of CP and clinical outcomes was evaluated using the Kaplan–Meier model.

Results

We identified rare pathogenic genotypes involving the SPINK1, PRSS1, CTRC, and/or CFTR genes in 535 (50.42%) CP patients but in only 71 (5.94%) controls (odds ratio = 16.12; P < 0.001). Mutation-positive patients had significantly earlier median ages at disease onset and at diagnosis of pancreatic stones, diabetes mellitus and steatorrhea than mutation-negative ICP patients. Pathogenic genotypes were present in 57.1, 39.8, and 32.1% of the ICP, ACP, and SCP patients, respectively, and influenced age at disease onset and clinical outcomes in all subgroups.

Conclusions

We provide evidence that rare pathogenic variants in the SPINK1, PRSS1, CTRC, and CFTR genes significantly influence the age of onset and clinical outcomes of CP. Extensive gene–environment interactions were also identified.

Chronic pancreatitis: an update on genetic risk factors

PURPOSE OF REVIEW

Genetic mutations in genes within and outside of the trypsin-dependent pathologic pathway have been found to be associated with chronic pancreatitis. This review highlights recent developments.

RECENT FINDINGS

CTRB1-CTRB2 has been identified as a new risk locus for chronic pancreatitis and the disease mechanism may involve trypsin degradation. Misfolding mutations in PRSS1, CPA1, and CEL, as well as environmental stress factors like tobacco and alcohol can trigger endoplasmic reticulum stress (ER-Stress).

SUMMARY

Protein misfolding as well as enzyme activity changes due to altered autoactivation, intracellular degradation, or enzyme inhibition represent the most important pathological mechanisms of chronic pancreatitis to date. Analysis of composite risk patterns by next-generation sequencing will help elucidate complex gene interactions and identify new potential therapeutic targets.