- Research article
- Open Access
Comprehensive analysis of lncRNA and mRNA based on expression microarray profiling reveals different characteristics of osteoarthritis between Tibetan and Han patients
Journal of Orthopaedic Surgery and Research volume 16, Article number: 133 (2021)
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.
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.
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.
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.
Osteoarthritis (OA) is the most common skeletal disease worldwide, causing pain, dyskinesia for patients, and a severe burden to society . The incidence of osteoarthritis is rising because of the aging population and the epidemic of obesity . Pain and loss of function are the main clinical characteristics that lead to treatment containing non-pharmacological, pharmacological, and surgical approaches . Clinicians recognize that the diagnosis of osteoarthritis may be too late to expect much help from disease-modifying drugs . Despite efforts over the past decades to excavate markers of disease, still-imaging procedures and biochemical marker analyses need to be improved and possibly extended with more specific and sensitive methods to reliably describe disease processes, to diagnose the disease earlier, and to follow the course of disease and treatment effectiveness [5, 6]. Hence, the investigation and exploration of the biological mechanism and biomarkers of OA are critical to understanding the OA pathophysiology. China is the country with the largest incidence of OA, such as Kashin Beck Disease (KBD) in the world, and the Tibetans are the people most affected by OA in China . In particular, in the Qing-Tibet Plateau, a selenium-deficient region, the prevalence of OA is serious and still increasing and continues to damage public health for Tibetans . More importantly, due to the specific climate and living environment and life habits in Tibet, it exhibits particular characteristics of OA between Tibetans and other populations, including the Han nationality . In Tibet, where severe selenium deficiency is endemic, iodine deficiency is a risk factor for OA . Besides, the long-term or severe cold stimulation and the mutual labor weaken the metabolism and immune defense of articular cartilage, leading to cartilage surface ulceration and damage inflammation of OA patients in Tibet [11, 12]. Thus, the investigation of the genetic and biological difference in OA progression between Tibetans and the Han nationality will benefit in understanding the different mechanisms of the development of OA between Tibetans and Chinese Han. However, the exploration in this field is still extremely limited.
Over the last decade, the development of high-throughput technologies has allowed an in-depth examination of the non-coding genome with unprecedented resolution and scale, in which long non-coding RNAs (lncRNAs) play crucial roles in many physiological and pathological processes . But their role in biological function and the signaling pathway involved in OA remain unknown. It has been reported that lncRNAs are expressed aberrantly in OA . Current evidence indicates that lncRNAs not only serve as positive or negative regulators of OA but also crosstalk with multiple potential targets to impact the critical events in the OA process [15, 16]. The identification of lncRNAs in OA patients reveals its diagnostic and therapeutic value [17,18,19]. However, the comprehensive landscape of lncRNA in OA development remains unclear. Besides, aberrant mRNA expression provides us with an opportunity to better diagnose and stratify the OA via established panels or discover new approaches . Analysis of mRNA expression identifies multiple markers in OA patients and may be the key contributing factor to OA . Gene expression microarray has been the primary biomarker platform ubiquitously applied in biomedical research, leading to enormous data, predictive models, and biomarkers accrued [22, 23]. Genetic profiling of the mRNA and lncRNA expression can be used for staging OA at the molecular level . However, the comprehensively integrated analysis of lncRNA and mRNA in OA of Tibetans and Chinese Han based on gene expression profiling remains unreported.
In the present study, we performed genome-wide screening for mRNA and lncRNA profiles in OA patients of Tibetans and Han nationality. Our study aimed to characterize the different features of lncRNA and mRNA expression between Tibetans and Han nationality. We also focused on determining the biological processes and pathways that exhibited significant changes in this event. Strikingly, we identified OTOA as a critical factor among the interaction of lncRNAs and mRNAs in the different lncRNA-mRNA networks of OA development between Tibetans and Chinese Han. Our comprehensive analysis will provide valuable information for understanding the variant molecular mechanisms in the OA progression of Tibetans and Han nationality.
Materials and methods
Clinical cartilage sample collection
A total of 8 human articular cartilage samples were obtained from consented patients undergoing joint replacement surgery due to OA, including 3 Han nationality patients and 5 Tibetan patients. All the samples were snapped frozen in liquid nitrogen within 2 h of surgery and stored at 80 °C before RNA extraction. All cartilage was washed extensively in phosphate-buffered saline (PBS) to avoid contamination. All study subjects provided informed consent before the inclusion of people in the study. The study protocol was approved by the Qinghai Provincial People’s Hospital, China.
The processes of cartilage harvest, sectioning, grinding, and RNA extraction were performed as previously described [25, 26]. Briefly, articular cartilage samples were sectioned and powdered under liquid nitrogen; 100 mg of articular cartilage powder was used for RNA isolation with 5 ml of Trizol (Invitrogen, CA). The RNA concentration and quality (RNA integrity number (RIN) and 28S/18S ratio) were determined by a NanoDrop (NanoDrop Technologies, DE) and the RNA 6000 Nano Assay on an Agilent 2100 Bioanalyzer (Agilent Technologies, CA), respectively.
LncRNA and mRNA microarray expression profiling
The lncRNA and mRNA microarray expression profiling was performed in articular cartilage samples from OA patients of Han nationality and Tibetans. Briefly, total RNA (100 ng) was labeled with the mRNA Complete Labeling and Hyb Kit (Agilent Technologies) and hybridized on the SurePrint G3 Human Gene Expression 8x60K v2 Microarray (Agilent Technologies). The microarray contains 50,599 probes for 32,776 human mRNA and 17,438 human lncRNAs, which are derived from authoritative databases, including RefSeq, Ensemble, GenBank, and the Broad Institute. After hybridization and washing, processed slides were scanned with the Agilent G2505C microarray scanner (Agilent Technologies). Raw data were extracted using Feature Extraction (version10.7.1.1; Agilent Technologies). Next, quantile normalization and subsequent data processing were done through the Genespring software (version 12.0; Agilent Technologies). The microarray profiling was conducted in Beijing Capitalbio Technology company.
The analysis of differential expression level of mRNA and lncRNA from microarray
The differential expression level of mRNA and lncRNA from microarray was analyzed by using a limma R package. After quantile normalization, raw signals from microarray were log2 transformed. Differential expression of an mRNA or lncRNA was defined by absolute value of fold change (FC) > 2 (|log2FC| > 1) and P value < 0.05 (Student’s t test).
The differentially expressed mRNAs were submitted to the DAVID database (https://david.ncifcrf.gov) to be classified into different Gene Ontology  including molecular function (MF), biological process (BP) and cellular component (CC), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation groups.
LncRNA-mRNA co-expression analysis
The differential expression level of mRNA and lncRNA from microarray was analyzed by using a limma R package. After quantile normalization, raw signals from microarray were log2 transformed. Differential expression of an mRNA or lncRNA was defined by absolute value of fold change (FC) > 2 (|log2FC| > 1) and P value < 0.05 (Student’s t test).
The relationship between the differently expressed lncRNAs and mRNAs in Tibetan OA samples was assessed by Pearson correlation analysis. We calculated the Pearson correlation coefficient (PCC) of differentially expressed lncRNAs and mRNAs to unearth the co-expressed lncRNA and mRNAs. The absolute value of PCC > 0.7, as well as P < 0.05, was considered statistically significant.
Based on the co-expression, we further searched for the mechanism of how these aberrant lncRNAs realize functions through cis- or trans-regulating mRNAs. When the correlation of lncRNAs and mRNAs was significant (PCC > 0.7) and the mRNA loci were within 300 k windows up- and downstream of the given lncRNAs, we identified them as cis-regulated mRNAs of the corresponding lncRNAs. LncRNAs involved in certain biological pathways may interact with transcription factors (TFs) . Thus, we performed the trans-regulation analysis to identify the co-expressed mRNAs significantly related to the dysregulated lncRNAs by using a threshold (PCC > 0.7, distance larger than 300 k) based on the Cytoscape software.
OA clinical samples of Tibetans demonstrate significantly altered lncRNA and mRNA expression patterns compared with that of Han nationality
It has been reported that the development of OA in Tibet is prevalent, and many risk factors are involved in the OA progression of Tibetans . Several lncRNAs have now been identified as being either differentially expressed in OA tissue or as candidate central regulators for OA progression . Previous studies highlighted the importance of aberrant expression of lncRNA and mRNA in both the pathogenesis and potential treatment of OA . Besides, the microarray analysis is an effective tool for the investigation of OA [31, 32]. Hence, we performed gene expression profiling based on the microarray analysis in the articular cartilage tissues of Tibetans and the Chinese Han. Interestingly, the volcano plot filtering has found total 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 in Tibetan samples (|log2 fold change | > 1, and p < 0.05) (Fig. 1 a and b). On the other hand, 68 lncRNAs and 139 mRNAs were found to be downregulated in the cartilage tissues of Tibetans (n = 5) comparing with those of Chinese Han (n = 3) (fold change > 2, and P < 0.05) (Fig. 1 a and b), suggesting that the different contributed factors are involved in the development of OA between Tibetans and the Han nationality.
Many key genes and pathways are involved in the variant pathogenesis of OA between Tibetans and Chinese Han
To further explore the potential function of the differently expressed genes, we performed Gene Ontology  and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis based on the above 297 identified mRNAs. Surprisingly, GO analysis revealed that the functions of these mRNAs were related to many important processes such as muscle cell proliferation, cell-cell adhesion, and chemokine ligand 2 production (Fig. 2 a and Table S1), implying that the different cellular processes may play specific roles in OA progression of Tibetans and Chinese Han. Significantly, KEGG analysis demonstrated that these mRNAs participated in several key signaling pathways, including the Hippo signaling pathway, Wnt signaling pathway, and Ras signaling pathway (Fig. 2b and Table S2), indicating the particular landscape of signaling pathway of OA in Tibetan and Han nationality.
OTOA plays a crucial role in the different lncRNA-mRNA networks of OA progression between Tibetan and Chinese Han
The interaction of lncRNA and mRNA significantly affects the progression of many diseases, such as cancer and inflammation [33, 34]. It has been reported that the integrated network analysis of lncRNA-mRNA identifies the biomarkers in patients with OA of the knee . Hence, we further explored the co-expression status of the above variant lncRNAs and mRNAs of OA in Tibetan and Chinese Han. Fundamentally, the Pearson correlation analysis showed the significant relationship of expression between the lncRNAs and mRNAs in Tibetan and Chinese Han samples, implying that the lncRNA-mRNA interaction may contribute to different clinical OA progression of Tibetans and Chinese Han (Fig. 3a) As we all know, lncRNA is able to modulate its targeted genes in the manner of cis-regulation and trans-regulation . When the mRNA loci are within 300 k windows up- and downstream of the given lncRNAs, we can identify them as cis-regulated mRNAs of the corresponding lncRNAs . Besides, lncRNAs involved in certain biological pathways may interact with transcription factors (TFs) . The information of lncRNA-mediated cis- and trans-regulation is crucial for our understanding of RNA-mediated gene regulation and genome complexity . Therefore, we tried to explore the cis-regulation and trans-regulation networks of lncRNA-mRNA in the system. Interestingly, we identified the top 100 significantly interacted lncRNA-mRNA in the trans-regulation analysis, in which the OTOA was the most remarkable gene (Fig. 3b).
Then, we further analyzed the lncRNA-mRNA cis-regulation networks based on the threshold of the Pearson correlation coefficient ≥ 0.7 and P value < 0.05. Surprisingly, we identified lnc-SDIM1 and SDIM1, lnc-CYB5R2 and CYB5R2, lnc-CCDC67 and CCDC67, and lnc-OTOA and OTOA with the Pearson correlation coefficients of 0.805, 0.833, 0.833, and 0.857, respectively (Fig. 4a, b). Moreover, lnc-OTOA and OTOA was the most significant lncRNA-mRNA relationship among these interactions (Fig. 4a, b). Taken together, we conclude that the OTOA, which was elevated in the cartilage tissues of Tibetans with OA compared to those in Chinese Han, was the most critical gene in the variant lncRNA-mRNA networks of OA progression in Tibetan and Han nationality. OTOA might be served as the potential biomarkers to characteristic the difference between the OA of Tibetans and Han nationality.
OA is the most frequent musculoskeletal disease, leading to functional decline and loss of quality of life . Clinically, the disease is featured by joint pain, tenderness, crepitus, stiffness, limitation of movement with occasional effusion, and variable degrees of local inflammation . The OA is mostly prevalent in the western regions of China, especially in Sichuan and Qinghai provinces and Tibet autonomous regions . Despite various preventive and control measures, new cases still have been emerging in Tibet due to the specific living environment and life habits of Tibetan [43,44,45]. Unfortunately, the biological investigation and different processes of the development of OA between Tibetans and Han nationality is extremely limited and unclear. It has been reported that the single-nucleotide polymorphism (SNPs) rs7033979 in the ASPN gene may play a protective role in susceptibility to OA, and genes COL10A1 and HABP2 may play a role in the risk of developing OA in the Tibetan population [46, 47]. LncRNAs are involved in many diseases, and the exquisite cell type specificity of lncRNA function has clear implications for targeted therapy . The differential expression of long noncoding RNAs by bioinformatics method in human knee OA may as novel dignosis biomarkers and therapeutic targets . In addition, lncRNAs such as HOTAIR, MEG2, MALAT1, and DILC are involved in the modification of OA development [50,51,52,53]. In addition, recent studies have examined the expression of extracellular genomic material, including lncRNA and mRNA, which are aberrantly expressed in the body fluids of OA patients . In this study, we identified a total of 117 lncRNAs and 297 mRNAs differently expressed in the cartilage tissues of Tibetans comparing with those of Chinese Han. We observed 49 lncRNAs and 158 mRNAs were upregulated, and 68 lncRNAs and 139 mRNAs were downregulated in the cartilage tissues of Tibetans. These data provide new evidence that the different modulation mechanisms may be involved in the Tibetan and Han nationality OA development. The differently expressed lncRNA and mRNA may be the discriminative biomarkers for the OA of Tibetan and Chinese Han.
Previous GO enrichment analysis showed that 119 upregulated genes were significantly enriched in blood vessel development, and KEGG pathway enrichment showed that genes were involved in the circadian rhythm pathway . It has been reported that some risk pathways such as TGF-β signaling pathway and inflammatory bowel disease and risk biological processes such as chondrocyte differentiation and regulation of cartilage development are involved in the OA pathology . Functional analysis of lncRNA in the synovial membrane of osteoarthritis patients revealed that lipopolysaccharide, angiogenesis, tumor necrosis factor (TNF) signaling, and mitogen-activated protein kinase (MAPK) signaling pathways might contribute to OA development . Interestingly, our GO and KEGG analysis showed that several unreported processes such as muscle cell proliferation, cell-cell adhesion, and chemokine ligand 2 production and signaling pathways including Hippo signaling pathway, Wnt signaling pathway, and Ras signaling pathway were particularly identified in Tibetan OA patients compared to the Chinese Han OA patients. Our data revealed some different biological changes which were not similar to previous studies. These data suggest that the development of Tibetan and Chinese Han OA may feature particularly direct or indirect biological progression, and the specific mechanism is needed to further investigate.
It has been reported that the lncRNA-mRNA interaction is crucial for many pathological processes including OA . CRNDE and LINC00152 play active roles in the development of OA . The deregulated genes including USP46, CPVL, FKBP5, FOSL2, GADD45B, PTGS1, ZNF423, ADAMTS1, and TFAM might be involved in the pathology of OA . Surprisingly, we demonstrated a significant relationship of expression between the lncRNAs and mRNAs in Tibetan OA samples compared to that of Chinese Han. It indicates that the c-expression of lncRNA and mRNA may be involved in the variant regulation mechanism of the OA progression of Tibetans and Han nationality. Strikingly, our trans-regulation analysis demonstrated the top 100 significantly interacted lncRNA-mRNA, in which the OTOA was the most remarkable gene. Moreover, the cis-regulation analysis further identified that lnc-OTOA and OTOA was the most significant lncRNA-mRNA relationship. These data imply that OTOA may play a critical role in the different lncRNA-mRNA networks between Tibetans and Han nationality OA progression. The effect of OTOA on the OA development of Tibetans and Han nationality and its molecular mechanism are needed to further explore. Moreover, it has been well-identified that microRNAs (miRNAs) play critical roles in the pathogenesis of OA by interacting with lncRNAs and various targeted genes from the aspects of diagnosis, physiopathology, and therapy . It will be crucial and beneficial to investigate the interplay of lncRNA, miRNA, and mRNA networks in the development of OA, which may need to be explored further in the future. Furthermore, the different lncRNAs and mRNAs involved in multiple critical cellular biological processes and signaling identified in this study may be helpful for the exploration of diagnosis biomarkers for Tibetan OA patients. The significant lncRNA-mRNA association of lnc-OTOA and OTOA may serve as the potential diagnostic or therapeutic for the Tibetan OA patients in the clinic.
In conclusion, this study demonstrated that 117 lncRNAs and 297 mRNAs are differently expressed in Tibetan OA patients compared to the Chinese Han, which are involved in many characteristic biological processes and signaling pathways. LncRNA-mRNA co-expression analysis identified OTOA as the critical gene in the variant lncRNA-mRNA networks of OA progression between Tibetans and Han nationality. Our finding provides new insights into the differential lncRNA-mRNA scenario of OA development in Tibetan and Han nationality, improving the understanding of the discriminative mechanism of OA progression between Tibetan and Han nationality. Consequently, this should be helpful in shedding novel light on the development of nationality-specific therapeutic for OA.
Availability of data and materials
The datasets used during the current study are available from the corresponding author on reasonable request. RNA was extracted from 3 Han nationality patients and 5 Tibetan patients with OA.
Kyoto Encyclopedia of Genes and Genomes
Kashin Beck disease
Pearson correlation coefficient
Tumor necrosis factor
Mitogen-activated protein kinase
Bijlsma JW, Berenbaum F, Lafeber FP. Osteoarthritis: an update with relevance for clinical practice. Lancet. 2011;377(9783):2115–26.
Barnett R. Osteoarthritis. Lancet. 2018;391(10134):1985.
He T, Wu D, He L, Wang X, Yang B, Li S, Chen Y, Wang K, Chen R, Liu B, et al. Casein kinase 1 epsilon facilitates cartilage destruction in osteoarthritis through JNK pathway. FASEB J. 2020;34(5):6466–78.
Anakwe RE, Middleton SD. Osteoarthritis at the base of the thumb. BMJ. 2011;343:d7122.
Appleton CT. Osteoarthritis year in review 2017: biology. Osteoarthritis Cart. 2018;26(3):296–303.
Huang Z, Chen J, Li B, Zeng B, Chou CH, Zheng X, Xie J, Li H, Hao Y, Chen G, et al. Faecal microbiota transplantation from metabolically compromised human donors accelerates osteoarthritis in mice. Ann Rheum Dis. 2020;79(5):646–56.
Huang Q, Zhou ZK, Ma J, Li Y, Yang X, Shen B, Yang J, Kang PD, Pei FX. The arthropathic and functional impairment features of adult Kashin-Beck disease patients in Aba Tibetan area in China. Osteoarthritis Cartil. 2015;23(4):601–6.
Zhao ZJ, Li Q, Yang PZ, Wang H, Kong LC, Wang LH, Sun LY. Selenium: a protective factor for Kaschin-Beck disease in Qing-Tibet Plateau. Biol Trace Elem Res. 2013;153(1-3):1–4.
Chasseur C, Suetens C, Nolard N, Begaux F, Haubruge E. Fungal contamination in barley and Kashin-Beck disease in Tibet. Lancet. 1997;350(9084):1074.
Moreno-Reyes R, Suetens C, Mathieu F, Begaux F, Zhu D, Rivera MT, Boelaert M, Neve J, Perlmutter N, Vanderpas J. Kashin-Beck osteoarthropathy in rural Tibet in relation to selenium and iodine status. N Engl J Med. 1998;339(16):1112–20.
Hinsenkamp M, Mathieu F, Claus W, Collard JF, de Maertelaer V. Effects of physical environment on the evolution of Kashin-Beck disease in Tibet. Int Orthop. 2009;33(4):1085–8.
Mathieu F, Suetens C, Begaux F, De Maertelaer V, Hinsenkamp M. Effects of physical therapy on patients with Kashin-Beck disease in Tibet. Int Orthop. 2001;25(3):191–3.
Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172(3):393–407.
Cen X, Huang XQ, Sun WT, Liu Q, Liu J. Long noncoding RNAs: a new regulatory code in osteoarthritis. Am J Transl Res. 2017;9(11):4747–55.
Fathollahi A, Aslani S, Jamshidi A, Mahmoudi M. Epigenetics in osteoarthritis: novel spotlight. J Cell Physiol. 2019;234(8):12309–24.
Chen WK, Yu XH, Yang W, Wang C, He WS, Yan YG, Zhang J, Wang WJ. lncRNAs: novel players in intervertebral disc degeneration and osteoarthritis. Cell Prolif. 2017;50(1):e12313.
Sun D, Li C, Liu J, Wang Z, Liu Y, Luo C, Chen Y, Wen S. Expression profile of microRNAs in hypertrophic cardiomyopathy and effects of microRNA-20 in inducing cardiomyocyte hypertrophy through regulating gene MFN2. DNA Cell Biol. 2019;38(8):796–807.
Wu Y, Lu X, Shen B, Zeng Y. The therapeutic potential and role of miRNA, lncRNA, and circRNA in osteoarthritis. Curr Gene Ther. 2019;19(4):255–63.
Jiang SD, Lu J, Deng ZH, Li YS, Lei GH. Long noncoding RNAs in osteoarthritis. Joint Bone Spine. 2017;84(5):553–6.
Watt FE. Osteoarthritis biomarkers: year in review. Osteoarthritis Cartil. 2018;26(3):312–8.
Zhang R, Guo H, Yang X, Li Z, Zhang D, Li B, Zhang D, Li Q, Xiong Y. Potential candidate biomarkers associated with osteoarthritis: evidence from a comprehensive network and pathway analysis. J Cell Physiol. 2019;234(10):17433–43.
Su Z, Fang H, Hong H, Shi L, Zhang W, Zhang W, Zhang Y, Dong Z, Lancashire LJ, Bessarabova M, et al. An investigation of biomarkers derived from legacy microarray data for their utility in the RNA-seq era. Genome Biol. 2014;15(12):523.
Nimse SB, Sonawane MD, Song KS, Kim T. Biomarker detection technologies and future directions. Analyst. 2016;141(3):740–55.
Zhong L, Huang X, Karperien M, Post JN. Correlation between Gene Expression and Osteoarthritis Progression in Human. Int J Mol Sci. 2016;17(7):1126.
Chou CH, Lee MT, Song IW, Lu LS, Shen HC, Lee CH, Wu JY, Chen YT, Kraus VB, Wu CC. Insights into osteoarthritis progression revealed by analyses of both knee tibiofemoral compartments. Osteoarthritis Cartil. 2015;23(4):571–80.
Xu Y, Barter MJ, Swan DC, Rankin KS, Rowan AD, Santibanez-Koref M, Loughlin J, Young DA. Identification of the pathogenic pathways in osteoarthritic hip cartilage: commonality and discord between hip and knee OA. Osteoarthritis Cartil. 2012;20(9):1029–38.
Chen S, Dong C, Qian X, Huang S, Feng Y, Ye X, Miao H, You Q, Lu Y, Ding D. Microarray analysis of long noncoding RNA expression patterns in diabetic nephropathy. J Diab Comp. 2017;31(3):569–76.
Moreno-Reyes R, Suetens C, Mathieu F, Begaux F, Zhu D, Rivera T, Boelaert M, Neve J, Perlmutter N, Vanderpas J. Kashin-Beck disease and iodine deficiency in Tibet. Int Orthop. 2001;25(3):164–6.
Pearson MJ, Jones SW. Review: long noncoding RNAs in the regulation of inflammatory pathways in rheumatoid arthritis and osteoarthritis. Arthritis Rheumatol. 2016;68(11):2575–83.
Jeffries MA. Osteoarthritis year in review 2018: genetics and epigenetics. Osteoarthritis Cartil. 2019;27(3):371–7.
Kuttapitiya A, Assi L, Laing K, Hing C, Mitchell P, Whitley G, Harrison A, Howe FA, Ejindu V, Heron C, et al. Microarray analysis of bone marrow lesions in osteoarthritis demonstrates upregulation of genes implicated in osteochondral turnover, neurogenesis and inflammation. Ann Rheum Dis. 2017;76(10):1764–73.
de Kroon LMG, van den Akker GGH, Brachvogel B, Narcisi R, Belluoccio D, Jenner F, Bateman JF, Little CB, Brama PAJ, Blaney Davidson EN, et al. Identification of TGFbeta-related genes regulated in murine osteoarthritis and chondrocyte hypertrophy by comparison of multiple microarray datasets. Bone. 2018;116:67–77.
Chen Y, Lin Y, Bai Y, Cheng D, Bi Z. A long noncoding RNA (lncRNA)-associated competing endogenous RNA (ceRNA) network identifies eight lncRNA biomarkers in patients with osteoarthritis of the knee. Med Sci Monit. 2019;25:2058–65.
Han BW, Ye H, Wei PP, He B, Han C, Chen ZH, Chen YQ, Wang WT. Global identification and characterization of lncRNAs that control inflammation in malignant cholangiocytes. BMC Genomics. 2018;19(1):735.
Fang L, Wang H, Li P. Systematic analysis reveals a lncRNA-mRNA co-expression network associated with platinum resistance in high-grade serous ovarian cancer. Invest New Drugs. 2018;36(2):187–94.
Wang Y, Wu T, Zou L, Xiong L, Zhang T, Kong L, Xue Y, Tang M. Genome-wide identification and functional analysis of long non-coding RNAs in human endothelial cell line after incubation with PM2.5. Chemosphere. 2019;216:396–403.
Liu TL, Fan XC, Li YH, Yuan YJ, Yin YL, Wang XT, Zhang LX, Zhao GH. Expression profiles of mRNA and lncRNA in HCT-8 cells infected with Cryptosporidium parvum IId subtype. Front Microbiol. 2018;9:1409.
Chen R, Jiang T, She Y, Xie S, Zhou S, Li C, Ou J, Liu Y. Comprehensive analysis of lncRNAs and mRNAs with associated co-expression and ceRNA networks in C2C12 myoblasts and myotubes. Gene. 2018;647:164–73.
Yan P, Luo S, Lu JY, Shen X. Cis- and trans-acting lncRNAs in pluripotency and reprogramming. Curr Opin Genet Dev. 2017;46:170–8.
Pereira D, Ramos E, Branco J. Osteoarthritis. Acta Med Port. 2015;28(1):99–106.
Stoffer MA, Smolen JS, Woolf A, Ambrozic A, Berghea F, Boonen A, Bosworth A, Carmona L, Dougados M, de Wit M et al: Development of patient-centred standards of care for osteoarthritis in Europe: the eumusc.net-project. Ann Rheum Dis 2015, 74(6):1145-1149.
Sun LY, Yuan LJ, Fu Y, Deng JY, Wang LH. Prevalence of Kaschin-Beck disease among Tibetan children in Aba Tibetan and Qiang Autonomous Prefecture: a 3-year epidemiological survey. World J Pediatr. 2012;8(2):140–4.
Guo Y, Li H, Yang L, Li Y, Wei B, Wang W, Gong H, Guo M, Nima C, Zhao S, et al. Trace element levels in scalp hair of school children in Shigatse, Tibet, an endemic area for Kaschin-Beck disease (KBD). Biol Trace Elem Res. 2017;180(1):15–22.
Moreno-Reyes R, Mathieu F, Boelaert M, Begaux F, Suetens C, Rivera MT, Neve J, Perlmutter N, Vanderpas J. Selenium and iodine supplementation of rural Tibetan children affected by Kashin-Beck osteoarthropathy. Am J Clin Nutr. 2003;78(1):137–44.
Schepman K, Engelbert RH, Visser MM, Yu C, de Vos R. Kashin Beck disease: more than just osteoarthrosis: a cross-sectional study regarding the influence of body function-structures and activities on level of participation. Int Orthop. 2011;35(5):767–76.
Tai Z, Huang L, Lu F, Shi Y, Ma S, Cheng J, Lin H, Liu X, Li Y, Yang Z. Association study of candidate genes for susceptibility to Kashin-Beck disease in a Tibetan population. BMC Med Genet. 2017;18(1):69.
Yang Y, Zhao Q, Liu Y, Liu X, Chu Y, Yan H, Fan Y, Huo S, Wang L, Lou Q, et al. FRZB1 rs2242070 polymorphisms is associated with brick tea type skeletal fluorosis in Kazakhs, but not in Tibetans, China. Arch Toxicol. 2018;92(7):2217–25.
Liu SJ, Horlbeck MA, Cho SW, Birk HS, Malatesta M, He D, Attenello FJ, Villalta JE, Cho MY, Chen Y, et al. CRISPRi-based genome-scale identification of functional long noncoding RNA loci in human cells. Science. 2017;355(6320).
Xiao K, Yang Y, Bian Y, Feng B, Li Z, Wu Z, Qiu G, Weng X. Identification of differentially expressed long noncoding RNAs in human knee osteoarthritis. J Cell Biochem. 2019;120(3):4620–33.
Huang W, Ding X, Ye H, Wang J, Shao J, Huang T. Hypoxia enhances the migration and invasion of human glioblastoma U87 cells through PI3K/Akt/mTOR/HIF-1alpha pathway. Neuroreport. 2018;29(18):1578–85.
Zhang Y, Wang F, Chen G, He R, Yang L. LncRNA MALAT1 promotes osteoarthritis by modulating miR-150-5p/AKT3 axis. Cell Biosci. 2019;9:54.
Xu J, Xu Y. The lncRNA MEG3 downregulation leads to osteoarthritis progression via miR-16/SMAD7 axis. Cell Biosci. 2017;7:69.
Huang J, Liu L, Yang J, Ding J. Xu X: lncRNA DILC is downregulated in osteoarthritis and regulates IL-6 expression in chondrocytes. J Cell Biochem. 2019;120(9):16019–24.
Budd E, Nalesso G, Mobasheri A. Extracellular genomic biomarkers of osteoarthritis. Expert Rev Mol Diagn. 2018;18(1):55–74.
Rao ZT, Wang SQ, Wang JQ. Exploring the osteoarthritis-related genes by gene expression analysis. Eur Rev Med Pharmacol Sci. 2014;18(20):3056–62.
Li L, Dong L, Zhao D, Gao F, Yan J. Classical dendritic cells regulate acute lung inflammation and injury in mice with lipopolysaccharide induced acute respiratory distress syndrome. Int J Mol Med. 2019;44(2):617–29.
Shui X, Xie Q, Chen S, Zhou C, Kong J, Wang Y. Identification and functional analysis of long non-coding RNAs in the synovial membrane of osteoarthritis patients. Cell Biochem Funct. 2020;38(4):460–71.
Xiang S, Li Z, Bian Y, Weng X. Identification of changed expression of mRNAs and lncRNAs in osteoarthritic synovium by RNA-sequencing. Gene. 2019;685:55–61.
Cui J, Zhang Q, Song Q, Wang H, Dmitriev P, Sun MY, Cao X, Wang Y, Guo L, Indig IH, et al. Targeting hypoxia downstream signaling protein, CAIX, for CAR T-cell therapy against glioblastoma. Neuro Oncol. 2019;21(11):1436–46.
Zhang X, Bu Y, Zhu B, Zhao Q, Lv Z, Li B, Liu J. Global transcriptome analysis to identify critical genes involved in the pathology of osteoarthritis. Bone Joint Res. 2018;7(4):298–307.
Oliviero A, Della Porta G, Peretti GM, Maffulli N. MicroRNA in osteoarthritis: physiopathology, diagnosis and therapeutic challenge. Br Med Bull. 2019;130(1):137–47.
This work was supported by the Chinese Natural Scientific Grant #81460022, Qinghai Natural Scientific Grant 2016-ZJ-747, and Hunan Provincial Innovation Fund S2017SFYLJS0070 from Junming Luo; Chinese Natural Scientific Grant #81760396 from Zhongshu Shan.
Ethics approval and consent to participate
All procedures performed in studies involving human participants were approved by the Qinghai Provincial People’s Hospital, and in accordance with the ethical standards of the institutional research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards.
Consent for publication
Informed consent was obtained from all individual participants included in the study.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Luo, J., Luo, X., Duan, Z. et al. Comprehensive analysis of lncRNA and mRNA based on expression microarray profiling reveals different characteristics of osteoarthritis between Tibetan and Han patients. J Orthop Surg Res 16, 133 (2021). https://doi.org/10.1186/s13018-021-02213-y
- Microarray profiling