Abstract
Aim:
Osteoarthritis (OA) is increasingly understood as a condition influenced by the diversity of chondrocytes and immune system dysfunction; however, the specific molecular elements connecting distinct chondrocyte populations to the immune environment remain inadequately explored. The protein CORO1A, associated with the cytoskeleton and implicated in immune regulation, has not been thoroughly investigated in the context of OA or through single-cell analysis.
Approach:
We employed multi-omics techniques alongside single-cell RNA sequencing, supported by experimental validation using human OA cartilage samples, a murine model of medial meniscus destabilization (DMM), and in vitro inflammatory chondrocyte models to elucidate the expression patterns and functional significance of CORO1A in OA.
Findings:
The CORO1A gene is recognized as significantly expressed in osteoarthritic cartilage. Transcriptomic evaluations reveal that CORO1A is predominantly expressed in inflammatory and proliferative chondrocyte subpopulations, while showing lower levels in other chondrocyte types, underscoring the specificity of these groups. The upregulation of CORO1A is associated with increased infiltration of various immune cells, including M2 macrophages, plasma cells, and natural killer cells, and is linked to the activation of IL-6/JAK-STAT3 and TNF-α/NF-κB signaling pathways. Additionally, experimental data indicate that CORO1A is markedly elevated in the cartilage regions affected by osteoarthritis in both humans and mice. In terms of its function, the expression of CORO1A in primary chondrocytes is enhanced by inflammatory stimuli, while its suppression leads to reduced MMP13 expression and diminished chondrocyte migration. Conclusion: This study identifies CORO1A as a regulator specific to subpopulations within the immune microenvironment, connecting immune remodeling to both inflammatory and proliferative chondrocytes in OA. Considering the diversity of chondrocytes, the presence of immune cells, and catabolic processes, our results emphasize the significance of CORO1A as a mechanistic factor contributing to cartilage degradation, suggesting its potential as a therapeutic target for OA.
1 Introduction
Osteoarthritis (OA) is a long-lasting and degenerative joint condition marked by the gradual deterioration of articular cartilage (1). This ailment primarily affects the body’s weight-bearing joints, such as the knees and hips. As the condition advances, individuals experience significant pain, functional limitations, and reduced mobility, placing a considerable strain on global healthcare systems (2, 3). Knee osteoarthritis (KOA), the most common form of OA, is particularly alarming due to its status as a leading cause of disability worldwide (4). Despite its high prevalence and considerable socioeconomic impact, existing treatment methods largely offer only symptom relief. The main aim of these approaches is to alleviate symptoms and slow disease progression rather than address the underlying structural or functional joint issues (5–7). There is a significant gap in understanding the molecular mechanisms that contribute to the onset and advancement of osteoarthritis. Therefore, it is crucial to investigate the underlying factors involved, which could pave the way for the development of innovative and effective treatment options that go beyond mere symptom management.
The progression of osteoarthritis (OA) is influenced not only by cartilage deterioration but also by complex interactions among various chondrocyte populations and the immune microenvironment (8). Recent studies indicate that communication between immune cells—such as macrophages, T cells, and natural killer (NK) cells—and cartilage plays a significant role in both cartilage damage and repair, rather than merely participating in inflammatory responses (9, 10). Recent studies indicate that communication between immune cells—such as macrophages, T cells, and natural killer (NK) cells—and cartilage plays a significant role in both cartilage damage and repair, rather than merely participating in inflammatory responses. These immune cells are crucial in regulating the remodeling of the extracellular matrix, which involves matrix metalloproteinases (MMPs); their activity must be carefully controlled to maintain a balance between cartilage degradation and regeneration (11). Importantly, M2 macrophages release anti-inflammatory cytokines and growth factors like IGF1 and PDGF, which are also activated during the adaptive immune response in the context of OA (12, 13). Nevertheless, despite growing awareness of chondrocyte diversity and the immune microenvironment’s role in OA, the specific regulatory molecules that govern immune cell polarization related to cartilage repair have yet to be identified.
CORO1A is a crucial protein within the cytoskeleton and belongs to the Coronin protein family, which plays significant roles in various functions, including the movement of immune cells, development of the cytoskeleton, and signaling pathways. These processes are vital for maintaining immune readiness and initiating appropriate inflammatory responses, as highlighted in recent studies (14, 15). In addition to its fundamental cellular roles, CORO1A is increasingly recognized as a pivotal regulator of inflammation and tissue remodeling driven by immune activity. A notable aspect of CORO1A’s function is its dependence on β2 integrin for the adhesion, spreading, and migration of neutrophils under flow conditions, underscoring its essential role in the early stages of inflammatory responses (16).
Additionally, CORO1A has emerged as a key immune-related regulator in inflammatory joint disorders, with gene ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways associated with CORO1A showing significant enrichment in inflammation and immune signaling. This highlights its potential involvement in the dynamic remodeling of the pathological immune microenvironment (17). Beyond its inflammatory role, CORO1A has been recognized as a modulator of bone resorption, influencing lysosomal secretion and protease release in osteoclasts. This suggests a crucial link between immune cell activity and the mechanisms of bone remodeling (18). Furthermore, alterations in CORO1A expression have been observed in the early stages of osteoarthritis (OA) within both subchondral bone and synovial fluid. These changes are closely associated with inflammatory activity and disease progression, reinforcing the significance of CORO1A in the pathophysiology of OA (19–21).
A key area that warrants further investigation is the role of CORO1A in regulating immune cell activity within the osteoarthritis (OA) immune microenvironment, as well as its specific contributions to cartilage repair mechanisms and disease progression. To comprehensively explore this topic, we undertake a thorough examination of CORO1A’s involvement in OA, combining bioinformatics approaches with experimental validation. Initially, we pinpoint immune-related genes that show differential expression in OA, elucidating their relationship with CORO1A and significant immune-regulatory pathways through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Subsequently, we identify the immune cell populations most closely associated with CORO1A expression via immune infiltration analyses, with particular emphasis on macrophage polarization, the activity of plasma cells, and the role of natural killer (NK) cells in relation to this protein. Ultimately, we investigate the potential functional roles of CORO1A in altering the immune microenvironment of OA and its implications for cartilage repair processes. Collectively, our findings position CORO1A as a previously overlooked immunoregulatory element in the OA context, offering new insights into its mechanisms and paving the way for potential clinical strategies to modulate immune responses.
2 Results
2.1 Osteoarthritic cartilage exhibits a distinct and coordinated transcriptional signature
To investigate the relationship between osteoarthritis (OA) and the overall changes in gene expression, rather than focusing solely on individual transcriptional variations, we analyzed the cartilage transcriptomes from the GSE286154 dataset, comparing OA samples with control samples. Our differential expression analysis revealed 778 genes with significant dysregulation, comprising 557 that were upregulated and 221 that were downregulated. This pattern indicates a broad activation of transcription in cartilage affected by OA (Figure 1A). Additionally, the relative expression levels of the top 50 upregulated and 50 downregulated DEGs across the OA and control samples were visualized using a heatmap (Figure 1B). These findings underscore the unique and reproducible molecular characteristics of OA cartilage. Collectively, these results imply that OA is marked by a consistent reorganization of transcription, laying a strong foundation for further functional annotation and network-level studies aimed at identifying key regulators associated with the condition.
Figure 1
2.2 Transcriptional changes linked to OA focus on immune control and microenvironmental
Interaction Pathways To explore the biological themes that elucidate the consistent transcriptional alterations seen in osteoarthritic cartilage, we conducted an enrichment analysis of the differentially expressed genes (DEGs) utilizing Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG). Rather than being scattered across various unrelated cellular functions, these genes exhibited a pronounced tendency towards activities associated with immune responses and interactions within the microenvironment, such as chemotaxis, taxis, cell-cell adhesion, promotion of cell adhesion, and regulation of leukocyte growth (Figure 1C). At the level of cellular components, the enriched terms were primarily linked to the outer regions of the plasma membrane, collagen-rich extracellular matrix, and the membranes of secretory granules, underscoring their significant role in communication between cells and between cells and the matrix. Additionally, the terms related to molecular functions encompassed immune receptor activity, integrin binding, extracellular matrix interaction, and protein-glycan binding, thereby emphasizing the concept of modified immune-matrix interactions in cartilage affected by OA.
The results of the KEGG pathway analysis reveal a notable enrichment in signaling pathways crucial for the regulation of the immune system and the microenvironment. This encompasses pathways such as the PI3K-Akt signaling pathway, interactions between cytokines and receptors, as well as cell adhesion molecules (CAMs) (Figure 1D). Overall, these enrichment findings suggest that immune regulation, the relay of inflammatory signals, and interactions with the extracellular matrix are fundamental biological processes that initiate transcriptional changes associated with OA. This underscores the importance of a network-based approach to pinpoint essential regulatory genes that govern these mechanisms.
2.3 Analysis of weighted gene co-expression networks reveals a gene module linked to OA
In order to investigate the coordinated gene expression patterns relevant to osteoarthritis, we developed a weighted gene co-expression network utilizing the complete expression dataset from GSE286154. The initial evaluation involved clustering the samples to assess data quality and identify any potential outliers (Figure 2A). This analysis revealed that a topology could be established, achieving an approximate scale-free property with R2>0.85 at a soft-thresholding power of β=11, which facilitated the formation of a biologically meaningful network with sufficient gene connectivity (Figure 2B). We performed gene clustering based on topological overlap to create hierarchical co-expression modules, each represented by distinct colors (Figure 2C). To pinpoint the modules relevant to the disease, we analyzed the correlations between module eigengenes and the osteoarthritis phenotype. Among the modules examined, the turquoise module exhibited the strongest positive correlation with osteoarthritis (Figures 2D, E), suggesting that the genes within this module may be crucial in the molecular processes of OA. As a result, this module was chosen for further investigation of hub genes and their functional roles.
Figure 2
2.4 Identification of CORO1A as a potential gene linked to OA via network-centric hub gene evaluation
Recognizing that transcriptional alterations linked to osteoarthritis (OA) are associated with immune response pathways and interactions within the microenvironment, we aimed to pinpoint the crucial regulators that might orchestrate these processes at the network level. By integrating intramodular connectivity (kWithin) with gene significance (GS), we identified the top 100 hub genes and cross-referenced them with previously identified differentially expressed genes (DEGs), resulting in 97 common candidates (Figure 2F). This method of overlap ensured that our subsequent analyses concentrated on genes that were not only significant within the co-expression network but also exhibited transcriptional dysregulation during OA.
To clarify the functional relationships among the selected candidates, we developed a protein-protein interaction (PPI) network utilizing the STRING database. Analyzing the network’s topology showed a high degree of interconnectivity at its core, indicating the presence of essential regulatory nodes. By employing the CytoHubba plug-in, we ranked the genes based on two criteria: betweenness centrality and stress centrality, which assess the gene’s capacity to influence information flow within the network. The overlap of the top 40 genes identified through both metrics resulted in a collection of 39 reliable hub candidates (Figure 2H).
Among the potential hub candidates, CORO1A stood out due to its significant concentration of various analytical metrics, such as its position within a module, its connections to other nodes, and its association with diseases. Among the potential hub candidates, CORO1A stood out due to its significant concentration of various analytical metrics, such as its position within a module, its connections to other nodes, and its association with diseases. Notably, at the pathway level, these genes showed enrichment in immune and inflammatory pathways, encompassing phagosome, cell adhesion molecules (CAMs), systemic lupus erythematosus, rheumatoid arthritis, and graft-versus-host disease (Figure 3A). Furthermore, a functional enrichment analysis revealed that CORO1A and its co-expressed neighbors were significantly associated with immune-related biological processes, including leukocyte proliferation, modulation of immune effector processes, leukocyte-mediated immunity, monocyte proliferation, and phagocytosis (Figure 3B).
Figure 3
In summary, the findings from the network analysis highlight CORO1A as a crucial hub linking immune activation, inflammatory signaling, and cell interaction pathways that play a significant role in osteoarthritis progression. Due to its prominent centrality within the network, consistent differential expression, and potential role in immune-related joint disorders, CORO1A has been identified as the primary candidate gene for further experimental investigation.
2.5 CORO1A high expression correlates with pro-inflammatory and pro-fibrotic signaling programs
In order to explore the biological processes associated with elevated CORO1A expression, we performed a comparative analysis of global pathway activities in osteoarthritis (OA) samples exhibiting high versus low levels of CORO1A through Gene Set Enrichment Analysis (GSEA). Rather than identifying individual pathway alterations, our findings revealed that increased CORO1A expression was consistently associated with a synchronized upregulation of inflammatory responses, stress-related mechanisms, and tissue remodeling activities.
Among the gene sets that showed the highest enrichment in samples exhibiting elevated CORO1A levels, notable pathways included IL6-JAK-STAT3 signaling (NES = 2.83, FDR < 0.001), Inflammatory Response (NES = 2.49, FDR < 0.001), and TNF-alpha signaling through NF-κB (NES = 1.68, FDR < 0.001) (Figure 3C). Additionally, pathways associated with structural changes and fibrotic activation were also prominently enriched, such as Epithelial-Mesenchymal Transition (NES = 2.24, FDR < 0.001) and Apoptosis (NES = 1.97, FDR = 0.001). The simultaneous activation of inflammatory and mesenchymal-like pathways indicates that elevated CORO1A expression is associated not only with immune responses but also with catabolic and pro-fibrotic cellular conditions linked to cartilage damage and joint remodeling.
In summary, the findings across the entire system suggest that the expression of CORO1A is useful for categorizing OA samples within a distinct spectrum characterized by pro-inflammatory and pro-remodeling traits, providing a mechanistic understanding of its link to immune infiltration and the severe disease characteristics noted in subsequent evaluations.
2.6 Single-cell transcriptomics reveals preferential enrichment of CORO1A in inflammatory and proliferative chondrocyte subpopulations
In order to investigate whether the increased expression of CORO1A noted in large-scale transcriptomic studies of osteoarthritic cartilage indicates specific regulatory mechanisms at the single-cell level, we analyzed single-cell RNA sequencing (scRNA-seq) data from both normal and osteoarthritic cartilage samples (GSE169454). Following rigorous quality control measures, we identified 52,374 high-quality cells for subsequent analysis.
Through the application of unsupervised clustering techniques, we have successfully identified nine major groups of transcriptionally unique cell populations within the dataset, with chondrocytes emerging as the predominant cell type (Figures 4B–D). In the chondrocyte compartment, we recognized eight distinct subpopulations based on the expression of canonical marker genes and automated annotation methods. These subpopulations include inflammatory chondrocytes (InfC), proliferating chondrocytes (ProC), hypertrophic-like chondrocytes (HTC), homeostatic chondrocytes (HomC), precursor chondrocytes (PreC), fibrocartilage chondrocytes (FC), along with several other minor groups (Figures 4E, F).
Figure 4
Crucially, the expression of CORO1A demonstrated a notable selectivity for specific cell types, with a significant lack of enrichment in the chondrocyte population (Figure 4E). Upon examining various chondrocyte subtypes, it became evident that the expression levels of CORO1A were markedly elevated in these specific subtypes, which are known to contribute to inflammation associated with osteoarthritis (OA) and impaired cartilage remodeling. Additionally, our findings revealed that CORO1A was expressed at much higher levels in OA samples compared to control samples, both in the overall chondrocyte population and within the subset of cells that express CORO1A (Figures 4H, I).
Investigation on individual cells has revealed that the expression of CORO1A in osteoarthritis varies among chondrocytes, being particularly elevated in certain inflammatory and proliferative subgroups. This indicates that CORO1A plays a role in activating chondrocytes that contribute to the progression of OA.
2.7 CORO1A is elevated in degenerated cartilage found in human OA and experimental models
To confirm the cellular findings from single-cell analysis within the tissue context, we assessed the natural expression of CORO1A in the cartilage of individuals with osteoarthritis and in a mouse model of the disease. Cartilage samples were obtained from patients undergoing total knee replacement surgery, with areas of intact and damaged cartilage identified based on anatomical location and visual characteristics (Figure 5A). Western blot analysis of these samples revealed a marked increase in CORO1A protein levels in damaged cartilage (Figures 5B, C). Consistently, immunohistochemical staining demonstrated a notable increase in the number of CORO1A-positive chondrocytes in the damaged regions compared to less affected areas (Figure 5D).
Figure 5
To validate these findings in a controlled setting, we employed the destabilization of the medial meniscus (DMM) mouse model. As expected, the DMM procedure led to significant cartilage deterioration and a reduction in proteoglycans, as evidenced by elevated OARSI scores (Figure 5G). Importantly, both conditions of cartilage affected by DMM exhibited a considerable increase in the proportion of CORO1A-positive chondrocytes compared to control samples, indicating that the presence of CORO1A-positive chondrocytes in degenerated cartilage is a notable aspect of osteoarthritis pathology (Figures 5H, I). In summary, the results indicate that the upregulation of CORO1A is a consistent observation in both human osteoarthritis and the mouse model.
2.8 CORO1A expression relates to a pro-inflammatory immune microenvironment in OA
Given the notable connections between CORO1A levels and inflammatory states in chondrocytes, we proceeded to investigate if CORO1A levels are linked to overall alterations in the immune microenvironment within osteoarthritic joints. We analyzed the proportions of 22 immune cell subtypes in osteoarthritis samples through CIBERSORT deconvolution applied to bulk transcriptomic data (see Figures 6A, B). The immune structure of OA cartilage exhibited considerable diversity, predominantly featuring macrophages, T cells, and B cells.
Figure 6
When analyzing lineage cells, significant differences in immune cell composition emerged based on CORO1A expression levels. Notably, samples with elevated CORO1A expression exhibited a marked rise in the infiltration of M2 macrophages, plasma cells, and activated NK cells, while showing a decrease in naive B cells and inactive mast cells (Figures 6C, D). These immune cell subsets are crucial for chronic inflammation, antibody-driven immune responses, and tissue remodeling in osteoarthritis. This finding reinforces the connection between CORO1A and an enhanced immune environment that favors osteoarthritis, highlighting its role as a molecular link between chondrocyte dysfunction and immune-related pathology.
2.9 CORO1A enhances catabolic behavior and chondrocyte mobility in inflammatory settings
To investigate the function of CORO1A in chondrocyte behavior, we employed an inflammatory model using IL-1β in primary mouse chondrocyte cultures. Immunofluorescence analysis revealed a notable increase in CORO1A protein levels following IL-1β treatment (Figures 7A, B). Notably, the reduction of CORO1A resulted in a significant decline in MMP13 expression induced by IL-1β (Figure 7D). To confirm CORO1A’s active involvement in these processes, we conducted a knockdown using siRNA. Among the three siRNAs evaluated, si-673 proved to be the most effective in diminishing CORO1A expression and was chosen for subsequent experiments (Figure 7E).
Figure 7
The involvement of CORO1A in the catabolic processes of chondrocytes during inflammation was investigated. Given that CORO1A regulates the cytoskeletal structure, we assessed its impact on chondrocyte movement. Our wound-healing assays revealed a notable reduction in the migration ability of chondrocytes with CORO1A knockdown at both the 24 and 48-hour marks post-scratch injury (Figures 7F, G).
In summary, our findings suggest that CORO1A acts as a mediator of catabolism and aids in the movement of chondrocytes during inflammation, thereby reinforcing the theory that CORO1A promotes harmful cellular processes associated with the development of osteoarthritis.
3 Discussion
Osteoarthritis (OA) is acknowledged as a widespread joint disorder characterized by cartilage degradation, synovial inflammation, and subchondral bone alterations. To comprehend the reasons behind the inability of cartilage to sustain homeostasis as OA progresses, it is important to identify factors linking inflammatory signaling to changes in chondrocyte behavior (22, 23).
Our bioinformatics analysis identified immune-regulatory, cytokine-receptor, PI3K-Akt, and cell adhesion-related pathways as major features of osteoarthritic cartilage. These enrichment results suggest that OA-associated transcriptional changes are not driven solely by classical inflammatory cytokines such as IL-1β and TNF-α, but are also linked to cell survival, adhesion, extracellular matrix interactions, and potentially motility-related processes. Consistently, clinical trials of biologic anti-inflammatory agents in hand OA, including TNF inhibitors, have shown limited symptomatic efficacy, suggesting that OA-related inflammatory-catabolic responses cannot be fully explained by canonical cytokine-driven inflammation alone (24). In this context, CORO1A-related changes may represent part of a broader network linking inflammatory signaling, adhesion/cytoskeletal remodeling, and chondrocyte state changes, rather than merely a common inflammatory signature.
Our focus on CORO1A stems from the pivotal role that coronin proteins play in linking actin dynamics and signal transduction (14, 25). Although CORO1A has been extensively examined in immune contexts, particularly regarding its capacity to influence receptor signaling and cytoskeletal organization (26), its significance in the context of osteoarthritis (OA) cartilage remains inadequately defined. Our findings suggest that CORO1A may also be relevant to cartilage pathology beyond its established role in immune-related cytoskeletal regulation.
In both human osteoarthritic cartilage and the DMM model, a consistent increase in CORO1A levels was noted, indicating an association with a degenerative-inflammatory cartilage state. Notably, in certain DMM cases, a rise in CORO1A staining was observed alongside increased cartilage damage (see Supplementary Figure 2). While these instances were insufficient for subgroup analysis, these observations are consistent with an association between CORO1A expression and a more inflammatory or degenerative cartilage state. Supporting this notion, CORO1A has been detected in notably higher amounts in more systemically inflamed arthritic conditions, such as the synovial fluid of RA (20), implying that CORO1A might reflect an inflammatory response rather than being specific to a particular disease (12). This differentiation is crucial for its translational significance: CORO1A could serve as a more valuable indicator of the interactions between inflammatory and catabolic processes and cellular states, rather than merely as a diagnostic marker for a single condition.
Our scRNA-seq analysis further refined this observation by showing that CORO1A was not uniformly elevated across all chondrocytes. Instead, its up-regulation is specifically observed in two disease-relevant states: inflammatory chondrocytes (InfC) and proliferative chondrocytes (ProC). This pattern suggests that CORO1A is associated with specific disease-related chondrocyte states rather than generalized cartilage degeneration (27).
The presence of CORO1A in both InfC and ProC raises the possibility that persistent inflammatory signals may engage a cytoskeleton-linked signaling program associated with both catabolic gene expression and cell motility. As a member of the coronin family and an actin-binding signaling molecule, CORO1A may act at the interface of membrane dynamics and signal transduction, providing a plausible link to cytokine responsiveness and cell motility (25, 28).
In vitro, siRNA-mediated CORO1A knockdown attenuated IL-1β-induced MMP13 expression and reduced scratch closure, supporting a possible association between CORO1A and inflammatory-catabolic as well as motility-related responses in chondrocytes. Mechanistically, this may involve small GTPase signaling, as CORO1A has been reported to scaffold Rac-family signaling complexes and Rac signaling is known to regulate cytoskeletal rearrangement and cell migration (14). Further studies are needed to determine whether Rac activation and downstream actin dynamics mediate CORO1A-associated catabolic and migration-related responses under inflammatory stimulation.
Immune infiltration analysis showed that samples with high CORO1A expression exhibited increased infiltration of M2 macrophages. Although M2 macrophages are generally considered anti-inflammatory and repair-associated, their functions in chronically inflamed tissues may be highly context-dependent; persistent M2-like polarization does not necessarily indicate effective inflammation resolution or tissue repair, but may instead reflect dysregulated repair, immune remodeling, or unresolved inflammation. This association may involve both macrophage-intrinsic and microenvironment-dependent mechanisms. One testable explanation is that CORO1A may be linked to TGF-β/SMAD3-related signaling, as TGF-β is an important inducer of M2-like macrophage polarization and Coronin 1A has been reported to regulate TGF-β receptor/SMAD3 signaling in Th17 CD4+ T cells (26, 29). In addition, studies in other disease contexts have linked CORO1A to PI3K/AKT/mTOR-related signaling and myeloid cell regulation (30, 31), providing indirect support for a possible connection between CORO1A-associated signaling programs and macrophage polarization. Alternatively, CORO1A-high samples may reflect broader cytokine- and chemokine-associated inflammatory programs, including pathways related to NF-κB or IL-6/STAT3 signaling (32, 33), which could indirectly influence macrophage recruitment or polarization by reshaping the local inflammatory microenvironment. However, the current data cannot distinguish whether increased M2 macrophage infiltration reflects macrophage-intrinsic CORO1A activity, paracrine signals derived from CORO1A-high chondrocyte states, or both. Therefore, the association between CORO1A expression and M2 macrophage infiltration should be interpreted as hypothesis-generating rather than causal.
Several limitations should be acknowledged. First, the scRNA-seq dataset is based on a limited number of samples, restricting generalizability. Second, the in vitro model reflects acute cytokine stimulation and does not capture chronic disease dynamics. Third, the immune infiltration analysis is correlative and does not establish causality. Fourth, the absence of a matched isotype control limits definitive assessment of the specificity of immunofluorescence staining. Finally, the use of macroscopically unaffected cartilage as a control may not fully represent true healthy tissue, as such regions may already exhibit molecular alterations.
In addition, while this study focuses on ECM catabolic processes, impaired matrix synthesis is also a critical component of OA progression. Future studies are needed to determine whether CORO1A affects anabolic pathways, including COL2 and Aggrecan expression.
In conclusion, our study identifies CORO1A as a context-associated regulator linked to inflammatory signaling and chondrocyte state remodeling in OA. The selective enrichment of CORO1A in inflammatory and proliferative chondrocytes suggests that OA progression is associated with the activation of specific CORO1A-high pathogenic programs rather than uniform activation across all chondrocytes. Through multi-omics analysis and functional perturbation, we propose that CORO1A represents a testable component of a cartilage-immune network associated with inflammatory-catabolic responses, providing a basis for further mechanistic studies.
4 Methods
4.1 Data acquisition
Transcriptomic datasets related to osteoarthritis (OA) and cartilage were sourced from the Gene Expression Omnibus (GEO) database (https://www.ncbi.nlm.nih.gov/geo/). For this analysis of transcriptional changes in OA-affected tissues, two specific datasets were selected. One of these is the bulk cartilage gene expression dataset, GSE286154, which involves a paired-sample approach utilizing 16 cartilage specimens from eight OA patients. Each patient provided cartilage samples from both less affected areas, which were relatively smooth, and more severely damaged regions that bear greater weight. The classification of these samples was guided by the Mankin scoring system, a method that was subsequently employed for expression and network analysis differentiation.
After collecting and reviewing the data, we selected the single-cell RNA sequencing dataset GSE169454 (34) for our analysis of cellular diversity. Ultimately, we investigated the expression patterns and distribution of the CORO1A gene across different chondrocyte subpopulations in individual samples from both osteoarthritis patients and non-osteoarthritis controls, focusing on articular cartilage specimens.
4.2 Analysis of differential expression and functional enrichment
The analysis of differential gene expression was conducted utilizing the limma R package (35) which compared samples of damaged cartilage to those of intact cartilage from the same osteoarthritis patients in the GSE2861454 dataset. To address the paired nature of the study, patient identities were included as a blocking factor in the linear model. Genes were deemed significantly differentially expressed if they had an adjusted P-value (using the Benjamini-Hochberg method) below 0.05 and an absolute log2 fold change exceeding 1.896. Visualization of the differentially expressed genes (DEGs) was carried out using the pheatmap and ggplot2 packages. Following this, enrichment analyses for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were performed with the clusterProfiler package, considering a significance threshold of less than 0.05.
4.3 Gene set enrichment analysis
We categorized the gene expression profiles based on CORO1A levels into two groups: CORO1A-low and CORO1A-high, using the median as a threshold. Phenotypic data were organized in CLS format, while mRNA expression data were compiled in GCT format. Subsequently, we conducted Gene Set Enrichment Analysis (GSEA) (36) on the differentially expressed genes and associated signaling pathways, utilizing Hallmark gene sets from the MSigDB database. The gene sets of particular interest were those that exhibited a normalized enrichment score (NES) with an absolute value exceeding 1.5 and a false discovery rate (FDR) below 0.05.
4.4 Identification of OA-related modules with WGCNA
Gene expression data underwent normalization and correction through the NormalizeBetweenArrays algorithm. Subsequently, a weighted gene co-expression network was constructed utilizing the WGCNA package in R. The soft-thresholding power was set to 11, achieving a scale-free topology with the highest R2 value. Gene hierarchical clustering was performed using a topological overlap matrix (TOM) to assess pairwise dissimilarity, employing average linkage and dynamic tree cutting methods. Module definition was based on a minimum of 75 genes, and modules with similar expression patterns (height cut-off = 0.25) were merged. Genes within the same module exhibited high topological overlap, indicating co-expression. Significant correlations with osteoarthritis (OA) traits were identified, and the modules with the strongest associations were selected for further investigation. To visualize the relationship between module membership (inter-module connectivity) and gene significance in OA, scatter plots were created for the key modules of interest (Supplementary Table 3), with the connection quantified using the Pearson correlation coefficient. Hub genes within a module were identified as those with high intramodular connectivity and significant relevance to OA. Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed on the most relevant OA module, the CORO1A gene, along with its direct neighbors, using the cluster Profiler R package.
4.5 Construction of protein-protein interaction network and identification of hub genes
A network of protein-protein interactions (PPI) was created to elucidate the functional relationships among the differentially expressed genes (DEGs) and identify key hub genes. The DEGs were submitted to the STRING database (version 11.5), where only interactions with a minimum required score of 0.7 were retained to ensure high-quality data. This interaction data was then imported into Cytoscape software (version 3.9.1) for network visualization. The CytoHubba plug-in was employed to calculate centrality metrics for each node within the network. Gene rankings were primarily determined using two topological algorithms: Betweenness and Stress centrality. The top 40 genes from each ranking were selected, resulting in an intersection of 39 high-confidence candidate hub genes for further analysis. It is important to mention that the subsequent focused investigation of CORO1A was based on a hypothesis, as it has been previously linked to the pathophysiology of arthritis (20).
4.6 Immune infiltration analysis
The estimation of immune cell fractions was conducted using CIBERSORT (37). along with the LM22 signature matrix and involved 1,000 permutations. Only those samples that demonstrated a CIBERSORT deconvolution P-value of less than 0.05 were included in subsequent correlation analyses. The relationship between CORO1A expression and immune cell fractions was analyzed using Spearman correlation. To ascertain the immune cell fractions through CIBERSORT, the analysis utilized 1,000 permutations with the LM22 signature matrix. Samples meeting the criteria of a CIBERSORT deconvolution P-value below 0.05 and the 1,000 permutations were retained for further correlation studies. The connections between CORO1A expression and immune cell proportions were assessed via Spearman correlation.
4.7 Single-cell transcriptomics identifies OA-associated chondrocyte subpopulations enriched for CORO1A
The GSE169454 dataset underwent processing and analysis of single-cell RNA sequencing (scRNA-seq) utilizing the Seurat R package (version 5.0+). Initially, a quality control (QC) procedure was applied to filter the dataset, ensuring that only genes expressed in a minimum of three cells were retained. Cells were selected based on having between 300 and 4,500 unique genes, while those with mitochondrial gene content exceeding 20% of the total were excluded. Gene expression counts were normalized using the logNormalization method. The top 2,000 highly variable genes (HVGs) were identified through the variance-stabilizing transformation (vst) technique in FindVariableFeatures (Figure 4A). Data scaling was performed with the ScaleData function, which also accounted for variations related to mitochondrial gene expression. Following this, principal component analysis (PCA) was conducted, and the optimal number of principal components for subsequent analysis was determined via an elbow plot. To mitigate potential batch effects on individual sample data, the R package Harmony was employed for further correction. Cell clustering was achieved by constructing a common nearest neighbor graph with the FindNeighbors function, followed by community detection using the FindClusters function at a resolution of 0.4. The resulting clusters were visualized in two dimensions using Uniform Manifold Approximation and Projection (UMAP). Finally, cell type annotations were performed by examining the expression of cluster-specific markers (identified with the FindAllMarkers function) alongside well-documented markers from existing literature.
4.8 Human samples
Cartilage tissues from individuals with osteoarthritis (OA) who were receiving total knee replacement surgery were obtained for this study. The research was sanctioned by the Ethical Review Board of the Second Affiliated Hospital of Guangxi Medical University (Approval No. 2024-KYL (017)). Clinical details of the donors are summarized in Supplementary Table 2. To facilitate comparisons between affected and unaffected cartilage while minimizing variations among patients, samples were gathered from both the lateral and medial tibial plateau. Following collection, the samples were either preserved in liquid nitrogen or fixed in 4% paraformaldehyde for subsequent analysis as required. Full-thickness cartilage underwent histological examination and immunohistochemistry (IHC). Additionally, the cartilage was minced and lysed using RIPA buffer containing protease inhibitors. Western blotting was employed to analyze the expression patterns of CORO1As in both damaged and intact regions.
4.9 Animal model
The mice utilized in this study were kept in the SPF-grade laminar flow facilities at the Experimental Animal Center of Guangxi Medical University. The temperature was consistently regulated at 25 ± 2 °C, with humidity levels ranging from 55% to 75%. They had unrestricted access to standard pellet food and water, and their bedding was replaced daily while their cages were regularly disinfected to maintain hygiene. To induce osteoarthritis (OA), we performed destabilization of the medial meniscus (DMM) surgery on 12-week-old mice, which involved making an incision in the capsule of the right knee and severing the medial meniscus-tibial ligament. Mice in the sham group underwent the same surgical procedures but did not have their ligaments cut. After an eight-week period, the mice were euthanized, and their knee joints were extracted for histological and immunohistochemical analysis, focusing on cartilage degeneration. This research received approval from the Ethics Committee of the Animal Care Committee at Guangxi Medical University (Approval No. 202411045).
4.10 Processing of tissue samples, histological examination, and immunohistochemical staining
Biochemical and histological evaluations of articular cartilage were performed using collected samples. Specimens from the knee joints were stored in 4% paraformaldehyde for 24 hours, followed by a decalcification process in 10% EDTA (pH 7.4) for 21 days (or 7 days for mouse knee samples, with the solution changed every 3 days). Subsequently, the samples underwent dehydration and were embedded in paraffin. Cross-sections were prepared at a thickness of 4 μm, and the condition of the articular cartilage was assessed using Safranin O and Hematoxylin and Eosin (H&E) staining. Images were captured using a Cytation5 microscope, and cartilage degradation was quantified using the Osteoarthritis Research Society International (OARSI) scoring system along with ImageJ software. Antigen retrieval was performed with citrate buffer, and endogenous peroxidases were inhibited using 3% H2O2 in methanol. Sections were blocked with goat serum and incubated overnight at 4 °C with a primary antibody against CORO1A (1:300, Cell Signaling Technology). Following this, a horseradish peroxidase (HRP)-conjugated goat anti-rabbit IgG polymer was applied, and the samples were stained with diaminobenzidine (DAB) and hematoxylin. Finally, the sections were dried, fixed, and scanned.
4.11 Extraction of primary cells lines, cell culture, biological reagents and transfection
The primary chondrocytes were extracted from the knee joint cartilage of C57BL/6 wild type mouse pups aged 0 to 3 days postnatally. This extraction involved digesting the cartilage tissue with a solution of type II collagenase.
After digestion, the isolated primary chondrocytes were re-suspended in DMEM/F12 medium supplemented with 10% fetal bovine serum, 50 U/mL penicillin, and 0.05 mg/mL streptomycin. To promote optimal growth, the cells were cultured in sterile monolayers at 37 °C in a humidified environment with 5% CO2. To ensure the reliability of the experimental outcomes, the chondrocytes were not subcultured beyond the third generation, preserving their viability for future studies. To establish an in vitro model of cartilage inflammation, the cultured chondrocytes were treated with 10 ng/mL recombinant mouse IL-1β for 24 hours, simulating an inflammatory setting to observe the chondrocyte response. Additionally, to explore the function of the Corola gene, the cells were transfected with double-stranded siRNA targeting Corola. All siRNA sequences (refer to Supplementary Table 1) were obtained from Sangon Biotech (Sangon, Shanghai, China). This experiment aimed to assess the effects of reduced Corola expression on chondrocyte functionality and inflammatory responses under the specified conditions.
4.12 Immunofluorescence staining
Chondrocytes were cultured on glass coverslips and subsequently fixed using a 4% paraformaldehyde solution (prepared in 0.1 M PBS, pH 7.4) for 15 minutes at ambient temperature. Following this, the cells were rinsed with PBS and permeabilized with 0.5% Triton X-100 for 10 minutes. A blocking step was performed using a 3% BSA solution for one hour. The cells were then incubated overnight at 4 °C with a rabbit antibody against CORO1A (Cell Signaling Technology) at a dilution of 1:100. After additional washing, the samples were treated with goat anti-rabbit IgG conjugated to Alexa Fluor 594 (1:500; ZSGB-BIO) under light-protected conditions at room temperature. Finally, the nuclei were stained with DAPI (1:5000; Sigma-Aldrich) for 10 minutes. The coverslips were mounted using an anti-fade medium and examined with an OLYMPUS IX73 microscope equipped with a 60x oil immersion objective and appropriate filter sets for detecting CORO1A and DAPI-stained nucleic acids.
4.13 Cell scratch assay
The scratch assay procedure begins with the growth of primary cells in six-well plates. After the cells have attached and multiplied, various treatments, such as the introduction of interleukin-1 and siRNA transfection, are applied according to the designated experimental groups. Finally, a wound gap is created by using a 200μL pipette tip to scratch the cell monolayer. The width of this wound is recorded at three time points: at the start of the experiment (0 hours), and then at 24 and 48 hours. Wound closure was monitored using phase-contrast microscopy without additional staining.
4.14 Western blot analysis
Cell lysis was performed using a lysis buffer created by combining protease inhibitors, phosphatase inhibitors, phenylmethylsulfonyl fluoride (PMSF), and RIPA buffer in a ratio of 1:1:1:1000. The lysate was centrifuged to clarify it, and the protein concentration in the pellet was measured using the bicinchoninic acid (BCA) assay. The protein samples were then mixed uniformly with SDS-PAGE loading buffer and subjected to SDS-PAGE electrophoresis under denaturing conditions, followed by transfer to polyvinylidene difluoride (PVDF) membranes. After blocking, the membranes were incubated overnight at 4 °C with primary antibodies: anti-CORO1A (1:1,000; Cell Signaling Technology, USA), anti-MMP13 (1:5,000; Proteintech, USA), and GAPDH (1:50,000; Proteintech, USA) as the internal control. Finally, the membranes were treated with a fluorescent secondary antibody (goat anti-rabbit IgG (H+L), DyLight ® 800 4X PEG; Invitrogen, USA). The protein bands were thoroughly washed with TBST solution and imaged using the ImageQuant LAS-4000 system.
4.15 Statistical analysis
The quantification of immunohistochemical metrics, such as the percentage of stained areas and the proportion of positive cells, was conducted using QuPath software (version 0.6.0) (34). The analysis of the results was performed with R statistical software or GraphPad Prism (version 10.1.2). To determine statistically significant differences in mean expression values, either one-way ANOVA followed by Tukey’s post-hoc test or Student’s t-test was employed. A P value of less than 0.05 was deemed statistically significant.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.
Ethics statement
The studies involving humans were approved by the Ethical Review Committee of The Second Affiliated Hospital of Guangxi Medical University. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study. The animal study was approved by the Animal Care Committee of Guangxi Medical University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
QL: Methodology, Writing – original draft, Writing – review & editing. CN: Data curation, Funding acquisition, Methodology, Writing – review & editing. SG: Writing – original draft, Project administration. DL: Data curation, Supervision, Validation, Writing – review & editing. JZ: Formal analysis, Investigation, Project administration, Writing – review & editing. LH: Conceptualization, Resources, Visualization, Writing – review & editing. JH: Conceptualization, Formal analysis, Investigation, Writing – review & editing. ZT: Conceptualization, Funding acquisition, Investigation, Writing – review & editing.
Funding
The author(s) declared that financial support was received for this work and/or its publication. This project was funded by the National Natural Science Foundation of China (No. 8216140155), the Joint Project on Regional High-Incidence Diseases Research of Guangxi Natural Science Foundation (No. 2022JJB140140), the Guangxi Zhuang Autonomous Region Health Department (No. Z-A20230593), the Joint Project of Youth Science Foundation of Guangxi Medical University (No. GXMUYSFB2026103), and the Project for Enhancing the Basic Research Capacity of Young and Middle-aged Teachers in Guangxi Universities (No. 2023KY0109).
Acknowledgments
Thanking the participants for providing cartilage samples.
Conflict of interest
The author(s) declared that this work 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) declared that generative AI was not 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.2026.1820391/full#supplementary-material
References
1
SteinmetzJDCulbrethGTHaileLMRaffertyQLoJFukutakiKGet al. Global, regional, and national burden of osteoarthritis, 1990–2020 and projections to 2050: a systematic analysis for the global burden of disease study 2021. Lancet Rheumatol. (2023) 5:e508–22. doi: 10.1016/S2665-9913(23)00163-7
2
AllenKDThomaLMGolightlyYM. Epidemiology of osteoarthritis. Osteoarthritis Cartilage. (2022) 30:184–95. doi: 10.1016/j.joca.2021.04.020
3
JohnsonVLHunterDJ. The epidemiology of osteoarthritis. Best Pract Res Clin Rheumatol. (2014) 28:5–15. doi: 10.1016/j.berh.2014.01.004
4
JohnsonVLHunterDJ. Knee osteoarthritis. Ann Intern Med. (2024) 177:ITC129–44. doi: 10.7326/ANNALS-24-01249
5
DeyleGDAllenCSAllisonSCGillNWHandoBRPetersenEJet al. Physical therapy versus glucocorticoid injection for osteoarthritis of the knee. N Engl J Med. (2020) 382:1420–9. doi: 10.1056/NEJMoa1905877
6
O’NeillTWFelsonDT. Mechanisms of osteoarthritis (OA) pain. Curr Osteoporos Rep. (2018) 16:611–6. doi: 10.1007/s11914-018-0477-1
7
LitwicAEdwardsMHDennisonEMCooperC. Epidemiology and burden of osteoarthritis. Br Med Bull. (2013) 105(1):185–99. doi: 10.1093/bmb/lds038
8
WangWChuYLuYXuJZhaoWLiangZet al. Skatole alleviates osteoarthritis by reprogramming macrophage polarization and protecting chondrocytes. Research. (2025) 8:604. doi: 10.34133/research.0604
9
NongJLuGHuangYOthers. Identification of cuproptosis-related subtypes, characterization of immune microenvironment infiltration, and development of a prognosis model for osteoarthritis. Front Immunol. (2023) 14. doi: 10.3389/fimmu.2023.1178794
10
LiMYinHYanZLiHWuJWangYet al. The immune microenvironment in cartilage injury and repair. Acta Biomater. (2022) 140:23–42. doi: 10.1016/j.actbio.2021.12.006
11
SebastianAHumNRMcCoolJLWilsonSPMurugeshDKMartinKAet al. Single-cell RNA-seq reveals changes in immune landscape in post-traumatic osteoarthritis. Front Immunol. (2022) 13. doi: 10.3389/fimmu.2022.938075
12
FernandesTLGomollAHLattermannCHernandezAJBuenoDFAmanoMT. Macrophage: a potential target on cartilage regeneration. Front Immunol. (2020) 11:111. doi: 10.3389/fimmu.2020.00111
13
HousmansBACNeefjesMSurtelDAMVitíkMCremersAVan RhijnLWet al. Synovial fluid from end-stage osteoarthritis induces proliferation and fibrosis of articular chondrocytes via MAPK and RhoGTPase signaling. Osteoarthritis Cartilage. (2022) 30:862–74. doi: 10.1016/j.joca.2021.12.015
14
PickRBegandtDStockerTJSalvermoserMThomeSBöttcherRTet al. Coronin 1A, a novel player in integrin biology, controls neutrophil trafficking in innate immunity. Blood. (2017) 130:847–58. doi: 10.1182/blood-2016-11-749622
15
DengXLiuBJiangQLiGLiJXuK. CREBH promotes autophagy to ameliorate NASH by regulating Coro1a. Biochim Biophys Acta Mol Basis Dis. (2024) 1870:166914. doi: 10.1016/j.bbadis.2023.166914
16
ShaverskyiAHegermannJBrandKLeeKHFögerN. Coronin 1a-mediated F-actin disassembly controls effector function in murine neutrophils. Redox Biol. (2025) 82:103618. doi: 10.1016/j.redox.2025.103618
17
ChenYLinJShiDMiaoYXueFLiuKet al. Identification of WDR74 and TNFRSF12A as biomarkers for early osteoarthritis using machine learning and immunohistochemistry. Front. Immunol. Lausanne: Frontiers Media SA. (2025) 16:1517646. doi: 10.3389/fimmu.2025.1517646
18
OhmaeSNomaNToyomotoMShinoharaMTakeiriMFujiHet al. Actin-binding protein coronin 1A controls osteoclastic bone resorption by regulating lysosomal secretion of cathepsin K. Sci Rep. (2017) 7:41710. doi: 10.1038/srep41710
19
LiangWFengRLiXDuanXFengSChenJet al. A RANKL-UCHL1-sCD13 negative feedback loop limits osteoclastogenesis in subchondral bone to prevent osteoarthritis progression. Nat Commun. (2024) 15:8792. doi: 10.1038/s41467-024-53119-2
20
BalakrishnanLBhattacharjeeMAhmadSNirujogiRSRenuseSSubbannayyaYet al. Differential proteomic analysis of synovial fluid from rheumatoid arthritis and osteoarthritis patients. Clin Proteomics. (2014) 11:1. doi: 10.1186/1559-0275-11-1
21
KhorevaAButovKRNikolaevaEIMartyanovAKulakovskayaEPershinDet al. Novel hemizygous CORO1A variant leads to combined immunodeficiency with defective platelet calcium signaling and cell mobility. J Allergy Clin Immunol Glob. (2023) 3:100172. doi: 10.1016/j.jacig.2023.100172
22
XinRXuYLongDMaoGLiaoHZhangZet al. Mitochonic acid-5 inhibits reactive oxygen species production and improves human chondrocyte survival by upregulating SIRT3-mediated, parkin-dependent mitophagy. Front Pharmacol. (2022) 13:911716. doi: 10.3389/fphar.2022.911716
23
GanDTaoCJinXWuXYanQZhongYet al. Piezo1 activation accelerates osteoarthritis progression and the targeted therapy effect of artemisininPiezo1???????????????????? J Adv Res. (2024) 62:105–17. doi: 10.1016/j.jare.2023.09.040
24
RichettePLatourteA. Hand osteoarthritis: a fresh look. Joint Bone Spine. (2024) 91:105652. doi: 10.1016/j.jbspin.2023.105652
25
Tokarz-DeptułaBMalinowskaMAdamiakMDeptułaW. Coronins and their role in immunological phenomena. Cent-Eur J Immunol. (2016) 41:435–41. doi: 10.5114/ceji.2016.65143
26
KaminskiSHermann-KleiterNMeiselMThuilleNCroninSHaraHet al. Coronin 1A is an essential regulator of the TGFβ receptor/SMAD3 signaling pathway in Th17 CD4+ T cells. J Autoimmun. (2011) 37:198–208. doi: 10.1016/j.jaut.2011.05.018
27
HuangCZengBZhouBChenGZhangQHouWet al. Single-cell transcriptomic analysis of chondrocytes in cartilage and pathogenesis of osteoarthritis. Genes Dis. (2024) 12:101241. doi: 10.1016/j.gendis.2024.101241
28
OkuTKanekoYMurofushiKSeyamaYToyoshimaSTsujiT. Phorbol ester-dependent phosphorylation regulates the association of p57/coronin-1 with the actin cytoskeleton∗. J Biol Chem. (2008) 283:28918–25. doi: 10.1074/jbc.M709990200
29
GengKMaXJiangZGuJHuangWWangWet al. WDR74 facilitates TGF-β/smad pathway activation to promote M2 macrophage polarization and diabetic foot ulcer wound healing in mice. Cell Biol Toxicol. (2023) 39:1577–91. doi: 10.1007/s10565-022-09748-8
30
JinYHuangYZengGHuJLiMTianMet al. Advanced glycation end products regulate macrophage apoptosis and influence the healing of diabetic foot wound through miR-361-3p/CSF1R and PI3K/AKT pathway. Heliyon. (2024) 10:e24598. doi: 10.1016/j.heliyon.2024.e24598
31
DingZShenHXuKWuYWangSYiFet al. Comprehensive analysis of mTORC1 signaling pathway–related genes in the prognosis of HNSCC and the response to chemotherapy and immunotherapy. Front Mol Biosci. (2022) 9:792482. doi: 10.3389/fmolb.2022.792482
32
SiegmundKThuilleNPoschNFresserFBaierG. Novel protein kinase C θ: Coronin 1A complex in T lymphocytes. Cell Commun Signal CCS. (2015) 13:22. doi: 10.1186/s12964-015-0100-3
33
ZhaoCXieRQianQYanJWangHWangX. Triclosan induced zebrafish immunotoxicity by targeting miR-19a and its gene socs3b to activate IL-6/STAT3 signaling pathway. Sci Tot Environ. (2022) 815:152916. doi: 10.1016/j.scitotenv.2022.152916
34
FuWHettinghouseAChenYHuWDingXChenMet al. 14-3-3 epsilon is an intracellular component of TNFR2 receptor complex and its activation protects against osteoarthritis. Ann Rheum Dis. (2021) 80:1615–27. doi: 10.1136/annrheumdis-2021-220000
35
RitchieMEPhipsonBWuDHuYLawCWShiWet al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007
36
ChenYLiCWangNWuZZhangJYanJet al. Identification of LINC00654-NINL regulatory axis in diffuse large B-cell lymphoma in silico analysis. Front Oncol. (2022) 12:883301. doi: 10.3389/fonc.2022.883301
37
NewmanAMLiuCLGreenMRGentlesAJFengWXuYet al. Robust enumeration of cell subsets from tissue expression profiles. Nat Methods. (2015) 12:453–7. doi: 10.1038/nmeth.3337
Summary
Keywords
bioinformatics analysis, chondrocytes, CORO1A, immune microenvironment, osteoarthritis, therapeutic targets
Citation
Liang Q, Ning C, Gao S, Liu D, Zhao J, Huang L, Hou J and Tan Z (2026) CORO1A links inflammatory chondrocyte subpopulations to immune microenvironment alterations in osteoarthritis: an integrative multi-omics and single-cell study. Front. Immunol. 17:1820391. doi: 10.3389/fimmu.2026.1820391
Received
09 March 2026
Revised
09 May 2026
Accepted
26 May 2026
Published
11 June 2026
Volume
17 - 2026
Edited by
Verena Fischer, Universitaetsklinikum Ulm, Germany
Reviewed by
Stinne Ravn Greisen, Aarhus University Hospital, Denmark
Yunkun Qu, Huazhong University of Science and Technology, China
Updates
Copyright
© 2026 Liang, Ning, Gao, Liu, Zhao, Huang, Hou and Tan.
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: Zhen Tan, flextz@126.com; Jun Hou, hj@sr.gxmu.edu.cn; Linke Huang, huanglinke28@sr.gxmu.edu.cn
†These authors have contributed equally to this work
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.