- Research article
- Open Access
- Open Peer Review
Abnormal PITX1 gene methylation in adolescent idiopathic scoliosis: a pilot study
BMC Musculoskeletal Disorders volume 19, Article number: 138 (2018)
The gene of pituitary homeobox 1 (PITX1) has been reported to be down-regulated in adolescent idiopathic scoliosis (AIS), of which the cause has not been well addressed. The abnormal DNA methylation was recently assumed to be an important mechanism for the down-regulated genes expression. However, the association between PITX1 promoter methylation and the etiology of AIS was not clear.
The peripheral blood samples of 50 AIS patients and 50 healthy controls were collected and the genomic DNA was extracted. The pyrosequencing assay was used to assess the methylation status of PITX1 promoter and real-time quantitative polymerase chain reaction (PCR) was used to detect the PITX1 gene expression. Comparison analysis was performed using independent t test and Chi-square tests, while correlation analysis were performed with 2-tailed Pearson coefficients.
The mean methylation level was (3.52 ± 0.96)% in AIS and (1.40 ± 0.81)% in healthy controls (P < 0.0001). The PITX1 gene expression was 0.15 ± 0.08 in AIS and 0.80 ± 0.55 in healthy controls (P < 0.0001). The comparative analysis showed significant difference in age (P = 0.021) and Cobb angle of the main curve (P = 0.0001) between AIS groups with positive and negative methylation. The methylation level of 6 CpG sites in PITX1 promoters was significantly associated with Cobb angle of the main curve (P < 0.001) in AIS. No statistical relationship between PITX1 promoter methylation and gene expression was found in AIS (P = 0.842).
Significantly higher methylation level and lower PITX1 gene expression are found in AIS patients. PITX1 methylation is associated with Cobb angles of the main curves in AIS. DNA methylation thus plays an important role in the etiology and curve progression in AIS.
Adolescent idiopathic scoliosis (AIS) is a three-dimensional spinal deformity commonly seen in females with prevalence ranging from 1 to 3% in adolescents [1, 2]. The origin and cause of AIS remain obscure although there are several proposed etiological hypotheses such as genetic variations [2, 3], abnormal bone mineral density [4, 5], asymmetrical growth and biomechanical conditions [6, 7], abnormal changes of the nervous system [8, 9], and hormonal variations [10, 11]. In addition, Fendri et al. [12, 13] comparatively analyzed the gene expression patterns of AIS osteoblasts and healthy osteoblasts by microarray analysis, in which the gene of pituitary homeobox 1 (PITX1) was reported to be down-regulated more than 1.5-fold in AIS. This is an interesting observation of which may potentially shed lights on the etiopathogenesis of AIS.
Even though the significantly lower expression of PITX1 gene in AIS has been previously reported , the cause of this abnormality remains unclear. Recently, the DNA methylation has been proposed as an important epigenetic mechanism operating at the interface between genome and environment to regulate phenotypic plasticity with complex regulations across the genome during the first decade of life . DNA methylation may be the most suitable epigenetic mark for large-scale epidemiological studies since methyl groups are covalently bounded to CpG dinucleotides and are not lost during routine DNA extraction. This opens the possibility of exploiting existing DNA biomarkers in order to discover epigenetic risk factors for complex disease . It was hypothesized that the abnormal DNA methylation might play an important role in the down-regulation of PITX1 gene expression, serving as an epigenetic factor leading to the occurrence and development of AIS.
The purposes of the current study were: 1) to evaluate the DNA methylation of PITX1 promoter in AIS; 2) to investigate the relationship between abnormal DNA methylation and PITX1 gene expression; and 3) further to analyze the role of abnormal DNA methylation in the occurrence and development of AIS through genetic and epigenetic perspectives.
TIANamp Blood DNA kit (Tiangen),DNA isolation kit (Qiagen AB, Solna, Sweden), DNA Methylation Kit™ (Zymo Research, HiSS Diagnostics, CA, USA), ZymoTaq™ Pre mix (Zymo Research, HiSS Diagnostics, CA, USA), PyroMark® Q24 pyrosequencer (Qiagen, Germany), RNeasy Plus Mini kit (QIAGEN), ReverTra Ace qPCR RT kit (Toyobo), SYBR-Green PCR master mix (Toyobo).
New AIS patients seen in our scoliosis center from January 2014 to January 2015 were included in the study. The inclusion criteria were: 1) patients with standing whole spine x-rays showing single thoracic or thoracolumbar curve (Cobb angle ranging from 10 deg. to 50 deg.); 2) age between 10 to 17 years; 3) with MRI of whole spine showing no additional abnormality of spinal cord or vertebrae; and 4) with complete demographic data including sex, age, etc. The exclusion criteria were as follows: 1) patients with previous spinal surgery; and 2) with any signs of growth abnormalities, neurological abnormality, skeletal dysplasia, dwarfism, or recognized genetic syndrome, for example Marfan syndrome. Finally, a total of 50 AIS patients were included, and 50 sex matched and age matched healthy controls were recruited in the study. Considering the association between CpGs methylation of PITX1 gene and cancer [16, 17], the possibility of cancer in both AIS and control groups were excluded with physical examination and inquiry of family history. The peripheral blood samples were collected from all the subjects. This retrospective study was approved by the ethics review board of the affiliated Drum Tower Hospital of Nanjing University, and the methods were carried out in accordance with the approved guidelines. All participating subjects or their legal guardians signed the written informed consents.
DNA extraction and sodium bisulfite modification
Genomic DNA was extracted from blood samples using TIANamp Blood DNA kit (Tiangen). Bisulfite modification was done with 500 ng samples of DNA using DNA Methylation Kit™ (Zymo Research, HiSS Diagnostics, CA, USA). The samples were then eluted in 50 μL of elution buffer and stored at − 20°C.
Generation of PCR products and pyrosequencing
Pyrosequencing  was a sequence-by-synthesis technique based on the fact that DNA polymerization could be monitored by measuring pyrophosphate production, which was detected by light. The pyrosequencing assay was used to assess the methylation status of PITX1 promoter (chr5:135032889–135,034,242) (Fig. 1). This assay detected the level of methylation in the region Chr5:135033586–135,033,643 of the PITX1 gene (chr5:134,363,424–134,369,964). The polymerase chain reaction (PCR) and pyrosequencing primers were designed using PyroMark Assay Design 2.0. The PCR volume was 20 μL, and was incorporated with 0.5 μM forward and reverse primers respectively, 10 μl ZymoTaqTM Pre mix (Zymo Research, HiSS Diagnostics, CA, USA) and 2 μl bisulfite-modified DNA. The forward and reverse primer sequences were 5′- GAT TAA GAG GGG TTG TTA GTT TAA TT -3′ and 5′-biotin- CCC CAA AAA CCT CAA AAA CTT TCT TTT TC -3′, respectively. PCR testing was carried out at 95°C for 10 min, followed by 40 cycles of 95°C for 30s, 55°C for 30s, and 72°C for 1 min, and a final extension at 72°C for 10 min. PCR product quality was confirmed by electrophoresis with 1% agarose gels with ethidium bromide staining. The pyrosequencing primer was 5′- AGT TGT AGG TTA GTG TAT TTG TTA -3′. Pyrosequencing was subsequently carried out in the PyroMark® Q24 pyrosequencer (Qiagen, Germany) (Fig. 2). The degree of methylation of all 6 CpG sites were automatically analyzed by the PyroMark® Q24 software. The exact frequency of methylation could be determined by the resulting peak patterns, which differed between methylated and unmethylated samples .
Total RNA was extracted from blood samples using RNeasy Plus Mini kit (QIAGEN). A total of 1 μg RNA was reversely transcribed using ReverTra Ace qPCR RT kit (Toyobo) according to the manufacturer’s guideline. The resulting cDNA samples were amplified by real-time PCR using primers specific for PITX1 gene expression and GAPDH served as an internal control. Primers were as follows: PITX1: forward, 5′-GTG GCG TAA GCG CGA GCG TAA-3′; Reverse, 5′- GAC AGC GGG CTC ATG GAG TTG AAG-3′. GAPDH (housekeeping gene): forward, 5’-CGG ATT TGG TCG TAT TGG G-3′ and reverse, 5’-CTG GAA GAT GGT GAT GGG ATT-3′. SYBR-Green PCR master mix (Toyobo) containing 333 nM of each forward and reverse primer and 1 μl cDNA were transferred into a Eco plate (Illumina). PCR was performed on Eco™ real-time PCR system (Illumina). The amplification program was as follows: 95°C for 1 min, followed by 40 cycles of denaturation at 95°C for 15 s and annealing at 60°C for 30s. Melting curve was used to monitor specificity of PCR products, which was further quantitated by 2-△△Ct methods.
Clinical and biological variables were compared using independent t test and Chi-square tests. Correlation analysis were performed using the 2-tailed Pearson coefficients of correlation. Statistical significance was set at a level of p < 0.05. All statistical analysis was performed with SPSS Statistics 17.0 software (SPSS, Inc., Somers, NY, USA).
Comparison between healthy controls and AIS
A total of 50 AIS patients and 50 healthy controls were included in the study. Overall, the mean methylation level was (3.52 ± 0.96)% in AIS patients and (1.40 ± 0.81)% in healthy controls (P < 0.0001) (Fig. 3a). Besides, the PITX1 gene expression was 0.80 ± 0.55 in healthy controls and 0.15 ± 0.08 in AIS patients, respectively (P < 0.0001) (Fig. 3b).
Comparison between positive and negative methylation AIS groups
Two subgroups were defined for statistical analysis according to the average of CpG methylation status: the negative methylation group (defined as methylation level < 3.8%) and the positive methylation group (defined as methylation level ≥ 3.8%). A total of 22 AIS patients were divided into positive methylation group while 28 patients were in negative methylation group, respectively (Table 1). The comparative analysis between 2 groups showed no significant difference in sex (P = 0.308). However, significant difference was found in age (P = 0.021) and Cobb angle of main curve (P = 0.0001) between groups. The PITX1 gene expression was 0.07 ± 0.04 and 0.24 ± 0.11 in positive and negative methylation groups, revealing significant difference between 2 groups (P < 0.0001).
Correlation between methylation status and clinical characteristics in AIS
The methylation level of the 6 CpG sites in PITX1 promoters was significantly associated with Cobb angle of main curve (P < 0.001) in AIS patients. No significant correlation was found between methylation status and other clinical parameters including age, height or weight (P > 0.05 for all). Notably, no statistical linear relationship between PITX1 promoter methylation and gene expression was found (P = 0.842).
Several etiological hypotheses including both genetic and environmental factors have recently been proposed for the occurrence of AIS [12, 14]. Being reversible and heritable modifications, epigenetic markers such as DNA methylation are important to reflect the interactions between genetic factors and environmental exposures [20, 21]. We further assumed that the DNA methylation can also be up- or down-regulated in AIS patients, which serves as an important mechanism of this longitudinal and complex spinal deformity. Therefore, this study compared the DNA methylation between AIS patients and normal controls, aiming to explore the etiology of AIS through genetic and epigenetic aspects.
According to previous studies, the PITX1 is a member of the RIEG/PITX homeobox transcription factors and acts as a transcriptional regulator involved in basal and hormone-regulated activity of prolactin . Members of this family are involved in organ development and left-right asymmetry . The PITX1 abnormality is associated with many bone related diseases including congenital clubfoot, with or without deficiency of long bones and/or mirror-image polydactyly, and Liebenberg syndrome [24,25,26]. Besides, the abnormal PITX1 gene expression was also found in other pathological conditions such as lung cancer, cutaneous malignant melanoma and others [27, 28]. Hence, there is potential association between PITX abnormality and AIS. Recently, Fendri et al.  identified several genes involved in various bone regulatory and developmental pathways and showed many of them can be grouped into clusters to participate in a particular biological pathway . Fendri’s study proved that the PITX1 gene expression was significantly lower in AIS patients. However, the cause of reduced expression of PITX1 and its corresponding mechanism in AIS remained unknown.
The current study was conducted to elucidate the relationship between DNA methylation and PITX1 abnormality in AIS. The comparative analysis showed a significantly higher PITX1 methylation and lower gene expression in AIS patients as compared with healthy controls. These results further supported the results reported by Fendri et al.  Most importantly, our study clearly showed that the DNA methylation played an important role in the decreased PITX1 gene expression due to the relatively high promoter region methylation and this may be associated with the etiology of AIS.
To further investigate the specific effect of abnormal PITX1 methylation on the clinical and biological characteristics of AIS patients, the age, sex and Cobb angles were compared between AIS patients with positive and negative PITX1 methylation. The results revealed significant difference in age and Cobb angle of the main curves between the two groups. Bork et al.  reported the in-vivo hypo-methylated changes upon aging of CpGs from PITX1 promoter in mesenchymal stromal cells. The percentages of positive methylation were significantly higher in younger patients (10–13 years) than elder patients (14–16 years) in the current study, which supported the conclusion of Bork et al.  The rates of positive methylation in AIS patients with Cobb angle >30° were significantly higher, which showed that the DNA methylation, especially PITX methylation, might be correlated with the curve progression and curve severity. The significant correlation between methylation level of the 6 CpG sites in PITX1 promoters and Cobb angles further supported the findings. In addition, it seemed AIS boys exhibited more hyper-methylation pattern than girls even though the difference was not statistically significant (P = 0.308). If this phenomenon can be proven by further studies with a larger sample size, the influence of estrogen should not be ignored . The role of estrogen in the regulation of PITX1 expression has been illustrated in other diseases such as breast cancer , however, how estrogen regulates PITX1 expression in AIS has not been reported before.
Combining the above results, we further assumed that the PITX1 promoter methylation and gene expression should be well associated. However, linear correlation analysis indicated no statistically significant correlation between PITX1 promoter methylation and gene expression (P = 0.842). In addition, the results of the current study showed that only 44.0% (22 of 50) patients were identified with positive methylation of PITX1 promoter. We presumed 2 possible explanations for this significant deviation. Firstly, the methylation region of PITX1 gene involved in the pyrosequencing assay was only part of the whole gene, and other fragments might also be responsible for the regulation of gene expression. Secondly, there should be other regulation pathways apart from PITX1 promoter methylation, such as the cartilage oligomeric matrix protein (COMP), protocadherin10 (PCDH10) and others . Hence, further studies focusing on the methylation of other genes in AIS are indeed needed.
There were several limitations with this study. Firstly, the relatively small sample size was the largest limitation of our study. Secondly, only PITX1 promoter methylation was included in the analysis. Since it has been widely accepted that AIS resulted from combination of multiple genes, abnormalities of other regions of PITX1 gene and other potential genes was not assessed. The decision of cut-off point for positive and negative methylated subgroups based on the results of the current study was arbitrary. Nevertheless, this study put forward a novel hypothesis based on the association between genetics and epigenetics to explain the etiology and curve progression in AIS.
Recent data suggested that epigenetic responses including DNA methylation were involved not only in cellular differentiation but also in modulation of genome function in response to signals from various environments . Therefore, our study not only indicated the relationship between abnormal gene methylation and occurrence of AIS but also implied the link between the Cobb angle and the aberrant methylation of PITX1 gene promoter. In addition, the results suggested that PITX1 gene methylation may be a new molecular diagnostic factor and a new biomarker for curve progression.
AIS patients have significantly higher methylation level and lower PITX1 gene expression when compared to normal controls. In AIS patients, PITX1 methylation is significantly associated with the Cobb angles of main curves. This close association between genetics and epigenetics in AIS puts forward a novel hypothesis which suggests DNA methylation may play an important role in the etiology and curve progression of AIS.
Adolescent idiopathic scoliosis
Polymerase chain reaction
Pituitary homeobox 1
Weinstein SL, Dolan LA, Cheng JC, Danielsson A, Morcuende JA. Adolescent idiopathic scoliosis. Lancet. 2008;371(9623):1527–37.
Zhu Z, Tang NL, Xu L, Qin X, Mao S, Song Y, Liu L, Li F, Liu P, Yi L, et al. Genome-wide association study identifies new susceptibility loci for adolescent idiopathic scoliosis in Chinese girls. Nat Commun. 2015;6:8355.
Xu L, Qiu X, Sun X, Mao S, Liu Z, Qiao J, Qiu Y. Potential genetic markers predicting the outcome of brace treatment in patients with adolescent idiopathic scoliosis. Eur Spine J. 2011;20(10):1757–64.
Yu WS, Chan KY, Yu FW, Yeung HY, Ng BK, Lee KM, Lam TP, Cheng JC. Abnormal bone quality versus low bone mineral density in adolescent idiopathic scoliosis: a case-control study with in vivo high-resolution peripheral quantitative computed tomography. Spine J. 2013;13(11):1493–9.
Cheuk KY, Zhu TY, Yu FW, Hung VW, Lee KM, Qin L, Cheng JC, Lam TP. Abnormal bone mechanical and structural properties in adolescent idiopathic scoliosis: a study with finite element analysis and structural model index. Calcif Tissue Int. 2015;97(4):343–52.
Wang S, Qiu Y, Zhu Z, Ma Z, Xia C, Zhu F. Histomorphological study of the spinal growth plates from the convex side and the concave side in adolescent idiopathic scoliosis. J Orthop Surg Res. 2007;2:19.
Huynh AM, Aubin CE, Rajwani T, Bagnall KM, Villemure I. Pedicle growth asymmetry as a cause of adolescent idiopathic scoliosis: a biomechanical study. Eur Spine J. 2007;16(4):523–9.
Deng M, Hui SC, Yu FW, Lam TP, Qiu Y, Ng BK, Cheng JC, Chu WC. MRI-based morphological evidence of spinal cord tethering predicts curve progression in adolescent idiopathic scoliosis. Spine J. 2015;15(6):1391–401.
Kong Y, Shi L, Hui SC, Wang D, Deng M, Chu WC, Cheng JC. Variation in anisotropy and diffusivity along the medulla oblongata and the whole spinal cord in adolescent idiopathic scoliosis: a pilot study using diffusion tensor imaging. AJNR Am J Neuroradiol. 2014;35(8):1621–7.
Nissinen M, Heliovaara M, Seitsamo J, Poussa M. Trunk asymmetry, posture, growth, and risk of scoliosis. A three-year follow-up of Finnish prepubertal school children. Spine (Phila Pa 1976). 1993;18(1):8–13.
Tam EM, Liu Z, Lam TP, Ting T, Cheung G, Ng BK, Lee SK, Qiu Y, Cheng JC. Lower muscle mass and body fat in adolescent idiopathic scoliosis are associated with abnormal leptin bioavailability. Spine. 2015;
Fendri K, Patten SA, Kaufman GN, Zaouter C, Parent S, Grimard G, Edery P, Moldovan F. Microarray expression profiling identifies genes with altered expression in adolescent idiopathic scoliosis. Eur Spine J. 2013;22(6):1300–11.
Fendri K, Moldovan F, Parent S, Zaouter C, Moreau A, Patten KS. Microarray analysis of gene expression in primary human osteoblasts derived from spine. Osteoarthritis Cartilage. 2012;20(4):S197–8.
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.
Talens RP, Boomsma DI, Tobi EW, Kremer D, Jukema JW, Willemsen G, Putter H, Slagboom PE, Heijmans BT. Variation, patterns, and temporal stability of DNA methylation: considerations for epigenetic epidemiology. FASEB J. 2010;24(9):3135–44.
Wei JH, Haddad A, Wu KJ, Zhao HW, Kapur P, Zhang ZL, Zhao LY, Chen ZH, Zhou YY, Zhou JC, et al. A CpG-methylation-based assay to predict survival in clear cell renal cell carcinoma. Nat Commun. 2015;6:8699.
Bell A, Bell D, Weber RS, El-Naggar AK. CpG island methylation profiling in human salivary gland adenoid cystic carcinoma. Cancer. 2011;117(13):2898–909.
Harrington CT, Lin EI, Olson MT, Eshleman JR. Fundamentals of pyrosequencing. Arch Pathol Lab Med. 2013;137(9):1296–303.
Zainuddin N, Kanduri M, Berglund M, Lindell M, Amini RM, Roos G, Sundstrom C, Enblad G, Rosenquist R. Quantitative evaluation of p16(INK4a) promoter methylation using pyrosequencing in de novo diffuse large B-cell lymphoma. Leuk Res. 2011;35(4):438–43.
Wong CC, Caspi A, Williams B, Craig IW, Houts R, Ambler A, Moffitt TE, Mill J. A longitudinal study of epigenetic variation in twins. Epigenetics. 2010;5(6):516–26.
Cheng J, Wang Y, Zhou K, Wang L, Li J, Zhuang Q, Xu X, Xu L, Zhang K, Dai D, et al. Male-specific association between dopamine receptor D4 gene methylation and schizophrenia. PLoS One. 2014;9(2):e89128.
Logan M, Tabin CJ. Role of Pitx1 upstream of Tbx4 in specification of hindlimb identity. Science (New York, NY). 1999;283(5408):1736–9.
Gurnett CA, Alaee F, Kruse LM, Desruisseau DM, Hecht JT, Wise CA, Bowcock AM, Dobbs MB. Asymmetric lower-limb malformations in individuals with homeobox PITX1 gene mutation. Am J Hum Genet. 2008;83(5):616–22.
Alvarado DM, McCall K, Aferol H, Silva MJ, Garbow JR, Spees WM, Patel T, Siegel M, Dobbs MB, Gurnett CA. Pitx1 haploinsufficiency causes clubfoot in humans and a clubfoot-like phenotype in mice. Hum Mol Genet. 2011;20(20):3943–52.
Al-Qattan MM, Al-Thunayan A, Alabdulkareem I, Al Balwi M. Liebenberg syndrome is caused by a deletion upstream to the PITX1 gene resulting in transformation of the upper limbs to reflect lower limb characteristics. Gene. 2013;524(1):65–71.
Pandey SN, Cabotage J, Shi R, Dixit M, Sutherland M, Liu J, Muger S, Harper SQ, Nagaraju K, Chen YW. Conditional over-expression of PITX1 causes skeletal muscle dystrophy in mice. Biology Open. 2012;1(7):629–39.
Chen Y, Knosel T, Ye F, Pacyna-Gengelbach M, Deutschmann N, Petersen I. Decreased PITX1 homeobox gene expression in human lung cancer. Lung Cancer (Amsterdam, Netherlands). 2007;55(3):287–94.
Osaki M, Chinen H, Yoshida Y, Ohhira T, Sunamura N, Yamamoto O, Ito H, Oshimura M, Kugoh H. Decreased PITX1 gene expression in human cutaneous malignant melanoma and its clinicopathological significance. Eur Spine J. 2013;23(3):344–9.
Bork S, Pfister S, Witt H, Horn P, Korn B, Ho AD, Wagner W. DNA methylation pattern changes upon long-term culture and aging of human mesenchymal stromal cells. Aging Cell. 2010;9(1):54–63.
Leboeuf D, Letellier K, Alos N, Edery P, Moldovan F. Do estrogens impact adolescent idiopathic scoliosis? Trends Endocrinol Metab. 2009;20(4):147–52.
Stender JDSF, Funk CC, Charn TH, Barnett DH, Katzenellenbogen BS. The estrogen-regulated transcription factor PITX1 coordinates gene-specific regulation by estrogen receptor-alpha in breast cancer cells. Mol Endocrinol. 2011;25(10):1699–709.
Szyf M. The early life social environment and DNA methylation: DNA methylation mediating the long-term impact of social environments early in life. Epigenetics. 2011;6(8):971–8.
Thanks are due to Prof. Tsz-Ping Lam from Department of Orthopaedics and Traumatology, Chinese University of Hong Kong, for assistance with the revision of the manuscript.
This work was financially supported by the China Postdoctoral Science Foundation (2017 M610323) and the Natural Science Foundation of Jiangsu Province (BK20170126).
Availability of data and materials
The data and materials in current paper may be made available upon request through sending an e-mail to first author.
Ethics approval and consent to participate
The Drum Tower Hospital’s review board approved this retrospective study and all participating subjects or their legal guardians signed the written informed consents..
ZZ is a member of the Editorial Board of BMC Musculoskeletal Disorders. The other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Shi, B., Xu, L., Mao, S. et al. Abnormal PITX1 gene methylation in adolescent idiopathic scoliosis: a pilot study. BMC Musculoskelet Disord 19, 138 (2018) doi:10.1186/s12891-018-2054-2
- Adolescent idiopathic scoliosis
- Pituitary homeobox 1
- Gene methylation