Integrative Single-Cell RNA-Seq and ATAC-Seq Analysis of Peripheral Mononuclear Cells in Patients With Ankylosing Spondylitis

Advancing skeletal health and disease research with single-cell RNA sequencing

Orthopedic conditions have emerged as global health concerns, impacting approximately 1.7 billion individuals worldwide. However, the limited understanding of the underlying pathological processes at the cellular and molecular level has hindered the development of comprehensive treatment options for these disorders. The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized biomedical research by enabling detailed examination of cellular and molecular diversity. Nevertheless, investigating mechanisms at the single-cell level in highly mineralized skeletal tissue poses technical challenges. In this comprehensive review, we present a streamlined approach to obtaining high-quality single cells from skeletal tissue and provide an overview of existing scRNA-seq technologies employed in skeletal studies along with practical bioinformatic analysis pipelines. By utilizing these methodologies, crucial insights into the developmental dynamics, maintenance of homeostasis, and pathological processes involved in spine, joint, bone, muscle, and tendon disorders have been uncovered. Specifically focusing on the joint diseases of degenerative disc disease, osteoarthritis, and rheumatoid arthritis using scRNA-seq has provided novel insights and a more nuanced comprehension. These findings have paved the way for discovering novel therapeutic targets that offer potential benefits to patients suffering from diverse skeletal disorders.

Single-cell transcriptome revealed dysregulated RNA-binding protein expression patterns and functions in human ankylosing spondylitis

Objective

To explore the expression characteristics and regulatory patterns of RBPs in different immune cell types of AS, and to clarify the potential key role of RBPs in the occurrence and development of AS disease.

Methods

PBMC sample data from scRNA-seq (HC*29, AS*10) and bulk RNA-seq (NC*3, AS*5) were selected for correlation analysis.

Results

(1) Compared with the HC group, the numbers of B, DC (dendritic cells), CD14+ Mono and CD8+ T cells were increased in AS group, while the numbers of platelet (platelets), CD8+ NKT, CD16+ Mono (non-classical monocytes), Native CD4+ T and NK were decreased. (2) Through the analysis of RBP genes in B cells, some RBPs were found to play an important role in B cell differentiation and function, such as DDX3X, SFPQ, SRRM1, UPF2. (3) It may be related to B-cell receptor, IgA immunity, NOD-like receptor and other signaling pathways; Through the analysis of RBP genes in CD8+ T cells, some RBPs that play an important role in the immune regulation of CD8+ T were found, such as EIF2S3, EIF4B, HSPA5, MSL3, PABPC1 and SRSF7; It may be related to T cell receptor, TNF, IL17 and other signaling pathways. (4) Based on bulk RNA-seq, it was found that compared with HC and AS patients, differentially expressed variable splicing genes (RASGs) may play an important role in the occurrence and development of AS by participating in transcriptional regulation, protein phosphorylation and ubiquitination, DNA replication, angiogenesis, intracellular signal transduction and other related pathways.

Conclusion

RBPs has specific expression characteristics in different immune cell types of AS patients, and has important regulatory functions. Its abnormal expression and regulation may be closely related to the occurrence and development of AS.

Functional CRISPR screens in T cells reveal new opportunities for cancer immunotherapies

T cells are fundamental components in tumour immunity and cancer immunotherapies, which have made immense strides and revolutionized cancer treatment paradigm. However, recent studies delineate the predicament of T cell dysregulation in tumour microenvironment and the compromised efficacy of cancer immunotherapies. CRISPR screens enable unbiased interrogation of gene function in T cells and have revealed functional determinators, genetic regulatory networks, and intercellular interactions in T cell life cycle, thereby providing opportunities to revamp cancer immunotherapies. In this review, we briefly described the central roles of T cells in successful cancer immunotherapies, comprehensively summarised the studies of CRISPR screens in T cells, elaborated resultant master genes that control T cell activation, proliferation, fate determination, effector function, and exhaustion, and highlighted genes (BATF, PRDM1, and TOX) and signalling cascades (JAK-STAT and NF-κB pathways) that extensively engage in multiple branches of T cell responses. In conclusion, this review bridged the gap between discovering element genes to a specific process of T cell activities and apprehending these genes in the global T cell life cycle, deepened the understanding of T cell biology in tumour immunity, and outlined CRISPR screens resources that might facilitate the development and implementation of cancer immunotherapies in the clinic.

CD_99 G1 neutrophils modulate osteogenic differentiation of mesenchymal stem cells in the pathological process of ankylosing spondylitis

OBJECTIVES

This study aimed to identify the types and heterogeneity of cells within the spinal enthesis and investigate the underlying mechanisms of osteogenesis.

METHODS

Single-cell RNA sequencing was used to identify cell populations and their gene signatures in the spinal enthesis of five patients with ankylosing spondylitis (AS) and three healthy individuals. The transcriptomes of 40 065 single cells were profiled and divided into 7 clusters: neutrophils, monocytic cells, granulomonocytic progenitor_erythroblasts, T cells, B cells, plasma cells and stromal cells. Real-time quantitative PCR, immunofluorescence, flow cytometry, osteogenesis induction, alizarin red staining, immunohistochemistry, short hairpin RNA and H&E staining were applied to validate the bioinformatics analysis.

RESULTS

Pseudo-time analysis showed two differentiation directions of stromal cells from the mesenchymal stem cell subpopulation MSC-C2 to two Cxcl12-abundant-reticular (CAR) cell subsets, Osteo-CAR and Adipo-CAR, within which three transcription factors, C-JUN, C-FOS and CAVIN1, were highly expressed in AS and regulated the osteogenesis of mesenchymal stem cells. A novel subcluster of early-stage neutrophils, CD99_G1, was elevated in AS. The proinflammatory characteristics of monocyte dendritic cell progenitor-recombinant adiponectin receptor 2 monocytic cells were explored. Interactions between Adipo-CAR cells, CD99_G1 neutrophils and other cell types were mapped by identifying ligand-receptor pairs, revealing the recruitment characteristics of CD99_G1 neutrophils by Adipo-CAR cells and the pathogenesis of osteogenesis induced in AS.

CONCLUSIONS

Our results revealed the dynamics of cell subpopulations, gene expression and intercellular interactions during AS pathogenesis. These findings provide new insights into the cellular and molecular mechanisms of osteogenesis and will benefit the development of novel therapeutic strategies.

Identification of anterior cruciate ligament fibroblasts and their contribution to knee osteoarthritis progression using single-cell analyses

The mechanical properties of the anterior cruciate ligament (ACL) in the knee have been highlighted, but its role in the regulation of the joint microenvironment remains unclear, especially in the progression of Knee Osteoarthritis (KOA). Here, single-cell RNA sequencing (scRNA-seq) and single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq) data were integrated to reveal the transcriptional and epigenomic landscape of ACL in normal and OA states. We identified a novel subpopulation of fibroblasts in ACL, which provides new insights into the role of the ACL in knee homeostasis and disease. Degeneration of the ACL during OA mechanically alters the knee joint homeostasis and influences the microenvironment by regulating inflammatory- and osteogenic-related factors, thereby contributing to the progression of KOA. Additionally, the specific mechanism by which these Inflammation-associated Fibroblasts (IAFs) regulate KOA progression was uncovered, providing new foundation for the development of targeted treatments for KOA.

Network pharmacology combined with molecular docking and experimental validation to explore the potential mechanism of Cinnamomi ramulus against ankylosing spondylitis

BACKGROUND

Cinnamomi ramulus (C. ramulus) is frequently employed in the treatment of ankylosing spondylitis (AS). However, the primary constituents, drug targets, and mechanisms of action remain unidentified.

METHODS

In this study, various public databases and online tools were employed to gather information on the compounds of C. ramulus, drug targets, and disease targets associated with ankylosing spondylitis. The intersection of drug targets and disease targets was then determined to identify the common targets, which were subsequently used to construct a protein-protein interaction (PPI) network using the STRING database. Network analysis and the analysis of hub genes and major compounds were conducted using Cytoscape software. Furthermore, the Metascape platform was utilized for Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses. Molecular docking studies and immunohistochemical experiments were performed to validate the core targets.

RESULTS

The network analysis identified 2-Methoxycinnamaldehyde, cinnamaldehyde, and 2-Hydroxycinnamaldehyde as the major effective compounds present in C. ramulus. The PPI network analysis revealed PTGS2, MMP9, and TLR4 as the most highly correlated targets. GO and KEGG analyses indicated that C. ramulus exerts its therapeutic effects in ankylosing spondylitis through various biological processes, including the response to hormones and peptides, oxidative stress response, and inflammatory response. The main signaling pathways involved were IL-17, TNF, NF-kappa B, and Toll-like receptor pathways. Molecular docking analysis confirmed the strong affinity between the key compounds and the core targets. Additionally, immunohistochemical analysis demonstrated an up-regulation of PTGS2, MMP9, and TLR4 levels in ankylosing spondylitis.

CONCLUSIONS

This study provides insights into the effective compounds, core targets, and potential mechanisms of action of C. ramulus in the treatment of ankylosing spondylitis. These findings establish a solid groundwork for future fundamental research in this field.

Single-cell analysis of the miRNA activities in tuberculous meningitis (TBM) model mice injected with the BCG vaccine

BACKGROUND

Our previous study revealed the transcriptome atlas of specific cell types in tuberculous meningitis (TBM) model mice injected with the BCG vaccine via scRNA sequencing. However, the activities of miRNAs in TBM at single-cell resolution remain to be explored.

METHOD

Cell type-specific miRNA activities were investigated by using motif enrichment analyses (miReact) on the transcriptome data of 15 cell types. The target mRNAs of miRNAs were predicted and subjected to enrichment analysis. Furthermore, miRNAs and their target mRNAs with opposite expression trends were chosen to construct functional networks. Besides, qRT-PCR and RNA scope were performed to verify the expression level of representative miRNA.

RESULTS

The tSNE dimensionality reduction presented 15 cell types in TBM model mice, in which microglia and endothelial cells accounted for the majority. Target mRNAs of each cell type were predicted for verification or network construction. The immune and inflammation-related miRNA-mRNA networks of macrophages and microglia, oxidative phosphorylation-related miRNA-mRNA networks of neurons, ion and protein transport-related networks of epididymal cells, and angiogenesis-related miRNA-mRNA networks of VSMCs were constructed. The miRNA activity analysis revealed that miR-21a-3p activity was increased in microglia, macrophages, neurons and epididymal cells. The result of qRT-PCR and RNA scope indicate that miR-21a-3p was significantly higher-expressed in TBM brain tissue compared with normal brain tissue.

CONCLUSION

In our study, an in-depth exploration of the mRNA expression and miRNA activity of macrophages, microglia, epididymal cells, neurons and vascular smooth muscle cells during TBM progression was conducted using scRNA-Seq, which provided novel insights into the immune cell engagement in TBM patients.

Emerging Concepts in Precision Medicine in Axial Spondyloarthritis

PURPOSE OF REVIEW

Axial spondyloarthritis (AxSpA) is among the rheumatology's most heritable complex diseases, yet precision medicine at clinics still needs to be explored. We reviewed the emerging concepts and recent developments in polygenic risk scores, Mendelian randomization, pharmacogenomics, single-cell sequencing, and spatial transcriptomics.

RECENT FINDINGS

Polygenic risk score has resulted in encouraging results with potential diagnostic utility as it appears to outperform current diagnostic tools. Its performance and generalizability vary with ethnicity. Mendelian randomization has elucidated multiple causal associations, particularly between inflammatory bowel disease and AxSpA. Single-cell transcriptomics (particularly scRNA-seq and scATAC-seq) has identified numerous cell types, including synovial and blood immunological cells, to understand the contribution of both innate and adaptative immunity in AxSpA. Current molecular tools provide an exciting opportunity to advance precision medicine for AxSpA patients. However, extensive research and implementation strategies are still required before they can have an impact in the clinic.

Application of single-cell RNA sequencing methods to develop B cell targeted treatments for autoimmunity

The COVID-19 pandemic coincided with several transformative advances in single-cell analysis. These new methods along with decades of research and trials with antibody therapeutics and RNA based technologies allowed for highly effective vaccines and treatments to be produced at astonishing speeds. While these tools were initially focused on models of infection, they also show promise in an autoimmune setting. Self-reactive B cells play important roles as antigen-presenting cells and cytokine and autoantibody producers for many autoimmune diseases. Yet, current therapies to target autoreactive B cells deplete all B cells irrespective of their pathogenicity. Development of self-reactive B cell targeting therapies that would spare non-pathogenic B cells are needed to treat disease while allowing effective immune responses to other ailments. Single-cell RNA sequencing (scRNA-seq) approaches will aid in identification of the pathogenic self-reactive B cells operative in autoimmunity and help with development of more favorable precision targeted therapies.

Disease-specific expansion of CD29+IL-17RA+ T effector cells possessing multiple signalling pathways in spondyloarthritis

OBJECTIVES

T cells adhere to enthesis fibrocartilage via integrins and intrinsically require IL-17RA-mediated signals to maintain their effector function. We analysed CD29+IL-17RA+ T cells in inflamed lesions and peripheral blood in patients with SpA and investigated their association with disease activity and therapeutic response.

METHODS

Transcriptome analysis of synovial fluid T cells from PsA was performed using publicly available bulk cell RNA sequencing data. Blood samples were obtained from healthy controls (n = 37), RA (n = 12), IgG4-related disease (IgG4-RD; n = 12), large vessel vasculitis (LVV; n = 12) and SpA (n = 28) and were analysed by flow cytometry.

RESULTS

T cells in the inflamed joints of PsA showed CD29 and IL-17RA expression. CD29+IL-17RA+ T cells showed enriched CXCR3+CD45RA+ effector cells and activation of spleen tyrosine kinase (Syk), nuclear factor κB (NF-κB) and Janus kinase-signal transducer and activator of transcription (JAK-STAT) pathways. The proportion of peripheral blood CD29+IL-17RA+ T cells was significantly increased in patients with SpA compared with patients with RA, IgG4-RD or LVV and in healthy controls. Based on the ASDAS-CRP scores, the proportion of CD29+IL-17RA+ T cells was positively correlated with disease activity in treatment-naïve patients with active SpA. Anti-IL-17 but not anti-TNF monoclonal antibodies reduced CD29+IL-17RA+ T cells.

CONCLUSIONS

CD29+IL-17RA+ T effector cells with enhanced Syk, NF-κB and JAK-STAT pathways were specifically increased in SpA and were correlated with disease activity, implicating a role of this newly identified T cell population in the pathogenesis. Anti-IL-17 monoclonal antibodies may be effective for patients by reducing this pathogenic T cell population.

Single-Cell Genomics for Investigating Pathogenesis of Inflammatory Diseases

Recent technical advances have enabled unbiased transcriptomic and epigenetic analysis of each cell, known as "single-cell analysis". Single-cell analysis has a variety of technical approaches to investigate the state of each cell, including mRNA levels (transcriptome), the immune repertoire (immune repertoire analysis), cell surface proteins (surface proteome analysis), chromatin accessibility (epigenome), and accordance with genome variants (eQTLs; expression quantitative trait loci). As an effective tool for investigating robust immune responses in coronavirus disease 2019 (COVID-19), many researchers performed single-cell analysis to capture the diverse, unbiased immune cell activation and differentiation. Despite challenges elucidating the complicated immune microenvironments of chronic inflammatory diseases using existing experimental methods, it is now possible to capture the simultaneous immune features of different cell types across inflamed tissues using various single-cell tools. In this review, we introduce patient-based and experimental mouse model research utilizing single-cell analyses in the field of chronic inflammatory diseases, as well as multi-organ atlas targeting immune cells.

No Evidence to Support a Causal Relationship between Circulating Adiponectin Levels and Ankylosing Spondylitis: A Bidirectional Two-Sample Mendelian Randomization Study

Based on previous observational studies, the causal association between circulating adiponectin (CA) levels and ankylosing spondylitis (AS) risk remains unclear. Therefore, this study aims to investigate whether CA levels are related to the risk of AS. We carried out a bidirectional two-sample Mendelian randomization (MR) analysis to examine the causal correlation between CA levels and AS via published genome-wide association study (GWAS) datasets. Single-nucleotide polymorphisms (SNPs) related to CA levels were derived from a large GWAS that included 39,883 individuals of European descent. SNPs related to AS were obtained from the FinnGen consortium (2252 cases and 227,338 controls). The random-effects inverse variance weighted (IVW) method was the primary method utilized in our research. We also used four complementary approaches to improve the dependability of this study (MR-Egger regression, Weighted median, Weighted mode, and Simple mode). Random-effects IVW (odds ratio [OR], 1.00; 95% confidence interval [CI], 0.79-1.27, p = 0.984) and four complementary methods all indicated that genetically predicted CA levels were not causally related to the risk of AS. In reverse MR analysis, there is little evidence to support the genetic causality between the risk of AS and CA levels.

Single Cell Transcriptome and Surface Epitope Analysis of Ankylosing Spondylitis Facilitates Disease Classification by Machine Learning

Ankylosing spondylitis (AS) is an immune-mediated inflammatory disorder that primarily affects the axial skeleton, especially the sacroiliac joints and spine. This results in chronic back pain and, in extreme cases, ankylosis of the spine. Despite its debilitating effects, the pathogenesis of AS remains to be further elucidated. This study used single cell CITE-seq technology to analyze peripheral blood mononuclear cells (PBMCs) in AS and in healthy controls. We identified a number of molecular features associated with AS. CD52 was found to be overexpressed in both RNA and surface protein expression across several cell types in patients with AS. CD16⁺ monocytes overexpressed TNFSF10 and IL-18Rα in AS, while CD8⁺ T_EM cells and natural killer cells overexpressed genes linked with cytotoxicity, including GZMH, GZMB, and NKG7. Tregs underexpressed CD39 in AS, suggesting reduced functionality. We identified an overrepresented NK cell subset in AS that overexpressed CD16, CD161, and CD38, as well as cytotoxic genes and pathways. Finally, we developed machine learning models derived from CITE-seq data for the classification of AS and achieved an Area Under the Receiver Operating Characteristic (AUROC) curve of > 0.95. In summary, CITE-seq identification of AS-associated genes and surface proteins in specific cell subsets informs our understanding of pathogenesis and potential new therapeutic targets, while providing new approaches for diagnosis via machine learning.