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.

The p.E152K-STIM1 mutation deregulates Ca2+ signaling contributing to chronic pancreatitis.

Since deregulation of intracellular Ca2+ can lead to intracellular trypsin activation, and stromal interaction molecule-1 (STIM1) protein is the main regulator of Ca2+ homeostasis in pancreatic acinar cells, we explored the Ca2+ signaling in 37 STIM1 variants found in three pancreatitis patient cohorts. Extensive functional analysis of one particular variant, p.E152K, identified in three patients, provided a plausible link between dysregulated Ca2+ signaling within pancreatic acinar cells and chronic pancreatitis susceptibility. Specifically, p.E152K, located within the STIM1 EF-hand and sterile α-motif domain, increased the release of Ca2+ from the endoplasmic reticulum in patient-derived fibroblasts and transfected HEK293T cells. This event was mediated by altered STIM1-sarco/endoplasmic reticulum calcium transport ATPase (SERCA) conformational change and enhanced SERCA pump activity leading to increased store-operated Ca2+ entry (SOCE). In pancreatic AR42J cells expressing the p.E152K variant, Ca2+ signaling perturbations correlated with defects in trypsin activation and secretion, and increased cytotoxicity after cholecystokinin stimulation.This article has an associated First Person interview with the first author of the paper.

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.

Classification of PRSS1 variants responsible for chronic pancreatitis: An expert perspective from the Franco-Chinese GREPAN study group.

BACKGROUND: PRSS1 was the first reported chronic pancreatitis (CP) gene. The existence of both gain-of-function (GoF) and gain-of-proteotoxicity (GoP) pathological PRSS1 variants, together with the fact that PRSS1 variants have been identified in CP subtypes spanning the range from monogenic to multifactorial, has made the classification of PRSS1 variants very challenging. METHODS: All currently reported PRSS1 variants (derived primarily from two databases) were manually reviewed with respect to their clinical genetics, functional analysis and population allele frequency. They were classified by variant type and pathological mechanism within the framework of our recently proposed ACMG/AMP guidelines-based seven-category system. RESULTS: The total number of distinct germline PRSS1 variants included for analysis was 100, comprising 3 copy number variants (CNVs), 12 5' and 3' variants, 19 intronic variants, 5 nonsense variants, 1 frameshift deletion variant, 6 synonymous variants, 1 in-frame duplication, 3 gene conversions and 50 missense variants. Based upon a combination of clinical genetic and functional analysis, population data and in silico analysis, we classified 26 variants (all 3 CNVs, the in-frame duplication, all 3 gene conversions and 19 missense) as "pathogenic", 3 variants (missense) as "likely pathogenic", 5 variants (four missense and one promoter) as "predisposing", 13 variants (all missense) as "unknown significance", 2 variants (missense) as "likely benign", and all remaining 51 variants as "benign". CONCLUSIONS: We describe an expert classification of the 100 PRSS1 variants reported to date. The results have immediate implications for reclassifying many ClinVar-registered PRSS1 variants as well as providing optimal guidelines/standards for reporting PRSS1 variants.

Identification of protease-sensitive but not misfolding PNLIP variants in familial and hereditary pancreatitis.

Mutations in the PNLIP gene have recently been implicated in chronic pancreatitis. Several PNLIP missense variants have been reported to cause protein misfolding and endoplasmic reticulum stress although genetic evidence supporting their association with chronic pancreatitis is currently lacking. Protease-sensitive PNLIP missense variants have also been associated with early-onset chronic pancreatitis although the underlying pathological mechanism remains enigmatic. Herein, we provide new evidence to support the association of protease-sensitive PNLIP variants (but not misfolding PNLIP variants) with pancreatitis. Specifically, we identified protease-sensitive PNLIP variants in 5 of 373 probands (1.3%) with a positive family history of pancreatitis. The protease-sensitive variants, p.F300L and p.I265R, were found to segregate with the disease in three families, including one exhibiting a classical autosomal dominant inheritance pattern. Consistent with previous findings, protease-sensitive variant-positive patients were often characterized by early-onset disease and invariably experienced recurrent acute pancreatitis, although none has so far developed chronic pancreatitis.

Genetic and functional analysis of chymotrypsin-like protease (CTRL) in chronic pancreatitis.

BACKGROUND: Genetic predisposition is crucial in the pathogenesis of early-onset chronic pancreatitis (CP). So far, several genetic alterations have been identified as risk factors, predominantly in genes encoding digestive enzymes. However, many early-onset CP cases have no identified underlying cause. Chymotrypsins are a family of serine proteases that can cleave trypsinogen and lead to its degradation. Because genetic alterations in the chymotrypsins CTRC, CTRB1, and CTRB2 are associated with CP, we genetically and functionally investigated chymotrypsin-like protease (CTRL) as a potential risk factor. METHODS: We screened 1005 non-alcoholic CP patients and 1594 controls for CTRL variants by exome sequencing. We performed Western blots and activity assays to analyse secretion and proteolytic activity. We measured BiP mRNA expression to investigate the potential impact of identified alterations on endoplasmic reticulum (ER) stress. RESULTS: We identified 13 heterozygous non-synonymous CTRL variants: five exclusively in patients and three only in controls. Functionality was unchanged in 6/13 variants. Four alterations showed normal secretion but reduced (p.G20S, p.G56S, p.G61S) or abolished (p.S208F) activity. Another three variants (p.C201Y, p.G215R and p.C220G) were not secreted and already showed reduced or no activity intracellularly. However, intracellular retention did not lead to ER stress. CONCLUSION: We identified several CTRL variants, some showing potent effects on protease function and secretion. We observed these effects in variants found in patients and controls, and CTRL loss-of-function variants were not significantly more common in patients than controls. Therefore, CTRL is unlikely to play a relevant role in the development of CP.

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.

Frameshift coding sequence variants in the LPL gene: identification of two novel events and exploration of the genotype-phenotype relationship for variants reported to date.

BACKGROUND: Lipoprotein lipase (LPL) is the rate-limiting enzyme for triglyceride hydrolysis. Homozygous or compound heterozygous LPL variants cause autosomal recessive familial chylomicronemia syndrome (FCS), whereas simple heterozygous LPL variants are associated with hypertriglyceridemia (HTG) and HTG-related disorders. LPL frameshift coding sequence variants usually cause complete functional loss of the affected allele, thereby allowing exploration of the impact of different levels of LPL function in human disease. METHODS: All exons and flanking intronic regions of LPL were Sanger sequenced in patients with HTG-related acute pancreatitis (HTG-AP) or HTG-AP in pregnancy. Previously reported LPL frameshift coding sequence variants were collated from the Human Gene Mutation Database and through PubMed keyword searching. Original reports were manually evaluated for the following information: zygosity status of the variant, plasma LPL activity of the variant carrier, disease referred for genetic analysis, patient's age at genetic analysis, and patient's disease history. SpliceAI was employed to predict the potential impact of collated variants on splicing. RESULTS: Two novel rare variants were identified, and 53 known LPL frameshift coding sequence variants were collated. Of the 51 variants informative for zygosity, 30 were simple heterozygotes, 12 were homozygotes, and 9 were compound heterozygotes. Careful evaluation of the 55 variants with respect to their clinical and genetic data generated several interesting findings. First, we conclude that 6-7% residual LPL function could significantly delay the age of onset of FCS and reduce the prevalence of FCS-associated syndromes. Second, whereas a large majority of LPL frameshift coding sequence variants completely disrupt gene function through their "frameshift" nature, a small fraction of these variants may act wholly or partly as "in-frame" variants, leading to the generation of protein products with some residual LPL function. Third, we identified two candidate LPL frameshift coding sequence variants that may retain residual function based on genotype-phenotype correlation or SpliceAI-predicted data. CONCLUSIONS: This study reported two novel LPL variants and yielded new insights into the genotype-phenotype relationship as it pertains to LPL frameshift coding sequence variants.

Functionally deficient TRPV6 variants contribute to hereditary and familial chronic pancreatitis.

The recent discovery of TRPV6 as a pancreatitis susceptibility gene served to identify a novel mechanism of chronic pancreatitis (CP) due to Ca2+ dysregulation. Herein, we analyzed TRPV6 in 81 probands with hereditary CP (HCP), 204 probands with familial CP (FCP), and 462 patients with idiopathic CP (ICP) by targeted next-generation sequencing. We identified 25 rare nonsynonymous TRPV6 variants, 18 of which had not been previously reported. All 18 variants were characterized by a Ca2+ imaging assay, with 8 being identified as functionally deficient. Evaluation of functionally deficient variants in the three CP cohorts revealed two novel findings: (i) functionally deficient TRPV6 variants appear to occur more frequently in HCP/FCP patients than in ICP patients (3.2% vs. 1.5%) and (ii) functionally deficient TRPV6 variants found in HCP and FCP probands appear to be more frequently coinherited with known risk variants in SPINK1, CTRC, and/or CFTR than those found in ICP patients (66.7% vs 28.6%). Additionally, genetic analysis of available HCP and FCP family members revealed complex patterns of inheritance in some families. Our findings confirm that functionally deficient TRPV6 variants represent an important contributor to CP. Importantly, functionally deficient TRPV6 variants account for a significant proportion of cases of HCP/FCP.

Trypsinogen (PRSS1 and PRSS2) gene dosage correlates with pancreatitis risk across genetic and transgenic studies: a systematic review and re-analysis.

Trypsinogen (PRSS1, PRSS2) copy number gains and regulatory variants have both been proposed to elevate pancreatitis risk through a gene dosage effect (i.e., by increasing the expression of wild-type protein). However, to date, their impact on pancreatitis risk has not been thoroughly evaluated whilst the underlying pathogenic mechanisms remain to be explicitly investigated in mouse models. Genetic studies of the rare trypsinogen duplication and triplication copy number variants (CNVs), and the common rs10273639C variant, were collated from PubMed and/or ClinVar. Mouse studies that analyzed the influence of a transgenically expressed wild-type human PRSS1 or PRSS2 gene on the development of pancreatitis were identified from PubMed. The genetic effects of the different risk genotypes, in terms of odds ratios, were calculated wherever appropriate. The genetic effects of the rare trypsinogen duplication and triplication CNVs were also evaluated by reference to their associated disease subtypes. We demonstrate a positive correlation between increased trypsinogen gene dosage and pancreatitis risk in the context of the rare duplication and triplication CNVs, and between the level of trypsinogen expression and disease risk in the context of the heterozygous and homozygous rs10273639C-tagged genotypes. We retrospectively identify three mouse transgenic studies that are informative in relation to the pathogenic mechanism underlying the trypsinogen gene dosage effect in pancreatitis. Trypsinogen gene dosage correlates with pancreatitis risk across genetic and transgenic studies, highlighting the fundamental role of dysregulated expression of wild-type trypsinogen in the etiology of pancreatitis. Specifically downregulating trypsinogen expression in the pancreas may serve as a potential therapeutic and/or prevention strategy for pancreatitis.

Variants in the pancreatic CUB and zona pellucida-like domains 1 (CUZD1) gene in early-onset chronic pancreatitis - A possible new susceptibility gene.

OBJECTIVE: Non-alcoholic chronic pancreatitis (NACP) frequently develops in the setting of genetic susceptibility associated with alterations in genes that are highly expressed in the pancreas. However, the genetic basis of NACP remains unresolved in a significant number of patients warranting a search for further risk genes. DESIGN: We analyzed CUZD1, which encodes the CUB and zona pellucida-like domains 1 protein that is found in high levels in pancreatic acinar cells. We sequenced the coding region in 1163 European patients and 2018 European controls. In addition, we analyzed 297 patients and 1070 controls from Japan. We analyzed secretion of wild-type and mutant CUZD1 from transfected cells using Western blotting. RESULTS: In the European cohort, we detected 30 non-synonymous variants. Using different prediction tools (SIFT, CADD, PROVEAN, PredictSNP) or the combination of these tools, we found accumulation of predicted deleterious variants in patients (p-value range 0.002-0.013; OR range 3.1-5.2). No association was found in the Japanese cohort, in which 13 non-synonymous variants were detected. Functional studies revealed >50% reduced secretion of 7 variants, however, these variants were not significantly enriched in European CP patients. CONCLUSION: Our data indicate that CUZD1 might be a novel susceptibility gene for NACP. How these variants predispose to pancreatitis remains to be elucidated.

The reversion variant (p.Arg90Leu) at the evolutionarily adaptive p.Arg90 site in CELA3B predisposes to chronic pancreatitis.

A gain-of-function missense variant in the CELA3B gene, p.Arg90Cys (c.268C>T), has recently been reported to cause pancreatitis in an extended pedigree. Herein, we sequenced the CELA3B gene in 644 genetically unexplained French chronic pancreatitis (CP) patients (all unrelated) and 566 controls. No obvious loss-of-function variants were identified. None of the six low-frequency or common missense variants detected showed significant association with CP. Nor did the aggregate rare/very rare missense variants (n = 14) show any significant association with CP. However, p.Arg90Leu (c.269G>T), which was found in four patients but no controls, and affects the same amino acid as p.Arg90Cys, serves to revert p.Arg90 to the human elastase ancestral allele. As p.Arg90Leu has previously been shown to exert a similar functional effect to that of p.Arg90Cys, our findings not only confirm the involvement of CELA3B in the etiology of CP but also pinpoint a new evolutionarily adaptive site in the human genome.

Variants That Affect Function of Calcium Channel TRPV6 Are Associated With Early-Onset Chronic Pancreatitis.

BACKGROUND & AIMS: Changes in pancreatic calcium levels affect secretion and might be involved in development of chronic pancreatitis (CP). We investigated the association of CP with the transient receptor potential cation channel subfamily V member 6 gene (TRPV6), which encodes a Ca2+-selective ion channel, in an international cohort of patients and in mice. METHODS: We performed whole-exome DNA sequencing from a patient with idiopathic CP and from his parents, who did not have CP. We validated our findings by sequencing DNA from 300 patients with CP (not associated with alcohol consumption) and 1070 persons from the general population in Japan (control individuals). In replication studies, we sequenced DNA from patients with early-onset CP (20 years or younger) not associated with alcohol consumption from France (n = 470) and Germany (n = 410). We expressed TRPV6 variants in HEK293 cells and measured their activity using Ca2+ imaging assays. CP was induced by repeated injections of cerulein in TRPV6mut/mut mice. RESULTS: We identified the variants c.629C>T (p.A210V) and c.970G>A (p.D324N) in TRPV6 in the index patient. Variants that affected function of the TRPV6 product were found in 13 of 300 patients (4.3%) and 1 of 1070 control individuals (0.1%) from Japan (odds ratio [OR], 48.4; 95% confidence interval [CI], 6.3-371.7; P = 2.4 × 10-8). Twelve of 124 patients (9.7%) with early-onset CP had such variants. In the replication set from Europe, 18 patients with CP (2.0%) carried variants that affected the function of the TRPV6 product compared with 0 control individuals (P = 6.2 × 10-8). Variants that did not affect the function of the TRPV6 product (p.I223T and p.D324N) were overrepresented in Japanese patients vs control individuals (OR, 10.9; 95% CI, 4.5-25.9; P = 7.4 × 10-9 for p.I223T and P = .01 for p.D324N), whereas the p.L299Q was overrepresented in European patients vs control individuals (OR, 3.0; 95% CI, 1.9-4.8; P = 1.2 × 10-5). TRPV6mut/mut mice given cerulein developed more severe pancreatitis than control mice, as shown by increased levels of pancreatic enzymes, histologic alterations, and pancreatic fibrosis. CONCLUSIONS: We found that patients with early-onset CP not associated with alcohol consumption carry variants in TRPV6 that affect the function of its product, perhaps by altering Ca2+ balance in pancreatic cells. TRPV6 regulates Ca2+ homeostasis and pancreatic inflammation.

The three common polymorphisms p.A986S, p.R990G and p.Q1011E in the calcium sensing receptor (CASR) are not associated with chronic pancreatitis.

BACKGROUND: The calcium sensing receptor (CASR) is a G protein-coupled receptor that is responsible for assessing extracellular Ca2+ levels and thus plays a crucial role in calcium homeostasis. Hypercalcemia is a metabolic risk factor for pancreatitis and rare CASR variants have been described in patients with chronic pancreatitis. At the carboxy-terminal tail of CASR, there is a cluster of three common polymorphisms, p.A986S (rs1801725), p.R990G (rs1042636) and p.Q1011E (rs1801726), which have been associated with chronic pancreatitis in various studies, but with conflicting results. METHODS: We examined 542 German and 339 French patients with chronic pancreatitis as well as 1025 German controls for the 3 common CASR polymorphism by melting curve analysis. For comparison, we used genotype data from 583 French controls from a previous study. In addition, we functionally analyzed the three variants by NFAT and SRE luciferase reporter systems as well as Western blotting and verified cell surface expression by ELISA. RESULTS: In both cohorts, neither the genotype nor the allele frequencies differed significantly between patients and controls. In both luciferase assays, p.R990G showed a significant leftward shift, indicating an increased responsiveness of the receptor. p.A986S showed a leftward shift in the SRE but not in the NFAT reporter assay, while the responsiveness of p.Q1011E did not differ from the wild-type. These functional studies therefore do not support the contributions of variant CASR to increasing the risk of pancreatitis. CONCLUSIONS: The three frequent CASR polymorphisms are unlikely to increase the risk for chronic pancreatitis.

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.

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.

Pathogenic and likely pathogenic variants in at least five genes account for approximately 3% of mild isolated nonsyndromic thrombocytopenia.

BACKGROUND: Thrombocytopenia has a variety of different etiologies, both acquired and hereditary. Inherited thrombocytopenia may be associated with other symptoms (syndromic forms) or may be strictly isolated. To date, only about half of all the familial forms of thrombocytopenia have been accounted for in terms of well-defined genetic abnormalities. However, data are limited on the nature and frequency of the underlying causative genetic variants in individuals with mild isolated nonsyndromic thrombocytopenia. STUDY DESIGN AND METHODS: Thirteen known or candidate genes for isolated thrombocytopenia were included in a gene panel analysis in which targeted next-generation sequencing was performed on 448 French blood donors with mild isolated nonsyndromic thrombocytopenia. RESULTS: A total of 68 rare variants, including missense, splice site, frameshift, nonsense, and in-frame variants (all heterozygous) were identified in 11 of the 13 genes screened. Twenty-nine percent (N = 20) of the variants detected were absent from both the French Exome Project and gnomAD exome databases. Using stringent criteria and an unbiased approach, we classified seven predicted loss-of-function variants (three in ITGA2B and four in TUBB1) and four missense variants (one in GP1BA, two in ITGB3 and one in ACTN1) as being pathogenic or likely pathogenic. Altogether, they were found in 13 members (approx. 3%) of our studied cohort. CONCLUSION: We present the results of gene panel sequencing of known and candidate thrombocytopenia genes in mild isolated nonsyndromic thrombocytopenia. Pathogenic and likely pathogenic variants in five known thrombocytopenia genes were identified, accounting for approximately 3% of individuals with the condition.

Characterization of CEL-DUP2: Complete duplication of the carboxyl ester lipase gene is unlikely to influence risk of chronic pancreatitis.

BACKGROUND/OBJECTIVES: Carboxyl ester lipase is a pancreatic enzyme encoded by CEL, an extremely polymorphic human gene. Pathogenic variants of CEL either increases the risk for chronic pancreatitis (CP) or cause MODY8, a syndrome of pancreatic exocrine and endocrine dysfunction. Here, we aimed to characterize a novel duplication allele of CEL (CEL-DUP2) and to investigate whether it associates with CP or pancreatic cancer. METHODS: The structure of CEL-DUP2 was determined by a combination of Sanger sequencing, DNA fragment analysis, multiplex ligation-dependent probe amplification and whole-genome sequencing. We developed assays for screening of CEL-DUP2 and analyzed cohorts of idiopathic CP, alcoholic CP and pancreatic cancer. CEL protein expression was analyzed by immunohistochemistry. RESULTS: CEL-DUP2 consists of an extra copy of the complete CEL gene. The allele has probably arisen from non-allelic, homologous recombination involving the adjacent pseudogene of CEL. We found no association between CEL-DUP2 carrier frequency and CP in cohorts from France (cases/controls: 2.5%/2.4%; P = 1.0), China (10.3%/8.1%; P = 0.08) or Germany (1.6%/2.3%; P = 0.62). Similarly, no association with disease was observed in alcohol-induced pancreatitis (Germany: 3.2%/2.3%; P = 0.51) or pancreatic cancer (Norway; 2.5%/3.2%; P = 0.77). Notably, the carrier frequency of CEL-DUP2 was more than three-fold higher in Chinese compared with Europeans. CEL protein expression was similar in tissues from CEL-DUP2 carriers and controls. CONCLUSIONS: Our results support the contention that the number of CEL alleles does not influence the risk of pancreatic exocrine disease. Rather, the pathogenic CEL variants identified so far involve exon 11 sequence changes that substantially alter the protein's tail region.

Analysis of GPRC6A variants in different pancreatitis etiologies.

BACKGROUND: The G-protein-coupled receptor Class C Group 6 Member A (GPRC6A) is activated by multiple ligands and is important for the regulation of calcium homeostasis. Extracellular calcium is capable to increase NLRP3 inflammasome activity of the innate immune system and deletion of this proinflammatory pathway mitigated pancreatitis severity in vivo. As such this pathway and the GPRC6A receptor is a reasonable candidate gene for pancreatitis. Here we investigated the prevalence of sequence variants in the GPRC6A locus in different pancreatitis aetiologies. METHODS: We selected 6 tagging SNPs with the SNPinfo LD TAG SNP Selection tool and the functional relevant SNP rs6907580 for genotyping. Cohorts from Germany, further European countries and China with up to 1,124 patients and 1,999 controls were screened for single SNPs with melting curve analysis. RESULTS: We identified an association of rs1606365(G) with alcoholic chronic pancreatitis in a German (odds ratio (OR) 0.76, 95% confidence interval (CI) 0.65-0.89, p = 8 × 10-5) and a Chinese cohort (OR 0.78, 95% CI 0.64-0.96, p = 0.02). However, this association was not replicated in a combined cohort of European patients (OR 1.18, 95% CI 0.99-1.41, p = 0.07). Finally, no association was found with acute and non-alcoholic chronic pancreatitis. CONCLUSIONS: Our results support a potential role of calcium sensing receptors and inflammasome activation in alcoholic chronic pancreatitis development. As the functional consequence of the associated variant is unclear, further investigations might elucidate the relevant mechanisms.

The Experimentally Obtained Functional Impact Assessments of 5' Splice Site GT'GC Variants Differ Markedly from Those Predicted.

INTRODUCTION: 5' splice site GT>GC or +2T>C variants have been frequently reported to cause human genetic disease and are routinely scored as pathogenic splicing mutations. However, we have recently demonstrated that such variants in human disease genes may not invariably be pathogenic. Moreover, we found that no splicing prediction tools appear to be capable of reliably distinguishing those +2T>C variants that generate wild-type transcripts from those that do not. METHODOLOGY: Herein, we evaluated the performance of a novel deep learning-based tool, SpliceAI, in the context of three datasets of +2T>C variants, all of which had been characterized functionally in terms of their impact on pre-mRNA splicing. The first two datasets refer to our recently described "in vivo" dataset of 45 known disease-causing +2T>C variants and the "in vitro" dataset of 103 +2T>C substitutions subjected to full-length gene splicing assay. The third dataset comprised 12 BRCA1 +2T>C variants that were recently analyzed by saturation genome editing. RESULTS: Comparison of the SpliceAI-predicted and experimentally obtained functional impact assessments of these variants (and smaller datasets of +2T>A and +2T>G variants) revealed that although SpliceAI performed rather better than other prediction tools, it was still far from perfect. A key issue was that the impact of those +2T>C (and +2T>A) variants that generated wild-type transcripts represents a quantitative change that can vary from barely detectable to an almost full expression of wild-type transcripts, with wild-type transcripts often co-existing with aberrantly spliced transcripts. CONCLUSION: Our findings highlight the challenges that we still face in attempting to accurately identify splice-altering variants.