Article Text

Short report
Comparative whole transcriptome analysis of Parkinson’s disease focusing on the efficacy of zonisamide
  1. Tatsuhiko Naito1,2,
  2. Wataru Satake1,3,
  3. Pei-Chieng Cha3,4,
  4. Kazuhiro Kobayashi3,
  5. Miho Murata5,
  6. Tatsushi Toda1,3
  1. 1 Division of Neurology/Molecular Brain Science, Kobe University Graduate School of Medicine, Kobe, Japan
  2. 2 Department of Neurology, Graduate School of Medicine, The University of Tokyo, Tokyo, Japan
  3. 3 Department of Statistical Genetics, Osaka University Graduate School of Medicine, Osaka, Japan
  4. 4 Department of Genomic Medicine, Research Institute, National Cerebral and Cardiovascular Center, Suita, Japan
  5. 5 Department of Neurology, NCNP Hospital, Kodaira, Japan
  1. Correspondence to Dr Wataru Satake, Neurology, The University of Tokyo, Bunkyo-ku, Tokyo, Japan; satake{at}; Professor Tatsushi Toda; toda{at}


Objective Interindividual variations in responsiveness to zonisamide in patients with Parkinson’s disease (PD) have been observed in clinical settings. To decipher the molecular mechanisms determining the efficacy of zonisamide, we conducted whole transcriptome sequencing analysis of patients with PD.

Methods We selected 23 super-responders (SRs) and 25 non-responders (NRs) to zonisamide from patients with PD who had participated in a previous clinical trial for the approval of zonisamide for the treatment of ‘wearing-off’. Whole transcriptome analysis of peripheral blood was conducted on samples taken before and 12 weeks after zonisamide treatment. We performed differential gene expression analysis to compare between the SRs and NRs at each time point.

Results Differentially expressed genes in the pre-treatment samples were significantly enriched for glutamatergic synapses and insulin-like growth factor binding (P adj=7.8 × 10−3 and 0.029, respectively). The gene sets associated with these functions changed more dynamically by treatment in SRs than NRs (p=7.2 × 10−3 and 8.2 × 10−3, respectively).

Conclusions Our results suggest that the efficacy of zonisamide in PD patients is associated with glutamate-related synaptic modulation and p53-mediated dopaminergic neural loss. Their transcriptomic differences could be captured before treatment, which would lead to the realisation of future personalised treatment.

  • parkinson's disease
  • Neurogenetics

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See:

Statistics from


Patients with advanced Parkinson’s disease (PD) often experience shortening of beneficial effects of levodopa, termed as ‘wearing-off’. Zonisamide, which was originally developed as an antiepileptic drug, is an effective adjunct for PD patients that improves various symptoms, especially ‘wearing-off’.1 Interindividual variations in responsiveness to zonisamide in PD patients are often observed in clinical settings,1 and we previously conducted a genome-wide association study (GWAS) on PD patients to decipher genetic factors influencing their responsiveness to zonisamide. We identified a single-nucleotide variant (SNP) significantly associated with the effect of zonisamide to reduce ‘off’ time in patients with PD, which could alter the expression of MDM4 and downregulate p53 in neural cells.2

Despite the successful identification of the genetic basis, the mechanisms behind the interindividual variation in zonisamide responsiveness have not been fully elucidated. In addition, considering that a genetic testing is time and cost-intensive, a serological biomarker would be helpful in personalised treatment. To provide insight to these issues, we conducted a comparative whole transcriptome sequencing analysis of the peripheral blood samples of patients with PD, which were taken during a clinical trial for the approval of zonisamide in the treatment of PD.1

Materials and methods

Subjects and definition of zonisamide responsiveness

Participants were selected from 220 Japanese patients with PD who participated in a previous clinical trial for the approval of zonisamide as an adjunct to levodopa therapy for ‘wearing-off’.1 They received either 25 mg or 50 mg of zonisamide per day, and the plasma concentration of zonisamide was measured at the week 4 of their treatment. Blood samples for transcriptome sequencing were collected before and 12 weeks after the start of treatment. Responsiveness to zonisamide was evaluated by the changes in ‘off’ time and the total score of Unified Parkinson’s Disease Rating Scale (UPDRS) part III, from baseline (before treatment) to 12 weeks after the start of zonisamide treatment. A subject who showed a decrease of at least 1.5 hour from baseline in ‘off’ time as well as a decrease of at least five points from baseline in the UPDRS part III total score was defined as a super-responder (SR). On the other hand, a subject who showed no improvement or worsening in both ‘off’ time (change ≥0) and UPDRS part III total score (change ≥0) was defined as a non-responder (NR). Only those who met the criteria of SR or NR were considered for subsequent RNA sequencing analysis.

RNA sequencing

RNA was extracted from blood samples using the RNeasy Mini Kit (QIAGEN) following the manufacturer’s protocol. The quality and quantity of RNA samples were examined using an Agilent Technologies 2200 TapeStation. RNA Integrity Number equivalent (RINe) and rRNA ratio were determined and used to evaluate the quality of the samples. The sequencing library was prepared using the Library Protocol TruSeq RNA Sample Preparation v2 Guide (Illumina). RNA sequencing was performed using the Illumina platform to generate 101 bp paired-end sequences. Sequencing data were converted into raw data for subsequent analysis.

Genotyping, quality control and read counting

Gene expression data preparation was conducted using the following steps: (1) Alignment: the paired-end sequences were mapped to the human reference genome (hg38) by using STAR aligner3 with a two-pass method. PCR duplicates were removed using Picard ( (2) Read counting: the number of reads was obtained at the gene level mapped to Ensembl 70 annotations using RSEM.4

Differential gene expression analysis

Differential gene expression analysis and other statistical analyses were conducted using the Bioconductor package DESeq25 and R (3.6.1) (, respectively. In differential gene expression analysis, we fitted a generalised linear model of negative binomial distribution with a logarithmic link and conducted Wald’s test with age, sex and zonisamide dose as covariates. We set a study-wide significance threshold of p=1.9 × 10−6 based on Bonferroni correction of the total number of tested genes (=0.05/26,475). We analysed the gene expression differences between SRs and NRs before and after treatment. To test whether the expression levels of a gene set changed more dynamically in SRs than in NRs, we calculated the means of absolute logarithm of fold changes by treatment of individual genes, and evaluated the intergroup difference using the Wilcoxon signed-rank test.

Gene set enrichment analysis

Gene set enrichment analysis was performed for differentially expressed genes (DEGs) for gene ontology (GO)6 and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway,7 respectively, using g:Profiler with the ordered query option.8 The order was based on the statistical significance of DEGs.


Clinical characteristics

Among clinical trial participants with complete records, only 25 and 26 participants met the criteria for SR and NR, respectively. Additionally, of these 51 participants, three failed at RNA sample collection or RNA quality tests. Hence, a total of 48 samples (23 SRs and 25 NRs), were processed for subsequent analyses.

Demographic and clinical characteristics of the patients are summarised in table 1. In general, SRs and NRs have matched demographic and clinical information. SRs achieved a mean decrease of 2.8 hours in ‘off’ time whereas NRs showed a mean increase of 1.5 hour in ‘off’ time, after 12 weeks of zonisamide treatment.

Table 1

Clinical characteristics of super-responders and non-responders to zonisamide in this study

RNA-sequencing data profile

The 96 RNA samples from 48 subjects were successfully sequenced and passed quality control. The summary statistics for raw RNA sequencing data are shown in online supplemental table 1). We further removed one sample from the pre-treatment NR data that was strongly suspected of viral infection from the transcriptomic profile (online supplemental note1 and online supplemental figure 1).

Supplemental material

Supplemental material

Differential gene expression analysis

There were no individual DEGs between SRs and NRs that satisfied the study-wide significance threshold both before and after treatment (p>1.9 × 10−6; Figure 1A,B online supplemental data 1). However, we further performed enrichment analysis for genes that satisfied nominal significance (p<0.05) to capture the combined effects of weak but meaningful signals. We identified several significant enrichments at each time point (figure 1C,D). The 613 DEGs in the pre-treatment comparison were significantly enriched for GO terms of glutamatergic synapse, insulin-like growth factor (IGF) binding, and corticosteroid signalling pathway (P adj=7.8 × 10−3, 0.029, and 0.026, respectively; figure 1C). The 595 DEGs in the post-treatment comparison were significantly enriched for a GO term of positive regulation of T cell migration (P adj=0.041; figure 1D). There was no significant enrichment in the KEGG pathway.

Supplemental material

Figure 1

The results of differential gene expression analysis and GO term enrichment analysis of DEGs. (A, B) Volcano plots showing the fold change of transcripts between SRs and NRs (horizontal axis) and their p values (vertical axis) in pre-treatment (A) and post-treatment samples (B). Dashed grey lines represent the nominal significance threshold (p=0.05). There were no DEGs that satisfied the study-wide significance threshold (p>1.9×10−6). (C, D) The results of GO term enrichment analysis of DEGs in pre-treatment (C) and post-treatment (D) samples. Each bar represents the adjusted p values of each term along with the horizontal axis. The vertical dashed line indicates statistical significance threshold (Padj=0.05). BP, biological process; CC, cellular component; DEGs, differentially expressed genes; FC, fold change; MF, molecular function; NR, non-responder; SR, super-responder.

Among the ontologies enriched for DEGs in the pre-treatment, we focused on genes related to glutamatergic synapses because the glutamatergic-mediated pathway was implicated in the mechanism of zonisamide on PD,9 and our previous GWAS suggested its relevance to the responsiveness of zonisamide.2 We then revealed that genes belonging to glutamatergic synapses changed more dynamically in SRs than in NRs (p=7.2 × 10−3). Similarly, genes belonging to IGF binding changed more dynamically in SRs than in NRs (p=8.2 × 10−3), in contrast to the corticosteroid signalling pathway, which resulted in no significance changes (p=0.50).


In the current study, significant enrichment of DEGs between SRs and NRs to zonisamide in PD patients was identified for several biological functions, providing meaningful insights. First, glutamatergic synapses-associated genes were differentially expressed between SRs and NRs before treatment. Zonisamide could increase the release of metabotropic glutamate receptor agonists from astrocytes and inhibit neurotransmission, which may explain the efficacy and tolerability of zonisamide in PD.9 10 In the current study, zonisamide treatment resulted in increased expression levels of glutamatergic synapse-related genes in SRs and reduced differences in expression levels against NRs. This suggests that this pathway might be dysregulated in SRs and that zonisamide is more effective. Since the current study was conducted on blood samples, validation using brain samples is warranted. Second, genes related to IGF binding were also differentially expressed between SRs and NRs. Our GWAS suggested that an SNP located in MDM4 could determine the efficacy of zonisamide for ‘wearing-off’ through suppression of p53.2 Upregulation of p53 is reportedly associated with dopaminergic neuronal apoptosis in PD.11 Several studies supported that zonisamide exerts neuroprotective effects by inhibiting caspase-3, which is regulated by p53,12 against endoplasmic reticulum stress.13 Thus, the caspase-3-mediated effect of zonisamide may be reinforced by suppression of p53 via the overexpression of MDM4 in responders. Since p53 also regulates the IGF-1 pathway,14 the current results suggest that differential expression of genes involved in this pathway may also be present in the blood before treatment. Finally, considering the implications of the immune system in the aetiology of PD15 and the hypothesised antineuroinflammation effect of zonisamide,10 the enrichment in immune-related functions may also be significant.

Overall, our results indicate that the efficacy of zonisamide in PD patients is associated with glutamate-related synaptic modulation and p53-mediated dopaminergic neural loss. Their transcriptomic differences could be captured in blood sample before treatment; thus, more thorough investigations with larger samples would lead to the development of a molecular biomarker that predicts the responsiveness to zonisamide.

Ethics statements

Patient consent for publication

Ethics approval

This study was approved by the Research Ethics Committee of Medicine of the University of Tokyo with the approval number 2021124G. Participants gave informed consent to participate in the study before taking part.


The authors appreciate all investigators included in The Japan Zonisamide on PD Study Group who participated in the clinical study.


Supplementary materials


  • Contributors TN, WS and P-CC analysed the data and wrote the manuscript. WS and KK conducted RNA sequencing. MM and TT conceived the original idea.

  • Funding This study was funded and supported by Sumitomo Dainippon Pharma, of which the roles were the collection of DNA samples in the previous clinical trial and the approval of the study design; AMED, undergrant numbers JP17km0405206 to TT and 17ek0109207h0001 to WS; and JSPS under grant number 21H02836 to WS and 20J12189 to TN.

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

  • Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.

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.