Long non-coding RNA expression profiling of subchondral bone reveals AC005165.1 modifying FRZB expression during osteoarthritis

LncRNA AC005165.1 Alleviates IL-1β-Induced Osteoarthritis via miR-199a-3p/TXNIP Axis

Osteoarthritis (OA) is a chronic musculoskeletal disease and often causes impaired joint mobility and disability. Long noncoding RNAs (lncRNAs) play pivotal roles in OA development. This study was done to explore the role and mechanism of the lncRNA AC005165.1 in the cell model of interleukin-1β (IL)-1β-treated chondrocytes. This study recruited 20 surgically treated OA patients and 12 age- and gender-matched controls. Real-time reverse transcription quantitative polymerase chain reaction was used to examine the expression levels of AC005165.1, miR-199a-3p, and thioredoxin-interacting protein (TXNIP) in articular cartilage of patients and IL-1β-treated human chondrocytes. Cell viability and apoptosis were evaluated by cell counting kit-8 and flow cytometry assays, respectively. The protein levels of inflammatory cytokines were assessed by western blotting. Enzyme-linked immunosorbent assay was conducted to detect the concentrations of the inflammatory cytokines in chondrocytes. Luciferase reporter assay and Pearson's correlation analysis were used for analyzing the interaction and the correlation among AC005165.1, miR-199a-3p, and TXNIP. AC005165.1 expression was downregulated in cartilage of OA patients and chondrocytes treated with IL-1β, compared to that in the control groups. AC005165.1 knockdown increased apoptosis and aggravated inflammatory response in IL-1β-treated chondrocytes. AC005165.1 interacted with miR-199a-3p, and TXNIP was targeted by miR-199a-3p. In rescue assay, miR-199a-3p knockdown and TXNIP overexpression significantly reduced apoptosis and mitigated inflammatory response in IL-1β-treated chondrocytes with AC005165.1 knockdown. AC005165.1 knockdown promoted apoptosis and inflammatory response in IL-1β-treated chondrocytes via the miR-199a-3p/TXNIP axis.

Biological functions and applications of LncRNAs in the regulation of the extracellular matrix in osteoarthritis

Osteoarthritis (OA) is a major cause of disability, characterized by chronic pain, irreversible destruction, and loss of function of the articular cartilage. The integrity and arrangement of the composition and structure of the extracellular matrix (ECM) are essential for maintaining the elasticity, integrity, and mechanical support function of the cartilage tissue. Osteoarthritis causes substantial changes in the ECM, driving the progression of the disease. Recent studies have shown that the ECM plays a critical role in the development of cartilage tissue as well as the occurrence and development of osteoarthritis by directly or indirectly regulating chondrocyte proliferation, apoptosis, differentiation, and gene expression. Long non-coding RNAs (lncRNAs) are a class of non-coding RNAs derived from large transcripts. Mutations and disorders of lncRNAs are closely related to the development of osteoarthritis. Abnormal expression of lncRNAs in osteoarthritic cartilage regulates the synthesis and decomposition of the cartilaginous ECM. Therefore, the use of lncRNAs as nucleic acid drugs that regulate their targets may reduce ECM degradation, thereby delaying the pathological progression of osteoarthritis. In this review, the regulatory effects of lncRNAs on ECM in different cell behaviors related to OA are summarized. The roles of lncRNAs in the proliferation, apoptosis, differentiation, and ECM-related gene activity of chondrocytes, as well as the application of lncRNAs as potential gene therapy drugs for the repair and regeneration of osteoarthritic tissue, are also reviewed. A better understanding of the roles of lncRNAs in guiding chondrocyte behavior and ECM metabolism is critical for their future applications in osteoarthritis therapy and regenerative medicine.

Coding and noncoding RNA profile of human heterotopic ossifications - Risk factors and biomarkers

Heterotopic ossification (HO) means the formation of bone in muscles and soft tissues, such as ligaments or tendons. HO could have a genetic history or develop after a traumatic event, as a result of muscle injury, fractures, burns, surgery, or neurological disorders. Many lines of evidence suggest that the formation of HO is related to the pathological differentiation of stem or progenitor cells present within soft tissues or mobilized from the bone marrow. The cells responsible for the initiation and progression of HO are generally called HO precursor cells. The exact mechanisms behind the development of HO are not fully understood. However, several factors have been identified as potential contributors. For example, local tissue injury and inflammation disturb soft tissue homeostasis. Inflammatory cells release growth factors and cytokines that promote osteogenic or chondrogenic differentiation of HO precursor cells. The bone morphogenetic protein (BMP) is one of the main factors involved in the development of HO. In this study, next-generation sequencing (NGS) and RT-qPCR were performed to analyze the differences in mRNA, miRNA, and lncRNA expression profiles between muscles, control bone samples, and HO samples coming from patients who underwent total hip replacement (THR). As a result, crucial changes in the level of gene expression between HO and healthy tissues were identified. The bioinformatic analysis allowed to describe the processes most severely impacted, as well as genes which level differed the most significantly between HO and control samples. Our analysis showed that the level of transcripts involved in leukocyte migration, differentiation, and activation, as well as markers of chronic inflammatory diseases, that is, miR-148, increased in HO, as compared to muscle. Furthermore, the levels of miR-195 and miR-143, which are involved in angiogenesis, were up-regulated in HO, as compared to bone. Thus, we suggested that inflammation and angiogenesis play an important role in HO formation. Importantly, we noticed that HO is characterized by a higher level of TLR3 expression, compared to muscle and bone. Thus, we suggest that infection may also be a risk factor in HO development. Furthermore, an increased level of transcripts coding proteins involved in osteogenesis and signaling pathways, such as ALPL, SP7, BGLAP, BMP8A, BMP8B, SMPD3 was noticed in HO, as compared to muscles. Interestingly, miR-99b, miR-146, miR-204, and LINC00320 were up-regulated in HO, comparing to muscles and bone. Therefore, we suggested that these molecules could be important biomarkers of HO formation and a potential target for therapies.

Comprehensive analysis on subchondral bone marrow lesions of human osteoarthritis by integrating bulk and single-cell transcriptomes

OBJECTIVE

This study aims to demonstrate the cellular composition and underlying mechanisms in subchondral bone marrow lesions (BMLs) of knee osteoarthritis (OA).

METHODS

BMLs were assessed by MRI Osteoarthritis Knee Score (MOAKS)≥2. Bulk RNA-sequencing (bulk-seq) and BML-specific differentially expressed genes (DEGs) analysis were performed among subchondral bone samples (including OA-BML=3, paired OA-NBML=3; non-OA=3). The hub genes of BMLs were identified by verifying in independent datasets and multiple bioinformatic analyses. To further estimate cell-type composition of subchondral bone, we utilized two newly developed deconvolution algorithms (MuSiC, MCP-counter) in transcriptomic datasets, based on signatures from open-accessed single-cell RNA sequencing (scRNA-seq). Finally, competing endogenous RNA (ceRNA) and transcription factor (TF) networks were constructed through multiple predictive databases, and validated by public non-coding RNA profiles.

RESULTS

A total of 86 BML-specific DEGs (up 79, down 7) were identified. IL11 and VCAN were identified as core hub genes. The "has-miR-424-5p/lncRNA PVT1" was determined as crucial network, targeting IL11 and VCAN, respectively. More importantly, two deconvolution algorithms produced approximate estimations of cell-type composition, and the cluster of heterotopic-chondrocyte was discovered abundant in BMLs, and positively correlated with the expression of hub genes.

CONCLUSION

IL11 and VCAN were identified as the core hub genes of BMLs, and their molecular networks were determined as well. We profiled the characteristics of subchondral bone at single-cell level and determined that the heterotopic-chondrocyte was abundant in BMLs and was closely linked to IL11 and VCAN. Our study may provide new insights into the microenvironment and pathological molecular mechanism of BMLs, and could lead to novel therapeutic strategies.

Advances in Research on the Regulatory Roles of lncRNAs in Osteoarthritic Cartilage

Osteoarthritis (OA) is the most common degenerative bone and joint disease that can lead to disability and severely affect the quality of life of patients. However, its etiology and pathogenesis remain unclear. It is currently believed that articular cartilage lesions are an important marker of the onset and development of osteoarthritis. Long noncoding RNAs (lncRNAs) are a class of multifunctional regulatory RNAs that are involved in various physiological functions. There are many differentially expressed lncRNAs between osteoarthritic and normal cartilage tissues that play multiple roles in the pathogenesis of OA. Here, we reviewed lncRNAs that have been reported to play regulatory roles in the pathological changes associated with osteoarthritic cartilage and their potential as biomarkers and a therapeutic target in OA to further elucidate the pathogenesis of OA and provide insights for the diagnosis and treatment of OA.

RNA-seq identifies differentially expressed genes involved in csal1 overexpression in granulosa cells of prehierarchical follicles in Chinese Dagu hens

The transcription factor csal1 is an important molecule that plays a critical regulatory function in ovarian follicle development, as confirmed by our previous data. However, the candidate genes of csal1 and its regulatory mechanism remain poorly understood in the granulosa cells (GCs) of chicken prehierarchical follicles (PFs). Six transcriptomes of csal1 and empty vector were analyzed in Chinese Dagu hens by RNA sequencing. Six cDNA libraries were constructed, with more than 42 million clean reads and 16,779 unigenes. Of these 16,779 unigenes, 2,762 differentially expressed genes (DEGs) were found in GCs, including 1,605 upregulated and 1,157 downregulated unigenes. Fourteen genes, including BMP5, TACR2, AMH, PLAG1, MYOD1, BOP1, SIPA1, NOTCH1, BCL2L1, SOX9, ADGRA2, WNT5A, SLC7A11, and GATAD2B, were related to GC proliferation and differentiation, hormone production, ovarian follicular development, regulation of reproductive processes, and signaling pathways in the PFs. Further analysis demonstrated the DEGs in GCs of ovarian follicles were enriched in neuroactive ligand-receptor interaction, cell adhesion molecules, and pathways related to cytochrome P450, indicating a critical function for csal1 in the generation of egg-laying features by controlling ovarian follicle development. For the first time, the current study represents the transcriptome analysis with ectopic csal1 expression. These findings provide significant evidence for investigating the molecular mechanism by which csal1 controls PF development in the hen ovary.

Emerging role of lncRNAs in osteoarthritis: An updated review

Osteoarthritis (OA) is a prevalent joint disease, which is associated with progressive articular cartilage loss, synovial inflammation, subchondral sclerosis and meniscus injury. The molecular mechanism underlying OA pathogenesis is multifactorial. Long non-coding RNAs (lncRNAs) are non-protein coding RNAs with length more than 200 nucleotides. They have various functions such as modulating transcription and protein activity, as well as forming endogenous small interfering RNAs (siRNAs) and microRNA (miRNA) sponges. Emerging evidence suggests that lncRNAs might be involved in the pathogenesis of OA which opens up a new avenue for the development of new biomarkers and therapeutic strategies. The purpose of this review is to summarize the current clinical and basic experiments related to lncRNAs and OA with a focus on the extensively studied H19, GAS5, MALAT1, XIST and HOTAIR. The potential translational value of these lncRNAs as therapeutic targets for OA is also discussed.

Allelic expression imbalance in articular cartilage and subchondral bone refined genome-wide association signals in osteoarthritis

OBJECTIVES

To present an unbiased approach to identify positional transcript single nucleotide polymorphisms (SNPs) of osteoarthritis (OA) risk loci by allelic expression imbalance (AEI) analyses using RNA sequencing of articular cartilage and subchondral bone from OA patients.

METHODS

RNA sequencing from 65 articular cartilage and 24 subchondral bone from OA patients was used for AEI analysis. AEI was determined for all genes present in the 100 regions reported by the GWAS catalog that were also expressed in cartilage or bone. The count fraction of the alternative allele (φ) was calculated for each heterozygous individual with the risk-SNP or with the SNP in linkage disequilibrium (LD) with it (r2 > 0.6). Furthermore, a meta-analysis was performed to generate a meta-φ (null hypothesis median φ = 0.49) and P-value for each SNP.

RESULTS

We identified 30 transcript SNPs (28 in cartilage and 2 in subchondral bone) subject to AEI in 29 genes. Notably, 10 transcript SNPs were located in genes not previously reported in the GWAS catalog, including two long intergenic non-coding RNAs (lincRNAs), MALAT1 (meta-φ = 0.54, FDR = 1.7x10-4) and ILF3-DT (meta-φ = 0.6, FDR = 1.75x10-5). Moreover, 12 drugs were interacting with 7 genes displaying AEI, of which 7 drugs have been already approved.

CONCLUSIONS

By prioritizing proxy transcript SNPs that mark AEI in cartilage and/or subchondral bone at loci harboring GWAS signals, we present an unbiased approach to identify the most likely functional OA risk-SNP and gene. We identified 10 new potential OA risk genes ready for further, translation towards underlying biological mechanisms.

Long Non-Coding RNA PVT1 and Its Target miRNA-146a as Potential Prognostic Biomarkers in Rheumatoid Arthritis Patients

Objective: Long non-coding RNAs (lncRNAs) and their target microRNAs were documented in multiple studies to have a significant role in different joint disorders such as rheumatoid arthritis (RA) and osteoarthritis (OA). The current work aimed to determine the potential role of lnc-PVT1 and miR-146a as promising biomarkers to distinguish between RA, OA patients, and healthy individuals. Methods: The expression levels of lnc-PVT1 and its target miR-146a in the serum were measured for three different groups, including patients with RA (40), OA patients (40), and healthy controls (HCs) (40). Participating individuals were subjected to a full history investigation and clinical examination. Blood samples were tested for ESR, RF, CBC, as well as liver and renal functions. Serum was used to detect the relative expression levels of lnc-PVT1 and miR-146a and we correlated the levels with RA and OA activity and severity signs. Results: Lnc-PVT1 expression level was greater among patients with RA compared to that of OA patients, with a fold change median of 2.62 and 0.22, respectively (p = 0.001). The miR-146a fold change was significantly demonstrated between the RA, OA, and HCs groups. There was no correlation between both biomarkers with the disease activity scales (DAS28) of RA, the Knee injury Osteoarthritis Outcome Score (KOOS), or any sign of detection of the disease severity of OA. Conclusions: lnc-PVT1 and miR-146a could be considered as promising biomarkers for the diagnosis of RA and OA and may have an important role as therapeutic targets in the future.