- Open Access
Suv39h1 promotes facet joint chondrocyte proliferation by targeting miR-15a/Bcl2 in idiopathic scoliosis patients
Clinical Epigenetics volume 11, Article number: 107 (2019)
Idiopathic scoliosis (IS) is a complex disease with an unclear etiology, and the worldwide prevalence is approximately 2–3%. As an important link between environmental factors and phenotypic differences, epigenetic changes, such as lncRNA, miRNA, and DNA methylation, have recently been reported to be associated with the development of IS. However, the correlation between histone methylation, another classical epigenetic mechanism, and IS has not been determined. In this study, we investigated the morphological changes, alterations in the levels of histone methylation, and cell proliferation-related pathway in inferior facet joint cartilage in 11 IS patients and 10 comparable controls.
Compared with the control group, narrowed facet joint cartilage but increased proliferative chondrocytes and upregulated collagen type II (COL2A1) and B-cell lymphoma-2 (Bcl2) were observed in IS patients. Additionally, tri-methylation levels of H3K9 (H3K9me3) rather than other lysine sites were significantly increased in IS patients, coinciding with the upregulation of its specific enzyme, suppressor of variegation 3-9, drosophila homolog of 1 (SUV39H1). In addition, Bcl2-targeted miR-15a was downregulated in IS patients, and the level of H3K9me3 in the promoter region of the miR-15a host gene was remarkably increased in IS patients compared with the control group. Moreover, overexpressing SUV39H1 in ATDC5 cells with increased H3K9me3 levels led to similar changes, with increased expression of COL2A1 and Bcl2, decreased expression of miR-15a, and increased cell proliferation.
Thus, our study suggests that increased chondrocyte proliferation occurs in the facet joint cartilage of IS patients compared with the control group and may be promoted by the elevated levels of H3K9me3 and SUV39H1, which regulate the miR-15a/Bcl2 pathway. This dysregulation of chondrocyte proliferation could result in abnormal spinal growth and may additionally participate in the development and progression of IS.
Idiopathic scoliosis (IS) is a three-dimensional spinal deformity disease that mostly occurs during adolescence (adolescent idiopathic scoliosis, AIS). The worldwide prevalence of AIS ranges from 0.47–5.2% . Although the exact etiology of IS remains unknown, several theories related to the development of IS have been established, including genetics, endocrine metabolism, biochemics, biomechanics, and environmental factors .
Some studies have suggested that the fundamental problem in scoliosis is the imbalanced growth between the anterior and posterior components of the spine, manifested by longer anterior elements of the spine than the posterior elements in IS patients [3, 4]. Moreover, high spinal growth velocity is predisposed to a more rapid curve progression . In addition to growth-related factors, joint hypermobility may be susceptible to spinal instability and the development of IS, since researchers have found that girls with IS have significantly more laxity than that of healthy sex- and age-matched controls [6, 7]. Pal et al. found that scoliosis developed rapidly and was considerable after both transverse processes and facet joints were removed . Therefore, abnormal growth of the vertebral column indeed may be tightly associated with the development of IS.
Discordant phenotypes for monozygotic twins with AIS suggested that environmental factors are important to the development of IS, although what kind of environmental factors may be involved is unclear . The influence of environmental factors can be linked to phenotypic differences by epigenetic changes, mechanisms of which mainly contain DNA methylation, histone modification, and noncoding RNAs [10, 11]. Epigenetic changes can regulate phenotypic plasticity with complex regulation during the early stage of life but are also involved in the development of many diseases, such as cancer, diabetes, cardiovascular disease, and liver disease [9, 10, 12, 13].
Similarly, changes in epigenetics were also found in AIS patients compared with the control group. A total of 139 long noncoding RNAs (lncRNAs) were differentially expressed between AIS patients and healthy controls, and 4 lncRNAs were expressed differently between different patients when grouped according to age, height, classification, severity of scoliosis, and Risser grade . Garcia-Gimenez et al. observed several deregulated microRNAs (miRNAs) that participate in the differentiation mechanisms of osteoblasts/osteoclasts in AIS patients . Ogura et al. found that the risk allele of a functional variant in MIR4300HG, the host gene of miR4300, can lead to a decrease in miR4300 and be related to AIS progression . In addition, AIS patients had significantly higher COMP (cartilage oligomeric matrix protein) promoter methylation and lower gene expression compared with those of the control groups, which was correlated with young age and a high Cobb angle of the main curve . Although strong associations between rs12459350 in Disruptor of telomeric silencing 1-like (DOT1L) that catalyze methylation of histone H3 lysine 79 (H3K79) and susceptibility of AIS were found , the correlation between alteration of histone methylation and development of IS has not been discovered. In this study, we examined histone methylation levels in facet joint chondrocytes, which are important to spinal growth, in IS patients and the control group, and further investigation between changes in histone methylation and the development of IS was conducted.
General characteristics of IS patients and control groups
Overall, a population of 11 IS patients and 10 comparable controls were investigated in our study. The mean age of IS patients was 16.86 years, ranging from 15 to 21 years, while the mean age of the control group was 19.80 years, ranging from 16 to 28 years (Table 1 and Additional file 1: Table S1). However, there was no significant difference between the ages of the two groups (Table 1). Moreover, the gender difference between the two groups was not observed, since there were 6 males/4 females in the control group and 5 males/6 females in the IS patient group (Table 1). Additionally, the average Cobb angle of the major curve in IS patients was 45.59° (Table 1).
Significant changes in the cartilage layer of the facet joint in IS patients
The histopathological features of facet joint cartilage in IS patients and the control group were analyzed; however, severely damaged facet joints of IS patients that lacked an obvious cartilage layer were excluded from the analysis. In general, articular cartilage has a highly organized structure composed of 4 zones: the superficial zone, middle zone, deep zone, and calcified zone . The chondrocytes in the superficial region are smaller and flatter when compared to the middle and deep zones where cells are organized in a more columnar-shaped cartilage lacuna . As depicted in Fig. 1 a, hematoxylin-eosin (H&E) staining results suggested that the surface, middle, and deep zones of the cartilage were well preserved in both the control group and IS patients. However, a significantly narrowed cartilage layer of the facet joint was observed in IS patients compared with the control group (Fig. 1a, b). In addition, changes in the chondrocyte phenotype in the middle and deep zones of cartilage were observed between the two groups, showing increased smaller and fewer columnar chondrocytes in IS patients compared with the control group (Fig. 1a, b).
Furthermore, we examined the mRNA levels of extracted primary chondrocyte extracellular matrix components aggrecan (ACAN), collagen type II (COL2A1), and collagen type X (COL10A1), which is a reliable marker for new bone formation in articular cartilage [21, 22]. As shown in Fig. 1 d and e, the mRNA level of COL2A1 was significantly increased in IS patients compared with the control group, whereas COL10A1 markedly decreased. There was no significant change in ACAN mRNA levels between the two groups (Fig. 1c).
Alteration in histone methylation in the facet joint chondrocytes of IS patients compared with the control group
Given that histone methylation plays important roles in chondrocyte proliferation and hypertrophy , we detected the histone methylation levels of facet joint chondrocytes in both IS patients and the control group. According to the WB results, a significant increase in the H3K9me3 level, which is a marker of gene repression, was found in IS patients compared with the control group (Fig. 2a) . However, another marker of gene repression, H3K27me3, was not changed in IS patients compared with the control group (Fig. 2a) . Surprisingly, the markers of gene activation, such as H3K4me3, H3K36me3, and H3K79me3, did not change between the two groups (Fig. 2a) . Next, we investigated the expression levels of H3K9 tri-methyltransferases, including suppressor of variegation 3-9, drosophila homolog of 1 and 2 (SUV39H1, SUV39H2) . Although no changes in SUV39H2 expression were observed, the mRNA and protein levels of SUV39H1 were increased in IS patients compared with the control group (Fig. 2b–d). Meanwhile, no significant changes of the mRNA levels of another H3K9 methyltransferase SET domain bifurcated-1 (SETDB1) and H3K9 demethylases lysine (K)-specific demethylases 4 (KDM4A, KDM4B, and KDM4C) were observed between the two groups (Additional file 2: Figure S1).
Increased proliferation of chondrocytes in the facet joint cartilage of IS patients
We further measured chondrocyte proliferation in the facet joint cartilage of IS patients and the control group. In accordance with the H&E staining results, the immunofluorescent staining of proliferating cell nuclear antigen (PCNA), a marker for proliferation, suggested markedly increased proliferative cells of the facet joint chondrocytes in IS patients compared with the control group (Fig. 3a, b). Moreover, the expression of the anti-apoptosis protein B-cell lymphoma-2 (Bcl2) was increased in both the mRNA and protein levels of chondrocytes in IS patients, but no significant changes in either the mRNA or protein level of the pro-apoptosis protein Bcl2-associated x (Bax) were found in IS patients compared with the control group (Fig. 3c–e).
These data showed that the increased smaller chondrocytes we observed in IS patients may result from the enhanced proliferation of facet joint cartilage, which could be further induced by the increased expression of Bcl2. However, whether these changes were related to the alteration in H3K9me3 levels remains unclear, especially the gene inactivation role of H3K9me3 modification.
H3K9me3 modification caused a decrease in miRNAs targeting Bcl2 in the facet joint chondrocytes of IS patients
Based on the hypothesis that elevated Bcl2 expression may be caused indirectly by H3K9me3 modification, we searched for genes that downregulated the expression of Bcl2. By using miRTarBase (http://mirtarbase.mbc.nctu.edu.tw/php/search.php) and a review of the literature, we found four miRNAs that inhibit Bcl2 expression: the microRNA cluster miR-15a/miR-16, miR-29b, and miR-let7d [25,26,27,28,29]. As shown in Fig. 4 a–d, the expression level of hsa-miR-15a was significantly decreased in IS patients compared with the control group, while hsa-miR-16, hsa-miR-29b, and hsa-miR-let7d did not significantly change between the two groups. The analysis of the miR-15a/miR-16 host gene showed no changes in deleted in lymphocytic leukemia 2 (DLEU2) mRNA levels between IS patients and the control group (Fig. 4e) . We further tested the association between the downregulation of has-miR-15a and elevation in H3K9me3 levels by chromatin immunoprecipitation (ChIP) assay. In accordance with our hypothesis, the H3K9me3 level was significantly increased in the promoter region of DLEU2 (Fig. 4f and Additional file 2: Figure S2).
Changes caused by overexpressing SUV39H1 in ATDC5 cells
To verify the effect of elevated levels of H3K9me3 and SUV39H1 on chondrocytes, we overexpressed the highly conserved gene SUV39H1 in ATDC5 cells, which is regarded as a promising in vitro model to study the factors that influence cell behaviors during chondrogenesis (Additional file 2: Figure S3) . Successful overexpression of SUV39H1 was observed not only at the mRNA level but also at the protein level (Fig. 5a, b). Additionally, overexpression of SUV39H1 led to an increase in the level of H3K9me3 (Fig. 5a). Similar to the results obtained in facet joint chondrocytes, elevated levels of H3K9me3 and SUV39H1 increased both the mRNA and protein levels of Bcl2 rather than Bax (Fig. 6a, b). Although the mRNA level of Col10a1 did not change, the mRNA level of Col2a1 was significantly increased in SUV39H1-overexpressing ATDC5 cells compared with control cells (Fig. 5c). And the 5-ethynyl-2′ -deoxyuridine (EdU) staining results showed increased cell proliferation in SUV39H1-overexpressing ATDC5 cells compared with control cells (Fig. 5d). In addition, the expression level of miR-15a instead of miR-16 was significantly decreased in SUV39H1-overexpressing ATDC5 cells (Fig. 6c). In accordance with the changes observed in IS patients, the H3K9me3 level was also significantly increased in the promoter region of Dleu2 in SUV39H1-overexpressing ATDC5 cells, although no differences were found in the mRNA levels of Dleu2 between SUV39H1-overexpressing and control ATDC5 cells (Fig. 6d, e and Additional file 2: Figure S4).
Taken together, these results suggest that H3K9me3 levels elevated by SUV39H1 in the promoter of the miR-15a/miR-16 host gene may decrease the expression of miR-15a inconsistent with its host gene expression, which could further increase the expression of Bcl2. The increased expression of Bcl2 may promote chondrocyte proliferation and ultimately impact spinal growth.
Environmental and lifestyle factors are thought to play an important role in the onset of IS in light of discordant findings for monozygotic twins with AIS [9, 32]. In general, environmental and lifestyle factors could be hormones, nutrition, traffic pollution, noise, alcohol, smoking, viruses, drugs, medications, toxins, physical activity, playing a musical instrument, the type and weight of school bags, the time spent studying, sleep posture, and the number of hours spent sleeping [2, 9, 32]. Although no lifestyle-related factor was identified to be significantly associated with AIS, classical ballet training, a family history of scoliosis, and low BMI may be associated with AIS . Moreover, early exposure to indoor heated swimming pools has been found to increase the independent odds of AIS . Therefore, epigenetics, which is thought to represent the link between environmental factors and phenotypic differences, is a new approach to examine the etiopathogenesis of IS, given that epigenetic differences have arisen during the lifetime of phenotypic discordant monozygotic twins [9, 34].
So far, both DNA methylation and noncoding RNAs have been reported to be connected with IS in some studies [14,15,16,17]. However, little research has been conducted on the association between histone modification and IS. As one of the major histone modifications, histone methylation is deeply involved in endochondral ossification and articular cartilage development . Ito et al. found that the histone methylase protein arginine N-methyltransferase-4 (PRMT4) could regulate endochondral osteogenesis and chondrocyte proliferation through the methylation of sex-determining region Y-box 9 (Sox9) . Conditional inactivation of H3K27-specific methyltransferase enhancer of zeste homolog 2 (EZH2) prevented craniofacial bone and cartilage formation . H3K9 methyltransferases, such as SUV39H1 and SUV39H2, and H3K9me3 were detected in prehypertrophic and hypertrophic chondrocytes in the growth plate of fetal mice . Czvitkovich et al. observed that overexpression of SUV39H1 induces altered proliferation and differentiation in transgenic mice . In addition, DOT1L, a single nucleotide polymorphism, was reported to be strongly associated with IS and played a role in chondrogenic differentiation and adult articular cartilage [18, 39].
In our study, we examined the tri-methylation levels of H3K4, H3K9, H3K27, H3K36, and H3K79 in facet joint chondrocytes of both IS patients and the control group. Significant upregulation of the level of H3K9me3 was observed, as well as the specific methyltransferase SUV39H1 (Fig. 2). Combining these results with the histopathological findings of the narrowed facet joint cartilage layer but increased number of smaller chondrocytes and expression of COL2A1 (Fig. 1), we further investigated the proliferation of facet joint chondrocytes. Accordingly, PCNA-positive cells and the expression levels of Bcl2 were increased in IS patients (Fig. 3), which may suggest increased proliferative chondrocytes in the facet joint of IS patients.
Previous studies have revealed faster growth of AIS patients during puberty and overgrowth in the anterior compared to the posterior spinal column in AIS patients, which indicated the presence of abnormal regulation and modulation of skeletal growth and endochondral ossification in patients with AIS, although the underlying mechanism is unclear [4, 40, 41]. Since the inferior facet joint pertains to the posterior column of the spine and there was an increase in proliferative facet joint chondrocytes in the IS patients we observed, the abnormal growth of the facet joint may occur in IS patients. This abnormal growth may further play roles in the development and progression of IS. However, generally, the increased expression levels of H3K9me3 and SUV39H1 cannot upregulate the expression of Bcl2 directly because of its function of gene repression. Surprisingly, the expression of the negative regulator of Bcl2, miR-15a, was significantly decreased in the facet joint chondrocytes of IS patients compared with the control group (Fig. 4). Moreover, the increased H3K9me3 levels in the promoter region of the miR-15a host gene, as evidenced by the ChIP assay, suggested that the downregulation of miR-15a may result from the elevation in H3K9me3 levels and SUV39H1 expression (Fig. 4). In parallel, overexpression of SUV39H1 in ATDC5 cells led to the analogous results that the expression of miR-15a was decreased while the H3K9me3 levels in the promoter region of the miR-15a host gene, the expression levels of Bcl2 and Col2a1 were increased, and the cell proliferation was also increased (Figs. 5 and 6).
Therefore, our study showed that SUV39H1 and H3K9me3 can regulate chondrocyte proliferation via miR-15a/Bcl2 signaling, although the reason of the inconsistent expression of miR-15a and its host gene DLEU2 in both IS patients and SUV39H1-overexpressing ATDC5 cells were still unclear. Because studies had shown that intronic miRNA could feedback regulate its host gene, and Morenos et al. had observed that downregulation of DLEU2 by the increased DLEU2 promoter DNA methylation in pediatric acute myeloid leukemia was independent of embedded miR-15a/16-1 [42,43,44]. Thus, we considered that complex mechanisms and processes like feedback regulation, post-transcriptional effects of miRNAs or host genes may be involved, and a large amount of further research may be well needed to figure it out. In addition to chondrocyte proliferation, the progression of endochondral ossification may be related to abnormal spinal growth in IS patients, given the markedly decreased cartilage lacuna and lighter blue color as well as the downregulated expression of COL10A1 in the facet joint cartilage of IS patients compared with the control group (Fig. 1). Nevertheless, the role of aberrant endochondral ossification in the development and progression of IS and the exact underlying mechanism requires further investigation.
In conclusion, as shown in Fig. 7, the increase in H3K9me3, induced by SUV39H1 in chondrocytes of the inferior articular process, activated Bcl2 by inhibiting the expression of miR-15a. However, the dysregulation of the miR-15a/Bcl2 pathway could affect the proliferation of chondrocytes and further lead to abnormal spine growth, which results in the development and progression of IS (Fig. 7). Because of the unclear etiology of IS, prevention-aimed biomechanical intervention rather than effective causality-rooted treatment of IS has been widely applied. Since earlier intervention of IS is a better treatment, further knowledge of the pathogenesis of IS is still urgently needed. In this study, the association of increased H3K9me3 levels at the promoter of miR-15a host gene with the development of IS has been revealed in some way, and it may lay a theoretical foundation for better prevention and treatment of IS in the future.
Patients and controls
The inferior facet joint of 11 IS patients and 10 controls (non-IS patients with lumbar disk herniation and spine fracture), scheduled to be removed during the posterior approach surgery, were collected carefully and harmlessly (Table 1 and Additional file 1: Table S1). The patients and controls were identified based on their clinical manifestations, X-ray, CT, MRI results, etc. The exclusion criteria of the patients and controls were as follows: individuals with severe neuromuscular or genetic disease or routinely using hormones or immune inhibitors.
Histological analysis and pathological evaluation
Fixed and decalcified facet joint tissues were routinely embedded in paraffin and cut into 5 μm thickness. Facet joint sections were deparaffinized and rehydrated in decreasing concentrations of ethanol and then stained with H&E as previously reported . For analysis, pictures of 3–5 different fields per sample were taken under a Leica microscope equipped with a digital CCD (Leica Camera, Wetzlar, Germany), and semi-quantitative analyses of facet joint cartilage histology were performed in a blinded way.
Deparaffinized and rehydrated facet joint sections or fixed cells were blocked by 5% bovine serum albumin with 0.3% trion X-100 at room temperature (RT) for 30 min and then were incubated with primary antibody for overnight at 4 °C. After extensive washing, sections were incubated with fluorescein-labeled anti-rabbit immunoglobulin G (Servicebio, Wuhan, China) for 2 h at RT. Slides or cells were co-stained with DAPI (Servicebio). The primary antibody against PCNA (Catalogue GB11010-1) and Collagen II (Catalogue GB12021) was purchased from Servicebio. For analysis, pictures of 3–5 different fields per sample were taken under a Leica fluorescence microscope equipped with a digital CCD, and semi-quantitative analyses were performed.
Extraction of primary chondrocyte and cell culture
The facet joint tissues of IS patients and control group were washed with phosphate buffer solution (PBS) containing 1% penicillin-streptomycin, and then the cartilage layer was carefully dissected and chopped. The chopped cartilage was digested at 37 °C with 0.25% trypsin for 30 min and with 0.1% collagenase II for 4 h successively. The cell pellet obtained after screening and centrifuging were resuspended with the complete culture medium and cultivated in proper dishes. Toluidine blue staining and Collagen II immunofluorescent staining were used for chondrocytes identification (Additional file 2: Figure S5), and the expanded cultured P3 cells were used for subsequent experiments. Primary chondrocytes and ATDC5 cell lines were cultured in high glucose DMEM (Hyclone, Logan, USA) supplemented with 10% fetal bovine serum and 1% penicillin-streptomycin at 37 °C/5% CO2.
To analyze the ATDC5 proliferation after transfection, EdU staining was conducted using the BeyoClick™ EdU Cell Proliferation Kit with Alexa Fluor 594 (Beyotime Biotechnology, Jiangsu, China, Cat. No. C00788S) according to the manufacturer’s instructions. EdU was added into the medium until the final concentration was 10 μM, and the cells were incubated for 2 hours. After the incubation, the cells were washed with PBS to remove the DMEM and the free EdU probe and then fixed in 4% paraformaldehyde at RT for 15 min. After being co-stained with DAPI, the cells were observed under a Leica fluorescence microscope. EdU positive cells were counted at 5 different areas per sample and reported as a percentage of EdU-positive cells.
Primary chondrocytes were sonicated in RIPA lysis buffer (Beyotime Biotechnology, Jiangsu, China) to obtain whole-cell lysates. Western blot analysis was performed as previously reported and detected by ChemiDoc MP Imager (Bio-Rad Laboratories, Inc., Hercules, CA, USA) . Antibodies for glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (Catalogue 10494-1-AP), Bcl2 (Catalogue 12789-1-AP), and Bax (Catalogue 50599-2-Ig) were obtained from Proteintech Group, Inc, (Chicago, IL, USA), while the remaining antibodies were purchased from Cell Signaling Technologies (Danvers, MA, USA) (Additional file 1: Table S2). The expression levels of target proteins were quantified with the Quantity One 1-D Analysis Software (Bio-Rad) and normalized to the GAPDH in the same sample.
Total cell RNA was extracted using Trizon (KWBio, Beijing, China). For mRNA analysis, equal nanograms of RNA of each sample were reverse transcribed into cDNA using HiFiScript cDNA Synthesis Kit (KWBio). The PCR was performed by NovoStart® Probe qPCR Super Mix (Novoprotein, Shanghai, China) and the transcriptional levels were quantified using GAPDH in the same sample as internal control. For miRNA analysis, equal mass of RNA of each sample was reverse transcribed into cDNA using All-in-One miRNA First-Strand cDNA Synthesis Kit and the real-time PCR was performed with specific primers for the targeted miRNAs using miRNA qPCR Kit (GeneCopoeia, Rockville, MD, USA). The expression levels were quantified using RNU6-B in the same sample as internal control. And the relative fold change was calculated using the 2−ΔΔCT. All of the sequences of primers used are described in Tables 2 and 3 and Additional file 1: Table S3.
ChIP experiment was performed as previously reported  with some modifications. Briefly, primary chondrocytes or ATDC5 cells were fixed with 1% formaldehyde and stopped by 0.125 M glycine. Pellets were resuspended in WB/IP lysis buffer (Beyotime Biotechnology). DNA was sheared into 200–1000 bp in length, and an aliquot was saved as input. Ten micrograms of DNA from each sample were diluted in WB/IP lysis buffer containing H3K9me3 antibody (Catalogue 13969, CST, Danvers, MA, USA) or IgG (Beyotime Biotechnology) plus salmon sperm DNA-protein A+G agarose (Beyotime Biotechnology) for overnight incubation. The eluted DNA-protein cross-links were reversed by proteinase K incubating at 65 °C. DNA was extracted and PCR analyses were performed on miR-15a host gene DLEU2 promoter regions. Primer sequences of ChIP-PCR for human and mouse were shown in Tables 4 and 5, respectively, and the location of primers was indicated in Additional file 2: Figure S6. The input samples were used as an internal control for comparison between samples.
Plasmids and transfection
The CDS sequence of human SUVS9H1 was amplified by Phanta HS Super-Fidelity DNA Polymerase (Vazyme Biotech Co., Nanjing, China), and inserted into the pcDNA3.1+ vector by NheI/BamHI (New England Biolabs Inc., Ipswich, MA, USA) digestion. The sequenced positive plasmids were transfected into ATDC5 cells by lipofectamine 3000 (Invitrogen, Carlsbad, USA) following the user guide. The CDS amplifying primers were listed in Table 6.
All results were expressed as mean ± SD (standard deviation). Statistical significance was determined by analyzing the data of variance by the Student’s t test, except for that the gender difference was analyzed by x2 test. Differences were considered statistically significant when p < 0.05.
Availability of data and materials
All data generated during this study are included in this published article and its supplementary information files.
Adolescent idiopathic scoliosis
BCL2 associated x
Collagen type X
Collagen type II
Deleted in lymphocytic leukemia 2
Disruptor of telomeric silencing 1-like
Enhancer of zeste homolog 2
Histone H3 lysine 79
Lysine (K)-specific demethylases 4 A-C
Long noncoding RNAs
Phosphate buffer solution
Proliferating cell nuclear antigen
Protein arginine N-methyltransferase-4
SET domain bifurcated-1
Sex-determining region Y-Box 9
Suppressor of variegation 3-9, drosophila homolog of 1/2
Konieczny MR, Senyurt H, Krauspe R. Epidemiology of adolescent idiopathic scoliosis. Journal of children’s orthopaedics. 2013;7(1):3–9.
Fadzan M, Bettany-Saltikov J. Etiological theories of adolescent idiopathic scoliosis: past and present. The open orthopaedics journal. 2017;11:1466–89.
Roaf R. The basic anatomy of scoliosis. The Journal of bone and joint surgery British volume. 1966;48(4):786–92.
Guo X, Chau WW, Chan YL, Cheng JC. Relative anterior spinal overgrowth in adolescent idiopathic scoliosis. Results of disproportionate endochondral-membranous bone growth. The Journal of bone and joint surgery British volume. 2003;85(7):1026–31.
Shi B, Mao S, Liu Z, Sun X, Zhu Z, Zhu F, et al. Spinal growth velocity versus height velocity in predicting curve progression in peri-pubertal girls with idiopathic scoliosis. BMC musculoskeletal disorders. 2016;17(1):368.
Binns M. Joint laxity in idiopathic adolescent scoliosis. The Journal of bone and joint surgery British volume. 1988;70(3):420–2.
Czaprowski D, Kotwicki T, Pawlowska P, Stolinski L. Joint hypermobility in children with idiopathic scoliosis: SOSORT award 2011 winner. Scoliosis. 2011;6:22.
Pal GP, Bhatt RH, Patel VS. Mechanism of production of experimental scoliosis in rabbits. Spine. 1991;16(2):137–42.
Burwell RG, Dangerfield PH, Moulton A, Grivas TB. Adolescent idiopathic scoliosis (AIS), environment, exposome and epigenetics: a molecular perspective of postnatal normal spinal growth and the etiopathogenesis of AIS with consideration of a network approach and possible implications for medical therapy. Scoliosis. 2011;6(1):26.
Lafon-Hughes L, Di Tomaso MV, Mendez-Acuna L, Martinez-Lopez W. Chromatin-remodelling mechanisms in cancer. Mutation research. 2008;658(3):191–214.
Wolffe AP, Matzke MA. Epigenetics: regulation through repression. Science. 1999;286(5439):481–6.
Brunet A, Berger SL. Epigenetics of aging and aging-related disease. The journals of gerontology Series A, Biological sciences and medical sciences. 2014;69(Suppl 1):S17–20.
Hardy T, Mann DA. Epigenetics in liver disease: from biology to therapeutics. Gut. 2016;65(11):1895–905.
Liu XY, Wang L, Yu B, Zhuang QY, Wang YP. Expression signatures of long noncoding RNAs in adolescent idiopathic scoliosis. BioMed research international. 2015;2015:276049.
Garcia-Gimenez JL, Rubio-Belmar PA, Peiro-Chova L, Hervas D, Gonzalez-Rodriguez D, Ibanez-Cabellos JS, et al. Circulating miRNAs as diagnostic biomarkers for adolescent idiopathic scoliosis. Scientific reports. 2018;8(1):2646.
Ogura Y, Kou I, Takahashi Y, Takeda K, Minami S, Kawakami N, et al. A functional variant in MIR4300HG, the host gene of microRNA MIR4300 is associated with progression of adolescent idiopathic scoliosis. Human molecular genetics. 2017;26(20):4086–92.
Mao SH, Qian BP, Shi B, Zhu ZZ, Qiu Y. Quantitative evaluation of the relationship between COMP promoter methylation and the susceptibility and curve progression of adolescent idiopathic scoliosis. Eur Spine J. 2018;27(2):272–7.
Mao S, Xu L, Zhu Z, Qian B, Qiao J, Yi L, et al. Association between genetic determinants of peak height velocity during puberty and predisposition to adolescent idiopathic scoliosis. Spine. 2013;38(12):1034–9.
Carballo CB, Nakagawa Y, Sekiya I, Rodeo SA. Basic Science of Articular Cartilage. Clinics in sports medicine. 2017;36(3):413–25.
O’Leary SA, Link JM, Klineberg EO, Hu JC, Athanasiou KA. Characterization of facet joint cartilage properties in the human and interspecies comparisons. Acta biomaterialia. 2017;54:367–76.
Shen G. The role of type X collagen in facilitating and regulating endochondral ossification of articular cartilage. Orthodontics & craniofacial research. 2005;8(1):11–7.
Schneevoigt J, Fabian C, Leovsky C, Seeger J, Bahramsoltani M. In vitro expression of the extracellular matrix components aggrecan, collagen types I and II by articular aartilage-derived chondrocytes. Anatomia, histologia, embryologia. 2017;46(1):43–50.
Allas L, Boumediene K, Bauge C. Epigenetic dynamic during endochondral ossification and articular cartilage development. Bone. 2018.
Kouzarides T. Histone methylation in transcriptional control. Current opinion in genetics & development. 2002;12(2):198–209.
Pekarsky Y, Balatti V, Croce CM. BCL2 and miR-15/16: from gene discovery to treatment. Cell death and differentiation. 2018;25(1):21–6.
Chen HX, Xu XX, Tan BZ, Zhang Z, Zhou XD. MicroRNA-29b inhibits angiogenesis by targeting VEGFA through the MAPK/ERK and PI3K/Akt signaling pathways in endometrial carcinoma. Cellular physiology and biochemistry : international journal of experimental cellular physiology, biochemistry, and pharmacology. 2017;41(3):933–46.
Xia L, Zhang D, Du R, Pan Y, Zhao L, Sun S, et al. miR-15b and miR-16 modulate multidrug resistance by targeting BCL2 in human gastric cancer cells. International journal of cancer. 2008;123(2):372–9.
Inoue A, Mizushima T, Wu X, Okuzaki D, Kambara N, Ishikawa S, et al. A miR-29b byproduct sequence exhibits potent tumor-suppressive activities via inhibition of NF-kappaB signaling in KRAS-mutant colon cancer cells. Molecular cancer therapeutics. 2018;17(5):977–87.
Garzon R, Pichiorri F, Palumbo T, Visentini M, Aqeilan R, Cimmino A, et al. MicroRNA gene expression during retinoic acid-induced differentiation of human acute promyelocytic leukemia. Oncogene. 2007;26(28):4148–57.
Klein U, Lia M, Crespo M, Siegel R, Shen Q, Mo T, et al. The DLEU2/miR-15a/16-1 cluster controls B cell proliferation and its deletion leads to chronic lymphocytic leukemia. Cancer cell. 2010;17(1):28–40.
Yao Y, Wang Y. ATDC5: an excellent in vitro model cell line for skeletal development. Journal of cellular biochemistry. 2013;114(6):1223–9.
Watanabe K, Michikawa T, Yonezawa I, Takaso M, Minami S, Soshi S, et al. Physical activities and lifestyle factors related to adolescent idiopathic scoliosis. The Journal of bone and joint surgery American volume. 2017;99(4):284–94.
McMaster ME, Lee AJ, Burwell RG. Physical activities of patients with adolescent idiopathic scoliosis (AIS): preliminary longitudinal case-control study historical evaluation of possible risk factors. Scoliosis. 2015;10:6.
Fraga MF, Ballestar E, Paz MF, Ropero S, Setien F, Ballestar ML, et al. Epigenetic differences arise during the lifetime of monozygotic twins. Proceedings of the National Academy of Sciences of the United States of America. 2005;102(30):10604–9.
Ito T, Yadav N, Lee J, Furumatsu T, Yamashita S, Yoshida K, et al. Arginine methyltransferase CARM1/PRMT4 regulates endochondral ossification. BMC developmental biology. 2009;9:47.
Schwarz D, Varum S, Zemke M, Scholer A, Baggiolini A, Draganova K, et al. Ezh2 is required for neural crest-derived cartilage and bone formation. Development. 2014;141(4):867–77.
Ideno H, Shimada A, Imaizumi K, Kimura H, Abe M, Nakashima K, et al. Predominant expression of H3K9 methyltransferases in prehypertrophic and hypertrophic chondrocytes during mouse growth plate cartilage development. Gene expression patterns : GEP. 2013;13(3-4):84–90.
Czvitkovich S, Sauer S, Peters AH, Deiner E, Wolf A, Laible G, et al. Over-expression of the SUV39H1 histone methyltransferase induces altered proliferation and differentiation in transgenic mice. Mechanisms of development. 2001;107(1-2):141–53.
Castano Betancourt MC, Cailotto F, Kerkhof HJ, Cornelis FM, Doherty SA, Hart DJ, et al. Genome-wide association and functional studies identify the DOT1L gene to be involved in cartilage thickness and hip osteoarthritis. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(21):8218–23.
Little DG, Song KM, Katz D, Herring JA. Relationship of peak height velocity to other maturity indicators in idiopathic scoliosis in girls. The Journal of bone and joint surgery American volume. 2000;82(5):685–93.
Wang WW, Man GC, Wong JH, Ng TB, Lee KM, Ng BK, et al. Abnormal response of the proliferation and differentiation of growth plate chondrocytes to melatonin in adolescent idiopathic scoliosis. International journal of molecular sciences. 2014;15(9):17100–14.
Petrocca F, Visone R, Onelli MR, Shah MH, Nicoloso MS, de Martino I, et al. E2F1-regulated microRNAs impair TGFbeta-dependent cell-cycle arrest and apoptosis in gastric cancer. Cancer cell. 2008;13(3):272–86.
Cao G, Huang B, Liu Z, Zhang J, Xu H, Xia W, et al. Intronic miR-301 feedback regulates its host gene, ska2, in A549 cells by targeting MEOX2 to affect ERK/CREB pathways. Biochemical and biophysical research communications. 2010;396(4):978–82.
Morenos L, Chatterton Z, Ng JL, Halemba MS, Parkinson-Bates M, Mechinaud F, et al. Hypermethylation and down-regulation of DLEU2 in paediatric acute myeloid leukaemia independent of embedded tumour suppressor miR-15a/16-1. Molecular cancer. 2014;13:123.
Li J, Huang J, Li JS, Chen H, Huang K, Zheng L. Accumulation of endoplasmic reticulum stress and lipogenesis in the liver through generational effects of high fat diets. Journal of hepatology. 2012;56(4):900–7.
Wang YJ, Zhang HQ, Han HL, Zou YY, Gao QL, Yang GT. Taxifolin enhances osteogenic differentiation of human bone marrow mesenchymal stem cells partially via NF-kappaB pathway. Biochemical and biophysical research communications. 2017;490(1):36–43.
The authors thank all the staff of the Department of Spine Surgery, Xiangya Hospital, Central South University for their dedicated assistance in patient sample collection.
This work was supported by the National Natural Science Foundation of China (No. 81772298) and the National Natural Science Foundation of Hunan Province, China (No. 2018JJ3848).
Ethics approval and consent to participate
This study has been conducted adhering to the principles of the Declaration of Helsinki II and approved by the medical ethics committee of Xiangya Hospital, Central South University (ethical code, 201703358). Written informed consent was obtained from all the study participants or their parents or guardians.
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. Clinical data of study subjects. Table S2. List of antibodies from Cell Signaling Technologies. Table S3. Primers sequences for genes. (DOCX 20 kb)
Figure S1. No significant changes of other H3K9 methyltransferases and demethylases between facet joint chondrocytes of IS patients and controls. Figure S2. ChIP-PCR data of IS patients and controls. Figure S3. SUV39H1 is highly conserved in different species. Figure S4. ChIP-PCR data of SUV39H1-overexpressing and control ATDC5 cells. Figure S5. Identification of extracted human facet joint chondrocytes. Figure S6. Characterization of primers for ChIP-PCR. (DOCX 3899 kb)