- Open Access
Spatiotemporal expression patterns of genes coding for plasmalemmal chloride transporters and channels in neurological diseases
Molecular Brain volume 16, Article number: 30 (2023)
Neuronal voltage changes which are dependent on chloride transporters and channels are involved in forming neural functions during early development and maintaining their stability until adulthood. The intracellular chloride concentration maintains a steady state, which is delicately regulated by various genes coding for chloride transporters and channels (GClTC) on the plasmalemma; however, the synergistic effect of these genes in central nervous system disorders remains unclear. In this study, we first defined 10 gene clusters with similar temporal expression patterns, and identified 41 GClTC related to brain developmental process. Then, we found 4 clusters containing 22 GClTC were enriched for the neuronal functions. The GClTC from different clusters presented distinct cell type preferences and anatomical heterogeneity. We also observed strong correlations between clustered genes and diseases, most of which were nervous system disorders. Finally, we found that one of the most well-known GClTC, SLC12A2, had a more profound effect on glial cell-related diseases than on neuron-related diseases, which was in accordance with our observation that SLC12A2 was mainly expressed in oligodendrocytes during brain development. Our findings provide a more comprehensive understanding of the temporal and spatial expression characteristics of GClTC, which can help us understand the complex roles of GClTC in the development of the healthy human brain and the etiology of brain disorders.
Chloride, the most abundant physiological anion, is involved in various biophysical processes, including osmotic equilibrium, membrane potential modulation, and intracellular pH maintenance . The electrochemical gradient of chloride is determined by the extracellular and intracellular concentrations as dictated by the Nernst equation, and the corresponding equilibrium potential (ECl) and the membrane potential (Vm) set the driving force (Vm-ECl) of Cl− currents mediated by γ-aminobutyric acid sub-type A and glycine receptors (GABAARs and GlyRs) .
During brain development, the intracellular chloride concentration of neural cells changes dynamically. This dynamic variance plays a pivotal role in the regulation of neural function [3, 4]. For example, the action of γ-aminobutyric acid (GABA) critically depends on the changes of intracellular chloride concentration. In the early stages of development, the intracellular chloride concentration is relatively high. ECl is above Vm, and thus GABAAR activation results in chloride ions efflux and neuronal depolarization. Depolarizing response not only promotes developmental processes, but also governs developmental switching points, from proliferation to migration, from migration to differentiation, and finally from differentiation to synapse formation . We in the later stage of development, the level of intracellular chloride is relatively low. ECl is lower than Vm, which allows adult neurons to hyperpolarize in response to GABAAR activation by GABA. Hyperpolarizing response regulates the electrical activity and tightens neuronal plasticity to optimize neuronal networks [6, 7].
Although multiple factors can influence chloride regulation , the chloride channels and transporters are the direct executors to mediate chloride ions to enter or exit the cells to maintain intracellular chloride homeostasis, such as ligand-gated channels, the cystic fibrosis transmembrane conductance regulator (CFTR), and cation-chloride co-transporters [9, 10]. In this study, we use the abbreviation “GClTC” to refer to the genes coding for plasmalemmal chloride transporters and channels. Several studies have highlighted that the dysfunction of these GClTC leads to central nervous system (CNS) diseases, including autism spectrum disorder (ASD), neonatal epilepsy, and schizophrenia [11,12,13,14]. However, the underlying mechanisms are not fully understood.
Chloride homeostasis is affected by the transcriptional regulations of multiple GClTC. However, current studies have focused more on the roles of single genes in a restricted time and space, instead of the roles of a group of GClTC on a dynamic spatio-temporal scale . In this study, we first clustered the genes with similar temporal expression patterns and focused on 4 out of 10 clusters whose functions were related to the CNS. Then we investigated the biological roles of GClTC in those clusters via Gene Ontology (GO) database, which provides a specific biological information of protein functions. Finally, based on the cellular and spatial distribution characteristics of GClTC, we expounded on the relationship between the gene expression features and CNS diseases, which provided a more comprehensive view of the underlying mechanisms of CNS disorders induced by disrupted chloride homeostasis.
Gene expression profile acquisition
To explore the specific gene expression in the human brain at different developmental periods, RNA-seq data from different developmental periods were downloaded from the Atlas of the Developing Human Brain  (http://www.brainspan.org/static/download). Brain single-cell transcriptome data were obtained from the Multiple Cortical Areas-Smart-Seq  (https://portal.brain-map.org/atlases-and-data/rnaseq). The gene expression landscape, cell type heterogeneity, and temporal dynamics across the human brain were obtained using the Spatio-Temporal cell Atlas of the human Brain (STAB) . Heatmaps of the gene expression profiles and enrichment analyses were generated using TBtools . The workflow chart is shown in Fig. 1.
Gene expression pattern analysis
Gene expression patterns during development were analyzed using Short Time-Series Expression Miner (STEM) software . The options were set to default because STEM already provided optimal results for both biological and simulated data. The expression parameter used in this analysis was the reads per kilobase million (RPKM) value, and a p-value < 0.05 for the clustered profile was considered significant. The related R packages were used for data visualization.
Functions, pathways, and diseases associated with developmental gene expression patterns
Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of different clusters were performed using the R package “clusterprofiler”, and p-values for the representative GO and KEGG terms shown in this study were adjusted for multiple testing with the appropriate Benjamini–Hochberg correction. The functional categories and pathways with adjusted p-values < 0.05 were deemed to be significantly enriched. Gene-Disease Association (GDA) analysis of differentially expressed genes was performed according to the DisGeNET association-type ontology . In this analysis, a q-value < 0.05 was considered significant.
Spatial gene expression analysis
The spatial expression information of genes in the brain was obtained using the Allen Institute Human Brain Atlas . Whole-brain gene expression maps for GClTC were downloaded from the Neurosynth  (https://neurosynth.org/genes/). Functional connectivity was offered by Neurosynth, and we adopted this index to determine the relative strengths of the associations. This analysis yielded the Pearson correlation coefficient (r) between the whole-brain (uncorrected) reverse-inference map for each term and the functional connectivity map seeded at the current location. Gene maps and whole-brain maps were merged using Python.
Gene expression analysis in specific diseases
Specific disease gene expression data were acquired from the Gene Expression Omnibus (GEO) database GSE108000  and GSE122228 . Gene set enrichment analysis (GSEA) was conducted to identify differences in the enrichment of GClTC between the control and disease groups. Shapiro–Wilk test was applied to assess data normality. Bartlett's test was applied to estimate the variance homogeneity. The Student’s t-test or Wilcoxon test was used to analyze the statistical difference. These statistical analyses were performed by R package "stats". Differences were considered statistically significant if p < 0.05.
A developmental gene expression landscape of brain
Similar neural functions are often generated in a specific time window, and associative transcriptomic states present consistent neuronal development dynamics [26, 27]. Studying transcriptional regulation across brain development periods can facilitate the understanding of brain functions and cellular mechanisms that underlie brain diseases . Short time-series expression miner (STEM) software is widely used to study gene expression data of biological processes with short time series . It adopts a clustering algorithm to classify time series gene expression trends. In this study, we used STEM software to cluster genes with similar temporal expression patterns, whose expression changes synchronously during brain development.
The transcriptomic data for human brain genes were obtained from the Atlas of the Developing Human Brain database. We focused on transcriptomic changes during hippocampal development as it is a representative biological process in the development of brain regions . Ten time points (from 12 post-conception weeks to 18 years) representing different developmental periods were selected for further analysis (Additional file 1: Table S1). A total of 9543 genes were screened out by STEM, and we defined these genes with temporal expression features as “developmental-stage dependent genes”. They were further classified into different groups as clusters, according to their unique time sequence change characteristics and different expression kinetics. 10 clusters with statistical differences were selected. They were ordered by their statistical significance in Fig. 2a and Additional file 2: Table S2. The number of clusters is marked in a pie chart (Fig. 2b).
Considering the intracellular chloride ion concentration is directly mediated by chloride transporters and channels, we focus on their GClTC in our study and explore their roles during development. We first identified the target genes using the term “chloride transmembrane transport (GO:1902476): the process in which chloride is transported across a membrane” in the Gene Ontology database . 454 genes of “Homo sapiens” were included in this term. After ruling out the duplicated gene names, 104 GClTC were finally selected for further analysis (Additional file 3: Table S3). Next, we intersected these 104 GClTC with 9543 developmental-stage dependent genes. A total of 41 GClTC were identified (Fig. 2c and Table 1). They were evenly distributed among different clusters, and no obvious clustering was observed in one group.
Pathways associated with expression patterns
Genes with temporal expression features may be involved in time-dependent developmental functions. Therefore, we conducted GO and KEGG enrichment analyses of the genes from 10 clusters. Each cluster exhibited unique biological characteristics (Additional file 8: Fig. S1 and Additional file 4: Table S4).
Notably, we found genes from clusters 4, 6, 8, and 9 were mainly enriched for CNS-specific terms (Fig. 3). For example, in biological processes (BP) most terms of clusters 4 and 9 were related to signal regulation in the CNS. Cluster 6 contained terms for the regulation of cell morphogenesis. Cluster 8 included terms related to ensheathment formation, which belonged to non-neuronal cell components. Moreover, the enriched cellular component (CC) analysis corresponded with the BP-associated locations, such as synapse in clusters 4 and 9, and myelin sheath in cluster 8. KEGG analysis revealed that clusters 4, 6, and 9 tended to trigger pathways involved in regulating neural signals, such as phosphatidylinositol signaling pathway, oxytocin signaling pathway, and axon guidance. According to the above results, the genes in clusters 4, 6, 8, and 9 may be of great importance in neural functions. We observed that 22 GClTC were present in those 4 clusters, suggesting they may regulate the corresponding neuronal functions.
Heterogeneity of GClTC expression across cell subtypes
We then investigated the distribution of GClTC involved in the regulation of CNS functions in different cell types. The single-cell transcriptomics data were obtained from the Multiple Cortical Areas-Smart-Seq . As shown in Fig. 4, the genes in clusters 4 and 6 showed no preference for certain cell types; however, the genes in cluster 9 were expressed mainly in neurons, with minimal expression in non-neuronal cells. The genes in cluster 8 were mainly expressed in glial cells. To investigate whether these genes and the specific cell types were determined during early development, we selected some GClTC using the Spatio-Temporal cell Atlas of the human Brain (STAB) to evaluate their expression dynamics in different cell types. Representative genes selected from the different clusters are shown in Additional file 9: Fig. S2. The results indicated that remarkable changes in expression levels during development were detected in specific cell types, consistent with the distribution of these genes in adults.
Spatial specificity of GClTC
To explore the spatial expression features of GClTC clusters in the brain, we downloaded spatial gene expression maps of specific genes of interest from the Neurosynth . Neurosynth is a large database of mappings structured by text-mining, meta-analysis, and machine-learning techniques. We used the Neurosynth-Gene module to explore gene expression levels in a whole-brain map. The distribution of all genes of interest is shown in Additional file 10: Fig. S3 and Additional file 5: Table S5, and the representative genes from different clusters are shown in Fig. 5. We found that the members of cluster 4 were observed with relatively higher expression levels in the cortical area. The GClTC in cluster 6 were mainly expressed in the basal brainstem and ganglia. Moreover, the GClTC in cluster 8 were also highly expressed in the brainstem, basal ganglia, and thalamus. The gene expression levels in cluster 9 were higher in the visual cortex, temporal lobe, parietal lobe, and hippocampus. These results showed that GClTC present anatomical heterogeneity in different clusters.
GClTC enriched for different CNS disorders
The GClTC from different clusters have unique expression characteristics. As the genes in clusters 4, 6, 8, and 9 were enriched for some CNS-specific terms, we further investigated the connections between those genes and diseases. The DisGeNET database, which offers rich information of genes and variants associated with many kinds of human diseases, was used to enrich for the diseases with all genes in clusters 4, 6, 8, and 9. We observed that genes from these 4 clusters were enriched for distinct diseases, most of which were CNS disorders (Fig. 6 and Additional file 6: Table S6). Clusters 4 and 9 were positively and similarly enriched for epilepsy. Cluster 6 was enriched for neurodevelopmental disorders. “Demyelinating diseases” was glaringly enriched in cluster 8.
As shown in our findings (Figs. 3, 4, 5), GClTC in different clusters exhibited functional and anatomical heterogeneity. We speculated that the heterogeneity might account for their different roles in different CNS diseases. To test our hypothesis, we selected typical GClTC to illustrate the intrinsic connection between their expression features and CNS diseases.
The neuronal chloride ion gradient is primarily established by cation-chloride cotransporters (CCCs) from the SLC12A gene family, especially the chloride ion importer sodium potassium chloride cotransporter-1 (NKCC1) and the chloride ion exporter potassium chloride cotransporter-2 (KCC2). NKCC1 and KCC2 encoded by SLC12A2 and SLC12A5, respectively . The abnormal transcriptional regulations of these two genes lead to multiple neurological disorders resulting from disrupted chloride homeostasis. According to the previous analysis (Table 1 and Fig. 6), SLC12A2 was classified in cluster 8 related to demyelinating diseases, while SLC12A5 was classified in cluster 9 related to epilepsy. To further confirm their relationship with the enriched diseases, we examined their heterogeneity in two Gene Expression Omnibus (GEO) datasets, multiple sclerosis (MS) (GSE108000) (one of the most common demyelinating diseases) and epilepsy (GSE122228), both of which were strongly associated with chloride homeostasis disorder [32, 33]. We next defined 104 GClTC as a gene set (Additional file 2: Table S2). GSEA was conducted to determine whether those genes showed statistically significant, concordant differences in MS and epilepsy. The results showed that the defined gene set was significantly enriched in multiple sclerosis (p < 0.05) (Fig. 7a) and epilepsy (p < 0.05) (Fig. 7b), indicating that GClTC were involved in disease mechanisms. Then, we compared mRNA levels of SLC12A2 and SLC12A5 in MS and epilepsy separately (Fig. 7c, d and Additional file 7: Table S7). SLC12A2 showed no significant difference in models with epilepsy compared to the controls (control mean ± SD = 0.952 ± 0.943, KA mean ± SD = 0.412 ± 0.994, p = 0.532, Student's t-test), but it was significantly decreased in multiple sclerosis [control median (IQR) = 0.511 (0.392), MS median (IQR) = 0.200 (0.597), p = 0.035, Wilcoxon test]. In contrast, SLC12A5 expression levels were significantly decreased in epilepsy (control mean ± SD = − 0.086 ± 0.486, KA mean ± SD = − 1.428 ± 0.192, p = 0.028, Student's t-test) but not in demyelinating diseases [control median (IQR) = − 1.530 (1.246), MS median (IQR) = − 1.955 (0.451), p = 0.237, Wilcoxon test]. This result implied that the expressions of SLC12A5 and SLC12A2 had different effects on the pathological mechanisms of these two diseases.
Abnormal chloride homeostasis occurs in various CNS disorders; however, the underlying mechanisms are not yet fully understood. Chloride homeostasis is delicately regulated by various GClTC, so we explored the relationship between GClTC and neurological diseases according to their temporal and spatial expression patterns and provided a new perspective to illustrate the intrinsic connections between them.
As temporal demands drive the evolution of neural diversity, time is regarded as a key metric for all brain operations. With neural signals sending out specific programs to start or stop neurogenesis, distinct functions are generated in specific time windows during neuronal development, which is driven by multiple molecular events [34, 35]. Therefore, those gene expression that change synchronously during development may play a key role in the establishment of brain functions from simple to complex. Thus, understanding the transcriptional dynamics of crucial genes may help interpret the functional heterogeneity underlying neurodevelopmental processes. Time series analysis can help to understand related information surrounding biological processes, and it has been used in studies of a variety of species, including Drosophila, Arabidopsis thaliana, and so on [36, 37]. The STEM, a time series analysis software, allows us to screen genes with similar variation expression dynamics in mRNA levels. In this study, we used STEM to cluster genes whose expression changed synchronously during brain development. As chloride-mediated signals are involved in the formation of neural function during early development, we believe that taking a developmental perspective in investigation can help us better understand the neural regulatory mechanisms where GClTC are involved.
To illustrate the relevance between expression trends and biological functions, we clustered genes with similar expression dynamics and enriched their biological functions. Some clusters targeted enrichment terms for transcriptional regulation, division regulation, and energy metabolism, which belong to general cell functions (Additional file 8: Fig. S1). Interestingly, some enrichment terms were related to specific neural cell functions. For instance, genes in cluster 4 were enriched for the function of synaptic components, signal transmissions, and targeted neuron-related pathways such as the phosphatidylinositol signaling pathway. The phosphatidylinositol signaling system regulates cellular functions such as receptor signaling, secretion, and endocytosis. The disorder of phosphoinositide system leads to brain fail to react appropriate responsive to stimuli, which is associated with several neurological diseases [38, 39]. Therefore, we suggest that classifying genes based on their temporal expression may help screen for biological characteristics related to neurodevelopmental processes.
We further focused on the GClTC in these CNS-related clusters. Most GClTC have been reported to be involved in the corresponding functions, as shown in our results. These gene functions are consistent with the enriched term in cluster 4 (regulating synaptic signaling), Anoctamin 3 (ANO3) has been reported to control the excitability of synapses in the hippocampus in response to hyperthermia [40, 41]. Solute carrier family 12 member 6 (SLC12A6) and solute carrier family 12 member 6 (SLC12A7) help maintain cell volume and intracellular chloride levels , concordant with their functions in cluster 6. The gene tweety family member 2 (TTYH2) in cluster 8 constitutes the volume-regulated anion channel (VRAC) that decreases intracellular volumes during cell swelling [43, 44]; it executes glial-related functions. This suggested that GClTC from different clusters might play different and specific roles as opposed to all the GClTC playing the same role during development.
The similarities between gene in expression patterns may reflect cell type-specific expression . Thus, we suppose that this distinction of functions might result from differences in the cell types. As mentioned before, clusters 4, 6, and 9 were enriched for neural functions, consistent with our finding that the GClTC in those clusters were mainly detected in neurons. For example, the family of GABA receptor genes, which are mainly expressed in neurons , are found in cluster 9. Furthermore, the expression of GABRA5 and GABRD were mostly detected in excitatory neurons, indicating their important role in this cell type. It has been shown that GABRA5 contributes to excitatory synapse development [47, 48]. Moreover, GABRD codes for the delta subunit of the GABA receptor, which generates tonic inhibition [49, 50]; its effect on excitatory neurons warrants further study. In addition, the data from STAB (Additional file 9: Fig. S2) indicated that dramatic changes in gene expression levels are present in specific cells during development, which confirm the important role of GClTC in neurological function formation. However, it should be noted that the data from STAB only reported gene expression levels at a few time points after birth; therefore, detailed studies for gene expression levels in the postnatal stage are necessary.
The brain is functionally organized into distinct regions. We speculated that GClTC with similar functions might have similar spatial distributions in the brain. Our results showed that the GClTC in the same clusters presented higher expression in specific brain areas, and the differences in distribution were shown between clusters. Considering that abnormal connections between certain brain regions appear in different diseases, we further analyzed the correspondence between gene expression distribution and diseases. We noticed that the brain regions with high gene expression levels were in accordance with the pathological sites of these diseases. For example, the genes in clusters 4 and 9 were mainly expressed in the cortex, hippocampus, and temporal lobe, which constitute seizure circuits and regulate different stages of seizures [51, 52]. The pathogenesis of epilepsy involves the disruption of synaptic transmission and membrane potential , which also corresponds with the functional enrichment analysis, indicating that the GClTC in clusters 4 and 9 may play more critical roles in the mechanism of seizures than other GClTC. We also noticed that some of the genes in these two clusters encode GABA receptors and that these genes were associated with a broad spectrum of epileptic syndromes. Variants in GABA receptor genes can either be associated with monogenic epileptic syndromes or constitute risk factors for genetically complex generalized epilepsies . Interestingly, some genes in cluster 9 were expressed in the hippocampus, whereas genes in cluster 4 were less expressed in this region. This may provide a hint that genes in cluster 9 might have a higher therapeutic significance for epilepsy, whose pathology involves the hippocampus.
In previous studies, chloride has been identified as a signal that broadly regulates neural activity in the CNS. Thus far, we have suggested that GClTC might play specific roles in certain diseases because they have specific molecular characteristics, cell type distributions, and brain region distributions. Moreover, we suggest that disease-related chloride dyshomeostasis may not be caused by a single GClTC but rather by a group of genes with the same expression characteristics.
To validate abnormal chloride homeostasis in different diseases is associated with the abnormal transcriptional levels of GClTC from different clusters, we selected two widely known GClTC, SLC12A2 and SLC12A5, which encode proteins NKCC1 and KCC2, for further analysis [55, 56]. KCC2 is the major extruder of intracellular chloride, whereas NKCC1 mediates the influx of chloride ions . These proteins are essential for neural chloride homeostasis. Dysfunction of them could result in the pathogenesis of several brain disorders  or abnormal GABAergic signal during brain development . As reported in some research, SLC12A2 presents a relatively high expression level in early development stage in brain and decreases after birth . Interestingly, in temporal single-cell results (Additional file 9: Fig. S2), we found that SLC12A2 expression remained at a low level in neurons, and the slight upregulation was only observed mainly in oligodendrocytes after birth but was stable in neurons (Additional file 9: Fig. S2). Recent research also revealed that NKCC1 protein and mRNA are expressed at remarkably high levels in oligodendrocytes in mice . If dramatic changes in mRNA occur mainly in non-neuronal cells, how it regulates neuronal chloride homeostasis during development warrants further study. Moreover, based on our data, SLC12A2 and SLC12A5 had completely different expression patterns (Fig. 8). SLC12A5 has a pattern similar to that of GABA receptors; it is located in neurons and is mainly involved in regulating neural functions, consistent with previous reports [62,63,64]. In contrast, SLC12A2 mediates the function of non-neuronal cells. It is associated with myelin dysfunction, which has rarely been reported. Although, SLC12A2 is regarded as a target gene in epilepsy treatment [65, 66], we found no significant changes in SLC12A2 expression in epileptic models, but they seemed to act more sensitively in demyelinating diseases. Overall, we suggest SLC12A2 may play a more important role in oligodendrocytes than in neurons, and myelination disorders might be much more widespread than what has been thought so far.
In summary, although chloride homeostasis disorders occur in various neurological diseases, the GClTC are specific to certain types of disease. The findings of our current study are compelling, but they have some limitations. Not all the genes relevant to chloride homeostasis were included in our gene clusters. We only focused on the chloride transporters and channels, considering they are the direct executors to maintain intracellular chloride homeostasis. However, as reported by other study , genes like Na–K ATPase α2 (ATP1A2) and synaptophysin (SYP) also manifested similar expression pattern as SLC12A5, even though these genes are not chloride-regulatory genes. The lack of abundant postnatal data on critical periods (CPs) limits the interpretation of the results in those periods. It should not be ignored that the delicate regulation of post-translational modification during development can magnify the effects of protein functions even when the mRNA expression level is relatively low. Besides, even though we observed strong association between GClTC and diseases, it is still unclear whether the changes of gene expression are the cause of these diseases, or the changes are merely the adaptive responses due to the diseases. Therefore, we believe that more comprehensive studies on GClTC are required to elucidate the intricate mechanisms underlying brain diseases.
Availability of data and materials
Data will be made available on reasonable request.
- CNS :
Central nervous system
- ASD :
Autism spectrum disorder
- STAB :
Spatio-temporal cell atlas of the human brain
- MS :
- KA :
- IQR :
- CPs :
Pressey JC, de Saint-Rome M, Raveendran VA, Woodin MA. Chloride transporters controlling neuronal excitability. Physiol Rev. 2023;103(2):1095–135. https://doi.org/10.1152/physrev.00025.2021.
Kaila K, Price TJ, Payne JA, Puskarjov M, Voipio J. Cation-chloride cotransporters in neuronal development, plasticity and disease. Nat Rev Neurosci. 2014;15(10):637–54. https://doi.org/10.1038/nrn3819.
Blaesse P, Airaksinen MS, Rivera C, Kaila K. Cation-chloride cotransporters and neuronal function. Neuron. 2009;61(6):820–38. https://doi.org/10.1016/j.neuron.2009.03.003.
Virtanen MA, Uvarov P, Mavrovic M, Poncer JC, Kaila K. The multifaceted roles of KCC2 in cortical development. Trends Neurosci. 2021;44(5):378–92. https://doi.org/10.1016/j.tins.2021.01.004.
Peerboom C, Wierenga CJ. The postnatal GABA shift: a developmental perspective. Neurosci Biobehav Rev. 2021;124:179–92. https://doi.org/10.1016/j.neubiorev.2021.01.024.
Ge S, Pradhan DA, Ming G-L, Song H. GABA sets the tempo for activity-dependent adult neurogenesis. Trends Neurosci. 2007;30(1):1–8.
Giachino C, Barz M, Tchorz JS, Tome M, Gassmann M, Bischofberger J, et al. GABA suppresses neurogenesis in the adult hippocampus through GABAB receptors. Development. 2014;141(1):83–90. https://doi.org/10.1242/dev.102608.
Sedmak G, Jovanov-Milosevic N, Puskarjov M, Ulamec M, Kruslin B, Kaila K, et al. Developmental expression patterns of KCC2 and functionally associated molecules in the human brain. Cereb Cortex. 2016;26(12):4574–89. https://doi.org/10.1093/cercor/bhv218.
Medina I, Chudotvorova I. GABA neurotransmission and neural cation-chloride co-transporters: actions beyond ion transport. Crit Rev Neurobiol. 2006;18(1–2):105–12.
Singh AK, Venglarik CJ, Bridges RJ. Development of chloride channel modulators. Kidney Int. 1995;48(4):985–93.
De Koninck Y. Altered chloride homeostasis in neurological disorders: a new target. Curr Opin Pharmacol. 2007;7(1):93–9.
Raveendran VA, Pressey JC, Woodin MA. A novel small molecule targets NKCC1 To restore synaptic inhibition. Trends Pharmacol Sci. 2020;41(12):897–9. https://doi.org/10.1016/j.tips.2020.10.002.
Merlin LR. Chloride’s exciting role in neonatal seizures suggests novel therapeutic approach. Epilepsy Curr. 2011;11(3):92–3. https://doi.org/10.5698/1535-7511-11.3.92.
Kalkman HO. Alterations in the expression of neuronal chloride transporters may contribute to schizophrenia. Prog Neuropsychopharmacol Biol Psychiatry. 2011;35(2):410–4. https://doi.org/10.1016/j.pnpbp.2011.01.004.
Doyon N, Vinay L, Prescott SA, De Koninck Y. Chloride regulation: a dynamic equilibrium crucial for synaptic inhibition. Neuron. 2016;89(6):1157–72. https://doi.org/10.1016/j.neuron.2016.02.030.
Miller JA, Ding S-L, Sunkin SM, Smith KA, Ng L, Szafer A, et al. Transcriptional landscape of the prenatal human brain. Nature. 2014;508(7495):199–206. https://doi.org/10.1038/nature13185.
Hodge RD, Bakken TE, Miller JA, Smith KA, Barkan ER, Graybuck LT, et al. Conserved cell types with divergent features in human versus mouse cortex. Nature. 2019;573(7772):61–8. https://doi.org/10.1038/s41586-019-1506-7.
Song L, Pan S, Zhang Z, Jia L, Chen W-H, Zhao X-M. STAB: a spatio-temporal cell atlas of the human brain. Nucleic Acids Res. 2021;49(D1):D1029–37. https://doi.org/10.1093/nar/gkaa762.
Chen C, Chen H, Zhang Y, Thomas HR, Frank MH, He Y, et al. TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol Plant. 2020;13(8):1194–202. https://doi.org/10.1016/j.molp.2020.06.009.
Ernst J, Bar-Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinform. 2006;7:191.
Piñero J, Ramírez-Anguita JM, Saüch-Pitarch J, Ronzano F, Centeno E, Sanz F, et al. The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 2020;48(D1):D845–55. https://doi.org/10.1093/nar/gkz1021.
Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature. 2012;489(7416):391–9. https://doi.org/10.1038/nature11405.
Yarkoni T, Poldrack RA, Nichols TE, Van Essen DC, Wager TD. Large-scale automated synthesis of human functional neuroimaging data. Nat Methods. 2011;8(8):665–70. https://doi.org/10.1038/nmeth.1635.
Hendrickx DAE, van Scheppingen J, van der Poel M, Bossers K, Schuurman KG, van Eden CG, et al. Gene expression profiling of multiple sclerosis pathology identifies early patterns of demyelination surrounding chronic active lesions. Front Immunol. 2017;8:1810. https://doi.org/10.3389/fimmu.2017.01810.
Conte G, Parras A, Alves M, Ollà I, De Diego-Garcia L, Beamer E, et al. High concordance between hippocampal transcriptome of the mouse intra-amygdala kainic acid model and human temporal lobe epilepsy. Epilepsia. 2020;61(12):2795–810. https://doi.org/10.1111/epi.16714.
Kang HJ, Kawasawa YI, Cheng F, Zhu Y, Xu X, Li M, et al. Spatio-temporal transcriptome of the human brain. Nature. 2011;478(7370):483–9. https://doi.org/10.1038/nature10523.
Zhong S, Zhang S, Fan X, Wu Q, Yan L, Dong J, et al. A single-cell RNA-seq survey of the developmental landscape of the human prefrontal cortex. Nature. 2018;555(7697):524–8. https://doi.org/10.1038/nature25980.
Polioudakis D, de la Torre-Ubieta L, Langerman J, Elkins AG, Shi X, Stein JL, et al. A single-cell transcriptomic atlas of human neocortical development during mid-gestation. Neuron. 2019. https://doi.org/10.1016/j.neuron.2019.06.011.
Skutella T, Nitsch R. New molecules for hippocampal development. Trends Neurosci. 2001;24(2):107–13.
Gene OC. Gene Ontology Consortium: going forward. Nucleic Acids Res. 2015;43(Database issue):D1049–56. https://doi.org/10.1093/nar/gku1179.
Zhang S, Zhou J, Zhang Y, Liu T, Friedel P, Zhuo W, et al. The structural basis of function and regulation of neuronal cotransporters NKCC1 and KCC2. Commun Biol. 2021;4(1):226. https://doi.org/10.1038/s42003-021-01750-w.
Boscia F, Elkjaer ML, Illes Z, Kukley M. Altered expression of ion channels in white matter lesions of progressive multiple sclerosis: what do we know about their function? Front Cell Neurosci. 2021;15: 685703. https://doi.org/10.3389/fncel.2021.685703.
Auer T, Schreppel P, Erker T, Schwarzer C. Impaired chloride homeostasis in epilepsy: molecular basis, impact on treatment, and current treatment approaches. Pharmacol Ther. 2020;205: 107422. https://doi.org/10.1016/j.pharmthera.2019.107422.
Reh RK, Dias BG, Nelson CA, Kaufer D, Werker JF, Kolb B, et al. Critical period regulation across multiple timescales. Proc Natl Acad Sci U S A. 2020;117(38):23242–51. https://doi.org/10.1073/pnas.1820836117.
Silva CG, Peyre E, Adhikari MH, Tielens S, Tanco S, Van Damme P, et al. Cell-intrinsic control of interneuron migration drives cortical morphogenesis. Cell. 2018. https://doi.org/10.1016/j.cell.2018.01.031.
Ma D, Przybylski D, Abruzzi KC, Schlichting M, Li Q, Long X, et al. A transcriptomic taxonomy of circadian neurons around the clock. Elife. 2021. https://doi.org/10.7554/eLife.63056.
Monroe JG, Powell T, Price N, Mullen JL, Howard A, Evans K, et al. Drought adaptation in by extensive genetic loss-of-function. Elife. 2018. https://doi.org/10.7554/eLife.41038.
Di Paolo G, De Camilli P. Phosphoinositides in cell regulation and membrane dynamics. Nature. 2006;443(7112):651–7.
Pacheco MA, Jope RS. Phosphoinositide signaling in human brain. Prog Neurobiol. 1996;50(2–3):255–73.
Wang TA, Chen C, Huang F, Feng S, Tien J, Braz JM, et al. TMEM16C is involved in thermoregulation and protects rodent pups from febrile seizures. Proc Natl Acad Sci U S A. 2021. https://doi.org/10.1073/pnas.2023342118.
Feenstra B, Pasternak B, Geller F, Carstensen L, Wang T, Huang F, et al. Common variants associated with general and MMR vaccine-related febrile seizures. Nat Genet. 2014;46(12):1274–82. https://doi.org/10.1038/ng.3129.
Mount DB, Mercado A, Song L, Xu J, George AL, Delpire E, et al. Cloning and characterization of KCC3 and KCC4, new members of the cation-chloride cotransporter gene family. J Biol Chem. 1999;274(23):16355–62.
Concepcion AR, Wagner LE, Zhu J, Tao AY, Yang J, Khodadadi-Jamayran A, et al. The volume-regulated anion channel LRRC8C suppresses T cell function by regulating cyclic dinucleotide transport and STING-p53 signaling. Nat Immunol. 2022;23(2):287–302. https://doi.org/10.1038/s41590-021-01105-x.
Yang J, Vitery MDC, Chen J, Osei-Owusu J, Chu J, Qiu Z. Glutamate-releasing SWELL1 channel in astrocytes modulates synaptic transmission and promotes brain damage in stroke. Neuron. 2019. https://doi.org/10.1016/j.neuron.2019.03.029.
Huisman SMH, van Lew B, Mahfouz A, Pezzotti N, Höllt T, Michielsen L, et al. BrainScope: interactive visual exploration of the spatial and temporal human brain transcriptome. Nucleic Acids Res. 2017;45(10): e83. https://doi.org/10.1093/nar/gkx046.
Fritschy JM, Mohler H. GABAA-receptor heterogeneity in the adult rat brain: differential regional and cellular distribution of seven major subunits. J Comp Neurol. 1995;359(1):154–94.
Nuwer JL, Brady ML, Povysheva NV, Coyne A, Jacob TC. Sustained treatment with an α5 GABA A receptor negative allosteric modulator delays excitatory circuit development while maintaining GABAergic neurotransmission. Neuropharmacology. 2021;197: 108724. https://doi.org/10.1016/j.neuropharm.2021.108724.
Fischell J, Van Dyke AM, Kvarta MD, LeGates TA, Thompson SM. Rapid antidepressant action and restoration of excitatory synaptic strength after chronic stress by negative modulators of alpha5-containing GABAA receptors. Neuropsychopharmacology. 2015;40(11):2499–509. https://doi.org/10.1038/npp.2015.112.
Ahring PK, Liao VWY, Gardella E, Johannesen KM, Krey I, Selmer KK, et al. Gain-of-function variants in GABRD reveal a novel pathway for neurodevelopmental disorders and epilepsy. Brain. 2021. https://doi.org/10.1093/brain/awab391.
Whissell PD, Avramescu S, Wang D-S, Orser BA. δGABAA receptors are necessary for synaptic plasticity in the hippocampus: implications for memory behavior. Anesth Analg. 2016;123(5):1247–52.
Bertram EH. Neuronal circuits in epilepsy: do they matter? Exp Neurol. 2013;244:67–74. https://doi.org/10.1016/j.expneurol.2012.01.028.
Wang Y, Chen Z. An update for epilepsy research and antiepileptic drug development: toward precise circuit therapy. Pharmacol Ther. 2019;201:77–93. https://doi.org/10.1016/j.pharmthera.2019.05.010.
Keezer MR, Sisodiya SM, Sander JW. Comorbidities of epilepsy: current concepts and future perspectives. Lancet Neurol. 2016;15(1):106–15. https://doi.org/10.1016/S1474-4422(15)00225-2.
May P, Girard S, Harrer M, Bobbili DR, Schubert J, Wolking S, et al. Rare coding variants in genes encoding GABA receptors in genetic generalised epilepsies: an exome-based case-control study. Lancet Neurol. 2018;17(8):699–708. https://doi.org/10.1016/S1474-4422(18)30215-1.
Koumangoye R, Bastarache L, Delpire E. NKCC1: newly found as a human disease-causing ion transporter. Function (Oxf). 2021;2(1):zqaa028. https://doi.org/10.1093/function/zqaa028.
Stödberg T, McTague A, Ruiz AJ, Hirata H, Zhen J, Long P, et al. Mutations in SLC12A5 in epilepsy of infancy with migrating focal seizures. Nat Commun. 2015;6:8038. https://doi.org/10.1038/ncomms9038.
Nalamalapu RR, Yue M, Stone AR, Murphy S, Saha MS. The gene family: from embryo to disease. Front Mol Neurosci. 2021;14: 672511. https://doi.org/10.3389/fnmol.2021.672511.
Huang H, Song S, Banerjee S, Jiang T, Zhang J, Kahle KT, et al. The WNK-SPAK/OSR1 kinases and the cation-chloride cotransporters as therapeutic targets for neurological diseases. Aging Dis. 2019;10(3):626–36. https://doi.org/10.14336/AD.2018.0928.
Watanabe M, Fukuda A. Development and regulation of chloride homeostasis in the central nervous system. Front Cell Neurosci. 2015;9:371. https://doi.org/10.3389/fncel.2015.00371.
Ben-Ari Y. NKCC1 chloride importer antagonists attenuate many neurological and psychiatric disorders. Trends Neurosci. 2017;40(9):536–54. https://doi.org/10.1016/j.tins.2017.07.001.
Kurki SN, Uvarov P, Pospelov AS, Trontti K, Hübner AK, Srinivasan R, et al. Expression patterns of NKCC1 in neurons and non-neuronal cells during cortico-hippocampal development. Cereb Cortex. 2022. https://doi.org/10.1093/cercor/bhac470.
Sedmak G, Jovanov-Milošević N, Puskarjov M, Ulamec M, Krušlin B, Kaila K, et al. Developmental expression patterns of KCC2 and functionally associated molecules in the human brain. Cereb Cortex. 2016;26(12):4574–89. https://doi.org/10.1093/cercor/bhv218.
Vu TQ, Payne JA, Copenhagen DR. Localization and developmental expression patterns of the neuronal K-Cl cotransporter (KCC2) in the rat retina. J Neurosci. 2000;20(4):1414–23.
Coull JAM, Boudreau D, Bachand K, Prescott SA, Nault F, Sík A, et al. Trans-synaptic shift in anion gradient in spinal lamina I neurons as a mechanism of neuropathic pain. Nature. 2003;424(6951):938–42.
Dzhala VI, Kuchibhotla KV, Glykys JC, Kahle KT, Swiercz WB, Feng G, et al. Progressive NKCC1-dependent neuronal chloride accumulation during neonatal seizures. J Neurosci. 2010;30(35):11745–61. https://doi.org/10.1523/JNEUROSCI.1769-10.2010.
Hampel P, Johne M, Gailus B, Vogel A, Schidlitzki A, Gericke B, et al. Deletion of the Na-K-2Cl cotransporter NKCC1 results in a more severe epileptic phenotype in the intrahippocampal kainate mouse model of temporal lobe epilepsy. Neurobiol Dis. 2021;152: 105297. https://doi.org/10.1016/j.nbd.2021.105297.
This work was supported by the National Natural Science Foundation, People's Republic of China (Grant Nos. 81730031 to YW Wang, Grant Nos. 82001455 to WJ Du, Grant Nos. 32200777 to ZJ Hao) and the Foundation of Shanghai Municipal Key Clinical Specialty (shslczdzk06901 to YW Wang).
Ethics approval and consent to participate
Consent for publication
All authors have revised and have agreed with the content of the manuscript. Also, all authors have given explicit consent to submit.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Definitions of the periods of human development.
Additional file 2: Table S2.
The time series analysis of genes in cluster 1 ~ 10.
Additional file 3: Table S3.
List of 104 GClTC.
Additional file 4: Table S4.
The top 5 significant terms enriched from cluster 1 ~ 10.
Additional file 5: Table S5.
Distribution of GClTC in brain regions.
Additional file 6: Table S6.
List of clusters 4, 6, 8, and 9 annotated through gene-disease association analysis.
Additional file 7: Table S7.
The normality and variance homogeneity for data in Fig. 7.
Additional file 8: Fig. S1.
Molecular characteristics of the developmental expression patterns for 1–10 cluster.
Additional file 9: Fig. S2.
Gene expression dynamics across cell subtypes of the selected brain region.
Additional file 10: Fig. S3.
Distribution of neural GClTC in different brain regions.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.
About this article
Cite this article
Huang, Y., Wang, Q., Peng, Y. et al. Spatiotemporal expression patterns of genes coding for plasmalemmal chloride transporters and channels in neurological diseases. Mol Brain 16, 30 (2023). https://doi.org/10.1186/s13041-023-01018-w
- Gene expression pattern
- Chloride transporter
- Chloride channel