Integrated Analysis of lncRNA and mRNA Expression Profiles Indicates Age-Related Changes in Meniscus

Little has been known about the role of long non-coding RNA (lncRNA) involves in change of aged meniscus. Microarray analyses were performed to identify lncRNAs and mRNAs expression profiles of meniscus in young and aging adults and apple bioinformatics methods to analyse their potential roles. The differentially expressed (DE) lncRNAs and mRNAs were confirmed by qRT-PCR. A total of 1608 DE lncRNAs and 1809 DE mRNAs were identified. Functional and pathway enrichment analyses of all DE mRNAs showed that DE mRNAs were mainly involved in the TGF-beta, Wnt, Hippo, PI3K-Akt signaling pathway. The expressions of TNFRSF11B and BMP2 were significantly upregulated in aging group. LASSO logistic regression analysis of the DE lncRNAs revealed four lncRNAs (AC124312.5, HCG11, POC1B-AS1, and AP001011.1) that were associated with meniscus degradation. CNC analysis demonstrated that AP001011 inhibited the expression of TNFRSF11B and AC1243125 upregulated the expression of TNFRSF11B. CeRNA analysis suggested that POC1B-AS1 regulates the expression of BMP2 by sponging miR 130a-3p, miR136-5p, miR 18a-3p, and miR 608. Furthermore, subcellular localization and m6A modification sites prediction analysis of these four lncRNAs was performed. These data lay a foundation for extensive studies on the role of lncRNAs in change of aged meniscus.


INTRODUCTION
As a crucial part of the tibiofemoral articulation and an articular cartilage of the knee joint, the meniscus protects the underlying articular cartilage, contributing to normal function of the knee joint and protecting the articular surfaces (Noble and Turner, 1986;Huang et al., 2019). Meniscus injuries include traumatic tears caused by sports injury and age-related degenerative tears. Degeneration of the meniscus further increases the risk of traumatic tears. In recent years the incidence of meniscus injury is on rise probably due to the general aging of the population. Aging is a time-dependent decline of physiological abilities. Physiological degeneration is a major risk factor for age-related diseases such as type 2 diabetes, cancer, and dementia (Degirmenci and Lei, 2016). Many changes occur in articular cartilages with aging, including dysregulated expression of related genes, deterioration of mechanical properties and degenerative changes in the morphology (Sacitharan and Vincent, 2016). Age-related tissue changes, such as myxoid degeneration, collagen fiber thickening, and accumulation of advanced glycation end products, ultimately increases the tissue vulnerability of the meniscus (Pauli et al., 2011). Knee cartilage thickness is also negatively correlated with age. All these changes aggravate the course of osteoarthritis, increase the incidence of meniscal tears in the elderly population, and seriously affect people's quality of life. However, there are few studies on meniscus aging, and key genes involved in meniscal aging are still remaining to be identified.
Long non-coding RNAs (lncRNAs), which do not have protein-coding functions, have received increasing research attention in recent years (Quinn and Chang, 2016). LncRNAs play key roles in many biological processes. Studies have shown that lncRNAs play important roles in the aging process of tissues and organs. For example, the expression of lncRNA H19 is downregulated during endothelial cell senescence. LncRNA H19 controls senescence, angiogenic sprouting, proliferation, and inflammatory activation of endothelial cell via inhibiting STAT3 activation (Hofmann et al., 2019). Su et al. found 88 significantly upregulated lncRNAs and 46 significantly downregulated lncRNAs in the cochlea of aged C57BL/6 mice. LncRNA AW112010 alleviated age-related hearing loss by promoting mitochondrial biogenesis to maintain mitochondrial function and hair cell survival in cochlear hair cells (Su et al., 2019). Moreover, dysregulated lncRNAs were involved in meniscus degeneration. Interleukin-1 β (IL-1β) treatment of degenerative meniscus resulted in significant alterations in 375 mRNAs, 15 miRNAs, 56 lncRNAs and 90 circRNAs. LncRNA LOC107986251 is an important player in meniscal degeneration (Jiang et al., 2021). Zhao et al. performed an analysis on the GSE98918 dataset, which included 24 meniscus samples and related clinical data. The expressions of 208 lncRNAs and mRNAs were significantly differential between patients with knee osteoarthritis and non-knee osteoarthritis group. In the knee osteoarthritis group, a total of five lncRNAs were significantly associated with meniscus degeneration (Zhao et al., 2020a). However, until now, the role of lncRNAs in age-related change of meniscus is poorly known.
In order to investigate the molecular mechanisms of age-related change of meniscus and further explore potential biologic therapeutic targets, we performed microarray analysis of meniscus to determine lncRNAs and mRNAs expression profiles in young and aging adults. Then, we carried out bioinformatics analysis of differentially expressed molecules to identify key molecules and functional pathways in the aged meniscus tissues.

Human meniscus Sample Collection
The torn segment of the injured meniscus was resected from eight patients who underwent arthroscopic meniscectomy surgery with no significant cartilage degeneration, osteoarthritis, rheumatoid arthritis and other metabolic diseases. In the partial meniscectomy, in addition to the torn tissue, a small amount of normal tissue was removed in order to make the cutting edge smooth, and these tissues were used for this study. For the current study, these tissues were divided into two categories based on patients' age: young group (n = 4, 20-25 years old) and aging group (n = 4, 50-55 years old). The study was approved by Peking University Third Hospital Ethical Committee and a written informed consent was obtained by all study subjects (IRB00006761-M2019299).

Histology Staining
The meniscus samples were fixed for 24-48 h in 4% paraformaldehyde and dehydrated in graded ethanol prior to paraffin embedding. Sections were then cut in 5 μm using a microtome (Lecia, Bannockburn, IL, United States) and mounted on slides. Hematoxylin-eosin (HE) staining was performed to verify the histological differences between the young and aging groups (Figure 1). A 5x magnification setting was used to show the integral meniscus sections. A 10x magnification was utilized to highlight the remaining extracellular martix conditions of two groups. In addition, a Rodeo score system was used to assess tissue histological features of two groups from cellularity, predominant cell type, collagen organization, and matrix morphology (Rodeo et al., 2000;Longo et al., 2013). The total histological score ranges from 0 to 6 points. A lower final score indicates a significantly higher rate of fibrocartilage changes (Mesiha et al., 2007).

Microarray Hybridization and Analysis
Total RNA was extracted and purified from aging group (n = 4) and young group (n = 4) for lncRNA and mRNA, respectively. TRIzol Reagent (Invitrogen, Waltham, MA, United States) and RNasey Mini Kit (Qiagen, German) were used for RNA extract and purify. The quantity and purity of total RNA samples were measured by NanoDrop ND-1000 (ThermoFisher, Waltham, MA, United States). Human lncRNA & mRNA microarray v5.0 (Arraystar, 8 × 60K, Rockville, MD, United States) were used for detection of differential expression of lncRNAs and mRNAs. The threshold set for differentially expressed (DE) lncRNAs and mRNAs was p < 0.05, Fold change ≥2.0 or fold change ≤0.5. All the microarray hybridization and analysis were carried out by Aksomics (Shanghai, China).

Quantitative Real-Time PCR
Quantitative real-time PCR (qRT-PCR) experiments were performed using uperScriptTM III Reverse Transcriptase (Invitrogen, Waltham, MA, United States) and 2 × SYBR Green PCR Master Mix (Arraystar, Rockville, MD, United States) according to the manufacturer's instructions, with β-actin as the internal control. Details of the genes and primers were listed in Supplementary Tables S1, S2.

Protein and Protein Interaction Analysis
Protein and protein interaction (PPI) analysis of all DE mRNAs was performed based on STRING database (https://string-db.org/) with minimum required interaction score ≥0.9. Textmining, experiments, databases, coexpression, neighborhood, gene fusion and co-occurrenc were active interaction sources. The PPI sub network was constructed by Cytoscape (v3.8.2).

Least Absolute Shrinkage and Selection Operator Logistic Regression Analysis to Identify Target Differentially Expressed LncRNAs
The DE lncRNAs were narrowed down by using the least absolute shrinkage and selection operator (LASSO) regression. The LASSO algorithm was applied with the Stata software version 15.0 "lassopack" package. Lasso regression is used for eliminating automated variables and the selection of features. The Lasso penalty shrinks or reduces the coefficient value towards zero. The less contributing variable is therefore allowed to have a zero or near-zero coefficient (Chen and Chen, 2008). The top 20 up-regulated and top 20 down-regulated DE LncRNAs were selected with the largest foldchange for LASSO regression.

LncRNAs-mRNAs Co-Expressed Analysis
Nineteen DE mRNAs co-expressed with all DE lncRNAs were identifified to reveal the potential regulatory relationships between DE lncRNAs and DE mRNAs. The targeted DE mRNAs and DE lncRNA with Pearson correlation coefficients values |r|≥0.94 and p-values of no less than 0.001 were filtered out to construct coding and non-coding co-expression (CNC) network. Also, 3 DE lncRNAs from LASSO logistic regression analysis and 3 DE mRNAs were screened out to constructed sub CNC network. Both of them were constructed by Cytoscape (v3.8.2).

Bioinformatic Analysis of Targete Differentially Expressed LncRNAs
to predict the m6A modification sites with combined score ≥0.6 (High confidence and Very high confidence).

Statistical Analysis
Statistical analyses were carried out with the GraphPad Prism v5.0 package (GraphPad Sofware, Inc., La Jolla, CA, United States). Data were expressed as mean ± standard error of the mean (SEM), and Student's t test (Mann-Whitney U) was used to determine the diference between two groups. p < 0.05 was considered significant differences. Table 1 summarizes the basic clinical parameters of patients. HE staining showed that the young menisci exhibit a smooth and intact surface, while the aged menisci display a rough surface with fibrillation. Microscopically, young meniscus tissue showed a well-organized collagen fiber network and normal fibrous matrix. Besides, evenly distributed meniscus cells can be seen between the collagen bundles ( Figure 1A). However, extracellular matrix (ECM) mucoid degeneration, disruption of collagen fiber, disarrangement of native cells and decrease in cell density were observed in the aged meniscus tissue. This demonstrated ECM structural damage and cellular disarrangement. Rodeo score system was used to was used to assess tissue histological features of young and aging groups (n = 5). Aging group had a significantly lower (p < 0.05) Rodeo score compared with young group. This scoring data are in agreement with the HE staining tissue features.

LncRNA and mRNA Expression Profiles Between Young and Aging Group
To understand the effects of age on the expression of lncRNAs and mRNAs in meniscus tissues, Arraystar Human lncRNA and mRNA Array (V5.0) was used to detect meniscus tissues of young and aging groups (n = 4/group). Using fold change ≥2 or fold change ≤0.5 and p < 0.05 as the threshold of differential expression, a total of 1608 DE lncRNAs (1,183 upregulated lncRNAs, 425 downregulated lncRNAs) and 1809 DE mRNAs (1,592 upregulated mRNAs, 217 downregulated mRNAs) were detected. Violcano plot was conducted to explore signal intensity of DE lncRNAs and DE mRNAs. Hierarchical clustering was performed to display the distinguishable DE lncRNAs and DE mRNAs expression pattern among samples ( Figure 2). The data have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE191157 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc= GSE191157).

Functional and Pathway Enrichment
Analyses of all Differentially Expressed mRNAs GO analyses of 1809 DE mRNAs were performed. The top 10 BPs, CCs and MFs based on their p-values are shown in Figure 3A. The most enriched GO terms were Gene expression (BP), Nucleoplasm (CC), poly (A) RNA binding (MF).
Pathway analyses of all DE mRNAs were also performed to detect enriched pathways associated with cartilage damage in different databases. As a result, TGF-beta signaling pathway, FGF signaling pathway and Wnt signaling pathway were enriched from PANTHER database ( Figure 3B); TGF-beta signaling pathway, PI3K-Akt signaling pathway, Th17 cell differentiation, mTOR signaling pathway were enriched from KEGG-pathway database ( Figure 3C); Signaling by hedgehog, PTEN regulation, MAPK family signaling cascades, signaling by WNT were enriched from Reactome database ( Figure 3D). At last, a total of 18 DE mRNAs were obtained by intersecting DE mRNAs with pathways detected in the three databases ( Figure 3E). In the 18 DE mRNAs, FGF6 and PPP2R5C were significantly downregulated, and the others were significantly upregulated. Our pathway analyses showed that the 18 DE mRNAs were mainly involved in the TGF-beta signaling pathway, Wnt signaling pathway, Hippo signaling pathway, PI3K-Akt signaling pathway ( Figure 3F; Table 2).
Our analysis of all disease-related DE mRNAs using the David database and kobas V3.0 revealed several disorders related to articular injury and aging, including osteoarthritis, osteoporosis, BMD, musculoskeletal disorders, and skeletal disease. At the same time, six related molecules were screened out ( Figure 3G). The expressions of COL1A1, COL1A2, COL6A1 and MMP2 were significantly downregulated, and the expressions of TNFRSF11B and BMP2 were significantly upregulated ( Table 2).

qRT-PCR Validation of Target Differentially Expressed mRNAs
A total of 18 and 6 disease-associated DE mRNAs were uncovered by pathway analyses and enrichment analyses, respectively. At last we produced 23 disease-related DE mRNAs except for one   Figure 4A). The gene expression abundance of the 19 mRNAs is shown in Figure 4B.

Protein-Protein Interaction Analysis
PPI analysis of all DE mRNAs showed 1,724 nodes and 1,743 edges in the map (average node degree = 2.02, PPI enrichment p = 0.00043) ( Figure 5A). The results of qRT-PCR validation showed that the expression changes of 19 DE mRNAs in 23 disease-related DE mRNAs were consistent with the microarray results. After that, the 19 DE mRNAs that were consistent with the microarray results and their interacting proteins were used to construct the sub PPI network, generating 86 nodes, 115 edges (average node degree = 2.775) ( Figure 5B).

Least Absolute Shrinkage and Selection Operator Logistic Rregression Analysis to Identify Target Differentially Expressed LncRNAs
To further narrow down the varieties and identify the target DE lncRNAs, we selected the top 20 up-regulated and top 20 downregulated DE lncRNAs with the largest fold change for LASSO regression analysis. As the tuning parameter lambda changes, the corresponding coefficients of the DE lncRNAs were reduced to zero. Finally, 4 lncRNAs (3 upregulated DE lncRNAs including lncRNA AC124312.5, lncRNA HCG11 and lncRNA POC1B-AS1,1 downregulated DE lncRNA including lncRNA AP001011.1) were screened out for further analysis ( Figure 6A). Relative expression in microarray of 4 DE lncRNAs was showed in Figure 6B.

Coding and Noncoding Co-Expression Analysis
We used CNC analysis to found lncRNAs related to meniscus damage. CNC analysis was performed on 19 DE mRNAs that were consistent with the microarray results and all DE lncRNAs to find the lncRNAs that regulate these 19 DE mRNAs. Based on | r | >0.94, p < 0.001, FDR <0.05, a total of 632 nodes (19 DE mRNAs and 613 DE lncRNAs) and 1525 DE mRNA-DE lncRNA pairs (1,356 positive correlation, 169 negative correlation) were obtained ( Figure 7A). Based on the CNC network analysis, 3 target DE lncRNAs in LASSO regression analysis (lncRNA AC124312.5, lncRNA HCG11 and lncRNA AP001011.1) were filtered out. We further constructed a CNC sub network based on 3 DE lncRNAs and 3 DE mRNAs ( Figure 7B).

Competing Endogenous RNA Analysis
Based on the ceRNA mechanism, a total of 98 nodes (8 DE mRNAs, 44 miRNAs, 46 DE lncRNAs) and 229 edges were predicted using the 19 DE mRNAs that were consistent with microarray results and all DE lncRNAs. We constructed ceRNA network via Cytoscape V3.8.2. (Figure 8A). In the ceRNA network, we found the lncRNA POC1B-AS1 obtained by the cable regression analysis. We choose lncRNA POC1B-AS1 and its miRNA binding sites and corresponding target genes to construct a ceRNA sub network ( Figure 8B).

Bioinformatic Analysis of 4 Differentially Expressed LncRNAs
In this work, we used m6A modification site predictor SRAMP to predict the m6A modification site of 4 DE LncRNAs. The result showed that lncRNA HCG1 had a high confidence m6A modification site in position 556 with combined score 0.657 ( Figure 9A). lncRNA POCIB-AS1 had thirteen m6A modification sites including two very high confidence and eleven high confidence m6A modification sites ( Figure 9B). And lncRNA AP001011.1 had two high confidence m6A modification sites in position 314 and 368 with combined score 0.602 and 0.609, respectively. However, no m6A modification sites were predicted out on lncRNA AC124312.5 ( Figure 9C). Then, we used two sequence-based predictors called "lncLocater" and "iLoc-LncRNA" to predict the subcellular locations of 4 DE LncRNAs. The result showed that lncRNA HCG1 was located in nucleus and the probability score was 0.8159 and 0.7755 in two database, respectively. lncRNA POCIB-AS1 was located in cytoplasm with probability score 0.9034 and 0.8496, respectively. The localization of lncRNA AC124312.5 and lncRNA AP001011.1 was not only in nucleus by lncLocater with probability score 0.5861 and 0.5005, but also in cytoplasm by iLoc-LncRNA with probability score 0.6971 and 0.9848 ( Figure 9D).

DISCUSSION
Age-related degeneration of the meniscus is a pathological process that occurs at the tissue, cell, and molecular level of the meniscus. It is manifested as significant morphological and physiological changes of the meniscus, and finally manifested as a decrease in its load-bearing capacity. Since aging is an inevitable process, finding mechanisms to protect the meniscus tissue from avoidable damage has become the focus of research. In order to understand the influence of age factors on meniscus tissues, we analyzed the changes of lncRNA and mRNA expression profiles in meniscus tissues of young and aging people by microarray analysis. A total of 1,608 lncRNAs and 1809 mRNAs have undergone significant changes.
We first analyzed the function of DE mRNA. After verification by qRT-RCR, we found the expression of metalloproteinase-2 (MMP-2) and type I α 2 the expression of colagen (COL1A2) were significantly decreased in the aging group. ECM of meniscus plays a key role in regulating cell activity due to non-vascular nature of the adult meniscus. The ECM is composed of various macromolecules, mainly including types I and III collagens, as well as degradative components such as matrix MMPs (Fourniere et al., 2021). Meniscus ECM plays an important functional role in regulation of the physicochemicall and mechanical environments at multiple scales, which affects overall joint health, cartilage homeostasis, and chondrocyte metabolism (Wilusz et al., 2014). MMPs degrades structural components of ECM, such as collagen (Amaral et al., 2020). MMPs are directly related to almost every biological process that involves ECM remodeling. In the aging process, changes in collagen structure and disorder of MMP will affect ECM remodeling, leading to a variety of diseases, such as vascular diseases, rheumatoid arthritis, osteoarthritis, fibrosis and cancer (Panwar et al., 2018). In many animal models and human diseases, imbalances in MMP activity are associated with myocardial aging status, even in the absence of cardiovascular disease (Meschiari et al., 2017). In aging mouse myocardium and cartilage, increased MMP activity is associated with increased inflammation, ECM deposition, and decreased angiogenesis (Meschiari et al., 2017;Sebastian et al., 2020). Dysregulated ECM leads to myocardial fibrosis and cartilage ECM degradation, which results in dysfunction at the cellular, extracellular, and whole organ levels. MMP-mediated ECM degradation is the core process of intervertebral disc degeneration. In nucleus pulposus cells, miR-874-3p targets and inhibits the expression of MMP2 and MMP3 (Song et al., 2021), and PART1 controls nucleus pulposus cells degeneration through the miR-93/ MMP2 pathway (Gao et al., 2020). The meniscus is a fibrocartilage and contains a lot of type I collagen. Type I collagen is a heterotrimer protein consisting of one α1 chains (COL1A1) and two α2 chain (COL1A2) (Karsenty and Park, 1995). COL1A2 is a major component of ECM. Some studies have reported that COL1A2 is involved in regulation of the osteoarthritis, Ehlers-Danlos syndrome, osteogenesis imperfecta . Our research showed that the mRNA expression of MMP2 and COL1A2 significantly decreased in the aging group. We believed that abnormal expression of MMP2 and COL1A2 in the aging group may contribute to ECM disorder. This might be one of the reasons that affect the structural and mechanical integrity of the meniscus.
During aging, the level of apoptosis has changed in many tissues (Tower, 2015). For example, excessive proliferation and decreased apoptosis of vascular smooth muscle cells in aged humans and mice may contribute to plaque formation (Qian et al., 2011). However, in other tissues, the increase in apoptosis leads to a decrease in the number of tissue cells, and the loss of cells increases the destruction of tissue homeostasis and dysfunction. The process of intervertebral disc degeneration is associated with higher rates of apoptosis and senescence. This results in cell loss in the nucleus pulposus tissue led to significant changes in the disc morphology and physiology which ultimately resulting in its decreased capability to bear compressive loads (Cazzanelli and Wuertz-Kozak, 2020). In our research, we found that the 18 DE mRNAs identified in the three databases are associated with cartilage injury and are mainly involved in the Hippo signaling pathway. The Hippo signaling pathway is a highly conserved pathway involved in tissue development and regeneration, which controls organ size by regulating cell proliferation and apoptosis (Ma et al., 2019;Meng et al., 2021). Activation of Hippo signaling pathway by overexpression of mammalian Ste-20 like kinase 1 (Mst1) or downregulation of YAP inhibits the development of Natural Killer T-Cell Lymphoma in nude mice xenograft models (Chang et al., 2019). During ischaemia-reperfusion injury or myocardial infarction, the Hippo-YAP pathway is activated and then activates Mst1 to increase caspase activation and cardiomyocyte apoptosis. Since the human heart lacks the ability to self-repair, loss of cardiomyocytes after injury such as myocardial infarction can lead to heart failure and death . Our results implicated that the Hippo pathway is involve in the aged meniscus. The Hippo pathway might be involved in meniscal tissue cell loss through apoptosis, leading to meniscal morphological and physiological changes. Currently, the Hippo pathway has emerged as a potential therapeutic target for controlling pathological cardiac remodelling. For example, the new inhibitor XMU-MP-1 can be used to pharmacologically inhibit Mst1/2, the core component of the Hippo pathway, to control the adverse effects of hypertrophy caused by pressure overload (Triastuti et al., 2019). Therefore, amelioration of meniscal degeneration by controlling the Hippo pathway may become a new research direction that deserves further investigation. In addition to ECM degradation and cell loss, inflammation also plays an important role in the degenerative process. The 18 DE mRNAs that we identified in the three databases were involved in TGF-beta signaling pathway. It was suggested that inflammation occurred during meniscus degeneration. However, the regulatory mechanism of TGF-beta signaling pathway during aging is still unclear. In recent years, studies have found that a large number of lncRNAs play a regulatory role in inflammatory responses, including TGF-beta signaling pathway (Hodgson et al., 2019). Tang et al. (2019) found that lncRNA-ATB promotes TGF-β induced glioma cell invasion through NF-κB and P38/ MAPK pathways. LncRNA H19 expression was downregulated in vitreous humour of patients with proliferative diabetic retinopathy. LncRNA H19 via TGF-β Influences endothelial mesenchymal transition during diabetic retinopathy (Thomas et al., 2019). LncRNA SNHG3 promotes proliferation and migration of non-small-cell lung cancer cell by activating TGF-β pathway and IL-6/JAK2/STAT3 pathway (Shi et al., 2020). Through CNC analysis of DE lncRNAs, we found that lncRNA AP001011 1, which is involved in modulating inflammation, Inhibits the expression of tumor necrosis factor receptor superfamily member 11B (TNFRSF11B), and AC124312.5 upregulates the expression of TNFRSF11B. TNFRSF11B is a member of tumor necrosis factor receptor super family (TNFRSF). It is a secreted for tumor necrosis factor (TNF)-related apoptosis inducing ligand (TRAIL) (Luan et al., 2020). The expression of TNFRSF11B is significantly upregulated in patients with osteoarthritis. The miR-145 plays an important role in chondrocyte proliferation and fibrosis by directly targeting TNFRSF11B, which inhibits the proliferation and fibrosis of OA (Wang et al., 2017). Upregulation of TNFRSF11B expression in injured osteoarthritic cartilage is a factor driving the transition of chondrocytes into osteoblasts in 3D vitro chondrogenesis mode (Ruiz et al., 2021). TNFRSF11B gene encodes osteoprotegerin. Due to the important role of osteoprotegerin in bone biology, TNFRSF11B gene is considered as a candidate gene for osteoporosis. Some studies believe that TNFRSF11B gene t245g (rs3134069) polymorphism is a risk factor for osteoporotic fracture (Boronova et al., 2015). However, there is still a lack of research on the role of TNFRSF11B in aged meniscus, especially its regulatory mechanism. Through CNC analysis, we found that AP001011 inhibited the expression of TNFRSF11B and AC1243125 upregulated the expression of TNFRSF11B, which provided a direction for further study on the mechanism of regulating TNFRSF11B during aging. Our CNC analysis also showed that lncRNA human leukocyte antigen complex group 11 (lncRNA HCG11) upregulated the expression of Smd5. It has been confirmed that lncRNA HCG11 is involved in regulating the progression and prognosis of various diseases. HCG11 inhibits growth and invasion of cervical cancer cell by sponging miR-942-5p and targeting growth factorindependent transcription repressor 1 (Zhang et al., 2020). Moreover, HCG11 promotes the expression of MMP13 and exacerbates osteosarcoma through sponging miR-579. However, the relationship between HCG11 and meniscus degradation has not been studied. According to our research, we proposed that HCG11 expression is significantly upregulated in the meniscus tissues in aging group, and Smad5 may be the target molecule of HCG11.
LncRNA act as miRNA sponges alleviating the inhibitory effects of miRNAs on target genes and increasing target gene amounts, known as ceRNA. Our ceRNA analysis indicated that lncRNA POC1B-AS1 regulates the expression of bone morphogenetic protein 2 (BMP2), MAPK14, Smad3 and TGFbeta receptor 2 (TGFBR2) through sponging miRNA. BMP2 is a superfamily member of TGF-β, it is one of the main chondrogenic growth factors involved in cartilage regeneration (Balaskas et al., 2020). Subcutaneous stem cell implantation studies have shown that BMP2 induces cartilage formation but also promotes endochondral ossification during ectopic bone/ cartilage formation (Zhou et al., 2016). Melatonin improves agerelated bone loss and osteoporosis by upregulating the expression of BMP2. Our findings suggested that lncRNA POC1B-AS1 regulates the expression of BMP2 by sponging miR130a-3p, miR136-5p, mir18a-3p, and miR608. However, the involvement of these miRNAs in meniscus degradation should be studied further.
In this work, we also predict the m 6 A modification site and the subcellular locations of 4 DE lncRNAs. The m 6 A modification is the most abundant internal modification in mRNAs and occurs also in lncRNAs . Alterations in RNA m 6 A methylation status can lead to cell dysfunction and disease. At present, m 6 A-modified lncRNAs have received extensive attention. M 6 A modification may control gene expression by regulating the translation efficiency and stability of lncRNAs (Coker et al., 2019). Shen et al. (2021) found that decreased m 6 A levels of lncRNA-XR_343955 in rat aortic tissue affect the inflammatory response through the cell adhesion molecule pathway during lipopolysaccharide-induced sepsis. The decrease in m 6 A modification of lncRNA XR_346,771 may be related to cation import in rat cardiac tissue during sepsis . Moreover, the function of lncRNA is closely related to their subcellular localization. In the nucleus, lncRNAs regulate gene expression at the epigenetic and transcriptional levels, and in the cytoplasm, they regulate gene expression at the posttranscriptional and translational levels (Zhao et al., 2020b). Our results showed that lncRNA HCG1 was located in nucleus and lncRNA POCIB-AS1 was located in cytoplasm. This indicated that the lncRNA HCG1 and lncRNA POCIB-AS1 have different functions. M 6 A modification site and the subcellular locations analysis provided a new direction for indepth study of these lncRNAs. In summary, this study investigated the lncRNAs and mRNAs expression profile of meniscus between young and aging group, and further analyzed regulatory relationships between lncRNAs and mRNAs that may be associated with change of aged meniscus. ECM degeneration, apoptosis, and inflammation were the focus of this study, especially the abnormal expression of lncRNAs and its regulatory role. At present, people have a high demand for new and targeted therapies methods that can delay the degeneration of the meniscus and reduce the risk of meniscus tear. Several possible approaches are currently being studied and tested. Among them, cell therapy, endogenous repair strategies, and treatments based on biological factors such as no coding RNA (miRNA, lncRNA et al.) are the most promising. In order to prevent or limit the impact of meniscal degeneration on articulation, we analyzed the molecular background of meniscus in the aging people used bioinformatics methods, providing directions for the development of new targeted therapies.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories