Introduction

After birth, the brain undergoes developmental changes for a prolonged time that involve a series of complex and accurately orchestrated processes (Rakic 2009). The production and migration of neurons is largely complete at the beginning of postnatal development, and then the intrauterine developmental processes continue: gray and white matter thickening, myelination, synaptogenesis, pruning, and establishment of the basic anatomical architecture for initial neural pathway function. Subsequently, local connections within cortical circuits are fine-tuned, and increasingly complex, longer-term connections are established between circuits (Stiles 2010). After that, the changes do not end, but continue throughout human life. They are mostly driven by environmental influences and experiences and lead to changes in metabolic activities (Kuzawa et al. 2014), changes in functional connectivity patterns (Kelly et al. 2009) and, with the maturation of white matter (Beck et al 2021, Yeatman et al 2014) changes in the speed of long-distance transmission (van Blooijs et al. 2023). The final phase is aging, where it slowly declines with advancing age, leading to a decline in cognitive signal processing functions and often resulting in neurodegenerative diseases (Peters 2006). The cortical supragranular glutamatergic cell (or pyramidal cell) provides the excitatory synaptic inputs for local inhibitory circuitry and other pyramidal cells by which they create distinct subnetworks (Yoshimura 2005). The development and formation of dendrites (Petanjek 2008, Koenderink 1995) and synapses (Huttenlocher and Dabholkar 1997) of pyramidal cells in the human cerebral cortex has been documented to some extent by postmortem studies, besides much less is known about their biophysical maturation and electrical properties in the early stages of development and the subsequent change or maintenance in later ages. Numerous studies demonstrated in non-primate animal models that the electrical characteristics of neurons change prominently in the early postnatal stage (Molnár et al. 2020). Changes in the intrinsic membrane properties (Picken Bahrey and Moody 2003), the input resistance or the kinetics of the elicited action potentials were reported (Kroon 2019, Elston and Fujita 2014) in connection with maturation of macaque and rodent pyramidal cells. To date, however, no cross-age studies have been conducted on the electrophysiological parameters of human pyramidal neurons. We have studied in detail the postnatal lifetime profile of the physiological and morphological properties of supragranular (layer 2/3) neurons of human pyramidal cells from neurosurgical resections. To this end, we performed whole-cell patch-clamp recordings and 3D anatomical reconstructions of human cortical pyramidal cells from 109 patients aged 1 month to 85 years for comprehensive data analysis to obtain the morphoelectric lifetime profile of supragranular pyramidal cells.

Results

Age-dependent differences in intrinsic subthreshold membrane properties

To extract biophysical properties of excitatory cells of human brain specimens we performed whole cell patch-clamp recordings of pyramidal cells from neurosurgically removed human neocortical tissue sections. The samples were mainly from the frontal and temporal lobes (Figure 1D, Supplementary Figure 1D), mostly from patients with tumors or hydrocephalus (Supplementary Figure 1B). Data were collected from 109 patients aged 1 month to 85 years (Figure 1B-C) from 498 human cortical layer 2/3 (L2/3) pyramidal cells. To confirm that the studied pyramidal cells originate from the L2/3, we measured the distance between the cell body and the L1 border (Berg 2021). 36% of the cells recovered their soma, with a distance of 129.69 ± 130.77 μm from the L1 border (Fig 1E). The data set was divided into seven age groups (Kang 2011, Bethlehem 2022): infant: < 1 year, early childhood: 1-6 years, late childhood: 7-12 years, adolescence: 13-19 years, young adulthood: 20-39 years, middle adulthood: 40-59 years, late adulthood: ≥ 60 years (Figure 1A).

Figure 1.

We evaluated the voltage deflections induced by negative and positive current injections and extracted subthreshold membrane features such as resting membrane potential, input resistance, time constant (tau), and sag ratio in 457 cortical pyramidal cells from 99 patients. We found that the subthreshold features from samples of infant significantly different from those from other age groups (Figure 2A-D; resting membrane potential: P =3.53×10^-8, input resistance: P = 1.29×10^-16, tau: P = 1.31×10^-15, sag ratio: P = 5.2×10^-4, Kruskal-Wallis test) (Supplementary Figure 2). The resting membrane potential was significantly more positive in the first year of life than in the other age groups. Before adulthood a slight decrease was observed in the resting membrane potential across the groups (Table 1; Figure 2A). Input resistance, tau and sag ratio were measured on voltage deflections elicited by injecting negative (-100 pA) current steps into the cells. A significant decrease in input resistance (Table 1) was observed with the largest reduction after the first year of age (P = 1.88×10^-6, Kruskal–Wallis test with post-hoc Dunn test) (Figure 2B). The membrane time constant also decreased significantly in the older groups compared to the infant group, after infancy we found more conserved mean values of membrane time constant into older age (Table 1; Figure 2C). The ratio of the maximal deflection and the steady-state membrane potential during a negative current step (sag ratio) is significantly higher in late adulthood than in the early stages of life (Table 1; Figure 2D). Note that the high variance of the infant group data (e.g. resting membrane potential, input resistance, tau) are due to the dynamic change over the 0-1 year period (Supplementary Figure 3).

Figure 2.

Table 1.

Suprathreshold properties across age-groups

To initiate action potentials we injected positive current steps increased by 20 pA into the cells and recorded various types of input-output transformations. We extracted 25 features from different action potentials (AP) and firing patterns and assessed active membrane properties from recordings filtered for appropriate electrophysiological quality (see Methods). We found that the infant group differs most from the other age groups in several of the suprathreshold properties (Figure 3), but other trends are also apparent (rheobase current: P = 8.71×10^-12, AP half-width: P = 9.57×10^-25, AP up-stroke velocity: P = 1.63×10^-12, AP amplitude: P = 2.24×10^-11, Kruskal-Wallis test)(Supplementary Figure 4). For example, the average rheobase current, the minimum current that can trigger an action potential, is significantly lower in the early ages of life than in data collected from adolescence stage (Table 2) (infant vs. adolescence P = 4.15×10^-13,Kruskal–Wallis test with post-hoc Dunn test). Further on the age scale we found that rheobase current was increased with age reaching a maximum value at the adolescent ages then declining to a significantly lower level (adolescent vs. late adulthood, P = 5.23×10^-4, Kruskal–Wallis test with post-hoc Dunn test, Figure. 3A). The action potential (AP) half-width averages among a declining trend through the groups of age (Table 2) varied considerably across age groups forming significant differences between childhood and adulthood ages (Figure 3B). Action potential up-stroke velocities were significantly slower in the infant APs than in all the other ages (Table 2; Figure 3C). The amplitude of the elicited APs also showed age-dependent differences between age groups (Table 2); in the first year of life, the amplitude of APs were significantly lower than in other age groups. In adulthood, a significant decrease with age was observed (Figure 3D). Electrophysiological differences markedly separated the infant group from the older age groups shown on UMAP projection (Figure 3F, Supplementary Figure 5).

Figure 3.

Table 2.

Table 3.

Table 4.

Next, we investigated how somatic current inputs are transformed to action potential output by evaluating the firing patterns of cells evoked by injecting prolonged positive current steps (Figure 4A) (Supplementary Figure 6). The neurons (n = 331) were regular spiking cells with moderate adaptation. The slope of the firing frequency versus the current curve (f-I slope) showed no significant difference across the age groups (P = 0.055, Kruskal-Wallis test) (Table 2; Figure 4B). Age related difference was observed in the latency of the first AP during rheobase current injection (First AP latency), in the first year of life the latency of the first spike is significantly higher than all the other groups of age (P = 7.67×10^-4, Kruskal-Wallis test), which is the most prominent with the oldest age group (P = 8.41×10^-6, Kruskal–Wallis test with post-hoc Dunn test) (Table 2; Figure 4C). The adaptation of the AP frequency response to the same current injection stimulus also showed differences (P = 0.032, Kruskal-Wallis test) between the younger and the groups older than 13 years patients, we found the lowest adaptation values in early childhood (Table 2; Figure 4D).

Figure 4.

Previous research has shown that the biophysical properties of the human pyramidal cells show depth-related correlations throughout the L2/3 (Berg 2021, Kalmbach 2018, Chameh 2021). Although there are some deeper cells in our dataset, the majority comes from the upper region of the L2/3. We compared the electrophysiological characteristics according to their depth from the border of L1 and L2 to exclude the possibility that the biophysical differences we found were a result of depth dependence. We did not find any overall differences related to distance of the soma from the L1 border within the age groups with a few exceptions. For example, the values of input resistance (P = 0.02, Mann-Whitney test) and AP up-stroke velocity (P = 0.04, Mann-Whitney test) differ significantly in the middle adulthood group. We found a significant difference in AP amplitude (P = 0.02, Mann-Whitney test) and adaptation (P = 0.009, Mann-Whitney test) in the adolescence age group (Supplementary Figure 7-9).

Morphological features of layer 2/3 pyramidal cells in different stages of life

To investigate possible morphological differences between the age groups, we filled the pyramidal cells with biocytin during recordings. Only neurons with no signs of deterioration and with complete apical dendrites and no signs of truncated dendritic branches or tufts were considered for morphological analysis. 63 pyramidal cells (Supplementary Figure 10) were reconstructed in 3D at ages 0-73 years (infant n = 7, early childhood n = 8, late childhood n = 11, adolescence n = 11, young adulthood n = 9, middle adulthood n = 9, late adulthood n = 8)(Supplementary Figures 11-12). Figure 5A shows examples of the reconstructed pyramidal cells. We did not detect significant change in total dendritic length (P = 0.37, Kruskal-Wallis test) (Table 3; Figure 5B), apical dendritic length (P = 0.6, Kruskal-Wallis test) (Table 3; Figure 5C), or basal dendritic length (P = 0.28, Kruskal-Wallis test) (Table 3; Figure 5D) at different ages. To investigate dendritic complexity we measured the total number of dendritic branching and found no significant developmental changes (P = 0.18, Kruskal-Wallis test) (Table 3; Figure 5E). We also did not observe significant differences in the size of dendritic branching when we measured the maximum horizontal (P = 0.64, Kruskal-Wallis test) (Table 3; Figure 5F) and vertical (P = 0.51, Kruskal-Wallis test) (Table 3; Figure 5G) extent of the reconstructed cells. We found significant differences in the average length of the uncut terminal segments of apical (P = 0.033, Kruskal-Wallis test) (Table 3; Figure 5H) but not of the basal dendrites (P = 0.85, Kruskal-Wallis test) (Table 3; Figure 5I) of the cells across the age groups.

Figure 5.

To analyze the distribution of dendritic spines we identified and labeled each spine on n = 6 fully 3D-reconstructed cells (Figure 6). We compared the spine density of selected pyramidal cells of two age groups: infant (83 days old, n = 3 of one patient, parietal lobe) vs. late adulthood (64.3 ± 2.08 years old, n = 3 of 3 patients, frontal, temporal and parietal lobes). The investigated cells are located in L2 (infant: 144.43 ± 45.26 µm, late adulthood: 161.22 ± 66.22 µm). We found that the total spine density was higher (P = 7.57×10^-40, Mann-Whitney test) and also the spine density of both apical (P = 2.02×10^-31, Mann-Whitney test) and basal (P = 3.8×10^-12, Mann-Whitney test) dendrites was higher in the infant than in the late adult group (Figure 6B). To evaluate the age-dependence of spine morphology, we classified the spines into commonly used phenotypes based on their morphological characteristics (Li 2023), specifically distinguishing between mushroom-shaped, thin, filopodial, branched, and stubby spines (Figure 6C-G). Spines with large spine heads were classified as mushroom-shaped, those with small heads as thin, and long protrusions were distinguished as filopodial. Those that did not have peduncles were classified as stubby, spines with two heads emerging from the same spot were called branched spines (Figure 6C-G) (Luebke 2015). Only fully visible spines were included in the classification analysis. The composition of spine types varies between the two age groups, mushrooms are present at a higher percentage on apical branches in late adulthood cells (P = 4.4×10^-9, Mann-Whitney test), and on the basal processes (P = 9.04×10^-8, Mann-Whitney test) (Figure 6C,H). In contrast, thin spines and filopodia are present in significantly higher numbers on the apical (Figure 6D-E) (thin spines: P = 7.34×10^-14; filopodia: P = 1.11×10^-39, Mann-Whitney test) and basal (thin: P = 2.46×10^-8; filopodia: P = 2.14×10^-12, Mann-Whitney test) (Figure 6I-J) dendritic branches of the infant. Both apical and basal infant branches had a significantly higher percentage of branched spines (apical: P = 1.64×10^-11; basal: P = 8.9×10^-5, Mann-Whitney test) (Figure 6F,K), which were present in modest numbers on both. Stubby spines were also more prevalent on the elderly pyramidal cells, either on the apical (P = 7.19×10^-5, Mann-Whitney test) or basal (P = 6.97×10^-9, Mann-Whitney test) processes (Figure 6G,L). Comparing the spine density of the six individual cells, we found differences across the two age groups similar to those previously mentioned, alongside slight variation within age groups.(Supplementary Figure 13, Kruskal-Wallis test with post-hoc Dunn test).

Figure 6.

Cortical tissue was dissected during neurosurgical procedures for various pathologies, but not to the same extent in the different groups. Specimen collection in adult/elderly patients was mostly for tumor resection, whereas in children hydrocephalus was the most common reason for brain surgery (Supplementary Figure 1). To determine the difference in medical condition could contribute to age-related differences in cellular features, we compared the extracted electrophysiological and morphological features based on the medical condition in the different age groups. It should be noted that all circumstances of surgical steps, tissue dissection, transport (time, media, temperature etc.) and cutting procedure remained the same in the different conditions. When comparing passive electrophysiological properties, we found no significant differences between the different medical conditions in the age groups from infancy to middle adulthood, only time constants were found to be significantly lower in the hydrocephalus patients than in the tumor patients in the young adulthood (P = 0.048, Mann-Whitney test) and late adulthood groups (P = 0.01, Mann-Whitney test) (Supplementary Figure 14). Comparing the action potential kinetics and firing pattern-related parameters between the different pathology groups we found no overall differences with some sporadic exceptions. For example, the parameters of AP half-width in young adulthood (P = 0.01, Mann-Whitney test) and the values of rheobase in early childhood (P = 0.04, Mann-Whitney test) differed significantly between the different pathologies (Supplementary Figure 15). Regarding the firing pattern-related properties we found significant changes in the infant F-I slope (P = 0.02, Mann-Whitney test), first AP latency (P = 0.006, two sample t-test) and adaptation (P = 0.002, Mann-Whitney test) parameters. Also in the late adulthood group the F-I slope (P = 0.04, Mann-Whitney test), the first AP latency (P = 0.02, Mann-Whitney test) and the adaptation (P = 0.047, Mann-Whitney test) were significantly different between the hydrocephalus and tumor patients (Supplementary Figure 16). Further comparison of morphological features revealed no statistical difference between tumor and hydrocephalus groups (Supplementary Figure 17).

To ascertain whether age-related variations in cell properties may be influenced by the gender of the patient, we compaired the examined characteristics of the cells by their sex (Supplementary Figure. 18-21). Comparing the passive membrane properties of the cells in different age groups we found higher resting membrane potential values in adolescent females than in males (P = 0.04, Two-sample t-test) and input resistance was lower in female patients in the late childhood group (P = 0.03, Mann-Whitney test) (Supplementary Figure 18). We found differences when comparing the AP characteristics in AP half-with cells from male patient showed higher values in early childhood (P = 0.036) and adolescence (P = 0.02). AP up-stoke velocity was higer in young adulthood females than in males (P = 0.34). AP amplitude was significantly higher in male patients in early childhood (P = 0.003) and late adulthood (P = 0.035) age groups (Supplementary Figure 19). Under the comparison of the firing pattern features cells from the infancy showed significant differences in F-I slope (P = 0.2, Mann-Whitney test), first AP latency (P = 0.0096, Two-sample t-test) and adaptation (P = 0.002, Mann-whitney test). The latency of the first AP was higher in cells from male patient in the young (P = 0.03, Mann-Whitney test) and middle (P = 0.018, Mann-Whitney test) adulthood cells. The adaptation of the APs in males was higher during adolescence (P = 0.027, Mann-Whitney test) but lower in middle adulthood (P = 0.034, Mann-Whitney test)(Supplementary Figure 20). Regarding the morphological property of the reconstructed cell we found differences in the young adulthood age group between cells from male and female patients in total (P = 0.03, Mann-Whitney test) and apical (P = 0.03, Mann-Whitney test) dendritic length and the total number of bifurcations (P = 0.03, Mann-Whitney test), in the middle adulthood group we found significant difference in the average apical terminal length across males females(P = 0.03, Mann-Whitney test)(Supplementary Figure 21).

Discussion

In this study, using the whole-cell patch-clamp technique and 3D reconstructions, we have studied the differences of human cortical L2/3 pyramidal cells at various stages of life. We found that during the lifespan the most significant changes take place early in life during the first year in most of the biophysiological characteristics of human pyramidal cells (see Supplementary Figure 3). We recorded from n = 457 human cortical excitatory pyramidal cells from the supragranular layer from birth to 85 years. Most differences in sub- and suprathreshold features were found in the youngest age groups. There were particular differences in resting membrane potential, input resistance, time constant, rheobasic current, AP halfwidth and AP upstroke velocity. Other age groups that differed the most from other age groups were the oldest (60 - 85 years) with modest differences in sag ratios. In our morphological analysis, we found no significant changes in the overall apical and basal dendritic dimensions across different ages. When evaluating the number of spines in two anatomical reconstructions (from the infant and late adulthood groups) we found a higher overall spine density in the infant than in the elder sample consisting of mainly branched and thin spines and filopodia. In contrast, stubby and mushroom-shaped spines were more prevalent in the older pyramidal cell.

The nervous system experiences a multitude of plastic changes throughout an individual’s lifetime. These changes extend from the early stages of development and continue through the gradual degenerative processes that come with old age. Although precise intrinsic cellular modifications linked to aging remain somewhat elusive, there is an observable morphological transformation during early development. As the nervous system matures, dendritic branching (Kasper et al. 1994, McAllister 2000) and cell body size increases (Zhang 2004). Also, learning and experience increase the extent of dendritic branching (Holloway 1966, Greenough and Volkmar 1973) and dendritic spines (Globus 1973). An extensive body of research has shown that synaptic plasticity induces formation of new dendritic spines and the enlargement of existing spines (Kastellakis 2023). These modifications may affect the input-output relation of a neuron as model studies suggest that dendritic geometry influences electrical and firing properties (Mainen and Sejnowski 1996, van Elburg and van Ooyen 2010, Eyal 2014, Vetter et al 2001). Furthermore, in parallel to the changes in morphological complexity there is ample experimental evidence that the maturation of ion channels and the change in increasing channel density in the neuronal membrane also play a role in influencing the electrical properties of neurons (Huguenard et al. 1988, Connors 1994). These findings point to a close causal relationship between structure and electrophysiological properties at cellular level.

Lifespan changes in electrophysiological properties

The intrinsic electrophysiological properties of a neuron are determined by several factors: the resistivity of the cytoplasm and membrane, the membrane capacitance, and the shape of the soma and dendrites. The basic passive parameters we recorded, such as the resting membrane potential, the input resistance, and the membrane time constant show a dramatical change in the first year. This change is consistent with data found in previous studies on neocortical pyramidal cells from rodent models (Kasper 1994, Kroon 2019, Picken Bahrey and Moody 2003, Zhang 2004, Zhu 2000, McCormick 1987) and on xenotransplated human cortical neurons (Linaro 2019). We found a progressive shift of the resting membrane potential from -60 mV to -68 mV in the first part of the lifespan (<40 years), which could be due to the age-dependent change in ion channel composition (Picken Bahrey and Moody 2003). The input resistance and time constant showed a sudden and sharp decrease between the first and second age groups, indicating the most dramatic changes in the infant brain followed by a generalized and slower decrease. As pyramidal cells mature, the size and volume of the cell body increase (Petanjek 2019, Zhu 2000) which correlates with the input resistance, meaning that a small membrane area has a higher input resistance (Luhmann 2000). In addition to the increase in cell surface area, the increase in the expression of certain potassium leak channels, e.g. from the two-pore domain potassium channel (KCNK) family, can also influence the input resistance during cell maturation (Aller 2008, Goldstein 2001). Similar reasons can lead to a decrease in the membrane time constant. The membrane time constant is the product of the specific input resistance and the specific membrane capacitance. Assuming constant membrane capacitance (∼1 pF/cm², Oláh et al 2024, Beaulieu-Laroche 2018), the decrease is likely also due to a change in specific membrane resistance. The result of developmental changes in membrane passive parameters during early development improves spike timing precision and temporal responsiveness in the cortical microcircuit (Oswald and Reyes 2011, Doischer et al. 2008). However, in our human samples at later ages (>40 years) that we associate with the aging brain, the previous trends of changes in intrinsic characteristics are reversed. Input resistance and time constant increase, which may be associated with the decline in KCNK expression in the aging brain (Erraji-Benchenkroun 2005). These membrane parameters are also generally higher in studies in rhesus monkeys (Moore 2023), except for the resting membrane potential, which is not different from that of young cortical L2/3 (Chang 2005) or L5 pyramidal cells (Luebke 2007) in the brain. However, in our samples, we found a positive shift in the resting potential between the middle- and old-age groups by a few millivolts. It should be noted that this discrepancy may be due to the fact that the age groups in our work are not the same as the works mentioned.

Changes in resting membrane potential and input resistance will add up to produce a rising and then falling curve of the averaged rheobasic current, with the value three times higher in the adolescent group than in the youngest group and 1.5 times higher than in late adulthood. Our results on increased rheobase during early life are consistent with previous studies in rodent cortical L3 and 5 pyramidal cells (Kroon 2019, Perez-García 2021, Popescu 2021) and with findings regarding aging in rhesus monkey cortical neurons from rhesus monkeys (Moore 2023).

When we assessed the voltage sag ratio with age we found a consistent increase that peaked in late adulthood. The voltage sag that occurs during membrane hyperpolarization is the result of activation of the hyperpolarization-activated cyclic nucleotide-gated (HCN) current, also known as the h-current. Similar to our findings the age-dependent increase in sag deflection has been observed in human cortical pyramidal cells (Guet-McCreight 2023, Kalmbach 2018) which is consistent with studies by Wang et al. (2011, 2007) that have linked inhibited HCN activation to working memory deficits in the elderly.

The electrical output of neurons is driven by the composition of neuronal molecules and ion channels whose expressions and distributions change most dynamically in two periods, early development, and the aging brain. This dynamic essentially produced an inverted U-shaped curve in the distribution of the active electrical properties over the lifespan in our dataset, an increase in the first year, a relatively stable middle phase, and then a return in older age groups. The origin of these changes is mainly due to changes in the distribution of sodium and potassium channels. During development, there is evidence that the number of voltage-gated Na⁺ channels increases in rodents (Huguenard et al. 1988, Picken Bahrey and Moody 2003) and in humans (Erraji-Benchekroun 2005), which may accelerate the kinetics of the action potential (up-stroke) and, similar to our results, increase the amplitude of the action potential (Perez-García 2021, Kroon 2019, Zhang 2004, McCormick 1987, Etherington 2011). During aging, however, regulation is disrupted by the change in channel expression, which can lead to altered conduction rate of ionic currents. In brain senescence, the amplitude of the action potential and a slowing of the action potential (Luebke 2007) have been shown to be similar to our data, which may be caused by the decreased expression of voltage-gated Na⁺ channels (Erraji-Benchekroun 2005).

We examined the firing patterns of pyramidal cells at different ages by comparing the F-I slope, and we found that the patterns are fairly conserved across ages. In rodent developmental studies, a decrease in F-I slope was found only in the first two weeks after birth (Zhang 2004, McCormick 1987), which time window is not represented in our data. The latency of the first spike is influenced by the time-constant and specific ion channels, like the transient (or A-type) K⁺ current and the transient Ca²⁺ current (or T-type) (Molineux 2005). The same ion channels are involved in burst firing at the onset of stimulation. Note that in our sample only ∼10 % of pyramidal cells showed burst firing at the beginning of stimulation, but it is known that burst firing is more frequent from layer 3 onwards at depth (Berg 2021). The presence of A-current has been found in immature pyramidal neurons from sensorimotor areas in rodent studies (McCormick 1987), and the density of T-type currents also remains unchanged during the development of the visual cortex (Horibe 2014). It has been shown that the expression of the T-type calcium channel family (CACNG3) decreases during aging (Erraji-Benchekroun 2005). However, it must be mentioned that the differences in the expression of all subtypes of the Kv4 (A-type) and Ca_V3.x (T-type) channel families during brain development and aging have not been systematically described.

Properties of dendritic trees and spine distribution with age

Transcriptomic analyses of the human neocortex show relatively stable expression patterns of genes for dendritic and synaptic development after about 1 year of age (Kang 2011). In early postmortem studies human pyramidal neurons were studied for morphological changes as a function of age. In these studies, using Golgi staining methods, it was found that the dendritic trees of infant layer 2/3 pyramidal cells are well developed (Mrzljak et al. 1990). In the human prefrontal cortex after birth, the dendritic tree of layer 3 pyramidal neurons reaches its structurally mature form after about 3 months (Petanjek 2008). Thereafter, it is assumed that a smaller increase in dendritic length reaches its final form between about 7.5-12-months (Koenderink 1995, Petanjek 2008). In our study, we observed stable patterns of apical and basal dendrite lengths and dendritic tree complexity across the lifespan. Even between the early postnatal groups (infancy vs. early childhood), no significant difference in dendritic tree size was observed. This can be explained by the influence of various factors, such as the large difference between individual subjects or the subtype specific dendritic morphology of pyramidal cells (Berg et al 2021) or different maturation curves of neurons from different brain regions (Elston and Fujita 2014) or cortical layers (Petanjek 2008).

Neurons exhibit a phase of synaptogenesis overproducing synapses lasting months or years in the human cortex (Huttenlocher 1979, Huttenlocher et al. 1982, Mrzljak et al. 1990, Petanjek et al. 2011), followed by dendritic spine/synapse pruning that is reported to last more than a decade in the cortex (Petanjek et al. 2011, Jacobs 1997, Mavroudis 2015, Benavides-Piccione 2013, Coskren 2015). Our results from comparing two groups: infant and late adult pyramidal cells also show that the overall number of dendritic spines decreases with increasing age. When evaluating the distribution of spine shapes in young and old pyramidal cells, we found spines with larger heads (mushroom shape) in a greater number in the dendrite of old pyramidal cells, which are considered as mature synapses. The size of the head is an indicator of the size of the postsynaptic density, the number of glutamate receptors, and the strength of the synapse (Bourne and Harris 2008, Matsuzaki 2001, Nusser 1998). In the developing nervous system, however, the filopodia, thin dendritic protrusions without a head, are the most characteristic type of spines. The density of filopodial spines was higher on the dendrites of infant pyramidal cells. In the older cells, we still found filopodia, which are the silent precursor of active synapses (Vardalaki 2022). The change in the ratio of subtypes from young to old age is well established and is associated with the basis for changes in cognitive function during aging (Luebke 2015, Dumitriu 2010, Jacobs 1997).

The complexity of human brain activity and cognitive abilities increases with development and decreases with aging, resulting in an inverted U-shaped curve across the lifespan (Craik and Bialystok 2006). Cognitive functions depends on many age dependent factors, such as the density and specificity of synaptic connections formed by synaptic pruning and plasticity, or the degree of myelination and white matter maturation and there is also an age-related modulation of neurotransmitters, hormones, ion transporters and receptors (Luna 2015). Here we have shown that pyramidal cells become less excitable and temporarily more precise during development by changing their intrinsic functional properties and with aging these changes occur somewhat in opposite directions making the intrinsic parameters change symmetrically. In addition some of the changes are asymmetrical either occur with development or with aging, such as the resting membrane potential which mostly changes in young ages or the ratio of sag, which is shifted most in old age. These changes in the intrinsic properties of cells during the first and last stages of life also contribute to the input-output functions of a neuron and ultimately to the age-related development of cognitive abilities.

Materials and methods

Slice preparation

Experiments were performed according to the Declaration of Helsinki with the approval of the University of Szeged Ethical Committee and Regional Human Investigation Review Board (ref. 75/2014). Prior to surgery, the patients provided written consent for all tissue material. We used human cortical tissue adjacent to the pathological lesion that had to be surgically removed from patients (n = 64 female n = 45 male) as part of the treatment for tumors, hydrocephalus, apoplexy, cysts, and arteriovenous malformation. 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. The patients received 0.5 mg/kg rocuronium to facilitate endotracheal intubation. The trachea was intubated, and the patient was ventilated with a mixture of O₂-N₂O at a ratio of 1:2. Anesthesia was maintained with sevoflurane at a care volume of 1.2–1.5. During the surgical procedure tissue blocks were removed from parietal (n = 22), temporal (n = 35), frontal (n = 42), and occipital (n = 10) regions, the resected tissue blocks were immediately immersed in ice-cold solution. Slices were cut perpendicular to the pia mater at a thickness of 320 µm with a vibrating blade microtome (Microm HM 650 V) in ice-cold solution (in mM) 75 sucrose, 84 NaCl, 2.5 KCl, 1 NaH₂PO₄, 25 NaHCO₃, 0.5 CaCl₂, 4 MgSO₄, 25 D(+)-glucose, saturated with 95% O₂ and 5% CO₂. The slices were incubated in the same solution for 30 minutes at 36°C following that the solution was changed to (in mM) 130 NaCl, 3.5 KCl, 1 NaH₂PO₄, 24 NaHCO₃, 1 CaCl₂, 3 MgSO₄, 10 D(+)-glucose, saturated with 95% O₂ and 5% CO₂, the slices were kept in it until use.

In vitro electrophysiological recordings

Somatic whole-cell current-clamp recordings were obtained at ∼36 °C in solution containing (in mM) 130 NaCl, 3.5 KCl, 1 NaH₂PO₄, 24 NaHCO₃, 3 CaCl₂, 1.5 MgSO₄, 10 D(+)-glucose, from layer 2/3 pyramidal cells visualized by infrared differential interference contrast (DIC) video microscopy equipped with micromanipulators (Luigs and Neumann, 652 Ratingen, Germany) and HEKA EPC 9&10 patch clamp amplifier (HEKA Elektronik GmbH, Lambrecht, Germany). Micropipettes (3-5 MΩ) were filled with intracellular solution containing (in mM) 126 potassium-gluconate, 4 KCl, 4 ATP-Mg, 0.3 GTP-Na₂, 10 HEPES, 10 phosphocreatine, and 8 biocytin (pH 7.20; 300 mOsm). After whole-cell configuration was obtained stepwise currents were injected to measure the evoked sub- and suprathreshold membrane potential properties. For the analysis of the electrophysiological recordings n = 457 recordings with a series resistance (Rs) of 24.93 ± 11.18 MΩ (max: 63.77 MΩ) were used. For the analysis of fast parameters related to the action potential (AP half-width, AP upstroke velocity, AP amplitude and rheobase), higher quality requirements were set and cells with Rs > 30 MΩ were excluded. This reduced the data set to n = 331 cells with Rs 19.42 ± 6.2 MΩ.

Data analysis

Electrophysiological features were measured from voltage responses elicited by 800 ms long current steps increasing by 20 pA from -100 pA. We analyzed the electrophysiological data with Fitmaster software (HEKA Elektronik GmbH, Lambrecht), and custom MATLAB (The Math Works, Inc.) scripts.

Analysis of morphological features was made by NeuroExplorer software (MBF Bioscience, Williston, VT, USA) and Origin 9 (OriginLab, Northampton, MA).

Statistics

We used custom-written R scripts (R 4.1.2) and the gamm function from the mgcv R package (mgcv 1.8.38). Data presented as the mean ± s.d. Normality was tested with the Lilliefors test, for statistical analysis, ANOVA with posthoc Bonferroni test, Kurskal-Wallis with posthoc Dunn test, for pairwise comparison two-sample t-test or Mann-Whitney test was used. Differences were accepted as significant if p < 0.05. The data are shown on boxplots, boxes indicate 25th, 50th (median), and 75th percentiles, rectangle represents the mean value, and whiskers indicate s.d.

Histology and reconstruction

Slices were fixed in a fixative of 4% paraformaldehyde, 15% picric acid, and 1.25% glutaraldehyde in 0.1 M phosphate (PB, pH = 7.4) for at least 12 hours after electrophysiological recording. After multiple washes in 0.1 M PB, slices were cryoprotected in 10% then 20% sucrose solution in 0.1 M PB. The slices were frozen in liquid nitrogen and then thawed in PB. Slices were embedded in 10% gelatin and further sectioned into 70 µm thick sections. Sections were incubated in a solution containing conjugated avidin-biotin horseradish peroxidase (ABC; 1:100; Vector Labs) in Tris-buffered saline (TBS, pH = 7.4) overnight at 4°C. The enzyme reaction became visible by using 0.05% 3’3-diaminobenzidine tetrahydrochloride as a chromogen and 0.01% H₂O₂ as an oxidant. Sections were post-fixed with 1% OsO₄ in 0.1 M PB. Following several washes with distilled water, sections were stained in 1% uranyl acetate, and dehydrated in an ascending series of ethanol. The sections were infiltrated with epoxy resin (Durcupan, Sigma-Aldrich) overnight and embedded on glass slides. After the electrophysiologically recorded cells have been visualized by DAB staining, 3D light microscopic reconstructions were carried out using the Neurolucida system (MBF Bioscience, Williston, VT, USA) with a 100× objective. The length of terminal segments was measured between the terminal tip of the dendrites and the last branching point before the terminal tip. Dendritic spine density on the dendritic branches was calculated as spine/μm between two bifurcations. The percentage of spines on branches was calculated as the percentage of the given spine type from all the types of spines on that branch.

Supplementary figures and tables

Supplementary table 1.

Supplementary table 2.

Supplementary Figure 1.

Supplementary Figure 2.

Supplementary Figure 3.

Supplementary figure 4.

Supplementary Figure 5.

Supplementary Figure 6.

Supplementary Figure 7.

Supplementary Figure 8.

Supplementary Figure 9.

Supplementary Figure 10.

Supplementary Figure 11.

Supplementary Figure 12.

Supplementary Figure 13.

Supplementary Figure 14.

Supplementary Figure 15.

Supplementary Figure 16.

Supplementary Figure 17.

Supplementary Figure 18.

Supplementary Figure 19.

Supplementary Figure 20.

Supplementary Figure 21.

Acknowledgements

The authors thank Éva Tóth, Katalin Kocsis, Leona Mezei, Otília Kelemen and Bettina Lehóczki for assistance in anatomical experiments and Gergely Komlósi, Szabina Horváth-Furdan, Ildikó Piszár, Sándor Lovas, Balázs Kovács, Rajmund Lákovics, Márton Rózsa, Gáspár Oláh, Attila Ozsvár, Joanna Grace Sandle and Zoltán Péterfi for electrophysiological recordings.

Additional information

Funding

This work was supported by Eötvös Loránd Research Network grants HUN-REN-SZTE Agykérgi Neuronhálózatok Kutatócsoport and KÖ-36/2021 (G.T.)

Ministry of Human Capacities Hungary (20391-3/2018/FEKUSTRAT and NKP 16-3-VIII-3)(G.T.); National Research, Development and Innovation Office grants GINOP 2.3.2-15-2016-00018, Élvonal KKP 133807, ÚNKP-20-5 - SZTE-681, 2019-2.1.7-ERA-NET-2022-00038,

TKP2021-EGA-09, TKP-2021-EGA-28 (G.T.) ÚNKP-21-5-SZTE-580

National Institutes of Health awards U01MH114812 (G.T.) and UM1MH130981 (G.T.)