Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2019 Jul;18(1):373-383.
doi: 10.3892/etm.2019.7538. Epub 2019 Apr 30.

RNA-seq analysis of blood of valproic acid-responsive and non-responsive pediatric patients with epilepsy

Affiliations

RNA-seq analysis of blood of valproic acid-responsive and non-responsive pediatric patients with epilepsy

Yan Wang et al. Exp Ther Med. 2019 Jul.

Abstract

Epilepsy is the most common chronic neurological disorder, affecting ~70 million individuals worldwide. However, approximately one-third of the patients are refractory to epilepsy medication. Of note, 100% of patients with genetic epilepsy who are resistant to the traditional drug, valproic acid (VPA), are also refractory to the other anti-epileptic drugs. The aim of the present study was to compare the transcriptomes in VPA responders and non-responders, to explore the mechanism of action of VPA and identify possible biomarkers to predict VPA resistance. Thus, RNA-seq was employed for transcriptomic analysis, differentially expressed genes (DEGs) were analyzed using Cuffdiff software and the DAVID database was used to infer the functions of the DEGs. A protein-protein interaction network was obtained using STRING and visualized with Cytoscape. A total of 389 DEGs between VPA-responsive and non-responsive pediatric patients were identified. Of these genes, 227 were upregulated and 162 were downregulated. The upregulated DEGs were largely associated with cytokines, chemokines and chemokine receptor-binding factors, whereas the downregulated DEGs were associated with cation channels, iron ion binding proteins, and immunoglobulin E receptors. In the pathway analysis, the toll-like receptor signaling pathway, pathways in cancer, and cytokine-cytokine receptor interaction were mostly enriched by the DEGs. Furthermore, three modules were identified by protein-protein interaction analysis, and the potential hub genes, chemokine (C-C motif) ligand 3 and 4, chemokine (C-X-C motif) ligand 9, tumor necrosis factor-α and interleukin-1β, which are known to be closely associated with epilepsy, were identified. These specific chemokines may participate in processes associated with VPA resistance and may be potential biomarkers for monitoring the efficacy of VPA.

Keywords: RNA-seq; biomarkers; efficacy; valproic acid.

PubMed Disclaimer

Figures

Figure 1.
Figure 1.
Heat map of differentially expressed genes between VPA responsive and non-responsive patients. The samples were clustered according to the VPA response. Blue and red color represent low and high levels of expression, respectively. VPA, valproic acid; Resp., responsive; Resi, resistant.
Figure 2.
Figure 2.
Flow chart for the identification of differentially expressed genes. In total, 1,153 genes exhibited variations with P<0.05, and 389 of these had a fold change of >2. After P-value adjustment, 121 genes exhibited Q<0.05. VPA, valproic acid; DEGs, differentially expressed genes.
Figure 3.
Figure 3.
The 2−∆∆Cq values of CCL3 and FOS in the valproic acid non-responsive (n=11) and responsive patients (n=6). (A) CCL3; (B) FOS. Cq, quantification cycle; CCL3, chemokine (C-C motif) ligand 3.
Figure 4.
Figure 4.
Gene ontology enrichment analysis for the DEGs (Q<0.05). (A-C) Terms enriched by the upregulated genes in the categories (A) biological process, (B) cellular component and (C) molecular function. (D-F) Terms enriched by the downregulated genes in the categories (D) biological process, (E) cellular component and (F) molecular function. IgE, immunoglobulin E; MHC, major histocompatibility complex; DEG, differentially expressed gene.
Figure 5.
Figure 5.
Pathways enriched by the 121 DEGs (Q<0.05). DEG, differentially expressed gene; Ig, immunoglobulin; RIG, retinoic acid-inducible gene.
Figure 6.
Figure 6.
Protein-protein interaction network. Nodes with a degree of connectivity of 3–5, 5–10 and 10–13 were indicated in green, yellow and pink, respectively.
Figure 7.
Figure 7.
The modules of the 121 DEGs (Q<0.05). (A) Module 1; (B) module 2; (C) module 3. Pink and yellow nodes represent upregulated and downregulated genes, respectively.

Similar articles

Cited by

References

    1. Singh A, Trevick S. The epidemiology of global epilepsy. Neurol Clin. 2016;34:837–847. doi: 10.1016/j.ncl.2016.06.015. - DOI - PubMed
    1. Lv RJ, Shao XQ, Cui T, Wang Q. Significance of MDR1 gene C3435T polymorphism in predicting childhood refractory epilepsy. Epilepsy Res. 2017;132:21–28. doi: 10.1016/j.eplepsyres.2017.02.010. - DOI - PubMed
    1. Voll A, Hernández-Ronquillo L, Buckley S, Téllez-Zenteno JF. Predicting drug resistance in adult patients with generalized epilepsy: A case-control study. Epilepsy Behav. 2015;53:126–130. doi: 10.1016/j.yebeh.2015.09.027. - DOI - PubMed
    1. Rakitin A, Kõks S, Haldre S. Valproate modulates glucose metabolism in patients with epilepsy after first exposure. Epilepsia. 2015;56:e172–e175. doi: 10.1111/epi.13114. - DOI - PubMed
    1. Fernando-Dongas MC, Radtke RA, Vanlandingham KE, Husain AM. Characteristics of valproic acid resistant juvenile myoclonic epilepsy. Seizure. 2000;9:385–388. doi: 10.1053/seiz.2000.0432. - DOI - PubMed