Analysis of the characteristics and expression profiles of coding and noncoding RNAs of human dental pulp stem cells in hypoxic conditions
Stem Cell Research & Therapy volume 10, Article number: 89 (2019)
Human dental pulp stem cell (DPSC)-mediated regenerative endodontics is a promising therapy for damaged teeth; however, the hypoxic environment in root canals can affect tissue regeneration. In this study, we investigate the characteristics and possible regulatory mechanisms of DPSC function under hypoxic conditions.
Human DPSCs were cultured under normoxia (20% O2) and hypoxia (3% O2). DPSC proliferation and osteo/odontogenic differentiation potential were assessed by Cell Counting Kit-8 (CCK8) assay, carboxyfluorescein succinimidyl ester (CFSE) assay, alkaline phosphatase (ALP) activity, Alizarin red staining, real-time RT-PCR assays, and western blot analysis. Microarray and bioinformatic analyses were performed to investigate the differences in the mRNA, lncRNA, and miRNA expression profiles of DPSCs.
DPSCs exhibited a more powerful proliferation ability and lower osteo/odontogenic differentiation potential in hypoxic conditions. A total of 60 mRNAs (25 upregulated and 35 downregulated), 47 lncRNAs (20 upregulated and 27 downregulated), and 14 miRNAs (7 upregulated and 7 downregulated) in DPSCs were differentially expressed in the hypoxia group compared with the normoxia group. Bioinformatic analysis identified that 7 mRNAs (GRPR, ERO1L, ANPEP, EPHX1, PGD, ANGPT1, and NQO1) and 5 lncRNAs (AF085958, AX750575, uc002czn.2, RP3-413H6.2, and six-twelve leukemia (STL)) may be associated with DPSCs during hypoxia according to CNC network analysis, while 28 mRNAs (including GYS1, PRKACB, and NQO1) and 13 miRNAs (including hsa-miR-3916 and hsa-miR-192-5p) may be involved according to miRNA target gene network analysis. The depletion of one candidate lncRNA, STL, inhibited the osteo/odontogenic differentiation potentials of DPSCs.
Our results revealed that hypoxia could enhance the proliferation ability and impair the osteo/odontogenic differentiation potential of DPSCs in vitro. Furthermore, our results identified candidate coding and noncoding RNAs that could be potential targets for improving DPSC function in regenerative endodontics and lead to a better understanding of the mechanisms of hypoxia’s effects on DPSCs.
Pulpitis and periapical periodontitis are the most prevalent oral diseases, and root canal treatment is the current primary treatment. Although the treatment has been shown to be highly effective, after treatment, teeth can become devitalized, brittle, and susceptible to fracture, leading to a higher incidence of extraction [1, 2]. During the past decade, advanced discoveries in cell-based therapies have provided novel insights into tooth regeneration. Regenerative endodontics (RE) is a clinical procedure that aims to regenerate the dentin-pulp complex (DPC) inside root canals . The DPC can provide nutrition and biological barriers for teeth, prevent the invasion of bacteria, reduce the risk of tooth fracture, and extend the life of a tooth. RE is a more ideal treatment for teeth that are damaged from pulpitis and periapical periodontitis, which is garnering increased attention [4,5,6]. Human dental pulp stem cells (DPSCs), a type of mesenchymal stem cell (MSC), are isolated from dental pulp and are ideal seeding cells for RE. A previous study reported that pulp-like tissues were formed by transplanting 3D DPSC constructs into human root canals after subcutaneous implantation into immunodeficient mice . Another study revealed that hDPSCs combined with NF-SMS were used to fill the pulp cavity and form pulp-like tissues similar to native pulp in a rat model . However, true regeneration of DPC has not been achieved in RE [9,10,11,12,13,14]. To improve the outcomes of RE, a key issue is to improve DPSC function in dental clinic conditions; therefore, it is necessary to explore the mechanisms controlling the directed differentiation of DPSCs in these conditions.
The microenvironmental niche of MSCs, which maintains and regulates their function, is the key element for determining the characteristics of MSCs. Oxygen in the microenvironment is necessary for the normal growth and differentiation of cells. Oxygen tensions in human organs and tissues vary from 2 to 9% in physiological conditions. In pathologic conditions or in some relatively hypoxic areas, oxygen tensions can drop to 1% or, in some cases, lower . Oxygen is required for ATP generation, which is necessary to provide free energy for cell survival. Cells maintain the supply of ATP via oxidative phosphorylation in mitochondria. The survival of cells is threatened with the loss of oxygen. Cells can sense changes in oxygen tension, initiate or inhibit gene transcription via the stimulation of hypoxia-inducible factors, and activate adaptive processes that can enhance cell survival in the case of limited oxygen . The effective isolation, culture, and expansion of MSCs under hypoxic conditions will mimic the physical microenvironmental niche and is crucial for maintaining and protecting implanted MSCs and beneficial for tissue regeneration in the physical hypoxic environment. Currently, MSCs are principally isolated and cultured under normoxia (20–21% O2) in vitro, which is markedly different from their typical in vivo environment. Previous studies have indicated that after being cultured in normoxic conditions, dental MSCs recultured in hypoxic conditions possess altered proliferation capacities and differentiation potentials; however, the results are inconclusive due to different cell sources and experimental conditions. Most studies have shown that hypoxia can promote the proliferation of dental stem cells [17,18,19,20,21,22,23,24]. However, some studies have reported that hypoxia had no or negative effects on the proliferation and survival of dental stem cells [25,26,27,28]. Some studies have shown that hypoxia can enhance the osteo/odontogenic differentiation of dental stem cells [17, 21, 23, 24, 26, 29,30,31], while other studies reported the opposite result . In addition, it has been reported that cell death can occur under severe hypoxia (0.01% O2) or anoxia . MSC-mediated regenerative endodontics is a promising method for damaged teeth; however, root canals also provide a hypoxic environment for MSCs . There is little known about the effects of hypoxia on dental stem cells and how to improve their function in a hypoxic environment. Thus, there is a need to discover DPSC characteristics and regulatory mechanisms in the hypoxic environment of root canals.
Here, 3% O2 tension was used to mimic hypoxia, and the characteristics and gene expression profiles of DPSCs in hypoxia were investigated. In addition, bioinformatic analysis was conducted to predict the interactions of coding and noncoding RNA and identify core regulatory factors of DPSCs in hypoxia. Furthermore, loss-of-function assays were performed to detect the effect of one candidate noncoding RNA, six-twelve leukemia (STL), on DPSCs. These results may help reveal the molecular mechanisms of DPSCs in a hypoxic environment, which may contribute to the identification of the key regulators that are necessary for DPSC-mediated dental tissue regeneration.
DPSC isolation and culture
Human teeth were collected with informed patient consent and based on the protocol approved by Beijing Stomatological Hospital, School of Stomatology, Capital Medical University, (Ethics Committee Agreement, Beijing, Stomatological Hospital Ethics Review No. 2011-02). Human dental pulps were obtained from healthy teeth extracted for impaction or other orthodontic reasons. The tissues from the crown and superior 2/3 root pulp were gently separated, and DPSCs were isolated, cultured, and characterized as previously described . The culture media was changed every 3 days. DPSCs from passages 3–5 were used in further experiments. For hypoxic culture, cells were grown in a hypoxic chamber flushed with 3% O2 and 5% CO2, with a balance of 92% N2 at 37 °C by using the hypoxic incubator (Thermo Scientific™ Heracell VIOS 160i, USA) which can accurately adjust the O2 concentration of 1–21%.
Cell Counting Kit-8 (CCK8) assay
The proliferation of DPSCs was determined using the CCK8 kit, according to the manufacturer’s instructions (Dojindo Laboratories, Kumamoto, Japan). A volume of 100 μL of DPSC suspension (3000 cells/well) was dispensed into the wells of a 96-well plate. After preincubation for 2 or 3 days in a humidified incubator (37 °C, 5% CO2), the culture medium was replaced with 100 μL MEM containing 10 μL CCK8, then incubated for 2 h at 37 °C. Absorbance (OD) at 450 nm was measured in a multiwell spectrophotometer.
Carboxyfluorescein succinimidyl ester (CFSE) assays were performed according to the CellTrace™ CFSE Cell Proliferation Kit protocol (Invitrogen, Carlsbad, CA, USA) for labeling cells in suspension as previously described . Briefly, DPSCs were stained, then seeded at a density of 1.0 × 105 cells/plate in six well plates. After 3 days of culture, DPSCs were harvested and analyzed using a flow cytometer (FACSCalibur, BD Biosciences, USA) with 635 nm excitation and emission filter appropriate for far red. The proliferation index was calculated by Modfit LT (Verity Software House, Topsham, ME, USA).
Alkaline phosphatase (ALP) activity assay and Alizarin red staining
To detect mineralization, DPSCs were grown in mineralization-inducing medium consisting of α-MEM, 15% FBS, 2 mM/L glutamine, 100 U/mL penicillin, 100 μg/mL streptomycin, 100 μM/L ascorbic acid, 10 mM/L β-glycerophosphate (Sigma-Aldrich, St. Louis, MO, USA), 1.8 mM/L of KH2PO4, and 10 nM/L dexamethasone (Sigma-Aldrich). Cells were induced for 5 days, and an ALP activity assay was performed as previously described . DPSCs were induced for 2 or 3 weeks, and Alizarin red staining and quantitative calcium analysis was performed as previously described .
Real-time RT-PCR was used to detect the expression levels of osteo/odontogenic marker genes and to verify the changes in gene expression that were detected by microarray analysis. Total RNA was isolated from DPSCs with TRIzol reagent (Invitrogen). cDNA was synthesized using reverse transcriptase with random primers for mRNA and lncRNA and with a special stem-loop primer for miRNA. The QuantiTect SYBR Green PCR kit (Qiagen, Hilden, Germany) was used in real-time PCR reactions, and the IcycleriQ multicolor real-time PCR detection system was used (Bio-Rad, Hercules, CA, USA). The gene expression levels were evaluated by normalizing the PCR signal to that of GAPDH (mRNA and lncRNA) or U6 (miRNA) and applying the 2−△△CT calculation. Each reaction was run in triplicate, and the entire procedure was repeated three times. The primers were designed by applying the online D-LUX DesignerTM program (Invitrogen) and are shown in Additional file 1: Table S1.
Western blot analysis
To detect proteins of osteo/odontogenic differentiation genes, DPSCs were cultured in mineralization-inducing medium for 2 weeks, and then western blot analysis was performed as previously described . Primary antibodies were obtained from commercial sources as follows: polyclonal antibody against bone sialoprotein (BSP) (Cat. No. bs-2668R, Bioss, Beijing, China); monoclonal antibody against osteocalcin (OCN) (Cat. No. ab133612, Abcam, Cambridge, UK); polyclonal antibody against dentin sialophosphoprotein (DSPP) (Cat. No. bs-10316R Bioss, Beijing, China); polyclonal antibody against beta-actin (β-actin) (Cat. No. AC026, Abclonal, Wuhan, China).
RNA preparation and microarray analysis
To detect differential coding and noncoding transcript expression profiles, DPSCs from three different individuals were cultured in normoxia or hypoxia for 3 days. Total RNA was then extracted from DPSCs using TRIzol and the RNeasy mini kit (Qiagen). The quality and quantity of RNA were verified by multiImager and spectrophotometer (Meriton, Beijing, China) analysis. The RNA labeling, hybridization, and scanning of the chips (Affymetrix HTA2.0 Array for mRNA and lncRNA, Affymetrix GeneChip miRNA 4.0 Array for miRNA) were performed as outlined in the Affymetrix technical manual by Cnkingbio Corp (Beijing, China). The threshold set for up- or downregulated mRNA, lncRNA, or miRNA was a fold change > 1.5 and a p value < 0.05.
The microarray raw data were normalized for follow-up analysis with Expression Console software (Affymetrix, CA, USA) using the MAS 5.0 statistical algorithm. The bioinformatic analysis was performed after normalized data was compared and filtered. Hierarchical clusters were performed by EPCLUST. Gene ontology (GO) analysis was used to analyze the genes’ primary functions. Pathway analysis was used to determine the significantly changed pathways of the differential genes based on KEGG, Biocarta, and Reatome. A coding-noncoding gene coexpression (CNC) network and miRNA target gene network were constructed to identify the interactions among genes and to locate core regulatory factors that played an important role in these networks.
STL shRNA (STLsh) and control shRNA (Consh) were purchased from Genepharma (Suzhou, China). The lentiviruses were transfected into DPSCs in the presence of polybrene (6 μg/mL, Sigma) for 12 h. After 48 h, the transfected DPSCs were selected with 1 μg/mL puromycin for 3 days. The target sequences for the shRNA were STLsh, 5′-CGGCATGACTAAGAGATATCG-3′ and Consh, 5′-TTCTCCGAACGTGTCACGT-3′.
Oil Red O staining
To detect adipogenic differentiation, DPSCs were grown in adipogenic-inducing medium for 3 weeks under normoxic conditions, and then Oil Red O staining was performed as previously described . The OD at 500 nm was measured using 100% isopropanol as a blank.
The statistical analyses were performed by using SPSS 16.0 software (SPSS Inc., Chicago, IL, USA). Significance was determined by the Student’s t test or one-way ANOVA analysis. For all tests, p < 0.05 was considered statistically significant.
Hypoxia enhanced the proliferation ability and inhibited the osteo/odontogenic differentiation potential of DPSCs
To investigate the effects of hypoxia on DPSCs, DPSCs were cultured in both normoxic and hypoxic conditions. CCK8 assays showed that the cell numbers of DPSCs were increased in hypoxic conditions compared with normoxic conditions at days 2 and 3 (Fig. 1a). CFSE assays also confirmed that the cell proliferation of DPSCs was increased in the hypoxia group compared with the normoxia group at day 3 (Fig. 1b, c). We then tested whether hypoxia could affect the osteo/odontogenic differentiation abilities of DPSCs. DPSCs were cultured in mineralization-inducing medium in both normoxic and hypoxic conditions. ALP activity, an early indicator for osteo/odontogenic differentiation in DPSCs, was determined at 5 days postinduction. The results showed that hypoxia resulted in a significant decrease in the ALP activity of DPSCs (Fig. 2a). After 2 weeks of induction, Alizarin red staining and calcium quantitation revealed weakened mineralization of DPSCs in hypoxic conditions compared with the control group (Fig. 2b, c). Real-time RT-PCR showed that the osteo/odontogenic marker genes BSP, OCN, and DSPP were downregulated at 7 and 14 days postinduction in the hypoxia group compared with the normoxia group (Fig. 2d–f). Western blot results also showed that the expressions of BSP, OCN, and DSPP in DPSCs were downregulated at 2 weeks postinduction in the hypoxia group (Fig. 2g).
The profiling of differentially expressed mRNAs, lncRNAs, and miRNAs in DPSCs under hypoxic and normoxic conditions
To illuminate why hypoxia could affect the proliferation capacity and differentiation potential of DPSCs, microarray and bioinformatic analyses were performed to identify and locate core regulators of DPSCs in hypoxia. The RVM t test was applied to filter the genes that were differentially expressed, and the differentially expressed genes with 1.5-fold changes and p values < 0.05 were selected for subsequent analysis. From the microarray data, a total of 60 mRNAs with differential expression in DPSCs were identified between the hypoxia and normoxia groups, of which 25 were upregulated and 35 were downregulated in the hypoxia group compared with those in the normoxia group (Additional file 2: Table S2). A total of 47 lncRNAs in DPSCs were significantly differentially expressed between the hypoxia and normoxia groups, of which 20 were upregulated and 27 were downregulated in the hypoxia group compared with those in the normoxia group (Additional file 3: Table S3). In addition, a total of 14 miRNAs in DPSCs were differentially expressed between the hypoxia and normoxia groups, of which seven were upregulated and seven were downregulated in the hypoxia group compared with the normoxia group (Additional file 4: Table S4). To confirm the reliability of the microarray data, seven differentially expressed mRNAs (GRPR, CA12, GFRA2, ERO1L, EPAS1, TXNRD1, and NQO1), two differentially expressed lncRNAs (LINC00707 and STL), and two differentially expressed miRNAs (hsa-miR-3916 and hsa-miR-6744-5p) were selected, and real-time RT-PCR was performed to detect their expression levels. The results showed that the expression levels of the selected coding and noncoding RNAs were consistent with the microarray results and confirmed the reliability of the microarray data (Fig. 3).
GO analysis was performed to determine the significantly changed functions obtained from the differentially expressed mRNAs of DPSCs in hypoxic conditions compared with normoxic conditions. A total of 91 GO functions (p < 0.01, FDR < 0.05) were obtained, of which 28 were upregulated and 63 were downregulated in the hypoxia group. The negative logarithm of the p value (− LgP) was applied to show the correlation of gene expression levels and the relevant biological processes (Additional file 5: Figure S1 and Additional file 6: Figure S2). The GO functions enriched from the upregulated genes were related to glycogen biosynthetic process, positive regulation of DNA metabolic process, wound healing, negative regulation of mesoderm development, negative regulation of mitotic sister chromatid separation, cellular protein modification process, and several others. The GO functions enriched from downregulated genes included mesoderm formation, positive regulation of glutamate-cysteine ligase activity, response to hypoxia, oxidation-reduction process, pentose biosynthetic process, regulation of cell growth, and several others.
Pathway analysis using the KEGG database was then performed to enrich for the major pathways that were significantly associated with the differentially expressed mRNAs of DPSCs in hypoxia compared with the normoxia group. A total of 17 pathways (p < 0.05) were identified, among which 5 were upregulated and 12 were downregulated in the hypoxia group (Additional file 7: Figure S3 and Additional file 8: Figure S4). The upregulated pathways were starch and sucrose metabolism, insulin resistance, neurotrophin signaling pathway, oocyte meiosis, and nitrogen metabolism. The downregulated pathways included Wnt signaling pathway, pentose phosphate pathway, glutathione metabolism, ubiquinone and other terpenoid-quinone biosynthesis, hypoxia-inducible factor-1 (HIF-1) signaling pathway, carbon metabolism, and several others.
Next, a CNC network was constructed to identify the interactions and significance degrees among the differentially expressed mRNAs and lncRNAs based on clustering coefficients and degrees (Fig. 4). Twelve core regulatory genes were identified including seven mRNAs (GRPR, ERO1L, ANPEP, EPHX1, PGD, ANGPT1, and NQO1) and five lncRNAs (AF085958, AX750575, uc002czn.2, RP3-413H6.2, and STL). In addition, a miRNA target gene network was constructed to identify the interactions and significance degrees among the differentially expressed miRNAs and mRNAs (Fig. 5). A total of 41 core regulatory genes were identified including 28 mRNAs (GYS1, PRKACB, NQO1, and several others) and 13 miRNAs (hsa-miR-3916, hsa-miR-192-5p, and several others).
Depletion of STL inhibited the osteo/odontogenic differentiation potential of DPSCs
Then, we selected one candidate lncRNA, STL, which was identified as a possible core regulatory factor of DPSCs according to the CNC network, to investigate its function in DPSCs. The real-time RT-PCR results showed that the expression of STL increased with the osteo/odontogenic differentiation of DPSCs (Fig. 6a). To investigate the possible regulatory function of STL on DPSCs, STL in DPSCs was knocked down through STL lentiviral shRNA infection. Approximately 70% knockdown efficiency was achieved as determined by real-time RT-PCR (Fig. 6b). To test whether STL could affect the differentiation potential of DPSCs, DPSCs were cultured in osteogenic-inducing medium under normoxic conditions. ALP activity was determined 5 days postinduction, and the results showed a significant decrease in ALP activity after knockdown of STL compared with the Consh and Mock groups (Fig. 6c). After 2 weeks of induction, Alizarin red staining and calcium quantitation revealed weakened mineralization in vitro in STL-depleted DPSCs compared to the Consh and Mock groups (Fig. 6d, e). Western blot results showed that the expressions of BSP and DSPP were downregulated at 2 weeks postinduction in STL-depleted DPSCs compared to the Consh and Mock groups (Fig. 6f). Real-time RT-PCR showed that NQO1 was downregulated and ERO1L was upregulated in STL-depleted DPSCs compared to the control group (Fig. 6g, h). Interestingly, we also found a significant decrease of lipid deposits in STL-depleted DPSCs compared to the Consh and Mock groups at 3 weeks after adipogenic induction (Additional file 9: Figure S5); these indicated that depletion of STL decreased the adipogenic differentiation in DPSCs.
In tissue engineering, MSCs need to be implanted into the tissue of the host after in vitro amplification. During this process, changes in oxygen tension in the extracellular environment can affect the growth and differentiation of MSCs. In the present study, we found that DPSCs exhibited significantly higher proliferation under 3% oxygen. Our results are consistent with those of previous studies, which demonstrated that culturing under appropriate hypoxic conditions or pretreatment with 1–5% oxygen could enhance the proliferation potential of DPSCs [17,18,19,20]. In the present study, we also found that hypoxia reduced ALP activity, weakened the mineralization ability of DPSCs, and decreased the expression of BSP, OCN, and DSPP, indicating that the osteo/odontogenic differentiation of DPSCs was inhibited in 3% oxygen. Our results are in agreement with those of a previous study that found ALP activity and the expression of DMP1, DSPP, and OCN were suppressed when DPCs were cultured under hypoxic conditions . However, other studies showed that hypoxia could promote the mineralization of DPSCs and upregulate mineral-associated genes [17, 21, 29]. These different results may be caused by different oxygen tensions and action times. Our results confirm that oxygen tension is an important regulator and suggest that 3% oxygen promotes cell proliferation and inhibits DPSC differentiation.
For successful cell therapy, the molecular mechanisms that underlie the effects of hypoxia on human DPSCs need to be elucidated. Therefore, the expression profiles of mRNAs, lncRNAs, and miRNAs of DPSCs (in hypoxia vs normoxia) were determined, and the reliability of the microarray data was confirmed by real-time RT-PCR. In total, 60 differentially expressed mRNAs, 47 differentially expressed lncRNAs, and 14 differentially expressed miRNAs were identified. These differentially expressed coding and noncoding RNAs may be candidates for key regulators controlling the proliferation and differentiation of DPSCs in hypoxia and warrant further bioinformatic analysis.
By GO and pathway analysis, we screened out the significantly altered functions and pathways associated with DPSCs in hypoxia, including the Wnt and HIF-1 signaling pathways. The canonical Wnt/β-catenin signaling pathway has different functions in hard tissue development and regeneration, depending on cell type and stage of development. Previous reports found that the proliferation and differentiation potentials of DPSCs could be regulated via the Wnt/β-catenin signaling pathway [35, 36]. HIF-1, which consists of two subunits, HIF-1α and HIF-1β, is known as a key nuclear transcription factor related to hypoxia. HIF-1α is regulated by the cellular oxygen concentration and determines the transcriptional activity of HIF-1. HIF-1α can play an important role in regulating the proliferation, differentiation, and pluripotency of stem cells through signaling pathways such as Wnt, Notch, and Sonic hedgehog (Shh) [25, 37]. It was reported that HIF-1α protects MSCs via the induction of autophagy and the PI3K/AKT/mTOR pathway . A study in porcine dental pulp-derived cells under severe hypoxic conditions reported that the expression of BSP, DMP-1, and DSPP was repressed, while the expression of HIF-1α was increased at 6 and 12 h then downregulated after 18 h . From these data, we speculate that the Wnt and HIF-1 signaling pathways could contribute to DPSC characteristics in hypoxic conditions.
According to CNC network and miRNA target gene network analysis, some of the differentially expressed RNAs in the two groups were screened out as possible regulatory factors of DPSCs in hypoxia, including ERO1L, NQO1, STL, hsa-miR-192-5p, and several others. NQO1 can function as a direct superoxide anion radical scavenger or a major target of the Nrf2 pathway to maintain cellular redox homeostasis and attenuate cellular oxidative stress. In addition to oxidation-reduction process, NQO1 is related to obesity, hypertension, renal injury, bone metabolism, tumor growth, and other processes [39,40,41]. It was reported that NQO1 and HO-1 were significantly downregulated in human primary endothelial cells that had been adapted to physiological O2 (5%) according to microarray analysis and further validation, which is consistent with our results . NQO1 may act downstream of the Nrf2 pathway, and it was revealed that the mineralization of human periodontal ligament fibroblasts was improved by downregulating Nrf2 through activation of the Wnt/β-catenin pathway . ERO1L belongs to the family of classic oxygen-regulated genes, and its expression is augmented under hypoxia, which is consistent with our results . It was revealed that the expression of ERO1L was regulated via the HIF pathway during hypoxia [44, 45]. ERO1L controls several proteins via reoxidation of “client” protein disulfide isomerase. Furthermore, it can regulate cell proliferation and plays crucial roles in the progression of various tumors, cardiovascular diseases, diabetes, and other diseases [46,47,48,49,50,51]. As a miRNA, miR-192-5p plays an important role in cell proliferation and differentiation [52,53,54]. Emerging evidence has further suggested its role in the development of tumors, cardiovascular diseases, liver disease, and other diseases [55,56,57]. Consistent with our results, it has been reported that the release of miR-192-5p may be attributed to oxidative stress, and its expression was regulated by hypoxia [58, 59]. As an lncRNA, STL does not contain any long open reading frames. Emerging evidence revealed that lncRNAs play important roles in the proliferation and osteogenic differentiation capacity of MSCs through different regulatory mechanisms, including transcription factor binding and chromatin modification [60,61,62,63,64]. Some lncRNAs have been identified as key regulators of the proliferation capacity and differentiation potential of DPSCs [33, 65]. It has also been reported that STL was related to hematologic malignancies and uterine tumors resembling ovarian sex cord tumors [66,67,68]. According to the CNC network, STL showed positive and negative regulatory effects on NQO1 and ERO1L, respectively. Due to the effects of NQO1 and ERO1L on the functions of stem cells, we speculate that STL may be a key regulator of these processes.
Indeed, we found that the expression of STL in DPSCs was enhanced along with osteo/odontogenic differentiation of DPSCs, indicating that STL may play a positive role in this process. We therefore conducted loss-of-function assays to investigate the functions of STL in DPSCs. Our results showed that the knockdown of STL inhibited the osteo/dentinogenic differentiation potential of DPSCs, indicating that STL was involved in the positive regulation of osteo/dentinogenic differentiation ability of DPSCs. Consistent with the CNC network, NQO1 was downregulated while ERO1L was upregulated when STL was knocked down in DPSCs. In addition, STL and NQO1 were downregulated and ERO1L was upregulated in hypoxia, indicating that decreased STL might be one factor for the impaired osteo/dentinogenic differentiation ability of DPSCs in hypoxia and that this process might be downregulated via NQO1 and upregulated via ERO1L. It is possible that STL works on DPSCs via interactions with these two genes, and therefore, STL may be considered as a possible therapeutic target for dental tissue regeneration. However, further studies are needed to verify this hypothesis and the underlying mechanism.
In conclusion, our results revealed that 3% oxygen tension could enhance the proliferative ability and impair the osteo/odontogenic differentiation potential of DPSCs, indicating that hypoxia is an important factor affecting DPSC-mediated dental tissue regeneration. Furthermore, our results identified that candidate coding and noncoding RNA could be potential targets for improving the DPSC characteristics in hypoxia and lead to a better understanding of the mechanisms of hypoxia’s effects on DPSCs. In addition, our results highlighted the significant involvement of one lncRNA, STL, in the positive regulation of the osteo/odontogenic differentiation ability of DPSCs and indicated that STL could be a potential target in regenerative endodontics.
Carboxyfluorescein succinimidyl ester
Coding-noncoding gene coexpression
Dental pulp stem cells
Endoplasmic reticulum oxidoreductase 1-like
Kyoto Encyclopedia of Genes and Genomes
Long noncoding RNA
Mesenchymal stem cells
NAD(P)H dehydrogenase quinone 1
Reverse transcription-polymerase chain reaction
Wingless-type MMTV integration site family
Ng YL, Mann V, Gulabivala K. Tooth survival following non-surgical root canal treatment: a systematic review of the literature. Int Endod J. 2010;43:171–89.
Moule AJ, Kahler B. Diagnosis and management of teeth with vertical root fractures. Aust Dent J. 1999;44:75–87.
Gangolli RA, Devlin SM, Gerstenhaber JA, et al. A bilayered poly (lactic-co-glycolic acid) scaffold provides differential cues for the differentiation of dental pulp stem cells. Tissue Eng Part A. 2018; doi: 10.1089 /ten.TEA.2018.0041. [Epub ahead of print]
Iwaya SI, Ikawa M, Kubota M. Revascularization of an immature permanent tooth with apical periodontitis and sinus tract. Dent Traumatol. 2001;17:185–7.
Ding RY, Cheung GSP, Chen J, et al. Pulp revascularization of immature teeth with apical periodontitis: a clinical study. J Endod. 2009;35:745–9.
Sun HH, Jin T, Yu Q, et al. Bological approaches toward dental pulp regeneration by tissue engineering. J Tissue Eng Regen Med. 2011;5:e1–16.
Itoh Y, Sasaki JI, Hashimoto M, et al. Pulp regeneration by 3-dimensional dental pulp stem cell constructs. J Dent Res. 2018;97:1137–43.
Kuang R, Zhang Z, Jin X, et al. Nanofibrous spongy microspheres for the delivery of hypoxia-primed human dental pulp stem cells to regenerate vascularized dental pulp. Acta Biomater. 2016;33:225–34.
Yamauchi N, Yamauchi S, Nagaoka H, et al. Tissue engineering strategies for immature teeth with apical periodontitis. J Endod. 2011;37:390–7.
Tawfik H, Abu-Seida AM, Hashem AA, et al. Regenerative potential following revascularization of immature permanent teeth with necrotic pulps. Int Endod. 2013;46:910–22.
Wang XJ, Thibodeau B, Trope M, et al. Histologic characterization of regenerated tissues in canal space after the revitalization/revascularization procedure of immature dog teeth with apical periodontitis. J Endod. 2010;36:56–63.
Becerra P, Ricucci D, Loghin S, et al. Histologic study of a human immature permanent premolar with chronic apical abscess after revascularization/revitalization. J Endod. 2014;40:133–9.
Mangione F, EzEldeen M, Bardet C, et al. Implanted dental pulp cells fail to induce regeneration in partial pulpotomies. J Dent Res. 2017;96:1406–13.
Altaii M, Richards L, Rossi-Fedele G. Histological assessment of regenerative endodontic treatment in animal studies with different scaffolds: a systematic review. Dent Traumatol. 2017;33:235–44.
Mohyeldin A, Garzón-Muvdi T, Quiñones-Hinojosa A. Oxygen in stem cell biology: a critical component of the stem cell niche. Cell Stem Cell. 2010;7:150–61.
Guzy RD, Schumacker PT. Oxygen sensing by mitochondria at complex III:the paradox of increased reactive oxygen species during hypoxia. Exp Physiol. 2016;91:807–19.
Li L, Zhu YQ, Jiang L, et al. Hypoxia promotes mineralization of human dental pulp cells. J Endod. 2011;37:799–802.
Iida K, Takeda-Kawaguchi T, Tezuka Y, et al. Hypoxia enhances colony formation and proliferation but inhibits differentiation of human dental pulp cells. Arch Oral Biol. 2010;55:648–54.
Gong Q, Quan J, Jiang H, et al. Regulation of the stromal cell–derived factor-1α–CXCR4 axis in human dental pulp cells. J Endod. 2010;36(9):1499–503.
Sakdee JB, White RR, Pagonis TC, et al. Hypoxia-amplified proliferation of human dental pulp cells. J Endod. 2009;35:818–23.
Amemiya K, Kaneko Y, Muramatsu T, et al. Pulp cell responses during hypoxia and reoxygenation in vitro. Eur J Oral Sci. 2003;111:332–8.
Kanafi MM, Ramesh A, Gupta PK, et al. Influence of hypoxia, high glucose, and low serum on the growth kinetics of mesenchymal stem cells from deciduous and permanent teeth. Cells Tissues Organs. 2013;198:198–208.
Zhang QB, Zhang ZQ, Fang SL, et al. Effects of hypoxia on proliferation and osteogenic differentiation of periodontal ligament stem cells: an in vitro and in vivo study. Genet Mol Res. 2014;13:10204–14.
Dai Y, He H, Wise GE, et al. Hypoxia promotes growth of stem cells in dental follicle cell populations. J Biomed Sci Eng. 2011;4:454–61.
Zhou Y, Fan W, Xiao Y. The effect of hypoxia on the stemness and differentiation capacity of PDLC and DPC. Biomed Res Int. 2014;2014:890675.
Vanacker J, Viswanath A, De Berdt P, et al. Hypoxia modulates the differentiation potential of stem cells of the apical papilla. J Endod. 2014;40:1410–8.
Agata H, Kagami H, Watanabe N, et al. Effect of ischemic culture conditions on the survival and differentiation of porcine dental pulp-derived cells. Differentiation. 2008;76:981–93.
Wang J, Wei X, Ling J, et al. Side population increase after simulated transient ischemia in human dental pulp cell. J Endod. 2010;36:453–8.
Ito K, Matsuoka K, Matsuzaka K, et al. Hypoxic condition promotes differentiation and mineralization of dental pulp cells in vivo. Int Endod J. 2015;48:115–23.
Choi H, Jin H, Kim JY, et al. Hypoxia promotes CEMP1 expression and induces cementoblastic differentiation of human dental stem cells in an HIF-1-dependent manner. Tissue Eng Part A. 2014;20:410–23.
Wu Y, Yang Y, Yang P, et al. The osteogenic differentiation of PDLSCs is mediated through MEK/ERK and p38 MAPK signalling under hypoxia. Arch Oral Biol. 2013;58:1357–68.
Yu CY, Boyd NM, Cringle SJ, et al. Oxygen distribution and consumption in rat lower incisor pulp. Arch Oral Biol. 2002;47:529–36.
Yi Q, Liu O, Yan F, et al. Analysis of senescence-related differentiation potentials and gene expression profiles in human dental pulp stem cells. Cells Tissues Organs. 2017;203:1–11.
Li WZ, Lin X, Yang HQ, et al. Depletion of HOXA5 inhibits the osteogenic differentiation and proliferation potential of stem cells from the apical papilla. Cell Biol Int. 2018;42:45–52.
Rahman SU, Oh JH, Cho YD, et al. Fibrous topography-potentiated canonical Wnt signaling directs the odontoblastic differentiation of dental pulp-derived stem cells. ACS Appl Mater Interfaces. 2018;10:17526–41.
Liu N, Zhou M, Zhang Q, et al. Stiffness regulates the proliferation and osteogenic/odontogenic differentiation of human dental pulp stem cells via the WNT signalling pathway. Cell Prolif. 2018;51:e12435.
Mazumdar J, Dondeti V, Simon MC. Hypoxia-inducible factors in stem cells and cancer. J Cell Mol Med. 2009;13:4319–28.
Lv B, Hua T, Li F, et al. Hypoxia-inducible factor 1 α protects mesenchymal stem cells against oxygen-glucose deprivation-induced injury via autophagy induction and PI3K/AKT/mTOR signaling pathway. Am J Transl Res. 2017;9:2492–9.
Gu DR, Lee JN, Oh GS, et al. The inhibitory effect of beta-lapachone on RANKL-induced osteoclastogenesis. Biochem Biophys Res Commun. 2017;482:1073–9.
Al-Naqeb G, Rousová J, Kubátová A, et al. Pulicaria jaubertii E. Gamal-Eldin reduces triacylglyceride content and modifies cellular antioxidant pathways in 3T3-L1 adipocytes. Chem Biol Interact. 2016;253:48–59.
Tong YH, Zhang B, Yan YY, et al. Dual-negative expression of Nrf2 and NQO1 predicts superior outcomes in patients with non-small cell lung cancer. Oncotarget. 2017;8:45750–8.
Chapple SJ, Keeley TP, Mastronicola D, et al. Bach1 differentially regulates distinct Nrf2-dependent genes in human venous and coronary artery endothelial cells adapted to physiological oxygen levels. Free Radic Biol Med. 2016;92:152–62.
Kook SH, Lee D, Cho ES, et al. Activation of canonical Wnt/β-catenin signaling inhibits H2O2-induced decreases in proliferationand differentiation of human periodontal ligament fibroblasts. Mol Cell Biochem. 2016;411:83–94.
Gess B, Hofbauer KH, Wenger RH, et al. The cellular oxygen tension regulates expression of the endoplasmic oxidoreductase ERO1-Lalpha. Eur J Biochem. 2003;270:2228–35.
May D, Itin A, Gal O, et al. Ero1-L alpha plays a key role in a HIF-1-mediated pathway to improve disulfide bond formation and VEGF secretion under hypoxia: implication for cancer. Oncogene. 2005;24:1011–20.
Kutomi G, Tamura Y, Tanaka T, et al. Human endoplasmic reticulum oxidoreductin 1-α is a novel predictor for poor prognosis of breast cancer. Cancer Sci. 2013;104:1091–6.
Seol SY, Kim C, Lim JY, et al. Overexpression of endoplasmic reticulum oxidoreductin 1-α (ERO1L) is associated with poor prognosis of gastric cancer. Cancer Res Treat. 2016;48:1196–209.
Fujii M, Yoneda A, Takei N, et al. Endoplasmic reticulum oxidase 1α is critical for collagen secretion from and membrane type 1-matrix metalloproteinase levels in hepatic stellate cells. J Biol Chem. 2017;292:15649–60.
Han F, Xu Q, Zhao J, et al. ERO1L promotes pancreatic cancer cell progression through activating the Wnt/catenin pathway. J Cell Biochem. 2018;119:8996–9005.
Chin KT, Kang G, Qu J, et al. The sarcoplasmic reticulum luminal thiol oxidase ERO1 regulates cardiomyocyte excitation-coupled calcium release and response to hemodynamic load. FASEB J. 2011;25:2583–91.
Awazawa M, Futami T, Sakada M, et al. Deregulation of pancreas-specific oxidoreductin ERO1β in the pathogenesis of diabetes mellitus. Mol Cell Biol. 2014;34:1290–9.
Yan-Chun L, Hong-Mei Y, Zhi-Hong C, et al. MicroRNA-192-5p promote the proliferation and metastasis of hepatocellular carcinoma cell by targeting SEMA3A. Appl Immunohistochem Mol Morphol. 2017;25:251–60.
Jin H, Qiao F, Wang Y, et al. Curcumin inhibits cell proliferation and induces apoptosis of human non-small cell lung cancer cells through the upregulation of miR-192-5p and suppression of PI3K/Akt signaling pathway. Oncol Rep. 2015;34:2782–9.
Raut A, Khanna A. Enhanced expression of hepatocyte-specific microRNAs in valproic acid mediated hepatic trans-differentiation of human umbilical cord derived mesenchymal stem cells. Exp Cell Res. 2016;343:237–47.
Zhou S, Xiong M, Dai G, et al. MicroRNA-192-5p suppresses the initiation and progression of osteosarcoma by targeting USP1. Oncol Lett. 2018;15:6947–56.
Wu T, Chen Y, Du Y, et al. Serum exosomal MiR-92b-5p as a potential biomarker for acute heart failure caused by dilated cardiomyopathy. Cell Physiol Biochem. 2018;46:1939–50.
Liu XL, Cao HX, Wang BC, et al. miR-192-5p regulates lipid synthesis in non-alcoholic fatty liver disease through SCD-1. World J Gastroenterol. 2017;23:8140–51.
Matsuzaki J, Ochiya T. Extracellular microRNAs and oxidative stress in liver injury: a systematic mini review. J Clin Biochem Nutr. 2018;63:6–11.
Zhi F, Shao N, Xue L, et al. Characteristic microRNA expression induced by δ-opioid receptor activation in the rat liver under prolonged hypoxia. Cell Physiol Biochem. 2017;44:2296–309.
Yang Q, Jia L, Li X, et al. Long noncoding RNAs: new players in the osteogenic differentiation of bone marrow- and adipose-derived mesenchymal stem cells. Stem Cell Rev. 2018;14:297–308.
Zhang J, Tao Z, Wang Y. Long non-coding 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.
Li Z, Jin C, Chen S, et al. Long non-coding RNA MEG3 inhibits adipogenesis and promotes osteogenesis of human adipose-derived mesenchymal stem cells via miR-140-5p. Mol Cell Biochem. 2017;433:51–60.
Zhang W, Dong R, Diao S, et al. Differential long noncoding RNA/mRNA expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells. Stem Cell Res Ther. 2017;8:30.
Zhuang W, Ge X, Yang S, et al. Upregulation of lncRNA MEG3 promotes qsteogenic differentiation of mesenchymal stem cells from multiple myeloma patients by targeting BMP4 transcription. Stem Cells. 2015;33:1985–97.
Zeng L, Sun S, Han D, et al. Long non-coding RNA H19/SAHH axis epigenetically regulates odontogenic differentiation of human dental pulp stem cells. Cell Signal. 2018;52:65–73.
Wang J, Blakey GL, Zhang L, et al. Uterine tumor resembling ovarian sex cord tumor: report of a case with t(X;6)(p22.3;q23.1) and t (4;18)(q21.1;q21.3). Diagn Mol Pathol. 2003;12:174–80.
Suto Y, Sato Y, Smith SD, et al. A t(6,12) (q23;p13) results in the fusion of ETV6 to a novel gene, STL, in a B-cell ALL cell line. Genes Chromosomes Cancer. 1997;18:254–68.
Fontanari Krause LM, Japp AS, Krause A, et al. Identification and characterization of OSTL (RNF217) encoding a RING-IBG-RING protein adjacent to a translocation breakpoint involving ETV6 in childhood ALL. Sci Rep. 2014;4:6565.
This work was supported by grants from the National Natural Science Foundation of China (81625005 to Z.P.F., 81170952 to B.X.H.), and the Discipline Construction Fund from the Beijing Stomatological Hospital, School of Stomatology, Capital Medical University (18-09-11 to R.T.S.).
Availability of data and materials
The microarray datasets supporting the conclusions of this article are available. All other relevant datasets have been uploaded as part of additional files.
Ethics approval and consent to participate
All human stem cells involved in the present research followed the International Society for Stem Cell Research “Guidelines for the Conduct of Human Embryonic Stem Cell Research.” The study was approved by the Ethical Committee of Beijing Stomatological Hospital (Review No. 2011–02) before the study started. Written informed consent was obtained from each patient, under the approved guidelines set by Beijing Stomatological Hospital, School of Stomatology, Capital Medical University.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Primer sequences used in real-time RT-PCR analysis. (DOCX 16 kb)
Table S2. The differentially expressed mRNAs of hDPSCs in hypoxic and normoxic conditions. (DOCX 18 kb)
Table S3. The differentially expressed lncRNAs of hDPSCs in hypoxic and normoxic conditions. (DOCX 20 kb)
Table S4. The differentially expressed miRNAs of hDPSCs in hypoxic and normoxic conditions. (DOCX 14 kb)
Figure S1. Significant GOs of upregulated genes of DPSCs in hypoxia. The x axis shows − LgP and the y axis shows GO category. A larger – LgP indicates a smaller p value for the difference. (TIFF 1845 kb)
Figure S2. Significant GOs of downregulated genes of DPSCs in hypoxia. The x axis shows − LgP and the y axis shows GO category. A larger − LgP indicates a smaller p value for the difference. (TIFF 2370 kb)
Figure S3. Significant pathways of upregulated genes of DPSCs in hypoxia. The x axis shows − LgP and the y axis shows pathways category. A larger – LgP indicates a smaller p value for the difference. (TIFF 512 kb)
Figure S4. Significant pathways of downregulated genes of DPSCs in hypoxia. The x axis shows − LgP and the y axis shows pathways category. A larger – LgP indicates a smaller p value for the difference. (TIFF 1140 kb)
Figure S5. Knockdown of STL inhibited the adipogenic differentiation of DPSCs in vitro. (a, b) Oil red staining (a) and a quantitation analysis (b) revealed a significant decrease in lipid deposits in STL-depleted DPSCs compared to the control groups 3 weeks after adipogenic induction. One-way ANOVA was performed to determine statistical significance. All error bars represent SD (n = 3). NS; no significant difference, **p < 0.01. Scale bar 200 μm. (TIFF 1495 kb)
About this article
Cite this article
Shi, R., Yang, H., Lin, X. et al. Analysis of the characteristics and expression profiles of coding and noncoding RNAs of human dental pulp stem cells in hypoxic conditions. Stem Cell Res Ther 10, 89 (2019). https://doi.org/10.1186/s13287-019-1192-2
- Dental pulp stem cells (DPSCs)
- Coding RNA
- Long noncoding RNA (lncRNA)
- MicroRNA (miRNA)