- Research article
- Open Access
The enzymatic hydrolysates from deer sinew promote MC3T3-E1 cell proliferation and extracellular matrix synthesis by regulating multiple functional genes
BMC Complementary Medicine and Therapies volume 21, Article number: 59 (2021)
Deer Sinew serves as a medicinal food, and has been used for treating skeletal diseases, especially bone diseases in a long history. Thus, it could become an alternative option for the prevention and therapeutic remedy of bone-related diseases. In our previous study, we established an optimal extraction process of the enzymatic hydrolysates from Chinese Sika deer sinews (DSEH), and we demonstrated that DSEH significantly promoted the proliferation of MC3T3-E1 cells (an osteoblast-like cell line) with a certain dose-effect relationship. However, the precise molecular mechanism of deer sinew in regulating bone strength is still largely unknown. The aim of this study was to explore the underlying molecular mechanism of DSEH on MC3T3-E1 cells proliferation and extracellular matrix synthesis.
Preparation and quality control were performed as previously described. The effect of DSEH at different administrated concentrations on cell proliferation was measured using both CCK-8 and MTT assays, and the capacity of DSEH on extracellular matrix synthesis was detected by Alizarin red staining and quantification. The gene expression pattern change of MC3T3-E1 cells under the treatment of DSEH was investigated by RNA-seq analysis accompanied with validation methods.
We demonstrated that DSEH promoted MC3T3-E1 cell proliferation and extracellular matrix synthesis by regulating multiple functional genes. DSEH significantly increased the expression levels of genes that promoted cell proliferation such as Gstp1, Timp1, Serpine1, Cyr61, Crlf1, Thbs1, Ctgf, P4ha2, Sod3 and Nqo1. However, DSEH significantly decreased the expression levels of genes that inhibited cell proliferation such as Mt1, Cdc20, Gas1, Nrp2, Cmtm3, Dlk2, Sema3a, Rbm25 and Hspb6. Furthermore, DSEH mildly increased the expression levels of osteoblast gene markers.
Our findings suggest that DSEH facilitate MC3T3-E1 cell proliferation and extracellular matrix synthesis to consolidate bone formation and stability, but prevent MC3T3-E1 cells from oxidative stress-induced damage, apoptosis and further differentiation. These findings deepened the current understanding of DSEH on regulating bone development, and provided theoretical support for the discovery of optional prevention and treatment for bone-related diseases.
Deer Sinew, a precious animal production, has been widely used in traditional Chinese remedies to support bone health for over 20 centuries. According to the theory of traditional Chinese medicine, Deer sinews are considered to be more effective than cattle sinews in nourishing kidney yang (shen yang) and strengthening bones and tendons . Several studies have shown that deer sinew extract could prevent bone loss in osteoporosis model and reduce the injury risk of musculoskeletal system. By using an ovariectomized rat model of osteoporosis, the researchers demonstrated that deer sinew extract could significantly increase bone mineral density and serum hydroxyproline, and improve the histomorphometric parameters and mechanical indicators of bone . Consistent with the above results, according to the investigation of a retinoic acid-induced rat model of osteoporosis, the researchers demonstrated that deer sinew extract could increase bone mineral density and bone weight . In a clinical investigation of the effect of deer sinew supplementation on the exercise performance and risk training injury of a population of athletes, the researchers demonstrated that deer sinew extract could improve exercise performance and reduce the injury risk of musculoskeletal system . Furthermore, several studies have also shown that deer sinew extract has potential anti-inflammatory effects [5,6,7]. Sinews, also known as tendons, are fibrous connective tissues that localize between bones and muscles . The major components of sinews are collagens, especially type I collagen, which account for 65 to 80% of the dry mass . Meanwhile, collagens including types I, III and V are the major extracellular matrices in bone, and type I collagen accounts for about 95% of the total collagen components present in bone and approximately 80% of the entire proteins of bone. Those collagens thus play pivotal roles in maintaining bone strength and further prevent bone fragility related diseases, such as osteogenesis imperfecta and osteoporosis . However, the precise molecular mechanism of deer sinew in regulating bone strength is still largely unknown.
In recent years, high-throughput sequencing technology has led rapid increase in understanding the physiological and pathological mRNA dynamics, which further influence the protein expression that involved in the growth, development and related functions of cells as well as organ systems . It is particularly noteworthy that RNA sequencing (RNA-seq) has also been widely used in skeletal biology research to solve tough problems regarding the development and disease control of skeletal tissues, including bone, cartilage, tendon and ligament . During bone formation and repair, osteoblasts serve as bone-forming cells that are required for the synthesis of bone matrix, including collagen proteins (mainly type I collagen), noncollagen proteins (osteocalcin, osteonectin, bone sialoprotein II and osteopontin), and proteoglycans (decorin and biglycan) [13, 14].
In our previous study, we established an optimal extraction process of the enzymatic hydrolysates from Chinese Sika deer sinews (DSEH). The yield of DSEH was 58.56%, and the protein concentration was 48.4%. Cell proliferation assay showed that DSEH significantly promoted the proliferation of MC3T3-E1 cells (an osteoblast-like cell line) with a certain dose-effect relationship, and the proliferative effect was most significant at an administration concentration of 8 mg/ml . Therefore, in the present study, we performed RNA-seq analysis accompanied with validation methods to explore the underlying molecular mechanism of DSEH on MC3T3-E1 cells via deeply dissecting the gene expression patterns under DSEH treatment. We demonstrated that DSEH promoted MC3T3-E1 cell proliferation and extracellular matrix synthesis by regulating multiple functional genes. DSEH significantly increased the expression levels of genes that promoted cell proliferation such as Gstp1, Timp1, Serpine1, Cyr61, Crlf1, Thbs1, Ctgf, P4ha2, Sod3 and Nqo1. However, DSEH significantly decreased the expression levels of genes that inhibited cell proliferation such as Mt1, Cdc20, Gas1, Nrp2, Cmtm3, Dlk2, Sema3a, Rbm25 and Hspb6. Furthermore, DSEH mildly increased the expression levels of osteoblast gene markers. Our findings deepened the current understanding of DSEH on regulating bone development, and provided theoretical support for the discovery of optional prevention and treatment for bone-related diseases.
DSEH treatment and cell proliferation assay
Deer sinews were purchased from the Shuangyang deer farm in Changchun, China. The preparation of the enzymatic hydrolysates from Chinese Sika deer sinews (DSEH) was carried out as previously described . All procedures were submitted and approved by the ethics committee of Changchun University of Chinese Medicine in accordance with the guidelines of the ethical protocol (No. ccucm-2017-0015). Primary osteoblasts from individual calvaria of newborn C57BL/6 mice (n = 6) were isolated according to the previous protocol . Primary osteoblasts and MC3T3-E1 cells (ATCC, USA) were inoculated into a 96-well cell culture plate (Thermo, USA) at a cell density of 5 × 103 cells/well, respectively, and incubated in a humidified incubator (Thermo, USA) containing 5% CO2 at 37 °C for 12 h. The effects of DSEH at different administrated concentrations (0, 2 mg/ml, 4 mg/ml, 6 mg/ml, 8 mg/ml and 10 mg/ml) on cell proliferation were measured using a CCK-8 assay kit (Sigma, USA) and a MTT assay kit (Sigma, USA) according to the manufacturers’ protocols.
Alizarin red staining and quantification
Primary osteoblasts and MC3T3-E1 cells were inoculated into a 24-well cell culture plate at a cell density of 2 × 104 cells/well, and incubated in a humidified incubator (Thermo, USA) containing 5% CO2 at 37 °C for 24 h. The cells were either treated by DSEH at an optimum concentration based on the result of the cell proliferation assay or treated with plain culture medium (Thermo, USA), and incubated for 24 h. All culture mediums were removed, and cells were washed with phosphate-buffered saline buffer (Thermo, USA). Cell fixation was performed at room temperature by incubating with 4% paraformaldehyde for 15 min. After washing with Millipore purified water for three times, the fixed cells were stained with 1% Alizarin red S staining solution (Solarbio, China) for 0.5 h at room temperature, and washed with purified water for three times. Images were taken with an optical microscope (Olympus, Japan) accompanied with a digital camera. Quantitative analysis was carried out by dissolving the stained cells with 10% acetic acid (Thermo, USA) at room temperature for 0.5 h with gentle shaking, followed by vortexing, heating, centrifuging and neutralizing by ammonium hydroxide, and quantified using a plate reader (Life science, USA) at an optical density (OD) of 405 nm .
RNA purification and quality control
MC3T3-E1 cells were inoculated into a 6-well cell culture plate at a cell density of 1 × 106 cells/well, and incubated in a humidified incubator (Thermo, USA) containing 5% CO2 at 37 °C for 12 h. The cells were either treated with DSEH at an optimum concentration based on the result of the cell proliferation assay or treated with plain culture medium (Thermo, USA). All culture mediums were removed, and cells were gently rinsed with precooled phosphate-buffered saline buffer (Thermo, USA). Total RNA from experiments in triplicate was isolated with TRIzol reagent (Invitrogen, USA) according to the manufacturer’s instructions. The total RNA was qualified and quantified based on the RNA integrity number (RIN) measured by an Agilent Bioanalyzer 2100 (Agilent Technologies, USA) following the instrument’s protocols.
Library construction and sequencing
The RNA-seq libraries were generated using a TruSeq Stranded mRNA kit (Illumina, USA) in accordance with the company’s recommendations. Briefly, the messenger RNA containing poly (A) was purified from the total RNA using a magnetic poly (T) beads. Messenger RNA was broke into short fragments with a fragmentation buffer. Double stranded cDNA was synthesized by random hexamer primers, following by purification, repair and sequencing adapter ligation, and amplified by polymerase chain reaction to generate the libraries. The libraries were sequenced by performing high-throughput sequencing via the Illumina HiSeq 2500 platform with a paired-end read length of 150 bp (Illumina, USA).
Sequencing data analysis
After illumine sequencing, the image data was transformed via base calling into raw data in FASTQ format. The raw data were processed using Perl scripts to generate clean reads by removing the low-quality reads and adapter sequences, following by Q30 evaluation and GC content calculation. The data sets were submitted into the NCBI Sequence Read Archive (SRA) database with an accession number PRJNA612675. Subsequently, the clean reads were aligned with the mouse (Mus musculus) reference genome using the HISAT program . The gene expression levels of each transcript were estimated according to the FPKM algorithm . The BLAST program was used to perform annotation against the Non-Redundant (NR) and Swiss-Prot protein databases . The differential expression between the DSEH treated group and untreated group were analyzed by the DEGseq program . Differentially expressed genes (DEGs) were identified using the following criteria: log2 fold change ≥1 or ≤ − 1 and with a p value ≤0.001.
The enrichment analysis of Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) was conducted by R function phyper to determine the biological functions and signaling pathways of those identified DEGs. The significance of enriched functional clusters and signaling pathways were determined based on the adjusted p value (Q value) of less than 0.05, which was calculated according to the Hypergeometric test and Bonferroni correction .
Quantitative real-time PCR (qRT-PCR) verification
The RNA-seq results were further verified through examining the expression levels of DEGs by qRT-PCR assay. Briefly, total RNA was extracted using the TRIzol reagent (Invitrogen, USA) in accordance with the manufacturer’s protocols. Reverse-transcribed cDNA was synthesized using the iScript cDNA Synthesis kit (Bio-Rad, USA), following by qRT-PCR detection under standard amplification condition using a generic SYBR_ Green Supermix kit (Bio-Rad, USA) on a CFX Connect Real-Time PCR Detection System (Bio-Rad, USA). The gene of mouse glyceraldehyde 3-phosphate dehydrogenase (Gapdh) was selected as the internal reference gene. The normalized expression of validated genes was calculated based on the 2−ΔΔCT method .
Effects of DSEH on osteoporotic model caused by hydrogen peroxide (H2O2)-induced oxidative injury
Effects of DSEH on MC3T3-E1 cells under hydrogen peroxide (H2O2)-induced oxidative injury were evaluated to further validate the regulatory function of DSEH on osteoblasts . Briefly, MC3T3-E1 cells were either treated with 500 μM H2O2 for 6 h (H2O2 group), or treated with 500 μM H2O2 for 6 h followed by the addition of 8 mg/ml DSEH (H2O2 + DSEH group), untreated MC3T3-E1 cells serve as a negative control (Blank group). Cells were double stained with propidium iodide (PI) and Annexin V-FITC using a BD Apoptosis Detection Kit (Thermo, USA), and analyzed by a FlowSight® Imaging Flow Cytometer (Merck Millipore, USA). Analysis of cell apoptosis was performed using the IDEAS Application V6.1 software (Amnis, USA).
DSEH promoted MC3T3-E1 cell and primary osteoblast proliferation in a dose-dependent manner
We first investigated the effects of DSEH on the cell proliferation of MC3T3-E1 cells and primary osteoblasts. The results of CCK-8 assay showed that DSEH promoted the proliferation of MC3T3-E1 cells and primary osteoblasts in a dose-dependent manner, as shown in Fig. 1. Compared with the untreated group (0 mg/ml), the viabilities of both cells were significantly increased under the DSEH treatment at progressively increasing concentrations. Meanwhile, the results of MTT assay also showed similar results as those of CCK-8 assy. Since DSEH treatments at the concentration of 8 mg/ml had the highest cell viability, thus, in the following experiments, the concentration of DSEH treatment was selected as 8 mg/ml.
DSEH enhanced MC3T3-E1 cell extracellular matrix synthesis in a time-dependent manner
In order to investigate the effects of DSEH on osteoblast extracellular matrix synthesis, we performed Alizarin red staining and quantification at different time points as follows, 1 d, 3 d, 5 d, 7 d, 11 d and 13 d. As shown in Fig. 2, DSEH significantly enhanced extracellular matrix synthesis of MC3T3-E1 cells at day 1 and day 3 with a similar pattern, but inhibit extracellular matrix synthesis of MC3T3-E1 cells from day 7 to day 13. However, DSEH showed no significant effects on extracellular matrix synthesis of primary osteoblasts.
Gene expression analysis of MC3T3-E1 cells with or without DSEH treatment was performed by RNA sequencing, genome mapping, and transcript annotation
To explore the molecular mechanism underlying the effects of DSEH on MC3T3-E1 cells, we performed RNA-seq analysis on MC3T3-E1 cells with or without DSEH treatment. As shown in Table 1, after Illumina sequencing and data processing, 62,258,284 and 42,159,210 clean reads were obtained from MC3T3-E1 cells treated with plain medium (Blank group) and those treated with DSEH (DSEH group), respectively. The quality assessment showed that the Q30 percentages were more than 92%, and the GC content percentages were approximately 52%. For the Blank and DSEH-treated samples, 59,578,848 and 40,340,992 clean reads were mapped to the mouse genome, respectively. In total, 13,962 out of 14,295 (Blank group), and 14,057 out of 14,373 (DSEH group) transcripts were annotated by searching against the Non-Redundant (NR) NCBI protein database and Swiss-Prot protein database, respectively.
A series of DEGs were identified in MC3T3-E1 cells with or without DSEH treatment
Differential expression analysis was carried out to discover DEGs significantly changed under DSEH treatment. As shown in Table 2, in total, 841 genes with significant differences in the expression levels between DSEH group and Blank group were identified. Among those identified DEGs, 494 DEGs were significantly up-regulated, and 347 genes were significantly down-regulated (DSEH group vs. Blank group). The top 20 up-regulated DEGs mainly consisted of glutathione S-transferase P1 (Gstp1), metalloproteinase inhibitor 1 (Timp1), plasminogen activator inhibitor 1 (Serpine1), protein CYR61 (Cyr61), H-2 class I histocompatibility antigen D-B alpha chain (H2d1) and connective tissue growth factor (Ctgf), etc., as shown in Table 3, whereas the top 20 down-regulated DEGs mainly consisted of metallothionein-1 (Mt1), non-histone chromosomal protein HMG-17 (Hmgn2), prolactin-2C2 (Prl2c2), cell division cycle protein 20 homolog (Cdc20), prolactin-2C3 (Prl2c3) and transmembrane glycoprotein NMB (Gpnmb), etc., as shown in Table 4.
DEGs were classified into a variety of functional categories by functional enrichment analysis
GO enrichment analysis was performed to specify the distribution of these identified DEGs in cellular location, function and physiological processes. According to the functional enrichment analysis of DEGs, the identified DEGs were classified into the following GO categories: cellular component, molecular function and biological process, as shown in Fig. 3. Cellular component classification showed that most of the DEGs were located in the regions of extracellular matrix, intracellular part and organelle. Molecular function classification showed that the dominant functions of these DEGs were involved in binding activities, such as protein binding and ion binding. Biological process classification showed that these DEGs mainly participated in the biological processes including regulation of metabolic process, single-multicellular organism process and developmental process.
DEGs were mapped into multiple signaling pathways by pathway enrichment analysis
KEGG enrichment analysis was carried out to uncover the related signaling pathways by searching the identified DEGs against the KEGG database. As shown in Fig. 4, the DEGs were mainly mapped to the following nine signaling pathways, including thyroid hormone signaling pathway, regulation of actin cytoskeleton, PI3K-Akt signaling pathway, FoxO signaling pathway, focal adhesion, ECM-receptor interaction, cAMP signaling pathway and axon guidance.
DSEH increased the expression levels of genes positively regulating the proliferation of MC3T3-E1 cells
Next, we screened the DEGs involved in positively regulating cell proliferation to further dissect the molecular mechanisms of DSEH on regulating MC3T3-E1 cells. Among the identified DEGs, there were 10 up-regulated DEGs, including glutathione S-transferase P1 (Gstp1), metalloproteinase inhibitor 1 (Timp1), metalloproteinase inhibitor 1 (Serpine1), protein CYR61 (Cyr61), cytokine receptor-like factor 1 (Crlf1), thrombospondin-1 (Thbs1), connective tissue growth factor (Ctgf), prolyl 4-hydroxylase subunit alpha-2 (P4ha2), extracellular superoxide dismutase (Sod3) and NAD(P) H dehydrogenase [quinone] 1 (Nqo1), which positively regulated the proliferation of osteoblasts, as shown in Table 5.
DSEH decreased the expression levels of genes negatively regulating the proliferation of MC3T3-E1 cells
Then, we screened the DEGs involved in negatively regulating the proliferation of MC3T3-E1 cells. In consistent with the above results, among the identified DEGs, there were 9 down-regulated DEGs, including metallothionein-1 (Mt1), cell division cycle protein 20 homolog (Cdc20), growth arrest-specific protein 1 (Gas1), neuropilin-2 (Nrp2), CKLF-like MARVEL transmembrane domain-containing protein 3 (Cmtm3), semaphorin-3A (Sema3a), protein delta homolog 2 (Dlk2), RNA-binding protein 25 (Rbm25) and heat shock protein beta-6 (Hspb6), which negatively regulated the proliferation of osteoblasts, as shown in Table 6.
DSEH mildly increased the expression levels of osteoblast markers in MC3T3-E1 cells
Finally, we analyzed the expression levels of osteoblast markers under DSEH treatment in order to have a deeper recognition regarding the effects of DSEH on MC3T3-E1 cells. The results were consistent with our above findings that DSEH mildly increased the expression levels of multiple osteoblast marker genes including osteopontin (Spp1), collagen alpha-1(I) chain (Col1a1), fibronectin (Fn1), cyclic AMP-dependent transcription factor ATF-4 (Atf4), bone marrow stromal antigen 2 (Bst2) and bone morphogenetic protein 1 (Bmp1), etc., as shown in Table 7.
The accuracy of RNA-seq results was further consolidated by qRT-PCR assay
To further verify and consolidate the accuracy of the RNA-seq data, 12 DEGs including 6 up-regulated genes (Gstp1, Timp1, Serpine1, Cyr61, Crlf1 and Thbs1) and 6 down-regulated genes (Mt1, Cdc20, Gas1, Nrp2, Cmtm3 and Dlk2) were chosen to verify the accuracy of RNA-seq results by analyzing their expression levels using qRT-PCR assay. Gene specific primer sequences were shown in Table 8. The relative mRNA expression levels were presented in Fig. 5. The results showed that the qRT-PCR assay obtained consistent results with those obtained by RNA-seq analysis.
DSEH inhibited MC3T3-E1 cell apoptosis caused by H2O2-induced oxidative injury
In order to confirm the results obtained by cell proliferation assay, Alizarin red staining and RNA-seq analysis, effects of DSEH on MC3T3-E1 cells under H2O2-induced oxidative injury were analyzed using Annexin V/PI staining and flow cytometry following 24 h of incubation with DSEH. As shown in Fig. 6, compared with the negative control group (Blank), the percentage of apoptotic cells (R4: early apoptotic cells; R5: late apoptotic cells) were increased under H2O2 induction (H2O2). However, DSEH treatment (H2O2 + DSEH) decreased the percentage of apoptotic cells, especially the late apoptotic cells.
Sinew shares a series of types of collagen components with bone, including type I, III and V collagens, of which type I collagen is expressed in the whole process of bone growth and development, and type III collagen play pivotal role in regulating bone regeneration during the process of bone repair . Deer Sinew serves as a medicinal food, and has been used for treating skeletal diseases, especially bone diseases in a long history. Thus, it could become an alternative option for the prevention and therapeutic remedy of bone-related diseases. However, the exact regulation mechanism remains to be elucidated. In the present study, the results of CCK-8 assay showed that DSEH promoted the proliferation of both MC3T3-E1 cells and primary osteoblasts in a dose-dependent manner. Based on the analysis of Alizarin red staining and quantification, DSEH significantly enhanced extracellular matrix synthesis of MC3T3-E1 cells at day 1 and day 3 with a similar pattern, but inhibit extracellular matrix synthesis of MC3T3-E1 cells from day 7 to day 13. However, DSEH showed no significant effects on extracellular matrix synthesis of primary osteoblasts. Therefore, we further explored the molecular regulation of DSEH on the gene expression of MC3T3-E1 cells using RNA-seq technique together with qRT-PCR verification method, and demonstrated that DSEH promoted the proliferation of MC3T3-E1 cells and suppressed the apoptosis of those cells by regulating multiple functional genes.
According to our findings, we addressed that DSEH promoted osteoblast proliferation and extracellular matrix synthesis by increasing the expression levels of genes positively regulating the proliferation of MC3T3-E1 cells, including Gstp1, Timp1, Serpine1, Cyr61, Crlf1, Thbs1, Ctgf, P4ha2, Sod3 and Nqo1. Among these genes, Gstp1 is a glutathione transferase that acts as a regulator of signaling pathways that are involved in cell proliferation, and contributes in protection against osteoblast apoptosis . Timp1, a member of the tissue inhibitor of metalloproteinases family, increases osteoblast proliferation and decreased apoptosis . Serpine1 is a primary inhibitor of endogenous plasminogen activators, which serves as a regulator during cell growth and cell proliferation, and plays an important role in regulating resorption during bone fracture repair . Cyr61 is a cysteine-rich angiogenic inducer that is capable of regulating a broad range of cellular activities, such as cell adhesion and proliferation, and could stimulate osteoblast proliferation during bone remodeling . Crlf1 is a member of the ciliary neurotrophic factor receptor pathway, and is involved in the anabolic therapy of osteoporosis through regulating osteoblast proliferation and bone formation . Thbs1 is an adhesive glycoprotein that is directly deposited by bone-forming osteoblasts, and is present in the mineralized bone matrix . Ctgf is an extracellular matrix-associated growth factor that plays a pivotal role in osteoblast proliferation and bone formation . P4ha2 is a subunit of collagen prolyl-4-hydroxylase that is mainly expressed in osteoblasts, chondrocytes and capillary endothelial cells, and serves as an important enzyme during collagen biosynthesis . Sod3 is a member of the superoxide dismutase protein family, and has a protective effect on oxidative stress-induced damage in osteoblastic cells . Nqo1 is a pleiotropic enzyme that is expressed in calvarial osteoblasts, and protects the cells from oxidative damage . Thus, these findings suggest that DSEH facilitate osteoblast proliferation, maintain osteoblast homeostasis and protect oxidative stress-induced damage.
We next addressed that DSEH promoted osteoblast proliferation and extracellular matrix synthesis by decreasing the expression levels of genes negatively regulating the proliferation of MC3T3-E1 cells, including Mt1, Cdc20, Gas1, Nrp2, Cmtm3, Dlk2, Sema3a, Rbm25 and Hspb6. Among these genes, Mt1 is a family member of the cysteine-rich metallothioneins, and is mainly expressed during osteoblastic differentiation and mineralization . Cdc20 is an activator of the anaphase-promoting complex, and is weakly expressed in osteoblasts but highly expressed in the osteocytes . Gas1 is a cell growth repressor gene that serves as a suppressive regulator of cell growth and as a mediator of apoptosis . Nrp2 is a family member of the neuropilins, and one of its isoform Nrp2b inhibits cultured cell proliferation . Cmtm3, a member of the chemokine-like factor gene superfamily, is an inhibitor of cell growth, and an inducer of cell apoptosis . Dlk2 is a member of the EGF-like family of membrane proteins, and interacts with Notch signaling to regulate cell growth and apoptosis . Sema3a is a member of the semaphoring family, which suppresses cell proliferation and promotes cell apoptosis . Rbm25, a RNA-binding protein, negatively regulates cell proliferation and positively regulating cell apoptosis . Hspb6, also known as Hsp20, is a small heat shock protein that participates in reducing cell proliferation and inducing cell apoptosis . Thus, these findings suggest that DSEH prevent osteoblast apoptosis, maintain osteoblast homeostasis and prevent osteoblasts form further differentiation.
We further addressed that DSEH mildly increased the expression levels of a majority of the osteoblast markers, including Spp1, Col1a1, Fn1, Atf4, Bst2, Bmp1, Bmp4, Sox4, Sox11, Msx1, Znf25, Bglap2, Sp7, Bglap, Dmp1, Alpl, Bglap3, Ostn and Omd. Among these genes, Spp1, Col1a1, Fn1, Bglap, Bglap2, Bglap3 and Dmp1 are the major bone matrix protein genes [45,46,47]. Atf4, Sox4, Sox11, Msx1, Znf25 and Sp7 are key transcription factors involved in osteoblast commitment, proliferation, differentiation and extracellular matrix synthesis during bone formation [48,49,50,51,52]. Bst2, also known as tetherin, is a lipid raft associated protein that is highly expressed in osteoblasts, and is involved in the process of osteogenic differentiation . Bmp1 and Bmp4, two members of the bone morphogenetic protein family, are important stimulator of osteoblast proliferation and differentiation, bone formation and stability [54, 55]. Alpl is a well-known osteoblastic marker, and has been used as an evaluation indicator to estimate the capacity of bone formation . Ostn is a secreted bone-active protein that is highly expressed in osteoblasts, and plays key role in regulating osteoblast activity in developing bone and at the sites of bone remodeling . Omd, also known as osteoadherin, is mainly expressed in osteoblasts at early differentiated stage, and serves as an organizer of bone mineral formation . However, the expression level of Runx2 was mildly decreased under DSEH treatment, and the underlying molecular mechanism still need to be further studied. Thus, these findings suggest that DSEH facilitate osteoblast proliferation, differentiation and extracellular matrix synthesis to consolidate bone formation and stability.
Consistent with the above results, our enrichment analysis of GO function and KEGG pathway indicated that the DEGs were primarily classified into the cellular components of extracellular matrix, intracellular part and organelle with dominant molecular functions of binding activities, and predominantly participated in the biological processes including regulation of metabolic process, single-multicellular organism process and developmental process. Those DEGs were mainly mapped to the signaling pathways such as thyroid hormone signaling pathway, regulation of actin cytoskeleton, PI3K-Akt signaling pathway, FoxO signaling pathway, focal adhesion, ECM-receptor interaction, cAMP signaling pathway and axon guidance. Among these pathways, thyroid hormone signaling pathway plays a key role in the regulation of bone metabolism by modulating osteoblast activity during skeletal development . Actin cytoskeleton has been considered to be involved in various functions of osteoblast by discerning cell condition and regulating extracellular matrix mineralization and bone formation . PI3K-Akt signaling pathway is a central regulation hub for osteoblast homeostasis and function, and has the capacity to work together with other signaling networks to regulate osteoblast phenotype and bone formation . FoxO signaling is indispensable for bone mass homeostasis since this pathway could prevent oxidative stress and osteoblast apoptosis, and promote osteoblast proliferation to facilitate bone formation and stability . Focal adhesion facilitates the perpendicular arrangement of bone matrix orthogonal to osteoblast alignment during bone formation . ECM-receptor interaction plays a crucial role in regulating osteoblast adhesion, proliferation and survival during bone formation and repair . cAMP signaling pathway is capable of enhancing cell adhesion of osteoblasts to facilitate bone regeneration . Axon guidance signaling pathway plays a pivotal role in regulating bone formation and resorption by promoting osteoblast migration and suppressing osteoclast differentiation . Finally, the results of DSEH on osteoporotic model caused by hydrogen peroxide (H2O2)-induced oxidative injury further indicated that DSEH treatment decreased the percentage of apoptotic cells, especially the late apoptotic cells. Thus, these findings confirmed that DSEH facilitated osteoblast proliferation and extracellular matrix synthesis, and prevented osteoblast oxidative stress and apoptosis to consolidate bone formation and repair.
In summary, we performed cell proliferation assay, Alizarin red staining and quantification, RNA-seq analysis accompanied with verification methods to dissect the molecular mechanism of DSEH in the regulation of an osteoblast-like cell line MC3T3-E1. Our findings suggest that DSEH facilitate osteoblast proliferation and extracellular matrix synthesis to consolidate bone formation and stability, but prevent osteoblasts from oxidative stress-induced damage, apoptosis and further differentiation. These findings deepened the current understanding of DSEH on regulating bone development, and provided theoretical support for the discovery of optional prevention and treatment for bone-related diseases.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Deer sinews enzymatic hydrolysate
RNA integrity number
Basic Local Alignment Search Tool
Differentially expressed genes
Quantitative real-time PCR
Sequence Read Archive
Glyceraldehyde 3-phosphate dehydrogenase
Kyoto Encyclopedia of Genes and Genomes
Jiang LL, Liu CL, Wong YL, Nip CF, Shaw PC. Differentiation of deer tendons from cattle tendons by a loop-mediated isothermal amplification (LAMP) test and bone remodeling bioassays. Chin Med. 2015;10:33.
Zhang H, Dong Y, Qi B, et al. Preventive effects of collagen peptide from deer sinew on bone loss in ovariectomized rats. Evid Based Complement Alternat Med. 2014;2014:627285.
Zhang H, Zhao Y, Li YQ, Sun XD, Bai XY, Zhao DQ. Effects of deer tendons collagen on osteoporosis rats induced by retinoic acid. Zhong Yao Cai. 2010;33(3):411–4.
Wang IL, Hsiao CY, Shen J, Wang Y, Huang CC, Chen YM. The effects of Jilin sika Deer’s (Cervus dybowski) tendon liquid supplementation on endurance drop jumps performance, biochemistry profile of free boxing players. J Ethnopharmacol. 2019;245:112119.
Sun X, Li Y, Zhao Y, Qu Y, Zhang H, Qu X. Therapeutic effect of deer sinew collagens on rat osteoarthritis model. Zhongguo Zhong Yao Za Zhi. 2009;23:3135–8.
Zhao Y, Li Y, Zhang H, Bai X, Han Y. Study on anti-inflammatory and analgesic effects of deer sinew collagens. Liaoning J Trad Chin Med. 2010;09:1825–7.
Sun X, Li Y, Zhao Y, Xing N, Bai X, Qu X. Study on the effects of anti-inflammatory and immune regulation of deer sinew collagens. Lishizhen Med Mater Med Res. 2010;04:853–4.
Thorpe CT, Screen HR. Tendon structure and composition. Adv Exp Med Biol. 2016;920:3–10.
Kannus P. Structure of the tendon connective tissue. Scand J Med Sci Sports. 2000;10(6):312–20.
Viguet-Carrin S, Garnero P, Delmas PD. The role of collagen in bone strength. Osteoporos Int. 2006;17(3):319–36.
Xu Z, Asakawa S. Physiological RNA dynamics in RNA-Seq analysis. Brief Bioinform. 2019;20(5):1725–33.
Ayturk U. RNA-seq in skeletal biology. Curr Osteoporos Rep. 2019;17(4):178–85.
Shapiro F. Bone development and its relation to fracture repair. The role of mesenchymal osteoblasts and surface osteoblasts. Eur Cell Mater. 2008;15:53–76.
Florencio-Silva R, Sasso GR, Sasso-Cerri E, Simões MJ, Cerri PS. Biology of bone tissue: structure, function, and factors that influence bone cells. Biomed Res Int. 2015;2015:421746.
Zhang P, Yao B, Leng X. Optimization of extraction conditions of collagen from deer sinew. J Agric Sci Yanbian Univ. 2020;42(1):69–75.
Ducy P, Starbuck M, Priemel M, Shen J, Pinero G, Geoffroy V, Amling M, Karsenty G. A Cbfa1-dependent genetic pathway controls bone formation beyond embryonic development. Genes Dev. 1999;13(8):1025–36.
Gregory CA, Gunn WG, Peister A, Prockop DJ. An Alizarin red-based assay of mineralization by adherent cells in culture: comparison with cetylpyridinium chloride extraction. Anal Biochem. 2004;329(1):77–84.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Trapnell C, Williams BA, Pertea G, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.
Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.
Yao B, Liu J, Zhang M, Leng X, Zhao D. Deciphering the potential pharmaceutical mechanism of Guzhi Zengsheng Zhitongwan on rat bone and kidney based on the “kidney governing bone” theory. J Orthop Surg Res. 2020;15(1):146.
Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative C(T) method. Nat Protoc. 2008;3(6):1101–8.
Fatokun AA, Stone TW, Smith RA. Hydrogen peroxide-induced oxidative stress in MC3T3-E1 cells: the effects of glutamate and protection by purines. Bone. 2006;39(3):542–51.
Liu SH, Yang RS, al-Shaikh R, Lane JM. Collagen in tendon, ligament, and bone healing. A current review. Clin Orthop Relat Res. 1995;318:265–78.
Ragab SM, Badr EA, Ibrahim AS. Evaluation of glutathione-S-transferase P1 polymorphism and its relation to bone mineral density in Egyptian children and adolescents with beta-thalassemia major. Mediterr J Hematol Infect Dis. 2016;8(1):e2016004.
Merciris D, Schiltz C, Legoupil N, Marty-Morieux C, de Vernejoul MC, Geoffroy V. Over-expression of TIMP-1 in osteoblasts increases the anabolic response to PTH. Bone. 2007;40(1):75–83.
Rundle CH, Wang X, Wergedal JE, Mohan S, Lau KH. Fracture healing in mice deficient in plasminogen activator inhibitor-1. Calcif Tissue Int. 2008;83(4):276–84.
Liu H, Peng F, Liu Z, et al. CYR61/CCN1 stimulates proliferation and differentiation of osteoblasts in vitro and contributes to bone remodeling in vivo in myeloma bone disease. Int J Oncol. 2017;50(2):631–9.
Standal T, Johnson RW, McGregor NE, et al. gp130 in late osteoblasts and osteocytes is required for PTH-induced osteoblast differentiation. J Endocrinol. 2014;223(2):181–90.
Amend SR, Uluckan O, Hurchla M, et al. Thrombospondin-1 regulates bone homeostasis through effects on bone matrix integrity and nitric oxide signaling in osteoclasts. J Bone Miner Res. 2015;30(1):106–15.
Safadi FF, Xu J, Smock SL, et al. Expression of connective tissue growth factor in bone: its role in osteoblast proliferation and differentiation in vitro and bone formation in vivo. J Cell Physiol. 2003;196(1):51–62.
Xiong G, Deng L, Zhu J, Rychahou PG, Xu R. Prolyl-4-hydroxylase α subunit 2 promotes breast cancer progression and metastasis by regulating collagen deposition. BMC Cancer. 2014;14:1.
Jung WW. Protective effect of apigenin against oxidative stress-induced damage in osteoblastic cells. Int J Mol Med. 2014;33(5):1327–34.
Sun YX, Xu AH, Yang Y, Li J. Role of Nrf2 in bone metabolism. J Biomed Sci. 2015;22:101.
Dohi Y, Shimaoka H, Ikeuchi M, Ohgushi H, Yonemasu K, Minami T. Role of metallothionein isoforms in bone formation processes in rat marrow mesenchymal stem cells in culture. Biol Trace Elem Res. 2005;104(1):57–70.
Billiard J, Moran RA, Whitley MZ, et al. Transcriptional profiling of human osteoblast differentiation. J Cell Biochem. 2003;89(2):389–400.
Dong C, Wang X, Li N, et al. microRNA-mediated GAS1 downregulation promotes the proliferation of synovial fibroblasts by PI3K-Akt signaling in osteoarthritis. Exp Ther Med. 2019;18(6):4273–86.
Gemmill RM, Nasarre P, Nair-Menon J, et al. The neuropilin 2 isoform NRP2b uniquely supports TGFβ-mediated progression in lung cancer. Sci Signal. 2017;10(462):eaag0528.
Zhang H, Zhang J, Nan X, et al. CMTM3 inhibits cell growth and migration and predicts favorable survival in oral squamous cell carcinoma. Tumour Biol. 2015;36(10):7849–58.
Huang Q, Li J, Zheng J, Wei A. The carcinogenic role of the notch signaling pathway in the development of hepatocellular carcinoma. J Cancer. 2019;10(6):1570–9.
Capparuccia L, Tamagnone L. Semaphorin signaling in cancer cells and in cells of the tumor microenvironment--two sides of a coin. J Cell Sci. 2009;122(Pt 11):1723–36.
Ge Y, Schuster MB, Pundhir S, et al. The splicing factor RBM25 controls MYC activity in acute myeloid leukemia. Nat Commun. 2019;10(1):172.
Nagasawa T, Matsushima-Nishiwaki R, Toyoda H, Matsuura J, Kumada T, Kozawa O. Heat shock protein 20 (HSPB6) regulates apoptosis in human hepatocellular carcinoma cells: direct association with Bax. Oncol Rep. 2014;32(3):1291–5.
Komori T. Regulation of proliferation, differentiation and functions of osteoblasts by Runx2. Int J Mol Sci. 2019;20(7):1694.
Diegel CR, Hann S, Ayturk UM, et al. An osteocalcin-deficient mouse strain without endocrine abnormalities. PLoS Genet. 2020;16(5):e1008361.
Lu Y, Yuan B, Qin C, et al. The biological function of DMP-1 in osteocyte maturation is mediated by its 57-kDa C-terminal fragment. J Bone Miner Res. 2011;26(2):331–40.
Elefteriou F, Benson MD, Sowa H, et al. ATF4 mediation of NF1 functions in osteoblast reveals a nutritional basis for congenital skeletal dysplasiae. Cell Metab. 2006;4(6):441–51.
Lefebvre V, Bhattaram P. SOXC genes and the control of skeletogenesis. Curr Osteoporos Rep. 2016;14(1):32–8.
Komori T. Regulation of osteoblast differentiation by transcription factors. J Cell Biochem. 2006;99(5):1233–9.
Twine NA, Harkness L, Kassem M, Wilkins MR. Transcription factor ZNF25 is associated with osteoblast differentiation of human skeletal stem cells. BMC Genomics. 2016;17(1):872.
Hojo H, Ohba S, He X, Lai LP, McMahon AP. Sp7/Osterix is restricted to bone-forming vertebrates where it acts as a Dlx co-factor in osteoblast specification. Dev Cell. 2016;37(3):238–53.
Yoo SH, Kim JG, Kim BS, et al. BST2 mediates osteoblast differentiation via the BMP2 signaling pathway in human alveolar-derived bone marrow stromal cells. PLoS One. 2016;11(6):e0158481.
Asharani PV, Keupp K, Semler O, et al. Attenuated BMP1 function compromises osteogenesis, leading to bone fragility in humans and zebrafish. Am J Hum Genet. 2012;90(4):661–74.
Chang SF, Chang TK, Peng HH, et al. BMP-4 induction of arrest and differentiation of osteoblast-like cells via p21 CIP1 and p27 KIP1 regulation. Mol Endocrinol. 2009;23(11):1827–38.
Liu W, Zhang L, Xuan K, et al. Alpl prevents bone ageing sensitivity by specifically regulating senescence and differentiation in mesenchymal stem cells. Bone Res. 2018;6:27.
Bord S, Ireland DC, Moffatt P, Thomas GP, Compston JE. Characterization of osteocrin expression in human bone. J Histochem Cytochem. 2005;53(10):1181–7.
Rehn AP, Chalk AM, Wendel M. Differential regulation of osteoadherin (OSAD) by TGF-beta1 and BMP-2. Biochem Biophys Res Commun. 2006;349(3):1057–64.
Kim HY, Mohan S. Role and mechanisms of actions of thyroid hormone on the skeletal development. Bone Res. 2013;1(2):146–61.
Suzuki H, Tatei K, Ohshima N, Sato S, Izumi T. Regulation of MC3T3-E1 differentiation by actin cytoskeleton through lipid mediators reflecting the cell differentiation stage. Biochem Biophys Res Commun. 2019;514(2):393–400.
McGonnell IM, Grigoriadis AE, Lam EW, Price JS, Sunters A. A specific role for phosphoinositide 3-kinase and AKT in osteoblasts? Front Endocrinol (Lausanne). 2012;3:88.
Ambrogini E, Almeida M, Martin-Millan M, et al. FoxO-mediated defense against oxidative stress in osteoblasts is indispensable for skeletal homeostasis in mice. Cell Metab. 2010;11(2):136–46.
Nakanishi Y, Matsugaki A, Kawahara K, Ninomiya T, Sawada H, Nakano T. Unique arrangement of bone matrix orthogonal to osteoblast alignment controlled by Tspan11-mediated focal adhesion assembly. Biomaterials. 2019;209:103–10.
Shekaran A, García AJ. Extracellular matrix-mimetic adhesive biomaterials for bone repair. J Biomed Mater Res A. 2011;96(1):261–72.
Lo KW, Ashe KM, Kan HM, Lee DA, Laurencin CT. Activation of cyclic amp/protein kinase: a signaling pathway enhances osteoblast cell adhesion on biomaterials for regenerative engineering. J Orthop Res. 2011;29(4):602–8.
Iqbal J, Yuen T, Kim SM, Zaidi M. Opening windows for bone remodeling through a SLIT. J Clin Invest. 2018;128(4):1255–7.
This work was supported by the National Key Research and Development Program of China (Grant No. 2018YFC1706605), the TCM Clinical Research Center for Bone diseases of Jilin Province (Grant No. 20180623048TC) and the Science and Technology Development Plan Project of Jilin Province (Grant No. 20190304095YY). The funding bodies provided financial support, and the awardees performed the research. The founding sponsor had no role in the design of the study, the collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
All procedures were submitted and approved by the ethics committee of Changchun University of Chinese Medicine in accordance with the guidelines of the ethical protocol (No. ccucm-2017-0015).
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.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Zhou, Z., Zhao, D., Zhang, P. et al. The enzymatic hydrolysates from deer sinew promote MC3T3-E1 cell proliferation and extracellular matrix synthesis by regulating multiple functional genes. BMC Complement Med Ther 21, 59 (2021). https://doi.org/10.1186/s12906-021-03240-2
- Deer sinew
- Enzymatic hydrolysates
- MC3T3-E1 cells
- RNA sequencing
- Differentially expressed genes
- Molecular mechanism