- Open Access
Human neuronal changes in brain edema and increased intracranial pressure
- Nóra Faragó1, 2, 3,
- Ágnes Katalin Kocsis1,
- Csilla Braskó1,
- Sándor Lovas1,
- Márton Rózsa1,
- Judith Baka1,
- Balázs Kovács1,
- Katalin Mikite1,
- Viktor Szemenyei1,
- Gábor Molnár1,
- Attila Ozsvár1,
- Gáspár Oláh1,
- Ildikó Piszár1,
- Ágnes Zvara2,
- Attila Patócs4,
- Pál Barzó5,
- László G. Puskás2, 3 and
- Gábor Tamás1Email author
© The Author(s). 2016
Received: 13 June 2016
Accepted: 28 July 2016
Published: 4 August 2016
Functional and molecular changes associated with pathophysiological conditions are relatively easily detected based on tissue samples collected from patients. Population specific cellular responses to disease might remain undiscovered in samples taken from organs formed by a multitude of cell types. This is particularly apparent in the human cerebral cortex composed of a yet undefined number of neuron types with a potentially different involvement in disease processes. We combined cellular electrophysiology, anatomy and single cell digital PCR in human neurons identified in situ for the first time to assess mRNA expression and corresponding functional changes in response to edema and increased intracranial pressure. In single pyramidal cells, mRNA copy numbers of AQP1, AQP3, HMOX1, KCNN4, SCN3B and SOD2 increased, while CACNA1B, CRH decreased in edema. In addition, single pyramidal cells increased the copy number of AQP1, HTR5A and KCNS1 mRNAs in response to increased intracranial pressure. In contrast to pyramidal cells, AQP1, HMOX1and KCNN4 remained unchanged in single cell digital PCR performed on fast spiking cells in edema. Corroborating single cell digital PCR results, pharmacological and immunohistochemical results also suggested the presence of KCNN4 encoding the α-subunit of KCa3.1 channels in edema on pyramidal cells, but not on interneurons. We measured the frequency of spontaneous EPSPs on pyramidal cells in both pathophysiological conditions and on fast spiking interneurons in edema and found a significant decrease in each case, which was accompanied by an increase in input resistances on both cell types and by a drop in dendritic spine density on pyramidal cells consistent with a loss of excitatory synapses. Our results identify anatomical and/or physiological changes in human pyramidal and fast spiking cells in edema and increased intracranial pressure revealing cell type specific quantitative changes in gene expression. Some of the edema/increased intracranial pressure modulated and single human pyramidal cell verified gene products identified here might be considered as novel pharmacological targets in cell type specific neuroprotection.
Brain edema and an increase in intracranial pressure (ICP) might result from a number of insults including traumatic brain injury, cerebral ischemia, hypoxia, infection, brain tumors, and neuroinflammation . The localization and extent of brain damage correlates with the neurological outcome, but permanent deterioration and death are also causally linked to edema or infarct of the perilesional tissue. The mainstream treatment for edema and increased ICP has been unchanged during the last 35–40 years (cerebrospinal fluid drainage, sedation, hyperventilation, osmotherapy, steroids, barbiturate and decompressive craniectomy as an ultimum refugium) [36, 44] and no target designed pharmacological treatment specifically lowering brain edema is available to patients [36, 44].
Advances in unbiased large scale molecular approaches are proven to be successful in pinpointing individual or multiple genes associated with a disease . Omics based assessment of human pathological conditions is usually based on tissue samples containing a mixture of cell types. Blood samples are relatively easily sorted into different cell classes prior to molecular characterization, but similarly automatized methods are not readily available for biopsies taken from solid tissues [8, 28]. The cerebral cortex, composed of several distinct classes of neurons, glial cells and cell types forming the vasculature, is particularly challenging from this respect [22, 46] and pioneering laser-capture microdissection studies confirm cell type specific expression patterns in disease . Previous work showed that neuron classes recorded in acute brain slices made from biopsies of the cerebral cortex can be classified based on their electrophysiological features [21, 24, 34, 43]. Such recordings allow intracellular labeling and anatomical analysis of the recorded cells further refining the identification of cell types. Moreover, whole cell patch clamp recordings provide an opportunity for harvesting cytoplasm from electrophysiologically and morphologically identified neurons [25, 41]. Quantification of mRNA with single molecule precision from the cytoplasm collected from individual neurons [13, 23] could provide cell type specific information in addition to tissue level changes, and, in turn, electrophysiological recordings can further validate the function of therapeutically relevant genes in identified neurons . We present a workflow for detecting disease related mRNA expression changes in biopsies of the cerebral cortex which is capable of validating the alterations predicted at the tissue level with the precision of individual mRNA molecules in single neurons combined with functional validation in identified cells targeted in situ.
Materials and methods
We analysed neocortical tissue surgically removed from patients (n = 37, n = 22 female and n = 15 male, aged 49 ± 14 years) in a course of two years as part of the treatment protocol for aneurysms, brain tumors and/or increased intracranial pressure. Tissue blocks were removed from prefrontal (n = 22), temporal (n = 8) and parietal (n = 7) areas. Patients were grouped according to the presence of edema based on MRI scans showing hypointense signal in T1 weighted sequences and hyperintense signal in T2 weighted sequences  and according to measurements of ICP using Codman Microsensor™ ventricular and/or parenchymal kit . Patient group 1 (Control, n = 11, n = 5 female and n = 6 male, aged 46 ± 15 years) had no detectable edema or increase in ICP. These patients were diagnosed with unruptured aneurisms (n = 2) and benign (n = 56) or malignant (n = 3) tumors (WHO grade II) detected serendipitously at early stages. Patient group 2 (Edema, n = 18, n = 12 female and n = 6 male, aged 52 ± 11 years) showed edema on MRI scans and 0–14 Hgmm ICP. These patients were diagnosed with primary benign (n = 3) or malign brain tumors (n = 12), tumor metastases (n = 2) or abscess (n = 1). Patient group 3 (Pressure, n = 10, n = 7 female and n = 3 male, aged 42 ± 16 years) showed no edema on MRI scans and had 21–43 Hgmm ICP. These patients suffered from hydrocephalus caused by benign (n = 4) or malign (n = 4) brain tumors, by a colloid cyst (n = 1) or by Arnold-Chiari syndrome (n = 1).
All procedures were performed according to the Declaration of Helsinki with the approval of the University of Szeged Ethical Committee. Human slices were derived from material which had to be removed according to treatment protocols from prefrontal regions (gyrus frontalis superior and medialis of either hemispheres and right gyrus frontalis inferior, the left gyrus frontalis inferior was avoided in order to keep Broca’s area intact), from temporal regions of both hemispheres (gyrus temporalis superior, medialis and inferior) and from parietal regions with written informed consent patients prior to surgery over the last three years. Patients having a record of epileptic seizures and with drugs related to epilepsy were excluded from the study. Anesthesia was induced with intravenous midazolam and fentanyl (0.03 mg/kg, 1–2 μg/kg respectively). A bolus dose of propofol (1–2 mg/kg) was administered intravenously. To facilitate endotracheal intubation, the patient received 0.5 mg/kg rocuronium. After 120 s the trachea was intubated and the patient was ventilated with a mixture of O2 and N2O at a ratio of 1:2. Anesthesia was maintained with sevoflurane at minimum alveolar concentration volume of 1.2–1.5. Blocks of tissue were immersed into ice cold solution containing (in mM) 130 NaCl, 3.5 KCl, 1 NaH2PO4, 24 NaHCO3, 1 CaCl2, 3 MgSO4, 10 d (+)-glucose, saturated with 95 % O2 and 5 % CO2 in the operating theatre, sliced at a thickness of 350 μm with a vibrating blade microtome (Microm HM 650 V) and were incubated at room temperature for 1 h in the same solution. The solution used during recordings differed only in that it contained 2 mM CaCl2 and 1.5 mM MgSO4. Recordings were obtained at ~36 °C from up to four concomitantly recorded cells visualized in layer 2/3 by infrared differential interference contrast videomicroscopy at depths 60–130 μm from the surface of the slice. Micropipettes (5–7 MΩ) were filled with an intracellular solution containing 126 mM K-gluconate, 4 mM KCl, 10 mM HEPES, 10 mM creatine phosphate, 8 mM biocytin; pH 7.25; 300 mOsm supplemented with RNase Inhibitor (1U/μl, Life Technologies) to prevent any RNA degradation. Access resistance was monitored with −10 mV voltage steps in between experimental epochs and data collection was terminated if access resistance exceeded 30 MΩ. Signals were filtered at 8 kHz, digitized at 16 kHz and analyzed with PULSE software. Membrane properties of human neurons or polysynaptic events did not show significant changes for up to 20 h after slicing, but recordings included in the analysis were arbitrarily terminated 15 h after slice preparation. Traces shown are single sweeps (firing patterns) or averages or of 50–100 consecutive episodes (pharmacology). The effects during drug application were tested following 5 min wash-in periods. Visualization of biocytin and microscopy was performed as described earlier .
Single cell reverse transcription and digital PCR
At the end of electrophysiological recordings, the intracellular content was aspirated into the recording pipettes by application of a gentle negative pressure while maintaining the tight seal. Pipettes were then delicately removed to allow outside-out patch formation, and the content of the pipettes (~1.5 μl) was expelled into a low-adsorbtion test tube (Axygen) containing 0.5 μl SingleCellProtectTM (Avidin Ltd. Szeged, Hungary) solution in order to prevent nucleic acid degradation and to be compatible with direct reverse transcription reaction. Samples were snap-frozen in liquid nitrogen and stored or immediately used for reverse transcription. The reverse transcription of the harvested cytoplasm was carried out in two steps. The first step was done for 5 min at 65 °C in a total reaction volume of 5 μl containing 2 μl intracellular solution and SingleCellProtectTM mix with the cytoplasmic contents of the neuron, 0.3 μl TaqMan Assays, 0.3 μl 10 mM dNTPs, 1 μl 5× first-strand buffer, 0.3 μl 0.1 mol/L DTT, 0.3 μl RNase inhibitor (Life Technologies) and 100 U of reverse transcriptase (Superscript III, Invitrogen). The second step of the reaction was carried out at 55 °C for 1 h and then the reaction was stopped by heating at 75 °C for 15 min. The reverse transcription reaction mix was stored at −20 °C until PCR amplification.
For digital PCR analysis, half of the reverse transcription reaction mixture (2.5 μl), 2 μl TaqMan Assays (Life Technologies), 10 μl OpenArray Digital PCR Master Mix (Life Technologies) and nuclease free water (5.5 μl) were mixed in a total volume of 20 μl. The mixture was evenly distributed on an OpenArray plate. RT mixes were loaded into 4 wells of a 384-well plate from which the OpenArray autoloader transferred the cDNA master mix by capillary action into 256 nanocapillary holes (4 subarrays) on an OpenArray plate. Processing of the OpenArray slide, cycling in the OpenArray NT cycler and data analysis were done as previously described . For our dPCR protocol amplification, reactions with CT confidence values below 100 as well as reactions having CT values less than 23 or greater than 33 were considered primer dimers or background signals, respectively, and excluded from the data set.
RNA preparation, amplification and labeling
Total RNA was purified from each sample using an RNA purification kit (Macherey Nagel, Düren, Germany) according to the manufacturer’s instructions. At a final concentration of 0.8 U/μl, an RNase inhibitor (Fermentas, Lithuania) was added to the samples. RNA quantity was determined using Agilent Bioanalyzer 2100 NanoDrop 3.1.0. RNA samples were stored at −80 °C before used. An aliquot of the total RNA (1 μg) was amplified with the AminoAllyl MessageAmpTM II aRNA Amplification Kit (Ambion, USA) according to the manufacturer’s instructions. Six microgram aminoallyl-modified amplified RNA (aaRNA) was labeled with Cy3 dye according to the manufacturer’s instructions (Ambion, USA) then purified. The dye incorporation rate and the labeled aRNA concentration were detected using NanoDrop 3.1.0. The incorporation rate of the samples was 30–60 dye molecules per 1000 nucleotides.
Microarray hybridization, raw data extraction and analysis
The human oligonucleotide microarray (Whole Human Genome Microarray, 4 × 44 K, G4112F, Design ID 014850) from Agilent Technologies Inc. (Santa Clara, CA, USA) was used to determine gene expression changes. 825 ng of Cy3 labeled aaRNA, 11 μl 10× Blocking Agent and 2.2 μl 25× Fragmentation Buffer were mixed together in a final volume of 55 μl and incubated at 60 °C for 30 min then 55 μl 2× GEx Hybridization Buffer were added to each sample to stop the fragmentation reaction. All these steps were done using Gene Expression Hybridization Kit of Agilent Technologies according to the manufacturer’s instructions. 100 μl of these mixes were used to fill the 4-plex backing slides and the microarray was placed onto it. This “hybridization sandwich” was assembled in an Agilent microarray hybridization chambers. The chambers were then loaded into a hybridization rotator rack (~5 rpm) and incubated at 65 °C for 17 h. After hybridization the slides were washed in Wash buffer 1 from Agilent Technologies at room temperature for 1 min than in Wash buffer 2 at 37 °C for another 1 min before scanning. Each array was scanned at 543 nm (for Cy3 labeling) in Agilent Scanner (G2505B) using the extended dynamic range function with 5 μm resolution. Output image analysis and feature extraction was done using Feature Extraction 9.5.1 software using the single-color gene expression protocol (GE1_1100_Jul11).
All the raw data extracted by Feature Extraction software was statistically analyzed by GeneSpringGX 13.0 software from Agilent Technologies. Percentile Shift Normalization method was applied to globally normalize all the spot intensities. All the ratios were calculated from the average signal intensities in each groups (control n = 9, odema n = 12 and press n = 6) and one-way ANOVA with unequal variances (Welch) test with a twofold change cutoff was applied to determine which genes showed significantly altered expression levels. If more than a half of the data points in an experimental group had a “not detected” value - which corresponded to a not significant data point (not significantly above the background) - were excluded. Gene expression ratios with p value of < 0.05 and log2 ratio of < −1 or log2 ratio of >1 (2 fold change) were considered as repression or overexpression respectively in gene activity. The MIAME formatted microarray data are deposited to EMBL - European Bioinformatics Institute ArrayExpress system [http://www.ebi.ac.uk/microarray-as/ae/] with the Accession Number: E-MTAB-3678.
Quantitative real-time PCR (QRT-PCR)
In order to validate gene expression changes obtained by DNA microarray, QRT-PCR was performed on an Exicycler-96 instrument (Bioneer) with gene-specific TaqMan assays and TaqMan protocol to monitor gene expression. 1 μg of total RNA was reverse transcribed using the High-Capacity cDNA Archive Kit (Life Technologies) according to the manufacturer’s instructions in a final volume of 30 μL. After dilution with 60 μL of water, 1 μL of the diluted reaction mix was used as template in the QRT-PCR with FastStart TaqMan Probe Master mix (Roche Applied Science) with the following protocol: 10 min at 95 °C followed by 55 cycles of 95 °C for 10 s and 60 °C for 30 s. The fluorescence intensity of FAM dye was detected after each amplification step. Relative expression ratios were calculated as normalized ratios to rat GAPDH housekeeping gene. The final relative gene expression ratios were calculated as delta-delta Ct values. Fold change refers to 2-ΔΔCt (in the case of up-regulated genes) and – (1/2-ΔΔCt) (in the case of down-regulated genes).
Cortical sections were prepared as described above for electrophysiology at 320–500 μm thickness and immersed immediately in a fixative containing 1 % paraformaldehyde and 0.2 % (w/v) picric acid dissolved in 0.1 M phosphate buffer (pH 7.3–7.4) for 4 h. Slices were resectioned at 50 μm thickness and a proteolytic antigen retrieval method was used to localise membrane–bound epitopes. Briefly, the tissue sections were incubated at 37 °C for 10 min in 0.1 M phosphate buffer followed by 5 min in 0.2 M HCl containing preheated 0.2 mg/ml pepsin (Dako), then sections were washed in 0.1 M phosphate buffer for 2 × 10 and 1 × 30 min. For the comparison of KCNN4 expression in human control and cerebral edema cortex sections were kept in the same vial to ensure equal conditions throughout the procedure. Following several washes in tris-buffered saline, non-specific binding was blocked by incubating sections in 20 % dried milk and 0.05 % Tween 20 for 2 h at room temperature. The tissue sections were incubated with Rabbit-KCNN4 1:300 (Alomone, APC-064) or Rabbit-KCNN4 1:100 (Thermo PA5-33875) primary antibody diluted in tris-buffered saline, 2 days at 4 °C. After several washes with TBS, sections were incubated in Alexa 488 -donkey anti-rabbit (Jackson) diluted 1∶500, for 2 h at room temperature. Sections were washed 3 times in tris-buffered saline and then in 0.1 M phosphate buffer and mounted in Vectashield mounting medium (Vector Laboratories, Burlingame, CA).
Data are given as mean ± S.D., datasets were statistically compared using one-way ANOVA or Kruskal-Wallis test and Wilcoxon test was used for pharmacological experiments with SPSS software (IBM), differences were accepted as significant if p ≤ 0.05.
Comparison of microarray and QPCR measurements on homogenates of cortical gray matter derived from the Control versus Edema and Control versus Pressure groups of patients
Log2 change ([odema] vs [cont])
p value (OneWay Welch ANOVA)
Log2 change (ΔΔCt [odema] vs [cont])
p value (Student’s ttest)
Control versus Edema
Control versus Pressure
Post hoc anatomical recovery and analysis of dendritic and axonal arborizations further ensured the identity of pyramidal cells and classified interneurons as basket cells based on axonal terminals frequently observed adjacent to unlabeled somata (Fig. 1). To test potential morphological alterations in pathological conditions, we first compared the length of apical and basal dendrites and the length of axons remaining in the slices based on three-dimensional reconstructions human pyramidal cells recorded in the Control (n = 5; 3845 ± 2123, 3829 ± 1840, 17201 ± 9563 μm, respectively), Edema (n = 10; 3987 ± 1157, 4817 ± 1854, 16780 ± 7272 μm, respectively) and Pressure (n = 5; 4528 ± 1384, 4861 ± 1591, 1629 ± 3985 μm, respectively) groups, however, we detected no significant differences (p ≥ 0.12). Dendritic spines are dynamically altered in response to functional changes , but are remarkably stable in adulthood . Earlier experiments found that the density of dendritic spines differs between frontal and parietal areas in primates  and developmentally stabilizes in adulthood around the average age of patients involved in our study in humans . When measuring the density of spines on dendritic branches of basal and oblique dendrites at distances between 70 and 150 μm from the soma in samples taken from the prefrontal gyrus, we found that spine densities on dendrites of human pyramidal neurons in the Edema (n = 38, 0.38 ± 0.18/μm) and Pressure (n = 30, 0.43 ± 0.17/μm) groups were smaller of what we found on pyramidal cells of the Control (n = 31, 0.55 ± 0.23/μm, respectively, p < 0.028) group (Fig. 1). Spines receive the majority of glutamatergic input arriving to pyramidal neurons and the drop in spine densities indicate a loss of excitatory synapses in disease and, in turn, lower synapse numbers might contribute to higher input resistances measured in edema. The morphology of individual dendrites and length of the dendritic arborization of basket cells in the Control (n = 3, 4155 ± 804 μm) and Edema group (n = 3, 4497 ± 741 μm) appeared similar, axonal lenghts were not measured due to a high number of thick branches cut at the surface of slices.
The workflow presented above paves the way for detecting neuron type specific changes at molecular, electrophysiological and morphological level in various pathological conditions in the human brain. Cell type specific therapy targeted against key molecular components might form the basis in future treatment of psychiatric or neurological disorders [15, 19]. Emerging differences between the human cerebral cortex and its animal models widely used in experimentation [24, 43] might contribute to factors leading to failing clinical trials in the field  and argue for the use of human samples. In spite of recent treatment targets for cerebral edema and ICP including the Na+-K+-2Cl− co-transporter (NKCC1) and the SUR1-regulated NCCa-ATP (SUR1/TRPM4) channel , osmotherapy remains typically administered [36, 44]. Some of the edema/ICP modulated and single human pyramidal cell verified gene products identified here might be considered as novel pharmacological targets in cell type specific neuroprotection. For example, aquaporins were suggested to be crucial in extracellular space- and osmoregulation in the brain [18, 26] and are reportedly expressed in astrocytes (AQP4 and AQP9) [10, 29] and in the epithelium of the choroid plexus (AQP1) , but neuronal aquaporin expression is not clear. In agreement with earlier reports [10, 18, 29, 30], we detected no expression of aquaporins in pyramidal cells of the Control group, however, human pyramidal cells contained AQP1 and AQP3 mRNA in edema and AQP1 in patients with increased ICP, suggesting a redistribution of water (AQP1) and water plus glycerin permeable (AQP3) aquaporins in the cerebral cortex in response to disease. Furthermore, human glia-derived brain tumors (gliomas, glioblastomas, meningiomas and schwannomas) do not express AQP3  arguing for a potentially neuron specific expression. It remains to be elucidated whether potentially interneuron selective expression of AQP1 in the Control group contributes to neuron type selective water homeostasis. Differential sensitivity of different neuronal and glial populations in various types of head injury is well documented [7, 37] and cell type specific analysis under these conditions might identify underlying pathomechanisms leading to prevention or treatment.
Pharmacological and immunohistochemical results confirm our single cell digital PCR data and suggest a functionally detectable presence of KCNN4 and KCa3.1 channels in edema relative to control conditions in human pyramidal cells, however, these channels are absent in fast spiking interneurons in both groups of patients. Moreover, our results support earlier findings functionally linking KCNN4 and edema  and raise the possibility that the functional expression of Kcnn4 changes from purely glial to glial and neuronal in control and edemic conditions, respectively. A potential neuroprotective function of Kcnn4 channels on human pyramidal cells in edema might be the contribution to the suppression of EPSP summation and to the curtailment of action potentials leading to lower Ca2+ entry and excitotoxicity . Speculatively, this effect of Kcnn4 might counteract Ca2+ influx due to large amplitude EPSPs resulting from the increased input resistance and rising from a more depolarized resting membrane potential in edema. The cell type dependent absence of this mechanism suggests differential survival of neuron populations in edema making fast spiking interneurons especially vulnerable when responding to extremely powerful unitary EPSPs driving these cells to fire in human microcircuits . In addition, voltage-gated Ca2+ channels, like Cav2.2 conducting N-type Ca2+ currents encoded by the CACNA1B gene found here to be down-regulated in human pyramidal cells in edema, are essential in action potential to release coupling at presynaptic terminals [6, 40] and could effectively suppress excitatory output and counteract excitotoxicity. Indeed, blockade of N-type VGCCs decrease neuronal damage associated with ischemic brain injury in animal models  and the CACNA1B gene is associated with cerebral infarction in a human population . The voltage-gated Na+-channel β3 subunit encoding gene SNC3B found to be up-regulated in human pyramidal cells in edema might also be involved mechanisms governing excitability. Studies in Scn3b −/− mice revealed spontaneous cardiac arrhythmia and conduction abnormalities  and the recently uncovered structure of trimeric β3 subunits of Na+-channels together with their suggested interactions with α subunits  predict future testing of β subunit dependent μ-conotoxin sensitivity of Na+-channels  in disease or treatment. Pyramidal neuron protective regulation of KCNN4, CACNA1B, and SCN3B suggested here could be effectively boosted by other cytoprotective processes involving heme oxygenase 1  and mitochondrial superoxide dismutase 2 encoded by the HMOX1 and SOD2 genes which are up-regulated in pyramidal neurons in edema. The voltage-gated K+-channel α subunit Kv9.1 encoded by the KCNS1 gene up-regulated in human pyramidal cells in the Pressure group has been found to be associated with multiple chronic pain states , however, interpretation of our findings is difficult knowing that electrically silent Kv9.1 subunits are not capable of forming functional homo-multimeric channels  but can suppress currents mediated by Kv2 and Kv3 α subunit families .
We applied the combination of cellular electrophysiology, anatomy and single cell digital PCR to human neurons identified in brain slices in situ to assess mRNA expression and corresponding functional changes in response to edema and increased intracranial pressure. In material derived from patients with edema, mRNA copy numbers of AQP1, AQP3, HMOX1, KCNN4, SCN3B and SOD2 increased, while CACNA1B, CRH decreased in individual pyramidal cells, however, AQP1, HMOX1and KCNN4 remained unchanged in single cell digital PCR performed on fast spiking cells. In single pyramidal cells derived from patients with increased intracranial pressure, the copy number of AQP1, HTR5A and KCNS1 mRNAs was increased. Pharmacological and immunohistochemical results also suggested the presence of KCNN4 encoding the α-subunit of KCa3.1 channels in edema on pyramidal cells, but not on interneurons validating single cell digital PCR data. The frequency of spontaneous EPSPs decreased on pyramidal cells in both pathophysiological conditions and on fast spiking interneurons in edema. This was parallelled by an increase in input resistances on both cell types and by a drop in dendritic spine density on pyramidal cells consistent with a loss of excitatory synapses. The anatomical and/or physiological changes in individually identified human pyramidal and fast spiking cells in edema and increased intracranial pressure reveal cell type specific quantitative changes in gene expression and some of the modulated gene products identified here might be considered as novel pharmacological targets in cell type specific neuroprotection.
The authors thank Éva Tóth and Nelli Tóth for technical assistance. This work was supported by the ERC INTERIMPACT project (GT), the Hungarian Academy of Sciences (GT), the National Research, Development and Innovation Office of Hungary (VKSZ-14-1-2015-0155, NF, CB, GT), and the National Brain Research Program (GT, MG).
RNA preparation, amplification and labelling, single cell reverse transcription, digital PCR and quantitative real-time PCR (QRT-PCR) was performed by NF and KÁK, immunocytochemistry was performed by CB, microarray hybridization, raw data extraction and analysis was performed by ÁZ and AP, electrophysiology was performed by MR, SL, GM, BK, KM, IP, AO, GO, anatomical reconstructions were performed by JB, neurosurgery was performed by PB and LGP and GT designed experiments and wrote the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Altshuler D, Daly MJ, Lander ES. Genetic mapping in human disease. Science. 2008;322:881–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Bhatt DH, Zhang S, Gan W-B. Dendritic spine dynamics. Annu Rev Physiol. 2009;71:261–82.View ArticlePubMedGoogle Scholar
- Bocksteins E, Snyders DJ. Electrically silent Kv subunits: their molecular and functional characteristics. Physiology (Bethesda). 2012;27:73–84.View ArticleGoogle Scholar
- Bourne JN, Harris KM. Balancing structure and function at hippocampal dendritic spines. Annu Rev Neurosci. 2008;31:47–67.View ArticlePubMedPubMed CentralGoogle Scholar
- Bratton SL, Chestnut RM, Ghajar J, McConnell Hammond FF, Harris OA, Hartl R, Manley GT, Nemecek A, Newell DW, Rosenthal G, Schouten J, Shutter L, Timmons SD, Ullman JS, Videtta W, Wilberger JE, Wright DW. Guidelines for the management of severe traumatic brain injury. VII. Intracranial pressure monitoring technology. J Neurotrauma. 2007;24 Suppl 1:S45–54.PubMedGoogle Scholar
- Catterall WA. Voltage-gated calcium channels. Cold Spring Harb Perspect Biol. 2011;3:a003947.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen Y, Swanson RA. Astrocytes and brain injury. J Cereb Blood Flow Metab. 2003;23:137–49.View ArticlePubMedGoogle Scholar
- Cookson W, Liang L, Abecasis G, Moffatt M, Lathrop M. Mapping complex disease traits with global gene expression. Nat Rev Genet. 2009;10:184–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Costigan M, Belfer I, Griffin RS, Dai F, Barrett LB, Coppola G, Wu T, Kiselycznyk C, Poddar M, Lu Y, Diatchenko L, Smith S, Cobos EJ, Zaykin D, Allchorne A, Gershon E, Livneh J, Shen P-H, Nikolajsen L, Karppinen J, Männikkö M, Kelempisioti A, Goldman D, Maixner W, Geschwind DH, Max MB, Seltzer Z, Woolf CJ. Multiple chronic pain states are associated with a common amino acid-changing allele in KCNS1. Brain. 2010;133:2519–27.View ArticlePubMedPubMed CentralGoogle Scholar
- Elkjaer M, Vajda Z, Nejsum LN, Kwon T, Jensen UB, Amiry-Moghaddam M, Frøkiaer J, Nielsen S. Immunolocalization of AQP9 in liver, epididymis, testis, spleen, and brain. Biochem Biophys Res Commun. 2000;276:1118–28.View ArticlePubMedGoogle Scholar
- Elston GN. Pyramidal cells of the frontal lobe: all the more spinous to think with. J Neurosci. 2000;20:RC95.PubMedGoogle Scholar
- Engbers JDT, Anderson D, Asmara H, Rehak R, Mehaffey WH, Hameed S, McKay BE, Kruskic M, Zamponi GW, Turner RW. Intermediate conductance calcium-activated potassium channels modulate summation of parallel fiber input in cerebellar Purkinje cells. Proc Natl Acad Sci U S A. 2012;109:2601–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Faragó N, Kocsis ÁK, Lovas S, Molnár G, Boldog E, Rózsa M, Szemenyei V, Vámos E, Nagy LI, Tamás G, Puskás LG. Digital PCR to determine the number of transcripts from single neurons after patch-clamp recording. Biotechniques. 2013;54:327–36.View ArticlePubMedGoogle Scholar
- Hakim P, Brice N, Thresher R, Lawrence J, Zhang Y, Jackson AP, Grace AA, Huang CL-H. Scn3b knockout mice exhibit abnormal sino-atrial and cardiac conduction properties. Acta Physiol (Oxf). 2010;198:47–59.View ArticleGoogle Scholar
- Hyman SE. Revitalizing psychiatric therapeutics. Neuropsychopharmacology. 2014;39:220–9.View ArticlePubMedGoogle Scholar
- Jazwa A, Cuadrado A. Targeting heme oxygenase-1 for neuroprotection and neuroinflammation in neurodegenerative diseases. Curr Drug Targets. 2010;11:1517–31.View ArticlePubMedGoogle Scholar
- Khanna R, Roy L, Zhu X, Schlichter LC. K+ channels and the microglial respiratory burst. Am J Physiol Cell Physiol. 2001;280:C796–806.PubMedGoogle Scholar
- King LS, Kozono D, Agre P. From structure to disease: the evolving tale of aquaporin biology. Nat Rev Mol Cell Biol. 2004;5:687–98.View ArticlePubMedGoogle Scholar
- Krystal JH, State MW. Psychiatric disorders: diagnosis to therapy. Cell. 2014;157:201–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Mauler F, Hinz V, Horváth E, Schuhmacher J, Hofmann HA, Wirtz S, Hahn MG, Urbahns K. Selective intermediate-/small-conductance calcium-activated potassium channel (KCNN4) blockers are potent and effective therapeutics in experimental brain oedema and traumatic brain injury caused by acute subdural haematoma. Eur J Neurosci. 2004;20:1761–8.View ArticlePubMedGoogle Scholar
- McCormick DA, Williamson A. Convergence and divergence of neurotransmitter action in human cerebral cortex. Proc Natl Acad Sci U S A. 1989;86:8098–102.View ArticlePubMedPubMed CentralGoogle Scholar
- Mirnics K, Levitt P, Lewis DA. Critical appraisal of DNA microarrays in psychiatric genomics. Biol Psychiatry. 2006;60:163–76.View ArticlePubMedGoogle Scholar
- Molnár G, Faragó N, Kocsis ÁK, Rózsa M, Lovas S, Boldog E, Báldi R, Csajbók É, Gardi J, Puskás LG, Tamás G. GABAergic neurogliaform cells represent local sources of insulin in the cerebral cortex. J Neurosci. 2014;34:1133–7.View ArticlePubMedGoogle Scholar
- Molnár G, Oláh S, Komlósi G, Füle M, Szabadics J, Varga C, Barzó P, Tamás G. Complex events initiated by individual spikes in the human cerebral cortex. PLoS Biol. 2008;6:e222.View ArticlePubMedPubMed CentralGoogle Scholar
- Monyer H, Lambolez B. Molecular biology and physiology at the single-cell level. Curr Opin Neurobiol. 1995;5:382–7.View ArticlePubMedGoogle Scholar
- Nagelhus EA, Ottersen OP. Physiological roles of aquaporin-4 in brain. Physiol Rev. 2013;93:1543–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Namadurai S, Balasuriya D, Rajappa R, Wiemhöfer M, Stott K, Klingauf J, Edwardson JM, Chirgadze DY, Jackson AP. Crystal structure and molecular imaging of the Nav channel β3 subunit indicates a trimeric assembly. J Biol Chem. 2014;289:10797–811.View ArticlePubMedPubMed CentralGoogle Scholar
- Navin N, Kendall J, Troge J, Andrews P. Tumour evolution inferred by single-cell sequencing. Nature. 2011;472:90–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Nielsen S, Nagelhus EA, Amiry-Moghaddam M, Bourque C, Agre P, Ottersen OP. Specialized membrane domains for water transport in glial cells: high-resolution immunogold cytochemistry of aquaporin-4 in rat brain. J Neurosci. 1997;17:171–80.PubMedGoogle Scholar
- Nielsen S, Smith BL, Christensen EI, Agre P. Distribution of the aquaporin CHIP in secretory and resorptive epithelia and capillary endothelia. Proc Natl Acad Sci U S A. 1993;90:7275–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Niu D, Kondo T, Nakazawa T, Yamane T, Mochizuki K, Kawasaki T, Matsuzaki T, Takata K, Katoh R. Expression of Aquaporin3 in human neoplastic tissues. Histopathology. 2012;61:543–51.PubMedGoogle Scholar
- Petanjek Z, Judaš M, Šimic G, Rasin MR, Uylings HBM, Rakic P, Kostovic I. Extraordinary neoteny of synaptic spines in the human prefrontal cortex. Proc Natl Acad Sci U S A. 2011;108:13281–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Pietersen CY, Lim MP, Macey L, Woo T-UW, Sonntag KC. Neuronal type-specific gene expression profiling and laser-capture microdissection. Methods Mol Biol. 2011;755:327–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Prince DA, Wong RK. Human epileptic neurons studied in vitro. Brain Res. 1981;210:323–33.View ArticlePubMedGoogle Scholar
- Rabinstein A. Treatment of cerebral edema. Neurologist. 2006;12:59–73.View ArticlePubMedGoogle Scholar
- Ropper AH. Hyperosmolar therapy for raised intracranial pressure. N Engl J Med. 2012;367:746–52.View ArticlePubMedGoogle Scholar
- Scharfman HE, Myers CE. Hilar mossy cells of the dentate gyrus: a historical perspective. Front Neural Circuits. 2012;6:106.PubMedGoogle Scholar
- Schwab A, Oberleithner H. Plasticity of renal epithelial cells: the way a potassium channel supports migration. Pflugers Arch Eur J Physiol. 1996;432:R87–93.Google Scholar
- Stocker M, Hellwig M, Kerschensteiner D. Subunit assembly and domain analysis of electrically silent K+ channel alpha-subunits of the rat Kv9 subfamily. J Neurochem. 1999;72:1725–34.View ArticlePubMedGoogle Scholar
- Sudhof TC. The synaptic vesicle cycle. Annu Rev Neurosci. 2004;27:509–47.View ArticlePubMedGoogle Scholar
- Sugino K, Hempel CM, Miller MN, Hattox AM, Shapiro P, Wu C, Huang ZJ, Nelson SB. Molecular taxonomy of major neuronal classes in the adult mouse forebrain. Nat Neurosci. 2006;9:99–107.View ArticlePubMedGoogle Scholar
- Valentino K, Newcomb R, Gadbois T, Singh T, Bowersox S, Bitner S, Justice A, Yamashiro D, Hoffman BB, Ciaranello R. A selective N-type calcium channel antagonist protects against neuronal loss after global cerebral ischemia. Proc Natl Acad Sci U S A. 1993;90:7894–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Verhoog MB, Goriounova NA, Obermayer J, Stroeder J, Hjorth JJJ, Testa-Silva G, Baayen JC, de Kock CPJ, Meredith RM, Mansvelder HD. Mechanisms underlying the rules for associative plasticity at adult human neocortical synapses. J Neurosci. 2013;33:17197–208.View ArticlePubMedGoogle Scholar
- Walcott BP, Kahle KT, Simard JM. Novel treatment targets for cerebral edema. Neurotherapeutics. 2012;9:65–72.View ArticlePubMedGoogle Scholar
- Yamaguchi M, Nakayama T, Fu Z, Sato N, Soma M, Morita A, Hinohara S, Doba N, Mizutani T. The haplotype of the CACNA1B gene associated with cerebral infarction in a Japanese population. Hereditas. 2010;147:313–9.View ArticlePubMedGoogle Scholar
- Yao W, Wu C. Exploring the genetic underpinnings of brain and behavioral disorders. J Neurogenet. 2014;28:1–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang M-M, Wilson MJ, Azam L, Gajewiak J, Rivier JE, Bulaj G, Olivera BM, Yoshikami D. Co-expression of Na(V)β subunits alters the kinetics of inhibition of voltage-gated sodium channels by pore-blocking μ-conotoxins. Br J Pharmacol. 2013;168:1597–610.View ArticlePubMedPubMed CentralGoogle Scholar