1. Medicine

A single-cell atlas depicting the cellular and molecular features in human anterior cruciate ligamental degeneration: A single cell combined spatial transcriptomics study

  1. Runze Yang
  2. Tianhao Xu
  3. Lei Zhang
  4. Minghao Ge
  5. Liwei Yan
  6. Jian Li
  7. Weili Fu  Is a corresponding author
  1. Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, China
https://doi.org/10.7554/eLife.85700
Share this article
Cite this article
  1. Runze Yang
  2. Tianhao Xu
  3. Lei Zhang
  4. Minghao Ge
  5. Liwei Yan
  6. Jian Li
  7. Weili Fu
(2023)
A single-cell atlas depicting the cellular and molecular features in human anterior cruciate ligamental degeneration: A single cell combined spatial transcriptomics study
eLife 12:e85700.
https://doi.org/10.7554/eLife.85700
  2. Download BibTeX
  3. Download .RIS

Abstract

Background:

To systematically identify cell types in the human ligament, investigate how ligamental cell identities, functions, and interactions participated in the process of ligamental degeneration, and explore the changes of ligamental microenvironment homeostasis in the disease progression.

Methods:

Using single-cell RNA sequencing and spatial RNA sequencing of approximately 49,356 cells, we created a comprehensive cell atlas of healthy and degenerated human anterior cruciate ligaments. We explored the variations of the cell subtypes’ spatial distributions and the different processes involved in the disease progression, linked them with the ligamental degeneration process using computational analysis, and verified findings with immunohistochemical and immunofluorescent staining.

Results:

We identified new fibroblast subgroups that contributed to the disease, mapped out their spatial distribution in the tissue and revealed two dynamic trajectories in the process of the degenerative process. We compared the cellular interactions between different tissue states and identified important signaling pathways that may contribute to the disease.

Conclusions:

This cell atlas provides the molecular foundation for investigating how ligamental cell identities, biochemical functions, and interactions contributed to the ligamental degeneration process. The discoveries revealed the pathogenesis of ligamental degeneration at the single-cell and spatial level, which is characterized by extracellular matrix remodeling. Our results provide new insights into the control of ligamental degeneration and potential clues to developing novel diagnostic and therapeutic strategies.

Funding:

This study was funded by the National Natural Science Foundation of China (81972123, 82172508, 82372490) and 1.3.5 Project for Disciplines of Excellence of West China Hospital Sichuan University (ZYJC21030, ZY2017301).

Editor's evaluation

This is a valuable study of particular interest to those researchers working in the field. The authors provide a single-cell atlas of healthy and degenerative human anterior cruciate ligament using transcriptomic profiling. The quality of the study was considered valuable as it expands current knowledge.

https://doi.org/10.7554/eLife.85700.sa0

Introduction

Skeletal ligaments, an important part of the musculoskeletal system, are defined as the dense connective tissue that spans the joint and then is fixed at both ends of the bone (Hanhan et al., 2016). One of the main functions of ligaments is mechanical, as they passively maintain joint stability and assist in guiding those joints through their normal range of motion when a tensile loading is applied (Frank, 2004). The knee is the largest and most complicated hinge joint associated with weight bearing in the human body. The anterior cruciate ligament (ACL) is crucial for knee kinematics, especially in rotation, and functions as an anterior/posterior stabilizer (Fleming, 2003). ACL degeneration can gradually result in chronic knee pain, instability, and even poor life quality (Thompson et al., 2015). It has been reported that aging-related degenerative changes in the ACL might also contribute to OA occurrence and progression (Loeser, 2010), but mechanisms of ACL aging and degeneration remain unclear.

The extracellular matrix (ECM) is the physical basis for the biological functions of ligaments, and ACL degeneration is accompanied by alterations in the ECM. The ECM of ACL consists of collagen types I, II, III, and V, elastin, and proteoglycans (Laurencin and Freeman, 2005), and collagen type I is a major determinant of tensile strength (Corps et al., 2006). ACL degeneration is characterized by the disorganization of collagen fibers, cystic changes, mucoid degeneration, and chondroid metaplasia (Hasegawa et al., 2012). Different types of pathological changes correspond to different ECM alterations. For example, mucoid degeneration reflects the degradation of collagen and deposition of new glycosaminoglycans, and cystic changes are devoid of ECM in the diseased area (Hasegawa et al., 2012). It has been reported that some regions in the degenerated ACL have decreased type I collagen, whereas type II, III, and X collagen and aggrecan are significantly increased, and this abnormal ECM production can lead to a biomechanical fault (Hasegawa et al., 2013; Hayashi et al., 2003). In recent years, many scholars have paid attention to the degeneration mechanism of ligaments. TGFβ1, a member of the TGF superfamily, can be used to induce chondrogenic differentiation of ligament-derived stem cells in vitro (Schwarz et al., 2019) and is considered to be an essential molecule contributing to chondroid metaplasia. It also has been reported that mechanical stretching force can promote TGF-β1 production by ligament cells, contributing to ligamental hypertrophy (Nakatani et al., 2002). Some reports suggested that complement cascade could lead to ECM catabolism and contribute to ligament degeneration (Busch et al., 2013; Schulze-Tanzil, 2019). Akihiko Hasegawa et al. suggested that there is inflammation in the degenerative ligament, manifesting as leukocyte infiltration and neovascularization between collagen fibers within ACL substance, and once inflammation is initiated, its synergy with mechanical forces could facilitate ACL degeneration (Fleming et al., 2005; Hasegawa et al., 2012).

There are various types of cells in the ligament (Kharaz et al., 2018), and building a detailed ligamental cell landscape is essential to understanding ligamental characteristics and underlying pathogenesis of ACL degeneration. Currently, the reports on the ACL cellular changes during the process of degeneration mostly focus on chondroid metaplasia and progenitor cells or myofibroblasts recruitment and proliferation (Hasegawa et al., 2013; Mullaji et al., 2008). However, the cellular changes during ligamental degeneration are still unclear. Here, we performed single-cell RNA sequencing (scRNA-seq) and spatial RNA sequencing (spRNA-seq) to obtain an unbiased atlas of ACL cell clusters. Our findings provide a better understanding of the inherent heterogeneity, construct the classifications of fibroblasts, and profile the spatial information of identified cell clusters in the ACL. Notably, we also identified the cell interactions during the disease process and demonstrated the role of FGF and TGF-β signaling pathways in ligament degeneration. Thus, our study reveals the cellular landscape of the human ACL and provides insight that could help to identify therapeutic targets for human ligamental degeneration.

Methods

Human ligament cell sample preparation

The ligament specimens were collected from four joint replacement patients with osteoarthritis (degenerated ACLs) and four amputation patients with severe trauma or malignant bone tumor (healthy ACLs; Supplementary file 1).

All specimens removed were placed immediately in Dulbecco modified Eagle medium (DMEM) free of antibiotics and FBS under 4℃. Ligament samples were rinsed in precooled PBS and then were cut into 1 mm3 pieces. Containing Collagenase type I (2 mg/ml) and 0.25% (w/v) trypsin DMEM was used to digest the specimens for 1 hr at 37 °C using a Thermomixer at 1200 rpm (Eppendorf, Hamburg, Germany). After that, Ham’s F-12 media containing 10% FBS was added to stop the process of digestion, and then the digested tissue passed through a 100 µm cell strainer. Finally, after centrifugation, cells were collected for subsequent batch analysis.

Single-cell RNA-seq analysis

The raw single cell sequencing data was mapped and quantified with the 10×Genomics Inc software package CellRanger (v5.0.1) and the GRCh38 reference genome. Using the table of unique molecular identifiers produced by Cell Ranger, we identified droplets that contained cells using the call of functional droplets generated by Cell Ranger. After cell-containing droplets were identified, gene expression matrices were first filtered to remove cells with >5% mitochondrial genes, <250 or >8000 genes, and <500 UMI. Downstream analysis of Cellranger matrices was carried out using R (4.1.3) and the Seurat package (v 4.1.0, https://satijalab.org/seurat/). In total, 49,356 cells with an average of 2435 genes per cell were selected for ongoing analysis. Of these single cells, 24,721 were obtained from lesions ligament, including L2, L5, L6, and L8. The remainder were cells from healthy and included L1, L3, L4, and L10.

Single-cell trajectory analysis

We used the Monocle3 v. 2.8.0 R package to infer a hierarchical organization part of fibroblasts to organize these cells in pseudotime. We took these subpopulations from the Seurat data set from which we reperformed shared nearest neighbor clustering and differential expression analysis as described previously. We then selected the differentially expressed genes based on the fold-change expression for Monocle to order the cells using the DDRTree method and reverse graph embedding. We find gene co-expression modules according to the trend of gene expression and show the relationship between cells and modules in the form of heat maps.

Ligand-receptor interaction model

We used the ligand-receptor interaction database from the cellchat database to determine potential ligand-receptor interactions between fibroblast subpopulations and other cell types. The expression data were preprocessed for subsequent cell-cell communication analysis. The ligand or receptor that is overexpressed in a class of cells is first identified, and the gene expression data is then projected into a protein-protein interaction network. Whenever a ligand or receptor is overexpressed, the ligand-receptor interaction is recognized. CellphoneDB is based on the expression of a receptor in one cell type with a ligand in another cell type, resulting in rich receptor-ligand interactions between two cell populations. For the gene expressed by the cell population, the percentage of cells expressing the gene and the average gene expression were calculated. The gene was removed if it was expressed only in 10% or less of the cells in the population (the default value is 0.1).

Multiplex immunofluorescence (OPAL) staining

The tissue slices underwent deparaffinization and rehydration, and then heat-induced antigen retrieval was used. Following that, the paraffin slide containing the ligament tissue was subjected to continuous staining using the Opal Polaris Multiple-Color Manual IHC Kit (NEL861001KT), which involved incubating the tissue with different primary antibodies specific to various cell markers. We used different primary antibodies to simultaneously label pericyte (ACTA2), pericyte1 (MYH11), and pericyte2 (CD90) on the same tissue slide. APOE was used to label Fib.9 and TOP2A was used to label Fib.8. DAPI was used to stain cell nuclei. An automated staining system (BOND-RX, Leica Microsystems, Vista, CA) was used to perform the chromogen-based multiplex immunohistochemistry labeling. The Opal Polaris dyes were used to pair with these antibodies containing fluorophores for tyramide signal amplification (TSA) to enhance sensitivity. The following primary antibodies: FGF7 (160128, Abclonal), TGFβ1 (HA721143, HUABIO), APOE (YT0273, Immunoway), TOP2A (YM6914, Immunoway), ACTA2 (ab5694, Abcam), MYH11 (A4064, Abclonal), and CD90 (ab307736, Abcam) were used.

Spatial-transcriptome analysis

We used L1, and L8 samples for spatial transcriptome analysis. We performed the Seurat standard analysis. Machine prediction and marker gene methods were used to determine the cell type, and SCTransform was used for standard analysis of the data. First, we determine it is specific enough to find a topic profile (similar to a feature vector) for each cell type and then deconvolve the information by cell type and superimpose it onto the slice in pie chart form.

Deconvolution analysis

SCDC (v 0.0.0.9000) tool combined with single cell data was used to deconvolve bulk RNAseq data for cell type and content proportion analysis bulk RNAseq data. The bulk RNAseq data is also our own data. Histogram and heat map can be drawn for the predicted data, and the proportion of cell type content can be compared between groups. Then, we used GSVA (v 1.42.0) to calculate the same data and draw heat maps to aid the analysis.

Statistical analysis

A Nonparametric Wilcoxon rank sum test was used to analyze the differences between the two groups. All statistical analyses were performed in R or GraphPad Prism (version 5.0). Statistical significance was defined as *p<0.05, **p<0.01, ***p<0.001.

Results

Comprehensive scRNA-seq analyses reveal the major cell types in the human normal and degenerated ligament

To decipher cellular heterogeneity and construct the cell landscape for ligament degeneration, we performed single-cell transcriptomic profiling of cells from four healthy and four degenerated ACL samples (Figure 1—figure supplement 1). The characteristic pathological findings in the ACL of OA were identified by arthroscopy (Figure 1A and B). In the normal group, the ligaments are white stripes with a clear structure, and blood vessels can be seen on their surfaces. In the degenerated group, the ligaments were swollen, the structure was disorganized, and the color of the tissue was dark. From the images of knee MRI, we can also find that the signal of normal ACL on the T2 phase is uniform and the structure is clear, while the signal of the ACL in the degenerated group is mixed, and the clear structure cannot be distinguished (Figure 1—figure supplement 2A). After rigorous quality control, we obtained the transcriptomes of 49,356 cells (normal: 24635 degenerated: 24721) and then employed differential gene expression analysis to discern cluster-specific markers. Firstly, we used principal component analysis (PCA) to reduce the dimension, and then we adopted the uniform manifold approximation and projection (UMAP) method to conduct the next analysis. Unbiased clustering based on UMAP identified four cell classes and histrionic lineage-defining genes were detected, including endothelial cells: PECAM1, VWF, and PLVAP (3862, 7.83%), fibroblasts: DCN, LUM, and COL1A2 (32644, 66.14%), immune cells: PTPRC, SRGN, and CD163 (8021, 16.25%), and pericytes: ACTA2, MCAM, and RGS5 (4829, 9.78%) (Figure 1C–F and Figure 1—figure supplement 2B). According to Figure 1I, we could further confirm that the selected cell markers could well identify various cell types. UMAP plots and stacked bar plots in Figure 1G described the distributions of cells in the healthy and degenerated samples. We can conclude that, in general, fibroblasts are the main cell type in both normal and degenerated ligaments and the degenerated group had a higher immune cell ratio and lower pericyte ratio than the normal group. We next analyzed the number of differentially expressed genes (Degs) between healthy and degenerated ligament clusters. The results demonstrated that the fibroblast had the largest difference, implying that fibroblasts undergo significant changes during the degenerative process (Figure 1H). We also collected normal as well as degenerated ligament specimens for bulk sequencing and performed deconvolution analysis between bulk sequencing results and scRNA-seq results. The results illustrated that scRNA-seq samples were in good agreement with bulk sequencing samples (Figure 1J and K and Figure 1—figure supplement 2C). After that, through Degs analysis in the scRNA-seq results, we selected genes highly expressed in the degenerative group of the four-cell subsets (fibroblast, immune cell, endothelial cell, and pericyte). By examining the expression of these genes in bulk sequencing data, it can be found that these genes are also highly expressed in the diseased sample of bulk RNA-seq (Figure 1—figure supplement 2D–G). The results of the direct comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se also illustrated that scRNA-seq samples present high consistency with bulk RNA-seq samples (Figure 1L).

Figure 1 with 2 supplements see all
Download asset Open asset
Single-cell RNA-seq reveals major cell classes in human ligament.

(A and B) A photograph of typical normal (left) and degenerated (right) anterior cruciate ligaments (ACL) under arthroscopy. (C) UMAP visualization of all cell clusters in collected ligamental specimens. (D and E) UMAP visualization of the donor origins in normal/diseased samples. (F) Heatmap of selected marker genes in each cell cluster. (G) The percentages of the identified cell classes in normal/diseased ligament. (H) Number of differentially expressed genes (DEGs) in each cell type of normal/diseased status. (I) Feature plots of expression distribution for selected cluster-specific genes. Brighter colors indicate higher expression levels. (J) Gene expression profiles from bulk samples (n=6) and in single cell samples (n=8) were averaged and plotted on X and Y axes, respectively. Red lines indicate linear model fit and the diagonal. (K) Correlation heatmap shows Pearson’s correlation between all bulk and in single cell samples. (L) Consistency analysis of Degs between normal and diseased groups in single cell and bulk RNA sequencing data.

Characterization of fibroblast subpopulations in healthy and degenerative ACL

Because fibroblasts are known to play a vital role in ECM homeostasis and the degeneration of the ligament and undergo remarkable changes during the process of ligament degeneration, we next conducted further analysis of fibroblasts in normal and degenerated ligaments. We subclustered fibroblasts and identified 10 subsets (Fib.1-Fib.10; Figure 2A). The distribution shown in Figure 2B indicated that the batch effects are resolved well. As illustrated in Figure 2C, we obtained the cell proportions of the fibroblast subpopulations between normal and diseased groups. We can conclude that the ratio of fibroblast subclusters varied greatly between these two groups. The ratio of Fib.1, Fib.2, Fib.8, and Fib.9 in the diseased group were higher than that in the healthy group, and the ratio of Fib.3, Fib.4, and Fib.7 in the healthy group were higher compared to the diseased group. The proportion of remaining fibroblast subclusters did not differ significantly between the two groups.

Figure 2
Download asset Open asset
Characterization of fibroblast subclusters across different statuses in human ACL.

(A) UMAP visualization of the subclusters of fibroblasts. (B) UMAP visualization of the distribution of fibroblast subclusters at different samples. (C) The proportions of 10 fibroblast subpopulations in normal and diseased ligaments. (D) A cell-level heatmap reveals the normalized expression of DEGs for each fibroblast cluster defined in. (E) Volcano plot showing the DEGs between two statues of ligamental fibroblasts. The x axis represents highly expressed genes in normal cells, and the y axis represents highly expressed genes in diseased cells. (F) GO and KEGG enrichment analysis of DEGs in ligamental fibroblasts between normal and degenerated states. (G) GSEA enrichment plots for representative signaling pathways upregulated in fibroblasts of diseased samples, compared with normal samples. (H) Gene-gene interaction networks between DEGs in ligamental fibroblasts of normal group and fibroblasts of diseased group.

According to the analysis of Degs among these fibroblast subgroups, we tried to reveal the identity of each subpopulation (Figure 2D). Combined with the analysis of the proportion of cell subclusters, we inferred that Fib.1 and Fib.2 are ECM remodeling-associated fibroblasts as they highly expressed the genes related to ECM components and ECM decomposition such as COL1A1, CTHRC1, COL3A1, FBLN1, MMP2, MMP14, and TPPP3. From the perspective of ECM synthesis, Fib.1 and Fib.2 are slightly different. Fib.1 highly expressed the collagen-related genes, but Fib.2 highly expressed the genes related to the synthesis of connexins. Fib.3 expressing metallothionein family genes such as MT1X, MT1E, and MT1M is a population of homeostasis-associated fibroblasts with defensive functions. Fib.4, expressing genes associated with the formation and organization of the ligamental ECM such as ANGPTL7, MYOC, CILP2, and THBS1 is a group of resident fibroblasts that function normally. Fib.5 is a cluster of structural fibroblasts within the ligament as they highly expressed multiple collagen-related genes such as COL3A1, COL6A1, COL5A1, COL1A1, and TGFVBI. Fib.6 is a cluster of chondrocytes as they highly expressed cartilage-related genes such as COL2A1, FMOD, CHAD, ACAN, and CILP2. We cannot define Fib.7 well because they expressed Degs without specific tags. Fib.8 is a group of cycling cells as they highly expressed cell cycle-related genes such as CENPF, TOP2A, and MLI67. Fib.9 is a group of inflammation-related fibroblasts as they highly expressed CXCL14, CXCL12, C3, and C7. Fib.10 may have the function of damage repair as they highly expressed anti-inflammation and anti-ECM decomposing associated genes such as PRG4, DEFB1, TIMP3, and HBEGF. To verify the accuracy of the classification of fibroblast subclusters, the top 50 Degs of each subpopulation were deconvolved on the bulk RNA-seq data by the GSVA method. The results illustrated that fibroblast 3,4 scored higher in the normal samples and fibroblast 1,2,9 scored higher in the diseased samples, which were consistent with the changes in proportion of scRNA-seq.

For further in-depth investigation of changes in ligamental fibroblasts during the degeneration process, we performed the Degs analysis in the whole fibroblast population between the normal and degenerated ligaments. We observed genes upregulated in the degenerated group such as CFD, SPP1, COL1A2, and COL3A1, which were related to scar healing and ECM remolding (Figure 2E). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were conducted for further functional interpretation of these Degs. The results illustrated that the upregulated genes in the diseased group were enriched for the terms ‘extracellular region’, ‘TGF-beta signaling pathway’ and ‘Complement and coagulation cascades’ (Figure 2F). Gene Set Enrichment Analysis (GSEA) also demonstrated that ECM-associated signaling, ECM binding and ECM structural constituent, were activated in the degenerated group, which implied that the changes of ECM play an important role in the ligament degeneration (Figure 2G). And then, by building gene interaction networks, we identified that genes upregulated in degenerated groups were concentrated in two parts, ECM-receptor interaction and leukocyte trans-endothelial migration, which once again proved the significant function of ECM remolding and immune inflammation in the disease process (Figure 2H).

Dynamic transcriptional changes in ligamental degeneration

To gain insight into the cellular progression of fibroblast subclusters during the disease process, we performed cell trajectory analysis, revealing two significant routes in the process of degeneration. Five fibroblast subclusters, including Fib.2, Fib.4, Fib.5, Fib.9, and Fib.10, were involved in the reconstruction of the disease trajectories using Monocle 3, an algorithm for the reconstruction of lineage programs based on similarity at the transcriptional level (Cao et al., 2019). We set Fib.4 and Fib.5 as the starting point of the trajectories because they are structural fibroblasts with high expression of ECM component genes and represent normal ligamental functions well, and then computed pseudotime for cells along the inferred developmental axis (Figure 3A and B). More specifically, Fib.4 and Fib.5 were predicted to transform into two distinct cell fates, including the cell fate1, which includes Fib.2 and Fib.9, and the cell fate2, which includes Fib.10. With this in mind, we tried to explore the genic dynamics that distinguished these two cell fates. The expression profile of cell fate1 showed increasingly high expression of genes (CXCL12, MMP2, MMP14, C7) related to ‘leukocyte trans-endothelial migration’, ‘cellular response to chemokine’, and ‘inflammatory response pathway’. Along with the cell fate2, we observed gradually high expression of ECM organization and skeletal muscle tissue development-related genes (FOS, EGR1, FGF10, FN1; Figure 3C–H).

Figure 3
Download asset Open asset
Evolution trajectory and transcriptional fluctuation during ligamental degeneration progression.

(A and B) UMAP visualization of fibroblast 2, 4, 5, 9, 10. Developmental pseudotime for cells present along the trajectory inferred by Monocle 3, with cell fate1 and cell fate2 branches coming from fibroblast 4. (C) Heatmap showing the expression changes of the highly variable genes along the cell fate1 in normal and degenerated groups. (D and E) Representative gene expression levels along cell fate1 trajectory of normal and diseased statuses. (F) Heatmap showing the expression changes of the highly variable genes along the cell fate2 in normal and degenerated groups. (G and H) Representative gene expression levels along cell fate2 trajectory of normal and diseased statuses. (I) Heatmap showing the scaled mean expression of modules of coregulated genes grouped by Louvain community analysis across the subclusters. (J and K) Enrichment analysis results of Module 3 and 5.

To investigate gene expression dynamics along the trajectories, we group genes that varied between cell clusters into Module 10 using Louvain community analysis. We can observe the aggregated expression of each module in Figure 3I. We identified that the expression of genes in module 3 was gradually elevated along the cell fate1 and gradually declined along the cell fate2, which were enriched for genes related to rheumatoid arthritis, antigen processing and presentation, and complement activation. In contrast, the expression of wound repair and ECM assembly-related genes were gradually elevated along the cell fate2 but declined along the cell fate1, such as ADAM12, HS3ST3A1, FN1, and SERPINE2 in module 5 (Figure 3I–K). According to these pseudotime analysis results, we inferred that there are two opposite dynamic trajectories in the process of ligamental degeneration: one is the progressive damage process of complement inflammation leading to the continuous progression of the disease, and the other is the process of repairing the damage to delay or even reverse the disease process.

Characterization of stromal cells in coordinating ligament microenvironment during degeneration progression

Besides fibroblasts, various endothelial cells (ECs) and immune cells also play significant roles in the occurrence and development of ligament degeneration. In our ligament samples, we identified seven clusters of blood-vessel-derived cells, including five clusters of ECs (EC1-5) and two clusters of pericytes (pericyte1 and pericyte2). Figure 4A shown the distribution of EC1 expressing a high level of ACKR1 is a population of venous ECs; EC2 highly expressed COL4A1, COL4A2, H19 was related to anti-angiogenesis; EC3 highly expressed genes associated with ECM organization and we inferred this cluster related to damage repair. EC4 expressing a high level of CCL2, CCL8, STEAP4, and SYNPO2 is a cluster of ECs related to inflammatory chemotaxis. EC5 is a cluster of lymphatic ECs. As for pericytes, we identified two populations that expressed pericyte markers ACTA2, and MCAM. Pericyte1 highly expressed myocyte-related genes such as MYH11, MUSTN1, and LBH, which implied that this cluster may have myocyte-like properties. Pericyte2 highly expressed fibroblast markers such as COL1A1 and COL3A1, which implied that this cluster may have fibroblast-like properties (Figure 4G). We next compared the proportions of these 7 clusters between normal and diseased ligaments (Figure 4Cand D). The results illustrated that most ECs subsets were elevated in the diseased group, implying increased angiogenesis during ligamental degeneration. The ratio of pericyte1 in the diseased group was higher than that in the healthy group, but the ratio of pericyte2 was lower. An increased ratio of pericyte1 may change the tissue’s biological properties, as ligament was dominant by fibroblast. To gain more biological insights underlying ECs, we performed the Degs analysis between normal and degenerated groups (Figure 4E). Combined with the result of enrichment analysis and GSEA analysis, we identified that ligamental ECs upregulated inflammation and complement-related genes, such as C1QA, PLA2G2A, and ECM-related genes such as PRG4, VIM, and POSTN during degeneration procession (Figure 4H). The terms ‘inflammation mediated by chemokine and cytokine signaling pathway’ and ‘TGF-beta signaling pathway’ activated in the diseased group also implied that inflammation and ECM remodeling have significant roles in the disease process (Figure 4F). By immunofluorescence staining, we demonstrated the presence of two pericyte subsets in normal and diseased samples. In the normal ligament, there were more pericyte2 than pericyte1 and compared with the normal group, pericyte1 was increased in the diseased group (Figure 4I and J and Figure 4—figure supplement 1).

Figure 4 with 1 supplement see all
Download asset Open asset
Identification of blood vessel derived cell and immune cell subclusters in human ligament.

(A) UMAP visualization of the subclusters of endothelial cell and pericyte. (B) UMAP visualization of the distribution of endothelial cell and pericyte subclusters at different sample statuses. (C and D) Summarized subpopulations of endothelial cell and pericyte percentage changes. (E) Volcano plots displaying the DEGs in endothelial cell between normal group and diseased group. Each dot represented one gene. Red dots, differentially up-regulated genes; blue dots, differentially down-regulated genes; gray dot, non-differentially expressed genes. (F) GSEA enrichment plots for representative signaling pathways upregulated in endothelial cell of diseased samples, compared with normal samples. (G) Violin plots showing representative marker genes associated with different types of pericyte expressed in pericyte1 and pericyte2. (H) GO and KEGG enrichment analysis of DEGs between normal and diseased endothelial cells. (I and J) Immunofluorescence staining of pericyte related markers in normal and diseased groups. (K) UMAP visualization of the subclusters of immune cell and the distribution of immune cell subclusters at different sample statuses. (L) The proportion of each subcluster of immune cells in the lesioned and normal ligament. (M) Summarized subpopulations of immune cell percentage changes. (N) Violin plots showing representative genes of macrophages between normal and degenerated states. (O) GO and KEGG enrichment analysis of DEGs between normal and diseased macrophages.

As an indispensable cellular component in the ligament, we identified six subpopulations of immune cells, including basophils, CD8+ NKT-like cells, macrophages, memory CD4+ T cells, memory CD8+ T cells, and myeloid dendritic cells (Figure 4K). From Figure 4L and M, we can conclude that macrophages are the most abundant immune cells in both normal and diseased groups and compared with the normal group, only the proportion of macrophages was increased in the degenerated group. We explored the gene profiling change of macrophages in the process of degeneration. As is shown in the violin plot Figure 4N, we suggested that ECM remolding-related genes have higher expression levels in the diseased group, such as PRG4, FN1, and HTRA1. The enrichment analysis also implied that ECM organization and complement and coagulation cascades were activated in the degenerated ligament, which was consistent with the results of ECs analysis (Figure 4O).

Putative signaling network for the intercellular crosstalk regulating the microenvironment homeostasis during ligamental degeneration

Elucidating the explicit interaction among fibroblasts, ECs, and immune cells in the ligamental microenvironment will shed light on the mechanisms of the pathogenesis of ligamental degeneration. CellphoneDB and CellChat analysis were used to investigate the signaling network among the main cell clusters in the normal and diseased ligaments. The heatmap illustrated the intensity of interactions among different cell clusters. From the results, we can find that in addition to interactions within fibroblast subclusters, cell-cell cross-talks were mainly between fibroblast subpopulations and ECs and between fibroblast subpopulations and macrophages. At the same time, we can find that fibroblast subsets dominated by the diseased group had stronger interaction associations with other cells, which implied that there are stronger and more complex interactions in degenerated ligaments (Figure 5A).

Figure 5 with 1 supplement see all
Download asset Open asset
Cell–cell crosstalk during ligamental degeneration progression.

(A) Heatmap depicting the significant interactions among the identified subclusters of fibroblast, immune cell, and endothelial cell of normal and degenerated ligaments. (B and C) Circle plots showing the inferred TGF-β and FGF signaling networks. (D and E) Violin plots showing ligand-receptor interactions related to TGF-β and FGF signaling pathways among subclusters of fibroblast, endothelial cell and immune cell. (F) Immunofluorescent staining of TGFβ1, TOP2A, APOE, and FGF7 in normal and diseased groups.

To determine the important factors, we further analyzed the intercellular signaling networks of FGF and TGFB. Interestingly, Fib.9 was involved in FGF signaling, both autocrine and paracrine (Figure 5B). Fibroblasts were the leading receiver of FGF signals, as expected. Through the violin plot, we identified that Fib.9 acts on FGFR1 located in each fibroblast subpopulation by expressing FGF7 and FGF signaling pathways, which can promote pathological proliferation of cells and participate in scar repair (Figure 5D). We inferred that the FGF signaling network may disrupt the homeostasis of fibroblasts in the ligaments and promote the progression of ligamental degeneration. Moreover, the TGF-β pathway was involved in many cell-cell interactions among fibroblast subpopulations and macrophages via TGFB1-TGFBR1 or TGFB1-TGFBR2 (Figure 5C and E). As shown above, fibroblasts of the degenerated group highly expressed PRG4, COL3A1, and SSP1, and these genes are all the important downstream targets of TGF-β, which participate in the process of ECM remolding. Hence, we suggested that TGFB1 is a key gene and Fib.8 and macrophages promote the development of the disease through this molecule. According to the Degs, we selected APOE as the marker gene of Fib.9 and TOP2A as the marker gene of Fib.8. The results of immunofluorescent staining illustrated that the number of Fib.8 and Fib.9 in the disease group were more than those in the normal group, and Fib.8 highly expressed TGFB1 and Fib.9 highly expressed FGF7 in the disease group, which validated the above findings (Figure 5F and Figure 5—figure supplement 1).

SpRNA-seq of ligament deciphering the spatial interactions of cell subclusters

For further insight, and multi-angle interpretation of the cell composition changes that occurred in ligamental degeneration, we performed spRNA-seq on normal and lesioned ligaments. Firstly, by using the feature marker transfer method, we mapped the individual subsets identified by single-cell sequencing into spatial transcriptome data. The heatmap illustrated that Fib.4, as a representative of fibroblasts in the normal ligament tissue, was widely distributed in L1 and the amount was significantly higher than that in L8. The number of Fib.9, as a representative of fibroblasts in the degenerated ligament, in L8 was higher than that in L1 (Figure 6A–D). We also examined the immune cells and ECs. As expected, there were more immune cells and ECs in L8 than in L1, implying immune infiltration and vascular hyperplasia were presented in the degenerative specimens (Figure 6E–J). In addition, we also examined the FGF and TGFB signaling pathways related genes expression in spRNA-seq. From the results, we can find that TGFB1, TGFBR1, FGF7, and FGFR1 were highly expressed in the diseased sample (Figure 6—figure supplement 1). These results are consistent with those obtained by scRNA-seq. After that, we performed spotlight analysis in these two groups. All cell subsets are rendered in four colors, with red being the fibroblast subsets specific to the disease group. Finally, we hide the areas that are not red in both samples. We can find that the red region in L8 was widely distributed, and the area was significantly higher than that in L1. After careful observation, we found that almost all the red dots had yellow and black areas (Figure 6K and L). This indicates that some fibroblast subpopulations specific to the diseased group are adjacent to endothelial and immune cells in the degenerative state. The proximity of different cells in space makes it possible for them to interact with each other. Therefore, this further provided strong evidence for the results of our single-cell interaction analysis.

Figure 6 with 1 supplement see all
Download asset Open asset
Spatial transcriptome sequencing deciphers the microenvironment changes during ligamental degeneration progression.

(A–D) Spatial heatmaps showing the representative normal and degenerated fibroblast subcluster distribution in L1 and L8 samples. (E and F) Spatial heatmaps showing the immune cells distribution in normal and degenerated ligamental samples. (G) Quantitative analysis of immune cells in these two types of tissues. (H and I) Spatial heatmaps showing the endothelial cells distribution in normal and degenerated ligamental samples. (J) Quantitative analysis of endothelial cells in these two types of tissues. (K and L) Spot light showing fibroblast subpopulation unique to diseased samples proportion in L1 and L8 samples.

Discussion

ACL degeneration can contribute to cartilage injury and even OA onset and progression and may put a great burden on the normal life of many patients (Georgoulis et al., 2010; Hasegawa et al., 2012). However, the mechanisms of this disease are not well characterized and treatments to prevent or treat ACL degeneration are scarce and not effective (Shane Anderson and Loeser, 2010; Tozer and Duprez, 2005). Healthy and degenerated ACL tissues include multiple cell subpopulations with diverse genetic and phenotypic characteristics. How this heterogeneity emerges in development degeneration remains unclear. Herein, we built a single-cell atlas of normal and degenerated human ACL and explored the characteristics and key regulatory pathways of distinct fibroblast subtypes. These findings will help us understand the pathogenesis of ACL degeneration in-depth, and provide potential targets for clinical therapies for this disease.

Fibroblasts are increasingly considered dominating cell types in ligaments and central mediators of ligamental degeneration, and here, we identified 10 fibroblast subpopulations in normal and degenerated ACL samples by using scRNA-seq. Further cluster analysis profiled characteristics of each subcluster and combined with the cell proportion analysis, we found that pro-inflammation and ECM remodeling-related fibroblast subgroups were significantly increased in the degenerative ligamental tissues, which implied that inflammation and ECM remodeling are key events in the disease process. Through comparing the Degs of fibroblasts between degenerative and healthy conditions, we found that ECM-related genes upregulated in the diseased group, such as FN1, COL3A1, PRG4, and SPP1. The results of enrichment analysis and GSEA also illustrated that ECM-related pathways were activated in the degenerated group, which suggested that the modules of ECM play an important role in the process of ligamental degeneration. According to the gene interaction analysis, two gene modules were activated in the diseased group, leukocyte trans-endothelial migration and ECM-receptor interaction. These findings were consistent with previous studies which suggested that ligament inflammation and ECM changes contribute to ACL degeneration (Busch et al., 2013; Hasegawa et al., 2012). It is due to the changes and remodeling of the ECM in the degenerative ligament tissue that the biomechanical changes of the ligament are caused and eventually may result in the rupture of the ligament. From the pseudotime analysis, we identified two cell fates in the ligamental degeneration. Cell fate 1 represents a progressive process of the disease characterized by inflammatory damage and extracellular matrix degradation. Cell fate 2 implies a chronic repair process for the disease characterized by ECM remodeling, ligament development, and damage repair.

Immune cells are closely related to ligamental degeneration (Kim-Wang et al., 2021). We analyzed the heterogeneity of immune cells in ACL using single-cell RNA sequencing. The results illustrated that macrophage-dominating immune cells in ACL and inflammation-related genes PLA2G2A and ECM-related genes PRG4, FN1, and HTRA1 were upregulated in the macrophages of the diseased group. Enrichment analysis also revealed that complement and ECM-related pathways were activated in the degenerated group. This implied that macrophages highly expressed ECM and inflammation-associated genes in the ligamental degeneration progression. So, we suggested that macrophages may contribute to ligamental degeneration. We also identified two types of pericytes in these tissues. As we all know, pericytes in the tissue have some stem cell characteristics and can differentiate into fibroblasts, and chondrocytes (Armulik et al., 2011; Smyth et al., 2018). Pericyte 1 highly expressed stem cell-related genes MCAM and MUSTN1 and have the characteristics of myofibroblasts, and pericyte 2 tends to have the characteristics of fibroblasts. Pericyte 1 was increased and pericyte 2 decreased in the degenerated group, which may change the properties of ligamental cells and contribute to ligamental degeneration.

Cells can communicate via ligand-receptor interactions (Kumar et al., 2018), so targeting cell-cell interactions is frequently utilized in clinical treatment. The results of CellphoneDB analysis demonstrated that the high content of fibroblasts and immune cells in the degenerative group had significantly higher interaction intensity, such as fibroblast 1, 2, and 8, which means that the interaction between cells in the ligament was enhanced in the degenerative state. To identify the key factors regulating the disease process, CellChat analysis was used to dissect the intercellular crosstalk based on the signaling network in the human ACL. We found that the FGF signaling pathway and TGF-β signaling pathway were involved in the crosstalk network. TGFβ1 contributes to osteophyte formation by inducing endochondral ossification, is thought to be closely related to the onset of OA (Murata et al., 2019), and can be used to induce chondrogenic differentiation of ligament-derived stem cells in vitro (Schwarz et al., 2019). It has been reported that in myocardial fibrosis, liver fibrosis, glaucoma, and other diseases, the TGF-β signaling pathway is also involved in angiogenesis, ECM remodeling, pathological scar healing (Ilieș et al., 2021; Lee and Massagué, 2022; Ferrão et al., 2018; Penn et al., 2012). FGFs can induce the fibroblast to myofibroblast differentiation, promote tissue repair by regulating cell proliferation, survival, and angiogenesis and contribute to ECM remodeling (Kendall and Feghali-Bostwick, 2014; Ma et al., 2017). ECM remodeling by overexpression, degradation, and cross-linking of ECM proteins in the lesioned tissue are direct factors leading to fibrosis tissues (Sun et al., 2022). In this study, we observed that fibroblast9, inflammation-related fibroblast, highly expressed FGF7 and can act on FGFR1 mainly distributed in the fibroblast subgroups. Immune cells especially macrophages highly expressed TGFB1 and can act on TGFBR1 and TGFBR2 mainly distributed in the fibroblast and endothelial subgroups. These findings suggested that FGF and TGF-β signaling pathways may induce ECM remolding in the process of ligamental degeneration and FGF7-FGFR1 and TGFB1-TGFBR2 may be potential targets for the treatment of this disease. SpRNA-seq was used to detect the spatial information underlying the normal and degenerated ligaments and validate the findings acquired from the scRNA-seq. Through the results, we observed that the disease group specific fibroblasts and immune cells recognized by single cells were more numerous and more widely distributed in the disease group specimens. We also got some spatial information, that fibroblasts in the disease group were spatially closer to immune and endothelial cells, which was more conducive to cell interaction.

Conclusions

In conclusion, our study described the cell atlas of the human ligament, providing a valuable resource for further investigation of ACL homeostasis and the pathogenesis of ligamental degeneration. The cellular heterogeneity and signaling network we uncovered help to increase the understanding of the human ACL at a single-cell level and provide crucial clues for establishing new diagnostic and therapeutic strategies for this disease in the future.

Data availability

Data are available in a public, open-access repository. The single-cell RNA-seq data and cluster annotations are available at GSA for human (https://ngdc.cncb.ac.cn/gsa-human/) with the accession number PRJCA014157.

The following data sets were generated
    1. Fu W
    (2023) CNCB BioProject
    ID PRJCA014157. The single-cell RNA-seq data and cluster annotations are available at GSA for human.
    https://ngdc.cncb.ac.cn/bioproject/browse/PRJCA014157

References

    1. Fleming BC
    (2003)
    Biomechanics of the anterior cruciate ligament
    The Journal of Orthopaedic and Sports Physical Therapy 33:A13–A15.
    1. Frank CB
    (2004)
    Ligament structure, physiology and function
    Journal of Musculoskeletal & Neuronal Interactions 4:199–201.
    1. Penn JW
    2. Grobbelaar AO
    3. Rolfe KJ
    (2012)
    The role of the TGF-β family in wound healing, burns and scarring: a review
    International Journal of Burns and Trauma 2:18–28.

Decision letter

  1. Carlos Isales
    Senior and Reviewing Editor; Augusta University, United States

Our editorial process produces two outputs: i) public reviews designed to be posted alongside the preprint for the benefit of readers; ii) feedback on the manuscript for the authors, including requests for revisions, shown below. We also include an acceptance summary that explains what the editors found interesting or important about the work.

Decision letter after peer review:

Thank you for submitting your article "A single-cell atlas depicting the cellular and molecular features in human ligamental degeneration: a single cell combined spatial transcriptomics study" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Carlos Isales as the Senior Editor. The reviewers have opted to remain anonymous.

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The initial findings provided by the scRNA-seq technique need to be validated by in vitro and in vivo experiments.

2) The "Introduction" section could be improved substantially by removing irrelevant and redundant information. This section should focus on the biology of ACL, particularly current knowledge about the cellular and molecular mechanisms that maintain ACL homeostasis or cause degeneration.

3) The legend of Figure S1 D-G and the corresponding text needs to be revised to increase clarity. The current layouts are confusing, and those violin plots seem to present scRNA-seq data. There wasn't any comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se.

In addition, the authors need to address other comments raised by 3 reviewers.

Reviewer #1 (Recommendations for the authors):

In this study, the authors performed single-cell RNA sequencing using normal and diseased human anterior cruciate ligaments. However, the study is very preliminary, which needs to be confirmed by animal studies and in vitro studies.

1) The authors identified two types of pericyte cells with different characteristics. If this can be reproduced by animal study, it will strengthen the findings reported by single-cell RNA sequencing assay.

2) Roles of TGF-β and FGF signaling in the anterior cruciate ligament degeneration need to be further verified by animal studies and by in vitro studies.

Reviewer #2 (Recommendations for the authors):

1) Some of the descriptions in the text are inappropriate. Such as L25 "two stages of the degenerative process". L224: "obviously higher than that in the healthy group".

2) In figure3, the author only kept 4 clusters for the pseudotime analysis, what's the reason for excluding the other subpopulations?

3) The grammar throughout the manuscript needs significant improvement.

Reviewer #3 (Recommendations for the authors):

Concerns:

1. The "Introduction" section could be improved substantially by removing irrelevant and redundant information. This section should focus on the biology of ACL, particularly current knowledge about the cellular and molecular mechanisms that maintain ACL homeostasis or cause degeneration.

2. The authors should check and review extensively for improvements to the use of English.

3. The labels in Figure S1C-G used to define the normal and disease groups and cell clusters (fibroblasts, immune cells, pericytes, and endothelial cells) are too small to read. In Figure 1F, the genes listed on the left didn't match the corresponding heatmap, and their color codes were not defined.

4. Page 10, lines 204-206 and Figure S1 D-G: these figures, figure legends, and corresponding text need to be revised to increase clarity. The current layouts are confusing, and those violin plots seem to present scRNA-seq data. There wasn't any comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se.

5. Page 10, lines 208-209: Please consider revising it as "Characterization of fibroblast subpopulations in healthy and degenerative ACL". The rationale and criteria used for classifying subclusters of the ACL fibroblasts are ambiguous (Figures 2&3). This classification lacks the support of functional validation, therefore dampening its authenticity and significance.

6. The Data shown in Figure 2 F-H are computational models; they need to be validated by additional experiments. This applies to other claimed findings as well.

7. Considering the effect of biological variables and age on this study, information about the donors' gender and age should be provided.

[Editors' note: further revisions were suggested prior to acceptance, as described below.]

Thank you for resubmitting your work entitled "A single-cell atlas depicting the cellular and molecular features in human ligamental degeneration: a single cell combined spatial transcriptomics study" for further consideration by eLife. Your revised article has been evaluated by Carlos Isales (Senior Editor) and a Reviewing Editor.

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

1. The authors are suggested to run the proofreading carefully since many typos are found in the text. Below is a list of some, but not all.

P3 line 48: consider changing " the musculoskeletal" to "the musculoskeletal system".

P3 line 59: Please consider revising "ACL aging and degradation" to "ACL aging and degeneration". Also, "degradation" has been heavily used in the text. Consider replacing it with "degeneration" wherever it is appropriate.

P4 line 85: replace "synergize" with "synergy"; P10 line 204: please replace "resolve" with "reveal"?

2. The color codes for Figure 2B are confusing, please consider changing the color codes for L3 to L5 to increase visibility.

3. Markers used for pericyte immunostaining in Figure 4G were incorrect! ACTA2, MYH11, and CD90 (Thy1) are not specific markers for pericytes. ACTA2 (also called sSMA) and MYH11 (myosin11) are generally considered markers for myofibroblasts and vascular smooth muscle, not pericytes. Also, most of the mesenchymal progenitors express CD19. So pericyte classification in this ms is confusing and not supported by the immunostaining data.

4. The data in Figure 5F doesn't convincingly demonstrate there is an increase in the expression of FGF7 and TGFβ1 (where and which types of cells) in the specimens from the diseased groups. In addition, these data need to be quantified, and the "n" number should be provided.

https://doi.org/10.7554/eLife.85700.sa1

Author response

Essential revisions:

1) The initial findings provided by the scRNA-seq technique need to be validated by in vitro and in vivo experiments.

Thanks for the comments and suggestions. We followed them and mainly validated the changes of fibroblasts during ligamental degeneration and the changes of two pericyte subsets in normal and degenerative specimens by immunohistochemical staining, and immunofluorescence staining. In addition, two disease-related pathways identified by the scRNA-seq: TGF-β and FGF signaling were verified by immunohistochemical staining and validated at the spatial transcriptome level.

2) The "Introduction" section could be improved substantially by removing irrelevant and redundant information. This section should focus on the biology of ACL, particularly current knowledge about the cellular and molecular mechanisms that maintain ACL homeostasis or cause degeneration.

Thank you for your suggestions. We have reworded the “Introduction” section according to your requirement. The revised “Introduction” section focuses on the current knowledge about the cellular and molecular mechanisms that maintain ACL homeostasis or cause degeneration.

3) The legend of Figure S1 D-G and the corresponding text needs to be revised to increase clarity. The current layouts are confusing, and those violin plots seem to present scRNA-seq data. There wasn't any comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se.

Sorry for the inconvenience. We have now re-written the legend of Figure S1 D-G and the corresponding text. These violin plots do not represent scRNA-seq data. Here, we make a detailed explanation as follows. We first selected genes that were significant highly expressed in the diseased group than in the normal group in the four clusters (fibroblast, endothelial cell, immune cell, and pericyte) of our scRNA-seq data, and then we detected these genes in the bulk RNA-seq data to verify whether these genes were also highly expressed in the diseased group of the bulk RNA-seq. Thus, the consistency of our scRNA-seq data with the bulk RNA-seq data was demonstrated by these violin plots. In addition, according to your suggestions, we added the direct comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se.

Reviewer #1 (Recommendations for the authors):

In this study, the authors performed single-cell RNA sequencing using normal and diseased human anterior cruciate ligaments. However, the study is very preliminary, which needs to be confirmed by animal studies and in vitro studies.

1) The authors identified two types of pericyte cells with different characteristics. If this can be reproduced by animal study, it will strengthen the findings reported by single-cell RNA sequencing assay.

Thanks. According to your suggestions, we selected the marker genes of two types of pericytes and performed immunofluorescence staining in normal and diseased samples, which confirmed the findings reported by scRNA-seq.

2) Roles of TGF-β and FGF signaling in the anterior cruciate ligament degeneration need to be further verified by animal studies and by in vitro studies.

Thanks for your suggestions. We verified the expression of key genes in these two signaling pathways: FGF7 and TGFB1 in the normal and diseased groups by immunohistochemical staining, and also detected these two pathways at the spatial transcriptome level. As a single-cell atlas type article, we aim to construct single-cell maps of both normal and degenerative ligaments to provide tools and resources for future research in ligament-related diseases. So relatively little experimental verification has been done.

Reviewer #2 (Recommendations for the authors):

1) Some of the descriptions in the text are inappropriate. Such as L25 "two stages of the degenerative process". L224: "obviously higher than that in the healthy group".

Sorry for that. According to your suggestions, we have reworded the inappropriate descriptions in the manuscript.

2) In figure3, the author only kept 4 clusters for the pseudotime analysis, what's the reason for excluding the other subpopulations?

Thanks for your great questions. We kept 5 clusters (fibroblast2,4,5,9,10) for the pseodotime analysis. According to the Degs and the definition of different fibroblast subclusters, these five cell populations are the focus of our attention and are closely related to the disease progression. Besides, based on the current understanding of ligamental degeneration, we hypothesized that the process of ligamental degeneration could be reflected in these cell clusters by the pseodotime analysis. In addition, the reasons for excluding the other subpopulations are also described. Fibroblast8: a subpopulation that are in the cell cycle is not associated with disease progression. Fibroblast7: a subpopulation that has not yet been defined. Fibroblast6: a cluster of chondrocytes and not associated with ligamental degeneration. Fibroblast3: high expression of metalloenzyme related genes and low correlation with the disease. Fibroblast1: a subpopulation associated with ECM remodeling. It was initially included in the pseodotime analysis, but the results did not fit the biological process, so it was also excluded.

3) The grammar throughout the manuscript needs significant improvement.

Thanks for your comments and suggestions. We carefully checked and revised the grammar in the manuscript. Besides, we also asked for help from the native speaker on the grammar and syntax issues.

Reviewer #3 (Recommendations for the authors):

Concerns:

1. The "Introduction" section could be improved substantially by removing irrelevant and redundant information. This section should focus on the biology of ACL, particularly current knowledge about the cellular and molecular mechanisms that maintain ACL homeostasis or cause degeneration.

Thanks. According to your suggestions, we have reworded the “Introduction” section according to your requirement. The revised “Introduction” section focuses on the current knowledge about the cellular and molecular mechanisms that maintain ACL homeostasis or cause degeneration.

2. The authors should check and review extensively for improvements to the use of English.

Thanks for your reminding. We have carefully checked and reviewed the use of English of the manuscript and ask for help from the native speaker. Besides, we have extensively and rigorously revised the English use in the manuscript.

3. The labels in Figure S1C-G used to define the normal and disease groups and cell clusters (fibroblasts, immune cells, pericytes, and endothelial cells) are too small to read. In Figure 1F, the genes listed on the left didn't match the corresponding heatmap, and their color codes were not defined.

Thanks for your suggestions. We adjusted the labels in Figure S1C-G and redrew the Figure 1F.

4. Page 10, lines 204-206 and Figure S1 D-G: these figures, figure legends, and corresponding text need to be revised to increase clarity. The current layouts are confusing, and those violin plots seem to present scRNA-seq data. There wasn't any comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se.

Sorry for the inconvenience. We have now re-written the legend of Figure S1 D-G and the corresponding text. These violin plots do not represent scRNA-seq data. Here, we make a detailed explanation as follows. We first selected genes that were significant highly expressed in the diseased group than in the normal group in the four clusters (fibroblast, endothelial cell, immune cell, and pericyte) of our scRNA-seq data, and then we detected these genes in the bulk RNA-seq data to verify whether these genes were also highly expressed in the diseased group of the bulk RNA-seq. Thus, the consistency of our scRNA-seq data with the bulk RNA-seq data was demonstrated by these violin plots.

In addition, according to your suggestions, we added the direct comparison of gene expression between the bulk RNA-seq and scRNA-seq data per se and a GSVA analysis. Through the results of direct comparison and GSVA analysis (figure 1L and supplementary figure 1C), we further proved the reliability of our scRNA-seq results.

5. Page 10, lines 208-209: Please consider revising it as "Characterization of fibroblast subpopulations in healthy and degenerative ACL". The rationale and criteria used for classifying subclusters of the ACL fibroblasts are ambiguous (Figures 2&3). This classification lacks the support of functional validation, therefore dampening its authenticity and significance.

Thank you for your suggestions. For lines 208-209, We have revised it as “Characterization of fibroblast subpopulations in healthy and degenerative ACL”. We believed that our classification of fibroblast subclusters are objective and rational. We defined them based on the characteristic Degs of each subclusters and combined the function of multiple genes to predict the function of each fibroblast subcluster. At the same time, we added the GSVA method to validate our fibroblast classification criteria in bulk RNA-seq samples, and the results further proved the rationality of our classification.

6. The Data shown in Figure 2 F-H are computational models; they need to be validated by additional experiments. This applies to other claimed findings as well.

Thanks. Immunohistochemical and immunofluorescence staining of normal and disease ligaments were used to validate our results. At the same time, the results of spatial transcriptome sequencing also confirmed the results of scRNA-seq.

7. Considering the effect of biological variables and age on this study, information about the donors' gender and age should be provided.

Thank you for your reminding. We added the information about the donors’ gender and age in the supplementary table.

[Editors’ note: what follows is the authors’ response to the second round of review.]

The manuscript has been improved but there are some remaining issues that need to be addressed, as outlined below:

1. The authors are suggested to run the proofreading carefully since many typos are found in the text. Below is a list of some, but not all.

P3 line 48: consider changing " the musculoskeletal" to "the musculoskeletal system".

P3 line 59: Please consider revising "ACL aging and degradation" to "ACL aging and degeneration". Also, "degradation" has been heavily used in the text. Consider replacing it with "degeneration" wherever it is appropriate.

P4 line 85: replace "synergize" with "synergy"; P10 line 204: please replace "resolve" with "reveal"?

Thanks for your suggestions. We have carefully proofread and revised the whole manuscript.

2. The color codes for Figure 2B are confusing, please consider changing the color codes for L3 to L5 to increase visibility.

Thanks. According to your suggestions, we have changed the color codes in Figure 2B to make the figure clearer.

3. Markers used for pericyte immunostaining in Figure 4G were incorrect! ACTA2, MYH11, and CD90 (Thy1) are not specific markers for pericytes. ACTA2 (also called sSMA) and MYH11 (myosin11) are generally considered markers for myofibroblasts and vascular smooth muscle, not pericytes. Also, most of the mesenchymal progenitors express CD19. So pericyte classification in this ms is confusing and not supported by the immunostaining data.

Thank you for the suggestion. Our definition of pericytes is well-founded and markers we selected used for pericyte immunostaining are also reasonable. Firstly, pericytes and smooth muscle cells originate from the same mesenchymal lineage, share many marker expressions like PDGFRB, ACTA2, CD13, CD1461. Pericytes also present strong muscle contractile gene expression signatures, including ACTA2, MYH/MYL genes, TAGLN2. There are several literature demonstrating pericytes expressing ACTA2 and MYH/MYL genes across multiple human organs such as meniscus2, intervertebral disc3, heart4, brain5. Secondly, it is known that pericytes attach to the surface of vascular endothelial cells and maintain vasculature stability via the endothelium–pericyte crosstalk6. Our immunofluorescence imaging found that this population of cells formed a tube shape and it seems to surround the blood vessels. Therefore, according to the spatial localization, this group of cells also conforms to the characteristics of pericytes. Thirdly, it has been reported that there are two subsets of pericytes. One subset expresses more stem cell markers such as CD90, RGS5, CD146. Another subset expressed more markers associated with muscle contraction, such as MYH11, MYH91, 2. In our study, we also found similar two populations of pericytes and performed multiple immunofluorescent staining according to their markers (ACTA2: shared marker; MYH11: a marker of pericyte1; THY1/CD90: a marker of pericyte2). So, we believe that pericyte classification in our study is clear and can be supported by our immunostaining data.

4. The data in Figure 5F doesn't convincingly demonstrate there is an increase in the expression of FGF7 and TGFβ1 (where and which types of cells) in the specimens from the diseased groups. In addition, these data need to be quantified, and the "n" number should be provided.

Thank you for your suggestions. Fibroblasts are the main cells in the ligament and are the focus of our study. In order to better reflect which types of cells express TGFB1 and FGF7, we conducted immunofluorescence co-staining of Fib.9 marker gene APOE and FGF7 and Fib.8 marker gene TOP2A and TGFB1 according to the results of cell interaction analysis. According to the results, we can see that fib.8 and fib.9 groups in the disease group were significantly more than those in the normal group, and fib8 highly expressed TGFB1 and Fib9 highly expressed FGF7 in the disease group, which was consistent with the results of scRNA-seq.

1. Kumar, A.; D'Souza, S. S.; Moskvin, O. V.; Toh, H.; Wang, B.; Zhang, J.; Swanson, S.; Guo, L. W.; Thomson, J. A.; Slukvin, II, Specification and Diversification of Pericytes and Smooth Muscle Cells from Mesenchymoangioblasts. Cell reports 2017, 19 (9), 1902-1916.

2. Fu, W.; Chen, S.; Yang, R.; Li, C.; Gao, H.; Li, J.; Zhang, X., Cellular features of localized microenvironments in human meniscal degeneration: a single-cell transcriptomic study. eLife 2022, 11.

3. Gan, Y.; He, J.; Zhu, J.; Xu, Z.; Wang, Z.; Yan, J.; Hu, O.; Bai, Z.; Chen, L.; Xie, Y.; Jin, M.; Huang, S.; Liu, B.; Liu, P., Spatially defined single-cell transcriptional profiling characterizes diverse chondrocyte subtypes and nucleus pulposus progenitors in human intervertebral discs. Bone research 2021, 9 (1), 37.

4. Litviňuková, M.; Talavera-López, C.; Maatz, H.; Reichart, D.; Worth, C. L.; Lindberg, E. L.; Kanda, M.; Polanski, K.; Heinig, M.; Lee, M.; Nadelmann, E. R.; Roberts, K.; Tuck, L.; Fasouli, E. S.; DeLaughter, D. M.; McDonough, B.; Wakimoto, H.; Gorham, J. M.; Samari, S.; Mahbubani, K. T.; Saeb-Parsy, K.; Patone, G.; Boyle, J. J.; Zhang, H.; Zhang, H.; Viveiros, A.; Oudit, G. Y.; Bayraktar, O. A.; Seidman, J. G.; Seidman, C. E.; Noseda, M.; Hubner, N.; Teichmann, S. A., Cells of the adult human heart. Nature 2020, 588 (7838), 466-472.

5. Smyth, L. C. D.; Rustenhoven, J.; Scotter, E. L.; Schweder, P.; Faull, R. L. M.; Park, T. I. H.; Dragunow, M., Markers for human brain pericytes and smooth muscle cells. Journal of chemical neuroanatomy 2018, 92, 48-60.

6. Eilken, H. M.; Diéguez-Hurtado, R.; Schmidt, I.; Nakayama, M.; Jeong, H. W.; Arf, H.; Adams, S.; Ferrara, N.; Adams, R. H., Pericytes regulate VEGF-induced endothelial sprouting through VEGFR1. Nature communications 2017, 8 (1), 1574.

https://doi.org/10.7554/eLife.85700.sa2

Article and author information

Author details

  1. Runze Yang

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Resources, Software, Validation, Methodology, Writing - original draft, Project administration
    Competing interests
    No competing interests declared

  2. Tianhao Xu

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Software, Supervision, Validation, Visualization
    Competing interests
    No competing interests declared

  3. Lei Zhang

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Supervision, Validation, Visualization
    Competing interests
    No competing interests declared

  4. Minghao Ge

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Supervision, Validation, Visualization
    Competing interests
    No competing interests declared

  5. Liwei Yan

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Software, Supervision, Methodology
    Competing interests
    No competing interests declared

  6. Jian Li

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Conceptualization, Resources, Data curation, Formal analysis, Supervision
    Competing interests
    No competing interests declared

  7. Weili Fu

    Sports Medicine Center, Department of Orthopedic Surgery/ Orthopedic Research Institute, West China Hospital, Sichuan University, Chengdu, China
    Contribution
    Conceptualization, Data curation, Formal analysis, Supervision, Funding acquisition, Investigation, Writing - original draft, Project administration, Writing - review and editing
    For correspondence
    foxwin2008@163.com
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4438-2760

Funding

National Natural Science Foundation of China (81972123)

  • Weili Fu

National Natural Science Foundation of China (82172508)

  • Weili Fu

National Natural Science Foundation of China (82372490)

  • Weili Fu

1.3.5 Project for Disciplines of Excellence of West China Hospital Sichuan University (ZYJC21030, ZY2017301)

  • Weili Fu

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We appreciate all patients who participated in this study. The authors are grateful to Yi Zhang, Yue Li, and Wanli Zhang for the help in multiplex immunofluorescence staining. We would like to thank Genergy Bio-technology who provided scRNA-seq and spRNA-seq support for this work.

Ethics

This study was reviewed and approved by our University Ethics Committee (Ethics Committee on Biomedical Research, West China Hospital of Sichuan University No. 658 2020-(921)) and all procedures complied with the Helsinki Declaration. Participants gave informed consent to participate in the study.

Senior and Reviewing Editor

  1. Carlos Isales, Augusta University, United States

Version history

  1. Received: December 21, 2022
  2. Preprint posted: January 7, 2023 (view preprint)
  3. Accepted: November 15, 2023
  4. Accepted Manuscript published: November 16, 2023 (version 1)
  5. Version of Record published: November 30, 2023 (version 2)

Copyright

© 2023, Yang et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 205
    Page views
  • 69
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Runze Yang
  2. Tianhao Xu
  3. Lei Zhang
  4. Minghao Ge
  5. Liwei Yan
  6. Jian Li
  7. Weili Fu
(2023)
A single-cell atlas depicting the cellular and molecular features in human anterior cruciate ligamental degeneration: A single cell combined spatial transcriptomics study
eLife 12:e85700.
https://doi.org/10.7554/eLife.85700

Categories and tags

Research organism

Further reading

    1. Cell Biology
    2. Medicine

    Glucagon-like peptide-1 receptor activation stimulates PKA-mediated phosphorylation of Raptor and this contributes to the weight loss effect of liraglutide

    Thao DV Le, Dianxin Liu ... Julio E Ayala
    Research Article Updated

    The canonical target of the glucagon-like peptide-1 receptor (GLP-1R), Protein Kinase A (PKA), has been shown to stimulate mechanistic Target of Rapamycin Complex 1 (mTORC1) by phosphorylating the mTOR-regulating protein Raptor at Ser791 following β-adrenergic stimulation. The objective of these studies is to test whether GLP-1R agonists similarly stimulate mTORC1 via PKA phosphorylation of Raptor at Ser791 and whether this contributes to the weight loss effect of the therapeutic GLP-1R agonist liraglutide. We measured phosphorylation of the mTORC1 signaling target ribosomal protein S6 in Chinese Hamster Ovary cells expressing GLP-1R (CHO-Glp1r) treated with liraglutide in combination with PKA inhibitors. We also assessed liraglutide-mediated phosphorylation of the PKA substrate RRXS*/T* motif in CHO-Glp1r cells expressing Myc-tagged wild-type (WT) Raptor or a PKA-resistant (Ser791Ala) Raptor mutant. Finally, we measured the body weight response to liraglutide in WT mice and mice with a targeted knock-in of PKA-resistant Ser791Ala Raptor. Liraglutide increased phosphorylation of S6 and the PKA motif in WT Raptor in a PKA-dependent manner but failed to stimulate phosphorylation of the PKA motif in Ser791Ala Raptor in CHO-Glp1r cells. Lean Ser791Ala Raptor knock-in mice were resistant to liraglutide-induced weight loss but not setmelanotide-induced (melanocortin-4 receptor-dependent) weight loss. Diet-induced obese Ser791Ala Raptor knock-in mice were not resistant to liraglutide-induced weight loss; however, there was weight-dependent variation such that there was a tendency for obese Ser791Ala Raptor knock-in mice of lower relative body weight to be resistant to liraglutide-induced weight loss compared to weight-matched controls. Together, these findings suggest that PKA-mediated phosphorylation of Raptor at Ser791 contributes to liraglutide-induced weight loss.

    1. Epidemiology and Global Health
    2. Medicine

    Association between bisphosphonate use and COVID-19 related outcomes

    Jeffrey Thompson, Yidi Wang ... Ulrich H von Andrian
    Research Article Updated

    Background:

    Although there are several efficacious vaccines against COVID-19, vaccination rates in many regions around the world remain insufficient to prevent continued high disease burden and emergence of viral variants. Repurposing of existing therapeutics that prevent or mitigate severe COVID-19 could help to address these challenges. The objective of this study was to determine whether prior use of bisphosphonates is associated with reduced incidence and/or severity of COVID-19.

    Methods:

    A retrospective cohort study utilizing payer-complete health insurance claims data from 8,239,790 patients with continuous medical and prescription insurance January 1, 2019 to June 30, 2020 was performed. The primary exposure of interest was use of any bisphosphonate from January 1, 2019 to February 29, 2020. Bisphosphonate users were identified as patients having at least one bisphosphonate claim during this period, who were then 1:1 propensity score-matched to bisphosphonate non-users by age, gender, insurance type, primary-care-provider visit in 2019, and comorbidity burden. Main outcomes of interest included: (a) any testing for SARS-CoV-2 infection; (b) COVID-19 diagnosis; and (c) hospitalization with a COVID-19 diagnosis between March 1, 2020 and June 30, 2020. Multiple sensitivity analyses were also performed to assess core study outcomes amongst more restrictive matches between BP users/non-users, as well as assessing the relationship between BP-use and other respiratory infections (pneumonia, acute bronchitis) both during the same study period as well as before the COVID outbreak.

    Results:

    A total of 7,906,603 patients for whom continuous medical and prescription insurance information was available were selected. A total of 450,366 bisphosphonate users were identified and 1:1 propensity score-matched to bisphosphonate non-users. Bisphosphonate users had lower odds ratios (OR) of testing for SARS-CoV-2 infection (OR = 0.22; 95%CI:0.21–0.23; p<0.001), COVID-19 diagnosis (OR = 0.23; 95%CI:0.22–0.24; p<0.001), and COVID-19-related hospitalization (OR = 0.26; 95%CI:0.24–0.29; p<0.001). Sensitivity analyses yielded results consistent with the primary analysis. Bisphosphonate-use was also associated with decreased odds of acute bronchitis (OR = 0.23; 95%CI:0.22–0.23; p<0.001) or pneumonia (OR = 0.32; 95%CI:0.31–0.34; p<0.001) in 2019, suggesting that bisphosphonates may protect against respiratory infections by a variety of pathogens, including but not limited to SARS-CoV-2.

    Conclusions:

    Prior bisphosphonate-use was associated with dramatically reduced odds of SARS-CoV-2 testing, COVID-19 diagnosis, and COVID-19-related hospitalizations. Prospective clinical trials will be required to establish a causal role for bisphosphonate-use in COVID-19-related outcomes.

    Funding:

    This study was supported by NIH grants, AR068383 and AI155865, a grant from MassCPR (to UHvA) and a CRI Irvington postdoctoral fellowship, CRI2453 (to PH).

    1. Medicine

    Resistance exercise protects mice from protein-induced fat accretion

    Michaela E Trautman, Leah N Braucher ... Dudley W Lamming
    Research Article

    Low-protein (LP) diets extend the lifespan of diverse species and are associated with improved metabolic health in both rodents and humans. Paradoxically, many athletes and bodybuilders consume high-protein (HP) diets and protein supplements, yet are both fit and metabolically healthy. Here, we examine this paradox using weight pulling, a validated progressive resistance exercise training regimen, in mice fed either an LP diet or an isocaloric HP diet. We find that despite having lower food consumption than the LP group, HP-fed mice gain significantly more fat mass than LP-fed mice when not exercising, while weight pulling protected HP-fed mice from this excess fat accretion. The HP diet augmented exercise-induced hypertrophy of the forearm flexor complex, and weight pulling ability increased more rapidly in the exercised HP-fed mice. Surprisingly, exercise did not protect from HP-induced changes in glycemic control. Our results confirm that HP diets can augment muscle hypertrophy and accelerate strength gain induced by resistance exercise without negative effects on fat mass, and also demonstrate that LP diets may be advantageous in the sedentary. Our results highlight the need to consider both dietary composition and activity, not simply calories, when taking a precision nutrition approach to health.