Proteomic analysis of human epileptic neocortex predicts vascular and glial changes in epileptic regions

Abstract

Epilepsy is a common neurological disorder, which is not well understood at the molecular level. Exactly why some brain regions produce epileptic discharges and others do not is not known. Patients who fail to respond to antiseizure medication (refractory epilepsy) can benefit from surgical removal of brain regions to reduce seizure frequency. The tissue removed in these surgeries offers an invaluable resource to uncover the molecular and cellular basis of human epilepsy. Here, we report a proteomic study to determine whether there are common proteomic patterns in human brain regions that produce epileptic discharges. We analyzed human brain samples, as part of the Systems Biology of Epilepsy Project (SBEP). These brain pieces are in vivo electrophysiologically characterized human brain samples withdrawn from the neocortex of six patients with refractory epilepsy. This study is unique in that for each of these six patients the comparison of protein expression was made within the same patient: a more epileptic region was compared to a less epileptic brain region. The amount of epileptic activity was defined for each patient as the frequency of their interictal spikes (electric activity between seizures that is a parameter strongly linked to epilepsy). Proteins were resolved from three subcellular fractions, using a 2D differential gel electrophoresis (2D-DIGE), revealing 31 identified protein spots that changed significantly. Interestingly, glial fibrillary acidic protein (GFAP) was found to be consistently down regulated in high spiking brain tissue and showed a strong negative correlation with spike frequency. We also developed a two-step analysis method to select for protein species that changed frequently among the patients and identified these proteins. A total of 397 protein spots of interest (SOI) were clustered by protein expression patterns across all samples. These clusters were used as markers and this analysis predicted proteomic changes due to both histological differences and molecular pathways, revealed by examination of gene ontology clusters. Our experimental design and proteomic data analysis predicts novel glial changes, increased angiogenesis, and changes in cytoskeleton and neuronal projections between high and low interictal spiking regions. Quantitative histological staining of these same tissues for both the vascular and glial changes confirmed these findings, which provide new insights into the structural and functional basis of neocortical epilepsy.

About the Authors:

Gal Keren-Aviram

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Chemistry & Biochemistry, Montana State University, Bozeman, Montana, United States of America

Fabien Dachet

Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing – review & editing

Affiliations The Center for Molecular Medicine and Genetics, Wayne State University School of Medicine, Detroit, Michigan, United States of America, Department of Neurology and Rehabilitation, University of Illinois at Chicago, Chicago, Illinois, United States of America

ORCID http://orcid.org/0000-0002-5146-3449

Shruti Bagla

Roles Conceptualization, Data curation, Investigation, Methodology

Affiliation: The Center for Molecular Medicine and Genetics, Wayne State University School of Medicine, Detroit, Michigan, United States of America

Karina Balan

Roles Investigation, Methodology, Validation

Affiliation: The Center for Molecular Medicine and Genetics, Wayne State University School of Medicine, Detroit, Michigan, United States of America

Jeffrey A. Loeb

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations The Center for Molecular Medicine and Genetics, Wayne State University School of Medicine, Detroit, Michigan, United States of America, Department of Neurology and Rehabilitation, University of Illinois at Chicago, Chicago, Illinois, United States of America

Edward A. Dratz

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

Affiliation: Department of Chemistry & Biochemistry, Montana State University, Bozeman, Montana, United States of America

Introduction

Epilepsy is a neurological disorder that affects about 1% of the population worldwide. It is often called a spectrum disorder because of the large number of different etiologies or ‘brain insults’ that are associated with epilepsy. Similarly, symptomatic seizures can be clinically quite different, depending on where in the brain they arise. While seizures can often be treated symptomatically with anti-epileptic drugs, about 30% of the patients have seizures that are refractory to medications and many of these patients can benefit from a surgical procedure that removes electrically defined brain regions, while largely preserving normal regions [1,2]. The samples removed in these surgical procedures offer an opportunity to use proteomic and other high-throughput approaches on the resected human brain tissues to probe the molecular basis of human epilepsy [3]

The heterogeneity in resected tissues generated from such surgeries has made it hard to use these tissues to study the pathogenesis of epilepsy, in part due to the lack of suitable control samples. Previous studies have used specimens from non-epileptic, age and gender matched individuals as control samples [4,5], but this is far from optimal, since natural biological variation and differences in antiepileptic medication could well mask epilepsy-specific changes. We have developed a Systems Biology of Epilepsy Project to provide a new approach to study human epilepsy. Rather than studying histologically pathological tissues, the Systems Biology of Epilepsy Project focuses on human neocortical tissues that are located by in vivo recorded electrical activities, by quantifying the frequency of epileptic potentials between seizures (interictal activity). Omic and histological comparisons are made within each patient, between high- and nearby low-spiking control tissue [6]. This is a powerful experimental design, since high spike frequency in certain brain areas has been found to predict the ictal-onset locations, and resection of the interictal high spike tissue has been associated with good surgical outcomes for seizure frequency reduction [7]. In the present studies, the focus on tissues with abnormal electrophysiology allows for differentiation of proteomic changes specifically associated with epileptic activity and not due to less localized secondary responses, such as inflammation or gliosis.

Using this experimental approach, we have previously found highly consistent changes in gene transcription between regions that produce seizures [8,9] and high levels of interictal spiking [10,11], compared to nearby ‘control’ regions with little to no epileptic activities. Here, we have taken a proteomic approach, using 2-D differential in gel electrophoresis (2D-DIGE), to study six pairs of high and low spiking human neocortical tissues, to determine if there are common proteomic changes as function of abnormal epileptic activity across these patients. The 2D gels have the advantage that they directly display the relative abundances and patterns of protein isoforms that change in the diseased subjects. All of the tissues investigated here were also used to investigate transcriptional changes [10,11]. Our results show consistent changes in proteins that predict decreased glial cells and increased vascularity, as well as a hierarchy of ontological changes that could underlie the development of epilepsy. While some of these changes have been independently found through gene expression studies, many changes implicated here by proteomics are unique, providing independent and added value to our understanding of the epileptic human cortex.

Material and methods

Human tissue collection and processing

Human neocortex samples were acquired from six patients with refractory epilepsy, written informed consent was obtained from all patients prior to their surgeries through an approved IRB at Wayne State University (Detroit, MI). Only tissues not needed for their clinical care that would otherwise be discarded were used. For minors, signed consent was obtained by parents or guardians. Only tissue regions that were classified in vivo, using long-term subdural grid electrical recordings, were used after removal, as part of the patients medical treatment in a two-stage surgical procedure at the Comprehensive Epilepsy Center at Wayne State University, as described earlier [8–11]. Briefly, electrical recordings, taken over several days, were quantified for their intrinsic interictal spike frequency. Regions showing high spike frequency were compared to nearby ‘control’ regions, defined as regions showing no or the least amount of epileptic activity. After excision, the brain samples were further dissected, as described by Loeb [3] and Lipovich et al [10], to provide adequate samples for a range of analyses within the Systems Biology of Epilepsy Project, with plans for future integration of the data (as illustrated in the experimental flowchart in Fig 1). Samples underwent subcellular fractionation to produce three fractions: P1 –nuclear, P2 –membrane, and Supernatant–cytosol, as previously described by Beaumont et al [9]. The Bradford protein assay was used for protein quantification [12]. The fractionated samples were frozen and sent to Montana State University on dry ice for proteomic analysis.

[Figure omitted. See PDF.]

Fig 1. Experiment flowchart.

Pairs of tissues from low and high spiking neocortical brain regions are removed as part of the patient’s clinically-determined surgical procedure to treat refractory epilepsy. Each tissue sample is homogenized and separated by centrifugation into three subcellular fractions (Cytosolic, Nucleic and Membrane) and analyzed by 2D-DIGE. Spots found to vary were identified by LC-MS/MS. The proteomic information was used in pathway analyses.

https://doi.org/10.1371/journal.pone.0195639.g001

Protein labeling and 2D-DIGE

Frozen samples were lyophilized at Montana State University and suspended in denaturing buffer (7 M urea, 2 M thiourea, 3% (w/v) CHAPS, 1% (w/v) amidosulfobetaine-14, 10 mM Tris-HCl pH 8.4) to bring the protein concentration to 5 mg/mL. Samples were reduced with 5 mM tributylphosphine for 1 h, followed by alkylation with 10 mM 4-vinylpyridine for 2 h, all at 4°C. The fluorescent dye covalent labeling procedure was done as recommended by GE Healthcare with a few modifications. Briefly, all dyes, Cy2 (GE), Cy3 & Cy5 (synthesized by the Grieco group at MSU) were brought to 400 pMole/μL in dry dimethylformamide, and reacted with the proteins at the ratio 8 pMole dye/μg protein on ice in the dark for 30 min. The reactions were quenched with 1 μL of 100 mM lysine for 10 min. After labeling, reaction tubes were combined to contain 50 μg of Cy3 labeled proteins from a High Spike sample, 50 μg of Cy5 labeled proteins from a Low Spike sample from the same patient and 50 μg of Cy2 labeled, pooled internal standard. The pooled internal standard of each subcellular fraction contained an equal amount of protein extract from High and Low Spike samples from all patients. A technical replicate was prepared from the same unlabeled protein samples, using reciprocal Cy3 and Cy5 dye labeling between High and Low Spike samples. The Cy2, Cy3 and Cy5 labeled samples were combined and brought to 450 μL with re-hydration buffer (7 M urea, 2 M thiourea, 2% (w/v) CHAPS, a trace of bromophenol blue, and 0.5% (v/v) final concentration of 3–11 non-linear ampholines (GE Healthcare)). The solutions were loaded on 24 cm pH 3–11 non-linear IPG strips (GE Healthcare) for passive re-hydration over night at 4°C in the dark. An Ettan IPGphor II was used for iso-electric focusing, with the protocol recommended by the manufacturer. Focused IPGs were agitated for 30 min in equilibration buffer (6 M urea, 375 mM Tris-HCl pH 8.8, 20% (v/v) glycerol, 2% (w/v) SDS), loaded onto 11% SDS gels, and overlaid with agarose. The second dimension SDS-PAGE was carried out using Ettan DALT twelve System, and the protocol recommended by GE Healthcare.

Image acquisition and analysis

Fluorescent scanning of the gels was carried out using a Typhoon Trio (GE Healthcare). The scan resolution was 100μm/pixel and the manufacturer’s recommended settings were used. The image analysis was carried out using Progenesis SameSpots software (Nonlinear Dynamics, version 3.1). The three subcellular fractions were analyzed as independent data sets, each set containing 12 gels and 36 analytical gel images at each of three fluorescent colors. Standard settings were used for spot detection, background subtraction and ratiometric normalization. The initial statistical analysis utilized the ANOVA test with p-value < 0.05 and fold change (FC) ≥ 1.25 to detect protein spot expression differences between high and low spike samples of each patient. Spots that were found to have expression differences in 3 or more patients were designated “variable spots”. Ratiometrically normalized spot volumes for each spot were exported for further analysis.

In order to quantify relative concentrations of GFAP in the different patient samples, we added the intensities of all GFAP isoforms together. Using ratiometric normalization the changes in different spots cannot be related to each other, thus we carried out a second analysis, using normalization by the total spot volume on each color image. This second analysis allowed us to unite the data from all three sub cellular fractions for all subjects for analysis. The spot intensity values used for each spot were the average of the normalized volumes of all of the technical replicates.

Histochemical staining of blood vessels

To quantity the blood vessels we used 3,30-diaminobenzidine (DAB, Sigma) directly on the slides and stopped the staining as soon as a dark coloration appeared. This reaction gives a dark coloration of all the cells and allowed visualization of the blood vessels. by microscopy [11]. The length and density of the blood vessels were measured using Metamorph® 6.2 software (Molecular Devices Inc, CA, USA).

Protein identification by mass spectrometry

Preparative gels were run to resolve the proteins of each of the subcellular fractions for protein spot picking and mass spectral analysis. The preparative 2D gels contained a total of 500 μg protein, including 450 μg unlabeled fractionated epilepsy brain tissue, plus 50 μg of the Cy2 labeled internal standard, which was used for spot alignment with the analytical gels. After fluorescence imaging, preparative gels were stained with silver Coomassie [13]. Spots of interest were manually excised from the gels and In-gel digested with porcine trypsin (Promega) overnight at 37°C, as previously described [14]. Extracted peptides were analyzed, using one of two LC/MS/MS systems, previously described [14,15] with small modifications. One system was an Agilent 1100 HPLC coupled to an XCT-Ultra 6330 Ion Trap and an in-line Agilent microfluidic Chip LC (43 mm 300A C18 separation column with 40nL trapping column). Peptides were eluted with a 5–95% acetonitrile gradient over 22 min. The second system was an Agilent 1100 HPLC coupled with an 6520 Accurate-Mass Quadrupole Time-of-Flight LC/MS (Q-TOF) and an in-line microfluidic Chip LC, as described above (Agilent Technologies). Peptides were eluted with a 5–50% acetonitrile gradient over 16 min. MGF compound list files were generated, using the default algorithms in DataAnalysis 3.3 (Bruker Daltonics) and Qualitative Analysis B.04.00 (Agilent Technologies) software for XCT and Q-TOF, respectively. We used MASCOT in-house version 2.2 (Matrix Science, London, UK) to query the MS/MS spectra, using the following parameters: database: SwissProt 2012 (534242 sequences; 189454791 residues), taxonomy: human (20317 sequences), enzyme: trypsin, allowing up to 1 missed cleavage, fixed modification: pyridylethyl (C), variable modifications: deamidation (NQ) and oxidation (M). Peptide tolerance was ± 1.2 Da and 30 ppm, MS/MS tolerance ± 0.6 Da for the XCT and 30 mmu for the QTOF. Protein identification was considered significant using the default threshold of p < 0.05. For MS results see supplementary material.

Western blots

A mass of five μg of protein of the high and low spike samples from the nuclear fraction from each patient, were each labeled with Cy3, as previously described, and resolved on mini 10% SDS-PAGE gels, following the standard Bio-Rad procedure. The proteins were then blotted onto Hybond-P membrane (GE Healthcare), using an XCell II apparatus (Invitrogen), following the manufacturer’s recommended procedure. The membrane was blocked with 5% nonfat dry milk in TBST (20 mM Tris pH 7.6, 140 mM NaCl, 0.1% Tween 20) for 2 h, reacted with monoclonal mouse anti-GFAP (NeuroMab N206A/8) at 1:1000 dilution in TBST for 2 h, washed three times with TBST, reacted with Goat anti-mouse IgG Alexa 488 conjugated antibody (Invitrogen) at 1:2000 dilution in TBST for 2 h, and imaged with the Typhoon Trio (GE Healthcare). ImageJ [16] was used for quantification. Six high intensity Cy3-labeled nuclear fraction bands were used as loading controls.

Statistical and pathway analysis

We used the lme function in the R package nlme [17] to independently test for significant natural log fold-change (high/low) difference from 0 for each protein spot. A p-value < 0.05 for the two-sided t-test and FC ≥ 1.25 was deemed significant. A total of 34 spots passed the t-test cutoffs out of the 4400 measured spots. Thus, we have not corrected for multiple testing, but many of these spots were independently authenticated. GFAP was validated by Western blots, and excluding two spots with ambiguous MS identifications, all other statistically changing proteins had multiple isoforms, and each isoform independently passed the t-test cutoff or very close to it. Thus, the multiple independent identifications of the proteins in different isoforms strengthened the notion that these protein changes were not found by chance. By implementing the mixed model mechanism of lme, we accounted for the variance structure of technical repetitions that were nested within the values from the individual patients, which were imposed by the reciprocal dye-labeling design of the paired experiments.

Spot of interest (SOI) cluster analysis was carried out using Ward's Minimum Variance Method to perform hierarchical clustering in R [18]. Pairwise dissimilarities (distances) between observations were computed with Cluster Package [19], using the "Daisy" function and Gower's Metric [20].

Data were analyzed through the use of IPA (Ingenuity^® Systems, www.ingenuity.com, version 42012434, release date 2017-12-07). The IPA core analysis was applied to the list of SOI with fold change values ≥ 1.3 FC cutoff and default parameters. The enrichment within the nervous system was tested using Fisher’s exact test and predicted the involved pathways for each patient. A comparison analysis across the patient population was used to identify pathways for which involvement was frequently predicted (–log(p-value of pathway) ≥ 1.3, in at least 3 patients).

Results

The human epileptic neocortex proteome

Tissues from patients who underwent cortical resections for intractable seizures, consenting to be part of our System Biology of Epilepsy Project, were used for this proteomic study. The patients undergo a two stage procedure that begins with several days of subdural electroencephalography recording, followed by excision of the epileptic foci that often contain nearby, electrically quieter “normal” surrounding tissue, as part of their surgical resections. Only tissues that were part of the planned surgical procedure and resected by the neurosurgeon were used. Pairs of 1 cm³ blocks of tissue precisely localized to high and low spiking regions from each patient were further subdivided for histology, genomics, and proteomics, as described [3].

Differential proteomics analysis was carried out after subcellular fractionation and differential fluorescent labeling of High and Low Spiking sample pairs from each of 6 patients (Table 1), using the 2D DIGE method to provide improved sample comparison, protein resolution, sensitivity and increased dynamic range [21]. In total, more than 4,400 protein spots were analyzed in independent data sets from the three subcellular fractions. In a t-test (p<0.05, FC > 1.25), 34 spots were found to change between High and Low Spiking samples. 31 of these spots were identified by mass spectroscopy as originating from 18 gene products, some of which existed in multiple protein isoforms. In high spiking regions eight gene products were up regulated (SNCA, STMN1, UGP2, DSP, CA1, PRDX2, SYN2 and DPYSL2) and ten were down regulated (GFAP, HNRNPK, CPNE6, CRYAB, GNAO1, PHYHIP, HNRPDL, ALDH2, GAPDH and LASP1) (S1 Table). Given the wide range of spike frequencies and potential variance in absolute protein expression between patients, we also hypothesized that spots that change between high and low spiking regions within individual patients, but have a high variance between patients (“variable spots”) may be biologically significant and might indicate involvement in an active process or molecular pathway that is changing, but would not meet the criteria to pass a simple t-test for all the subjects. Therefore, we also selected these “variable” spots (FC > 1.25, which were up or down regulated in the high spiking tissue, in at least half the patients) in each sub-cellular fraction. A data set of 397 SOI resulted, including statistically changing (t-test) and “variable” spots from all three fractions. 2D maps of SOIs shown in S1 Fig. These SOI were further analyzed, using clustering methods.

[Figure omitted. See PDF.]

Table 1. Patient clinical information and tissue samples used for this study.

https://doi.org/10.1371/journal.pone.0195639.t001

In cases where a protein displayed a spot train of isoforms, in which each spot had a similar up or down regulated expression pattern, only the most intense isoform spots were picked for mass spectrometry-based protein identification, and the other isoforms in the similarly regulated spot train were assumed to be the same protein. Some SOIs were very low in abundance, while some were not ubiquitously expressed, leading to their apparent absence in the preparative gels. In total, more than 90% of the 397 SOIs were identified as arising from 146 gene products (S1 Table). We did not attempt to characterize the numerous isoform modifications present in the proteins, at this stage.

Predicted histological changes through hierarchical clustering of protein spot expression patterns

In human studies, biological variation is typically conceived of as an inherent limitation which might mask disease or treatment affects. However, in this part of our study, where our focus was to find proteins which change their level of expression in the same direction relative to the internal standard across all the samples, biological variation has been beneficial. The use of a pooled internal standard control, containing proteins from all samples, facilitated the comparison of spot expression levels across the whole dataset. To provide more expression data points and higher pattern specificity (similarities in relative expression), samples were treated independently, regardless of patient and spike frequency, and variation in expression between and within patients were weighted equally. The added variation revealed unique protein expression patterns and enabled better separation of co-expressing protein groups (Fig 2). The SOIs were clustered and visualized on a dendrogram, according to how the level of their expression varied across the 12 samples (6 patients x 2 spike frequencies). SOIs that showed similar levels of variation for all samples were clustered together (Fig 2A and 2B). Some spots clustered more tightly together, sharing short vertical branches on the dendrogram. These tight clustering observations were also validated by Leave One Out Cross Validation.

[Figure omitted. See PDF.]

Fig 2. Hierarchical clustering by protein spot expression patterns.

(A) Spots of interest (SOIs) from all three subcellular fractions were clustered by expression patterns (level and direction of variation relative to the internal standard) across all samples (dendrogram) and then aligned with the heat map based on their fold changes in SOI (High/Low spiking) by patients. (B) Examples of different expression patterns, which represent the colored spots in the 2D gel image to the right (Fig 2C), grouped together into two distinct clusters. Each of the lines connects a normalized spot volume of technical replicates of high and low spike frequency regions of all 6 patients (the high and low patient sample identifications are in the colored boxes at the top of Fig 2B). (C) 2D gel location of soluble fraction erythrocyte protein spots (red) and extracellular proteins (green) represented in clusters * and ** respectively in Fig 2A. Note: horizontal spot trains are largely due to several isoforms of single proteins. (D) Representative histochemical staining comparing the blood vessel density in high and low spike frequency samples from the same patient demonstrate the predicted increase in vascularity. The scale bar in D indicates 0.5 mm. For detailed method and analysis see [11]. (E) Box and Whisker plots [22] of the changes in vascularity predicted by the ratio (High/Low spiking) of 23 protein spots in the erythrocyte cluster * (in Fig 2A) correlated with blood vessel length ratio (r = 0.75) The ep122 high spiking tissue sample was not available for staining, leading to the exclusion of patient ep122 from this analysis.

https://doi.org/10.1371/journal.pone.0195639.g002

It is important to emphasize that the similarities in expression pattern of a group of spots, across patients, indicates that these changes should be treated as a single finding, allowing for reduction in complexity of the dataset. For example, groups of similarly expressed isoforms, of a single protein, which are changing in expression across the samples, will have essentially the same expression patterns (Fig 2B and 2C). Another example is a protein complex that is expressed at different levels in different samples, since the subunits of the protein complex will have the same expression patterns (change in the same direction). It is likely that the tightness of the cluster will decrease with increasing numbers of gene products in the complex, due to multiple possible regulations and involvement of proteins which are ubiquitously expressed. However, due to the high biological variation between individuals, and the resulting enhanced clustering power of this approach, we were able to apply the same logic to identify clusters that we hypothesize represent changes in the phenotype of a given cell type, within the neocortical tissue. Particular cell types were predicted, based on known cell-specific protein spots within a given group, that were identified using mass spectrometry. Detailed information on protein identification by fraction and spot is given in S3 Table, in the SOIs cluster assignment in S2 Table, and on MS identification of SOIs by Mascot, as provided in S4 Table.

Several histological changes were predicted from the groups of proteins found by hierarchical clustering. Following the right most branch of the four main dendrogram branches in Fig 2A (see S2 Fig for enlarged dendrogram and cluster number designations), we found several indicators of blood, including platelets—indicated by a fibrinogen cluster (cluster 37), that includes multiple isoforms of two of the fibrinogen subunits (FGB, FGG). Evidence for fibroblasts, smooth muscle or possibly pericytes, was indicated in this branch by co-expression of TPM1, TPM2 and TAGLN (cluster 34). The presence of plasma in clusters 28 & 33 was indicated by albumin and components of high-density lipoprotein (APOA1, APOD). Finally, the presence of erythrocytes was indicated in cluster 27 primarily by hemoglobin, but also by a cluster of proteins, all of which are present in high abundance in red blood cells (HBA, HBB, catalase, CA1, BLVRB) [23–25]. Since all tissue sample surfaces were rinsed to remove blood before homogenization, we attributed changes in blood levels to changes in vascularity in the tissues. The relative change in vascularity in the brain tissue, predicted by the erythrocyte cluster, was validated by histochemical staining (Fig 2D and 2E).

The next dendrogram branch, second from the right in Fig 2A (cluster number 24 in S2 Fig), contains mostly albumin isoforms, as well as other extracellular proteins. The distinct protein expression patterns between the extracellular and erythrocyte proteins (Fig 2B and cluster numbers 24 and 27 in S2 Fig) indicate that the changes in albumin levels are not solely due to changes in vascularity (circulating serum albumin), and suggested the presence of additional albumin in the extracellular space.

Many of the identified SOIs were related to cytoskeleton proteins, suggesting structural or migrational changes in affected cells. Dominating this group of proteins, with statistically significant down regulation in the high spiking regions in all fractions for all patients was GFAP (~50kDa): a commonly used astrocytic marker. As one of the major astrocytic intermediate filaments, the cytoskeletal GFAP mostly pellets in the nuclear fraction. The total nuclear fraction GFAP SOIs (the sum of 6 spots of ~50 kDa + 5 spots of 38–40 kDa) normalized spot volumes (Fig 3A), confirmed the decrease of GFAP in the high spiking regions (paired t-test p-value 0.002, average FC 1.53, Fig 3B). Western blot analysis validated the decrease of total nuclear GFAP in high spiking tissues (p-value 0.047, Fig 3C). The spot group of lower molecular weight (LMW, 38–40 kDa) GFAP presented a different trend of expression across the patients. These LMW GFAP spots were mostly up-regulated in high spiking regions and notably clustered with blood proteins, as shown in cluster 32–33. The six 50 kDa GFAP isoforms in the P1 (nuclear fraction) added together demonstrated a stronger down regulation in the high spiking regions (p-value 0.002, FC 1.79), if the LMW GFAPs are excluded from the analysis. Interestingly, the increase in total GFAP fold change (L/H spike) of the 50kDa isoforms was strongly correlated with increased spike frequency (frequency difference High-Low = Δspike), as shown in Fig 3D. The linear correlation coefficient was r = 0.96 with a p-value 0.002, using the Pearson product-moment correlation coefficient. Alternatively, we evaluated the non-parametric association and found a Kendall tau coefficient of r = 1, p-value = 0.003 for the data plotted in Fig 3D.

[Figure omitted. See PDF.]

Fig 3. GFAP quantification studies.

(A) 11 GFAP SOIs are the blue marked spots in the 2D SDS-PAGE image of the P1 fraction, which were used for quantification of total GFAP. (B) Total GFAP as the sum of total image intensity normalized volume (non-ratiometric) of nuclear GFAP spots shown by patient and spike frequency (L-low, H-High). GFAP is down regulated in high spiking samples (paired t-test for all 11 spots or upper higher MW (~50kDa) 6 spots, p-value 0.002, average FC 1.53 or 1.79, respectively). (C) Western blot analysis for all nuclear fraction samples with a monoclonal anti-GFAPα antibody, validated the down regulation of GFAP in high spiking samples (p-value 0.047). (D) Correlation of fold change of nuclear fraction 50kDa GFAP (upper 6 spots) (L/H spike ratios, circles are sized by the weighted total intensity) with differences in spike frequency for individual patients, Δspike = H-L.Y = 200X – 251, r = 0.96, p-value 0.002. (E) Sum normalized volume of cytosolic GFAP spots by patient and spike frequency. (F) Cytosol to nuclear total GFAP ratio (11 spots). Note—for GFAP quantification, fluorescence spot intensities were normalized to total fluorescent protein signals in each gel image. Values are given as the sum of relevant spot normalized volumes in millions of volume units.

https://doi.org/10.1371/journal.pone.0195639.g003

The total GFAP changed between high and low spike regions only in the nuclear pellet fraction and did not significantly change in the P2 and Cytosol fractions. The total GFAP in the Cytosol fraction (Fig 3E), although still mostly higher in low spiking regions, demonstrated a very different expression pattern between the patients, when compared with the nuclear fraction (Fig 3B & 3E). The levels of cytosolic GFAP were more in agreement with the level of gliosis, which was assigned to the patients by the pathological analysis (Table 1), but was not a good predictor for astrogliosis status since much of the gliosis seen was subcortical in white matter, not in the gray matter where our proteomic analysis was performed. The GFAP soluble pool (in the cytosolic fraction) is presumably related to the level of GFAP filaments (nuclear fraction) and might be affected by pathological conditions [26]. Therefore, we evaluated the ratio of total cytosolic to total nuclear GFAP concentration and found the ratio to be a good predictor for the patient’s astrogliosis level, as can be seen by comparing the patient histopathology in Table 1 to the data in Fig 3F. This analysis found that tissues with no apparent histopathology had a soluble GFAP concentration between one-quarter to one-fifth of the concentration in the nuclear fraction. Diffuse and mild levels of gliosis were associated with increased shift to soluble GFAP, until the soluble and insoluble GFAP concentrations were about equal for the patient with the most pronounced gliosis.

Bioinformatic analysis of proteins of interest

The 397 SOIs originate from 146 gene products. Multiple spots for the same gene product might arise from alternative splicing and/or posttranslational modification, as well as a SOI occurring in more than one subcellular fraction. In this exploratory phase of the project, we have identified candidate proteins for further study and did not investigate specific modifications of gene products.

Description of the 146 proteins of interest is given in S1 Table and their FCs are shown with color codes by patient. Since multiple isoforms were often summed to a single gene product, the protein fold change represents the average expression of sets of isoforms for each patient. Therefore, the absence of a significant, net fold change in at least three patients of a given protein with multiple SOIs is likely to indicate approximate cancelation of FC by differential post-translational modifications. An exception was made with regard to analysis of GFAP, which had two very distinct groups of multiple isoforms. First, the spot train having the expected MW (~50kDa), including the expected pI for full length GFAP. Second, several isoforms, identified as GFAP by MS, are modified with lower than the predicted MW and different levels of acidic shift in isoelectric point. These were analyzed separately as two GFAP products, GFAP ~50kDa and GFAP LMW (low molecular weight).

Ingenuity Pathway Analysis (IPA) was used to predict involvement of Canonical Pathways for each patient. The differential expression and biological involvement of the erythrocyte proteins in relation to changes in vascularity, was already addressed. Therefore, erythrocyte proteins were removed from the analyzed genes (unless the proteins were also expressed outside the erythrocyte cluster). After considering the differences in expression patterns of the blood components, as well as their size and likelihood of leaking out of the vessels under pathological conditions, the rest of the blood components were included in the analysis, since they might have interactions or cellular effects exterior to the blood vessels. Common involved pathways across the patient population were found using IPA comparison analysis. Pathways predicted to be involved in at least three patients were considered. The top 25 ranking canonical pathways involved in human neocortical epilepsy, and the relative contribution of each patient to the score, ranked by p value, are shown in Fig 4.

[Figure omitted. See PDF.]

Fig 4. Top 25 canonical pathways involved in epileptic brain regions are predicted by Ingenuity Pathway Analysis ranked by significance.

The sum of the logarithm of the p-value for all the patients and for each canonical pathway found is ranked from the most significant to the less significant. The contribution that each patient has on the sum is represented by color codes.

https://doi.org/10.1371/journal.pone.0195639.g004

Discussion

Surgical resection of electrophysiologically characterized brain tissue from patients with refractory partial epilepsy can be highly beneficial for reduction or elimination of seizures [27,28], and provides a unique opportunity to study the molecular and histological changes associated with abnormal electrophysiological properties of epileptic human neocortical brain regions. The Systems Biology of Epilepsy Project is taking advantage of this opportunity to generate a clinical and multi-omics database to further our knowledge of the epileptogenic process and advance the search for better treatments [3,6]. In the present study we used differential fluorescent labeling in 2D-DIGE to resolve and compare protein expression and protein modifications between high and low spike frequency tissue sample pairs from each patient. In the initial differential analysis we identified a set of proteins of interest, suggesting that these proteins might be involved in the cellular changes that lead to epilepsy. Subsequent analysis focused on attempts to relate this set of proteins to epilepsy-related biological changes. We took advantage of the natural biological variation between patients as well as natural, or epilepsy induced variation within the samples of each patient, using hierarchical clustering, to allow the separation of groups of co-expressing protein spots that predicted histological changes.

High interictal spike frequency is associated with vascular modification

Many of the changing proteins, when clustered by their expression patterns, suggest significant increases in blood vessels in high spiking neocortex. Although all but one (CA1) of the individual protein changes did not quite reach statistical significance alone, the group of proteins in the blood cluster suggested increased vascularity in high spiking samples in four of the six patients. Patient ep122 showed no change and patient ep158 showed an opposing trend, for this cluster, as well as in most of the changes in the dataset. Since patient ep158 was only one year old, his brain is in a very different developmental stage than the other patients, which might be expected to affect many of its “normal” and “epileptic” features.

The process of angiogenesis is associated with brain injury [29,30] and has previously been described in temporal lobe epilepsy [31], with an association to increased blood brain barrier (BBB) permeability [32,33]. Here, for the first time using precisely mapped cortical regions with more and less epileptic activity, we predicted increased vascularity by both proteomic and previous transcriptomic and metabolomic studies [11,34] and these findings were validated using histochemical staining [11]. Together, this work suggests that angiogenesis is closely associated with interictal high spiking neocortical tissues and may play a role in metabolic changes linked to increased cortical activity. Given the correlative nature of these human studies, we cannot say whether the vascular changes cause the epileptic activity or result from increased activity. The consistent finding with all of our transcriptomic, proteomic and metaobolomic studies, however, suggests a clear correlation between epileptic spiking and increased vascularity.

These findings also demonstrate the power of protein expression clustering to interpret proteomic studies. Although carbonic anhydrase 1 (CA1) was statistically up-regulated in high spiking regions, its correlation with hemoglobin levels and other proteins highly abundant in erythrocytes, led us to interpret CA1 up-regulation and the erythrocyte protein cluster as part of increased vascularity. You et al [35] suggested four proteins (hemoglobin, catalase, peroxiredoxin and carbonic anhydrase I) as markers for blood contamination in proteomic analysis of cerebrospinal fluid. It is important to note that, because we use a paired analysis from high and low spiking tissues that were resected in the same manner, the increase in erythrocyte proteins were not due to blood contamination of the tissues, but rather to the observed increased in vascularity in the high spiking regions. We also suggest that the similar expression pattern of most of these proteins, with the addition of biliverdin reductase B (BLVRB), defines them as erythrocytes markers. It is important in gel-based proteomic studies, to use several blood marker proteins since hemoglobin isoforms run close to the electrophoretic front and may be excluded from the analysis. In a gel-free proteomic study, Visanji et al [36] also found increased expression of catalase and CA1 in human epileptic brain, compared with controls, and localized both proteins to astrocytes. In addition to CA1 and catalase, they found seven hemoglobin isoforms, albumin, and alpha-2-macroglobulin to be up-regulated in epileptic brain [36]. They suggested catalase and CA1 had epileptogenic roles in response to oxidative stress (catalase) and support to carbonic anhydrase 2 pro-apoptotic/epileptogenic role (CA1). Although these mechanisms cannot be ruled out, our proposed erythrocyte protein cluster marker and correlation with histological measurements of increased vascularity provides an alternative explanation. We suggest that the protein changes found by Visanji et al [36] are fully consistent with increased vascularity in the epileptic brain.

The astrocytic foot is a structural component of brain vasculature, which can explain the localization of catalase and CA1 to astrocytes. The hemoglobin isoforms, catalase and CA1 are the most abundant proteins in erythrocytes [25]. Furthermore, CA1 expression is limited to erythrocytes and the gastrointestinal tract [37]. Additionally, alpha-2-macroglobulin and albumin are two major carriers of proteins in the serum [38,39]. Taken together, our alternative explanation of increased vascularity in the epileptic tissues, accounts for changes in levels of eleven of the twelve proteins found to change by Visanji et al [36].

High interictal spike frequency is associated with decreased GFAPα

Reactive astrocytes over-expressing GFAP (gliosis) is a hallmark of posttraumatic brain injury and can be see in hippocampal sclerosis, associated with epilepsy [40]. GFAP is the main intermediate filament in astrocytes, its expression increases as a function of age and is affected by multiple pathological conditions [41]. GFAP is normally found primarily in the water-insoluble subcellular pellet fraction [42], as a polymer, with a small amount in the soluble protein fraction. However, astrogliosis increases the proportion of soluble GFAP [26]. We evaluated the relative GFAP concentration in both the cytosolic and nuclear fractions and found that the ratio of soluble to insoluble GFAP correlated with levels of gliosis from the pathological assessment of the patients (Table 1), regardless of spike frequency.

Surprisingly, the 50kDa GFAP, as well as total nuclear GFAP, was down regulated in high spiking tissue of all the patients, compared with low spiking tissue from the same patient. Furthermore, the fold change decrease in GFAP was strongly related to the level of spiking, suggesting a causal relationship between the two phenomena. In addition to the unexpected decrease in GFAP, α-B-crystallin (CRYAB), another protein whose over-expression is a marker of epileptic foci [43], was found in our study to be down regulated in the P2 subcellular fraction of high spiking tissue. CRYAB is small heat shock protein, known to be strongly induced in reactive astrocytes [44]. CRYAB has a role in the anti-inflammatory response, and an anti-apoptotic role in astrocytes [45]. CRYAB is also known to regulate GFAP assembly and to co-aggregate with it in astrocytic inclusions [46]. It is unclear whether the changes in expression of the two proteins are related. It is possible that the co-reduction is due to a CRYAB & GFAP complex, to a reduction in reactive astrocytes, or to astrocytic apoptosis that may be mediated by the lack of protective CRYAB. The decrease in these two epilepsy foci markers in the high spike tissue, suggests that the cellular changes associated with abnormal electrophysiology are separate from the less localized cellular changes, such as astrogliosis, which are likely to be associated with general trauma or stress. One possible explanation is that more severely damaged cortex with more gliosis is less able to generate epileptic discharges than nearby areas without extensive gliosis and more normal neuronal connectivity.

Several GFAP splice variants have been described in the literature [41], some of which have been found to be differentially expressed in epileptic lesions [47,48]. However, our MS analysis of 2D gel electrophoresis spots, in agreement with 1D western analysis using a specific GFAPα monoclonal antibody, indicates that the canonical 50kDa GFAPα is the variant in the SOIs examined here. While the 50kDa GFAPα was down regulated in high spiking tissues, the insoluble low molecular weight (38–40 kDa), modified GFAPα products had a mixed expression pattern and clustered with the blood proteins. Although the exact nature of this GFAP cleavage process was not studied here, the co-expression with the blood proteins might indicate that insoluble 38–40 kDa GFAPα is an indicator for perivascular astrocytes, or specific GFAP modifications that are associated with the astrocytic foot process, which is an integral component of the blood brain barrier [49].

An increasing body of evidence implicating astrocytic dysfunction in epilepsy has been described [50–52]. Witcher & Ellis [53] addressed a paradox in temporal lobe epilepsy, of increased GFAP and decreased glutamate transporter EAAT2, which is predominantly expressed in astroglial cells and throughout the astrocytic membrane. They consolidated these co-occurring changes with the morphological interpretation of “…decrease in non-GFAP containing perisynaptic astrocytic processes…” Indeed, GFAP might not be an ideal marker for non-reactive astrocytes, since it is not expressed throughout the astrocyte cytoplasm and is not detectable in fine astrocytic processes [54]. In epilepsy models, expression studies of key astrocytic proteins such as aquaporin-4 and inwardly rectifying potassium channel Kir4.1, that are suspected to be epileptogenic, revealed reduction in their immunoreactivity, specifically in cortical (and not hippocampus) astrocytic processes [55,56]. Rakhade & Loeb [57] found that the glutamate transporter EAAT2 was down-regulated in human neocortical epileptic foci in experimental settings comparing high and low spiking regions, similar to our study. Thus, the decrease in GFAP in high spiking regions in our study supports a hypothesis that GFAP positive astrocytes are decreased. In addition, ALDH1L1, an alternative candidate for an astrocytic marker [54], also exhibits a decrease in our dataset (p-value = 0.013, FC = 1.18). One of the challenges of this type of study is that it is difficult to deduce whether a particular change in epileptic brain regions contributes towards hyperexcitability or may be compensatory or protective. Since all of the patients in this study failed medical management (were pharmacoresistant), one can argue that if compensatory, the changes were still not sufficient to prevent seizures. Additional studies are needed to more fully characterize the nature of the changes in GFAP positive astrocytes, whether they are due to a reduction in the number of cells or to a change in their morphology and whether differences in their expression are potentially protective against epileptic activities.

Canonical pathways involved in refractory neocortical epilepsy point to changes in the blood brain barrier and synaptic plasticity

Knowledge-based software, which predicts the involved pathways from lists of changing proteins, is a powerful tool in the complex field of neuroproteomics [58]. We queried Ingenuity Knowledge-Base Ingenuity Pathway Analysis (IPA) to predict the pathways involved for each of the patients and to identify the common, probable pathways that are associated with biological changes in the high spike frequency tissues.

Many of the identified common pathways partially overlap at the molecular level and share some of the same protein networks. For example, we have identified changes in several protein components of the adherens junction (AJ), which pointed to several of the proposed signaling pathways. IPA identified “Remodeling of epithelial adherens junctions” and “Epithelial adherens junction signaling”, as shown in Fig 4, that can also indicate changes in AJs of endothelial cells, and in our study appears to indicate changes involved in blood vessels or BBB. This notion is supported by other signaling pathways identified by Ingenuity analysis in our study, “germ cell–sertoli cell junctions signaling” and “sertoli cell–sertoli cell junction signaling”. Sertoli cells are obviously not part of this neurological system, but the blood testis barrier and the BBB share some structural and protein component similarities. For example, components of tight junctions, AJs and gap junctions, which together generate the non-fenestrated barrier between compartments [59,60]. AJ changes in the BBB can affect its permeability [61] and are also related to endothelial cells motility and angiogenesis [62,63], re-enforcing our focus on increased vascularity and BBB dysfunction, as was also suggested by the hierarchical clustering of protein expression patterns. An alternative interpretation to AJ changes in the nervous system can be due to re-modeling of synaptic junctions, where the synapses are held together by AJ’s [64], and could change as part of epilepsy related plasticity changes.

“Acute phase response signaling” was another pathway suggested by our study. Although this pathway might have relevance to epilepsy, our study, due to its experimental design, lacks the potential of finding systemic changes (in the blood), such as the acute phase response. However, this pathway, which is represented by a group of modulated plasma proteins, might be mistakenly indicated in cases where there is increased vascularity or BBB dysfunction. Several studies address dysfunctional BBB in epilepsy [65–69]. Increase in extracellular space albumin is assumed to be the result of a leaky BBB [69,70]. Therefore, extracellular proteins changes (clusters 24–26) might be interpreted as further evidence to a leaky BBB in human epileptic neocortex.

“Semaphorin signaling in neurons” and “axonal guidance signaling” found in our Ingenuity analysis, have a high degree of component overlap and suggest fundamental changes in synaptic plasticity in human epileptic neocortex. In a meta-analysis of refractory epilepsy transcription, Mirza et al [71] also found these pathways to be up-regulated. We have identified several proteins involved in these pathways, including multiple isoforms of the collapsin response mediator protein (CRMP) family. The CRMP family is homologous to dihydropyrimidinase, an enzyme responsible for uracil and thymine degradation, consistent with the identification of these degradation pathways in our analysis. Of potentially great interest is CRMP2 (formally, dihydropyrimidinase-like protein-2 or DPYSL2), which we found to be up-regulated in high spiking tissues. CRMP2 is the most abundant and studied member of the CRMPs. CRMP2 is known to have multiple neurological roles, and to be involved in several neuropathological conditions [72]. The CRMP2 antagonist, Lacosamide, is specifically used to treat refractory partial epilepsy seizures, by mechanisms that might involve stabilization of neuronal firing through regulation of voltage-gated sodium channels, or inhibition of abnormal axon sprouting [73–75].

Biological significance of statistically changing proteins

Understanding the biological significance of the change in a single protein in a proteomics study can be challenging. Therefore, we combined the statistically changing protein spots and the commonly changing spots ('variable spots') for hierarchical clustering and pathway analyses. In addition to these histological modifications and cellular pathways involved in high spiking tissues there were a number of statistically significant individual proteins that warrant further discussion. These include the guanine nucleotide-binding protein G(o) subunit alpha (GNAO1), and LIM and SH3 protein 1 (LASP1). Both proteins are involved in multiple cellular pathways. Our study did not elucidate a single pathway which clearly indicates the biological significance of the down regulation of GNAO1 and LASP1 observed in high spiking tissue. However, we demonstrated that the expression patterns of GNAO1 in the P2 fraction and LASP1 in the P1 fraction have parallel patterns of expression as the insoluble 50 kDa GFAP and CRYAB isoforms (see clusters 1 & 13 respectively in S1 Table). Since GFAP and CRYAB changes were attributed to changes in GFAP positive astrocytes, the changes in GNAO1 and LASP1 expressions may also be of astrocytic origin.

Two neuronal phosphoproteins, alpha synuclein (SNCA) and synapsin II (SYN2), were up-regulated in high spiking tissues. Both proteins might be involved in multiple pathways and processes. Neuritogenesis and synaptic transmission are of special interest, since both of these proteins localize to the presynaptic terminal and interact with cytoskeletal proteins in the nerve growth cone [76,77]. SYN2 has been implicated previously in epilepsy, since certain genetic variants of this protein contribute to the predisposition to epilepsy and its differential effect on inhibitory and excitatory synapses [76]. Another phosphoprotein stathmin1 (STMN1), which is involved in microtubule regulation, was found here to be up-regulated in high spiking tissues. Zhao et al [78] found over-expression of STMN1 in the neocortex of patients with intractable temporal lobe epilepsy, when compared with neocortex of normal temporal lobes. They used immunoassays to localize the changes in STMN1 expression to neurons and suggested that STMN1 participates in neurite sprouting and synapse remodeling. Taken together, our results suggest a heightened level of neurite sprouting and synapse remodeling in the high spiking brain areas that could be linked to the development of hypersynchrony. Our recent studies support this interpretation, showing quantitative, superficial layer-specific increases in presynaptic nerve terminals in epileptic brain regions [9], however, at the same time in another study, we discovered focal reductions of synaptic complexity in deeper layers of high spiking regions called microlesions [11].

Remarkably, transcriptional profiling of the exact same brain regions that we used for the proteomic analysis described here led us to overlapping as well as completely different pathways [9–11]. In both analyses, an increase in vasculature was in common, however the most salient transcriptional pathways centered around MAPK/CREB signaling, which did not come up in the proteomic analysis. One likely explanation is that our proteomic approach using 2D gels works best for higher abundance structural proteins, such as those in astrocytes and blood vessels, rather than lower abundance signaling and transcription factors that may be detected in transcriptional profiles. Given the complexity of the 6-layered neocortex and possible alterations in epilepsy, further studies will be needed to understand how and where posttranscriptional and posttranslational modifications are linked to changes in the cytoarchitecture, that produces epileptic brain regions.

Conclusions

This study is the first to correlate proteomic changes directly with electrical brain activities recorded in vivo, allowing for the unbiased identification of protein changes specifically associated with human epileptic activity. We combined proteomic, statistical and bioinformatic tools to predict and then validate histological changes and cellular pathways associated with the abnormal electrophysiology. We found a decrease in astrocyte markers and an increase in vascularity in the high spiking regions to be a common histological change, that we also previously found in both transcriptomic and metabolomic studies [11,34]. Furthermore, the increase in spike frequency in different patients was strongly related to the level of decreased astrocyte marker GFAP, suggesting causality between astrocytic reduction and abnormal electrophysiology.

Contrary to what is thought for forms of hippocampal epilepsy, where there is a diffuse increase in astrogliosis, our study suggests brain regions with increased astrocytes actually have less epileptic firing. The strong correlation of the reduction of reactive astrocytosis (marked by GFAP) in higher spiking regions, suggests that these reactive astrocytes may protect the neocortex from epileptic discharges, rather than cause them. While this type of analysis has not yet been done for the hippocampus, our findings suggest that local regions of the hippocampus with extensive gliosis may not in fact be more epileptically active as nearby hippocampal areas with less gliosis.

The ongoing integration of the proteomic data with the clinical, histological and other high-throughput studies of Systems Biology of Epilepsy Project is leading to the discovery of cellular and molecular interactomes of human epileptic neocortex. Such an interactome can deepen our understanding and derive new, knowledge based, hypotheses for epileptogenic mechanisms and therefore possibly lead to new therapeutic opportunities.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. 2D maps of all spots of interest (SOIs).

Each image represents 2D-PAGE of proteins from subcellular fractions using pH range of 3–11 (non-linear) and MW range 200-15kDa. SOIs are marked by blue perimeter and spot number. Note—the SOIs numbers in the cytosolic fraction are split across two images (CytA & CytB) in order to allow their visualization. P1 is the nuclear fraction and P2 is the membrane fraction. CytA and CytB are identical images of the cytosolic fraction used to accommodate different SOI numbers.

https://doi.org/10.1371/journal.pone.0195639.s001

(TIF)

S2 Fig. Spots of interest (SOIs) dendrogram.

Visualization of hierarchical clustering of 397 SOIs by their expression patterns. The dendrogram was divided into two main branches (A & B) and an enlarged image of each branch is shown. An arbitrary distance cutoff was used to separate the dataset to 37 groups of spots.

https://doi.org/10.1371/journal.pone.0195639.s002

(TIF)

S1 Table. 146 Proteins of interest, representing 397 protein isoform SOIs.

Notes: Number indicates the number of isoforms (spots) of the identified protein, using bin-averages for fold change values. (‡)—indicate a spot with an ambiguous MS identification (see supplementary material for details and definitions). (†)–indicates a spot where the observed MW does not agree with expected MW, implying a gene product modification. ↑ and ↓ indicate a spot with statistically significant up and down regulation, respectively. t-test p-value <0.05 (a) and < 0.005 (b). Patient FC color: dark green ≤ (-1.5), light green ≤ (-1.25), (-1.25) < white < 1.25, orange ≥ 1.25, red ≥ 1.5.

https://doi.org/10.1371/journal.pone.0195639.s003

(XLSX)

S2 Table. Spots of interest by cluster.

SOIs ordered by their dendrogram position. Include the protein assignment and the spot expression ratios (H/L) for each patient. Four levels of identification are used and each one coded by a label: i) '$': identification by MS of a single protein with multiple peptides or the spot has been identified by Western blot. ii) '*': identification of the spot inferred from MS and/or 2D-PAGE information indicate only one potential protein. iii) '/': identification of the spot inferred from MS and/or 2D-PAGE information indicate more than one potential protein. iv) the protein in the spot cannot be identified or predicted by our techniques.

https://doi.org/10.1371/journal.pone.0195639.s004

(XLSX)

S3 Table. Spots of interest by fraction and spot number.

Spots are ordered by fractions and spot number, including protein assignment and description. Four levels of identification are used and each one coded by a label: i) '$': identification by MS of a single protein with multiple peptides or the spot has been identified by Western blot. ii) '*': identification of the spot inferred from MS and/or 2D-PAGE information indicate only one potential protein. iii) '/': identification of the spot inferred from MS and/or 2D-PAGE information indicate more than one potential protein. iv) the protein in the spot cannot be identified or predicted by our techniques.

https://doi.org/10.1371/journal.pone.0195639.s005

(XLSX)

S4 Table. Spots of interest Mascot search information.

This table provides Mascot search information of identified proteins. Four levels of identification are used and each one coded by a label: i) '$': identification by MS of a single protein with multiple peptides or the spot has been identified by Western blot. ii) '*': identification of the spot inferred from MS and/or 2D-PAGE information indicate only one potential protein. iii) '/': identification of the spot inferred from MS and/or 2D-PAGE information indicate more than one potential protein. iv) the protein in the spot cannot be identified or predicted by our techniques.

https://doi.org/10.1371/journal.pone.0195639.s006

(XLSX)

Acknowledgments

We would like to thank Dr. Ilai Keren for R programming and statistical advice, Dr. Paul Grieco and Dr. Kevin Spicka for synthesizing Cy3 and Cy5 for our study, and Mr. Duane Mooney for meaningful discussion and critical reading.

Citation: Keren-Aviram G, Dachet F, Bagla S, Balan K, Loeb JA, Dratz EA (2018) Proteomic analysis of human epileptic neocortex predicts vascular and glial changes in epileptic regions. PLoS ONE 13(4): e0195639. https://doi.org/10.1371/journal.pone.0195639

References

1. Yun C-H, Lee SK, Lee SY, Kim KK, Jeong SW, et al. (2006) Prognostic Factors in Neocortical Epilepsy Surgery: Multivariate Analysis. Epilepsia 47: 574–579. pmid:16529624

2. Kwan P, Brodie MJ (2000) Early Identification of Refractory Epilepsy. New England Journal of Medicine 342: 314–319. pmid:10660394

3. Loeb JA (2010) A human systems biology approach to discover new drug targets in epilepsy. Epilepsia 51: 171–177. pmid:20618426

4. Eun J-P, Choi H-Y, Kwak Y-G (2004) Proteomic analysis of human cerebral cortex in epileptic patients. Experimental & Molecular Medicine 36: 185–191.

5. Yang J-W, Czech T, Gelpi E, Lubec G (2005) Extravasation of plasma proteins can confound interpretation of proteomic studies of brain: A lesson from apo A-I in mesial temporal lobe epilepsy. Molecular Brain Research 139: 348–356. pmid:16095751

6. Loeb JA (2011) Identifying targets for preventing epilepsy using systems biology. Neuroscience Letters 497: 205–212. pmid:21382442

7. Asano E, Muzik O, Shah A, Juhász C, Chugani DC, et al. (2003) Quantitative Interictal Subdural EEG Analyses in Children with Neocortical Epilepsy. Epilepsia 44: 425–434. pmid:12614399

8. Rakhade SN, Yao B, Ahmed S, Asano E, Beaumont TL, et al. (2005) A common pattern of persistent gene activation in human neocortical epileptic foci. Annals of Neurology 58: 736–747. pmid:16240350

9. Beaumont TL, Yao B, Shah A, Kapatos G, Loeb JA (2012) Layer-Specific CREB Target Gene Induction in Human Neocortical Epilepsy. Journal of Neuroscience 32: 14389–14401a. pmid:23055509

10. Lipovich L, Dachet F, Cai J, Bagla S, Balan K, et al. (2012) Activity-Dependent Human Brain Coding/Noncoding Gene Regulatory Networks. Genetics 192: 1133–1148. pmid:22960213

11. Dachet F, Bagla S, Keren-Aviram G, Morton A, Balan K, et al. (2015) Predicting novel histopathological microlesions in human epileptic brain through transcriptional clustering. Brain.

12. Bradford M (1976) A Rapid and Sensitive Method for the Quantitation of Microgram Quantities of Protein Utilizing the Principle of Protein-Dye Binding. Analytical Biochemistry 72: 248–254. pmid:942051

13. Candiano G, Bruschi M, Musante L, Santucci L, Ghiggeri GM, et al. (2004) Blue silver: A very sensitive colloidal Coomassie G-250 staining for proteome analysis. ELECTROPHORESIS 25: 1327–1333. pmid:15174055

14. Maaty WS, Steffens JD, Heinemann J, Ortmann AC, Reeves BD, et al. (2012) Global Analysis of Viral Infection in an Archaeal Model System. Frontiers in Microbiology 3.

15. Maaty WS, Lord CI, Gripentrog JM, Riesselman M, Keren-Aviram G, et al. (2013) Identification of C-terminal Phosphorylation Sites ofN-Formyl Peptide Receptor-1 (FPR1) in Human Blood Neutrophils. Journal of Biological Chemistry 288: 27042–27058. pmid:23873933

16. Schneider CA, Rasband WS, Eliceiri KW (2012) NIH Image to ImageJ: 25 years of image analysis. Nature Methods 9: 671–675. pmid:22930834

17. Ezzet F, Pinheiro JC (2006) Linear, Generalized Linear, and Nonlinear Mixed Effects Models. Pharmacometrics: John Wiley & Sons, Inc. pp. 103–135.

18. R Development Core Team (2010) R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.

19. Maechler M RP, Struyf A, Hubert M, Hornik K. (2012) cluster: Cluster Analysis Basics and Extensions. R package version 1.14.2.

20. Gower JC (1971) A general coefficient of similarity and some of its properties. Biometrics 27: 857–874.

21. Lilley KS, Friedman DB (2004) All about DIGE: quantification technology for differential-display 2D-gel proteomics. Expert Review of Proteomics 1: 401–409. pmid:15966837

22. Fisher NI, Chambers JM, Cleveland WS, Kleiner B, Tukey PA (1984) Graphical Methods for Data Analysis. Biometrics 40: 567.

23. Roux-Dalvai F, Gonzalez de Peredo A, Simó C, Guerrier L, Bouyssié D, et al. (2008) Extensive Analysis of the Cytoplasmic Proteome of Human Erythrocytes Using the Peptide Ligand Library Technology and Advanced Mass Spectrometry. Molecular & Cellular Proteomics 7: 2254–2269.

24. Walpurgis K, Kohler M, Thomas A, Wenzel F, Geyer H, et al. (2012) Validated hemoglobin-depletion approach for red blood cell lysate proteome analysis by means of 2D PAGE and Orbitrap MS. ELECTROPHORESIS 33: 2537–2545. pmid:22899261

25. D'Amici GM, Rinalducci S, Zolla L (2011) Depletion of hemoglobin and carbonic anhydrase from erythrocyte cytosolic samples by preparative clear native electrophoresis. Nature Protocols 7: 36–44. pmid:22157974

26. Newcombe J, Woodroofe MN, Cuzner ML (1986) Distribution of Glial Fibrillary Acidic Protein in Gliosed Human White Matter. Journal of Neurochemistry 47: 1713–1719. pmid:3772373

27. Ziemba KS, Wellik KE, Hoffman-Snyder C, Noe KH, Demaerschalk BM, et al. (2011) Timing of Antiepileptic Drug Withdrawal in Adult Epilepsy Patients After Neocortical Surgical Resection. The Neurologist 17: 176–178. pmid:21532391

28. Spencer S, Huh L (2008) Outcomes of epilepsy surgery in adults and children. The Lancet Neurology 7: 525–537. pmid:18485316

29. Zhang Y, Xiong Y, Mahmood A, Zhang ZG, Chopp M (2013) Angiogenesis and Functional Recovery After Traumatic Brain Injury. Vascular Mechanisms in CNS Trauma: Springer New York. pp. 141–156.

30. Greenberg DA, Jin K (2005) From angiogenesis to neuropathology. Nature 438: 954–959. pmid:16355213

31. Papageorgiou IE, Gabriel S, Fetani AF, Kann O, Heinemann U (2011) Redistribution of astrocytic glutamine synthetase in the hippocampus of chronic epileptic rats. Glia 59: 1706–1718. pmid:21780187

32. Morin-Brureau M, Lebrun A, Rousset MC, Fagni L, Bockaert J, et al. (2011) Epileptiform Activity Induces Vascular Remodeling and Zonula Occludens 1 Downregulation in Organotypic Hippocampal Cultures: Role of VEGF Signaling Pathways. Journal of Neuroscience 31: 10677–10688. pmid:21775611

33. Rigau V, Morin M, Rousset MC, de Bock F, Lebrun A, et al. (2007) Angiogenesis is associated with blood-brain barrier permeability in temporal lobe epilepsy. Brain 130: 1942–1956. pmid:17533168

34. Wu HC, Dachet F, Ghoddoussi F, Bagla S, Fuerst D, et al. (2017) Altered metabolomic-genomic signature: A potential noninvasive biomarker of epilepsy. Epilepsia 58: 1626–1636. pmid:28714074

35. You J-S, Gelfanova V, Knierman MD, Witzmann FA, Wang M, et al. (2005) The impact of blood contamination on the proteome of cerebrospinal fluid. PROTEOMICS 5: 290–296. pmid:15672452

36. Visanji NP, Wong JC, Wang SX, Cappel B, Kleinschmidt-DeMasters BK, et al. (2012) A proteomic analysis of pediatric seizure cases associated with astrocytic inclusions. Epilepsia 53: e50–e54. pmid:22220588

37. Boulpaep EL, Boron WF (2009) FOUNDATIONS OF PHYSIOLOGY. Medical Physiology: Elsevier. pp. 3–6.

38. Mantuano E, Mukandala G, Li X, Campana WM, Gonias SL (2008) Molecular Dissection of the Human α2-Macroglobulin Subunit Reveals Domains with Antagonistic Activities in Cell Signaling. Journal of Biological Chemistry 283: 19904–19911. pmid:18499670

39. Kragh-Hansen U, Chuang VTG, Otagiri M (2002) Practical Aspects of the Ligand-Binding and Enzymatic Properties of Human Serum Albumin. Biological & Pharmaceutical Bulletin 25: 695–704.

40. Oberheim NA, Tian GF, Han X, Peng W, Takano T, et al. (2008) Loss of Astrocytic Domain Organization in the Epileptic Brain. Journal of Neuroscience 28: 3264–3276. pmid:18367594

41. Middeldorp J, Hol EM (2011) GFAP in health and disease. Progress in Neurobiology 93: 421–443. pmid:21219963

42. Eng LF (1985) Glial fibrillary acidic protein (GFAP): the major protein of glial intermediate filaments in differentiated astrocytes. Journal of Neuroimmunology 8: 203–214. pmid:2409105

43. Sarnat HB, Flores-Sarnat L (2009) α-B-Crystallin as a Tissue Marker of Epileptic Foci in Paediatric Resections. Canadian Journal of Neurological Sciences / Journal Canadien des Sciences Neurologiques 36: 566–574. pmid:19831124

44. Che Y, Piao CS, Han P-L, Lee J-K (2001) Delayed induction of α-B-crystallin in activated glia cells of hippocampus in kainic acid-treated mouse brain. Journal of Neuroscience Research 65: 425–431. pmid:11536326

45. Ousman SS, Tomooka BH, van Noort JM, Wawrousek EF, O’Conner K, et al. (2007) Protective and therapeutic role for αB-crystallin in autoimmune demyelination. Nature 448: 474–479. pmid:17568699

46. Hagemann TL, Boelens WC, Wawrousek EF, Messing A (2009) Suppression of GFAP toxicity by αB-crystallin in mouse models of Alexander disease. Human Molecular Genetics 18: 1190–1199. pmid:19129171

47. Boer K, Middeldorp J, Spliet WGM, Razavi F, van Rijen PC, et al. (2010) Immunohistochemical characterization of the out-of frame splice variants GFAP Δ164/Δexon 6 in focal lesions associated with chronic epilepsy. Epilepsy Research 90: 99–109. pmid:20430588

48. Martinian L, Boer K, Middeldorp J, Hol EM, Sisodiya SM, et al. (2009) Expression patterns of glial fibrillary acidic protein (GFAP)-delta in epilepsy-associated lesional pathologies. Neuropathology and Applied Neurobiology 35: 394–405. pmid:19508443

49. Abbott NJ, Patabendige AAK, Dolman DEM, Yusof SR, Begley DJ (2010) Structure and function of the blood–brain barrier. Neurobiology of Disease 37: 13–25. pmid:19664713

50. Seifert G, Steinhäuser C (2013) Neuron–astrocyte signaling and epilepsy. Experimental Neurology 244: 4–10. pmid:21925173

51. Heinemann U, Kaufer D, Friedman A (2012) Blood-brain barrier dysfunction, TGFβ signaling, and astrocyte dysfunction in epilepsy. Glia 60: 1251–1257. pmid:22378298

52. Aronica E, Ravizza T, Zurolo E, Vezzani A (2012) Astrocyte immune responses in epilepsy. Glia 60: 1258–1268. pmid:22331574

53. Witcher MR, Ellis TL (2012) Astroglial Networks and Implications for Therapeutic Neuromodulation of Epilepsy. Frontiers in Computational Neuroscience 6.

54. Sofroniew MV, Vinters HV (2009) Astrocytes: biology and pathology. Acta Neuropathologica 119: 7–35. pmid:20012068

55. Stewart TH, Eastman CL, Groblewski PA, Fender JS, Verley DR, et al. (2010) Chronic Dysfunction of Astrocytic Inwardly Rectifying K+ Channels Specific to the Neocortical Epileptic Focus After Fluid Percussion Injury in the Rat. Journal of Neurophysiology 104: 3345–3360. pmid:20861444

56. Kim J-E, Yeo S-I, Ryu HJ, Kim M-J, Kim D-S, et al. (2010) Astroglial loss and edema formation in the rat piriform cortex and hippocampus following pilocarpine-induced status epilepticus. The Journal of Comparative Neurology 518: 4612–4628. pmid:20886625

57. Rakhade SN, Loeb JA (2008) Focal reduction of neuronal glutamate transporters in human neocortical epilepsy. Epilepsia 49: 226–236. pmid:17868051

58. Liu X, Wen F, Yang J, Chen L, Wei Y-Q (2009) A review of current applications of mass spectrometry for neuroproteomics in epilepsy. Mass Spectrometry Reviews: n/a-n/a.

59. Abbott NJ, Rönnbäck L, Hansson E (2006) Astrocyte–endothelial interactions at the blood–brain barrier. Nature Reviews Neuroscience 7: 41–53. pmid:16371949

60. Cheng CY, Mruk DD (2009) An intracellular trafficking pathway in the seminiferous epithelium regulating spermatogenesis: a biochemical and molecular perspective. Critical Reviews in Biochemistry and Molecular Biology 44: 245–263. pmid:19622063

61. Marchi N, Angelov L, Masaryk T, Fazio V, Granata T, et al. (2007) Seizure-Promoting Effect of Blood?Brain Barrier Disruption. Epilepsia 48: 732–742. pmid:17319915

62. Carmeliet P, Lampugnani M-G, Moons L, Breviario F, Compernolle V, et al. (1999) Targeted Deficiency or Cytosolic Truncation of the VE-cadherin Gene in Mice Impairs VEGF-Mediated Endothelial Survival and Angiogenesis. Cell 98: 147–157. pmid:10428027

63. Guo R, Sakamoto H, Sugiura S, Ogawa M (2007) Endothelial cell motility is compatible with junctional integrity. Journal of Cellular Physiology 211: 327–335. pmid:17167782

64. Yap AS, Brieher WM, Gumbiner BM (1997) MOLECULAR AND FUNCTIONAL ANALYSIS OF CADHERIN-BASED ADHERENS JUNCTIONS. Annual Review of Cell and Developmental Biology 13: 119–146. pmid:9442870

65. Oby E, Janigro D (2006) The Blood?Brain Barrier and Epilepsy. Epilepsia 47: 1761–1774. pmid:17116015

66. Weissberg I, Reichert A, Heinemann U, Friedman A (2011) Blood-Brain Barrier Dysfunction in Epileptogenesis of the Temporal Lobe. Epilepsy Research and Treatment 2011: 1–10.

67. Liu JYW, Thom M, Catarino CB, Martinian L, Figarella-Branger D, et al. (2012) Neuropathology of the blood–brain barrier and pharmaco-resistance in human epilepsy. Brain 135: 3115–3133. pmid:22750659

68. Kovács R, Heinemann U, Steinhäuser C (2012) Mechanisms underlying blood-brain barrier dysfunction in brain pathology and epileptogenesis: Role of astroglia. Epilepsia 53: 53–59.

69. Marchi N, Granata T, Ghosh C, Janigro D (2012) Blood-brain barrier dysfunction and epilepsy: Pathophysiologic role and therapeutic approaches. Epilepsia 53: 1877–1886. pmid:22905812

70. Ivens S, Kaufer D, Flores LP, Bechmann I, Zumsteg D, et al. (2007) TGF- receptor-mediated albumin uptake into astrocytes is involved in neocortical epileptogenesis. Brain 130: 535–547. pmid:17121744

71. Mirza N, Vasieva O, Marson AG, Pirmohamed M (2011) Exploring the genomic basis of pharmacoresistance in epilepsy: an integrative analysis of large-scale gene expression profiling studies on brain tissue from epilepsy surgery. Human Molecular Genetics 20: 4381–4394. pmid:21852245

72. Hensley K, Venkova K, Christov A, Gunning W, Park J (2011) Collapsin Response Mediator Protein-2: An Emerging Pathologic Feature and Therapeutic Target for Neurodisease Indications. Molecular Neurobiology 43: 180–191. pmid:21271304

73. Wilson SM, Xiong W, Wang Y, Ping X, Head JD, et al. (2012) Prevention of posttraumatic axon sprouting by blocking collapsin response mediator protein 2-mediated neurite outgrowth and tubulin polymerization. Neuroscience 210: 451–466. pmid:22433297

74. Beyreuther BK, Freitag J, Heers C, Krebsfänger N, Scharfenecker U, et al. (2007) Lacosamide: A Review of Preclinical Properties. CNS Drug Reviews 13: 21–42. pmid:17461888

75. Kelemen , Halász P (2010) Lacosamide for the prevention of partial onset seizures in epileptic adults. Neuropsychiatric Disease and Treatment: 465.

76. Fornasiero EF, Bonanomi D, Benfenati F, Valtorta F (2009) The role of synapsins in neuronal development. Cellular and Molecular Life Sciences 67: 1383–1396. pmid:20035364

77. Hashimoto M, Hsu LJ, Rockenstein E, Takenouchi T, Mallory M, et al. (2002) α-Synuclein Protects against Oxidative Stress via Inactivation of the c-Jun N-terminal Kinase Stress-signaling Pathway in Neuronal Cells. Journal of Biological Chemistry 277: 11465–11472. pmid:11790792

78. Zhao F, Hu Y, Zhang Y, Zhu Q, Zhang X, et al. (2012) Abnormal expression of stathmin 1 in brain tissue of patients with intractable temporal lobe epilepsy and a rat model. Synapse 66: 781–791. pmid:22535533