Abstract
Blockade of Angiotensin type 1 receptor (AT1R) has potential therapeutic utility in the treatment of numerous detrimental consequences of epileptogenesis, including oxidative stress, neuroinflammation, and blood–brain barrier (BBB) dysfunction. We have recently shown that many of these pathological processes play a critical role in seizure onset and propagation in the Scn8a-N1768D mouse model. Here we investigate the efficacy and potential mechanism(s) of action of candesartan (CND), an FDA-approved angiotensin receptor blocker (ARB) indicated for hypertension, in improving outcomes in this model of pediatric epilepsy. We compared length of lifespan, seizure frequency, and BBB permeability in juvenile (D/D) and adult (D/+) mice treated with CND at times after seizure onset. We performed RNAseq on hippocampal tissue to quantify differences in genome-wide patterns of transcript abundance and inferred beneficial and detrimental effects of canonical pathways identified by enrichment methods in untreated and treated mice. Our results demonstrate that treatment with CND gives rise to increased survival, longer periods of seizure freedom, and diminished BBB permeability. CND treatment also partially reversed or ‘normalized’ disease-induced genome-wide gene expression profiles associated with inhibition of NF-κB, TNFα, IL-6, and TGF-β signaling in juvenile and adult mice. Pathway analyses reveal that efficacy of CND is due to its known dual mechanism of action as both an AT1R antagonist and a PPARγ agonist. The robust effectiveness of CND across ages, sexes and mouse strains is a positive indication for its translation to humans and its suitability of use for clinical trials in children with SCN8A epilepsy.
Introduction
Epilepsy is a common neurological disorder that can have devastating consequences for children [1]. Advances in genomic research have implicated many individually rare genetic variants as causal for a significant proportion of early-onset epilepsy [2]. Pathogenic gain-of-function variants in the voltage-gated sodium channel gene, SCN8A, underlie a severe early onset seizure syndrome that is accompanied by global developmental delays, intellectual disabilities, and a host of comorbidities [3]. These variants alter the function of the NaV1.6 channel, the most abundant sodium channel in the human brain and a master regulator of neuronal excitability [4]. Unfortunately, more than 75% of SCN8A patients treated with currently marketed antiseizure medications (ASM) do not achieve seizure freedom, and many suffer adverse side effects. Indeed, there are no commercially available drugs that specifically target SCN8A or that prevent progression of the disease [5].
We previously investigated a mouse model with a pathogenic Scn8a knockin variant (N1768D) and found that prior to seizure onset there was increased blood–brain barrier (BBB) permeability and genome-wide alterations in transcript abundance, both of which were exacerbated once seizures were established [6]. These patterns of disease progression fit a model of epileptogenesis (i.e., the process by which neuronal networks are altered toward increased seizur e susceptibility and generation of chronic spontaneous seizures), similar to that seen in post-traumatic epilepsy in which there is an initial injury, a latent phase and a chronic epilepsy phase [7]. We suggested that pathological conditions associated with SCN8A-related epilepsy include the same pathogenic ‘tetrad’ of glutamate excitotoxicity, oxidative stress (OXS), inflammation, and BBB dysfunction (BBBD) that characterizes the neurobiology of other brain disorders [6,8].
Angiotensin receptor blockers (ARBs) represent a class of drugs that has shown promise in treating the ‘injury cascade’ associated with brain insults by effectively blocking the physiological effects of angiotensin II, the main active factor of the renin–angiotensin system (RAS) [9]. After brain injury, and in various brain disorders, increased Angiotensin II (Ang II) activates the Angiotensin type 1 receptor (AT1R) leading to glutamate excitotoxicity, OXS, BBB breakdown, brain inflammation, cerebrovascular remodeling, as well as other risk factors for neuronal injury, cognitive decline, post-traumatic epilepsy, and neurodegenerative diseases [10,11]. Blocking AT1R activation with ARBs has therapeutic potential in many chronic disease states [11,12] as they are known to have powerful anti-inflammatory, anti-fibrotic, and antioxidant activities [11,13,14]. In addition, some studies have shown promising effects in various epilepsy models [15–19]. Therefore, blocking AT1R activation with ARBs has therapeutic potential in many chronic disease states [11,12].
In this study, we test the efficacy of candesartan (CND), an FDA-approved ARB indicated for hypertension, for its ability to prevent BBB dysfunction and mitigate epileptogenesis in our Scn8a mouse model, which includes the N1768D variant in heterozygous (D/+) and homozygous (D/D) form on two strain backgrounds (C57BL/6J and C3H/HeJ). In particular, we examine survival and seizure frequency in untreated and treated juveniles, as well as in untreated and treated adults of both sexes and on both strain backgrounds. We carry out BBB permeability studies and RNAseq on hippocampal tissues in order to determine (1) what distinguishes gene expression signatures of untreated and treated juveniles and adults, (2) which pathways are altered in a manner consistent with overstimulation of AT1R, and (3) by what mechanism(s) CND treatment leads to improved outcomes.
Materials and methods
Mouse strains
The C57BL/6J-N1768D congenic mouse (B6-D/+) was used for phenotyping, RNAseq and BBB studies [20]. In addition, a hybrid strain that was produced by backcross breeding to transfer the Scn8a-N1768D allele from the B6 to the C3H/HeJ strain background was used to further assess the effect of CND on survival and seizure frequency. After constructing an F1 by crossing male B6-D/+ heterozygotes to female C3H/HeJ, five rounds of backcrossing (N6) produced a 2.5% B6 and 98.5% C3H/HeJ-N1768D line (C3H-D/+). Sister-brother mating maintained the 98.5% C3H background. This line (referred to here as C3H) was observed to have an altered phenotype in terms of survival of the homozygote (C3H-D/D), which typically only survive for ∼25 days after a series of ∼60–80 tonic-clonic seizures (TCs) beginning at ∼15–20 days of age (Figure 1A). C3H-D/D mice have a delayed seizure onset (∼40–45 days) and an extended survival (65–70 days) with intermittent bouts of TCs (Figure 1A).
Mouse lines and dosing regimen for B6-D/D juvenile and B6-D/+ adult mice
Phenotyping
Female and male mice on both the B6 and C3H backgrounds were housed in sex-specific groups of 3-4 per cage in a pathogen-free mouse facility with a 14 h light/10 h dark cycle (lights turned on at 5am). As described previously [6], a 24/7 video monitoring system was utilized to collect seizure data, with infrared illumination to monitor behavior during the dark period. Seizures were counted as individual TCs as described in [6,20]. Adult mice were followed from the age of 6 weeks, a period of time well before seizure onset, in order to identify the day of seizure onset. We define a seizure bout as a cluster of seizures on consecutive days, and a gap as a seizure-free period of 3 days or more [6]. This research received formal approval from the University of Arizona Institutional Animal Care and Use Committee Program (IACUC #16-160). All experiments were designed in accordance with the Animal Research: Reporting In Vivo Experiments (ARRIVE) guidelines [21] and efforts were made to minimize animal stress and suffering and to reduce the number of mice used. Genotyping at the Scn8a-N1768D site was carried out as previously described [6,20].
Drug administration
Preparation of CND suspension for injections and oral dosing
Juvenile B6-D/D and adult B6-D/+ mice were given subcutaneous (s.c.) injections with either CND (2–4 mg/kg/day; doses that fall within the FDA-approved range for CND in both adults and children) or vehicle (VEH) (Sigma-Aldrich) [22,23] (Figure 1B). Oral dosing of adults was performed by thoroughly mixing 160 μl of suspension containing 1.28 mg of CND or VEH with 5 g of pure (organic) peanut butter in a petri dish [20]. Mice were given 1 g pellet per mouse each day. Pellets were weighed each day to record amount consumed and 24/7 video used to verify that each mouse was consuming pellet. An average of 0.79 g (±0.36 g) was consumed by an adult mouse (∼20–25 g), corresponding to 9.1 ± 4.2 mg/kg/day. No correlations were found for amount of CND consumed and lifetime survival or number of TCs (Supplementary Figure S1).
Survival and seizure frequency studies
B6-D/D mice were injected with CND (n=10) or VEH (n=10) at 15 days of age (P15) and followed to monitor TCs by 24/7 video recording until death (Figure 1B). A separate set of 10 B6-D/D mice was administered a daily s.c. injection of 25 mg/kg/day phenytoin (PHT) suspension or VEH (US Pharmacopeia, Sigma-Aldrich) [24]. Oral administration of CND or VEH was initiated at the time of seizure onset (range P75-P85) for adult B6-D/+ females (N=10) and males (N=10) (Figure 1B). Mice were monitored for seizures by 24/7 video recording until time of death.
Transcriptome studies
B6-D/D and wild-type (+/+) (n=3) were administered CND or VEH by s.c. injection for a period of 5 days starting at P15 and sacrificed for removal of hippocampal tissue (n=12) (Figure 1B). Adult B6-D/+ male (n=3) mice were administered CND or VEH pellets as described above at the time of seizure onset and monitored for seizures for an extended period before sacrifice and removal of hippocampal tissue (100-160 d) (i.e., to match the age at which a VEH-treated male B6-D/+ experience ∼20 TCs) (Figure 1B).
Blood–brain barrier studies
B6-D/+ females (n=6) and males (n=6) were administered CND (2–4 mg/kg/day) or VEH by s.c. injection at the time of seizure onset for a period of 10 days (Figure 1B). On the 11th day these two groups of mice were subjected to BBB analysis as described in the next section.
Blood–brain barrier analysis
RNA sequencing
Hippocampal tissues from 21 mice were obtained according to a sampling strategy that compared (1) untreated (VEH only) and treated B6-D/D (n=3) with age-matched B6-+/+ controls (n = 3 per group), (2) untreated (VEH only) male B6-D/+ mice that had experienced ∼20 TCs (n = 3 per group) with age-matched B6-+/+ controls (n = 3 per group), and (3) male B6-D/+ mice that were treated with CND for an extended period (n = 3 per group). Brains of treated and untreated mice were dissected to yield tissue samples from the hippocampus [28]. Bulk tissue was stored in RNALater (Qiagen, Valencia, CA) at −80 degrees centigrade. The technique for analyzing hippocampal gene expression was performed as previously described [6]. Briefly, RNA was isolated from hippocampal tissue and initial QC performed. Libraries were constructed using a stranded mRNA-Seq Kit and average fragment size was assessed. After concentrations were determined with an adaptor-specific qPCR kit, equimolar samples were pooled and clustered for sequencing on a Novaseq instrument (Illumina). Sample data were demultiplexed, trimmed and quality filtered, and Fastq files were splice aligned against the GRCh37 reference genome using STAR aligner version 2.5.2b [29]. Gene expression counts were obtained using htseq-count version 0.6.1 [30]. Both splice alignment and counting were performed with Ensembl Annotation of the NCBI reference genome and raw counts analyzed with edgeR version 3.16.5 [31].
Data and pathway analysis
Results of statistical summaries were generally expressed as mean ± SD. Kaplan–Meier survival curves were used to test for differences in survival. In cases where groups did not have the same variance, we performed two-sample t-tests. We used a ‘perturbation signature’ approach to identify genome-wide differences in transcript abundance between transgenic mice (i.e., D/+, D/D or +/+) that were untreated (VEH) or treated with CND [32]. Differential expression analysis was performed as previously described [33,34]. Briefly, we utilized the exactTest function in edgeR and gene expression counts were first normalized using the calcNormFactors function. DEG analysis was performed on treatment versus control groups at two different ages using three biological replicates per group, which has been shown in power analyses to be sufficient to yield a true positive rate greater than 80% under the conditions used here [35]. Multidimensional scaling (MDS) plots were constructed using the ‘plotMDS’ function in edgeR, which plots samples on a 2D scatterplot so that distances on the plot approximate the typical log2 fold changes between samples. All significant differentially expressed genes (DEGs) (false discovery rate [FDR] < 0.05) were analyzed with Ingenuity® Pathway Analysis (IPA) to identify biological pathways that were significantly activated or deactivated as compared with controls and to identify putative upstream transcriptional regulators (Qiagen, Hilden, Germany). Rather than focusing on any single gene, bioinformatic analyses of our RNAseq data identified the most statistically significant biological pathways that were enriched given the set of DEGs in each experiment.
Pathway literature searches and feature selection
ChatGTP was used to aid in literature searches on the role of canonical pathways identified by IPA in the context of brain disorders. Initially, the question was posed: ‘what does “canonical pathway name” have to do with brain injury?’ Several pathway ‘effects’ commonly reported in pathological conditions (e.g., neuroinflammation, BBB, etc.) were identified in the untreated and treated B6-D/D juvenile and B6-D/+ male adults. A second round of searches was performed including the term ‘activation’ or ‘deactivation’; for example, ‘what does activation/deactivation of the “canonical pathway” name have to do with brain injury?’ For each canonical pathway, a score of 1 or 2 was assigned to each effect depending on whether it was ‘detrimental’ or ‘beneficial’, respectively, in association with an activated or deactivated state as appropriate. For example, in the case of the pathway effect ‘neuroinflammation’, a ‘1’ or ‘2’ was assigned if pro- or anti-inflammatory cytokines were inferred to be produced early and/or late in the injury process, respectively. A score of 1/2 was assigned if the effect was beneficial in one context and detrimental in another. If an effect was not involved in a given pathway it was assigned ‘0’. An underlying assumption in the interpretation of whether an effect was beneficial or detrimental was that mice were sampled in the later stages of epileptogenesis (i.e., after the establishment of spontaneous seizures) [6]. Finally, a list of cited references in each search was compiled and inspected to verify results for the pathways that were deemed to be most relevant in the study.
Results
Differential survival and seizure frequency for untreated and treated mice
Adults
Supplementary Table S1 displays survival and seizure statistics for female and male untreated (N=17) and treated (n=10) B6-D/+ monitored 24/7 by video beginning at P30 and continuing for the remainder of their entire life span. Data for the untreated mice were previously reported [6]. Supplementary Table S1 also provides C3H-D/D summary statistics for female untreated (N=15) and treated (N=15) and male untreated (N=27) and treated (N=16) monitored until time of death. Age at seizure onset was not significantly different for untreated and treated mice within each line, except in the case of male C3H-D/D (44.0 ± 6.5 days vs 48.1 ± 5.1 days, t-test P-value = 0.018).
In all cases, treated mice lived longer and experienced a lower post-onset seizure frequency. For example, there was a significant increase in survival of B6-D/+ treated versus untreated females (167.6 ± 19.1 vs 132.9 ± 21.4 days, t-test P-value = 1.38 × 10−4) and B6-D/+ treated versus untreated males (135.8 ± 28.4 vs 91.3 ± 15.2 days, t-test P-value = <1.0 × 10−5), reflecting a 34.7% and 48.9% increase, respectively (Figure 2A and Supplementary Table S1). Similarly, C3H-D/D treated versus untreated females (69.9 ± 12.9 vs 117.4 ± 28.0 days, t-test P-value = <1.0 × 10−5) and C3H-D/D treated versus untreated males (67.9 ± 9.9 vs 117.1 ± 20.4 days, t-test P-value = <1.0 × 10−5), reflecting a 47.5% and 49.2% increase, respectively (Supplementary Figure S2 and Table S1).
Lifetime survival, seizure frequency and seizure gaps for adult mice
Kaplan–Meier survival curves are shown for all adults in Supplementary Figure S3 and for C3H-D/D in Figure 2B. Unequal survival of treated versus untreated B6-D/+ females and males was strongly supported in a goodness of fit test using the χ2 distribution (right-tailed) (P-value = 2.1 × 10−3 and 3.5 × 10−4, respectively), which also indicated a medium and large observed standard effect sizes of 0.59 and 0.69, respectively (Supplementary Figure S3A and S3B). C3H-D/D treated females and males showed similar patterns, with unequal survival strongly supported (P-values = <1.0 × 10−5 in both cases) with large observed standard effect sizes of 0.83 (supplementary Figure S3C) and 0.79, respectively (Figure 2B).
Treated adult mice also had a lower post-onset seizure frequency, with B6-D/+ females and males experiencing 88.9% and 59.1% reductions, respectively (t-test P-value = <1.0 × 10−5 and 9.4 × 10−4, respectively) (Figure 2C); and C3H-D/D females and males experiencing 78.4% and 48.0% reductions, respectively (t-test P-value = 4.2 × 10−2 and = 1.2 × 10−4, respectively) (Supplementary Figure S2D and Table S1). These reductions are reflected in a greater number of post-onset seizure gaps (i.e., ≥3 seizure-free days) [6] for B6-D/+ males (t-test P-value = 3.9 × 10−2) and C3H-D/D females (t-test P-value = <1.0 × 10−5) and males (t-test P-value = <1.0 × 10−5) (Supplementary Figure S4 and Table S1). In addition, there was a greater gap length in B6-D/+ females (t-test P-value = 1.6 × 10−4) and males (t-test P-value = 3.6 × 10−5), as well as in C3H-D/D males (t-test P-value = 3.3 × 10−2) (Supplementary Figure S3 and Table S1), resulting in a greater number of ‘gap days’ (shown for B6-D/+ in Figure 2D).
Juveniles
Supplementary Table S2 displays survival and seizure statistics for two groups of untreated and treated B6-D/D juvenile mice: a set treated with CND (N=10) and VEH (N=10) and a second set treated with PHT (N=10) and VEH (N=10). There was a significant increase (27%) in survival of juveniles treated with CND (31.8 ± 2.7 vs 25.5 ± 2.1 days, t-test P-value = <1.0 × 10−5) (Figure 3A; Supplementary Table S2). Juveniles treated with PHT experienced a minor increase in survival that did not reach statistical significance (30.1 ± 5.1 days vs 27.6 ± 3.2) (t-test, P-value = 0.083). Kaplan–Meier survival curves similarly show unequal survival of CND treated versus untreated juveniles (P-value = 4.2 × 10−5) with a large observed standard effect sizes of 0.79 (Figure 3B), while PHT treated juveniles did not show a statistically significant increase in survival (P-value = 0.132) (Figure 3C). For juveniles treated with CND there was a slight decrease in post-onset seizure frequency (15.4 ± 10.3 vs 28.4 ± 17.7 TCs, t-test, P-value = 2.9 × 10−2), and a more robust reduction in post-onset seizure frequency for juveniles treated with PHT (10.2 ± 7.1 vs 24.8 ± 11.0 TCs, t-test, P-value = 7.4 × 10−4) (Figure 3D; Supplementary Table S2).
Lifetime survival and seizure frequency for juveniles treated with CND and PHT
Effect of CND on blood–brain barrier permeability in female and male adults
Our results previously demonstrated that BBB paracellular permeability (i.e., ‘leak’) to [14C]sucrose, a small molecule tracer that does not cross the intact BBB [36], increased in both pre-TC B6-D/+ females and males relative to wild-type controls [6]. The magnitude of the sucrose permeability increase was shown to further increase over wild-type and pre-seizure levels in untreated post-TC B6-D/+ females and males. Treatment with CND before seizure onset prevents BBB paracellular permeability from increasing above physiological levels in both B6-D/+ females and males (data not shown). After seizure onset, treatment with CND results in a statistically significantly reduced BBB paracellular permeability (i.e., ‘leak’) relative to untreated post-seizure females (P<0.05) and untreated post-seizure males (P<0.001) (Figure 4).
BBB paracellular permeability to sucrose is increased in untreated B6-D/+ female and male mice relative to those treated with CND after seizure activity
Effects of CND on hippocampal gene expression in juvenile and adult mice
Figure 5A shows an MDS plot of the 500 most variably expressed genes for the 21 B6-D/D and B6-D/+ individuals submitted to RNAseq analysis. The plot shows untreated juveniles and adults on the left side of the plot and CND-treated juveniles and adults clustering on the right side with age-matched wild-type individuals. The co-clustering of juvenile and adult mice on both sides of the plot suggests that treatment with CND is the chief explanatory factor (i.e., more important than age or genotype) underlying inter-group genome-wide gene expression differences.
Number of differentially expressed genes and pathways identified by RNAseq and pathway enrichment analysis for untreated and treated B6-D/D juvenile and B6-D/+ adult mice
To further investigate the effects of CND treatment on gene expression in mice experiencing TCs, we characterized the number of DEGs that are up-regulated and down-regulated in the untreated and treated groups. For untreated juveniles, 1518 transcripts were identified with an FDR-adjusted P-value <0.05, 956 of which were up-regulated and 562 of which were down-regulated (Figure 5B, upper panel). CND-treated B6-D/D mice had a much smaller number of DEGs: 50 up-regulated and 25 down-regulated. Only 6 of the up-regulated and 2 of the down-regulated transcripts were shared with untreated mice. Untreated B6-D/+ adult males had a total of 602 DEGs (FDR-adjusted P-value <0.05), 454 of which were up-regulated and 148 were down-regulated. CND-treatedf B6-D/+ males had 86 DEGs: 76 up-regulated and 10 down-regulated. Four upregulated transcripts were shared between untreated and treated males (Figure 5B, upper panel).
Canonical pathways altered in untreated and treated juvenile mice
For untreated B6-D/D juveniles, pathway enrichment procedures identified 52 canonical pathways with P-values ≤0.05 and z-scores with absolute values ≥2.0 (Supplementary Table S3), including 10 predicted to be deactivated and 42 predicted to be activated (Figure 5B, lower panel). There were no enriched canonical pathways for CND-treated B6-D/D under strict FDR. After relaxing the FDR cutoff from 0.05 to 0.25 a set of related pathways became statistically significant in treated juveniles: Oxidative Phosphorylation (-log(P-value) = 9.9, z-score = 5.0), EIF2 Signaling (-log(P-value) = 4.2, z-score = 3.2), Mitochondrial Dysfunction (-log(P-value) = 7.6, z-score = -3.8), and Granzyme A Signaling (-log(P-value) = 5.5, z-score = -3.4). The enrichment results reflect a common set of 19 up-regulated electron transport chain (ETC) subunit genes that alternately result in positive and negative activation z-scores for the oxidative phosphorylation (OXPHOS) and mitochondrial dysfunction (MD) pathways, respectively (another 24 non-ETC genes distinguish the MD from OXPHOS pathways). These genes include 12 within ETC complex I (Ndufa2, Ndufa4, Ndufa6, Ndufb4, Ndufb6, Ndufs3, Ndufs4, Ndufs6, Ndufs7, Ndufs8, Ndufv1, and Ndufv2), 3 within ETC complex III (Uqcrc1, Uqcr2 and Uqcr11), and 4 within ETC complex IV (Cox6c, Cox7a2, Cox7b, and Surf1). Only a single transcript out of 96 ETC genes was differentially expressed in untreated juveniles; the up-regulation of UCP2 (+3.32-fold, P-value = 1.39 × 10−8, z-score = 1.06 × 10−5) has been shown to shift metabolism away from OXPHOS toward glycolysis [37]. A Fisher exact test comparing the number of differentially expressed ETC genes (FDR ≤0.25) in treated versus untreated yields a P-value of 1.64 × 10−5, and results in an odds ratio representing a 23-fold enhancement of ETC gene activation in treated juveniles (95%CI = 3.5–979).
Canonical pathways altered in untreated and treated adult mice
For untreated B6-D/+ adults, pathway enrichment procedures identified 79 canonical pathways with P-values ≤0.05 and z-scores with absolute values ≥2.0 (Supplementary Table S3), including 5 predicted to be deactivated and 74 predicted to be activated (Figure 5B, lower panel). Treatment with CND resulted in an enrichment of 14 canonical pathways: 3 predicted to be deactivated and 11 predicted to activated. Only three activated pathways were shared between untreated and treated B6-D/+ males (Supplementary Table S3). Eight pathways uniquely activated in treated adults are involved in T-helper 1 cell-mediated immune responses, 5 of which involve Regulatory T cells (Tregs) (Dendritic Cell Maturation, Calcium-induced T Lymphocyte Apoptosis, ICOS-ICOSL Signaling in T Helper Cells, PKCθ Signaling in T Lymphocytes, T Cell Receptor Signaling).
Altered canonical pathways shared between untreated juvenile and adult mice
Figure 5C shows a bubble chart of the 30 canonical pathways shared between untreated B6-D/D juveniles and B6-D/+ adults. Of the 30 shared enriched pathways, 2 are predicted to be deactivated (PPAR and RHOGDI Signaling) and 28 are predicted to be activated (Figure 5D and Supplementary Table S3). The top enriched pathways (i.e., -log(P-value) ≥7.0) that are shared between juvenile and adults include Pulmonary Fibrosis Idiopathic Signaling, Osteoarthritis Pathway, Hepatic Fibrosis Signaling, S100 Family Signaling, Wound Healing Signaling, Role of JAK family kinases in IL-6-type Cytokine Signaling, and IL-6 Signaling.
Of the pathways uniquely altered in untreated juveniles, 8 were predicted to be deactivated and 11 were predicted to be activated (Figure 5D). Of those uniquely altered in untreated adults, 3 are predicted to be deactivated and 43 predicted to be activated (Figure 5D). None of the abovementioned pathways was found to be activated in juveniles or adult males treated with CND (Supplementary Table S3).
Predicted upstream regulator molecules
To identify potential drivers of the differential expression pattern observed within each dataset we used the upstream regulator function in IPA. Lipopolysaccharide (LPS), TNF, IL1B, and transforming growth factor-β (TGF-β) were the top predicted upstream activators in both untreated juveniles and adult males (including their rank order), with activation z-score ranging from 7.1 to 9.6 and 6.8 to 8.6, respectively (Supplementary Table S4). None of the z-scores was statistically significant for these upstream activators in either juvenile or adults treated with CND (Supplementary Table S4). Similarly, neither of the top predicted upstream inhibitors (also shared between) untreated juvenile and adults – APOE (z-scores: -4.0 and -4.6, respectively) and PPARGC1A (z-scores: -2.4 and -3.8, respectively) – reached statistical significance in treated juveniles and adults (Supplementary Table S4). While not appearing on the list of upstream regulators in juveniles, the top predicted upstream inhibitor for treated adults (PTGS2 or cyclooxygenase-2, z-score = -2.8) was oppositely predicted to be an upstream activator in the case of untreated adults (z-score = 2.5).
Beneficial and detrimental pathway effects
Literature searches querying cellular processes associated with canonical pathways in Supplementary Table S3 resulted in many that are commonly reported in brain injury or neurological disorders. A total of 12 such effects was compiled with tallies of the number of times each effect was involved for all 109 pathways. Many of these pathway effects are considered hallmarks of neurological disease or known to play prominent roles in disease pathogenesis or response to brain insults [38,39]. The most common pathway effect was ‘Neuroinflammation’ (INF), which was involved in 80 of the 109 (73.4%) pathways. This was followed by BBB (51.4%), ‘neurogenesis/repair’ (NER, 51.4%), ‘apoptosis/clearance’ (APO, 48.6%), OXS (36.7%), ‘glial activation’ (GLA, 35.8%), ‘synaptic plasticity’ (SYN, 29.4%), ‘angiogenesis’ (ANG, 27.5%), ‘secondary injury’ (SIC, 23.9%), ‘extracellular matrix’ (ECM, 17.4%), ‘neurodegeneration’ (NDG, 15.6%), and ‘fibrosis/scar formation’ (FBS, 13.8%). Juvenile versus adult effect frequencies differed slightly, with some higher in adults: APO (36.5% vs 55.7%), BBB (50.0% vs 64.6%), SIC (19.2% vs 30.4%), SYN (23.1% vs 29.1%), GLA (34.6% vs 48.1%), respectively, and OXS higher in juveniles (46.2% vs 35.4%), respectively (Figure 6A). The pathway effects most often inferred to be beneficial are ANG (86.7%), NER (85.7%), SYN (78.1%), and APO (67.9%), while the effects most often inferred to be detrimental are NDG (94.1%), FBS (93.3%), SIC (73.1%), GLA (66.7%), and BBB (64.3%).
Detrimental and beneficial pathway effects in untreated and treated B6-D/D juvenile and B6-D/+ adult male mice
The mean number of effects associated with a pathway is 4.2 ± 2.0. Comparing the total number of detrimental and beneficial effects associated with untreated juvenile and adult pathways, we find that juveniles have a slightly higher ratio of detrimental to beneficial effects (62.3%) compared with adults (52.2%) (Figure 6B). The mean number of detrimental effects for juvenile pathways (2.6 ± 2.6) is statistically significantly higher than beneficial effects (1.6 ± 1.7) (t-test, P-value = 9.4 × 10−4), which is not the case for adult pathways (2.5 ± 2.7 vs 2.3 ± 1.3, respectively; t-test, P-value = 0.293). Interestingly, there is a statistically significant shift toward increased detrimental effects (73.8%) among the 33 pathways shared between untreated juvenile and adults (Figure 6B): the mean number of detrimental and beneficial effects is 3.7 ± 2.6 and 1.3 ± 1.8, respectively (t-test, P-value = 3.4 × 10−6). This shift toward detrimental effects is clearest for OXS (+41.7%), APO (+39.4%), NER (+32.4%), and SIC (+26.9%), with only ANG (13.3%) shifting toward beneficial. On the contrary, we find that the mean number of beneficial effects is statistically significantly greater (1.7 ± 1.7) than that for detrimental effects (0.4 ± 0.9) (t-test, P-value = 0.014) for pathways identified in treated mice (Figure 6C). Detrimental tallies were cut by more than half for nearly all of the effects (Figure 6D) other than for INF and NER. In addition, a new pathway effect appeared in treated adults: Activation of T cells (Tregs).
When we tally the number of detrimental and beneficial effects for each pathway with ≥3 effects, we find that 26 are inferred to be associated solely with detrimental effects and 36 solely with beneficial effects (Supplementary Table S3). The RAS and PPAR Signaling pathways are associated with the largest number of detrimental effects, each with 9 detrimental and 0 beneficial effects (Figure 6C), while the PI3K/AKT, HGF and Prolactin Signaling pathways are associated with the largest number of beneficial effects, each with 0 detrimental and 7 beneficial effects (Supplementary Table S3).
Discussion
This is the first study to test CND administration in a mouse model with a natural onset of seizures, and to investigate the cellular and molecular mechanisms involved in improved outcomes. The results indicate significant efficacy of CND from three perspectives: survival and seizure frequency, BBB function, and genome-wide hippocampal gene expression. The findings of extended lifespan, longer seizure-free periods, and reduced seizure frequency in Scn8a-N1768D mice treated with CND are robust to age, sex, and strain background (Figures 2 and 3). We previously tested cannabidiol (CBD) on B6-D/+ and B6-D/D mice and found that males (but not females) benefited with increased survival and reduced seizure frequencies: B6-D/+ adult and B6-D/D juvenile males survived for a mean of 17.4 (18.2%, t-test, P-value = 4.5 × 10−3) 2.5 (10.2%, t-test, P-value = 0.032) days longer, while each experienced reduced seizure frequencies of 31.3% and 44.0%, respectively [20]. Treatment with CND yielded greater increases in survival, with B6-D/+ males surviving for a mean of 44.5 days longer (48.9%, t-test p-value = <1.0 × 10−5) and B6-D/D males surviving for 6.3 days longer (27%, t-test, P-value = <1.0 × 10−5). These increases represent larger standard effect sizes of 0.69 and 0.79, respectively, compared with 0.50 and 0.35 for CBD-treated B6-D/+ and B6-D/D males [20]. In addition, CND treatment resulted in greater reductions in seizure frequency for B6-D/+ males (56.5%) and B6-D/D juveniles (45.8%). The present study found CND also outperformed phenytoin – an antiseizure medication shown to be efficacious for patients with SCN8A [40], in terms of survival in juvenile Scn8a-D/D mice. While phenytoin showed a greater reduction in post-onset seizure frequency, CND provided significant improvements in both survival and seizure control. The same mouse model was used to test chronic administration with Prax-330 (GS967) – a novel sodium channel modulator [41]. While we couldn’t make direct comparisons with testing on B6-D/+ adults due to differences in study design, Prax-330 treatment of B6-D/D mice yielded an increased survival of 5.7 days (29.8%), representing a slightly smaller standard effect size (0.69) as we found for CND.
We also found that adult B6-D/+ females and males both exhibit significantly reduced BBB ‘leak’ when treated with CND after the establishment of seizures (Figure 4). Untreated B6-D/D juvenile and B6-D/+ adult males exhibited significant ‘disease-induced’ genome-wide changes in hippocampal transcript abundance, many of which were returned to physiological levels in treated mice. Characterization of enriched pathways in untreated mice revealed a host of cellular and molecular processes representing both adaptive and maladaptive responses to increased neuronal excitability (Supplementary Table S3).
In the following sections we discuss physiological effects of Scn8a-induced neuronal hyperexctability in the context of the juvenile and adult brain, and then ask whether the known mechanism(s) of action of ARBs can explain our findings of improved outcomes in CND-treated mice. We used a ‘two-prong approach’ to determining how administration of CND impacts the injury cascade following seizure onset in this model. First, we compare disease-induced pathway alterations in untreated mice with those that are expected to occur after increased AT1R signaling in various disease models. Then we turn to the drug-induced pathway alterations and ask whether they are consistent with the known mechanism of action (MoA) of CND. ARBs are known to have a dual MoA resulting from AT1R antagonism and partial PPARγ agonism. Our comparative analysis evaluates whether observed pathway alterations reflect predicted effects of these MoAs in untreated and treated mice.
Hippocampal gene expression signatures of untreated juvenile and adult mice
At baseline, the excitability of the juvenile or neonatal brain differs from that of the adult brain. In neonates and juveniles, the brain is undergoing rapid growth and maturation, with ongoing synaptogenesis, dendritic arborization, and myelination. These developmental processes affect the balance between excitation and inhibition, making the immature brain more prone to hyperexcitability compared with the fully matured adult brain [42–44]. The more detrimental pattern in juveniles (Figure 6B) along with the increased number of DEGs (Figure 5B), may indicate increased vulnerability to the harmful effects of seizures. Despite these differences there are many similarities in the ways juvenile and adult brains respond to epileptic conditions. In this section we compare genome-wide patterns of transcript abundance in untreated juveniles and adult males, first on shared pathways and then on unique aspects of each age group. For a more complete discussion of genome-wide changes in gene expression in untreated adult B6-D/+ females and males, see [6].
Shared pathways
Figure 5C shows the 32 pathways shared by untreated juvenile and adults.
Identification in independent experiments lends credence to the veracity of the enrichment results and the importance of these pathways in the epileptic brain. The topmost significantly activated pathways (mean -log(P-values) ≥7.0) include pulmonary fibrosis idiopathic signaling, S100 family signaling, wound healing, hepatic fibrosis signaling, role of JAK family kinases in IL-6 type cytokine signaling, IL-6 signaling, and osteoarthritis signaling (Figure 5). In the context of the CNS, activation of the two fibrosis pathways implicates the process of astrogliosis and processes involving TGF-β signaling (also activated in both age groups), extracellular matrix (ECM) remodeling, and collagen deposition [45]. In addition to TGF-β, these pathways share a subset of genes with wound healing, including collagens and laminins that comprise the basement membrane component of the BBB [6]. While initially beneficial, these processes can be exacerbated after seizures are established, leading to BBB dysfunction, astrogliosis and scar formation [45]. Indeed, direct activation of the TGF-β pathway by TGF-β1 – one of the top predicted activators in untreated mice – results in epileptiform activity similar to that after exposure to albumin [46]. S100 family signaling proteins are known to function in stress responses and a variety of Ca2+-dependent intracellular functions and represent biomarkers of neuronal injury [47,48]. The JAK-STAT signaling pathway, including the JAK family of kinases and IL-6 cytokine signaling are known to play a significant role in neuroinflammation following brain injury [49], including in the development of epilepsy [50].
Deactivation of RHODGI and PPAR signaling (discussed below) were also found for both age groups (Figure 5C). The deactivation of RHOGDI, a negative regulator of the Rho family of GTPases, can result in the constitutive activation of multiple Rho GTPases, including RhoA2 and its downstream effector Rho kinase (ROCK) [51–53]. This in turn can promote actin cytoskeleton rearrangements in brain microvascular endothelial cells, leading to increased permeability of the BBB [54,55]. Activating the RhoA/ROCK1 pathways also has been linked to neuronal death and subsequent production of tumor necrosis factor-α (TNF-α) and IL-1β after traumatic brain injury (TBI) [56].
We note a clear enhancement in detrimental effects associated with shared pathways (Figure 6B). Pathways with the most detrimental effects (n≥6) include both deactivated pathways (PPAR and RHOGDI signaling), and the six activated pathways with asterisks in Figure 5C: LPS-stimulated MAPK signaling, IL-6 Signaling, iNOS signaling, Neuroinflammation Signaling, TGF-β signaling, and role of JAK family kinases in IL-6-type cytokine signaling. The detrimental nature of the effects of these activated pathways may represent the signal of an advanced neuroinflammatory process, possibly mediated by activated M1 pro-inflammatory microglia [57]. The pro-inflammatory cytokines produced through the LPS-stimulated MAPK and HMGB1 signaling pathways can lead to a cytokine storming [58] and the propagation of seizure activity [59,60]. Despite this increased detrimental signal among shared pathways, we note activation of some pathways with potential beneficial effects, such as RAC, Paxillin and FAK signaling – which may represent compensatory mechanisms to protect against BBB dysfunction and mitigate inflammatory responses during chronic injury [61–63] (Figure 6C).
Pathways unique to juveniles
There were 19 dysregulated pathways unique to juveniles (Figure 5D; Supplementary Table S3). Activation of CSDE1 signaling could potentially influence multiple pathways as it results in increases and decreases in the abundance of RNAs [64]. Activation of the production of nitric oxide and reactive oxygen species (ROS) by activated macrophages contributes to OXS and neuroinflammation, which can further exacerbate neuronal damage [57,65,66]. This inflammatory milieu may also activate signaling pathways such as TNFR1 and TNFR2 signaling, which mediate the effects of TNF-α on immune responses and cell survival [67]. The deactivation of eicosanoid signaling may result in prolonged inflammation and impaired resolution of neuroinflammation, exacerbating neuronal damage and tissue injury [68,69]. Additionally, OXS potentially can impact pathways such as white adipose tissue browning and urate biosynthesis/inosine 5'-phosphate degradation, which are sensitive to cellular redox status and metabolic changes [70].
Eleven of the 19 altered juvenile pathways are involved in cellular metabolism, energy homeostasis and neurotransmitter metabolism (Supplementary Table S3). While there isn’t a single common pathway that regulates all these pathways, metabolic pathways involved in nucleotide degradation (adenosine, purine, and guanosine), urea biosynthesis, and ethanol degradation are interconnected and regulated by common metabolic enzymes and signaling pathways.
For example, the enzymes involved in the degradation pathways of ethanol, leucine, histamine, noradrenaline, and dopamine (many of which localized in mitochondria) may be sensitive to OXS [8,71]. In response to seizures, the antioxidant defense mechanisms may be activated to counteract oxidative damage and potentially modulate the activity of these enzymes [72,73].
Effects of CND treatment on gene expression signatures in juvenile and adult mice
The greatly reduced number of DEGs, and the lack of significant pathway enrichment in juveniles treated with CND (Figure 5B) reflects transcript abundance patterns similar to those under physiological conditions. In treated juveniles, a relaxed FDR revealed a significant enrichment signal for OXPHOS and EIF2 signaling pathways. This pattern is very similar to that observed in untreated pre-seizure females, which was discussed previously [6]. Briefly, we suggested that the increase in OXPHOS gene expression was a beneficial response to cellular stress and the depletion of ATP in the hyperexcitable state before seizure onset [74,75]. The coordinated activation of EIF2 signaling may have represented a compensatory mechanism providing for the production of proteins involved in the ETC [76–78]. CND treatment may induce a similar process, boosting mitochondrial function in a manner similar to the mitohormesis effect observed in exercise-trained mice [75,79]. In so doing, CND treatment may have shifted the metabolic and redox balance toward physiological conditions.
Similar to juveniles, treated adults exhibited genome-wide transcript abundance profiles approximatlng those under physiological conditions; however, there was a set of 8 canonical pathways exclusively activated in treated adult males (Supplementary Table S3), 5 of which are involved in modulating the immune response [80–82] – specifically in enhancing the activation of Tregs. Tregs are a subset of T cells that play a crucial role in immune tolerance (i.e., dampening excessive immune responses and promoting immune homeostasis) [82,83]. Tregs have been shown to protect against pathological neuroinflammation by polarizing microglia toward a M2-like phenotype [84] in several injury models [84–86]. In epilepsy, brain Treg depletion reduces neuronal survival and exacerbates hippocampal neuroinflammation and OXS and after SE, while Treg enhancement exerts an antiepileptic effect [87]. The anti-inflammatory milieu so produced may protect neurons from secondary damage and contribute to the maintenance of seizure control by reducing excitability of neuronal networks and minimizing risk of seizure recurrence [87]. It is not possible to infer whether treatment with CND directly or indirectly led to the activation of Tregs and the long-term stabilization of the immune response.
Another interesting finding is the activation of the macrophage alternative activation signaling pathway (Supplementary Table S3). This pathway (also known as the ‘M2’ macrophage activation pathway), plays a significant role in the resolution of inflammation, tissue repair, and remodeling. In the early phases of CNS injury, skewing of microglia and macrophages toward pro-inflammatory M1-like phenotype is adaptive (e.g., aids in clearing dead cells and debris). Over time, switching to a M2-like anti-inflammatory phenotype promotes improved tissue repair and remodeling, BBB integrity, angiogenesis, and neuronal survival [65,88].
Are the effects of Scn8a-induced gene expression patterns consistent with the pathological effects of RAS as mediated by AT1Rs?
Many of the pathological conditions that are known to occur during epileptogenesis also occur in the Scn8a-N1768D mouse (Supplementary Table S3). Figure 6D presents an interactive model in which neuronal hyperexcitability associated with persistent NaV1.6 current triggers a sequential chain including glutamate excitotoxicity, calcium overload, OXS/MD, inflammatory responses, glial activation, BBB dysfunction, fibrosis, and seizures. The model includes several feedback loops and highlights the importance of neuroinflammatory processes [89].
Likewise, in the context of various brain disorders and injuries, Ang II is increased and AT1R signaling activated, triggering release of pro-inflammatory cytokines via activation of NF-kβ expression and leading to an injury cascade characterized by many of the abovementioned pathological processes [10,11,90–93]. Indeed, increased Ang II has been shown to play a pathophysiological role in epilepsy [16,17,94–97]. Experiments with LPS, which also activates AT1R signaling, point to positive feedback between Ang II and inflammation [98]. Additionally, there is crosstalk between the RAS and PPAR signaling pathways as AT1R overstimulation inhibits PPARγ [92] (see below).
Activation of AT1Rs involves stimulating kinase pathways, such as mitogen-activated protein kinases (MAPK, p38) and JAK/STAT pathways [11]. MAPK signaling was inferred to be activated in both untreated juvenile and adult mice as evidenced by activation of LPS-stimulated MAPK signaling and p38 MAPK signaling (Supplementary Table S3). Neuroinflammation can also be mediated by the release of damage associated molecular pattern mediators (DAMPs), which trigger inflammatory responses from microglia, astrocytes, and neurons [11,89]. Indeed, HMGB1 signaling was inferred to be activated in both untreated juvenile and adult mice, as were IL-6 signaling, neuroinflammation signaling, pathogen induced cytokine storm signaling, and JAK family kinases in IL-6-type cytokine storm signaling pathways (Supplementary Table S3).
Our results reveal that activation of the renin–angiotensin signaling and deactivation of PPAR signaling pathway are associated with the largest number of detrimental effects of any pathway in Supplementary Table S3 (Figure 6C). The top predicted upstream activators for both juveniles and adult males (LPS, TNF-α, IL-1B, and TGFB1) are key mediators of inflammation and inducible by Ang II via microglia activation [92,99] (Supplementary Table S4). The top predicted upstream inhibitor was PPARGC1A (PGC-1α), a coactivator recruited during PPAR stimulation. Among the targets of PGC-1α are several mitochondria-related genes that play important roles in mitochondrial oxidative metabolism and biogenesis, as well as in the maintenance of intracellular calcium levels [100]. We infer from these results that Scn8a pathophysiology closely mimics (and may indirectly result from or converge on) the pathological processes associated with overstimulation of AT1R and deactivation of PPAR signaling.
Are the effects of CND treatment consistent with the known mechanism of action of ARBs?
Nearly all of the above pathway alterations identified in untreated were ‘normalized’ (i.e., no longer significantly enriched by DEGs) in both treated juvenile and adult male mice (Supplementary Table S3). Here we evaluate whether the known MOAs of CND can explain these results. ARBs are known to have powerful anti-inflammatory, anti-fibrotic, and antioxidant activities [11,13,14]. The anti-inflammatory action of ARBs may involve both their AT1R antagonist and partial PPARγ agonist activities. One of the known anti-inflammatory effects of ARBs derives from their ability to inhibit HMGB1 signaling (see above) [11]. Activation of widely expressed PPAR receptors leads to additional anti-inflammatory, neuroprotective, and metabolic effects [14]. PPARγ agonism can block NF-κB signaling in the brain, reduce the production of ROS, mitigate microgliosis, and promote mitochondrial function [14,101]. It is important to note that AT1R blockade and PPARγ agonism can act as independent anti-inflammatory mechanisms; however, as stated above, AT1R signaling suppresses PPARγ induction and PPARγ signaling reduces AT1R activation [11].
In animal models of TBI, AT1R blockade has been shown to reduce OXS and stabilize the BBB. CND in particular has consistently demonstrated beneficial effects by reducing expression of genes involved in the inflammatory response and microglia activation, and by inhibiting IL-6, IL-1β, iNOS, TNF-α, and TGF-β1 signaling. By inhibiting signaling of the pro-fibrotic cytokine TGFβ, ARBs can modulate ECM remodeling, reduce fibrosis and gliotic scar formation, and mitigate the development of recurrent seizures [11]. Most of the studies that have explored the use of ARBs to control vulnerability to seizures [16] have focused on losartan. This ARB was shown to mitigate BBB breakdown and status epilepticus (SE)-induced epileptogenesis in rats by blocking TGFβ signaling [15,102], and to inhibit production of ROS by preventing NOX2 activation [15,17,103–105]. BBB dysfunction and gliosis were improved in a rat pilocarpine model [103]. Treatment with another ARB, telmisartan, reduced seizure frequency in dogs [106].
We were not able to find another published study that examined the effects of CND directly in an epilepsy model; however, a number of studies have examined CND effects in related models. For example, CND was able to suppress acute hypoxia-induced seizures and reduce later-life seizure susceptibility in a neonatal mouse hypoxia model [107]. In rats treated with LPS, CND administration ameliorated microglial activation, reduced NF-kβ signaling, and lowered the level of COX-2 and the generation of ROS [108,109]. Treatment of neuronal cells with CND normalized the expression of hundreds of genes associated with glutamate-induced inflammation, diabetes signaling, and amyloid β metabolism [9]. Similar to this study, the top predicted upstream activators for neuronal cells stimulated with excitotoxic levels of glutamate (and down-regulated by CND) included LPS, TNF-α, IL-1B, and TGFB1 (Supplementary Table S4).
While we cannot tease apart the relative roles of AT1R and PPAR signaling in adults, PPAR agonism may play a larger role in the juvenile case. RAS signaling is not activated in untreated juveniles, and there is indirect evidence for increased detrimental effects of OXS (Figure 6A). The beneficial effects of CND treatment may thus relate to its ability to induce PPARγ activity, which is known to mitigate OXS and tissue damage. Indeed, the PPARγ agonist pioglitazone has been shown to increase the expression of ETC subunit proteins [101]. In the case of the adult brain, CND treatment produced an ‘anti-similar’ result in which cyclooxygenase-2 (COX-2) was predicted as an upstream activator in untreated adults and as an upstream inhibitor in treated mice (Supplementary Table S4). Both AT1R blockade and PPAR activation mediate down-regulation of a wide array of pro-inflammatory mediators in microglia, including TNF-α, IL-1β, IL-6, iNOS, and COX-2 [13,14] (Supplementary Figure S5). Under pathologic conditions, COX-2 activity can produce ROS and toxic prostaglandin metabolites that can exacerbate brain injury [110]. A similar anti-similar pattern held for ERBB2 and NUPR1, the latter is activated in response to OXS [111]. Inhibition of ErbB receptor activity, which is elevated in reactive astrogliosis, suppresses hypertrophic remodeling via reduced Src/FAK activity [112].
Conclusions and challenges
These results support the hypothesis that the efficacy of CND in our novel model of epileptogenesis is related to its known dual MoA [11] (Figure 6 D). The robust efficacy of CND across ages, sexes and mouse strains is a positive indication for successful translation and drug repurposing in humans. Evidence for its efficacy in in the context of the juvenile brain is critical given the early onset of seizures in children with SCN8A-related epilepsy. Interestingly, results of recent large cohort studies of adults with hypertension indicated that use of ARBs was associated with a significantly decreased incidence of epilepsy compared with other classes of antihypertensive drugs [113,114]. Importantly, the use of CND has been approved to treat hypertension in adolescents and children aged one and older, reflecting its excellent safety profile and suitability for use in clinical trials for children with SCN8A epilepsy [115]. In terms of lifetime survival in Scn8a-D/D juveniles, CND outperformed PHT – an often-prescribed ASM for epilepsy [40]. In contrast with many ASMs that primarily target ion channels or neurotransmitter systems to reduce seizure occurrence, ARBs appear to address underlying pathological processes involved in epileptogenesis, such as inflammation, oxidative stress, and BBB dysfunction. Importantly, the effect of CND on gene expression and multiple pathological processes suggest it may offer broader neuroprotective benefits beyond seizure control. Interestingly, removal of administration of Prax-330 from B6-D/+ mice resulted in rapid loss of protection, suggesting that it acts primarily as an anticonvulsant rather than a disease-modifying therapy [41]. Future studies of CND and gene-based treatments [116] will help to assess to what extent these different therapeutic approaches are truly disease-modifying. Alternatively, treatment strategies incorporating CND as an adjuvant with sodium channel blockers may also yield improved outcomes for seizures and co-morbid conditions [16]. Finally, many cellular and molecular processes identified here are shared with other forms of epilepsy [6], brain injury [117], and inflammatory brain disorders [6,10,12,39]. In that light, the translational data set presented in this study supports the hypothesis that ARBs may be a useful therapeutic strategy for other neurodevelopmental disorders with similar primary hallmarks [12,39].
Challenges
The response to brain insult is highly complex, involving a multitude of signaling pathways that interact in a dynamic and context-dependent manner. The balance between these signals can influence whether the outcome leans more towards regeneration and recovery or towards scar formation and chronic impairment [118]. Many pathways represent a double-edged sword in that they are associated with both beneficial and detrimental effects, depending on the temporality, duration, and magnitude of pathway (de)activation, as well as on the specific cellular and molecular context in which the pathway(s) operate [66,119,120]. In general, we interpreted the effects of pathway alterations under the assumption that they were occurring in the later stages of epileptogenesis (i.e., after the establishment of spontaneous recurrent seizures). One advantage of the systems approach taken here is that comparison of drug-induced gene expression profiles can identify potential off-targets [32] – of which few, if any, were identified. We also faced the previously discussed methodological challenge [6] of false positives called by IPA [121]. The approach taken here helped to mitigate this issue to some extent by focusing on a set of the most informative and commonly reported pathway effects (Supplementary Table S3). This form of ‘feature selection’ reduces the impact of false positives. None of the 12 pathway effects were highly correlated suggesting that this set of effects captures distinct features of clinical presentation. However, we note that alternative inferences regarding beneficial or detrimental effects are possible and that additional features/effects may have been excluded.
Clinical perspectives
Epilepsy is a common neurological disease; however, few if any of the currently marketed antiseizure medications prevent or cure epilepsy. Angiotensin receptor blockers (ARBs) represent a promising strategy to treat the numerous detrimental consequences of epileptogenesis.
Scn8a pathology mimics the known effects of angiotensin type 1 receptor (AT1R) overstimulation and thus, blocking AT1R activation has potential therapeutic value in this model and other neurological disorders with similar pathological features.
Results of testing candesartan, an FDA-approved ARB indicated for hypertension in pediatric patients older than 1 year, on this transgenic Scn8a mouse model yield strong evidence for its efficacy and support its use in clinical trials for children with SCN8A-related epilepsy.
Data Availability
The data that support the findings of the present study are available from the corresponding author on reasonable request. Raw data of RNA sequencing have been submitted to the Figshare repository at https://portlandpress.figshare.com/articles/dataset/RNAseq_readcount_files_for_untreated_and_treated_juveniles/26417218
Competing Interests
The authors declare that there are no competing interests associated with the manuscript.
Funding
P.T.R. was funded by a Translational Research Grant from the American Heart Association [grant number 19TPA34910113]. This work was also supported by funding from the Shay Emma Hammer Research Foundation.
Open Access
Open access for this article was enabled by the participation of University of Arizona in an all-inclusive Read & Publish agreement with Portland Press and the Biochemical Society under a transformative agreement with Individual.
CRediT Author Contribution
Michael F. Hammer: Conceptualization, Resources, Data curation, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Visualization, Methodology, Writing—original draft, Project administration, Writing—review & editing. Erfan Bahramnejad: Data curation, Investigation, Methodology. Joseph C. Watkins: Formal analysis, Validation. Patrick T. Ronaldson: Data curation, Formal analysis, Funding acquisition, Writing—original draft.
Ethics Statement
Ethical approval for animal work was obtained from the University of Arizona Institutional Animal Care and Use Committee Program (IACUC #16-160). All experiments were designed in accordance with the Animal Research: Reporting In Vivo Experiments (ARRIVE) guidelines. Animal experiments were conducted in approved facilities in the laboratories of MH and PTR at the University of Arizona. At the conclusion of each experiment or when necessary due to morbidity, animals were euthanized using ketamine/xylazine followed by cervical dislocation.
We confirm that we have read the Journal's position on issues involved in ethical publication and affirm that this report is consistent with those guidelines.
Acknowledgements
We thank Aurora Hurtado, Emily Barnes, Joshua Hack and Sarah Lester for technical support. We also thank the Shay Emma Hammer Research Foundation for funding.
Abbreviations
- ANG
angiogenesis
- AngII
Angiotensin II
- APO
apoptosis/clearance
- ARB
angiotensin receptor blocker
- ARRIVE
Animal Research: Reporting In Vivo Experiments
- ASM
antiseizure medication
- AT1R
angiotensin type 1 receptor
- ATP
adenosine triphosphate
- B6
C57BL/6J
- BBB
blood–brain barrier
- BBBD
blood–brain barrier dysfunction
- C3H
C3H/HeJ
- CND
candesartan
- COX-2
cyclooxygenase-2
- CSDE1
Cold Shock Domain Containing E1
- DAMP
damage-associated molecular pattern mediator
- DEG
differentially expressed genes
- ECM
extracellular matrix
- EIF2
eukaryotic initiation factor 2
- ERBB2
erythroblastic oncogene B
- ETC
electron transport chain
- FAK
focal adhesion kinase
- FBS
fibrosis/scar formation
- FDA
Food and Drug Administration
- FDR
false discovery rate
- GLA
glial activation
- HGF
hepatocyte growth factor
- HMGB1
high-mobility group box 1
- IACUC
Institutional Animal Care and Use Committee Program
- INF
neuroinflammation
- iNOS
inducible nitric oxide synthase
- IPA
Ingenuity® Pathway Analysis
- LPS
lipopolysaccharide
- MAPK
mitogen-activated protein kinases
- MD
mitochondrial dysfunction
- MDS
multidimensional scaling
- MoA
mechanism of action
- NDG
neurodegeneration
- NER
neurogenesis/repair
- Nox2
NADPH oxidase 2
- NUPR1
nuclear protein 1
- OXPHOS
oxidative phosphorylation
- OXS
oxidative stress
- PGC-1α
peroxisome proliferator-activated receptor-γ coactivator
- PHT
phenytoin
- PI3K/Akt
phosphoinositide-3-kinase/protein kinase B
- PPAR
peroxisome proliferator-activated receptor
- RAS
renin–angiotensin system
- RHODGI
Rho GDP-dissociation inhibitor
- ROCK
Rho-associated coiled-coil containing protein kinase
- ROS
reactive oxygen species
- s.c.
subcutaneous
- SIC
secondary injury cascade
- Src
proto-oncogene tyrosine-protein kinase
- STAT
Janus kinase
- SYN
synaptic plasticity
- TBI
traumatic brain injury
- TC
tonic-clonic seizure
- TGF-β
transforming growth factor-β
- TNF-α
tumor necrosis factor-α
- TNFR
tumor necrosis factor receptor
- Tregs
regulatory T cells
- VEH
vehicle