Theranostics 2016; 6(5):638-649. doi:10.7150/thno.14596

Research Paper

Spectrum and Classification of ATP7B Variants in a Large Cohort of Chinese Patients with Wilson's Disease Guides Genetic Diagnosis

Yi Dong1*, Wang Ni1*, Wan-Jin Chen2*, Bo Wan3*, Gui-Xian Zhao4, Zhu-Qing Shi5, Yue Zhang4, Ning Wang2, Long Yu3, Jian-Feng Xu5, Zhi-Ying Wu1 Corresponding address

1. Department of Neurology and Research Center of Neurology in Second Affiliated Hospital, and the Collaborative Innovation Center for Brain Science, Zhejiang University School of Medicine, Hangzhou, China;
2. Department of Neurology and Institute of Neurology, First Affiliated Hospital, Fujian Medical University, Fuzhou, China;
3. State Key Laboratory of Genetic Engineering, Institutes of Biomedical Sciences, School of Life Science, Fudan University, Shanghai, China;
4. Department of Neurology and Institute of Neurology, Huashan Hospital, Shanghai Medical College, Fudan University, Shanghai, China;
5. Center for Genomic Translation and Health Intervention, School of Public Health, Fudan University, Shanghai, China.
* These authors contributed equally to this work.

This is an open access article distributed under the terms of the Creative Commons Attribution (CC BY-NC) License. See http://ivyspring.com/terms for full terms and conditions.
How to cite this article:
Dong Y, Ni W, Chen WJ, Wan B, Zhao GX, Shi ZQ, Zhang Y, Wang N, Yu L, Xu JF, Wu ZY. Spectrum and Classification of ATP7B Variants in a Large Cohort of Chinese Patients with Wilson's Disease Guides Genetic Diagnosis. Theranostics 2016; 6(5):638-649. doi:10.7150/thno.14596. Available from http://www.thno.org/v06p0638.htm

Abstract

Background: Wilson's disease (WD) is an autosomal recessive disorder of copper metabolism caused by ATP7B pathogenic mutations. The symptoms of WD can be effectively prevented if the affected individuals are identified and intervened early. However, clinical utility of this molecular analysis is challenging due to hundreds of variants with various clinical effects in the gene. Here, we aim to describe the spectrum of ATP7B variants and assess their clinical effects in the Han Chinese population.

Methods: The ATP7B gene was directly sequenced in 632 unrelated WD patients and 503 unrelated healthy individuals. The effects of identified variants were classified according to the American College of Medical Genetics and Genomics (ACMG) standards and guidelines. Different frequency of variants observed in both cases and controls were tested using Chi-square or Fisher's exact tests.

Results: We detected 161 non-synonymous variants in these 632 WD patients, 58 of which were novel. Among these variants, 78, 64, 8, 4, and 7 were classified as 'pathogenic variants', 'likely pathogenic variants', 'variants with uncertain significance', 'likely benign variants', and 'benign variants', respectively. Ninety percent (569/632) of these WD patients can be genetically diagnosed with two or more 'pathogenic' or 'likely pathogenic' variants. The 14 most common disease-causing variants were found at least once in 94% (537/569) of genetically diagnosed patients.

Conclusions: These data expand the spectrum of ATP7B variants and facilitate effective screening for ATP7B variants for early diagnosis of WD and development of individualized treatment regimens.

Keywords: Wilson's disease, ATP7B, variants classification.

Introduction

Wilson's disease (WD, OMIM #277900), or progressive hepatolenticular degeneration, is an autosomal recessive disease caused by pathogenic mutations within the ATP7B gene. 1-3ATP7B encodes a P-type ATPase that is involved in the transport of copper into the plasma protein ceruloplasmin (Cp) and in the excretion of copper from the liver. ATP7B protein malfunction leads to massive accumulation of copper in the liver, brain, and other tissues. The accumulation of copper in these areas can present a wide spectrum of symptoms, including liver cirrhosis; neuronal degeneration of the brain, particularly in the basal ganglia; Kayser-Fleischer (K-F) rings at the corneal limbus; and kidney damage.4 The incidence of WD is estimated to be 1 in 30,000 individuals, and the heterozygous carrier rate is about 1 in 90 among many ethnic groups.5 However, a recent study identified the genetic prevalence of WD as 1:7,026 in the United Kingdom by sequencing ATP7B in 1000 control subjects,6 much higher than the typically reported prevalence of WD of 1:30,000. In addition, WD is thought to be more frequent in Asian populations. The presumed prevalence of WD is approximately 1 in 3,000 in the Korean population7 and the expected frequency of WD is 1 in 5,400 in the Hong Kong Han Chinese population.8

WD is among a limited number of inherited diseases for which symptoms can be prevented if the affected individuals can be identified and intervened early. However, accurate and early clinical diagnoses are often difficult to make due to the wide array of phenotypic variations.9 While genetic analysis is feasible, its clinical utility has been limited due to more than 700 ATP7B reported variants with various clinical effects (https://portal.biobase-international.com/hgmd/pro/gene.php?gene=ATP7B, Human Gene Mutation Database Professional, access date: 20 October, 2015). Therefore, better understanding of the spectrum of ATP7B variants and the clinical effects of these variants are necessary for clinical application.

Here, we reported variants identified from a study that directly sequenced the ATP7B gene in 632 consecutively treated WD patients at three Chinese academic medical centers between 2004 and 2015, and 503 unrelated phenotypically normal individuals. The clinical effects of identified variants were classified according to the American College of Medical Genetics and Genomics (ACMG) standards and guidelines.

Subjects and Methods

Participants

Consecutive patients who sought for diagnosis and treatment of WD between January 15, 2004 and April 30, 2015 at the Departments of Neurology at three leading Chinese academic medical centers were recruited, including Second Affiliated Hospital, Zhejiang University School of Medicine, Hangzhou; First Affiliated Hospital, Fujian Medical University, Fuzhou; Huashan Hospital, Fudan University, Shanghai. Each patient was examined by at least two senior neurologists, including medical history review and physical exams. The patients were clinically diagnosed with WD according to the diagnostic criteria described previously.4 Blood samples of 632 unrelated patients (346 males and 286 females, aged 9 to 56) who have average onset age of 18.5 ± 0.6 ranging from 3 to 48 years and 503 unrelated phenotypically normal controls with no known family history of WD were included in the study. Family history were detected in 16 out of 632 WD patients. All participants are of Han Chinese descent. The probands, or their legal guardians, provided informed consents. This study was approved by three institutional review boards.

Direct DNA Sequencing for Mutation Analysis of ATP7B

Genomic DNA was isolated from peripheral blood lymphocytes using a DNA isolation kit (Qiagen Inc, Valencia, CA). All 21 exons and their flanking regions as well as the 5'-untranslated region (UTR) and promoter regions of the ATP7B gene were sequenced. Exon 2, which is 1,234 base pair (bp) in length, was covered by three overlapping PCR fragments. The primer sequences, annealing temperatures, and sizes of the PCR products of exons 1, 2, 4, 16, and 18 were listed in Table S1. Information related to the remaining primers has been described previously.10 PCR amplification and direct sequencing were performed as previously reported.11 Obtained sequences were compared, and nucleotide changes were numbered according to their position in ATP7B mRNA (NM_000053). Novel variants were further confirmed by sequencing both the forward and reverse strands. The sequence variants were interpreted and classified according to ACMG Standards and Guidelines.12

Statistical Analysis

Different frequency of variants observed in both WD patients and controls was tested using Chi-square test or the Fisher's exact test (for variants with the expected number of subjects below 5). The Bonferroni correction was used to declare statistical significance.

Results

Variants identified in the ATP7B Gene

A total of 173 variants in the coding region and the adjacent splice sites of the ATP7B gene were identified among 632 unrelated WD patients and 503 controls, including 161 non-synonymous and 12 synonymous variants. Among the 161 non-synonymous variants, 150 were detected only in WD patients and 11 were detected in both WD patients and normal controls. To assess link of the latter group of variants with WD, we compared allele frequency for 10 of these 11 variants between WD patients and controls (one of the 11 variants, p.R778L, was excluded because it is a well-established pathogenic variant13). Three variants (p.I390V, p.T935M, and p.V1106I) were considerably more common in WD patients than controls and the differences reached a statistical significance level (p<0.005) after adjusting for 10 multiple tests (0.05/10) (Table 1). The estimated odds ratios (ORs) of these three variants greater than 10 (p.I390V, OR=11.26; p.T935M, OR=77.04; p.V1106I, OR=10.44). For the remaining 7 variants, the allele frequency was similar between cases and controls (p>0.005), with ORs close to 1.00. Homozygotes in controls were found for 5 of these 7 variants (Table 1).

Among 12 synonymous variants, 5 were novel, including c.2145c>T (p.Y715Y), c.3243G>A (p.E1081E), c.3861C>G (p.G1287G), c.4014T>A (p.I1338I), c.4194c>T (p.S1398S). The remaining variants have been previously identified (http://asia.ensembl.org/Homo_sapiens/Transcript/Sequence_cDNA?db=core;g=ENSG00000123191;r=13:51932673-52011494;t=ENST00000242839, access date: 20, October, 2015), including c.747G>A (p.L249L, rs554554415), c.1620c>T (p.L540L,rs145798966), c.2292c>T (p.F764F, rs372979339), c.2310C>G (p.L770L, rs398123136), c.2973G>A (p.T991T, rs1801246), c.3009G>A (p.A1003A, rs1801247) and c.4251A>G (p.T1417T, rs546721020).

Classification of Variants

Among 161 non-synonymous variants, 7 were similar between cases and controls (Table 1), thus classified as benign variants according to the ACMG Standards and Guidelines.12 The remaining 154 non-synonymous variants were only observed in WD patients or significantly more frequent in WD patients (Table 2). They are distributed throughout the ATP7B gene exons 1 to 21 (Figure 1). Among these 154 non-synonymous variants, 18 (11.7%) are nonsense variants, 25 (16.2%) are small deletions or insertions, 11 (7.1%) are splice site variants, and 100 (64.9%) are missense variants. According to the ACMG Standards and Guidelines,12 18 nonsense variants can be classified as 'pathogenic variants'. Among 25 small deletions/insertions, three shift variants including c.2316_2317insCTCTTTGTG (p.V772insLFV), c.2790_2792del(p.I930del) and c.4005_4006insTTATAATGGGTTGCG (p.G1335insLXWVA) can only be classified as 'likely pathogenic variants' due to its in-frame characteristics, other 22 ones can be classified as 'pathogenic variants' as well. For the 11 splice site mutations, six (c.51+1g>a, c.1543+1g>t, c.1708-1g>c, c.1870-2a>g, c.2356-1g>c, and c.3557-2a>g) are predicted to result in exon skipping and lead to the production of a defective protein, therefore can also be classified as 'pathogenic variants'. For the other five splice site variants (c.1543+4a>g, c.1708-5t>g, c.1946+5g>a, c.2447+5g>t, and c.3903+5g>a), because additional functional analysis is required to assess their impact on RNA and protein, they can be classified as 'variants with uncertain significance' at this stage. For the 100 missense variants (Table 3), 32, 61, 4, and 3 can be classified as 'pathogenic variants', 'likely pathogenic variants'', 'likely benign variants'', and 'variants with uncertain significance', respectively, based on results of SIFT, PolyPhen-2, 1000 Genomes Project, Exome Aggregation Consortium and the ACMG Standards and Guidelines. In total, among 161 non-synonymous variants, 78 are classified as 'pathogenic variants', 64 as 'likely pathogenic variants', 8 as 'variants with uncertain significance', 4 as 'likely benign variants' and 7 as 'benign variants'. Therefore, 142 variants could be considered as potential disease-causing variants ('pathogenic variants' and 'likely pathogenic variants') at this stage.

 Table 1 

The Allele Frequency of 10 Non-Synonymous Variants in 632 WD Patients and 503 Controls.

VariantsControlsWD PatientsProbabilities
OR(95%CI)
WildtypeHeterozygoteHomozygoteAllele
Frequency
WildtypeHeterozygoteHomozygoteAllele
Frequency
p.I390V502100.0016191210.0110.003*11.256(1.478-85.743)
p.S406A1462231340.4881532712080.5440.0441.249(1.058-1.474)
p.L456V1502401130.4631573041710.5110.0851.211(1.026-1.430)
p.R832K339591050.267441471440.2650.0450.988(0.819-1.192)
p.I929V4881500.015624800.0060.0420.421(0.178-0.997)
p.T935M502100.0015438810.0711.40E-22*77.044(10.716-553.910)
p.R952K164261780.415206335910.4090.8620.978(0.826-1.157)
p.V1106I502100.0016191300.0100.0049*10.444(1.364-79.969)
p.V1140A195240680.374258286880.3660.7060.965(0.813-1.146)
p.V1297I498500.005628400.0030.5200.636(0.17-2.373)

*Based on Fisher's Exact test due to small observed number of variants.

 Figure 1 

(A) Schematic representation of the ATP7B gene with 21 exons. All here identified variants in WD patients with Chinese Han descent are visualized in the corresponding ATP7B protein regions. Fifty-eight novel variants are indicated in red characters. Cu 1~6, six copper binding domains; Tm 1~8, eight transmembrane domains; Td, a transduction domain converting the ATP hydrolysis energy to the impetus of transporting copper cation; Ch, the transmembrane cation channel; Ph, the phosphorylation domain; ATP loop/binding/hinge, the ATP binding domain. (B) Distribution of ATP7B variants is shown within 21 exons of ATP7B gene.

Theranostics Image (Click on the image to enlarge.)
 Table 2 

One Hundred and Fifty-Four Non-Synonymous Variants were Only Observed in WD Patients or Significantly More Frequent in WD Patients.

Mutation analysisDomainFrequency of MU (%)No. of patients
Nucleotide mutationProtein alterationExonMU/MUWT/MU
c.51+1g>aNa1before Cu10.0801
c.268_271delp.K90FfsX102Cu20.0801
c.314C>Ap.S105X2Cu20.0801
c.367delGp. A123PfsX302Cu20.0801
c.525dupAp.V176SfsX282Cu21.50117
c.588C>Ap.D196E2Cu20.4005
c.592A>Gp.R198G2Cu20.0801
c.695delCp.P232QfsX302Cu30.0801
c.994G>Tp.E332X2Cu3/Cu40.5507
c.1168A>Gp.I390V2Cu3/Cu41.11112
c.1403_1416delp.A468GfsX333Cu4/Cu50.0801
c.1426G>Ap.A476T3Cu4/Cu50.0801
c.1449_1456delp.R483SfsX203Cu50.1602
c.1470C>Ap.C490X3Cu50.2403
c.1531C>Tp.Q511X3Cu52.37226
c.1543+1g>tNa3Cu50.4706
c.1543+4a>gNa3Cu50.0801
c.1544G>Tp.G515V4Cu50.0801
c.1545delTp.G515GfsX94Cu50.0801
c.1552_1553delTCp.S518RfsX154Cu50.0801
c.1639C>Tp.Q547X4Cu50.0801
c.1708-5t>gNa5Cu61.03111
c.1708-1g>cNa5Cu61.50215
c.1760C>Tp.T587M5Cu60.0801
c.1782T>Ap.Y594X5Cu60.0801
c.1803delCp.S602AfsX465Cu60.0801
c.1817T>Gp.V606G5Cu60.1602
c.1820dupAp.F608VfsX25Cu60.1602
c.1846C>Tp.R616Y5Cu60.0801
c.1870-2a>gna6Cu60.0801
c.1925A>Gp.D642G6Cu60.0801
c.1946+5g>ana6Cu6/TM10.0801
c.1950G>Ap.W650Term7Cu6/TM10.0801
c.2038C>Tp.Q680X7Cu6/TM10.0801
c.2043delCp.S681SfsX157TM20.0801
c.2078C>Gp.S693C7TM20.1602
c.2128G>Ap.G710S8TM2/TM30.0801
c.2156A>Gp.Y719C8TM2/TM30.0801
c.2157C>Ap.Y719X8TM2/TM30.2403
c.2192T>Ap.V731E8TM30.0810
c.2195T>Cp.L732P8TM30.0801
c.2223T>Ap.Y741X8TM30.0801
c.2251G>Tp.A751S8TM30.0801
c.2261A>Gp.E754G8TM3/TM40.0801
c.2267C>Gp.A756G8TM3/TM40.0801
c.2294A>Gp.D765G8TM40.2403
c.2297C>Tp.T766M8TM40.1602
c.2304dupCp.M769HfsX268TM40.87011
c.2305A>Gp.M769V8TM40.0801
c.2308C>Tp.L770F8TM40.0801
c.2316_2317ins
CTCTTTGTG
p.V772insLFV8TM40.0801
c.2332C>Tp.R778W8TM40.1602
c.2333G>Tp.R778L8TM429.6764247
c.2333G>Ap.R778Q8TM41.98319
c.2336G>Ap.W779X8TM40.0801
c.2341G>Ap.E781K8TM40.0801
c.2356-1g>cna9TM4/Td0.0801
c.2383C>Tp.L795F9TM4/Td0.0801
c.2390C>Tp.S797F9TM4/Td0.0801
c.2447+5g>tNa9TM4/Td0.0801
c.2455C>Tp.Q819X10TM4/Td0.0801
c.2525A>Gp.D842G10Td0.0801
c.2587C>Tp.P863S11Td0.0801
c.2605G>Ap.G869R11Td0.1602
c.2620G>Cp.A874P11Td/TM50.0801
c.2621C>Tp.A874V11Td/TM53.56339
c.2648_2649delp.V883AfsX311Td/TM50.0801
c.2662A>Cp.T888P11Td/TM50.3204
c.2668G>Ap.V890M11Td/TM50.1602
c.2740C>Tp.Q914X12Td/TM50.0801
c.2755C>Gp.R919G12Td/TM51.98025
c.2790_2792delp.I930del12TM50.4706
c.2794_2795insGTp.S932CfsX412TM50.0801
c.2804C>Tp.T935M12TM57.12188
c.2810delTp.V937GfsX512TM50.6308
c.2827G>Ap.G943S12TM50.5515
c.2828G>Ap.G943D12TM52.21126
c.2848G>Tp.V950F12TM5/TM60.0801
c.2853_2856delp.Q951HfsX1512TM5/TM60.0801
c.2905C>Tp.R969W13TM60.0801
c.2924C>Ap.S975Y13TM60.79010
c.2930C>Tp.T977M13TM60.0801
c.2957C>Tp.S986F13TM60.0801
c.2975C>Tp.P992L13TM6/Ph14.5626132
c.3007G>Ap.A1003T13TM6/Ph0.1602
c.3008C>Tp.A1003V13TM6/Ph0.0801
c.3010C>Tp.Q1004X13TM6/Ph0.0801
c.3029A>Cp.K1010T13TM6/Ph0.1602
c.3044T>Cp.L1015P13TM6/Ph0.0801
c.3053C>Tp.A1018V13TM6/Ph0.2403
c.3083A>Gp.K1028R14Ph0.0801
c.3087delTp.G1030AfsX9114Ph0.0801
c.3089G>Ap.G1030D14Ph0.4005
c.3095T>Cp.I1032T14Ph0.0801
c.3098C>Tp.T1033I14Ph0.0801
c.3104G>Tp.G1035V14Ph0.0801
c.3122G>Cp.R1041P14ATP loop0.1610
c.3121C>Tp.R1041W14ATP loop0.0801
c.3140A>Tp.D1047V14ATP loop0.2411
c.3155C>Tp.P1052L14ATP loop0.1602
c.3209C>Gp.P1070R14ATP loop0.0801
c.3221C>Tp.A1074V14ATP loop0.0801
c.3263T>Cp.L1088S15ATP loop0.3204
c.3271T>Cp.C1091R15ATP loop0.0801
c.3274A>Cp.T1092P15ATP loop0.0810
c.3311G>Ap.C1104Y15ATP loop0.0801
c.3316G>Ap.V1106I15ATP loop1.03013
c.3368C>Tp.P1123L15ATP loop0.0801
c.3377_3378delACp.H1126PfsX315ATP loop0.0801
c.3426G>Cp.Q1142H16ATP loop1.58118
c.3443T>Cp.I1148T16ATP loop3.32140
c.3446G>Ap.G1149E16ATP loop0.0801
c.3451C>Tp.R1151C16ATP loop0.0801
c.3451C>Gp.R1151G16ATP loop0.0801
c.3452G>Ap.R1151H16ATP loop0.1602
c.3459G>Tp.W1153C16ATP loop0.0801
c.3517G>Ap.E1173K16ATP loop0.7918
c.3532A>Gp.T1178A16ATP loop0.1602
c.3557-2a>gna17ATP loop0.0801
c.3563T>Gp.L1188R17ATP loop0.0801
c.3587A>Gp.D1196G17ATP loop0.0801
c.3605C>Gp.A1202G17ATP loop0.0801
c.3646G>Ap.V1216M17ATP bind1.98221
c.3653T>Cp.L1218P17ATP bind0.0801
c.3700delGp.V1234LfsX9618ATP bind0.3204
c.3715G>Tp.V1239F18ATP bind0.2403
c.3733C>Gp.P1245A18ATP bind0.0801
c.3741C>Gp.H1247Q18ATP bind0.0801
c.3766_3767dupCAp.Q1256PfsX7518ATP bind0.0801
c.3776G>Tp.G1259V18ATP hinge0.0811
c.3791T>Cp.M1264T18ATP hinge0.0801
c.3802G>Ap.G1268R18ATP hinge0.0801
c.3809A>Gp.N1270S18ATP hinge2.22028
c.3818C>Tp.P1273L18ATP hinge0.0801
c.3818C>Ap.P1273Q18ATP hinge0.1602
c.3836A>Gp.D1279G18ATP hinge0.2403
c.3848C>Tp.A1283V18ATP hinge0.0801
c.3859G>Ap.G1287S18ATP hinge0.0801
c.3877G>Ap.E1293K18ATP hinge0.0801
c.3884C>Tp.A1295V18ATPhinge/TM70.4005
c.3896T>Gp.L1299R18ATP hinge/TM70.0801
c.3901_3902insAp.R1301KfsX318ATP hinge/TM70.0801
c.3903+5g>ana18ATP hinge/TM70.0801
c.3955C>Tp.R1319X19ATP hinge/TM70.2403
c.3982G>Ap.A1328T19TM70.1602
c.4003G>Cp.G1335R19TM70.2403
c.4005_4006insTTATAATGGGTTGCGp.G1335insLXWVA19TM70.1602
c.4057T>Cp.W1353R20TM80.1602
c.4059G>Ap.W1353X20TM80.0801
c.4064G>Ap.G1355D20TM80.0801
c.4112T>Cp.L1371P20TM80.4005
c.4114C>Tp.Q1372X20TM80.6308
c.4162delGp.A1388RfsX521after TM80.0801
c.4272T>Gp.Y1424X21after TM80.0801

MU: mutant; WT: wild type; Cu: copper binding domain; Td: transduction domain converting energy from ATP hydrolysis to cation transportation; Tm: transmembrane domain; Ch: ion channel; Ph: phosphorylation loop.

 Table 3 

Classified mutations within 100 missense mutations.

Nucleotide mutationProtein
alteration
ExonSiftPolyPhen 21000gExAcClassificationEvidence of
pathogenicity
ScorePredictionScorePrediction
c.588C>Ap.D196E21Tolerated0.924Probably damaging05US4*PP
c.592A>Gp.R198G20.05Damaging0.994Probably damaging00LP1*PM,4*PP
c.1168A>Gp.I390V20.57Tolerated0.001Benign02LB1*BS,1*BP
c.1426G>Ap.A476T30.63Tolerated0.503Possibly damaging253US1*BS
c.1544G>Tp.G515V40Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.1760C>Tp.T587M50.19Tolerated0.099Benign00LB1*BS,1*BP
c.1817T>Gp.V606G50Damaging1Probably damaging00LP1*PM,5*PP
c.1846C>Tp.R616Y50Damaging1Probably damaging04LP2*PM,5*PP
c.1925A>Gp.D642G60.01Damaging1Probably damaging00LP1*PM,4*PP
c.2078C>Gp.S693C70.01Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.2128G>Ap.G710S80.02Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.2156A>Gp.Y719C80Damaging0.991Probably damaging00LP2*PM,4*PP
c.2192T>Ap.V731E80Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.2195T>Cp.L732P80Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.2251G>Tp.A751S80.16Tolerated0.999Probably damaging00LP2*PM,3*PP
c.2261A>Gp.E754G80.34Tolerated0.999Probably damaging00LP2*PM,3*PP
c.2267C>Gp.A756G80.21Tolerated1Probably damaging00P1*PS,2*PM,5*PP
c.2294A>Gp.D765G80Damaging1Probably damaging00LP3*PM,5*PP
c.2297C>Tp.T766M80Damaging1Probably damaging01P1*PS,2*PM,5*PP
c.2305A>G
c.2308C>T
p.M769V
p.L770F
8
8
00Damaging
Damaging
1
1
Probably damaging
Probably damaging
008
0
LP
LP
3*PM,5*PP
2*PM,4*PP
c.2332C>Tp.R778W80Damaging1Probably damaging05P2*PS,3*PM
c.2333G>Ap.R778Q80Damaging1Probably damaging05P2*PS,3*PM
c.2333G>Tp.R778L80Damaging1Probably damaging05P2*PS,3*PM
c.2341G>Ap.E781K80Damaging1Probably damaging00P1*PS,2*PM,4*PP
c.2383C>Tp.L795F90Damaging1Probably damaging04P1*PS,3*PM
c.2390C>Tp.S797F90Damaging1Probably damaging00LP1*PM,4*PP
c.2525A>Gp.D842G100Damaging1Probably damaging00LP1*PM,4*PP
c.2587C>Tp.P863S110Damaging1Probably damaging00LP1*PM,4*PP
c.2605G>Ap.G869R110Damaging1Probably damaging582LP2*PM,5*PP
c.2620G>Cp.A874P110.01Damaging1Probably damaging00LP2*PM,5*PP
c.2621C>Tp.A874V110.01Damaging1Probably damaging09LP2*PM,5*PP
c.2662A>Cp.T888P110Damaging1Probably damaging00LP2*PM,5*PP
c.2668G>Ap.V890M110.01Damaging1Probably damaging00LP2*PM,5*PP
c.2755C>Gp.R919G120.01Damaging0.988Probably damaging02P1*PS,2*PM,5*PP
c.2804C>Tp.T935M120Damaging1Probably damaging021P1*PS,1*PM,5*PP
c.2827G>Ap.G943S120.16Tolerated1Probably damaging02P1*PS,1*PM,4*PP
c.2828G>Ap.G943D120Damaging1Probably damaging02P1*PS,2*PM,5*PP
c.2848G>Tp.V950F120.01Damaging0.994Probably damaging00LP2*PM,4*PP
c.2905C>Tp.R969W130.02Damaging1Probably damaging05LP2*PM,5*PP
c.2924C>Ap.S975Y130.01Damaging1Probably damaging01LP2*PM,5*PP
c.2930C>Tp.T977M130Damaging1Probably damaging010P1*PS,2*PM,5*PP
c.2957C>Tp.S986F130Damaging1Probably damaging00LP2*PM,5*PP
c.2975C>Tp.P992L130Damaging1Probably damaging05P2*PS,3*PM
c.3007G>Ap.A1003T130Damaging1Probably damaging12LP2*PM,5*PP
c.3008C>Tp.A1003V130Damaging1Probably damaging07LP2*PM,5*PP
c.3029A>Cp.K1010T130Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.3044T>Cp.L1015P130Damaging1Probably damaging00LP2*PM,4*PP
c.3053C>Tp.A1018V130.07Tolerated1Probably damaging03LP2*PM,5*PP
c.3083A>Gp.K1028R140.02Damaging0.998Probably damaging00LP2*PM,4*PP
c.3089G>Ap.G1030D140Damaging1Probably damaging00LP2*PM,5*PP
c.3095T>Cp.I1032T140Damaging0.998Probably damaging00LP2*PM,4*PP
c.3098C>Tp.T1033I140.01Damaging1Probably damaging00P1*PS,2*PM,4*PP
c.3104G>Tp.G1035V140Damaging1Probably damaging00LP2*PM,5*PP
c.3121C>Tp.R1041W140Damaging1Probably damaging05LP2*PM,5*PP
c.3122G>Cp.R1041P140.01Damaging1Probably damaging00LP2*PM,5*PP
c.3140A>Tp.D1047V140.13Tolerated0.997Probably damaging00LP2*PM,4*PP
c.3155C>Tp.P1052L140.06Tolerated0.998Probably damaging01LP2*PM,4*PP
c.3209C>Gp.P1070R140Damaging1Probably damaging00LP2*PM,4*PP
c.3221C>Tp.A1074V140Damaging1Probably damaging00LP2*PM,5*PP
c.3263T>Cp.L1088S150.2Tolerated1Probably damaging00LP1*PS,1*PM,3*PP
c.3271T>Cp.C1091R150Damaging0.96Probably damaging00LP1*PM,4*PP
c.3274A>Cp.T1092P150.11Tolerated0.832Possibly damaging00US1*PM,3*PP
c.3311G>Ap.C1104Y150Damaging1Probably damaging01P1*PS,2*PM,5*PP
c.3316G>Ap.V1106I150.15Tolerated0.984Probably damaging216P1*PS,2*PM,4*PP
c.3368C>Tp.P1123L150.31Tolerated0.025Benign025LB1*BS,1*BP
c.3426G>Cp.Q1142H160.16Tolerated0.007Benign03LB1*BS,1*BP
c.3443T>Cp.I1148T160Damaging0.999Probably damaging05LP1*PM,5*PP
c.3446G>Ap.G1149E160Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.3451C>Gp.R1151G160Damaging1Probably damaging00P1*PS,1*PM,4*PP
c.3451C>Tp.R1151C160Damaging1Probably damaging05P1*PS,1*PM,4*PP
c.3452G>Ap.R1151H160.01Damaging1Probably damaging02LP3*PM,5*PP
c.3459G>Tp.W1153C160Damaging1Probably damaging00LP3*PM,5*PP
c.3517G>Ap.E1173K160Damaging1Probably damaging01LP3*PM,5*PP
c.3532A>Gp.T1178A160Damaging0.988Probably damaging00LP1*PM,5*PP
c.3563T>Gp.L1188R170Damaging0.998Probably damaging00LP1*PM,4*PP
c.3587A>Gp.D1196G170Damaging1Probably damaging00LP1*PM,4*PP
c.3605C>Gp.A1202G170Damaging1Probably damaging00LP1*PM,4*PP
c.3646G>Ap.V1216M170Damaging1Probably damaging05LP1*PM,4*PP
c.3653T>Cp.L1218P170Damaging0.999Probably damaging00LP1*PM,4*PP
c.3715G>Tp.V1239F180Damaging0.998Probably damaging02P1*PS,1*PM,4*PP
c.3733C>Gp.P1245A180Damaging1Probably damaging00LP1*PM,4*PP
c.3741C>Gp.H1247Q180.17Tolerated1Probably damaging01LP1*PM,4*PP
c.3776G>Tp.G1259V180Damaging0.997Probably damaging00LP1*PM,4*PP
c.3791T>Cp.M1264T180Damaging0.998Probably damaging00LP1*PM,4*PP
c.3802G>Ap.G1268R180Damaging1Probably damaging00LP2*PM,5*PP
c.3809A>Gp.N1270S180Damaging1Probably damaging018P1*PS,2*PM,5*PP
c.3818C>Ap.P1273Q180Damaging1Probably damaging00P1*PS,2*PM,5*PP
c.3818C>Tp.P1273L180Damaging1Probably damaging04P1*PS,2*PM,5*PP
c.3836A>Gp.D1279G180.01Damaging1Probably damaging02LP2*PM,5*PP
c.3848C>Tp.A1283V180Damaging1Probably damaging00LP2*PM,5*PP
c.3859G>Ap.G1287S180Damaging1Probably damaging03P1*PS,2*PM,5*PP
c.3877G>Ap.E1293K180Damaging1Probably damaging01LP2*PM,5*PP
c.3884C>Tp.A1295V180Damaging1Probably damaging00P1*PS,2*PM,4*PP
c.3896T>Gp.L1299R180Damaging0.999Probably damaging00P1*PS,2*PM,4*PP
c.3982G>Ap.A1328T190Damaging1Probably damaging00LP1*PM,5*PP
c.4003G>Cp.G1335R190Damaging1Probably damaging00LP1*PM,5*PP
c.4057T>Cp.W1353R200Damaging1Probably damaging00LP2*PM,5*PP
c.4064G>Ap.G1355D200Damaging1Probably damaging00LP2*PM,5*PP
c.4112T>Cp.L1371P200.01Damaging1Probably damaging00LP1*PM,5*PP

PolyPhen 2/ SIFT: Software prediction programs used for sequence variant effect explanation; 1000g: 1000 Genomes Project; ExAc: Exome Aggregation Consortium.

Novel Variants

Among the 161 non-synonymous variants, 58 are novel and the chromatograms of these variants were illustrated in Figures 2A, 2B, 2C and 2D, grouped by the type of variation. None of which was found in the 503 normal individuals. Of these 58 novel variants, 16 are small deletions or insertions, 4 are nonsense variants, 6 are splice site variants, and 32 are missense variants. The other 103 variants have been recorded in the professional version of HGMD (Human Gene Mutation Database Professional, access date: 20 October, 2015).

Among the 16 small deletions or insertions, 13 alter the reading frame of ATP7B and result in production of a truncated dysfunctional protein. They are c.268_271del (p.K90FfsX10), c.367delG (p.A123PfsX30), c.695delC (p.P232QfsX30), c.1403_1416del (p.A468GfsX33), c.1545delT (p.G515GfsX9), c.1552_1553delTC (p.S518RfsX15), c.2043delC (p.S681SfsX15), c.2794_2795insGT (p.S932CfsX4), c.2853_2856delGAGA (p.Q951HfsX15), c.3377_3378delAC (p.H1126PfsX3), c.3766_3767dupCA (p.Q1256PfsX75), c.3901_3902insA (p.R1301KfsX3) and c.4162delG (p.A1388RfsX5). The other 3 small deletions or insertions, c.2316_2317insCTCTTTGTG (p.V772insLFV),c.2790_2792delCAT (p.I930del) and c.4005_4006insTTATAATGGGTTGCG (p.G1335insLXWVA) do not change the reading frame of ATP7B, but cause prolonged (p.V772insLFV and p.G1335insLXWVA) or shortened (p.I930del) dysfunctional ATP7B protein. Similarly, the 4 nonsense variants could result in production of a shortened, dysfunctional protein.

 Figure 2 

Chromatograms of 58 novel ATP7B variants identified in the present study. The lower chromatogram in each frame represents the variant, while the upper one represents the normal sequence. The c.4162delG variant is shown in reverse sequence, while the other 57 variants are illustrated in forward sequence. 2A, 2B, 2C and 2D respectively illustrates 32 missense changes, 16 small deletions or insertions, 6 splicing site variants as well as 4 nonsense variants.

Theranostics Image (Click on the image to enlarge.)
 Figure 3 

Homology comparisons of ATP7B protein sequences. The highlighted zones respectively indicate 32 novel missense variant sites among 8 species.

Theranostics Image (Click on the image to enlarge.)

For six splice site variants, three (c.51+1g>a, c.2356-1g>c and c.3557-2a>g) are classified as 'pathogenic variants', the other three (c.1543+4a>g, c.1946+5g>a, 3903+5g>a) are 'variants with uncertain significance'. In addition, among the 32 novel missense variants, 2 (p.T587M, p.P1123L) are classified as 'likely benign variants', 2 (p.A476T, p.T1092P) as 'variants with uncertain significance', 5 (p.E781K, p.T1033I, p.R1151G, p.V1239F, p.L1299R) as 'pathogenic variants' and the other 23 as 'likely pathogenic variants' (Table 3). A homology search of the ATP7B protein in different species demonstrated that 32 missense variants occur within highly conserved regions of ATP7B protein (Figure 3). In total, among 58 novel variants, 25 are classified as 'pathogenic variants', 26 as 'likely pathogenic variants', 5 as 'variants with uncertain significance' and 2 as 'likely benign variants'.

Most Common Variants

Among the 632 WD patients, 569 patients (90%) were identified with homozygous or compound heterozygous potential disease-causing variants, 58 patients (9%) with one heterozygous variant and other 5 patients (1%) did not have any potential disease-causing variant. Therefore, 90% (569/632) of patients can be genetically diagnosed with WD. Among the 142 potential disease-causing variants including 78 'pathogenic variants' and 64 'likely pathogenic variants', 14 were relatively common in the WD patient cohort, each with an allelic frequency 1% or higher. These 14 most common disease-causing variants were found in 94% (537/569) of genetically diagnosed WD patients with two or more 'pathogenic' or 'likely pathogenic' variants. The allelic frequencies and numbers of patients with each of these variant are presented in Table 4. Notably, the three most prevalent variants, p.R778L, p.P992L and p.T935M, were detected at least once in 78% (445/569) of genetically diagnosed WD patients. The allelic frequencies of p.R778L, p.P992L and p.T935M are 0.319, 0.155 and 0.077, respectively; all of which are 'pathogenic variants'.

Patients with More than Two Variants

Six patients carried three disease-causing variants and their 3-variant genotypes are described in Figure 4. As two patients had the same genotype, there are 5 unique 3-variant genotypes. Four of these 3-variant genotype include at least one 'likely pathogenic' variant, while one genotype (T935M+G943D+V1106I) consists of three 'pathogenic' variants.

 Table 4 

Common disease-causing variants within ATP7B among 569 WD patients.

MutationDomain affectedNumber of patientsAllelic
Frequencies
Classification
WWWMMM
p.R778LTM4270235640.319Pathogenic
p.P992LTM6/Ph419124260.155Pathogenic
p.T935MTM54828610.077Pathogenic
p.A874VTd/TM55293730.038Likely pathogenic
p.I1148TATP loop5303810.035Likely pathogenic
p.Q511XCu55412620.026Pathogenic
p.G943DTM55432510.024Pathogenic
p.N1270SATP hinge5442500.022Pathogenic
p.R778QTM45481830.021Pathogenic
p.R919GTd/TM55452400.021Pathogenic
p.V1216MATP bind5472020.021Likely pathogenic
p.V176SfsX28Cu25511710.017Pathogenic
c.1708-1g>cCu65531420.016Pathogenic
p.V1106IATP loop5561300.011Pathogenic

WW: neither of chromosome carries mutation; WM: one chromosome carries mutation;

MM: both chromosomes carry mutations.

 Figure 4 

Three distinct ATP7B disease-causing variants co-occur in 6 WD patients.

Theranostics Image (Click on the image to enlarge.)

Discussion

In the present study, the prevalence of ATP7B variants was systematically investigated in the largest Chinese WD cohort to date, with 632 patients and 503 normal controls. One hundred and sixty-one non-synonymous variants within the ATP7B gene were found in WD patients, including 58 novel variants. This study catalogs ATP7B variants in Han Chinese WD patients and expands the spectrum of ATP7B variants.

Another major contribution of the study is the classification of ATP7B variants. Based on the type of alterations, their predicted impact, and frequency between WD cases and controls, these 161 non-synonymous variants are classified as 'pathogenic variants' (N=78), 'likely pathogenic variants' (N=64), 'variants with uncertain significance' (N=8), 'likely benign variants' (N=4), and 'benign variants' (N=7). The observation that 90% (569/632) of these WD patients had two or more 'pathogenic' or 'likely pathogenic' variants demonstrates the clinical utility of the catalog and classification of ATP7B variants.

The 14 most common disease-causing variants were found at least once in 94% (537/569) of genetically diagnosed patients. Notably, the three most prevalent pathogenic variants, p.R778L, p.P992L and p.T935M, were detected in 78% (445/569) of the patients. These results demonstrate the feasibility of developing a rapid and cost-effective genetic test such as multiplex allele-specific PCR to screen for WD.

We found 6 unrelated WD patients (1.1%) carrying three disease-causing variants, and which was reported previously in Caucasian populations.6 This observation is important because it suggests that two disease-causing variants may reside in a chromosome. The implication is that if an individual carries two disease-causing variants but his/her clinical features do not support the diagnosis of WD, it is necessary to test these variants in parents to identify the specific location of them. If the two variants reside in a chromosome, the individual is a heterozygous carrier and cannot be genetically diagnosed. If two variants reside in different chromosomes, then a genetic diagnosis of WD can be made.

According to ACMG Standards and Guidelines,12 we classified 4 variants as 'likely benign' variants, one of which (p.Q1142H) was previously considered as a pathogenic variant.14 Tsai et al. found that Q1142 was located at the ATP pocket of ATP7B protein, and the amino acid replacement of the site would directly disrupt ATP7B function. However, the amino acid change in p.Q1142H was predicted to be tolerated by SIFT (score: 0.16) and benign by PolyPhen-2 (score: 0) software programs, this variant should be considered as 'likely benign'. This finding is important because the variant is relatively frequent and found in 19 WD patients (18 heterozygous and one homozygous). Caution should be made in interpreting the carrier of the Q1142 variant to avoid misdiagnosis and unnecessary treatment.

Sixty-three clinically diagnosed WD patients cannot be genetically diagnosed with WD because they were detected to carry only one disease-causing variant or none. Several factors may induce the observation. First, a subset of these 63 patients may have been misdiagnosed. For example, these patients may have been included in the patient cohort on the basis of false positive K-F rings. Second, large hemizygous deletions may occur in a subset of these patients and these alterations are difficult to be detected using the Sanger sequencing method. In fact, a hemizygous large deletion spanning the exon 20 and adjacent introns has previously been reported.15 Other methods such as multiplex ligation-dependent probe amplification (MLPA) may be used to detect large deletions. Third, other genetic alterations outside the ATP7B coding region and adjacent splice sites as well as other cellular factors associated with WD may contribute to the clinical development of WD in these patients.16

In conclusion, the current study considerably expands the spectrum of ATP7B variants and provides classification of their clinical effects. These results improve the genetic diagnosis of suspected WD patients and facilitate genetic screening for WD among asymptomatic children in the general Chinese population.

Abbreviations

Cp: ceruloplasmin; K-F: Kayser-Fleischer; ACMG: American College of Medical Genetics and Genomics; HGMD: Human Gene Mutation Database; MLPA: multiplex ligation-dependent probe amplification.

Supplementary Material

Attachment

Supplementary Table S1.

Acknowledgements

This work was supported by a grant from the National Natural Science Foundation of China to Zhi-Ying Wu (81125009). We gratefully acknowledge all participants for their help and willingness to participate this study.

Conflicts of Interest

All authors reported no biomedical financial interests or potential conflicts of interest.

References

1. Bull PC, Thomas GR, Rommens JM, Forbes JR, Cox DW. The Wilson disease gene is a putative copper transporting P-type ATPase similar to the Menkes gene. Nat Genet. 1993;5:327-337

2. Tanzi RE, Petrukhin K, Chernov I, Pellequer JL, Wasco W, Ross B. et al. The Wilson disease gene is a copper transporting ATPase with homology to the Menkes disease gene. Nat Genet. 1993;5:344-350

3. Yamaguchi Y, Heiny ME, Gitlin JD. Isolation and characterization of a human liver cDNA as a candidate gene for Wilson disease. Biochem Biophys Res Commun. 1993;197:271-277

4. Scheinberg IH, Sternlieb I. Wilson's disease. In: (ed.) Smith LH Jr. Major Problems in Internal Medicine. Vol 23. Philadelphia, Pa: WB Saunders Co. 1984

5. Ala A, Walker AP, Ashkan K, Dooley JS, Schilsky ML. Wilson's disease. Lancet. 2007;369:397-408

6. Coffey AJ, Durkie M, Hague S, Mclay K, Emmerson J, Lo C. et al. A genetic study of Wilson's disease in the United Kingdom. Brain. 2013;136:1476-1487

7. Park HD, Ki CS, Lee SY, Kim JW. Carrier frequency of the R778L, A874V, and N1270S mutations in the ATP7B gene in a Korean population. Clin Genet. 2009;75:405-407

8. Mak CM, Lam CW, Tam S, Lai CL, Chan LY, Fan ST. et al. Mutational analysis of 65 Wilson disease patients in Hong Kong Chinese: Identification of 17 novel mutations and its genetic heterogeneity. J Hum Genet. 2008;53:55-63

9. Saito T. An assessment of efficiency in potential screening in Wilson disease. J Epidemiol Community Health. 1981;35:274-280

10. Thomas GR, Forbes JR, Roberts EA, Walshe JM, Cox DW. The Wilson disease gene: spectrum of mutations and their consequences. Nat Genet. 1995;9:210-217

11. Wu ZY, Zhao GX, Chen WJ, Wang N, Wan B, Lin MT. et al. Mutation analysis of 218 Chinese patients with Wilson disease revealed no correlation between the canine copper toxicosis gene MURR1 and Wilson disease. J Mol Med. 2006;84:438-442

12. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al.Standards, guidelines for the interpretation of sequence variants. a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med. 2015;17:405-424

13. Wu ZY, Wang N, Lin MT, Fang L, Murong SX, Yu L. Mutation analysis and the correlation between genotype and phenotype of Arg778Leu mutation in Chinese patients with Wilson disease. Arch Neurol. 2001;58:971-976

14. Tsai CH, Tsai FJ, Wu JY, Chang JG, Lee CC, Lin SP. et al. Mutation analysis of Wilson disease in Taiwan and description of six new mutations. Hum Mutat. 1998;12:370-376

15. Moller LB, Ott P, Lund C, Horn N. Homozygosity for a gross partial gene deletion of the C-terminal end of ATP7B in a Wilson patient with hepatic and no neurological manifestations. Am J Med Genet A. 2005;138:340-343

16. Mufti AR, Burstein E, Csomos RA. XIAP Is a copper binding protein deregulated in Wilson's disease and other copper toxicosis disorders. Mol Cell. 2006;21:775-785

Author contact

Corresponding address Corresponding author: Prof. Zhi-Ying Wu, Department of Neurology and Research Center of Neurology in Second Affiliated Hospital, and the Collaborative Innovation Center for Brain Science, Zhejiang University School of Medicine, 88 Jiefang Rd, Hangzhou, 310009, China. E-mail: zhiyingwuedu.cn.


Received 2015-12-3
Accepted 2016-1-27
Published 2016-3-3