Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Immunol., 10 October 2025

Sec. Nutritional Immunology

Volume 16 - 2025 | https://doi.org/10.3389/fimmu.2025.1655321

Immunomodulatory inhibition of osteoclastogenesis by a marine microalgal ethanol fraction targeting T-cells, antigen presentation, and macrophage fate

Alessio Carletti,&#x;Alessio Carletti1,2†Katia Pes&#x;Katia Pes1†Marco Tarasco&#x;Marco Tarasco1†Joana T. Rosa,Joana T. Rosa1,3Sunil Poudel,&#x;Sunil Poudel1,2†Hugo Galvo PereiraHugo Galvão Pereira4Bruno LouroBruno Louro1M. Leonor Cancela,,M. Leonor Cancela1,2,5Vincent Laiz,Vincent Laizé1,3Paulo J. Gavaia,,*Paulo J. Gavaia1,2,4*
  • 1Centro de Ciências do Mar do Algarve (CCMAR/CIMAR LA), Campus de Gambelas, Universidade do Algarve, Faro, Portugal
  • 2Faculty of Medicine and Biomedical Sciences, University of Algarve, Faro, Portugal
  • 3S2AQUA - Collaborative Laboratory, Association for a Sustainable and Smart Aquaculture, Olhão, Portugal
  • 4GreenCoLab – Associação Oceano Verde, University of Algarve, Faro, Portugal
  • 5Algarve Biomedical Center, University of Algarve, Faro, Portugal

Background: Targeting immune pathways to prevent bone loss represents a promising, yet underexplored therapeutic strategy.

Methods: An ethanol-soluble fraction derived from the freeze-dried biomass of the marine microalga Skeletonema costatum (SKLT) was tested for its ability to modulate immune responses and inhibit osteoclastogenesis. Its effects were evaluated in a zebrafish model of bone regeneration, a medaka model of RANKL-induced osteoporosis, and in vitro using murine RAW 264.7 macrophages. Transcriptomic profiling of regenerating fin blastemas at 24 hours post-amputation was performed to identify the affected molecular pathways.

Results: In zebrafish, SKLT treatment suppressed the recruitment of osteoclast precursors and altered mineralization dynamics. Transcriptomic profiling revealed downregulation of genes involved in inflammation, antigen presentation, T-cell activation, and macrophage commitment towards osteoclastogenesis, accompanied by reduced expression of chemokines and cytokines that promote osteoclast precursor recruitment and fusion. In medaka, SKLT significantly reduced vertebral bone loss and enhanced neural arch mineralization in larvae with high RANKL expression. In vitro, SKLT inhibited proliferation and osteoclastic differentiation of murine RAW 264.7 macrophages exposed to RANKL without inducing cytotoxicity.

Conclusion: These findings identify S. costatum as a source of bioactive immunomodulatory compounds capable of interfering with key osteoimmune mechanisms. Beyond providing proof of concept for their therapeutic potential in bone erosive disorders, this work opens avenues for isolating and characterizing the active molecules, optimizing their delivery, and evaluating their efficacy in preclinical mammalian models. Such strategies could expand the repertoire of safe, nutraceutical-based or adjuvant therapies for osteoporosis and other inflammation-driven skeletal diseases, complementing and potentially enhancing current antiresorptive and anabolic treatments.

Introduction

Bone erosive pathologies such as osteoporosis represent urgent medical challenges requiring novel therapeutic approaches (1, 2). These disorders arise from the disruption of the equilibrium between bone formation and resorption, ultimately resulting in mineral loss and deterioration of bone structure (3, 4). Current therapeutic strategies rely on a limited number of pharmaceutical agents that either stimulate osteoblast-mediated bone formation (i.e. osteoanabolic compounds) or inhibit osteoclast-driven bone resorption (i.e. antiresorptive drugs). Examples of osteoanabolic drugs are parathyroid hormone analogs (e.g., teriparatide and abaloparatide) (57), whereas antiresorptive agents include bisphosphonates and RANKL-monoclonal antibodies (e.g., Denosumab) (810). More recently, dual-action drugs such as the anti-sclerostin antibody Romosozumab have emerged as promising options for treating bone loss (11).

Although these drugs demonstrate short-term efficacy, they are not without limitations: their effectiveness diminishes over time (12), and they carry the risk of adverse events—ranging from rare complications such as osteonecrosis of the jaw for antiresorptives (13, 14), to potential associations with osteosarcoma for anabolic agents (57).

This therapeutic gap has fueled demand for novel approaches that offer long-lasting, safer alternatives. Accordingly, research has intensified into the mechanisms driving bone erosive diseases, as well as the identification of novel osteoactive compounds with improved safety.

A growing body of research underscores the essential role of the immune system in regulating bone metabolism, a concept that has crystallized into the field of osteoimmunology (15).

Immune cells, particularly macrophages and T-cells, exert critical regulatory control over bone homeostasis (16, 17). Macrophages share a common monocytic lineage with osteoclasts (18), and represent a conserved pool of osteoclasts progenitors in many vertebrates, including fish, mice and humans (19, 20). Together with macrophages, T-cells are the main cellular regulators of inflammation, which is a potent inducer of bone resorption. Inflammation is required for the receptor activator of the nuclear factor kappa-Β ligand (RANKL) signaling pathway, the central molecular driver of osteoclast differentiation (2125). Chronic inflammation is widely recognized as a key root cause of bone loss in various erosive pathologies (26). For example, postmenopausal osteoporosis is increasingly acknowledged as an inflammatory disease and was (re)named immunoporosis (2730), where activated T-cells establish chronically increased bone resorption by stimulating osteoclast differentiation through RANKL-dependent and RANKL-independent pathways (31, 32).

The leading role of the immune system in the pathophysiology of osteoporosis and other bone disorders raises the question of whether immunomodulatory therapies could effectively prevent bone loss in patients suffering from such disorders. Surprisingly, few investigations have focused on the potential of targeting immune cells to control bone loss.

The search for novel therapeutic approaches has recently refreshed the role of natural compounds, and marine-derived molecules have become promising assets in drug discovery (33, 34).

This is particularly true for metabolic bone disorders, as an increasing number of studies have identified molecules and extracts obtained from marine organisms with immunomodulatory and osteoprotective bioactivities (35). Among those, microalgae are attracting significant interest owing to a combination of their capacity to synthesize bioactive compounds, technologically advanced cultivation technologies (3638), and amenability to genetic engineering (39, 40).

Among the recently emerging microalgal species, Skeletonema costatum, a marine diatom commonly cultivated in Europe and approved for human consumption (41), has shown promising potential as a source of compounds with biological activity. In detail, the ethanol-soluble compounds extracted from dry biomass of this microalga were characterized by a potent anti-inflammatory activity emerged from in vitro studies (42, 43). In the mentioned works, authors performed chemical profiling of ethanolic extract from Skeletonema spp. and found them rich in polyunsaturated fatty acids, saturated fatty acids, monounsaturated fatty acids, long-chain aliphatic hydrocarbons, long-chain alcohols, esters, ethers and sterols. Additionally, in a recent work published from our laboratory, ethanolic extract of S. costatum had a pro-osteogenic effect in zebrafish larvae and juveniles, as well as fish osteoblastic cell lines (44).

The aim of this study was to evaluate the immunomodulatory and anti-osteoclastogenic potential of the ethanol-soluble fraction of S. costatum, and investigate the underlying molecular mechanisms via tissue-specific transcriptomics. To assess the fraction’s osteoactivity, a zebrafish model of bone regeneration (45), a medaka genetic model of RANKL-induced osteoporosis (46), a murine macrophage cell line, were implemented.

Results

Skeletonema suppresses osteoclast recruitment and alters bone regeneration dynamics in zebrafish

The impact of the ethanolic fraction of Skeletonema costatum (hereafter referred to as SKLT) on bone regeneration was initially evaluated using the zebrafish caudal fin regeneration assay. This model was selected due to its simplicity, high reproducibility, and suitability for studying the coordinated actions of multiple cell types during de novo bone formation. It also allows for precise in vivo cell tracking (45, 47, 48). Following caudal fin amputation, adult zebrafish were exposed to SKLT by immersion in water supplemented with the fraction. Regenerative outgrowth, mineralization, and morphological patterning of the regenerating bony rays were assessed at a single time point: 5 days post-amputation (5 dpa, Figures 1A, B). Quantitative analysis showed a trend toward reduced regenerative area in SKLT-treated animals, though this reduction did not reach statistical significance (Figure 1C). In contrast, SKLT exposure resulted in a significant increase in the area of mineralized tissue (Figure 1D), suggesting enhanced mineral deposition. While the width of regenerated rays remained unchanged (Figure 1E), a marked distalization of ray bifurcation points along the proximal–distal axis was observed in SKLT-treated fish (Figure 1F).

Figure 1
An illustration showing a comparison of fish fin regeneration. Panel A displays a fin schematic with color-coded segments representing different regions. Panel B highlights a section with labels indicating ray, cut, and bifurcation areas. Panels C-F present graphs comparing control and SKLT groups for regenerated area, mineralized area, ray width, and ray bifurcation, with SKLT showing significant differences.

Figure 1. SKLT shifts the mineral equilibrium of regenerating rays toward pro-mineralogenic outcomes. Caudal fin regeneration and mineralization were assessed in young adults exposed to the ethanolic fraction of S. costatum (SKLT) or 0.1% ethanol (control). (A) Representative image of a caudal fin illustrating the morphometric measurements used to assess the regenerative and mineralogenic effects of SKLT, i.e., ray width (RAY), stump length (STU), regenerated area (REG) and mineralized area (MIN). (B) Representative image of an alizarin red S (AR-S) stained fin ray illustrating the morphometric measurements used to assess the patterning effect of SKLT, i.e., amputation plan (CUT), ray width before amputation (RAY), average ray width calculated as the average width of the first two segments after the amputation plan (W), length of bifurcation (BIF) and total length of regenerated ray (TOT). (C–F) Effects of SKLT on the regenerated area (C), mineralized area (D), average ray width (E) and ray bifurcation (F). Normality was tested through an Anderson–Darling test (p < 0.05). Statistical differences were tested through an unpaired t test (p < 0.05) or a nonparametric Mann–Whitney test (p < 0.05) whenever the data distribution resulted non-normal. Asterisks indicate values significantly different at p < 0.0005 (***) and p < 0.0001 (****). Scale bars are 1 mm in (A) and 250 µm in (B).

To investigate whether the observed effects on mineralization were associated with antiresorptive activity, the dynamics of osteoclast precursors—identified as rounded cathepsin K (ctsk)-expressing cells – in Tg(Ola.ctsk:FRT-DsRed-FRT-Cre, myl7:EGFP)mh201 transgenic zebrafish (49), hereafter referred to as Tg(ctsk:DsRed), was monitored in fish continuously exposed to SKLT for 10 days following caudal fin amputation (Figure 2A).

Figure 2
Panel A shows a fluorescent microscopy image highlighting ctsk expression in a sample. Panels B and C compare control and SKLT conditions at 24, 48, and 72 hours post-amputation (hpa), showing red fluorescent areas. Panel D presents a violin plot comparing normalized ctsk positive areas between control and SKLT at multiple time points from 24 to 240 hpa, with significant differences marked at certain time points. Scale bars are present.

Figure 2. SKLT inhibits the recruitment of osteoclast precursors to the caudal fin blastema. Osteoclastic recruitment was determined from DsRed fluorescence signals observed in the regenerating caudal fin of Tg(ctsk:DsRed) transgenic fish. (A) Representative image illustrating the regenerated area (dotted line) and ctsk+ area (ctsk, red signal). Time-course of the ctsk+ area in regenerating fins of SKLT-treated (B) and control fish (C). (D) Quantification of the ctsk+ area in regenerating fins of SKLT-treated and control fish from 24 to 240 hours post-amputation (hpa). At each timepoint, significant differences were tested through an unpaired t test (p < 0.05) or a non-parametric Mann-Whitney test (p < 0.05) whenever the data distribution resulted non-normal. Asterisks indicate values significantly different at p < 0.0001 (****). Scale bars are 500 µm.

In vehicle-treated controls, the number of ctsk+ cells peaked between 24- and 48 hours post-amputation (hpa) (Figures 2B, D), followed by a progressive decline in cell density within the regenerated area. In contrast, SKLT-treated fish showed significantly fewer ctsk+ cells recruited to the blastema at both 24 and 48 hpa (Figures 2C, D). A delayed peak was observed at 72 hpa, after which ctsk+ cell counts were comparable between SKLT-treated and control animals.

Notably, elongated, tubular ctsk-expressing osteoclasts resembling osteolytic tubules (OLTs), previously described in regenerating fin tissue (48), appeared adjacent to the regenerated rays by 96 hpa in control fish, but were delayed until 120 hpa in SKLT-treated fish (Supplementary Figure S1).

This reduction in osteoclast precursor recruitment correlated with a significant decline in early osteoclastic activity. TRAP (tartrate-resistant acid phosphatase) staining, a marker of mature osteoclast function, was nearly absent at 24 hpa in SKLT-exposed animals but returned to levels comparable to controls by 240 hpa (Figures 3A, B). Furthermore, in the 10-day exposure experiment, SKLT-treated fish exhibited a sustained reduction in the area of regenerative outgrowth beginning at 72 hpa and persisting through 240 hpa (Figure 4). This suggests that the ethanolic fraction may interfere with molecular pathways governing blastema expansion during the regeneration process.

Figure 3
Panel A shows two histological images comparing a control sample and SKLT sample, both highlighting differences in staining patterns and color intensity. Panel B is a violin plot showing the normalized TRAP positive area for control and SKLT at 24 and 240 hours post application (hpa). Control displays wider distribution at both time points, with significant differences indicated at 24 hpa. Control is denoted by black, SKLT by blue.

Figure 3. SKLT suppresses bone resorption activity at 24 hours post-amputation. Osteoclast activity was assessed through TRAP (tartrate-resistant acid phosphatase) staining of the regenerating caudal fin of adult zebrafish exposed to the ethanolic fraction of S. costatum (SKLT) or 0.1% ethanol (control) at 24 and 240 hours post-amputation (hpa). (A) Representative images showing TRAP staining at 24 hpa in control and SKLT-treated fish. (B) Quantification of TRAP+ area at the early regeneration (24 hpa) and late regeneration (240 hpa) stages. Statistical differences were tested through an unpaired t test (p < 0.05). Asterisks indicate values significantly different at p < 0.0001 (****). Scale bars are 500 µm.

Figure 4
Panel A and B show images of fish fins at different time points: 24, 48, 72, 96, 120, and 240 hours post-amputation (hpa). Panel A is labeled “Control” and Panel B “SKLT.” Each image displays regeneration progression with clear markings indicating the regrowth area. Panel C presents a graph comparing regenerated area percentages over time between Control (black) and SKLT (blue). The x-axis represents time in hpa, and the y-axis shows regenerated area percentage. The graph indicates statistical significance at various time points, marked by asterisks.

Figure 4. SKLT reduces the overall regenerative performance. Regeneration of the caudal fin was assessed in adult zebrafish exposed to the ethanolic fraction of S. costatum (SKLT) or 0.1% ethanol (control) from 24 to 240 hours post amputation (hpa). (A, B) Representative images depicting the time-course of caudal fin regeneration (dotted line indicates regenerate) in control (A) and SKLT-treated (B) fish. (C) Quantification of the regenerated area at relevant time points. Statistical differences were tested through an unpaired t test (p < 0.05) or a non-parametric Mann-Whitney test (p < 0.05) whenever the data distribution resulted non-normal. Asterisks indicate values significantly different at p < 0.05 (*), p < 0.01 (**). Scale bars are 500 µm. Picture illumination was artificially adjusted for better visualization.

SKLT suppresses inflammation, T-cell activation, and antigen presentation in the regenerating fin blastema

To investigate whether immunomodulatory mechanism underlies the effects observed in SKLT-treated zebrafish, tissue-specific transcriptomic analysis was performed on regenerating caudal fin blastemas collected at 24 hours post-amputation — the time point at which SKLT exerted the strongest suppression of osteoclast precursor recruitment.

Bulk RNA sequencing identified 720 differentially expressed genes (false discovery rate, FDR< 0.01), of which 361 were upregulated and 359 were downregulated (Supplementary Figure S2). Gene ontology (GO) analysis was performed using the DAVID platform, and significantly enriched Biological Process (BP) terms were filtered for uniqueness and plotted using REVIGO (Figure 5A). Heatmaps were then generated to visualize the expression patterns of genes associated with the most relevant enriched processes.

Figure 5
A multi-panel scientific image includes:  A) A bubble chart showing Gene Ontology biological processes in a semantic space, with bubble size indicating significance and color representing log value.  B) A heatmap titled “Immune system and chemokines” with hierarchical clustering. Colors represent expression levels from high (red) to low (blue).  C) A heatmap titled “Cell cycle and cytoskeleton organization,” also with hierarchical clustering and a similar color scheme.  D) A heatmap titled “Cell metabolism” with a similar layout and color coding. Each panel compares control and SKLT samples.

Figure 5. Analysis of the transcriptome of regenerating caudal fin blastema revealed that SKLT modulates immune response, cell cycle and metabolism. RNA-Seq analysis identified differentially expressed genes (DEGs) and enriched gene ontology terms for caudal fin blastemas at 24 hours post-amputation (hpa) in control and SKLT-exposed fish. (A) REVIGO plot of biological processes filtered for redundancy and dispensability. Cluster analysis and heatmap of differentially expressed genes associated with (B) immune system (GO:0002376), eicosanoid metabolism (GO:0006690), inflammatory response (GO:0006954), T-cell activation (GO:0042110), and cytokine-mediated signaling pathway (GO:0019221). (C) Cell cycle (GO:0022402) and cytoskeleton organization (GO:0007010); (D) glycolysis (GO:0006096), fatty acid metabolism (GO:0006099), tricarboxylic acid cycle (GO:0006631), and carboxylic acid metabolic process (GO:0019752).

Among the most significantly affected biological processes were those associated with both innate and adaptive immune responses. Genes involved in inflammatory signaling—including irg1l, nos2a, thbs1, and il20ra—as well as those linked to the interferon system (ifi30, ifng1r, irf4b, irf10) and the complement cascade (c4, c7a) were consistently downregulated in SKLT-treated blastemas (Figure 5B). Since interferon gamma (IFN-γ), primarily produced by Th1 cells (50), promotes antigen-dependent T-cell activation and the expression of pro-osteoclastogenic cytokines (50), genes involved in T-cell signaling and antigen presentation were examined. Indeed, multiple genes central to T-cell activation and MHC class II antigen presentation—including socs1a, itk, il7r, batf, cd4, cd74a, cd74b, and jak1 (30, 31)—were downregulated following SKLT exposure (Figure 5B).

Conversely, several immunoregulatory genes such as fn1a, fn1b (fibrin isoforms), and il13ra2 (interleukin 13 receptor alpha 2) were upregulated, consistent with a shift toward anti-inflammatory immune signaling. Chemokine-related paracrine communication appeared differentially modulated: genes encoding Ccl34 (ccl34a.3) and Cxcl8 (cxcl8a) were upregulated, whereas Ccl19 (ccl19a.1) and the chemokine receptor Cxcr3 (cxcr3.1) were downregulated.

Together, these transcriptomic data support the hypothesis that S. costatum induces complex immunomodulatory effects in the regenerating fin blastema, suppressing both inflammation and T-cell–mediated antigen presentation. This shift in immune signaling likely contributes to reduced recruitment and differentiation of osteoclast precursors.

SKLT inhibits retinoic acid synthesis and cell proliferation while upregulating metabolic gene expression

Cluster analysis of genes involved in cell cycle regulation and cytoskeletal organization revealed that SKLT exposure may suppress cell proliferation in fin blastemas at 24 hpa (Figure 5C).

Downregulated genes included those associated with proliferation (ttk, cdca8, dlgap5, tacc3, nuf2, ndc80), mitosis (bub1, ccnb1, ccnb3, ccna2, cdk1, tmpoa, mad2l1, foxm1), and cytokinesis (ect2, aurka, aurkb, kif23, racgap1, anln), indicating a broad inhibitory effect of SKLT on cell cycle progression. In contrast, transcriptomic analysis revealed increased expression of genes involved in energy metabolism (Figure 5D), suggesting a metabolic shift. Glycolysis-related genes—including eno3, eno1a, gapdhs, aldocb, tpi1b, and gpia—as well as fatty acid oxidation genes (gcdh, asah1a, sdhb, sdhaf2, idh1, sdha, idh3a) were significantly upregulated in SKLT-treated blastemas.

Interestingly, aldh1a2, encoding aldehyde dehydrogenase 1 family member A2, was downregulated in SKLT-treated fish. These molecular signatures are consistent with the reduced regenerative outgrowth observed at later stages in SKLT-treated animals (Figure 4).

SKLT protects against RANKL-induced bone loss in a medaka model of osteoporosis

Whether the immunomodulatory properties of SKLT could prevent bone loss in a disease setting was investigated in the inducible osteoporosis medaka model Tg(rankl:HSE: CFP)TG1135, in which bone resorption is triggered by conditional overexpression of RANKL (19). Successful induction is monitored via cyan fluorescent protein (CFP) fluorescence, which is co-expressed with RANKL under the control of a bidirectional heat-shock promoter. Higher CFP intensity corresponds to greater RANKL expression and increased osteoclastic activity, resulting in bone loss.

To evaluate the protective effects of SKLT, medaka larvae were heat-shocked and then exposed to SKLT or vehicle for 6 days. RANKL expression levels were estimated by measuring CFP fluorescence intensity (Figures 6A, B). Because the distribution of CFP fluorescence intensity was bimodal or right-skewed, larvae were stratified into two subgroups based on CFP intensity (Figure 6C). Bone loss was assessed by measuring the mineralized area of the vertebral column following alizarin red staining. In the subgroup with low CFP fluorescence—indicative of low RANKL expression—no difference in vertebral mineralization was observed between control and SKLT-treated fish (Figure 6D). However, in fish with high CFP (and thus high RANKL) expression, exposure to SKLT significantly increased vertebral mineralization (Figure 6E).

Figure 6
Composite image showing various analyses of fish samples. Panels A and B display fish with areas labeled AR-S and CFP, comparing control and SKLT specimens with red and blue fluorescence. Panel C presents a violin plot of Low and High CFP between control and SKLT groups. Panels D and E show bar graphs of mineralized area percentages with statistical significance. Panels F and G depict magnified sections illustrating mineralization differences in control and SKLT samples. Panel H is a stacked bar chart comparing mineralized neural arches, with key indicating quantities zero to three or more.

Figure 6. SKLT protects against bone loss in a medaka model of RANKL-induced osteoporosis. Mineralization of the vertebral column was assessed in adult medaka induced for osteoporosis and exposed to the ethanolic fraction of S. costatum (SKLT) or 0.1% ethanol (control). (A, B) Representative images illustrating the total body area (dotted line) and rankl-CFP+ area in SKLT-treated (A) and control (B) fish. (C) Frequency distribution of the CFP mean intensity. (D, E) Quantification of the mineralized area of the vertebral column in the different fish clusters. (F) Representative images depicting alizarin red stained vertebral columns of control (upper image) and SKLT-exposed (lower image) fish in the High CFP intensity group. (G) Representative images depicting the number of mineralized neural arches in control (upper image) and SKLT-exposed (lower image) fish in the High CFP intensity group. (H) Number of mineralized neural arches in SKLT-exposed and control fish in the High CFP intensity group. Statistical differences were tested through an unpaired t test (p < 0.05). Scale bars are 360 µm in (A, B) and 730 µm in (J, K).

Additionally, these fish exhibited a greater number of mineralized neural arches (Figures 6F–H). These findings suggest that SKLT provides protection against RANKL-induced bone loss in cases of strong RANKL induction, but has no effect when RANKL expression is low. Thus, SKLT appears to counteract the osteoporotic phenotype in a severity-dependent manner.

SKLT inhibits osteoclastic differentiation in murine macrophages without inducing cytotoxicity

The ability of SKLT to inhibit osteoclast differentiation was evaluated using the murine macrophage cell line RAW 264.7, which undergoes osteoclastic differentiation upon exposure to RANKL. To assess the effects of SKLT on cell viability and proliferation, RAW 264.7 cells were exposed to various concentrations of SKLT (50, 100, and 200 µg/mL) for 24 and 48 hours. No cytotoxic effects were observed at any concentration (Figure 7A). However, significant inhibition of cell proliferation was detected after 1, 3, 5, and 7 days of treatment, even at the lowest concentration (50 µg/mL) (Figure 7B). To evaluate osteoclastogenesis, cells were co-treated for 6 days with RANKL (50 ng/mL) and SKLT (50 µg/mL), and stained with DAPI and TRAP (Figures 7C, D). No large multinucleated (>2 nuclei) TRAP+ cells were detected in the negative control group (no RANKL), while the RANKL-treated group produced an average of 14.3 ± 2.6 osteoclasts per 5.0x field. Co-treatment with SKLT significantly suppressed osteoclast differentiation, reducing the number of multinucleated TRAP+ osteoclasts to 0.78 ± 0.67 per 5.0x field (Figures 7C, E).

Figure 7
Graphs and images analyzing the effects of SKLT on cells. (A) Bar graphs showing relative cytotoxicity at 24 and 48 hours, indicating no significant toxicity at different SKLT concentrations. (B) Line graph illustrating cell proliferation over seven days, with significant differences between SKLT concentrations and control. (C) Bar graph displaying the number of multinucleated osteoclasts, showing a significant reduction with SKLT treatment. (D) Microscopic images with DAPI and TRAP staining highlighting cellular morphology and osteoclast activity. (E) Microscopic images comparing untreated, RANKL-treated, and RANKL with SKLT treatment, illustrating changes in cellular structure.

Figure 7. SKLT inhibits cell proliferation and osteoclastic differentiation in murine macrophages. The cytotoxicity, proliferation and osteoclastic differentiation of RAW 264.7 cells were assessed upon exposure to the ethanolic fraction of S. costatum (SKLT) or 0.1% ethanol (control). (A) Cytotoxicity (LDH assay) in cells exposed for 24 h (left) or 48 h (right) to different concentrations of SKLT. (B) Proliferation (XTT assay) of cells exposed to different concentrations of SKLT. (C) Number of multinucleated osteoclasts per field (magnification 5.0X). (D) Representative fluorescence image of DAPI-stained multinucleated osteoclasts (left) and brightfield images of TRAP+ osteoclasts (right). (E) Representative phase-contrast images of the 3 experimental groups. Arrowheads indicate multinucleated osteoclasts. Differences between each treatment and the control were tested via one-way ANOVA followed by Dunnett’s multiple comparison test (p < 0.05). Asterisks indicate values significantly different at p < 0.0001 (****). Scale bars are 150 µm.

These results confirm that SKLT exerts anti-proliferative and anti-osteoclastogenic effects on murine macrophages without causing cytotoxicity, supporting its cross-species activity and potential as an immunomodulatory inhibitor of osteoclast formation.

Discussion

The immune system is increasingly recognized as a key driver of bone erosive pathologies, including primary osteoporosis, where dysregulation of innate and adaptive immune cells contributes to sustained bone resorption (30, 51).

Although no currently available therapeutics act via immunomodulation, growing evidence supporting immune involvement in bone homeostasis has laid the foundation for novel therapeutic approaches targeting immune pathways (52, 53).

Here, the immunomodulatory potential of an ethanol-soluble fraction from the marine microalga Skeletonema costatum in treating bone loss was explored. In a zebrafish fin regeneration model (45), exposure to this fraction (SKLT) led to reduced ray bifurcation. During fin regeneration, the exact positioning of the bifurcation point in regenerating rays depends on the conjunct activity of re-differentiated osteoblasts (51, 54, 55), and osteolytic tubules (OLTs) (48). A shift in branching position indicates that SKLT modulates the morphogenetic pattern of ray formation, potentially by altering the balance between osteogenic and osteoclastic activity, as previously observed for other osteo-modulatory drugs (48). Previous work has demonstrated that osteoanabolic and antiresorptive agents alter ray bifurcation patterns, making this model a suitable in vivo screening platform for compounds that affect bone remodeling (48). Therefore, the distalization of ray bifurcations in SKLT-treated fins reflects either enhanced osteogenesis or suppressed resorption.

Given prior reports of anti-inflammatory (42, 43), and pro-mineralogenic (44) activities in S. costatum extracts, and the central role of inflammation in osteoclastogenesis (16), we hypothesized that SKLT impairs osteoclast function through immunomodulation.

This was supported by reduced expression of cathepsin K (an osteoclast precursor marker) and weaker TRAP activity (a marker of mature osteoclasts), consistent with previous findings (48).

These data suggest that SKLT delays osteoclast precursor recruitment or differentiation during early fin regeneration. Interestingly, SKLT also impaired late-stage regenerative outgrowth.

Since inflammation is essential for initiating regeneration across multiple tissues—including fin (51, 55), heart (56, 57), and spinal cord (58)—its suppression could have impaired tissue regeneration (59).

This hypothesis is supported by studies showing that pharmacologically suppressing inflammation in zebrafish leads to impaired fin regeneration (60, 61). Similarly, macrophage ablation in adult zebrafish reduces regeneration and abolishes ray bifurcation (62).

Transcriptomic data confirmed that SKLT downregulates genes associated with early inflammatory responses in blastema tissue, including il20ra, a receptor for the pro-inflammatory and pro-osteoclastogenic cytokine in interleukin 20 (IL-20) (63, 64).

SKLT also affected adaptive immune signaling. Downregulation of genes involved in T-cell activation, interferon gamma (IFNγ) signaling, and antigen presentation was observed in fish treated with SKLT. The contribution of T-cells to estrogen deficiency-induced bone loss has been known since the early 2000s (65), when the field of osteoimmunology was formally proposed (66). Numerous T-cell–related genes were suppressed in SKLT-treated fish, including cd4, cd74a/b, ccl19a.1, and jak1, all of which are associated with osteoclastogenesis (30). Similarly, components of the interferon pathway (ifng1r, irf4b, irf10) (6769), and itk, encoding the IL2-inducible T-cell kinase, a regulator of natural killer T-cells (NKT) differentiation (70), were downregulated.

Notably, il7ra, linked to pro-osteoclastogenic cytokine production (71), and bone resorption–formation uncoupling (72), was also suppressed.

In contrast, il13ra2, encoding a receptor for the osteoprotective cytokine interleukin 13 (IL-13) (73, 74), was upregulated, suggesting a shift toward anti-inflammatory signaling.

Modulation of genes involved in osteoimmune paracrine communication was also observed. In medaka, Cxcl9l, produced by osteoblasts under osteoporotic conditions, directs a subset of macrophages toward osteoclast fate via Cxcr3.2 (75).

In zebrafish, cxcr3.1, one of the three Cxcr3 paralogs, was downregulated, possibly indicating reduced responsiveness of osteoclast progenitors to chemotactic signals. Moreover, the matrix metalloproteinase mmp13b, recently implicated in osteoclast maturation and resorptive function (76), was also downregulated in SKLT-exposed fish.

Together, these data indicate that SKLT induces broad suppression of immune-related molecular programs at 24 hpa, including those regulating T-cell activation, antigen presentation, and macrophage–osteoblast cross-talk. These immunological shifts correlate with reduced recruitment and maturation of osteoclast precursors in vivo.

Additionally, transcriptomic data indicate that SKLT reduced the expression of aldh1a2, the main enzyme responsible for the synthesis of retinoic acid (RA), a molecule whose expression is normally elevated in early stages of fin regeneration (7779). Given RA’s role in promoting blastema proliferation and growth, and in coordinating osteoblast dedifferentiation, proliferation, and redifferentiation during fin regeneration (7779), its suppression by SKLT likely contributes to the observed downregulation of cell cycle genes, and is consistent with our in vitro data showing an anti-proliferative effect on murine macrophages.

Subsequently, the potential of SKLT to prevent bone loss in a disease context was validated in a medaka model of RANKL-induced osteoporosis (19), which showed prevention of bone loss, consistent with the suppression of osteoclast activation found in zebrafish and by our transcriptomic analysis. These effects were reproduced in a murine macrophage line, where SKLT reduced the formation of multinucleated osteoclasts under RANKL stimulation, supporting its cross-species activity and translational potential (8082).

Taken together, our findings establish SKLT as a source of immunomodulatory and anti-osteoclastogenic activity, and based on the known phytochemical repertoire of Skeletonema costatum, including carotenoids such as zeaxanthin (83), long-chain polyunsaturated fatty acids like palmitoleic acid, and phenolic derivatives (83, 84). These classes of compounds emerge as plausible mediators of the observed anti-osteoclastogenic effects, as previous works have shown (8588), warranting future bioassay-guided fractionation to pinpoint their contribution to bone-protective activity.

Conclusions

The findings presented here highlight the therapeutic potential of targeting both innate and adaptive immune responses to treat bone erosive pathologies. We demonstrate that Skeletonema costatum contains ethanol-soluble immunomodulatory compounds capable of downregulating T-cell activation, antigen presentation, and macrophage differentiation toward the osteoclastic lineage. Although the precise identity of the active component(s) remains to be determined and represents a limitation of the present study, our data establish a foundation for future pharmacological exploration of this microalga. While further work is needed to isolate and characterize individual bioactives, the consistent immunomodulatory and anti-osteoclastogenic activity observed across zebrafish, medaka, and murine models supports the broader therapeutic promise of this ethanol-soluble fraction.

These results align with growing calls for holistic, nutraceutical-based approaches to managing chronic diseases involving immune dysregulation and bone loss.

Our study provides a proof of concept for the use of marine-derived immunomodulators to prevent bone erosion by regulating immune-driven osteoclastogenesis.

Materials and methods

Preparation of microalgal ethanolic fractions

The freeze-dried biomass of Skeletonema costatum (Necton S.A., Olhão, Portugal) was macerated with 96% ethanol (Laborspirit Lda, Lisbon, Portugal) at a biomass–solvent ratio of 1 g:40 mL (M/V) by gently stirring at 24 °C for 18 h. The macerate was centrifuged for 5 min at 1,000 × g via an Allegra 6R centrifuge (Beckman Coulter Inc., Brea, USA), and the supernatant was collected. The pellet was washed twice with 96% ethanol, and all the supernatants were pooled and then vacuum filtered sequentially through 0.45 µm and 0.22 µm nylon membranes (Labbox Labware S.L., Barcelona, Spain). The filtrate was concentrated with an RV 10 digital rotary evaporator (IKA-Werke GmbH & Co. KG, Staufen im Breisgau, Germany), with the temperature set at 40 °C and pressure at 178 mbar, until a dense, paste-like extract was obtained. The extraction yield (37.9 ± 2.4%) was calculated from 2 mL aliquots (n = 3) placed under a gentle flow of 99.8% nitrogen until complete evaporation of the solvent.

Fish maintenance

Zebrafish wild-type line AB and transgenic line Tg(Ola.ctsk:FRT-DsRed-FRT-Cre, myl7:EGFP)mh201,48, referred to as Tg(ctsk:DsRed) throughout the manuscript, were maintained in a water recirculating system ZebTEC (Tecniplast, Buguggiate, Italy) with the following conditions: temperature 28 ± 0.1 °C, pH 7.5 ± 0.1, conductivity 700 ± 50 μS, ammonia and nitrites at levels below 0.1 mg/L, nitrates lower than 50 mg/L, and a photoperiod of 14:10 h light-dark. The medaka transgenic line Tg(rankl:HSE: CFP)TG1135, hereafter referred to as Tg(rankl:HSE: CFP), was purchased from the National BioResource Project Medaka (NBRP Medaka) (89) and maintained in a water recirculating system with the following conditions: temperature 27 ± 0.1 °C, pH 7.0 ± 0.1, conductivity 300 ± 100 μS, ammonia and nitrites at levels below 0.1 mg/L, nitrates lower than 50 mg/L and a photoperiod of 14:10 h light–dark. For zebrafish and medaka, system water was prepared by supplementing reverse osmosis water with a salt mixture (Instant Ocean, City, USA) and sodium bicarbonate (Sigma–Aldrich, St. Louis, USA). All the fish were fed commercial dry food Zebrafeed (Sparos Lda, Olhão, Portugal) daily.

Zebrafish caudal fin regeneration assay

The caudal fin of wild-type or transgenic adult zebrafish aged 3–4 months was amputated 1–2 segments anterior to the bifurcation of the most peripheral branching lepidotrichia, as described by Cardeira et al. (2016) (45). For the first experiment, AB wild type zebrafish were sampled at a single time-point, 5 days post-amputation (dpa). After finectomy, fish (n = 14 for SKLT, and n = 32 for CTRL) were placed at 33 ± 1 °C in 3 L plastic containers at a density of 5 fish/L. The fish were exposed to the S. costatum ethanolic fraction (SKLT) at 56 µg/mL or ethanol (CTRL) at 0.1% supplemented in system water for 5 days, to evaluate their mineralogenic performance.

For the second experiments to track osteoclastic cells, transgenic Tg(ctsk:DsRed) fish (n = 18) were amputated and then exposed to 56 µg/mL (SKLT) or 0.1% ethanol (CTRL).

In all experiments, moderate water dynamics and air-water exchanges were facilitated by bubbling the fish tanks with an air pump. Water quality parameters were monitored daily and maintained stable for the duration of the experiment as follows: dissolved oxygen, 7.0 ± 0.5 mg/L; pH, 7.4 ± 0.2; and conductivity, 680 ± 20 μS. Wild-type fish were sacrificed at 120 hours post-amputation (hpa) via lethal anesthesia with 0.6 mM tricaine methanesulfonate (MS-222, pH 7.0; Sigma–Aldrich), immersed for 30 min in 0.03% alizarin red S (AR-S, pH 7.4; Sigma–Aldrich) and washed two times for 5 min with system water. The stained fish were imaged for fin morphometric analysis. Transgenic For the osteoclast-tracking experiment, At 24 hours post-amputation, half of the transgenic Tg(ctsk:DsRed) (n = 9) for each group were given a lethal anesthesia with MS-222, imaged for ctsk signal, then caudal fins were dissected and collected for TRAP staining. The other half (n = 9) were continuously kept in 56 µg/mL (SKLT) or 0.1% ethanol (CTRL), and fluorescence was at different time points (48, 72, 96, 120 and 240 hpa). Fish were immersed for 30 min in 0.2% calcein (pH 7.4; Fluorexon, Sigma–Aldrich), and washed twice in system water for 10 min. After staining, the fish were anesthetized for 5 min in tricaine and imaged.

For bulk the RNA sequencing (RNA-seq) analysis, fin blastemas were collected at 24 hpa, pooled (n = 5, 10 blastemas per pool), and stored at -80 °C until further processing.

TRAP activity in caudal fins

Caudal fins were amputated at the level of the caudal peduncle, washed once with 1X phosphate-buffered saline (PBS, pH 7.4) and fixed for 4 h in 4% paraformaldehyde solution (PFA, solubilized in PBS, pH 7.4) at 24 °C. Tartrate-resistant acid phosphatase (TRAP) staining was performed as previously described by Blum & Begemann (2015) (78), and the fins were imaged as described below.

Morphometric analysis of caudal fins

Fins were imaged under an MZ10F fluorescence stereomicroscope (Leica, Wetzlar, Germany) coupled to a DFC7000T color camera (Leica). Bright-field images were collected to assess the progression of fin regeneration. Fluorescence images were collected to assess (i) de novo bone formation in wild-type fish stained with AR-S or calcein and (ii) the involvement of ctsk-expressing cells in transgenic fish labeled with DsRed. Bright-field images were acquired with an exposure time of 1 ms. Fluorescence images were acquired with the filter set ET560/40x - ET630/75m and an exposure time of 600 ms for mCherry, and the filter set ET470/40x - ET525/50m and an exposure time of 80 ms for GFP. Other image parameters were as follows: gamma, 1.00; image format, 1920×1440 pixels; and binning, 1×1. The fluorescence images were analyzed using ImageJ software version 2.0.0-rc-69/1.52p and processed with the ZFBONE toolset for caudal fin morphometrics (90). Fin regeneration and mineralization were assessed following the method described by Cardeira et al. (2016) (45) by calculating the regenerated area (REG), the stump width (STU), the mineralized area (MIN), and the average width of the rays before the amputation (RAYs). Fin ray patterning was also assessed by calculating the average ray width ratio (W¯) and the average bifurcation ratio (Bt¯), as shown in Supplementary Equations a, b. For the quantification of the ctsk signal in transgenic fish, the DsRed-positive area was measured using a color threshold on fluorescence images and normalized with REG/STU. To quantify the TRAP signal, TRAP-positive areas were measured using a color threshold on bright-field images and subsequently normalized with REG/STU.

RNA preparation

Total RNA was extracted from pools of blastemas at 24 hpa (n = 5) using NZYol (NZYTech, Lisbon, Portugal) and quantified using a NanoDrop OneC spectrophotometer (Thermo Fisher Scientific, Waltham, USA). The RNA integrity was confirmed using an Experion Automated Electrophoresis system (Bio-Rad, Hercules, USA). Only RNA with an RNA integrity number (RIN) greater than 7 was used.

Bulk RNA sequencing and analysis of differentially expressed genes

Bulk RNA sequencing was outsourced to STABVIDA Lda (Caparica, Portugal). DNA libraries were constructed using a Stranded mRNA Library Preparation Kit (STABVIDA) and sequenced on a NovaSeq platform (Illumina, San Diego, USA) to generate 150 bp paired-end sequencing reads. The raw sequence data were processed using CLC Genomics Workbench 12.0.3. Trimming was performed in 3 steps: quality trimming on the basis of quality scores (error probability of 0.01), ambiguity trimming (ambiguous limit of 2 nucleotides), and length trimming to discard reads shorter than 30 nucleotides. The quality-checked sequencing reads were mapped against the zebrafish reference genome GRCz11 (GCF_000002035.6) using length fraction and similarity fraction equal to 0.8. The TPM (transcripts per million) and RPKM (reads per kilobase of transcript per million mapped reads) values were then determined from the mapped data. Differential expression analysis was performed with the multifactorial EdgeR method in R. Overall, the gene expression results were represented with a clustering heatmap of gene expression, with a principal component analysis (PCA) plotting the amount of variance explained by the two principal components and a volcano plot representing overall gene expression with a fold change on the x-axis and significance of expression on the y-axis (Supplementary Figure S2). These overall representations of gene expression were performed to assess variation patterns in the gene expression dataset and identify outlier samples for quality control. For functional annotation enrichment analysis, the lists of upregulated and downregulated genes were analyzed with the online resource Database for Annotation, Visualization and Integrated Discovery (DAVID) (91, 92) for Gene Ontology–Biological Processes (GO: BP) with FDR< 0.01 and enriched biological processes.

Medaka model of RANKL overexpression-induced osteoporosis

Eggs of the medaka line Tg(rankl:HSE: CFP) were produced following an in-house breeding program and maintained in Petri dishes with 40 mL of system water supplemented with 0.0002% (w/v) methylene blue until 8 days post-fertilization (dpf). At 8 dpf, the hatched larvae were incubated at 39 °C for 2 h (heat shock) and screened for CFP (cyan fluorescent protein) signal using a fluorescence microscopy (see parameters below). CFP-positive fish were distributed into a 6-well plate at a density of 5 larvae/well. Each well was filled with 10 mL of system water supplemented with SKLT at 56 µg/mL or 0.1% ethanol (vehicle). The treatment was renewed 100% daily. Six days after heat shock, the larvae (n = 17 for CTRL, and n = 12 for SKLT) were sacrificed with a lethal dose of tricaine (see above), stained with AR-S and imaged as described above. Bright-field and fluorescence images were collected for each fish. Bright-field images were acquired with an exposure time of 2 ms. AR-S fluorescence images were acquired using the filter set ET560/40x - ET630/75m and an exposure time of 700 ms. CFP fluorescence images were acquired using the filter set ET436/20x - ET480/40m and an exposure time of 200 ms. Other image parameters were as follows: gamma 1.00, image format 1920×1440 pixels, and binning 1×1. Images were analyzed using ImageJ using two ad hoc macros (Supplementary Files 1, 2). CFP fluorescence images were transformed into 8-bit images, and the CFP-positive area was measured using a color threshold (min intensity of 7 and max intensity of 255). The pixel mean intensity inside the CFP-positive area was calculated and used as a proxy for CPF fluorescence intensity. AR-S fluorescence images were transformed into 8-bit images, and AR-S positive areas were measured using a color threshold (min intensity 5 and max intensity 255). The AR-S positive area was normalized using the total body area (determined manually from bright-field images) to correct for inter-specimen size variation. The AR-S positive area in the abdominal and caudal vertebrae was used as a proxy for mineralization of the vertebral column. The nomenclature proposed by Di Biagio et al. (2022) (83) was used to identify vertebrae from 4-29. The number of mineralized neural arches was manually counted in each fish from the AR-S images. Following the collection of the data, the frequency distribution of the CFP intensity and tested data normality was determined through an Anderson–Darling test (p< 0.05). Normality was not met in the control group (p = 0.0035), but CFP intensity was normally distributed in SKLT-exposed fish (p = 0.220). The data distribution for CFP intensity was right skewed (median< mean) for the control group and bimodal for the SKLT group. These findings implied that clustering was necessary to proceed with the statistical analysis. The fish were therefore grouped according to CFP intensity into “High CFP” (mean pixel intensity > 50) and “Low CFP” (mean pixel intensity< 50) for further analysis.

Culture of mouse RAW 264.7 macrophages

RAW 264.7 cells were cultured in 10 cm cell culture dishes with Dulbecco’s Modified Eagle Medium (DMEM) supplemented with 10% fetal bovine serum (FBS; Sigma–Aldrich), 1% of penicillin–streptomycin (100x), 1% of 200 mM L-glutamine, and 0.2% of 250 µg/mL amphotericin B at 37 °C in a humidified atmosphere containing 5% CO2. PreconfluenT-cell cultures were sub-cultured at a ratio of 1:4 every other day using trypsin-EDTA solution (0.2% trypsin, 1.1 mM EDTA: prepared in 1x PBS, pH 7.4). All cell culture reagents were obtained from GIBCO-Thermo Fisher Scientific, unless otherwise stated.

Cytotoxicity and cell proliferation

An LDH cytotoxicity assay kit and a XTT-cell proliferation assay kit (Canvax Biotech, Córdoba, Spain) were used to evaluate the effects of SKLT on cellular toxicity and proliferation, respectively. For cytotoxicity, RAW 264.7 cells were seeded in a 96-well plate at 1.0 × 104 cells/well (n = 4) in 100 µL of culture medium supplemented with either the SKLT at 50, 75, 100, 200 µg/mL or 0.1% ethanol (vehicle) and cultured for 24 and 48 h. Relative cytotoxicity was calculated as a percentage of LDH lysis control. For cell proliferation, RAW 264.7 cells were seeded in a 96-well plate at 500 cells/well (n = 6) in 100 µL of culture medium supplemented with either the SKLT at 50, 75, 100, 200 µg/mL or vehicle and cultured for 1, 3, 5 or 7 days. The culture medium was renewed every other day. Relative cell proliferation was calculated as a percentage of the negative control.

Osteoclast differentiation

RAW 264.7 cells were seeded in 12-well plates at 1.5 × 104 cells/well in 1 mL of culture medium (as described above, n = 9) supplemented with either SKLT at 50 µg/ml or 0.1% ethanol (vehicle) and treated for 6 days with 50 ng/mL RANKL (Preprotech, London, UK) prepared in 0.1% bovine serum albumin (BSA) to induce their differentiation into multinucleated osteoclasts. The experimental conditions were as follows: undifferentiated control (0.1% ethanol); differentiated control (50 ng/mL RANKL + 0.1% ethanol); differentiated cells exposed to SKLT (50 ng/mL RANKL + 50 µg/mL SKLT in 0.1% ethanol). The culture medium was freshly prepared and replaced daily. Osteoclast differentiation was assessed by counting the number of tartrate-resistant acid phosphatase (TRAP)-positive cells and the number of nuclei following 4′,6-diamidino-2-phenylindole (DAPI) staining. For this purpose, the cells were fixed in 4% PFA at 24 °C (RT) for 10 min and stained according to the protocol of Blum & Begemann (2015) (78). The cells were subsequently stained with DAPI for nuclear detection and counting. The cells were imaged using an Axio Vert. A1 inverted microscope (ZEISS, Jena, Germany) coupled with (i) an Axiocam 202 monocolor camera (ZEISS) for DAPI images, or (ii) a VWR VisiCam 5 Plus (VWR, Radnor, USA) for bright-field TRAP images. The number of multinucleated osteoclasts per field (5.0X magnification; n = 9) was calculated considering only TRAP+ cells with at least two nuclei.

Statistical analysis

For all the experiments, normality was tested with a D’Agostino–Pearson omnibus normality test or with an Anderson–Darling test (p< 0.05). Homoscedasticity was tested through the Brown–Forsythe test (p< 0.05). When the distribution of the data of all the experimental groups resulted in normal and homogeneous, statistical differences between the control and the treated groups were tested with either an unpaired t test or one-way ANOVA, followed by Dunnett’s multiple comparison test (p< 0.05). If the distribution of the data of any of the experimental conditions resulted in non-normal or nonhomogeneous, significant differences between the control and the treatments were tested with a Mann–Whitney test or a nonparametric test followed by Dunn’s multiple comparison test (p< 0.05). Statistical analyses were performed using Prism version 9.00 (GraphPad Software Inc., La Jolla, United States).

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below:https://www.ebi.ac.uk/biostudies/, RNA sequencing data have been uploaded onto the free repository BioStudies EMBL-EBI with accession number E-MTAB-14677, and will be made publicly available upon article acceptance. All other raw data generated in this study will be made available upon request.

Ethics statement

Procedures involving animals were performed following the EU and Portuguese legislation for animal experimentation and welfare (Directives 86/609/CEE and 2010/63/EU; Portaria 1005/92, 466/95 and 1131/97; Decreto-Lei 113/2013). Animal handling and experimentation were performed by qualified operators accredited by the Portuguese Direção-Geral de Alimentação e Veterinária under authorization no. 012769/2021. All efforts were made to minimize pain, distress, and discomfort. The experiments were terminated (fish were returned to normal conditions or euthanized) whenever adverse effects were observed. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

AC: Funding acquisition, Conceptualization, Methodology, Writing – review & editing, Formal Analysis, Writing – original draft, Project administration, Validation, Visualization, Investigation, Data curation. KP: Formal Analysis, Writing – review & editing, Methodology, Investigation, Visualization. MT: Writing – review & editing, Software, Visualization. JR: Resources, Writing – review & editing, Visualization, Investigation. SP: Data curation, Methodology, Investigation, Writing – review & editing, Visualization, Formal Analysis. HP: Writing – review & editing, Resources. BL: Visualization, Software, Formal Analysis, Data curation, Methodology, Investigation, Writing – review & editing. MC: Funding acquisition, Writing – review & editing, Resources, Supervision, Project administration. VL: Writing – review & editing, Investigation, Funding acquisition, Supervision, Conceptualization, Visualization, Validation, Resources, Project administration. PG: Funding acquisition, Visualization, Resources, Validation, Conceptualization, Supervision, Project administration, Writing – review & editing, Investigation.

Funding

The author(s) declare financial support was received for the research and/or publication of this article. This work was financed by the European Maritime and Fisheries Fund (EMFF/FEAMP) through the National Operational Programme MAR2020 (grant 16-02-01-FMP-0057/OSTEOMAR), by the European Regional Development Fund (ERDF/FEDER) through the Transnational Cooperation Programme Atlantic Area (grant EAPA/151/2016/ BLUEHUMAN), by the Marie Skłodowska-Curie innovative training network BIOMEDAQU (grant H2020-MSCA-ITN/ 766347), This study received Portuguese national funds from FCT - Foundation for Science and Technology through contracts UID/04326/2025, UID/PRR/04326/2025 and LA/P/0101/2020 (DOI:10.54499/LA/P/0101/2020), and from the operational programmes CRESC Algarve 2020 and COMPETE 2020 through contract EMBRC.PT ALG-01-0145-FEDER-022121, and through the doctoral fellowship 2021.05406.BD, and by the operational programs CRESC Algarve 2020 and COMPETE 2020 through project EMBRC.PT ALG-01-0145-FEDER-02212.

Acknowledgments

The authors would also like to thank João Navalho and Necton S.A. (Olhão, Portugal) for kindly providing the microalgal biomass used to prepare the ethanolic extracts evaluated in this work.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2025.1655321/full#supplementary-material

References

1. Wu A, Bisignano C, James SL, Abady GG, Abedi A, Abu-Gharbieh E, et al. Global, regional, and national burden of bone fractures in 204 countries and territories, 1990–2019: a systematic analysis from the Global Burden of Disease Study 2019. Lancet Healthy Longevity. (2021) 2:e580–92. doi: 10.1016/S2666-7568(21)00172-0

PubMed Abstract | Crossref Full Text | Google Scholar

2. Williamson S, Landeiro F, McConnell T, Fulford-Smith L, Javaid MK, Judge A, et al. Costs of fragility hip fractures globally: a systematic review and meta-regression analysis. Osteoporos Int. (2017) 28:2791–800. doi: 10.1007/s00198-017-4153-6

PubMed Abstract | Crossref Full Text | Google Scholar

3. Feng X and McDonald JM. Disorders of bone remodeling. Annu Rev Pathol. (2011) 6:121–45. doi: 10.1146/annurev-pathol-011110-130203

PubMed Abstract | Crossref Full Text | Google Scholar

4. Mays S. Metabolic bone disease. In: Advances in Human Palaeopathology (2007). Chichester, UK: John Wiley & Sons, Inc. p. 215–51.

Google Scholar

5. Jilka RL. Molecular and cellular mechanisms of the anabolic effect of intermittent PTH. Bone. (2007) 40:1434–46. doi: 10.1016/j.bone.2007.03.017

PubMed Abstract | Crossref Full Text | Google Scholar

6. Haas AV and LeBoff MS. Osteoanabolic agents for osteoporosis. J Endocr Soc. (2018) 2:922–32. doi: 10.1210/js.2018-00121

PubMed Abstract | Crossref Full Text | Google Scholar

7. Minisola S, Cipriani C, Grotta GD, Colangelo L, Occhiuto M, Biondi P, et al. Update on the safety and efficacy of teriparatide in the treatment of osteoporosis. Ther Adv Musculoskelet Dis. (2019) 11. doi: 10.1177/1759720X19877994

PubMed Abstract | Crossref Full Text | Google Scholar

8. Seeman E and Martin TJ. Antiresorptive and anabolic agents in the prevention and reversal of bone fragility. Nat Rev Rheumatol. (2019) 15:225–36. doi: 10.1038/s41584-019-0172-3

PubMed Abstract | Crossref Full Text | Google Scholar

9. Kenkre JS and Bassett JHD. The bone remodelling cycle. Ann Clin Biochem. (2018) 55:308–27. doi: 10.1177/0004563218759371

PubMed Abstract | Crossref Full Text | Google Scholar

10. Dahiya N, Khadka A, Sharma AK, Gupta AK, Singh N, and Brashier DBS. Denosumab: A bone antiresorptive drug. Med J Armed Forces India. (2015) 71:71–5. doi: 10.1016/j.mjafi.2014.02.001

PubMed Abstract | Crossref Full Text | Google Scholar

11. Miller SA, St. Onge EL, and Whalen KL. Romosozumab: a novel agent in the treatment for postmenopausal osteoporosis. J Pharm Technol. (2021) 37:45–52. doi: 10.1177/8755122520967632

PubMed Abstract | Crossref Full Text | Google Scholar

12. Langdahl B. Treatment of postmenopausal osteoporosis with bone-forming and antiresorptive treatments: combined and sequential approaches. Bone. (2020) 139:115516. doi: 10.1016/j.bone.2020.115516

PubMed Abstract | Crossref Full Text | Google Scholar

13. Skjødt MK, Frost M, and Abrahamsen B. Side effects of drugs for osteoporosis and metastatic bone disease. Br J Clin Pharmacol. (2019) 85:1063–71. doi: 10.1111/bcp.13759

PubMed Abstract | Crossref Full Text | Google Scholar

14. Japanese Allied Committee on Osteonecrosis of the Jaw. Antiresorptive agent-related osteonecrosis of the jaw: Position Paper 2017. J Bone Miner Metab. (2017) 35:6–19. doi: 10.1007/s00774-016-0810-7

PubMed Abstract | Crossref Full Text | Google Scholar

15. Lorenzo J and Choi Y. Osteoimmunology. Immunol Rev. (2005) 208:5–6. doi: 10.1111/j.0105-2896.2005.00340.x

PubMed Abstract | Crossref Full Text | Google Scholar

16. Weitzmann MN. Bone and the immune system. In: Bone Toxicology (2017). Cham, Switzerland: Springer International. p. 363–98.

Google Scholar

17. Mori G, D'Amelio P, Faccio R, and Brunetti G. The interplay between the bone and the immune system. J Immunol Res. (2013) 2013:720504. doi: 10.1155/2013/720504

PubMed Abstract | Crossref Full Text | Google Scholar

18. Teitelbaum SL. Bone resorption by osteoclasts. Science. (2000) 289:1504–8. doi: 10.1126/science.289.5484.1504

PubMed Abstract | Crossref Full Text | Google Scholar

19. Phan QT, Liu R, Tan WH, Imangali N, Cheong B, Schartl M, et al. Macrophages switch to an osteo-modulatory profile upon RANKL induction in a Medaka (Oryzias latipes) osteoporosis model. JBMR Plus. (2020) 4:e10409. doi: 10.1002/jbm4.10409

PubMed Abstract | Crossref Full Text | Google Scholar

20. Lampiasi N, Russo R, and Zito F. The alternative faces of macrophage generate osteoclasts. BioMed Res Int. (2016) 2016:9089610. doi: 10.1155/2016/9089610

PubMed Abstract | Crossref Full Text | Google Scholar

21. Hofbauer LC, Lacey DL, Dunstan CR, Spelsberg TC, Riggs BL, and Khosla S. IL-1β and TNF-α, but not IL-6, stimulate osteoprotegerin ligand gene expression in human osteoblastic cells. Bone. (1999) 25:255–9. doi: 10.1016/S8756-3282(99)00162-3

PubMed Abstract | Crossref Full Text | Google Scholar

22. Zhang YH, Heulsmann A, Tondravi MM, Mukherjee A, and Abu-Amer Y. TNF-α stimulates RANKL-induced osteoclastogenesis via coupling of TNF type 1 receptor and RANK signaling pathways. J Biol Chem. (2001) 276:563–8. doi: 10.1074/jbc.M008198200

PubMed Abstract | Crossref Full Text | Google Scholar

23. Zha L, He L, Liang Y, Qin H, Yu B, Chang L, et al. TNF-α contributes to postmenopausal osteoporosis by synergistically promoting RANKL-induced osteoclast formation. BioMed Pharmacother. (2018) 102:369–74. doi: 10.1016/j.biopha.2018.03.080

PubMed Abstract | Crossref Full Text | Google Scholar

24. Fuller K, Murphy C, Kirstein B, Fox SW, and Chambers TJ. TNFα potently activates osteoclasts through a direct action independent of and strongly synergistic with RANKL. Endocrinology. (2002) 143:1108–18. doi: 10.1210/endo.143.3.8701

PubMed Abstract | Crossref Full Text | Google Scholar

25. Kanazawa K, Azuma Y, Nakano H, and Kudo A. TRAF5 functions in both RANKL- and TNFα-induced osteoclastogenesis. J Bone Miner Res. (2003) 18:443–50. doi: 10.1359/jbmr.2003.18.3.443

PubMed Abstract | Crossref Full Text | Google Scholar

26. Briot K, Geusens P, Em Bultink I, Lems WF, and Roux C. Inflammatory diseases and bone fragility. Osteoporos Int. (2017) 28:3301–14. doi: 10.1007/s00198-017-4189-7

PubMed Abstract | Crossref Full Text | Google Scholar

27. Weitzmann MN and Pacifici R. Estrogen deficiency and bone loss: an inflammatory tale. J Clin Invest. (2006) 116:1186–94. doi: 10.1172/JCI28550

PubMed Abstract | Crossref Full Text | Google Scholar

28. Cenci S, Weitzmann MN, Roggia C, Namba N, Novack D, Woodring J, et al. Estrogen deficiency induces bone loss by enhancing T-cell production of TNF-α. J Clin Invest. (2000) 106:1229–37. doi: 10.1172/JCI11066

PubMed Abstract | Crossref Full Text | Google Scholar

29. Lencel P and Magne D. Inflammaging: the driving force in osteoporosis? Med Hypotheses. (2011) 76:317–21. doi: 10.1016/j.mehy.2010.07.029

PubMed Abstract | Crossref Full Text | Google Scholar

30. Srivastava RK, Dar HY, and Mishra PK. Immunoporosis: immunology of osteoporosis—role of T-cells. Front Immunol. (2018) 9:657. doi: 10.3389/fimmu.2018.00657

PubMed Abstract | Crossref Full Text | Google Scholar

31. Weitzmann MN. T-cell activation induces human osteoclast formation via RANKL-dependent and-independent mechanisms. J Bone Miner Res. (2001) 16:328–37. doi: 10.1359/jbmr.2001.16.2.328

PubMed Abstract | Crossref Full Text | Google Scholar

32. Rifas L and Weitzmann MN. A novel T-cell cytokine, secreted osteoclastogenic factor of activated T-cells, induces osteoclast formation in a RANKL-independent manner. Arthritis Rheumatol. (2009) 60:3324–35. doi: 10.1002/art.24877

PubMed Abstract | Crossref Full Text | Google Scholar

33. Hu Y, Chen J, Hu G, Yu J, Zhu X, Lin Y, et al. Statistical research on the bioactivity of new marine natural products discovered during the 28 years from 1985 to 2012. Mar Drugs. (2015) 13:202–21. doi: 10.3390/md13010202

PubMed Abstract | Crossref Full Text | Google Scholar

34. Karthikeyan A, Joseph A, and Nair BG. Promising bioactive compounds from the marine environment and their potential effects on various diseases. J Genet Eng Biotechnol. (2022) 20:14. doi: 10.1186/s43141-021-00290-4

PubMed Abstract | Crossref Full Text | Google Scholar

35. Carletti A, Gavaia PJ, Cancela ML, and Laizé V. Metabolic bone disorders and the promise of marine osteoactive compounds. Cell Mol Life Sci. (2024) 81:11. doi: 10.1007/s00018-023-05033-x

PubMed Abstract | Crossref Full Text | Google Scholar

36. Sathasivam R, Radhakrishnan R, Hashem A, and Abd Allah EF. Microalgae metabolites: A rich source for food and medicine. Saudi J Biol Sci. (2019) 26:709–22. doi: 10.1016/j.sjbs.2017.11.003

PubMed Abstract | Crossref Full Text | Google Scholar

37. Jha D, Jain V, Sharma B, Kant A, and Garlapati VK. Microalgae-based pharmaceuticals and nutraceuticals: an emerging field with immense market potential. ChemBioEng Rev. (2017) 4:257–72. doi: 10.1002/cben.201600023

Crossref Full Text | Google Scholar

38. Saide A, Martínez KA, Ianora A, and Lauritano C. Unlocking the health potential of microalgae as sustainable sources of bioactive compounds. Int J Mol Sci. (2021) 22:4383. doi: 10.3390/ijms22094383

PubMed Abstract | Crossref Full Text | Google Scholar

39. Sproles AE, Fields FJ, Smalley TN, Le CH, Badary A, and Mayfield SP. Recent advancements in the genetic engineering of microalgae. Algal Res. (2021) 53:102158. doi: 10.1016/j.algal.2020.102158

Crossref Full Text | Google Scholar

40. Grama SB, Liu Z, and Li J. Emerging trends in genetic engineering of microalgae for commercial applications. Mar Drugs. (2022) 20:285. doi: 10.3390/md20050285

PubMed Abstract | Crossref Full Text | Google Scholar

41. Enzing C, Ploeg M, Barbosa M, and Sijtsma L. Microalgae-based products for the food and feed sector: an outlook for Europe. Inst Prospect Technol Stud. (2014), 1017–24.

Google Scholar

42. Bhattacharjya R, Marella TK, Tiwari A, Saxena A, Singh PK, and Mishra B. Bioprospecting of marine diatoms Thalassiosira, Skeletonema, and Chaetoceros for lipids and other value-added products. Mar Drugs. (2023) 21:32. doi: 10.3390/md21010032

PubMed Abstract | Crossref Full Text | Google Scholar

43. Silva M, Kamberovic F, Uota ST, Kovan IM, Viegas CS, Simes DC, et al. Microalgae as potential sources of bioactive compounds for functional foods and pharmaceuticals. Appl Sci. (2022) 12:5877. doi: 10.3390/app12125877

Crossref Full Text | Google Scholar

44. Carletti A, Rosa JT, Pes K, Borges I, Santos T, Barreira L, et al. The osteogenic and mineralogenic potential of the microalgae Skeletonema costatum and Tetraselmis striata CTP4 in fish models. Cell Mol Life Sci. (2023) 80:310. doi: 10.1007/s00018-023-04953-y

PubMed Abstract | Crossref Full Text | Google Scholar

45. Cardeira J, Gavaia PJ, Fernández I, Cengiz IF, Moreira-Silva J, Oliveira JM, et al. Quantitative assessment of the regenerative and mineralogenic performances of the zebrafish caudal fin. Sci Rep. (2016) 6:39191. doi: 10.1038/srep39191

PubMed Abstract | Crossref Full Text | Google Scholar

46. To TT, Witten PE, Renn J, Bhattacharya D, Huysseune A, and Winkler C. Rankl-induced osteoclastogenesis leads to loss of mineralization in a medaka osteoporosis model. Development. (2012) 139:141–50. doi: 10.1242/dev.071035

PubMed Abstract | Crossref Full Text | Google Scholar

47. Rosa JT, Tarasco M, Gavaia PJ, Cancela ML, and Laizé V. Screening of mineralogenic and osteogenic compounds in zebrafish—Tools to improve assay throughput and data accuracy. Pharmaceutics. (2022) 15:983. doi: 10.3390/ph15080983

PubMed Abstract | Crossref Full Text | Google Scholar

48. Cardeira-da-Silva J, Bensimon-Brito A, Tarasco M, Brandão AS, Rosa JT, Borbinha J, et al. Fin ray branching is defined by TRAP+ osteolytic tubules in zebrafish. Proc Natl Acad Sci U.S.A. (2022) 119. doi: 10.1073/pnas.2209231119

PubMed Abstract | Crossref Full Text | Google Scholar

49. Caetano-Lopes J, Henke K, Urso K, Duryea J, Charles JF, Warman ML, et al. Correction: Unique and non-redundant function of csf1r paralogues in regulation and evolution of post-embryonic development of the zebrafish. Development. (2020) 147:dev192211. doi: 10.1242/dev.192211

PubMed Abstract | Crossref Full Text | Google Scholar

50. Gao Y, Grassi F, Ryan MR, Terauchi M, Page K, Yang X, et al. IFN-γ stimulates osteoclast formation and bone loss in vivo via antigen-driven T-cell activation. J Clin Invest. (2007) 117:122–32. doi: 10.1172/JCI30074

PubMed Abstract | Crossref Full Text | Google Scholar

51. Sehring IM and Weidinger G. Recent advancements in understanding fin regeneration in zebrafish. Wiley Interdiscip Rev Dev Biol. (2020) 9:e367. doi: 10.1002/wdev.367

PubMed Abstract | Crossref Full Text | Google Scholar

52. Miron RJ, Bohner M, Zhang Y, and Bosshardt DD. Osteoinduction and osteoimmunology: Emerging concepts. Periodontology. (2024), 9–26. doi: 10.1111/prd.12519

PubMed Abstract | Crossref Full Text | Google Scholar

53. Wu D, Cline-Smith A, Shashkova E, Perla A, Katyal A, and Aurora R. T-cell mediated inflammation in postmenopausal osteoporosis. Front Immunol. (2021) 12:687551. doi: 10.3389/fimmu.2021.687551

PubMed Abstract | Crossref Full Text | Google Scholar

54. Laforest L, Brown CW, Poleo G, Géraudie J, Tada M, Ekker M, et al. Involvement of the sonic hedgehog, patched 1 and bmp2 genes in patterning of the zebrafish dermal fin rays. Development. (1998) 125:4175–84. doi: 10.1242/dev.125.21.4175

PubMed Abstract | Crossref Full Text | Google Scholar

55. Pfefferli C and Jaźwińska A. The art of fin regeneration in zebrafish. Regeneration. (2015) 2:72–83. doi: 10.1002/reg2.33

PubMed Abstract | Crossref Full Text | Google Scholar

56. González-Rosa JM, Burns CE, and Burns CG. Zebrafish heart regeneration: 15 years of discoveries. Regeneration. (2017) 4:105–23. doi: 10.1002/reg2.83

PubMed Abstract | Crossref Full Text | Google Scholar

57. Gao Y, Qian N, Xu J, and Wang Y. The roles of macrophages in heart regeneration and repair after injury. Front Cardiovasc Med. (2021) 8:744615. doi: 10.3389/fcvm.2021.744615

PubMed Abstract | Crossref Full Text | Google Scholar

58. Tsarouchas TM, Wehner D, Cavone L, Munir T, Keatinge M, Lambertus M, et al. Dynamic control of proinflammatory cytokines Il-1β and Tnf-α by macrophages in zebrafish spinal cord regeneration. Nat Commun. (2018) 9:4670. doi: 10.1038/s41467-018-07036-w

PubMed Abstract | Crossref Full Text | Google Scholar

59. Iribarne M. Inflammation induces zebrafish regeneration. Neural Regener Res. (2021) 16:1693–701. doi: 10.4103/1673-5374.306059

PubMed Abstract | Crossref Full Text | Google Scholar

60. Nguyen-Chi M, Laplace-Builhé B, Travnickova J, Luz-Crawford P, Tejedor G, Lutfalla G, et al. TNF signaling and macrophages govern fin regeneration in zebrafish larvae. Cell Death Dis. (2017) 8:e2979. doi: 10.1038/cddis.2017.374

PubMed Abstract | Crossref Full Text | Google Scholar

61. Hasegawa T, Hall CJ, Crosier PS, Abe G, Kawakami K, Kudo A, et al. Transient inflammatory response mediated by interleukin-1β is required for proper regeneration in zebrafish fin fold. Elife. (2017) 6:e22716. doi: 10.7554/eLife.22716.021

PubMed Abstract | Crossref Full Text | Google Scholar

62. Petrie TA, Strand NS, Tsung-Yang C, Rabinowitz JS, and Moon RT. Macrophages modulate adult zebrafish tail fin regeneration. Development. (2015) 142:406. doi: 10.1242/dev.120642

PubMed Abstract | Crossref Full Text | Google Scholar

63. Hsu YH, Chen WY, Chan CH, Wu CH, Sun ZJ, and Chang MS. Anti–IL-20 monoclonal antibody inhibits the differentiation of osteoclasts and protects against osteoporotic bone loss. J Exp Med. (2011) 208:1849–61. doi: 10.1084/jem.20102234

PubMed Abstract | Crossref Full Text | Google Scholar

64. Meng B, Yang B, Qu Y, Liu Y, Wu D, Fu C, et al. Dual role of interleukin-20 in different stages of osteoclast differentiation and its osteoimmune regulation during alveolar bone remodeling. Int J Mol Sci. (2023) 24:3810. doi: 10.3390/ijms24043810

PubMed Abstract | Crossref Full Text | Google Scholar

65. Garcia-Perez MA, Noguera I, Hermenegildo C, Martínez-Romero A, Tarin JJ, and Cano A. Alterations in the phenotype and function of immune cells in ovariectomy-induced osteopenic mice. Hum Reprod. (2006) 21:880–7. doi: 10.1093/humrep/dei413

PubMed Abstract | Crossref Full Text | Google Scholar

66. Arron JR and Choi Y. Bone versus immune system. Nature. (2000) 408:535–6. doi: 10.1038/35046196

PubMed Abstract | Crossref Full Text | Google Scholar

67. Hu CM, Jang SY, Fanzo JC, and Pernis AB. Modulation of T-cell cytokine production by interferon regulatory factor-4. J Biol Chem. (2002) 277:49238–46. doi: 10.1074/jbc.M205895200

PubMed Abstract | Crossref Full Text | Google Scholar

68. Wang S, He Q, Ma D, Xue Y, and Liu F. Irf4 regulates the choice between T lymphoid-primed progenitor and myeloid lineage fates during embryogenesis. Dev Cell. (2015) 34:621–31. doi: 10.1016/j.devcel.2015.07.011

PubMed Abstract | Crossref Full Text | Google Scholar

69. Biburger M and Tiegs G. α-Galactosylceramide-induced liver injury in mice is mediated by TNF-α but independent of Kupffer cells. J Immunol. (2005) 175:1540–50. doi: 10.4049/jimmunol.175.3.1540

PubMed Abstract | Crossref Full Text | Google Scholar

70. Gadue P and Stein PL. NK T-cell precursors exhibit differential cytokine regulation and require Itk for efficient maturation. J Immunol. (2002) 169:2397–406. doi: 10.4049/jimmunol.169.5.2397

PubMed Abstract | Crossref Full Text | Google Scholar

71. Toraldo G, Roggia C, Qian WP, Pacifici R, and Weitzmann MN. IL-7 induces bone loss in vivo by induction of receptor activator of nuclear factor κB ligand and tumor necrosis factor α from T-cells. Proc Natl Acad Sci U S A. (2003) 100:125–30. doi: 10.1073/pnas.0136772100

PubMed Abstract | Crossref Full Text | Google Scholar

72. Weitzmann MN, Roggia C, Toraldo G, Weitzmann L, and Pacifici R. Increased production of IL-7 uncouples bone formation from bone resorption during estrogen deficiency. J Clin Invest. (2002) 110:1643–50. doi: 10.1172/JCI0215687

PubMed Abstract | Crossref Full Text | Google Scholar

73. Raphael I, Nalawade S, Eagar TN, and Forsthuber TG. T-cell subsets and their signature cytokines in autoimmune and inflammatory diseases. Cytokine. (2015) 74:5–17. doi: 10.1016/j.cyto.2014.09.011

PubMed Abstract | Crossref Full Text | Google Scholar

74. Palmqvist P, Lundberg P, Persson E, Johansson A, Lundgren I, Lie A, et al. Inhibition of hormone and cytokine-stimulated osteoclastogenesis and bone resorption by interleukin-4 and interleukin-13 is associated with increased osteoprotegerin and decreased RANKL and RANK in a STAT6-dependent pathway. J Biol Chem. (2006) 281:2414–29. doi: 10.1074/jbc.M510160200

PubMed Abstract | Crossref Full Text | Google Scholar

75. Phan QT, Tan WH, Liu R, Sundaram S, Buettner A, Kneitz S, et al. Cxcl9l and Cxcr3.2 regulate recruitment of osteoclast progenitors to bone matrix in a medaka osteoporosis model. Proc Natl Acad Sci U.S.A. (2020) 117:19276–86. doi: 10.1073/pnas.2008422117

PubMed Abstract | Crossref Full Text | Google Scholar

76. Liu R, Imangali N, Ethiraj LP, Carney TJ, and Winkler C. Transcriptome profiling of osteoblasts in a Medaka (Oryzias latipes) osteoporosis model identifies Mmp13b as crucial for osteoclast activation. FronT-cell Dev Biol. (2022) 10:775512. doi: 10.3389/fcell.2022.775512

PubMed Abstract | Crossref Full Text | Google Scholar

77. Blum N and Begemann G. Retinoic acid signaling controls the formation, proliferation and survival of the blastema during adult zebrafish fin regeneration. Development. (2012) 139:107–16. doi: 10.1242/dev.065391

PubMed Abstract | Crossref Full Text | Google Scholar

78. Blum N and Begemann G. Osteoblast de- and redifferentiation are controlled by a dynamic response to retinoic acid during zebrafish fin regeneration. Development. (2015) 142:2894–903. doi: 10.1242/dev.120204

PubMed Abstract | Crossref Full Text | Google Scholar

79. Blum N and Begemann G. Retinoic acid signaling spatially restricts osteoblasts and controls ray-interray organization during zebrafish fin regeneration. Development. (2015) 142:2888–93. doi: 10.1242/dev.120212

PubMed Abstract | Crossref Full Text | Google Scholar

80. Pandey MK, Gupta SC, Karelia D, Gilhooley PJ, Shakibaei M, and Aggarwal BB. Dietary nutraceuticals as backbone for bone health. Biotechnol Adv. (2018) 36:1633–48. doi: 10.1016/j.biotechadv.2018.03.014

PubMed Abstract | Crossref Full Text | Google Scholar

81. Raut N, Wicks SM, Lawal TO, and Mahady GB. Epigenetic regulation of bone remodeling by natural compounds. Pharmacol Res. (2019) 147:104350. doi: 10.1016/j.phrs.2019.104350

PubMed Abstract | Crossref Full Text | Google Scholar

82. Roseti L, Borciani G, Grassi F, Desando G, Gambari L, and Grigolo B. Nutraceutical benefits for bone health: Role in osteoporosis prevention. Basel, Switzerland: MDPI AG (2023).

Google Scholar

83. Di Biagio C, Dellacqua Z, Martini A, Huysseune A, Scardi M, Witten PE, et al. A baseline for skeletal investigations in Medaka (Oryzias latipes): The effects of rearing density on the postcranial phenotype. Front Endocrinology. (2022) 13:893699. doi: 10.3389/fendo.2022.893699

PubMed Abstract | Crossref Full Text | Google Scholar

84. Ashour M, Ringø E, Mabrouk MM, Mansour AI, Naiel MA, Mansour AT, et al. Marine Diatom Skeletonema costatum Dietary Supplementation Improves Growth, Immunological Responses, Antioxidant Activities, Gene Expressions, and Digestive Enzymes of Shrimp Litopenaeus vannamei. Aquaculture Nutr. (2025), 5013608. doi: 10.1155/anu/5013608

PubMed Abstract | Crossref Full Text | Google Scholar

85. Bastos CR, Maia IB, Pereira H, and Navalho J. & Varela JC: Optimisation of biomass production and nutritional value of two marine diatoms (Bacillariophyceae), Skeletonema costatum and Chaetoceros calcitrans. Biology. (2022) 11:594. doi: 10.3390/biology11040594

PubMed Abstract | Crossref Full Text | Google Scholar

86. Poudel S, Martins G, Cancela ML, and Gavaia PJ. Resveratrol-mediated reversal of doxorubicin-induced osteoclast differentiation. Int J Mol Sci. (2022) 23:15160. doi: 10.3390/ijms232315160

PubMed Abstract | Crossref Full Text | Google Scholar

87. Odes-Barth S, Khanin M, Linnewiel-Hermoni K, Miller Y, Abramov K, Levy J, et al. Inhibition of osteoclast differentiation by carotenoid derivatives through inhibition of the NF-κB pathway. Antioxidants. (2020) 9:1167. doi: 10.3390/antiox9111167

PubMed Abstract | Crossref Full Text | Google Scholar

88. Van Heerden B, Kasonga A, Kruger MC, and Coetzee M. Palmitoleic acid inhibits RANKL-induced osteoclastogenesis and bone resorption by suppressing NF-κB and MAPK signalling pathways. Nutrients. (2017) 9:441. doi: 10.3390/nu9050441

PubMed Abstract | Crossref Full Text | Google Scholar

89. National BioResource Project Medaka. (2022). Available online at: https://shigen.nig.ac.jp/medaka/.

Google Scholar

90. Tarasco M, Cordelières FP, Cancela ML, and Laizé V. ZFBONE: An ImageJ toolset for semi-automatic analysis of zebrafish bone structures. Bone. (2020) 138:115480. doi: 10.1016/j.bone.2020.115480

PubMed Abstract | Crossref Full Text | Google Scholar

91. Huang DW, Sherman BT, and Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. (2009) 4(1):44-57. doi: 10.1038/nprot.2008.211

PubMed Abstract | Crossref Full Text | Google Scholar

92. Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. (2022) 50(W1):W216-W221. doi: 10.1093/nar/gkac194

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: osteoimmunology, osteoporosis, T-cells, antigen presentation, macrophages, osteoclast differentiation, microalgae, Skeletonema costatum

Citation: Carletti A, Pes K, Tarasco M, Rosa JT, Poudel S, Pereira HG, Louro B, Cancela ML, Laizé V and Gavaia PJ (2025) Immunomodulatory inhibition of osteoclastogenesis by a marine microalgal ethanol fraction targeting T-cells, antigen presentation, and macrophage fate. Front. Immunol. 16:1655321. doi: 10.3389/fimmu.2025.1655321

Received: 27 June 2025; Accepted: 25 September 2025;
Published: 10 October 2025.

Edited by:

Kosuke Nishi, Ehime University, Japan

Reviewed by:

Jen-Tsung Chen, National University of Kaohsiung, Taiwan
Chalermpong Saenjum, Chiang Mai University, Thailand

Copyright © 2025 Carletti, Pes, Tarasco, Rosa, Poudel, Pereira, Louro, Cancela, Laizé and Gavaia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Paulo J. Gavaia, cGdhdmFpYUB1YWxnLnB0

Present addresses: Alessio Carletti, Max Plank Institute for Hearth and Lung Research – Department of Developmental Genetics, Bad Nauheim, Germany
Marco Tarasco, Max Plank Institute for Hearth and Lung Research – Department of Developmental Genetics, Bad Nauheim, Germany
Sunil Poudel, Center for Transplantation Sciences, Massachusetts General Hospital, Harvard Medical School, Boston, MA, United States
Katia Pes, Excellence Cluster Cardio-Pulmonary Institute CPI, Justus-Liebig University of Giessen, Giessen, Germany

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.