- Research article
- Open Access
Integrative analysis of genome-wide DNA methylation and single-nucleotide polymorphism identified ACSM5 as a suppressor of lumbar ligamentum flavum hypertrophy
Arthritis Research & Therapy volume 23, Article number: 251 (2021)
Hypertrophy of ligamentum flavum (HLF) is a common lumbar degeneration disease (LDD) with typical symptoms of low back pain and limb numbness owing to an abnormal pressure on spinal nerves. Previous studies revealed HLF might be caused by fibrosis, inflammatory, and other bio-pathways. However, a global analysis of HLF is needed severely.
A genome-wide DNA methylation and single-nucleotide polymorphism analysis were performed from five LDD patients with HLF and five LDD patients without HLF. Comprehensive integrated analysis was performed using bioinformatics analysis and the validated experiments including Sanger sequencing, methylation-specific PCR, qPCR and ROC analysis. Furthermore, the function of novel genes in ligamentum flavum cells (LFCs) was detected to explore the molecular mechanism in HLF through knock down experiment, overexpression experiment, CCK8 assay, apoptosis assay, and so on.
We identified 69 SNP genes and 735 661 differentially methylated sites that were enriched in extracellular matrix, inflammatory, and cell proliferation. A comprehensive analysis demonstrated key genes in regulating the development of HLF including ACSM5. Furthermore, the hypermethylation of ACSM5 that was mediated by DNMT1 led to downregulation of ACSM5 expression, promoted the proliferation and fibrosis, and inhibited the apoptosis of LFCs.
This study revealed that DNMT1/ACSM5 signaling could enhance HLF properties in vitro as a potential therapeutic strategy for HLF.
Lumbar degeneration diseases (LDDs) are common spinal lesions with typical symptoms of low back pain and limb numbness owing to an abnormal pressure on spinal nerves . LDDs limit the patients’ mobility immensely and decrease their life quality significantly. It is worth noting that LDDs mostly occur in the elderly population, which brings high pressure on the whole society, especially with the rapid population aging in the major countries. LDDs are mainly composed of lumbar spinal canal stenosis (LSCS), degenerative lumbar spondylolisthesis (DLS), lumbar disc herniation (LDH), of which LSCS is most familiar with a prevalence of 47% in patients with legs pain and numbness. Among the causative factors, ligamentum flavum hypertrophy (HLF) is thought to be the most critical factor in the pathogenesis of LSCS [2, 3]. However, the ligamenta flava in LSCS patients are often 4-8 mm, compared with less than 4 mm in normal ligamenta flava .
Ligamenta flava consist of 80% elastic fibers and 20% collagen fibers , which makes them able to connect the adjacent vertebrae, keep vertebrae’ mobility, and protect the spinal cord . As we know, lumbar suffers almost the whole weight of the upper body. Similar to that, lumbar ligamenta flava also bear the highest pressure among all the ligamenta flava in the vertebral column. Because of the high loads, lumbar ligamenta flava have a significantly higher frequency of becoming injured, especially the lowest L4/5 lumbar ligamentum flavum . The injuries will further lead to fibrotic changes with a loss of elastic fibers and a cumulation of collagen fibers . And it has been reported that fibrosis may be the primary pathology of HLF .
Previous studies have demonstrated some molecular mechanisms underlying HLF, for example, inflammatory pathways including TGF-β1, COX-2/PGE2, and iNOS; JNK pathway; p38 pathway; and angiogenesis processes including Angptl2 and VEGF [10,11,12,13]. However, the mechanisms behind HLF are still poorly understood and lack global screens. As we know, mechanical stress may be the primary causative factor for HLF, and mechanical stress has a great relationship with the habits and the environment. Furthermore, the critical influence environment on the epigenetics has been well established and thus affects the gene expression profoundly . So, we hypothesize epigenetics may play an essential role in HLF.
Here, we investigated the epigenetic changes in the LDD patients with or without HLF (5 patients in each group) by the use of genome-wide DNA methylation and single-nucleotide polymorphism (SNP) analysis to widen our understanding of HLF in a global field. We identify the methylation levels of COL4A3, HSPA6, PIGC, and ACSM5 had positive correlation with HLF. Then, we found DNMT1-mediated hypermethylation of ACSM5 promoter leads to downregulation of ACSM5 expression in LFCs from HLF patients, and ACSM5 inhibits the proliferation and fibrosis and promotes the apoptosis of LFCs from HLF patients. Finally, we demonstrated that inhibition of DNMT1 suppressed HLF properties via regulation of ACSM5 expression.
Materials and methods
Human ligamentum flavum sample collection
Human ligamentum flavum samples in this study were collected according to protocols approved by the Ethical Committee of the Zhujiang Hospital of Southern Medical University (Guangzhou, China). Informed consents were obtained from all subjects. Ligamentum flavum samples were isolated from the anatomical region (L4/5) of LDD patients. Ligamentum flavum samples acquired from LDD patients with ligamentum flavum thickening (≤ 4 mm) were considered to be the non-ligamentum flavum hypertrophy (non-HLF) group, whereas the pathological ligamentum flavum samples harvested from LDD patients with ligamentum flavum thickening (> 4 mm) were assigned to the HLF group. The patients’ clinical characters are shown in Table 1.
Ligamentum flavum samples and cultured ligamentum flavum cells were used for DNA extraction with a QIAamp DNA Midi Kit (Qiagen). Extracted DNA was digested with RNase (Thermo Scientific) and purified with a Genomic DNA Clean Kit (Zymo Research).
Whole-exome sequencing (WES)
WES was performed as previously described . Briefly, Agilent SureSelect All Exon V6 (Agilent Technologies) with reagents were used for the synthesis of the library. Sequencing was performed on an Illumina HiSeq 2500 platform (Illumina).
Single-nucleotide polymorphism (SNP) detection and analysis
The raw data of WES were performed QC process with fastqc and removed adapters with Cutadapter 2.10. Then, the clean data were mapped to the human reference genome GRCh37/hg19 by BWA-MEM. We used SAMtools for SNP detection and ANNOVAR for SNP annotation.
Infinium Methylation EPIC BeadChip assay
The methylation assays were performed the same as described previously . Genomic DNA was sulfurated with the EZ DNA Methylation kit (Zymo Research). Genome-wide DNA methylation profiles were generated by Illumina’s Infinium MethylationEPIC BeadChip assay (EPIC array) (Illumina), determining > 850,000 methylated sites. The assay was scanned on an Illumina HiScan.
EPIC array data processing
Raw EPIC array data were performed with CHAMP, an R package. Filtered data were used for evaluating the methylation level of each methylated site with a β value. A site with a β value > 0.6 was considered an up-methylated site, whereas a site with a β value < 0.2 was considered a down-methylated site. Then, the data were normalized using Illumina’s default normalization method. Normalized data were used for the detection of differentially methylated sites between HLF group and CT group with an R package “limma.” In detail, a linear model was established to calculate the P value of each site, and the P value was subsequently adjusted with a Benjamini & Hochberg method. A site with a P value < 0.05 was considered as a differentially methylated site.
The Volcano plot, heatmap, Manhattan plot, and pie chart to illustrate the differentially methylated sites were drawn with an R package “ggplot2”.
GO and KEGG pathway enrichment analysis
Differential genes were mapped to the GO database or KEGG database, and the enrichment analysis of GO term or KEGG pathway was performed using the function Clusterprofiler in R package . The results are shown in figures drawn with ggplot2 or Cytoscape 3.7.2.
Protein-protein interaction (PPI) analysis
PPI analysis was performed using STRING database (https://string-db.org/). The network was established with a parameter of score > 900. Established gene network was plotted using Cytoscape 3.7.2.
Validation of the SNPs
The SNPs to be validated were chosen based on the WES data. The primers used for PCR were listed in Table 2. The PCR products were sequenced using Sanger method by Invitrogen.
Methylation-specific analysis (MSP)
The methylation level of ACSM5, COL4A3, PIGC, and HSPA6 were validated using methylation-specific PCR [21, 22]. A methylation-independent assay with non-CpG including primers for the β-actin gene (ACTB) was used in order to verify DNA quality and to normalize results. Specificity and cross-reactivity of methylated and unmethylated primers (Table 3) were evaluated by using unconverted gDNA, SB-converted methylated and non-methylated DNA standards. Analytical sensitivity of qMSP assays was evaluated by using serial dilutions of SB-converted methylated and nonmethylated DNA standards and was found to be 0.1%. The assay efficiency (expressed as E = 10− 1/slope − 1) was evaluated by using serial dilutions of the SB-converted methylated DNA standards in H2O and was in the range of 91–105%. The analysis was performed according to the RQ sample (Relative Quantification) = 2−ΔΔCT method. Specifically, ΔΔCT values were generated for each target after normalization by ACTB values and using 1% methylation as calibrator and then were multiplied by 100 (RQ = 2−ΔΔCT× 100).
Correlation analysis and ROC analysis
The correlations of the key genes hypermethylation and flavonoid ligament hypertrophy was analyzed by Pearson correlation. The coefficient of association and P value were performed to confirm their association. The ROC curves were analyzed by pROC package in R based on the gene methylation level and the patients’ clinical characters.
Isolation and culture of ligamentum flavum cells
To extract of human ligamentum flavum cell, ligamentum flavum samples were isolated from LDD patients with or without ligamentum flavum hypertrophy. Then, the samples were minced and digested with 0.2% type I collagenase 2 for 1 h at 37 °C. Digested tissues were seeded in 6-cm dishes and cultured in DMEM medium (Lonza) containing 10% fetal bovine serum (Gibco) and 1% penicillin/streptomycin (Invitrogen) for 7 days. Migrated cells were digested and pass on for subsequent experiments.
The expression levels of collagen I (Abcam, ab260043), collagen III (Abcam, ab184993), vimentin (Abcam, ab16700), and fibronectin (Abcam, ab2413) were detected by immunofluorescence staining and observed under a Leica Microsystems DMI4000 fluorescence microscope with a DFC365FX camera (Leica Microsystems) to identify the cultured ligamentum flavum cells (LFCs).
Transfection and cell treatment
The pCMV-DNMT1, pCMV-ACSM5, and pCMV were purchased from Hanbio Biotechnology. sh-DNMT1, sh-ACSM5, and sh-NC were purchased from Thermofisher Biotechnology. Transfection of plasmid and shRNAs was performed using Lipo2000 (Thermofisher) according to the protocol of the manufacture. The sequences for sh-DNMTs and sh-ACSM5 are shown in Table 4.
For the treatment of 5-AzaC, 0 μM, 1 μM, 5 μM, and 10 μM 5-AzaC were added to LFCs from HLF tissues for 24 h, respectively. Then, the cells were used for other experiments.
Cells (5000 cells/well in 96-well plates) were adjusted to different treatments for 24–72 h at 37 °C. Then, 10 μl CCK-8 reagent (Abcam) was added to each well and cultured for another 1 h at 37 °C. The absorption values were measured by a microplate reader (Thermofisher) and were used to calculate the cell viability.
Apoptosis assays and the test of caspase 3 activity
Cells (5 × 105 cells/well in 6-well plates) were adjusted to different treatments for 24 h at 37 °C. Then, the cells were digested and stained with an Annexin V-FITC/propidine iodide (PI) double-staining kit (Abcam) according to the manufacturer’s instructions. The cells were then analyzed with FACS using CytoFlex (Beckman).
Caspase 3 activity was measured with a caspase 3 assay kit (Abcam) according to the manufacturer’s instructions.
Total RNA was extracted with TRIzol (Invitrogen) according to the manufacturer’s protocol. The synthesis of cDNA was used for the mRNA Reverse Transcription Kit (Qiagen). Amplification reactions was performed with SYBR Green PCR Master Mix (Applied Biosystems) using a CFX96 Real-Time System thermocycler (Bio-Rad) and normalized to GAPDH. The relative expression of genes were calculated using the 2−△△Ct method. All the primers’ sequences were listed in Table 4.
Western Blot analysis
The procedure of western blot was the same as the previous paper . In brief, the total protein from the cells was lysed with radio immunoprecipitation assay buffer (Solarbio Biotech). Protein concentrations were measured using BCA protein assay kit (Beyotime) according to the manufacturer’s instructions. Then, equivalent amount of protein were separated by SDS-PAGE and transferred into polyvinylidene fluoride membranes (Bio-Rad). After blocking with 5% non-fat milk for 1 h at room temperature, the membranes were incubated with primary antibodies at 4 °C overnight and then incubated with secondary antibodies for 2 h at room temperature. Bands were visualized with an ECL detection kit (Thermofisher) using ChemiDoc XRS Plus (BioRad) and analyzed using the Image J software. All primary antibodies and secondary antibodies in the study are listed as follows: anti-ACSM5 antibody (1:500, Proteintech), anti-collagen I (1:1000, Abcam), anti-collagen III (1:1000, Abcam), anti-MMP2 (1:1000, Proteintech), anti-MMP13 (1:1000, Proteintech), and anti-β-actin antibody (1:5000, Abcam).
Data are shown as the mean ± S.E.M. Statistical analyses of data were performed using GraphPad 6.0 Software. Analyses of significant differences between groups were performed using Student’s t tests or ANOVA with Tukey’s test.
SNP and DNA methylation profiles in HLF
We collected five L4/5 ligamentum flavum samples from LDD patients with HLF (HLF group) and five samples in the same location from LDD patients without HLF (CT group). The detailed clinical characters of the patients are shown in Table 1. Each sample was separated into two and performed whole-exome sequencing (WES) and Illumina Infinium MethylationEPIC BeadChips assays (EPIC arrays).
We used SAMtools to identify the SNPs and ANNOVAR to annotate the SNPs from the WES data. Then, we evaluated the SNPs using SIFT (> 0.05), CADD (> 15), and Plyphen2_HVAR (< 0.909), and we got 69 SNP genes with significance compared to the reference genome sequence (Fig. 1A). We found broad missense SNPs among all the 69 genes, indicating they might probably participate in the development of HLF. After quality control and filtering, a total of 735,661 differentially methylated sites were identified, of which 415,237 were up-methylated sites and 320,424 were down-methylated (Supplementary Fig. 1). Heatmap of the top 50 differentially methylated sites shown a significant preference with up-methylated sites (45 up-methylated sites, contrary to 5 down-methylated sites; Fig. 1B). The differentially methylated sites were equally distributed among the autosomes (Fig. 1C). However, most of the differentially methylated sites were located at the gene body, with a prevalence of 40.88%. Furthermore, we have noticed that about 20% of the differentially methylated sites were at the promoter region, of which 9.87% were at the area 1500 bp up to the transcription start site (TSS) and 7.93% were at the 5' UTR (Fig. 1D). Methylation changes at either the gene body or promoter may influence the gene expression immensely. Therefore, the majority of the differentially methylated sites between the HLF group and the CT group might play an important role in regulating the signal pathways involved in the HLF.
To verify which bio-functions or pathways might be regulated by the differentially methylated sites, we performed GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis. The GO terms consist of three themes: biology process (BP), cell component (CC), and molecular function (MF). The top 15 enriched terms in each analysis are shown in Supplementary Fig. 2. In summary, the bio-functions or pathways that were influenced by the differentially methylated sites were enriched in ECM-related terms, inflammation-related terms, and cell proliferation-related pathways including wnt pathway, Akt pathway, Ras pathway, and MAPK pathway. These findings are consistent with previous studies, like ECM, inflammation, and cell proliferation, which are thought to be the main causative factors for HLF [10,11,12,13].
Integrated analysis of SNP and methylation data
To investigate the HLF-causative genes, we combined the SNP and methylation data for integrated analysis. The 7334 up-methylated genes or the 2864 down-methylated genes were taken intersection with the 69 SNP genes, respectively. Therefore, we got 30 up-methylated SNP genes (Fig. 2A, C) and eight down-methylated SNP genes (Fig. 2B, D). Further functional categories enrichment analysis revealed that the 38 intersection genes were mainly enriched to sequence variant and Fibronectin type-III related pathways (Fig. 2E).
To ulteriorly verify the key genes that participate in the process of HLF in the 38 filtered genes, we performed protein-protein interaction (PPI) analysis using STRING database with a filtered threshold of score > 900. The result indicated a network of following genes: CLSPN, COL12A1, GPRIN2, IDO2, MUC16, MYH4, NEB, OBSCN, MKI67, COL4A3, HYDIN, INMT, and TTN (Fig. 2F). Meanwhile, we mapped the 38 candidates to the KEGG database and filtered out five genes: COL4A3, MACC1, HSPA6, PIGC, and ACSM5. These five genes were all up-methylated genes and mainly participated in glucose and lipid metabolism, PI3K pathway, collagen related pathways, and so on (Fig. 2C, D, G).
Combined those results above, we propose that COL4A3, MACC1, HSPA6, PIGC, and ACSM5 might be the key genes in the process of HLF. And the relationship of those genes with HLF has never been reported.
Methylation and SNP verification of the key genes in clinical samples
We collected ten L4/5 ligamentum flavum samples from LDD patients with HLF (HLF group) and ten samples in the same location from LDD patients without HLF (non-HLF group). The clinical characters of the patients were shown in Table 1.
To verify the SNP sites of the five genes, we extracted DNA from the patients' samples and performed DNA sequencing with the Sanger method. For COL4A3 (rs10205042), MACC1 (rs975263), and PIGC (rs2230471 and rs8062344), the verification data were the same with the WES data (Fig. 3C–F). For ACSM5 (rs 8062344) and COL4A3 (rs6436669), only the HLF samples had the same SNPs with our WES data, while the non-HLF samples were the same with the reference genome (Fig. 3A, B). Those data revealed that the majority of verified results were consistent with our WES data.
To validate the methylation level, we performed methylation-specific PCR (MSP). ACSM5, COL4A3, PIGC, and HSPA6 showed immensely increased methylation levels in the HLF group compared to the non-HLF group, consistent with our genome-wide DNA methylation data (Figs. 1B, 2c, 4). However, MSP results showed no methylation difference of MACC1 between the HLF group and the non-HLF group (Fig. 4).
Then, we explored the relationship between the genes’ methylation level with the patients’ clinical characters. We investigated the correlation between the methylation levels of COL4A3, HSPA6, PIGC, or ACSM5 with the patients’ ligamentum flavum thickness and performed ROC analysis to evaluate the clinical efficiency and sensitivity of the indicated genes. We found that the methylation levels of COL4A3, HSPA6, PIGC, and ACSM5 were all positively correlated with the patients’ ligamentum flavum thickness (correlation index arrange from − 0.69 to − 0.87; Fig. 5A). And it is expected that the four genes all shown perfect efficiency and sensitivity for HLF with AUC values of 1.00 (Fig. 5B).
Therefore, we validated the SNP status and methylation levels of COL4A3, MACC1, HSPA6, PIGC, and ACSM5. The results were mostly consistent with the data of our genome-wide DNA methylation and SNP analysis. Furthermore, we found that the methylation level of COL4A3, HSPA6, PIGC, and ACSM5 had a perfect sensitivity for HLF, highlighting their therapeutic potential for HLF.
DNA methylation downregulated the expression of ACSM5 in a DNMT1-dependent manner
To explore the therapeutic potential of COL4A3, HSPA6, PIGC, and ACSM5, we performed functional experiments in vitro. We isolated ligamentum flavum cells (LFCs) from ligamentum flavum samples belonged to LDD patients with HLF (HLF group) or without HLF (non-HLF group). LFCs of the 3rd passage were identified through immunofluorescence staining (Fig. 6A, B and Fig.S4) and used for subsequent experiments. The cultured cells had a typical LFCs phenotype and uniformly expressed of LFCs markers (collagen I, collagen III, vimentin or fibronectin) in each cell, suggesting the high purity of the LFCs obtained. The expression of COL4A3, HSPA6, PIGC, and ACSM5 in HLF cells and non-HLF cells was first examined by qPCR. The results showed that HLF, ACSM5 was downregulated, HSPA6 was upregulated, and COL4A3 and PIGC had no change in HLF cells (Fig. 6C). Commonly, a higher level of DNA methylation brings a decrease in the gene expression; however, other factors may diminish or even converse this phenomenon [24, 25], indicating the regulation of COL4A3, HSPA6, and PIGC’s expression might be complex while ACSM5’s expression was mainly influenced by DNA methylation. Then, we validated that ACSM5 had a higher DNA methylation level in HLF group (Fig. 6D). Furthermore, decreasing DNA methylation levels were observed with the treatment of increasing doses of 5-AzaC (Fig. 6E) and liberated the expression of ACSM5 in the HLF group (Fig. 6F), proving the solid relationship between DNA methylation level and the expression of ACSM5.
As we know, DNA methylation is mainly regulated by the DNMT family. To explore the upstream regulators of ACSM5’s DNA methylation, we knocked down DNMT1, DNMT3A, and DNMT3B in HLF cells (Fig. 6G) and found that only knocking down of DNMT1 could rescue the expression level of ACSM5 (Fig. 6H), indicating the methylation of ACSM5 was mainly contributed by DNMT1. Nevertheless, overexpression of DNMT1 could decrease the expression of ACSM5 (Supplementary Fig. 3). The above results indicate that the downregulation of ACSM5 in HLF is related to its promoter hypermethylation, which is regulated by DNMT1.
ACSM5 regulates the proliferation, fibrosis, and apoptosis of LFCs
We over-expressed ACSM5 in LFCs from HLF patients (Fig. 7A) and found that the over-expression of ACSM5 could inhibit the proliferation of LFCs from HLF patients (Fig. 7B). At the same time, ACSM5 could also promote the apoptosis and enhance the activity of caspase3 in LFCs from HLF patients (Fig. 7C, D). Moreover, over-expression of ACSM5 could decrease the fibro features of LFCs from HLF patients, manifesting as decreasing expression of COL1, COL3, MMP13, and MMP2 (Fig. 7E, F). Collectively, these results demonstrated that ACSM5 could inhibit the proliferation and fibrosis and promote the apoptosis of LFCs from HLF patients.
DNMT1 regulates the proliferation, fibrosis, and apoptosis of LFCs via downregulating ACSM5 expression
We have proved that knocking down of DNMT1 could increase the expression of ACSM5; however, knocking down of ACSM5 at the same time could diminish this influence (Fig. 8A). Then, we observed that knocking down DNMT1 could inhibit the proliferation, promote apoptosis, and decrease the fibro features of LFCs from HLF patients (Fig. 8B–F), which were similar to the influences of ACSM5 over-expression. Furthermore, knocking down of ACSM5 could promote the proliferation, inhibit apoptosis, and enhance the fibro features of LFCs from HLF patients; however, knocking down of DNMT1 at the same time could rescue the influence of knocking down of ACSM5 (Fig. 8B–F). Taken together, these data proved that the silencing DNMT1 suppresses the proliferation and fibrosis and promote the apoptosis of LFCs from HLF patients via regulating ACSM5 expression.
With the acceleration of global aging and changes in lifestyle, the incidence of LSCS increases year by year and tends to be younger, which has seriously threatened the health and life of patients . Currently, HLF is thought to be a critical pathological leading to LSCS , but the precise molecular mechanisms of HLF remain unknown so far. Herein, we used the genome-wide DNA methylation sequencing and SNP sequencing for the first time to analyze the epigenetic changes in HLF tissues and non-HLF tissues. Integration analysis of the above two sequencing results would help us to identify the key molecules that regulate the occurrence and development of HLF. Fortunately, we found that five key genes may be involved in the occurrence and development of HLF, among which ACSM5 was hypermethylated and low expressed in HLF tissues and cells.
Previous studies have raised the importance of SNP to LDDs , but a genome-wide screen is lacking in HLF tissues in patients with LSCS. Our research revealed that the HLF tissues shared most SNPs with the non-HLF tissues, but there are also some different SNPs, which might be involved in the occurrence and development of HLF. In addition, DNA methylation is an important epigenetic mechanism. Studies have reported abnormal DNA methylation in ossified spinal LF tissues and calcified thoracic LF tissues [29, 30], but the differential expression profile of DNA methylation between HLF tissues and non-HLF tissues has not been reported. In this study, the whole genome epigenetic landscape of HLF and non-HLF was studied for the first time. We found that there were a large number of differential DNA methylation genes between HLF tissues and non-HLF tissues, and these differential DNA methylation genes were mainly enriched in MAPK, WNT, and PI3K/AKT signaling pathways, suggesting that these signaling pathways played an important role in the occurrence and development of HLF. Our findings were consistent with those of previous studies . For example, Cao et al. reported that p38 MAPK signaling pathway promoted the homeostasis of the extracellular matrix and the hypertrophy of LF . Zheng et al. indicated that the inhibitor of MAPK/ERK signaling pathway restrained TGF-1β-induced LF fibrosis, suggesting a crucial role of MAPK/ERK signaling pathway in HLF . The involvement of WNT signaling pathway in HLF has been reported by Taiki et al. . Moreover, Chen et al. demonstrated that estradiol upregulated the expression of MMP-13 via activating PI3K pathway and contributed to the homeostasis of ECM in the LF . Zhou et al. revealed that LPA bind to the receptor LPAR1 to induce LF cells proliferation and anti-apoptosis through activating AKT signaling pathway .
We further integration analysis of DNA methylation sequencing data and SNP sequencing data to identify five key genes: COL4A3, HSPA6, PIGC, ACSM5, and MACC1. Validation of the SNP sites indicated that only ACSM5 rs80623344 and COL4A3 rs6436669 mutated in HLF tissues compared with non-HLF tissues. COL4A3 gene variants have been reported to be involved in the occurrence and development of various diseases, such as diabetic nephropathy , focal segmental glomerulosclerosis , and Alport syndrome [38, 39]. Here, it is the first time that COL4A3 gene mutation has been identified in HLF. In addition, ACSM5, as a member of acyl-CoA synthetase family, is located in mitochondrial matrix and involved in the first step of fatty acid metabolism, but physiological characteristics of ACSM5 have been seriously neglected by human beings . Recently, the role of ACSM5 in human diseases has gradually attracted attention, and preliminary exploratory studies have been conducted in some fields. For example, studies have shown that ACSM5 is low expressed in tumor tissues of breast cancer and lung cancer and its expression level is related to the poor prognosis of patients, which can be used as a potential marker for disease prognosis and diagnosis [41, 42]. Qi et al. have found that ACSM5 is closely related to attention-deficit/hyperactivity disorder (ADHD), suggesting that ACSM5 may be a key target in the prevention and treatment of ADHD . Cao et al. reported that ACSM5 was low expressed in renal tissue of rats with chronic glomerulonephritis (CGN), while Qiteng Xiaozhao granules significantly increased the low expression of ACSM5 in the renal tissue of rats with CGN, suggesting that ACSM5 was involved in the occurrence and development of CGN . Cao et al. found that ACSM5 gene was hypermethylated in the peripheral blood of patients with myocardial infarction, and its methylation site and level might be closely related to the occurrence of myocardial infarction . However, the role of ACSM5 in the HLF has not been reported yet. In this study, we identified for the first time that the low expression and hypermethylation of ACSM5 gene in tissues and cells of HLF, and overexpression of ACSM5 restrained proliferation, anti-apoptosis, and fibrosis of HLF cells in vitro. DNMT1, a type of DNA methyltransferases (DNMTs), can maintain and catalyze abnormal DNA hypermethylation of CpG islands to cause the loss expression of genes specific, thereby regulating the occurrence and development of various fibrosis-related diseases [46,47,48]. Our data demonstrated that the low expression of ACSM5 in HLF cells was negatively regulated by DNMT1, and suppression of DNMT1 restrained the proliferation, anti-apoptosis, and fibrosis of HLF cells in vitro, which can be reversed by silencing ACSM5. Taken together, our findings not only shed a new light on the mechanism of HLF, but also revealed that interference of DNMT1/ACSM5 signaling pathway might be a potential preventive and therapeutic strategy against patients with HLF.
However, there are certain limitations in this study. Firstly, the role of DNMT1/ACSM5 signaling pathway in the development of HLF has not been further verified in animal model with HLF. Secondly, the downstream signaling of ACSM5 in HLF remained largely unknown, which would be the focus of our future research.
The present study revealed that positive correlations between the methylation level of COL4A3, HSPA6, PIGC, and ACSM5 and the development of HLF by the use of genome-wide DNA methylation and SNP analysis. Most importantly, we proved that downregulation of ACSM5 by DNMT1 could inhibit the proliferation and fibrosis, and promote the apoptosis of LFCs from HLF patients, shedding light on the therapeutic potential of DNMT1/ACSM5 signaling pathway.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ligamentum flavum hypertrophy
Lumbar degeneration disease
Lumbar spinal canal stenosis
Quantitative methylation-specific analysis
Kyoto Encyclopedia of Genes and Genomes
Frost BA, Camarero-Espinosa S, Foster EJ. Materials for the spine: anatomy, problems, and solutions. Materials (Basel). 2019;12(2). https://doi.org/10.3390/ma12020253.
Munns JJ, Lee JYB, Espinoza Orías AA, Takatori R, Andersson GBJ, An HS, et al. Ligamentum flavum hypertrophy in asymptomatic and chronic low back pain subjects. PLoS One. 2015;10(5):e0128321. https://doi.org/10.1371/journal.pone.0128321.
Sakai Y, Ito S, Hida T, Ito K, Harada A, Watanabe K. Clinical outcome of lumbar spinal stenosis based on new classification according to hypertrophied ligamentum flavum. J Orthop Sci. 2017;22(1):27–33. https://doi.org/10.1016/j.jos.2016.08.007.
Botwin KP, Gruber RD. Lumbar spinal stenosis: anatomy and pathogenesis. Phys Med Rehabil Clin N Am. 2003;14(1):–15. https://doi.org/10.1016/S1047-9651(02)00063-3.
Viejo-Fuertes D, Liguoro D, Rivel J, Midy D, Guerin J. Morphologic and histologic study of the ligamentum flavum in the thoraco-lumbar region. Surg Radiol Anat. 1998;20(3):171–6. https://doi.org/10.1007/BF01628891.
Schroeder GD, Kurd MF, Vaccaro AR. Lumbar spinal stenosis: how is it classified? J Am Acad Orthop Surg. 2016;24(12):843–52. https://doi.org/10.5435/JAAOS-D-15-00034.
Adams MA, Roughley PJ. What is intervertebral disc degeneration, and what causes it? Spine. 2006;31(18):2151–61. https://doi.org/10.1097/01.brs.0000231761.73859.2c.
Zhong Z-M, Zha D-S, Xiao W-D, Wu S-H, Wu Q, Zhang Y, et al. Hypertrophy of ligamentum flavum in lumbar spine stenosis associated with the increased expression of connective tissue growth factor. J Orthop Res. 2011;29(10):1592–7. https://doi.org/10.1002/jor.21431.
Sairyo K, Biyani A, Goel VK, Leaman DW, Booth R, Thomas J, et al. Lumbar ligamentum flavum hypertrophy is due to accumulation of inflammation-related scar tissue. Spine. 2007;32(11):E340–7. https://doi.org/10.1097/01.brs.0000263407.25009.6e.
Chao Y-H, Tsuang Y-H, Sun J-S, Sun M-G, Chen M-H. Centrifugal force induces human ligamentum flavum fibroblasts inflammation through activation of JNK and p38 pathways. Connect Tissue Res. 2012;53(5):422–9. https://doi.org/10.3109/03008207.2012.685132.
Hur JW, Kim B-J, Park J-H, Kim J-H, Park Y-K, Kwon T-H, et al. The mechanism of ligamentum flavum hypertrophy: introducing angiogenesis as a critical link that couples mechanical stress and hypertrophy. Neurosurgery. 2015;77(2):–282. https://doi.org/10.1227/NEU.0000000000000755.
Nakatani T, Marui T, Hitora T, Doita M, Nishida K, Kurosaka M. Mechanical stretching force promotes collagen synthesis by cultured cells from human ligamentum flavum via transforming growth factor-beta1. J Orthop Res. 2002;20(6):1380–6. https://doi.org/10.1016/S0736-0266(02)00046-3.
Tyack PL, Calambokidis J, Friedlaender A, Goldbogen J, Southall B, Formal Comment on, et al. First long-term behavioral records from Cuvier’s beaked whales (Ziphius cavirostris) reveal record-breaking dives. PLoS One. 2014;9(3):e92633. https://doi.org/10.1371/journal.pone.0092633 PLoS ONE 2015, 10(12):e0142287.
Soll JM, Sobol RW, Mosammaparast N. Regulation of DNA alkylation damage repair: lessons and therapeutic opportunities. Trends Biochem Sci. 2017;42(3):206–18. https://doi.org/10.1016/j.tibs.2016.10.001.
Choi YJ, Ohn JH, Kim N, Kim W, Park K, Won S, et al. Family-based exome sequencing combined with linkage analyses identifies rare susceptibility variants of MUC4 for gastric cancer. PLoS One. 2020;15(7):e0236197. https://doi.org/10.1371/journal.pone.0236197.
Biancari F, Dalén M, Ruggieri VG, Demal T, Gatti G, Onorati F, et al. Prognostic impact of multiple prior percutaneous coronary interventions in patients undergoing coronary artery bypass grafting. J Am Heart Assoc. 2018;7(20):e010089. https://doi.org/10.1161/JAHA.118.010089.
Ng PC, Henikoff S. SIFT: Predicting amino acid changes that affect protein function. Nucleic Acids Res. 2003;31(13):3812–4. https://doi.org/10.1093/nar/gkg509.
Ramensky V, Bork P, Sunyaev S. Human non-synonymous SNPs: server and survey. Nucleic Acids Res. 2002;30(17):3894–900. https://doi.org/10.1093/nar/gkf493.
Souren NY, Gerdes LA, Lutsik P, Gasparoni G, Beltrán E, Salhab A, et al. DNA methylation signatures of monozygotic twins clinically discordant for multiple sclerosis. Nat Commun. 2019;10(1):2094. https://doi.org/10.1038/s41467-019-09984-3.
Yu GC, Wang LG, Han YY. He QY: clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7. https://doi.org/10.1089/omi.2011.0118.
Ramalho-Carvalho J, Henrique R, Jerónimo C. Methylation-specific PCR. Methods Mol Biol. 1708;2018:447–72. https://doi.org/10.1007/978-1-4939-7481-8_23.
Panagopoulou M, Karaglani M, Balgkouranidou I, Biziota E, Koukaki T, Karamitrousis E, et al. Circulating cell-free DNA in breast cancer: size profiling, levels, and methylation patterns lead to prognostic and predictive classifiers. Oncogene. 2019;38(18):3387–401. https://doi.org/10.1038/s41388-018-0660-y.
Mahmood T, Yang P-C. Western blot: technique, theory, and trouble shooting. N Am J Med Sci. 2012;4(9):429–34. https://doi.org/10.4103/1947-2714.100998.
Unnikrishnan A, Freeman WM, Jackson J, Wren JD, Porter H, Richardson A. The role of DNA methylation in epigenetics of aging. Pharmacol Ther. 2019;195:172–85. https://doi.org/10.1016/j.pharmthera.2018.11.001.
Mehdipour P, Murphy T, De Carvalho DD. The role of DNA-demethylating agents in cancer therapy. Pharmacol Ther. 2020;205:107416. https://doi.org/10.1016/j.pharmthera.2019.107416.
Wessberg P, Frennered K. Central lumbar spinal stenosis: natural history of non-surgical patients. European spine journal : official publication of the European Spine Society, the European Spinal Deformity Society, and the European Section of the Cervical Spine Research Society. 2017;26(10):2536–42. https://doi.org/10.1007/s00586-017-5075-x.
Sun C, Zhang H, Wang X, Liu X. Ligamentum flavum fibrosis and hypertrophy: molecular pathways, cellular mechanisms, and future directions. FASEB J. 2020;34(8):9854–68. https://doi.org/10.1096/fj.202000635R.
Kawaguchi Y. Genetic background of degenerative disc disease in the lumbar spine. Spine Surg Relat Res. 2018;2(2):–112. https://doi.org/10.22603/ssrr.2017-0007.
Chiba N, Furukawa K, Takayama S, Asari T, Chin S, Harada Y, et al. Decreased DNA methylation in the promoter region of the WNT5A and GDNF genes may promote the osteogenicity of mesenchymal stem cells from patients with ossified spinal ligaments. J Pharmacol Sci. 2015;127(4):467–73. https://doi.org/10.1016/j.jphs.2015.03.008.
Fan T, Meng X, Sun C, Yang X, Chen G, Li W, et al. Genome-wide DNA methylation profile analysis in thoracic ossification of the ligamentum flavum. J Cell Mol Med. 2020;24(15):8753–62. https://doi.org/10.1111/jcmm.15509.
Cao YL, Duan Y, Zhu LX, Zhan YN, Min SX, Jin AM. TGF-β1, in association with the increased expression of connective tissue growth factor, induce the hypertrophy of the ligamentum flavum through the p38 MAPK pathway. Int J Mol Med. 2016;38(2):391–8. https://doi.org/10.3892/ijmm.2016.2631.
Zheng Z, Ao X, Li P, Lian Z, Jiang T, Zhang Z, et al. CRLF1 is a key regulator in the ligamentum flavum hypertrophy. Front Cell Dev Biol. 2020;8:858. https://doi.org/10.3389/fcell.2020.00858.
Mori T, Sakai Y, Kayano M, Matsuda A, Oboki K, Matsumoto K, et al. MicroRNA transcriptome analysis on hypertrophy of ligamentum flavum in patients with lumbar spinal stenosis. Spine Surg Related Res. 2017;1(4):211–7. https://doi.org/10.22603/ssrr.1.2017-0023.
Chen MH, Hu CK, Chen PR, Chen YS, Sun JS, Chen MH. Dose-dependent regulation of cell proliferation and collagen degradation by estradiol on ligamentum flavum. BMC Musculoskelet Disord. 2014;15(1):238. https://doi.org/10.1186/1471-2474-15-238.
Zhou T, Du L, Chen C, Han C, Li X, Qin A, et al. Lysophosphatidic acid induces ligamentum flavum hypertrophy through the LPAR1/Akt pathway. Cell Physiol Biochem. 2018;45(4):1472–86. https://doi.org/10.1159/000487574.
Wang Y, Zhang J, Zhao Y, Wang S, Zhang J, Han Q, et al. COL4A3 gene variants and diabetic kidney disease in MODY. Clin J Am Soc Nephrol. 2018;13(8):1162–71. https://doi.org/10.2215/CJN.09100817.
Malone AF, Phelan PJ, Hall G, Cetincelik U, Homstad A, Alonso AS, et al. Rare hereditary COL4A3/COL4A4 variants may be mistaken for familial focal segmental glomerulosclerosis. Kidney Int. 2014;86(6):1253–9. https://doi.org/10.1038/ki.2014.305.
Storey H, Savige J, Sivakumar V, Abbs S, Flinter FA. COL4A3/COL4A4 mutations and features in individuals with autosomal recessive Alport syndrome. J Am Soc Nephrol. 2013;24(12):1945–54. https://doi.org/10.1681/ASN.2012100985.
Ozdemir G, Gulhan B, Atayar E, Saygılı S, Soylemezoglu O, Ozcakar ZB, et al. COL4A3 mutation is an independent risk factor for poor prognosis in children with Alport syndrome. Pediatr Nephrol (Berlin, Germany). 2020;35(10):1941–52.
van der Sluis R, Erasmus E. Xenobiotic/medium chain fatty acid: CoA ligase - a critical review on its role in fatty acid metabolism and the detoxification of benzoic acid and aspirin. Expert Opin Drug Metab Toxicol. 2016;12(10):1169–79. https://doi.org/10.1080/17425255.2016.1206888.
Behnaz Y, Mahnoosh J, Negin J, Marjan K, Samira R, Mansoureh A, et al. Altered expression level of ACSM5 in breast cancer: an integrative analysis of tissue biomarkers with diagnostic potential. Gene Rep. 2021;22:100992. https://doi.org/10.1016/j.genrep.2020.100992.
Ma C, Luo H, Cao J, Zheng X, Zhang J, Zhang Y, et al. Identification of a novel tumor microenvironment-associated eight-gene signature for prognosis prediction in lung adenocarcinoma. Front Mol Biosci. 2020;7:571641. https://doi.org/10.3389/fmolb.2020.571641.
Qi X, Wang S, Zhang L, Liu L, Wen Y, Ma M, et al. An integrative analysis of transcriptome-wide association study and mRNA expression profile identified candidate genes for attention-deficit/hyperactivity disorder. Psychiatry Res. 2019;282:112639. https://doi.org/10.1016/j.psychres.2019.112639.
Gao JR, Qin XJ, Jiang H, Wang T, Song JM, Xu SZ. The effects of Qi Teng Xiao Zhuo granules, traditional Chinese medicine, on the expression of genes in chronic glomerulonephritis rats. J Ethnopharmacol. 2016;193:140–9. https://doi.org/10.1016/j.jep.2016.08.011.
Koseler A, Ma F, Kilic ID, Morselli M, Kilic O, Pellegrini M. Genome-wide DNA methylation profiling of blood from monozygotic twins discordant for myocardial infarction. In Vivo (Athens, Greece). 2020;34(1):361–7.
Zhu H, He C, Zhao H, Jiang W, Xu S, Li J, et al. Sennoside A prevents liver fibrosis by binding DNMT1 and suppressing DNMT1-mediated PTEN hypermethylation in HSC activation and proliferation. FASEB J. 2020;34(11):14558–71. https://doi.org/10.1096/fj.202000494RR.
Pan XY, You HM, Wang L, Bi YH, Yang Y, Meng HW, et al. Methylation of RCAN1.4 mediated by DNMT1 and DNMT3b enhances hepatic stellate cell activation and liver fibrogenesis through Calcineurin/NFAT3 signaling. Theranostics. 2019;9(15):4308–23. https://doi.org/10.7150/thno.32710.
Tian L, Wu D, Dasgupta A, Chen KH, Mewburn J, Potus F, et al. Epigenetic metabolic reprogramming of right ventricular fibroblasts in pulmonary arterial hypertension: a pyruvate dehydrogenase kinase-dependent shift in mitochondrial metabolism promotes right ventricular fibrosis. Circ Res. 2020;126(12):1723–45. https://doi.org/10.1161/CIRCRESAHA.120.316443.
This work was supported by grants from the National Natural Science Foundation of China (81470121).
Ethics approval and consent to participate
This study was approved by Zhujiang Hospital, Southern Medical University. All subjects gave their informed written consent by signing a document that had been carefully reviewed by Zhujiang Hospital, Southern Medical University.
Consent for publication
All authors have read the final manuscript and approved the manuscript for publication.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Volcano map showing the up-methylated sites and the down-methylated sites. The red dots indicate the methylation upregulation and the blue dots indicate the methylation downregulation in ligamentum flavum hypertrophy.
GO and KEGG enrichment analysis of the differentially methylated sites. Top 15 bio-functions or pathways enriched of the differentially methylated sites in GO_BP, GO_CC, GO_MF, and KEGG were shown, respectively.
qPCR and Western Blot results showing the expression level of ACSM5 in LFCs from HLF patients with the treatment of pCMV-DNMT1 or empty plasmid. Data were presented as the mean ± SME of three independent experiments. Two-tailed paired t-test, #P <0.05, ##P <0.01.
Identification of the phenotype of cultured ligamentum flavum (LF) cells using immunofluorescence staining. Immunofluorescence staining for collagen III and fibronectin in cultured cells, scale bar=100 μm, n=3 experiments.
About this article
Cite this article
Cao, Y., Zhan, Y., Qiu, S. et al. Integrative analysis of genome-wide DNA methylation and single-nucleotide polymorphism identified ACSM5 as a suppressor of lumbar ligamentum flavum hypertrophy. Arthritis Res Ther 23, 251 (2021). https://doi.org/10.1186/s13075-021-02625-5
- Genome-wide DNA methylation
- Single-nucleotide polymorphism
- Hypertrophy of ligamentum flavum
- Cell fibrosis