- 1State Key Laboratory of Oral Diseases and National Clinical Research Center for Oral Diseases, West China Hospital of Stomatology, Sichuan University, Chengdu, China
- 2Department of Oral Implantology, Shanghai Engineering Research Center of Tooth Restoration and Regeneration, School and Hospital of Stomatology, Tongji University, Shanghai, China
- 3Institute of Engineering Medicine, Beijing Institute of Technology, Beijing, China
- 4School of Basic Medical Sciences, Chengdu University, Chengdu, China
- 5Department of Orthodontics, China Shanghai Ninth People’s Hospital, College of Stomatology, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 6National Center for Stomatology, National Clinical Research Center for Oral Diseases, Shanghai Key Laboratory of Stomatology, Shanghai Jiao Tong University, Shanghai, China
Mechanical force, being so ubiquitous that it is often taken for granted and overlooked, is now gaining the spotlight for reams of evidence corroborating their crucial roles in the living body. The bone, particularly, experiences manifold extraneous force like strain and compression, as well as intrinsic cues like fluid shear stress and physical properties of the microenvironment. Though sparkled in diversified background, long noncoding RNAs (lncRNAs) concerning the mechanotransduction process that bone undergoes are not yet detailed in a systematic way. Our principal goal in this research is to highlight the potential lncRNA-focused mechanical signaling systems which may be adapted by bone-related cells for biophysical environment response. Based on credible lists of force-sensitive mRNAs and miRNAs, we constructed a force-responsive competing endogenous RNA network for lncRNA identification. To elucidate the underlying mechanism, we then illustrated the possible crosstalk between lncRNAs and mRNAs as well as transcriptional factors and mapped lncRNAs to known signaling pathways involved in bone remodeling and mechanotransduction. Last, we developed combinative analysis between predicted and established lncRNAs, constructing a pathway–lncRNA network which suggests interactive relationships and new roles of known factors such as H19. In conclusion, our work provided a systematic quartet network analysis, uncovered candidate force-related lncRNAs, and highlighted both the upstream and downstream processes that are possibly involved. A new mode of bioinformatic analysis integrating sequencing data, literature retrieval, and computational algorithm was also introduced. Hopefully, our work would provide a moment of clarity against the multiplicity and complexity of the lncRNA world confronting mechanical input.
Introduction
Mechanic force is recognized as a crucial factor regulating cellular physiological properties and directing cell fate (Thompson et al., 2012; Hao et al., 2015). Particularly, as one of the most active organs of human body undergoing consecutive remodeling, bone is exposed to a dynamic mechanic environment and shows dependence on mechanical cues for homeostasis maintenance (Thompson et al., 2012; Petridou et al., 2017). Forces like strain, fluid shear stress (FSS), compression, and microgravity (MG) have been reported to regulate bone reconstruction (Thompson et al., 2012). Remarkably, signaling pathways responding to force as well as the final impacts on cellular fate after force application differ among types, duration, intensity, and other parameters of the mechanic stimuli, which has already been detailed by us before (Hao et al., 2015). Though the whole picture of how components of bone function during the process is far from clear, what can be sure is the involvement of bone-related cells such as osteoblasts, osteocytes, osteoclasts, chondrocytes, and mesenchymal stem cells (MSCs) of various origins, in forming an entwined concert both spatially and temporally. Notably, periodontal ligament cells (PDLCs) and periodontal ligament stem cells (PDLSCs), as the specific cells seeding in jaws, turn out to be highly sensitive to force and responsible for bone remodeling during orthodontic tooth movement, which make them ideal model for mechanical study (Huang et al., 2018).
Yet with an obscure landscape so far, the ways cells sense and react to those fickle and perplexing mechanical cues are extensively explored. To begin with, an organism would depend on the so-called “mechanosensors” to perceive and convert the information into cellular ones. Focal adhesions (FAs) and cytoskeleton, as well as membrane channels which contribute to force sensing, are mostly explored. FAs play as the crucial initiator in detecting and transducing extracellular matrix (ECM)-related biomechanical properties. The FA complex functions by having the integrin heterodimer at the cell surface to interact with the extracellular information and subsequently putting docking proteins to connect F-actin bundles (Qin et al., 2020). Hence, the cues induce the contraction of cytoskeleton, contributing to the opening of nuclear pores by direct contact with nucleus via the linker of nucleoskeleton and cytoskeleton (LINC) (Bouzid et al., 2019). Piezo1/2, as the calcium (Ca2+)-permeable channels, was proven to be another mechanosensor essential for bone homeostasis. Via the critical intracellular Ca2+ sensor Ppp3ca, Piezo1/2 would contribute to the dephosphorylation and formation of the NFAT-Yap1-Ctnnb1 transcription factor (TF) network. The TFs were then activated and showed nuclear translocation for downstream gene regulation (Wang L et al., 2020; Zhou et al., 2020).
Once mechanical cues are transmitted into cells, cascade signaling pathways and diversified molecules would be involved in the response. In recent years, noncoding RNAs (ncRNAs) have gained the limelight for their potent regulation effects. NcRNAs hold a share of almost 98% of all genomic outputs in humans, contributing to an extensive manipulation of the transcriptional activities. MicroRNAs (miRNAs), as the most extensively studied type of ncRNAs, have gained wide recognition for their mechanosensitive roles, which were detailed before (Chen et al., 2017; Yuan et al., 2017). Nonetheless, a force-related study on long noncoding RNA (lncRNA) was far from completed. As transcribed RNA molecules longer than 200 nucleotides in length, lncRNAs could transform into unique secondary conformation, enabling its direct binding with DNA, RNA, and protein and indicating its significant roles (Mercer and Mattick, 2013). Meanwhile, lncRNA has also been attested to be a ubiquitous regulator in diversified biological processes via multiple function modes (Ørom and Shiekhattar, 2013; Cech and Steitz, 2014). Significantly, the competing endogenous RNA (ceRNA) mechanism accounts for a great part of lncRNA-oriented research where lncRNAs would act as the “sponge” to hijack miRNAs from their binding with mRNAs, thus relieving the degradation or decoy effects of miRNAs on the targets (Kartha and Subramanian, 2014; Tay et al., 2014). Alongside lncRNA’s well-known function modes, the upstream regulation of its transcription controlled by TFs is another eye-catching project. Intriguingly, lncRNAs and TFs could form entwined crosstalk and even contribute to a feedback loop, considering that TFs, as one particular type of mRNAs, can be regulated by ceRNA interactions in the posttranscriptional level as well as regulate gene transcription after its translocation to the nuclear (Ji et al., 2020). Thus, networks illustrating ceRNAs and TF–lncRNA interaction could provide us with an overview of the lncRNA world from both upstream regulation and downstream effects (Alam et al., 2014; Ji et al., 2020).
Moreover, there are well-documented compelling proofs that lncRNAs participate in the chondrogenic and osteogenic differentiation of MSCs and beyond, indicating the significance of lncRNAs in biophysical processes of bones (Huang et al., 2019; Ju et al., 2019; Sun et al., 2019; Wang H et al., 2020). Emerging evidence postulated that lncRNA could react to force condition and H19 was already proven to sponge miR-138 to regulate hMSC activity via FAK in a tension-induced way (Wu et al., 2018; Liu H. et al., 2019). A new TMSB4 pseudogene lncRNA related to mechanical stress (lncRNA-MSR) would hijack miR-152 to control TMSB4 expression in cyclic tensile strain (CTS)-induced regulation of chondrocytes and promote cartilage degradation (Liu et al., 2016). However, the systemic roles of lncRNAs in a mechanosensitive way remain obscure. Thus, aiming at highlighting an lncRNA-focused network, in which bone cells and their precursors would possibly activate to respond to the mechanic stimuli, we first explored a Gene Expression Omnibus (GEO) dataset for force-sensitive mRNA (FS mRNA) identification. Gene Ontology (GO) analysis was exploited for deeper knowledge of both the general and unique mechanisms that cells utilized for force responses. Instead of utilizing miRNA sequencing data which lack potent evidences, we screened out the creditable force-sensitive miRNAs (CFS miRNAs) based on reported research including qPCR and function experiment data to ensure reliability. Force-related lncRNA (FR lncRNAs) and ceRNA (FR ceRNA) networks were then predicted and constructed. Furthermore, to explicate the functional mode of FR lncRNAs, we displayed the direct interaction as well as the co-expression among FS mRNAs and FR lncRNAs. Crosstalk among TFs and lncRNAs by TF binding site (TFBS) prediction and identified several lncRNA-TF loops, demonstrating the potential upstream regulation on lncRNAs. Last, based on the analysis of our networks, we deciphered and summarized the potential mechanotransduction signaling pathways that FR lncRNAs may participate in, which were later integrated with the established data to build a pathway–lncRNA network. The network-centered strategy could highlight potential interaction and new roles of identified factors while suggesting a promising field worth exploring (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g001.jpg)
FIGURE 1. Workflow of this research. Orange squares stand for the database and evidences on which our research is based. Light and dark blue colors represent the crucial factors and networks predicted though bioinformatic strategy, respectively. Circles are used to highlight the quartet network of lncRNAs.
Methods and Materials
High-Throughput Sequencing Data Acquisition and Differentially Expressed Gene Analysis
Systematic searching of datasets concerning mechanical conditions was carried out in the GEO database (http://www.ncbi.nlm.nih.gov/geo/) for inclusion of FS mRNA. GSE112122 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE112122) was finally chosen based on its ideal cell type, experimental strategy, and data quality. Raw data of the gene expression profiles of GSE112122 were downloaded for analysis, and several grouping strategies were developed. Namely, intermittent force group analysis applied differentiated expression comparison based on data from control (GSM3058133-3058135) vs. intermittent force (GSM3058136-3058138). Similarly, static force group analysis utilized strategy as control (GSM3058139-3058141) vs. static force (GSM3058142-3058144). To explore the deviation of mRNA profiles caused by force types, intermittent force (GSM3058136-3058138) and static force (GSM3058142-3058144) are set as grouping modes. The Empirical Analysis of Digital Gene Expression Data in R (edgeR) package (Robinson et al., 2010; McCarthy et al., 2012) (http://www.bioconductor.org/packages/release/bioc/html/edgeR.html) was used to screen out differentially expressed mRNAs (DE mRNAs). Gene ID was converted to gene name and noted for gene type based on the Ensembl database (http://asia.ensembl.org/index.html). Only genes with type annotation as “protein coding” were selected for further analysis, which ruled out several noncoding RNA expression profiles with low confidence. The p-value was calculated with t-test and adjusted for multiple testing using the false discovery rate (FDR) method. A threshold was set as FDR < 0.05 and |log2foldChange (FC)| > 1 for DE mRNA identification. The Pretty Heatmaps (pheatmap) package was exploited to visualize the differences. Top 50 genes with the highest values of |log2foldChange (FC)| were chosen for heatmap creation. All analyses were performed with the R 3.5.0 framework. Intersection of DE mRNAs from intermittent force group and static force group analyses contributed to the list of FS mRNA. DE mRNAs gained from the comparison between profiles of intermittent force and static force were named as force-type-sensitive mRNAs (FTS mRNAs). VENNY 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html) was explored for visualization. Functional annotation for DE mRNAs was conducted based on GO analyses. All functional clustering was assessed by the clusterProfiler package (Yu et al., 2012) (http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html). The Enrichplot package (http://www.bioconductor.org/packages/release/bioc/html/enrichplot.html) and ggplot2 package were implemented to visualize the enrichment results. Significant GO terms were identified with a cutoff of p < 0.05.
Exploration and Pathway Analysis of CFS miRNAs
To identify reliable miRNAs which are sensitive to force among bone cells, we explored and collected the miRNAs which were verified to act in a mechanosensitive way via functional experiments or displayed certified different expression profiles in a force-related situation by qPCR. Cell types included in our work are confined to MSCs and their derivations in bone tissues. Research made before 2010 was excluded, as well as those without highly reliable evidence. CFS miRNAs identified in species other than Homo sapiens were further analyzed in miRBase (Kozomara et al., 2019) (http://www.mirbase.org/) for validation of their conservation among species, especially the similarity of sequences compared to their counterpart hsa-miRNAs. CFS miRNAs were uploaded into DIANA-mirPath v3.0 (Vlachos et al., 2015) (http://snf-515788.vm.okeanos.grnet.gr/) for KEGG analysis.
Construction of FR CeRNA Network and Predication of FR LncRNAs
Aiming at constructing the potential FR ceRNA network, we first screened out the lncRNA targets of CFS miRNAs, utilizing ENCORI (the Encyclopedia of RNA Interactomes) (Li et al., 2014) (http://starbase.sysu.edu.cn/) for prediction. The analysis was based on human genome and hg19 assembly, and the cutoff value was set as CLIP Data ≥ 3 (high stringency). Then, CFS-miRNA-targeting mRNAs were also predicted based on the intersection of datasets from miRTarBase (Chou et al., 2018) (http://mirtarbase.mbc.nctu.edu.tw), miRDB (Chen and Wang, 2020) (http://www.mirdb.org/), and TargetScan (Lewis et al., 2005) (http://www.targetscan.org/mamm_31/). MiRNA–mRNA pairs that gained validation from all three datasets were included. The predicted pairs were then intersected with the FS mRNA. Only mRNAs which were both force sensitive and predicted as miRNA targets were included to get CFS miRNA–FS mRNA pairs. Finally, the lncRNA–CFS miRNA pairs and CFS miRNA–FS mRNA pairs were imported and analyzed in Cytoscape 3.8.1 (Shannon et al., 2003) (https://cytoscape.org/). Pairs which could not form an lncRNA–miRNA–mRNA triplet were excluded. Alongside that, lncRNAs with a degree lower than three were also deleted. The ceRNA network was finally visualized by Cytoscape 3.8.1.
Validation of the Expression Changes of FR lncRNA
C3H/10T1/2 cells were purchased from Procell Life Science and Technology Co., Ltd. and were then cultured at 37°C under 5% carbon dioxide using modified Eagle’s medium (MEM, Invitrogen) supplemented with 10% fetal bovine serum (FBS, Gibco), 100 U/ml penicillin, and 100 μg/ml streptomycin (Invitrogen). The Flexcell FX-5000TM Tension System (Flexcell International Corporation) was used to apply mechanical cyclical stretch. A total of 5 × 104 cells were seeded onto the collagen type-I-coated BioFlex silicon plates per well for 24 h. Cells were then subjected to cyclical mechanical stimulation with 10% elongation, 0.5 Hz for 24 h. Control cells were cultured under the same conditions without stretch force. Once the programmed stimuli finished, the total RNA of cells was extracted by RNA-Quick Purification Kit (Esunbio, China). Then, 1 µg of RNA was transcribed to cDNA using a HiScript® III RT SuperMix for qPCR (+gDNA wiper) (Vazyme, China) according to the manufacturer’s protocol. The qRT-PCR assay was carried out on an ABI QuantStudio 6 system using Taq Pro universal SYBR qPCR Master Mix (Vazyme, China) according to the manufacturer’s protocol. GAPDH was used as endogenous references. The 2−ΔΔCt method was used to determine the relative quantification of gene expression levels. Primers for all the genes are listed in Table 1.
Identifying the FS mRNA–FR lncRNA Interaction and Co-expression
Interaction relationships between FS mRNA and FR lncRNA were screened out based on the database of ENCORI, which recorded the pairing relationship supported by at least one experiment. Loop relationships were extracted from the ceRNA network and a Ce-loop network was constructed. For co-expression network construction, the expression profile of GSE112122 was first normalized to TPM (transcripts per kilobase million) value, and the Pearson correlation coefficient was then calculated among FS mRNAs and FR lncRNAs. A threshold was set as a COR of >0.8 and a p-value of <0.05.
TF–LncRNA Regulation and Potential Feedback Circuit
Transcription activity is tightly regulated by TFs via their sequence-specific binding with the TFBSs on genome. We managed to map out the possible regulation of TFs on our FR lncRNAs for elucidation of the potential upstream mechanism. The prediction strategy was based on the position weight matrix (PWM) method (D’haeseleer, 2006). The upstream sequences of genes were downloaded from the University of California, Santa Cruz (UCSC) Genome Browser, and 1,000 bp were used for TFBS prediction of lncRNA upstream regulatory sites. The R 3.5.0 framework and TFBSTools package (Tan and Lenhard, 2016) (http://www.bioconductor.org/packages/release/bioc/html/TFBSTools.html) as well as JASPAR2016 package (Mathelier et al., 2016) (http://www.bioconductor.org/packages/release/data/experiment/html/JASPAR2016.html) were used for calculation, and the cutoff score was set as 95%. The TF–FR lncRNA regulation was later analyzed and visualized by Cytoscape 3.8.1. For possible bidirectional regulation between TFs and lncRNAs, we intersected the FR ceRNA as well as FTR ceRNA network with the TF list gained by TFBS prediction. TF–FR lncRNA pairs predicted to have a ceRNA relationship via miRNAs were then chosen. Triplets which could form the circle as TF–FR lncRNA–CFS miRNA–TF were finally included and visualized.
Paper-Based Pathway Mapping of FR lncRNAs
Based on the extensively reported mechanosensitive pathways as well as those essential for bone development, we matched lncRNAs in our networks to possible mechanisms based on their known effects. Specifically, compelling masters of mechanotransduction such as YAP/TAZ, Wnt/β-catenin, BMP/Smad, and MAPK are the main focus. Possible mechanisms of lncRNAs were analyzed based on the reported research in bone and beyond.
Identification of FS LncRNA and Pathway–LncRNA Network Construction
To corroborate our bioinformatics prediction, we reviewed the published articles concerning lncRNAs in bone cells undergoing simulated mechanical condition. Our searching strategy aimed to find research related to all kinds of mechanical condition, including tension, compression, fluid shear stress, vibration, microgravity, stiffness, topography, and so on. Besides, research irrelevant to the skeletal system was excluded. All data were updated to 2021/3/1. We named lncRNAs reported before as force-sensitive lncRNAs (FS lncRNAs). Furthermore, we built the pathway–lncRNA network integrating prediction results of FR lncRNA with verified data of FS lncRNA, which was then visualized by Cytoscape.
Results
Identification of FS and FTS mRNAs
As for the intermittent force group, 1,989 DE mRNAs were screened out, including 886 upregulated and 1,103 downregulated mRNAs. Additionally, 451 DE mRNAs were validated in the static force group, with 281 and 170 DE mRNAs found to be upregulated and downregulated, respectively. The discrepancy of the DE mRNA volumes between two force types conformed to the conclusion of the original study (Manokawinchoke et al., 2019). The intersection of DE mRNAs from intermittent and static force groups contributed to 404 genes in common, which were identified as FS mRNAs. Among them, 252 were promoted and 152 were inhibited after the mechanical loading regardless of types. Data analysis concerning force types contributed to 1,086 DE mRNAs, with 521 upregulated and 565 downregulated, which were defined as FTS mRNAs. Heatmap plots were used to assess differentially expressed genes, and a VENNY diagram was displayed for FS mRNA screening (Figures 2A–D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g002.jpg)
FIGURE 2. Differentially expressed profiles and functional analysis of mRNAs in a force-related way. Heatmaps display the top 50 genes with the highest values of |log2foldChange (FC)| in the intermittent force group (A), static force group (B), and different force type group (C). Intersection of results gained from two force groups contributes to the list of force-sensitive mRNA (FS mRNA), which is visualized by the VENNY diagram (D). Top 10 GO annotations grouped by BP, CC, and MF of FS mRNAs (E) and FTS mRNAs (F). GO terms are ranked by the GeneRatio scores. The size and color of the dots represent the gene count and q-value score of the enrichment analysis, respectively.
GO annotation and enrichment analysis of the DE mRNAs in different groups provided us with a better understanding of the physiological process of PDLCs responding to mechanical force. Under the cutoff limitation of p < 0.05, 396 GO terms were screened out, with 359 in the biological process (BP), 21 in the cellular component (CC), and 13 in the molecular function (MF) based on the FS mRNAs. As for the FTS mRNA group, total terms of GO are counted as 582 with 535, 32, 15 for BP, CC, and MF, respectively. BP terms enriched in FS mRNAs are most related to development processes like cellular developments, bone mineralization, and ossification, verifying that the cellular differentiation processes were impacted by applying force. As for CC terms, extracellular matrix, contractile fibers, and communication elements are most enriched among FS mRNAs, suggesting that the force would significantly induce the reaction modes of cells. The mechanically responsive cellular activity could also be indicated by the MF terms enriched in stimulus reception and communication such as receptor ligand activity. Interestingly, when it comes to FTS mRNAs, the GO terms are largely connected with nuclear and chromosome activity, implying the possibility for transcriptional differences when reacting to different force types (Figures 2E,F).
Validation of CFS miRNAs
In total, research carried out in MSCs, bone marrow stem cells (BMSCs), adipose-derived stem cells (ADSCs), PDLSCs, PDLCs, tendon-derived stem cells (TDSCs), alveolar bone cells (ABCs), primary rat osteoblasts (prOB), osteoblasts, MC3T3-E1 cells, MLO-Y4 osteocytes, and C2C12 cells was included and is listed in Table 2. Apart from hsa-miRNA, miRNAs that are conservative among species by sequence matching were also included. Finally, 38 CFS miRNAs, with 22 negatively and 16 positively related to osteogenic differentiation, were collected for further analysis. Based on CFS miRNAs, KEGG enrichment was carried out, and 59 pathways were screened out by mirPath. To focus on force/osteogenesis-related processes, we ranked pathways according to p-value and visualized the top 22 pathways after filtering those relating to cancer and chemical synthesis (Figure 3). Notably, TGF-β, Hippo, PI3K-Akt, Wnt, FoxO, and Ras signaling pathways are all recognized FS pathways, corroborating the credibility of the bioinformatics strategy. Furthermore, KEGG terms such as ECM–receptor interaction, adherens junction, regulation of actin cytoskeleton, and focal adhesion also ranked among the top 20, suggesting the involvement of ncRNA in controlling cell morphology to mechanical input. It is explainable for the cancerous pathway like proteoglycans in cancer, pathways in cancer, and melanoma (nos. 5, 6, and 21, respectively) to be enriched in our analysis, considering cancer represents the most common model in miRNA studies. Overall, KEGG results of CFS miRNAs highly conform to the established data, ensuring the confidence of the following analysis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g003.jpg)
FIGURE 3. KEGG enrichment of CFS miRNAs corroborating the reliability of bioinformatics. Pathways are ranked according to p-value, and the top 22 are displayed after filtering those relating to cancer and chemical synthesis. KEGG terms included are highly consistent with verified FR pathways at present. Pathways excluded are proteoglycans in cancer, pathways in cancer, fatty acid biosynthesis, mucin-type O-glycan biosynthesis, morphine addiction, and melanoma (nos. 5, 6, 8, 11, 18, and 21, respectively).
FR CeRNA Network Visualization and FR LncRNA Exploration
One hundred seventy-seven lncRNAs and 442 pairs of lncRNA–miRNA, as well as 1,445 target mRNAs and 2,463 miRNA–mRNA pairs, were predicted based on CFS miRNAs. The intersection of predicted mRNA targets and FS mRNAs was carried out, and 47 FS mRNAs and 79 mRNA-miRNA pairs were then included. Based on the 24 miRNAs identified by FS mRNA–CFS miRNA pairs, 147 lncRNAs and 408 CFS miRNA–lncRNA pairs were chosen for FR ceRNA network construction. Last, lncRNAs with a degree lower than three were excluded based on the network analysis carried out by Cytoscape. After filtration, 122 nodes as well as 340 edges were finally reserved, which include 51 lncRNAs, 24 miRNAs, and 47 mRNAs (Figure 4A). Strikingly, the network analysis displayed by Cytoscape indicated that all CFS miRNAs show a value of degree ≥5, and 14 out of the 24 CFS miRNAs have their degrees higher than 10, highlighting the core role as well as the credibility of our CFS miRNAs. Notably, the intricate and entwined edges around miR-20a, miR-195, and miR-424, the top three miRNAs with the highest degrees of 20, 29, and 31, respectively, suggest possible dominating roles of these miRNAs in mechanobiology, which is worthy of future exploration. As for lncRNA prediction, the exclusion of lncRNAs with a degree lower than three contributed to a more confidential outcome as 51 lncRNAs were finally included, named the FR lncRNAs. Interestingly, when ranked by degree, betweenness, or closeness, 13 lncRNAs, namely, NEAT1, XIST, OIP5-AS1, MALAT1, AC021078.1, TUG1, SNHG16, AC005261.1, AC021092.1, FGD5-AS1, H19, SNHG1, and AL035425.3, are always listed among the top 15 regardless of the indicator of ranking, highlighting their potential key roles in the mechanoresponsive process. For in vitro validation of our prediction, expression changes of 10 FR lncRNAs were checked by qPCR (Figure 4B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g004.jpg)
FIGURE 4. Identification of FR lncRNA through the ceRNA network and RT-qPCR validation. In the ceRNA network, rhombuses, squares, and circles stand for lncRNAs, miRNAs, and mRNAs, respectively. Upregulated RNAs are shown in pink, while downregulated ones are shown in blue. Gray lines represent regulatory relationships (A). The validation of lncRNA expression changes via RT-qPCR under mechanical stretch was normalized to Gapdh (n = 3) (B).
FR LncRNA Exploration Based on FS mRNA Interaction and Co-expression Network
Grounded on the list of FR lncRNAs gained from the strong evidence provided by FS mRNA–CFS miRNAs, the next goal is to explore the functionality of FR lncRNAs. To start with, lncRNA, as an intricate type, could not only impact mRNA indirectly though miRNA as a ceRNA mechanism, it could directly interact with mRNAs as well. Sixty-eight FS mRNAs and nine FR lncRNAs were verified to have direct interaction relationships. XIST, MALAT1, NORAD, and SNHG1/12/16 showed highly intertwined interaction with FS mRNAs (Figures 5A,B). Second, based on the guilt-by-association method, it is supposed that highly co-expressed genes during the biological process hold great possibility to share a similar biological function. Thus, novel lncRNAs’ function could largely be predicted by co-expressed protein-coding mRNAs (Liao et al., 2011). In our work, 910 pairs of FS mRNA–FR lncRNA showed Pearson’s correlation coefficients ≥0.8 (p < 0.05), which further verified the credibility of our prediction. For a BP-based explanation, we screened out several force-responsive pathways which were highly enriched among these pairs (p < 0.05), such as the PI3K-Akt signaling pathway, focal adhesion, regulation of actin cytoskeleton, calcium signaling pathway, and Wnt and Hippo signaling pathways. Connecting the co-expression network with a pathway could further suggest the potential roles of FR lncRNAs (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g005.jpg)
FIGURE 5. FS mRNA and FR lncRNA interaction. The direct interaction relationships were visualized (A) and the Ce-loop RNA network was also extracted (B). Rhombuses, squares, and circles stand for lncRNAs, miRNAs, and mRNAs, respectively. Upregulated RNAs are shown in pink, while downregulated ones are shown in blue. Gray lines represent regulatory relationships.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g006.jpg)
FIGURE 6. The co-expression network among FS mRNA and FR lncRNA. Rhombuses and squares stand for lncRNAs and mRNAs, respectively. Upregulated RNAs are shown in pink, while downregulated ones are shown in blue. Red and blue lines represent positive and negative correlations, respectively. The widths of lines are mapped to correlation scores.
Regulation of FR LncRNAs Realized by TFs and the Possible Bidirectional Interaction Between TF and LncRNAs
Based on the prediction through PWMs, 24 out of 51 FR lncRNAs were screened out, and 259 corresponding TFs were enrolled in the network, suggesting another potential mechanism of force-sensitive response. TFs with a degree lower than five were excluded for demonstration, and 121 TFs were finally reserved. RHOXF1, MZF1, ZNF354C, NFIX, KLF5, MEIS1, NFIC, HIC2, SP1, FOXD2, and FOXP3 are the top 11 TFs with a degree higher than 10, highlighting their key roles of potential regulation on lncRNAs. What bears our attention is that a large proportion of these TFs were explored to function in bone-related biological processes as well as mechanosensitive regulation, which could serve as validation for the plausibility of our network from one perspective. In turn, our work provided a systematic view and suggested the potential intricate pathways these TFs may involve in (Figure 7A) Interestingly, the shuttling of TFs from the cytoplasm to nuclear area provides the possibility for their bidirectional interaction with lncRNAs. From one perspective, as the mRNA target of miRNA sponged by lncRNAs in the cytoplasm, TFs are exposed to the epigenetic modification by lncRNAs. However, after TFs translocated to the nuclear area, it is the TFs which bind to the upstream transcripts of lncRNAs that master the transcription regulation. Hence, we explored the possibility of a triplet circular network in both the FS mRNA group and FTS group. Fortunately, we found NFIA-(MIR4458HG/H19/MIR497HG)-miR-29a-NFIA based on the FS group (Figure 7B), and the regulation effect between miR-29 and NFIA in bone remodeling was already confirmed in a previous study (Franceschetti et al., 2013). Meanwhile, the regulation loops were also identified in the FTS group as FOS-(ARHGAP27P1-BPTFP1-KPNA2P3)-miR-221-FOS, FOS-(MIR4458HG/XIST/OIP5-AS1)-miR-29a/b/c-FOS, and E2F1-(FGD5-AS1)-miR-20a-E2F1 (Figure 7C). In conclusion, NFIA, FOS, and E2F1 showed significant alternation of expression profiles when exposed to the force-related environment while embodying the potential to form the force-sensitive circuit with lncRNAs and miRNAs, which could be studied in the future.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g007.jpg)
FIGURE 7. The FR TF–lncRNA regulatory network. The TF–FR lncRNA network is visualized with lncRNAs forming the middle circle while TFs with degree < 10 standing on the left and those with degree ≥ 10 on the right. Node colors are mapped to degree scores. (A) The feedback loops based on the triplet of TF (FS mRNA)–CFS miRNA–FR lncRNA (B) as well as the triplet of TF (FTS mRNA)–CFS miRNA–FR lncRNA (C) are displayed.
Elucidation of the Downstream Pathway Map of FR lncRNAs
Literature retrieval provided evidences of lncRNAs’ effects on mediating the cellular response to mechanical cues. Here, we mapped out a systematic evaluation of the pathways that our FR lncRNAs may involve in (Table 3). NEAT1 and H19, with high scores in our network, have gained wide recognition in mechanobiology by previous study, which confirmed the reliability of our work. Remarkably, studies of NEAT1 or H19 were carried out under different force conditions in various cells, which indicated that these lncRNAs may be universal mechanoresponsive factors like our study aimed to define (Hirose et al., 2014; Heo et al., 2015; Le et al., 2016; Li et al., 2017, 2018; Izadpanahi et al., 2018; Wu et al., 2018; An et al., 2019; Wang S. et al., 2019; Zhang Y. et al., 2019; Huang et al., 2019; Todorovski et al., 2020).
Though research of lncRNAs’ involvement in mechanotransduction of the bone remains relatively void, studies of other cellular systems exploring essential pathways of force sensing could also fuel our understanding. As aforementioned, once the mechanical cues are transformed as intracellular signals, several intersected pathways are then involved to conduct the following cytoplasm-to-nuclear information delivery. RhoA/ROCK, Hippo/LATS1/2 kinase cascade, YAP/TAZ, AKT, JUN/p38 MAPK, BMP/Smad1/5/8, TGF-β/Smad2/3, and canonical and noncanonical Wnt pathways are all crucial factors involved in mechanobiology (Azzolin et al., 2012; Teramura et al., 2012; Furumatsu et al., 2013; Ji et al., 2014; Zhang et al., 2015; Gao et al., 2016; Panciera et al., 2017; Stanley et al., 2019). Therefore, lncRNAs predicted in our network, connecting these crucial factors and pathways, hold the promise to be key players in the process and are listed in Table 3, which could support our prediction as well boost downstream research.
Nevertheless, the project of lncRNA in the mechanical field is far from complete. For a deeper understanding, we also explored the possible mechanoresponsive and osteogenic pathways that these lncRNAs may participate in. H19, KCNQ1OT1, TUG1, SNHG1, and MALAT1, highlighted in our network, are closely correlated with osteogenesis via the ceRNA mechanism and beyond, which mainly involves BMP/Smad, RUNX2, MAPK, and Wnt/β-catenin signaling pathways (Yu C. et al., 2018; Gong et al., 2018; He et al., 2018; Gu H. et al., 2019; Liu H. et al., 2019; Wang C.-G. et al., 2019; Wu et al., 2019a; Jiang et al., 2019; Wu D. et al., 2020; Bi et al., 2020; Wang J. L et al., 2020; Xiang et al., 2020; Xiaoling et al., 2020; Yu et al., 2020), and detailed relationships are concluded in Table 4.
Pathway–LncRNA Network Building Combining Predicted and Established Data
Genome-wide sequencing and bioinformatics analyses, as potent tools in ncRNA research, facilitate an overall understanding of lncRNA and were utilized by researchers to target ideal genes and molecules. To begin with, we conclude sequencing outcomes and computational strategies concerning FS lncRNAs published in recent years in Table 4, which could, on one hand, provide hints to future researches and, on the other hand, attenuate the improving quality and popularity of bioinformatics analysis.
Based on these published articles, we integrated the predicted pathways detailed in the last section with the proven ones, to illustrate the intricate pathway network of osteogenesis provoked by mechanical input modulated by lncRNAs (Figure 8). An overview of the pathways could stand as a valid proof of the conformity between predicted data with confirmed ones. More importantly, a pathway-centered network could provide insightful views of the interaction between lncRNAs. For instance, the Wnt/β-catenin signaling pathway displays a connection with lncRNAs as H19, SNHG1, SNHG16, TUG1, KCNQ1OT1, and LINC-00662, suggesting the possible interactive relationships of ncRNA modulation. The same goes for pathways like YAP/TAZ, BMP, Smad, and TGF-β. Though the function mode of ncRNAs as working via the network has long been recognized, previous studies mostly concentrated on the connections predicted by molecular bonding. Our work introduced a novel strategy to explore the interaction though pathway connection, which could be more specific when applying to a set of biological process like mechanical stimulation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g008.jpg)
FIGURE 8. Intricate pathway network provoked by mechanical input is modulated by lncRNAs. Blue ellipses and green squares stand for FR lncRNAs and verified pathways, respectively.
Moreover, a combinative strategy based on our prediction and verified background would also provide new insights into a particular molecule. For instance, considering the centric role of H19 in force response, revealed by both our network analysis and established data, we integrated and visualized the ceRNA modulation and pathways H19 stimulates by confronting the mechanical input (Figure 9). Certificated roles of H19 could echo our prediction, while its complex interaction with CFS–miRNA and FS mRNA would suggest a promising field of research.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/780211/fbioe-10-780211-HTML/image_m/fbioe-10-780211-g009.jpg)
FIGURE 9. ceRNA network and downstream pathways revealed H19 as a key FS lncRNA confronting the mechanical input. Orange squares and dark yellow ellipses in the left semicircle stands for FS mRNA and CFS miRNA related to H19, respectively. White ellipses in the right semicircle represent the verified pathway H19 involved in FS condition and osteogenic processes.
Discussion
Mechanical force, which used to be so prevailing that it was often taken for granted, is now gaining the spotlight as a crucial factor for homeostasis maintenance (Wang et al., 2018). Particularly, as the supporting organ receiving and confronting complex force stimuli, the bone represents one of the most active tissue undergoing frequent remodeling, responding to mechanic cues (Yong et al., 2020). Outstandingly, though ncRNAs largely outnumber coding ones and master various properties, lncRNAs, as one type of ncRNA verified to have a diversified function mode and conformation, are less explored for their function in a mechanosensitive way. Aiming at providing a glimpse into universal lncRNA responding to mechanical stress, we first constructed a ceRNA network based on FS mRNAs and CFS miRNAs, which were individually screened out from the intersection of DE mRNA under different mechanical environments and from paper-based experimental data of high confidence. Consequently, lncRNAs involved in our ceRNA network with high degree were identified as FR lncRNAs. We carried out the functional analysis of FR lncRNAs from three aspects, namely, the FS mRNA interaction, TF regulation, and pathway network. To start, the direct interactions among FS mRNA and FR lncRNAs were displayed as well as the co-expression networks. Second, we explored the potential upstream regulation mediated by TFs on lncRNAs and put forward the underlying crosstalk loop between them, demonstrating the possible regulatory transcriptional activities. As for the downstream exploration, pathways concerning osteogenesis and mechanotransduction were mainly focused for the overview of FR lncRNAs. Later on, we developed combinative strategy to illustrate our predicted data against established ones of FS lncRNA, constructing an lncRNA-pathway network and identifying H19 as a critical factor in mechanobiology with more potential roles. In conclusion, our work contributed to a systematic prediction of a quartet network with lncRNAs at its core, elucidating the list of potential FR lncRNAs as well as their functional modes based on creditable evidence, suggesting a promising field worthy of further exploration and a new mode of bioinformatic analysis combining sequencing data, paper-based evidence, and computational prediction.
As the core component of our analysis, miRNAs and lncRNAs are no more than the so-called “noise” of transcriptome. Speculation as well as evidences are growing that ncRNAs may function in a more extensive way in the bone and beyond (Chen et al., 2017; Fang et al., 2019; Zhu et al., 2019). The past few years have seen the topic of miRNA’s role in the mechanical condition rush into the forefront (Chen et al., 2017; Wang et al., 2018; Fang et al., 2019), which was also concluded in our work. The FR ceRNA network demonstrated that miR-195, miR-424, and miR-20a are the top three miRNAs with the highest degree. Coherently, presented studies confirmed that miR-195 plays as an inhibitor of PDLC osteogenesis in both mice and human, targeting essential bone formation genes such as WNT family member 3A (WNT3A), fibroblast growth factor 2 (FGF2), and bone morphogenetic protein receptor-1A (BMPR1A) (Chang et al., 2017). Meanwhile, miR-424 was identified as one of the core miRNAs of tension-force-induced bone formation of PDLC in a former study (Chang et al., 2015). Additionally, in vivo study based on loaded and unloaded tibia of mice indicated that miR-20a was the only member of the miR17/92 cluster, which displayed significant changes responding to force (Mohan et al., 2015). Alongside these reported functions, our network indicated a more entwined relationship of these three miRNAs with mechanotransduction, which is open for further exploration. What also merits our attention are both the consistency and the discrepancy displayed by miRNAs under different stimuli or in diversified cellular environments. For instance, miR-33-5p promotes osteoblast differentiation regardless of the types of mechanical cues as MG or fluid shear stress (Wang et al., 2016). MiR-494 is a negative regulator of osteogenesis under both compressive and unloading conditions in MC3T3-E1 cells, C2C12 cells, and osteoblasts (Iwawaki et al., 2015; Qin et al., 2019). In hPDLCs, miR-29b shows altered expression level to both cyclic stretch and compression forces (Chen et al., 2015). However, it is only responsive to the mechanical tensile strain in osteocytes, but not osteoblasts of mouse (Zeng et al., 2019). MiR-138 was proven to damper the osteogenesis of MSCs regardless of their origins as bone, adipose, or tendon under the extracorporeal shockwave or cyclic strain (Hu et al., 2016; Wu et al., 2018). However, Yun Wu et al. (2019b) found miR-138-5p to be upregulated in qPCR and high sequencing under the cyclic strain of hPDLCs, which lacks reasonable explanation for the difference. Generally, though deeper exploration and validation should be carried out for the specific variance, the mechanosensitive properties of miRNAs remain coherent among cells and forces, supporting our systematic evaluation of the force-sensitive network.
Notably, another crucial prediction network in our study is carried out to evaluate the TF–lncRNA pairs, which was highlighted as an essential regulator under mechanical cues and beyond (Pylayeva-Gupta, 2011; Alam et al., 2014; Ji et al., 2020). Strikingly, TFs predicted to be the regulators of FR lncRNAs in our network are grounded for their roles in mechanotransduction and bone formation. Among the top 11 TFs with the highest degree, KLF5 mediates the response of hPDLCs to cyclic tensile stress and contributes to the progenitor maintenance in bone marrow via Rab5-β1/β2 integrin trafficking (Taniguchi Ishikawa et al., 2013; Gao J. Y. et al., 2018). Sensing the tensile forces, SP1 would enhance its interaction with p38 with elevated phosphorylation and activate its regulation on filamin-A (D’Addario et al., 2002). One orthodontic tooth movement model also proposed a mechanotransduction mode involving FAK, Ras-GTPase, and ERK1/2 regulating ephrin B2 via SP1 in periodontal ligament fibroblasts (PDLFs) (Diercke et al., 2011). The LINC00511/miR-150-5p/SP1 feedback loop in osteoarthritis (OA) was proven to master the extracellular matrix synthesis of chondrocyte (Zhang L. et al., 2020). Meanwhile, SP1 and MZF1 contribute to, whereas SP3 antagonizes, N-cadherin promoter activity in osteoblasts (Le Mée et al., 2005). MZF1, on the other hand, would regulate the bone formation of hMSCs via the osteopontin–MZF1–TGF-β1 linked pathway (Driver et al., 2015). NFIC was valued as a transcriptional switch of BMSCs to osteogenic or adipogenic differentiation for its interaction with diversified factors such as TGFβ1, KLF4, Wnt/β-catenin, and Runx2 (Lee et al., 2014; Roh and Park, 2017; Zhou et al., 2017; Huang et al., 2020). NFIX, regulated by miR-25-3p, also turned out to be a positive regulator of osteoclast proliferation and bone ossification (Driller et al., 2007). Elucidating TFs and their interactive relationship with FR lncRNAs, on one hand, contributes to a better understanding of the regulation network in mechanical transduction and, on the other hand, corroborates our prediction work as a qualified one worthy of further study.
As the main part of our studies, FR lncRNAs also show coincidence with a previous study. Remarkably, NEAT1, as the FR lncRNA with the highest score in our network, may exert an extensive epigenetic role in mechanotransduction. Elevated expression of NEAT1 in PDLSCs confronting compressive force (CF) was verified (Huang et al., 2019), and it was proven to connect the mechanosensitive epigenetic silencing via PRC2-H3K27me3 (Hirose et al., 2014; Le et al., 2016; An et al., 2019; Todorovski et al., 2020). In the myogenic process of C2C12 cells, this role was suggested to depend on the interaction between NEAT1 and EZH2 (Wang S. et al., 2019). Meanwhile, in PDLSCs, EZH2 is responsible for the global change of H3K27me3 and H3K27ac activity to encode the mechanical information onto chromatin when CF was applied both in vivo and in vitro (Heo et al., 2015). Thus, NEAT1 holds the promise of being one epigenetic modification strategy utilized by cells to memorize the mechanical cues after the withdrawal of physical stimuli. Considering its verified osteogenic promotion role in hBMSCs (Zhang Y. et al., 2019), a deeper investigation of NEAT1 functionality in connecting force and bone remodeling is worth exploring. H19, another lncRNA involved in our network, can also be a case in point. Under mechanical tension, H19 was elevated as a miRNA sponge of miR-138 to release PTK2, which codes FAK, inducing osteogenesis of hBMMSC via the FAK-ERK1/2-Runx2 pathway (Wu et al., 2018). H19 is also involved in controlling MSCs’ fate, confronting nanotopographical cues of niches as the alignment of fibers via the BMP signaling pathway (Izadpanahi et al., 2018). During mechanical unloading, alternation of H19 may be credited to elevated DNMT1 expression. Three to four weeks of stimulated unloading resulted in an abnormally increased profile and higher nucleic-to-cytoplastic location ratio of DNMT1, causing 5 mC enrichment at CpGs in the promoter region of H19. Hypermethylation of the promoter results in lower levels of H19, followed by inhibited phosphorylation and activation of ERK and impact on TGF-β, WNT, and JAK-STAT pathways, resulting in microgravity (MG)-induced bone loss (Li et al., 2018). Another study confirmed that H19 may serve as the initiation of DOP via Dkk4, controlling the GSK phosphorylation ratio, which impacts the Wnt signaling pathway (Li et al., 2017). Our work, promisingly, could suggest possible network they function by and boost further studies. Besides these two well-studied lncRNAs, various creditable FR pathways were also highlighted for lncRNA identification. For instance, the RhoA/ROCK pathway is activated as a linker for the extracellular matrix to affect cellular processes, concerning the rearrangement of cytoskeleton fibers and following cell morphology changes. As an upper-stream regulator, RhoA would facilitate Smad1/5/8 phosphorylation (Ji et al., 2014) and activate YAP/TAZ (Stanley et al., 2019) while inhibiting AKT phosphorylation (Teramura et al., 2012). Significantly, YAP/TAZ stands out as the master integrator conducting the mechanic information. YAP/TAZ would shuttle from the cytoplasm to the nucleus for its binding with the transcriptional enhanced associate domain (TEAD) once sensing the signals from FA, actomyosin cytoskeleton, membrane channels, and the classic Hippo/LATS1/2 kinase cascade (Panciera et al., 2017). Furthermore, pathways such as the JUN/p38 MAPK, BMP/Smad1/5/8, TGF-β/Smad2/3, and canonical and noncanonical Wnt pathways, which are main factors involved in osteogenesis, embody the intersection with YAP/TAZ for mechanical force conduction (Azzolin et al., 2012; Gao et al., 2016). Taken together, the pathways themselves contribute to a network. Thus, we dug into these pathways to find the evidence concerning the roles of lncRNAs in Table 3. Specifically, lncRNAs like MALAT1 and NORAD, are found to have a perplexing interaction with YAP in diversified processes (Ventura, 2016; Gu L. et al., 2019; Jin et al., 2019; Tan et al., 2019; Yao et al., 2019). Besides, NORAD is highly involved with the TGF-β/Smad pathway direct mRNA or miRNA targeting (Kawasaki et al., 2018; Yang et al., 2019; Zhou et al., 2019). LINC00662, SNHG16, and SNHG1 are entwined with the Wnt pathway (Christensen et al., 2016; Wu W. et al., 2020; Tian et al., 2020). SNHG1, binding to EZH2, would also impact the PI3K/AKT/mTOR (Chen et al., 2020) as well as ROCK1 pathway via miR-101-3p (Deng et al., 2019). Though these models represent a close interaction of lncRNAs and key pathways in cancerous background, which is the most generalized model for lncRNA research, it could also stand as a testimony for their promising relationship in FR environments.
Clinically, essential roles of FS lncRNAs in skeletal mechanobiology as well as their therapeutic application potential could not be overlooked. Generally, MG, tensile, and compressive forces are most utilized as mechanical stimuli, mimicking varied physiological and pathological conditions. Firstly, MG implies, for the study of disuse osteoporosis (DOP), a common bone disorder when patients are in a long-term-immobilization condition like fracture, therapeutic bed rest, and space flight. H19 was found to be downregulated in the MG condition and responsible for disuse-induced bone loss (Li et al., 2017; Li et al., 2018). Its overexpression turns out to be a remedy for postmenopausal osteoporosis considering its potent effect on miR-19b-3p inhibition (Xiaoling et al., 2020). Alongside that, in a bone-related disease such as delayed fracture healing, which is characterized by the impaired ability of bone regeneration, reduced expression of lncRNA KCNQ1OT1 was suggested to be a potential biomarker while the overexpression of this lncRNA would be a possible therapeutic strategy (Gao X. et al., 2018). Novel lncRNAs such as ODSM and OGRU are also recognized as MG-sensitive ones, whose promising therapeutic roles in osteoporosis were confirmed both in vitro and in vivo (Wang M et al., 2020; Wang Y et al., 2020). While mechanical disuse would impair bone quality, muscle atrophy also represents a concern under MG circumstances. Besides, muscle–bone crosstalk against disuse condition was verified, and muscle loss would also cause bone loss (Bettis et al., 2018). Inspiringly, enforced expression of skeletal muscle-specific and MG-sensitive lncRNA, such as lncMUMA, would not only prevent muscle atrophic development but also ameliorate and even reverse the established atrophy both in vivo and in vitro (Zhang et al., 2018). Other than mechanical unloading, compressive and tensile forces contribute to bone remodeling during orthodontic tooth movement (OTM) progress, where bone resorption and formation are kept in a delicate balance. Axiomatically, minor changes induced by ncRNAs would greatly change the inclination, indicating the possibility of ncRNA modification for complication avoidance and improving efficacy of therapy. DANCR and H19 were both proved to favor osteogenesis when cells were exposed to strain force (Wu et al., 2018; Zhang X. et al., 2019). Depletion of DANCR could hinder osteoclast formation and ameliorate root resorption induced by CF both in vitro and in vivo (Zhang X. et al., 2019). Upregulated exosomal H19 of umbilical MSCs (U-MSCs), induced by an appropriate rotary cell culture condition mimicking a tensile force, exert significant remedial effects in promoting chondral regeneration and pain relieving (Yan et al., 2021). However, though adequate loading serves as a positive factor for bone homeostasis, overuse or hypertension would cause undesirable outcomes like osteoarthritis (OA). LncRNA-MSR and lncRNA-CIR, activated by excessive mechanical stress, would mediate consequent cytoskeleton disorganization and ECM degradation (Liu et al., 2016). Nonetheless, as implant and bioengineering gain popularity, the topography of the matrix and physical properties of the cell culturing environment will begin to attract researchers’ sights. MEG3 and H19 are involved in controlling MSCs’ fate in confronting nanotopographical cues of niches as a way of fiber alignment (Izadpanahi et al., 2018). Upregulated lncRNA PWRN1-209 in hBMSCs exposed to a Ti surface with microtopography would promote osteogenesis via integrin-FAK-ALP signaling (Wang M et al., 2020). Extrinsic force, like MG, CF, and strain, and intrinsic force like topography and matrix stiffness serve as essential factors for bone maintenance, and illustrating ncRNAs’ role in modulation would favor the development of tissue engineering and clinical therapy. More importantly, network analysis could greatly promote the study of skeletal diseases, considering their dynamic properties and intricate pathogenesis.
Bioinformatics has gained popularity in recent years. Though in silico analysis could pale in comparison with biological experiments, its strength in data processing and network analysis is impeccable. Meanwhile, growing evidences showed that predicted data of lncRNAs concerning mechanobiology could be observed in biological systems. For instance, based on the expression profile of mouse cementoblasts under static compressive force gained from sequencing, H. Liu et al. (2019) selected lncRNAs such as Prkcz2, Hklos, Trp53cor1, Gdap10, and Ak312-ps for validation, and results of RT-qPCR highly match those of sequencing data. The research group also applied compressive force to hPDLSCs (Huang et al., 2019), where RT-qPCR also validated the microarray data of FS lncRNAs such as FER1L4, HIF1A-AS2, MIAT, NEAT1, ADAMTS9-AS2, LUCAT1, MIR31HG, and DHFRP1. Moreover, upregulated FER1L4 was later proven to facilitate the compression-induced autophagy via AKT/FOXO3 as predicted (Huang et al., 2021). However, previous studies concerning lncRNAs’ role in a mechanosensitive way mainly depended on the expression profile by microarray or high-throughput sequencing (Zhang et al., 2017; Lv et al., 2018; Huang et al., 2019; Wang H et al., 2020). A limited scope and confined force type should be considered among these studies. While creating a general work of force condition to provide an overview of FR lncRNAs, we also adopted diversified strategies to guarantee and emphasize the reliability of our bioinformatic research. First, we developed a strict inclusion criteria of data and adopted prediction strategies of high stringency. The intersection of DE mRNAs under different mechanical forces contributed to the reliable identification of the FS mRNAs, and paper-based inclusion of CFS miRNAs ruled out false-positive results which usually appeared in miRNA sequencing (Iwawaki et al., 2015; Qin et al., 2019). Moreover, the predicted pairs of mRNA–miRNA were verified in all three datasets based on different computational algorithms, while lncRNA was only included when it met the standard of acquiring a highly stringent pair relationship with CFS miRNA and a degree of no fewer than three in the ceRNA network. Last, systematic analysis of lncRNAs, including the upstream factors such as TFs, the downstream pathway-mapping strategy, and the combinative analysis against FS lncRNAs, would in turn corroborate our results. Axiomatically, drawbacks exist in our work. Despite all the efforts to ensure confidence level, the flaw of bioinformatic analysis is that uncertainty persists. Besides, when our scope was expanded to all mechanical force types and included diversified bone-related cells, mechanisms specialized by cells or mechanical cues were overlooked. Certainly, though the intricate network could provide us with a blueprint, the whole network may not be activated in all mechanical situation, and factors involved under specific conditions should be further confirmed. Undeniably, concerted efforts with studies in the future would be helpful to provide further validation to our work.
Unequivocally, identification of lncRNAs needs a case-by-case basis, yet a systematic illustration is just as important, if not more so, as to provide a stepping stone to uncovering the panoramic view of this specific type of ncRNA. Though a question mark still hangs over the world of lncRNAs against mechanic cues, our work, hopefully, would produce a moment of clarity to break this impasse.
Data Availability Statement
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding authors.
Author Contributions
Conceptualization, JC and DJ; Data acquisition, JC, CL, and SL; Data analysis, JC, CL, SL, and KY; Writing—Original Draft, JC, KY, XXF, JW, and PY; Visualization, JC, KY, and XG; Writing—Review and Editing, JY, JW, DJ, and ZZ; Funding Acquisition and Supervision, DJ and ZZ.
Funding
This study was supported in part by grants from the National Natural Science Foundation of China (Nos. 81901041, 81901035, 81771048, 81801018, and 82071127), Sichuan Science and Technology Program (2020YFS0170 and 2019YJ0097), China Postdoctoral Science Foundation funding (No. 2020M673264), and Research and Development Foundation of West China Hospital of Stomatology, Sichuan University (No. RD-02-201911).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
CTS, cyclic tensile strain; CFS miRNA, creditable force-sensitive miRNA; FS mRNA, force-sensitive mRNA; FTS mRNAs, force type sensitive mRNAs; FR lncRNAs, force-related lncRNAs; FR ceRNA, force-related ceRNA; FS lncRNA, force-sensitive lncRNA; TFBSs, transcriptional factor binding sites.
References
Alam, T., Medvedeva, Y. A., Jia, H., Brown, J. B., Lipovich, L., and Bajic, V. B. (2014). Promoter Analysis Reveals Globally Differential Regulation of Human Long Non-coding RNA and Protein-Coding Genes. PLoS One 9, e109443. doi:10.1371/journal.pone.0109443
An, H., Tan, J. T., and Shelkovnikova, T. A. (2019). Stress Granules Regulate Stress-Induced Paraspeckle Assembly. J. Cel Biol. 218, 4127–4140. doi:10.1083/jcb.201904098
Azzolin, L., Zanconato, F., Bresolin, S., Forcato, M., Basso, G., Bicciato, S., et al. (2012). Role of TAZ as Mediator of Wnt Signaling. Cell 151, 1443–1456. doi:10.1016/j.cell.2012.11.027
Bettis, T., Kim, B.-J., and Hamrick, M. W. (2018). Impact of Muscle Atrophy on Bone Metabolism and Bone Strength: Implications for Muscle-Bone Crosstalk with Aging and Disuse. Osteoporos. Int. 29, 1713–1720. doi:10.1007/s00198-018-4570-1
Bi, H., Wang, D., Liu, X., Wang, G., and Wu, X. (2020). Long Non-coding RNA H19 Promotes Osteogenic Differentiation of Human Bone Marrow-Derived Mesenchymal Stem Cells by Regulating microRNA-140-5p/SATB2 axis. J. Biosci. 45. doi:10.1007/s12038-020-0024-y
Bouzid, T., Kim, E., Riehl, B. D., Esfahani, A. M., Rosenbohm, J., Yang, R., et al. (2019). The LINC Complex, Mechanotransduction, and Mesenchymal Stem Cell Function and Fate. J. Biol. Eng. 13. doi:10.1186/s13036-019-0197-9
Cai, X., Liu, Y., Yang, W., Xia, Y., Yang, C., Yang, S., et al. (2016). Long Noncoding RNA MALAT1 as a Potential Therapeutic Target in Osteosarcoma. J. Orthop. Res. 34, 932–941. doi:10.1002/jor.23105
Cech, T. R., and Steitz, J. A. (2014). The Noncoding RNA Revolution-Trashing Old Rules to Forge New Ones. Cell 157, 77–94. doi:10.1016/j.cell.2014.03.008
Chang, M., Lin, H., Fu, H., Wang, B., Han, G., and Fan, M. (2017). MicroRNA-195-5p Regulates Osteogenic Differentiation of Periodontal Ligament Cells under Mechanical Loading. J. Cel. Physiol. 232, 3762–3774. doi:10.1002/jcp.25856
Chang, M., Lin, H., Luo, M., Wang, J., and Han, G. (2015). Integrated miRNA and mRNA Expression Profiling of Tension Force-Induced Bone Formation in Periodontal Ligament Cells. In Vitro Cell.Dev.Biol.-Animal 51, 797–807. doi:10.1007/s11626-015-9892-0
Chen, J., Wang, F., Xu, H., Xu, L., Chen, D., Wang, J., et al. (2020). Long Non-coding RNA SNHG1 Regulates the Wnt/β-Catenin and PI3K/AKT/mTOR Signaling Pathways via EZH2 to Affect the Proliferation, Apoptosis, and Autophagy of Prostate Cancer Cell. Front. Oncol. 10, 552907. doi:10.3389/fonc.2020.552907
Chen, N., Sui, B. D., Hu, C. H., Cao, J., Zheng, C. X., Hou, R., et al. (2016). MicroRNA-21 Contributes to Orthodontic Tooth Movement. J. Dent. Res. 95, 1425–1433. doi:10.1177/0022034516657043
Chen, Y., Mohammed, A., Oubaidin, M., Evans, C. A., Zhou, X., Luan, X., et al. (2015). Cyclic Stretch and Compression Forces Alter microRNA-29 Expression of Human Periodontal Ligament Cells. Gene 566, 13–17. doi:10.1016/j.gene.2015.03.055
Chen, Y., and Wang, X. (2020). miRDB: an Online Database for Prediction of Functional microRNA Targets. Nucleic Acids Res. 48, D127–D131. doi:10.1093/nar/gkz757
Chen, Z., Zhang, Y., Liang, C., Chen, L., Zhang, G., and Qian, A. (2017). Mechanosensitive miRNAs and Bone Formation. Ijms 18, 1–16. doi:10.3390/ijms18081684
Chou, C.-H., Shrestha, S., Yang, C.-D., Chang, N.-W., Lin, Y.-L., Liao, K.-W., et al. (2018). miRTarBase Update 2018: a Resource for Experimentally Validated microRNA-Target Interactions. Nucleic Acids Res. 46, D296–D302. doi:10.1093/nar/gkx1067
Christensen, L. L., True, K., Hamilton, M. P., Nielsen, M. M., Damas, N. D., Damgaard, C. K., et al. (2016). SNHG16 Is Regulated by the Wnt Pathway in Colorectal Cancer and Affects Genes Involved in Lipid Metabolism. Mol. Oncol. 10, 1266–1282. doi:10.1016/j.molonc.2016.06.003
D'Addario, M., Arora, P. D., Ellen, R. P., and McCulloch, C. A. G. (2002). Interaction of P38 and Sp1 in a Mechanical Force-Induced, β1 Integrin-Mediated Transcriptional Circuit that Regulates the Actin-Binding Protein Filamin-A. J. Biol. Chem. 277, 47541–47550. doi:10.1074/jbc.M207681200
D'haeseleer, P. (2006). What Are DNA Sequence Motifs? Nat. Biotechnol. 24, 423–425. doi:10.1038/nbt0406-423
Deng, R., Zhang, J., and Chen, J. (2019). lncRNA SNHG1 Negatively Regulates miRNA‑101‑3p to Enhance the Expression of ROCK1 and Promote Cell Proliferation, Migration and Invasion in Osteosarcoma. Int. J. Mol. Med. 43, 1157–1166. doi:10.3892/ijmm.2018.4039
Diercke, K., Kohl, A., Lux, C. J., and Erber, R. (2011). Strain-dependent Up-Regulation of Ephrin-B2 Protein in Periodontal Ligament Fibroblasts Contributes to Osteogenesis during Tooth Movement. J. Biol. Chem. 286, 37651–37664. doi:10.1074/jbc.M110.166900
Driller, K., Pagenstecher, A., Uhl, M., Omran, H., Berlis, A., Grunder, A., et al. (2007). Nuclear Factor I X Deficiency Causes Brain Malformation and Severe Skeletal Defects. Mol. Cel. Biol. 27, 3855–3867. doi:10.1128/MCB.02293-06
Driver, J., Weber, C. E., Callaci, J. J., Kothari, A. N., Zapf, M. A., Roper, P. M., et al. (2015). Alcohol Inhibits Osteopontin-dependent Transforming Growth Factor-Β1 Expression in Human Mesenchymal Stem Cells. J. Biol. Chem. 290, 9959–9973. doi:10.1074/jbc.M114.616888
Fang, F., Zhang, K., Chen, Z., and Wu, B. (2019). Noncoding RNAs: New Insights into the Odontogenic Differentiation of Dental Tissue-Derived Mesenchymal Stem Cells. Stem Cel Res. Ther. 10. doi:10.1186/s13287-019-1411-x
Franceschetti, T., Kessler, C. B., Lee, S.-K., and Delany, A. M. (2013). MiR-29 Promotes Murine Osteoclastogenesis by Regulating Osteoclast Commitment and Migration. J. Biol. Chem. 288, 33347–33360. doi:10.1074/jbc.M113.484568
Furumatsu, T., Matsumoto, E., Kanazawa, T., Fujii, M., Lu, Z., Kajiki, R., et al. (2013). Tensile Strain Increases Expression of CCN2 and COL2A1 by Activating TGF-β-Smad2/3 Pathway in Chondrocytic Cells. J. Biomech. 46, 1508–1515. doi:10.1016/j.jbiomech.2013.03.028
Gao, J. Y., Yu, X. Q., and Wang, J. Q. (2018). KLF5 Modulates Proliferation and Osteogenic Differentiation of Human Periodontal Ligament Cells Subjected to Cyclic Tensile Stress. Shanghai Kou Qiang Yi Xue 27, 28–33.
Gao, Q., Walmsley, A. D., Cooper, P. R., and Scheven, B. A. (2016). Ultrasound Stimulation of Different Dental Stem Cell Populations: Role of Mitogen-Activated Protein Kinase Signaling. J. Endodontics 42, 425–431. doi:10.1016/j.joen.2015.12.019
Gao, X., Ge, J., Li, W., Zhou, W., and Xu, L. (2018). LncRNA KCNQ1OT1 Promotes Osteogenic Differentiation to Relieve Osteolysis via Wnt/β-Catenin Activation. Cell Biosci. 8, 19. doi:10.1186/s13578-018-0216-4
Gao, Y., Zhao, H., and Li, Y. (2019). LncRNA MCM3AP-AS1 Regulates miR-142-3p/HMGB1 to Promote LPS-Induced Chondrocyte Apoptosis. BMC Musculoskelet. Disord. 20, 1–7. doi:10.1186/s12891-019-2967-4
Gong, Y. Y., Peng, M. Y., Yin, D. Q., and Yang, Y. F. (2018). Long Non-coding RNA H19 Promotes the Osteogenic Differentiation of Rat Ectomesenchymal Stem Cells via Wnt/β-Catenin Signaling Pathway. Eur. Rev. Med. Pharmacol. Sci. 22, 8805–8813. doi:10.26355/eurrev_201812_16648
Gu, H., Li, Z., Lv, X. F., Zhao, A. B., Zhu, M. Y., and Zhang, Y. (2019). LncRNA KCNQ1OT1 Delayed Fracture Healing through the Wnt/β-Catenin Pathway. Eur. Rev. Med. Pharmacol. Sci. 23, 4575–4583. doi:10.26355/eurrev_201906_18034
Gu, L., Liu, J., Xu, D., and Lu, Y. (2019). Reciprocal Feedback Loop of the MALAT1-MicroRNA-194-YAP1 Pathway Regulates Progression of Acute Pancreatitis. Med. Sci. Monit. 25, 6894–6904. doi:10.12659/MSM.915598
Guo, Y., Wang, Y., Liu, Y., Liu, Y., Zeng, Q., Zhao, Y., et al. (2015). MicroRNA-218, microRNA-191*, microRNA-3070a and microRNA-33 Are Responsive to Mechanical Strain Exerted on Osteoblastic Cells. Mol. Med. Rep. 12, 3033–3038. doi:10.3892/mmr.2015.3705
Hao, J., Zhang, Y., Jing, D., Shen, Y., Tang, G., Huang, S., et al. (2015). Mechanobiology of Mesenchymal Stem Cells: Perspective into Mechanical Induction of MSC Fate. Acta Biomater. 20, 1–9. doi:10.1016/j.actbio.2015.04.008
He, Q., Yang, S., Gu, X., Li, M., Wang, C., and Wei, F. (2018). Long Noncoding RNA TUG1 Facilitates Osteogenic Differentiation of Periodontal Ligament Stem Cells via Interacting with Lin28A. Cell Death Dis. 9. doi:10.1038/s41419-018-0484-2
Heo, S.-J., Thorpe, S. D., Driscoll, T. P., Duncan, R. L., Lee, D. A., and Mauck, R. L. (2015). Biophysical Regulation of Chromatin Architecture Instills a Mechanical Memory in Mesenchymal Stem Cells. Sci. Rep. 5, 1–14. doi:10.1038/srep16895
Hirose, T., Virnicchi, G., Tanigawa, A., Naganuma, T., Li, R., Kimura, H., et al. (2014). NEAT1 Long Noncoding RNA Regulates Transcription via Protein Sequestration within Subnuclear Bodies. MBoC 25, 169–183. doi:10.1091/mbc.E13-09-0558
Hu, J., Liao, H., Ma, Z., Chen, H., Huang, Z., Zhang, Y., et al. (2016). Focal Adhesion Kinase Signaling Mediated the Enhancement of Osteogenesis of Human Mesenchymal Stem Cells Induced by Extracorporeal Shockwave. Sci. Rep. 6, 1–12. doi:10.1038/srep20875
Hu, Z., Wang, H., Wang, Y., Zhou, H., Shi, F., Zhao, J., et al. (2017). Genome-wide Analysis and Prediction of Functional Long Noncoding RNAs in Osteoblast Differentiation under Simulated Microgravity. Mol. Med. Rep. 16, 8180–8188. doi:10.3892/mmr.2017.7671
Hu, Z., Wang, Y., Sun, Z., Wang, H., Zhou, H., Zhang, L., et al. (2015). MiRNA-132-3p Inhibits Osteoblast Differentiation by Targeting Ep300 in Simulated Microgravity. Sci. Rep. 5, 1–13. doi:10.1038/srep18655
Huang, H., Yang, R., and Zhou, Y.-H. (2018). Mechanobiology of Periodontal Ligament Stem Cells in Orthodontic Tooth Movement. Stem Cell Int. 2018, 1, 7. doi:10.1155/2018/6531216
Huang, Y., Liu, H., Guo, R., Han, Y., Yang, Y., Zhao, Y., et al. (2021). Long Non-coding RNA FER1L4 Mediates the Autophagy of Periodontal Ligament Stem Cells under Orthodontic Compressive Force via AKT/FOXO3 Pathway. Front. Cell Dev. Biol. 9, 631181. doi:10.3389/fcell.2021.631181
Huang, Y., Ren, K., Yao, T., Zhu, H., Xu, Y., Ye, H., et al. (2020). MicroRNA-25-3p Regulates Osteoclasts through Nuclear Factor I X. Biochem. Biophysical Res. Commun. 522, 74–80. doi:10.1016/j.bbrc.2019.11.043
Huang, Y., Zhang, Y., Li, X., Liu, H., Yang, Q., Jia, L., et al. (2019). The Long Non-coding RNA Landscape of Periodontal Ligament Stem Cells Subjected to Compressive Force. Eur. J. Orthod. 41, 333–342. doi:10.1093/ejo/cjy057
Iwawaki, Y., Mizusawa, N., Iwata, T., Higaki, N., Goto, T., Watanabe, M., et al. (2015). MiR-494-3p Induced by Compressive Force Inhibits Cell Proliferation in MC3T3-E1 Cells. J. Biosci. Bioeng. 120, 456–462. doi:10.1016/j.jbiosc.2015.02.006
Izadpanahi, M., Seyedjafari, E., Arefian, E., Hamta, A., Hosseinzadeh, S., Kehtari, M., et al. (2018). Nanotopographical Cues of Electrospun PLLA Efficiently Modulate Non-coding RNA Network to Osteogenic Differentiation of Mesenchymal Stem Cells during BMP Signaling Pathway. Mater. Sci. Eng. C 93, 686–703. doi:10.1016/j.msec.2018.08.023
Ji, H., Tang, H., Lin, H., Mao, J., Gao, L., Liu, J., et al. (2014). Rho/Rock Cross-Talks with Transforming Growth Factor-β/Smad Pathway Participates in Lung Fibroblast-Myofibroblast Differentiation. Biomed. Reports 2, 787–792. doi:10.3892/br.2014.323
Ji, Y., Gu, Y., Hong, S., Yu, B., Zhang, J. H., and Liu, J. N. (2020). Comprehensive Analysis of lncRNA‐TF Crosstalks and Identification of Prognostic Regulatory Feedback Loops of Glioblastoma Using lncRNA/TF‐mediated ceRNA Network. J. Cel. Biochem. 121, 755–767. doi:10.1002/jcb.29321
Jiang, Y., Wu, W., Jiao, G., Chen, Y., and Liu, H. (2019). LncRNA SNHG1 Modulates P38 MAPK Pathway through Nedd4 and Thus Inhibits Osteogenic Differentiation of Bone Marrow Mesenchymal Stem Cells. Life Sci. 228, 208–214. doi:10.1016/j.lfs.2019.05.002
Jin, D., Guo, J., Wu, Y., Du, J., Yang, L., Wang, X., et al. (2019). M6A mRNA Methylation Initiated by METTL3 Directly Promotes YAP Translation and Increases YAP Activity by Regulating the MALAT1-miR-1914-3p-YAP axis to Induce NSCLC Drug Resistance and Metastasis. J. Hematol. Oncol. 12. doi:10.1186/s13045-019-0830-6
Ju, C., Liu, R., Zhang, Y.-W., Zhang, Y., Zhou, R., Sun, J., et al. (2019). Mesenchymal Stem Cell-Associated lncRNA in Osteogenic Differentiation. Biomed. Pharmacother. 115, 108912. doi:10.1016/j.biopha.2019.108912
Kanzaki, H., Wada, S., Yamaguchi, Y., Katsumata, Y., Itohiya, K., Fukaya, S., et al. (2019). Compression and Tension Variably Alter Osteoprotegerin Expression via miR-3198 in Periodontal Ligament Cells. BMC Mol. Cel Biol 20. doi:10.1186/s12860-019-0187-2
Kartha, R. V., and Subramanian, S. (2014). Competing Endogenous RNAs (ceRNAs): New Entrants to the Intricacies of Gene Regulation. Front. Genet. 5, 8. doi:10.3389/fgene.2014.00008
Kawasaki, N., Miwa, T., Hokari, S., Sakurai, T., Ohmori, K., Miyauchi, K., et al. (2018). Long Noncoding RNANORADregulates Transforming Growth Factor-β Signaling and Epithelial-To-Mesenchymal Transition-like Phenotype. Cancer Sci. 109, 2211–2220. doi:10.1111/cas.13626
Kozomara, A., Birgaoanu, M., and Griffiths-Jones, S. (2019). miRBase: from microRNA Sequences to Function. Nucleic Acids Res. 47, D155–D162. doi:10.1093/nar/gky1141
Le, H. Q., Ghatak, S., Yeung, C.-Y. C., Tellkamp, F., Günschmann, C., Dieterich, C., et al. (2016). Mechanical Regulation of Transcription Controls Polycomb-Mediated Gene Silencing during Lineage Commitment. Nat. Cel Biol. 18, 864–875. doi:10.1038/ncb3387
Le Mée, S., Fromigué, O., and Marie, P. J. (2005). Sp1/Sp3 and the Myeloid Zinc finger Gene MZF1 Regulate the Human N-Cadherin Promoter in Osteoblasts. Exp. Cel Res. 302, 129–142. doi:10.1016/j.yexcr.2004.08.028
Lee, D.-S., Choung, H.-W., Kim, H.-J., Gronostajski, R. M., Yang, Y.-I., Ryoo, H.-M., et al. (2014). NFI-C Regulates Osteoblast Differentiation via Control of Osterix Expression. Stem Cells 32, 2467–2479. doi:10.1002/stem.1733
Lewis, B. P., Burge, C. B., and Bartel, D. P. (2005). Conserved Seed Pairing, Often Flanked by Adenosines, Indicates that Thousands of Human Genes Are microRNA Targets. Cell 120, 15–20. doi:10.1016/j.cell.2004.12.035
Li, B., Liu, J., Zhao, J., Ma, J.-X., Jia, H.-B., Zhang, Y., et al. (2017). LncRNA-H19 Modulates Wnt/β-Catenin Signaling by Targeting Dkk4 in Hindlimb Unloaded Rat. Orthop. Surg. 9, 319–327. doi:10.1111/os.12321
Li, B., Zhao, J., Ma, J.-x., Li, G.-m., Zhang, Y., Xing, G.-s., et al. (2018). Overexpression of DNMT1 Leads to Hypermethylation of H19 Promoter and Inhibition of Erk Signaling Pathway in Disuse Osteoporosis. Bone 111, 82–91. doi:10.1016/j.bone.2018.03.017
Li, J.-H., Liu, S., Zhou, H., Qu, L.-H., and Yang, J.-H. (2014). starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and Protein-RNA Interaction Networks from Large-Scale CLIP-Seq Data. Nucl. Acids Res.42, D92–D97. doi:10.1093/nar/gkt1248
Li, J., Hu, C., Han, L., Liu, L., Jing, W., Tang, W., et al. (2015). MiR-154-5p Regulates Osteogenic Differentiation of Adipose-Derived Mesenchymal Stem Cells under Tensile Stress through the Wnt/PCP Pathway by Targeting Wnt11. Bone 78, 130–141. doi:10.1016/j.bone.2015.05.003
Liao, Q., Liu, C., Yuan, X., Kang, S., Miao, R., Xiao, H., et al. 2011. Large-scale Prediction of Long Non-coding RNA Functions in a Coding-Non-Coding Gene Co-expression Network. Nucleic Acids Res. 39(9):3864–3878. doi:10.1093/nar/gkq1348
Liu, H., Huang, Y., Zhang, Y., Han, Y., Zhang, Y., Jia, L., et al. (2019). Long Noncoding RNA Expression Profile of Mouse Cementoblasts under Compressive Force. Angle Orthod. 89, 455–463. doi:10.2319/061118-438.1
Liu, L., Liu, M., Li, R., Liu, H., Du, L., Chen, H., et al. (2017). MicroRNA-503-5p Inhibits Stretch-Induced Osteogenic Differentiation and Bone Formation. Cell Biol. Int. 41, 112–123. doi:10.1002/cbin.10704
Liu, Q., Hu, X., Zhang, X., Dai, L., Duan, X., Zhou, C., et al. (2016). The TMSB4 Pseudogene LncRNA Functions as a Competing Endogenous RNA to Promote Cartilage Degradation in Human Osteoarthritis. Mol. Ther. 24, 1726–1733. doi:10.1038/mt.2016.151
Liu, S. C., Sun, Q. Z., Qiao, X. F., Li, X. G., Yang, J. H., Wang, T. Q., et al. (2019). LncRNA TUG1 Influences Osteoblast Proliferation and Differentiation through the Wnt/β-Catenin Signaling Pathway. Eur. Rev. Med. Pharmacol. Sci. 23, 4584–4590. doi:10.26355/eurrev_201906_18035
Luo, Y., Ge, R., Wu, H., Ding, X., Song, H., Ji, H., et al. (2019). The Osteogenic Differentiation of Human Adipose-Derived Stem Cells Is Regulated through the Let-7i-3p/LEF1/β-Catenin axis under Cyclic Strain. Stem Cel Res. Ther. 10. doi:10.1186/s13287-019-1470-z
Lv, L., Liu, Y., Zhang, P., Bai, X., Ma, X., Wang, Y., et al. (2018). The Epigenetic Mechanisms of Nanotopography-Guided Osteogenic Differentiation of Mesenchymal Stem Cells via High-Throughput Transcriptome Sequencing. Ijn Vol. 13, 5605–5623. doi:10.2147/IJN.S168928
Manokawinchoke, J., Pavasant, P., Sawangmake, C., Limjeerajarus, N., Limjeerajarus, C. N., Egusa, H., et al. (2019). Intermittent Compressive Force Promotes Osteogenic Differentiation in Human Periodontal Ligament Cells by Regulating the Transforming Growth Factor-β Pathway. Cel Death Dis. 10. doi:10.1038/s41419-019-1992-4
Mathelier, A., Fornes, O., Arenillas, D. J., Chen, C.-Y., Denay, G., Lee, J., et al. (2016). JASPAR 2016: a Major Expansion and Update of the Open-Access Database of Transcription Factor Binding Profiles. Nucleic Acids Res. 44, D110–D115. doi:10.1093/nar/gkv1176
McCarthy, D. J., Chen, Y., and Smyth, G. K. (2012). Differential Expression Analysis of Multifactor RNA-Seq Experiments with Respect to Biological Variation. Nucleic Acids Res. 40, 4288–4297. doi:10.1093/nar/gks042
Mercer, T. R., and Mattick, J. S. (2013). Structure and Function of Long Noncoding RNAs in Epigenetic Regulation. Nat. Struct. Mol. Biol. 20, 300–307. doi:10.1038/nsmb.2480
Mohan, S., Wergedal, J. E., Das, S., and Kesavan, C. (2015). Conditional Disruption of miR17-92 Cluster in Collagen Type I-Producing Osteoblasts Results in Reduced Periosteal Bone Formation and Bone Anabolic Response to Exercise. Physiol. Genomics 47, 33–43. doi:10.1152/physiolgenomics.00107.2014
Ørom, U. A., and Shiekhattar, R. (2013). Long Noncoding RNAs Usher in a new era in the Biology of Enhancers. Cell 154, 1190–1193. doi:10.1016/j.cell.2013.08.028
Panciera, T., Azzolin, L., Cordenonsi, M., and Piccolo, S. (2017). Mechanobiology of YAP and TAZ in Physiology and Disease. Nat. Rev. Mol. Cel Biol. 18, 758–770. doi:10.1038/nrm.2017.87
Petridou, N. I., Spiró, Z., and Heisenberg, C.-P. (2017). Multiscale Force Sensing in Development. Nat. Cel Biol. 19, 581–588. doi:10.1038/ncb3524
Ponting, C. P., Oliver, P. L., and Reik, W. (2009). Evolution and Functions of Long Noncoding RNAs. Cell 136, 629–641. doi:10.1016/j.cell.2009.02.006
Pylayeva-Gupta, Y. (2011). Transcription Factor Regulation by Mechanical Stress. Bone 23, 1–7. doi:10.1038/jid.2014.371
Qi, L., and Zhang, Y. (2014). The microRNA 132 Regulates Fluid Shear Stress-Induced Differentiation in Periodontal Ligament Cells through mTOR Signaling Pathway. Cell. Physiol. Biochem. 33, 433–445. doi:10.1159/000358624
Qin, L., Liu, W., Cao, H., and Xiao, G. (2020). Molecular Mechanosensors in Osteocytes. Bone Res. 8. doi:10.1038/s41413-020-0099-y
Qin, W., Liu, L., Wang, Y., Wang, Z., Yang, A., and Wang, T. (2019). Mir-494 Inhibits Osteoblast Differentiation by Regulating BMP Signaling in Simulated Microgravity. Endocrine 65, 426–439. doi:10.1007/s12020-019-01952-7
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26, 139–140. doi:10.1093/bioinformatics/btp616
Roh, S. Y., and Park, J.-C. (2017). The Role of Nuclear Factor I-C in Tooth and Bone Development. J. Korean Assoc. Oral Maxillofac. Surg. 43, 63–69. doi:10.5125/jkaoms.2017.43.2.63
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Stanley, A., Heo, S. j., Mauck, R. L., Mourkioti, F., and Shore, E. M. (2019). Elevated BMP and Mechanical Signaling through YAP1/RhoA Poises FOP Mesenchymal Progenitors for Osteogenesis. J. Bone Miner Res. 34, 1894, 1909. doi:10.1002/jbmr.3760
Sun, H., Peng, G., Ning, X., Wang, J., Yang, H., and Deng, J. (2019). Emerging Roles of Long Noncoding RNA in Chondrogenesis, Osteogenesis, and Osteoarthritis. Am. J. Transl. Res. 11, 16–30. Available at: www.ajtr.org [Accessed August 24, 2020].
Sun, Z., Cao, X., Hu, Z., Zhang, L., Wang, H., Zhou, H., et al. (2015). MiR-103 Inhibits Osteoblast Proliferation Mainly through Suppressing Cav1.2 Expression in Simulated Microgravity. Bone 76, 121–128. doi:10.1016/j.bone.2015.04.006
Tan, B.-S., Yang, M.-C., Singh, S., Chou, Y.-C., Chen, H.-Y., Wang, M.-Y., et al. (2019). LncRNA NORAD Is Repressed by the YAP Pathway and Suppresses Lung and Breast Cancer Metastasis by Sequestering S100P. Oncogene 38, 5612–5626. doi:10.1038/s41388-019-0812-8
Tan, G., and Lenhard, B. (2016). TFBSTools: an R/bioconductor Package for Transcription Factor Binding Site Analysis. Bioinformatics 32, 1555–1556. doi:10.1093/bioinformatics/btw024
Taniguchi Ishikawa, E., Chang, K. H., Nayak, R., Olsson, H. A., Ficker, A. M., Dunn, S. K., et al. (2013). Klf5 Controls Bone Marrow Homing of Stem Cells and Progenitors through Rab5-Mediated β1/β2-integrin Trafficking. Nat. Commun. 4, 1660. doi:10.1038/ncomms2645
Tay, Y., Rinn, J., and Pandolfi, P. P. (2014). The Multilayered Complexity of ceRNA Crosstalk and Competition. Nature 505, 344–352. doi:10.1038/nature12986
Teramura, T., Takehara, T., Onodera, Y., Nakagawa, K., Hamanishi, C., and Fukuda, K. (2012). Mechanical Stimulation of Cyclic Tensile Strain Induces Reduction of Pluripotent Related Gene Expressions via Activation of Rho/ROCK and Subsequent Decreasing of AKT Phosphorylation in Human Induced Pluripotent Stem Cells. Biochem. Biophysical Res. Commun. 417, 836–841. doi:10.1016/j.bbrc.2011.12.052
Thompson, W. R., Rubin, C. T., and Rubin, J. (2012). Mechanical Regulation of Signaling Pathways in Bone. Gene 503, 179–193. doi:10.1016/j.gene.2012.04.076
Tian, X., Wu, Y., Yang, Y., Wang, J., Niu, M., Gao, S., et al. (2020). Long Noncoding RNA LINC00662 Promotes M2 Macrophage Polarization and Hepatocellular Carcinoma Progression via Activating Wnt/β‐catenin Signaling. Mol. Oncol. 14, 462–483. doi:10.1002/1878-0261.12606
Todorovski, V., Fox, A. H., and Choi, Y. S. (2020). Matrix Stiffness-Sensitive Long Noncoding RNA NEAT1 Seeded Paraspeckles in Cancer Cells. MBoC 31, 1654–1662. doi:10.1091/mbc.E20-02-0097
Ventura, A. (2016). NORAD: Defender of the Genome. Trends Genet. 32, 390–392. doi:10.1016/j.tig.2016.04.002
Vlachos, I. S., Zagganas, K., Paraskevopoulou, M. D., Georgakilas, G., Karagkouni, D., Vergoulis, T., et al. (2015). DIANA-miRPath v3.0: Deciphering microRNA Function with Experimental Support. Nucleic Acids Res. 43, W460–W466. doi:10.1093/nar/gkv403
Wang, C.-G., Liao, Z., Xiao, H., Liu, H., Hu, Y.-H., Liao, Q.-D., et al. (2019). LncRNA KCNQ1OT1 Promoted BMP2 Expression to Regulate Osteogenic Differentiation by Sponging miRNA-214. Exp. Mol. Pathol. 107, 77–84. doi:10.1016/j.yexmp.2019.01.012
Wang, H., Feng, C., Li, M., Zhang, Z., Liu, J., and Wei, F. (2020). Analysis of lncRNAs‐miRNAs‐mRNAs Networks in Periodontal Ligament Stem Cells under Mechanical Force. Oral Dis., 27., 325–337. doi:10.1111/odi.13530
Wang, H., Sun, Z., Wang, Y., Hu, Z., Zhou, H., Zhang, L., et al. (2016). miR-33-5p, a Novel Mechano-Sensitive microRNA Promotes Osteoblast Differentiation by Targeting Hmga2. Sci. Rep. 6, 1–15. doi:10.1038/srep23170
Wang, J. L., Wei, X., Wang, A. G., Bai, Y., and Wu, X. J. (2020). KCNQ1OT1 Regulates Osteogenic Differentiation of hBMSC by miR-320a/Smad5 axis. Eur. Rev. Med. Pharmacol. Sci. 24, 2843–2854. doi:10.26355/eurrev_202003_20648
Wang, K., Wang, Y., Hu, Z., Zhang, L., Li, G., Dang, L., et al. (2020). Bone-targeted lncRNA OGRU Alleviates Unloading-Induced Bone Loss via miR-320-3p/Hoxa10 axis. Cel Death Dis. 11. doi:10.1038/s41419-020-2574-1
Wang, L., You, X., Lotinun, S., Zhang, L., Wu, N., and Zou, W. (2020). Mechanical Sensing Protein PIEZO1 Regulates Bone Homeostasis via Osteoblast-Osteoclast Crosstalk. Nat. Commun. 11, 282. doi:10.1038/s41467-019-14146-6
Wang, M., Ge, X., Zheng, Y., Wang, C., Zhang, Y., and Lin, Y. (2020). Microarray Analysis Reveals that lncRNA PWRN1 ‐209 Promotes Human Bone Marrow Mesenchymal Stem Cell Osteogenic Differentiation on Microtopography Titanium Surface In Vitro. J. Biomed. Mater. Res. 108, 2889–2902. doi:10.1002/jbm.b.34620
Wang, S., Zuo, H., Jin, J., Lv, W., Xu, Z., Fan, Y., et al. (2019). Long Noncoding RNA Neat1 Modulates Myogenesis by Recruiting Ezh2. Cel Death Dis. 10. doi:10.1038/s41419-019-1742-7
Wang, Y., Wang, K., Zhang, L., Tan, Y., Hu, Z., Dang, L., et al. (2020). Targeted Overexpression of the Long Noncoding RNA ODSM Can Regulate Osteoblast Function In Vitro and In Vivo. Cel Death Dis. 11, 133. doi:10.1038/s41419-020-2325-3
Wang, Y., Jia, L., Zheng, Y., and Li, W. (2018). Bone Remodeling Induced by Mechanical Forces Is Regulated by miRNAs. Biosci. Rep. 38. doi:10.1042/BSR20180448
Wei, F., Liu, D., Feng, C., Zhang, F., Yang, S., Hu, Y., et al. (2015). MicroRNA-21 Mediates Stretch-Induced Osteogenic Differentiation in Human Periodontal Ligament Stem Cells. Stem Cell Develop. 24, 312–319. doi:10.1089/scd.2014.0191
Wu, D., Yin, L., Sun, D., Wang, F., Wu, Q., Xu, Q., et al. (2020). Long Noncoding RNA TUG1 Promotes Osteogenic Differentiation of Human Periodontal Ligament Stem Cell through Sponging microRNA-222-3p to Negatively Regulate Smad2/7. Arch. Oral Biol. 117, 104814. doi:10.1016/j.archoralbio.2020.104814
Wu, J., Zhao, J., Sun, L., Pan, Y., Wang, H., and Zhang, W.-B. (2018). Long Non-coding RNA H19 Mediates Mechanical Tension-Induced Osteogenesis of Bone Marrow Mesenchymal Stem Cells via FAK by Sponging miR-138. Bone 108, 62–70. doi:10.1016/j.bone.2017.12.013
Wu, W., Guo, L., Liang, Z., Liu, Y., and Yao, Z. (2020). Lnc-SNHG16/miR-128 axis Modulates Malignant Phenotype through WNT/β-catenin Pathway in Cervical Cancer Cells. J. Cancer 11, 2201–2212. doi:10.7150/jca.40319
Wu, Y., Ou, Y., Liao, C., Liang, S., and Wang, Y. (2019b). High-throughput Sequencing Analysis of the Expression Profile of microRNAs and Target Genes in Mechanical Force-Induced Osteoblastic/cementoblastic Differentiation of Human Periodontal Ligament Cells. Am. J. Transl. Res. 11, 3398–3411.
Wu, Y., Jiang, Y., Liu, Q., and Liu, C.-Z. (2019a). LncRNA H19 Promotes Matrix Mineralization through Up-Regulating IGF1 by Sponging miR-185-5p in Osteoblasts. BMC Mol. Cel Biol 20. doi:10.1186/s12860-019-0230-3
Xiang, J., Fu, H. Q., Xu, Z., Fan, W. J., Liu, F., and Chen, B. (2020). lncRNA SNHG1 Attenuates Osteogenic Differentiation via the miR‑101/DKK1 axis in Bone Marrow Mesenchymal Stem Cells. Mol. Med. Rep. 22, 3715–3722. doi:10.3892/mmr.2020.11489
Xiao, X., Zhou, T., Guo, S., Guo, C., Zhang, Q., Dong, N., et al. (2017). LncRNA MALAT1 Sponges miR-204 to Promote Osteoblast Differentiation of Human Aortic Valve Interstitial Cells through Up-Regulating Smad4. Int. J. Cardiol. 243, 404–412. doi:10.1016/j.ijcard.2017.05.037
Xiaoling, G., Shuaibin, L., and Kailu, L. (2020). MicroRNA-19b-3p Promotes Cell Proliferation and Osteogenic Differentiation of BMSCs by Interacting with lncRNA H19. BMC Med. Genet. 21. doi:10.1186/s12881-020-0948-y
Yan, L., Liu, G., and Wu, X. (2021). Exosomes Derived from Umbilical Cord Mesenchymal Stem Cells in Mechanical Environment Show Improved Osteochondral Activity via Upregulation of LncRNA H19. J. Orthopaedic Translation 26, 111–120. doi:10.1016/j.jot.2020.03.005
Yang, X., Cai, J. B., Peng, R., Wei, C. Y., Lu, J. C., Gao, C., et al. (2019). The Long Noncoding RNA NORAD Enhances the TGF‐β Pathway to Promote Hepatocellular Carcinoma Progression by Targeting miR‐202‐5p. J. Cel. Physiol. 234, 12051–12060. doi:10.1002/jcp.27869
Yao, Q., Yang, J., Liu, T., Zhang, J., and Zheng, Y. (2019). Long Noncoding RNA MALAT1 Promotes the Stemness of Esophageal Squamous Cell Carcinoma by Enhancing YAP Transcriptional Activity. FEBS Open Bio 9, 1392–1402. doi:10.1002/2211-5463.12676
Yong, K. W., Choi, J. R., Choi, J. Y., and Cowie, A. C. (2020). Recent Advances in Mechanically Loaded Human Mesenchymal Stem Cells for Bone Tissue Engineering. Ijms 21, 5816–5820. doi:10.3390/ijms21165816
Yu, C., Li, L., Xie, F., Guo, S., Liu, F., Dong, N., et al. (2018). LncRNA TUG1 Sponges miR-204-5p to Promote Osteoblast Differentiation through Upregulating Runx2 in Aortic Valve Calcification. Cardiovasc. Res. 114, 168–179. doi:10.1093/cvr/cvx180
Yu, G., Wang, L.-G., Han, Y., and He, Q.-Y. (2012). clusterProfiler: an R Package for Comparing Biological Themes Among Gene Clusters. OMICS: A J. Integr. Biol. 16, 284–287. doi:10.1089/omi.2011.0118
Yu, Y., Chen, Y., Zhang, X., Lu, X., Hong, J., Guo, X., et al. (2018b). Knockdown of lncRNA KCNQ1OT1 Suppresses the Adipogenic and Osteogenic Differentiation of Tendon Stem Cell via Downregulating miR-138 Target Genes PPARγ and RUNX2. Cell Cycle 17, 2374–2385. doi:10.1080/15384101.2018.1534510
Yu, Y., Chen, Y., Zheng, Y.-j., Weng, Q., Zhu, S., Weng, Q.-h., et al. 2020). LncRNA TUG1 Promoted Osteogenic Differentiation through Promoting bFGF Ubiquitination. In Vitro Cell.Dev.Biol.-Animal 56, 42–48. doi:10.1007/s11626-019-00410-y
Yuan, Y., Guo, J., Zhang, L., Tong, X., Zhang, S., Zhou, X., et al. (2019). MiR-214 Attenuates the Osteogenic Effects of Mechanical Loading on Osteoblasts. Int. J. Sports Med. 40, 931–940. doi:10.1055/a-1015-0285
Yuan, Y., Zhang, L., Tong, X., Zhang, M., Zhao, Y., Guo, J., et al. (2017). Mechanical Stress Regulates Bone Metabolism through MicroRNAs. J. Cel. Physiol. 232, 1239–1245. doi:10.1002/jcp.25688
Zeng, Q., Wang, Y., Gao, J., Yan, Z., Li, Z., Zou, X., et al. (2019). miR-29b-3p Regulated Osteoblast Differentiation via Regulating IGF-1 Secretion of Mechanically Stimulated Osteocytes. Cell. Mol. Biol. Lett. 24. doi:10.1186/s11658-019-0136-2
Zhang, L., Li, G., Wang, K., Wang, Y., Dong, J., Wang, H., et al. (2020). MiR-30 Family Members Inhibit Osteoblast Differentiation by Suppressing Runx2 under Unloading Conditions in MC3T3-E1 Cells. Biochem. Biophysical Res. Commun. 522, 164–170. doi:10.1016/j.bbrc.2019.11.057
Zhang, T., Wen, F., Wu, Y., Goh, G. S. H., Ge, Z., Tan, L. P., et al. (2015). Cross-talk between TGF-Beta/SMAD and Integrin Signaling Pathways in Regulating Hypertrophy of Mesenchymal Stem Cell Chondrogenesis under Deferral Dynamic Compression. Biomaterials 38, 72–85. doi:10.1016/j.biomaterials.2014.10.010
Zhang, W., Dong, R., Diao, S., Du, J., Fan, Z., and Wang, F. (2017). Differential Long Noncoding RNA/mRNA Expression Profiling and Functional Network Analysis during Osteogenic Differentiation of Human Bone Marrow Mesenchymal Stem Cells. Stem Cel Res. Ther. 8, 1–13. doi:10.1186/s13287-017-0485-6
Zhang, X., Zhao, Y., Zhao, Z., Han, X., and Chen, Y. (2019). Knockdown of DANCR Reduces Osteoclastogenesis and Root Resorption Induced by Compression Force via Jagged1. Cell Cycle 18, 1759–1769. doi:10.1080/15384101.2019.1632638
Zhang Y, Y., Chen, B., Li, D., Zhou, X., and Chen, Z. (2019). LncRNA NEAT1/miR-29b-3p/BMP1 axis Promotes Osteogenic Differentiation in Human Bone Marrow-Derived Mesenchymal Stem Cells. Pathol. - Res. Pract. 215, 525–531. doi:10.1016/j.prp.2018.12.034
Zhang, Y., Dong, Q., and Sun, X. (2020b). Positive Feedback Loop LINC00511/miR-150-5p/SP1 Modulates Chondrocyte Apoptosis and Proliferation in Osteoarthritis. DNA Cel Biol. 39, 1506–1512. doi:10.1089/dna.2020.5718
Zhang, Z.-K., Li, J., Guan, D., Liang, C., Zhuo, Z., Liu, J., et al. (2018). Long Noncoding RNA lncMUMA Reverses Established Skeletal Muscle Atrophy Following Mechanical Unloading. Mol. Ther. 26, 2669–2680. doi:10.1016/j.ymthe.2018.09.014
Zheng, D., Wang, B., Zhu, X., Hu, J., Sun, J., Xuan, J., et al. (2019a). LncRNA OIP5-AS1 Inhibits Osteoblast Differentiation of Valve Interstitial Cells via miR-137/TWIST11 axis. Biochem. Biophysical Res. Commun. 511, 826–832. doi:10.1016/j.bbrc.2019.02.109
Zheng, S., Wang, Y. B., Yang, Y. L., Chen, B. P., Wang, C. X., Li, R. H., et al. (2019b). LncRNA MALAT1 Inhibits Osteogenic Differentiation of Mesenchymal Stem Cells in Osteoporosis Rats through MAPK Signaling Pathway. Eur. Rev. Med. Pharmacol. Sci. 23, 4609–4617. doi:10.26355/eurrev_201906_18038
Zhou, J., Wang, S., Qi, Q., Yang, X., Zhu, E., Yuan, H., et al. (2017). Nuclear Factor I‐C Reciprocally Regulates Adipocyte and Osteoblast Differentiation via Control of Canonical Wnt Signaling. FASEB j. 31, 1939–1952. doi:10.1096/fj.201600975RR
Zhou, K., Ou, Q., Wang, G., Zhang, W., Hao, Y., and Li, W. (2019). High Long Non-coding RNA NORAD Expression Predicts Poor Prognosis and Promotes Breast Cancer Progression by Regulating TGF-β Pathway. Cancer Cel Int. 19, 63. doi:10.1186/s12935-019-0781-6
Zhou, T., Gao, B., Fan, Y., Liu, Y., Feng, S., Cong, Q., et al. (2020). Piezo1/2 Mediate Mechanotransduction Essential for Bone Formation through Concerted Activation of NFAT-YAP1-SS-Catenin. Elife 9, 1–38. doi:10.7554/eLife.52779
Zhu, J., Yu, W., Wang, Y., Xia, K., Huang, Y., Xu, A., et al. (2019). LncRNAs: Function and Mechanism in Cartilage Development, Degeneration, and Regeneration. Stem Cel Res. Ther. 10, 1–12. doi:10.1186/s13287-019-1458-8
Keywords: mechanotransduction, long noncoding RNA (lncRNA), competing endogenous RNA (ceRNA), transcription factor (TF), bone, mesenchymal stem cells (MSCs), microRNA (miRNAs), force
Citation: Cai J, Li C, Li S, Yi J, Wang J, Yao K, Gan X, Shen Y, Yang P, Jing D and Zhao Z (2022) A Quartet Network Analysis Identifying Mechanically Responsive Long Noncoding RNAs in Bone Remodeling. Front. Bioeng. Biotechnol. 10:780211. doi: 10.3389/fbioe.2022.780211
Received: 20 September 2021; Accepted: 20 January 2022;
Published: 09 March 2022.
Edited by:
Maurizio Pesce, Monzino Cardiology Center (IRCCS), ItalyReviewed by:
Mar Cóndor, Imec, BelgiumGiada Andrea Cassanmagnago, Humanitas Research Hospital, Italy
Copyright © 2022 Cai, Li, Li, Yi, Wang, Yao, Gan, Shen, Yang, Jing and Zhao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dian Jing, janetjingdian@hotmail.com; Zhihe Zhao, zhzhao@scu.edu.cn