Background. HTRA2 has already been nominated as PARK13 which may cause Parkinson’s disease, though there are still discrepancies among these results. Recently, Gulsuner et al.’s study found that HTRA2 p.G399S is responsible for hereditary essential tremor and homozygotes of this allele develop Parkinson’s disease by examining a six-generation family segregating essential tremor and essential tremor coexisting with Parkinson’s disease. We performed this study to validate the condition of HTRA2 gene in Chinese familial essential tremor and familial Parkinson’s disease patients, especially essential tremor. Methods. We directly sequenced all eight exons, exon-intron boundaries, and part of the introns in 101 familial essential tremor patients, 105 familial Parkinson’s disease patients, and 100 healthy controls. Results. No exonic variant was identified, while one exon-intron boundary variant (rs2241028) and one intron variant (rs2241027) were detected, both with no clinical significance and uncertain function. There was no difference in allele, genotype, and haplotype between groups. Conclusions. HTRA2 exonic variant might be rare among Chinese Parkinson’s disease and essential tremor patients with family history, and HTRA2 may not be the cause of familial Parkinson’s disease and essential tremor in China.
National Natural Science Fund8143002291332107813714071. Introduction
As two of the most prevalent tremor disorders, essential tremor (ET) and Parkinson’s disease (PD), which are estimated to constitute 0.9% and 0.3% of worldwide population, respectively, are considered as distinctively different entities formerly [1, 2]. Several lines of evidence showed that there are remarkable overlaps in clinical features, epidemiology, imaging, genetics, and pathology between PD and ET, including a fourfold increase of risk developing Parkinson’s disease in essential tremor cases [3, 4].
ET is widely regarded as caused by genetic with no disease-causing gene ever been focused; Contrarily, though PD is mainly sporadic, up to now 22 PARK loci have been identified [5, 6]. To be specific, 50% of ET patients demonstrate familial aggregation, while less than 15% of PD patients have affected first-degree relatives [7–9]. Due to the overlap phenomena between ET and PD, investigations into the relationship between PD risk variants and ET patients have been done, involving LINGO1, LINGO2, LRRK2, SLC1A2, and HTRA2 genes [3, 10–12].
HTRA2 has already been nominated as PARK13 which may cause Parkinson’s disease, though there are still discrepancies among these results. Recently, a research by Gulsuner and colleagues examining a six-generation family segregating ET and ET coexisting with PD revealed that HTRA2 p.G399S is responsible for hereditary essential tremor and homozygotes for this allele develop Parkinson’s disease [13]. Replications conducted in Western Norway and Asian population to address the association between p.G399S and ET, PD, ET/PD, and tremulous cervical dystonia failed to reach a consensus [14, 15]. In addition, report from a small sample (29 FETs) in Germany adopting coding exon Sanger sequencing did not reconfirm it either [16]. To validate the condition in Chinese familial essential tremor (FET) and familial Parkinson’s disease (FPD) patients, we performed a Sanger sequencing of eight exons and exon-intron boundaries of HTRA2 instead of just one variant (p.G399S).
2. Methods2.1. Patients
This study enrolled 221 unrelated Chinese patients, including 105 PD patients with autosomal dominant inheritance (2 or more affected relatives in 2 consecutive generations), 101 ET patients with family history, and 15 patients of ET coexisting with PD. All patients were from the Movement Disorder Clinic of Department of Neurology at Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine. PD and ET patients were diagnosed by senior movement disorder specialist on the basis of MDS clinical diagnostic criteria for Parkinson’s disease and Consensus Statement on Tremor of the Movement Disorders Society, respectively [17, 18]. Patients presenting secondary Parkinsonism, Parkinson-plus syndrome, or hyperthyroidism were excluded from the study. We also included 100 healthy controls without any symptom of movement disorders. The demographic information of patients is shown in Table 1. We received approval from the Ethics Committee of Ruijin Hospital affiliated to Shanghai Jiao Tong University School of Medicine. Written informed consents were obtained from all patients and controls participating in the study as well.
Demographics of participants.
Details
FET
FPD
ET-PD
Control
Total
101
105
15
100
Agea (range, pb value)
61.24 ± 12.62 (28–90, 0.12)
59.28 ± 11.21 (36–84, 0.86)
67.80 ± 8.65 (56–79, NA)
59.06 ± 6.21 (49–74, NA)
Male/female, pb value
51/50, 0.52
53/52, 0.52
12/3, N/A
46/54, N/A
N/A: not applicable; adata are mean ± SD; bdata are compared with control; FPD: familial PD; FET: familial ET; ET-PD: ET coexisting with PD.
2.2. DNA Sequencing and Mutation Analysis
Genomic DNA was extracted from venous blood applying standardized phenol/chloroform extraction method from patients and controls. The 8 coding sequences, exon-intron boundaries, and part of introns were sequenced by Sanger sequencing in 4 products of PCR (polymerase chain reaction) amplification using 4 pairs of primers (Table 2). DNASTAR Lasergene MegAlign (v7.1.0) and Chromas (v2.33) were used to conduct sequence alignment, and the chromatograms were double checked to avoid missing any variants. Variants detected were searched in NCBI to get access to their clinical significance and MAF in ExAC and 1000 Genomes Projects database.
Primers of HTRA2.
Name
Forward
Reverse
Products
1
GTC TCA CAA CTC GCG TCC G
GCC TGA AAT GGA GGG AAA GCA
Exon 1 and boundaries
2
TCG AGA TCC TGG ACC GGT AA
GGC CAC ATT TTT GCA GCC TAA
Exons 2, 3 and intron 2; boundaries
3
GCA GCT ATT GAT GTG CGT CC
TGA AGG GAG ACA GCT CTT GTG
Exons 4, 5, 6 and introns 4, 5; boundaries
4
ACT CAG CCA ACC TGA TTT CCT AC
TTC AGA GCC CAG GAG TCA GT
Exons 7, 8 and intron 7; boundaries
2.3. Statistical Analysis
Statistical analysis was performed with Statistical Analysis System V8 (SAS V8). Difference of age was assessed applying t-test or t′-test. Hardy-Weinberg equilibrium (HWE) was calculated by Chi-square analysis. Chi-square or Fisher’s exact test was used to test the differences in genotype and gender between groups. Odds ratios (ORs) and 95% confidence intervals (95% CI) were evaluated by Mantel-Haenszel Chi-squared test to verify the association between variants and FPD or FET. The evaluation of the association was also conducted using logistic regression under different genetic models adjusted for age and gender. Online SHEsis program was used to conduct haplotype analysis [19]. Two-tailed p value < 0.05 was considered significant. The statistical power was performed using Quanto.
3. Results
The patients and controls in the study are well matched for mean age (p=0.12 for FET and p=0.86 for FPD) and sex distribution (p=0.52 for FET and p=0.52 for FPD) (Table 1). By sequencing all the four products in all 221 patients (FET, FPD, and ET-PD) and 100 controls, no exonic variant was identified, while one exon-intron boundary variant (rs2241028) and one intron variant (rs2241027) were detected. In NCBI SNP database, MAF of rs2241027 and rs2241028 were 0.05/0.10, 0.06/0.07, respectively, from ExAC/1000 Genomes Project, both with no clinical significance. The function of both variants was defined as uncertain by MyGenostics. The variants distribution was within the range of Hardy-Weinberg equilibrium in controls (p=0.82, 0.71 resp., Table 3). Given the present sample sizes, we have 80% power to detect an odds ratio of 1.83 in both PD and ET for rs2241027 adopting an additive model and OR of 1.91 in both PD and ET for rs2241028 adopting an additive model. What is worth noting is that there are big differences in MAFs between our control and database in both two variants, which may be caused by ethnical diversity, so we calculate the power considering MAFs of 0.28 and 0.21, respectively, in our control, which is higher than in database; otherwise, it would require much bigger sample sizes. Additionally, we only have 34% power to detect an OR of 1.44 (the OR in Krüger et al.’s study) for rs2241028.
Statistics of minor allele frequency (MAF).
RS number
Position
Function
MAF
ExAC/1000
HWE (pa)
OR (95% CI), p value
FET
FPD
ET-PD
Control
Genomes MAF
FET versus control
FPD versus control
rs2241027
Intron
Uncertain
0.33
0.29
0.33
0.28
0.05/0.10
0.82
1.28 (0.83–1.96), 0.26
0.80 (0.47–1.35), 0.40
rs2241028
Exon-intron boundary
Uncertain
0.14
0.16
0.23
0.21
0.06/0.07
0.71
1.08 (0.70–1.66), 0.73
0.70 (0.42–1.16), 0.17
aHWE for controls.
As for allele and genotype distribution of both variants, we failed to detect any significant differences either in FET versus controls or in FPD versus controls (Tables 3 and 4). No significant difference was observed in the logistic regression either (data not shown). Moreover, haplotypes of two variants showed hardly any association with the risk of FET or FPD (Table 5). Regarding ET-PD, in which we attempted to investigate the situation of HTRA2 in case there were some dramatic mutations, owing to the limitation of sample size, we quitted further statistical analysis.
Statistics of genotype.
Participants
Genotype (rs2241027/rs2241028)
pa value
GG
GA
AA
rs2241027
rs2241028
FET
44/75
48/23
9/3
0.50
0.29
FPD
53/73
43/31
9/1
0.77
0.14
ET-PD
8/9
4/5
3/1
N/A
N/A
Control
51/64
43/30
6/6
N/A
N/A
pa value compared with controls.
Haplotype analysis.
Haplotype
FET (%)
FPD (%)
Control (%)
χ2 valuea/b
Fisher’s pa/b
Pearson’s pa/b
OR (95% CI)a/b
A-A
0
0
0
—
—
—
—
A-G
33
29
28
1.28/0.12
0.26/0.73
0.26/0.73
1.28 (0.83–1.96)/1.08 (0.70–1.66)
G-A
14
16
21
3.05/1.92
0.08/0.17
0.08/0.17
0.63 (0.38–1.06)/0.70 (0.42–1.16)
G-G
53
55
52
0.09/0.58
0.77/0.45
0.77/0.45
1.06 (0.72–1.57)/1.16 (0.79–1.71)
a/bvalue for FET versus controls/FPD versus controls.
4. Discussion
The high temperature requirement A2 (HTRA2), known as a mitochondria protein, plays distinct different roles in mitochondria homeostasis and cellular apoptosis regulation [20]. As one study indicated, deficiency of HTRA2 can cause damage and mutation of mitochondria DNA [21]. Another study revealed that HTRA2 was regulated by PINK1, which might contribute to early-onset PD, in the proteolytic activity [22].
Many researches concerning the association of PD with HTRA2 variants have been done. The earliest mutation screening of HTRA2 in PD patients was done in a German population after the finding that targeted disruption of HTRA2 can cause neurodegeneration and a Parkinsonian phenotype in mice, which resulted in the identification of two mutations (G399S and A141S) related to the risk of PD [23, 24]. Later on, replications with contradictory consequences have been conducted [25–31], and one large scale genetic association study is worth noting, which showed no evidence for an overall association of common variants in HTRA2 with PD [32], while Gulsuner et al.’s study of a six-generation family provides further evidence for the probability of HTRA2 acting as a cause for PD and ET, especially those with family history [13]. So the aim of our study is to investigate the situation of HTRA2 by Sanger sequencing of the whole coding sequence in FET, FPD, and ET-PD in Chinese population, especially FET and ET-PD.
Our study detected two variants (rs2241028, rs2241027). Variant rs2241028 has been reported in several studies with similar negative results except for Krüger et al.’s study, in which rs2241028 was considered as a susceptible factor for PD in the Scandinavian population and their descent from USA [32], while there is no report of this variant in studies about ET. Since rs2241028 is near the splicing region, it may affect the transcript efficiency of HTRA2 to some extent or influence the expression of HTRA2 in some other way, so it would be promising to do some research into the function of this variant and the relationship with PD. Variant rs2241027 has never been mentioned in the previous study no matter about PD or ET. Our study showed that neither of two variants was related to the risk of developing ET or PD, and two variants were defined as no clinical significance in database. Meanwhile, we have not detected mutations (G399S and A141S) mentioned in other studies. So we provided no evidence of association of HTRA2 with FET and FPD. As for ET-PD, the result of our study was not so convincing due to the sample size though we found nothing significant as well. Admittedly, there are some limitations in our study. On the one hand, the sample sizes were only able to detect a moderate correlation with enough power and not for a relatively weaker correlation, which may cause false negative error. On the other hand, it would be more persuasive if the promoter of HTRA2 gene has been sequenced as well.
In conclusion, HTRA2 might not be a cause of familial ET or PD in China. Studies with larger sample size are needed to investigate thoroughly the role of HTRA2 in ET and ET-PD in China and other places in the world.
Competing Interests
The authors report no competing interests.
Authors’ Contributions
Ya-Chao He and Pei Huang contributed equally to this work as first authors.
Acknowledgments
The authors are grateful to all the participants recruited in this study for their donation of blood samples without which the research could not be conducted. The study was supported by National Natural Science Fund (81430022, 91332107, and 81371407). The Sanger sequencing was performed by MyGenostics and Biosune.
LouisE. D.FerreiraJ. J.How common is the most common adult movement disorder? Update on the worldwide prevalence of essential tremor201025553454110.1002/mds.228382-s2.0-77951230101SchmouthJ.-F.DionP. A.RouleauG. A.Genetics of essential tremor: from phenotype to genes, insights from both human and mouse studies2014119-12011910.1016/j.pneurobio.2014.05.0012-s2.0-84905400654ThenganattM. A.JankovicJ.The relationship between essential tremor and Parkinson's disease201622, supplement 1S162S16510.1016/j.parkreldis.2015.09.0322-s2.0-84947703838Benito-LeónJ.LouisE. D.Bermejo-ParejaF.Risk of incident Parkinson's disease and parkinsonism in essential tremor: A Population Based Study200980442342510.1136/jnnp.2008.1472232-s2.0-64749083422KalinderiK.BostantjopoulouS.FidaniL.The genetic background of Parkinson's disease: current progress and future prospects2016134531432610.1111/ane.125632-s2.0-84958665773FunayamaM.OheK.AmoT.FuruyaN.YamaguchiJ.SaikiS.LiY.OgakiK.AndoM.YoshinoH.TomiyamaH.NishiokaK.HasegawaK.SaikiH.SatakeW.MogushiK.SasakiR.KokuboY.KuzuharaS.TodaT.MizunoY.UchiyamaY.OhnoK.HattoriN.CHCHD2 mutations in autosomal dominant late-onset Parkinson's disease: a genome-wide linkage and sequencing study201514327428210.1016/s1474-4422(14)70266-22-s2.0-84923226964McDonnellS. K.SchaidD. J.ElbazA.StrainK. J.BowerJ. H.AhlskogJ. E.MaraganoreD. M.RoccaW. A.Complex segregation analysis of Parkinson's disease: the Mayo Clinic Family Study200659578879510.1002/ana.208442-s2.0-33646384570WhaleyN. R.PutzkeJ. D.BabaY.WszolekZ. K.UittiR. J.Essential tremor: phenotypic expression in a clinical cohort200713633333910.1016/j.parkreldis.2006.12.0042-s2.0-34547683269WiderC.RossO. A.WszolekZ. K.Genetics of Parkinson disease and essential tremor201023438839310.1097/WCO.0b013e32833b1f4c2-s2.0-77954622198ChaoY. X.NgE. Y.TanL.PrakashK. M.AuW.-L.ZhaoY.TanE.-K.Lrrk2 R1628P variant is a risk factor for essential tremor20155, article 902910.1038/srep090292-s2.0-84924565514YuS.-W.ChenC.-M.ChenY.-C.ChangC. W.ChangH.-S.LyuR.-K.RoL.-S.WuY.-R.SLC1A2 variant is associated with essential tremor in taiwanese population201388, article e7191910.1371/journal.pone.00719192-s2.0-84881341410García-MartínE.MartínezC.Alonso-NavarroH.Benito-LeónJ.Lorenzo-BetancorO.PastorP.Ortega-CuberoS.López-AlburquerqueT.AgúndezJ. A. G.Jiménez-JiménezF. J.No association of the SLC1A2 rs3794087 allele with risk for essential tremor in the Spanish population2013231158759010.1097/FPC.0b013e328364db9d2-s2.0-84886599936GulsunerH. U.GulsunerS.MercanF. N.OnatO. E.WalshT.ShahinH.LeeM. K.DoguO.KansuT.TopalogluH.ElibolB.AkbostanciC.KingM.-C.OzcelikT.TekinayA. B.Mitochondrial serine protease HTRA2 p.G399S in a kindred with essential tremor and Parkinson disease201411151182851829010.1073/pnas.14195811112-s2.0-84919898161TzoulisC.ZayatsT.KnappskogP. M.MüllerB.LarsenJ. P.TysnesO.-B.BindoffL. A.JohanssonS.HaugarvollK.HTRA2 p.G399S in Parkinson disease, essential tremor, and tremulous cervical dystonia201511218E226810.1073/pnas.15031051122-s2.0-84928974641ChaoY. X.NgE. Y.FooJ. N.LiuJ.ZhaoY.TanE.-K.Mitochondrial serine protease HTRA2 gene mutation in Asians with coexistent essential tremor and Parkinson disease201516324124210.1007/s10048-015-0443-32-s2.0-84931568294HopfnerF.MüllerS. H.LorenzD.AppenzellerS.KlebeS.DeuschlG.KuhlenbäumerG.Mutations in HTRA2 are not a common cause of familial classic ET20153081149115010.1002/mds.262522-s2.0-84929439026PostumaR. B.BergD.SternM.PoeweW.OlanowC. W.OertelW.ObesoJ.MarekK.LitvanI.LangA. E.HallidayG.GoetzC. G.GasserT.DuboisB.ChanP.BloemB. R.AdlerC. H.DeuschlG.MDS clinical diagnostic criteria for Parkinson's disease201530121591160110.1002/mds.264242-s2.0-84944452191DeuschlG.BainP.BrinM.Consensus statement of the Movement Disorder Society on Tremor. Ad Hoc Scientific Committee199813supplement 3223LiZ.ZhangZ.HeZ.TangW.LiT.ZengZ.HeL.ShiY.A partition-ligation-combination-subdivision em algorithm for haplotype inference with multiallelic markers: update of the SHEsis (http://analysis.bio- x.cn)200919451952310.1038/cr.2009.332-s2.0-64149129845Vande WalleL.LamkanfiM.VandenabeeleP.The mitochondrial serine protease HtrA2/Omi: an overview200815345346010.1038/sj.cdd.44022912-s2.0-39449115415GooH.-G.JungM. K.HanS. S.RhimH.KangS.HtrA2/Omi deficiency causes damage and mutation of mitochondrial DNA2013183381866187510.1016/j.bbamcr.2013.03.0162-s2.0-84877796195Plun-FavreauH.KlupschK.MoisoiN.GandhiS.KjaerS.FrithD.HarveyK.DeasE.HarveyR. J.McDonaldN.WoodN. W.MartinsM. L.DownwardJ.The mitochondrial protease HtrA2 is regulated by Parkinson's disease-associated kinase PINK120079111243125210.1038/ncb16442-s2.0-35748935851StraussK. M.MartinsL. M.Plun-FavreauH.MarxF. P.KautzmannS.BergD.GasserT.WszolekZ.MüllerT.BornemannA.WolburgH.DownwardJ.RiessO.SchulzJ. B.KrügerR.Loss of function mutations in the gene encoding Omi/HtrA2 in Parkinson's disease200514152099211110.1093/hmg/ddi2152-s2.0-25444498785JonesJ. M.DattaP.SrinivasulaS. M.JiW.GuptaS.ZhangZ.DaviesE.HajnóczkyG.SaundersT. L.Van KeurenM. L.Fernandes-AlnemriT.MeislerM. H.AlnemriE. S.Loss of Omi mitochondrial protease activity causes the neuromuscular disorder of mnd2 mutant mice2003425695972172710.1038/nature020522-s2.0-0142246441BogaertsV.NuytemansK.ReumersJ.PalsP.EngelborghsS.PickutB.CorsmitE.PeetersK.SchymkowitzJ.De DeynP. P.CrasP.RousseauF.TheunsJ.Van BroeckhovenC.Genetic variability in the mitochondrial serine protease HTRA2 contributes to risk for Parkinson disease200829683284010.1002/humu.207132-s2.0-44849115171RossO. A.SotoA. I.Vilariño-GüellC.HeckmanM. G.DiehlN. N.HulihanM. M.AaslyJ. O.SandoS.GibsonJ. M.LynchT.Krygowska-WajsA.OpalaG.BarcikowskaM.CzyzewskiK.UittiR. J.WszolekZ. K.FarrerM. J.Genetic variation of Omi/HtrA2 and Parkinson's disease200814753954310.1016/j.parkreldis.2008.08.0032-s2.0-53749093969Simón-SánchezJ.SingletonA. B.Sequencing analysis of OMI/HTRA2 shows previously reported pathogenic mutations in neurologically normal controls200817131988199310.1093/hmg/ddn0962-s2.0-45749135862LinC.-H.ChenM.-L.ChenG. S.TaiC.-H.WuR.-M.Novel variant Pro143Ala in HTRA2 contributes to Parkinson's disease by inducing hyperphosphorylation of HTRA2 protein in mitochondria2011130681782710.1007/s00439-011-1041-62-s2.0-82555192473WangC.-Y.XuQ.WengL.ZhangQ.ZhangH.-N.GuoJ.-F.TanL.-M.TangJ.-G.YanX.-X.TangB.-S.Genetic variations of Omi/HTRA2 in Chinese patients with Parkinson's disease2011138529329710.1016/j.brainres.2011.02.0372-s2.0-79953191414TianJ.-Y.GuoJ.-F.WangL.SunQ.-Y.YaoL.-Y.LuoL.-Z.ShiC.-H.HuY.-C.YanX.-X.TangB.-S.Mutation analysis of LRRK2, SCNA, UCHL1, HtrA2 and GIGYF2 genes in Chinese patients with autosomal dorminant Parkinson's disease2012516220721110.1016/j.neulet.2012.03.0862-s2.0-84860470675ChenC.-M.WuC.-H.HsiehC.-H.LinC.-H.ChenI.-C.ChenY.-C.LeeL.-C.LeeC.-M.TsengY.-C.Lee-ChenG.-J.WuY.-R.HTRA2 variations in Taiwanese Parkinson's disease2014121549149810.1007/s00702-013-1131-92-s2.0-84901262991KrügerR.SharmaM.RiessO.GasserT.Van BroeckhovenC.TheunsJ.AaslyJ.AnnesiG.BentivoglioA. R.BriceA.DjarmatiA.ElbazA.FarrerM.FerrareseC.GibsonJ. M.HadjigeorgiouG. M.HattoriN.IoannidisJ. P. A.Jasinska-MygaB.KleinC.LambertJ.-C.LesageS.LinJ.-J.LynchT.MellickG. D.de NigrisF.OpalaG.PrigioneA.QuattroneA.RossO. A.SatakeW.SilburnP. A.TanE. K.TodaT.TomiyamaH.WirdefeldtK.WszolekZ.XiromerisiouG.MaraganoreD. M.A large-scale genetic association study to evaluate the contribution of Omi/HtrA2 (PARK13) to Parkinson's disease2011323548.e9548.e1810.1016/j.neurobiolaging.2009.11.0212-s2.0-79952900158