Developmental expression of neural and synaptic proteins
We selected proteins which, in rodent, are developmentally regulated and/or related to neuronal plasticity. Blot by blot we aimed to compare expression profiles in VC to SC beginning at midgestation E65, just after the begin of gyrification (Ernst et al., 2018), at late gestation E80, E95, E100, and at P90 (adult). Levels at P90 had been set to 1 because we expected proteins to increase towards adulthood.
Caspase-3 expression indicates cell death which is high early in development. As expected, expression was high early on in VC declining after E95. Decline in SC commenced already after E65. Lowest level was seen at P90 (Fig. 1A). Translocase of the outer membrane, TOM70, works in a complex to transport proteins into mitochondria. It is a general metabolic marker useful for assessing biogenesis of mitochondria (Kreimendahl and Rassow, 2020; Ma et al., 2023). Given the dramatic volume increase of the cortex from midgestation to postnatal and the rapid differentiation of energy-hungry neurons as well as the addition of microglia, myelin-forming oligodendrocytes and blood vessels (Ernst et al., 2018; Sobierajski et al., 2022; Sobierajski et al., 2023; Sobierajski et al., 2024) we expected an increase. Indeed, TOM70 expression increased substantially reaching the adult level at E100 in SC, albeit not yet in VC (Fig. 1B). Looking more to neurons, the neuron-specific RNA splice factor NeuN comes in two isoforms of 46 and 48 kDa. The larger isoform is faintly expressed at E65 (Fig. 1C). Both isoforms increased from E80 onwards, and the adult level is reached at E100 in SC, albeit not yet in VC. The E65 levels were higher than levels at E80 and E95 resulting, unexpectedly, in a U-shaped expression profile (Fig. 1C). The microtubule-associated Ring finger E3 ligase protein Trim46 is involved in neuronal polarity and axonal outgrowth (Curcio and Bradke, 2015). Expression was highest at E65 remaining at a constant level thereafter (Fig. 1D). The enrichment of Trim46 in the axon initial segment maintains neuronal polarity and it is in line with the presence of βIV-spectrin-positive axon initial segments in pig cortex as early at E70 (Ernst et al., 2018).
Growing axons secrete transmitters, and since synaptic vesicle proteins become selectively delivered to axons (Watson et al., 2023) their expression could be a proxy for development of axonal connectivity. The vesicle docking and fusion protein Munc-18 has a crucial role for exocytosis. It is at near-adult levels already from E65 onwards in VC and in SC (Fig. 1E). Munc-18 is involved in organizing the presynapse rather than being a marker for synaptogenesis, since it was already highly expressed when the vesicle-associated proteins p38 and p70 were still very low. Yet, transmitter release from nascent axons can occur without molecularly mature vesicles (Andreae and Burrone, 2014). P38/synaptophysin is a marker for presynapse development. It was present weakly in SC at E65 but not yet in VC at E65. Protein amounts increased slowly until E100 and then doubled to P90 (Fig. 1F). Similarly, p70/synapsin-1a/b was detectable at E65 with higher levels in SC than in VC and increasing in both areas (Fig. 1G). Another presynaptic active zone protein is RIM1a. The expression profile was U-shaped in VC and in SC. Levels at E100 were not yet near P90 levels in both areas (Fig. 1H).
We also compared P90 pig VC to rat VC (Supplemental Fig. 1A-D) by using blot images not included in Figs. 1–4. The aim was to confirm the quantification with independent blots and to explicitly focus on the delayed development of the VC. The comparison revealed that NeuN, Trim46, Munc-18, p70/synapsin-1 and RIM1 expression were fairly similar to rat in terms of postnatal band intensity and molecular weight. Pig caspase-3 ran lower than rat caspase-3 and the smaller isoform seen in pig up to E100 has not been reported in rodent. TOM70 had a doublet band in rat lysates but never in pig lysates. Rat p38/synaptophysin ran slightly lower than pig p38 (Supplemental Fig. 1A).
A. Cell-death-related caspase-3. B. Mitochondrial marker translocator of outer membrane, TOM70. C. Neuron-specific RNA splice factors, NeuN. D. Neuronal polarity and axon initial segment marker Trim46. E. Presynaptic organizer Munc-18. F. Presynaptic vesicle protein p38/synaptophysin. G. Presynaptic vesicle protein p70/synapsin-1a/b. H. Presynaptic organizer RIM1. β-actin at 42 kDa was used for normalization.
In Fig. 1, 2, 3, 4, VC is left, SC is right, each with 5 stages from embryonic day (E) 65, E80, E95, E100 to postnatal day P90 the average of which has been set to 1. Numbers in the bars indicate the number of independent lysates which delivered measurable bands. n.m., not measurable; the protein of interest was not detectable in these lysates. In case a lysate was run twice for a given marker, the average has been included to avoid technical repeats. Bars represent the mean with S.E.M.. Representative gel stripes are arranged together with the housekeeping protein from the same lanes. Molecular weight [kDa] is given.
Before continuing describing the other proteins, we shall briefly address the expression profiles. Seeing more intense bands with age, as with TOM70, p38/synaptophysin, p70/synapsin-1a/b, is commonly interpreted as an increase of protein amounts. This has been also the case for the expression of myelin-related proteins and GFAP detected with blots from the very same material (Sobierajski et al., 2023). Seeing less intense bands, as with caspase-3 or Trim46, is commonly interpreted as a decrease. This has been also the case PDGFRα, a marker for immature oligodendrocytes (Sobierajski et al., 2023). In contrast, U-shaped profiles as revealed for NeuN, RIM1 or PSD95 (Fig. 3) look peculiar. It might suggest a transient expression, or transient interaction with binding partners masking the antigenic epitope, up to harvest at a wrong circadian time point or simply postmortem protein degradation. Myelin proteins are particularly susceptible to the latter, however, within our postmortem time periods (Ernst et al., 2018) we could not detect any degradation (Sobierajski et al., 2023). Moreover, such interpretations make little sense in case of NeuN and PSD95 which are continuously increasing in neurons with ongoing differentiation, and both are not regulated by circadian mechanisms. Rather, we have to consider the dramatic volume increase of the brain from ~ 5 g at E45 to ~ 90 g at P90 (Ernst et al., 2018; Sobierajski et al., 2023). At E65, the future cortical gray matter consists mostly of densely packed immature neurons and the still migrating neurons of the cortical plate, whereas microglia (Sobierajski et al., 2022) and immature oligodendrocytes are still scarce, and astrocytic GFAP protein expression is barely detectable (Sobierajski et al., 2023). Therefore, at E65 neuronal proteins expressed by zillions of neurons will be the dominant proteins in a lysate of 30 µg total protein. With time, these early expressed neuronal proteins become diluted because proliferating, invading and/or maturating microglia, macroglia and endothelial cells contribute more and more non-neuronal proteins to the lysate. Accordingly, the band intensity drops transiently only to recover towards the postnatal stage with ongoing neuronal differentiation.
We therefore suggest an alternative interpretation. U-shaped profiles are seen with proteins which are expressed at substantial amounts early on and which then remain rather on a plateau. A constant protein level over the ages should be rather regarded as a continuous increase of expression with the “increase” being compensated for by the “dilution”, be it via neurons or glia which for instance also express many transmitter receptors. An increasing intensity of a band should be regarded as an even steeper increase of protein expression, again, in neurons or in glia cells. The caspase-3 expression then would initially increase before declining. Reports on U-shaped expression profiles are rare presumably because midgestation brain samples are rarely included. For instance, expression of protein tyrosine phosphatase 1D in fetal superior colliculus displays a U-shaped profile which has been interpreted as reflecting a switch in function (Reinhard et al., 2009). Increasing, decreasing and U-shaped profiles have also been shown for claudin-3, Jam-A and claudin-12, respectively, in fetal sheep from 0.60 gestational age onwards (Sadowska et al., 2015).
Developmental expression of plasticity-related proteins
PSD95 anchors NMDA receptors. The mRNA encoding the obligatory subunit GluN1 increases in rat cortex from birth to peak levels during the fourth postnatal week followed by a decline towards adulthood (Hofer et al., 1994) and the protein increases about 3-fold from birth to P40 (Luo et al., 1996). In pig, however, the strongest expression was at E65 in VC and less so also in SC with levels being 3-fold higher than at P90 (Fig. 2A) confirming the view that this subunit is not limiting the assembly of NMDA receptors. At the following fetal ages, the GluN1 band became extremely faint and was not always measurable by densitometry in all lysates tested. The band was again well detectable at P90. With all caution towards measuring such a faint expression, the bar graphs did suggest a small increase from E80 to E95. In contrast, the synaptic plasticity-related subunit GluN2B was well detectable from E65 onwards in VC and SC. A steep increase to the adult levels has occurred at E100 in SC, albeit not yet in VC (Fig. 2B). Of note, the expression of glutamate and GABA receptors may not solely mirror neuronal maturation since AMPA, NMDA and GABAAR are also expressed in macroglia (Gundersen et al., 2015, for review).
The 55 kDa CamKIIα subunit is expressed in pyramidal cells and upregulates in rodent VC around eye opening at the end of the second postnatal week. CamKIIα autophosphorylation is essential for NMDAR-dependent plasticity (Gordon et al., 1996; Glazewski et al., 2000). Highly interesting, CamKIIα protein was expressed and steadily increasing in pig VC and SC from E65 onwards (Fig. 2C) suggesting that CamKIIα mirrored the functional differentiation of excitatory neurons. Also important, the T286-phosphorylated form was detectable at E100, a bit stronger already in SC than in VC. Substantially higher CamKIIα levels was seen at P90 (Fig. 2D) suggesting a massive increase after E100. The CamKIIα interaction with α-actinin-2 is required for structural long-term potentiation (Curtis et al., 2023) and spine growth and stability. Indeed, α-actinin-2 was well expressed in pig VC and SC and steadily increasing from E65 onwards (Fig. 2E). Also essential for spine maturation and maintenance is the spine apparatus with synaptopodin via its interaction with α-actinin-2. Synaptopodin comes as a long and a short transcript, and PCR detects both in brain and in kidney (Asanuma et al., 2005). At the protein level, only the smaller, ~ 100 kDa isoform is expressed in brain. In pig, expression was U-shaped with synaptopodin present at E65, not measurable at E85 and E95 and again detectable at near-adult levels at E100 in SC and less so in VC (Fig. 2F). Comparing pig VC to rat VC revealed the bands of GluN1, GluN2B, CamKIIα, α-actinin-2 and synaptopodin at the same molecular weight as in rat VC. The intensity of the P90 GluN1 and synaptopodin bands was lower in pig VC than in rat VC, which might reflect a difference in antibody binding (Supplemental Fig. 1B).
A. NMDAR subunit GluN1. B. NMDAR subunit GluN2B. C. Postsynaptic enzyme CamKIIα, total protein. D. T286 phosphorylated CamKIIα. E. Structural plasticity-related actin bundling protein α-actinin-2. F. Spine apparatus and cisternal organelle-associated synaptopodin. β-actin at 42 kDa was used for normalization. n.m., the protein was not measurable at this age in this cortex.
Development of pyramidal cell dendritic spines and spine types
The density and morphology of dendritic spines are important indicators of the development of neuronal maturation. Here, we acquired confocal scans of DiI-stained neurons from the VC and SC from E70 to P90 and reconstructed the basal dendrites and classified the dendritic protrusions (addressed as ‘spines’ below). Only second and higher order dendritic segments were included. Dendritic end segments were excluded as the thin distal tips of dendrites often have fewer spines. Thick proximal dendritic shafts were excluded because short spines may be masked by the thick shaft (Xu et al., 2020). At E70, neurons had small somata (Fig. 5A). They form spines with immature morphologies which made classification difficult, and we did not determine spine types at this stage (Fig. 5A1). At E85, neurons had a much more mature morphology (Fig. 5B). Now, spine types could be determined (Fig. 5B1). A substantial increase in density proceeded between E85 and E110, and pyramidal cells developed larger somata with longer and more branched dendrites (Fig. 5C). Pyramidal neurons of the P5 domestic piglet did not obviously differ from E110 boar cells in morphology (Fig. 5D) but had a lower spine density (Fig. 5D1). Interestingly, volume, branching, and segment length of blood vessels of the P5 domestic piglet were somewhat lower than values obtained in E100/E110 pig marginal zone and gray matter suggesting that also vascular development tends to proceed slower in domestic pig than in wild boar (Sobierajski et al., 2024). Until P90, the complexity of pyramidal neurons (Fig. 5E) and the density of spines had increased further (Fig. 5E1).
The data base of the spine analysis is given in Fig. 6A. From E70 to E85, spine density moderately increased on average from 8 to 18 spines per 100 µm (Fig. 6B) followed by a roughly threefold increase until E110 to a median of 60 spines per 100 µm. As expected from the qualitative data, the P5 domestic piglet had a spine density of on average 43 spines per 100 µm which was in the lower range of what is present in boar at E110. Until P90, spine density increased to 79 spines per 100 µm (Fig. 6B). At E85, the majority of spines were classified as thin and filopodial, mushroom and stubby spines were rare (Fig. 6C). Thin spines were still the dominant fraction at E110, but the density of mushroom and stubby spines had substantially increased whereas filopodia became rare (Fig. 6D). This trend continued at P90 with stubby spines remaining at the level seen at E110, while thin and mushroom spines continued to increase (Fig. 6E).
When looking at Fig. 6B, the almost three-fold variability between individual cells at E110 and at P90 might indicate large cell-individual differences in spine densities. This is known to occur in vivo but could potentially suggest a sampling bias. The latter can be excluded. The variability decreased substantially when separating the cells by area. At E110 and at P90, SC neurons had significantly higher spine densities than VC neurons. Moreover, spine density increased significantly between E110 and P90 in VC neurons, whereas the spine density in SC neurons at E110 was statistically not different from the density at P90 (Fig. 6F). The developmental advance of SC pyramidal neurons and the much slower development of VC pyramidal cells, respectively, was mirrored by the proportion of mushroom spines. In SC, the proportion seen at P90 was reached already at E110, whereas in VC the proportion of mushroom spines was still low at E110, increasing to P90 albeit not yet to the level reached in P90 SC (Fig. 6G).