Based on the results of a previous study (
Partevian et al., 2023), the initial phase of the current work involved a detailed phenotypic analysis of
D. rerio following the suppression of
anxa2a gene expression by injecting 3 ng of morpholino at 2 days post-fertilization (dpf).
Upon analyzing the phenotype of fish in the Mo-Anxa2a group, the following changes were observed in multiple experimental replicates only in this group: hydrocephalus (enlargement of the hindbrain ventricle (HBV)) and ocular and olfactory pit malformations (Table 4). It should be noted that some fish (30.14%) with specific morphogenesis disruptions in the head region also exhibited developmental abnormalities of the tail and notochord. Furthermore, deformations of the body, tail and notochord (without concomitant specific abnormalities in the head region) were also observed in control and experimental groups, indicating the nonspecificity of the observed phenotype. Therefore, for subsequent expression analysis, we selected individuals with specific head region deformations from the Mo-Anxa2a group and individuals with a normal phenotype from the Co-Mo group.
Next, a detailed analysis of the severity of hydrocephalus was performed in fish with specific deformations using the previously developed approach
(Partevian et al., 2025). Two subgroups of fish were identified: those with pronounced (or severely expressed) enlargement of the HBV area (2.49±0.43 fold increase, 42.55% of individuals) and those with mild HBV enlargement area (1.33±0.23 fold increase, 17.73% of individuals) compared to the control group (Fig 1).
Western blot analysis was performed to determine the efficiency of suppressing Anxa2a expression at the translational level using the Morpholino oligonucleotide (MO). The analysis demonstrated that Anxa2a protein expression was reduced twofold in the Mo-Anxa2a group (Me = 0.17±0.04) compared to the Co-Mo group (Me = 0.36±0.08) at 2 dpf (Fig 2). This reduction confirms the effectiveness of the 3 ng MO injection in suppressing Anxa2a protein expression. Crucially, the expression of this gene did not change statistically at the mRNA level, thereby confirming that the suppression of Anxa2a protein occurs at the translational level.
Then, an analysis of the change in the relative mRNA expression level of the main Anxa2a partner genes (
egf, egfra, ptena, ptenb, gfap, psen1, psen2) was conducted. The analysis was performed both for the entire Mo-Anxa2a group of fish with specific deformations and separately for fish with pronounced and mild hydrocephalus. The results are presented in Table 5.
As can be seen from Table 5, all groups of fish with suppressed Anxa2a expression exhibited pronounced, statistically significant decrease in the relative expression level of the following genes:
ptena, ptenb, gfap, psen1 and
psen2. In contrast, no statistically significant changes in the expression of the
ahnak gene were observed across any of the experimental groups. Furthermore, fish with severe hydrocephalus also showed a decline in the expression of
egf and
egfra, suggesting a significant contribution of these genes to the development of severe hydrocephalus.
Currently, there is a growing interest in investigating the role of ANXA2 in nervous system development under both physiological and pathological conditions. However, most studies have been conducted using cell cultures and murine models. In the
D. rerio model, the roles of the ANXA2 orthologs - Anxa2a and Anxa2b - in nervous system function remain poorly understood
(Quoseena et al., 2020). Therefore, the aim of our study was to evaluate the impact of Anxa2a knockdown on nervous system functioning in
D. rerio.
In the present study, it was demonstrated that suppression of Anxa2a expression at the protein level leads to the development of morphological changes of varying severity in the nervous system of
D. rerio. Such phenotypic alterations may be associated with the involvement of ANXA2 in a range of biological processes potentially significant for nervous system function. ANXA2 is involved in these processes both independently and through interaction with its primary partners
(Partevian et al., 2025).
To gain a more complete understanding of the effect of altered Anxa2a expression on the nervous system of
D. rerio, we carried out a selection of the main partners (proteins and their encoding genes) of Anxa2a using the Pathway Studio v.12.4.0.3 program. This analysis enabled the identification of 31 ANXA2 partners. These partners were analyzed for compliance with the following criteria: interaction with ANXA2 shown in at least two separate publications, association with neurological diseases and/or expression in the nervous system and presence of orthologs in
D. rerio. Ultimately, we identified six partners that satisfied all three conditions, namely, GFAP (Gfap), EGFR (Egfra), PTEN (Ptena and Ptenb), EGF (Egf) and S100A10 (S100a10a and S100a10b) and AHNAK (Ahnak). PSEN1 (Psen1) and PSEN2 (Psen2) satisfy two of the three criteria: they are associated with Alzheimer’s disease and have orthologs in fish. Furthermore, it has been previously shown that suppression of
psen1 and
psen2 gene expression by MO in
D. rerio results in the development of hydrocephalus (
Nornes et al., 2009). In addition, an interaction between PSEN1 and ANXA2 was demonstrated in the study by
Bustos et al. (2017). Consequently, these genes were also included in our subsequent expression analysis.
Further analysis of the Anxa2a partners using the STRING resource (v. 12.0) indicates close regulatory interactions among the selected partners.
As shown in Fig 3, Egf and Egfra are targets regulated by Anxa2a. Based on the data presented in Table 5, a statistically significant decrease in
egf and
egfra expression was observed in the overall Mo-Anxa2a group, as well as in
D. rerio with severe hydrocephalus. The most pronounced drop in expression was observed for the
egfra gene. Based on this finding, it can be hypothesized that
egfra makes a direct contribution to the development of the morphological changes we observed, which are associated with nervous system development. These changes may be linked to the fact that ANXA2 is involved in the clathrin-dependent internalization of EGFR at the cell membrane and promotes the activation of its downstream signaling pathways (
Grewal and Enrich, 2009;
Chaudhary et al., 2014).
It has also been shown in
D. rerio that the activation of the MAPK/ERK signaling pathway
via HB-EGF/EGFR is essential for the proliferation of basal cells and the restoration of olfactory epithelium neurogenesis after injury
(Sireci et al., 2024). These data indirectly support our findings. It can be hypothesized that the olfactory pit deformations we observed may be linked to the disruption of the PI3K/Akt and MAPK/ERK signaling pathways, where EGFR is an important participant. We hypothesize that the decrease in Anxa2a expression might impair EGFR internalization and its subsequent intracellular vesicular transport within the cell. The accumulation of EGFR on the cell surface then leads to a reduction in its mRNA level.
Furthermore, a slight but significant decrease in the expression of the main ligand of
egfra -
egf was also observed in the group with severe hydrocephalus. It is possible that such a decrease in expression may lead to the disruption of metabolic processes associated with cell division, thereby contributing to the formation of severe hydrocephalus.
In
D. rerio developing hydrocephalus, a decrease in the expression of the
ptena and
ptenb genes was observed. The PTEN protein (orthologs in
D. rerio: Ptena and Ptenb) is part of the same signaling pathway as the EGF, EGFR and ANXA2 proteins
(Castaldo et al., 2019). Previously,
Croushore et al., (2005) investigated the expression patterns of
ptena and
ptenb in
D. rerio under normal and pathological conditions. It was demonstrated that
ptena and
ptenb expression in
D. rerio is most pronounced in the central nervous system, eyes, pharyngeal arches and pectoral fins by 48 hours post fertilization
(Croushore et al., 2005). Despite similar distribution patterns within the organism, distinct phenotypes were observed upon the suppression of these proteins. Specifically, Ptena knockdown led to deformations in the head and notochord regions, as well as impaired vasculogenesis and otolith development. Conversely, Ptenb knockdown was associated with brain ventricle enlargement (domed head, hydrocephalus), as well as tail and yolk deformations
(Croushore et al., 2005). It is possible that the suppression of
anxa2a gene expression leads to the formation of a specific phenotype via the downregulation of
ptena and
ptenb and the activation of the PI3K/Akt pathway, followed by the development of oxidative stress. For instance, reduced
ptena expression may result in otolith and tail deformations, while reduced
ptenb expression may lead to ventricular enlargement, a finding confirmed by previous independent studies. However, this hypothesis requires further experimental validation.
A decrease in the expression of
psen1 and
psen2 genes was also observed in both groups. It is noteworthy that PSEN1, in conjunction with ANXA2, participates in autophagosome formation and the subsequent cleavage of amyloid-beta
(Bustos et al., 2017). Furthermore, PSEN1, as a component of the γ-secretase complex, is involved in the processing of the Notch protein
(Hurley et al., 2023). A previously published independent study demonstrated that the suppression of Psen1 and Psen2 expression using MOs in
D. rerio resulted in the development of hydrocephalus and a reduction in melanocyte number
(Nornes et al., 2009). Based on these findings, it can be hypothesized that the downregulation of Anxa2a may indirectly lead to the development of hydrocephalus
via the reduction of
psen1 and
psen2 expression. Such changes in expression may lead to further disruption of autophagy, the Notch signaling pathway and the formation of a specific phenotype.
A statistically significant decrease in the relative mRNA levels of the
gfap gene, which encodes the glial fibrillary acidic protein (Gfap), was identified in all groups. This protein is a constituent of intermediate filaments in astrocytes. It is known that ANXA2 is also expressed in astrocytes and can regulate their proliferation
(Chen et al., 2017). Previous studies have shown that ANXA2 can interact with GFAP and promote its assembly into glial filaments in a Ca
2+-dependent manner
(Garbuglia et al., 1995). We hypothesize that the reduction in
gfap expression, coupled with the inhibition of Anxa2a translation, impairs the assembly of glial filaments. This disruption of astrocyte function may serve as a primary mechanism underlying the development of inflammatory hydrocephalus.
Also, no statistically significant decrease in the expression of the
ahnak gene, whose protein product forms a complex with Anxa2a, was observed in either group
(Jin et al., 2020). Possibly, this is due to the fact that the metabolic pathways jointly involving Anxa2a and Ahnak are not implicated in the development of hydrocephalus.
Thus, the data obtained suggest that an altered Anxa2a level may lead to pathological processes in the nervous system and the development of hydrocephalus through the modification of metabolic pathways associated with its main partners:
egf and
egfra (cell division),
ptena and
ptenb (reactive oxygen species formation),
psen1 and
psen2 (autophagy) and
gfap (regulation of inflammation in astrocytes). However, this remains a hypothesis and requires further research to be fully validated.