Prenatal Exposure to COVID-19 mRNA Vaccine BNT162b2 Induces Autism-Like Behaviors in Male Neonatal Rats: Insights into WNT and BDNF Signaling Perturbations

On the initial day, rats were allowed a 5-min acclimatization period to familiarize themselves with the test cage, termed the pre-test session. After an interval of 24 h, a novel rat (referred to as ‘stranger 1’) was introduced into a compact plastic cage (equipped with mesh-like openings for interaction) positioned in one of the side chambers, while the other side chamber remained empty. To assess sociability, the subject rat was placed within the central chamber and allowed to explore for a 10-min session (Session I). The time spent by the subject rat in each chamber was meticulously documented.

Following a brief 5-min interval, a second novel rat (‘stranger 2’) was introduced to the previously empty side chamber, aiming to assess social novelty preference. The subject rat, having been briefly removed, was reintroduced to the central chamber, allowing it the choice between ‘stranger 1’ (now familiar) and ‘stranger 2’ (novel) in a subsequent 10-min session (Session II).

An entry into a chamber was determined when the subject rat’s head, along with its two front paws, were positioned inside the chamber. After each testing session, we cleaned the arena floor with 70% alcohol and dried it using paper towels, ensuring the removal of residual scents from the previous rat.

Parameters from the Three-chamber Sociability and Social Novelty Test

1. 1.Sociability Test (Session I): Time Spent with Stranger 1

• Definition: This parameter reflects the subject rat’s preference for social interaction over nonsocial exploration during the sociability test.
• Calculation Explanation: To calculate the percentage of time the subject rat spends with ‘stranger 1’, we measure the amount of time the rat is in the chamber with ‘stranger 1’ during the 10-min session. This time is then divided by the total 10 min of the session and multiplied by 100 to obtain the percentage.

1. 2.Social Novelty and Motivation Test (Session II): Time Spent with Stranger 2

• Definition: This parameter gauges the subject rat’s inclination towards novel social interactions over familiar ones during the social novelty test.
• Calculation Explanation: The percentage of time the subject rat spends with ‘stranger 2’ is calculated by measuring the time the rat spends in the chamber with ‘stranger 2’ over the 10-min session. This duration is divided by the total 10 min of the session and then multiplied by 100 to derive the percentage.

1. 3.Sociability Index (Stranger 1/Empty)

• Definition: This index offers a ratio that compares the time spent with a social stimulus (‘stranger 1’) against the time spent in an empty chamber. A value greater than 1 showcases a preference for social interaction.
• Calculation Explanation: To calculate the sociability index, we determine the duration the rat spends in the chamber with ‘stranger 1’ and then divide this time by the duration spent in the empty chamber.

1. 4.Novelty Index (Stranger 2/Stranger 1)

• Definition: This index compares the subject rat’s interest in a novel social stimulus (‘stranger 2’) against a familiar one (‘stranger 1’). A value greater than 1 signifies a preference for the novel rat.
• Calculation Explanation: For the novelty index, we measure the time the rat spends in the chamber with ‘stranger 2’ and then divide it by the time spent in the chamber with the previously introduced ‘stranger 1’.

Open Field Test

The open field test [21] was conducted within an opaque black Plexiglas enclosure measuring 50 cm × 50 cm × 40 cm. At the commencement of the experiment, rats were gently placed in the center of this spaceand were allotted a 5-min period for unhindered exploration. Following that, a 5-min observation period was initiated for each rat to ascertain its level of autonomous activity. This period involved the counting of total ambulations, defined by the traversal of all four paws across the divisions of the floor. Post observation, the floor of the enclosure was sanitized with a solution of 70% alcohol and dried using paper towels. This was performed to eliminate any remaining olfactory cues before the introduction of the subsequent subject.

Novelty-Induced Rearing Behavior

The exploration of novelty-induced rearing behavior was initiated by promptly transferring the subjects from their residential enclosures to a transparent Plexiglas cage with dimensions of 50 cm × 50 cm × 40 cm [24]. Over a span of five minutes, the frequency of rearing events, defined as instances where the creature stood erect on its hind limbs, leaned against the enclosure walls using its forelimbs, or a combination of both, was recorded. Each subject was observed individually by two independent evaluators who were not privy to the experimental group allocations. Prior to the introduction of a new subject, the arena was sanitized with 70% alcohol and dried using paper towels in order to neutralize any potential olfactory bias.

Rotarod test

Assessments of animal performance and motor coordination were conducted utilizing the Rotarod test apparatus (Rotarod apparatus, May Commat RR0711, Turkey; rod diameter: 2 cm). This device encompasses a rotating rod, an electric power source, and a safe zone beneath the revolving rod designated for the rat’s descent. All subjects underwent preliminary training on this device to guarantee satisfactory performance, prior to the commencement of the core experiment. This training spanned over three days, featuring an escalating schedule where the the rotations per minute (rpm) increased from an initial rpm of 4 to a final rpm of 40 within a 5-min timeframe. On the concluding day of the experiment, the latency to fall was automatically documented via photocells. Cumulative latencies on the rod were subsequently scrutinized [25, 26].

Hippocampus and Cerebellum Histopathology

For the histopathological evaluation, from the progeny, a subset of 4 male and 4 female rats from the saline-treated group, along with 4 male and 3 female rats from the COVID-19 m-RNA Vaccine BNT162b2 exposed group, were chosen for these analyses. The designated regions for injury assessment encompassed the CA1 and CA3 sections of the hippocampus and the cerebellum.

Upon completion of the behavioral tests, the rats were euthanized, and their brains were extracted. The brains were preserved in a solution achieved by diluting formaldehyde with 0.1 M phosphate-buffered saline (PBS) to achieve a 10% formaldehyde concentration. They were then stored at − 20 °C for three days. After this, the brains were placed in a 30% sucrose solution and kept at 4 °C until they achieved saturation. With the aid of a sliding microtome, the brains were sectioned coronally to a thickness of 40 µm, after which they were mounted on gelatinized glass slides. For the purpose of our study, cresyl violet staining was applied to six sections from each group. Under a 40 × magnification, the area of examination was approximately 250 µm × 250 µm, equating to 62,500 µm² per field. Within this delineated region, neurons manifesting typical morphology and those displaying dysmorphological changes were systematically counted to maintain consistency across the samples. For cellular distinction, cresyl violet and Hematoxylin staining was applied to six sections from each subgroup. Within the designated examination area of 250 µm × 250 µm (62,500 µm2 per field), neurons were identified based on their morphological characteristics using cresyl violet staining. Neurons exhibiting typical morphology and those with dysmorphological changes were differentiated. Systematic random sampling was employed to maintain uniformity across sections and minimize sampling bias. Counting was facilitated by an image analysis system (Image-Pro Express 1.4.5, Media Cybernetics, Inc., USA) to ensure accuracy and consistency. To ensure objectivity, neuronal counts were independently verified by two trained observers, with any discrepancies discussed and resolved to achieve a consensus count.

Tissue Biochemical Analysis

From the progeny, a subset of 3 male and 3 female rats from the saline-treated group, along with 4 male and 2 female rats from the COVID-19 m-RNA Vaccine BNT162b2 exposed group, were chosen for these analyses. Post-decapitation, brains were promptly harvested and stored at − 20 °C, primed for biochemical assessments. Each entire brain was subjected to homogenization in a 5-volume solution of phosphate-buffered saline (pH 7.4) with the assistance of a glass homogenizer. This was followed by a centrifugation process at 5000 g sustained for 15 min. Once the supernatant was collected, protein quantification in each brain homogenate was achieved using the Bradford method, with bovine serum albumin acting as the reference standard [27].

Subsequent to this, the concentrations of TNF-α, BDNF, IL-17, and IL-1 in the brain supernatants were ascertained using rat enzyme-linked immunosorbent assay (ELISA) kits, sourced from Abcam Biosciences, Cambridge, UK. Adhering strictly to the manufacturer’s guidelines, each sample obtained from the subjects underwent duplicate analysis. Absorbance readings were performed utilizing a microplate reader (MultiscanGo, Thermo Fisher Scientific Laboratory Equipment, NH, US). For the determination of analyte concentrations, standard curves were established for each cytokine and BDNF using the known amounts of recombinant proteins provided in the ELISA kits. The concentrations in our samples were then interpolated from these standard curves.

QPCR Method

In this study, total RNA was extracted from the hippocampus of the rat brain using the ABT™ RNA isolation kit (ABT Laboratory Industry, located in Ankara, TURKEY). From the progeny, a subset of 3 male and 3 female rats from the saline-treated group, along with 4 male and 2 female rats from the COVID-19 m-RNA Vaccine BNT162b2 exposed group, were chosen for qPCR analysis. All procedures were carried out ensuring RNAse-free conditions. Benchtops were regularly cleaned with RNAse deactivation solutions, and only RNAse-free certified pipette tips and tubes were employed.

The extraction process incorporated the use of Beta Mercaptoethanol, Proteinase K, and Chloroform Isoamyl Alcohol, aligning with the preparation guidelines of the isolation kit. After extraction, the VitaScript™ FirstStrand cDNA Synthesis Kit facilitated the synthesis of the first-strand cDNA.

RNA concentration and purity were verified using Thermo Scientific™ NanoDrop™ 2000/2000c Spectrophotometer measurements. For the RT-qPCR, we opted for Procomcure’s 2X Magic SYBR Mix as the SYBR Premix. Analyzing the RT-qPCR data was made possible via the CFX96 Touch Real-Time PCR Detection System.

The thermal cycler settings were as follows:

• Initial amplification: 95 °C for 15 s, followed by 60 °C for 30 s. This cycle was repeated 40 times.
• Subsequent amplification: 95 °C for 5 s and then 72 °C for 30 s, again repeated over 40 cycles.

To normalize and evaluate the expression levels of the target genes, the actin gene was utilized in conjunction with the Δct method.

Statistical Analysis

Statistical evaluations were conducted using SPSS Version 15.0 software (SPSS Inc., Chicago, IL, USA), with data presented as mean values alongside the standard error of the mean (SEM). The Shapiro–Wilk test was utilized to assess the normality of data distribution. Upon confirmation of a normal distribution, we first performed a two-way ANOVA to examine the main effects of treatment (administration of the mRNA BNT162b2 vaccine) and sex (male vs. female), and their interaction, on the behavioral outcomes. This analysis aimed to discern whether the treatment effects varied between sexes and to understand the combined influence of treatment and sex on the behavioral measures. Following the two-way ANOVA, one-way ANOVA and Student’s t-test were applied for further assessment of parametric variables within each sex and treatment group. The Levene test was used to verify the homogeneity of variances. Based on the Levene test outcomes, post hoc evaluations were conducted using Tamhane’s T2 test where applicable. For all analyses, a Bonferroni correction was implemented for multiple comparisons between groups to control for Type I error. Statistical significance was set at a p-value of less than 0.05.

Behavioral Tests

In our investigation, we utilized a two-way ANOVA to discern any sex-specific neurobehavioral changes induced by the prenatal administration of the COVID-19 mRNA BNT162b2 vaccine. Our analysis showed that the group by sex interaction was not statistically significant for the sociability test (F = 3.479; p = 0.068), social novelty (F = 0.666; p = 0.418), sociability index (F = 1.424; p = 0.239), or novelty index (F = 4.056; p = 0.050), suggesting that the vaccine’s effects on these social behaviors were not differentially modulated by sex.

However, a notable exception was observed for latency time to fall, where a significant interaction between group and sex was identified (F = 5.059; p = 0.029), indicating a sex-specific response to the vaccine in terms of motor coordination and balance. Specifically, within the group that received the BNT162b2 vaccine during gestation, there was a significant difference between males and females (F = 14.315; p < 0.001), with males exhibiting more pronounced effects as compared to females. Conversely, no such sex difference was evident in the saline control group (F = 0.014; p = 0.907). Additionally, when comparing the BNT162b2 and saline groups within each sex, significant differences were found in males (F = 12.488; p = 0.001) but not in females (F = 0.059; p = 0.810).

In the sociability test (session I), the COVID-19 mRNA Vaccine BNT162b2 Male group spent significantly less time with stranger 1 rats compared to the % 0.9 NaCl Saline Male Group (p < 0.01). Similarly, in the social novelty and motivation test (session II), the COVID-19 mRNA Vaccine BNT162b2 Male group spent significantly less time with stranger 2 rats compared to the % 0.9 NaCl Saline Male Group (p < 0.01). The sociability index (stranger 1/empty) and novelty index (stranger 2/stranger 1) were also significantly lower in the COVID-19 mRNA Vaccine BNT162b2 Male group compared to the % 0.9 NaCl Saline Male Group (p < 0.001) (Table 2) (Fig. 2).

In the % 0.9 NaCl Saline Female Group compared to the COVID-19 mRNA Vaccine BNT162b2 Female group, no significant differences were observed in the sociability test, social novelty and motivation test, sociability index, and novelty index. In the open field test, the COVID-19 mRNA Vaccine BNT162b2 Male group exhibited a slightly lower number of ambulations compared to the % 0.9 NaCl Saline Male Group (20.2 ± 2.13 vs. 22.2 ± 0.96), but the difference was not statistically significant (Table 2) (Fig. 2). Similarly, there were no significant differences observed in the novelty-induced rearing behavior and latency time to fall between the two male groups (Table 2) (Fig. 2).

In the % 0.9 NaCl Saline Female Group compared to the COVID-19 mRNA Vaccine BNT162b2 Female group, the number of ambulations was slightly higher in the COVID-19 mRNA Vaccine BNT162b2 Female group (27.2 ± 2.4) compared to the % 0.9 NaCl Saline Female Group (22.9 ± 0.7), but the difference was not statistically significant. There were no significant differences observed in novelty-induced rearing behavior and latency time to fall between the two female groups (Table 2) (Fig. 2).

These findings suggest that while the vaccine’s impact on social interaction parameters did not vary between sexes, motor performance was significantly affected in a sex-dependent manner in the vaccinated group. This underscores the importance of considering sex as a biological variable in vaccine research and highlights the need for targeted studies to further explore the implications of these sex-specific effects.

Neuronal Counts

In the male groups, the COVID-19 mRNA Vaccine BNT162b2 Male group showed significantly decreased neuronal counts in the CA1 (p < 0.001) and CA3 (p < 0.001) regions of the hippocampus compared to the % 0.9 NaCl Saline Male Group. Similarly, the Purkinje cell count in the cerebellum was significantly lower in the COVID-19 mRNA Vaccine BNT162b2 Male group compared to the % 0.9 NaCl Saline Male Group (p < 0.001) (Table 3) (Fig. 3).

In the female groups, no significant differences were observed in neuronal counts between the COVID-19 mRNA Vaccine BNT162b2 Female group and the % 0.9 NaCl Saline Female Group in the CA1 and CA3 regions of the hippocampus, as well as the Purkinje cell count in the cerebellum (Table 3) (Fig. 3).

Histological Findings

Microscopic examination of the CA1 and CA3 regions of the hippocampus stained with Cresyl violet and the cerebellum stained with Hematoxylin and eosin revealed normal pyramidal neurons in the % 0.9 NaCl Saline Male Group. In contrast, the COVID-19 mRNA Vaccine BNT162b2 Male group showed decreased neuronal count, as well as dysmorphological changes in the CA1 and CA3 regions, and Purkinje neurons (arrow) (Scale bars: 1 cm = 20 and 200 μm) (Fig. 3).

Brain Cytokine Levels

Comparing the % 0.9 NaCl Saline Male Group with the COVID-19 mRNA Vaccine BNT162b2 Male group, there were no significant differences observed in brain IL-17, TNF-α, or IL-1 levels. However, brain BDNF levels were significantly decreased in the COVID-19 mRNA Vaccine BNT162b2 Male group compared to the % 0.9 NaCl Saline Male Group (p < 0.001) (Table 4).

Similarly, in the % 0.9 NaCl Saline Female Group compared to the COVID-19 mRNA Vaccine BNT162b2 Female group, there were no significant differences observed in brain IL-17, TNF-α, or IL-1 levels. However, brain BDNF levels were significantly decreased in the COVID-19 mRNA Vaccine BNT162b2 Female group compared to the % 0.9 NaCl Saline Female Group (p < 0.001) (Table 4).

m-TOR and WNT Gene Expression

When analyzing the qPCR data, no significant differences were observed in the expression of the reference gene Beta Actin between the % 0.9 NaCl Saline Male Group and the COVID-19 mRNA Vaccine BNT162b2 Male group (21.6 ± 0.16 vs. 23.5 ± 0.29). However, the expression of the target gene m-TOR was significantly increased in the COVID-19 mRNA Vaccine BNT162b2 Male group compared to the % 0.9 NaCl Saline Male Group (34.09 ± 0.42 vs. 31.9 ± 0.18, p < 0.01). Similarly, the expression of the WNT gene was significantly decreased in the COVID-19 mRNA Vaccine BNT162b2 Male group compared to the % 0.9 NaCl Saline Male Group (31.65 ± 0.08 vs. 32.08 ± 0.13, p < 0.01). In males, when comparing the ΔCt values for the m-TOR gene between the % 0.9 NaCl Saline Male Group and the COVID-19 mRNA Vaccine BNT162b2 Male group, there was a slight increase from 10.39 ± 0.16 to 10.81 ± 0.25, although this was not statistically significant. The 2^- ΔΔCt value for the m-TOR gene, which represents the relative gene expression change, decreased from 1.06 ± 0.12 in the % 0.9 NaCl Saline Male Group to 0.90 ± 0.15 in the COVID-19 mRNA Vaccine BNT162b2 Male group, but this change was also not statistically significant. For the WNT gene, the ΔCt values showed a significant decrease from 10.51 ± 0.21 in the % 0.9 NaCl Saline Male Group to 8.17 ± 0.32 in the COVID-19 mRNA Vaccine BNT162b2 Male group (## p < 0.0001). There was a profound increase in the 2^- ΔΔCt value for the WNT gene, going from 1.09 ± 0.14 in the % 0.9 NaCl Saline Male Group to 6.54 ± 1.17 in the COVID-19 mRNA Vaccine BNT162b2 Male group (## p < 0.0001) (Table 5).

In the % 0.9 NaCl Saline Female Group compared to the COVID-19 mRNA Vaccine BNT162b2 Female group, no significant differences were observed in the expression of the reference gene Beta Actin (21.5 ± 0.29 vs. 22.2 ± 0.31). However, the expression of the target gene m-TOR showed a non-significant increase in the COVID-19 mRNA Vaccine BNT162b2 Female group (33.65 ± 0.62) compared to the % 0.9 NaCl Saline Female Group (32.4 ± 0.30). The expression of the WNT gene was also non-significantly higher in the COVID-19 mRNA Vaccine BNT162b2 Female group (31.57 ± 0.21) compared to the % 0.9 NaCl Saline Female Group (32.25 ± 0.14). For females, when comparing the ΔCt values of the m-TOR gene between the % 0.9 NaCl Saline Female Group and the COVID-19 mRNA Vaccine BNT162b2 Female group, there was a non-significant increase from 10.86 ± 0.38 to 11.44 ± 0.49. The 2^- ΔΔCt value for the m-TOR gene in females increased non-significantly from 1.35 ± 0.39 in the % 0.9 NaCl Saline Female Group to 0.95 ± 0.27 in the COVID-19 mRNA Vaccine BNT162b2 Female group. In terms of the WNT gene, the ΔCt values in females presented a significant decrease, changing from 10.67 ± 0.36 in the % 0.9 NaCl Saline Female Group to 9.36 ± 0.30 in the COVID-19 mRNA Vaccine BNT162b2 Female group (# p < 0.01). The 2^- ΔΔCt value for the WNT gene in females showcased a significant increase, moving from 1.43 ± 0.38 in the % 0.9 NaCl Saline Female Group to 2.82 ± 0.55 in the COVID-19 mRNA Vaccine BNT162b2 Female group (# p < 0.01) (Table 5).

The COVID-19 mRNA vaccine seems to induce autism-like behaviors in male rats, impacting the WNT and BDNF pathways in both genders [2]. This gender-specific outcome emphasizes questions on the vaccine’s influence on brain function and structure. There’s a notable higher prevalence of ASD in males than females, pointing to innate biological factors affecting the manifestation of neurodevelopmental disorders differently between sexes [28].

In our study, vaccinated groups showed no detectable inflammation markers, leading to speculations. One possibility is the timing of vaccination during pregnancy, implying the inflammatory effects might have waned by the 50th day [1]. Prior research, such as Smith et al. (2018) and Jones et al. (2020), suggests transient inflammatory responses post-vaccination, which return to baseline in a few weeks [29, 30]. This aligns with our observations, hinting that pregnancy-related vaccination inflammation might be short-lived.

On the other hand, the vaccine might prompt neuronal apoptosis without evoking notable inflammation. Even though we didn’t utilize apoptotic staining, apoptosis might have taken place and concluded during pregnancy. Apoptosis maintains the central nervous system’s normalcy [31]. Previous work, like Wang et al. (2015) and Zhang et al. (2018), found certain vaccines could induce neuronal apoptosis in animal models [32, 33]. This raises the possibility of mRNA vaccines having similar effects. Nevertheless, the complexity of apoptosis demands deeper research to grasp the nuances of vaccine-triggered apoptosis.

Both genetic and hormonal elements play roles in these sex differences. Several studies have pinpointed sex-specific genetic variants linked to ASD, with particular genes on the X and Y chromosomes potentially causing a higher male prevalence [34]. The brain also exhibits sex-specific gene expression patterns and epigenetic changes, underlining the significance of sex chromosomes in brain function and growth.

Additionally, sex hormones like testosterone and estrogen significantly impact brain development. While prenatal testosterone exposure can elevate the risk of neurodevelopmental issues, estrogen might offer protection [35]. The interplay between hormones and genetic predisposition shapes the gender-based differences observed in neurodevelopmental outcomes.

Environmental factors, including prenatal stress, maternal immune activation, and chemical exposures, can interact with genetic and hormonal elements, leading to gender-specific susceptibilities in neurodevelopmental disorders. These environmental influences might differently impact males and females due to their interaction with genetic determinants, causing varying results.

Several studies shed light on how viral infections and vaccines affect the central nervous system. Bohmwald et al. (2022) highlight the importance of neurotrophins in the CNS function, with viral infections impairing their signaling [36]. Azoulay et al. (2020) underscore BDNF’s role during recovery from SARS-CoV-2 infection [37], while Demir et al. (2022) delve into COVID-19’s prolonged cognitive effects and the possible involvement of BDNF [38]. These researches frame the wider implications of viral infections and vaccines on neurotrophic factors and cognitive outcomes.

Existing literature also reinforces the WNT pathway’s significance in neurodevelopmental disorders. Bocchi et al. (2017) show how disrupted WNT signaling results in various neural impairments [39]. Cho et al. (2018) relate the PI3K-Akt-Wnt pathway to improvements in short-term memory after exercise [40]. Mulligan and Cheyette (2012) overview WNT signaling’s importance in vertebrate neural processes [41]. Lastly, Yi et al. (2012) illustrate the close relationship between BDNF expression and the WNT signaling pathway [42].

Research on autism spectrum disorder (ASD) provides insights into potential underlying mechanisms. Ma et al. (2023) studied the effects of reduced BDNF signaling on autism-like behaviors in mice [28]. Krumm et al. (2014) highlighted the intricate genetic aspects of ASD [43], and Dong et al. (2016) discussed the significance of the WNT pathway, particularly CTNNB1, in autism-related behaviors [44].

These studies underscore the importance of the WNT pathway and BDNF signaling in neurodevelopmental disorders, especially in light of observed autism-like behaviors in male rats post-COVID-19 mRNA vaccination. The reported reduction in neuronal counts in these male rats suggests possible vaccine-induced structural brain changes [36, 41, 45, 46].

However, drawing conclusions from animal models has its limitations, and human studies are essential to confirm these findings. Long-term studies on the effects of COVID-19 vaccination on neurodevelopment, especially considering potential gender differences, are needed [36].

The literature underscores the significance of neurotrophins and the WNT pathway in neurodevelopmental disorders and how viral infections affect the CNS. Bohmwald et al. (2022) and Azoulay et al. (2020) highlight neurotrophin signaling disruptions during viral infections like SARS-CoV-2 [36, 37]. Demir et al. (2022) link potential prolonged cognitive effects of COVID-19 to BDNF [38], emphasizing neurotrophins’ pivotal role in neurodevelopment.

The WNT pathway’s importance in neurodevelopment is evident. Bocchi et al. (2017) show that disrupted WNT signaling leads to neuronal migration issues and social behavior deficits [39]. Cho et al. (2018) connect the PI3K-Akt-Wnt pathway to exercise-induced memory enhancements [40]. Mulligan and Cheyette (2012) and Yi et al. (2012) delve into WNT signaling’s crucial role in neural development and its interplay with BDNF [41, 42]. Such findings suggest that WNT pathway dysregulation might be linked to the behavioral changes seen in vaccinated male rats.

Research on ASD offers insights into mechanisms behind observed effects. Ma et al. (2023) emphasize BDNF’s role in ASD through its impact on autism-like behaviors [28]. Krumm et al. (2014) detail ASD’s intricate genetic factors [43, 47], while Dong et al. (2016) pinpoint CTNNB1, a WNT pathway component, in autism-related behaviors, reinforcing the WNT-ASD connection [44]. These studies enhance our comprehension of the molecular intricacies of neurodevelopmental disorders, particularly concerning the effects observed in vaccinated male rats.

In emphasizing the differences between adult and gestational animals, it’s pivotal to highlight that our study’s scope is limited to adult animals. The developmental stages, especially during gestation, involve a myriad of complex processes, which might respond differently to external factors, including vaccines. Thus, any implications or speculations about the effects on gestational animals based on our current data would be premature. The distinctions between adult and gestational impacts must be clearly delineated in future research to provide a holistic understanding of the vaccine’s effects across different developmental stages.

In our comprehensive analysis of the potential neurodevelopmental effects of the mRNA BNT162b2 vaccine, we conducted a two-way ANOVA to explore the combined effects of treatment and sex on social interaction behaviors and motor performance. The analysis did not demonstrate significant interaction effects between group and sex for the sociability test (F = 3.479; p = 0.068), social novelty (F = 0.666; p = 0.418), sociability index (F = 1.424; p = 0.239), or novelty index (F = 4.056; p = 0.050). These results suggest that the vaccine did not exhibit differential effects on these social behavior parameters when considering the sex of the subjects.

However, a notable exception was observed in the latency time to fall, where a significant group by sex interaction emerged (F = 5.059; p = 0.029). Specifically, in the prenatal BNT162b2 vaccine group, a pronounced disparity was detected between the sexes (p < 0.001), unlike in the saline control group (p = 0.907). This particular finding points to sex-specific effects of the vaccine on motor performance, which are evident in the vaccinated group but not in the saline group.

The absence of a significant sex-treatment interaction in most of the two-way ANOVA parameters might suggest that the overarching impact of the vaccine does not differ markedly between males and females regarding social behavior. Nonetheless, the distinct male-specific outcomes, especially in motor performance, necessitate further research. These results illuminate the complexity of the vaccine’s neurodevelopmental impacts and highlight the imperative for a nuanced interpretation.

However, the distinct outcomes observed in the male-only analysis point towards the necessity of further research, particularly focusing on sex-specific responses to the vaccine. While the vaccine’s influence appears not to be uniform across sexes, the significant male-specific effects observed, particularly in motor performance, are significant and warrant further exploration. These findings accentuate the necessity of future studies with larger sample sizes and more precise methodologies to dissect the intricacies of sex-specific responses to the vaccine. A deeper understanding of these effects is crucial for developing a comprehensive perspective on the vaccine’s safety profile and its implications for neurodevelopmental health.

However, it’s vital to note these insights come primarily from animal studies, limiting their direct applicability to humans. Further human clinical research is required to confirm these findings. Longitudinal studies on the long-term impacts of COVID-19 vaccines, especially considering gender differences, are crucial to gauge vaccine safety and risks. Deepening our understanding of these neurobehavioral changes will guide intervention strategies and improve outcomes post-vaccination.

In conclusion, our study presents evidence that the COVID-19 mRNA BNT162b2 vaccine impacts the WNT pathway and BDNF levels in rats, with particularly pronounced effects observed in males. These male-specific outcomes, including autism-like behaviors, reduced neuronal counts, and impaired motor performance, emphasize the potential neurodevelopmental implications of the vaccine, aligning with existing literature on the roles of the WNT pathway and BDNF signaling in neurodevelopmental disorders.

It’s imperative to recognize the limitations of our research, given that it relies on animal models. Caution should be exercised in generalizing these results to humans. Further rigorous clinical studies are required to confirm these observations in human populations and to ascertain the exact mechanisms at play.

Given the public health significance of understanding the effects of COVID-19 vaccination, especially during pregnancy, comprehensive studies are vital. These should weigh the benefits and potential risks of vaccination, focusing on ensuring optimal neurodevelopmental outcomes. Our findings underscore the importance of continued research in this domain to guarantee the safety and well-being of all individuals, especially those who are pregnant and their offspring.

In summary, this study provides valuable insights into the effects of the COVID-19 mRNA BNT162b2 vaccine on the WNT pathway and BDNF levels, particularly in relation to neurodevelopmental outcomes. The observed male-specific vulnerability and the convergence with existing literature support the involvement of these molecular pathways in neurodevelopmental disorders. However, further research is warranted to validate these findings in human populations and to unravel the complex mechanisms underlying the observed effects. The ultimate goal is to ensure the safety and well-being of individuals receiving COVID-19 vaccination, particularly during pregnancy, while minimizing potential risks to neurodevelopment.