Long Non-coding RNAs in Traumatic Brain Injury Accelerated Fracture Healing

It is commonly observed that patients with bone fracture concomitant with traumatic brain injury (TBI) had significantly increased fracture healing, but the underlying mechanisms were not fully revealed. Long non-coding RNAs (lncRNAs) are known to play complicated roles in bone homeostasis, but their role in TBI accelerated fracture was rarely reported. The present study was designed to determine the role of lncRNAs in TBI accelerated fracture via transcriptome sequencing and further bioinformatics analyses. Blood samples from three fracture-only patients, three fracture concomitant with TBI patients, and three healthy controls were harvested and were subsequently subjected to transcriptome lncRNA sequencing. Differentially expressed genes were identified, and pathway enrichment was performed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. High-dimensional data visualization by self-organizing map (SOM) machine learning was applied to further interpret the data. An xCell method was then used to predict cellular behavior in all samples based on gene expression profiles, and an lncRNA–cell interaction network was generated. A total of 874 differentially expressed genes were identified, of which about 26% were lncRNAs. Those identified lncRNAs were mainly enriched on TBI-related and damage repair-related pathways. SOM analyses revealed that those differentially expressed lncRNAs could be divided into three major module implications and were mainly enriched on transcriptional regulation and immune-related signal pathways, which promote us to further explore cellular behaviors based on differentially expressed lncRNAs. We have predicted that basophils, CD8+ T effector memory cells, B cells, and naïve B cells were significantly downregulated, while microvascular endothelial cells were predicted to be significantly upregulated in the Fr/TBI group, was the lowest and highest, respectively. ENSG00000278905, ENSG00000240980, ENSG00000255670, and ENSG00000196634 were the most differentially expressed lncRNAs related to all changes of cellular behavior. The present study has revealed for the first time that several critical lncRNAs may participate in TBI accelerated fracture potentially via regulating cellular behaviors of basophils, cytotoxic T cells, B cells, and endothelial cells.

Introduction

Fractures, the most common injuries of the musculoskeletal system of humans, are a major cause of global mortality and disability (1). Epidemiological investigation revealed that roughly 4.39 million Chinese have a traumatic fracture in 2014 (2). Fracture healing is a complex process that requires multiple cell involvement and signal transduction. Generally, as a mineralized tissue, bone was intended to resist tension, compression, and torsional stresses (1, 3). A fracture occurs when the deforming energy exceeds the maximum threshold of bone structure (4, 5). The initial response of bone injury was vascular disruption-induced hematoma and clot formation at the fracture site. This danger signal subsequently induces the infiltration of immune cells and the release of inflammatory cytokine, chemokines, and growth factors. This bone injury further activates stem/progenitor cells involved in bone regeneration at the periosteum and in the bone marrow. The stem/progenitor cells were then differentiated into chondrocyte lineage or osteoblast lineage to promote the formation of cartilage for repairing the injured bones. Cartilage callus tissue further undergoes calcification, and blood vessels invade so that bone could be formatted and remodeled secondarily (1, 5, 6). Although massive progress has been made for fracture treatment, ~10% of fractures do not heal as expected, termed as non-union. Other compromised fracture healing like delayed union and inappropriate tissue regeneration was also commonly seen in the clinic (5, 6). Therefore, the development of novel molecules for predicting or treating fracture that could be utilized in the clinic is still of great importance.

In clinical practice, we and other orthopedists found that accelerated fracture healing could be observed in patients with traumatic brain injury (TBI) (7–9). As previously reported, significantly accelerated callus formation was found in cases concomitant with TBI when compared with those in patients with fracture alone (7–9). The possible working mechanism of this clinical phenomenon was also little studied (8, 9). Currently, TBI accelerated fracture healing was observed to be highly associated with changes in mesenchymal stem cells, cytokines and growth factors, hormones, and other substances like calcium and phosphate (9). However, most of the abovementioned mechanisms were preliminary ones that originated from preclinical studies. There is still little known by which cytokines, growth factors, and hormones that mediated the signal pathway finally promote fracture healing.

Long non-coding RNAs (lncRNAs), a family of transcripts with strand length exceeding 200 bp that do not encode protein, were demonstrated to play critical roles in regulating gene transcription, epigenetic regulation, cell proliferation, and cell differentiation (10–12). According to GENCODE annotation for humans, over 16,000 lncRNAs were identified, and those lncRNAs were reported to participate in numerous diseases and were not considered as novel targets/markers for disease diagnosis, treatment, or prognosis prediction. As previously reported, lncRNAs were vital in maintaining bone homeostasis since they could regulate the expression of key enzymes and proteins in the process of bone formation and resorption. For example, H19, an upregulated lncRNA during osteogenesis in bone mesenchymal stem cells (BMSCs), together with forkhead box C2, was proved to have the ability to promote osteogenic differentiation of BMSCs (13). Another lncRNA, differentiation antagonizing non-protein coding RNA (DANCR), was also demonstrated to inhibit osteogenic differentiation of human BMSCs (14). In animal models of fracture, LOC103691336 was reported to modulate osteogenic differentiation in a rat femoral fracture model (15); downregulated lncRNA-ANCR was found to promote tibial fracture healing in rabbits by regulating the expression of runt-related transcription factor 2 (RUNX2) (16); in mouse tibia fracture models, silencing or inhibiting of maternally expressed gene 3 (MEG3) (17) or TSIX (18) was proved to accelerate fracture healing. However, the role of lncRNAs in TBI accelerated fracture healing has never been reported before. The present study was designed to identify differentially expressed lncRNA via high-throughput sequencing and to further determine their roles in TBI accelerated fracture healing by bioinformatics analysis. The present study might form the experimental basis of potential targets of lncRNA for fracture treatment.

Patients and Methods

Demographic Data

The research protocol was approved by the ethical committee of Affiliated Hospital of Zunyi Medical University and was performed in compliance with the Declaration of Helsinki. Approximately 30 TBI patients, 100 fracture patients, and 50 TBI concomitant with fracture patients are admitted into our hospital monthly. Due to funding limitations, the sample size for this study was limited to three TBI concomitant with fracture patients. These were randomly selected in accordance with the method of random number table based on the inclusion and exclusion criteria. Likewise, three patients matched for age, gender, and fracture sites without concomitant TBI were also selected in accordance with the method of random number table. To avoid potential selection bias, randomization of selection was performed by a member of the study team who is not on the staff of our hospital. Three fracture-only patients (two male and one female, average age was 45.33 ± 6.51 years, defined as the FrO group) and three fracture concomitant with TBI patients (two male and one female, average age was 46.67 ± 2.88 years, defined as the Fr/TBI group) admitted in the Department of Emergency, Affiliated Hospital of Zunyi Medical University from January 2019 to March 2019, were involved in the current research (as presented in Table 1). Three age- and sex-matched healthy controls (two male and one female, average age was 47.33 ± 6.11 years, defined as HC) were all medical staff from our department. Patients in the Fr/TBI group have significantly accelerated fracture healing time than those patients in the FrO group (Table 1). Participants involved in the current study have signed an informed consent form.

TABLE 1

Table 1. Demographic data.

Fractures were diagnosed by X-ray photography and computed tomography (CT). All fractures were treated by surgery with internal fixation as performed by the same surgeon. TBI was diagnosed by CT, and the severity of TBI was evaluated by Glasgow Coma Scale (GCS). All TBI patients were treated with individualized protocol based on patients' conditions. Fracture healing time was recorded when bone union emerged as determined by X-ray photography.

Inclusion and Exclusion Criteria

Patients of the FrO or Fr/TBI group were included or excluded by the following criteria. Inclusion criteria: (1) clear fracture determined by X-ray photography and CT; (2) normal functions in critical organs like the heart (normal cardiac function as determined by electrocardiography and echocardiography), liver (indicated by serum biochemical tests), and kidney (indicated by serum and urine biochemical tests); and (3) GCS score of 9 to 12 for TBI patients. Exclusion criteria: (1) with severe infection complications; (2) immunodeficiency; (3) complications with other acute or chronic diseases, like cardiovascular diseases and metabolic diseases; (4) severe TBI (GCS score of 8 or less); and 5) severe open fractures (type II and III Gustilo open fracture), or transverse fractures.

Long Non-coding RNA High-Throughput Sequencing

Blood samples were harvested 7 days after admission. Total RNA from 4 ml of peripheral whole blood of each sample was extracted by 12 ml of TRIzol reagent according to manufacturer's protocols. Purity, concentration, and integrity of isolated RNA were further determined or verified to make sure the isolated total RNA was suitable for further analyses. Three micrograms of total RNA per sample was used to generate sequencing libraries after removal of ribosomal RNA. After cluster generation, the harvested libraries were sequenced on an Illumina HiSeq platform. Raw data of sequenced RNA in FASTAQ format were then subjected to quality control, and clean data were generated for further downstream analyses. Assembled transcripts were then annotated to reference genomics.

Gene Ontology and Kyoto Encyclopedia of Genes and Genomes Analyses

Differentially expressed genes or lncRNA target genes were subjected to Gene Ontology (GO) enrichment analysis by using clusterProfiler (19) R package. GO terms with a corrected P-value < 0.05 were considered to be significantly enriched.

Kyoto Encyclopedia of Genes and Genomes (KEGG) is a database resource for understanding high-level functions and utilities of a biological system, at the cell, organism, and ecosystem levels, using molecular-level information, especially large-scale molecular datasets generated by genome sequencing and other high-throughput experimental technologies (http://www.genome.jp/kegg/). An R package of clusterProfiler was used to determine the enrichment of differentially expressed genes or lncRNA target genes in KEGG pathways.

High-Dimensional Data Portraying by Self-Organizing Map Machine Learning

To better interpret harvested data, a machine learning-based self-organizing map (SOM) was used by utilizing R package of oposSOM (20, 21) (v2.8.0). The oposSOM was able to translate RNA-seq expression data into metadata of a reduced dimension. This package can provide various sample-centered and group-centered visualizations, sample similarity analyses, and functional enrichment analyses. In the current study, we adopted a sample-centered visualization of harvested lncRNAs, and map size was selected in accordance with recommendations from the developer. Pre-processed expression matrix was presented in a log scale and then subjected to SOM training from different gene expression matrices. Generally, a cluster of genes with similar expression profiles was taken as a metagene, which was selected by SOM training. The training process was terminated until all expression patterns in the data matrix metagenes were captured by metagenes. In the current study, a total of 10 expression modules were harvested and labeled from A to J based on K-means cluster. Neighbor-joining was used to obtain a dendrogram.

Inference of Infiltrating Cells From Gene Expression Profiles

To determine the fraction of infiltrated cells, we used xCell (22), which was a computational method based on ssGSEA and provides immune scores for 64 cell subtypes, spanning multiple innate and adaptive immune cells. It employs a spillover compensation technique to reduce dependencies between closely related cell types.

Long Non-coding RNA–Cell Interaction Network and Long Non-coding RNA–mRNA Co-expression Network Construction

Spearman's correlation analysis was applied to calculate the cor-relationship between the expression of lncRNAs and the changes of cellular behavior. Pearson's correlation analysis was used to calculate the cor-relationship between the expression of lncRNAs and mRNAs. |Correlation coefficient| > 0.5 and P < 0.05 was considered to be interdependent. Cytoscape (version 3.6.2) was used to present the network.

Statistical Analysis

All computational and statistical analyses were performed using R (version 3.5.1). An unpaired Student's t-test or Wilcoxon's test was used to compare the differences between two groups. One-way ANOVA was used to analyze the differences among three or more groups. P < 0.05 was considered statistically significant.

Results

Expression Profile of Sequence Analysis

To determine whether expressions of lncRNAs were changed in TBI accelerated fracture, an lncRNA-based transcriptome sequencing was performed. As presented in Figure 1, we identified a total of 874 genes that were differentially expressed among the three groups (Supplementary Table 1). Of all identified genes, only about 26% were lncRNAs. Among all identified lncRNAs (Figure 1B, Supplementary Table 1), the majority were large intergenic non-coding RNAs (lincRNAs) and antisense lncRNAs (Figure 1C). These differentially expressed mRNAs and lncRNAs were mainly involved in TBI-related and damage repair-related pathways (like homologous recombination and phagosome pathway analyzed by KEGG; as well as regulation of hair follicle development and telencephalon glial cell migration biological process analyzed by GO, Figures 1D,E).

FIGURE 1

Figure 1. Expression profile of sequence analysis. (A) Heatmap showed all differentially expressed genes from healthy volunteers and patients. (B) Composition of identified differentially expressed genes from long non-coding RNA (lncRNA) sequencing. (C) Counts of each lncRNA subtype of identified lncRNAs. lincRNA, large intergenic non-coding RNA; TEC, to be experimentally confirmed lncRNA. (D) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of differentially expressed genes. (E) Gene Ontology (GO) enrichment analysis of differentially expressed genes.

Self-Organizing Map–Portrayal Expression Maps

As presented in the heatmap, a significantly changed expression pattern could be observed. To better interpret the sequencing results, we adopted a SOM machine learning-based data portraying method, which generated in a sample-centered manner. By using the neighbor-joining method based on the result of SOM, a dendrogram was obtained, which showed the state of development among the three groups, as expected (Figure 2A). To better identify the key different genes among the three groups, a total of 10 modules were generated based on similarities of expression profiles by K-means cluster (Figure 2B, Supplementary Table 1). Further analyses showed the expression pattern of each module (Figure 2C).

FIGURE 2

Figure 2. Self-organizing map (SOM)–portrayal expression maps. (A) Representative image of SOM analyses of the expression pattern of each sample. Dendrogram summarized the similarity of samples. (B) Representative image of gene expression modules. (C) Heatmap of expressions of modules in each sample.

Annotation of Harvested Expression Modules

According to the expression pattern of each module, the harvested modules were further divided into three groups (modules H and I were defined as module implications of C1; modules D, F, and G were defined as module implications of C2; modules A and B were defined as module implications of C3; while no expression changes were observed in modules C and E). In all subdivided module implications, normalized gene expressions were significantly downregulated or upregulated among the investigated groups. In module implications of C1, normalized gene expressions in FrO group were the lowest, while in the HC group, they were the highest (Figures 3A,B). In module implications of C2, normalized gene expressions in FrO group were the highest, while in the Fr/TBI group, they were the lowest (Figures 3A,B). In module implications of C3, normalized gene expressions in the HC group were the lowest, while in Fr/TBI group, they were the highest (Figures 3A,B). Further pathway enrichment analyses by KEGG revealed that top enriched pathways were transcriptional regulation and immune-related signal pathways (like the Fanconi anemia pathway, amoebiasis, and mismatch repair) in all subdivided module implications (Figure 3C). GO analysis also revealed the same tendency; as proved by the top enriched biological process, there was leukocyte differentiation in all subdivided module implications (Figure 3D).

FIGURE 3

Figure 3. Annotation of harvested expression modules. (A) Representative images of module implications subdivision. (B) Boxplot showed the normalized gene expressions of subdivided module implications in each group. P-values were obtained via Wilcoxon's test. ** and *** indicate P-value < 0.01 and P-value < 0.001, respectively. (C) Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of genes from each module implication. (D) Gene Ontology (GO) biological process enrichment analyses of genes from each module implication. The vertical and horizontal axes represent the KEGG pathways and different modules, respectively. The size and the color intensity of a circle represent gene number and –log10 (P-value), respectively.

Predicted Cellular Behavior

The abovementioned enrichment analysis led us further to determine the cellular response in TBI accelerated fracture healing. We have predicted several changed cellular behaviors, as basophils, CD8⁺ T effector memory cells, B cells, and naïve B cells were predicted to be significantly downregulated, while microvascular endothelial cells were predicted to be significantly upregulated in the Fr/TBI group, which were the lowest and highest, respectively (Figure 4A and Supplementary Table 2). Correlation analyses were utilized to determine the relationship between differentially expressed lncRNA and known cellular behavior, and we created an lncRNA–cell interaction network based on lncRNA from C3 (Figure 4B). The result revealed that most expressions of lncRNAs from C3 were correlated with these changes of cellular behavior (Supplementary Table 3). In addition, ENSG00000278905, ENSG00000240980, ENSG00000255670, and ENSG00000196634 were the most differentially expressed lncRNAs that related to changes of cellular behavior. As indicated in Genotype-Tissue Expression Project (GTEx, https://gtexportal.org/home/), those four lncRNAs were expressed in various human tissues (Supplementary Figure 1). To predict the potential function and targets of these four lncRNAs, we constructed lncRNA–mRNA co-expression network, and then GO analysis (Figure 4C) of this network revealed that the top four cellular behavior-related lncRNAs were enriched mainly on transcriptional regulation and immune-related biological process, suggesting that those lncRNAs may participate in TBI accelerated fracture healing via regulating immune-related signaling potentially through transcriptional regulation.

FIGURE 4

Figure 4. Predicted cellular behavior. (A) Representative heatmap is shown to predict significantly changed cellular behaviors in each group. P-values were obtained via one-way ANOVA; * and ** indicate P-value < 0.05 and P-value < 0.01, respectively. CD8+ Tem, CD8+ effector memory T-cells; mv endothelial cells, microvascular endothelial cells. (B) Representative image of interactions between identified long non-coding RNA (lncRNA) and changed cellular behavior. The larger size of lncRNA node indicates that more types of changed cellular behavior correlated with it. The thickness of a line represents the correlation coefficient. Solid and dotted lines showed the positive and negative relationships between the expression and changed cellular behavior, respectively. (C) Gene Ontology (GO) enrichment analysis of changed cellular behavior related to lncRNAs. The outer circle shows a scatter plot for each term of the correlation coefficient of the assigned target genes, and the z-score shows the number of positive genes minus the number of negative genes divided by the square root of the count.

Discussion

The presented results have revealed for the first time that lncRNAs may be important in regulating TBI accelerated fracture healing. It was reported that lncRNAs could be involved in the maintenance of bone homeostasis via regulating chromatin remodeling, transcriptional activation/repression, post-transcriptional modification, or interaction with microRNAs (10–12). Currently, the relationship of lncRNAs and fracture healing was not well-acknowledged, and lncRNAs were reported to participate in fracture healing through interaction with microRNAs or by regulating protein expression without fully clarified mechanisms (13–18). In the present study, we have identified four major lncRNAs, namely, ENSG00000278905 (AC106818.1), ENSG00000240980 (AC093904.1), ENSG00000255670 (AC007619.1), and ENSG00000196634 (LUADT1), which may critically participate in TBI accelerated fracture healing. Except for LUADT1, the other three were newly identified ones that may be involved in fracture healing and have not been reported before. Also, little is known about the specific role of LUADT1 in this process since its role was mainly studied in tumorigenesis (23). Current reports reveal that the role of LUADT1 in bone homeostasis was performed to uncover its role in the behavior of primary chondrocytes (24). Those researchers have found that LUADT1 could participate in chondrocyte apoptosis by regulating miR-34a/SIRT1 signal pathway (24). Therefore, the specific role of those four newly identified lncRNAs in fracture healing remains to be explored.

Although fracture healing is a complex process, the course of fracture could be divided into three partially overlapping phases, namely, inflammation, repair, and remodeling (3). Fracture-induced blood vessel rupture would result in hematoma (characterized by hypoxia and low pH), which further initiates immune cascades and promotes the immune cell infiltration to injury sites. According to published literature, the first accumulated immune cells are polymorphonuclear neutrophils (PMNs) (25) that attracted debris and signals from cell death. The short-lived PMNs were reported to secrete massive cytokines and chemokines [C-C motif chemokine 2 (CCL2), IL6, TNF-α, etc.] to subsequently promote the infiltration of other immune cells like macrophages and other lymphocytes (26, 27). In a rodent model of fracture, systematic PMN activation could impair fracture healing (28). Also, a report has found that enhanced fracture healing could result in neutropenia in animal models of fracture (29). In the present study, we have found that basophils were predicted to be significantly decreased in TBI accelerated fracture patients, which is consistent with previous findings. The macrophages infiltrated from local environment (osteomacs) or from circulation (inflammatory macrophage) promote intramembranous bone formation and regulate the process of endochondral ossification (30). In the current study, we have not found any predicted quality/quantity changes associated with macrophages, potentially suggesting that macrophages participate in TBI accelerated fracture at the current time setting.

Behaviors of T cells and B cells also alter the process of fracture healing. In a mouse model of fracture, recombination activating gene 1 knockout mouse (Rag-1^−/−, and T- and B-cell deficiency) showed delayed cartilage maturation and less mineralized tissue (31). However, another group has found that deficiency in the adaptive immune system could result in accelerated fracture healing (32). In a sheep fracture model, cytotoxic T cells (CD8⁺ T cells) were accumulated in a considerably increased abundance in delayed fracture healing group accompanied with enhanced inflammation (33), while in a rat model of fracture, depletion of cytotoxic T cells by neutralizing antibody did not affect bone formation by reducing the implant pull-out force and stiffness (34). It was reported that B cells can secrete osteoprotegerin to inhibit osteoclastogenesis to maintain bone homeostasis. The frequency of naïve B cells was decreased while plasmablasts were increased immediately after fracture, which suggested an inflammatory status of the body after fracture (35). At the remodeling phase of fracture healing, B cells with regulatory function critically affect tibia fracture healing as loss of B-cell regulatory function was observed to be associated with delayed healing in patients (36). In our current study, we have found that cytotoxic T cells, B cells, and naïve B cells were predicted to significantly decrease in patients of TBI accelerated fracture, which was not quite consistent with the previous findings. This paradoxical phenomenon could be explained by the following researches: on the one hand, cyclooxygenase-2 (COX-2) inhibitors (non-steroidal anti-inflammatory drugs) could significantly impair fracture healing (37); on the other hand, patients with diabetes mellitus (systematic inflammation) also have impaired fracture healing (38). It seems like a certain degree of inflammation is important in promoting fracture healing, while to what extent the inflammation is beneficial remains to be answered.

Another critical step of fracture healing is revascularization, which provides a structural basis for nutrient supply and osteoblast lineage fate specification (5). Endothelial cells were the major component of the neovascular. The interaction between those endothelial cells and osteoprogenitors would provide signals for osteoblast generation (39). Especially, type H endothelial cells not only were able to promote angiogenesis (40) but also could enhance osteogenesis through Notch signaling (41). In our current study, we have found that microvascular endothelial cells were predicted to be significantly enriched in patients of TBI accelerated fracture when compared with other groups, which was parallel with the positive role of endothelial cells in fracture healing.

This study has several limitations. Firstly, due to sample availability and funding limitation, the present study only includes three patients in each group, and sequencing was not performed on patients with TBI alone. Secondly, the diagnostic or prognostic value of identified lncRNAs was not validated in the current study. Thirdly, the source of identified lncRNAs was not determined. Fourthly, identified lncRNAs in the current study were not fully investigated to determine whether they were consequences or reasons of TBI accelerated fracture. Finally, whether severity of TBI affects fracture healing in a different manner was also not included in the current study.

Conclusions

The present study has identified several lncRNAs that may participate in TBI accelerated fracture via regulating behaviors of basophils, cytotoxic T cells, B cells, and endothelial cells. However, massive experimental validations were strongly encouraged to promote the reliability of the current study.

Data Availability Statement

The data presented in the study are deposited in the Gene Expression Omnibus repository (https://www.ncbi.nlm.nih.gov/geo/), accession number GSE171718.

Ethics Statement

The studies involving human participants were reviewed and approved by Ethical Committee of Affiliated Hospital of Zunyi Medical University. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

GG, TW, GZ, and AY designed this research project and drafted and revised the manuscript. GG, YG, XJ, SW, GZ, and RW carried out the experiments, analyzed the data and drafted the manuscript. CL, GZ, and GY checked the analysis and revised the manuscript. All authors reviewed the manuscript and approved the manuscript that was submitted for publication.

Funding

This study was funded in part by the National Natural Science Foundation of China, grant number #81471862, and by the Science and Technology Bureau of Zunyi, grant number #HZ-2020-259.

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.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fsurg.2021.663377/full#supplementary-material

References

1. Einhorn TA, Gerstenfeld LC. Fracture healing: mechanisms and interventions. Nat Rev Rheumatol. (2015) 11:45–54. doi: 10.1038/nrrheum.2014164

CrossRef Full Text | Google Scholar

2. Chen W, Lv H, Liu S, Liu B, Zhu Y, Chen X, et al. National incidence of traumatic fractures in China: a retrospective survey of 512 187 individuals. Lancet Glob Health. (2017) 5:e807–17. doi: 10.1016/S2214-109X(17)30222-X

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Claes L, Recknagel S, Ignatius A. Fracture healing under healthy and inflammatory conditions. Nat Rev Rheumatol. (2012) 8:133–43. doi: 10.1038/nrrheum.20121

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Long F. Building strong bones: molecular regulation of the osteoblast lineage. Nat Rev Mol Cell Biol. (2011) 13:27–38. doi: 10.1038/nrm3254

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Salhotra A, Shah HN, Levi B, Longaker MT. Mechanisms of bone development and repair. Nat Rev Mol Cell Biol. (2020) 21:696–711. doi: 10.1038/s41580-020-00279-w

CrossRef Full Text | Google Scholar

6. Hadjiargyrou M, O'keefe RJ. The convergence of fracture repair and stem cells: interplay of genes, aging, environmental factors and disease. J Bone Miner Res. (2014) 29:2307–22. doi: 10.1002/jbmr2373

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Morley J, Marsh S, Drakoulakis E, Pape HC, Giannoudis PV. Does traumatic brain injury result in accelerated fracture healing? Injury. (2005) 36:363–8. doi: 10.1016/j.injury.2004.08028

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Zhao XG, Zhao GF, Ma YF, Jiang GY. Research progress in mechanism of traumatic brain injury affecting speed of fracture healing. Chin J Traumatol. (2007) 10:376–80.

PubMed Abstract | Google Scholar

9. Hofman M, Koopmans G, Kobbe P, Poeze M, Andruszkow H, Brink PR, et al. Improved fracture healing in patients with concomitant traumatic brain injury: proven or not? Mediators Inflamm. (2015) 2015:204842. doi: 10.1155/2015/204842

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Hassan MQ, Tye CE, Stein GS, Lian JB. Non-coding RNAs: epigenetic regulators of bone development and homeostasis. Bone. (2015) 81:746–56. doi: 10.1016/j.bone.2015.05026

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Sikora M, Marycz K, Smieszek A. Small and long non-coding RNAs as functional regulators of bone homeostasis, acting alone or cooperatively. Mol Ther Nucleic Acids. (2020) 21:792–803. doi: 10.1016/j.omtn.2020.07017

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Song W, Xie J, Li J, Bao C, Xiao Y. The emerging roles of long noncoding RNAs in bone homeostasis and their potential application in bone-related diseases. DNA Cell Biol. (2020) 39:926–37. doi: 10.1089/dna.20205391

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Zhou P, Li Y, Di R, Yang Y, Meng S, Song F, et al. H19 and Foxc2 synergistically promotes osteogenic differentiation of BMSCs via Wnt-beta-catenin pathway. J Cell Physiol. (2019) 234:13799–806. doi: 10.1002/jcp28060

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang J, Tao Z, Wang Y. Long noncoding RNA DANCR regulates the proliferation and osteogenic differentiation of human bone-derived marrow mesenchymal stem cells via the p38 MAPK pathway. Int J Mol Med. (2018) 41:213–9. doi: 10.3892/ijmm.20173215

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Li D, Yu K, Xiao T, Dai Y, Liu L, Li H, et al. LOC103691336/miR-138-5p/BMPR2 axis modulates Mg-mediated osteogenic differentiation in rat femoral fracture model and rat primary bone marrow stromal cells. J Cell Physiol. (2019) 234:21316–30. doi: 10.1002/jcp28736

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Liu ZC, Xu YL, Jiang Y, Liu Y, Wei ZC, Liu SG, et al. Low-expression of lncRNA-ANCR promotes tibial fracture healing via targeting RUNX2. Eur Rev Med Pharmacol Sci. (2019) 23:60–6. doi: 10.26355/eurrev_201908_18629

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Liu YB, Lin LP, Zou R, Zhao QH, Lin FQ. Silencing long non-coding RNA MEG3 accelerates tibia fraction healing by regulating the Wnt/beta-catenin signalling pathway. J Cell Mol Med. (2019) 23:3855–66. doi: 10.1111/jcmm14229

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Xu WW, Xu Y, Ji F, Ji Y, Wang QG. Inhibition of long non-coding RNA TSIX accelerates tibia fraction healing via binding and positively regulating the SOX6 expression. Eur Rev Med Pharmacol Sci. (2020) 24:4070–9. doi: 10.26355/eurrev_202004_20983

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. (2012) 16:284–7. doi: 10.1089/omi.20110118

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Loffler-Wirth H, Kalcher M, Binder H. oposSOM: R-package for high-dimensional portraying of genome-wide expression landscapes on bioconductor. Bioinformatics. (2015) 31:3225–7. doi: 10.1093/bioinformatics/btv342

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Loeffler-Wirth H, Reikowski J, Hakobyan S, Wagner J, Binder H. oposSOM-Browser: an interactive tool to explore omics data landscapes in health science. BMC Bioinform. (2020) 21:465. doi: 10.1186/s12859-020-03806-w

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Aran D, Hu Z, Butte AJ. xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. (2017) 18:220. doi: 10.1186/s13059-017-1349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Qiu M, Xu Y, Wang J, Zhang E, Sun M, Zheng Y, et al. A novel lncRNA, LUADT1, promotes lung adenocarcinoma proliferation via the epigenetic suppression of p27. Cell Death Dis. (2015) 6:e1858. doi: 10.1038/cddis.2015203

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Ni S, Xu C, Zhuang C, Zhao G, Li C, Wang Y, et al. LncRNA LUADT1 regulates miR-34a/SIRT1 to participate in chondrocyte apoptosis. J Cell Biochem. (2020) 1–6. doi: 10.1002/jcb29637

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Chung R, Cool JC, Scherer MA, Foster BK, Xian CJ. Roles of neutrophil-mediated inflammatory response in the bony repair of injured growth plate cartilage in young rats. J Leukoc Biol. (2006) 80:1272–80. doi: 10.1189/jlb0606365

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Andrew JG, Andrew SM, Freemont AJ, Marsh DR. Inflammatory cells in normal human fracture healing. Acta Orthop Scand. (1994) 65:462–6. doi: 10.3109/17453679408995493

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Xing Z, Lu C, Hu D, Yu YY, Wang X, Colnot C, et al. Multiple roles for CCR2 during fracture healing. Dis Model Mech. (2010) 3:451–8. doi: 10.1242/dmm003186

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Göktürk E, Turgut A, Bayçu C, Günal I, Seber S, Gülbas Z. Oxygen-free radicals impair fracture healing in rats. Acta Orthop Scand. (1995) 66:473–5. doi: 10.3109/17453679508995590

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Recknagel S, Bindl R, Kurz J, Wehner T, Ehrnthaller C, Knöferl MW, et al. Experimental blunt chest trauma impairs fracture healing in rats. J Orthop Res. (2011) 29:734–9. doi: 10.1002/jor21299

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Alexander KA, Chang MK, Maylin ER, Kohler T, Müller R, Wu AC, et al. Osteal macrophages promote in vivo intramembranous bone healing in a mouse tibial injury model. J Bone Miner Res. (2011) 26:1517–32. doi: 10.1002/jbmr354

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Nam D, Mau E, Wang Y, Wright D, Silkstone D, Whetstone H, et al. T-lymphocytes enable osteoblast maturation via IL-17F during the early phase of fracture repair. PLoS ONE. (2012) 7:e40044. doi: 10.1371/journal.pone0040044

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Toben D, Schroeder I, El Khassawna T, Mehta M, Hoffmann JE, Frisch JT, et al. Fracture healing is accelerated in the absence of the adaptive immune system. J Bone Miner Res. (2011) 26:113–24. doi: 10.1002/jbmr185

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Schmidt-Bleek K, Schell H, Schulz N, Hoff P, Perka C, Buttgereit F, et al. Inflammatory phase of bone healing initiates the regenerative healing cascade. Cell Tissue Res. (2012) 347:567–73. doi: 10.1007/s00441-011-1205-7

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Bernhardsson M, Dietrich-Zagonel F, Tätting L, Eliasson P, Aspenberg P. Depletion of cytotoxic (CD8+) T cells impairs implant fixation in rat cancellous bone. J Orthop Res. (2019) 37:805–11. doi: 10.1002/jor24246

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Weitzmann MN. Bone and the immune system. Toxicol Pathol. (2017) 45:911–24. doi: 10.1177/0192623317735316

CrossRef Full Text | Google Scholar

36. Yang S, Ding W, Feng D, Gong H, Zhu D, Chen B, et al. Loss of B cell regulatory function is associated with delayed healing in patients with tibia fracture. Apmis. (2015) 123:975–85. doi: 10.1111/apm12439

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Cottrell J, O'connor JP. Effect of non-steroidal anti-inflammatory drugs on bone healing. Pharmaceuticals (Basel). (2010) 3:1668–93. doi: 10.3390/ph3051668

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Loder RT. The influence of diabetes mellitus on the healing of closed fractures. Clin Orthop Relat Res. (1988) 232:210–6. doi: 10.1097/00003086-198807000-00028

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Kusumbe AP, Ramasamy SK, Adams RH. Coupling of angiogenesis and osteogenesis by a specific vessel subtype in bone. Nature. (2014) 507:323–8. doi: 10.1038/nature13145

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Ramasamy SK, Kusumbe AP, Wang L, Adams RH. Endothelial Notch activity promotes angiogenesis and osteogenesis in bone. Nature. (2014) 507:376–80. doi: 10.1038/nature13146

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Goerke SM, Obermeyer J, Plaha J, Stark GB, Finkenzeller G. Endothelial progenitor cells from peripheral blood support bone regeneration by provoking an angiogenic response. Microvasc Res. (2015) 98:40–7. doi: 10.1016/j.mvr.2014.12001

PubMed Abstract | CrossRef Full Text | Google Scholar