Article Text

Original article
CFTR, SPINK1, CTRC and PRSS1 variants in chronic pancreatitis: is the role of mutated CFTR overestimated?
  1. Jonas Rosendahl1,2,
  2. Olfert Landt3,
  3. Jana Bernadova1,
  4. Peter Kovacs4,
  5. Niels Teich5,
  6. Hans Bödeker6,
  7. Volker Keim7,
  8. Claudia Ruffert7,
  9. Joachim Mössner7,
  10. Andreas Kage8,
  11. Michael Stumvoll9,
  12. David Groneberg10,
  13. Renate Krüger4,
  14. Werner Luck1,
  15. Matthias Treiber11,
  16. Michael Becker1,
  17. Heiko Witt1,12
  1. 1Department for Paediatrics, Charité, Berlin, Germany
  2. 2Department for Gastroenterology, University of Leipzig, Leipzig, Germany
  3. 3TIB Molbiol, Berlin, Germany
  4. 4Interdisciplinary Centre for Clinical Research Leipzig, University of Leipzig, Leipzig, Germany
  5. 5Internistische Gemeinschaftspraxis für Verdauungs- und Stoffwechselerkrankungen, Leipzig, Germany
  6. 6Department for Internal Medicine, Hospital Freiberg, Freiberg, Germany
  7. 7Department for Gastroenterology, University of Leipzig, Leipzig, Germany
  8. 8Department for Laboratory Medicine, Charité, Berlin, Germany
  9. 9Department for Endocrinology, University of Leipzig, Leipzig, Germany
  10. 10Department of Occupational Medicine, University of Frankfurt, Frankfurt, Germany
  11. 11Department for Gastroenterology, Technische Universität München (TUM), Munich, Germany
  12. 12Department for Paediatrics & Else Kröner-Fresenius-Zentrum für Ernährungsmedizin (EKFZ), Technische Universität München (TUM), Munich, Germany
  1. Correspondence to Professor Heiko Witt, Department for Paediatrics & Else Kröner-Fresenius-Zentrum (EKFZ), Technical University Munich (TUM), Gregor-Mendel-Str. 2, 85350 Freising, Germany; heiko.witt{at}lrz.tu-muenchen.de

Abstract

Objective In chronic pancreatitis (CP), alterations in several genes have so far been described, but only small cohorts have been extensively investigated for all predisposing genes.

Design 660 patients with idiopathic or hereditary CP and up to 1758 controls were enrolled. PRSS1, SPINK1 and CTRC were analysed by DNA sequencing, and cystic fibrosis transmembrane conductance regulator (CFTR) by melting curve analysis.

Results Frequencies of CFTR variants p.R75Q, p.I148T, 5T-allele and p.E528E were comparable in patients and controls. We identified 103 CFTR variants, which represents a 2.7-fold risk increase (p<0.0001). Severe cystic fibrosis (CF)-causing variants increased the risk of developing CP 2.9-fold, and mild CF-causing variants 4.5-fold (p<0.0001 for both). Combined CF-causing variants increased CP risk 3.4-fold (p<0.0001), while non-CF-causing variants displayed a 1.5-fold over-representation in patients (p=0.14). CFTR compound heterozygous status with variant classes CF-causing severe and mild represented an OR of 16.1 (p<0.0001). Notably, only 9/660 (1.4%) patients were compound heterozygotes in this category. Trans-heterozygosity increased CP risk, with an OR of 38.7, with 43/660 (6.5%) patients and 3/1667 (0.2%) controls being trans-heterozygous (p<0.0001).

Conclusions Accumulation of CFTR variants in CP is less pronounced than reported previously, with ORs between 2.7 and 4.5. Only CF-causing variants reached statistical significance. Compound and trans-heterozygosity is an overt risk factor for the development of CP, but the number of CFTR compound heterozygotes in particular is rather low. In summary, the study demonstrates the complexity of genetic interactions in CP and a minor influence of CFTR alterations in CP development.

  • Chronic pancreatitis
  • genetics
  • CFTR
  • SPINK1
  • CTRC
  • PRSS1
  • genotype
  • general practice
  • gene mutation
  • genetic testing
  • pancreatic disorders
  • pancreatic enzymes
  • Crohn's disease
  • 6-mercaptopurine
  • IBD
  • azathioprine
  • IBD clinical
  • chronic ulcerative colitis
  • pancreas
  • pancreatitis
  • genetic polymorphisms
  • genetic testing

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Significance of this study

What is already known on this subject?

  • In addition to variants of PRSS1, SPINK1 and CTRC, cystic fibrosis transmembrane conductance regulator (CFTR) variants contribute to the pathogenesis of chronic pancreatitis (CP).

  • Patient and control cohorts in studies that have screened for variants in all known genes were small, and CFTR variants in most cases have not been analysed comprehensively.

  • Data obtained for cystic fibrosis (CF) indicate that predominantly mild CFTR variants with residual CFTR function can be expected in CP.

What are the new findings?

  • CF-causing CFTR variants increase the risk of CP development (CF-causing severe OR 2.9; CF-causing mild OR 4.5), while non-CF-causing variants are more common in patients (1.5-fold), but this does not reach statistical significance; the overall risk contribution of CFTR variants to the pathogenesis of CP seems to be lower than previously reported (OR 2.7).

  • All CFTR compound heterozygotes carry at least one mild CFTR allele, and compound heterozygotes with two CF-causing variants are rare (1.4%), which shows that atypical CF is present in a small minority of CP patients only.

  • Complex genotypes such as homozygous, compound heterozygous and trans-heterozygous carriers are significantly more common in patients with ORs for distinct complex genotypes of 16.1–135.4, indicating that CP in many cases is a complex inherited disease.

How might it impact on clinical practice in the foreseeable future?

  • In patients with early-onset CP and no identifiable aetiological factor, genetic testing for distinct variants (PRSS1, SPINK1, CTRC and CFTR) should be offered after comprehensive education.

  • It seems that, for the pathogenesis of CP, the influence of CFTR variants has so far been overestimated, but a diagnosis of CF has to be ruled out, especially in patients with early-onset CP.

Introduction

Chronic pancreatitis (CP) is a relapsing or continuing inflammatory disease characterised by recurrent or persisting abdominal pain leading to permanent impairment of pancreatic function in some patients.1 ,2 A genetic background for CP was first described in a pedigree with an obviously autosomal dominant inheritance pattern in 1952.3 Depiction of the underlying defects in different genes established a model in which premature trypsin activation plays a central pathogenic role.

Genetic variations in cationic trypsinogen (PRSS1) (OMIM *276000) were predominantly found in families with a positive history of CP.4 Functional analyses indicate that gain-of-function—that is, premature trypsinogen activation—might be the underlying mechanism of mutated PRSS1.5 This concept of gain-of-function is further supported by a description of duplication and triplication of the trypsinogen locus in some CP families originating from France.6

In addition, alterations in the serine protease inhibitor, Kazal type 1 (SPINK1) (OMIM *167790), a potent trypsin inhibitor, have been described in patients with different forms of CP.7–10 In contrast with PRSS1, SPINK1 variants seem to induce a loss of SPINK1 inhibitor function. Even though there is a strong genetic association of the most commonly found SPINK1 alteration, p.N34S, with CP, functional data revealed no reduced inhibitory capacity of this variant.11 Following the pathophysiological concept, a loss-of-function variant (p.G191R) in the anionic trypsinogen (PRSS2) (OMIM *601564), the second major trypsin isoform was detected more commonly in controls.12 Variations in chymotrypsinogen C (CTRC) (OMIM *601405), an enzyme that promotes degradation of trypsin, have been associated with idiopathic and alcoholic CP. On a functional level, CTRC mutations exhibited diminished secretion and/or activity.13

Cystic fibrosis (CF) (OMIM 219700) is an autosomal recessive inherited disorder characterised by chronic endobronchial infection, pancreatic insufficiency and male infertility. Approximately 1–2% of patients with CF have recurrent attacks of pancreatitis. Similarly to CF, some patients with CP show raised sweat chloride levels. Moreover, in both disorders, obstruction of the pancreatic duct due to dense secretion can be found, suggesting a possible link between these different disease entities.14–18 Subsequently, variants of the cystic fibrosis transmembrane conductance regulator (CFTR) (OMIM *602421) were found to be associated with idiopathic and alcoholic CP.19 ,20

The complexity of the role of CFTR variants in CP and other diseases is disclosed by the large number of about 1900 variants that have so far been described (http://www.genet.sickkids.on.ca/cftr). Analysis of CFTR is further complicated by its size, comprising 27 exons. Moreover, the functional effects of CFTR variants are poorly understood in most cases.

Analysis of chloride channel properties has revealed abnormal sweat chloride levels and nasal potential differences in patients with CP compared with healthy controls when at least one CFTR variant was reported. The degree of impairment correlated with the number and type of CFTR variants detected, and was less pronounced than in pancreatic-sufficient or -insufficient CF patients.21 In CF, genotype–phenotype correlation has been described for exocrine pancreatic function. Patients with exocrine insufficiency in most cases carry two ‘severe’ CFTR variants, while patients with sufficient exocrine function bear at least one ‘mild’ CFTR variant. Pancreatic insufficiency prevalence scores in CF patients with pancreatitis were most likely ‘mild’, in contrast with patients without pancreatitis who displayed ‘moderate–severe’ scores.22 A possible conclusion is that, for the development of CP, residual CFTR function and pancreatic acinar cells are essential, which might be missing in CF patients with ‘moderate–severe’ pancreatic insufficiency prevalence scores and/or two ‘severe’ CFTR mutations.

These facts raise the question whether predominantly ‘mild’ CFTR variants can be found in CP and whether CFTR variants, like variants of other known pancreatitis-causing genes, are more common in different age groups. Owing to the limitations mentioned above, the number of patients and controls, if investigated at all, were rather small in previous studies, and a complete analysis of CFTR was performed in small cohorts only. To investigate the complex genetic basis of CP, we examined CFTR as well as PRSS1, SPINK1 and CTRC in a large cohort of patients with CP and in controls.

Patients and methods

Study population

The study protocol was approved by the medical ethics review committees of the University of Berlin (Charité) and the University of Leipzig. All participants gave informed consent. If patients were under the age of 18 years, parents gave informed consent.

Diagnosis of CP was based on two or more of the following findings as published previously13: briefly, a history of recurrent pancreatitis or recurrent abdominal pain typical of CP, pancreatic calcifications and/or pancreatic ductal irregularities revealed by endoscopic retrograde pancreaticography or by MRI of the pancreas, and/or pathological sonographic findings. Hereditary CP (HP) was diagnosed when one first-degree relative or two or more second-degree relatives had recurrent acute or chronic pancreatitis without any apparent precipitating factor. Idiopathic CP (ICP) was diagnosed in the absence of a positive family history or possible precipitating factors, such as alcohol abuse, trauma, medication, infection and metabolic disorders.

In total, 660 unrelated patients with HP or ICP from Berlin and Leipzig were enrolled (347 female; 107 HP, 553 ICP; median age 16 years, mean age 20.1 years, age range 1–81 years). The study population from Berlin comprised 439 patients (244 female; 101 HP, 338 ICP; median age 13 years, mean age 17.1 years, age range 1–70 years). A total of 221 German patients were recruited from Leipzig (118 female; six HP, 215 ICP; median age 23 years, mean age 26 years, age range 5–81 years).

Healthy unrelated German controls recruited for genetic association studies (medical students/staff, healthy retirees and healthy blood donors; n=1758) originated from Berlin and Leipzig (961 female, 104 of unknown gender; median age 32 years, mean age 38.7 years, age range 17–101 years). The control group from Berlin included 1644 individuals (927 female, 52 of unknown gender; median age 31 years, mean age 38.5 years, age range 17–101 years), whereas the Leipzig controls consisted of 114 individuals (34 female, 52 of unknown gender; median age 44.5 years, mean age 43.6 years, age range 20–64 years).

Methods

PCR

Genomic DNA was extracted from peripheral blood leucocytes. PCR was performed using 0.75 U AmpliTaq Gold polymerase (Applied Biosystems, Darmstadt), 400 μM dNTPs, 1.5 mM MgCl2 and 0.1 μM each primer in a total volume of 25 μl. Cycle conditions were an initial denaturation for 12 min at 95°C followed by 48 cycles of 20 s denaturation at 95°C, 40 s annealing at specific temperatures, 90 s primer extension at 72°C, and a final extension for 2 min at 72°C in an automated thermal cycler.

Primers were synthesised according to the published nucleotide sequences (PRSS1: GenBank: NM_002769.3; SPINK1: GenBank: NM_003122.3; CTRC: GenBank: NM_007272.2; CFTR: GenBank: NM_000492.3). Primers for the analysis of PRSS1, SPINK1 and CTRC were used as published previously.7 ,13 ,23 Oligonucleotide sequences and annealing temperatures are available upon request.

DNA sequencing of SPINK1, CTRC and PRSS1

PCR products were digested with shrimp alkaline phosphatase (USB) and exonuclease I (GE Healthcare, Munich). Cycle sequencing was performed using BigDye terminator mix (Applied Biosystems). Reaction products were purified by ethanol precipitation or on a Sephadex G-50 column (GE Healthcare, Munich) and loaded on to an ABI 3100-Avant or an ABI 3730 fluorescence sequencer (Applied Biosystems). DNA mutation numbering is based on cDNA sequences that use the A of the ATG start codon as nucleotide +1. Mutations are described according to the nomenclature recommended by the Human Genome Variation Society (http://www.hgvs.org/mutnomen).

Melting curve analysis of CFTR

Melting curve analysis for 37 CFTR variants was performed in the LightCycler instrument (Roche Diagnostics, Mannheim) using a pair of fluorescent resonance energy transfer (FRET) probes. The following CFTR variants were analysed with specific FRET probes: p.E60X, p.R75Q, p.G85E, p.R117H, p.I148T, c.621+1G>T (IVS4+1G>T), c.711+1G>T (IVS5+1G>T), c.1078delT, p.R334W, p.R347P, 9-13TG, 5-9T, p.A455E, p.M470V, p.F508del, c.1716G>A (p.E528E), c.1717-1G>A (IVS10-1G>A), p.G542X, p.S549N, p.R553X, p.R560T, c.1898+1G>A (IVS12+1G>A), c.2143delT, c.2183AA>G, c.2562T>G, c.2657+5G>A (IVS14B+5G>A), p.L997F, p.I1005R, p.Y1092X, p.D1152H, p.R1162X, c.3659delC, p.S1235R, p.S1251N, p.W1282X, p.N1303K, and c.4389G>A. Since FRET probes were designed in nearly all cases complementary to the mutant, the probes were able to detect more than 100 additional variants located adjacent to the specific target nucleotide (unpublished data). All probes were devised and synthesised by TIB Molbiol (Berlin, Germany). In complex target regions, locked nucleic acids were integrated for better discrimination. Sequences of the probes are available upon request. If melting curves differed from the expected curves for the wild-type or the specific variant, DNA sequencing was performed to identify the underlying alteration.

Statistical analysis

The significance of the differences between variant frequencies in affected individuals and controls was tested by two-tailed Fisher's exact test and was calculated using GraphPad Prism (v4.03). p<0.05 was considered to be significant. We excluded CFTR variants p.R75Q, p.I148T, the poly-T tract and p.E528E, PRSS1 variant p.S124F, SPINK1 variants c.1-52G>T and p.P55S, and CTRC variants p.R37Q, p.K151N and p.K172E from all calculations, because of missing functional data or a missing genetic association.

Results

Genetic analysis of SPINK1, CTRC and PRSS1

We found PRSS1 variants in 8.3% of the patients (55/660), but in none of 1758 controls (p<0.0001) (table 1). The following single variants were significantly more common in patients: p.A16V, p.N29I, p.R122C and p.R122H with ORs ranging from 29.5 to 141.1 (p=0.002 for p.R122C; for all other variants p<0.0001).

Table 1

Distribution of SPINK1, CTRC and PRSS1 variants in patients and controls

In SPINK1, 12 different variants were detected. In the patient group, c.(1-215G>A;194+2T>C) (p<0.0001, OR 38.1) and p.N34S (homozygous/heterozygous: p<0.0001, OR 95.6 and 10.5) were significantly accumulated. In summary, after exclusion of c.1-52G>T and p.P55S, 19.6% of the patients (129/660) carried a SPINK1 variant compared with 1.5% of the controls (27/1758) (p<0.0001) accounting for an OR of 15.6 (95% CI 10.2 to 23.9).

A CTRC variant was determined in 4% of patients (22/546) and in 0.8% of controls (13/1667). This represents a significant accumulation in patients (p<0.0001; OR 5.3, 95% CI 2.7 to 10.7). Two single variants, p.R254W and p.K247_R254del, were more common in patients (p=0.0005, OR 4; p=0.003, OR 28.2, respectively). To obtain stringent results, calculations were based on the smallest number of controls that were analysed for CTRC exons 2, 3 and 7 (table 1).

Genetic analysis of CFTR

Melting curve analysis detected 103 variants and characterised frequencies of five common alterations (c.1210-34TG(9_13), c.1210-12T(5_9), p.M470V, c.2562T>G and c.4389G>A). In total, CFTR variants were found in 103/660 (15.6%) patients compared with 112/1758 (6.4%) controls (OR 2.7, 95% CI 2 to 3.6). In subgroups, OR was 2.9 (95% CI 2 to 4.3) for CF-causing severe variants, 4.5 (95% CI 2.3 to 8.8) for CF-causing mild variants, and 3.4 (95% CI 2.4 to 4.8) for all CF-causing variants, whereas, for all non-CF-causing variants, there was a 1.5-fold over-representation in patients only (p=0.14) (table 2). Notably, p.I148T was never found in cis with c.3199del6, and p.S1235R was always associated with p.G628 (wild-type).

Table 2

CFTR variants detected by melting curve analysis

Three common variants, c.1210-34TG(9_13), p.M470V and c.4389G>A, were differently distributed in patients and controls (table 3). TG10-allele of c.1210-34TG(9_13) was found more commonly in patients (435/1320, 33%) than in controls (983/3516, 28%) (p=0.0008, OR 1.3), whereas allele TG11 was found more often in controls (2137/3516, 60.8%) than in patients (750/1320, 56.8%) (p=0.014, OR 0.9). With regard to p.M470V, the 470M allele accumulated in patients (607/1320, 46%) in comparison with controls (1480/3516, 42.1%) (p=0.02, OR 1.2). Analysis of c.4389G>A displayed an over-representation of the 4389A allele in patients (366/1320, 27.7%; controls: 544/2330, 23.3%; p=0.004, OR 0.8). Genotype data of these variants confirmed results for TG10/10, p.470MM and c.4389AA (p<0.0001, p<0.009 and p<0.0005, respectively).

Table 3

Allele frequencies of CFTR variants c.1210-34TG(9_13), c.1210-12T(5_9), p.M470V, c.2562T>G and c.4389G>A

To confirm that association of these variants is not due to co-segregation with p.F508del, which is found in 6.7% of patients and in 2.7% of control subjects, we calculated data for the variants after exclusion of p.F508del carriers in both patients and controls. Association of c.1210-34TG(9_13) (TG10-allele and TG10/10 genotype) and c.4389G>A (A-allele and c.4389AA genotype) persisted after correction (p=0.001, p=0.002; p=0.0005, p=0.0003). In contrast, association of p.M470 did not withstand correction with regard to the genotype (p=0.1), while the 470M allele was still over-represented in patients (p=0.01) (data not shown). The allele frequency for p.M470V did not withstand Bonferroni correction.

Homozygous and compound heterozygous carriers of CFTR, SPINK1 and CTRC variants

The SPINK1 combination p.N34S/p.N34S was found in 17/660 (2.6%) and p.N34S/c.(1-215G>A;194+2T>C) in 7/660 patients (1.1%) (total 24/660, 3.6%), but never in controls (p<0.0001; OR 135.4, 95% CI 8.2 to 2231). One patient was compound heterozygous for CTRC p.V235I/p.R254W and one homozygous for p.R254W (1/546; 0.2% each). None of the 1700 controls was homozygous or compound heterozygous for CTRC (p=0.2).

Compound heterozygous CFTR carriers with severe and mild CF-causing variants were found in 9/660 (1.4%) patients and in 1/1758 (0.06%) controls (p=0.002; OR 16.1, 95% CI 1.9 to 134.2). Combinations of CF-causing severe or CF-causing mild or non-CF-causing with non-CF-causing variants were present in 17/660 (2.6%) patients and in 14/1758 (0.8%) controls (p=0.001, OR 3.3, 95% CI 1.6 to 6.7). Except one (p.R117H (7T/7T)/p.S1235), these compound heterozygotes were excluded in the overall computations, because CFTR variant p.R75Q, p.I148T, 5T or p.E528E was present in at least one allele. In total, 10/660 (1.5%) patients and 1/1758 (0.06%) controls were compound heterozygous CFTR carriers (p<0.0001, OR 27, 95% CI 3.5 to 211.7). Individually, none of the combinations of two CFTR variants was significantly over-represented in patients or controls (table 4).

Table 4

Homozygous and compound heterozygous patients and controls with at least two CFTR, SPINK1 or CTRC variants

Trans-heterozygous carriers of CFTR, SPINK1, CTRC and PRSS1 variants

Trans-heterozygotes were found in 43/660 (6.5%) patients and 3/1667 (0.2%) controls (p<0.0001, OR 38.7, 95% CI 12 to 125.1). After exclusion of carriers in which a non-CF-causing variant was essential for trans-heterozygous status, 33/660 (5%) patients and 2/1667 (0.1%) controls were trans-heterozygous (p<0.0001, OR 43.8, 95% CI 10.5 to 183.2). In addition, 24/660 (3.6%) patients and 2/1700 (0.1%) controls were trans-heterozygotes in that CFTR variants p.R75Q, p.I148T, 5T and p.E528E were essential for trans-heterozygous status (p<0.0001, OR 30.8, 95% CI 7.3 to 131). These samples were not integrated into our computations, because the over-representation of these combinations is explained by the accumulation of the concomitant variant (eg, p.N34S, see below) and the mentioned CFTR variants were found with similar frequencies in patients and controls. It is noteworthy that nine of the patients carried a severe PRSS1 variant (p.N29I, p.R122C and p.R122H). As such in these cases the CFTR variants should be of minor or no importance for the development of CP. Calculations in patients were referred to the total number of 660 patients, although 114 patients were not analysed for CTRC variants, and in controls to 1667 individuals who were investigated completely for variants in all four genes to obtain more rigorous results (table 5).

Table 5

Trans-heterozygous carriers of CFTR, SPINK1, CTRC and PRSS1 variants

Individually, we found a significant strengthening for the following combinations in patients: p.A16V(PRSS1)/p.N34S(SPINK1), p.R122H(PRSS1)/p.E528E(CFTR), p.N34S(SPINK1)/p.R75Q(CFTR), p.N34S(SPINK1)/5T(CFTR), p.N34S(SPINK1)/p.F508del(CFTR) and p.N34S(SPINK1)/p.S1235R(CFTR) (p=0.007, p=0.03, p=0.03, p=0.002, p<0.0001 and p=0.03; range of OR 10.7–45.8). In these cases, the over-representation is also most probably due to the strengthening of p.R122H (PRSS1) and p.N34S (SPINK1) as shown for CFTR p.R75Q (see below). One patient (p.N34S/p.R254W/p.R117H(7T/7T); 1/546, 0.2%) carried variants in three genes (p value 0.6).

SPINK1 variant p.N34S (heterozygous and homozygous) was found in combination with CFTR variants p.R75Q (heterozygous and homozygous) in 6/660 patients (0.9%), with 5T in 8/660 patients (1.2%), and with p.E528E (heterozygous and homozygous) in 6/660 patients (0.9%) compared with 1/1758 (0.06%) controls for combination with p.R75Q and 5T, and 2/1758 controls (0.1%) for combination with p.E528E (p.R75Q: p=0.002, OR 16.1, 95% CI 1.92 to 134; 5T: p=0.0002, OR 21.6, 95% CI 2.7 to 172.8; p.E528E: p=0.008, OR 8, 95% CI 1.6 to 40). Since this might be explained by the over-representation of p.N34S in patients, we calculated frequencies of p.N34S carriers in patients (n=109) and controls (n=26). Thereafter, no significant association was obtained (p.R75Q: 6/109 patients, 5.5%, 1/26 controls, 3.9%, p=1.0; 5T: 8/109 patients, 7.3%, 1/26 controls, 3.9%, p=1.0; p.E528E: 6/109 patients, 5.5%, 2/26 controls, 7.7%, p=0.7). In contrast, 6/31 (19.4%) patients with p.R75Q also carried p.N34S, whereas only 1/60 (1.7%) of controls with p.R75Q showed the p.N34S variant, representing the data expected in patients and controls. On analysis of about 200 parents for CFTR, p.R75Q was found to pass from one parent to the child in 6/13 cases (46%) only, which further underlines that this variant most likely does not contribute to disease pathogenesis.

Table 6

Distribution of PRSS1, SPINK1, CTRC and CFTR variants in different age groups (age >20 or <20 years)

Variant distribution in patients with HP and ICP

In the online supplementary tables 1–5, we have additionally divided the patient group into patients with HP and ICP and calculated p values and ORs for these groups. As expected, PRSS1 variants were predominantly found in HP patients (30/107, 28%) compared with ICP patients (25/553, 4.5%). Frequencies of SPINK1, CTRC and CFTR variants did not differ significantly between HP and ICP. CF-causing severe variants were over-represented in HP patients (6/107, 5.6%), but this did not reach statistical significance in contrast with ICP patients (49/553, 8.9%). Compound heterozygous CFTR carriers were more common among HP (4/107, 3.7%) than ICP (5/553, 0.9%) patients.

Variant distribution in patients aged >20 and <20 years

In younger patients, overall PRSS1 variants were 2.9-fold more common (>20 years: 9/239, 3.8%; <20 years: 46/421, 10.9%; p=0.001, OR 3.1, 95% CI 1.5 to 6.5), whereas overall SPINK1 variants were similarly distributed (56/239, 23.4%; 73/421, 17.3%; p=0.07). CTRC variants were also found in similar frequencies in both age groups (>20 years: 6/170, 3.5%; <20 years: 16/376, 4.3%; p=0.8). All CFTR variants combined were more common in younger patients (74/421, 17.6%) than in older patients (29/239, 12.1%), but this difference was not significant (p=0.07).

Complex genotypes, which includes homozygous carriers of p.N34S and p.R254W and compound heterozygous and trans-heterozygous carriers, were found in similar distribution in younger (55/421, 13.1%) and older (24/239 (10%) patients (p=0.3), when summarised together (table 6).

Discussion

This is the largest study to date that presents genetic data for CFTR, SPINK1, CTRC and PRSS1 in patients with ICP and HP. Although alcohol abuse is the predominant cause of CP, we decided to exclude patients with alcoholic CP because ICP/HP and alcoholic CP seem only share the same genetic risks to some extent.9 ,13 ,25 Therefore, for alcoholic CP with expected modest genetic risk factors, a different approach such as genome-wide association studies is warranted. Extensive screening for all known CP predisposing genes allows us to display the complex genetic interactions in ICP and HP. Association of CFTR variants with CP is well established, but classification of some variants with regard to their influence on the development of CP and the aligned risk increase for CFTR carriers had not unambiguously been defined. In a recent consensus statement, the clinical consequences of some CFTR variants and their association with exocrine pancreatic function have been classified.24 Accordingly, several common variants should be excluded from the computation, which will provide a more conservative and realistic calculation.

PRSS1, SPINK1 and CTRC variants were over-represented in patients, as reported previously. PRSS1 variants were predominantly found in younger patients (<20 years) (p=0.001), whereas SPINK1 variants did not differ in the two age groups (p=0.07). In line with former studies, we detected the SPINK1 p.P55S variant in similar frequencies in patients and controls.7 ,8 ,26–28 Moreover, functional studies failed to show any effect of p.P55S on SPINK1 inhibitory capacity.29 Consequently, we excluded p.P55S from all calculations.

In contrast with other authors, we excluded CFTR variants p.R75Q, p.I148T, 5T and p.E528E from our calculations, because they were distributed similarly in patients and controls.30 ,31 Thereby, overall, CFTR variants displayed a 2.7-fold risk increase for CP development. We found the strongest association with CF-causing mild CFTR variants (OR 4.5); however, these variants were found in only 3.5% of patients, while CF-causing severe variants were found in 8.3% of our patients and increased CP risk 2.9-fold. The importance of non-CF-causing variants is questionable. Although an association has been described in former studies, the overall effect of these variants is expected to be rather low, and, in our study, these variants did not reach significance (p=0.14). In total, CFTR variants displayed a 3.1-fold risk increase (p<0.0001) in young patients under 20 years of age, whereas in older patients the risk increase was only twofold (p=0.003) (comparison of both groups: p=0.07).

Some common variants such as the poly-T tract in intron 9 influence CFTR transcript levels.32 In addition, common CFTR haplotypes seem to modulate susceptibility to CP.33 A large number of our patients were part of a study that identified a CFTR haplotype (TG10-7T-M470) that increases the risk of CP. Thus, it is not surprising that our analysis of these common variants is in line with the former finding. TG10 and c.4389A accumulated in CP patients even after exclusion of p.F508del carriers, highlighting the potential importance of common variants in CP.

Compound heterozygotes for CFTR and SPINK1 variants were found more commonly in patients than in controls. Overall, CFTR compound heterozygotes displayed an OR of 27 (95% CI 3.5 to 211.7), but this accounts for 10 patients (1.5%) only. In contrast, 24/660 (3.6%) patients were homozygous or compound heterozygous for a SPINK1 variant (p<0.0001; OR 135.4, 95% CI 8.2 to 2231), but only 2/546 (0.4%) for a CTRC variant (p=0.2). In this regard, it is justifiable to compare ORs of PRSS1 variant carriers (OR 322.4), homozygous and compound heterozygous SPINK1 carriers (OR 135.4) and compound heterozygous CFTR carriers (OR 27), while overall carriers of CTRC variants (OR 5.3) can be compared with overall SPINK1 (OR 15.7) and CFTR carriers (OR: CF-causing severe, 2.9; CF-causing mild, 4.5; all CFTR variants, 2.7) to classify their influence on the pathogenesis of CP. These data highlight the minor importance of CFTR variants in CP development.

Recent publications partly presented a stronger accumulation of CFTR variants and compound heterozygotes, but comprehensive screening of a circumscribed German group of patients and controls revealed comparable results.34 Different results are attributable to smaller sample sizes analysed, a lack of screening of controls, and inclusion of common but innocuous variants, whereas we chose a very stringent and, in our eyes, reliable calculation by excluding these variants.

The complexity of CP pathogenesis is demonstrated when interactions of all so far known pancreatitis genes are investigated. Trans-heterozygotes accumulated in patients, and 6.5% of our patients carried variants in at least two genes (OR 38.7, 95% CI 12 to 125.1). We could not replicate the data of a former study that reported a 900-fold increase in risk for patients with two CFTR variants and SPINK1 variant p.N34S, which is probably because we analysed control subjects, unlike the former publication.30 Nevertheless, trans-heterozygotes in particular demonstrate that carriage of different variants in different pancreatitis-associated genes increases the CP risk substantially and may in some individuals explain why the disease developed.

However, there is some dispute about how to interpret complex genetic results. In a recent study, trans-heterozygosity for SPINK1 p.N34S with CFTR p.R75Q was reported to increase CP risk.31 We also demonstrate a significant accumulation of p.N34S/p.R75Q trans-heterozygotes in patients (6/660, 0.9% vs controls 1/1758, 0.06%), and the data seem to portend an association for the combination of p.N34S with CFTR 5T-allele and p.E528E also. Of note, the calculation of trans-heterozygotes in the overall cohorts, which are all patients in comparison with all controls, may be incorrect. From our point of view, the data have to be compared between p.N34S carriers in patients and p.N34S carriers in controls. When we performed the calculation in these analogous groups, the described association cannot be sustained. Frequencies of p.R75Q in patients and controls carrying p.N34S were 6/109 (5.5%) and 1/26 (3.9%), respectively (p=1.0). This is also true for the combinations of p.N34S with 5T or p.E528E and supports the notion that the observed unequal distribution of trans-heterozygotes in the overall cohorts is due to the over-representation of p.N34S in the patient group.

In summary, we provide further evidence that CP is a complex genetic disorder characterised by different genetic alterations in different genes and a complex interplay of these alterations. CFTR variants found in CP patients are in most cases ‘mild’ variants with residual CFTR function, and in all CFTR compound heterozygotes at least one ‘mild’ CFTR allele was present. Frequencies of CFTR variants and CFTR compound heterozygotes in CP are lower than previously reported. With regard to the associated genes in CP, it seems conclusive to propose two terminal positions, which are CF on the one side and hereditary CP on the other (figure 1). Both forms will be characterised by ‘severe’ genetic changes sufficient for disease development in nearly all cases, while the influence of environmental factors may be minor. Pancreatitis will be rare in patients with two ‘severe’ CF-causing variants, and the numbers of patients with pancreatitis will increase when at least one ‘mild’ CFTR variant is present. In between, complex inherited CP including CFTR-related disorders and pancreatitis seems to be distributed, and there may be a floating transition to patients that carry variants in different pancreatitis-associated genes.

Figure 1

Chronic pancreatitis (CP) is a complex genetic disorder that develops through the interaction of environmental and genetic factors. Cystic fibrosis and hereditary CP seem to require no, or little, environmental influence apart from their disease-causing genetic basis. Patients with CP compound-heterozygous for CFTR variants may in many cases be designated as CFTR-related disorders. CFTR function decreases, which is displayed in positive sweat chloride measurements, when CFTR variants accumulate in patients, whereas exocrine pancreatic insufficiency increases. CF, CFTR variant; PRSS1, cationic trypsinogen; SPINK1, serine protease inhibitor, Kazal type 1; CTRC, chymotrypsinogen C; CFsev, CF-causing severe variant; CFmild, CF-causing mild variant; CFPoly, non-CF-causing variant.

Acknowledgments

We thank the individuals who participated in this study. We also thank Claudia Güldner, who has since died, and Markus Braun (both Berlin), Susanne Kistner (Leipzig), Knut Krohn, Kathleen Schön and Birgit Oelzner (IZKF core unit DNA technologies, Leipzig) for excellent technical assistance. We dedicate this article to the valued people who unfortunately died during this project.

References

View Abstract

Supplementary materials

  • Supplementary Data

    This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.

    Files in this Data Supplement:

Footnotes

  • Funding This work was supported by the Deutsche Forschungsgemeinschaft WI 2036/2-1, WI 2036/2-2 and WI 2036/2-3 (to HW) and RO 3929/1-1 and RO 3939/2-1 (to JR), by the Sonnenfeld-Stiftung (to HW and to RN), and by a grant from the Colora Stiftung GmbH (to JR).

  • Competing interests None.

  • Patient consent Obtained.

  • Ethics approval Ethics approval was provided by the University of Berlin and the University of Leipzig.

  • Provenance and peer review Not commissioned; externally peer reviewed.