Identification of changed expression of mRNAs and lncRNAs in osteoarthritic synovium by RNA-sequencing

LncRNA and mRNA expression characteristic and bioinformatic analysis in myocardium of diabetic cardiomyopathy mice

BACKGROUND

Diabetic cardiomyopathy (DCM) is becoming a very well-known clinical entity and leads to increased heart failure in diabetic patients. Long non-coding RNAs (LncRNAs) play an important role in the pathogenesis of DCM. In the present study, the expression profiles of lncRNAs and mRNAs were illuminated in myocardium from DCM mice, with purpose of exploring probable pathological processes of DCM involved by differentially expressed genes in order to provide a new direction for the future researches of DCM.

RESULTS

The results showed that a total of 93 differentially expressed lncRNA transcripts and 881 mRNA transcripts were aberrantly expressed in db/db mice compared with the controls. The top 6 differentially expressed lncRNAs like up-regulated Hmga1b, Gm8909, Gm50252 and down-regulated Msantd4, 4933413J09Rik, Gm41414 have not yet been reported in DCM. The lncRNAs-mRNAs co-expression network analysis showed that LncRNA 2610507I01Rik, 2310015A16Rik, Gm10503, A930015D03Rik and Gm48483 were the most relevant to differentially expressed mRNAs.

CONCLUSION

Our results showed that db/db DCM mice exist differentially expressed lncRNAs and mRNAs in hearts. These differentially expressed lncRNAs may be involved in the pathological process of cardiomyocyte apoptosis and fibrosis in DCM.

The emerging role of lncRNAs in osteoarthritis development and potential therapy

Osteoarthritis impairs the functions of various joints, such as knees, hips, hands and spine, which causes pain, swelling, stiffness and reduced mobility in joints. Multiple factors, including age, joint injuries, obesity, and mechanical stress, could contribute to osteoarthritis development and progression. Evidence has demonstrated that genetics and epigenetics play a critical role in osteoarthritis initiation and progression. Noncoding RNAs (ncRNAs) have been revealed to participate in osteoarthritis development. In this review, we describe the pivotal functions and molecular mechanisms of numerous lncRNAs in osteoarthritis progression. We mention that long noncoding RNAs (lncRNAs) could be biomarkers for osteoarthritis diagnosis, prognosis and therapeutic targets. Moreover, we highlight the several compounds that alleviate osteoarthritis progression in part via targeting lncRNAs. Furthermore, we provide the future perspectives regarding the potential application of lncRNAs in diagnosis, treatment and prognosis of osteoarthritis.

Construction and analysis of a lncRNA–miRNA–mRNA competing endogenous RNA network from inflamed and normal synovial tissues after anterior cruciate ligament and/or meniscus injuries

Background: Despite ample evidence demonstrating that anterior cruciate ligament (ACL) and meniscus tears are associated with posttraumatic osteoarthritis (PTOA) development, the contributing factors remain unknown. Synovial inflammation has recently been recognized as a pivotal factor in the pathogenesis of OA. However, there is a lack of data on synovial profiles after ACL or meniscus injuries, which may contribute to PTOA.

Methods: Twelve patients with ACL tears and/or meniscus injuries were recruited. During surgery, synovial tissues were obtained from the injured knees. The inflammation status of the synovium was characterized according to macroscopic criteria and histological synovitis grades. Then the synovial tissues were classified as control group or inflamed group. High-throughput RNA sequencing of the synovial samples (3 vs. 3) was conducted to identify differentially expressed (DE) RNAs. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, and protein–protein interaction (PPI) analyses were performed to investigate DE mRNAs. Next, competing endogenous RNA (ceRNA) networks were constructed based on bioinformatics analyses. Associations of the identified DE genes (DEGs) with infiltrating immune cells were explored using Pearson correlation analysis.

Results: The results showed that 2793 mRNAs, 3392 lncRNAs and 211 miRNAs were significantly DE between two groups. The top 3 significantly upregulated GO terms and KEGG pathways were immune response, adaptive immune response and immune system process, systemic lupus erythematosus, haematopoietic cell lineage and cytokine–cytokine receptor interaction, respectively. In PPI networks, the top 10 hub genes were IL6, CCR7, C3, CCR5, CXCR3, CXCL8, IL2, CCR3, CCR2 and CXCL1. Seven mRNAs (EPHA5, GSN, ORC1, TLN2, SOX6, NKD2 and ADAMTS19), 4 lncRNAs (MIR4435-2HG, TNXA, CEROX1 and TMEM92-AS1) and 3 miRNAs (miR-486-5p, miR-199a-3p and miR-21-3p) were validated by quantitative real-time polymerase chain reaction and sub-networks were constructed. In correlation analysis, MMP9 correlated positively with M0 macrophages and plasma cells, NKD2 positively with CD8 T cells, and CCR7 and IL2RB positively with naive B cells.

Conclusion: Our study provides foundational synovial inflammation profiles following knee trauma. The ceRNA and PPI networks provide new insight into the biological processes and underlying mechanisms of PTOA. The differential infiltration profiles of immune cells in synovium may contribute to PTOA development. This study also highlights immune-related DEGs as potential PTOA treatment biomarkers.

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.

Long Non-coding RNAs in Rheumatology

The last decade has seen an enormous increase in long non-coding RNA (lncRNA) research within rheumatology. LncRNAs are arbitrarily classed as non-protein encoding RNA transcripts that exceed 200 nucleotides in length. These transcripts have tissue and cell specific patterns of expression and are implicated in a variety of biological processes. Unsurprisingly, numerous lncRNAs are dysregulated in rheumatoid conditions, correlating with disease activity and cited as potential biomarkers and targets for therapeutic intervention. In this chapter, following an introduction into each condition, we discuss the lncRNAs involved in rheumatoid arthritis, osteoarthritis and systemic lupus erythematosus. These inflammatory joint conditions share several inflammatory signalling pathways and therefore not surprisingly many commonly dysregulated lncRNAs are shared across these conditions. In the interest of translational research only those lncRNAs which are strongly conserved have been addressed. The lncRNAs discussed here have diverse roles in regulating inflammation, proliferation, migration, invasion and apoptosis. Understanding the molecular basis of lncRNA function in rheumatology will be crucial in fully determining the inflammatory mechanisms that drive these conditions.

Role of long non-coding RNAs in adipogenesis: State of the art and implications in obesity and obesity-associated diseases

Obesity is an evolutionary, chronic, and relapsing disease that consists of a pathological accumulation of adipose tissue able to increase morbidity for high blood pressure, type 2 diabetes, metabolic syndrome, and obstructive sleep apnea in adults, children, and adolescents. Despite intense research over the last 20 years, obesity remains today a disease with a complex and multifactorial etiology. Recently, long non-coding RNAs (lncRNAs) are emerging as interesting new regulators as different lncRNAs have been found to play a role in early and late phases of adipogenesis and to be implicated in obesity-associated complications onset. In this review, we discuss the most recent advances on the role of lncRNAs in adipocyte biology and in obesity-associated complications. Indeed, more and more researchers are focusing on investigating the underlying roles that these molecular modulators could play. Even if a significant number of evidence is correlation-based, with lncRNAs being differentially expressed in a specific disease, recent works are now focused on deeply analyzing how lncRNAs can effectively modulate the disease pathogenesis onset and progression. LncRNAs possibly represent new molecular markers useful in the future for both the early diagnosis and a prompt clinical management of patients with obesity.

Roles of long non‑coding RNA in osteoarthritis (Review)

Osteoarthritis (OA) is a chronic bone and joint disease characterized by articular cartilage degeneration and joint inflammation and is the most common form of arthritis. The clinical manifestations of OA are chronic pain and joint activity disorder, which severely affect the patient quality of life. Long non-coding RNA (lncRNA) is a class of RNA molecules >200 nucleotides long that are expressed in animals, plants, yeast, prokaryotes and viruses. lncRNA molecules lack an open reading frame and are not translated into protein. The present review collated the results of recent studies on the role of lncRNA in the pathogenesis of OA to provide information for the prevention, diagnosis and treatment of OA.

Involvements of long noncoding RNAs in obesity-associated inflammatory diseases

Obesity is associated with chronic low-grade inflammation that affects the phenotype of multiple tissues and therefore is implicated in the development and progression of several age-related chronic inflammatory disorders. Importantly, a new family of noncoding RNAs, termed long noncoding RNAs (lncRNAs), have been identified as key regulators of inflammatory signalling pathways that can mediate both pretranscriptional and posttranscriptional gene regulation. Furthermore, several lncRNAs have been identified, which are differentially expressed in multiple tissue types in individuals who are obese or in preclinical models of obesity. In this review, we examine the evidence for the role of several of the most well-studied lncRNAs in the regulation of inflammatory pathways associated with obesity. We highlight the evidence for their differential expression in the obese state and in age-related conditions including insulin resistance, type 2 diabetes (T2D), sarcopenia, osteoarthritis and rheumatoid arthritis, where obesity plays a significant role. Determining the expression and functional role of lncRNAs in mediating obesity-associated chronic inflammation will advance our understanding of the epigenetic regulatory pathways that underlie age-related inflammatory diseases and may also ultimately identify new targets for therapeutic intervention.

Differentially-expressed mRNAs, microRNAs and long noncoding RNAs in intervertebral disc degeneration identified by RNA-sequencing

The underlying molecular mechanisms of intervertebral disc degeneration (IDD) remain unclear. This study aimed to identify the crucial molecules and explore the function of noncoding RNAs and related pathways in IDD. We randomly selected three samples each from an IDD and a spinal cord injury group (control) for RNA-sequencing. We identified 463 differentially-expressed long noncoding RNAs (lncRNAs), 47 differentially-expressed microRNAs (miRNAs), and 1,334 differentially-expressed mRNAs in IDD. Three hundred fifty-eight lncRNAs as cis-regulators could potentially target 865 genes. Protein–protein interaction (PPI) network analysis confirmed that IL-6, VEGFA, IGF1, MMP9, CXCL8, FGF2, IL1B, CCND1, ITGAM, PTPRC, FOS and PTGS2 were hub genes. We built a competing endogenous RNA (ceRNA) network and identified lncRNA XIST–hsa-miR-4775–PLA2G7 and lncRNA XIST–hsa-miR-424-5p–AMOT/TGFBR3 ceRNA axes. Quantitative real-time PCR (qRT-PCR) was implemented in 15 IDD samples and 15 controls to validate differentially-expressed genes in ceRNA axes. From the ceRNA network, gene ontology (GO) enrichment analysis indicated that noncoding RNAs were associated with several biological processes, including extracellular matrix organization, extracellular structure organization, leukocyte migration, and mesenchyme development. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis revealed that noncoding RNAs were associated with several pathways including the AGE-RAGE signaling pathway, PI3K-Akt signaling pathway, axon guidance, and osteoclast differentiation. These results indicate that some specific noncoding RNAs and ceRNA axes may be vital during the development of IDD, and may have potential as alternative diagnostic biomarkers as well as novel therapeutic strategies for IDD.

Comprehensive analysis of lncRNA and mRNA based on expression microarray profiling reveals different characteristics of osteoarthritis between Tibetan and Han patients

Background

Osteoarthritis (OA) is thought to be the most prevalent chronic joint disease, especially in Tibet of China. Here, we aimed to explore the integrative lncRNA and mRNA landscape between the OA patients of Tibet and Han.

Methods

The lncRNA and mRNA expression microarray profiling was performed by SurePrint G3 Human Gene Expression 8x60K v2 Microarray in articular cartilage samples from OA patients of Han nationality and Tibetans, followed by GO, KEGG, and trans-regulation and cis-regulation analysis of lncRNA and mRNA.

Results

We found a total of 117 lncRNAs and 297 mRNAs differently expressed in the cartilage tissues of Tibetans (n = 5) comparing with those of Chinese Han (n = 3), in which 49 lncRNAs and 158 mRNAs were upregulated, and 68 lncRNAs and 139 mRNAs were downregulated. GO and KEGG analysis showed that several unreported biological processes and signaling pathways were particularly identified. LncRNA-mRNA co-expression analysis revealed a remarkable lncRNA-mRNA relationship, in which OTOA may play a critical role in the different mechanisms of the OA progression between Tibetans and Chinese Han.

Conclusion

This study identified different lncRNA/mRNA expression profiling between OA patients of Tibetans and Han, which were involved in many characteristic biological processes and signaling pathways.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13018-021-02213-y.

Understanding osteoarthritis pathogenesis: a multiomics system-based approach

PURPOSE OF REVIEW

Osteoarthritis is a heterogeneous, multifactorial condition regulated by complex biological interactions at multiple levels. Comprehensive understanding of these regulatory interactions is required to develop feasible advances to improve patient outcomes. Improvements in technology have made extensive genomic, transcriptomic, epigenomic, proteomic, and metabolomic profiling possible. This review summarizes findings over the past 20 months related to omics technologies in osteoarthritis and examines how using a multiomics approach is necessary for advancing our understanding of osteoarthritis as a disease to improve precision osteoarthritis treatments.

RECENT FINDINGS

Using the search terms 'genomics' or 'transcriptomics' or 'epigenomics' or 'proteomics' or 'metabolomics' and 'osteoarthritis' from January 1, 2018 to August 31, 2019, we identified advances in omics approaches applied to osteoarthritis. Trends include untargeted whole genome, transcriptome, proteome, and metabolome analyses leading to identification of novel molecular signatures, cell subpopulations and multiomics validation approaches.

SUMMARY

To address the complexity of osteoarthritis, integration of multitissue analyses by multiomics approaches with the inclusion of longitudinal clinical data is necessary for a comprehensive understanding of the disease process, and for appropriate development of efficacious diagnostics, prognostics, and biotherapeutics.

Cell cycle-related lncRNAs and mRNAs in osteoarthritis chondrocytes in a Northwest Chinese Han Population

BACKGROUND

A group of differentially expressed long non-coding RNAs (lncRNAs) have been shown to play key roles in osteoarthritis (OA), although they represented only a small proportion of lncRNAs that may be biologically and physiologically relevant. Since our knowledge of regulatory functions of non-coding RNAs is still limited, it is important to gain better understanding of their relation to the pathogenesis of OA.

METHODS

We performed mRNA and lncRNA microarray analysis to detect differentially expressed RNAs in chondrocytes from three OA patients compared with four healthy controls. Then, enrichment analysis of the differentially expressed mRNAs was carried out to define disease molecular networks, pathways and gene ontology (GO) function. Furthermore, target gene prediction based on the co-expression network was performed to reveal the potential relationships between lncRNAs and mRNAs, contributing an exploration of a role of lncRNAs in OA mechanism. Quantitative RT-PCR analyses were used to demonstrate the reliability of the experimental results.

FINDINGS

Altogether 990 lncRNAs (666 up-regulated and 324 down-regulated) and 546 mRNAs (419 up-regulated and 127 down-regulated) were differentially expressed in OA samples compared with the normal ones. The enrichment analysis revealed a set of genes involved in cell cycle. In total, 854 pairs of mRNA and lncRNA were highly linked, and further target prediction appointed 12 genes specifically for their corresponding lncRNAs. The lncRNAs lncRNA-CTD-2184D3.4, ENST00000564198.1, and ENST00000520562.1 were predicted to regulate SPC24, GALM, and ZNF345 mRNA expressions in OA.

INTERPRETATION

This study uncovered several novel genes potentially important in pathogenesis of OA, and forecast the potential function of lnc-CTD-2184D3.4, especially for the cell cycle in the chondrocytes. These findings may promote additional aspects in studies of OA.

Interplay between genetics and epigenetics in osteoarthritis

Research into the molecular genetics of osteoarthritis (OA) has been substantially bolstered in the past few years by the implementation of powerful genome-wide scans that have revealed a large number of novel risk loci associated with the disease. This refreshing wave of discovery has occurred concurrently with epigenetic studies of joint tissues that have examined DNA methylation, histone modifications and regulatory RNAs. These epigenetic analyses have involved investigations of joint development, homeostasis and disease and have used both human samples and animal models. What has become apparent from a comparison of these two complementary approaches is that many OA genetic risk signals interact with, map to or correlate with epigenetic mediators. This discovery implies that epigenetic mechanisms, and their effect on gene expression, are a major conduit through which OA genetic risk polymorphisms exert their functional effects. This observation is particularly exciting as it provides mechanistic insight into OA susceptibility. Furthermore, this knowledge reveals avenues for attenuating the negative effect of risk-conferring alleles by exposing the epigenome as an exploitable target for therapeutic intervention in OA.

RNA sequencing reveals the circular RNA expression profiles of osteoarthritic synovium

The role of circular RNAs (circRNAs) in regulating cartilage homeostasis in osteoarthritis (OA) has been reported. However, the regulatory mechanisms of circRNAs in OA synovium remains basically unidentified. The current study intended to divulge the expression profile of circRNAs in OA synovium and investigate the possible molecular mechanisms of circRNAs in synovitis in OA through an integrated bioinformatics analysis. A total of 35 synovium samples were collected, including 17 from patients with knee OA and 18 from controls. circRNA sequencing was then carried out on five OA synovium samples as well as five controls to explore the expression pattern of the circRNAs. Real-time quantitative reverse transcription-polymerase chain reaction (qRT-PCR) was done to confirm the manifestation levels of six differentially expressed circRNAs. Gene Ontology (GO) as well as Kyoto Encyclopedia of Genes and Genomes pathway (KEGG) enrichment analyses were done for differentially expressed circRNAs using the DAVID database to annotate the functions. The circRNA-miRNA coexpression network was then created to estimate the probable molecular regulatory mechanisms of specifically expressed circRNAs in OA synovium. Total of 122 circRNAs were found to be differentially expressed in OA synovium through RNA sequencing. The expressions of five downregulated circRNAs as well as an upregulated circRNA were confirmed through the use of qRT-PCR. The circRNA-miRNA network was created to annotate the probable molecular regulatory mechanisms of specifically expressed circRNAs. Our outcomes revealed that circRNAs might be incorporated in the initiation as well as development of OA synovitis and might have prospective importance in OA diagnosis and therapy.

Identify CRNDE and LINC00152 as the key lncRNAs in age-related degeneration of articular cartilage through comprehensive and integrative analysis

Background

Osteoarthritis (OA) is one of the most important age-related degenerative diseases, and the leading cause of disability and chronic pain in the aging population. Recent studies have identified several lncRNA-associated functions involved in the development of OA. Because age is a key risk factor for OA, we investigated the differential expression of age-related lncRNAs in each stage of OA.

Methods

Two gene expression profiles were downloaded from the GEO database and differentially expressed genes (DEGs) were identified across each of the different developmental stages of OA. Next, gene ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to annotate the function of the DEGs. Finally, a lncRNA-targeted DEG network was used to identify hub-lncRNAs.

Results

A total of 174 age-related DEGs were identified. GO analyses confirmed that age-related degradation was strongly associated with cell adhesion, endodermal cell differentiation and collagen fibril organization. Significantly enriched KEGG pathways associated with these DEGs included the PI3K–Akt signaling pathway, focal adhesion, and ECM–receptor interaction. Further analyses via a protein–protein interaction (PPI) network identified two hub lncRNAs, CRNDE and LINC00152, involved in the process of age-related degeneration of articular cartilage. Our findings suggest that lncRNAs may play active roles in the development of OA. Investigation of the gene expression profiles in different development stages may supply a new target for OA treatment.