key: cord-0715492-ttg6l7cm
authors: Lu, Dinggui; Wei, Jihua; Chen, Jian; Zhao, Jingjie; Wang, Jiajia; Gong, Yuanxun; Wei, Liuzhi; Wei, Qiuju; Ban, Huadeng; Li, Yueyong; Wang, Zechen; Luo, Changtai; Zhou, Haidong; Shen, Jiajia; Liao, Qiujiao; He, Siyuan; Zhang, Weiyang; Luo, Qunqiang; Xie, Kegong; Song, Jian; Meng, Lingzhang
title: Apelin Alleviates Meniscus Endothelial Cell Apoptosis in Osteoarthritis
date: 2022-01-13
journal: Dis Markers
DOI: 10.1155/2022/3556372
sha: 05ec4939401aa74a3779bb3dfeac541714700ccc
doc_id: 715492
cord_uid: ttg6l7cm

Osteoarthritis (OA) is a degenerative disease characterized by articular cartilage and/or chondrocyte destruction, and although it has long been considered as a primary disease, the importance of meniscus endothelial cell modulation in the subchondral microenvironment has recently drawn attention. Previous studies have shown that apelin could potentially inhibit cellular apoptosis; however, it remains unclear whether apelin could play a protective role in protecting the endothelium in the OA meniscus. In this study, with the advantages of single-cell RNA sequencing (scRNA-seq) data, in combination with flow cytometry, we identified two endothelial subclusters in the meniscus, featured by high expression of Homeobox A13 (HOXA13) and Ras Protein-Specific Guanine Nucleotide Releasing Factor 2 (RASGRF2), respectively. Compared with control patients, both subclusters decreased in absolute cell numbers and exhibited downregulated APJ endogenous ligand (APLN, coding for apelin) and upregulated apelin receptor (APLNR, coding apelin receptor). Furthermore, we confirmed that in OA, decreased endothelial cell numbers, including both subclusters, were related to intrinsic apoptosis factors: one more relevant to caspase 3 (CASP3) and the other to BH3-Interacting Domain Death agonist (BID). In vitro culturing of meniscal endothelial cells purified from patients proved that apelin could significantly inhibit apoptosis by downregulating these two factors in endothelial cell subclusters, suggesting that apelin could potentially serve as a therapeutic target for patients with OA.

Osteoarthritis (OA), a degenerative joint disease, is generally associated with obesity and ageing. Although overexercise can elevate the risk of OA, it is usually considered a genetically driven disease [1] . In recent years, the pathological variations of the cellular components, especially the pathogenic meniscus and vasculature in OA, have garnered increasing attention [2] [3] [4] . Although the vasculature in the joint was observed to be abnormal in patients with OA [1, 5] , its molecular mechanism is still difficult to understand, as research is impacted by a technique limitation in acquiring sufficient endothelial cells. The rise of the scRNA-seq technique allows for exploring transcriptome in a highthroughput manner at a single-cell level so that the endothelial cell heterogeneity can be studied more conveniently [6, 7] . Notably, endothelial cells in different organs (e.g., brain, lung, and heart) exhibit tissue-specific heterogeneity, characterized by molecular signatures such as those related to pathways and cell-surface markers [6] [7] [8] . Particularly, research on tissue-specific endothelial cell heterogeneity will help to understand the development of various diseases.

In this study, scRNA-seq data of meniscus biopsies isolated from patients with OA were adopted, and endothelial cells were isolated to explore heterogeneity by complying with genetic features. On this basis, the degenerative meniscus vasculature in endothelial cell subclusters of patients with OA was confirmed. Furthermore, intrinsic apoptosis factors (i.e., CASP3 and BID) were reported to be related to two meniscus endothelial cell subclusters in patients with OA.

Apelin is a pleiotropic peptide known as an adipokine [9] . It has been shown to possess various biological activities, such as attenuation of endoplasmic reticulum stress and apoptosis in ischemic stroke [10] , influencing mitosis and apoptosis in intestinal cell lines [11] , and inhibiting bone marrow mesenchymal stem cell apoptosis [11] . However, whether this molecule could protect the meniscus endothelial cells in patients with OA is still unknown. In this study, the discovery of upregulation of the apelin receptor (APLNR), a marker for endothelial cell injury [12] , prompted us to investigate the role of apelin as a potential therapeutic method to alleviate meniscal endothelium injury in patients with OA.

This study reported that apelin alleviates meniscal endothelial cell apoptosis in patients with OA by inhibiting CASP3 and BID in different EC clusters.

2.1. Data Availability. The data and code generated or analyzed in this study are available from the corresponding author upon reasonable request.

2.2. Acquisition of Meniscus scRNA-seq Data. The human meniscus scRNA-seq data were provided by Sun et al. [13] , retrieved from the NCBI Gene Expression Omnibus database accessed at GSE133449. The dataset consisted of three meniscus biopsies isolated from patients during supra-knee amputation but without any history of OA, rheumatoid arthritis, or other diseases related to the meniscus. Therefore, the samples were used as controls. In addition, four meniscal biopsies were isolated from patients with OA during knee arthroplasty surgery. In the present study, cells with at least 50 genes, and genes expressed in at least three cells were analyzed.

Overall, control meniscus biopsies were isolated from 52 patients who underwent arthroscopic surgery after acute injury to the meniscus. A total of 84 patients with OA having undergone arthroscopic surgery provided biopsies. Each patient underwent a single biopsy. The biopsies were taken for the pathological examination. Furthermore, none of the patients had a history of knee-related disease.

2.4. Antibodies. In this study, the antibodies used for flow cytometry were as follows: 7-AAD (Invitrogen, #A1310), anti-human CD45 (clone HI30, eBiosciences, # 48-0459-42), anti-PECAM1 (clone WM59, Biolegend, #303106), antihuman/mouse/rat HOXA13 (Biorbyt, #orb461791), antihuman cleaved caspase 3 (clone Asp175, Cell Signalling, #9602), anti-human/mouse/rat RASGRF2-Biotin-conjugated (Biocompare, #MBS544685), anti-Bid Cleavage Site antibody (abcam, ab10640), streptavidin-APC (Invitrogen, #17-4317-82), and anti-human APJ/APLNR (BIOSS, #bs-2430R-Cy7).

Cytometry. The meniscus biopsies were cut into small pieces and digested with collagenase IV (Gibco, #17017-019) at 37°C for 30 min. In addition, PBS/ 0.5%BSA/2 mM EDTA was added 10 min before the completion of digestion. The single-cell suspension was produced by filtering through a 100 μm cell strainer. The nonspecific binding was blocked with an Fc receptor binding inhibitor (Invitrogen, #14-9161-73) for 15 min. Subsequently, the cells were washed with cold PBS/0.5%BSA and incubated with antibodies for 15 min on ice. The cells were fixed and permeabilized for intracellular staining with a commercial Cytofix/Cytoperm kit (BD, #554717). After washing twice with PBS/0.5%BSA, the cells were incubated with the antibodies for 20 min on ice. Specific to biotinylated anti-RASGFR2, the cells were washed twice and then incubated with streptavidin-APC for another 15 min on ice. After washing twice, the cells were resuspended in PBS/ 0.5%BSA/2 mM EDTA. Biopsies were measured using a flow cytometer (Thermal Fisher Attune NxT).

2.6. Fluorescence-Activated Cell Sorting. After the surface staining of fluorescent-labeled antibodies, the cells from patients with OA were placed into a FACS cell sorter (BD FACSMelody) to produce purified meniscus endothelial cells. 3 Disease Markers penicillin-streptomycin-neomycin (Sigma-Aldrich, #P4083) at a density of 2 × 10 4 cells/well. Subsequently, the cells were cultured overnight in 48-well plates. After adding apelin (100 nM, Sigma-Aldrich, #A6469) and/or ML221 (blocking apelin receptor, 100 nM, Sigma-Aldrich, #SML0919) for 8 h, the cells were harvested for indepth analysis.

After the meniscus endothelial cells purified from patients with OA were cultured for 24 h, mRNA was extracted using the TurboCapure mRNA kit (QIAGEN, #72251). cDNA was synthesized using TaqMan Reverse Transcription Reagents (Invitrogen, #4304134), and qPCR was performed using the Brilliant III Ultra-Fast QPCR Kit (Agilent, #60083) with a LightCycler 96 machine (ROCHE). Primers used for qPCR: APLN forward 5 ′ -GTCT CCTCCATAGATTGGTCTGC-3 ′ and reverse 5 ′ -GGAATC ATCCAAACTACAGCCAG-3 ′ ; CASP3 forward 5 ′ -GAAA TTGTGGAATTGATGCGTGA-3 ′ and reverse 5 ′ -CTAC AACGATCCCCTCTGAAAAA-3 ′ ; BID forward 5 ′ -ATGG ACCGTAGCATCC CTCC-3 ′ and reverse 5 ′ -GTAGGT GCGTAGGTTCTGGT-3′; GAPDH forward 5′-GGAGCG AGATCCCTCCAAAAT-3′ and reverse 5′-GGCTGTTGT CATACTTCTCATGG-3′ (internal control).

2.9. Statistical Analysis. Statistical analysis was performed using GraphPad Prism v5 software (La Jolla, USA). A nonparametric test was used for the data analysis. Results are presented as the mean ± SD in dot plots. NS: no significance (p > 0:05), * p < 0:05, * * p < 0:01, * * * p < 0:001.

It is reported that in meniscus, the immune cells and vasculature mainly distributed in the red zone [14] . Meniscus scRNA-seq data were analyzed using the R package Seurat [15] , and the cells were clustered at a resolution of 0.6. Overall, nine clusters were separated from the control meniscus biopsies, while 11 clusters were separated from the OA biopsies ( Figure 1(a) ). To identify endothelial cells, three endothelial cell feature genes (i.e., CD93, CDH5, and PECAM1) [16, 17] were employed (Figure 1 (b) and Supplementary Figure 2 ). In the control meniscus, cluster 6 was clearly identified as the endothelial cells for the high expression of the three feature genes observed. In OA biopsies, cluster 10 was identified as endothelial cells because of the high expression of CD93, CDH5, and PECAM1. Although, cluster 5 also exhibited the obvious expression of CD93, a high expression of monocyte/macrophage markers, CD14, and MS4A7 [13, 18] was observed (Figure 1 (b) and Supplementary Figure 2 ). Cluster 5 was excluded from the subsequent analysis to diminish ambiguous results because it was reported that certain monocyte subsets could also express CD93 [19] .

As revealed from the scRNA-seq measurements, endothelial cells comprised approximately 2% (cluster 6) of total meniscus cells in the controls, whereas they decreased to about 1% (cluster 10) in patients with OA. Flow cytometry experiments were performed to compare the meniscus endothelial cells between the controls and patients with OA to confirm these phenomena. Consistently, endothelial cells decreased by 36% in patients with OA compared with the control biopsies (Figures 1(d) and 1(e) ). The decreased endothelial cells could result from facilitated apoptosis, especially in chronic inflammation [20, 21] . and patients with OA. The unbiased subclustering analysis identified two subpopulations in the total meniscus endothelial cells (Figure 2(a) ). Notably, HOXA13 was predominantly expressed in the first subcluster, whereas RASGRF2 was preferentially expressed in the second subcluster; therefore, these two subclusters were termed as EC (HOXA13) and EC (RASGRF2), respectively (Figures 2(a) and 2(b) ). The genetic profiles of the two subclusters were compared (Supplementary Table 1 ) to confirm the heterogeneity of the meniscus endothelial cells. The heat map shows the disparity of the top 10 genes between the two subclusters ( Figure 2(c) ). Based on the differences in genetic profiles between the meniscus endothelial cell subclusters, the gene set enrichment analysis showed that EC (HOXA13) plays a pivotal role in the cellular response to organic substances, responding to chemical stimulus, regulation of response to stimulus, and modulation of extracellular region. In addition, another EC subcluster EC (RASGRF2) could potentially participate in the neutrophil apoptotic process and modulation of type V collagen (Figure 2(d) ). To compare the above two subclusters between the controls and patients with OA, split UMAPs showed that both subclusters decreased in a number of patients with OA (Figure 2(e) ). Moreover, flow cytometry analysis of isolated meniscus biopsies from patients showed decreased numbers of endothelial subclusters (Figures 2(f) and 2(g) ). 

Control OA 9 Disease Markers numbers in patients with OA, differentially expressed genes (DEGs) were calculated from the scRNA-seq data. Overall, 2831 DEGs were identified in EC (HOXA13), and 5297 DEGs were identified in EC (RASGRF2) (Figure 3(a) , Supplemental Tables 2 and 3 ). The most significant DEGs were labeled in volcano plots by cluster. Specifically, drastically decreased APLN was identified in both the subclusters (Figure 3(a), Supplementary Figure 3) since it was reported that apelin significantly inhibits apoptosis [9, [22] [23] [24] [25] . Interestingly, both meniscal endothelial subclusters exhibited upregulated expression of APLNR, coding the apelin receptors (Figure 3(a), Supplementary Figure 3) . Consistently, meniscus endothelial cells purified from patients showed downregulated APLN expression in patients with OA (Figure 3(b) ). Moreover, flow cytometry analysis showed that APLNR was upregulated in both meniscal endothelial subclusters in patients with OA (Figures 3(c) and 3(d) ), demonstrating that meniscus endothelial cells in patients with OA were severely injured and undergoing apoptosis [12] . This observation revealed that apelin could function as a therapeutic target against apoptosis in meniscal endothelial cells in both subclusters.

Apoptosis-related genes were extracted from DEGs (Supplemental Tables 2 and 3 ) and plotted by subclusters (Figure 3(e) ) to determine the genes potentially involved in the apoptosis of meniscus endothelial cells. It is noteworthy that CASP3 showed upregulated expression in EC (HOXA13) of patients with OA, while BID was upregulated in EC (RASGRF2) of patients with OA. In addition, other apoptosis-related genes selected from DEGs (e.g., CASP7, FAS, TNFRSF1A, and TNFSF10) were downregulated in patients with OA. Notably, CASP3 and BID could function 10 Disease Markers as intrinsic apoptosis genes [26] [27] [28] , while FAS, TNFRSF1A, and TNFSF10 are extrinsic apoptosis genes [29] . This observation suggests that decreased meniscus endothelial cell numbers in patients with OA could be attributed to intrinsic apoptosis rather than extrinsic apoptosis.

Fluorescent-activated cell sorting experiments were performed to acquire highly purified endothelial cells from meniscus biopsies to verify whether apelin can inhibit intrinsic apoptosis in the meniscus endothelial cells of patients with OA. The isolated endothelial cells were treated with BSA, APLN, or APLN and its receptor blocker (ML221) for 8 h. As indicated by flow cytometry analysis, early apoptosis (annexin V + PI-cells) of meniscus endothelial cells isolated from patients with OA was significantly inhibited by apelin, whereas it was abolished by ML221, a blocker of the apelin receptor [30] (Figures 4(a) and 4(b) ). To clarify whether it could be attributed to apelin's effect on intrinsic apoptosis genes, RNA was extracted from cultured meniscus endothelial cells isolated from patients with OA. As indicated by the quantitative PCR analysis, CASP3 and BID expression levels were inhibited by apelin treatment and restored by blocking the apelin receptor (treated with ML221). Furthermore, flow cytometry was performed to analyze cleaved caspase 3 (activated caspase 3) in EC (HOXA13) and cleavage site Bid in EC (RASGRF2). Consistently, apelin significantly inhibited these two activated intrinsic apoptosis factors in endothelial subclusters, and treatment with ML221 restored their expression.

Blood vessels modulate the tissue microenvironment and provide nutrients and immune factors affecting other cells [3, 31] , in which endothelial cells play a significant role in maintaining tissue homeostasis and in the development of various diseases. Therefore, extensive studies have been conducted to explore the biological and pathological features of endothelial cells in the brain, lung, heart, skin, and other tissues [6] [7] [8] 12] . Existing studies indicate that tissue-specific endothelial cells exhibit distinguished subclusters, genetic profiles, and biomarkers among diseases. However, biological/pathological features of meniscus endothelial cells have rarely been studied due to extremely poor vascularity in the human meniscus [32] , resulting in meniscus endothelial cell heterogeneity still considered an enigma. Accordingly, it is more urgent to illustrate the molecular mechanisms of meniscus endothelial cells under pathological conditions (e.g., in OA). This study used data from scRNA-seq of meniscal biopsies isolated from patients with OA [16] . This method effectively maintained cellularity, especially for rare cells, including endothelial cells, and allowed exploring the transcriptome of meniscus endothelial cells in a highthroughput manner. Indeed, clarifying the genetic features of meniscus endothelial cells will help gain insights into the development of OA [33] and cultivate new strategies against this disease. Normally, two feature genes were adopted to identify specific cell types; however, in this study, three feature genes (CD93, CDH5, and PECAM1) were used to isolate meniscus endothelial cells for downstream analysis, making it more reliable. The scRNA-seq analysis confirmed that the frequency of meniscus endothelial cells dropped to 1% in patients with OA, and flow cytometry analysis of freshly isolated meniscus biopsies confirmed this tendency. Furthermore, this phenomenon revealed that the cells underwent degeneration and injury.

Although OA has long been considered a primary degenerative disease characterized by an injured joint, the abnormal microenvironment of the subchondral bone has attracted attention in recent years [3] . This study confirmed that the blood vessels in the meniscus of patients with OA underwent pathological variations. Based on a highthroughput genetic profile, this study initially reported that meniscus endothelial cells contain two subclusters, each of which could be identified by feature genes: HOXA13 and RASGRF2. Interestingly, HOXA13 has been reported to be critical for vascular patterning [34, 35] , and RASGRF2 is involved in the response to endoplasmic reticulum stress in endothelial cells [36] . The drastic upregulation of the expression of these two genes in the meniscus endothelial cells of patients with OA indicated disordered cellular biological features, which inevitably contributed to the pathological meniscus microenvironment in patients with OA.

In addition, in patients with OA, both endothelial cell subclusters showed vigorous expression of APLNR (coding the receptor for apelin), a marker for injured and apoptotic endothelial cells [12] . Decreased endothelial cell numbers in the meniscus of patients with OA support this hypothesis. In contrast, a higher expression of APLNR could serve therapeutic purposes since its ligand, apelin, has an antiapoptotic effect [9] . Consistently, the downregulated expression of APLN has been observed in both meniscal endothelial cells in patients with OA, which could result in enhanced apoptosis of endothelial cells. Unexpectedly, two intrinsic apoptosis factors, CASP3 and BID, were upregulated in two endothelial cell subclusters in patients with OA. Downstreaming cell culturing of purified meniscus endothelial cells from patients with OA showed a pronounced outcome. This study demonstrated that apelin could function as a therapeutic agent against OA by protecting the meniscus endothelial cells.

Moreover, apelin also possesses the capability of controlling blood pressure [37] , promoting angiogenesis [38] , and interacting with ACE2, the main receptor for novel coronavirus COVID-19 [39, 40] . Apelin has been reported to rescue heart failure potentially and pulmonary arterial hypertension [41] , two syndromes frequently reported from COVID-19 infected patients [42, 43] . Therefore, evaluating the mentioned perspectives with biopsies from COVID-19 infected patients will help develop therapeutic strategies against the COVID-19 pandemic.

Vascular mechanisms in osteoarthritis

Vascular pathology and osteoarthritis: a systematic review

Microenvironment in subchondral bone: predominant regulator for the treatment of osteoarthritis

Functional MRI mapping of human meniscus functionality and its relation to degeneration

Vascular pathology and osteoarthritis

A molecular atlas of cell types and zonation in the brain vasculature

Single-cell RNA sequencing of mouse brain and lung vascular and vessel-associated cell types

Endothelial heterogeneity across distinct vascular beds during homeostasis and inflammation

Anti-apoptotic effect of apelin in human placenta: studies on BeWo cells and villous explants from third-trimester human pregnancy

Low dose of Apelin-36 attenuates ER stress-associated apoptosis in rats with ischemic stroke

Effect of apelin on mitosis, apoptosis and 12 Disease Markers DNA repair enzyme OGG 1/2 expression in intestinal cell lines IEC-6 and Caco-2

Single cell RNA sequencing identifies HSPG2 and APLNR as markers of endothelial cell injury in systemic sclerosis skin

Single-cell transcriptome analysis profile of meniscal tissue macrophages in human osteoarthritis

3D-bioprinting ready-toimplant anisotropic menisci recapitulate healthy meniscus phenotype and prevent secondary joint degeneration

Comprehensive integration of single-cell data

Single-cell RNA-seq analysis identifies meniscus progenitors and reveals the progression of meniscus degeneration

Endothelial rap 1 (Ras-association proximate 1) restricts inflammatory signaling to protect from the progression of atherosclerosis

Deciphering human macrophage development at single-cell resolution

Human monocyte subsets and phenotypes in major chronic inflammatory diseases

New markers of inflammation and endothelial cell activation: part I

New paradigms in inflammatory signaling in vascular endothelial cells

Lateral intracerebroventricular injection of Apelin-13 inhibits apoptosis after cerebral ischemia/reperfusion injury

Novel pathogenesis: regulation of apoptosis by apelin/APJ system

The effects of apoptosis and apelin on lymph node metastasis in invasive breast carcinomas

Protective effect of apelin on cultured rat bone marrow mesenchymal stem cells against apoptosis

Caspase-9, caspase-3 and caspase-7 have distinct roles during intrinsic apoptosis

BCL-2 family isoforms in apoptosis and cancer

The BCL-2 protein family, BH3-mimetics and cancer therapy

Apoptosis through death receptors in temporal lobe epilepsy-associated hippocampal sclerosis

Apelin (APLN) and apelin receptor (APLNR) in human ovary: expression, signaling, and regulation of steroidogenesis in primary human luteinized granulosa cells

Endothelial cells as tools to model tissue microenvironment in hypoxia-dependent pathologies

Vascularity for Healing of Meniscus Repairs

Vascular endothelial cell-secreted exosomes facilitate osteoarthritis pathogenesis by promoting chondrocyte apoptosis

HOXA13 is essential for placental vascular patterning and labyrinth endothelial specification

Analysis of homeobox gene action may reveal novel angiogenic pathways in normal placental vasculature and in clinical pregnancy disorders associated with abnormal placental angiogenesis

RasGRF couples Nox4-Dependent endoplasmic reticulum signaling to Ras

The role of apelin in cardiovascular diseases, obesity and cancer

Apelin signaling drives vascular endothelial cells toward a pro-angiogenic state

Structural basis for the recognition of SARS-CoV-2 by full-length human ACE2

Molecular interaction and inhibition of SARS-CoV-2 binding to the ACE2 receptor

Interaction between the apelinergic system and ACE2 in the cardiovascular system: therapeutic implications

COVID-19 and cardiovascular disease: from basic mechanisms to clinical perspectives

COVID-19 and pulmonary arterial hypertension: early data and many questions

The data and code generated or analyzed in this study are available from the corresponding author upon reasonable request.

All of the research protocols have been reviewed and approved by the Ethical Committee of Youjiang Medical University for Nationalities. Patients were informed and consented on this study. All the experiments including data analysis were performed at Youjiang Medical University for Nationalities.

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