Gene Expression in Temporal Lobe Epilepsy is Consistent with Increased Release of Glutamate by Astrocytes

Abstract

Patients with temporal lobe epilepsy (TLE) often have a shrunken hippocampus that is known to be the location in which seizures originate. The role of the sclerotic hippocampus in the causation and maintenance of seizures in temporal lobe epilepsy (TLE) has remained incompletely understood despite extensive neuropathological investigations of this substrate. To gain new insights and develop new testable hypotheses on the role of sclerosis in the pathophysiology of TLE, the differential gene expression profile was studied. To this end, DNA microarray analysis was used to compare gene expression profiles in sclerotic and non-sclerotic hippocampi surgically removed from TLE patients. Sclerotic hippocampi had transcriptional signatures that were different from non-sclerotic hippocampi. The differentially expressed gene set in sclerotic hippocampi revealed changes in several molecular signaling pathways, which included the increased expression of genes associated with astrocyte structure (glial fibrillary acidic protein, ezrin-moesin-radixin, palladin), calcium regulation (S100 calcium binding protein beta, chemokine (C-X-C motif) receptor 4) and blood-brain barrier function (Aquaaporin 4, Chemokine (C-C-motif) ligand 2, Chemokine (C-C-motif) ligand 3, Plectin 1, intermediate filament binding protein 55kDa) and inflammatory responses. Immunohistochemical localization studies show that there is altered distribution of the gene-associated proteins in astrocytes from sclerotic foci compared with non-sclerotic foci. It is hypothesized that the astrocytes in sclerotic tissue have activated molecular pathways that could lead to enhanced release of glutamate by these cells. Such glutamate release may excite surrounding neurons and elicit seizure activity.

Introduction

Medically refractory temporal lobe epilepsy (TLE) is a major seizure disorder characterized by a unique pathology first described in 1880 as ‘hippocampal sclerosis’ or ‘Ammon’s horn sclerosis’. This is a distinct neuropathological condition characterized by hippocampal atrophy, induration, distinctive neuronal loss, and astroglial proliferation (1,2) seen very prominently in area CA1 or the Sommer sector of the hippocampus.

A careful study of the pathology and electrophysiology of 151 hippocampi removed in the Yale surgical series has been undertaken (3). This study revealed that about 60% of the hippocampi had Ammon’s horn sclerosis and were described as mesial temporal lobe epilepsy (MTLE). Another 28% was associated with an extrahippocampal mass lesion and did not show sclerosis; this condition was described as mass associated temporal lobe epilepsy (MaTLE). A further 12% had no discernible pathology and their hippocampi also did not have sclerosis; a group described as paradoxical TLE (PTLE).

Several lines of evidence point to the sclerotic hippocampus, rather than the amygdala or other temporal lobe regions, as the most likely place of origination and maintenance of chronic seizures in TLE patients (4). Intracranial EEG recordings have shown that the seizures in MTLE originate in the sclerotic hippocampus — severe neuronal loss (sclerosis) in the anterior hippocampus correlates with anterior focal onset, whereas sclerosis in both anterior and posterior hippocampus was associated with seizure onset from a wider area of the hippocampus (5). It has been argued that hippocampal sclerosis is merely a consequence of prolonged seizures. If this were the case, then all epilepsies involving the temporal lobe should be associated with hippocampal sclerosis (6). However, PTLE and MaTLE are not associated with hippocampal sclerosis. The clinical presentation of PTLE is virtually indistinguishable from MTLE but its focus of origin is yet unclear (3). Additionally, hippocampal volume measurements by MRI as an indicator of sclerosis have shown that the degree of hippocampal atrophy is not correlated with the duration and severity of seizures (7). Hippocampal sclerosis was also reported to be absent in some patients with multiple daily seizures since infancy, again suggesting that seizures are not invariably a cause of sclerosis (8). A postmortem study of hippocampi from patients with poorly controlled seizures has confirmed a subgroup with absence of significant neuronal loss despite decades of seizures (9). Thus hippocampal sclerosis appears to be a distinctive pathological entity that is correlated with and possibly causally related to MTLE.

Strong evidence for the role of the sclerotic hippocampus in the maintenance of seizures comes from the success of resective surgery for the control of intractable epilepsy. Surgical removal of the sclerotic hippocampus, rather than adjacent regions, results in the best seizure free outcome (evidence reviewed in 10). In a more recent seven-center prospective observational study of resective epilepsy surgery, the authors examined the probability and predictors of patients entering two-year remission. They observed that about 70% of patients with medial temporal resections experienced two-year remission and found that only a history of absence of tonic-clonic seizures and the presence of hippocampal atrophy, were significantly and independently associated with remission (11).

Despite the above evidence, there is continuing debate as to the relationship of hippocampal sclerosis to the causation and maintenance of seizures in temporal lobe epilepsy. This debate has continued in the absence, until recently, of much information on the role of astrocytes in the pathology of disease or the organization and function of the sclerotic substrate. Astrocytes in sclerotic regions of the hippocampus have many unique properties, and may contribute to the hyperexcitability of the hippocampal seizure focus. The astrocytes have an increase in Na⁺ channels on their membranes, a more depolarized resting membrane potential, and a reduction of inward rectifying potassium (Kir) channel function (12,13); have an elevated flip-to-flop mRNA ratio of GluR1 that make them more responsive to glutamate (14); have a significant down regulation of the enzyme glutamine synthetase (15); and show increased expression of the transcription factor Nuclear Factor kappa B (NFκB), which is thought to play a crucial role in many acute inflammatory reactions (16).

This paper further explores the molecular organization of the sclerotic substrate in TLE and proposes the hypothesis that the molecular changes in astrocytes from sclerotic hippocampi are associated with calcium dependent glutamate release into the extracellular space, and that several important initial precipitating events in the causation of TLE may activate this mechanism.

Materials and Methods

Samples, Extract Preparation, and Labeling

The hippocampi used in this study were surgically removed from patients with medically intractable TLE who underwent anteromedial temporal lobectomy with hippocampectomy for seizure control, by procedures previously described (17). Tissue used in these studies was obtained after informed consent and with approval of the human investigations committees (HIC) of participating institutions. The patient characteristics are shown in Table 1. Twenty hippocampal specimens were subcategorized into mesial temporal lobe epilepsy with sclerosis (MTLE, n = 8), paradoxical temporal lobe epilepsy (PTLE, n = 6), mass associated temporal lobe epilepsy (MaTLE, n = 6); the latter two groups had no sclerosis. In the MaTLE patients used in this study, the mass lesion was outside the hippocampus. Subcategorization (3) was done on adjacent tissue sections from each patient that were immunostained for neuropeptides (18) and on which neuronal counts were performed (19).

Table 1 Clinical Features Related to Patient Tissue Used in Microarray Analysis

Hippocampal tissue slices, 5 mm or less in thickness, were collected in RNAlater (Ambion, Austin, TX) within 30 min of surgical removal of the hippocampus and stored overnight at 4° C to allow for penetration of the RNAlater, then frozen on dry ice and stored at −80° C. RNAlater has been shown to preserve RNA in tissue samples at least as well, if not better than, rapid freezing in liquid nitrogen (20).

The hippocampus was microdissected into subregions and total RNA extracted from the CA1 subregion using the RNeasy kit for lipid tissue (Qiagen, Valencia, CA.). The quality of the RNA was assessed by the Agilent Bioanalyzer (Agilent Technologies, Palo Alto, CA.) to visually verify the absence of genomic DNA contamination and integrity of 28S and 18S bands. The CA1 region was chosen, as this region, known as the Sommer sector (2), characterizes the sclerotic hippocampus.

The preparation of labeled cRNA and the subsequent hybridization to human U133A arrays was performed following guidelines provided by the Affymetrix Resource facility at Yale University (https://doi.org/keck.med.yale.edu/affymetrix). The U133A chip was chosen as it contained all the well-characterized genes at the time. Patient RNA was not pooled for this study. The RNA from each patient was profiled on a separate GeneChip.

Measurement of Data and Analytical Tools

The hybridized arrays were scanned using a confocal laser fluorescence scanner (Affymetrix GeneChip Scanner 3000). The two-way ANOVA model was applied to identify differentially expressed genes using the PM-only model-based expression values obtained from dChip, which implements the Li-Wong model (21).

To look at the relatedness of patient gene expression profiles, Multidimensional Scaling (22) analysis was conducted on representations of the expression of all 22,215 probe sets (14,500 genes). This approach projects the high-dimensional gene expression data onto two viewable lower dimensions. In this approach the distance between data points provides a measure of their differences. Another unsupervised hierarchical cluster analysis using Pearson correlation to measure the similarity of expression between samples (23,24) was also carried out on the microarray data set.

To study hippocampal sclerosis in MTLE, the non-sclerotic PTLE and MaTLE hippocampi were used for the comparison. The PTLE and MaTLE samples were used as they closely resembled the normal hippocampus in anatomy, their exposure profile to antiepileptic drugs was similar to the MTLE samples (Table 1) and they were subjected to the same surgical manipulations in removal. Further, while the MTLE hippocampus is an epileptogenic focus, the hippocampus of MaTLE and PTLE are not. Evidence in support of this comes from electrophysiological investigation of dentate granule cell excitability, which shows that the MTLE granule cells are hyperexcitable while those in MaTLE and PTLE are not (25). A study of the long term surgical outcome of patients with temporal lobe tumoral epilepsy (MaTLE) showed significant post operative seizure freedom was related only to gross complete tumor resection and unrelated to its pathology (26). This, too, suggests that in MaTLE the hippocampus itself is non-epileptogenic. The poor seizure free outcome (44%) in PTLE following hippocampectomy suggests that the hippocampus is unlikely to be epileptogenic in this group as well (3). Hippocampal tissue from autopsies of “neurologically normal” subjects was not used as a comparison or “control” group for this study as this tissue is subject to longer and variable post-mortem delay before tissue harvest, uncertain periagonal conditions, and other variable factors, which are shown themselves to influence gene expression patterns (27). Autopsy specimens also lack exposure to AEDs and anesthetics or exposure to seizure activity that surgical tissue has. Thus they are not a good tissue source for controlling the experimental variables in an experiment with surgically obtained tissue from epileptic patients.

Biological Pathways Analytical Tools

An extensive literature search has to be done to establish the functional roles, interactions, and relationships between the differentially expressed genes. This is considerably assisted by proprietary software, of which we chose Ingenuity Pathways Analysis (Ingenuity Systems, Mountain View, CA). The Ingenuity Pathways Knowledge Base is the largest curated database of previously published findings on mammalian biology from the public literature (https://doi.org/www.ingenuity.com). The Ingenuity Knowledge Base comprises a genomic database of over 10,000 human genes that were manually curated by reviewing over 200,000 full-text scientific articles supplemented with curated information parsed from MED-LINE abstracts. In addition, the database comprises biological networks created from millions of individually modeled relationships derived from findings reported in peer-reviewed scientific publications and constantly updated. Content and modeling experts systematically encoded these findings to derive binding relationships, post-translational modifications, expression regulation, transport, and other types of interaction. Hence, a molecular network of direct, physical, transcriptional, and enzymatic interactions between orthologues with high connectivity can be computed using this knowledge base. Each gene interaction in the biological network is supported by annotated literature (28)

From data generated by the ANOVA analysis, the added criteria of false discovery rate (FDR) of < 0.05, and fold-change cut-off of 1.5 were used. Because the brain is histologically heterogeneous, with often small numbers of different cell types, gene expression levels may not always be large, and it has been proposed that a 1.5 fold change is reasonable to use (29,30). Each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. These genes, called focus genes, were then used as the starting point for generating biological networks. These genes were overlaid onto a global molecular network developed from the information contained in the Ingenuity Pathways Knowledge Base. Networks of these focus genes were then algorithmically generated based on their connectivity and interactions. Ingenuity Pathway Analysis then computed a score for each network according to the fit of the set of the significant genes, which was derived from a P value and indicates the likelihood of the focus genes in a network being found together due to random chance. Only networks with at least a 99% confidence of not being generated by random chance alone were used for this study.

In addition, Ingenuity Pathways Analysis can generate Global Canonical Pathways, which associates genes from the dataset file with pathways of known metabolic and signaling pathways to examine their role in the context of these pathways. The significance values (P values) are measures for the likelihood that genes from the dataset file under investigation participate in that function and are calculated using the right-tailed Fisher’s Exact Test by comparing the number of genes from the dataset that participate in a given function or pathway, relative to the total number of occurrences of these genes in all functional/pathway annotations stored in the Ingenuity Pathways Knowledge Base. The significance of the association between the data set and the canonical pathway was measured in two ways: 1. The ratio of the number of genes from the data set that map to the pathway divided by the total number of genes that map to that canonical pathway; 2. Fisher’s exact test was used to calculate a P value determining the probability that the association between the genes in the dataset and the canonical pathway is explained by chance alone.

Validation of Microarray Data by Quantitative RT- PCR (QRT-PCR)

QRT-PCR was performed on five genes along with a housekeeping gene, β-actin or Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), for each sample as a means of validating the microarray data. The genes were chosen on the basis of previous reports of their presence in astrocytes. Each QRT-PCR reaction was performed with RNA from the same preparation used for the microarray analysis. Total RNA (2 µg) from each sample was subjected to reverse transcription using High Capacity cDNA Archive Kit (Applied Biosystems, Foster City, CA). QRT-PCR analysis was performed using ABI 7900 Sequence Detection System using Assays-on-Demand products for AQP4, CD44, CXCR4, GFAP, and PLEC1. All reactions were performed in triplicate in 15 µL final volumes. The following PCR conditions were used: 50° C for 2 min, then 95° C for 10 min, followed by 40 cycles at 95° C/0.15 min and 60° C /1 min. Standard curve was generated for each gene using cDNA obtained by pooling equal amounts from each sample. The expression level of target genes was normalized to internal β-actin or GAPDH. Data was analyzed using Microsoft Excel following ABI instructions, and calculated using the relative standard curve method (ABI, User Bulletin #2).

Immunohistochemistry

An approximately 3 mm thick slice of tissue adjacent to the one used for DNA microarray studies was fixed by immersion in a 4% paraformaldehyde solution in phosphate buffer (pH 7.4) for 6–8 h Vibratome sections (50 µm thick) were immunostained according to the Avidin-Biotin Complex method (31) with primary antibodies against AQP4 (1:10,000, Chemicon, Temecula, CA), DMN (1:5000, Novocastra, Newcastle, UK), CD44 (1:500, RDI, Flanders, NJ), PLEC1 (1:500, BD Transduction Laboratories) and CXCR4 (1:1000, ProSci Inc., Poway, CA) using procedures previously described (18).

Results

Patient Characteristics

The types of the hippocampal specimens used in this study were identified on the basis of immunohistochemistry for neuropeptides and neuronal density assessments (3). The MTLE specimens (n = 8) all had hippocampal sclerosis, i.e. significant neuronal loss with accompanying gliosis in area CA1 (Table 1). None of the MaTLE (n = 6) or PTLE (n = 6) specimens showed hippocampal sclerosis. The MTLE specimens had about a four-fold loss in neurons and a three-fold increase in glia in area CA1 compared with the MaTLE and PTLE specimens (Table 2). The MaTLE patients had an extrahippocampal temporal mass lesion, whereas the PTLE had no discernible pathology or identifiable etiology (Table 1). These features correspond to the previously described pathology for larger samples of similar specimens (3).

Table 2 Neuronal and Glial Counts in Surgically Removed Hippocampi used in this Study

Whereas there was a spread in the age of patients at the time of surgery, they were similar across the three groups (mean ± SD, MTLE 41.0 ± 11.9, MaTLE 26.2 ± 23, PTLE 31.1 ± 16, Table 1). The age at the first unprovoked seizure in MTLE was significantly younger (mean ± SD, 6.7 ± 3.3 years) than in either MaTLE (20 ± 25 years) or PTLE (13.8 ± 14.3 years); the age in the latter two groups being similar. The number of years to surgery (age at surgery minus age at first seizure) also showed differences. This period was similar in the MTLE (34.3 ± 11 years) and PTLE (17.3 ± 15.6 years) groups and significantly (Mann-Whitney U, two tailed, P ≤ 0.05) longer than in the MaTLE group (6.4 ± 5.5 years), and was consistent with previous observations on larger patient populations (3). All of the patients in this study had seizures that were uncontrolled by antiepileptic medications (AEDs). Examination of the lists of AEDs being used at the time of surgery or in AED use history showed no obvious difference between patient groups.

Patient Categorization by Gene Expression Patterns

One objective of this study was to determine if there were distinct molecular signatures in the transcriptional profiles of the different hippocampal samples studied. To this end, two approaches, with somewhat different underlying mathematical assumptions, were adopted. Multidimensional scaling, an unsupervised method of data reduction, was conducted on the representations of the expression of all the genes (22,215 probe sets representing 14,500 genes) and demonstrated that gene expression profiles of the eight MTLE hippocampi with sclerosis were grouped together and separable from the other hippocampi without sclerosis (data not shown).

Unsupervised hierarchical cluster analysis was also performed on all samples. This unbiased approach also organized the hippocampal specimens into two major groups. The MTLE hippocampi formed one tight cluster while the non-sclerotic hippocampi (MaTLE, PTLE) formed another group (P = 0.000002, Fisher Exact test). The hippocampal transcriptional profile of the MaTLE and PTLE patients were not distinguishable from each other, at least not from the number of samples used in our study (Figure 1).

Figure 1

The comparison of the gene expression patterns in 20 hippocampi from patients with TLE. The dendrogram was obtained by unsupervised hierarchical cluster analysis using Gene Spring software. The analysis included all genes expressed in a microarray analysis using the Affymetrix GeneChip U133A. It shows a primary branching pattern, which separates the sclerotic MTLE (dark green) hippocampi as one cluster. The PTLE (blue) and MaTLE (Red), both being non-sclerotic hippocampi, form a separate group. On the dendrogram, each horizontal color bar denotes the intensity of expression of a particular gene relative to the mean expression value of all samples. The color codes are shown in the color bar where green represents transcript levels below the median and red represents transcript levels above. Yellow is at median. The more intense the color, the greater the level of trust it represents. MTLE = mesial temporal lobe epilepsy; MaTLE = mass associated temporal lobe epilepsy; PTLE = paradoxical temporal lobe epilepsy.

Reproducibility and Distinctiveness of Transcription Profile in Sclerotic (MTLE) Hippocampi

A second objective of our analysis was to identify the population of genes whose expression was uniquely changed in the MTLE hippocampi. The assessment by ANOVA was used to discover genes that were differentially expressed across MTLE, PTLE, and MaTLE. The mean of the expression of each transcript in the eight MTLE hippocampi was compared with the mean expressions of the six PTLE and the mean expression of the six MaTLE hippocampi separately, generating two sets of comparisons (MTLE vs. PTLE, and MTLE vs. MaTLE). Only the gene transcripts that were differentially expressed in both comparisons were listed as significant genes for further evaluation, based on a false discovery rate (FDR) cut-off at < 0.05. The FDR is a multiple-hypothesis testing error measure, and is the expected proportion of false positive findings among all the rejected hypotheses (32). It is more appropriate to control for false positive results than the Bonferroni correction. Using the fold-change convention, we found 149 gene transcripts that were changed two-fold and greater in either direction for both comparisons — MTLE vs. PTLE and MTLE vs. MaTLE. Of these, 67 were up-regulated and 82 were down-regulated. If the fold change criteria were expanded to include 1.5-fold and greater, then there were a total of 674 gene transcripts, of which 366 were upregulated and 308 were down-regulated. (The complete list of genes is given in supplement 1, available on the Molecular Medicine Website (https://doi.org/www.molmed.org).)

We also conducted a separate analysis whereby MaTLE and PTLE were considered replicates of non-sclerotic hippocampi. This comparison of sclerotic vs. non-sclerotic specimens showed 138 genes changed two-fold (76 up-regulated and 66 down-regulated), while 947 genes were changed 1.5-fold or more (535 up-regulated and 412 down-regulated). (The complete list of genes is found in supplement two, available on the Molecular Medicine Website (https://doi.org/www.molmed.org).)

To assess the degree of variability in transcription patterns between individual MTLE specimens, correlation coefficients were calculated. The correlations ranged from 0.95 to 0.98, and demonstrated a high degree of similarity between MTLE replicates (33).

Identity of Genes Expressed in MTLE Specimens

One of the challenges in DNA microarray analyses is to find biological meaning in the lists of differentially expressed genes. To facilitate this process, several approaches were adopted. In one approach the list of genes and ESTs showing changes in expression of 1.5-fold and above in either direction was submitted to Ingenuity Pathway Analysis Software to identify the genes and ESTs belonging to canonical pathways. Among the genes associated with major canonical pathways, those involved in glutamate and GABA receptor signaling were down-regulated. Likewise a number of genes associated with G-protein-coupled receptor signaling were also down-regulated. Such reductions were not surprising in light of the anatomical observations of reduced neuronal populations in the CA1 area of the MTLE hippocampi. In contrast, many genes associated with immune and inflammatory functions showed increased expression. Among these were those associated with antigen presenting pathways, chemokine signaling, and interleukin signaling pathways (Table 3).

Table 3 List of genes up or down-regulated 1.5 fold in the 12 most significant canonical pathways according to Ingenuity Pathways Analysis

In another approach, the list of genes that were differentially expressed in MTLE were queried using Ingenuity Pathway Analysis Software to discover biological networks and pathways that may link genes and ESTs on the basis of known functional or structural relationships. Through such an analysis, several interesting observations were made in regards to the genes showing changed expression in the MTLE hippocampi.

Among the several networks of interrelated genes discerned by Ingenuity software, was one of special interest (Figure 2). In this network, several of the up-regulated genes were those associated with astrocytes (Figure 2). These included GFAP, CD44, PSEN1, VIL2, RDX, MSN, PLEC1, TNC, CSPG2, S100A6, S100B, AHNAK, and KIAA0992, CRYAB, CAPN3, CAPN2, and TTN. In addition to those in this network, there were other genes associated with astrocytes such as ARHGAP5, GSN (Gelsolin), AQP4, DMN (Desmuslin), and DTNA (Dystrobrevin alpha) which were also up-regulated, and FBXW7, which was down-regulated.

Figure 2

A pathway diagram developed using the Ingenuity Pathways Analysis and Knowledge Base in which all genes changed by more than 1.5 fold were queried. This diagram shows one of the relational networks in which all up-regulated genes (green) were those related to astrocyte structure or function. This network is very significant with a score of 40, which means that the association by chance alone is less than 99.99%. Those genes with a significant increase in fold change in our study are shaded green, while those with a decrease are shaded in red. In the diagram, the vertical diamond shape represents an enzyme, a horizontal diamond is a peptidase, downward pointing triangle is a kinase, horizontal rectangle is a nuclear receptor, the oval a transcription factor, and the round other genes. Lines with arrow head = acts on; line with short perpendicular line at end = inhibits; plain line = binds only.

Abbreviations of gene names: ADAMTS3, Disintegrin-like and metalloprotease with thrombospondin type 1 motif, 3; AHNAK, AHNAK nucleoprotein (desmoyokin); CAPN2, Calpain 2; CAPN3, Calpain 3; CAST, Calpastatin; CD44, CD44 antigen (homing function and Indian blood group system); CRYAB, Crystallin α B; CSPG2, Chondroitin sulfate proteoglycan 2 (Versican); CSPG3, Chondroitin sulfate proteoglycan 3 (neurocan); FBXW7, F-box and WD-40 domain protein 7; FGF1, Fibroblast growth factor 1 (acidic); FGFR2, Fibroblast growth factor receptor 2; G3PB2, Ras-GTPase activating protein SH3 domain-binding protein 2; GFAP, Glial fibrillary acidic protein; KIAA0992, Palladin; MET, Met proto-oncogene (hepatocyte growth factor receptor); MSN, Moesin; NET1, Neuroepithelial cell transforming gene 1; NFKBIA, Nuclear factor kappa light polypeptide gene enhancer in B-cells inhibitor, alpha; PCTK3, PCTAIRE protein kinase 3; PDE8A, Phosphodiesterase 8A; PLEC1, Plectin 1, intermediate filament binding protein 500kDa; PPFIA1, Protein tyrosine phosphatase, receptor type f, interacting protein (liprin), α 1; PSEN1, Presenilin 1 (Alzheimer disease 3); PTGS2, Prostaglandin-endoperoxide synthase 2 (G/H synthase and cyclooxygenase); RDX, Radixin; S100A6, S100 calcium binding protein A6 (calcyclin); S100B, S100 calcium binding protein, β (neural); TNC, Tenascin-C; TTN, Titin; VCAM1, Vascular cell adhesion molecule 1; VIL2, Villin 2 (ezrin); YWHAZ, Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta.

A group of immune and inflammatory response genes was also up-regulated. These genes included chemokines and their receptors (CCL2, CCL3, CCL4, CXCR4), cytokines and their receptors (FGF1, FGF2, FGF13, TNFSF7), signal transduction proteins (CALM1, CALM3, PPP3CA, PPP3R1, PTPRD, PTPRG, PTPRN, PTPRO), transcription factors (FKBP1B, FKPB1A) and transcription factor-related genes (AGT, COL1A1, COL21A1, NCAM1, VCAM1, CD44, IL11RA, IL13RA, IL15). Other genes associated with immune response were TNC, ID3, MOG, ITPKB, HLA-E, FGF1, HLA-DRB1, HLA-DRA, C3, C4, PMP2, CLESF12, CD74, and TYROBP.

In addition to the glutamate and GABA receptor signaling pathway genes that were down-regulated, some of the other genes whose expression were down-regulated also seem to be related to neuronal loss in CA1. These included, for example, synaptic vesicle cycle associated genes (34) SNAP 25, SYN1, RAB3A, RAB 6A, and RAB15. However, a few presumed synaptic vesicle associated genes were up-regulated (SYN 16, SNAP 23, HSP70, HSP40).

Validation of Candidate Genes

A small number of genes that were differentially expressed in microarray analysis were chosen for validation by QRT-PCR and comparison to microarray data. They were AQP4, CXCR4, CD44, PLEC1, GFAP, and NKKBIA. QRT-PCR confirmed the changes observed by microarray expression analysis, and suggested that microarray analysis tended to underestimate expression levels (Table 4).

Table 4 Gene Expression Fold Changes*

Immunohistochemical localization validated the expression of the corresponding proteins of GFAP, DMD, AQP4, PLEC1, CD44, and CXCR4 in tissue sections of the hippocampus (Figure 3). GFAP, AQP4, DMD, CD44, and PLEC1 proteins were localized in astrocytes, whereas CXCR4 was predominantly localized in microglia and a small population of astrocytes. The localization of these molecules showed different patterns in the sclerotic and non-sclerotic hippocampi in preliminary immunohistochemical studies. The localization of AQP4 (Figure 3A & B) and dystrophin (DMD) (Figure 3C&D) was as described previously on astrocyte end feet. In MTLE AQP 4 expression was increased on the end feet while DMD was decreased (35,36). In the non-sclerotic hippocampi, CD44 was localized in fine fiber-like extensions, which had globular profiles along their length and often ended in cone-like structures. These processes extended into the CA1 pyramidal cell body region from the stratum oriens border. Intermingled among them were immunoreactive astrocytes, whose processes also extended for 100 µm or greater distances within the pyramidal layer (Figure 3E). It was difficult, only on the basis of the light microscopic appearances of these processes to determine definitively whether they were exclusively astrocyte extensions or some were axonal processes. In the sclerotic CA1 region, CD44 immunoreactivity was seen very easily to be throughout the astrocytes and their processes and their perivascular staining strongly outlined the microvasculature (Figure 3F). The astrocytes formed an intricate network. Elongated fiber-like processes seen in the non-sclerotic specimens were not seen in sclerotic regions of CA1. PLEC1, which was only weakly expressed in astrocytes in the non-sclerotic CA1 region, was more strongly expressed in the astrocytes and their perivascular endfeet in the sclerotic hippocampus (Figure 3G&H). CXCR4 immunoreactivity, which was barely detectable in the CA1 region of the non-sclerotic hippocampus, was prominent in microglia and a small population of astrocytes in the sclerotic hippocampus (Figure 3 I&J).

Figure 3

This figure shows photomicrographs of the localization in area CA1 pyramidal layer of non-sclerotic (non-MTLE, A, C. E. G, I) and sclerotic (MTLE, B, D, F, H, J) hippocampi of some proteins whose genes showed changes in expression in this study. A & B show the immunolocalization of AQP4 (Aquaporin 4). Note the reduction in immunoreactivity in astrocytic foot processes around blood vessels in MTLE and increase in astrocytic processes throughout the pyramidal layer. C & D show reduction in DMD (Dystrophin) around blood vessels in MTLE. E & F show localization of CD44 (CD44 antigen). In the non-sclerotic hippocampus (E) immunoreactivity is seen on fiber-like processes (arrow) extending into the pyramidal layer from the stratum oriens. Scattered among these fibers are CD44 positive astrocyte cell bodies (not shown). In MTLE, many more astrocytes are immunoreactive and form a dense network throughout the region with more intense staining in perivascular foot processes (arrow in F). G & H show PLEC1 (Pectin 1, intermediate filament binding protein 500kDa) immunoreactivity, which is increased in expression in astrocyte cell bodies and perivascular foot processes (arrow). I & J show an increase in CXCR4 (Chemokine (C-X-C motif) receptor 4) expression on microglia in the sclerotic hippocampus (arrow). CXCR4 is increased on a small proportion of astrocytes as well.

Discussion

Gene Expression and Astrocytes

The identification of a set of genes showing differential expression only in the MTLE hippocampi opens the door to unraveling more fully the molecular mechanisms underlying hippocampal sclerosis. Area CA1 of the MTLE hippocampus differs in its cellular composition from the MaTLE and PTLE hippocampi in distinct ways. The MTLE hippocampi are characterized by a significant loss of neurons and an increase in astrocytes and microvasculature. It was thus not surprising to find prominent among the differentially expressed genes several related to many aspects of astrocyte function. However, the up-regulated genes are not merely a reflection of increased astrocyte numbers. Immunohistochemical localization of the proteins produced by the differentially expressed genes shows distinct differences in their cellular compartmentalization (this paper and 35). Several of the molecules observed to be up-regulated in this study have been shown previously to be expressed only in “reactive astrocytes” compared with resting astrocytes (37).

Among the astrocyte-related genes upregulated in MTLE were those associated with changes in astrocyte morphology. GFAP is associated with astrocytic process extension, Palladin (KIAA0992) regulates the shape of astrocytes (38) and the plectin gene (PLEC1), and its protein serves as a structural protein (39), the latter was also previously reported elevated in sclerotic hippocampi (40). The genes Ezrin (VIL2), Radixin (RDX), and Moesin (MSN) control three closely related proteins (ERM proteins), which constitute the machinery for the association of actin filament with the plasma membrane (41).

The amino terminal region of the ERM proteins is directly associated with the cytoplasmic domain of CD44 (42). The CD44 gene and protein were also increased in our study and localized in astrocytes. CD44 is a family of surface glycoproteins and is a receptor for hyaluronan (HA), the latter of which is also elevated in the hippocampus of patients with mesial temporal lobe epilepsy (43).

The up-regulation of the proteins of two extracellular matrix protein genes — the glycoprotein Tenascin-C (TNC) and the glycosoaminoglycan CSPG2 (Veriscan) were also previously reported in the sclerotic hippocampus (43, 44). Both of these extracellular matrix proteins are secreted by astrocytes (45).

The S100 proteins are a family of proteins involved in Ca²⁺ regulation in a variety of intracellular functions. In this study, the gene expression for two of these, S100A6 and S100B, was increased. A major and specific target protein for S100B is the giant phosphoprotein AHNAK (46) whose gene expression was up-regulated in the MTLE hippocampi in the present study.

It was reported as early as 1899 (47) that there is a proliferation of the microvasculature in the sclerotic Ammon’s horn. Our microarray studies point to a quite complex molecular reorganization of the sclerotic microvasculature. Gene expression of several molecules associated with the BBB was changed. They include AQP4, AHNAK, CD44, dystrobrevin A, and PLEC1. It is shown that human brain microvessels have chemokine receptors CCR1 and CCR2 on the parenchymal surface of the vessels (48). The genes for the ligands of these receptors, CCL3 and CCL2, were increased in expression in MTLE in the present study.

Astrocytes and microglia may play a central role in immune processes in the brain. Astrocytes are known to contribute to the inflammatory environment of the central nervous system by producing a wide range of immunologically relevant molecules (49,50). Among these are class II major histocompatibility complex antigens (49) HLA-DPA1, HLA-DQA1,HLA-DRA, HLA-DRB1, and HLA-DRB3 which were up-regulated in this study. Astrocytes also produce a variety of chemokines and cytokines and their receptors (49,50), in addition to those demonstrated in this paper. Additionally, we have demonstrated the presence of activated microglia that express CXCR4 in the sclerotic hippocampus. Such activated microglia may be an important source of TNF and IL-1β (51).

An obvious question that arises is whether such gene responses were due to some intervention other than sclerosis, such as surgical manipulation of tissue or injury due to depth electrode insertion. This appears not to be the case, as none of the sclerotic patients used in this study had depth electrodes placed in them, whereas four patients with non-sclerotic hippocampi had.

An Integrating Hypothesis

One of the functions of microarray analyses is the generation of new or alternative hypotheses to guide future research, and one such hypothesis can be proposed based on these data.

High levels of glutamate release by neurons at a seizure focus as the mechanism for the initiation and maintenance of seizures has been a popular hypothesis with much experimental data in support. In favor of such a mechanism, in vivo dialysis studies of hippocampal seizure foci in TLE patients show higher than baseline levels of extracellular glutamate and excitotoxic levels during ictus (52) as well as high cellular levels (53) in the sclerotic hippocampus compared with the non-sclerotic hippocampus. However, the source of such elevated extracellular levels of glutamate in sclerotic hippocampi have been a puzzle because the sclerotic hippocampus is conspicuously depleted of glutamatergic neurons, long thought to be the source of extracellular glutamate.

However, recent studies have shown that astrocytes can be a non-neuronal source of glutamate (54). Ca²⁺ elevations in astrocytes induce the exocytotic release of glutamate from astrocytes (55). More recently it has been shown in hippocampal slice models of epilepsy that astrocytic glutamate can generate paroxysmal depolarization shifts (PDSs) in neurons (56) suggesting that they may play an important role in interictal events related to seizures. Fellin et al. (57) confirmed the results of Tian et al. (56) in showing that astrocytes are activated during a period of epileptiform activity and their activation results in increased glutamate release. However, Fellin et al. conclude from their studies (57) that the astrocytic glutamate is not necessary for the initiation of epileptiform activity, although it may modulate and enhance the duration of ictal seizurelike events.

Whereas the experimental slice models of epilepsy discussed above provide clear support for the possibility of Ca²⁺-induced glutamate release by astrocytes, these models differ significantly from what prevails at a sclerotic hippocampal seizure focus in TLE. At such sclerotic foci, in addition to a lack of neurons, the astrocytes are “reactive” and behave differently (58). Indeed, we even found previously that astrocytes in primary cultures from human sclerotic hippocampi compared with those from non-sclerotic hippocampi show an increased oscillating frequency of Ca²⁺ waves, suggesting increased intracellular Ca²⁺ release (55) and perhaps increased glutamate release. Cortical astrocytes isolated from epileptic foci in childhood epilepsies (59) and epileptogenic cortical mass lesion in TLE (60) also exhibit elevated Ca²⁺ oscillation frequency).

The hypothesis suggested by our microarray expression data and some previously reported complementary molecular changes in sclerotic hippocampi by others (16,43) is that in the astrocytes in sclerotic foci are present the molecular mechanisms for calcium dependent glutamate release (Figure 4). In particular, that intracellular Ca²⁺ release in the sclerotic hippocampus is possible through stimulation by immune factors such as interleukin-1β and chemokines such as stromal derived factor 1α (SDF-1α). IL-1β levels are elevated in MTLE patients, and IL-1β acting through interleukin receptors can induce cyclooxygenase-2 (COX-2) and prostaglandin-E₂ (PGE₂) through the activation of the transcription factor NFκB (61). Several other molecules produced by the activation of the NFκB pathway can also directly or indirectly result in intracellular Ca²⁺ release. These include S100B (62); hyaluronan, which, acting on its receptor CD44, could trigger Ca²⁺ release through NFκB and the PGE2 pathway (63); and the chemokine stromal derived factor 1 (SDF-1), which, by binding to the CXCR4 receptor, can lead to calcium-dependent release of glutamate through a complex signaling pathway also involving PGE2 (64) (Figure 4). An increase in CXCR4 receptors on activated microglia (shown in this paper) could allow for increased SDF-1 binding, thereby inducing microglia to release TNFα, which potentiates prostaglandin-dependent Ca²⁺ activation and glutamate release. Up-regulation of CD44, S100B and CXCR4 gene expression in sclerotic tissue are demonstrated in this paper, while NFκB and hyaluronan increases have been shown in previous studies reported in the literature (16,43).

Figure 4

Hypothesized molecular pathways for glutamate release from astrocytes in sclerotic hippocampi. This diagram is constructed on the basis of molecules produced by genes whose expression was found to increase (outlined in red) in this study and those reported in other studies (outlined in blue; also see text) of human TLE. The figure also shows that many of the pathways could be activated through known initial precipitating factors of TLE. Abbreviations: CCL2, CCL3, and CCL4, CC family chemokines; CD44, CD44 antigen (homing function and Indian blood group system); CX3CL1, fractalkine; CXCR4, α-chemokine receptor; COX-2, cyclooxygenase-2; EP1,3, calcium mobilizing receptors for PGE2; ERK1,2, extracellular signal-related kinases; HA, hyaluronan; IL-1β, interleukin-1β; HSV, herpes simplex virus; HHV6, human herpes virus 6; IL-1R, interleukin 1 receptor; PGE2, prostaglandin E₂; SDF1, stromal derived factor 1, also known as CXCL12; TLR, Toll-like receptor; TNFR, tumor necrosis factor receptor; TNF-α, tumor necrosis factor α.

Ca²⁺ dependent glutamate release from astrocytes is shown to be a SNARE protein-dependent process that requires the presence of functional vesicle-associated proteins, and the glutamate being released from vesicles through an exocytotic pathway (65). Querying the gene expression levels of vesicle-associated proteins in our study reveals that SNAP23 (+ 1.57 fold) and Syntaxin-16 (+ 1.7 fold) were up-regulated. SNAP23 is an analogue of SNAP25 and is expressed in astrocytes (66). This provides evidence that astrocytes in sclerotic regions may possess vesicular mechanisms involving SNARE proteins for glutamate release.

In the context of a sclerotic hippocampal seizure focus elevated levels of glutamate, if produced by astrocytes, could initiate seizure activity in the brain by functioning as a paracrine agent, exciting neurons in the subicular region, which is in very near proximity and has an intact population of pyramidal neurons. The subiculum is the main source of output pathways from the hippocampus to the rest of the brain.

Furthermore, the immunological responsiveness of astrocytes may provide clues as to how a variety of initial precipitating factors contributing to TLE (4) may be associated with the common pathological substrate of hippocampal sclerosis. It is established that an early childhood febrile seizure can be a predisposing factor for TLE (67). Febrile seizures have been shown to be associated with increased IL-1β levels in CNS and blood plasma (68,69), and such IL-1β levels associated with a genetic polymorphism (70,71). IL-1β can act on the interleukin receptors on astrocytes and via the transcription factor NFκB induce a variety of molecules including the chemokines, cytokines and ERM proteins that are up-regulated in the sclerotic hippocampus.

Infections of the brain such as meningitis (72), Human Herpes Virus 6, and Herpes Simplex Virus (73) are shown to be associated with TLE. It is possible that the causal link in this association may also be via the astrocytes. Human astrocytes are shown to express Toll-like receptors (TLR) (74). TLRs are receptors that recognize pathogen-specific molecular patterns (75) and a wide range of TLRs is expressed by astrocytes, including those that respond to viruses and bacteria (76). Both IL-R1 and TLRs induce signal transduction pathways leading mainly to activation of the transcription factor NFκB (77) and presumably of molecular pathways that could lead to enhanced release of glutamate by the sclerotic hippocampus.

The proposed hypothesis, while providing an integrative framework, provides the basis for future experimental exploration of the pathophysiology of TLE.