Abstract
Aims:
Wilson disease (WD) is an autosomal recessive disorder of copper metabolism, characterized by its accumulation in tissues which results in hepatic, neurological, and/or psychiatric symptoms. The aim of this study was to investigate the genetics of WD in Croatian patients.
Methods:
Correlation of the clinical presentation subtype and the age at onset of the diagnosis of WD with the ATP7B genotype was investigated in a group of Croatian WD patients. DNA from peripheral blood samples was tested for the p.His1069Gln by direct mutational analysis and other polymorphisms were identified by sequence analysis of coding and flanking intronic regions of ATP7B gene.
Results:
In the group of 75 WD patients of Croatian origin, 18 different mutations in ATP7B gene were detected, three of which were novel. The p.His1069Gln mutation was most frequent, being detected in 44 Croatian WD patients (58.7%). Most ATP7B mutations (90.4%) were located in exons 5, 8, 13, 14, and 15.
Conclusions:
Clinical diagnosis of WD was confirmed in 59 patients by detecting mutations on both ATP7B alleles. The age at onset of WD and the type of WD clinical presentation showed no significant correlation with the ATP7B genotype.
Introduction
W
The ATP7B gene, coding the copper transporting P-type ATPase (Bull et al., 1993), has been reported to have 380 variants that were proven to affect the cellular pathways of copper metabolism, causing WD (EASL, 2012). The most frequent mutation found in populations of European origin is the p.His1069Gln mutation, located in the ATP loop of the ATP7B protein (Caca et al., 2001). The low incidence of WD, in addition to its genetic heterogeneity, makes the molecular characterization and comparison between genotype and phenotype difficult.
According to a meta-analysis conducted by Stapelbroek et al. (2004), the homozygous p.His1069Gln mutation genotype gives rise to neurological symptoms and later onset of the disease more frequently than both the heterozygous and the non-p.His1069Gln WD genotypes. Gromadzka et al. (2006) presented similar results although Vrabelova et al. (2005) did not find such correlations. It is hypothesized that missense mutations cause later onset of disease and less severe clinical presentation and disturbances of copper metabolism compared to mutations which have a profound effect on ATP7B protein product (Merle et al., 2010; Okada et al., 2010; Lee et al., 2011).
As no data about the genetic basis of WD in the Croatian population have been published thus far, the aim of this study is to investigate mutations in the ATP7B gene in Croatian WD patients, as well as to determine possible correlations of the type of clinical presentation and the age at onset of WD with the ATP7B genotype.
Materials and Methods
Subjects
DNA samples were collected from a Croatian cohort of 75 unrelated WD patients (40 males, 35 females). The diagnosis of WD was based on hepatic, neurological, and/or psychiatric symptoms, the presence of Kayser-Fleischer rings, and copper test results such as high hepatic copper concentration (>250 μg/g dry weight), low serum “free” copper, serum ceruloplasmin levels (<0.17 g/L), and increased 24-h urinary copper excretion (>1.7 μmol/24 h). WD presentation was classified according to Ferenci et al. (2003). Solely hepatic symptoms were present in 46 patients, including 43 patients with chronic liver symptoms (H2 type of WD) and 3 with acute liver symptoms (H1 type). Solely neurological presentation of WD (N2) was found in 14 patients, while 15 patients presented with both hepatic and neurological symptoms (N1) (Table 1).
KF, Kayser-Fleischer; CP, ceruloplasmin.
The median age (interquartile range [IQR]) at onset of WD symptoms in this cohort of patients was 20 years (IQR 16-29 years), median age (IQR) at diagnosis was 24 years (18-32 years), and the median difference between these two variables was 1 year (IQR 1-5 years). Patients were clinically characterized at the Department of Internal Medicine, Department of Neurology, and Department of Pediatrics, University Hospital Centre Zagreb. This study was approved by the Ethics Committee of the University Hospital Center Zagreb and an informed consent was obtained from all patients.
ATP7B gene mutation analysis
Genomic DNA was extracted from peripheral blood leukocytes by the salting-out method. Patients were first screened for the p.His1069Gln mutation, using primers and probes according to Witt and Landt (2001). DNA samples heterozygous for or without p.His1069Gln mutation were further analyzed by sequencing 21 exons and flanking intronic regions of the ATP7B gene according to Loudianos et al. (1998). PCR products were separated by capillary electrophoresis using the AB Genetic Analyzer 3130xl (Applied Biosystems; Life Technologies).
DNA variants were characterized using Wilson disease mutation database (Bugbee et al., 2001), and the effect of an amino acid substitution on a protein function for novel mutations was assessed using the following software tools: PolyPhen2 version 2.2.2 (Adzhubei et al., 2010), Provean version 1.1.3. (Choi and Chan, 2015), and SIFT (Kumar et al., 2009). The ATP7B gene mutations were named according to the GenBank NCBI reference sequences NM_000053.3 and NG_008806.1.
Statistical analysis
The Kolmogorov-Smirnov test was used to assess the data distribution and, according to its findings, appropriate nonparametric tests were used in the following analyses. Categorical variables were shown through frequencies and percentages, while medians and IQRs were used to describe the quantitative variables. Differences in quantitative variables were analyzed with the Kruskal-Wallis test. Chi-square test was used to analyze differences between investigated groups in categorical variables. All p-values below 0.05 were considered significant. Statistical package IBM SPSS Statistics version 19.0.0.1 (www.spss.com) was used for all statistical procedures.
Results
The analysis of ATP7B gene coding region and flanking intronic regions by sequencing method and/or direct mutation analysis was performed in a Croatian cohort of WD patients. We detected 18 different mutations, 11 of which were missense, 5 frameshift, 1 nonsense, and 1 splice-site mutation (Table 2). In most WD patients, segregation of variants in the ATP7B gene was tested and confirmed by analyzing their parents' DNA. The p.His1069Gln mutation proved to be the most prevalent in the Croatian population with a frequency of 44.7% among the total number of tested alleles, followed by p.Ala1003Thr (8.7%), c.2304dupC (6%), p.Arg616Gln (4.7%), and c.3402delC (4.7%).
Novel mutations are stated in bold letters.
WD, Wilson disease.
Three novel missense mutations, namely the p.Asp1027Asn, the p.Gly1030Ser, and the p.Ser1432Phe in the ATP7B gene were detected in the Croatian population. In patients in whom novel mutations were determined, all 21 exons of the ATP7B gene were analyzed to exclude the presence of other ATP7B mutations.
Thirty-two different ATP7B genotypes were determined in this study (Table 1). Clinical diagnosis of WD was confirmed in 59 patients (78.7%) as mutations were identified in both ATP7B alleles. A total of 23 WD patients (30.7%) were homozygous for p.His1069Gln mutation, 6 patients (8%) were homozygous for other mutations, while 30 patients (40%) were compound heterozygotes for ATP7B mutations. In seven patients (9.3%), mutations in only one ATP7B allele were detected, while no mutations were detected in nine patients (12%). Using the sequencing method combined with direct mutation analysis for p.His1069Gln, we were able to identify disease-causing mutations on 125 (83.3%) out of 150 tested ATP7B alleles.
Polymorphisms, including one silent mutation unreported thus far, that were found in Croatian WD patients are stated in Table 3.
Novel polymorphism is stated in bold letters.
Possible deleterious effect.
Genotype-phenotype correlation in Croatian WD patients
First, WD patients with mutations detected in both ATP7B alleles were classified into three groups according to the effect of mutation: the first group was composed of patients with missense mutations on both ATP7B alleles; the second group comprised patients with one missense mutation with an addition of one other mutation that has significant effect on the protein, such as a frameshift, splice-site, or a nonsense mutation; the third group included patients with both such mutations. WD phenotypes were classified into four groups according to Ferenci et al. (2003): N1, N2, H1, and H2. Statistical analysis of the correlation between the same genotype classification groups and the age at onset of the disease was also performed. No statistically significant correlation between the ATP7B mutation genotype and type of disease presentation, or between genotype and age at onset of the disease was found (p = 0.124 and p = 0.531, respectively).
Second, WD patients were classified into three groups based on the p.His1069Gln mutation, as follows: the first group included patients homozygous for p.His1069Gln, the second group was composed of compound heterozygotes for p.His1069Gln and some other mutation, while patients with other mutations on both alleles comprised the third group. WD phenotype was divided into two groups, hepatic (including H1 and H2 patients) and neurological (including N1 and N2 type of presentation). These genotype groups were also analyzed to investigate the possible correlation with age at onset of the disease. No statistically significant correlation between p.His1069Gln-based mutation genotype classification and type of disease presentation, or between genotype and age at onset of the disease was found (p = 0.195 and p = 0.453, respectively).
Discussion
Results obtained from this study reveal the genetic background in the Croatian cohort of WD patients for the first time. Detailed ATP7B gene analysis results showed a significant number of different mutations and confirmed the genetic heterogeneity among WD patients in Croatia. As expected, p.His1069Gln showed to be the most frequent mutation, as is the case in central and eastern European populations. The following three most frequent mutations are similar to those present in the Serbian population (Loudianos et al., 2003), which is not surprising considering the geographic proximity of the two populations. Frameshift mutation c.4374_4375del in exon 21, which was detected in four patients in this study, was previously detected only in the Serbian population thus far, further facilitating the similarity of WD genetic basis between these two populations (Lepori et al., 2007).
Pathogenicity estimation of novel missense mutations was carried out by using the following software tools: PolyPhen2 version 2.2.2 (Adzhubei et al., 2010), Provean version 1.1.3. (Choi and Chan, 2015), and SIFT (Kumar et al., 2009), and it indicated that all three novel mutations are likely pathogenic. Both p.Asp1027Asn and p.Gly1030Ser are located in exon 14, in the phosphorylation domain of the ATP7B protein, and the highly conserved amino acid residue motif Asp-Lys-Thr-Gly-Thr, which is supposed to have a significant impact on the normal function of the protein by disturbing catalytic phosphorylation. Mutation p.Gly1030Ser was found in two unrelated patients, which indicates that this variant could be more frequent in the Croatian population. Exon 14 appears to be a hot spot for mutations in the Croatian cohort of WD patients.
The third novel mutation, detected in exon 21, was p.Ser1432Phe and it is located in the C-terminal domain of the protein which is not involved in the catalytic function, but might have a role in maintaining the protein's stability or localization during its normal transfer between cellular compartments.
The newly described p.Asp1027Asn mutation was detected in a 26-year-old male patient who had previously undergone liver transplantation due to decompensated liver cirrhosis. The p.Gly1030Ser mutation was detected in a 17-year-old female patient with decompensated liver cirrhosis and in a 30-year-old male patient with neurological manifestations. The p.Ser1432Phe mutation was found in a male patient who presented with early neurological manifestations followed by liver cirrhosis and decompensation due to noncompliance to therapy. The functional effect of novel mutations on the ATP7B protein should also be investigated using further studies.
As the frequency of p.His1069Gln mutation alone is 44.7%, diagnostic sensitivity of WD genetic testing by sequencing the entire coding region of ATP7B compared to direct mutation analysis method for p.His1069Gln was increased by 38.6%, which allowed for 36 more WD patients (48%) to establish their molecular diagnosis of WD. According to the results of this study, the exons of the ATP7B gene that should be analyzed first for mutations in Croatian patients with WD are 5, 8, 13, 14, and 15 where 113 of 125 detected molecular defects are located (90.4%). By analyzing these five exons only, WD diagnosis was confirmed for 50 patients (66.7%) on the molecular level.
Finally, the five most frequent mutations in the Croatian population were present on 103 ATP7B alleles (68.7%), which indicate that a simpler and faster screening method would also be useful for the diagnosis of WD in this population. A new polymorphism c.2994C>T, p.Gly998Gly, which has also been detected in this study, is probably not a pathogenic variant since it causes a silent mutation and is also sufficiently distant from the exon-intron border not to cause irregularities in mRNA splicing.
The lack of statistically significant correlation between WD genotype and the type of disease presentation, as well as correlation between WD genotype and the age at onset of the disease, can be explained by a too small number of subjects, as well as by other genetic and environmental factors, which have a significant impact on age at onset and disease presentation even among siblings with the same WD genetic background (Kegley et al., 2010; Ferenci, 2014).
In conclusion, new mutations in ATP7B gene were discovered, which corroborate its genetic heterogeneity. This study contributes to the enrichment of knowledge of the genetics of WD in Europe, adds to the present understanding of its genotype/phenotype relationship, and provides data for a genetic screening approach for the most frequent mutations in the Croatian population.
Footnotes
Author Disclosure Statement
No competing financial interests exist.
